高振动激发态吡嗪碰撞传能的QCT计算研究

用准经典轨线(QCT)方法计算了高振动激发态吡嗪(C4N2H4)与N2、O2、NH3、基态吡嗪之间的碰撞传能. C4N2H4通过计算发现, 高振动激发态C4N2H4与N2、O2碰撞发生的主要是V-V传能, 与NH3碰撞发生的主要是V-R传能, 与基态C4N2H4碰撞发生的主要是V-V(R)传能. 通过比较高振动激发态C4N2H4、C6F6、C6H6与其基态分子的碰撞传能， 发现此类碰撞传能中, 若分子的对称性高， 则V-V传能更容易实现.     

Hybrid approach for ab initio molecular dynamics simulation combining energy density analysis and short-time Fourier transform: energy transfer spectrogram

We propose a new analysis technique for specifying molecular vibrational modes related with intramolecular and/or intermolecular energy transfer in ab initio molecular dynamics simulation of chemical reaction. The technique combines the short-time Fourier transform method with energy density analysis, which partitions the quantum chemical potential energy in the system into atomic contributions. The image obtained by the combined scheme, termed an energy transfer spectrogram (ETS), enables us to understand the dynamics of energy transfer by time-frequency representation. The time change of the local energy is quite important in chemical reactions. In order to assess the performance of the ETS, its application to the collision reaction between two carbon dioxide molecules is shown.     

Fragmentation cross sections of protonated water clusters

We have measured fragmentation cross sections of protonated water cluster cations (H(2)O)(n=30-50)H(+) by collision with water molecules. The clusters have well-defined sizes and internal energies. The collision energy has been varied from 0.5 to 300 eV. We also performed the same measurements on deuterated water clusters (D(2)O)(n=5-45)D(+) colliding with deuterated water molecules. The main fragmentation channel is shown to be a sequential thermal evaporation of single molecules following an initial transfer of relative kinetic energy into internal energy of the cluster. Unexpectedly, that initial transfer is very low on average, of the order of 1% of collision energy. We evaluate that for direct collisions (i.e., within the hard sphere radius), the probability for observing no fragmentation at all is more than 35%, independently of cluster size and collision energy, over our range of study. Such an effect is well known at higher energies, where it is attributed to electronic effects, but has been reported only in a theoretical study of the collision of helium atoms with sodium clusters in that energy range, where only vibrational excitation occurs.     

Ab initio QM/MM molecular dynamics study on the excited-state hydrogen transfer of 7-azaindole in water solution

Ab initio molecular dynamics (AIMD) simulations for the excited-state hydrogen transfer (ESHT) reaction of 7-azaindole (7AI-(H2O)n; n = 1, 2) clusters in the gas phase and in water are presented. The effective fragment potential (EFP) is employed to model the surrounding water molecules. The AIMD simulations for 7AI-H2O and 7AI-(H2O)2 clusters show an asynchronous hydrogen transfer at t approximately 50 fs after the photoexcitation. While the ESHT mechanism for 7AI-H2O in water does not change appreciably compared with that in the gas phase, the AIMD simulations on 7AI-(H2O)2 in water solution exhibit two different mechanisms. Since the tautomer form is lower in energy compared to the normal form in the S1 state, 7AI and (H2O) n fragments separate from each other after the ESHT. With the use of the results of the AIMD trajectories, the minimum energy conical intersection point in the tautomer region has also been located.     

Intra- and intermolecular energy transfer in highly excited ozone complexes

The energy transfer of highly excited ozone molecules is investigated by means of classical trajectories. Both intramolecular energy redistribution and the intermolecular energy transfer in collisions with argon atoms are considered. The sign and magnitude of the intramolecular energy flow between the vibrational and the rotational degrees of freedom crucially depend on the projection K(a) of the total angular momentum of ozone on the body-fixed a axis. The intermolecular energy transfer in single collisions between O(3) and Ar is dominated by transfer of the rotational energy. In accordance with previous theoretical predictions, the direct vibrational de-excitation is exceedingly small. Vibration-rotation relaxation in multiple Ar+O(3) collisions is also studied. It is found that the relaxation proceeds in two clearly distinguishable steps: (1) During the time between collisions, the vibrational degrees of freedom are "cooled" by transfer of energy to rotation; even at low pressure equilibration of the internal energy is slow compared to the time between collisions. (2) In collisions, mainly the rotational modes are "cool" by energy transfer to argon.     

Long-range solvent effects on the orbital interaction mechanism of water acidity enhancement in metal ion solutions: a comparative study of the electronic structure of aqueous Mg and Zn dications

We study the dissociation of water coordinated to a divalent metal ion center, M2+ = Mg2+, Zn2+ using density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations. First, the proton affinity of a coordinated OH- group is computed from gas-phase Mg2+(H2O)5(OH-), which yields a relative higher gas-phase acidity for a Zn2+-coordinated as compared to a Mg2+-coordinated water molecule, DeltapKa(gp) = 5.3. We explain this difference on the basis of a gain in stabilization energy of the Zn2+(H2O)5(OH-) system arising from direct orbital interaction between the coordinated OH- and the empty 4s state of the cation. Next, we compute the acidity of coordinated water molecules in solution using free-energy thermodynamic integration with constrained AIMD. This approach yields pKa Mg2+ = 11.2 and pKa Zn2+ = 8.4, which compare favorably to experimental data. Finally, we examine the factors responsible for the apparent decrease in the relative Zn2+-coordinated water acidity in going from the gas-phase (DeltapKa(gp) = 5.3) to the solvated (DeltapKa = 2.8) regime. We propose two simultaneously occurring solvation-induced processes affecting the relative stability of Zn2+(H2O)5(OH-), namely: (a) reduction of the Zn 4s character in solution states near the bottom of the conduction band; (b) hybridization between OH- orbitals and valence-band states of the solvent. Both effects contribute to hindering the OH- --> Zn2+ charge transfer, either by making it energetically unfavorable or by delocalizing the ligand charge density over several water molecules.     

Effect of moisture on the field dependence of mobility for gas-phase ions of organophosphorus compounds at atmospheric pressure with field asymmetric ion mobility spectrometry

The electric field dependence of the mobilities of gas-phase protonated monomers [(MH+(H2O)n] and proton-bound dimers [M2H+(H2O)n] of organophosphorus compounds was determined at E/N values between 0 and 140 Td at ambient pressure in air with moisture between 0.1 and 15 000 ppm. Field dependence was described as alpha (E/N) and was obtained from the measurements of compensation voltage versus field amplitude in a planar high-field asymmetric waveform ion mobility spectrometer. The alpha function for protonated monomers to 140 Td was constant from 0.1 to 10 ppm moisture in air with onset of effect at approximately 50 ppm. The value of alpha increased 2-fold from 100 to 1000 ppm at all E/N values. At moisture values between 1000 and 10 000 ppm, a 2-fold or more increase in alpha (E/N) was observed. In a model proposed here, field dependence for mobility through changes in collision cross sections is governed by the degree of solvation of the protonated molecule by neutral molecules. The process of ion declustering at high E/N values was consistent with the kinetics of ion-neutral collisional periods, and the duty cycle of the waveform applied to the drift tube. Water was the principal neutral above 50 ppm moisture in air, and nitrogen was proposed as the principal neutral below 50 ppm.     

Photodissociation spectroscopy of Zn+(H2O) and Zn+(D2O)

We report on a study of the photodissociation spectroscopy of weakly bound Zn+(H2O) and Zn+(D2O) complexes. The work is supported by ab initio electronic structure calculations of the ground and low-lying excited energy surfaces. We assign two molecular absorption bands in the near UV correlating to Zn+ (4s-4p)-based transitions, and identify vibrational progressions associated with both intermolecular and intramolecular vibrational modes of the cluster. Partially resolved rotational structure is consistent with a C(2V) equilibrium complex geometry. Experimental spectroscopic constants are in very good agreement with ab initio theoretical predictions. Results are compared with previous work on main group and transition metal ion-H2O clusters.     

Long range influence of an excess proton on the architecture of the hydrogen bond network in large-sized water clusters

Infrared spectra of completely size-selected protonated water clusters H+(H2O)n are reported for clusters ranging from n=15 to 100. The behavior of the dangling OH stretch bands shows that the hydrogen bond structure in H+(H2O)n is uniquely different to that of (H2O)n up to the size of n=100, at least. This finding indicates that the presence of an excess proton creates a characteristic morphology in the hydrogen bond network architecture of more than 100 surrounding water molecules.     

HCl hydrates as model systems for protonated water

Ab initio molecular dynamics simulations are presented of vibrational dynamics and spectra of crystal HCl hydrates. Depending on the composition, the hydrates include distinct protonated water forms, which in their equilibrium structures approximate either the Eigen ion H3O+(H2O)3 (in the hexahydrate) or the Zundel H2O...H+...OH2 ion (in the di- and trihydrate). Thus, the hydrates offer the opportunity to study spectra and dynamics of distinct species of protonated water trapped in a semirigid solvating environment. The experimentally measured spectra are reproduced quite well by BLYP/DZVP-level calculations employing Fourier transform of the system dipole. The large overall width (800-1000 cm-1) of structured proton bands reflects a broad range of solvating environments generated by crystal vibrations. The aqueous HCl solution was also examined in search of an objective criterion for separating the contributions of "Zundel-like" and "Eigen-like" protonated forms. It is suggested that no such criterion exists since distributions of proton-related structural properties appear continuous and unimodal. Dipole derivatives with respect to OH and O...H+ stretches in water and protonated water were also investigated to advance the understanding of the corresponding IR intensities. The effects of H bonding and solvation on the intensities were analyzed with the help of the Wannier centers' representation of electron density.     

Molecular dynamics with quantum transitions study of the vibrational relaxation of the HOD bend fundamental in liquid D2O

The molecular dynamics with quantum transitions method is used to study the vibrational relaxation of the HOD bend fundamental in liquid D(2)O. All of the vibrational bending degrees of freedom of the HOD and D(2)O molecules are described by quantum mechanics, while the remaining translational and rotational degrees of freedom are described classically. The effect of the coupling between the rotational and vibrational degrees of freedom of the deuterated water molecules is analyzed. A kinetic mechanism based on three steps is proposed in order to interpret the dynamics of the system. It is shown that intermolecular vibrational energy transfer plays an important role in the relaxation process and also that the transfer of energy into the rotational degrees of freedom is favored over the transfer of energy into the translational motions. The thermalization of the system after the relaxation is reached in a shorter time scale than that of the recovery of the hydrogen bond network. The relaxation and equilibration times obtained compare well with experimental and previous theoretical results.     

Theoretical investigations of proton-bound cluster ions

Several proton-bound cluster ions have been studied by means of coupled cluster calculations with large basis sets. Among these are complexes of a krypton or xenon atom with the cations HCO+, HN2+ and HNCH+. Various spectroscopic properties have been calculated in all cases. Effects of vibrational anharmonicity are particularly pronounced for the intramolecular stretching vibrations of Kr...HN2+ and Xe...HN2+. The proton stretching vibration of (N2)H+(N2) is predicted around 800 cm-1, with a large transition dipole moment of 1.15 D. Both (N2)H+(N2) and (HCN)H+(NCH) have linear centrosymmetric equilibrium structures. Those of (OC)H+(CO) and (HCC-)H+(CCH-) are asymmetric with barrier heights to the centrosymmetric saddle points of 382 and 2323 cm-1, respectively. The dissociation energy of the anionic complex Cl-...HCCH is calculated to be Do = 3665 cm-1, 650 cm-1 larger than the corresponding value for Br-...HCCH. The complex between a fluoride ion and acetylene is more strongly bound and shows strongly anharmonic behaviour, similar to the bihalides FHF- or ClHCl-. Strong Fermi resonance interaction is predicted between nu 3 (approximately proton stretch) and 2 nu 4 (first overtone of intermolecular stretch).     

Short-time Fourier transform analysis of ab initio molecular dynamics simulation: collision reaction between CN and C4H6

Collision reactions between cyano radical (CN) and dimethylacetylene (C4H6) are thought to occur in the atmosphere of Saturn's moon Titan. However, it is difficult to reproduce reactions occurring in unique environments to study their dynamical processes. In this study, collision reactions between CN and C4H6 were investigated using ab initio molecular dynamics (AIMD) simulations. The simulation results were categorized into three kinds: nonreactive collision, incorporation, and substitution. Short-time Fourier transform analysis of velocity autocorrelation functions obtained by the AIMD simulations, which has been recently developed by our research group, was performed to examine the nonequilibrium condition of the vibrational states. Spectrograms, which correspond to the time evolution of power spectra, clarify the relationship between the three reaction channels and the dynamical changes of the vibrational states.     

Reaction dynamics simulations of the identity S(N)2 reaction H(2)O + HOOH(2)(+)--> H(2)OOH(+)+ H(2)O. Requirements for reaction and competition with proton transfer

Despite the fact that the transition structure of the gas phase S(N)2 reaction H(2)O + HOOH(2)(+)--> HOOH(2)(+)+ H(2)O is well below the reactants in potential energy, the reaction has not yet been observed by experiment. Variational transition state RRKM theory reveals a strong preference for the competing proton transfer reaction H(2)O + HOOH(2)(+)--> H(3)O(+)+ HOOH due to entropy factors. Born-Oppenheimer reaction dynamics simulations confirm these results. However, by increasing the collision energy to around 7.5 eV the probability for nucleophilic substitution increases relative to proton transfer. These observations are explained by the presence of the key common intermediate HOO(H)[dot dot dot]H-OH(2)(+) which leads to effective proton transfer, but can be avoided with increasing collision energy. However, the S(N)2 probability remains below 0.2 since successful passage through the TS requires optimum initial orientation of the reactants, excitation of the relative translational motion and good phase correlation between the O-O vibration and the motion of the incoming water.     

Vibrational spectroscopy of protonated imidazole and its complexes with water molecules: ab initio anharmonic calculations and experiments

The results of anharmonic frequency calculations on neutral imidazole (C3N2H4, Im), protonated imidazole (ImH+), and its complexes with water (ImH+)(H2O)n, are presented and compared to gas phase infrared photodissociation spectroscopy (IRPD) data. Anharmonic frequencies are obtained via ab initio vibrational self-consistent field (VSCF) calculations taking into account pairwise interactions between the normal modes. The key results are: (1) Prediction of anharmonic vibrational frequencies on an MP2 ab initio potential energy surface show excellent agreement with experiment and outstanding improvement over the harmonic frequencies. For example, the ab initio calculated anharmonic frequency for (ImH+)(H2O)N2 exhibits an overall average percentage error of 0.6% from experiment. (2) Anharmonic vibrational frequencies calculated on a semiempirical potential energy surface fitted to ab initio harmonic data represents spectroscopy well, particularly for water complexes. As an example, anharmonic frequencies for (ImH+)H2O and (ImH+)(H2O)2 show an overall average deviation of 1.02% and 1.05% from experiment, respectively. This agreement between theory and experiment also supports the validity and use of the pairwise approximation used in the calculations. (3) Anharmonic coupling due to hydration effects is found to significantly reduce the vibrational frequencies for the NH stretch modes. The frequency of the NH stretch is observed to increase with the removal of a water molecule or replacement of water with N2. This result also indicates the ability of the VSCF method to predict accurate frequencies in a matrix environment. The calculation provides insights into the nature of anharmonic effects in the potential surface. Analysis of percentage anharmoncity in neutral Im and ImH+ shows a higher percentage anharmonicity in the NH and CH stretch modes of neutral Im. Also, we observe that anharmonicity in the NH stretch modes of ImH+ have some contribution from coupling effects, while that of neutral Im has no contribution whatsoever from mode-mode coupling. It is concluded that the incorporation of anharmonic effects in the calculation brings theory and experiment into much closer agreement for these systems.     

Time-domain calculations of the polarized Raman spectra, the transient infrared absorption anisotropy, and the extent of delocalization of the OH stretching mode of liquid water

The polarized Raman spectrum and the time dependence of the transient infrared (TRIR) absorption anisotropy are calculated for the OH stretching mode of liquid water (neat liquid H2O) by using time-domain formulations, which include the effects of both the diagonal frequency modulations (of individual oscillators) induced by the interactions between the dipole derivatives and the intermolecular electric field, and the off-diagonal (intermolecular) vibrational coupling described by the transition dipole coupling (TDC) mechanism. The IR spectrum of neat liquid H2O and the TRIR anisotropy of a liquid mixture of H2O/HDO/D2O are also calculated. It is shown that the calculated features of these optical signals, including the temperature dependence of the polarized Raman and IR spectra, are in reasonable agreement with the experimental results, indicating that the frequency separation between the isotropic and anisotropic components of the polarized Raman spectrum and the rapid decay (approximately 0.1 ps) of the TRIR anisotropy of the OH stretching mode of neat liquid H2O are mainly controlled by the resonant intermolecular vibrational coupling described by the TDC mechanism. Comparing with the time evolution of vibrational excitations, it is suggested that the TRIR anisotropy decays in the time needed for the initially localized vibrational excitations to delocalize over a few oscillators. It is also shown that the enhancement of the dipole derivatives by the interactions with surrounding molecules is an important factor in generating the spectral profiles of the OH stretching Raman band. The time-domain behavior of the molecular motions that affect the spectroscopic features is discussed.     

Solvent dependence of OH bend vibrational relaxation of monomeric water molecules in liquids

The vibrational relaxation rates of the OH bending mode of monomeric H(2)O molecules diluted in various liquid halogenated methane and ethane derivates have been determined by a picosecond infrared pump-probe study. Relaxation time constants between 4.8 and 40.5 ps have been obtained. The discussion of the general solvent dependence suggests that in all cases the solvent fundamental with the smallest energy mismatch is favorably populated by this intermolecular energy transfer process.     

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.     

Ab initio potential energy surfaces and nonadiabatic collision dynamics in H(+)+O(2) system

The adiabatic potential energy surfaces for the lowest five electronic states of (3)A" symmetry for the H(+)+O(2) collision system have been obtained at the multireference configuration interaction level of accuracy using Dunning's correlation consistent polarized valence triple zeta basis set. The radial nonadiabatic coupling terms and the mixing angle between the lowest two electronic states (1 (3)A" and 2 (3)A"), which adiabatically correlate in the asymptotic limit to H((2)S)+O(2) (+)(X (2)Pi(g)) and H(+)+O(2)(X (3)Sigma(g)(-)), respectively, have been computed using ab initio procedures at the same level of accuracy to yield the corresponding quasidiabatic potential energy matrix. The computed strengths of the vibrational coupling matrix elements reflect the trend observed for inelastic vibrational excitations of O(2) in the experiments at collision energy of 9.5 eV. The quantum dynamics has been preformed on the newly obtained coupled quasidiabatic potential energy surfaces under the vibrational close-coupling rotational infinite-order sudden framework at the experimental collision energy of 9.5 eV. The present theoretical results for vibrational elastic/inelastic excitations of O(2) are in overall good agreement with the available experimental data obtained from the proton energy-loss spectra in molecular beam experiments [F. A. Gianturco et al., J. Phys. B 14, 667 (1981)]. The results for the complementary charge transfer processes are also presented at this collision energy.