The calculated infrared spectrum of Cl- H2O using a full dimensional ab initio potential surface and dipole moment surface

We report a full dimensional ab initio-based global potential energy surface (PES) and dipole moment surface (DMS) for Cl(-)H(2)O. Both surfaces are symmetric with respect to interchange of the H atoms. The PES is a fit to thousands of electronic energies calculated using the coupled-cluster method (CCSD(T)) with a moderately large basis (aug-cc-pVTZ). The infrared spectrum and vibrational dynamics are reported and compared to experiment. These results are analyzed by examination of wave function and the dipole surface.     

New ab initio potential energy surface and the vibration-rotation-tunneling levels of (H2O)2 and (D2O)2

We report a new full-dimensional potential energy surface (PES) for the water dimer, based on fitting energies at roughly 30,000 configurations obtained with the coupled-cluster single and double, and perturbative treatment of triple excitations method using an augmented, correlation consistent, polarized triple zeta basis set. A global dipole moment surface based on Moller-Plesset perturbation theory results at these configurations is also reported. The PES is used in rigorous quantum calculations of intermolecular vibrational frequencies, tunneling splittings, and rotational constants for (H2O)2 and (D2O)2, using the rigid monomer approximation. Agreement with experiment is excellent and is at the highest level reported to date. The validity of this approximation is examined by comparing tunneling barriers within that model with those from fully relaxed calculations.     

The infrared spectrum of cyclic-N3: theoretical prediction

We have carried out the first calculations of the infrared absorption spectrum of cyclic-N(3). Accurate vibrational energies and wave functions computed with incorporation of the geometric phase effect (via gauge theory) and using an ab initio potential energy surface were employed in this work. A sophisticated fully dimensional dipole moment function was constructed using accurate ab initio calculations and a three-dimensional-spline interpolation. Transformation of the dipole moment vector function from the reference frame associated with instantaneous principal axes of inertia to the laboratory-fixed reference frame was carried out using hyperspherical coordinates. We found that the permanent dipole moment of cyclic-N(3) in the ground vibrational state is relatively small (170 mD). The excited vibrational states show permanent dipole moments in the 10-25 mD range. The most intense part of the infrared absorption spectrum is observed in the deep infrared part of spectrum, 75-275?cm(-1), where five lines exhibit absolute absorption intensities in the range between 0.5 and 1.2 km/mol. These transitions correspond to excitation of the pseudorotational progression of states. Several unique spectroscopic features discussed in the paper should help to identify cyclic-N(3) in the laboratory.     

Ab initio potential energy and dipole moment surfaces for H5O2 +

Full-dimensional ab initio potential energy surface (PES) and dipole moment surface (DMS) are reported for H(5)O(2) (+). Tens of thousands of coupled-cluster [CCSD(T)] and second-order Moller-Plesset (MP2) calculations of electronic energies, using aug-cc-pVTZ basis, were done. The energies were fit very precisely in terms of all the internuclear distances, using standard least-square procedures, however, with a fitting basis that satisfies permutational symmetry with respect to like atoms. The H(5)O(2) (+) PES is a fit to 48 189 CCSD(T) energies, containing 7962 polynomial coefficients. The PES has a rms fitting error of 34.9 cm(-1) for the entire data set up to 110 000 cm(-1). This surface can describe various internal floppy motions, including the H atom exchanges, monomer inversions, and monomer torsions. First- and higher-order saddle points have been located on the surface and compared with available previous theoretical work. In addition, the PES dissociates correctly (and symmetrically) to H(2)O+H(3)O(+), with D(e)=11 923.8 cm(-1). Geometrical and vibrational properties of the monomer fragments are presented. The corresponding global DMS fit (MP2 based) involves 3844 polynomial coefficients and also dissociates correctly.     

Flexible, ab initio potential, and dipole moment surfaces for water. I. Tests and applications for clusters up to the 22-mer

We report full-dimensional, ab initio potential energy and dipole moment surfaces, denoted PES and DMS, respectively, for arbitrary numbers of water monomers. The PES is a sum of 1-, 2-, and 3-body potentials which can also be augmented by semiempirical long-range higher-body interactions. The 1-body potential is a spectroscopically accurate monomer potential, and the 2- and 3-body potentials are permutationally invariant fits to tens of thousands of CCSD(T)/aug-cc-pVTZ and MP2/aug-cc-pVTZ electronic energies, respectively. The DMS is a sum of 1- and 2-body DMS, which are covariant fits to tens of thousands MP2/aug-cc-pVTZ dipole moment data. We present the details of these new 2- and 3-body potentials and then extensive applications and tests of this PES are made to the structures, classical binding energies, and harmonic frequencies of water clusters up to the 22-mer. In addition, we report the dipole moment for these clusters at various minima and compare the results against available and new ab initio calculations.     

Ab initio potential energy and dipole moment surfaces of (H2O)2

A full-dimensional ab initio potential energy surface (PES) and dipole moment surface (DMS) are reported for the water dimer, (H2O)2. The CCSD(T)-PES is a very precise fit to 19,805 ab initio energies obtained with the coupled-cluster (CCSD(T)) method, using an aug-cc-pVTZ basis. The standard counterpoise correction was applied to approximately eliminate basis set superposition errors. The fit is based on an approach that incorporates the permutational symmetry of identical atoms [Huang, X.; Braams, B.; Bowman, J. M. J. Chem.Phys. 2005, 122, 044308]. The DMS is a fit to the dipole moment obtained with M?ller-Plesset (MP2) theory, using an aug-cc-pVTZ basis. The PES has an RMS fitting error of 31 cm(-1) for energies below 20,000 cm(-1), relative to the global minimum. This surface can describe various internal floppy motions, including various monomer inversions, and isomerization pathways. Ten characteristic stationary points have been located on the surface, four of which are transition structures and the rest are higher order saddle points. Their geometrical and vibrational properties are presented and compared with available previous theoretical work. The CCSD(T)-PES and MP2-DMS dissociate correctly (and symmetrically) to two H2O monomers, with D(e) = 1665.7 cm(-1) (19.93 kJ/mol). Accurate quantum calculations of the zero-point energy of the dimer (using diffusion Monte Carlo) and the monomers (using a vibrational configuration interaction approach) are reported, and these together with D(e) give a value of D0 of 1042 cm(-1) (12.44 kJ/mol). A best estimated value is 1130 cm(-1) (13.5 kJ/mol).     

Predicting the infrared transition intensities in the Ar-HF complex: the key role of the dipole moment surface accuracy

The method proposed earlier for the generation of the full-dimensional energy surface for van der Waals complexes [P. Jankowski, J. Chem. Phys. 121, 1655 (2004)] is used to obtain a fulldimensional dipole moment surface for the atom-diatom complex in calculations based on the coupled-cluster with single, double, and noniterative triple excitation approach and the aug-cc-pVQZ basis sets. This surface has been employed to calculate transition intensities of the infrared spectra of Ar-HF. Special attention has been paid to study the problem of relative intensities of the different bands which have not been properly predicted within the long-range models of the dipole moment [A. E. Thornley and J. M. Hutson, J. Chem. Phys. 101, 5578 (1994)]. The intensities calculated with the present dipole moment surface agree very well with the experimental data, which indicate that the short-range interactions significantly affect the dipole moment surface and the calculated intensities. To investigate the role of the accuracy of the dipole moment surface on infrared transition intensities in atom-diatom complexes, four models of increasing complexity are studied. Their performance is shown to strongly depend on the region of the interaction energy surface probed by the initial and final states of the individual transitions.     

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) (+).     

Six-dimensional potential energy surface of the dissociative chemisorption of HCl on Au(111) using neural networks

We constructed a six-dimensional potential energy surface(PES)for the dissociative chemisorption of HCl on Au(111)using the neural networks method based on roughly 70000 energies obtained from extensive density functional theory(DFT)calculations.The resulting PES is accurate and smooth,based on the small fitting errors and good agreement between the fitted PES and the direct DFT calculations.Time-dependent wave packet calculations show that the potential energy surface is very well converged with respect to the number of DFT data points,as well as to the fitting process.The dissociation probabilities of HCl initially in the ground rovibrational state from six-dimensional quantum dynamical calculations are quite diferent from the four-dimensional fixed-site calculations,indicating it is essential to perform full-dimensional quantum dynamical studies for the title molecule-surface interaction system.     

Investigating the relationship between infrared spectra of shared protons in different chemical environments: a comparison of protonated diglyme and protonated water dimer

The energetics, dynamics, and infrared spectroscopy of the shared proton in different chemical environments is investigated using molecular dynamics simulations. A three-dimensional potential energy surface (PES) suitable for describing proton transfer between an acceptor and a donor oxygen atom is combined with an all-atom force field to carry out reactive molecular dynamics simulations. The construction of the fully dimensional PES is inspired from the established mixed quantum mechanics/molecular mechanics treatment of larger systems. The "morphing potential" method is used to transform the generic PES for proton transfer along an O...H+...O motif into a three-dimensional PES for proton transfer in protonated diglyme. Using molecular dynamics simulations at finite temperature, the gas phase infrared spectra are calculated for both species from the Fourier transform of the dipole moment autocorrelation function. For protonated diglyme the modes involving the H+ motion are strongly mixed with other degrees of freedom. At low temperature, the O...H+...O asymmetric stretching vibration is found at 870 cm-1, whereas for H5O2+ this band is at 724 cm-1. As expected, the vibrational bands of protonated diglyme show no temperature dependence whereas for H5O2+ at T = 100 K the proton transfer mode is found at 830 cm-1, in good agreement with 861 cm-1 from very recent molecular dynamics simulations.     

Quantum vibrational analysis and infrared spectra of microhydrated sodium ions using an ab initio potential

We present a full-dimensional potential energy surface and a dipole moment surface (DMS) for hydrated sodium ion. These surfaces are based on an n-body expansion for both the potential energy and the dipole moment, truncated at the two-body level for the H(2)O-Na(+) interaction and also for the DMS. The water-water interaction is truncated at the three-body level. The new full-dimensional two-body H(2)O-Na(+) potential is a fit to roughly 20,000 coupled-cluster single double (triple)/aug-cc-pVTZ energies. Properties of this two-body potential and the potential describing (H(2)O)(n)Na(+) clusters, with n up to 4 are given. We then report anharmonic, coupled vibrational calculations with the "local-monomer model" to obtain infrared spectra and also 0 K radial distribution functions for these clusters. Some comparisons are made with the recent infrared predissociation spectroscopy experiments of Miller and Lisy [J. Am. Chem. Soc. 130, 15381 (2008).].     

An ab initio potential surface describing abstraction and exchange for H+CH4

We report an ab initio-based global potential energy surface for H+CH4 that describes the abstraction and exchange reactions. The PES, which is invariant with respect to any permutation of five H atoms, is a fit to 20,728 electronic energies calculated using the partially spin-restricted coupled-cluster method (RCCSD(T)) with a moderately large basis (aug-cc-pVTZ). A first set of quasiclassical trajectory calculations using this PES are reported for the H+CD4-->HD+CD3 reaction at collision energies of 0.65 and 1.52 eV and are compared to experiment and recent direct dynamics calculations done with density functional theory.     

Accurate ab initio structure, dissociation energy, and vibrational spectroscopy of the F(-)-CH4 anion complex

Accurate equilibrium structure, dissociation energy, global potential energy surface (PES), dipole moment surface (DMS), and the infrared vibrational spectrum in the 0-3000 cm(-1) range of the F(-)-CH4 anion complex have been obtained. The equilibrium electronic structure calculations employed second-order M?ller-Plesset perturbation theory (MP2) and coupled-cluster (CC) method up to single, double, triple, and perturbative quadruple excitations using the aug-cc-p(C)VXZ [X = 2(D), 3(T), 4(Q), and 5] correlation-consistent basis sets. The best equilibrium geometry has been obtained at the all-electron CCSD(T)/aug-cc-pCVQZ level of theory. The dissociation energy has been determined based on basis set extrapolation techniques within the focal-point analysis (FPA) approach considering (a) electron correlation beyond the all-electron CCSD(T) level, (b) relativistic effects, (c) diagonal Born-Oppenheimer corrections (DBOC), and (d) variationally computed zero-point vibrational energies. The final D(e) and D0 values are 2398 +/- 12 and 2280 +/- 20 cm(-1), respectively. The global PES and DMS have been computed at the frozen-core CCSD(T)/aug-cc-pVTZ and MP2/aug-cc-pVTZ levels of theory, respectively. Variational vibrational calculations have been performed for CH4 and F(-)-CH4 employing the vibrational configuration interaction (VCI) method as implemented in Multimode.     

Quasiclassical trajectory calculations of the OH+NO2 association reaction on a global potential energy surface

We report a full-dimensional potential energy surface (PES) for the OH+NO(2) reaction based on fitting more than 55,000 energies obtained with density functional theory-B3LYP6-311G(d,p) calculations. The PES is invariant with respect to permutation of like nuclei and describes all isomers of HOONO, HONO(2), and the fragments OH+NO(2) and HO(2)+NO. Detailed comparison of the structures, energies, and harmonic frequencies of various stationary points on the PES are made with previous and present high-level ab initio calculations. Two hydrogen-bond complexes are found on the PES and confirmed by new ab initio CASPT2 calculations. Quasiclassical trajectory calculations of the cross sections for ground rovibrational OH+NO(2) association reactions to form HOONO and HONO(2) are done using this PES. The cross section to form HOONO is larger than the one to form HONO(2) at low collision energies but the reverse is found at higher energies. The enhancement of the HOONO complex at low collision energies is shown to be due, in large part, to the transient formation of a H-bond complex, which decays preferentially to HOONO. The association cross sections are used to obtain rate constants for formation of HOONO and HONO(2) for the ground rovibrational states in the high-pressure limit.     

Structures,harmonic frequencies and infrared intensities of the dimers of H2O and H2S

The structures of the dimers (H₂O)₂, HOHSH₂, HSHOH₂ and (H₂S)₂ are obtained at the SCF level. Vibrational frequencies, dipole moment derivatives and infrared intensities are calculated for the inter- and intra-molecular modes using analytic derivative methods.     

Quantum control of molecular chirality: optical isomerization of difluorobenzo[c]phenanthrene

The results of a theoretical study are presented on quantum control of a chiral exchange reaction of a polyatomic molecule by using infrared laser pulses. Difluorobenzo[c]phenanthrene was chosen to be the simplest model for its helical chirality exchange reaction. This molecule has two stable configurations: M and P forms. From the viewpoint of chemical reaction dynamics, isomerization is regarded as the movement of one of the two representative points that initially correspond to the two forms to the position of the other representative point, while the other representative point remains in its initial position. The ground-state potential energy surface and dipole moment functions required to control this reaction were evaluated at the MP2/6-31+G(d,p) and MP2/TZV+(d,p) levels of molecular orbital (MO) theory. An effective potential energy surface (PES) that is a function of twisting motion of the benzene rings and wagging motion of the CF(2) group was constructed on the basis of the MO results. An analytical expression for the effective PES and that for the dipole moment functions were prepared to make the isomerization control tractable. A quantum control method in a classical way was applied to the isomerization of preoriented difluorobenzo[c]phenanthrene in low temperature limits. The time evolution of the representative point of the M form and that of the P form are separately evaluated to determine the optimal laser fields. This means that the laser control produces pure helical enantiomers from a racemic mixture. Representative points are replaced by the corresponding nuclear wave packets in this treatment. The derived control laser field consists of two linearly polarized E(x)() and E(z)() components that are perpendicular to each other. These components are pi-phase-shifted when the representative point is in the transition-state regions. Under the irradiation of this laser pulse, one of the two representative points of the isomerization is transferred to the target position along the intrinsic reaction path between the enantiomers, while the other representative point remains in its initial potential well. This results in one-way isomerization control, that is, the M(P) to P(M) form. The isomerization is completed with yields of ca. 70% within a few picoseconds. Temporal behaviors of the nuclear wave packet whose center corresponds to the representative point are drawn to see how the desired chiral exchange reaction proceeds in the presence of the control field, while its reverse process is suppressed.     

The predicted spectrum of the hypermetallic molecule MgOMg

The present study of MgOMg is a continuation of our theoretical work on Group 2 M(2)O hypermetallic oxides. Previous ab initio calculations have shown that MgOMg has a linear (1)Σ(g)+ ground electronic state and a very low lying first excited triplet electronic state that is also linear; the triplet state has (3)Σ(u)+ symmetry. No gas phase spectrum of this molecule has been assigned, and here we simulate the infrared absorption spectrum for both states. We calculate the three-dimensional potential energy surface, and the electric dipole moment surfaces, of each of the two states using a multireference configuration interaction (MRCISD) approach based on full-valence complete active space self-consistent field (FV-CASSCF) wavefunctions with a cc-pCVQZ basis set. A variational MORBID calculation using our potential energy and dipole moment surfaces is performed to determine rovibrational term values and to simulate the infrared absorption spectrum of the two states. We also calculate the dipole polarizability of both states at their equilibrium geometry in order to assist in the interpretation of future beam deflection experiments. Finally, in order to assist in the analysis of the electronic spectrum, we calculate the vertical excitation energies, and electric dipole transition matrix elements, for six excited singlet states and five excited triplet states using the state-average full valence CASSCF-MRCISD/aug-cc-pCVQZ procedure.     

Theory of the photodissociation of ozone in the Hartley continuum: potential energy surfaces, conical intersections, and photodissociation dynamics

Ab initio potential energy and transition dipole moment surfaces are presented for the five lowest singlet even symmetry electronic states of ozone. The surfaces are calculated using the complete active space self consistent field method followed by contracted multireference configuration interaction (MRCI) calculations. A slightly reduced augmented correlation consistent valence triple-zeta orbital basis set is used. The ground and excited state energies of the molecule have been computed at 9282 separate nuclear geometries. Cuts through the potential energy surfaces, which pass through the geometry of the minimum of the ground electronic state, show several closely avoided crossings. Close examination, and higher level calculations, very strongly suggests that some of these seemingly avoided crossings are in fact associated with non-symmetry related conical intersections. Diabatic potential energy and transition dipole moment surfaces are created from the computed ab initio adiabatic MRCI energies and transition dipole moments. The transition dipole moment connecting the ground electronic state to the diabatic B state surface is by far the strongest. Vibrational-rotational wavefunctions and energies are computed using the ground electronic state. The energy level separations compare well with experimentally determined values. The ground vibrational state wavefunction is then used, together with the diabatic B<--X transition dipole moment surface, to form an initial wavepacket. The analysis of the time-dependent quantum dynamics of this wavepacket provides the total and partial photodissociation cross sections for the system. Both the total absorption cross section and the predicted product quantum state distributions compare well with experimental observations. A discussion is also given as to how the observed alternation in product diatom rotational state populations might be explained.     

Negative surface potential produced by self-assembled monolayers of helix peptides oriented vertically to a surface

Self-assembled monolayers (SAMs) of helix peptides oriented vertically to a gold surface were prepared and the surface potential measured using the Kelvin technique up to 140°C. Negative surface potentials of a few hundred millivolts were observed for the helix peptide SAMs, indicating the occurrence of the large dipole moment of the helices directing toward the surface. The longer the helix peptide, the larger was the negative surface potential obtained. The absolute value of the surface potential decreased with increase in temperature due to thermal perturbation in the helical structure. However, Fourier transform infrared reflection–absorption spectroscopy revealed that perturbation is not significant and the α-helical conformation is stable even at 140°C.     

Dipole active vibrations and dipole moments of N2 and O2 physisorbed on a metal surface

We have, in infrared reflection absorption measurements, observed narrow dipole active absorption lines associated with the fundamental internal vibrational transitions of N(2) and O(2) physisorbed at 30 K on the chemically inert Pt(111)(1 x 1)H surface. Such transitions are forbidden for free homonuclear molecules and become dipole active at a metal surface due to polarization induced surface dipole moments. The measurements show that the internal stretch vibration frequencies are lowered by 7-8 cm(-1) relative to the gas phase values. The measured static and dynamic dipole moments are in the ranges of 0.06-0.07 and 0.001-0.002 D, respectively. We find that good estimates of the induced dynamic as well as the static dipole moments can in general be obtained from a van der Waals model but that the ratios of the measured static and dynamic moments indicates a need for a refinement of the dipole moment function.