Instantaneous Pair Theory for High-Frequency Vibrational Energy Relaxation in Fluids

Notwithstanding the long and distinguished history of studies of vibrational energy relaxation, exactly how it is that high frequency vibrations manage to relax in a liquid remains somewhat of a mystery. Both experimental and theoretical approaches seem to say that there is a natural frequency range associated with intermolecular motions in liquids, typically spanning no more than a few hundred cm^{-1}. Landau-Teller-like theories explain how a solvent can absorb any vibrational energy within this "band", but how is it that molecules can rid themselves of superfluous vibrational energies significantly in excess of these values? We develop a theory for such processes based on the idea that the crucial liquid motions are those that most rapidly modulate the force on the vibrating coordinate -- and that by far the most important of these motions are those involving what we have called the mutual nearest neighbors of the vibrating solute. Specifically, we suggest that whenever there is a single solvent molecule sufficiently close to the solute that the solvent and solute are each other's nearest neighbors, then the instantaneous scattering dynamics of the solute-solvent pair alone suffices to explain the high frequency relaxation. The many-body features of the liquid only appear in the guise of a purely equilibrium problem, that of finding the likelihood of particularly effective solvent arrangements around the solute. These results are tested numerically on model diatomic solutes dissolved in atomic fluids (including the experimentally and theoretically interesting case of I_2 in Xe). The instantaneous pair theory leads to results in quantitative agreement with those obtained from far more laborious exact molecular dynamics simulations.Comment: 55 pages, 6 figures Scheduled to appear in J. Chem. Phys., Jan, 199

Relationship between quantum decoherence times and solvation dynamics in condensed phase chemical systems

A relationship between the time scales of quantum coherence loss and short-time solvent response for a solute/bath system is derived for a Gaussian wave packet approximation for the bath. Decoherence and solvent response times are shown to be directly proportional to each other, with the proportionality coefficient given by the ratio of the thermal energy fluctuations to the fluctuations in the system-bath coupling. The relationship allows the prediction of decoherence times for condensed phase chemical systems from well developed experimental methods.Comment: 10 pages, no figures, late

Entropy, Dynamics and Instantaneous Normal Modes in a Random Energy Model

It is shown that the fraction f of imaginary frequency instantaneous normal modes (INM) may be defined and calculated in a random energy model(REM) of liquids. The configurational entropy S and the averaged hopping rate among the states R are also obtained and related to f, with the results R~f and S=a+b*ln(f). The proportionality between R and f is the basis of existing INM theories of diffusion, so the REM further confirms their validity. A link to S opens new avenues for introducing INM into dynamical theories. Liquid 'states' are usually defined by assigning a configuration to the minimum to which it will drain, but the REM naturally treats saddle-barriers on the same footing as minima, which may be a better mapping of the continuum of configurations to discrete states. Requirements of a detailed REM description of liquids are discussed

Nodal domains on quantum graphs

We consider the real eigenfunctions of the Schr\"odinger operator on graphs, and count their nodal domains. The number of nodal domains fluctuates within an interval whose size equals the number of bonds $B$. For well connected graphs, with incommensurate bond lengths, the distribution of the number of nodal domains in the interval mentioned above approaches a Gaussian distribution in the limit when the number of vertices is large. The approach to this limit is not simple, and we discuss it in detail. At the same time we define a random wave model for graphs, and compare the predictions of this model with analytic and numerical computations.Comment: 19 pages, uses IOP journal style file

Spectral Statistics of Instantaneous Normal Modes in Liquids and Random Matrices

We study the statistical properties of eigenvalues of the Hessian matrix ${\cal H}$ (matrix of second derivatives of the potential energy) for a classical atomic liquid, and compare these properties with predictions for random matrix models (RMM). The eigenvalue spectra (the Instantaneous Normal Mode or INM spectra) are evaluated numerically for configurations generated by molecular dynamics simulations. We find that distribution of spacings between nearest neighbor eigenvalues, s, obeys quite well the Wigner prediction $s exp(-s^2)$, with the agreement being better for higher densities at fixed temperature. The deviations display a correlation with the number of localized eigenstates (normal modes) in the liquid; there are fewer localized states at higher densities which we quantify by calculating the participation ratios of the normal modes. We confirm this observation by calculating the spacing distribution for parts of the INM spectra with high participation ratios, obtaining greater conformity with the Wigner form. We also calculate the spectral rigidity and find a substantial dependence on the density of the liquid.Comment: To appear in Phys. Rev. E; 10 pages, 6 figure

Local Variational Principle

A generalization of the Gibbs-Bogoliubov-Feynman inequality for spinless particles is proven and then illustrated for the simple model of a symmetric double-well quartic potential. The method gives a pointwise lower bound for the finite-temperature density matrix and it can be systematically improved by the Trotter composition rule. It is also shown to produce groundstate energies better than the ones given by the Rayleigh-Ritz principle as applied to the groundstate eigenfunctions of the reference potentials. Based on this observation, it is argued that the Local Variational Principle performs better than the equivalent methods based on the centroid path idea and on the Gibbs-Bogoliubov-Feynman variational principle, especially in the range of low temperatures.Comment: 15 pages, 5 figures, one more section adde

Quantum and classical vibrational relaxation dynamics of N-methylacetamide on ab initio potential energy surfaces

Employing extensive quantum-chemical calculations at the DFT/B3LYP and MP2 level, a quartic force field of isolated N-methylacetamide is constructed. Taking into account 24 vibrational degrees of freedom, the model is employed to perform numerically exact vibrational configuration interaction calculations of the vibrational energy relaxation of the amide I mode. It is found that the energy transfer pathways may sensitively depend on details of the theoretical description. Moreover, the exact reference calculations were used to study the applicability and accuracy of (i) the quasiclassical trajectory method, (ii) time-dependent second-order perturbation theory, and (iii) the instantaneous normal mode description of frequency fluctuations. Based on the results, several strategies to describe vibrational energy relaxation in biomolecular systems are discussed.Comment: 18 pages, 6 figures, submitted to J. Phys. Chem.

Mean-atom-trajectory model for the velocity autocorrelation function of monatomic liquids

We present a model for the motion of an average atom in a liquid or supercooled liquid state and apply it to calculations of the velocity autocorrelation function $Z(t)$ and diffusion coefficient $D$. The model trajectory consists of oscillations at a distribution of frequencies characteristic of the normal modes of a single potential valley, interspersed with position- and velocity-conserving transits to similar adjacent valleys. The resulting predictions for $Z(t)$ and $D$ agree remarkably well with MD simulations of Na at up to almost three times its melting temperature. Two independent processes in the model relax velocity autocorrelations: (a) dephasing due to the presence of many frequency components, which operates at all temperatures but which produces no diffusion, and (b) the transit process, which increases with increasing temperature and which produces diffusion. Because the model provides a single-atom trajectory in real space and time, including transits, it may be used to calculate all single-atom correlation functions.Comment: LaTeX, 8 figs. This is an updated version of cond-mat/0002057 and cond-mat/0002058 combined Minor changes made to coincide with published versio