Second-order Brillouin-Wigner perturbation theory: size-extensivity correction

A convergence study of intensities of transitions from vibrational ground state to the lower lying states is done with respect to the rank of the dipole moment surface (DMS) in the Taylor series expansion of the DMS. The relative roles of the mechanical and electrical anharmonicity are analyzed in the calculation of the intensities of vibrational transitions from ground state. We find that at least a quadratic expansion of the dipole moment functional is necessary to predict the intensities of vibrational transitions. The mechanical anharmonicity becomes important when the resonances between the vibrational states are significant.     

Abnormal inversion splitting in NH2D: rotational analysis of the nu5 bending vibrational band system

A high-resolution Fourier transform infrared spectrum of the nu(5) bending vibrational band system region of the partially deuterated ammonia molecule NH(2)D has been measured and rotationally analyzed. The spectrum consists of strong a-type transitions between the states of same vibrational symmetry and weaker c-type transitions between the states of different vibrational symmetry. The Hamiltonian model used includes interaction terms between the rotational states of both upper and lower inversion doublets. The vibrational term values for the symmetric and the antisymmetric component of the upper-inversion doublet are 1,605.637 965(620) cm(-1) and 1,590.993 82(100) cm(-1), respectively, where the numbers in parentheses are one-standard deviations in the least significant digit. These figures are close to the corresponding values 1,605.62 cm(-1) and 1,590.72 cm(-1) obtained recently from results based on high-level ab initio calculations. The order of the vibrational term values is abnormal in the ammonia family, as typically the symmetric state is lower in wavenumber than the antisymmetric one.     

Molecular vibrations of methane molecules in the structure I clathrate hydrate from ab initio molecular dynamics simulation

Vibrational frequencies of guest molecules in clathrate hydrates reflect the molecular environment and dynamical behavior of molecules. A detailed understanding of the mechanism for the vibrational frequency changes of the guest molecules in the clathrate hydrate cages is still incomplete. In this study, molecular vibrations of methane molecules in a structure I clathrate hydrate are calculated from ab initio molecular dynamics simulation. The vibrational spectra of methane are computed by Fourier transform of autocorrelation functions, which reveal distinct separation of each vibrational mode. Calculated symmetric and asymmetric stretching vibrational frequencies of methane molecules are lower in the large cages than in the small cages (8 and 16 cm(-1) for symmetric and asymmetric stretching, respectively). These changes are closely linked with the C-H bond length. The vibrational frequencies for the bending and rocking vibrational modes nearly overlap in each of the cages.     

Optical control of excited-state vibrational coherences of a molecule in solution: The influence of the excitation pulse spectrum and phase in LD690

Spectral and phase shaping of femtosecond laser pulses is used to selectively excite vibrational wave packets on the ground (S0) and excited (S1) electronic states in the laser dye LD690. The transient absorption signals observed following excitation near the peak of the ground-state absorption spectrum are characterized by a dominant 586 cm(-1) vibrational mode. This vibration is assigned to a wave packet on the S0 potential energy surface. When the excitation pulse is tuned to the blue wing of the absorption spectrum, a lower frequency 568 cm(-1) vibration dominates the response. This lower frequency mode is assigned to a vibrational wave packet on the S1 electronic state. The spectrum and phase of the excitation pulse also influence both the dephasing of the vibrational wave packet and the amplitude profiles of the oscillations as a function of probe wavelength. Excitation by blue-tuned, positively chirped pulses slows the apparent dephasing of the vibrational coherences compared with a transform-limited pulse having the same spectrum. Blue-tuned negatively chirped excitation pulses suppress the observation of coherent oscillations in the ground state.     

Internal Energy State Distribution of Product CaBr in the Collision Reactions Ca+C2H5Br和Ca+nC3H7Br

The internal energy distributions of product CaBr in the collision reactions Ca+C2H5Br and Ca+nC3H7Br are studied by using the quasiclassical trajectory method. The average vibrational, rotational and translational energies and total available energies of the product CaBr molecules are calculated. The results indicate that when the collision energy is equal to 7.54 kJ/mol the energy of product CaBr is mainly the vibrational energy. As the reactant collision energy increases, the average translational and rotational energies of the product CaBr increase, the average vibrational energy decreases slightly, and the most probable vibrational state shifts to lower vibrational energy levels. The internal states of reagents have little influence on the internal energy distribution of the product. The bigger the radical group is, the higher ratio of the vibrational energy to the available energy of the product is. There exist two competitive reaction paths for the collision reactions Ca+C2H5Br and Ca+nC3H7Br, the migratory encounter and direct reaction paths. The former produces high vibrational excited state product CaBr and the latter causes C-Br bond to break. When the collision energy increases, the reactions tend to the latter path.     

Vibrational spectra and ab initio analysis of tert-butyl, trimethylsilyl, trimethylgermyl, trimethylstannyl, and trimethylplumbyl derivatives of 3,3-dimethylcyclopropene X. Some regularities

The changes in the vibrational frequencies of 1-tert-butyl and 1,2-di-tert-butyl derivatives of 3,3-dimethylcyclopropene brought about by substitution of the central carbon atom (X) of the tert-butyl moieties by Si, Ge, Sn, or Pb atoms are examined. The most important decrease in the vibrational frequencies implicating the X(CH(3))(3) moieties is noted for substitution of X=C by Si. The substitutions of Si by Ge or Ge by Sn or Sn by Pb are not accompanied by the pronounced frequency shifts observed for the C-->Si transition. An explanation is given for trends in these vibrational frequencies for the transitions X=C-->Si-->Ge-->Sn-->Pb. It is concluded that there are lower limiting values of the vibrational frequencies of a molecular moiety which are approached when the mass of its isovalent atom is increased. This leads to the formation of cluster regions in the vibrational spectra for the frequencies of the SnC(3) and PbC(3) moieties.     

Intramolecular vibrational coupling contribution to temperature dependence of vibrational mode frequencies

High-frequency vibrational modes in molecules in solution are sensitive to temperature and shift either to lower or higher frequencies with the temperature increase. These frequency shifts are often attributed to specific interactions of the molecule and to the solvent polarization effect. We found that a substantial and often dominant contribution to sensitivity of vibrational high-frequency modes to temperature originates from anharmonic interactions with other modes in the molecule. The temperature dependencies were measured for several modes in ortho-, meta-, and para-isomers of acetylbenzonitrile in solution and in a solid matrix and compared to the theoretical predictions originated from the intramolecular vibrational coupling (IVC) evaluated using anharmonic density functional theory calculations. It is found that the IVC contribution is essential for temperature dependencies of all high-frequency vibrational modes and is dominant for many modes. As such, the IVC contribution alone permits predicting the main trend in the temperature dependencies, especially for vibrational modes with smaller transition dipoles. In addition, an Onsager reaction field theory was used to describe the solvent contribution to the temperature dependencies.     

Langevin model of the temperature and hydration dependence of protein vibrational dynamics

The modification of internal vibrational modes in a protein due to intraprotein anharmonicity and solvation effects is determined by performing molecular dynamics (MD) simulations of myoglobin, analyzing them using a Langevin model of the vibrational dynamics and comparing the Langevin results to a harmonic, normal mode model of the protein in vacuum. The diagonal and off-diagonal Langevin friction matrix elements, which model the roughness of the vibrational potential energy surfaces, are determined together with the vibrational potentials of mean force from the MD trajectories at 120 K and 300 K in vacuum and in solution. The frictional properties are found to be describable using simple phenomenological functions of the mode frequency, the accessible surface area, and the intraprotein interaction (the displacement vector overlap of any given mode with the other modes in the protein). The frictional damping of a vibrational mode in vacuum is found to be directly proportional to the intraprotein interaction of the mode, whereas in solution, the friction is proportional to the accessible surface area of the mode. In vacuum, the MD frequencies are lower than those of the normal modes, indicating intramolecular anharmonic broadening of the associated potential energy surfaces. Solvation has the opposite effect, increasing the large-amplitude vibrational frequencies relative to in vacuum and thus vibrationally confining the protein atoms. Frictional damping of the low-frequency modes is highly frequency dependent. In contrast to the damping effect of the solvent, the vibrational frequency increase due to solvation is relatively temperature independent, indicating that it is primarily a structural effect. The MD-derived vibrational dynamic structure factor and density of states are well reproduced by a model in which the Langevin friction and potential of mean force parameters are applied to the harmonic normal modes.     

RDX分子内振动能量重分配的动力学研究

基于从头算分子动力学(Born-oppenheimer molecular dynamics, BOMD)模拟, 构建了环硝胺六氢-1,3,5-三硝基-1,3,5-三嗪(RDX)单分子不同振动模式之间的耦合矩阵, 并计算了在不同加载能量下从低频振动模式到高频振动模式的最优能量传输路径. 结果表明, RDX单分子中—NNO₂基团更有利于能量局域化, 振动模式v₃和v₄在从低频振动模式到高频振动模式的能量传输过程中扮演着重要角色. 通过对v₃和v₄两个振动模式的进一步分析发现, 加载能量的不同会导致RDX单分子能量传输路径的不同. 当加载能量较低时, RDX单分子倾向于从低频振动模式到中频振动模式再到高频振动模式的能量传输路径; 当加载能量较高时, 能量更倾向于从低频振动模式直接传输到高频振动模式上. 揭示了RDX分子内振动耦合能量转移的微观机制, 为进一步探索RDX将“机械能”转化为“化学能”的微观过程提供了理论基础.     

Vibrational frequencies and structural determinations of hexamethylenetetraamine

The normal mode frequencies and corresponding vibrational assignments of hexamethylenetetraamine (HMTA) in Td symmetry are examined theoretically using the Gaussian 98 set of quantum chemistry codes. All normal modes were successfully assigned to one of eight types of motion predicted by a group theoretical analysis. The vibrational modes of the deuterated form of HMTA (HMTA d-12) were also calculated and compared against experimental data. The normal mode vibrational frequencies were shifted to lower frequencies as on deuteration as expected. However, in some cases the dominant motion type changed on deuteration leading to an apparent 'blue shift' of some of the N-C stretching modes. It is possible that the observed frequency shifts are the result of a Fermi resonance condition.     

Theoretical study on the vibrational structure of S0 thiophosgene from the origin to dissociation

In this work, using our vibrational variational calculation method and a recently derived refined quartic potential energy surface for S₀ thiophosgene, we have carried out large scale vibrational calculations to analyze the vibrational structure of this electronic state in the whole range of vibrational excitation energies down from the origin and up to the dissociation limit (at ↼20,000 cm^↙1). In the lower excited vibrational range we have achieved satisfactory coincidence of calculated to experimentally measured frequencies, while at the higher vibrational excitations our main objective has been to estimate what part of the available vibrational level density is effectively involved into the vibrational mixing and IVR. The results from our calculations have been compared to the available experimentally obtained dataset (obtained from synchrotron IR, SEP and LIF spectra) as well as to conclusions from the analyses by other authors using local coupling models.     

Vibrational spectroscopy of nitroalkane chains using electron autodetachment and ar predissociation

If the binding energy of an excess electron is lower than some of the vibrational levels of its host anion, vibrational excitation can lead to autodetachment. We use excitation of CH stretching modes in nitroalkane anions (2700-3000 cm(-1)), where the excess electron is localized predominantly on the NO2 group. We present data on nitroalkane anions of various chain lengths, showing that this technique is a valid approach to the vibrational spectroscopy of such systems extending to nitroalkane anions at least the size of nitropentane. We compare spectra taken by using vibrational autodetachment with spectra obtained by monitoring Ar evaporation from Ar solvated nitroalkane anions. The spectra of nitromethane and nitroethane are assigned on the basis of ab initio calculations with a detailed analysis of Fermi resonances of CH stretching fundamentals with overtones and combination bands of HCH bending modes.     

Vibrational calculations in formaldehyde: the CH stretch system

An alternative procedure for the calculation of highly excited vibrational levels in S0 formaldehyde was developed to apply to larger molecules. It is based on a new set of symmetrized vibrational valence coordinates. The fully symmetrized vibrational kinetic energy operator is derived in these coordinates using the Handy expression [Molec. Phys. 61, 207 (1987)]. The potential energy surface is expressed as a fully symmetrized quartic expansion in the coordinates. We have performed ab initio electronic computations using GAMESS to obtain all force constants of the S0 formaldehyde quartic force field. Our large scale vibrational calculations are based on a fully symmetrized vibrational basis set, in product form. The vibrational levels are calculated one by one using an artificial intelligence search/selection procedure and subsequent Lanczos iteration, providing access to extremely high vibrational energies. In this work special attention has been given to the CH stretch system by calculating the energies up to the fifth CH stretch overtone at ∼16000 cm⁻¹, but the method has also been tested on two highly excited combination levels including other lower frequency modes.      

Spin�Corbit ab initio study of two low-lying states of chloroiodomethane cation

Spin?Corbit coupling plays a crucial role in the determination of molecular structure and calculation of vibrational frequencies of CH₂ClI⁺. We performed the geometry optimizations and vibrational frequency calculations of both the lower and upper spin?Corbit (SO) states using an ab initio SO method based on multiconfigurational wave function. The multistate complete active space perturbation theory second-order SO (MS-CASPT2-SO) method reasonably describes the structures of the lower SO state, yielding the C?CI distance and the I?CC?CCl angle close to the experimental values. The geometrical parameters of the upper SO state is quite similar to that of the lower SO state, whereas structures of two states differ substantially in calculations prior to the introduction of SO coupling. The MS-CASPT2-SO method reproduces the difference between the lower and upper SO states for the I?CC?CCl bending frequency. The vibrational frequencies calculated by MS-CASPT2-SO generally overestimate in comparison with the experiments. The energy gap between the two SO states calculated by MS-CASPT2-SO is reasonably close to the experimental value. To the best of our knowledge, this is the first attempt to calculate vibrational frequencies of two SO states of CH₂ClI⁺, and the first time to apply MS-CASPT2-SO method to the geometry optimization and vibrational frequency calculation of polyatomic molecule.     

Non-adiabatic effects on oxygen atom fine structure populations in the predissociation of the A2Σ+ state of OH

Multichannel scattering calculations are used to evaluate the role of non-adiabatic interactions and interference effect on the predissociation linewidths of several vibrational levels of the A²Σ⁺ state of OH due to coupling to three different repulsive electronic states. Oxygen atom fine structure branching ratios are calculated as a possible probe of the predissociation pathway. Non-adiabatic effects are found to be very strong for the lower vibrational predissociative A²Σ⁺ levels and to remain non-negligible even for the higher vibrational states of A ²Σ⁺.     

Three-dimensional quantum dynamics study of vibrational predissociation of HeI_2 van der Waals molecule for low vibrational excitation using the time-dependent wave packet method

Three-dimensional quantum mechanical calculations for vibrational predissociation of HeI2(B) van der Waals molecules are presented using the time-dependent wave packet technique within the golden rule approxima tion.The total and partial decay widths,lifetimes,rates and their dependence on initial vibrational states were obtained for HeI2 at low initial vibrational excited levels.Our calculations show that the calculated tota decay widths,lifetimes and rates agree well with those extrapolated from experimental data available The predicted total decay widths as a function of initial vibrational states exhibit highly nonlinear behavior.The very short propagation time (less.than 1 ps) required in the golden rule wave packet calculation is determined by the duration time of the final state inter-action between the fragments on the vibrationally deexcited adiabatic potential surface.The final state interaction between the fragments is shown to play an important role in determining the final rotational distri     

Time-dependent treatment of intramolecular energy transfer for HeICl complex

The time-dependent golden wave packet method has been used for calculating the decay widths of vibrational predissociation for HeICl complex in the B state with total angular momentum J=0. This is a good example of intramolecular energy transfer. We examine the dependence of the final rotational distribution (partial decay width) of ICl fragment on the stretching excitation. It is found that computed final rotational distributions are weakly dependent on the vibrational level being excited. Unlike the smoothly varying rotational distribution for lower initial vibrational levels, for higher initial vibrational levels the rotational distribution indicates the very pronounced oscillatory structure. The analysis of the rotational distribution as a function of propagation time reveals the predominant role of the final states interaction in determining the final rotational distribution.     

Electronic effects on photochemistry: the diverse reaction dynamics of highly excited stilbenes and azobenzene

Ultrafast time-resolved mass spectrometry and structural dynamics experiments on trans-stilbene, cis-stilbene, and azobenzene, with excitation to high-lying electronic states, reveal a rich diversity of photochemical reaction dynamics. All processes are found to be quite unlike the well-known photochemistry on lower electronic surfaces. While in trans-stilbene, excitation at 6 eV induces a phenyl twisting motion, in cis-stilbene it leads to an ultrafast ring-closing to form 4a,4b-dihydrophenanthrene. Azobenzene dissociates on an ultrafast time scale, rather than isomerizing as it does on a lower surface. The photochemical dynamics of the sample molecules proceed along steep potential energy surfaces and conical intersections. Because of that, the dynamics are much faster than vibrational relaxation, the randomizing effects from vibrational energy scrambling are avoided, and excitation-energy specific reaction dynamics results.     

一维精确量子散射研究——H+ClH(u)在反应势能面上的振动去激

HCl化学激光中存在振动激发的HCl及游离的H和Cl,故HCl在H原子和Cl原子碰撞下振动弛豫速率过程的研究很重要。不久前我们报道了Cl原子对HCl碰撞去激的一维精确量子散射研究,本文用类似方法,讨论H原子对激发态的HCl的碰撞去激。     

Three-dimensional quantum dynamics study of vibrational predissociation of HeI2 van der Waals molecule for low vibrational excitation using the time-dependent wave packet method

Three-dimensional quantum mechanical calculations for vibrational predissociation of He1₂(B) van der Waals molecules are presented using the time-dependent wave packet technique within the golden rule approximation. The total and partial decay widths, lifetimes, rates and their dependence on initial vibrational states were obtained for HeI₂ at low initial vibrational excited levels. Our calculations show that the calculated total decay widths, lifetimes and rates agree well with those extrapolated from experimental data available. The predicted total decay widths as a function of initial vibrational states exhibit highly nonlinear behavior. The very short propagation time (less than 1 ps) required in the golden rule wave packet calculation is determined by the duration time of the final state interaction between the fragments on the vibrationally deexcited adiabatic potential surface. The final state interaction between the fragments is shown to play an important role in determining the final rotational distribution. This interpretation clearly explains the dynamical effect that the final rotational distribution shifts to the lower rotational energy levels as the initial vibrational quantum numberu increases.