High energy three -body breakup of three -nucleon systems

The main goal of this dissertation was to study two- and three-nucleon Short Range Correlations (SRCs) in high energy three-body breakup of 3He nucleus in 3He(e, e'NN) N reaction. SRCs are characterized by quantum fluctuations in nuclei during which constituent nucleons partially overlap with each other. ^ A theoretical framework is developed within the Generalized Eikonal Approximation (GEA) which upgrades existing medium-energy methods that are inapplicable for high momentum and energy transfer reactions. High momentum and energy transfer is required to provide sufficient resolution for probing SRCs. GEA is a covariant theory which is formulated through the effective Feynman diagrammatic rules. It allows self-consistent calculation of single and double re-scatterings amplitudes which are present in three-body breakup processes. The calculations were carried out in detail and the analytical result for the differential cross section of 3He(e, e'NN)N reaction was derived in a form applicable for programming and numerical calculations. The corresponding computer code has been developed and the results of computation were compared to the published experimental data, showing satisfactory agreement for a wide range of values of missing momenta. ^ In addition to the high energy approximation this study exploited the exclusive nature of the process under investigation to gain more information about the SRCs. The description of the exclusive 3He( e, e'NN)N reaction has been done using the formalism of the nuclear decay function, which is a practically unexplored quantity and is related to the conventional spectral function through the integration of the phase space of the recoil nucleons. Detailed investigation showed that the decay function clearly exhibits the main features of two- and three-nucleon correlations. Four highly practical types of SRCs in 3He nucleus were discussed in great detail for different orders of the final state re-interactions using the decay function as an unique identifying tool. ^ The overall conclusion in this dissertation suggests that the investigation of the decay function opens up a completely new venue in studies of short range nuclear properties. ^

High energy three-body breakup of three-nucleon systems

The main goal of this dissertation was to study two- and three-nucleon Short Range Correlations (SRCs) in high energy three-body breakup of 3He nucleus in 3He(e, e'NN)N reaction. SRCs are characterized by quantum fluctuations in nuclei during which constituent nucleons partially overlap with each other. A theoretical framework is developed within the Generalized Eikonal Approximation (GEA) which upgrades existing medium-energy methods that are inapplicable for high momentum and energy transfer reactions. High momentum and energy transfer is required to provide sufficient resolution for probing SRCs. GEA is a covariant theory which is formulated through the effective Feynman diagrammatic rules. It allows self-consistent calculation of single and double re-scatterings amplitudes which are present in three-body breakup processes. The calculations were carried out in detail and the analytical result for the differential cross section of 3He(e, e'NN)Nreaction was derived in a form applicable for programming and numerical calculations. The corresponding computer code has been developed and the results of computation were compared to the published experimental data, showing satisfactory agreement for a wide range of values of missing momenta. In addition to the high energy approximation this study exploited the exclusive nature of the process under investigation to gain more information about the SRCs. The description of the exclusive 3He(e, e'NN)N reaction has been done using the formalism of the nuclear decay function, which is a practically unexplored quantity and is related to the conventional spectral function through the integration of the phase space of the recoil nucleons. Detailed investigation showed that the decay function clearly exhibits the main features of two- and three-nucleon correlations. Four highly practical types of SRCs in 3He nucleus were discussed in great detail for different orders of the final state re-interactions using the decay function as an unique identifying tool. The overall conclusion in this dissertation suggests that the investigation of the decay function opens up a completely new venue in studies of short range nuclear properties.

MULTI-BAND BOSE HUBBARD MODELS AND EFFECTIVE THREE-BODY INTERACTIONS

Experiments with ultracold atoms in optical lattice have become a versatile testing ground to study diverse quantum many-body Hamiltonians. A single-band Bose-Hubbard (BH) Hamiltonian was first proposed to describe these systems in 1998 and its associated quantum phase-transition was subsequently observed in 2002. Over the years, there has been a rapid progress in experimental realizations of more complex lattice geometries, leading to more exotic BH Hamiltonians with contributions from excited bands, and modified tunneling and interaction energies. There has also been interesting theoretical insights and experimental studies on “un- conventional” Bose-Einstein condensates in optical lattices and predictions of rich orbital physics in higher bands. In this thesis, I present our results on several multi- band BH models and emergent quantum phenomena. In particular, I study optical lattices with two local minima per unit cell and show that the low energy states of a multi-band BH Hamiltonian with only pairwise interactions is equivalent to an effec- tive single-band Hamiltonian with strong three-body interactions. I also propose a second method to create three-body interactions in ultracold gases of bosonic atoms in a optical lattice. In this case, this is achieved by a careful cancellation of two contributions in the pair-wise interaction between the atoms, one proportional to the zero-energy scattering length and a second proportional to the effective range. I subsequently study the physics of Bose-Einstein condensation in the second band of a double-well 2D lattice and show that the collision aided decay rate of the con- densate to the ground band is smaller than the tunneling rate between neighboring unit cells. Finally, I propose a numerical method using the discrete variable repre- sentation for constructing real-valued Wannier functions localized in a unit cell for optical lattices. The developed numerical method is general and can be applied to a wide array of optical lattice geometries in one, two or three dimensions.

Dynamic tunnelling ionization of H-+(2) in intense fields

Intense-field ionization of the hydrogen molecular ion by linearly polarized light is modelled by direct solution of the fixed-nuclei time-dependent Schrodinger equation and compared with recent experiments. Parallel transitions are calculated using algorithms which exploit massively parallel computers. We identify and calculate dynamic tunnelling ionization resonances that depend on laser wavelength and intensity, and molecular bond length. Results for lambda similar to 1064 nm are consistent with static tunnelling ionization. At shorter wavelengths lambda similar to 790 nm large dynamic corrections are observed. The results agree very well with recent experimental measurements of the ion spectra. Our results reproduce the single peak resonance and provide accurate ionization rate estimates at high intensities. At lower intensities our results confirm a double peak in the ionization rate as the bond length varies.

A discrete time-dependent method for metastable atoms and molecules in intense fields

The full-dimensional time-dependent Schrodinger equation for the electronic dynamics of single-electron systems in intense external fields is solved directly using a discrete method. Our approach combines the finite-difference and Lagrange mesh methods. The method is applied to calculate the quasienergies and ionization probabilities of atomic and molecular systems in intense static and dynamic electric fields. The gauge invariance and accuracy of the method is established. Applications to multiphoton ionization of positronium, the hydrogen atom and the hydrogen molecular ion are presented. At very high laser intensity, above the saturation threshold, we extend the method using a scaling technique to estimate the quasienergies of metastable states of the hydrogen molecular ion. The results are in good agreement with recent experiments. (C) 2004 American Institute of Physics.

Efficient grid treatment of the ionization dynamics of laser-driven H-2(+)

We implement a parallel, time-dependent hybrid finite-difference Lagrange mesh code to model the electron dynamics of the fixed-nuclei hydrogen molecular ion subjected to intense ultrashort laser Pulses, Ionization rates are calculated and compared with results from a previous finite-difference approach and also with published Floquet results. The sensitivity of the results to the gauge describing the electron-field interaction is studied. Visualizations of the evolving wave packets are also presented in which the formation of a stable bound-state resonance is observed.

Ionization dynamics of laser-driven H-2(+)

We set out aspects of a numerical algorithm used in solving the full-dimensionality time-dependent Schrodinger equation describing the electronic motion of the hydrogen molecular ion driven by an intense, linearly polarized laser pulse aligned along the molecular axis. This algorithm has been implemented within the fixed inter-nuclear separation approximation in a parallel computer code, a brief summary of which is given. Ionization rates are calculated and compared with results from other methods, notably the time-independent Floquet method. Our results compare very favourably with the precise predictions of the Floquet method, although there is some disagreement with other wavepacket calculations. Visualizations of the electron dynamics are also presented in which electron rescattering is observed.

Dissociation of H2+ from a short, intense, infrared laser pulse: proton emission spectra and pulse calibration

The proton energy spectrum from photodissociation of the hydrogen molecular ion by short intense pulses of infrared light is calculated. The time-dependent Schrödinger equation is discretized and integrated. For few-cycle pulses one can resolve vibrational structure, arising from the experimental preparation of the molecular ion. We calculate the corresponding energy spectrum and analyse the dependence on the pulse time delay, pulse length and intensity of the laser for ? ~ 790 nm. We conclude that the proton spectrum is a sensitive probe of both the vibrational populations and phases, and allows us to distinguish between adiabatic and nonadiabatic dissociation. Furthermore, the sensitivity of the proton spectrum from H2+ is a practical means of calibrating the pulse. Our results are compared with recent measurements of the proton spectrum for 65 fs pulses using a Ti:Sapphire laser (? ~ 790 nm) including molecular orientation and focal-volume averaging. Integrating over the laser focal volume, for the intensity I ~ 3 × 1015 W cm-2, we find our results are in excellent agreement with these experiments.

Many-body energy decomposition of hydrogen-bonded glycine clusters in gas-phase

A detailed analysis of the many-body contribution to the interaction energies of the gas-phase hydrogen-bonded glycine clusters, (Gly)(N), N = 1-4 is presented. The energetics of the hydrogen-bonded dimer, trimer and tetramer complexes have been analyzed using density-functional theory. The magnitude of the two-through four-body energy terms have been calculated and compared. The relaxation energy and the two-body energy terms are the principal contributors to the total binding energy. Four-body contribution is negligible. However, the three-body contribution is found to be sizable and the formation of the cyclic glycine trimer presents geometric strains that make it less favorable. (C) 2010 Elsevier B.V. All rights reserved.

Dynamical calculation of direct muon-transfer rates from thermalized muonic hydrogen to C6+ and O8+

Moun-transfer reactions from muonic hydrogen to carbon and oxygen nuclei employing a full quantum-mechanical few-body description of rearrangement scattering were studied by solving the Faddeev-Hahn-type equations using close-coupling approximation. The application of a close-coupling-type ansatz led to satisfactory results for direct muon-transfer reactions from muonic hydrogen to C6+ and O8+.

Applying many-body expansion aiming to make ab in initio molecular dynamics more efficient

In this present work we present a methodology that aims to apply the many-body expansion to decrease the computational cost of ab initio molecular dynamics, keeping acceptable accuracy on the results. We implemented this methodology in a program which we called ManBo. In the many-body expansion approach, we partitioned the total energy E of the system in contributions of one body, two bodies, three bodies, etc., until the contribution of the Nth body [1-3]: E = E1 + E2 + E3 + …EN. The E1 term is the sum of the internal energy of the molecules; the term E2 is the energy due to interaction between all pairs of molecules; E3 is the energy due to interaction between all trios of molecules; and so on. In Manbo we chose to truncate the expansion in the contribution of two or three bodies, both for the calculation of the energy and for the calculation of the atomic forces. In order to partially include the many-body interactions neglected when we truncate the expansion, we can include an electrostatic embedding in the electronic structure calculations, instead of considering the monomers, pairs and trios as isolated molecules in space. In simulations we made we chose to simulate water molecules, and use the Gaussian 09 as external program to calculate the atomic forces and energy of the system, as well as reference program for analyzing the accuracy of the results obtained with the ManBo. The results show that the use of the many-body expansion seems to be an interesting approach for reducing the still prohibitive computational cost of ab initio molecular dynamics. The errors introduced on atomic forces in applying such methodology are very small. The inclusion of an embedding electrostatic seems to be a good solution for improving the results with only a small increase in simulation time. As we increase the level of calculation, the simulation time of ManBo tends to largely decrease in relation to a conventional BOMD simulation of Gaussian, due to better scalability of the methodology presented. References [1] E. E. Dahlke and D. G. Truhlar; J. Chem. Theory Comput., 3, 46 (2007). [2] E. E. Dahlke and D. G. Truhlar; J. Chem. Theory Comput., 4, 1 (2008). [3] R. Rivelino, P. Chaudhuri and S. Canuto; J. Chem. Phys., 118, 10593 (2003).

Tailoring of ion species composition in complex plasmas with charge exchange collisions

A generic approach towards tailoring of ion species composition in reactive plasmas used for nanofabrication of various functional nanofilms and nanoassemblies, based on a simplified model of a parallel-plate rf discharge, is proposed. The model includes an idealized reactive plasma containing two neutral and two ionic species interacting via charge exchange collisions in the presence of a microdispersed solid component. It is shown that the number densities of the desired ionic species can be efficiently managed by adjusting the dilution of the working gas in a buffer gas, rates of electron impact ionization, losses of plasma species on the discharge walls, and surfaces of fine particles, charge exchange rates, and efficiency of three-body recombination processes in the plasma bulk. The results are relevant to the plasma-aided nanomanufacturing of ordered patterns of carbon nanotip and nanopyramid microemitters.

The 310 Helical Conformation of the Amino Terminal Decapeptide of Suzukacillin. 270 MHz 'H NMR Evidence for Eight Intramolecular Hydrogen Bonds

The amino terminal suzukacillin decapeptide fragment, Boc-Aib-Pro-Val-Aib-Val-Ala-Aib-Ala-Aib-Aitbh-eO Me, two pentapeptides Boc-AibPrc-Val-AibVal-OMe and Boc-Ala-AibAla-AibAibOMe, and the tripeptide Boc-Ala-AibAibOMe have been studied by 270-MHz 'H NMR spectroscopy. By use of solvent dependence of chemical shifts in a CDC13-(CD3),S0 system and temperature dependence of amide NH chemical shifts in (CD3),S0, the intramolecularly hydrogen bonded NH groups in these peptides have been identified. The tripeptide possesses one hydrogen bond, both pentapeptides show evidence for three intramolecular hydrogen bonds, and the decapeptide has eight NH groups participating in hydrogen bonding. An Ala( 1)-Aib(2) @ turn is proposed for the tripeptide. Both pentapeptides favor 310 helical conformations composed of three consecutive B turns. The decapeptide adopts a 310 helical conformation with some flexibility at the Va1(5)-Ala(6) segment. The proposed conformations are consistent with the known stereochemical preferences of Aib residues.

Charge Density Analysis of Ferulic Acid: Robustness of a Trifurcated C-H center dot center dot center dot O Hydrogen Bond

The toplogical features of a sporadic trifurcated C-H center dot center dot center dot O interaction region, where an oxygen atom acts as an acceptor of three weak hydrogen bonds, has been investigated by experimental and theoretical charge density analysis of ferulic acid. The interaction energy of the asymmetric molecular dimer formed by the trifurcated C-H center dot center dot center dot O motif, based on the multipolar model, is shown to be greater than the corresponding asymmetric O-H center dot center dot center dot O dimer in this crystal structure. Further, the hydrogen bond energies associated with these interaction motifs have been estimated from the local kinetic and potential energy densities at the bond critical points. The trends suggest that the interaction energy of the trifurcated C-H center dot center dot center dot O region is comparable to that of a single O-H center dot center dot center dot O hydrogen bond.