Diffusion Anomaly in Silicalite and VPI-5 from Molecular Dynamics Simulations

Molecular dynamics (MD) simulations on rigid and flexible framework models of silicalite and a rigid framework model of the aluminophosphate VPI-5 for different sorbate diameters are reported. The sorbate-host interactions are modeled in terms of simple atom-atom Lennard-Jones interactions. The results suggest that the diffusion coefficient exhibits an anomaly as gamma approaches unity. The MD results confirm the existence of a linear regime for sorbate diameters significantly smaller than the channel diameter and an anomalous regime observed for sorbate diameters comparable to the channel diameter. The power spectra obtained by Fourier transformation of the velocity autocorrelation function indicate that there is an increase in the intensity of the low-frequency component for the velocity component parallel to the direction of motion for the sorbate diameter in the anomalous regime. The present results suggest that the diffusion anomaly is observed irrespective of (1) the geometry and topology of the pore structure and (2) the nature of the host material. The results are compared with the work of Derouane and co-workers, who have suggested the existence of ''floating molecules'' on the basis of earlier theoretical and computational approaches.

Ultrafast solvation dynamics in water: Isotope effects and comparison with experimental results

A detailed theoretical study of solvation dynamics in water is presented. The motivation of the present study comes from the recent experimental observation that the dynamics of solvation of an ion in water is ultrafast and the solvation time correlation function decays with a time constant of about 55 fs. The slower decay in the long time can be described by a sum of two exponentials with time constants equal to 126 and 880 fs. The molecular theory (developed earlier) predicts a time constant equal to 52 fs for the initial Gaussian decay and time constants equal to 134 and 886 fs for the two exponential components at the long time. This nearly perfect agreement is obtained by using the most detailed dynamical information available in the literature. The present study emphasizes the importance of the intermolecular vibrational band originating from the O...O stretching mode of the O�H...O units in the initial dynamics and raises several interesting questions regarding the nature of the decay of this mode. We have also studied the effects of isotope substitution on solvation dynamics. It is predicted that a significant isotope effect may be observed in the long time. The experimental results have also been compared with the prediction of the dynamic mean spherical approximation (DMSA); the agreement is not satisfactory at the long time. It is further found that the molecular theory and the DMSA lead to virtually identical results if the translational modes of the solvent molecules are neglected in the former. DMSA has also been used to investigate the dynamics of solvation of a dipolar solute in water. It is found that the dynamics of dipolar solvation exhibit features rather different from those of ion solvation. © 1995 American Institute of Physics.

Ultrafast Solvation Dynamics Of An Ion In A Restricted Environment

The solvation dynamics of an excited coumarin dye molecule (C-480) enclosed within a restricted space have been studied using molecular hydrodynamic theory (MHT) and compared with the recent experimental findings. The solvation dynamics of the dye molecule within the cavity of a toroidal gamma-cyclodextrin molecule have been shown to be explained only in terms of the freezing of the solvent translational modes using MHT. The results of the theoretical calculation are in good agreement with the experimental results. The inertial components of the solvation time correlation function remain the same in both the restricted environment and in the free space. These results are interesting in the light of the simulation studies of Maroncelli and Fleming [J chem Phys, 89 (1988) 5044] which concludes that the participation of the different solvation shells in controlling the dynamics are much different. The earlier studies have been reviewed and the recent findings are discussed.

Monte Carlo and molecular dynamics simulation of argon clusters andn-alkanes in the confined regions of zeolites

Geometry and energy of argon clusters confined in zeolite NaCaA are compared with those of free clusters. Results indicate the possible existence of magic numbers among the confined clusters. Spectra obtained from instantaneous normal mode analysis of free and confined clusters give a larger percentage of imaginary frequencies for the latter indicating that the confined cluster atoms populate the saddle points of the potential energy surface significantly. The variation of the percentage of imaginary frequencies with temperature during melting is akin to the variation of other properties. It is shown that confined clusters might exhibit inverse surface melting, unlike medium-to-large-sized free clusters that exhibit surface melting. Configurational-bias Monte Carte (CBMC) simulations of n-alkanes in zeolites Y and A are reported. CBMC method gives reliable estimates of the properties relating to the conformation of molecules. Changes in the conformational properties of n-butane and other longer n-alkanes such as n-hexane and n-heptane when they are confined in different zeolites are presented. The changes in the conformational properties of n-butane and n-hexane with temperature and concentration is discussed. In general, in zeolite Y as well as A, there is significant enhancement of the gauche population as compared to the pure unconfined fluid.

A molecular dynamics study and molecular level explanation of pressure dependence of ionic conductivity of potassium chloride in water

Experimental ionic conductivity of different alkali ions in water shows markedly different dependences on pressure. Existing theories such as that of Hubbard-Onsager are unable to explain these dependences on pressure of the ionic conductivity for all ions. We report molecular dynamics investigation of potassium chloride solution at low dilution in water at several pressures between 1 bar and 2 kbar. Two different potential models have been employed. One of the models successfully reproduces the experimentally observed trend in ionic conductivity of K+ ions in water over the 0.001-2 kbar range. We also propose a theoretical explanation, albeit at a qualitative level, to account for the dependence of ionic conductivity on pressure in terms of the previously studied Levitation Effect. It also provides a microscopic picture in terms of the pore network in liquid water.

Solvation dynamics in monohydroxy alcohols: Agreement between theory and different experiments

Recently three different experimental studies on ultrafast solvation dynamics in monohydroxy straight-chain alcohols (C-1-C-4) have been carried out, with an aim to quantify the time constant (and the amplitude) of the ultrafast component. The results reported are, however, rather different from different experiments. In order to understand the reason for these differences, we have carried out a detailed theoretical study to investigate the time dependent progress of solvation of both an ionic and a dipolar solute probe in these alcohols. For methanol, the agreement between the theoretical predictions and the experimental results [Bingemann and Ernsting J. Chem. Phys. 1995, 102, 2691 and Horng et al. J: Phys, Chern, 1995, 99, 17311] is excellent. For ethanol, propanol, and butanol, we find no ultrafast component of the time constant of 70 fs or so. For these three liquids, the theoretical results are in almost complete agreement with the experimental results of Horng et al. For ethanol and propanol, the theoretical prediction for ionic solvation is not significantly different from that of dipolar solvation. Thus, the theory suggests that the experiments of Bingemann and Ernsting and those of Horng et al. studied essentially the polar solvation dynamics. The theoretical studies also suggest that the experimental investigations of Joo et al. which report a much faster and larger ultrafast component in the same series of solvents (J. Chem. Phys. 1996, 104, 6089) might have been more sensitive to the nonpolar part of solvation dynamics than the polar part. In addition, a discussion on the validity of the present theoretical approach is presented. In this theory the ultrafast component arises from almost frictionless inertial motion of the individual solvent molecules in the force field of its neighbors.

Solvation dynamics of a charge bubble in water

A microscopic theory is used to calculate the solvation-time correlation function, (S(t)), of a light, non-stationary charge bubble in water. The calculated correlation function is found to be similar to the energy-time correlation function of a solvated electron. The ionic mobility of a charge bubble of the size of the hydrated electron is also calculated. It is found that the mobility of the charge plays a very important role in its own solvation.

Structure and silver ion transport in AgI-Ag2MoO4 glasses: A molecular dynamics study

The structure and dynamics of silver ion conducting AgI-Ag2MoO4 glasses have been simulated by molecular dynamics simulation over a wide range of compositions. Formation of silver iodide like aggregates have been identified only in the AgI rich glasses. Increase in silver ion conductivity with an increase in AgI content in the glass is seen as in experiments. The dynamics of ion transport suggests that Ag+ ion transport occurs largely through paths connected by silver ion sites of mixed iodide-oxide coordination. The Van Hove correlation functions indicate that Ag+ ions prefer migration along the pathways formed with connected sites of similar coordination.

Isomerization dynamics in viscous liquids: Microscopic investigation of the coupling and decoupling of the rate to and from solvent viscosity and dependence on the intermolecular potential

A detailed investigation of viscosity dependence of the isomerization rate is carried out for continuous potentials by using a fully microscopic, self-consistent mode-coupling theory calculation of both the friction on the reactant and the viscosity of the medium. In this calculation we avoid approximating the short time response by the Enskog limit, which overestimates the friction at high frequencies. The isomerization rate is obtained by using the Grote-Hynes formula. The viscosity dependence of the rate has been investigated for a large number of thermodynamic state points. Since the activated barrier crossing dynamics probes the high-frequency frictional response of the liquid, the barrier crossing rate is found to be sensitive to the nature of the reactant-solvent interaction potential. When the solute-solvent interaction is modeled by a 6-12 Lennard-Jones potential, we find that over a large variation of viscosity (eta), the rate (k) can indeed be fitted very well to a fractional viscosity dependence: (k similar to eta(-alpha)), with the exponent alpha in the range 1 greater than or equal to alpha >0. The calculated values of the exponent appear to be in very good agreement with many experimental results. In particular, the theory, for the first time, explains the experimentally observed high value of alpha even at the barrier frequency, omega(b). similar or equal to 9 X 10(12) s(-1) for the isomerization reaction of 2-(2'-propenyl)anthracene in liquid eta-alkanes. The present study can also explain the reason for the very low value of vb observed in another study for the isomerization reaction of trans-stilbene in liquid n-alkanes. For omega(b) greater than or equal to 2.0 X 10(13) s(-1), we obtain alpha similar or equal to 0, which implies that the barrier crossing rate becomes identical to the transition-state theory predictions. A careful analysis of isomerization reaction dynamics involving large amplitude motion suggests that the barrier crossing dynamics itself may become irrelevant in highly viscous liquids and the rate might again be coupled directly to the viscosity. This crossover is predicted to be strongly temperature dependent and could be studied by changing the solvent viscosity by the application of pressure. (C) 1999 American Institute of Physics. [S0021-9606(9950514-X].

Polar and nonpolar solvation dynamics, ion diffusion, and vibrational relaxation: Role of biphasic solvent response in chemical dynamics

ChemInform is a weekly Abstracting Service, delivering concise information at a glance that was extracted from about 100 leading journals. To access a ChemInform Abstract of an article which was published elsewhere, please select a “Full Text” option. The original article is trackable via the “References” option.

Solvation dynamics in nonassociated polar solvents

Ultrafast solvation dynamics in three nonassociated polar solvents, namely, acetonitrile, dimethyl sulfoxide, and acetone, have been studied by using the molecular hydrodynamic theory. For solvation in acetonitrile, the solvent memory function required for this study has been obtained from recent dielectric relaxation measurements of Venabales and Schuttenmaer; earlier theoretical studies used only the Kerr relaxation data. As the latter provides only an indirect information regarding the polar dynamical response of the dipolar liquid, it fails to provide a fully quantitative description of the solvation time correlation function, S(t). The present study with full dielectric data, on the other hand, gives excellent agreement with the experimental results. The theory shows that the ultrafast part of the solvation dynamics originates almost entirely from the high-frequency component of dielectric relaxation (with time constant 0.177 ps), although the latter represents only a small part of the latter. For DMSO and acetone, however, the present theory predicts a decay slower than the experimental observation. It is proposed that for these two solvents specific chromophore-solvent interactions might be responsible for the-large discrepancy. On the basis of the theory, two experimental studies have also been proposed.

Solvation dynamics of a quadrupolar solute in dipolar liquids

A microscopic calculation of solvation dynamics of dipolar and quadrupolar solutes in liquid water and acetonitrile is presented. The solvation is found to he biphasic. The calculated solvation time correlation function of ionic quadrupolar solute (K+) in water is in good agreement with re cent computer simulation results. Present study reveals some interesting aspects of quadrupolar solvation dynamics which differ significantly from that of ionic and dipolar solvation.

Subquadratic quantum number dependence and other anomalies of vibrational dephasing in liquid nitrogen: Molecular dynamics simulation study from the triple point to the critical point and beyond

Vibrational phase relaxation near gas-liquid and liquid-solid phase coexistence has been studied by molecular dynamics simulations of N-N stretch in N-2. Experimentally observed pronounced insensitivity of phase relaxation from the triple point to beyond the boiling point is found to originate from a competition between density relaxation and resonant-energy transfer terms. The sharp rise in relaxation rate near the critical point (CP) can be attributed at least partly to the sharp, rise in vibration-rotation coupling contribution. Substantial subquadratic quantum number dependence of overtone dephasing rate is found near the CP and in supercritical fluids. [S0031-9007 (99)09318-7].

The molecular dynamics study of lithium ion conduction in phosphate glasses and the role of non-bridging oxygen

Molecular dynamics (MD) simulation of lithium phosphate (Li2O-P2O5) glasses with varying Li2O content has been carried out. Two different P-O distances corresponding to phosphorus coordination with bridging oxygen (BO) and non-bridging oxygen (NBO) were identified in the simulated glasses. NBO-BO interconversion or bond switching was noted, which results in a dynamic equilibration of the tetrahedral phosphate units (P-n, n = 1,3 indicates the number of bridging oxygen atoms in the coordination of phosphorus). The NBO-BO bond switching is mildly activated with an effective activation barrier of 0.03-0.05 eV. Lithium ion jumps do not appear to be strongly coupled to bond switching. But the number of Li+ ions coordinated to an optimum number of NBOs and the number of Li+ ions jumping out of their sites appear to be correlated. Detailed analysis was made of the dynamics of P-n species and new insights have been obtained regarding ion migration in network-modified phosphate glasses.

Biological water: Femtosecond dynamics of macromolecular hydration

The unique features of a macromolecule and water as a solvent make the issue of solvation unconventional, with questions about the static versus dynamic nature of hydration and the, physics of orientational and translational diffusion at the boundary. For proteins, the hydration shell that covers the surface is critical to the stability of its structure and function. Dynamically speaking, the residence time of water at the surface is a signature of its mobility and binding. With femtosecond time resolution it is possible to unravel the shortest residence times which are key for the description of the hydration layer, static or dynamic. In this article we review these issues guided by experimental studies, from this laboratory, of polar hydration dynamics at the surfaces of two proteins (Subtilisin Carlsberg (SC) and Monellin). The natural probe tryptophan amino acid was used for the interrogation of the dynamics, and for direct comparison we also studied the behavior in bulk water - a complete hydration in 1 ps. We develop a theoretical description of solvation and relate the theory to the experimental observations. In this - theoretical approach, we consider the dynamical equilibrium in the hydration shell, defining the rate processes for breaking and making the transient hydrogen bonds, and the effective friction in the layer which is defined by the translational and orientational motions of water molecules. The relationship between the residence time of water molecules and the observed slow component in solvation dynamics is a direct one. For the two proteins studied, we observed a "bimodal decay" for the hydration correlation function, with two primary relaxation times: ultrafast, typically 1 ps or less, and longer, typically 15-40 ps, and both are related to the residence time at the protein surface, depending on the binding energies. We end by making extensions to studies of the denatured state of the protein, random coils, and the biomimetic micelles, and conclude with our thoughts on the relevance of the dynamics of native structures to their functions.