Kinetic energy distribution of multiply charged ions in Coulomb explosion of Xe clusters.

We report on the calculations of kinetic energy distribution (KED) functions of multiply charged, high-energy ions in Coulomb explosion (CE) of an assembly of elemental Xe(n) clusters (average size (n) = 200-2171) driven by ultra-intense, near-infrared, Gaussian laser fields (peak intensities 10(15) - 4 × 10(16) W cm(-2), pulse lengths 65-230 fs). In this cluster size and pulse parameter domain, outer ionization is incomplete∕vertical, incomplete∕nonvertical, or complete∕nonvertical, with CE occurring in the presence of nanoplasma electrons. The KEDs were obtained from double averaging of single-trajectory molecular dynamics simulation ion kinetic energies. The KEDs were doubly averaged over a log-normal cluster size distribution and over the laser intensity distribution of a spatial Gaussian beam, which constitutes either a two-dimensional (2D) or a three-dimensional (3D) profile, with the 3D profile (when the cluster beam radius is larger than the Rayleigh length) usually being experimentally realized. The general features of the doubly averaged KEDs manifest the smearing out of the structure corresponding to the distribution of ion charges, a marked increase of the KEDs at very low energies due to the contribution from the persistent nanoplasma, a distortion of the KEDs and of the average energies toward lower energy values, and the appearance of long low-intensity high-energy tails caused by the admixture of contributions from large clusters by size averaging. The doubly averaged simulation results account reasonably well (within 30%) for the experimental data for the cluster-size dependence of the CE energetics and for its dependence on the laser pulse parameters, as well as for the anisotropy in the angular distribution of the energies of the Xe(q+) ions. Possible applications of this computational study include a control of the ion kinetic energies by the choice of the laser intensity profile (2D∕3D) in the laser-cluster interaction volume.

Time-resolved femtosecond photoelectron spectroscopy by field-induced surface hopping.

We present the extension of our field-induced surface hopping method for the description of the photoionization process and the simulation of time-resolved photoelectron spectra (TRPES). This is based on the combination of nonadiabatic molecular dynamics "on the fly" in the framework of TDDFT generalized for open shell systems under the influence of laser fields with the approximate quantum mechanical description of the photoionization process. Since arbitrary pulse shapes can be employed, this method can be also combined with the optimal control theory in order to steer the photoionization or to shape the outgoing electronic wavepackets. We illustrate our method for the simulation of TRPES on the prototype system of Ag(3), which involves excitation from the equilibrium triangular geometry, as well as excitation from the linear transition state, where in both cases nonadiabatic relaxation takes place in a complex manifold of electronic states. Our approach represents a generally applicable method for the prediction of time-resolved photoelectron spectra and their analysis in systems with complex electronic structure as well as many nuclear degrees freedom. This theoretical development should serve to stimulate new ultrafast experiments.

Extreme dynamics and energetics of Coulomb explosion of Xe clusters.

Unique features of Coulomb explosion (CE) of many-electron elemental Xe(n) (n = 13-2171) clusters driven by ultraintense and ultrashort near-infrared laser pulses (peak intensities 10(15)-10(20) W cm(-2) and pulse lengths of 10-100 fs) manifest ion dynamics and energetics in the extreme, with ultrafast (5-15 A fs(-1)) velocities and ultrahigh (keV-1 MeV) energies. Relations were established between the CE attributes, obtained from molecular dynamics simulations and from electrostatic models, and the extreme cluster inner ionization levels (5-36 per ion), in conjunction with the laser parameters required for the attainment of complete outer ionization, which was approximated by cluster vertical ionization (CVI) initial conditions. Interrelationship between electron dynamics and nuclear dynamics stems from the effects of the laser pulse length on the energetics and from the characterization of the border radius for complete outer ionization. Our computational-theoretical analysis semi-quantitatively describes CE dynamics and energetics in the CVI limit and also the energetics in the presence of a persistent nanoplasma, which is in accord with experiment.

Long-range electronic-to-vibrational energy transfer from nanocrystals to their surrounding matrix environment.

A radiationless transition process of long-range, resonance interconversion of electronic-to-vibrational energy transfer (EVET) is discovered between the band-gap excitation of nanocrystal quantum dots to matrix vibrational overtone modes using fluorescence lifetime measurements. A theoretical analysis based on long-range dipole-dipole nonadiabatic couplings, being distinct from the traditional adiabatic or "static-coupling" pictures, is given and is in qualitative agreement with experiments. EVET should be considered in matrix choices for near-infrared optoelectronic applications of nanocrystals.

Extreme ionization of Xe clusters driven by ultraintense laser fields.

We applied theoretical models and molecular dynamics simulations to explore extreme multielectron ionization in Xe(n) clusters (n=2-2171, initial cluster radius R(0)=2.16-31.0 A) driven by ultraintense infrared Gaussian laser fields (peak intensity I(M)=10(15)-10(20) W cm(-2), temporal pulse length tau=10-100 fs, and frequency nu=0.35 fs(-1)). Cluster compound ionization was described by three processes of inner ionization, nanoplasma formation, and outer ionization. Inner ionization gives rise to high ionization levels (with the formation of [Xe(q+)](n) with q=2-36), which are amenable to experimental observation. The cluster size and laser intensity dependence of the inner ionization levels are induced by a superposition of barrier suppression ionization (BSI) and electron impact ionization (EII). The BSI was induced by a composite field involving the laser field and an inner field of the ions and electrons, which manifests ignition enhancement and screening retardation effects. EII was treated using experimental cross sections, with a proper account of sequential impact ionization. At the highest intensities (I(M)=10(18)-10(20) W cm(-2)) inner ionization is dominated by BSI. At lower intensities (I(M)=10(15)-10(16) W cm(-2)), where the nanoplasma is persistent, the EII contribution to the inner ionization yield is substantial. It increases with increasing the cluster size, exerts a marked effect on the increase of the [Xe(q+)](n) ionization level, is most pronounced in the cluster center, and manifests a marked increase with increasing the pulse length (i.e., becoming the dominant ionization channel (56%) for Xe(2171) at tau=100 fs). The EII yield and the ionization level enhancement decrease with increasing the laser intensity. The pulse length dependence of the EII yield at I(M)=10(15)-10(16) W cm(-2) establishes an ultraintense laser pulse length control mechanism of extreme ionization products.

Tabletop nucleosynthesis driven by cluster Coulomb explosion.

Coulomb explosion of completely ionized (CH4)n, (NH3)n, and (H2O)n clusters will drive tabletop nuclear reactions of protons with 12C6+, 14N7+, and 16O8+ nuclei, extending the realm of nuclear reactions driven by ultraintense laser-heterocluster interaction. The realization for nucleosynthesis in exploding cluster beams requires complete electron stripping from the clusters (at laser intensities I(M) > or = 10(19) W cm(-2)), the utilization of nanodroplets of radius 300-700 A for vertical ionization, and the attainment of the highest energies for the nuclei (i.e., approximately 30 MeV for heavy nuclei and approximately 3 MeV for protons).

Spatial delocalization in para-H2 clusters.

We applied the quantum path integral Monte Carlo method for the study of (para-H)N (N = 5-33) clusters at T = 2 K, exploring static and dynamic order, which originates from the effects of zero-point energy, kinetic energy, and thermal fluctuations in quantum clusters. Information on dynamic structure was inferred from the asymptotic tails of the cage correlation function calculated from the centroid Monte Carlo trajectory. The centroid cage correlation function decays to zero for large clusters (N = 15-33), manifesting the interchange of molecules between different solvation shells, with statistically diminishing back interchange. Further evidence for the floppiness of para-hydrogen clusters emerges from the Monte Carlo evolution of the centroid of a tagged molecule, which exhibits significant changes in the list of its first and second solvation shells due to the interchange of molecules between these shells.

Conditions for the emergence of life on the early Earth: summary and reflections.

This review attempts to situate the emergence of life on the early Earth within the scientific issues of the operational and mechanistic description of life, the conditions and constraints of prebiotic chemistry, together with bottom-up molecular fabrication and biomolecular nanofabrication and top-down miniaturization approaches to the origin of terrestrial life.

Cluster dynamics transcending chemical dynamics toward nuclear fusion.

Ultrafast cluster dynamics encompasses femtosecond nuclear dynamics, attosecond electron dynamics, and electron-nuclear dynamics in ultraintense laser fields (peak intensities 10(15)-10(20) W.cm(-2)). Extreme cluster multielectron ionization produces highly charged cluster ions, e.g., (C(4+)(D(+))(4))(n) and (D(+)I(22+))(n) at I(M) = 10(18) W.cm(-2), that undergo Coulomb explosion (CE) with the production of high-energy (5 keV to 1 MeV) ions, which can trigger nuclear reactions in an assembly of exploding clusters. The laser intensity and the cluster size dependence of the dynamics and energetics of CE of (D(2))(n), (HT)(n), (CD(4))(n), (DI)(n), (CD(3)I)(n), and (CH(3)I)(n) clusters were explored by electrostatic models and molecular dynamics simulations, quantifying energetic driving effects, and kinematic run-over effects. The optimization of table-top dd nuclear fusion driven by CE of deuterium containing heteroclusters is realized for light-heavy heteroclusters of the largest size, which allows for the prevalence of cluster vertical ionization at the highest intensity of the laser field. We demonstrate a 7-orders-of-magnitude enhancement of the yield of dd nuclear fusion driven by CE of light-heavy heteroclusters as compared with (D(2))(n) clusters of the same size. Prospective applications for the attainment of table-top nucleosynthesis reactions, e.g., (12)C(P,gamma)(13)N driven by CE of (CH(3)I)(n) clusters, were explored.

Electron bubbles in helium clusters. I. Structure and energetics.

In this paper we present a theoretical study of the structure, energetics, potential energy surfaces, and energetic stability of excess electron bubbles in ((4)He)(N) (N=6500-10(6)) clusters. The subsystem of the helium atoms was treated by the density functional method. The density profile was specified by a void (i.e., an empty bubble) at the cluster center, a rising profile towards a constant interior value (described by a power exponential), and a decreasing profile near the cluster surface (described in terms of a Gudermannian function). The cluster surface density profile width (approximately 6 A) weakly depends on the bubble radius R(b), while the interior surface profile widths (approximately 4-8 A) increase with increasing R(b). The cluster deformation energy E(d) accompanying the bubble formation originates from the bubble surface energy, the exterior cluster surface energy change, and the energy increase due to intracluster density changes, with the latter term providing the dominant contribution for N=6500-2 x 10(5). The excess electron energy E(e) was calculated at a fixed nuclear configuration using a pseudopotential method, with an effective (nonlocal) potential, which incorporates repulsion and polarization effects. Concurrently, the energy V(0) of the quasi-free-electron within the deformed cluster was calculated. The total electron bubble energies E(t)=E(e)+E(d), which represent the energetic configurational diagrams of E(t) vs R(b) (at fixed N), provide the equilibrium bubble radii R(b) (c) and the corresponding total equilibrium energies E(t) (e), with E(t) (e)(R(e)) decreasing (increasing) with increasing N (i.e., at N=6500, R(e)=13.5 A and E(t) (e)=0.86 eV, while at N=1.8 x 10(5), R(e)=16.6 A and E(t) (e)=0.39 eV). The cluster size dependence of the energy gap (V(0)-E(t) (e)) allows for the estimate of the minimal ((4)He)(N) cluster size of N approximately 5200 for which the electron bubble is energetically stable.

Electron bubbles in helium clusters. II. Probing superfluidity.

In this paper we present calculations of electron tunneling times from the ground electronic state of excess electron bubbles in ((4)He)(N) clusters (N=6500-10(7), cluster radius R=41.5-478 A), where the equilibrium bubble radius varies in the range R(b)=13.5-17.0 A. For the bubble center located at a radial distance d from the cluster surface, the tunneling transition probability was expressed as A(0)phi(d,R)exp(-betad), where beta approximately 1 A(-1) is the exponential parameter, A(0) is the preexponential factor for the bubble located at the cluster center, and phi(d,R) is a correction factor which accounts for cluster curvature effects. Electron tunneling dynamics is grossly affected by the distinct mode of motion of the electron bubble in the image potential within the cluster, which is dissipative (i.e., tau(D)>tau(0)) in superfluid ((4)He)(N) clusters, where tau(D) is the bubble motional damping time (tau(D) approximately 4 x 10(-12) s for normal fluid clusters and tau(D) approximately 10 s for superfluid clusters), while tau(0) approximately 10(-9)-10(-10) s is the bubble oscillatory time. Exceedingly long tunneling lifetimes, which cannot be experimentally observed, are manifested from bubbles damped to the center of the normal fluid cluster, while for superfluid clusters electron tunneling occurs from bubbles located in the vicinity of the initial distance d near the cluster boundary. Model calculations of the cluster size dependence of the electron tunneling time (for a fixed value of d=38-39 A), with lifetimes increasing in the range of 10(-3)-0.3 s for N=10(4)-10(7), account well for the experimental data [M. Farnik and J. P. Toennies, J. Chem. Phys. 118, 4176 (2003)], manifesting cluster curvature effects on electron tunneling dynamics. The minimal cluster size for the dynamic stability of the bubble was estimated to be N=3800, which represents the threshold cluster size for which the excess electron bubble in ((4)He)(N) (-) clusters is amenable to experimental observation.

Reflections on physical chemistry: Science and scientists.

This is the story of a young person who grew up in Tel-Aviv during the period of the establishment of the State of Israel and was inspired to become a physical chemist by the cultural environment, by the excellent high-school education, and by having been trained by some outstanding scientists at the Hebrew University of Jerusalem and, subsequently, by the intellectual environment and high-quality scientific endeavor at the University of Chicago. Since serving as the first chairman of the Chemistry Department of the newly formed Tel-Aviv University he has been immersed in research, in the training of young scientists, and in intensive and extensive international scientific collaboration. Together with the members of his "scientific family" he has explored the phenomena of energy acquisition, storage and disposal and structure-dynamics-function relations in large molecules, condensed phase, clusters and biomolecules, and is looking forward to many future adventures in physical chemistry.

Fragmentation channels of large multicharged clusters.

We address unifying features of fragmentation channels driven by long-range Coulomb or pseudo-Coulomb forces in clusters, nuclei, droplets, and optical molasses. We studied the energetics, fragmentation patterns, and dynamics of multicharged (A+)n (n=55, 135, 321) clusters. In Morse clusters the variation of the range of the pair-potential induced changes in the cluster surface energy and in the fissibility parameter X=E(Coulomb)2E(surface). X was varied in the range of X=1-8 for short-range interactions and of X=0.1-1.0 for long-range interactions. Metastable cluster configurations were prepared by vertical ionization of the neutral clusters and by subsequent structural equilibration. The energetics of these metastable ionic clusters was described in terms of the liquid drop model, with the coefficients of the volume and surface energies depending linearly on the Morse band dissociation energy. Molecular-dynamics simulations established two distinct fragmentation patterns of multicharged clusters that involve cluster fission into a small number of large, multicharged clusters for X<1 and Coulomb explosion into a large number of individual ions and small ionic fragments for X>1. The Rayleigh instability limit X=1 separates between spatially anisotropic fission and spatially isotropic Coulomb explosion. Distinct features of the fragmentation energetics and dynamics were unveiled. For fission of n=55 clusters, large kinetic and internal energies of the large fragments are exhibited and the characteristic fragmentation time is approximately 700 fs, while for Coulomb explosion the major energy content of the small fragments involves kinetic energy and the characteristic fragmentation time of approximately 300 fs is shorter. The Rayleigh (X=1) limit, leading to isotropic Coulomb explosion, is transcended by a marked enhancement of the Coulomb energy, which is realized for extremely ionized clusters in ultraintense laser fields, or by a dramatic reduction of the surface energy as is the case for the expansion of optical molasses.

Femtosecond time-resolved geometry relaxation and ultrafast intramolecular energy redistribution in Ag2Au.

The ultrafast dynamics of the bimetallic cluster Ag2Au is investigated by pump-probe negative ion-to-neutral-to-positive ion (NeNePo) spectroscopy. Preparation of the neutral cluster in a highly nonequilibrium state by electron detachment from the mass-selected anion, and subsequent probing of the neutral nuclear dynamics through two-photon ionization to the cationic state, leads to strongly probe-energy-dependent transient cation-abundance signals. The origin of this pronounced time and wavelength dependence of the ionization probability on the femtosecond scale is revealed by ab initio theoretical simulations of the transient spectra. Based on the analysis of underlying dynamics, two fundamental processes involving geometry relaxation from linear to triangular structure followed by ultrafast intramolecular vibrational energy redistribution (IVR) have been identified and for the first time experimentally observed in the frame of NeNePo spectroscopy under conditions close to zero electron kinetic energy.

Regular multicharged transient soft matter in Coulomb explosion of heteroclusters.

Nanointerfaces of mobile, thin spherical shells of light ions that expand on the femtosecond time scale, can be produced by Coulomb explosion of extremely ionized molecular heteroclusters consisting of light and heavy ions, e.g., (D+Iq+)n (q = 7-35), which are generated in ultraintense laser fields (intensity, I, = 10(16) to 10(20) W.cm(-2)). Modeling, together with molecular dynamics simulations, reveals the expansion of 2D monolayers with high energies and narrow energy distributions [e.g., Eav approximately 23 keV and DeltaE/Eav = 0.16 for D+ from (D+I25+)(2171)] arising from kinematic run-over effects. The expanding regular, monoionic, spherical nanointerfaces manifest the attainment of transient self-organization in complex systems driven by repulsive Coulomb interactions.

Electron and nuclear dynamics of molecular clusters in ultraintense laser fields. IV. Coulomb explosion of molecular heteroclusters.

In this paper we present a theoretical and computational study of the temporal dynamics and energetics of Coulomb explosion of (CD4)(n) and (CH4)(n) (n=55-4213) molecular heteroclusters in ultraintense (I=10(16)-10(19) W cm(-2)) laser fields, addressing the manifestation of electron dynamics, together with nuclear energetic and kinematic effects on the heterocluster Coulomb instability. The manifestations of the coupling between electron and nuclear dynamics were explored by molecular dynamics simulations for these heteroclusters coupled to Gaussian laser fields (pulse width tau=25 fs), elucidating outer ionization dynamics, nanoplasma screening effects (being significant for I< or =10(17) W cm(-2)), and the attainment of cluster vertical ionization (CVI) (at I=10(17) W cm(-2) for cluster radius R(0)< or =31 A). Nuclear kinematic effects on heterocluster Coulomb explosion are governed by the kinematic parameter eta=q(C)m(A)/q(A)m(C) for (CA(4))(n) clusters (A=H,D), where q(j) and m(j) (j=A,C) are the ionic charges and masses. Nonuniform heterocluster Coulomb explosion (eta >1) manifests an overrun effect of the light ions relative to the heavy ions, exhibiting the expansion of two spatially separated subclusters, with the light ions forming the outer subcluster at the outer edge of the spatial distribution. Important features of the energetics of heterocluster Coulomb explosion originate from energetic triggering effects of the driving of the light ions by the heavy ions (C(4+) for I=10(17)-10(18) W cm(-2) and C(6+) for I=10(19) W cm(-2)), as well as for kinematic effects. Based on the CVI assumption, scaling laws for the cluster size (radius R(0)) dependence of the energetics of uniform Coulomb explosion of heteroclusters (eta=1) were derived, with the size dependence of the average (E(j,av)) and maximal (E(j,M)) ion energies being E(j,av)=aR(0) (2) and E(j,M)=(5a/3)R(0) (2), as well as for the ion energy distributions P(E(j)) proportional to E(j) (1/2); E(j)< or =E(j,M). These results for uniform Coulomb explosion serve as benchmark reference data for the assessment of the effects of nonuniform explosion, where the CVI scaling law for the energetics still holds, with deviations of the a coefficient, which increase with increasing eta. Kinematic effects (for eta>1) result in an isotope effect, predicting the enhancement (by 9%-11%) of E(H,av) for Coulomb explosion of (C(4+)H(4) (+))(eta) (eta=3) relative to E(D,av) for Coulomb explosion of (C(4+)D(4) (+))(eta) (eta=1.5), with the isotope effect being determined by the ratio of the kinematic parameters for the pair of Coulomb exploding clusters. Kinematic effects for nonuniform explosion also result in a narrow isotope dependent energy distribution (of width DeltaE) of the light ions (with DeltaE/E(H,av) approximately 0.3 and DeltaE/E(D,av) approximately 0.4), with the distribution peaking at the high energy edge, in marked contrast with the uniform explosion case. Features of laser-heterocluster interactions were inferred from the analyses of the intensity dependent boundary radii (R(0))(I) and the corresponding average D+ ion energies (E(D,av))(I), which provide a measure for optimization of the cluster size at intensity I for the neutron yield from dd nuclear fusion driven by Coulomb explosion (NFDCE) of these heteroclusters. We infer on the advantage of deuterium containing heteronuclear clusters, e.g., (CD4)(n) in comparison to homonuclear clusters, e.g., (D2)(n/2), for dd NFDCE, where the highly charged heavy ions (e.g., C4+ or C6+) serve as energetic and kinematic triggers driving the D+ ions to a high (10-200 keV) energy domain.

Electron and nuclear dynamics of molecular clusters in ultraintense laser 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 (D2)n/2 (n = 480 - 7.6 x 10(4), cluster radius R0 = 13.1 - 70 A) 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 EM and the average energy Eav 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 Eav, EM proportional to R0(2) and that the energy distribution is P(E) proportional to E1/2. For a rectangular laser pulse the CVI conditions are satisfied when complete outer ionization is obtained, with the outer ionization time toi being shorter than both the pulse width and the cluster radius doubling time tau2. By increasing toi, due to the increase of R0 or the decrease of I, we have shown that the deviation of Eav from the corresponding CVI value (Eav(CVI)) is (Eav(CVI) - Eav)/Eav(CVI) approximately (toi/2.91tau2)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 R0 < or = (R0)I, with (R0)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 (R0)I and the corresponding ion average energy (Eav)I were inferred from simulations and described in terms of an electrostatic model. Two independent estimates of (R0)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 (Eav)I proportional to (R0)I(2) 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 (Eav = 5 - 7 keV at I = 2 x 10(18) W cm(-2)), together with our simulation data, lead to the estimates of R0 = 51 - 60 A, which exceed the experimental estimate of R0 = 45 A. The predicted anisotropy of the D+ ion energies in the Coulomb explosion at I = 10(18) W cm(-2) is in accord with experiment. We also explored the laser frequency dependence of the energetics of Coulomb explosion in the range nu = 0.1 - 2.1 fs(-1) (lambda = 3000 - 140 nm), which can be rationalized in terms of the electrostatic model.

Electron and nuclear dynamics of molecular clusters in ultraintense laser fields. I. Extreme multielectron ionization.

In this paper we present a theoretical and computational study of extreme multielectron ionization (involving the stripping of all the electrons from light, first-row atoms, and the production of heavily charged ions, e.g., Xe(+q) (q< or =36) from heavy atoms) in elemental and molecular clusters of Xe(n),(D(2))(n), and (CD(4))(n) (n=55-1061) in ultraintense (intensity I=10(15)-10(19) W cm(-2)) laser fields. Single atom or molecule multielectron ionization can be adequately described by the semiclassical barrier suppression ionization (BSI) mechanism. Extreme cluster multielectron ionization is distinct from that of a single atomic or molecular species in terms of the mechanisms, the ionization level and the time scales for electron dynamics and for nuclear motion. The novel compound mechanism of cluster multielectron ionization, which applies when the cluster size (radius R(0)) considerably exceeds the barrier distance for the BSI of a single constituent, involves a sequential-parallel, inner-outer ionization. The cluster inner ionization driven by the BSI for the constituents is induced by a composite field consisting of the laser field and inner fields. The energetics and dynamics of the system consisting of high energy (< or =3 keV) electrons and of less, similar 100 keV ions in the laser field was treated by molecular dynamics simulations, which incorporate electron-electron, electron-ion, ion-ion, and charge-laser interactions. High-energy electron dynamics also incorporates relativistic effects and includes magnetic field effects. We treat inner ionization considering inner field ignition, screening and fluctuation contributions as well as small [(< or =13%)] impact ionization contributions. Subsequent to inner ionization a charged nanoplasma is contained within the cluster, whose response to the composite (laser+inner) field results in outer ionization, which can be approximately described by an entire cluster barrier suppression ionization mechanism.

Electron and nuclear dynamics of molecular clusters in ultraintense laser fields. II. Electron dynamics and outer ionization of the nanoplasma.

We explore electron dynamics in molecular (CD4)(1061) clusters and elemental Xen (n=249-2171) clusters, responding to ultraintense (intensity I=10(16)-10(19) W cm(-2)) laser fields. Molecular dynamics simulations (including magnetic field and relativistic effects) and analyses of high-energy electron dynamics and nuclear ion dynamics in a cluster interacting with a Gaussian shaped laser field (frequency 0.35 fs(-1), photon energy 1.44 eV, phase 0, temporal width 25 fs) elucidated the time dependence of inner ionization, the formation of a nanoplasma of unbound electrons within the cluster or its vicinity, and of outer ionization. We determined the cluster size and the laser intensity dependence of these three sequential-parallel electronic processes. The characteristic times for cluster inner ionization (tau(ii)) and for outer ionization (tau(oi)) fall in the femtosecond time domain, i.e., tau(ii)=2-9 fs and tau(oi)=4-15 fs for (CD4)(1061), tau(ii)=7-30 fs and tau(oi)=5-13 fs for Xe(n) (n=479,1061), with both tau(ii) and tau(oi) decreasing with increasing I, in accord with the barrier suppression ionization mechanism for inner ionization of the constituents and the cluster barrier suppression ionization mechanism for outer ionization. The positive delay times Deltatau(OI) between outer and inner ionization (e.g., Deltatau(OI)=6.5 fs for Xen at I=10(16) W cm(-2) and Deltatau(OI)=0.2 fs for (CD4)(1061) at I=10(19) W cm(-2)) demonstrate that the outer/inner ionization processes are sequential. For (CD4)(1061), tau(ii)tau(oi), reflecting on the energetic hierarchy in the ionization of the Xe atoms. Quasiresonance contributions to the outer ionization of the nanoplasma were established, as manifested in the temporal oscillations in the inner/outer ionization levels, and in the center of mass of the nanoplasma electrons. The formation characteristics, dynamics, and response of the nanoplasma in molecular or elemental clusters were addressed. The nanoplasma is positively charged, with a high-average electron density [rho(P)=(2-3)10(22) cm(-3)], being characterized by high-average electron energies epsilon(av) (e.g., in Xe(1061) clusters epsilon(av)=54 eV at I=10(16) W cm(-2) and epsilon(av)=0.56-0.37 keV at I=10(18) W cm(-2), with epsilon(av) proportional, variant I(1/2)). Beyond the cluster boundary the average electron energy markedly increases, reaching electron energies in the range of 1.2-40 keV for outer ionization of Xe(n) (n=249-2171) clusters. The nanoplasma exhibits spatial inhomogeneity and angular anisotropy induced by the laser field. Femtosecond time scales are predicted for the nanoplasma production (rise times 7-3 fs), for the decay (decay times approximately 5 fs), and for the persistence time (30-10 fs) of a transient nanoplasma at I=10(17)-10(18) W cm(-2). At lower intensities of I=10(16) W cm(-2) a persistent nanoplasma with a "long" lifetime of > 50 fs will prevail.