Effects of zero point vibration on the reaction dynamics of water dimer cations following ionization

Reactions of water dimer cation following ionization have been investigated by means of a direct ab initio molecular dynamics method. In particular, the effects of zero point vibration and zero point energy (ZPE) on the reaction mechanism were considered in this work. Trajectories were run on two electronic potential energy surfaces (PESs) of : ground state (²A″‐like state) and the first excited state (²A′ ‐ like state). All trajectories on the ground‐state PES lead to the proton‐transferred product: H₂O⁺(Wd)‐H₂O(Wa) → OH(Wd)‐H₃O⁺(Wa), where Wd and Wa refer to the proton donor and acceptor water molecules, respectively. Time of proton transfer (PT) varied widely from 15 to 40 fs (average time of PT = 30.9 fs). The trajectories on the excited‐state PES gave two products: an intermediate complex with a face‐to‐face structure (H₂O‐OH₂)⁺ and a PT product. However, the proton was transferred to the opposite direction, and the reverse PT was found on the excited‐state PES: H₂O(Wd)‐H₂O⁺ (Wa) → H₃O⁺(Wd)‐OH(Wa). This difference occurred because the ionizing water molecule in the dimer switched between the ground and excited states. The reaction mechanism of and the effects of ZPE are discussed on the basis of the results. © 2017 Wiley Periodicals, Inc.     

On the origin of the dynamical threshold for collision-induced dissociation processes

Collision-induced dissociation from the ground vibration-rotation state of the reactant diatom is studied in the systems He + H₂, H + H₂, and He + H⁺₂ by quasiclassical trajectory calculations, using ab initio potential energy surfaces. The dependence of the dynamical threshold values on the shape of the potential energy surface is discussed.     

High‐level ab initio study of the N+(3P)+SH2 reactions in the gas phase: Role of spin‐forbidden pathways

The singlet and triplet potential energy surfaces involved in N⁺+SH₂ reactions have been explored using high‐level ab initio techniques. The geometries of the stationary points were optimized at the QCISD/6‐311G(df,p) level. The final energies were obtained in CCSD(T)/6‐311+G(3df,2p) single‐point calculations. The results obtained show that, although the N⁺(¹D)+SH₂ entrance channel is higher in energy than the N⁺(³P)+SH₂ one, most of the [H₂, S, N]⁺ singlet state cations are lower in energy than the corresponding triplets, due to their different bonding characteristics. Both singlet and triplet potential energy surfaces are quite close each other, and crossover between them can occur. The minimum energy crossing points were located by means of CASSCF(6,5) calculations. The spin‐orbit couplings show that the transition probability from the triplet to the singlet potential energy surface is significantly large. One of the most important consequences is that some of the products of the reaction, such as SH⁺, can be formed in typical spin‐forbidden processes. Since all the relevant structures along these pathways are much lower in energy than the reactants, this mechanism should be accessible even at low impact energies and therefore could be important in processes taking place in interstellar media. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Variational transition‐state theory study of the atmospheric reaction OH + O3 → HO2 + O2

We report variational transition‐state theory calculations for the OH + O₃→ HO₂ + O₂ reaction based on the recently reported double many‐body expansion potential energy surface for ground‐state HO₄ [Chem Phys Lett 2000, 331, 474]. The barrier height of 1.884 kcal mol^?1 is comparable to the value of 1.77–2.0 kcal mol^?1 suggested by experimental measurements, both much smaller than the value of 2.16–5.11 kcal mol^?1 predicted by previous ab initio calculations. The calculated rate constant shows good agreement with available experimental results and a previous theoretical dynamics prediction, thus implying that the previous ab initio calculations will significantly underestimate the rate constant. Variational and tunneling effects are found to be negligible over the temperature range 100–2000 K. The O₁? O₂ bond is shown to be spectator like during the reactive process, which confirms a previous theoretical dynamics prediction. © 2007 Wiley Periodicals, Inc. 39: 148–153, 2007     

Ab initio potential energy surface of the lowest triplet BH+2 collision complex

Important parts of the lowest triplet potential energy surface of the BH⁺₂ molecular system are studied at the ab initio level. Obtained topological parameters of this surface and regions of its crossing with excited surfaces are compared with results of the DIM method and with previous ab initio studies.     

On the ionic states of vinylidene and acetylene

Electronic states of C₂H⁺₂ cation (with acetylenic and vinylidenic structure) are investigated by ab initio SCF and PNO CEPA calculations. The lowest excited state of the acetylene cation is calculated to be ⁴A₂ in a strongly bent cis conformation, and a qualitative explanation is given for the lack of detectable emission from theò∑_g⁺ state. Only the radical anion with the vinylidenic structure is bound.     

Thermal and state-selected rate constant calculations for O(3p) + H2 → OH + H and isotopic analogs

We present a new parametrization (based on ab initio calculations) of the bending potentials for the two lowest potential energy surfaces of the reaction O(³P) + H₂, and we use it for rate constant calculations by variational transition-state theory with multidimensional semiclassical tunneling corrections. We present results for the temperature range 250–2400 K for both the rate constants and the intermolecular kinetic isotope effects for the reactions of O(³P) with D₂ and HD. In general, the calculated rate constants for the thermal reactions are in excellent agreement with available experiments. We also calculate the enhancement effect for exciting H₂ to the first excited vibrational state. The calculations also provide information on which aspects of the potential energy surfaces are important for determining the predicted rate constants.     

Assessment of semiempirical methods for the photoisomerisation of a protonated Schiff base

The potential energy surfaces and non-adiabatic dynamics of the C₅H₆NH ₂ ⁺ protonated Schiff base (PSB3) have been investigated using the OM2 semiempirical Hamiltonian with GUGA configuration interaction. Three approaches to selecting the GUGA-CI active space are evaluated using closed-shell and open-shell molecular orbitals. Energy minima and minimum energy crossing points (MECPs) have been compared with ab initio CASSCF and CASPT2 results. Only the open-shell calculations give a qualitatively correct MECP. Minimum energy path (MEP) calculations demonstrate that a minimal active space gives a barrierless path from the planar S₁ minimum to the ground state, whereas larger active spaces result in a small barrier to torsional motion. Surface hopping dynamics calculations indicate that this barrier induces bi-exponential dynamics. The comparable CASSCF S₁ energy surface is barrierless, but the CASPT2 surface features an energy plateau, which may also lead to more complex dynamics.     

Noncollinear spin density functional theory for spin‐frustrated and spin‐degenerate systems

We have implemented ab initio linear combinations of Gaussian‐type orbital calculations with generalized localized spin density approximation (GLSDA) for a dimer of equilateral H₃ as a model of the noncollinear magnetic clusters. It has been found that the GLSDA solution with the three‐dimensional noncollinear spin structure is, contrary to prior band calculations by other groups, the ground state near the O_h conformation. Further computational results are compared to that of ab initio generalized Hartree–Fock. The difference between them and the influence of the correlation correction were discussed. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

A quantum chemical study of the infrared absorption intensities of the isolectronic C3v systems NH3, H3O+ and CH−3

The dipole moment derivatives and the infrared absorption intensities for the isoelectronic, isostructural species NH₃, H₃O⁺ and CH⁻₃, calculated by ab initio quantum methods within the double harmonic approximation, are reported. The calculations were performed at the SCF, CI and CPA″ levels of theory using basis sets of triple zeta+two polarization functions quality. For the ions H₃O⁺ and CH⁻₃, in the absence of adequate experimental information, the calculations are fully ab initio, since the equilibrium geometries as well as the force constants had to be computed. The applicability of the harmonic treatment to systems with inversion potentials is discussed, especially with regard to H₃O⁺. The dipole moment derivatives of the three systems show interesting, regular trends in accordance with the amount of electronic charge associated with the hydrogen atoms.     

Theoretical studies on structures and electronic spectra of linear carbon chains C2nH+ (n = 1−5)

The density functional theory (DFT) and the complete active space self‐consistent‐field (CASSCF) method have been used for full geometry optimization of carbon chains C_2nH⁺ (n = 1–5) in their ground states and selected excited states, respectively. Calculations show that C_2nH⁺ (n = 1–5) have stable linear structures with the ground state of X³Π for C₂H⁺ or X³Σ^? for other species. The excited‐state properties of C_2nH⁺ have been investigated by the multiconfigurational second‐order perturbation theory (CASPT2), and predicted vertical excitation energies show good agreement with the available experimental values. On the basis of our calculations, the unsolved observed bands in previous experiments have been interpreted. CASSCF/CASPT2 calculations also have been used to explore the vertical emission energy of selected low‐lying states in C_2nH⁺ (n = 1–5). Present results indicate that the predicted vertical excitation and emission energies of C_2nH⁺ have similar size dependences, and they gradually decrease as the chain size increases. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

On the Reaction Mechanism of the Complete Intermolecular O2 Transfer between Mononuclear Nickel and Manganese Complexes with Macrocyclic Ligands

The recently described intermolecular O₂ transfer between the side‐on Ni‐O₂ complex [(12‐TMC)Ni‐O₂]⁺ and the manganese complex [(14‐TMC)Mn]²⁺, where 12‐TMC and 14‐TMC are 12‐ and 14‐membered macrocyclic ligands, 12‐TMC=1,4,7,10‐tetramethyl‐1,4,7,10‐tetraazacyclododecane and 14‐TMC=1,4,8,11‐tetramethyl‐1,4,8,11‐tetraazacyclotetradecane, is studied by means of DFT methods. B3LYP calculations including long‐range corrections and solvent effects are performed to elucidate the mechanism. The potential energy surfaces (PESs) compatible with different electronic states of the reactants have been analyzed. The calculations confirm a two‐step reaction, with a first rate‐determining bimolecular step and predict the exothermic character of the global process. The relative stability of the products and the reverse barrier are in line with the fact that no reverse reaction is experimentally observed. An intermediate with a μ‐η¹:η¹‐O₂ coordination and two transition states are identified on the triplet PES, slightly below the corresponding stationary points of the quintet PES, suggesting an intersystem crossing before the first transition state. The calculated activation parameters and the relative energies of the two transition sates and the products are in very good agreement with the experimental data. The calculations suggest that a superoxide anion is transferred during the reaction.     

Nonempirical Investigation of Equilibrium Geometry and Electronic Structure of Kynurenine and 3-Hydroxykynurenine Molecule

By Hartree-Fock-Roothaan method with complete geometry optimization in the basis 6-31G* ab initio calculations of equilibrium geometry and electronic structure were performed for kynurenine C₁₀H₁₂N₂O₃ and 3-hydroxykynurenine C₁₀H₁₂N₂O₃ molecules in the singlet ground state and the first triplet excited state. The molecules in the triplet state can react at the oxygen of the carbonyl group adjacent to the aromatic ring by quite different pathway compared to the molecules in the ground singlet state.     

Diatomics-in-molecules calculation of the potential energy surfaces relevant to penning ionization in the He(2 3S)H2 system

Potential energy surfaces and the autoionization width for the Penning ionization transition He(2 ³S) + H₂ → He + H⁺₂ + e^? have been calculated using the DIM method. The surfaces compare favourably with the existing ab initio calculations, and the approximation to the autoioinization width appear to be reasonable.     

Magnetic Properties of a Terbium–[1]Ferrocenophane Complex: Analogies between Lanthanide–Ferrocenophane and Lanthanide–Bis‐phthalocyanine Complexes

A rare example of an organometallic terbium single‐ion magnet is reported. A Tb³⁺–[1]ferrocenophane complex displays a larger barrier to magnetization reversal than its isostructural Dy³⁺ analogue, which is reminiscent of trends observed for lanthanide–bis‐phthalocyanine complexes. Detailed ab initio calculations support the experimental observations and suggest a significantly larger ground‐state stabilization for the non‐Kramers ion Tb³⁺ in the Tb complex than for the Kramers‐ion Dy³⁺ in the Dy complex.     

Potential Energy Surfaces for the Na(3p)N2, Na(3p)HCN and Na(3p)C2H2 Complexes: Ab Initio Calculations and Model Considerations

Highly correlated ab initio calculations for the potential‐energy surfaces of the systems Na−N₂, Na−HCN, and Na−C₂H₂ have been performed, where the Na atom is in the 3s ground state or in the 3p excited state, and the molecules are kept fixed in their equilibrium configurations. The purpose of these calculations was to enable computer simulation of optical‐scattering experiments. In this paper, model considerations are presented that allow easy qualitative understanding of the shapes of the calculated surfaces.     

Origin of the Magnetic Anisotropy in Heptacoordinate NiII and CoII Complexes

The nature and magnitude of the magnetic anisotropy of heptacoordinate mononuclear Ni^II and Co^II complexes were investigated by a combination of experiment and ab initio calculations. The zero‐field splitting (ZFS) parameters D of [Ni(H₂DAPBH)(H₂O)₂](NO₃)₂ ? 2 H₂O ( 1 ) and [Co(H₂DAPBH)(H₂O)(NO₃)](NO₃) [ 2 ; H₂DAPBH=2,6‐diacetylpyridine bis‐ (benzoyl hydrazone)] were determined by means of magnetization measurements and high‐field high‐frequency EPR spectroscopy. The negative D value, and hence an easy axis of magnetization, found for the Ni^II complex indicates stabilization of the highest M_S value of the S=1 ground spin state, while a large and positive D value, and hence an easy plane of magnetization, found for Co^II indicates stabilization of the M_S=±1/2 sublevels of the S=3/2 spin state. Ab initio calculations were performed to rationalize the magnitude and the sign of D, by elucidating the chemical parameters that govern the magnitude of the anisotropy in these complexes. The negative D value for the Ni^II complex is due largely to a first excited triplet state that is close in energy to the ground state. This relatively small energy gap between the ground and the first excited state is the result of a small energy difference between the d_xy and ${{\rm{d}}_{x^2 - y^2 } }$ orbitals owing to the pseudo‐pentagonal‐bipyramidal symmetry of the complex. For Co^II, all of the excited states contribute to a positive D value, which accounts for the large magnitude of the anisotropy for this complex.     

An ab initio SCF and CI study of ketene imine

The electronic structures of the ground and excited states of ketene imine (_H^HCCNH) have been studied by ab initio SCF and CI calculations. The nucleophilic nature of the β carbon with respect to nitrogen has been discussed using calculated electrostatic potentials and by calculated energy differences between the parent and protonated species. The electronically excited ¹A″ and ³A″ states are found to be almost degenerate.     

Fe and Ni AB initio effective potentials for use in molecular calculations

We present effective potentials to replace the Ar core electrons of Fe and Ni. These effective potentials are obtained from ab initio ground state wavefunctions of Fe and Ni and are tested by comparing with ab initio SCF calculations for excited states of Fe, Fe⁺, Fe²⁺, Fe³⁺, Ni, Ni⁺, Ni²⁺, and the FeH⁺ molecule.