Nuclear Fusion Driven by Coulomb Explosion of Methane Clusters

Time-resolved Rayleigh scattering measurements of methane clusters for laser-cluster fusion experiments

We present a time-resolved analysis of Rayleigh scattering measurements to determine the average size of methane clusters and find the optimum timing for laser-cluster fusion experiments. We measure Rayleigh scattering and determine the average size of methane clusters varying the backing pressure (P₀) from 11 bar to 69 bar. Regarding the onset of clustering, we estimate that the average size of methane clusters at the onset of clustering is N_c0≅20 at 11 bar. According to our measurements, the average cluster radius r follows the power law of r∝P₀^1.86. Our ion time-of-flight measurements indicate that we have produced energetic deuterium ions with kT = 52±2 keV after laser-cluster interaction using CD₄ gas at 50 bar. We find that this ion temperature agrees with the predicted temperature from CD₄ clusters at 50 bar with r = 14 nm assuming the Coulomb explosion model.

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).

Angular distribution of neutrons from deuterated cluster explosions driven by femtosecond laser pulses

We have studied experimentally the angular distributions of fusion neutrons from plasmas of multi-keV ion temperature, created by 40 fs, multi-TW laser pulses in dense plumes of D2 and CD4 clusters. A slight anisotropy in the neutron emission is observed. We attribute this anisotropy to the fact that the differential cross section for DD fusion is anisotropic even at low collision energies, and this, coupled with the geometry of the gas jet target, leads to beam-target neutrons that are slightly directed. The qualitative features of this anisotropy are confirmed by Monte Carlo simulations.

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.

REFLECTIONS ON PHYSICAL CHEMISTRY: Science and Scientists

▪ Abstract  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.

“What to leave out and what to put in? That's the problem.” Hugh Lofting, Doctor Dolittle's Zoo, 1925

Dynamic acceleration effects in explosions of laser-irradiated heteronuclear clusters

Intense, femtosecond irradiation of atomic and molecular clusters can initiate Coulomb explosions, generating particle energies sufficient to drive nuclear fusion. Last and Jortner have proposed, based on particle dynamics simulations, that heteronuclear clusters with a mixture of heavy and light ions will not explode by the simple, equilibrium Coulomb model but that dynamic effects can lead to a boosting of energy of the lighter ejected ions [Phys. Rev. Lett. 87, 033401 (2001)]. We present experimental confirmation of this theoretically predicted ion energy enhancement in methane clusters.

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), as appropriate for sequential outer/inner ionization dynamics, while for Xen clusters 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.

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.

Shock shells in Coulomb explosions of nanoclusters

We predict that Coulomb explosion of a nanoscale cluster, which is ionized by high-intensity laser radiation and has a naturally occurring spatial density profile, will invariably produce shock waves. In most typical situations, two shocks, a leading and a trailing one, form a shock shell that eventually encompasses the entire cluster. Being the first example of shock waves on the nanometer scale, this phenomenon promises interesting effects and applications, including high-rate nuclear reactions inside each individual cluster.