Coulomb expansion of deuterated methane clusters irradiated by an ultrashort intense laser pulse

The simulations of three-dimensional particle dynamics show that when irradiated by an ultrashort intense laser pulse, the deuterated methane cluster expands and the majority of deuterons overrun the more slowly expanding carbon ions, resulting in the creation of two separated subclusters. The enhanced deuteron kinetic energy and a narrow peak around the energy maximum in the deuteron energy distribution make a considerable contribution to the efficiency of nuclear fusion compared with the case of homonuclear deuterium clusters. With the intense laser irradiation, the nuclear fusion yield increases with the increase of the cluster size, so that deuterated heteronuclear clusters with larger sizes are required to achieve a greater neutron yield.

Direct acceleration of solid-density plasma bunch by ultraintense laser

The interaction of a petawatt laser with a small solid-density plasma bunch is studied by particle-in-cell simulation. It is shown that when irradiated by a laser of intensity >10(21) W/cm(2), a dense plasma bunch of micrometer size can be efficiently accelerated. The kinetic energy of the ions in the high-density region of the plasma bunch can exceed ten MeV at a density in the 10(23)-cm(-3) level. Having a flux density orders of magnitude higher than that of the traditional charged-particle pulses, the laser-accelerated plasma bunch can have a wide range of applications. In particular, such a dense energetic plasma bunch impinging on the compressed fuel in inertial fusion can significantly enhance the nuclear-reaction cross section and is thus a promising alternative for fast ignition.

Pondermotive shift of photoelectron spectra in intense laser field

In a recent experimental work on the excess photon detachment (EPD) of H- ions [Phys. Rev. Lett. 87 (2001) 243001] it has been found that the ponderomotive shift of each EPD peak increases with the order of the EPD channel. By using a nonperturbative quantum scattering theory, we obtain the kinetic energy spectra for the differential detachment rate along the laser polarization for several laser intensities. It is demonstrated that higher order EPD peaks are produced mainly at relatively higher laser intensities. By calculating the overall EPD spectra with varying laser intensities, it is found that the ponderomotive shift of each EPD peak increases with the order of the EPD channel. Our calculations are in good agreement with the experimental observation. It is found that different EPD channels occur mainly when the laser field reaches some values, thus the intensity distribution of the laser field is responsible for the varying ponderomotive shifts.

Dynamical stability of nascent neutron stars

This thesis presents a study of the dynamical stability of nascent neutron stars resulting from the accretion induced collapse of rapidly rotating white dwarfs.

Chapter 2 and part of Chapter 3 study the equilibrium models for these neutron stars. They are constructed by assuming that the neutron stars have the same masses, angular momenta, and specific angular momentum distributions as the pre-collapse white dwarfs. If the pre-collapse white dwarf is rapidly rotating, the collapsed object will contain a high density central core of size about 20 km, surrounded by a massive accretion torus extending to hundreds of kilometers from the rotation axis. The ratio of the rotational kinetic energy to gravitational binding energy, β, of these neutron stars is all found to be less than 0.27.

Chapter 3 studies the dynamical stability of these neutron stars by numerically evolving the linearized hydrodynamical equations. A dynamical bar-mode instability is observed when the β of the star is greater than the critical value β_d ≈ 0.25. It is expected that the unstable mode will persist until a substantial amount of angular momentum is carried away by gravitational radiation. The detectability of these sources is studied and it is estimated that LIGO II is unlikely to detect them unless the event rate is greater than 10^-6/year/galaxy.

All the calculations on the structure and stability of the neutron stars in Chapters 2 and 3 are carried out using Newtonian hydrodynamics and gravity. Chapter 4 studies the relativistic effects on the structure of these neutron stars. New techniques are developed and used to construct neutron star models to the first post-Newtonian (1PN) order. The structures of the 1PN models are qualitatively similar to the corresponding Newtonian models, but the values of β are somewhat smaller. The maximum β for these 1PN neutron stars is found to be 0.24, which is 8% smaller than the Newtonian result (0.26). However, relativistic effects will also change the critical value β_d. A detailed post-Newtonian stability analysis has yet to be carried out to study the relativistic effects on the dynamical stability of these neutron stars.

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.

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.

A scaling law of photoionization in few-cycle regime

In this paper, a scaling law of photoionization of atoms irradiated by intense, few- cycle laser pulses is established. The scaling law sets a relation to the phase- dependent ionization with the kinetic energy of photoelectrons, the duration and peak intensity of short pulses, and the ionization potential of the target atoms. We find that it will be advantageous to manifest the phase- dependent photoionization by choosing the target atoms with larger ionization potential, using laser with smaller carrier- frequency, and increasing the pulse intensity. (c) 2007 Optical Society of America.

Towards understanding the mixing characteristics of turbulent buoyant flows

This work proposes a new simulation methodology in which variable density turbulent flows can be studied in the context of a mixing layer with or without the presence of gravity. Specifically, this methodology is developed to probe the nature of non-buoyantly-driven (i.e. isotropically-driven) or buoyantly-driven mixing deep inside a mixing layer. Numerical forcing methods are incorporated into both the velocity and scalar fields, which extends the length of time over which mixing physics can be studied. The simulation framework is designed to allow for independent variation of four non-dimensional parameters, including the Reynolds, Richardson, Atwood, and Schmidt numbers. Additionally, the governing equations are integrated in such a way to allow for the relative magnitude of buoyant energy production and non-buoyant energy production to be varied.

The computational requirements needed to implement the proposed configuration are presented. They are justified in terms of grid resolution, order of accuracy, and transport scheme. Canonical features of turbulent buoyant flows are reproduced as validation of the proposed methodology. These features include the recovery of isotropic Kolmogorov scales under buoyant and non-buoyant conditions, the recovery of anisotropic one-dimensional energy spectra under buoyant conditions, and the preservation of known statistical distributions in the scalar field, as found in other DNS studies.

This simulation methodology is used to perform a parametric study of turbulent buoyant flows to discern the effects of varying the Reynolds, Richardson, and Atwood numbers on the resulting state of mixing. The effects of the Reynolds and Atwood numbers are isolated by looking at two energy dissipation rate conditions under non-buoyant (variable density) and constant density conditions. The effects of Richardson number are isolated by varying the ratio of buoyant energy production to total energy production from zero (non-buoyant) to one (entirely buoyant) under constant Atwood number, Schmidt number, and energy dissipation rate conditions. It is found that the major differences between non-buoyant and buoyant turbulent flows are contained in the transfer spectrum and longitudinal structure functions, while all other metrics are largely similar (e.g. energy spectra, alignment characteristics of the strain-rate tensor). Also, despite the differences noted between fully buoyant and non-buoyant turbulent fields, the scalar field, in all cases, is unchanged by these. The mixing dynamics in the scalar field are found to be insensitive to the source of turbulent kinetic energy production (non-buoyant vs. buoyant).

Generation of an extreme ultraviolet supercontinuum in a two-color laser field

We theoretically investigate the high-order harmonic generation in a helium atom with a two-color optical field synthesized by an intense 6 fs pulse at 800 nm and a relatively weak 21.3 fs pulse at 400 nm. When the frequency-doubled pulse is properly time shifted with respect to the fundamental pulse, an ultrabroad extreme ultraviolet supercontinuum spectrum with a 148 eV spectral width can be generated which directly creates an isolated 65 as pulse even without phase compensation. We explain this extraordinary phenomenon by analyzing maximum electron kinetic energies at different return times.

甲烷团簇纯库仑爆炸的离子平均动能

利用一简单的经典静电模型研究了甲烷团簇纯库仑爆炸情况下产生的离子平均动能。研究表明，甲烷团簇爆炸后离子获得的平均动能和离子的初始平均静电势能的比值，与团簇的尺寸大小无关。这意味着在甲烷团簇纯库仑爆炸近似下，不必使用分子动力学模拟，离子获得的动能可以通过这一比值以及离子的初始静电势能进行估算。给出了不同碳离子价态下的离子平均动能和其初始平均静电势能的比值。

Impact of laser-accelerated micron-size projectile on dense plasma

The impact of a laser-accelerated micron-size projectile on a dense plasma target is studied using two-dimensional particle-in-cell simulations. The projectile is first accelerated by an ultraintense laser. It then impinges on the dense plasma target and merges with the latter. Part of the kinetic energy of the laser-accelerated ions in the projectile is deposited in the fused target, and an extremely high concentration of plasma ions with a mean kinetic energy needed for fusion reaction is induced. The interaction is thus useful for laser-driven impact fusion and as a compact neutron source.

Density functional theory embedding for correlated wavefunctions

Methods that exploit the intrinsic locality of molecular interactions show significant promise in making tractable the electronic structure calculation of large-scale systems. In particular, embedded density functional theory (e-DFT) offers a formally exact approach to electronic structure calculations in which the interactions between subsystems are evaluated in terms of their electronic density. In the following dissertation, methodological advances of embedded density functional theory are described, numerically tested, and applied to real chemical systems.

First, we describe an e-DFT protocol in which the non-additive kinetic energy component of the embedding potential is treated exactly. Then, we present a general implementation of the exact calculation of the non-additive kinetic potential (NAKP) and apply it to molecular systems. We demonstrate that the implementation using the exact NAKP is in excellent agreement with reference Kohn-Sham calculations, whereas the approximate functionals lead to qualitative failures in the calculated energies and equilibrium structures.

Next, we introduce density-embedding techniques to enable the accurate and stable calculation of correlated wavefunction (CW) in complex environments. Embedding potentials calculated using e-DFT introduce the effect of the environment on a subsystem for CW calculations (WFT-in-DFT). We demonstrate that WFT-in-DFT calculations are in good agreement with CW calculations performed on the full complex.

We significantly improve the numerics of the algorithm by enforcing orthogonality between subsystems by introduction of a projection operator. Utilizing the projection-based embedding scheme, we rigorously analyze the sources of error in quantum embedding calculations in which an active subsystem is treated using CWs, and the remainder using density functional theory. We show that the embedding potential felt by the electrons in the active subsystem makes only a small contribution to the error of the method, whereas the error in the nonadditive exchange-correlation energy dominates. We develop an algorithm which corrects this term and demonstrate the accuracy of this corrected embedding scheme.

Light-induced evaporative cooling in a magneto-optical trap

This paper presents an experimental demonstration of light-induced evaporative cooling in a magneto-optical trap. An additional laser is used to interact with atoms at the edge of the atomic cloud in the trap. These atoms get an additional force and evaporated away from the trap by both the magnetic field and laser fields. There remaining atoms have lower kinetic energy and thus are cooled. It reports the measurements on the temperature and atomic number after the evaporative cooling with different parameters including the distance between the laser and the centre of the atomic cloud, the detuning, the intensity. The results show that the light-induced evaporative cooling is a way to generate an ultra-cold atom source.

Development and Application of Embedding Methods for the Simulation of Large Chemical Systems

The high computational cost of correlated wavefunction theory (WFT) calculations has motivated the development of numerous methods to partition the description of large chemical systems into smaller subsystem calculations. For example, WFT-in-DFT embedding methods facilitate the partitioning of a system into two subsystems: a subsystem A that is treated using an accurate WFT method, and a subsystem B that is treated using a more efficient Kohn-Sham density functional theory (KS-DFT) method. Representation of the interactions between subsystems is non-trivial, and often requires the use of approximate kinetic energy functionals or computationally challenging optimized effective potential calculations; however, it has recently been shown that these challenges can be eliminated through the use of a projection operator. This dissertation describes the development and application of embedding methods that enable accurate and efficient calculation of the properties of large chemical systems.

Chapter 1 introduces a method for efficiently performing projection-based WFT-in-DFT embedding calculations on large systems. This is accomplished by using a truncated basis set representation of the subsystem A wavefunction. We show that naive truncation of the basis set associated with subsystem A can lead to large numerical artifacts, and present an approach for systematically controlling these artifacts.

Chapter 2 describes the application of the projection-based embedding method to investigate the oxidative stability of lithium-ion batteries. We study the oxidation potentials of mixtures of ethylene carbonate (EC) and dimethyl carbonate (DMC) by using the projection-based embedding method to calculate the vertical ionization energy (IE) of individual molecules at the CCSD(T) level of theory, while explicitly accounting for the solvent using DFT. Interestingly, we reveal that large contributions to the solvation properties of DMC originate from quadrupolar interactions, resulting in a much larger solvent reorganization energy than that predicted using simple dielectric continuum models. Demonstration that the solvation properties of EC and DMC are governed by fundamentally different intermolecular interactions provides insight into key aspects of lithium-ion batteries, with relevance to electrolyte decomposition processes, solid-electrolyte interphase formation, and the local solvation environment of lithium cations.

Storm track response to perturbations in climate

This thesis advances our understanding of midlatitude storm tracks and how they respond to perturbations in the climate system. The midlatitude storm tracks are regions of maximal turbulent kinetic energy in the atmosphere. Through them, the bulk of the atmospheric transport of energy, water vapor, and angular momentum occurs in midlatitudes. Therefore, they are important regulators of climate, controlling basic features such as the distribution of surface temperatures, precipitation, and winds in midlatitudes. Storm tracks are robustly projected to shift poleward in global-warming simulations with current climate models. Yet the reasons for this shift have remained unclear. Here we show that this shift occurs even in extremely idealized (but still three-dimensional) simulations of dry atmospheres. We use these simulations to develop an understanding of the processes responsible for the shift and develop a conceptual model that accounts for it.

We demonstrate that changes in the convective static stability in the deep tropics alone can drive remote shifts in the midlatitude storm tracks. Through simulations with a dry idealized general circulation model (GCM), midlatitude storm tracks are shown to be located where the mean available potential energy (MAPE, a measure of the potential energy available to be converted into kinetic energy) is maximal. As the climate varies, even if only driven by tropical static stability changes, the MAPE maximum shifts primarily because of shifts of the maximum of near-surface meridional temperature gradients. The temperature gradients shift in response to changes in the width of the tropical Hadley circulation, whose width is affected by the tropical static stability. Storm tracks generally shift in tandem with shifts of the subtropical terminus of the Hadley circulation.

We develop a one-dimensional diffusive energy-balance model that links changes in the Hadley circulation to midlatitude temperature gradients and so to the storm tracks. It is the first conceptual model to incorporate a dynamical coupling between the tropical Hadley circulation and midlatitude turbulent energy transport. Numerical and analytical solutions of the model elucidate the circumstances of when and how the storm tracks shift in tandem with the terminus of the Hadley circulation. They illustrate how an increase of only the convective static stability in the deep tropics can lead to an expansion of the Hadley circulation and a poleward shift of storm tracks.

The simulations with the idealized GCM and the conceptual energy-balance model demonstrate a clear link between Hadley circulation dynamics and midlatitude storm track position. With the help of the hierarchy of models presented in this thesis, we obtain a closed theory of storm track shifts in dry climates. The relevance of this theory for more realistic moist climates is discussed.