Vibrational spectrum of BiH3: Six-dimensional variational calculations on high-level ab initio potential energy surfaces

We report a theoretical study of the ground electronic state of BiH₃. The potential energy surface (PES) is obtained from coupled cluster CCSD(T) calculations with a large basis set (289 contracted Gaussian functions). The previously available quartic force field (P4) is extended by adding the dominant quintic and sextic stretching terms to yield improved potential functions in symmetry coordinates (P6) and Morse-type coordinates (M4). Second-order rovibrational perturbation calculations on the P4-PES and full variational calculations on the P6-PES and M4-PES yield almost identical vibrational term values which is rationalized by considering the local mode behavior of BiH₃ and the Morse-type character of the M4-PES. The remaining deviations between the computed and observed vibrational term values must thus be caused by imperfections in the CCSD(T) surface. A refinement of this ab initio surface by a restrained fit to experimental data allows an essentially perfect reproduction of the observed vibrational term values. Variational calculations on this refined surface provide predictions for several overtone and combination bands that have not yet been observed. Dedicated to Hermann Stoll on the occasion of his 60th birthday     

Ab initio based potential energy surfaces and Franck-Condon analysis of ionization thresholds of cyclic-C3H and linear-C3H

We report a Franck-Condon analysis in reduced dimensionality of the ionization thresholds of linear(l)-C3H and cyclic(c)-C3H using MP2-based potential energy surfaces and CCSD(T)/aug-cc-pVTZ calculations of electronic energies at selected geometries. The potential energy surfaces are fits to tens of thousands of MP2/aug-cc-pVTZ energies for the neutral and cation systems. These fits properly describe the invariance of the potential with respect to all permutations of the three C atoms. The realism of the potential surfaces is assessed by comparing stationary-point structures, energies, and normal-mode frequencies with previous high-level ab initio calculations. Several key vibrational modes in this ionization process are located at saddle points and so a numerical approach to obtain the Franck-Condon factors for those modes is done. On the basis of this analysis combined with a simple harmonic treatment of the energies of the remaining modes and key electronic energy differences obtained with CCSD(T)/aug-cc-pVTZ calculations, we find the threshold ionization energy of l-C3H to be 9.06 eV and for c-C3H we estimate the threshold to be in the range 9.70-9.76 eV. We estimate these values are accurate to within +/-0.05 eV.     

Calibration-quality adiabatic potential energy surfaces for H3(+) and its isotopologues

Calibration-quality ab initio adiabatic potential energy surfaces (PES) have been determined for all isotopologues of the molecular ion H(3)(+). The underlying Born-Oppenheimer electronic structure computations used optimized explicitly correlated shifted Gaussian functions. The surfaces include diagonal Born-Oppenheimer corrections computed from the accurate electronic wave functions. A fit to the 41,655 ab initio points is presented which gives a standard deviation better than 0.1 cm(-1) when restricted to the points up to 6000 cm(-1) above the first dissociation asymptote. Nuclear motion calculations utilizing this PES, called GLH3P, and an exact kinetic energy operator given in orthogonal internal coordinates are presented. The ro-vibrational transition frequencies for H(3)(+), H(2)D(+), and HD(2)(+) are compared with high resolution measurements. The most sophisticated and complete procedure employed to compute ro-vibrational energy levels, which makes explicit allowance for the inclusion of non-adiabatic effects, reproduces all the known ro-vibrational levels of the H(3)(+) isotopologues considered to better than 0.2 cm(-1). This represents a significant (order-of-magnitude) improvement compared to previous studies of transitions in the visible. Careful treatment of linear geometries is important for high frequency transitions and leads to new assignments for some of the previously observed lines. Prospects for further investigations of non-adiabatic effects in the H(3)(+) isotopologues are discussed. In short, the paper presents (a) an extremely accurate global potential energy surface of H(3)(+) resulting from high accuracy ab initio computations and global fit, (b) very accurate nuclear motion calculations of all available experimental line data up to 16,000 cm(-1), and (c) results suggest that we can predict accurately the lines of H(3)(+) towards dissociation and thus facilitate their experimental observation.     

Quantum and classical studies of vibrational motion of CH5+ on a global potential energy surface obtained from a novel ab initio direct dynamics approach

We report a full dimensional, ab initio based potential energy surface for CH(5) (+). The ab initio electronic energies and gradients are obtained in direct-dynamics calculations using second-order M?ller-Plesset perturbation theory with the correlation consistent polarized valence triple zeta basis. The potential energy and the dipole moment surfaces are fit using novel procedures that ensure the full permutational symmetry of the system. The fitted potential energy surface is tested by comparing it against additional electronic energy calculations and by comparing normal mode frequencies at the three lowest-lying stationary points obtained from the fit against ab initio ones. Well-converged diffusion Monte Carlo zero-point energies, rotational constants, and projections along the CH and HH bond lengths and the tunneling coordinates are presented and compared with the corresponding harmonic oscillator and standard classical molecular dynamics ones. The delocalization of the wave function is analyzed through comparison of the CH(5) (+) distributions with those obtained when all of the hydrogen atoms are replaced by (2)H and (3)H. The classical dipole correlation function is examined as a function of the total energy. This provides a further probe of the delocalization of CH(5) (+).     

Collision dynamics of O(3P) + DMMP using a specific reaction parameters potential form

Starting from previous benchmark CBS-QB3 electronic structure calculations (Conforti, P. F.; Braunstein, M.; Dodd, J. A. J. Phys. Chem. A 2009, 113, 13752), we develop two global potential energy surfaces for O((3)P) + DMMP collisions, using the specific reaction parameters approach. Each surface is simultaneously fit along the three major reaction pathways: hydrogen abstraction, hydrogen elimination, and methyl elimination. We then use these surfaces in classical dynamics simulations and compute reactive cross sections from 4 to 10 km s(-1) collision velocity. We examine the energy disposal and angular distributions of the reactive and nonreactive products. We find that for reactive collisions, an unusually large amount of the initial collision energy is transformed into internal energy. We analyze the nonreactive and reactive product internal energy distributions, many of which fit Boltzmann temperatures up to ~2000 K.     

An Anderson impurity model for efficient sampling of adiabatic potential energy surfaces of transition metal complexes

We present a model intended for rapid sampling of ground and excited state potential energy surfaces for first-row transition metal active sites. The method is computationally inexpensive and is suited for dynamics simulations where (1) adiabatic states are required "on-the-fly" and (2) the primary source of the electronic coupling between the diabatic states is the perturbative spin-orbit interaction among the 3d electrons. The model Hamiltonian we develop is a variant of the Anderson impurity model and achieves efficiency through a physically motivated basis set reduction based on the large value of the d-d Coulomb interaction U(d) and a Lanczos matrix diagonalization routine to solve for eigenvalues. The model parameters are constrained by fits to the partial density of states obtained from ab initio density functional theory calculations. For a particular application of our model we focus on electron transfer occurring between cobalt ions solvated by ammonium, incorporating configuration interaction between multiplet states for both metal ions. We demonstrate the capability of the method to efficiently calculate adiabatic potential energy surfaces and the electronic coupling factor we have calculated compares well to previous calculations and experiment. (     

Density functional theory study of the geometry, energetics, and reconstruction process of Si111 surfaces

We report the structures and energies from first principles density functional calculations of 12 different reconstructed (111) surfaces of silicon, including the 3x3 to 9x9 dimer-adatom-stacking fault (DAS) structures. These calculations used the Perdew-Burke-Ernzerhof generalized gradient approximation of density functional theory and Gaussian basis functions. We considered fully periodic slabs of various thicknesses. We find that the most stable surface is the DAS 7x7 structure, with a surface energy of 1.044 eV/1x1 cell (1310 dyn/cm). To analyze the origins of the stability of these systems and to predict energetics for more complex, less-ordered systems, we develop a model in which the surface energy is partitioned into contributions from seven different types of atom environments. This analysis is used to predict the surface energy of larger DAS structures (including their asymptotic behavior for very large unit cells) and to study the energetics of the sequential size change (SSC) model proposed by Shimada and Tochihara for the observed dynamical reconstruction of the Si(111) 1x1 structure. We obtain an energy barrier at the 2x2 cell size and confirm that the 7x7 regular stage of the SSC model (corresponding to the DAS 7x7 reconstruction) provides the highest energy reduction per unit cell with respect to the unreconstructed Si111 1x1 surface.     

A realistic double many-body expansion potential energy surface for SO2(X1A') from a multiproperty fit to accurate ab initio energies and vibrational levels

A single-valued double many-body expansion potential energy surface (DMBE I) recently obtained for the ground electronic state of the sulfur dioxide molecule by fitting correlated ab initio energies suitably corrected by scaling the dynamical correlation energy is now refined by fitting simultaneously available spectroscopic levels up to 6886 cm(-1) above the minimum. The topographical features of the novel potential energy surface (DMBE II) are examined in detail, and the method is emphasized as a robust route to fit together state-of-the-art theoretical calculations and spectroscopic measurements using a single fully dimensional potential form.     

A six-dimensional H(2)-H(2) potential energy surface for bound state spectroscopy

We present a six-dimensional potential energy surface for the (H(2))(2) dimer based on coupled-cluster electronic structure calculations employing large atom-centered Gaussian basis sets and a small set of midbond functions at the dimer's center of mass. The surface is intended to describe accurately the bound and quasibound states of the dimers (H(2))(2), (D(2))(2), and H(2)-D(2) that correlate with H(2) or D(2) monomers in the rovibrational levels (v,j)=(0,0), (0,2), (1,0), and (1,2). We employ a close-coupled approach to compute the energies of these bound and quasibound dimer states using our potential energy surface, and compare the computed energies for infrared and Raman transitions involving these states with experimentally measured transition energies. We use four of the experimentally measured dimer transition energies to make two empirical adjustments to the ab initio potential energy surface; the adjusted surface gives computed transition energies for 56 experimentally observed transitions that agree with experiment to within 0.036 cm(-1). For 26 of the 56 transitions, the agreement between the computed and measured transition energies is within the quoted experimental uncertainty. Finally, we use our potential energy surface to predict the energies of another 34 not-yet-observed infrared and Raman transitions for the three dimers.     

Rotationally resolved S1-S0 electronic spectra of 2,6-diaminopyridine: a four-fold barrier problem

A comparison of the electronic properties of the nitrogen-containing rings aniline, 2-aminopyridine, and 2,6-diaminopyridine (26DAP) shows that the potential energy surface of the molecule is significantly affected as more nitrogen atoms are added to the system. High resolution, rotationally resolved spectra of four vibrational bands in the S(1)-S(0) electronic transition of 26DAP were obtained in order to explain these changes. The zig-zagging inertial defects point to a double minimum excited state potential energy surface along the coupled amino group inversion vibrational mode, which becomes a four-fold well (and barrier) problem when the existence of two nearly degenerate isomers is taken into account. Assuming that the molecules are in the lower energy, opposite-side configuration, ab initio calculations were performed using the MP2/6-31G** level of theory to create a potential energy surface modeling the simultaneous antisymmetric NH(2)-inversion mode. The calculated potential energy surface shows a ground electronic state barrier to simultaneous NH(2) inversion of ~220 cm(-1), and a fit to experimental vibrational energy level spacings and relative intensities produces an excited electronic state barrier of ~400 cm(-1). The ground state barrier is less than that in aniline, but the excited state barrier is larger.     

Theoretical investigation of the He-I2(E3Πg) ion-pair state: ab initio intermolecular potential and vibrational levels

We present a theoretical study on the potential energy surface and vibrational bound states of the E electronic excited state of the HeI(2) van der Waals system. The interaction energies are computed using accurate ab initio methods and large basis sets. Relativistic small-core effective core potentials in conjunction with a quintuple-zeta quality basis set are employed for the heavy iodine atoms in multireference configuration interaction calculations for the (3)A' and (3)A" states. For the representation of the potential energy surface we used a general interpolation technique for constructing potential surfaces from ab initio data based on the reproducing kernel Hilbert space method. The surface presents global and local minima for T-shaped configurations with well-depths of 33.2 and 4.6 cm(-1), respectively. Vibrational energies and states are computed through variational quantum mechanical calculations. We found that the binding energy of the HeI(2)(E) T-shaped isomer is 16.85 cm(-1), in excellent agreement with recent experimental measurements. In lieu of more experimental data we also report our predictions on higher vibrational levels and we analyze the influence of the underlying surface on them. This is the first attempt to represent the potential surface of such a highly excited electronic state of a van der Waals complex, and it demonstrates the capability of the ab initio technology to provide accurate results for carrying out reliable studies to model experimental data.     

Analysis of local anharmonicity using Gaussian model potentials and cartesian oscillator basis sets: example, HCN

This paper examines local anharmonic vibrations in molecules using an analysis that starts with an ab initio potential energy surface, fits a model potential constructed of Gaussian basis functions, and proceeds to a quantum mechanical analysis of the anharmonic modes using Cartesian harmonic oscillator basis functions in a variational calculation. The objective of this work is to suggest methods, with origins in nuclear and molecular (electronic) quantum mechanics, that should be useful for the accurate analysis of the local anharmonic motions of hydrogen, and perhaps other atoms or small molecular fragments, residing in molecularly complicated but otherwise harmonic environments.     

High-dimensional quantum dynamical study of the dissociation of H2 on Pd110

We report the first six-dimensional quantum dynamical study of the dissociative adsorption of H(2) on a (110) surface. We have performed quantum coupled-channel calculations for the system H(2)/Pd(110) based on a potential energy surface (PES) that was derived from ab initio electronic structure calculations. In particular, we have focused on the effects of the corrugation and anisotropy of the PES on the H(2) dissociation probability. Our results agree well with the available experimental data for the sticking probability as a function of the initial kinetic energy and the angle of incidence. Because of the coupling between the anisotropy and corrugation of the potential energy surface our calculations predict an unusual rotational heating and a rather small rotational alignment in desorption.     

Deactivation mechanism of the green fluorescent chromophore

We report time-resolved fluorescence data for the anion of p-hydroxybenzylidene dimethylimidazolinone (p-HBDI), a model chromophore of the green fluorescence protein, in viscous glycerol-water mixtures over a range of temperatures, T. The markedly nonexponential decay of the excited electronic state is interpreted with the aid of an inhomogeneous model possessing a Gaussian coordinate-dependent sink term. A nonlinear least-squares fitting routine enables us to achieve quantitative fits by adjusting a single activation parameter, which is found to depend linearly on 1/T. We derive an analytic expression for the absolute quantum yield, which is compared with the integrated steady-state fluorescence spectra. The microscopic origins of the model are discussed in terms of two-dimensional dynamics, coupling the phenyl-ring rotation to a swinging mode that brings this flexible molecule to the proximity of a conical intersection on its multidimensional potential energy surface.     

Conformational energy landscapes of liquid crystal molecules

The conformational energy landscape of the prototypical nematic liquid crystal 4-pentyl-4cyanobiphenyl (5CB) is studied using first principles computer modelling. It is found that the most favourable conformation occurs when the two constituent phenyl rings are inclined at an angle of 31 with respect to each other. Also, the orientation of the alkyl chain is found to have an important influence on the ring-ring torsional potential. We fit the energy surface of these coupled torsions to yield an accurate intramolecular potential for use in empirical modelling. To test the strength of the coupling between the alkyl tail and the phenyl rings and the cyano group, we also calculate potentials for the relative orientation of the phenyl rings in biphenyl and cyanobiphenyl (0CB). Our calculations are performed using density functional theory using pseudo-potentials and the generalized gradient approximation to exchange and correlation. The molecular electronic wavefunction is expanded in terms of a plane wave basis set. We compare our results with recent NMR and Gaussian-based quantum chemistry calculations where available.     

Cis-->trans, trans-->cis isomerizations and N-O bond dissociation of nitrous acid (HONO) on an ab initio potential surface obtained by novelty sampling and feed-forward neural network fitting

The isomerization and dissociation dynamics of HONO are investigated on an ab initio potential surface obtained by fitting the results of electronic structure calculations at 21 584 configurations by using previously described novelty sampling and feed-forward neural network (NN) methods. The electronic structure calculations are executed by using GAUSSIAN 98 with a 6-311G(d) basis set at the MP4(SDQ) level of accuracy. The average absolute error of the NN fits varies from 0.012 eV (1.22 kJ mol(-1)) to 0.017 eV (1.64 kJ mol(-1)). The average computation time for a HONO trajectory using a single NN surface is approximately 4.8 s. These computation times compare very favorably with those required by other methods primarily because the NN fitting needs to be executed only one time rather than at every integration point. If the average result obtained from a committee of NNs is employed at each point rather than a single NN, increased fitting accuracy can be achieved at the expense of increased computational requirements. In the present investigation, we find that a committee comprising five NN potentials reduces the average absolute interpolation error to 0.0111 eV (1.07 kJ mol(-1)). Cis-trans isomerization rates with total energy of 1.70 eV (including zero point energy) have been computed for a variety of different initial distributions of the internal energy. In contrast to results previously reported by using an empirical potential, where cis-->trans to trans-->cis rate coefficient ratios at 1.70 eV total energy were found to lie in the range of 2.0-12.9 depending on the vibration mode excited, these ratios on the ab initio NN potential lie in the range of 0.63-1.94. It is suggested that this result is a reflection of much larger intramode coupling terms present in the ab initio potential surface. A direct consequence of this increased coupling is a significant decrease in the mode specific rate enhancement when compared to results obtained by using empirical surfaces. All isomerizations are found to be first order in accordance with the results reported by using empirical potentials. The dissociation rate to NO+OH has been investigated at internal HONO energies of 3.10 and 3.30 eV for different distributions of this energy among the six vibrational modes of HONO. These dissociations are also found to be first order. The computed dissociation rate coefficients exhibit only modest mode specific rate enhancement that is significantly smaller than that obtained on an empirical surface because of the much larger mode couplings present on the ab initio surface.     

The molecular potential energy surface and vibrational energy levels of methyl fluoride. Part II

New analytical bending and stretching, ground electronic state, potential energy surfaces for CH(3)F are reported. The surfaces are expressed in bond-length, bond-angle internal coordinates. The four-dimensional stretching surface is an accurate, least squares fit to over 2000 symmetrically unique ab initio points calculated at the CCSD(T) level. Similarly, the five-dimensional bending surface is a fit to over 1200 symmetrically unique ab initio points. This is an important first stage towards a full nine-dimensional potential energy surface for the prototype CH(3)F molecule. Using these surfaces, highly excited stretching and (separately) bending vibrational energy levels of CH(3)F are calculated variationally using a finite basis representation method. The method uses the exact vibrational kinetic energy operator derived for XY(3)Z systems by Manson and Law (preceding paper, Part I, Phys. Chem. Chem. Phys., 2006, 8, DOI: 10.1039/b603106d). We use the full C(3v) symmetry and the computer codes are designed to use an arbitrary potential energy function. Ultimately, these results will be used to design a compact basis for fully coupled stretch-bend calculations of the vibrational energy levels of the CH(3)F system.     

Molecular dynamics investigation of ferrous-ferric electron transfer in a hydrolyzing aqueous solution: calculation of the pH dependence of the diabatic transfer barrier and the potential of mean force

We present a molecular model for ferrous-ferric electron transfer in an aqueous solution that accounts for electronic polarizability and exhibits spontaneous cation hydrolysis. An extended Lagrangian technique is introduced for carrying out calculations of electron-transfer barriers in polarizable systems. The model predicts that the diabatic barrier to electron transfer increases with increasing pH, due to stabilization of the Fe3+ by fluctuations in the number of hydroxide ions in its first coordination sphere, in much the same way as the barrier would increase with increasing dielectric constant in the Marcus theory. We have also calculated the effect of pH on the potential of mean force between two hydrolyzing ions in aqueous solution. As expected, increasing pH reduces the potential of mean force between the ferrous and ferric ions in the model system. The magnitudes of the predicted increase in diabatic transfer barrier and the predicted decrease in the potential of mean force nearly cancel each other at the canonical transfer distance of 0.55 nm. Even though hydrolysis is allowed in our calculations, the distribution of reorganization energies has only one maximum and is Gaussian to an excellent approximation, giving a harmonic free energy surface in the reorganization energy F(DeltaE) with a single minimum. There is thus a surprising amount of overlap in electron-transfer reorganization energies for Fe(2+)-Fe(H2O)6(3+), Fe(2+)-Fe(OH)(H2O)5(2+), and Fe(2+)-Fe(OH)2(H2O)+ couples, indicating that fluctuations in hydrolysis state can be viewed on a continuum with other solvent contributions to the reorganization energy. There appears to be little justification for thinking of the transfer rate as arising from the contributions of different hydrolysis states. Electronic structure calculations indicate that Fe(H2O)6(2+)-Fe(OH)n(H2O)(6-n)(3-n)+ complexes interacting through H3O2- bridges do not have large electronic couplings.     

Potential energy surface determinations for nonrigid molecules: Application to acetone

In the present work, the problem of the determination of the potential energy surface for nonrigid molecules is examined in the case of the double rotation of the methyl groups in acetone. From the symmetry adapted functional form for the potential, the minimum number of configurations to be calculated is deduced in order to have a reliable surface. With this consideration in mind, the potential energy surface of acetone is determined in some Hartree–Fock semiempirical (CNDO /2) and ab initio procedures with different standard basis sets. In addition, ab initio calculations are performed using different sets of floating Gaussian functions in order to introduce some polarization effects in the wave function. Finally, the influence of the electronic correlation effects in the barrier height, and the role of the possible relaxation of the structure during the rotation is discussed.     

The photodissociation of ozone in the Hartley band: a theoretical analysis

Three-dimensional diabatic potential energy surfaces for the lowest four electronic states of ozone with 1A' symmetry-termed X, A, B, and R-are constructed from electronic structure calculations. The diabatization is performed by reassigning corresponding energy points. Although approximate, these diabatic potential energy surfaces allow one to study the uv photodissociation of ozone on a level of theory not possible before. In the present work photoexcitation in the Hartley band and subsequent dissociation into the singlet channel, O3X+hnu-->O(1D)+O2(a 1Deltag), are investigated by means of quantum mechanical and classical trajectory calculations using the diabatic potential energy surface of the B state. The calculated low-resolution absorption spectrum as well as the vibrational and rotational state distributions of O2(a 1Deltag) are in good agreement with available experimental results.