Electron and nuclear dynamics of molecular clusters in ultraintense lase fields.III.Coulomb explosion of deuterium clusters

摘要

In this paper we present a theoretical and computational study of the energetics and temporal dynamics of Coulomb explosion of molecular clusters of deuterium (D_2)_(n/2) n=480-7.6X10~4,cluster radius R_0=13.1-70A) in ultraintense laser fields (laser peak intensity I=10~(15)-10~(20) W cm~~2).The energetics of Coulomb explosion was inferred from the dependence of the maximal energy E_M and the average energy E_(av_ of the product D~+ ions on the laser intensity,the laser pulse shape,the cluster radius,and the laser frequency.Electron dynamics of outer cluster ionization and nuclear dynamics of Coulomb explosion were investigated by molecular dynamics simulations.Several distinct laser pulse shape envelopes,involving a rectangular field,a Gaussian field,and a truncated Gaussian field,were employed to determine the validity range of the cluster vertical ionization (CVI) approximation.The CVI predicts that E_M~(implied by)R_0~2 and that the energy distribution is P(E) implied by E~(1/2).For a rectangular laser pulse the CVT conditions are satisfied when complete outer ionization is obtained,with the outer ionization time t_(oi) being shorter than both the pulse width and the cluster radius doubling time tau_2.By increasing t_(oi),due to the increase of R_0 or the decrease of I,we have shown that the deviation of E_(av) from the corresponding CVI value (E_(av)~(CVI) ) is (E_(av)~(CVI) -E_(av))/E_(av)~(CVI) approx= (t_(oi)/2.9l tau_2)~2.The Gaussian pulses trigger outer ionization induced by adiabatic following of the laser field and of the cluster size,providing a pseudo-CVI behavior at sufficiently large laser fields.The energetics manifest the existence of a finite range of CVI size dependence,with the validity range for the applicability of the CVI being R_0<=(R_0)_I,with (R_0)_I representing an intensity dependent boundary radius.Relating electron dynamics of outer ionization to nuclear dynamics for Coulomb explosion induced by a Gaussian pulse,the boundary radius (R_0)_I and the corresponding ion average energy (E_(av))_I were inferred from simulations and described in terms of an electrostatic model.Two independent estimates of (R_0)_I which involve the cluster size where the CVI relation breaks down and the cluster size for the attainment of complete outer ionization,are in good agreement with each other,as well as with the electrostatic model for cluster barrier suppression.The relation (E_(av))_I implied by (R_0)~2_I provides the validity range of the pseudo-CVI domain for the cluster sizes and laser intensities,where the energetics of D~+ ions produced by Coulomb explosion of (D)_n clusters is optimized.The currently available experimental data [Madison et al.,Phys.Plasmas 11,1 (2004)] for the energetics of Coulomb explosion of (D)_n clusters (E_(av)=5-7keV at I=2X10~(18)Wcm~(-2)),together with our simulation data,lead to the estimates of R_0=51-60 A,which exceed the experimental estimate of R_0=45 A.The predicted anisotropy of the D~+ ion energies in the Coulomb explosion at I=10~(18) Wcm~(-2) is in accord with experiment.We also explored the laser frequency dependence of the energetics of Coulomb explosion in the range v=0.1 -2.1 fs~(-1) (gamma=3000-140 nm),which can be rationalized in terms of the electrostatic model.