Doctor of Philosophy / Department of Physics / Itzhak Ben-Itzhak / Present laser technology allows for the production of ultra short (&7 fs) intense (.1016 W/cm2)pulses, which are comparable in duration and interaction strength to the vibrational period and the interaction that binds the electron in molecules, respectively. In this intense-field ultra short-pulse regime one can both measure and manipulate dynamics on the femtosecond timescale. To probe the dynamics of laser-matter interactions in this regime, we have chosen to start from the simplest possible molecule - H+ 2 , which can either dissociate into H + p or ionize into p + p + e. We have designed and employ a coincidence three-dimensional momentum imaging technique which allows us to measure ionization and dissociation of a molecular ion beam target simultaneously, while completely separating the two channels from each other. By varying the laser intensity and the pulse duration, we measure the intensity and pulse length dependent momentum distributions for laser induced fragmentation of H+ 2 at 790 nm. These dissociation measurements are in agreement with the phenomena predicted using the adiabatic Floquet picture, e.g. bond softening, in addition to more sophisticated calculations done by solving the time-dependent Schrodinger equation in the Born-Oppenheimer representation. Furthermore, the structure seen in ionization in our measurements and soon after by others is explained via a unified diabatic Floquet picture, which includes both ionization and dissociation in a single intensity and wavelength dependent picture that includes nuclear motion. Additionally, we use the same experimental techniques and apparatus to probe the laser-induced dynamics of multi-electron diatomic molecules, e.g. O+2, N+2, and ND+. The most probable dissociation and ionization pathways producing the features seen in these measurements are discerned using the angular and kinetic-energy-release distributions in conjunction with the diabatic Floquet picture. Finally, we extend these experimental techniques and interpretive models to the simplest polyatomic molecule - H+ 3 , whose fragmentation presents challenges both in our first-of-their-kind experiments and in physical interpretation.

Doctor of Philosophy / Department of Physics / Itzik Ben-Itzhak / One ultimate goal of ultrafast, strong- field laser science is to coherently control chemical reactions. Present laser technology allows for the production of intense (>10[superscript]13 W/cm[superscript]2), ultrashort ( 5 fs), carrier-envelope phase-stabilized pulses. By knowing the electric field waveform, sub-cycle resolution on the order of 100's of attoseconds (1 as=10[superscript]-18 s) can be reached -- the timescale for electron motion. Meanwhile, the laser field strengths are comparable to that which binds electrons to atoms or molecules. In this intense-field ultrashort-pulse regime one can both measure and manipulate dynamics of strong-field, quantum-mechanical processes in atoms and molecules. Despite much progress in the technology, typical durations for which lasers can be reliably locked to a specific carrier-envelope phase ranges from a few minutes to a few hours. Experiments investigating carrier-envelope phase effects that have necessarily long data acquisition times, such as those requiring coincidence between fragments originating from the same atom or molecule, are thus challenging and uncommon. Therefore, we combined the new technology for measuring the carrier-envelope phase of each and every laser shot with other single-shot coincidence three-dimensional momentum imaging techniques to alleviate the need for carrier-envelope phase stabilized laser pulses. Using phase-tagged coincidence techniques, several targets and laser-induced processes were studied. One particular highlight uses this method to study the recollision process of non-sequential double ionization of argon. By measuring the momentum of the two electrons emitted in the process, we could study their energy sharing. Furthermore, by selecting certain carrier-envelope phase values, and therefore laser pulses with a particular waveform, events with single recollision could be isolated and further analyzed. Another highlight is our studies of carrier-envelope phase effects in the dissociation of the benchmark H[subscript]2[superscript[+] ion beam. Aided by near-exact quantum mechanical calculations, we could identify interfering pathways which lead to the observed spatial asymmetry. These and other similar experiments are described in this thesis as significant steps toward their ultimate control.

Doctor of Philosophy / Department of Physics / Uwe Thumm / Although the latest and most powerful supercomputer today, Tianhe-2 in China, can finish 33.86 quadrillion floating-point operations per second (www.top500.org), it is still a big challenge to simulate the simplest few-electron system - the helium atom - a threebody system with one nucleus and two electrons. Within the fixed-nucleus approximation and time-dependent close coupling (TDCC) approach, we developed software to solve the time-dependent Schrödinger equation (TDSE) accurately, implementing the finite-element discrete-variable representation (FE-DVR) scheme. The general idea of the method is to expand the wave functions in the eigenvectors of the angular momentum operator, which further transform the six-dimensional TDSE to a set of infinite two-dimensional coupled equations. Although there are infinitely many coupled equations, they can be truncated to a finite number of equations by applying selection rules and physical requirements, and solved with our current computational resources. By numerically solving the TDSE in full dimensionality, we investigate the double photoionization of helium atoms in external fields. In co-planar emission geometry with and without the presence of a comparatively weak infrared (IR) laser pulse, we discuss the double ionization (DI) dynamics of helium atoms irradiated by ultrashort pulses of extreme ultraviolet (XUV) laser light. We first investigate the degree of electronic correlation by correlated photoelectron angular distributions for two-photon double ionization (TPDI) of helium atoms in the sequential and non-sequential DI regime. We quantify sequential and non-sequential contributions to TPDI driven by an XUV pulse with central photon energy hw[subscript]xuv near the sequential DI threshold. If the spectral width of the XUV pulse is broad enough, both the sequntial (hw[subscript]xuv > 54.4 eV) and non-sequential (hw[subscript]xuv < 54.4 eV) channels are open. Therefore, the sequential and non-sequential DI mechanisms are difficult to distinguish. By tracking the DI asymmetry in joint photoelectron angular distributions, we introduce the forward-backward-emission asymmetry as a measure that allows the distinction of sequential and non-sequential contributions. Specifically, for hw[subscript]xuv = 50 eV pulses with a sine-squared temporal profile, we find that the sequential DI contribution is the largest at a pulse length of 650 as (1 as = 10[superscript]−18 s), due to competing temporal and spectral constraints. In addition, we validate a simple heuristic expression for the sequential DI contribution in comparison with ab initio calculations. We then investigate the influence of the laser field on the DI of helium by a single XUV pulse. For IR-laser-assisted single-XUV-photon DI our joint angular distributions show that the IR-laser field enhances back-to-back electron emission and induces a characteristic splitting in the angular distribution for electrons that are emitted symmetrically relative to the identical linear polarization directions of the XUV and IR pulse. These IR-pulse-induced changes in photoelectron angular distributions are (i) imposed by different symmetry constraints for XUV-pulse-only and laser-assisted XUV-photon DI, (ii) robust over a large range of energy sharing between the emitted electrons, and (iii) consistent with the transfer of discrete IR-photon momenta to both photoelectrons from the assisting IR-laser field. While selection-rule forbidden at equal energy sharing, for increasingly unequal energy sharing we find back-to-back emission to become more likely and to compete with symmetric emission. To obtain a high level of accuracy, accurate quantum-mechanical calculations of three Coulomb interacting particales exposed to an intense XUV and weak IR field are at the limit of current computational power. Any direct extension (such as strong laser-field intensity, elliptically-polarized field, and laser-induced DI) of our approach to more complicated systems appears to be currently out of reach. At the end of this thesis, we give suggestions on how to improve the efficiency of TDSE calculations for simulations of these complicated many-photon processes.

Doctor of Philosophy / Department of Physics / Brett D. Esry / The Efimov effect plays a central role in few-body systems at ultracold temperature and has thus accelerated a lot of studies on its manifestation in the collisional stability of the quantum degenerate gases. Near broad Feshbach resonances, Efimov physics has been studied both theoretically and experimentally through the zero-energy scattering observables. We have extended the theoretical studies of Efimov physics to a much broader extent. In particular, we have investigated the three-body Efimov physics near narrow Feshbach resonances and have also identified the Efimov features beyond the zero temperature limit. We have found, near a narrow Feshbach resonance, the non-trivial contribution from both of the resonance width and the short-range physics to the three-body recombination and vibrational dimer relaxation. Remarkably, the collisional stability of the Feshbach molecules are found to be opposite to that near the broad resonances: an increased stability for molecules made by bosons and a decreased stability for those made by fermions. The universal physics observed near the narrow Feshbach resonances is further found not to be limited to the zero temperature observables. We have found that the general features of Efimov physics and those pertaining to a narrow resonance are manifested in different energy ranges above zero temperature. This opens the opportunity to observe Efimov physics by changing the collisional energy while keeping the atomic interaction fixed. The landscape of the universal Efimov physics is thus delineated in both of the interaction and the energy domain. We have also investigated Efimov physics in heteronuclear four-body systems where the complexity can be reduced by approximations. In particular, we have proposed ways for controllable production of the Efimov tri-atomic molecules by three-body or four-body recombinations involving four atoms. We have also confirmed the existence of four-body Efimov effect in a system of three heavy particles and one light particle, which has resolved a decade-long controversy on this topic. Finally, we have studied the collisional properties of four identical bosons in 1D, which is important to the experiments on the quantum gases confined in the 1D optical lattices.

Doctor of Philosophy / Department of Physics / Uwe Thumm / We study the dynamics of diatomic molecules (dimers) in intense IR and XUV laser fields theoretically and compare the results with measured data in collaboration with different experimental groups worldwide. The first three chapters of the thesis cover the introduction and the background on solving time-independent and time-dependent Schrödinger equation. The numerical results in this thesis are presented in four chapters, three of which are focused on diatomic molecules in IR fields. The last one concentrates on diatomic molecules in XUV pulses. The study of nuclear dynamics of H[subscript]2 or D[subscript]2 molecules in IR pulses is given in Chapter 4. First, we investigate the optimal laser parameters for observing field-induced bond softening and bond hardening in D[subscript]2[superscript]+. Next, the nuclear dynamics of H[subscript]2[superscript]+ molecular ions in intense laser fields are investigated by analyzing their fragment kinetic-energy release (KER) spectra as a function of the pump-probe delay τ. Lastly, the electron localization is studied for long circularly polarized laser pulses. Chapter 5 covers the dissociation dynamics of O[subscript]2[superscript]+ in an IR laser field. The fragment KER spectra are analyzed as a function of the pump-probe delay τ. Within the Born-Oppenheimer approximation, we calculate ab-initio adiabatic potential-energy curves and their electric dipole couplings, using the quantum chemistry code GAMESS. In Chapter 6, the dissociation dynamics of the noble gas dimer ions He[subscript]2[superscript]+, Ne[subscript]2[superscript]+, Ar[subscript]2[superscript]+, Kr[subscript]2[superscript]+, and Xe[subscript]2[superscript]+ is investigated in ultrashort pump and probe laser pulses of different wavelengths. We observe a striking ‘‘delay gap’’ in the pump-probe-delay-dependent KER spectrum only if the probe-pulse wavelength exceeds the pump-pulse wavelength. Comparing pump-probe-pulse-delay dependent KER spectra for different noble gas dimer cations, we quantitatively discuss quantum-mechanical versus classical aspects of the nuclear vibrational motion as a function of the nuclear mass. Chapter 7 focuses on diatomic molecules in XUV laser pulses. We trace the femtosecond nuclear-wave-packet dynamics in ionic states of oxygen and nitrogen diatomic molecules by comparing measured kinetic-energy-release spectra with classical and quantum-mechanical simulations. Experiments were done at the free-electron laser in Hamburg (FLASH) using 38-eV XUV-pump–XUV-probe. The summary and outlook of the work is discussed in Chapter 8.

The mechanisms of thermal transport in defect-free silicon nanostructures are examined using a combination of lattice dynamics (LD) calculations and the Boltzmann transport equation (BTE). To begin, the thermal conductivity reduction in thin films is examined using a hierarchical method that first predicts phonon transport properties using LD calculations, and then solves the phonon BTE using the lattice Boltzmann method. This approach, which considers all of the phonons in the first Brillouin-zone, is used to assess the suitability of common assumptions used to reduce the computational effort. Specifically, we assess the validity of: (i) neglecting the contributions of optical modes, (ii) the isotropic approximation, (iii) assuming an averaged bulk mean-free path (i.e., the Gray approximation), and (iv) using the Matthiessen rule to combine the effect of different scattering mechanisms. Because the frequency-dependent contributions to thermal conductivity change as the film thickness is reduced, assumptions that are valid for bulk are not necessarily valid for thin films. Using knowledge gained from this study, an analytical model for the length-dependence of thin film thermal conductivity is presented and compared to the predictions of the LD-based calculations. The model contains no fitting parameters and only requires the bulk lattice constant, bulk thermal conductivity, and an acoustic phonon speed as inputs. By including the mode-dependence of the phonon lifetimes resulting from phonon-phonon and phonon-boundary scattering, the model predictions capture the approach to the bulk thermal conductivity better than predictions made using Gray models based on a single lifetime. Both the model and the LD-based method are used to assess a procedure commonly used to extract bulk thermal conductivities from length-dependent molecular dynamics simulation data. Because the mode-dependence of thermal conductivity is not included in the derivation of this extrapolation procedure, using it can result in significant error. Finally, phonon transport across a silicon/vacuum-gap/silicon structure is modelled using lattice dynamics and Landauer theory. The phonons transmit thermal energy across the vacuum gap via atomic interactions between the leads. Because the incident phonons do not encounter a classically impenetrable potential barrier, this mechanism is not a tunneling phenomenon. The heat flux due to phonon transport can be 4 orders of magnitude larger than that due to photon transport predicted from near-field radiation theory.

The mechanisms of thermal transport in defect-free silicon nanostructures are examined using a combination of lattice dynamics (LD) calculations and the Boltzmann transport equation (BTE). To begin, the thermal conductivity reduction in thin films is examined using a hierarchical method that first predicts phonon transport properties using LD calculations, and then solves the phonon BTE using the lattice Boltzmann method. This approach, which considers all of the phonons in the first Brillouin-zone, is used to assess the suitability of common assumptions used to reduce the computational effort. Specifically, we assess the validity of: (i) neglecting the contributions of optical modes, (ii) the isotropic approximation, (iii) assuming an averaged bulk mean-free path (i.e., the Gray approximation), and (iv) using the Matthiessen rule to combine the effect of different scattering mechanisms. Because the frequency-dependent contributions to thermal conductivity change as the film thickness is reduced, assumptions that are valid for bulk are not necessarily valid for thin films. Using knowledge gained from this study, an analytical model for the length-dependence of thin film thermal conductivity is presented and compared to the predictions of the LD-based calculations. The model contains no fitting parameters and only requires the bulk lattice constant, bulk thermal conductivity, and an acoustic phonon speed as inputs. By including the mode-dependence of the phonon lifetimes resulting from phonon-phonon and phonon-boundary scattering, the model predictions capture the approach to the bulk thermal conductivity better than predictions made using Gray models based on a single lifetime. Both the model and the LD-based method are used to assess a procedure commonly used to extract bulk thermal conductivities from length-dependent molecular dynamics simulation data. Because the mode-dependence of thermal conductivity is not included in the derivation of this extrapolation procedure, using it can result in significant error. Finally, phonon transport across a silicon/vacuum-gap/silicon structure is modelled using lattice dynamics and Landauer theory. The phonons transmit thermal energy across the vacuum gap via atomic interactions between the leads. Because the incident phonons do not encounter a classically impenetrable potential barrier, this mechanism is not a tunneling phenomenon. The heat flux due to phonon transport can be 4 orders of magnitude larger than that due to photon transport predicted from near-field radiation theory.

Doctor of Philosophy / Department of Physics / B.D. DePaola / In recent years, coherent excitation techniques have focused on the ability to efficiently prepare atomic or molecular systems into a selected state. Such population control plays a key role in cutting-edge research taking place today, such as in the areas of quantum information and laser-controlled chemical reactions. Stimulated Raman adiabatic passage (STIRAP) is a widely-used coherent excitation technique that provides a relatively robust control mechanism for efficiently exciting a target population into a desired state. While the technique is well proven, current experimental techniques yield little information on the population dynamics taking place throughout the excitation process, and experimentalists rely solely on final excited-state measurements to determine the efficiency of population transfer. This dissertation presents a unique diagnostic tool to measure multilevel coherent population transfer on a short (nanosecond) timescale. The technique described here uses magneto-optical trap recoil ion momentum spectroscopy (MOTRIMS) as a noninvasive probe of a coherently-controlled system. It provides extremely detailed information about the excitation process, and highlights some important characteristics seen in excited populations that would otherwise be misleading or completely overlooked if one were to use more traditional diagnostic techniques. This dissertation discusses both the theoretical and experimental results applied to three-level coherently excited target populations of Rb-87.

Doctor of Philosophy / Department of Physics / Uwe Thumm / The time-dependent Schrödinger equation is integrated on a numerical lattice for up to three-dimensional problems. The wave packet propagation technique has been applied to ion – atom collisions in a strong laser field, the vibrational nuclear motion in small homonuclear diatomic molecular ions, and for the scattering of an ion in front of a metallic surface. For laser-assisted proton – hydrogen collisions it is shown, that strong circularly polarized radiation significantly alters the capture and ionization probabilities and results in a dichroism with respect to the helicity. In a pump – control – probe scheme, “stroboscopic” exposure of a nuclear wave packet of the deuterium molecular ion by a single or a series of short and intense laser control pulses may be used to produce an almost stationary distribution of a single vibrational level, where the nodal structure can be tested using the Coulomb explosion imaging technique. Using a pump – probe setup with variable probe delays it is proposed to use Fourier analysis of the time dependence of the Coulomb explosion kinetic energy release spectrum to reveal insight into the initial vibrational state distribution for small diatomic molecules. A last application demonstrates, that resonant charge transfer for scattering of a negative hydrogen anion on a metal surface depends crucially on the position of surface and image states relative to the conduction and valence band, thereby implying different reaction mechanisms for different surface cuts of a metal.

Doctor of Philosophy / Department of Physics / Brett D. Esry / The nuclear rotational period of the simplest molecule H[subscript]2[superscript]+ is about 550 fs, which is more than 35 times longer than its vibrational period of 15 fs. The rotational time scale is also much longer than widely available ultra short laser pulses which have 10 fs or less duration. The large difference in rotational period and ultra short laser pulse duration raises questions about the importance of nuclear rotation in theoretical studies of H[subscript]2[superscript]+ dissociation by these pulses. In most studies, reduced-dimensionality calculations are performed by freezing the molecular axis in one direction, referred to as the aligned model. We have systematically compared the aligned model with our full-dimensionality results for total dissociation probability and field-free dynamics of the dissociating fragments. The agreement between the two is only qualitative even for ultra short 10 fs pulses. Post-pulse dynamics of the bound wave function show rotational revivals. Significant alignment of H[subscript]2[superscript]+ occurs at these revivals. Our theoretical formulation to solve the time-dependent Schrodinger equation is an important step forward to make quantitative comparison between theory and experiment. We accurately calculate observables such as kinetic energy, angular, and momentum distributions. Reduced-dimensionality calculations cannot predict momentum distributions. Our theoretical approach presents the first momentum distribution of H[subscript]2[superscript]+ dissociation by few cycle laser pulses. These observables can be directly compared to the experiment. After taking into account averaging steps over the experimental conditions, we find remarkable agreement between the theory and experiment. Thus, our theoretical formulation can make predictions. In H[subscript]2[superscript]+ dissociation by pulses less than 10 fs, an asymmetry in the momentum distribution occurs by the interference of different pathways contributing to the same energy. The asymmetry, however, becomes negligible after averaging over experimental conditions. In a proposed pump-probe scheme, we predict an order of magnitude enhancement in the asymmetry and are optimistic that it can be observed.