Search

Subsequent to our spectroscopic detection of the HBX (X=F, Cl, Br) free radicals (S.-G. He, F. X. Sunahori, and D. J. Clouthier, J. Am. Chem. Soc. 127, 10814 (2005)), the electronic spectrum of the à 2A″Π-X 2A′ system of the fluoroborane (HBF) radical in the 600–745 nm region has been studied in detail using the pulsed discharge jet technique. The band system involves a linear-bent transition between the two Renner–Teller components of what would be a 2Π state at linearity. Using the results of our theoretical study of the ground and excited state vibrational energy levels and 11B–10B isotope shifts (see the companion paper), the vibrational quantum numbers of the bands in the laser-induced fluorescence (LIF) spectra have been assigned. Rotational and vibrational analyses of the LIF and wavelength resolved emission spectra have been carried out, from which the linear excited state and the bent ground state equilibrium configurations have been confirmed. The ground state molecular geometry of HBF has been determined as r0(BH)=1.214(2) Å, r0(BF)=1.303 4(5) Å, and θ=120.7(1)°. Based on high-level ab initio calculations and symmetry considerations, predissociation of the excited state into H(2S)+BF(1Σ+) on the ground state potential energy surface is identified as the cause of the breaking off of fluorescence in the LIF spectra. [ABSTRACT FROM AUTHOR]

Zeeman quantum beat spectroscopy has been used to measure the 300 K rate constants for the angular momentum depolarization of OH(A 2Σ+) in the presence of Ar. We show that the beat amplitude at short times, in the absence of collisions, is well described by previously developed line strength theory for (1+1) laser induced fluorescence. The subsequent pressure dependent decay of the beat amplitude is used to extract depolarization rate constants and estimates of collisional depolarization cross sections. Depolarization accompanies both inelastic collisions, giving rise to rotational energy transfer, and elastic collisions, which change mj but conserve j. Previous experimental studies, as well as classical theory, suggest that elastic scattering contributes around 20% to the observed total depolarization rate at low j. Simulation of the experimental beat amplitudes, using theoretical calculations presented in the preceding paper, reveals that depolarization of OH(A) by Ar has a rate constant comparable to, if not larger than, that for energy transfer. This is consistent with a significant tilting or realignment of j′ away from j on collision. The experimental data are used to provide a detailed test of quantum mechanical and quasiclassical trajectory scattering calculations performed on a recently developed ab initio potential energy surface of Kłos et al. [J. Chem. Phys. 129, 054301 (2008)]. The calculations and simulations account well for the observed cross sections at high N, but underestimate the experimental results by between 10% and 20% at low N, possibly due to remaining inaccuracies in the potential energy surface or perhaps to limitations in the dynamical approximations made, particularly the freezing of the OH(A) bond. [ABSTRACT FROM AUTHOR]

Classical and quantum mechanical expressions for the j-j′ vector correlation (also referred to as the rotational tilt) are presented for the situation in which the initial and final relative velocity directions are unresolved. The quantum mechanical expressions are compared with previous descriptions in the literature. It is shown that in the case of 2S+1Σ radicals in collision with closed shell species, a tensor opacity formalism can be employed in quasiclassical trajectory calculations to provide classical estimates of both open shell spin-rotation state and nuclear hyperfine state changing (or conserving) cross sections. Polarization parameters are also obtained from the same formalism. The method is applied to calculations on the OH(A 2Σ+)–Ar system using a recently developed potential energy surface. The results of both the closed and open shell quasiclassical trajectory calculations are found to compare favorably with those from close-coupled closed and open shell quantum mechanical scattering calculations. The accompanying paper provides an experimental test of these calculations and of the potential energy surface they employ. [ABSTRACT FROM AUTHOR]

In this paper, we have computed the rovibrational spectrum of the H3+ molecule using a new global potential energy surface, invariant under all permutations of the nuclei, that includes the long range electrostatic interactions analytically. The energy levels are obtained by a variational calculation using hyperspherical coordinates. From the comparison with available experimental results for low lying levels, we conclude that our accuracy is of the order of 0.1 cm-1 for states localized in the vicinity of equilateral triangular configurations of the nuclei, and changes to the order of 1 cm-1 when the system is distorted away from equilateral configurations. Full rovibrational spectra up to the H++H2 dissociation energy limit have been computed. The statistical properties of this spectrum (nearest neighbor distribution and spectral rigidity) show the quantum signature of classical chaos and are consistent with random matrix theory. On the other hand, the correlation function, even when convoluted with a smoothing function, exhibits oscillations which are not described by random matrix theory. We discuss a possible similarity between these oscillations and the ones observed experimentally. [ABSTRACT FROM AUTHOR]

The multistate vibronic dynamics in the X-D electronic states of all three difluorobenzene radical cations are investigated theoretically by an ab initio quantum dynamical approach. The vibronic coupling scheme and the ab initio values of the system parameters are adopted from Paper I [S. Faraji and H. Köppel, J. Chem. Phys. 129, 074310 (2008)]. Extensive calculations by wave-packet propagation have been performed with the aid of the multiconfiguration time-dependent Hartree method. Five coupled electronic potential energy surfaces and 10 (11 in the case of the orthoisomer) vibrational degrees of freedom have been included in these calculations. The nonadiabatic interactions lead to the restructuring of the photoelectron spectral envelopes. Ultrafast internal conversion processes within the electronic manifolds in question demonstrate the strength of the nonadiabatic coupling effects and complement the analogous findings for the electronic spectra. The internal conversion dynamics is characterized by a stepwise transfer of the electronic population to the lowest electronic state on a time scale of femtoseconds to picoseconds. A difference between the three isomers is found to be related to the weaker interaction between the sets of X-Ã and B-C-D states (with high-energy conical intersections) in the meta isomer, as compared to the other isomers. The implications of these findings for the qualitative understanding of the fluorescence dynamics of fluorinated benzene radical cations are discussed. [ABSTRACT FROM AUTHOR]

In this paper, we present the results of a theoretical investigation on the dynamics of the title reaction at collision energies below 1.2 kcal/mol using rigorous quantum reactive scattering calculations. Vibrationally resolved integral and differential cross sections, as well as product rotational distributions, have been calculated using two electronically adiabatic potential energy surfaces, developed by us on the basis of semiempirical modifications of the entrance channel. In particular, we focus our attention on the role of the exothermicity and of the exit channel region of the interaction on the experimental observables. From the comparison between the theoretical results, insight about the main mechanisms governing the reaction is extracted, especially regarding the bimodal structure of the HF(v=2) nascent rotational state distributions. A good overall agreement with molecular beam scattering experiments has been obtained. [ABSTRACT FROM AUTHOR]

A new semiempirical exchange-Coulomb model potential energy surface for the N2–He interaction was reported recently [A. K. Dham et al., J. Chem. Phys. 127, 054302 (2007)] and, using it, the temperature dependence of bulk gas properties of N2–He mixtures, such as the second virial coefficient and traditional transport phenomena, most of which depend primarily on the isotropic component of the interaction potential energy surface, was determined. Values of these properties, along with values calculated using two high-quality ab initio potential energy surfaces [C.-H. Hu and A. J. Thakkar, J. Chem. Phys. 104, 2541 (1996); K. Patel et al., ibid 119, 909 (2003)] were compared critically to available experimental data. The present paper reports on the ability of the same three potential energy surfaces to predict state-to-state and total differential cross sections, total integral cross sections, and the temperature dependence of bulk gas relaxation phenomena (including magnetic field effects on transport coefficients). While all three potential energy surfaces give total differential and higher speed integral scattering results that fall within the experimental uncertainties, integral scattering results and state-to-state differential cross section measurements consistently exceed the calculated values. All three surfaces give similar agreement with the relaxation properties of N2–He binary mixtures, with the semiempirical exchange-Coulomb model potential energy surface giving slightly better overall agreement with experiment than the two ab initio potential energy surfaces. [ABSTRACT FROM AUTHOR]

The frozen natural orbital (FNO) coupled-cluster method increases the speed of coupled-cluster (CC) calculations by an order of magnitude with no consequential error along a potential energy surface. This method allows the virtual space of a correlated calculation to be reduced by about half, significantly reducing the time spent performing the coupled-cluster (CC) calculation. This paper reports the derivation and implementation of analytical gradients for FNO-CC, including all orbital relaxation for both noncanonical and semicanonical perturbed orbitals. These derivatives introduce several new orbital relaxation contributions to the CC density matrices. FNO-CCSD(T) and FNO-ΛCCSD(T) are applied to a test set of equilibrium structures, verifying that these methods are capable of reproducing geometries and vibrational frequencies accurately, as well as energies. Several decomposition pathways of nitroethane are investigated using CCSD(T) and ΛCCSD(T) with 60% of the FNO virtual orbitals in a cc-pVTZ basis, and find differences on the order of 5 kcal/mol with reordering of the transition state energies when compared to B3LYP 6-311+G(3df,2p). [ABSTRACT FROM AUTHOR]

Several features and the performance of the recently published [P. Jankowski and M. Ziólkowski, Mol. Phys. 104, 2293 (2006)] three-dimensional intermolecular potential energy surface for the Ar–HF complex have been investigated. This full-dimensional surface has been obtained using the method of the local expansion of the exact interaction energy surface [P. Jankowski, J. Chem. Phys. 121, 1655 (2004)] in the Taylor series with respect to intramolecular coordinates. The interaction energies have been calculated with the coupled-cluster supermolecular method with single, double, and noniterative triple excitations. The convergence of the interaction energy with respect to the size of the basis set is discussed. The two-dimensional surfaces resulting from averaging of the full-dimensional surface over the intramolecular vibration of HF have been obtained and directly compared to the empirical H6(4,3,2) set of surfaces proposed by Hutson [J. Chem. Phys. 96, 6752 (1992)]. A very good agreement has been observed. The averaged potentials have been used to calculate the rovibrational energy levels of the Ar–HF and Ar–DF complexes and compared to the experimental data. The accuracy of rovibrational calculations achieved with the new surface is much better than with any of the ab initio surfaces available so far. Predictions of the rovibrational energy levels and spectroscopic constants have also been done for Ar–HF with HF in the v=4,5 vibrational states, and for Ar–DF and DF in the v=3,4 states. The full-dimensional surface studied in this paper is the first ab initio surface which is fully compatible with the empirical H6(4,3,2) surface proposed by Hutson. [ABSTRACT FROM AUTHOR]

Local coupled-cluster theory provides an algorithm for measuring electronic correlation quickly, using only the spatial locality of localized electronic orbitals. Previously, we showed [J. Subotnik et al., J. Chem. Phys. 125, 074116 (2006)] that one may construct a local coupled-cluster singles-doubles theory which (i) yields smooth potential energy surfaces and (ii) achieves near linear scaling. That theory selected which orbitals to correlate based only on the distances between the centers of different, localized orbitals, and the approximate potential energy surfaces were characterized as smooth using only visual identification. This paper now extends our previous algorithm in three important ways. First, locality is now based on both the distances between the centers of orbitals as well as the spatial extent of the orbitals. We find that, by accounting for the spatial extent of a delocalized orbital, one can account for electronic correlation in systems with some electronic delocalization using fast correlation methods designed around orbital locality. Second, we now enforce locality on not just the amplitudes (which measure the exact electron-electron correlation), but also on the two-electron integrals themselves (which measure the bare electron-electron interaction). Our conclusion is that we can bump integrals as well as amplitudes, thereby gaining a tremendous increase in speed and paradoxically increasing the accuracy of our LCCSD approach. Third and finally, we now make a rigorous definition of chemical smoothness as requiring that potential energy surfaces not support artificial maxima, minima, or inflection points. By looking at first and second derivatives from finite difference techniques, we demonstrate complete chemical smoothness of our potential energy surfaces (bumping both amplitudes and integrals). These results are significant both from a theoretical and from a computationally practical point of view. [ABSTRACT FROM AUTHOR]

How useful it is to think about the potential energy landscape of a complex many-body system depends in large measure on how direct the connection is to the system’s dynamics. In this paper we show that, within what we call the potential-energy-landscape ensemble, it is possible to make direct connections between the geometry of the landscape and the long-time dynamical behaviors of systems such as supercooled liquids. We show, in particular, that the onset of slow dynamics in such systems is governed directly by the lengths of their geodesics—the shortest paths through their landscapes within the special ensemble. The more convoluted and labyrinthine these geodesics are, the slower that dynamics is. Geodesics in the landscape ensemble have sufficiently well-defined characteristics that it is straightforward to search for them numerically, a point we illustrate by computing the geodesic lengths for an ordinary atomic liquid and a binary glass-forming atomic mixture. We find that the temperature dependence of the diffusion constants of these systems, including the precipitous drop as the glass-forming system approaches its empirical mode-coupling transition, is predicted quantitatively by the growth of the geodesic path lengths. [ABSTRACT FROM AUTHOR]

The nonrelativistic, valence-shell-only-correlated ab initio potential energy curve of the F2 molecule, which was reported in the preceding paper, is complemented by determining the energy contributions that arise from the electron correlations that involve the core electrons as well as the contributions that are due to spin-orbit coupling and scalar relativistic effects. The dissociation curve rises rather steeply toward the energy of the dissociated atoms because, at larger distances, the atomic quadrupole-quadrupole repulsion and spin-orbit coupling counteract the attractive contributions from incipient covalent binding and correlation forces including dispersion. [ABSTRACT FROM AUTHOR]

This paper discusses practical scheme for correcting the linear response coupled cluster with singles and doubles (CCSD) equations by shifting the poles corresponding to the equation-of-motion CCSD excitation energies by adding noniterative corrections due to triples. A simple criterion is derived for the excited states to be corrected in the spectral resolution of similarity transformed Hamiltonian on the CCSD level. Benchmark calculations were performed to compare the accuracies of static and dynamic polarizabilities obtained in this way with the CC3 and CCSDT counterparts. [ABSTRACT FROM AUTHOR]

The present paper is devoted to a full quantum mechanical study of the cis→trans isomerization of HONO. In contrast to our previous study [Richter et al., J. Chem. Phys. 120, 6072 (2004)], the dynamics is now performed in the presence of an external time-dependent field in order to be closer to experimental conditions. A six-dimensional dipole surface is computed. Using a previously developed potential energy surface [Richter et al., J. Chem. Phys. 120, 1306 (2004)], all eigenstates up to 4000 cm-1 are calculated. We simulate the dynamics during and after excitation by an electromagnetic pulse whose parameters are chosen to efficiently trigger the isomerization. Our investigations show that there is a selective isomerization pathway. [ABSTRACT FROM AUTHOR]

The photodissociation dynamics of a HX (X=Cl,F) molecule adsorbed on a hexagonal ice surface at T=0 K is studied using time-dependent quantum wave packets and quasiclassical trajectories. The relevant potential energy surfaces are calculated using high-level ab initio methods. We present here two dimensional calculations for the dynamics of the hydrogen photofragment for both HCl and HF molecules. The purpose of this paper is to compare the photodissociation dynamics of the two molecules which are adsorbed on the ice surface with different equilibrium geometries. The total photodissociation cross section and the angular distribution are calculated. The comparison with classical trajectory calculations provides evidence for typical quantum effects and reveals rainbow structures. [ABSTRACT FROM AUTHOR]

In this paper, the structure and spectroscopic parameters of the C5 cluster are determined using multiconfigurational quantum chemical methods as implemented in the MOLCAS software. A number of spectroscopic properties (band center positions, l-doubling parameters, and rotational constants) have been characterized. From the new results, the assignments of previous astrophysical observations [J. Goicoechea et al., Astrophys. J. 609, 225 (2004)] are discussed. A detailed exploration of the global potential energy surface confirms that C5 has a X 1Σg+ linear isomer of prominent stability and, at least, three minimum energy structures showing singlet electronic ground states. Two of them are cyclic and one has a nonplanar geometry. Vertical and adiabatic electronic transitions and vibrational spectroscopic parameters are determined for the most stable linear isomer using multiconfigurational second order perturbation theory (CASPT2) using an active space containing 12 valence orbitals with 12 active electrons and extended ANO-type basis sets. The infrared spectrum has been analyzed from an anharmonic force field derived form the local surface, determined from the energies of a grid of 1350 geometries. The force field includes four coupling terms. The CASPT2 band center position of the ν7(πu) anharmonic fundamental has been calculated to be at 102 cm-1, which validates the assignment to C5 of the pattern of bands centered at 102 cm-1 observed with the ISO telescope. [ABSTRACT FROM AUTHOR]

For the charge transfer and vibrational and electronic deexcitations between O2/O2++CO+/CO, O/O++CO2+/CO2, and C/C++O3+/O3, multistep reaction pathways are discussed in light of the theoretical data of this and previous paper together with close comparison with the experimental observations. Our calculations show that these pathways involve both the long range and molecular region ranges of the potential energy surfaces of the electronic states of the stable isomers of OOCO+ and mostly those of the weakly bound charge transfer complex OOCO+. The couplings between these electronic states such as vibronic, Renner-Teller, Jahn-Teller, and spin orbit are viewed to play crucial roles here. Moreover, the initial orientation of the reactants, in the entrance channels, strongly influences the reaction mechanisms undertaken. We propose for the first time a mechanism for the widely experimentally studied spin-forbidden exothermic O+(4Su)+CO2(X 1Σg+)→O2+(X 2Πg)+CO(X 1Σ+) reaction where the O turns around the OCO molecule. [ABSTRACT FROM AUTHOR]

In this paper, we present a calculation for the bound states of A1 symmetry on the spin-aligned Li3(1 4A′) potential energy surface. We apply a mixture of discrete variable representation and distributed approximating functional methods to discretize the Hamiltonian. We also introduce a new method that significantly reduces the computational effort needed to determine the lowest eigenvalues and eigenvectors (bound state energies and wave functions of the full Hamiltonian). In our study, we have found the lowest 150 energy bound states converged to less than 0.005% error, and most of the excited energy bound states converged to less than 2.0% error. Furthermore, we have estimated the total number of the A1 bound states of Li3 on the spin-aligned Li3(1 4A′) potential surface to be 601. [ABSTRACT FROM AUTHOR]

In this paper the authors use the optimal internal vibrational coordinates previously determined for the electronic ground state of the ozone molecule to study the vibrational spectrum of the molecule employing the second empirical potential energy surface calculated by Tyuterev et al. [Chem. Phys. Lett. 316, 271 (2000)]. First, the authors compute variationally all the bound vibrational energy levels of the molecule up to the dissociation limit and state the usefulness of the optimal coordinates in this respect, which allows us to converge all the bound levels using relatively small anharmonic basis sets. By analyzing the expansion coefficients of the wave functions, they show then that a large portion of the vibrational spectrum of O3 can be structured in nearly separable polyadic groups characterized by the polyad quantum number N=n1+n2+nθ corresponding to the optimal internal coordinates. Accordingly, they determine an internal effective vibrational Hamiltonian for O3 by fitting the effective Hamiltonian parameters to the experimental vibrational frequencies, using as input parameters in the fit those extracted from an analytical second-order Van Vleck perturbation theory calculation. It is finally shown that the internal effective Hamiltonian thus obtained accurately describes the vibrational spectrum of ozone in the low and medium energy regimes. [ABSTRACT FROM AUTHOR]

Ab initio calculations on the metal (groups 1 and 11) cyanide complexes show two stable configurations for the ground state geometry, a linear cyanide (MCN) and a triangular (MNC) form with an obtuse M–N–C angle. Lithium complex may exist in a linear isocyanide (MNC) form, but it cannot be differentiated from the triangular configuration because of the flatness of the potential energy surface connecting the two isomers. The metal atom and cyano radical are bonded through a strongly ionic configuration (M+CN-) in both geometrical forms. The MNC triangular form is a very floppy structure having one low frequency for the bending mode, whereas the MCN linear form is more rigid. The CN complexes of the alkali atoms have a triangular geometry as the lowest energy conformer, while the noble metal atoms prefer the linear cyanide one. The relative stability of the two isomers, dipole moments, and effective charges are reported in this paper. The essential aspects of the potential energy surfaces for the ground and the first excited states exhibiting a closely avoided crossing are also explained. [ABSTRACT FROM AUTHOR]

We report in this paper a quantum dynamics study for the reaction H+NH3→NH2+H2 on the potential energy surface of Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The quantum dynamics calculation employs the semirigid vibrating rotor target model [J. Z. H. Zhang, J. Chem. Phys. 111, 3929 (1999)] and time-dependent wave packet method to propagate the wave function. Initial state-specific reaction probabilities are obtained, and an energy correction scheme is employed to account for zero point energy changes for the neglected degrees of freedom in the dynamics treatment. Tunneling effect is observed in the energy dependency of reaction probability, similar to those found in H+CH4 reaction. The influence of rovibrational excitation on reaction probability and stereodynamical effect are investigated. Reaction rate constants from the initial ground state are calculated and are compared to those from the transition state theory and experimental measurement. [ABSTRACT FROM AUTHOR]

The thermodynamics and kinetics of a many-body system can be described in terms of a potential energy landscape in multidimensional configuration space. The partition function of such a landscape can be written in terms of a density of states, which can be computed using a variety of Monte Carlo techniques. In this paper, a new self-consistent Monte Carlo method for computing density of states is described that uses importance sampling and a multiplicative update factor to achieve rapid convergence. The technique is then applied to compute the equilibrium quench probability of the various inherent structures (minima) in the landscape. The quench probability depends on both the potential energy of the inherent structure and the volume of its corresponding basin in configuration space. Finally, the methodology is extended to the isothermal-isobaric ensemble in order to compute inherent structure quench probabilities in an enthalpy landscape. [ABSTRACT FROM AUTHOR]

Cis, cis-peroxynitrous acid is known to be an intermediate in atmospheric reactions between OH and NO2 as well as HOO and NO. The infrared absorption spectra of matrix-isolated cc-HOONO and cc-DOONO in argon have been observed in the range of 500–8000 cm-1. Besides the seven fundamental vibrational modes that have been assigned earlier for this molecule [Zhang et al., J. Chem. Phys. 124, 084305 (2006)], more than 50 of the overtone and combination bands have been observed for cc-HOONO and cc-DOONO. Ab initio CCSD(T)/atomic natural orbital anharmonic force field calculations were used to help guide the assignments. Based on this study of the vibrational overtone transitions of cis, cis-HOONO that go as high as 8000 cm-1 and the earlier paper on the vibrational fundamentals, we conclude that the CCSD(T)/ANO anharmonic frequencies seem to correct to ±35 cm-1. The success of the theoretically predicted anharmonic frequencies {υ} in assigning overtone spectra of HOONO up to 8000 cm-1 suggests that the CCSD(T)/ANO method is producing a reliable potential energy surface for this reactive molecule. [ABSTRACT FROM AUTHOR]

This paper reports state-of-the-art ab initio calculations at the second order of Mo\ller-Plesset perturbation theory of molecular hydrogen interaction with polycyclic aromatic hydrocarbons (PAHs) up to coronene (C24H12). The effects of both H2 orientation with respect to the PAH plane and of PAH corrugation have been carefully investigated. It was found that the energetic preference for the perpendicular H2 orientation over the parallel one decreases rapidly as the PAH grows. The detailed study of coronene-H2 potential energy surface has made it possible to estimate from first principles the graphite-H2 binding energy. This ab initio estimate is shown to be in fair agreement with experiment. As revealed by the example of benzene-H2 dimer, the electrostatic energy component plays an important role in PAH-H2 interactions, that stems largely from the charge penetration effect and, therefore, cannot be reproduced by the simple quadrupole-quadrupole interaction model. [ABSTRACT FROM AUTHOR]

In this paper we compare semiclassical initial value representation, conventional transition state theory with Wigner and Eckart tunneling correction, quantum reduced dimensionality, and quasiclassical thermal rate coefficients for N+N2 exchange reaction. [ABSTRACT FROM AUTHOR]

One of the outstanding issues in the quantum instanton (QI) theory (or any transition-state-type theory) for thermal rate constants of chemical reactions is the choice of an appropriate “dividing surface” (DS) that separates reactants and products. (In the general version of the QI theory, there are actually two dividing surfaces involved.) This paper shows one simple and general way for choosing DSs for use in QI theory, namely, using the family of (hyper) planes normal to the minimum energy path on the potential energy surface at various distances s along it. Here the reaction coordinate is not one of the dynamical coordinates of the system (which will in general be the Cartesian coordinates of the atoms), but rather simply a parameter which specifies the DS. It is also shown how this idea can be implemented for an N atom system in three-dimensional space in a way that preserves overall translational and rotational invariance. Numerical application to a simple system (the collinear H+H2 reaction) is presented to illustrate the procedure. [ABSTRACT FROM AUTHOR]

This paper presents a new ground state (1 2A′) electronic potential energy surface for the F(2P)+HCl→HF+Cl(2P) reaction. The ab initio calculations are done at the multireference configuration interaction+Davidson correction (MRCI+Q) level of theory by complete basis set extrapolation of the aug-cc-pVnZ (n=2,3,4) energies. Due to low-lying charge transfer states in the transition state region, the molecular orbitals are obtained by six-state dynamically weighted multichannel self-consistent field methods. Additional perturbative refinement of the energies is achieved by implementing simple one-parameter correlation energy scaling to reproduce the experimental exothermicity (ΔE=-33.06 kcal/mol) for the reaction. Ab initio points are fitted to an analytical function based on sum of two- and three-body contributions, yielding a rms deviation of <0.3 kcal/mol for all geometries below 10 kcal/mol above the barrier. Of particular relevance to nonadiabatic dynamics, the calculations show significant multireference character in the transition state region, which is located 3.8 kcal/mol with respect to F+HCl reactants and features a strongly bent F–H–Cl transition state geometry (θ≈123.5°). Finally, the surface also exhibits two conical intersection seams that are energetically accessible at low collision energies. These seams arise naturally from allowed crossings in the C∞v linear configuration that become avoided in Cs bent configurations of both the reactant and product, and should be a hallmark of all X–H–Y atom transfer reaction dynamics between (2P) halogen atoms. [ABSTRACT FROM AUTHOR]

When a minimum on the potential energy surface is surrounded by multiple saddle points with similar energy barriers, the transition pathways with greater prefactors are more important than those that have similar energy barriers but smaller prefactors. In this paper, we present a theoretical formulation for the prefactors, computing the probabilities for transition paths from a minimum to its surrounding saddle points. We apply this formulation to a system of 2 degrees of freedom and a system of 14 degrees of freedom. The first is Brownian motion in a two-dimensional potential whose global anharmonicities play a dominant role in determining the transition rates. The second is a Lennard-Jones (LJ) cluster of seven particles in two dimensions. Low lying transition states of the LJ cluster, which can be reached directly from a minimum without passing through another minimum, are identified without any presumption of their characteristics nor of the product states they lead to. The probabilities are computed for paths going from an equilibrium ensemble of states near a given minimum to the surrounding transition states. These probabilities are directly related to the prefactors in the rate formula. This determination of the rate prefactors includes all anharmonicities, near or far from transition states, which are pertinent in the very sophisticated energy landscape of LJ clusters and in many other complex systems. [ABSTRACT FROM AUTHOR]

We present several numerical applications based upon the effective-mode formulation for the short-time dynamics through conical intersections in macrosystems, as detailed in the preceding paper and first proposed by Cederbaum et al. [Phys. Rev. Lett. 94, 113003 (2005)]. The macrosystem, containing a vast number of nuclear degrees of freedom (modes), is decomposed into a system part and an environment part. Only three effective environmental modes are needed—together with the system’s modes—to accurately calculate the low resolution spectra and the short-time dynamics of the entire macrosystem. For the systems discussed here, results are compared to those of a full quantum wave-packet propagation. Some rules are extracted to provide general tendencies; these rules allow one to understand and predict the dynamical properties in more general situations where the exact quantum dynamics of the macrosystem is out of reach. [ABSTRACT FROM AUTHOR]

The short-time dynamics through a conical intersection of a macrosystem comprising a large number of nuclear degrees of freedom (modes) is investigated. The macrosystem is decomposed into a “system” part carrying a limited number of modes, and an “environment” part. An orthogonal transformation in the environment’s space is introduced, as a result of which a subset of three effective modes can be identified which couple directly to the electronic subsystem. Together with the system’s modes, these govern the short-time dynamics of the overall macrosystem. The remaining environmental modes couple, in turn, to the effective modes and become relevant at longer times. In this paper, we present the derivation of the effective Hamiltonian, first introduced by Cederbaum et al. [Phys. Rev. Lett. 94, 113003 (2005)], and analyze its properties in some detail. Several special cases and topological aspects are discussed. [ABSTRACT FROM AUTHOR]

In this paper we present a theoretical study using time-dependent nonadiabatic reactant-product decoupling method for the state-to-state reactive scattering calculation of F(2P1/2)+H2 (ν=j=0) reaction on the Alexander-Stark-Werner potential energy surface. In this nonadiabatic state-to-state calculation, the full wave function is partitioned into reactant component and a sum of all product components. The reactant and product components of the wave function are solved independently. For the excited state reaction, the state-to-state reaction probabilities for J=0.5 are calculated. Comparing the state-to-state reaction probabilities, it is found that the vibrational population of the HF product is dominated by vibrational levels ν=2 and 3. The rotation specific reaction probabilities of HF product in j=1 and 2 are larger than those in other rotational levels. As the rotation quantum number j increases, the positions of the peak in the rotational reaction probability of HF product in ν=3 shift to higher collision energy. [ABSTRACT FROM AUTHOR]

The neural network (NN) procedure to interpolate ab initio data for the purpose of molecular dynamics (MD) simulations has been tested on the SiO2 system. Unlike other similar NN studies, here, we studied the dissociation of SiO2 without the initial use of any empirical potential. During the dissociation of SiO2 into Si+O or Si+O2, the spin multiplicity of the system changes from singlet to triplet in the first reaction and from singlet to pentet in the second. This paper employs four potential surfaces. The first is a NN fit [NN(STP)] to a database comprising the lowest of the singlet, triplet, and pentet energies obtained from density functional calculations in 6673 nuclear configurations. The other three potential surfaces are obtained from NN fits to the singlet, triplet, and pentet-state energies. The dissociation dynamics on the singlet-state and NN(STP) surfaces are reported. The results obtained using the singlet surface correspond to those expected if the reaction were to occur adiabatically. The dynamics on the NN(STP) surface represent those expected if the reaction follows a minimum-energy pathway. This study on a small system demonstrates the application of NNs for MD studies using ab initio data when the spin multiplicity of the system changes during the dissociation process. [ABSTRACT FROM AUTHOR]

Dissociations after the à 1B1→B 1A1 photoexcitation of SiH2, SiHD, and SiD2 were studied to investigate excited-state dynamics and effects of the initial vibrational state. The cross section (σ) for the photodissociation relative to SiH2(B)→Si(1D)+H2 and the rovibrational population of the H2 fragment were computed using the wave packet propagation technique based on the three-dimensional potential energy surfaces (PESs) of the à and B electronic states and the transition dipole surfaces, which were reported in our previous paper [J. Chem. Phys. 122, 144307 (2005)]. The photodissociation spectrum consists of a broadband and a number of sharp peaks. For SiH2 and SiD2, the sharp peaks correspond to the resonance structure of the vibrational levels of the B state and the broadbands are nearly independent of the photon energy. The broadband for SiHD increases steeply with the photon energy above 30 000 cm-1. The flux leaving the computational grid for SiH2 and SiD2 consists of at least two components, whereas that for SiHD consists of only a faster component. These large isotope effects were discussed based on the valley to the dissociation channel on PES and the difference in the position of the initial wave packet for three isotopomers. [ABSTRACT FROM AUTHOR]

A method for constructing diabatic potential-energy matrices from ab initio quantum chemistry data is described and tested for use in exact quantum reactive scattering. The method is a refinement of that presented in a previous paper, in that it accounts for the presence of the nonremovable derivative coupling. The accuracy of quantum dynamics on this type of diabatic potential is tested by comparison with an analytic model and for an ab initio description of the two lowest-energy states of H3. [ABSTRACT FROM AUTHOR]

We have calculated new potential-energy surfaces for the lowest two spin-aligned 4A′ states of the Li3 trimer. This calculation shows a seam of conical intersections between these states resulting from the extra symmetry of the system when the atoms are in a collinear arrangement. This seam is especially important because of its proximity to the three-body dissociation limit of the system; ultracold scattering calculations and the bound-state energies of the system will be affected by the presence of this conical intersection. In this paper we discuss the calculation of the potential-energy surface and the location of the conical intersection seam. [ABSTRACT FROM AUTHOR]

All the pure rotational transitions reported in the previous studies [J. Chem. Phys. 113, 10121 (2000); J. Mol. Spectrosc. 222, 22 (2003)] and newly observed rotation-vibration transitions, P=1/2←3/2, for Ar–SH and Ar–SD [J. Chem. Phys. (2005), the preceding paper] have been simultaneously analyzed to determine a new intermolecular potential-energy surface of Ar–SH in the ground state. A Schrödinger equation considering the three-dimensional freedom of motion for an atom-diatom complex in the Jacobi coordinate, R, θ, and r, was numerically solved to obtain energies of the rovibrational levels using the discrete variable representation method. A three-dimensional potential-energy surface is determined by a least-squares fitting with initial values of the parameters for the potential obtained by ab initio calculations at the RCCSD(T)/aug-cc-pVTZ level of theory. The potential well reproduces all the observed data in the microwave and millimeter wave regions with parity doublings and hyperfine splittings. Several low-lying rovibrational energies are calculated using the new potential-energy surface. The dependence of the interaction energy between Ar and SH(2Πi) on the bond length of the SH monomer is discussed. [ABSTRACT FROM AUTHOR]

Fundamental issues associated with the application of the nuclear-electronic orbital (NEO) approach to hydrogen transfer systems are addressed. In the NEO approach, specified nuclei are treated quantum mechanically on the same level as the electrons, and mixed nuclear-electronic wavefunctions are calculated with molecular orbital methods. The positions of the nuclear basis function centers are optimized variationally. In the application of the NEO approach to hydrogen transfer systems, the hydrogen nuclei and all electrons are treated quantum mechanically. Within the NEO framework, the transferring hydrogen atom can be represented by two basis function centers to allow delocalization of the proton vibrational wavefunction. In this paper, the NEO approach is applied to the [He-H-He]+ and [He-H-He]++ model systems. Analyses of technical issues pertaining to flexibility of the basis set to describe both single and double well proton potential energy surfaces, linear dependency of the hydrogen basis functions, multiple minima in the basis function center optimization, convergence of the number of hydrogen basis function centers, and basis set superposition error are presented. The accuracy of the NEO approach is tested by comparison to grid calculations for these model systems. [ABSTRACT FROM AUTHOR]

In this paper we report a combined experimental and theoretical study on the dynamics of the insertion reaction C(1D)+D2 at 15.5 kJ mol-1 collision energy. Product angular and velocity distributions have been obtained in crossed beam experiments and quasiclassical trajectory (QCT) and rigorous statistical calculations have been performed on the recent and accurate ab initio potential energy surface of Bussery-Honvault, Honvault, and Launay at the energy of the experiment. The molecular-beam results have been simulated using the theoretical calculations. Good agreement between experiment and both QCT and statistical predictions is found. [ABSTRACT FROM AUTHOR]

The paper presents a theoretical study of the low-energy dynamics of radiative association processes in the He++H2 collision system. Formation of the triatomic HeH2+ ion in its bound rotation-vibration states on the potential-energy surfaces of the ground and of the first excited electronic states is investigated. Close-coupling calculations are performed to determine detailed state-to-state characteristics (bound←free transition rates, radiative and dissociative widths of resonances) as well as temperature-average characteristics (rate constants, photon emission spectra) of the two-state (X←A) reaction He+(2S)+H2(X1Σg+)→HeH2+(X2A′)+hν and of the single-state (A←A) reaction He+(2S)+H2(X1Σg+)→HeH2+(A2A′)+hν. The potential-energy surfaces of the X- and A-electronic states of HeH2+ and the dipole moment surfaces determined ab initio in an earlier work [Kraemer, Sˇpirko, and Bludsky, Chem. Phys. 276, 225 (2002)] are used in the calculations. The rate constants k(T) as functions of temperature are calculated for the temperature interval 1≤T≤100 K. The maximum k(T) values are predicted as 3.3×10-15 s-1 cm3 for the X←A reaction and 2.3×10-20 s-1 cm3 for the A←A reaction at temperatures around 2 K. Rotationally predissociating states of the He+-H2 complex, correlating with the υ=0, j=2 state of free H2, are found to play a crucial role in the dynamics of the association reactions at low temperatures; their contribution to the k(T) function of the X←A reaction at T<30 K is estimated as larger than 80%. The calculated partial rate constants and emission spectra show that in the X←A reaction the HeH2+(X) ion is formed in its highly excited vibrational states. This is in contrast with the vibrational state population of the ion when formed via the (X←X) reaction He(1S)+H2+(X2Σg+)→HeH2+(X2A′)+hν. [ABSTRACT FROM AUTHOR]

A method for incorporating strong electric field polarization effects into optimal control calculations is presented. A Born–Oppenheimer-type separation, referred to as the electric-nuclear Born–Oppenheimer (ENBO) approximation, is introduced in which variations of both the nuclear geometry and the external electric field are assumed to be slow compared with the speed at which the electronic degrees of freedom respond to these changes. This assumption permits the generation of a potential energy surface that depends not only on the relative geometry of the nuclei but also on the electric field strength and on the orientation of the molecule with respect to the electric field. The range of validity of the ENBO approximation is discussed in the paper. A two-stage toolkit implementation is presented to incorporate the polarization effects and reduce the cost of the optimal control dynamics calculations. As an illustration of the method, it is applied to optimal control of vibrational excitation in a hydrogen molecule aligned along the field direction. Ab initio configuration interaction calculations with a large orbital basis set are used to compute the H–H interaction potential in the presence of the electric field. The significant computational cost reduction afforded by the toolkit implementation is demonstrated. [ABSTRACT FROM AUTHOR]

Infrared spectra of weakly bound hydrogen–krypton complexes have been studied at high spectral resolution (0.04 cm-1) using a long-path (154 m) low temperature (100 K) absorption cell and a Fourier transform spectrometer. In addition to spectra from the regions of the H2 and D2 fundamental vibrational bands in the midinfrared, the results also include the region of the pure rotational S0(0) transition of H2 in the far infrared. A total of 219 measured line positions from these spectra have been fully assigned to specific quantum transitions and form the basis for determining a greatly improved semiempirical three-dimensional intermolecular potential energy surface for hydrogen–krypton in an accompanying paper. [ABSTRACT FROM AUTHOR]

The recently developed completely renormalized (CR) coupled-cluster (CC) methods with singles, doubles, and noniterative triples or triples and quadruples [CR-CCSD(T) or CR-CCSD(TQ), respectively], which are based on the method of moments of CC equations (MMCC) [K. Kowalski and P. Piecuch, J. Chem. Phys. 113, 18 (2000)], eliminate the failures of the standard CCSD(T) and CCSD(TQ) methods at larger internuclear separations, but they are not rigorously size extensive. Although the departure from strict size extensivity of the CR-CCSD(T) and CR-CCSD(TQ) methods is small, it is important to examine the possibility of formulating the improved CR-CC methods, which are as effective in breaking chemical bonds as the existing CR-CCSD(T) and CR-CCSD(TQ) approaches, which are as easy to use as the CR-CCSD(T) and CR-CCSD(TQ) methods, and which can be made rigorously size extensive. This may be particularly useful for the applications of CR-CC methods and other MMCC approaches in calculations of potential energy surfaces of large many-electron systems and van der Waals molecules, where the additive separability of energies in the noninteracting limit is very important. In this paper, we propose different types of CR-CC approximations, termed the locally renormalized (LR) CCSD(T) and CCSD(TQ) methods, which become rigorously size extensive if the orbitals are localized on nointeracting fragments. The LR-CCSD(T) and LR-CCSD(TQ) methods rely on the form of the energy expression in terms of the generalized moments of CC equations, derived in this work, termed the numerator-denominator-connected MMCC expansion. The size extensivity and excellent performance of the LR-CCSD(T) and LR-CCSD(TQ) methods are illustrated numerically by showing the results for the dimers of stretched HF and LiH molecules and bond breaking in HF and H2O. [ABSTRACT FROM AUTHOR]

A potential energy surface is retrieved for the Ar-CO complex by carrying out a global analysis of its high-resolution spectroscopic data. The data set consists of already published microwave and infrared data and of new microwave transitions which are presented in the paper. The theoretical approach used to reproduce the spectrum is based on a model Hamiltonian which accounts simultaneously for the two large amplitude van der Waals modes and for the overall rotation of the complex. Only the vCO=0 state is considered. The root-mean-square deviation of the analysis is 18 MHz for the microwave data and 1.4×10-3 cm-1 for the infrared energy difference data. Fifteen parameters corresponding to the potential energy function are determined in addition to two kinetic energy parameters and two distortion-type parameters. The potential energy surface derived is in good agreement with the one obtained by Shin, Shin, and Tao [J. Chem. Phys. 104, 183 (1996)]. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Two recent papers presented calculations of the highly excited vibrational states of ozone [J. Chem. Phys. 119, 6512 (2003); 119, 6554 (2003)]. The nature and energies of these states may hold the key to the anomalous isotopic distribution of ozone in the atmosphere. Even though the same potential energy surface of Babikov et al. [J. Chem. Phys. 119, 2577 (2003)] was used in both calculations, the number of bound van der Waals states reported below dissociation differed significantly. In order to resolve the issue we present here the results of an independent computation of all the bound vibrational states of [sup 16]O[sup 16]O[sup 16]O and [sup 16]O[sup 16]O[sup 18]O up to dissociation. Our methods differ from both earlier calculations since we use hyperspherical coordinates and a direct product discrete variable representation of the Hamiltonian. The results of present work support the existence of several van der Waals states for J=0 on this potential energy surface. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The traditional multireference (MR) coupled-cluster (CC) methods based on the effective Hamiltonian are often beset by the problem of intruder states, and are not suitable for studying potential energy surface (PES) involving real or avoided curve crossing. State-specific MR-based approaches obviate this limitation. The state-specific MRCC (SS-MRCC) method developed some years ago [Mahapatra et al., J. Chem. Phys. 110, 6171 (1999)] can handle quasidegeneracy of varying degrees over a wide range of PES, including regions of real or avoided curve-crossing. Motivated by its success, we have suggested and explored in this paper a suite of physically motivated coupled electron-pair approximations (SS-MRCEPA) like methods, which are designed to capture the essential strength of the parent SS-MRCC method without significant sacrificing its accuracy. These SS-MRCEPA theories, like their CC counterparts, are based on complete active space, treat all the reference functions on the same footing and provide a description of potentially uniform precision of PES of states with varying MR character. The combining coefficients of the reference functions are self-consistently determined along with the cluster amplitudes themselves. The newly developed SS-MRCEPA methods are size-extensive, and are also size-consistent with localized orbitals. Among the various versions, there are two which are invariant with respect to the restricted rotations among doubly occupied and active orbitals separately. Similarity of performance of this latter and the noninvariant versions at the crossing points of the degenerate orbitals imply that the all the methods presented are rather robust with respect to the rotations among degenerate orbitals. Illustrative numerical applications are presented for PES of the ground state of a number of difficult test cases such as the model H4, H[sub 8] problems, the insertion of Be into H[sub 2], and Li[sub 2], where intruders exist and for a state of a molecule such as CH[sub 2], with pronounced MR character. Results obtained with SS-MRCEPA methods are found to be comparable in accuracy to the parent SS-MRCC and FCI/large scale CI results throughout the PES, which indicates the efficacy of our SS-MRCEPA methods over a wide range of geometries, despite their neglect of a host of complicated nonlinear terms, even when the traditional MR-based methods based on effective Hamiltonians fail due to intruders. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

In this paper, we employ the time-dependent quantum wave packet method to study the reaction of F([sup 2]P[sub 3/2], [sup 2]P[sub 1/2]) with HD on the Alexander–Stark–Werner potential energy surface. The reaction probabilities and total integral cross sections of the spin–orbit ground and excited states for the two possible products of the system are calculated. Because the reaction channel of the excited spin–orbit state is closed at the resonance energy, the resonance feature does not appear in the reaction probabilities and cross section for the F([sup 2]P[sub 1/2])+HD(v=j=0)→HF+D reaction, in contrast with that found for the ground spin–orbit state. We also compare the average cross sections of the two possible products with the experimental measurement. The resonance peak in the present average cross section for the HF+D product is slightly larger than the experimental result, but much smaller than that of the single-state calculations on the potential energy surface of Stark and Werner. It seems that the spin–orbit coupling would play a relatively important role in this reaction. Moreover, the isotope effects of the ground and excited spin–orbit states and the reactivity of the two product channels from the excited spin–orbit state are presented. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

In this paper effects of higher order Jahn-Teller coupling terms on the nonadiabatic dynamics are studied. Of particular interest is the case when the potential energy surfaces of the degenerate state show pronounced anharmonicity. In order to demonstrate the effects a two-dimensional Exe Jahn-Teller model system is treated which is based on the e[sup ′] stretching vibration of the photoactive [sup 2]E[sup ′] state of NO[sub 3] as a realistic example. The sixth order Exe Jahn-Teller Hamiltonian is derived in the diabatic representation which is valid for any system with a C[sub 3] rotation axis. This diabatization scheme is compared to lower-order Jahn-Teller Hamiltonians and to symmetry adapted as well as ad hoc approximations. Lower-order potentials result in pronounced quantitative and qualitative differences in the dynamics, including differences in the evolution of mean values, the autocorrelation functions (and thus the corresponding spectra), and the electronic population evolution. In the particular example treated, the results of fourth and fifth order potentials are very similar to the sixth order reference system. In contrast, the approximate sixth order Hamiltonians, though the corresponding adiabatic surfaces seem to be nearly identical, results in pronounced differences. The possible consequences for the dynamics of realistic systems with higher dimensionality are briefly discussed. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

This paper compares harmonic and anharmonic zero-point energies and thermodynamic functions for a number of molecules of small and medium size. Anharmonic corrections cannot be neglected for quantitative studies, but can be obtained quite effectively by a perturbative treatment including cubic force constants to the second order and semidiagonal quartic constants to the first order. Simple finite difference equations provide all the necessary terms by at most 6N-11 Hessian evaluations, where N is the number of atoms in the system. Accurate values are obtained by this method using the Becke three parameter Lee–Yang–Parr functional, medium size basis sets, and, when needed, proper treatment of internal rotations. The whole model has been completely automated in the Gaussian package. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

This paper examines the resonances in H+HLi scattering. The signature of these resonances is obtained from the oscillations in its reaction probability versus energy curves. They are identified here from a set of pseudospectra calculated for different initial locations of a stationary Gaussian wave packet on the ab initio potential energy surface (PES) reported by Dunne, Murrel, and Jemmer [Chem. Phys. Lett. 336, 1 (2001)]. The nuclear motion on this PES is monitored with the aid of a time-dependent wave packet method and the pseudospectrum are calculated by Fourier transforming the time autocorrelation function of the initial wave packet. The resonances are further examined and assigned by computing their eigenfunctions through spectral quantization algorithm. Both the linewidth as well as decay lifetimes of the resonances are reported. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The bound rovibronic levels of the He–HF[sup +] complex were calculated for total angular momentum J=1/2, 3/2, 5/2, 7/2, and 9/2 with the use of ab initio diabatic intermolecular potentials presented in Paper I and the inclusion of spin–orbit coupling. The character of the rovibronic states was interpreted by a series of calculations with the intermolecular distance R fixed at values ranging from 1.5 to 8.5 Å and by analysis of the wave functions. In this analysis we used approximate angular momentum quantum numbers defined with respect to a dimer body-fixed (BF) frame with its z axis parallel to the intermolecular vector R and with respect to a molecule-fixed (MF) frame with its z axis parallel to the HF[sup +] bond. The linear equilibrium geometry makes the He–HF[sup +] complex a Renner–Teller system. We found both sets of quantum numbers, BF and MF, useful to understand the characteristics of the Renner–Teller effect in this system. In addition to the properties of a “normal” semirigid molecule Renner–Teller system it shows typical features caused by large-amplitude internal (bending) motion. We also present spectroscopic data: stretch and bend frequencies, spin–orbit splittings, parity splittings, and rotational constants. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]