A study of molecular vibrational relaxation mechanism in condensed phase based upon mixed quantum-classical molecular dynamics. I. A test of IBC model for the relaxation of a nonpolar solute in nonpolar solvent at high density

In order to investigate vibrational relaxation mechanism in condensed phase, a series of mixed quantum-classical molecular dynamics calculations have been executed for nonpolar solute in nonpolar solvent and polar solute in polar solvent. In the first paper (Paper I), relaxation mechanism of I2 in Ar, where Lennard-Jones force is predominant in the interaction, is investigated as a function of density and temperature, focusing our attention on the isolated binary collision (IBC) model. The model was originally established for the relaxation in gas phase. A key question, here, is "can we apply the IBC model to the relaxation in the high-density fluid?" Analyzing the trajectory of solvent molecule as well as its interaction with the solute, we found that collisions between them may be defined clearly even in the high-density fluid. Change of the survival probability of the vibrationally first excited state on collision was traced. The change caused by collisions with a particular solvent molecule was also traced together with the interaction between them. Each collision makes a contribution to the relaxation by a stepwise change in the probability. The analysis clearly shows that the relaxation is caused by collisions even in the high-density fluid. The difference between stepwise relaxation and the continuous one found for the total relaxation in the low-density fluid and in the high-density one, respectively, was clarified to come from just the difference in frequency of the collision. The stronger the intensity of the collision is, the greater the relaxation caused by the collision is. Further, the shorter the collision time is, the greater the resultant relaxation is. The discussion is followed by the succeeding paper (Paper II), where we report that molecular mechanism of the relaxation of a polar molecule in supercritical water is significantly different from that assumed in the IBC model despite that the density dependence of the relaxation rate showed a linear correlation with the local density of water around the solute, the linear correlation being apparently in good accordance with the IBC model. The puzzle will be solved in Paper II.     

Umbrella sampling of solute vibrational line shifts in mixed quantum-classical molecular dynamics simulations

An umbrella sampling approach for vibrational frequency line shifts is presented. The technique allows for efficient sampling of the solvent configurations corresponding to frequency shifts of a solute in mixed quantum-classical simulations. The approach is generally applicable and can also be used within traditional perturbation theory calculations of frequency shifts. It is particularly useful in the extraction of detailed mechanistic information about the solute-solvent interactions giving rise to the frequency shifts. The method is illustrated by application to the simple I2 in a liquid Xe system, and the advantages are discussed.     

Mixed quantum-classical molecular dynamics simulation of vibrational relaxation of ions in an electrostatic field

The vibrational relaxation of ions in low-density gases under the action of an electrostatic field is reproduced through a molecular dynamics simulation method. The vibration is treated though quantum mechanics and the remaining degrees of freedom are considered classical. The procedure is tested through comparison against analytic results for a two-dimensional quantum model and by studying energy exchange during binary ion-atom collisions. Finally, the method has been applied successfully to the calculation of the mobility and the vibrational relaxation rate of O2+ in Kr as a function of the mean collision energy using a model interaction potential that reproduces the potential minimum of a previously known ab initio potential surface. The calculation of the steady mean vibrational motion of the ions in (flow) drift tubes seems straightforward, though at the expense of large amounts of computer time.     

Femtosecond quantum control of molecular dynamics in the condensed phase

We review the progress in controlling quantum dynamical processes in the condensed phase with femtosecond laser pulses. Due to its high particle density the condensed phase has both high relevance and appeal for chemical synthesis. Thus, in recent years different methods have been developed to manipulate the dynamics of condensed-phase systems by changing one or multiple laser pulse parameters. Single-parameter control is often achieved by variation of the excitation pulse's wavelength, its linear chirp or its temporal subpulse separation in case of pulse sequences. Multiparameter control schemes are more flexible and provide a much larger parameter space for an optimal solution. This is realized in adaptive femtosecond quantum control, in which the optimal solution is iteratively obtained through the combination of an experimental feedback signal and an automated learning algorithm. Several experiments are presented that illustrate the different control concepts and highlight their broad applicability. These fascinating achievements show the continuous progress on the way towards the control of complex quantum reactions in the condensed phase.     

A novel method for analyzing energy relaxation in condensed phases using nonequilibrium molecular dynamics simulations: application to the energy relaxation of intermolecular motions in liquid water

We present a novel method to investigate energy relaxation processes in condensed phases using nonequilibrium molecular dynamics simulations. This method can reveal details of the time evolution of energy relaxation like two-color third-order IR spectroscopy. Nonetheless, the computational cost of this method is significantly lower than that of third-order response functions. We apply this method to the energy relaxation of intermolecular motions in liquid water. We show that the intermolecular energy relaxation in water is characterized by four energy transfer processes. The structural changes of the liquid associated with the energy relaxation are also analyzed by the nonequilibrium molecular dynamics technique.     

Comparison between the Landau-Teller and flux-flux methods for computing vibrational energy relaxation rate constants in the condensed phase

The calculation of vibrational energy relaxation (VER) rate constants in the condensed phase is usually based on the Landau-Teller formula, which puts them in terms of the Fourier transform, at the vibrational frequency, of the autocorrelation function of the force exerted on the relaxing mode by the bath modes. An alternative expression for the VER rate constant puts it in terms of the autocorrelation function of the vibrational energy flux. In this paper, we compare the predictions obtained via those two methods in the case of iodine in liquid xenon. We find that the computational cost underlying both methods is comparable and that they predict similar VER rates. However, while the calculation of the VER rate via the Landau-Teller formula is somewhat more direct, the predictions obtained via the flux-flux formula are in somewhat better agreement with the VER rates obtained from nonequilibrium molecular dynamics simulations.     

Vibrational relaxation of OH and CH fundamentals of polar and nonpolar molecules in the condensed phase

Studies of vibrational energy flow in various polar and nonpolar molecules that follows the ultrafast excitation of the CH and OH stretch fundamentals, modeled using semiclassical methods, are reviewed. Relaxation rates are calculated using Landau-Teller theory and a time-dependent method, both of which consider a quantum mechanical solute molecule coupled to a classical bath of solvent molecules. A wide range of decay rates are observed, ranging from 1 ps for neat methanol to 50 ps for neat bromoform. In order to understand the flow rates, it is argued that an understanding of the subtle mixing between the solute eigenstates is needed and that solute anharmonicities are critical to facilitating condensed phase vibrational relaxation. The solvent-assisted shifts of the solute vibrational energy levels are seen to play a critical role of enhancing or decreasing lifetimes.     

Low-frequency vibrational spectrum of water in the hydration layer of a protein: a molecular dynamics simulation study

An atomistic molecular dynamics simulation has been carried out to understand the low-frequency intermolecular vibrational spectrum of water present in the hydration layer of the protein villin headpiece subdomain or HP-36. An attempt is made to explore how the heterogeneous rigidity of the hydration layers of different segments (three alpha helices) of the protein, strength of the protein-water hydrogen bonds, and their differential relaxation behavior influence the distribution of the intermolecular vibrational density of states of water in the hydration layers. The calculations revealed that compared to bulk water these bands are nonuniformly blue-shifted for water near the helices, the extent of shifts being more pronounced for water molecules hydrogen bonded to the protein residues. It is further noticed that the larger blue shift observed for the water molecules hydrogen bonded to helix 2 residues correlates excellently with the slowest structural relaxation of these hydrogen bonds. These results can be verified by suitable experimental measurements.     

A quantum equation of motion for chemical reaction systems on an adiabatic double-well potential surface in solution based on the framework of mixed quantum-classical molecular dynamics

We present a quantum equation of motion for chemical reaction systems on an adiabatic double-well potential surface in solution in the framework of mixed quantum-classical molecular dynamics, where the reactant and product states are explicitly defined by dividing the double-well potential into the reactant and product wells. The equation can describe quantum reaction processes such as tunneling and thermal excitation and relaxation assisted by the solvent. Fluctuations of the zero-point energy level, the height of the barrier, and the curvature of the well are all included in the equation. Here, the equation was combined with the surface hopping technique in order to describe the motion of the classical solvent. Applying the present method to model systems, we show two numerical examples in order to demonstrate the potential power of the present method. The first example is a proton transfer by tunneling where the high-energy product state was stabilized very rapidly by solvation. The second example shows a thermal activation mechanism, i.e., the initial vibrational excitation in the reactant well followed by the reacting transition above the barrier and the final vibrational relaxation in the product well.     

A comparative study of the properties of polar and nonpolar solvent/solute/polystyrene solutions in microwave fields via molecular dynamics

The influence of an applied microwave field on the dynamics of methylamine-dichloromethane (DCM) mixtures bound within atactic polystyrene (a-PS) over a range of polymer densities from 30 to 94 wt % polymer was examined using atomistic molecular dynamics simulations. This study is an extension of previous studies on methylamine transport in relatively polar polystyrene solutions of methanol and dimethylformamide [M. J. Purdue et al., J. Chem. Phys. 124, 204904 (2006)]. A direct comparison is made across the three types of polystyrene solutions. Consideration is given to both solvent and reagent transport within the polymer solutions under zero-field conditions and in an external electromagnetic field in the canonical ensemble (NVT) at 298.0 K. Various frequencies up to 10(4) GHz and a rms electric field intensity of 0.1 VA were applied. The simulation studies were validated by comparison of the simulated zero-field self-diffusion coefficients of DCM in a-PS with those obtained using pulsed-gradient spin-echo NMR spectrometry. Athermal effects of microwave fields on solute transport behavior within polymer solutions are discussed.     

Energy relaxation of intermolecular motions in supercooled water and ice: a molecular dynamics study

We investigate the energy relaxation of intermolecular motions in liquid water at temperatures ranging from 220 K to 300 K and in ice at 220 K using molecular dynamics simulations. We employ the recently developed frequency resolved transient kinetic energy analysis, which provides detailed information on energy relaxation in condensed phases like two-color pump-probe spectroscopy. It is shown that the energy cascading in liquid water is characterized by four processes. The temperature dependences of the earlier three processes, the rotational-rotational, rotational-translational, and translational-translational energy transfers, are explained in terms of the density of states of the intermolecular motions. The last process is the slow energy transfer arising from the transitions between potential energy basins caused by the excitation of the low frequency translational motion. This process is absent in ice because the hydrogen bond network rearrangement, which accompanies the interbasin transitions in liquid water, cannot take place in the solid phase. We find that the last process in supercooled water is well approximated by a stretched exponential function. The stretching parameter, β, decreases from 1 to 0.72 with decreasing temperature. This result indicates that the dynamics of liquid water becomes heterogeneous at lower temperatures.     

Density dependence of hydrogen bonding and the translational-orientational structural order in supercritical water: a molecular dynamics study

Molecular dynamics simulation have been performed with a wide range of densities along a near critical isotherm of supercritical water (SCW) in order to study the density dependence of the structure order and hydrogen bonding (HB). It is revealed that the translational structure order is nearly invariant while the orientational tetrahedral structure order is very sensitive to the bulk density under supercritical conditions. Meanwhile, some energetically unfavorable intermediate water dimer structures are found to appear under supercritical conditions due to the reduced energy difference and the enhanced energy fluctuation. As a consequence, a general geometrical criterion or the inclusion of a energy-based criterion instead of currently widely adopted pure r(OH)-based geometric criterion is suggested to be used in the HB statistics under supercritical conditions. It is found that the average HB number per H(2)O molecule (n(HB)) reduces with the decreasing SCW bulk density although a given pair of H(2)O molecules are shown to have a stronger ability to form a hydrogen bond under lower SCW bulk densities. Accordingly, the orientational tetrahedral structure order q decreases with the reducing bulk density under supercritical conditions. However, when the fluid is dilute with ρ ≤ 0.19ρ(c) (ρ(c) = 0.322 g/cm(3)), the energy fluctuation increases sharply and the short-range order is destroyed, signifying the supercritical fluid (SCF)-gas state transition. Accordingly, the orientational tetrahedral structure order q gets reversal around ρ = 0.19ρ(c) and approaches zero under very dilute conditions. The sensitivity of the orientational order to the density implies the microscopic origin of the significant dependence of SCF's physicochemical properties on the pressure.     

A molecular dynamics study of structural relaxation in tetrahedrally coordinated nanocrystals

The reorganisation of nanocrystals in order to reduce their surface energies has been examined in computer simulations. The relaxation takes a qualitatively different path for sphalerite- and wurtzite-structured particles. The surfaces of the sphalerite particles reconstruct into hexagonal nets, but the interior remains identifiable as sphalerite-like, whereas wurtzite particles form facetted, hexagonal nanorods by virtue of a reorganisation of the whole particle which involves the creation of a low energy internal interface between oppositely oriented domains. Despite the reorganisation, the diffraction patterns remain compatible with a wurtzite structure with some internal strain. The dipole moments of thermalized wurtzite particles are compared with experimental results for CdSe.     

Dependence of diffusivity on density and solute diameter in liquid phase: a molecular dynamics study of Lennard-Jones system

Investigations into the variation of self-diffusivity with solute radius, density, and degree of disorder of the host medium is explored. The system consists of a binary mixture of a relatively smaller sized solute, whose size is varied and a larger sized solvent interacting via Lennard-Jones potential. Calculations have been performed at three different reduced densities of 0.7, 0.8, and 0.933. These simulations show that diffusivity exhibits a maximum for some intermediate size of the solute when the solute diameter is varied. The maximum is found at the same size of the solute at all densities which is at variance with the prediction of the levitation effect. In order to understand this anomaly, additional simulations were carried out in which the degree of disorder has been varied while keeping the density constant. The results show that the diffusivity maximum gradually disappears with increase in disorder. Disorder has been characterized by means of the minimal spanning tree. Simulations have also been carried out in which the degree of disorder is constant and only the density is altered. The results from these simulations show that the maximum in diffusivity now shifts to larger distances with decrease in density. This is in agreement with the changes in void and neck distribution with density of the host medium. These results are in excellent agreement with the predictions of the levitation effect. They suggest that the effect of disorder is to shift the maximum in diffusivity towards smaller solute radius while that of the decrease in density is to shift it towards larger solute radius. Thus, in real systems where the degree of disorder is lower at higher density and vice versa, the effect due to density and disorder have opposing influences. These are confirmed by the changes seen in the velocity autocorrelation function, self part of the intermediate scattering function and activation energy.     

Ultrafast solvent response upon a change of the solute size in non-polar supercritical fluids

Non-polar solvation dynamics has been investigated using steady-state absorption and emission spectroscopy of the NO A ²Σ⁺(3sσ) Rydberg state in fluid Ar over a wide range of densities spanning the supercritical regime. Equilibrium molecular dynamics simulations were implemented to derive a new isotropic NO A(3sσ)–Ar pair potential which was further used to investigate the role of local density enhancements on the solvation process by non-equilibrium molecular dynamics simulations. These density inhomogeneities were found to have no influence on the solvation dynamics. Furthermore, the latter was shown to take place in a strongly non-linear regime, especially at low temperatures. This process results from the dramatic change of solute–solvent short range interaction associated with the large solute size change upon excitation to the Rydberg state.