Effect of reagent vibrational excitation and isotope substitution on the stereo-dynamics of the Ba+HF→BaF+H reaction

Based on an extended London-Eyring-Polanyi-Sato (LEPS) potential energy surface (PES), the Ba + HF reaction has been studied by the quasi-classical trajectory (QCT) method. The reaction integral cross section as a function of collision energy for the Ba + HF → BaF + H reaction is presented and the influence of isotope substitution on the differential cross sections (DCSs) and alignments of the product's rotational angular momentum have also been studied. The results suggest that the integral cross sections increase with increasing collision energy, and the vibrational excitation of the reagent has great influence on the DCS. In addition, the product's rotational polarization is very strong as a result of heavy-heavy-light (HHL) mass combination, and the distinct effect of isotope substitution on the stereodynamics is also revealed.     

Energy eigenvalues from an analytical transfer matrix method

A detailed procedure based on an analytical transfer matrix method is presented to solve bound-state problems. The derivation is strict and complete. The energy eigenvalues for an arbitrary one-dimensional potential can be obtained by the method. The anharmonic oscillator potential and the rational potential are two important examples. Checked by numerical techniques, the results for the two potentials by the present method are proven to be exact and reliable.     

Stereodynamics of the O（3P） with H2（D2） （ν=0,j=0） reaction

Quasi-classical trajectory theory is used to study the reaction of O(3 P) with H 2 (D 2) based on the ground 3 A″ potential energy surface (PES).The reaction cross section of the reaction O+H 2 →OH+H is in excellent agreement with the previous result.Vector correlations,product rotational alignment parameters P 2 (j · k) and several polarizeddependent differential cross sections are further calculated for the reaction.The product polarization distribution exhibits different characteristics that can be ascribed to different motion paths on the PES,arising from various collision energies or mass factors.     

A CLASSICAL TRAJECTORY STUDY BASED ON MBE POTENTIAL ENERGY FUNCTIONS—REACTIVE COLLISIONS OF C＋CH

Classical trajectories have been calculated for the reaction C CH→C2 H using many-body expansion(MBE) potential energy function of C2H.The present work gives the reactive cross sections at relative kinetic energy of 1.0,7.5,13.0,20.0and 30.0 kcal,From which it is found that this is a reaction without energy thres-hold.     

Analytical Potential Energy Function for Doubly Charged Ions O2+2, Be2+2 and He2+2

A new analytical potential function is proposed in our preceding paper as which can be used to describe the potential curves for doubly charged diatomic ions with both potential minimum and maximum where ρ=R-Rmin, if R=Rmin, ρ=0, if R=Rmax, ρ=Rmax-Rmin. The parameters a1～a4 are obtained by a least-squares fit based on the ab initio calculations, the values of parameters Rmin and Rmax are fixed from the same calculations. The application of this function is quite flexibly, for it can be also used to describe the potential curves without any potential minimum and maximum, i.e., no stationary point, where ρ=R. The present paper is to examine the proposed potential function form further by describing potential curves of O2+2, Be2+2 and He2+2. The calculations using the theoretical method CID with basis set 6-311++G* have showed the Rmin() and Rmax() to be 1.045 and 1.60, 2.141 and 3.50, and 0.736 and 1.05 for the ground state X1Σ+g of O2+2, Be2+2 and He2+2, respectively, which are in good agreement with reference data. The analytical potential function for the ground state of these doubly charged diatomic ions have been successfully derived for the first time     

Potential energy curves and spectroscopic properties of X2∑ and A2П states of 13C14N

The potential energy curves(PECs) of X2+Σand A2Π states of the CN molecule have been calculated with the multireference configuration interaction method and the aug-cc-pwCV5Z basis set. Based on the PECs, all of the vibrational and rotational levels of the13C14N molecule are obtained by solving the Schro¨dinger equation of the molecular nuclear motion.The spectroscopic parameters are determined by fitting the Dunham coefficients with the levels. Both the levels and the spectroscopic parameters are in good qualitative agreement with the experimental data available. The analytical potential energy functions are also deduced from the calculated PECs. The present results can provide a helpful reference for future spectroscopy experiments or dynamical calculations of the molecule.     

Accurate calculation of the potential energy curve and spectroscopic parameters of X~2Σ~+ state of ~(12)Mg~1H

High level calculations on the ground state of12Mg1 H molecule have been performed using multi-reference configuration interaction(MRCI) method with the Davidson modification. The core–valence correlation and scalar relativistic corrections are included into the present calculations at the same time. The potential energy curve(PEC) of the ground state, all of the vibrational levels and spectroscopic parameters are fitted. The results show that the levels and spectroscopic parameters are in good agreement with the available experimental data. The analytical potential energy function(APEF) is also deduced from the calculated PEC using the Murrell–Sorbie(M–S) potential function. The present results can provide a helpful reference for the future spectroscopic experiments or dynamical calculations of the molecule.     

Dynamics of the CH4+ O(~3 P) → CH_3(ν= 0) + OH(■= 0) reaction

The dynamics of the ground-state reaction of CH_4+ O(~3P) → CH_3(ν = 0) + OH( ν= 0) have attracted a great deal of attention both theoretically and experimentally. This rapid communication represents extensive quasi-classical trajectory calculations of the vibrational distributions on a unique full-dimensional ab initio potential energy surface for the title reaction, at the collision energy of relevance to previous crossed molecular beam experiments. The surface is constructed using the all electrons coupled-cluster singles and doubles approach plus quasi-perturbative triple excitations with optimized basis sets. A modified Shepard interpolation method is also employed for the construction. Good agreement between our calculations and the available experimental results has been achieved, opening the door for accurate dynamics on this surface.     

A new many-body potential with the second-moment approximation of tight-binding scheme for Hafnium

In this work,we develop a new many-body potential for alpha-hafnium(α-Hf)based on the second moment approximation of tight-binding(TB-SMA)theory by introducing an additional Heaviside step function into the potential model and a new analytical scheme of density function.All the parameters of the new potential have been systematically evaluated by fitting to ground-state properties including cohesive energy,lattice constants,elastic constants,vacancy formation energy,structure stability and equation of state.By using the present model,the melting point,melt heat,thermal expansion coefficient,point defects,and low-index surface energies ofα-Hf were calculated through molecular dynamics simulations.Comparing with experiment observations from others,it is shown that these properties can be reproduced reasonably by the present model,some results being more consistent to the experimental data than those by previous suggested models.This indicates that this work is sutiable in TB-SMA potential for hexagonal close packed metals.     

Quantum stereodynamics study for the reaction F + HD

This paper studies the quantum stereodynamics of the F + HD(υ = 0, j = 0) → HD + F/HF + D reaction at the collision energies of 0.52 and 0.87~kcal/mol. The quantum scattering calculations, based on Stark-Werner potential energy surfaces, show that the differential cross sections for the HF(υ' = 2) + D and DF(υ' = 3) + H channels are consistent with the recent theoretical results. Furthermore, the product rotational angular momentum orientation and alignment have been determined for some selected rovibrational states of the HF + D and DF + H channels.     

The equivalent potential of water molecules for electronic structure of lysine

In order to get more reliable electronic structures of proteins in aqueous solution, it is necessary to construct a potential of water molecules for protein’s electronic structure calculation. The lysine is a hydrophilic amino acid. It is positively charged (Lys+) in neutral water solution. The first-principles, all-electron, ab initio calcula-tions, based on the density functional theory, have been performed to construct such an equivalent potential of water molecules for lysine (Lys+). The process consists of three parts. First, the electronic structure of the cluster containing Lys+ and water molecules is calculated. By adjusting the positions of water molecules, the geometric structure of the cluster having minimum total energy is determined. Then, based on the structure, the electronic structure of Lys+ with the potential of water molecules is calculated using the self-consistent cluster-embedding (SCCE) method. Finally, the electronic structure of Lys+ with the potential of dipoles is calculated. The dipoles are adjusted so that the electronic structure of Lys+ with the potential of dipoles is close to that of water molecules. Thus the equivalent potential of water molecules for the electronic structure of lysine is obtained. The major effect of water molecules on lysine’s electronic structure is raising the occupied eigenvalues about 0.5032 eV, and broadening energy gap 89%. The effect of water molecules on the electronic structure of lysine can be simulated by dipoles potential.     

The characteristics of O＋HD （v = 0, j = 0） reaction dynamics

The analytical potential energy function of HDO is constructed at first using the many-body expansion method.The reaction dynamics of O+HD(v = 0,j = 0) in five product channels are all studied by quasi-classical trajectory(QCT) method.The results show that the long-lived complex compound HDO is the dominant product at low collision energy.With increasing collision energy,O+HD → OH+D and O+HD → OD+H exchange reactions will occur with remarkable characteristics,such as near threshold energies,different reaction probabilities,and different reaction cross sections,implying the isotopic effect between H and D.With further increasing collision energy(e.g.,up to 502.08 kJ/mol),O+HD → O+H+D will occur and induce the complete dissociation into single O,H,and D atoms.     

The structure and the potential energy function of AlSO(C_S,X~2A″)

By using the B3P86/aug-cc-pvtz method,the accurate equilibrium geometry of the AlSO(CS,X2A″) molecule has been calculated and compared with available theoretical values.The obtained results show that the AlSO molecule has a most stable structure with bond lengths of R OAl = 0.1864 nm,R OS = 0.1623 nm,R AlS = 0.2450 nm,together with a dissociation energy of 13.88 eV.The possible electronic states and their reasonable dissociation limits for the ground state of the AlSO molecule were determined based on the principle of atomic and molecular reaction statics.The analytic potential energy function of the AlSO molecule was derived by the many-body expansion theory and the contour lines were constructed for the first time,which show the internal information of the AlSO molecule,including the equilibrium structure and stable point.The analysis demonstrates that the obtained potential energy function of AlSO is reasonable and successful and the present investigations provide important insights for further study on molecular reaction dynamics.     

Structure and analytical potential energy function for the ground state of the BCx (x=0, －1)

In this paper, the electronic states of the ground states and dissociation limits of BC and BC- are correctly determined based on group theory and atomic and molecular reaction statics. The equilibrium geometries, harmonic frequencies and dissociation energies of the ground state of BC and BC- are calculated by using density function theory and quadratic CI method including single and double substitutions. The analytical potential energy functions of these states have been fitted with Murrell-Sorbie potential energy function from our ab initio calculation results. The spectroscopic data （αe, ωe and ωeχe） of each state is calculated via the relation between analytical potential energy function and spectroscopic data. All the calculations are in good agreement with the experimental data.     

Simulation of water potential for the electronic structureof serine

First-principles, all-electron, \textit{ab initio} calculations have been performed to construct an equivalent water potential for the electronic structure of serine (Ser) in solution. The calculation is composed of three steps. The first step is to search for the configuration of the Ser + NH₂O system with a minimum energy. The second step is to calculate the electronic structure of Ser with the water molecule potential via the self-consistent cluster-embedding method (SCCE), based on the result obtained in the first step. The last step is to calculate the electronic structure of Ser with the dipole potential after replacing the water molecules with dipoles. The results show that the occupied states of Ser are raised by about 0.017~Ry on average due to the effect of water. The water effect can be successfully simulated by using the dipole potential. The obtained equivalent potential can be applied directly to the electronic structure calculation of protein in solution by using the SCCE method.     

Analytical Potential Energy Function for the Ground State X1∑+ of Lanthanum Monofluoride

The equilibrium geometry, harmonic frequency and bond dissociation energy of lanthanum monofluoride have been calculated using Density-Functional Theory (DFT), post-HF methods MP2 and CCSD(T) with the energyconsistent relativistic effective core potentials. The possible electronic state and reasonable dissociation limit of the ground state of LaF are determined based on atomic and molecular reaction statics. Potential energy curve scans for the ground state X 1∑+ have been performed at B3LYP and CCSD(T) levels, due to their better results of harmonic frequency and bond dissociation energy. We find that the potential energy calculated with CCSD(T) is about 0.6 eV larger than the bond dissociation energy, when the internuclear distance is as large as 0.8 nm. The problem that single-reference ab initio methods do not meet dissociation limit during calculations of lanthanide heavy-metal elements is analyzed. We propose the calculation scheme to derive the analytical Murrell-Sorbie potential energy function. Vibrotational spectroscopic constants Be, ωe, ωeχe, αe, βe, De and He obtained by the standard Dunham treatment coincide well with the results of rotational analyses on spectroscopic experiments.     

Time-Dependent Quantum Dynamics of T＋CH4 → CH3＋HT Reaction

Reaction probability, cross section and rate constant are studied for polyatomic reaction T CH4 → CH3 HT using the semirigid vibrating rotor target (SVRT) model. The numerical calculation for the reaction system is carried out using the time-dependent wavepacket method, and the wavepacket is propagated by the splitoperator method. The calculation exhibits a variety of features that can be used for comparison with future experimental investigations. The reaction probability as a function of the translational energy shows slight oscillatory structures, similar to those observed in H abstraction reactions H H2 and H CH4. The comparisons with the H CH4 reaction are described.     

The ternary Ni–Al–Co embedded-atom-method potential for γ/γ Ni-based single-crystal superalloys: Construction and application

An Ni–Al–Co system embedded-atom-method potential is constructed for the γ(Ni)/γ(Ni3Al) superalloy based on experiments and first-principles calculations. The stacking fault energies(SFEs) of the Ni(Co, Al) random solid solutions are calculated as a function of the concentrations of Co and Al. The calculated SFEs decrease with increasing concentrations of Co and Al, which is consistent with the experimental results. The embedding energy term in the present potential has an important influence on the SFEs of the random solid solutions. The cross-slip processes of a screw dislocation in homogenous Ni(Co) solid solutions are simulated using the present potential and the nudged elastic band method. The cross-slip activation energies increase with increasing Co concentration, which implies that the creep resistance of γ(Ni)may be improved by the addition of Co.     

A typical slow reaction H（2S） ＋ S2（X3∑g） → SH（X2П） ＋S（3P） ona new surface： Quantum dynamics calculations

Quantum dynamics calculations for the title reaction H(2S) + S2(X3-Σg) → SH(X2Π) + S(3P) are performed by using a globally accurate double many-body expansion potential energy surface [J. Phys. Chem. A 115 5274(2011)].The Chebyshev real wave packet propagation method is employed to obtain the dynamical information, such as reaction probability, initial state-specified integral cross section, and thermal rate constant. It is found not only that there is a reaction threshold near 0.7 eV in both reaction probabilities and integral cross section curves, but also that both the probability and cross section increase firstly and then decrease as the collision energy increases. The existence of the resonance structure in both the probability and cross section curves is ascribed to the deep potential well. The calculation of the rate constant reveals that the reaction occurring on the potential energy surface of the ground-state HS2is slow to take place.     

Analytical potential energy function and spectroscopy parameters for B~1Π state of KH

Multi-reference configuration interaction is used to produce potential energy curves(PECs) for the excited B1Π state of KH molecule.To investigate the correlation effect of core-valence electrons,five schemes are employed which include the different correlated electrons and different active spaces.The PECs are fitted into analytical potential energy functions(APEFs).The spectroscopic parameters,ro-vibrational levels,and transition frequencies are determined based on the APEFs and compared with available experimental and theoretical data.The molecular properties for B1Π obtained in this letter,which are better than those available in literature,can be reproduced with calculations using the suitable correlated electrons and active space of orbitals.