Combined molecular dynamics-direct simulation Monte Carlo computational study of laser ablation plume evolution

Michael I. Zeifman, Barbara Jane Garrison, Leonid V. Zhigilei

Research output: Contribution to journalArticle

63 Citations (Scopus)

Abstract

A two-stage computational model of evolution of a plume generated by laser ablation of an organic solid is proposed and developed. The first stage of the laser ablation, which involves laser coupling to the target and ejection of molecules and clusters, is described by the molecular dynamics (MD) method. The second stage of a long-term expansion of the ejected plume is modeled by the direct simulation Monte Carlo (DSMC) method. The presence of clusters, which comprise a major part of the overall plume at laser fluences above the ablation threshold, presents the main computational challenge in the development of the combined model. An extremely low proportion of large-sized clusters hinders both the statistical estimation of their characteristics from the results of the MD model and the following representation of each cluster size as a separate species, as required in the conventional DSMC. A number of analytical models are proposed and verified for the statistical distributions of translational and internal energies of monomers and clusters as well as for the distribution of the cluster sizes, required for the information transfer from the MD to the DSMC parts of the model. The developed model is applied to simulate the expansion of the ablation plume ejected in the stress-confinement irradiation regime. The presence of the directly ejected clusters drastically changes the evolution of the plume as compared to the desorption regime. A one-dimensional self-similar flow in the direction normal to the ablated surface is developed within the entire plume at the MD stage. A self-similar two-dimensional flow of monomers forms in the major part of the plume by about 40 ns, while its counterpart for large clusters forms much later, leading to the plume sharpening effect. The expansion of the entire plume becomes self-similar by about 500 ns, when interparticle interactions vanish. The velocity distribution of particles cannot be characterized by a single translational temperature; rather, it is characterized by a spatially and direction dependent statistical scatter about the flow velocity. The cluster size dependence of the internal temperature is mainly defined by the size dependence of the unimolecular dissociation energy of a cluster.

Original languageEnglish (US)
Pages (from-to)2181-2193
Number of pages13
JournalJournal of Applied Physics
Volume92
Issue number4
DOIs
StatePublished - Aug 15 2002

Fingerprint

laser ablation
plumes
molecular dynamics
simulation
ablation
expansion
monomers
information transfer
organic solids
two dimensional flow
statistical distributions
internal energy
ejection
dynamic models
lasers
Monte Carlo method
proportion
fluence
flow velocity
velocity distribution

All Science Journal Classification (ASJC) codes

  • Physics and Astronomy(all)

Cite this

Zeifman, Michael I. ; Garrison, Barbara Jane ; Zhigilei, Leonid V. / Combined molecular dynamics-direct simulation Monte Carlo computational study of laser ablation plume evolution. In: Journal of Applied Physics. 2002 ; Vol. 92, No. 4. pp. 2181-2193.
@article{86b09733c61348cbbac638afea4259e2,
title = "Combined molecular dynamics-direct simulation Monte Carlo computational study of laser ablation plume evolution",
abstract = "A two-stage computational model of evolution of a plume generated by laser ablation of an organic solid is proposed and developed. The first stage of the laser ablation, which involves laser coupling to the target and ejection of molecules and clusters, is described by the molecular dynamics (MD) method. The second stage of a long-term expansion of the ejected plume is modeled by the direct simulation Monte Carlo (DSMC) method. The presence of clusters, which comprise a major part of the overall plume at laser fluences above the ablation threshold, presents the main computational challenge in the development of the combined model. An extremely low proportion of large-sized clusters hinders both the statistical estimation of their characteristics from the results of the MD model and the following representation of each cluster size as a separate species, as required in the conventional DSMC. A number of analytical models are proposed and verified for the statistical distributions of translational and internal energies of monomers and clusters as well as for the distribution of the cluster sizes, required for the information transfer from the MD to the DSMC parts of the model. The developed model is applied to simulate the expansion of the ablation plume ejected in the stress-confinement irradiation regime. The presence of the directly ejected clusters drastically changes the evolution of the plume as compared to the desorption regime. A one-dimensional self-similar flow in the direction normal to the ablated surface is developed within the entire plume at the MD stage. A self-similar two-dimensional flow of monomers forms in the major part of the plume by about 40 ns, while its counterpart for large clusters forms much later, leading to the plume sharpening effect. The expansion of the entire plume becomes self-similar by about 500 ns, when interparticle interactions vanish. The velocity distribution of particles cannot be characterized by a single translational temperature; rather, it is characterized by a spatially and direction dependent statistical scatter about the flow velocity. The cluster size dependence of the internal temperature is mainly defined by the size dependence of the unimolecular dissociation energy of a cluster.",
author = "Zeifman, {Michael I.} and Garrison, {Barbara Jane} and Zhigilei, {Leonid V.}",
year = "2002",
month = "8",
day = "15",
doi = "10.1063/1.1494129",
language = "English (US)",
volume = "92",
pages = "2181--2193",
journal = "Journal of Applied Physics",
issn = "0021-8979",
publisher = "American Institute of Physics Publising LLC",
number = "4",

}

Combined molecular dynamics-direct simulation Monte Carlo computational study of laser ablation plume evolution. / Zeifman, Michael I.; Garrison, Barbara Jane; Zhigilei, Leonid V.

In: Journal of Applied Physics, Vol. 92, No. 4, 15.08.2002, p. 2181-2193.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Combined molecular dynamics-direct simulation Monte Carlo computational study of laser ablation plume evolution

AU - Zeifman, Michael I.

AU - Garrison, Barbara Jane

AU - Zhigilei, Leonid V.

PY - 2002/8/15

Y1 - 2002/8/15

N2 - A two-stage computational model of evolution of a plume generated by laser ablation of an organic solid is proposed and developed. The first stage of the laser ablation, which involves laser coupling to the target and ejection of molecules and clusters, is described by the molecular dynamics (MD) method. The second stage of a long-term expansion of the ejected plume is modeled by the direct simulation Monte Carlo (DSMC) method. The presence of clusters, which comprise a major part of the overall plume at laser fluences above the ablation threshold, presents the main computational challenge in the development of the combined model. An extremely low proportion of large-sized clusters hinders both the statistical estimation of their characteristics from the results of the MD model and the following representation of each cluster size as a separate species, as required in the conventional DSMC. A number of analytical models are proposed and verified for the statistical distributions of translational and internal energies of monomers and clusters as well as for the distribution of the cluster sizes, required for the information transfer from the MD to the DSMC parts of the model. The developed model is applied to simulate the expansion of the ablation plume ejected in the stress-confinement irradiation regime. The presence of the directly ejected clusters drastically changes the evolution of the plume as compared to the desorption regime. A one-dimensional self-similar flow in the direction normal to the ablated surface is developed within the entire plume at the MD stage. A self-similar two-dimensional flow of monomers forms in the major part of the plume by about 40 ns, while its counterpart for large clusters forms much later, leading to the plume sharpening effect. The expansion of the entire plume becomes self-similar by about 500 ns, when interparticle interactions vanish. The velocity distribution of particles cannot be characterized by a single translational temperature; rather, it is characterized by a spatially and direction dependent statistical scatter about the flow velocity. The cluster size dependence of the internal temperature is mainly defined by the size dependence of the unimolecular dissociation energy of a cluster.

AB - A two-stage computational model of evolution of a plume generated by laser ablation of an organic solid is proposed and developed. The first stage of the laser ablation, which involves laser coupling to the target and ejection of molecules and clusters, is described by the molecular dynamics (MD) method. The second stage of a long-term expansion of the ejected plume is modeled by the direct simulation Monte Carlo (DSMC) method. The presence of clusters, which comprise a major part of the overall plume at laser fluences above the ablation threshold, presents the main computational challenge in the development of the combined model. An extremely low proportion of large-sized clusters hinders both the statistical estimation of their characteristics from the results of the MD model and the following representation of each cluster size as a separate species, as required in the conventional DSMC. A number of analytical models are proposed and verified for the statistical distributions of translational and internal energies of monomers and clusters as well as for the distribution of the cluster sizes, required for the information transfer from the MD to the DSMC parts of the model. The developed model is applied to simulate the expansion of the ablation plume ejected in the stress-confinement irradiation regime. The presence of the directly ejected clusters drastically changes the evolution of the plume as compared to the desorption regime. A one-dimensional self-similar flow in the direction normal to the ablated surface is developed within the entire plume at the MD stage. A self-similar two-dimensional flow of monomers forms in the major part of the plume by about 40 ns, while its counterpart for large clusters forms much later, leading to the plume sharpening effect. The expansion of the entire plume becomes self-similar by about 500 ns, when interparticle interactions vanish. The velocity distribution of particles cannot be characterized by a single translational temperature; rather, it is characterized by a spatially and direction dependent statistical scatter about the flow velocity. The cluster size dependence of the internal temperature is mainly defined by the size dependence of the unimolecular dissociation energy of a cluster.

UR - http://www.scopus.com/inward/record.url?scp=0037103676&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0037103676&partnerID=8YFLogxK

U2 - 10.1063/1.1494129

DO - 10.1063/1.1494129

M3 - Article

VL - 92

SP - 2181

EP - 2193

JO - Journal of Applied Physics

JF - Journal of Applied Physics

SN - 0021-8979

IS - 4

ER -