Generation of periodic ultrashort electron bunches and strongly asymmetric ion Coulomb explosion in nanometer foils interacting with ultra-intense laser pulse

The interaction of a linearly polarized intense laser pulse with an ultrathin nanometer plasma layer is investigated to understand the physics of the ion acceleration. It is shown by the computer simulation that the plasma response to the laser pulse comprises two steps. First, due to the vxB effect, electrons in the plasma layer are extracted and periodic ultrashort relativistic electron bunches are generated every half of a laser period. Second, strongly asymmetric Coulomb explosion of ions in the foil occurs due to the strong electrostatic charge separation, once the foil is burnt through. Followed by the laser accelerated electron bunch, the ion expansion in the forward direction occurs along the laser beam that is much stronger as compared to the backward direction. (c) 2008 American Institute of Physics.

Influence of Coulomb potential for photoionization of H atoms in an elliptically polarized laser field: Velocity gauge versus length gauge

We theoretically study the influence of Coulomb potential for photoionization of hydrogen atoms in an intense laser field with elliptical polarization. The total ionization rates, photoelectron energy spectra, and photoelectron angular distributions are calculated with the Coulomb-Volkov wave functions in the velocity gauge and compared with those calculated in the length gauge as well as those calculated with the Volkov wave functions. By comparing the results obtained by the Coulomb-Volkov and Volkov wave functions, we find that for linear polarization the influence of Coulomb potential is obvious for low-energy photoelectrons, and as the photoelectron energy and/or the laser intensity increase, its influence becomes smaller. This trend, however, is not so clear for the case of elliptical polarization. We also find that the twofold symmetry in the photoelectron angular distributions for elliptical polarization is caused by the cooperation of Coulomb potential and interference of multiple transition channels. About the gauge issue, we show that the difference in the photoelectron angular distributions obtained by the velocity and length gauges becomes rather obvious for elliptical polarization, while the difference is generally smaller for linear polarization.

Crystal Coulomb energies: I. Organic donor-acceptor complexes--a review. II. Classical Ewald calculations of the Coulomb binding energy of four donor-acceptor crystals and Wurster's blue perchlorate. III. A rapid-convergence quantum-mechanical formalism for crystal electronic energies

Chapter I

Theories for organic donor-acceptor (DA) complexes in solution and in the solid state are reviewed, and compared with the available experimental data. As shown by McConnell et al. (Proc. Natl. Acad. Sci. U.S., 53, 46-50 (1965)), the DA crystals fall into two classes, the holoionic class with a fully or almost fully ionic ground state, and the nonionic class with little or no ionic character. If the total lattice binding energy 2ε₁ (per DA pair) gained in ionizing a DA lattice exceeds the cost 2ε_o of ionizing each DA pair, ε₁ + ε_o less than 0, then the lattice is holoionic. The charge-transfer (CT) band in crystals and in solution can be explained, following Mulliken, by a second-order mixing of states, or by any theory that makes the CT transition strongly allowed, and yet due to a small change in the ground state of the non-interacting components D and A (or D⁺ and A^-). The magnetic properties of the DA crystals are discussed.

Chapter II

A computer program, EWALD, was written to calculate by the Ewald fast-convergence method the crystal Coulomb binding energy E_C due to classical monopole-monopole interactions for crystals of any symmetry. The precision of E_C values obtained is high: the uncertainties, estimated by the effect on E_C of changing the Ewald convergence parameter η, ranged from ± 0.00002 eV to ± 0.01 eV in the worst case. The charge distribution for organic ions was idealized as fractional point charges localized at the crystallographic atomic positions: these charges were chosen from available theoretical and experimental estimates. The uncertainty in E_C due to different charge distribution models is typically ± 0.1 eV (± 3%): thus, even the simple Hückel model can give decent results.

E_C for Wurster's Blue Perchl orate is -4.1 eV/molecule: the crystal is stable under the binding provided by direct Coulomb interactions. E_C for N-Methylphenazinium Tetracyanoquino- dimethanide is 0.1 eV: exchange Coulomb interactions, which cannot be estimated classically, must provide the necessary binding.

EWALD was also used to test the McConnell classification of DA crystals. For the holoionic (1:1)-(N,N,N',N'-Tetramethyl-para- phenylenediamine: 7,7,8,8-Tetracyanoquinodimethan) E_C = -4.0 eV while 2ε_o = 4.6₅ eV: clearly, exchange forces must provide the balance. For the holoionic (1:1)-(N,N,N',N'-Tetramethyl-para- phenylenediamine:para-Chloranil) E_C = -4.4 eV, while 2ε_o = 5.0 eV: again EC falls short of 2ε₁. As a Gedankenexperiment, two nonionic crystals were assumed to be ionized: for (1:1)-(Hexamethyl- benzene:para-Chloranil) E_C = -4.5 eV, 2ε_o = 6.6 eV; for (1:1)- (Napthalene:Tetracyanoethylene) E_C = -4.3 eV, 2ε_o = 6.5 eV. Thus, exchange energies in these nonionic crystals must not exceed 1 eV.

Chapter III

A rapid-convergence quantum-mechanical formalism is derived to calculate the electronic energy of an arbitrary molecular (or molecular-ion) crystal: this provides estimates of crystal binding energies which include the exchange Coulomb inter- actions. Previously obtained LCAO-MO wavefunctions for the isolated molecule(s) ("unit cell spin-orbitals") provide the starting-point. Bloch's theorem is used to construct "crystal spin-orbitals". Overlap between the unit cell orbitals localized in different unit cells is neglected, or is eliminated by Löwdin orthogonalization. Then simple formulas for the total kinetic energy Q^(XT)_λ, nuclear attraction [λ/λ]^XT, direct Coulomb [λλ/λ'λ']^XT and exchange Coulomb [λλ'/λ'λ]^XT integrals are obtained, and direct-space brute-force expansions in atomic wavefunctions are given. Fourier series are obtained for [λ/λ]^XT, [λλ/λ'λ']^XT, and [λλ/λ'λ]^XT with the help of the convolution theorem; the Fourier coefficients require the evaluation of Silverstone's two-center Fourier transform integrals. If the short-range interactions are calculated by brute-force integrations in direct space, and the long-range effects are summed in Fourier space, then rapid convergence is possible for [λ/λ]^XT, [λλ/λ'λ']^XT and [λλ'/λ'λ]^XT. This is achieved, as in the Ewald method, by modifying each atomic wavefunction by a "Gaussian convergence acceleration factor", and evaluating separately in direct and in Fourier space appropriate portions of [λ/λ]^XT, etc., where some of the portions contain the Gaussian factor.

Coulomb explosion of hydrogen clusters irradiated by an ultrashort intense laser pulse

The explosion dynamics of hydrogen clusters driven by an ultrashort intense laser pulse has been analyzed analytically and numerically by employing a simplified Coulomb explosion model. The dependence of average and maximum proton kinetic energy on cluster size, pulse duration, and laser intensity has been investigated respectively. The existence of an optimum cluster size allows the proton energy to reach the maximum when the cluster size matches with the intensity and the duration of the laser pulse. In order to explain our experimental results such as the measured proton energy spectrum and the saturation effect of proton energy, the effects of cluster size distribution as well as the laser intensity distribution on the focus spot should be considered. A good agreement between them is obtained.

Effects on the proton production from explosions of hydrogen clusters subjected to intense femtosecond laser fields

The dependence of the maximum and average energies of protons, which were produced in the interaction of an intense laser pulse (similar to 1 x 10(16) W cm(-2), 65 fs) with hydrogen clusters in a gas jet backed up to 80 bar at liquid nitrogen temperature (similar to 80 K), on the backing pressure has been studied. The general trend of the proton energy dependence on the square of the average cluster radius, which is determined by a calibrated Rayleigh scattering measurement, is similar to that described by theory under the single size approximation. Calculations are made to fit the experimental results under a simplified model by taking into account both a log-normal cluster size distribution and the laser intensity attenuation in the interaction volume. A very good agreement between the experimental proton energy spectra and the calculations is obtained in the high- energy part of the proton energy distributions, but a discrepancy of the fits is revealed in the low-energy part at higher backing pressures which are associated with denser flows. A possible mechanism which would be responsible for this discrepancy is discussed. Finally, from the fits, a variation of the cluster size distributions was revealed to be dependent on the gas backing pressure as well as on the evolving time of the gas flow of clusters.

Two overrun phenomena and their effects on fusion yield in the Coulomb explosion of heteronuclear clusters

Two overrun effects in the Coulomb explosion dynamics of heteronuclear clusters have been investigated theoretically by the use of a simplified electrostatic model. When the charge-to-mass ratio of light ions is higher than that of heavy ions, the light ions can overtake the heavy ions inside the cluster and acquire a higher kinetic energy. Further, if the charge density of the heavy ions is twice as high as that of the light ions, i.e. a proposed competitive parameter xi = rho BqB/rho AqA > 2, the inner light ions can overtake those light ions on the surface of the cluster and form a shock shell during the explosion, which might drive the intracluster collision and fusion of the light ions. Different regimes of nuclear fusion are discussed and the corresponding neutron yields are estimated. Our analysis indicates that the probability of intracluster fusion is quite low even if deuterated heteronuclear clusters such as (DI)(n) with large size and high competitive parameter are employed. However, heteronuclear clusters are still a better candidate compared with homonuclear clusters for enhancing the total intercluster fusion yield because both a higher energy region and a higher proportion of deuterons distributing in the energy region can be created in the deuterated heteronuclear clusters.

Finite differences and a coupled analytic technique with applications to explosions and earthquakes

An analytic technique is developed that couples to finite difference calculations to extend the results to arbitrary distance. Finite differences and the analytic result, a boundary integral called two-dimensional Kirchhoff, are applied to simple models and three seismological problems dealing with data. The simple models include a thorough investigation of the seismologic effects of a deep continental basin. The first problem is explosions at Yucca Flat, in the Nevada test site. By modeling both near-field strong-motion records and teleseismic P-waves simultaneously, it is shown that scattered surface waves are responsible for teleseismic complexity. The second problem deals with explosions at Amchitka Island, Alaska. The near-field seismograms are investigated using a variety of complex structures and sources. The third problem involves regional seismograms of Imperial Valley, California earthquakes recorded at Pasadena, California. The data are shown to contain evidence of deterministic structure, but lack of more direct measurements of the structure and possible three-dimensional effects make two-dimensional modeling of these data difficult.

Coulomb drag and tunneling studies in quantum Hall bilayers

The bilayer quantum Hall state at total filling factor ν_T=1, where the total electron density matches the degeneracy of the lowest Landau level, is a prominent example of Bose-Einstein condensation of excitons. A macroscopically ordered state is realized where an electron in one layer is tightly bound to a "hole" in the other layer. If exciton transport were the only bulk transportmechanism, a current driven in one layer would spontaneously generate a current of equal magnitude and opposite sign in the other layer. The Corbino Coulomb drag measurements presented in this thesis demonstrate precisely this phenomenon.

Excitonic superfluidity has been long sought in the ν_T=1 state. The tunneling between the two electron gas layers exihibit a dc Josephson-like effect. A simple model of an overdamped voltage biased Josephson junction is in reasonable agreement with the observed tunneling I-V. At small tunneling biases, it exhibits a tunneling "supercurrent". The dissipation is carefully studied in this tunneling "supercurrent" and found to remain small but finite.