Instantaneous pair theory for high-frequency vibrational energy relaxation in fluids

Vibrational energy relaxation of a diatomic molecule in a superfluid helium nanodroplet: influence of the nanodroplet size, interaction energy and energy gap

The influence of the nanodroplet size, molecule-helium interaction potential energy and ν = 1 - ν = 0 vibrational energy gap on the vibrational energy relaxation (VER) of a diatomic molecule (X₂) in a superfluid helium nanodroplet [HeND or (⁴He)_N; finite quantum solvent at T = 0.37 K] has been studied using a hybrid quantum approach recently proposed by us and taking as a reference the VER results on the I₂@(⁴He)₁₀₀ doped nanodroplet (Vilà et al., Phys. Chem. Chem. Phys., 2018, 20, 118, which corresponds to the first theoretical study on the VER of molecules embedded in a HeND). This has allowed us to obtain a deeper insight into the vibrational relaxation dynamics. The nanodroplet size has a very small effect on the VER, as this process mainly depends on the interaction between the molecule and the nanodroplet first solvation shell. Regarding the interaction potential energy and the energy gap, both factors play an important and comparable role in the VER time properties (global relaxation time, lifetime and transition time). As the former becomes stronger the relaxation time properties decrease in a significant way (their inverse follows a linear dependence with respect to the ν = 1 - ν = 0 coupling term) and they also decrease in a significant manner when the energy gap diminishes (linear dependence on the ν = 1 - ν = 0 energy difference). We expect that this study will motivate further work on the vibrational relaxation process in HeNDs.

Vibrational energy relaxation dynamics of diatomic molecules inside superfluid helium nanodroplets. The case of the I2 molecule

The vibrational relaxation (VER) of a X₂ molecule in a ⁴He superfluid nanodroplet (HeND; 0.37 K) was studied adapting a quantum approach recently proposed by us. In the first theoretical study on the VER of molecules inside HeND the I₂ molecule was examined [cascade mechanism (ν → ν − 1; ν − 1 → ν − 2; …) and time scale of ns].

Applicability of the Caldeira-Leggett Model to Vibrational Spectroscopy in Solution

Formulating a rigorous system-bath partitioning approach remains an open issue. In this context, the famous Caldeira-Leggett model that enables quantum and classical treatment of Brownian motion on equal footing has enjoyed popularity. Although this model is by any means a useful theoretical tool, its ability to describe anharmonic dynamics of real systems is often taken for granted. In this Letter, we show that the mapping between a molecular system under study and the model cannot be established in a self-consistent way, unless the system part of the potential is taken effectively harmonic. Mathematically, this implies that the mapping is not invertible. This "invertibility problem" is not dependent on the peculiarities of particular molecular systems and is rooted in the anharmonicity of the system part of the Caldeira-Leggett model potential.

Vibrational mode assignment of finite temperature infrared spectra using the AMOEBA polarizable force field

The Driven Molecular Dynamics approach has been adapted and associated with the AMOEBA polarizable force field to assign and visualize vibrational modes in infrared spectra obtained by molecular dynamics simulations.

Vibrational energy relaxation of large-amplitude vibrations in liquids

Given the limited intermolecular spaces available in dense liquids, the large amplitudes of highly excited, low frequency vibrational modes pose an interesting dilemma for large molecules in solution. We carry out molecular dynamics calculations of the lowest frequency ("warping") mode of perylene dissolved in liquid argon, and demonstrate that vibrational excitation of this mode should cause identifiable changes in local solvation shell structure. But while the same kinds of solvent structural rearrangements can cause the non-equilibrium relaxation dynamics of highly excited diatomic rotors in liquids to differ substantially from equilibrium dynamics, our simulations also indicate that the non-equilibrium vibrational energy relaxation of large-amplitude vibrational overtones in liquids should show no such deviations from linear response. This observation seems to be a generic feature of large-moment-arm vibrational degrees of freedom and is therefore probably not specific to our choice of model system: The lowest frequency (largest amplitude) cases probably dissipate energy too quickly and the higher frequency (more slowly relaxing) cases most likely have solvent displacements too small to generate significant nonlinearities in simple nonpolar solvents. Vibrational kinetic energy relaxation, in particular, seems to be especially and surprisingly linear.

Dynamics of relaxation and fragmentation in size-selected icosahedral Ar(n)[NO(-)(v = 1)] clusters

We study the vibrational relaxation and solvation dynamics in size-selected icosahedral Ar(n)(NO(-)) at 300 K, where NO(-)(X(3)Σ(-)) is in v = 1 and n = 1-12, using a classical dynamics method and an interaction model consisting of detailed host-guest and host-host interactions. Two relaxation time scales are found: (i) the short-time (<200 ps), in which rate is nearly independent of cluster size, and (ii) the ns scale, in which a slow energy transfer process occurs between NO(-) vibration and argon modes at a rate (~10(8) s(-1)) decreasing slightly from n = 12 to 6 and rapidly from n = 5 to 1 (~10(6) s(-1)). In Ar(12)(NO(-)), less than one-quarter of the host atoms sampled evaporate, nearly 60% of evaporation occurring within 200 ps caused by rapid energy transfer from NO(-) at short time. The fraction of evaporation decreases nearly exponentially with increasing evaporation time, but ~16% of evaporation still occurs on a time scale longer than 1 ns. Evaporation from one hemisphere of Ar(12)(NO(-)) dominates the rest. Final cluster sizes commonly produced from the fragmentation of Ar(12)(NO(-)) are n = 6-11 (evaporation of 6-1 atoms) and n = 12 (no evaporation).

Comparison of two simple models for high frequency friction: exponential versus Gaussian wings

Liquid phase vibrational energy relaxation (VER) times T1 typically depend critically on the relaxing mode's high frequency friction or wing function. The wing function may, in principle, be found from the mode's normalized force autocorrelation function (faf) C(t), since it is proportional to lim formula see text C(t)dt. However, the full form of C(t) is never available. Thus, the wing function is typically estimated from a model faf CM(t) which duplicates the known part of C(t) and which (hopefully) approximates its unknown part with enough realism to yield meaningful formula see text behavior. Unfortunately, apparently realistic CM(t)'s can predict unphysical wing functions, and T1's in error by tens of orders of magnitude. Thus, a condition is needed to discriminate between CM(t)'s which yield meaningful and unphysical forms for the high frequency friction. This condition is shown to be that model faf's CM(t) yield physical wing functions if and only if these functions derive from the short time "heads" of the faf's. This test is applied to the model faf's Cga(t) triple bond exp[-1/2(t/tau)2] and Cse(t) triple bond sech(t/tau). These faf's cannot both be physical, since they yield incompatible Gaussian and exponential wing functions. The test accepts Cga(t) as physical. It, however, rejects Cse(t), since its "tail" lim formula see text Cse(t) = 2 exp(-t/tau) (because of its long range) dominates the wing function.

Extracting effective normal modes from equilibrium dynamics at finite temperature

A general method for obtaining effective normal modes of a molecular system from molecular dynamics simulations is presented. The method is based on a localization criterion for the Fourier transformed velocity time-correlation functions of the effective modes. For a given choice of the localization function used, the method becomes equivalent to the principal mode analysis (PMA) based on covariance matrix diagonalization. On the other hand, a proper choice of the localization function leads to a novel method with a strong analogy with the usual normal mode analysis of equilibrium structures, where the Hessian system at the minimum energy structure is replaced by the thermal averaged Hessian, although the Hessian itself is never actually calculated. This method does not introduce any extra numerical cost during the simulation and bears the same simplicity as PMA itself. It can thus be readily applied to ab initio molecular dynamics simulations. Three such examples are provided here. First we recover effective normal modes of an isolated formaldehyde molecule computed at 20 K in very good agreement with the results of a normal mode analysis performed at its equilibrium structure. We then illustrate the applicability of the method for liquid phase studies. The effective normal modes of a water molecule in liquid water and of a uracil molecule in aqueous solution can be extracted from ab initio molecular dynamics simulations of these two systems at 300 K.

The molecular origins of nonlinear response in solute energy relaxation: The example of high-energy rotational relaxation

A key step in solution-phase chemical reactions is often the removal of excess internal energy from the product. Yet, the way one typically studies this process is to follow the relaxation of a solute that has been excited into some distribution of excited states quite different from that produced by any reaction of interest. That the effects of these different excitations can frequently be ignored is a consequence of the near universality of linear-response behavior, the idea that relaxation dynamics is determined by the solvent fluctuations (which may not be all that different for different kinds of solute excitation). Nonetheless, there are some clear examples of linear-response breakdowns seen in solute relaxation, including a recent theoretical and experimental study of rapidly rotating diatomics in liquids. In this paper we use this rotational relaxation example to carry out a theoretical exploration of the conditions that lead to linear-response failure. Some features common to all of the linear-response breakdowns studied to date, including our example, are that the initial solute preparation is far from equilibrium, that the subsequent relaxation promotes a significant rearrangement of the liquid structure, and that the nonequilibrium response is nonstationary. However, we show that none of these phenomena is enough to guarantee a nonlinear response. One also needs a sufficient separation between the solute time scale and that of the solvent geometry evolution. We illustrate these points by demonstrating precisely how our relaxation rate is tied to our liquid-structural evolution, how we can quantitatively account for the initial nonstationarity of our effective rotational friction, and how one can tune our rotational relaxation into and out of linear response.

Host-assisted intramolecular vibrational relaxation at low temperatures: OH in an argon cage

The vibrational relaxation of hydroxyl radicals in the A (2)Sigma(+) (v=1) state has been studied using the semiclassical perturbation treatment at cryogenic temperatures. The radical is considered to be trapped in a closest packed cage composed of the 12 nearest argon atoms and undergoes local translation and hindered rotation around the cage center. The primary relaxation pathway is towards local translation, followed by energy transfer to rotation through hindered-to-free rotational transitions. Free-to-free rotational transitions are found to be unimportant. All pathways are accompanied by the propagation of energy to argon phonon modes. The deexcitation probability of OH(v=1) is 1.3 x 10(-7) and the rate constant is 4.7 x 10(5) s(-1) between 4 and 10 K. The negligible temperature dependence is attributed to the presence of intermolecular attraction (>>kT) in the guest-host encounter, which counteracts the T(2) dependence resulting from local translation. Calculated relaxation time scales are much shorter than those of homonuclear molecules, suggesting the importance of the hindered and free motions of OH and strong guest-host interactions.

Dissipative dynamics of laser induced nonadiabatic molecular alignment

Nonadiabatic alignment induced by short, moderately intense laser pulses in molecules coupled to dissipative environments is studied within a nonperturbative density matrix theory. We focus primarily on exploring and extending a recently proposed approach [Phys. Rev. Lett. 95, 113001 (2005)], wherein nonadiabatic laser alignment is used as a coherence spectroscopy that probes the dissipative properties of the solvent. To that end we apply the method to several molecular collision systems that exhibit sufficiently varied behavior to represent a broad variety of chemical environments. These include molecules in low temperature gas jets, in room temperature gas cells, and in dense liquids. We examine also the possibility of prolonging the duration of the field free (post-pulse) alignment in dissipative media by a proper choice of the system parameters.

Conformational Properties of and a Reorientation Triggered by Sugar−Water Vibrational Resonance in the Hydroxymethyl Group in Hydrated β-Glucopyranose

In this paper, we discuss the conformational properties of the hydroxymethyl group of beta-glucopyranose in aqueous solution and its reorientation mechanism. First, using the values for the hydroxymethyl torsion (O5-C5-C6-O6) angle obtained by our ab initio simulations, we reestimate the experimental ratio of the hydroxymethyl rotamer populations. The reestimated ratio is found to be in agreement with those previously reported in several computational studies, which probably partly explains the discrepancies between theoretical and experimental studies that have been discussed in the literature. Second, our time-frequency analysis on a reorientation in the hydroxymethyl group in an ab initio molecular dynamics trajectory suggests that, before the reorientation, the O6-H6 stretching mode is vibrationally coupled with a proton-accepting first-hydration-shell water molecule, whereas the C6-O6 stretching mode is vibrationally coupled with a proton-donating one. The amount of the total vibrational energy induced by these vibrational couplings is estimated to be comparable to typical values for the potential barriers between hydroxymethyl rotamers. To elucidate the vibrational couplings, we investigate the hydrogen-bonding properties around the hydroxymethyl group during the pretransition period. The implications, validity, and limitation of a possible reorientation mechanism based on these findings are also discussed.

Mechanism for singular behavior in vibrational spectra of topologically disordered systems: short-range attractions

At low-enough fluid densities, we have found some naive singular behavior, like the van Hove singularities in the phonon spectra of lattices, appearing in the instantaneous normal mode spectra of the Lennard-Jones (LJ) 2n-n fluids, which serve as a prototype of topologically disordered systems. The singular behavior cannot be predicted by the mean-field theory, but interpreted by the perturbed binary modes of some special pairs, called the mutual nearest neighbor pairs, at separations corresponding to the extreme binary frequencies, which are solely determined by the attractive part of the LJ 2n-n pair potential. By reducing the range of attraction in the pair potential under the conditions of the same particle diameter and well depth, the tendency for the appearance of the singular behavior shifts to higher fluid densities. From this study, we conclude that pair potential with a short-range attraction can be a mechanism to produce a counterpart of the van Hove singularity in the vibrational spectra of disordered systems without a reference lattice.

Surface hopping simulation of the vibrational relaxation of I2 in liquid xenon using the collective probabilities algorithm

A surface hopping simulation of the vibrational relaxation of highly excited I(2) in liquid xenon is presented. The simulation is performed by using the collective probabilities algorithm which assures the coincidence of the classical and quantum populations. The agreement between the surface hopping simulation results and the experimental measurements for the vibrational energy decay curves at different solvent densities and temperatures is shown to be good. The overlap of the decay curves when the time axis is linearly scaled is explained in terms of the perturbative theory for the rate constants. The contribution of each solvent atom to the change of the quantum populations of the solute molecule is used to analyze the mechanism of the relaxation process

The workings of a molecular thermometer: The vibrational excitation of carbon tetrachloride by a solvent

An intriguing energy-transfer experiment was recently carried out in methanol/carbon tetrachloride solutions. It turned out to be possible to watch vibrational energy accumulating in three of carbon tetrachloride's modes following initial excitation of O-H and C-H stretches in methanol, in effect making those CCl(4) modes "molecular thermometers" reporting on methanol's relaxation. In this paper, we use the example of a CCl(4) molecule dissolved in liquid argon to examine, on a microscopic level, just how this kind of thermal activation occurs in liquid solutions. The fact that even the lowest CCl(4) mode has a relatively high frequency compared to the intermolecular vibrational band of the solvent means that the only solute-solvent dynamics relevant to the vibrational energy transfer will be extraordinarily local, so much so that it is only the force between the instantaneously most prominent Cl and solvent atoms that will significantly contribute to the vibrational friction. We use this observation, within the context of a classical instantaneous-pair Landau-Teller calculation, to show that energy flows into CCl(4) primarily via one component of the nominally degenerate, lowest frequency, E mode and does so fast enough to make CCl(4) an excellent choice for monitoring methanol relaxation. Remarkably, within this theory, the different symmetries and appearances of the different CCl(4) modes have little bearing on how well they take up energy from their surroundings--it is only how high their vibrational frequencies are relative to the solvent intermolecular vibrational band edge that substantially favors one mode over another.

Correlation functions in classical gases at high frequency

A general procedure is outlined to determine the exact asymptotic form of spectra in classical gases at high frequency. Examples are the force-force correlation or the velocity autocorrelation function of a tagged particle. For purely repulsive potentials of the form Ar(-n), the asymptotic spectra are proportional to omega(sigma)exp[-(omegatau)(nu)]. Exponent nu and time constant tau depend only on the interparticle potential, while exponent sigma depends in addition on the correlation studied. The analysis makes use of the fact that the high-frequency spectra are dominated by high-energy binary collisions. It is argued that for arbitrary potentials the spectra decay slower than exp(-constxomega(2/3)) and that the results are also relevant for dense fluids. The frequency range is estimated where quantum effects become important.