Strong parallel magnetic field effects on the hydrogen molecular ion

Equilibrium distances, binding energies and dissociation energies for the ground and low-lying states of the hydrogen molecular ion in a strong magnetic field parallel to the internuclear axis are calculated and refined, by using the two- dimensional pseudospectral method. High-precision results are presented for the binding energies over a wider field regime than already given in the literature (Kravchenko and Liberman 1997 Phys. Rev. A 55 2701). The present work removes a long- standing discrepancy for the R-eq value in the 1sigma(u) state at a field strength of 1.0 x 10(6) T. The dissociation energies of the antibonding 1pi(g) state induced by magnetic fields are determined accurately. We have also observed that the antibonding 1pi(g) potential energy curve develops a minimum if the field is sufficiently strong. Some unreliable results in the literature are pointed out and discussed. A way to efficiently treat vibrational processes and coupling between the nuclear and the electronic motions in magnetic fields is also suggested within a three-dimensional pseudospectral scheme.

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.

The closed orbits and the photo-excitation scaled spectrum of the hydrogen atom in crossed fields

In a recent Letter to the Editor (J Rao, D Delande and K T Taylor 2001 J. Phys. B: At. Mol. Opt. Phys. 34 L391-9) we made a brief first report of our quantal and classical calculations for the hydrogen atom in crossed electric and magnetic fields at constant scaled energy and constant scaled electric field strength. A principal point of that communication was our statement that each and every peak in the Fourier transform of the scaled quantum photo-excitation spectrum for scaled energy value epsilon = -0.586 538 871028 43 and scaled electric value (f) over tilde = 0.068 537 846 207 618 71 could be identified with a scaled action value of a found and mapped-out closed orbit up to a scaled action of 20. In this follow-up paper, besides presenting full details of our quantum and classical methods, we set out the scaled action values of all 317 closed orbits involved, together with the geometries of many.

Single- and multiphoton detachment of K

Single- and multiphoton detachment rates have been calculated for K- using the R-matrix Floquet approach. Single-photon detachment rates, obtained at a laser field peak intensity of 10(9) W cm(-2), are discussed and compared with other theoretical work. Two-photon detachment rates at the same intensity have also been obtained, and similarities with results from earlier calculations for Li- and Na- are discussed. Three-photon rates are also presented at this laser intensity, and are compared and contrasted with those arising in the single-photon case, since both involve resonance structure with P-1(o) symmetry. The influence of resonances such as the 5s(2) S-1(e) doubly excited state and excitations of the residual atom are also considered.

Multi-ionization of helium and lithium using the independent electron and independent event models with intrinsic CDW

Cross sections for the multi-ionization of He and Li are presented for impact energies in the range of 50 to 1000 keV/amu. These are calculated using the eikonal initial state approximation to represent the input and exit channels of the active electrons. The ionization process is simulated in a variety of ways, most notably an attempt to account for the effects of electron correlation via the inclusion of a continuum density of states (CDS) term. Inadequacies, of the CDW formulation at small impact parameters, and of the models themselves, are discussed and conclusions are drawn on what repercussions this has for the cross sections calculated.

Total and single differential cross sections for simple resonant collisions using a fully orthonormal continuum-distorted-wave basis

Simple electron capture processes are studied using an orthonormal two state continuum-distorted-wave (CDW) basis. The suitability of the basis set is tested by comparing predictions for total and differential cross sections with available experimental data. Overall good agreement is obtained and the authors conclude that a relatively small CDW basis set may be suitable to model a wide variety of low-energy collisions if the members of this extended set are astutely chosen.

Electron-impact ionization of atomic hydrogen close to threshold

A systematic study of the ionization of atomic hydrogen by electron impact from 0.3 eV to a few eV above the ionization threshold has been carried out using a semiclassical-quantal calculation. Differential and integrated cross sections are presented at 0.3 eV above the energy threshold. Triple- differential cross sections (TDCS) are presented at constant theta(12) geometry where theta(12)=180degrees and 150degrees. Good agreement is achieved with the measurement [Roder, Phys. Rev. Lett. 79, 1666 (1997)] and calculations based on exterior complex scaling at 2 eV and 4 eV above threshold. Results of triple-differential cross sections are also presented at 0.3, 0.5, and 1.0 eV above threshold at both theta(12)=180degrees and 150degrees. At theta(12)=180degrees the small local maximum in the TDCS around theta(1)=90degrees reported by Pan and Starace [Phys. Rev. A 45, 4588 (1992)] at 0.5 eV above threshold is not observed in our calculation at energies down to 0.3 eV above threshold. The shape of our double differential cross sections seems to disagree qualitatively with the available calculations as we found two local maxima around 15degrees and 165degrees in our calculation. Single differential cross sections in our formulation appear naturally as a function of total excess energy E and, therefore, constant for all combinations of individual electron energies E-1 and E- 2 with E=E-1+E-2. Total ionization cross sections are also compared with measurement and available theoretical calculations and found to be in reasonably good agreement up to 10 eV above ionization threshold.

Angular distributions of electrons in photoionization of 3s subshell of atomic chlorine

R-matrix calculated photoelectron angular distribution asymmetry parameters, beta for Cl+ 3s3p(5) P-3(o) and 3s(2)3p(3) (D-2(o))3d P-1(o) final ionic states in photoionization of the ground state of atomic Cl are presented in the photon energy range from threshold to 80 eV. The results, characterized by prominent autoionization structures which are sensitive to multielectron correlations, are compared with those recently measured by Whitfield et al (Whitfield S B, Kehoe K, Krause M 0 and Caldwell C D 2000 Phys. Rev. Lett. 84 4818). Contrary to experiment and previous theoretical calculations, our detailed CIV3 structure calculation (Deb N C, Crothers D S F, Felfli Z and Msezane A Z 2002 J. Phys. B: At. Mol. Opt. Phys. submitted) has identified the lowest P-1(o) level of Cl+ as 3S(2)3p(3)(D-2(o))3d P-1(o) rather than 3s3p(5) P-1(o). The implications and consequences of the measured data for the 3s P-1(o) level are also discussed in the context of our calculated energies for Cl+ and beta for 3d P-1(o).

Restricted rotation dynamics and far-infrared spectra of liquid water governed by elastic interactions

A semi-phenomenological model describing wideband dielectric and far-infrared spectra of liquid water was proposed recently by the same authors [J. Mol. Struct. 606 (2002) 9], where a small dipole-moment component changing harmonically with time determines a weak absorption band (termed here the R-band) centred at the wavenumber v similar to 200 cm(-1). In the present work, a rough molecular theory of the R-band based on the concept of elastic interactions is given. Stretching and bending of hydrogen bonds cause restricted rotation (RR) of a polar water molecule in terms of a dimer comprising the H- bonded molecules. Analytical expression for the RR frequency nu(str) is derived as a function of the RR amplitude, geometrical parameters and force constants. The density g(nu(str)) of frequency distribution is shown to be centred in the R-band. The spectrum of the dipolar auto-correlation function calculated for this structural-dynamical model is found. A composite model comprising two intermolecular potentials is proposed, which yields for water a good description of the experimental wideband (from 0 to 1000 cm(- 1)) spectra of complex permittivity and of absorption coefficient. The presented interpretation of these spectra is based on a concept that water presents a two-component solution, with components differing by the types of molecular rotation. (C) 2003 Elsevier B.V. All rights reserved.

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.

Specific and unspecific interactions in polar fluids in view of wideband dielectric far-infrared spectra

A simple molecular analytical theory of dielectric relaxation in strongly polar fluids is considered in terms of a semi- phenomenological approach. Theoretical spectra epsilon(v), a(v) of complex permittivity and absorption coefficient are fully determined by a form of intermolecular potential well, in which a dipole reorients. In a recent publication by VI. Gaiduk, O.F. Nielsen, and T.S. Perova [J. Molliq 95 (1002) 1-25] the wideband spectra of liquid H2O and D2O were described in terms of a composite model comprising the rectangular and the cosine squared potential wells. Much better results are achieved in this work, where the rectangular well is replaced by a well with a rounded bottom termed the hat-curved well. The spectrum of the auto-correlation function (ACF) is calculated for such a potential. The proposed theory of a composite model, comprising hat-curved and parabolic wells, is applied for liquid water. This model is capable for describing the Debye relaxation region, the second relaxation region in the submillimeter wavelength range, and the far infra-red (FIR) e(v), a(v) spectra, where an intense librational band and an additional weak band are placed, respectively, near 700 cm(-1) and 200 cm(-1). The latter band reflects the features of so-called specific (viz. directly related to H-bonds) interactions and the former band reflects the features of unspecific interactions. The physical mechanisms connected with these types of interactions are discussed in terms of two relevant types of water structure (types of molecular rotation). The proposed theory is also applied to a non-associated liquid in terms of one hat-curved potential well. (C) 2004 Elsevier B.V. All rights reserved.

Shifting of the electron-capture-to-the-continuum peak in proton-helium collisions at 10 and 20 keV

A refined theoretical approach has been developed to study the double-differential cross sections (DDCS's) in proton-helium collisions as a function of the ratio of ionized electron velocity to the incident proton velocity. The refinement is done in the present coupled-channel calculation by introducing a continuum distorted wave in the final state coupled with discrete states including direct as well as charge transfer channels. It is confirmed that the electron-capture-to-the-continuum (ECC) peak is slightly shifted to a lower electron velocity than the equivelocity position. Comparing measurements and classical trajectory Monte Carlo (CTMC) calculations at 10 and 20 keV proton energies, excellent agreement of the ECC peak heights is achieved at both energies. However, a minor disagreement in the peak positions between the present calculation and the CTMC results is noted. A smooth behavior of the DDCS is found in the present calculation on both sides of the peak whereas the CTMC results show some oscillatory behavior particularly to the left of the peak, associated with the statistical nature of CTMC calculations.

Nonharmonic transverse vibration of the H-bonded molecules and the THz spectra in ice and water

A nonlinear equation of motion is found for the dimer comprising two charged H2O molecules. The THz dielectric response to nonharmonic vibration of a nonrigid dipole, forming the hydrogen bond (HB), is found in the direction transverse to this bond. An explicit expression is derived for the autocorrelator that governs the spectrum generated by transverse vibration (TV) of such a dipole. This expression is obtained by analytical solution of the truncated set of recurrence equations. The far infrared (FIR) spectra of ice at the temperature - 7 degrees C are calculated. The wideband, in the wavenumber (frequency) v range 0... 100.0 cm(-1), spectra are obtained for liquid water at room temperature and for supercooled water at -5.6 degrees C. All spectra are represented in terms of the complex permittivity epsilon(v) and the absorption coefficient alpha(v). The obtained analytical formula for epsilon comprises the term epsilon(perpendicular to) pertinent to the studied TV mechanism with three additional terms Delta epsilon(q), Delta epsilon(mu), and epsilon(or) arising, respectively, from: elastic harmonic vibration of charged molecules along the H-bond; elastic reorientation of HB permanent dipoles; and rather free libration of permanent dipoles in 'defects' of water/ice structure. The suggested TV-dielectric relaxation mechanism allows us: (a) to remove the THz 'deficit' of loss epsilon" inherent in previous theoretical studies; (b) to explain the THz loss and absorption spectra in supercooled (SC) water; and (c) to describe, in agreement with the experiment, the low- and high-frequency tails of the two bands of ice H2O located in the range 10...300 cm(-1). Specific THz dielectric properties of SC water are ascribed to association of water molecules, revealed in our study by transverse vibration of HB charged molecules. (C) 2006 Published by Elsevier B.V.

Density functional and Monte Carlo studies of sulfur. I. Structure and bonding in S-n rings and chains (n=2-18) (vol 118, pg 9257, 2003)

Density functional calculations have been performed for ring isomers of sulfur with up to 18 atoms, and for chains with up to ten atoms. There are many isomers of both types, and the calculations predict the existence of new forms. Larger rings and chains are very flexible, with numerous local energy minima. Apart from a small, but consistent overestimate in the bond lengths, the results reproduce experimental structures where known. Calculations are also performed on the energy surfaces of S8 rings, on the interaction between a pair of such rings, and the reaction between one S8 ring and the triplet diradical S8 chain. The results for potential energies, vibrational frequencies, and reaction mechanisms in sulfur rings and chains provide essential ingredients for Monte Carlo simulations of the liquid–liquid phase transition. The results of these simulations will be presented in Part II.

Density functional and Monte Carlo studies of sulfur. II. Equilibrium polymerization of the liquid phase

The equilibrium polymerization of sulfur is investigated by Monte Carlo simulations. The potential energy model is based on density functional results for the cohesive energy, structural, and vibrational properties as well as reactivity of sulfur rings and chains [Part I, J. Chem. Phys. 118, 9257 (2003)]. Liquid samples of 2048 atoms are simulated at temperatures 450less than or equal toTless than or equal to850 K and P=0 starting from monodisperse S-8 molecular compositions. Thermally activated bond breaking processes lead to an equilibrium population of unsaturated atoms that can change the local pattern of covalent bonds and allow the system to approach equilibrium. The concentration of unsaturated atoms and the kinetics of bond interchanges is determined by the energy DeltaE(b) required to break a covalent bond. Equilibrium with respect to the bond distribution is achieved for 15less than or equal toDeltaE(b)less than or equal to21 kcal/mol over a wide temperature range (Tgreater than or equal to450 K), within which polymerization occurs readily, with entropy from the bond distribution overcompensating the increase in enthalpy. There is a maximum in the polymerized fraction at temperature T-max that depends on DeltaE(b). This fraction decreases at higher temperature because broken bonds and short chains proliferate and, for Tless than or equal toT(max), because entropy is less important than enthalpy. The molecular size distribution is described well by a Zimm-Schulz function, plus an isolated peak for S-8. Large molecules are almost exclusively open chains. Rings tend to have fewer than 24 atoms, and only S-8 is present in significant concentrations at all T. The T dependence of the density and the dependence of polymerization fraction and degree on DeltaE(b) give estimates of the polymerization temperature T-f=450+/-20 K. (C) 2003 American Institute of Physics.