We advance a scaling procedure for molecular dynamics (MD) simulations which allows for the computation of extreme ionization levels, multielectron dynamics, and Coulomb explosion (CE) energetics of very large elemental and molecular clusters driven by ultraintense (peak intensities $I_M={10}^{15}–{10}^{21}\mathrm{W}{\mathrm{cm}}^{-2}$) laser fields. The procedure of scaled particle dynamics uses a standard MD code for a scaled cluster with a reduced number of composite ions (nuclei) (heavy pseudoparticles) and of composite electrons (light pseudoparticles). A single scaling parameter $s$ was used for the composition, charge, and mass of all the pseudoparticles, while the initial distances between pseudoparticles and the short-range potential parameters were scaled by $s^{1∕3}$. The scaling procedure was applied for the computation of the ion energies from CE of ${\left({\mathrm{D}}_2\right)}_{n∕2}$ and ${\left(\mathrm{C}{\mathrm{H}}_4\right)}_n$ clusters over the size domain of $n={10}^3–{10}^7$, bridging between clusters and nanodroplets. Deuterons with a maximal energy of $110\mathrm{keV}$ can be effectively produced and utilized for nuclear fusion driven by CE of ${\left(\mathrm{D}\right)}_n$ clusters with $n=1.7\times {10}^6$. The intensity dependence of the ionization levels of ${\left(\mathrm{C}{\mathrm{H}}_4\right)}_n$ clusters $(n=8\times {10}^4–{10}^6)$ manifests ignition and nanoplasma screening effects, establishing the conditions for the production of bare ${\mathrm{C}}^{6+}+{\mathrm{H}}^{+}$ nuclear matter for nucleosynthesis. The scaling procedure will be applicable for large finite systems where the particle motion is mainly governed by long-range (e.g., Coulomb) interactions.

• Received 14 February 2007

DOI:https://doi.org/10.1103/PhysRevA.75.042507