Dear All,
I tried to simulate the ablation of the laser with the fix ehex command, but I see in the literature that it takes the laser energy density as a variable Like 3 × 10^9 W/cm2. I wonder how to calculate energy density.Here is the code I wrote about ablation, how should I calculate its laser energy density?
compute 1 all ave/sphere/atom
dump 1 all custom 100 dump.xyz id type x y z c_1[2]
fix 1 all nve
variable dy equal -15+stepdt2
region hot sphere -10 v_dy 65 10 units box
fix 2 all heat 1 400 region hot
run 33000
Best regards,
Hui Fu
That is a question that is very specific about the kind of research that you are doing. The chances to find somebody in a general forum like this with the specific expertise (and willing to help you) are very, very low. You should thus look for somebody with the particular expertise and discuss with that person.
From my limited experience, I am surprised that you are using fix ehex and specifically with the groups all and hot. That does not look to me like a good model for your purpose. Most people here that post inputs about materials interacting with laser radiation use either fix ttm or fix ttm/mod.
Thank you for your answer, I will try to use it