Rotational energy relaxation of individual rotational states in liquids

Rotational energy relaxation quantum dynamics of a diatomic molecule in a superfluid helium nanodroplet and study of the hydrogen isotopes case

The rotational energy relaxation (RER) of a molecule X2(j,mj) in a 4He superfluid nanodroplet [HeND or (4He)N; T = 0.37 K] has been investigated using a hybrid quantum dynamics approach recently proposed by us. As far as we know, this is the first theoretical study about rotational relaxation inside HeNDs, and here several (real and hypothetical) isotopes of H2 have been examined, in order to analyze the influence of the rotational constant Be of these fast rotors on the dynamics. The structure of the nanodroplet practically does not change during the RER process, which approximately takes place according to a cascade mechanism j → j - 2; j - 2 → j - 4; …; 2 → 0, and mj is conserved. The results are consistent with the very scarce estimated experimental data available. The lifetime of an excited rotational state (≈1.0-7.6 ns) increases when: (a) Be increases; (b) j increases; and (c) N decreases (above N = 100 there is a small influence of N on the lifetime). This also applies to the global relaxation time and transition time. The analysis of the influence of the coupling between the j and j - 2 rotational states (due to the X2-helium interaction) and the X2 angular velocity on the lifetime and related properties has been helpful to better understand the dynamics. In contrast to the RER results, for the vibrational energy relaxation (VER) in HeNDs, when the quantum number v increases a decrease is observed in the lifetime of the excited vibrational state. This difference can be interpreted taking into account that RER and VER are associated with very different types of motion. Besides, in VER the intermediate excited states show metastability, differing from the RER case. We hope that the present study will encourage more studies to be developed on the RER dynamics in HeNDs, a basic, interesting and difficult to study physical phenomenon about which we still know very little.

Local structural effects on orientational relaxation of OH-bond in liquid water over short to intermediate timescales

By simulating the rigid simple point charge extended model at temperature T = 300 K, the orientational relaxation of the OH-bond in water was investigated over short to intermediate timescales, within which molecules undergo inertial rotation and libration and then enter the rotational diffusion regime. According to the second-cumulant approximation, the orientational time correlation function (TCF) of each axis that is parallel or perpendicular to an OH-bond is related to an effective rotational density of states (DOS), which is determined using the power spectra of angular velocity autocorrelation functions (AVAFs) of the other two axes. In addition, the AVAF power spectrum of an axis was approximated as the rotational stable instantaneous normal mode (INM) spectrum of the axis. As described in a previous study [S. L. Chang, T. M. Wu, and C. Y. Mou, J. Chem. Phys. 121, 3605 (2004)], simulated molecules were classified into subensembles, according to either the local structures or the H-bond configurations of the molecules. For global molecules and the classified subensembles, the simulation results for the first- and second-rank orientational TCFs were compared with the second-cumulant predictions obtained using the effective rotational DOSs and the rotational stable-INM spectra. On short timescales, the OH-bond in water behaves similar to an inertial rotor and its anisotropy is lower than that of a water molecule. For molecules with three or more H-bonds, the OH-bond orientational TCFs are characterized by a recurrence, which is an indication for libration of the OH-bond. The recurrence can generally be described by the second-cumulant prediction obtained using the rotational stable-INM spectra; however, the orientational TCFs after the recurrence switch to a behavior similar to that predicted using the AVAF power spectra. By contrast, the OH-bond orientational TCFs of molecules initially connected with one or two H-bonds decay monotonically or exhibit a weak recurrence, indicating rapid relaxation into the rotational diffusion regime after the initial Gaussian decay. In addition to accurately describing the Gaussian decay, the second-cumulant predictions formulated using the rotational stable-INM spectra and the AVAF power spectra serve as the upper and lower limits, respectively, for the OH-bond orientational TCFs of these molecules after the Gaussian decay.

Vibrational Energy Relaxation Rates of H2 and D2 in Liquid Argon via the Linearized Semiclassical Method

The vibrational energy relaxation (VER) rates for H2 and D2 in liquid argon (T=152 K, rho=1.45x1022 cm-3) are calculated using the linearized semiclassical (LSC) method (J. Phys. Chem. 2003, 107, 9059, 9070). The calculation is based on Fermi's golden rule. The VER rate constant is expressed in terms of the quantum-mechanical force-force correlation function, which is then estimated using the LSC method. A local harmonic approximation (LHA) is employed in order to compute the multidimensional Wigner integrals underlying the LSC approximation. The H2-Ar and D2-Ar interactions are described by the three-body potential of Bissonette et al. (J. Phys. Chem. A 1996, 105, 2639). The LHA-LSC-based VER rate constants for both D2 and H2 are found to be about 2-3 orders of magnitude slower than those obtained experimentally. However, their ratio agrees quantitatively with the corresponding experimental result. In contrast, the classical VER rate constants are found to be 8-9 orders of magnitude slower than those obtained experimentally, and their ratio is found to be qualitatively different from the corresponding experimental result.

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.

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.

Intense laser alignment in dissipative media as a route to solvent dynamics

We extend the concept of alignment by short intense pulses to dissipative environments within a density matrix formalism and illustrate the application of this method as a probe of the dissipative properties of dense media. In particular, we propose a means of disentangling rotational population relaxation from decoherence effects via strong laser alignment. We illustrate also the possibility of suppressing rotational relaxation to prolong the alignment lifetime through choice of the field parameters. Implications to several disciplines and a number of potential applications are proposed.

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.