A justification for a nonadiabatic surface hopping Herman-Kluk semiclassical initial value representation of the time evolution operator

A justification is given for the validity of a nonadiabatic surface hopping Herman-Kluk (HK) semiclassical initial value representation (SC-IVR) method. The method is based on a propagator that combines the single surface HK SC-IVR method [J. Chem. Phys. 84, 326 (1986)] and Herman's nonadiabatic semiclassical surface hopping theory [J. Chem. Phys. 103, 8081 (1995)], which was originally developed using the primitive semiclassical Van Vleck propagator. We show that the nonadiabatic HK SC-IVR propagator satisfies the time-dependent Schrodinger equation to the first order of variant Planck's over 2pi and the error is O(variant Planck's over 2pi(2)). As a required lemma, we show that the stationary phase approximation, under current assumptions, has an error term variant Planck's over 2pi(1) order higher than the leading term. Our derivation suggests some changes to the previous development, and it is shown that the numerical accuracy in applications to Tully's three model systems in low energies is improved.     

Methane-water cross second virial coefficient with quantum corrections from an ab initio potential

We present our calculations of the cross second virial coefficient (B12) and of a related quantity, phi 12 = B12-TdB12/dT, for the methane-water system in the temperature range T = 200-1000 K. These calculations were performed using one of the ab initio potentials developed in previous work [Akin-Ojo and Szalewicz, J. Chem. Phys. 123, 134311 (2005)]. Quantum corrections of order variant Planck's over 2pi(2) were added to the computed classical values. We have estimated the uncertainties in our computed B12 and phi 12(T). This allowed evaluation of the quality of the experimental data to which we compare our results. We also provide an analytical expression for B12(T) as a function of the temperature T obtained by fitting the calculated values. This formula also predicts values of phi12(T) consistent with the directly calculated values.     

Speed dependent rotational angular momentum polarization of the O2 (a 1Deltag) fragment following ozone photolysis in the wavelength range 248-265 nm

The translational anisotropy and rotational angular momentum polarization of a selection of rotational states of the O2 (a 1Deltag; v=0) photofragment formed from ozone photolysis at 248, 260, and 265 nm have been determined using the technique of resonance enhanced multiphoton ionization in combination with time of flight mass spectrometry. At 248 nm, the dissociation is well described as impulsive in nature with all rotational states exhibiting similarly large, near-limiting values for the bipolar moments describing their angular momentum alignment and orientation. At 265 nm, however, the angular momentum polarization parameters determined for consecutive odd and even rotational states exhibit clear differences. Studies at the intermediate wavelength of 260 nm strongly suggest that such a difference in the angular momentum polarization is speed dependent and this proposal is consistent with the angular momentum polarization parameters extracted and reported previously for longer photolysis wavelengths [G. Hancock et al., Phys. Chem. Chem. Phys. 5, 5386 (2003); S. J. Horrocks et al., J. Chem. Phys. 126, 044308 (2007)]. The alternation of angular momentum polarization for successive odd and even J states may be a consequence of the different mechanisms leading to the formation of the two O2 (a 1Deltag) Lambda doublets. Specifically, the involvement of out of plane parent rotational motion is proposed as the origin for the observed depolarization for the Delta- relative to the Delta+ state.     

On the role of dynamical barriers in barrierless reactions at low energies: S(1D) + H2

Reaction probabilities as a function of total angular momentum (opacity functions) and the resulting reaction cross sections for the collision of open shell S((1)D) atoms with para-hydrogen have been calculated in the kinetic energy range 0.09-10 meV (1-120 K). The quantum mechanical hyperspherical reactive scattering method and quasi-classical trajectory and statistical quasi-classical trajectory approaches were used. Two different ab initio potential energy surfaces (PESs) have been considered. The widely used reproducing kernel Hilbert space (RKHS) PES by Ho et al. [T.-S. Ho, T. Hollebeek, H. Rabitz, S. D. Chao, R. T. Skodje, A. S. Zyubin, and A. M. Mebel, J. Chem. Phys 116, 4124 (2002)] and the recently published accurate double many-body expansion (DMBE)/complete basis set (CBS) PES by Song and Varandas [Y. Z. Song and A. J. C. Varandas, J. Chem. Phys. 130, 134317 (2009)]. The calculations at low collision energies reveal very different dynamical behaviors on the two PESs. The reactivity on the RKHS PES is found to be considerably larger than that on the DMBE/CBS PES as a result of larger reaction probabilities at low total (here also orbital) angular momentum values and to opacity functions which extend to significantly larger total angular momentum values. The observed differences have their origin in two major distinct topographic features. Although both PESs are essentially barrierless for equilibrium H-H distances, when the H-H bond is compressed the DMBE/CBS PES gives rise to a dynamical barrier which limits the reactivity of the system. This barrier is completely absent in the RHKS PES. In addition, the latter PES exhibits a van der Walls well in the entrance channel which reduces the height of the centrifugal barrier and is able to support resonances. As a result, a significant larger cross section is found on this PES, with marked oscillations attributable to shape resonances and/or to the opening of partial wave contributions. The comparison of the results on both PESs is illustrative of the wealth of the dynamics at low collision energy. It is also illuminating about the difficulties encountered in modeling an all-purpose global potential energy surface.     

Torsion-vibration coupling in methanol: the adiabatic approximation and intramolecular vibrational redistribution scaling

The four-dimensional model Hamiltonian of Wang and Perry [J. Chem. Phys. 109, 10795 (1998)] is used to compare the approximate adiabatic separation of the torsion and CH stretches in methanol to an exact solution of the same Hamiltonian. The adiabatic approximation accounts for the pattern of the energy levels in the lowest torsional states, including the inverted tunneling splittings, but does not account for the pattern of systematic two- and four-fold near degeneracies at high torsional excitation. In the adiabatic basis, the nonadiabatic couplings mix the torsional and vibrational degrees of freedom and hence are a source for intramolecular vibrational redistribution (IVR). These IVR matrix elements are found to decrease by only a factor of 2 or 3 with each higher coupling order, in agreement with the results of Pearman and Gruebele [Z. Phys. Chem. Munich 214, 1439 (2000)]. This gentle scaling behavior, which contrasts with a steeper falloff with coupling order in more rigid molecules, points to a more important role for direct high-order couplings in torsional molecules. In this model, the scaling behavior derives from a single coupling term that is low order in the torsional angular momentum in combination with one-dimensional torsional functions that include contributions from many torsional angular momenta.     

Quantum theory of molecular collisions in a magnetic field: efficient calculations based on the total angular momentum representation

An efficient method is presented for rigorous quantum calculations of atom-molecule and molecule-molecule collisions in a magnetic field. The method is based on the expansion of the wave function of the collision complex in basis functions with well-defined total angular momentum in the body-fixed coordinate frame. We outline the general theory of the method for collisions of diatomic molecules in the (2)Σ and (3)Σ electronic states with structureless atoms and with unlike (2)Σ and (3)Σ molecules. The cross sections for elastic scattering and Zeeman relaxation in low-temperature collisions of CaH((2)Σ(+)) and NH((3)Σ(-)) molecules with (3)He atoms converge quickly with respect to the number of total angular momentum states included in the basis set, leading to a dramatic (>10-fold) enhancement in computational efficiency compared to the previously used methods [A. Volpi and J. L. Bohn, Phys. Rev. A 65, 052712 (2002); R. V. Krems and A. Dalgarno, J. Chem. Phys. 120, 2296 (2004)]. Our approach is thus well suited for theoretical studies of strongly anisotropic molecular collisions in the presence of external electromagnetic fields.     

Resolutions of the Coulomb operator. VI. Computation of auxiliary integrals

We discuss the efficient computation of the auxiliary integrals that arise when resolutions of two-electron operators (specifically, the Coulomb operator [T. Limpanuparb, A. T. B. Gilbert, and P. M. W. Gill, J. Chem. Theory Comput. 7, 830 (2011)] and the long-range Ewald operator [T. Limpanuparb and P. M. W. Gill, J. Chem. Theory Comput. 7, 2353 (2011)]) are employed in quantum chemical calculations. We derive a recurrence relation that facilitates the generation of auxiliary integrals for Gaussian basis functions of arbitrary angular momentum and propose a near-optimal algorithm for its use.     

Inelastic scattering matrix elements for the nonadiabatic collision B(2P1/2)+H2(1Sigmag+,j)<-->B(2P3/2)+H2(1Sigmag+,j')

Inelastic scattering matrix elements for the nonadiabatic collision B(2P1/2)+H2(1Sigmag+,j)<-->B(2P3/2)+H2(1Sigmag+,j') are calculated using the time dependent channel packet method (CPM). The calculation employs 1 2A', 2 2A', and 1 2A" adiabatic electronic potential energy surfaces determined by numerical computation at the multireference configuration-interaction level [M. H. Alexander, J. Chem. Phys. 99, 6041 (1993)]. The 1 2A' and 2 2A', adiabatic electronic potential energy surfaces are transformed to yield diabatic electronic potential energy surfaces that, when combined with the total B+H2 rotational kinetic energy, yield a set of effective potential energy surfaces [M. H. Alexander et al., J. Chem. Phys. 103, 7956 (1995)]. Within the framework of the CPM, the number of effective potential energy surfaces used for the scattering matrix calculation is then determined by the size of the angular momentum basis used as a representation. Twenty basis vectors are employed for these calculations, and the corresponding effective potential energy surfaces are identified in the asymptotic limit by the H2 rotor quantum numbers j=0, 2, 4, 6 and B electronic states 2Pja, ja=1/2, 3/2. Scattering matrix elements are obtained from the Fourier transform of the correlation function between channel packets evolving in time on these effective potential energy surfaces. For these calculations the H2 bond length is constrained to a constant value of req=1.402 a.u. and state to state scattering matrix elements corresponding to a total angular momentum of J=1/2 are discussed for j=0<-->j'=0,2,4 and 2P1/2<-->2P1/2, 2P3/2 over a range of total energy between 0.0 and 0.01 a.u.     

High resolution electron momentum spectroscopy of dichlorodifluoromethane: unambiguous assignments of outer valence molecular orbitals

On account of controversial orbital assignment that appeared in previous works, [J. Chem. Phys. 120, 7933 (2004), and references therein] high resolution electron momentum spectroscopy (EMS) measurements on dichlorodifluoromethane has been carried out using a newly developed high resolution energy-momentum dispersive multichannel spectrometer employing asymmetric noncoplanar geometry at an impact energy of 2500 eV plus binding energy. Four resolved structures and two shoulders were obviously observed in high resolution binding energy spectrum in energy range covering eight outermost valence orbitals, whereas only two broad lobes were resolved in previous EMS studies [J. Chem. Phys. 120, 7933 (2004); Chin. Phys. 14, 2467 (2005)]. The ordering of these orbitals was reassigned unambiguously by simple comparison of experimental momentum distributions with theoretical ones.     

Study of the molecular structure, ionization spectrum, and electronic wave function of 1,3-butadiene using electron momentum spectroscopy and benchmark Dyson orbital theories

The scope of the present work is to reconcile electron momentum spectroscopy with elementary thermodynamics, and refute conclusions drawn by Saha et al. in J. Chem. Phys. 123, 124315 (2005) regarding fingerprints of the gauche conformational isomer of 1,3-butadiene in electron momentum distributions that were experimentally inferred from gas phase (e,2e) measurements on this compound [M. J. Brunger et al., J. Chem. Phys. 108, 1859 (1998)]. Our analysis is based on thorough calculations of one-electron and shake-up ionization spectra employing one-particle Green's function theory along with the benchmark third-order algebraic diagrammatic construction [ADC(3)] scheme. Accurate spherically averaged electron momentum distributions are correspondingly computed from the related Dyson orbitals. The ionization spectra and Dyson orbital momentum distributions that were computed for the trans-conformer of 1,3-butadiene alone are amply sufficient to quantitatively unravel the shape of all available experimental (e,2e) electron momentum distributions. A comparison of theoretical ADC(3) spectra for the s-trans and gauche energy minima with inner- and outer-valence high-resolution photoelectron measurements employing a synchrotron radiation beam [D. M. P. Holland et al., J. Phys. B 29, 3091 (1996)] demonstrates that the gauche structure is incompatible with ionization experiments in high-vacuum conditions and at standard temperatures. On the other hand, outer-valence Green's function calculations on the s-trans energy minimum form and approaching basis set completeness provide highly quantitative insights, within approximately 0.2 eV accuracy, into the available experimental one-electron ionization energies. At last, analysis of the angular dependence of relative (e,2e) ionization intensities nicely confirms the presence of one rather intense pi(-2) pi(*+1) satellite at approximately 13.1 eV in the ionization spectrum of the s-trans conformer.     

Recovering the Crooks equation for dynamical systems in the isothermal-isobaric ensemble: a strategy based on the equations of motion

The Crooks equation [Eq. (10) in J. Stat. Phys. 90, 1481 (1998)] relates the work done on a system during a nonequilibrium transformation to the free energy difference between the final and the initial state of the transformation. Recently, the authors have derived the Crooks equation for systems in the canonical ensemble thermostatted by the Nose-Hoover or Nose-Hoover chain method [P. Procacci et al., J. Chem. Phys. 125, 164101 (2006)]. That proof is essentially based on the fluctuation theorem by Evans and Searles [Adv. Phys. 51, 1529 (2002)] and on the equations of motion. Following an analogous approach, the authors derive here the Crooks equation in the context of molecular dynamics simulations of systems in the isothermal-isobaric (NPT) ensemble, whose dynamics is regulated by the Martyna-Tobias-Klein algorithm [J. Chem. Phys. 101, 4177 (1994)]. Their present derivation of the Crooks equation correlates to the demonstration of the Jarzynski identity for NPT systems recently proposed by Cuendet [J. Chem. Phys. 125, 144109 (2006)].     

Rotational excitations in para-H2+para-H2 collisions: full- and reduced-dimensional quantum wave packet studies comparing different potential energy surfaces

We study the process of rotational excitation in the collisions of para-H(2) with para-H(2) by propagating wave packets with the multiconfiguration time-dependent Hartree (MCTDH) algorithm. Transition probabilities are then calculated by the method of Tannor and Weeks based on time-correlation functions. Calculations were carried out up to a total angular momentum of J=70 to compute integral cross sections up to 1.2 eV in collision energy and thermal rate coefficients from 100 to 3000 K. The process is studied on the full-dimensional potential energy surface of Boothroyd-Martin-Keogh-Peterson (BMKP) as well as on the rigid rotor surface of Diep and Johnson. We test the validity of the rigid rotor approximation by also considering two rigid rotor restrictions of the BMKP potential energy surface (PES). Additionally, we investigate a variant of the BMKP PES suggested by Pogrebnya and Clary [Chem. Phys. Lett. 363, 523 (2002)] with reduced anisotropy. We compare our results with previous theoretical data for the cross sections and with experimental data for the rate coefficients at low temperatures.     

Electronic spectra of GdF reanalyzed by decomposing state functions according to f-shell angular momentum

The electronic structure of the GdF molecule was studied by means of four-component relativistic configuration interaction (CI) calculations [S. Yamamoto, H. Tatewaki, and T. Saue, J. Chem. Phys. 129, 244505 (2008)]. To analyze the electronic spectra more accurately, the CI wave function is decomposed according to the angular momentum (Ω(f)) generated from the (4f)(7) electrons. The weight of a specified Ω(f) is referred to as the "f-shell Omega component weight." This Ω(f) plays a crucial role in classifying the strong electronic transitions from the upper states (0.7 eV-3.0 eV) to the lower states (~0.55 eV). For these transitions, the upper and lower states have almost identical Ω(f) weights. This appears to be a necessary condition for a transition to be strong. The same condition is expected to hold for other lanthanide linear molecules. A point charge model is also studied, acting as a simplified model of GdF; it successfully reproduces the spectra of GdF, justifying studies based on ligand field theory.     

Complexes of type C6H7(+)·L (L = N2 and CO2) studied by explicitly correlated coupled cluster theory

Complexes of the benzenium ion (C(6)H(7)(+)) with N(2) or CO(2) have been studied by explicitly correlated coupled cluster theory at the CCSD(T)-F12x (x = a, b) level [T. B. Adler et al., J. Chem. Phys. 127, 221106 (2007)] and the double-hybrid density functional B2PLYP-D [T. Schwabe and S. Grimme, Phys. Chem. Chem. Phys. 9, 3397 (2007)]. Improved harmonic vibrational wavenumbers for C(6)H(7)(+) have been obtained by CCSD(T?)-F12a calculations with the VTZ-F12 basis set. Combining them with previous B2PLYP-D anharmonic contributions we arrive at anharmonic wavenumbers which are in excellent agreement with recent experimental data from p-H(2) matrix isolation IR spectroscopy [M. Bahou et al., J. Chem. Phys. 136, 154304 (2012)]. The energetically most favourable conformer of C(6)H(7)(+)·N(2) shows a π-bonded structure similar to C(6)H(7)(+)·Rg (Rg = Ne, Ar) [P. Botschwina and R. Oswald, J. Phys. Chem. A 115, 13664 (2011)] with D(e) ≈ 870 cm(-1). For C(6)H(7)(+)·CO(2), a slightly lower energy is calculated for a conformer with the CO(2) ligand lying in the ring-plane of the C(6)H(7)(+) moiety (D(e) ≈ 1508 cm(-1)). It may be discriminated from other conformers through a strong band predicted at 1218 cm(-1), red-shifted by 21 cm(-1) from the corresponding band of free C(6)H(7)(+).     

Efficient treatment of the Hartree interaction in the relativistic Kohn-Sham problem

We elaborate the two-component Douglas-Kroll reduction of the Dirac-Kohn-Sham problem of relativistic density-functional theory as introduced by Matveev and Rosch [J. Chem. Phys. 118, 3997 (2003)]. That method retains corrections to the Coulomb self-interaction (or Hartree) term of the energy functional that are due to the picture change. Using analytic expressions for the matrix elements, one is able to abandon the resolution of the identity approach for a crucial step of the relativistic transformation. Thus, a major source of uncertainties of the method is eliminated because basis sets no longer have to be extended by functions of higher angular momentum, previously required to ensure kinetic balance. This approach also relies on the electron charge-density fitting scheme via an auxiliary basis set. An efficient approximate implementation results if one restricts the relativistic transformation to the spherically symmetric atom-centered auxiliary functions. It provides accurate results while simplifying greatly the expressions for the matrix elements of the relativistically transformed operators and significantly reducing the computational effort. We demonstrate the performance of the method for the fine structure of one-electron levels of the Hg atom, the g-tensor shifts of NO2, and the properties of the diatomic molecules Bi2, Pb2, PbO, and TlH.     

A site-renormalized molecular fluid theory

The orientation-dependent pair distribution function for molecular fluids on site-site potentials is expanded in a topological analog of the diagrammatically proper site-site theory of liquids [D. Chandler et al., Mol. Phys. 46, 1335 (1982)]. The resulting functions are then used to diagrammatically renormalize the molecular fluid theory. A result is that the diagrammatically proper interaction site model theory is shown to be a linearized, minimal angular basis set approximation to this site-renormalized molecular theory. This framework is used to propose a new, exact, and proper closure to the diagrammatically proper interaction site model theory. The resulting equation system contains a bridge function expansion in the proper site-site theory. In addition, the construction of the theory is such that the molecular pair distribution function, in full dimensionality, is intrinsic to the theory. Furthermore, the theory is equivalent to the molecular Ornstein-Zernike treatment of site-site molecules in the basis set expansion of Blum and Torruella [J. Chem. Phys. 56, 303 (1971)]. A significant formal result of the theory is the demonstration that certain classes of diagrams which would otherwise be considered improper in the interaction site model formalism are included in the angular expansion of molecular interactions. Numerical results for several apolar homonuclear models and an apolar heteronuclear model are shown to quantitatively improve upon those of reference interaction site model and our recent proper variant with respect to simulation. Significant numerical results are that the various thermodynamic quantities obey the exact symmetries and sum rules within numerical error for the different sites in the heteronuclear case, even for the low order approximation used in this work, and the theory is independent of the so-called auxiliary site problem common to previous site-site theories.     

Polarizable and flexible model for ethanol

A polarizable, flexible model for ethanol is obtained based on an extensive series of B3LYP/6-311++G(d,p) calculations and molecular dynamics simulations. The ethanol model includes electric-field dependence in both the atomic charges and the intramolecular degrees of freedom. Field-dependent intramolecular potentials have been attempted only once previously, for OH and HH stretches in water [P. Cicu et al., J. Chem. Phys. 112, 8267 (2000)]. The torsional potential involving the hydrogen-bonding hydrogen in ethanol is found to be particularly field sensitive. The methodology for developing field-dependent potentials can be readily generalized to other molecules and is discussed in detail. Molecular dynamics simulations of bulk ethanol are performed and the results are assessed based on comparisons with the self-diffusion coefficient [N. Karger et al., J. Chem. Phys. 93, 3437 (1990)], dielectric constant [J. T. Kindt and C. A. Schmuttenmaer, J. Phys. Chem. 100, 10373 (1996)], enthalpy of vaporization [R. C. Wilhoit and B. J. Zwolinski, J. Phys. Chem. Ref. Data, Suppl. 2, 2 (1973)], and experimental interatomic distributions [C. J. Benmore and Y. L. Loh, J. Chem. Phys. 112, 5877 (2000)]. The simultaneous variation of the atomic charges and the intramolecular potentials requires modified equations of motion and a multiple time step algorithm has been implemented to solve these equations. The article concludes with a discussion of the bulk structure and properties with an emphasis on the hydrogen bonding network.     

Projection of two-dimensional diffusion in a curved midline and narrow varying width channel onto the longitudinal dimension

This study focuses on the derivation of a general effective diffusion coefficient to describe the two-dimensional (2D) diffusion in a narrow and smoothly asymmetric channel of varying width, in the simple diffusional motion of noninteracting pointlike particles under no external field. We present a generalization to the case of an asymmetric channel using the projection method introduced earlier by Kalinay and Percus [J. Chem. Phys. 122, 204701 (2005); and Phys. Rev. E 74, 041203 (2006)] to project the 2D diffusion equation into an effective one-dimensional generalized Fick-Jacobs equation. The expression for the diffusion coefficient given in Eq. (23) is our main result. This expression is a more general effective diffusion coefficient for narrow channels in 2D, which contains the well-known previous results as special cases, namely, those obtained by Bradley [Phys. Rev. E 80, 061142 (2009)], and more recently by Berezhkovskii and Szabo [J. Chem. Phys. 135, 074108 (2011)]. Finally, we study some specific 2D asymmetric channel configurations to test and show the broader applicability of Eq. (23).     

Comment on "MMM1D: A method for calculating electrostatic interactions in one-dimensional periodic geometries" [J. Chem. Phys. 123, 144103 (2005)

This comment gives an alternative derivation of the MMM1D method reported by Arnold and Holm [J. Chem. Phys.123, 144103 (2005)]. Moreover, several errors in expressions presented in the cited paper are identified.