Electronic states and spin-forbidden cooling transitions of AlH and AlF

The feasibility of laser cooling AlH and AlF is investigated using ab initio quantum chemistry. All the electronic states corresponding to the ground and lowest two excited states of the Al atom are calculated using multi-reference configuration interaction (MRCI) and the large AV6Z basis set for AlH. The smaller AVQZ basis set is used to calculate the valence electronic states of AlF. Theoretical Franck-Condon factors are determined for the A(1)Pi -> X(1)Sigma(+) transitions in both radicals and found to agree with the highly diagonal factors found experimentally, suggesting computational chemistry is an effective method for screening suitable laser cooling candidates. AlH does not appear to have a transition quite as diagonal as that in SrF (which has been laser cooled) but the A(1)Pi -> X(1)Sigma(+) transition transition of AlF is a strong candidate for cooling with just a single laser, though the cooling frequency is deep in the UV. Furthermore, the a (3)Pi -> X(1)Sigma(+) transitions are also strongly diagonal and in AlF is a practical method for obtaining very low final temperatures around 3 mu K.

Prospects for ultracold carbon via charge exchange reactions and laser cooled carbides

Strategies to produce an ultracold sample of carbon atoms are explored and assessed with the help of quantum chemistry. After a brief discussion of the experimental difficulties using conventional methods, two strategies are investigated. The first attempts to exploit charge exchange reactions between ultracold metal atoms and sympathetically cooled C+ ions. Ab initio calculations including electron correlation have been conducted on the molecular ions [LiC]+ and [BeC]+ to determine whether alkali or alkaline earth metals are a suitable buffer gas for the formation of C atoms but strong spontaneous radiative charge exchange ensure they are not ideal. The second technique involves the stimulated production of ultracold C atoms from a gas of laser cooled carbides. Calculations on LiC suggest that the alkali carbides are not suitable but the CH radical is a possible laser cooling candidate thanks to very favourable Frank-Condon factors. A scheme based on a four pulse STIRAP excitation pathway to a Feshbach resonance is outlined for the production of atomic fragments with near zero centre of mass velocity.

Ultracold radiative charge transfer in hybrid ion-atom traps

We present ab initio quantum chemistry calculations for elastic scattering and the radiative charge transfer reaction process and collision rates for trapped ytterbium ions immersed in a quantum degenerate rubidium vapor.
The collision of the ion (or ions) with the quasiatom is the key mechanism to transfer quantum coherences between the systems. We use first-principles
quantum chemistry codes to obtain the potential surfaces and coupling terms for the two-body interaction of Yb^+ with Rb. We find that the low energy collision has an inelastic radiative charge transfer process in agreement with recent experiments.
The charge transfer cross section agrees well with the semiclassical Langevin model at higher energies but is dominated by resonances at submillikelvin temperatures.

Strong electronic correlation in the hydrogen chain: A variational Monte Carlo study

In this paper, we report a fully ab initio variational Monte Carlo study of the linear and periodic chain of hydrogen atoms, a prototype system providing the simplest example of strong electronic correlation in low dimensions. In particular, we prove that numerical accuracy comparable to that of benchmark density-matrix renormalization-group calculations can be achieved by using a highly correlated Jastrow-antisymmetrized geminal power variational wave function. Furthermore, by using the so-called "modern theory of polarization" and by studying the spin-spin and dimer-dimer correlations functions, we have characterized in detail the crossover between the weakly and strongly correlated regimes of this atomic chain. Our results show that variational Monte Carlo provides an accurate and flexible alternative to highly correlated methods of quantum chemistry which, at variance with these methods, can be also applied to a strongly correlated solid in low dimensions close to a crossover or a phase transition.

Universal tight binding model for chemical reactions in solution and at surfaces. II. Water

A revised water model intended for use in condensed phase simulations in the framework of the self consistent polarizable ion tight binding theory is constructed. The model is applied to water monomer, dimer, hexamers, ice, and liquid, where it demonstrates good agreement with theoretical results obtained by more accurate methods, such as DFT and CCSD(T), and with experiment. In particular, the temperature dependence of the self diffusion coefficient in liquid water predicted by the model, closely reproduces experimental curves in the temperature interval between 230 K and 350 K. In addition, and in contrast to standard DFT, the model properly orders the relative densities of liquid water and ice. A notable, but inevitable, shortcoming of the model is underestimation of the static dielectric constant by a factor of two. We demonstrate that the description of inter and intramolecular forces embodied in the tight binding approximation in quantum mechanics leads to a number of valuable insights which can be missing from ab initio quantum chemistry and classical force fields. These include a discussion of the origin of the enhanced molecular electric dipole moment in the condensed phases, and a detailed explanation for the increase of coordination number in liquid water as a function of temperature and compared with ice-leading to insights into the anomalous expansion on freezing. The theory holds out the prospect of an understanding of the currently unexplained density maximum of water near the freezing point. 

Formation and characterization of carbon-radical precursors in char steam gasification

Highly reactive radicals play an important role in high-temperature gasification processes. However, the effect of radicals on gasification has not been systematically investigated. In the present study, the formation of carbon-radical precursors using atomic radicals such as OH, O, and H and molecules such as H₂ and O₂ was characterized, and the effect of the precursors on the adsorption step of steam char gasification was studied using quantum chemistry methods. The results revealed that the radicals can be chemisorbed exothermically on char active sites, and the following order of reactivity was observed: O > H₂ > H > OH > O ₂. Moreover, hydrogen bonds are formed between steam molecules and carbon-radical complexes. Steam molecule adsorption onto carbon-O and carbon-OH complexes is easier than adsorption onto clean carbon surfaces. Alternatively, adsorption on carbon-O₂, carbon-H₂, and carbon-H complexes is at the same level with that of clean carbon surfaces; thus, OH and O radicals accelerate the physical adsorption of steam onto the char surface, H radical and O₂ and H₂ molecules do not have a significant effect on adsorption. © 2010 American Chemical Society.

Radiative charge transfer in cold and ultracold sulphur atoms colliding with protons

Radiative decay processes at cold and ultra cold temperatures for sulfur atoms colliding with protons are investigated. The MOLPRO quantum chemistry suite of codes was used to obtain accurate potential energies and transition dipole moments, as a function of internuclear distance, between low-lying states of the SH+ molecular cation. A multi-reference configuration-interaction approximation together with the Davidson correction is used to determine the potential energy curves and transition dipole moments, between the states of interest, where the molecular orbitals are obtained from state-averaged multi-configuration-self-consistent field calculations. The collision problem is solved approximately using an optical potential method to obtain radiative loss, and a fully two-channel quantum approach for radiative charge transfer. Cross sections and rate coefficients are determined for the first time for temperatures ranging from 10 μK up to 10 000 K. Results are obtained for all isotopes of sulfur, colliding with H+ and D+ ions and comparison is made to a number of other collision systems.

Theoretical Investigation of NH3-SCR Processes over Zeolites: A Review

The selective catalytic reduction (SCR) of NOx compounds with NH3 is a hot topic in recent years. Among various catalysts, zeolites are proved to be efficient and promising for NH3-SCR, yet the whole processes and intrinsic mechanism are still not well understood due to the structural complexity of zeolites. With the improvement of theoretical chemistry techniques, quantum-chemical calculations are now capable of modeling the structure, acidity, adsorption, and ultimately reaction pathways over zeolites to some extent. In this review, a brief summary of relevant concepts of NH3-SCR is presented. Cluster approaches, embedded techniques, and periodic treatments are described as three main methods. Details of quantum-chemical investigations toward the key issues such as, the structure of active sites, the adsorption of small molecules, and the reaction mechanism of NH3-SCR over zeolites are discussed. Finally, a perspective for future theoretical research is given. 

Production of ultracold hydrogen and deuterium via Doppler-cooled Feshbach molecules

A counterintuitive scheme to produce ultracold hydrogen via fragmentation of laser cooled diatomic hydrides is presented where the final atomic H temperature is inversely proportional to the mass of the molecular parent. In addition, the critical density for formation of a Bose-Einstein condensate (BEC) at a fixed temperature is reduced by a factor (mH/mMH)3/2 over directly cooled hydrogen atoms. The narrow Feshbach resonances between a S01 atom and hydrogen are well suited to a tiny center of mass energy release necessary during fragmentation. With the support of ab initio quantum chemistry, it is demonstrated that BaH is an ideal diatomic precursor that can be laser cooled to a Doppler temperature of ∼26μK with just two rovibronic transitions, the simplest molecular cooling scheme identified to date. Preparation of a hydrogen atom gas below the critical BEC temperature Tc is feasible with present cooling technology, with optical pulse control of the condensation process.

Quantum manifestations of scattering orbits in the scaled spectrum of a non-hydrogenic atom in crossed fields

Evidence for scattering closed orbits for the Rydberg electron of the singly excited helium atom in crossed electric and magnetic fields at constant scaled energy and constant scaled electric field strength has been found through a quantum calculation of the photo-excitation spectrum. A particular 3D scattering orbit in a mixed regular and chaotic region has been investigated and the hydrogenic 3D closed orbits composing it identified. To the best of our knowledge, this letter reports the first quantum calculation of the scaled spectrum of a non- hydrogenic atom in crossed fields.

Analytical calculation of far infrared spectra of ice in terms of a molecular model

Wideband far infrared (FIR) spectra of complex permittivity e(p) of ice are calculated in terms of a simple analytical theory based on the method of dipolar autocorrelation functions. The molecular model represents a revision of the model recently presented for liquid water in Adv. Chem. Phys. 127 (2003) 65. A composite two-fractional model is proposed. The model is characterised by three phenomenological potential wells corresponding to the three FIR bands observed in ice. The first fraction comprises dipoles reorienting in a rather narrow and deep hat-like well; these dipoles generate the librational band centred at the frequency approximate to 880 cm(-1). The second fraction comprises elastically interacting particles; they generate two nearby bands placed around frequency 200 cm(-1). For description of one of these bands the harmonic oscillator (HO) model is used, in which translational oscillations of two charged molecules along the H-bond are considered. The other band is produced by the H-bond stretch, which governs hindered rotation of a rigid dipole. Such a motion and its dielectric response are described in terms of a new cut parabolic (CP) model applicable for any vibration amplitude. The composite hat-HO-CP model results in a smooth epsilon(nu) ice spectrum, which does not resemble the noise-like spectra of ice met in the known literature. The proposed theory satisfactorily agrees with the experimental ice spectrum measured at - 7 degrees C. The calculated longitudinal optic-transverse optic (LO-TO) splitting occurring at approximate to 250 cm(-1) qualitatively agrees with the measured data. (c) 2004 Elsevier B.V. All rights reserved.

Ion-Atom/Argon - Calculation of ionization cross sections by fast ion impact for neutral target atoms ranging from hydrogen to argon

A FORTRAN 90 program is presented which calculates the total cross sections, and the electron energy spectra of the singly and doubly differential cross sections for the single target ionization of neutral atoms ranging from hydrogen up to and including argon. The code is applicable for the case of both high and low Z projectile impact in fast ion-atom collisions. The theoretical models provided for the program user are based on two quantum mechanical approximations which have proved to be very successful in the study of ionization in ion-atom collisions. These are the continuum-distorted-wave (CDW) and continuum-distorted-wave eikonal-initial-state (CDW-EIS) approximations. The codes presented here extend previously published. codes for single ionization of. target hydrogen [Crothers and McCartney, Comput. Phys. Commun. 72 (1992) 288], target helium [Nesbitt, O'Rourke and Crothers, Comput. Phys. Commun. 114 (1998) 385] and target atoms ranging from lithium to neon [O'Rourke, McSherry and Crothers, Comput. Phys. Commun. 131 (2000) 129]. Cross sections for all of these target atoms may be obtained as limiting cases from the present code. Title of program: ARGON Catalogue identifier: ADSE Program summary URL: http://cpc.cs.qub.ac.uk/cpc/summaries/ADSE Program obtainable from: CPC Program Library Queen's University of Belfast, N. Ireland Licensing provisions: none Computer for which the program is designed and others on which it is operable: Computers: Four by 200 MHz Pro Pentium Linux server, DEC Alpha 21164; Four by 400 MHz Pentium 2 Xeon 450 Linux server, IBM SP2 and SUN Enterprise 3500 Installations: Queen's University, Belfast Operating systems under which the program has been tested: Red-hat Linux 5.2, Digital UNIX Version 4.0d, AIX, Solaris SunOS 5.7 Compilers: PGI workstations, DEC CAMPUS Programming language used: FORTRAN 90 with MPI directives No. of bits in a word: 64, except on Linux servers 32 Number of processors used: any number Has the code been vectorized or parallelized? Parallelized using MPI No. of bytes in distributed program, including test data, etc.: 32 189 Distribution format: tar gzip file Keywords: Single ionization, cross sections, continuum-distorted-wave model, continuum- distorted-wave eikonal-initial-state model, target atoms, wave treatment Nature of physical problem: The code calculates total, and differential cross sections for the single ionization of target atoms ranging from hydrogen up to and including argon by both light and heavy ion impact. Method of solution: ARGON allows the user to calculate the cross sections using either the CDW or CDW-EIS [J. Phys. B 16 (1983) 3229] models within the wave treatment. Restrictions on the complexity of the program: Both the CDW and CDW-EIS models are two-state perturbative approximations. Typical running time: Times vary according to input data and number of processors. For one processor the test input data for double differential cross sections (40 points) took less than one second, whereas the test input for total cross sections (20 points) took 32 minutes. Unusual features of the program: none (C) 2003 Elsevier B.V All rights reserved.

Determining the electronic structure and chemical potentials of molecules in solution

A simulation scheme is proposed for determining the excess chemical potential of a substance in solution. First, a Monte Carlo simulation is performed with classical models for solute and solvent molecules. A representative sample of these configurations is then used in a hybrid quantum/classical (QM/MM) calculation, where the solute is treated quantum-mechanically, and the average electronic structure is used to construct an improved classical model. This procedure is iterated to self-consistency in the classical model, which in practice is attained in one or two steps, depending on the quality of the initial guess. The excess free energy of the molecule within the QM/MM approach is determined relative to the classical model using thermodynamic perturbation theory with a cumulant expansion. The procedure provides a method of constructing classical point charge models appropriate for the solution and gives a measure of the importance of solvent fluctuations.

Classical computation with quantum systems

As semiconductor electronic devices scale to the nanometer range and quantum structures (molecules, fullerenes, quantum dots, nanotubes) are investigated for use in information processing and storage, it, becomes useful to explore the limits imposed by quantum mechanics on classical computing. To formulate the problem of a quantum mechanical description of classical computing, electronic device and logic gates are described as quantum sub-systems with inputs treated as boundary conditions, outputs expressed.is operator expectation values, and transfer characteristics and logic operations expressed through the sub-system Hamiltonian. with constraints appropriate to the boundary conditions. This approach, naturally, leads to a description of the subsystem.,, in terms of density matrices. Application of the maximum entropy principle subject to the boundary conditions (inputs) allows for the determination of the density matrix (logic operation), and for calculation of expectation values of operators over a finite region (outputs). The method allows for in analysis of the static properties of quantum sub-systems.