Time evolution of density fluctuation in the supercritical region. 2. Comparison of hydrogen- and non-hydrogen-bonded fluids

The time evolution of the density fluctuation of molecules is investigated by dynamic light scattering in six neat fluids in supercritical states. This study is the first to compare the dynamics of density inhomogeneity between hydrogen- and non-hydrogen-bonded fluids. Supercritical methanol and ethanol are used as hydrogen-bonded fluids, whereas four non-hydrogen-bonded fluids were used: CHF(3), C(2)H(4), CO(2), and Xe. We measure the time correlation function of the density fluctuation of each fluid at the same reduced temperatures and densities and investigate the relationship between the dynamic and static density inhomogeneities of those supercritical fluids. In all cases, the profile of the time correlation function of the density fluctuation is characterized by a single-exponential function, whose decay is responsible for the dynamics characterized by hydrodynamic conditions. We obtain correlation times from the time correlation function and discuss dynamic and static inhomogeneity using the Kawasaki theory and the Landau-Placzek theory. While the correlation times in the six fluids show noncoincidence, those values agree well with each other except for the supercritical alcohols when scaled to a dimensionless parameter. Although the principle of corresponding state is observed in the non-hydrogen-bonded fluids, both the supercritical methanol and ethanol deviate from that principle. This deviation is attributed to the presence of hydrogen bonding among alcohol molecules at high temperature and low density. The average cluster size of each fluid is estimated under the same thermodynamic conditions, and it is shown that the clusters of supercritical alcohols are on average 1.5-1.7 times larger than those of the four non-hydrogen-bonded fluids. Moreover, the thermal diffusivity of each neat fluid is obtained over wide ranges of density and temperature.     

Polar solvation dynamics in supercritical fluids: a mode-coupling treatment

A mode-coupling treatment of polar solvation dynamics in supercritical fluids is presented. The equilibrium solvation time correlation function for the solute fluctuating transition frequency is obtained from the mode-coupling theory method and from molecular-dynamics simulations. The theory is shown to be in good agreement with the simulation. The solvation time correlation function exhibits three distinct time scales, with rapid initial decay, followed by a recurrence at intermediate times, and a slowly decaying long-time tail. Our theoretical analysis shows that the short-time decay arises from the coupling of the solute energy gap to the solvent polarization modes, the recurrence at intermediate times is due to the energy modes, while the slow long-time decay reflects the coupling to the number density modes.     

Theoretical Investigation on Triplet Excitation Energy Transfer in Fluorene Dimer

Triplet-triplet energy transfer in fluorene dimer is investigated by combining rate theories with electronic structure calculations. The two key parameters for the control of energy transfer, electronic coupling and reorganization energy, are calculated based on the diabatic states constructed by the constrained density functional theory. The fluctuation of the electronic coupling is further revealed by molecular dynamics simulation. Succeedingly, the diagonal and off-diagonal fluctuations of the Hamiltonian are mapped from the correlation functions of those parameters, and the rate is then estimated both from the perturbationtheory and wavepacket diffusion method. The results manifest that both the static and dynamic fluctuations enhance the rate significantly, but the rate from the dynamic fluctuation is smaller than that from the static fluctuation.     

Mean field kinetic theory for a lattice gas model of fluids confined in porous materials

We consider the mean field kinetic equations describing the relaxation dynamics of a lattice model of a fluid confined in a porous material. The dynamical theory embodied in these equations can be viewed as a mean field approximation to a Kawasaki dynamics Monte Carlo simulation of the system, as a theory of diffusion, or as a dynamical density functional theory. The solutions of the kinetic equations for long times coincide with the solutions of the static mean field equations for the inhomogeneous lattice gas. The approach is applied to a lattice gas model of a fluid confined in a finite length slit pore open at both ends and is in contact with the bulk fluid at a temperature where capillary condensation and hysteresis occur. The states emerging dynamically during irreversible changes in the chemical potential are compared with those obtained from the static mean field equations for states associated with a quasistatic progression up and down the adsorption/desorption isotherm. In the capillary transition region, the dynamics involves the appearance of undulates (adsorption) and liquid bridges (adsorption and desorption) which are unstable in the static mean field theory in the grand ensemble for the open pore but which are stable in the static mean field theory in the canonical ensemble for an infinite pore.     

The landau-placzek ratio for multicomponent fluids

Abstract

Under the assumption that the coupling between the sound modes and modes associated with heat and mass diffusion can be neglected, an expression for the Landau-Placzek ratio for multicomponent fluids is derived using thermodynamic fluctuation theory. Applications of the general formula to ternary systems are discussed briefly     

Coarse-grained molecular dynamics modeling of strongly associating fluids: thermodynamics, liquid structure, and dynamics of symmetric binary mixture fluids

Symmetric binary mixtures capable of strong association via a highly directional and saturable specific interaction between unlike molecules are investigated by canonical molecular dynamics simulations. The specific interaction of the molecules is defined in a new coarse-grained pair potential that can be applied in continuous molecular dynamics as well as in Monte Carlo simulations. The thermodynamic, structural, and dynamic properties of the associating mixture fluids are investigated as a function of density, temperature, and association strength of the specific interaction. Detailed analysis of the simulation data confirms a two-stage mechanism in the formation of specific bonds with increasing interaction strength, including a fast dimerization process and a subsequent stage of perfecting the bonds. A large heat capacity peak is found during the formation or breaking of the bonds, reflecting the large energy fluctuation introduced by the strong association. The fractions of nonbonded molecules obtained from the simulations as a function of density, temperature, and interaction strength are in excellent agreement with the predictions of Wertheim's thermodynamic perturbation theory. The translational and rotational dynamics of the Tmer mixture are effectively retarded with increasing association strength and are analyzed in terms of autocorrelation functions and a non-Gaussian parameter for the translational dynamics. The lifetimes of molecules in bonded and nonbonded states provide detailed information about the transformation of molecules between the bonded state and the nonbonded state. Finally, simulation sampling problems inherent to strongly interacting systems are easily overcome using the parallel tempering simulation technique. This latter result confirms that with the new continuous coarse-grained simulation potential we have a versatile and flexible interaction potential that can be used with many available molecular dynamics and Monte Carlo algorithms under various ensembles.     

Attractive and repulsive interactions among methanol molecules in supercritical state investigated by Raman spectroscopy and perturbed hard-sphere theory

The short-range structure of supercritical methanol (CH(3)OH) is investigated by measuring the spontaneous Raman spectra of the C-O stretching mode. The spectra are obtained at a reduced temperature, T(r)=T/T(c)=1.02 (522.9 K), which permits the neat fluid to be studied isothermally as a function of density. As the density increases, the spectral peaks shift toward the lower energy side and the spectra broaden. In the supercritical region, the amount of shifting shows nonlinear density dependence and the width becomes anomalously large. We use the perturbed hard-sphere model to analyze these density dependencies along the vibrational coordinate. The amount of shifting is decomposed into attractive and repulsive components, and the changes in attractive and repulsive energies are evaluated as functions of density and packing fraction, both of which are continuously varied by a factor of 120. Here we show that the shift amount consists principally of the attractive component at all densities, since the attractive energy is about eight times the repulsive energy. The density dependence of the widths is analyzed by calculating homogeneous and inhomogeneous widths as a function of density. The results show that, although vibrational dephasing and density inhomogeneity contribute similarly to the width at low and middle densities, at high density the main contributor turns out to be the vibrational dephasing. We estimate the local density enhancements of supercritical CH(3)OH as function of bulk density by two methods. The results of these analyses show common features, and both the estimated local density enhancements of CH(3)OH are considerably larger than the local density enhancements of simple fluids, i.e., those having nonhydrogen bonding. It is revealed that the local density of supercritical CH(3)OH is 40%-60% greater than the local densities of the simple fluids. We also estimate the local density fluctuation using the obtained values of attractive shift, inhomogeneous width, and local density. The density fluctuation in the vicinity of a vibrating molecule is compared to the fluctuation of bulk density, which is obtained from the thermodynamic calculation.     

Model for dynamics of inhomogeneous and bulk fluids

An accurate model for the density of states (DOS) for strongly inhomogeneous and bulk fluids has been proposed based on gamma distributions. The contribution to the density of states from the collective dynamics is modeled as an incomplete gamma distribution and the high frequency region is obtained from the solution of the memory equation using a sech memory kernel. Using only the frequency moments as input, the model parameters for the collective dynamics are obtained by matching moments of the resulting distribution. The model results in an analytical expression for the self-diffusivity of the fluid. We present results for soft sphere fluids confined in slit-shaped pores as well as bulk soft sphere liquids. Comparisons of the DOS, velocity autocorrelation functions, and memory kernels with molecular dynamics simulations reveal that the model predicts features in the DOS over the entire frequency range and is able to capture changes in the DOS as a function of fluid density and temperature. As a result the predicted VACFs, memory kernels, and self-diffusivities are accurately predicted over a wide range of conditions. Since the frequency moments for bulk liquids can be obtained from pair correlation functions, our method provides a direct route from fluid structure to dynamics. For fluids confined in slit-shaped pores, where the frequency moments are obtained from molecular dynamics simulations, the predicted self-diffusivities capture the resulting oscillations due to variations in the solvation pressure, and in the case of smooth walled pores, the predictions are superior to those obtained using kinetic theory.     

Local density augmentation and dynamic properties of hydrogen-and non-hydrogen-bonded supercritical fluids: a molecular dynamics study

The local density inhomogeneities in neat supercritical fluids were investigated via canonical molecular dynamics simulations. The selected systems under investigation were the polar and hydrogen-bonded fluid methanol as well as the quadrupolar non-hydrogen-bonded carbon dioxide one. Effective local densities, local density augmentation, and enhancement factors were calculated at state points along an isotherm close to the critical temperature of each system (T(r)=1.03). The results obtained reveal strong influence of the polarity and hydrogen bonding upon the intensity of the local density augmentation. It is found that this effect is sufficiently larger in the case of the polar and associated methanol in comparison to those predicted for carbon dioxide. For both fluids the local density augmentation values are maximized in the bulk density region near 0.7rho(c), a result that is in agreement with experiment. In addition, the local density dynamics of each fluid were investigated in terms of the appropriate time correlation functions. The behavior of these functions reveals that the bulk density dependence of the local density reorganization times is very sensitive to the specific intermolecular interactions and to the size of the local region. Also, the estimated local density reorganization time as a function of bulk density of each fluid was further analyzed and successfully related to two different time-scale relaxation mechanisms. Finally, the results obtained indicate a possible relationship between the single-molecule reorientational dynamics and the local density reorganization ones.     

Pair correlation functions and the self-diffusion coefficient of Lennard-Jones liquid in the modified free volume theory of diffusion

In this paper, we apply the Matteoli-Mansoori empirical formula for the pair correlation function of simple fluids obeying the Lennard-Jones potential to calculate reduced self-diffusion coefficients on the basis of the modified free volume theory. The self-diffusion coefficient thus computed as functions of temperature and density is compared with the molecular dynamics simulation data and the self-diffusion coefficient obtained by the modified free volume theory implemented with the Monte Carlo simulation method for the pair correlation function. We show that the Matteoli-Mansoori empirical formula yields sufficiently accurate self-diffusion coefficients in the supercritical regime, provided that the minimum free volume activating diffusion is estimated with the classical turning point of binary collision at the mean relative kinetic energy 3k(B)T/2, where k(B) is the Boltzmann constant and T is the temperature. In the subcritical regime, the empirical formula yields qualitatively correct, but lower values for the self-diffusion coefficients compared with computer simulation values and those from the modified free volume theory implemented with the Monte Carlo simulations for the pair correlation function. However, with a slightly modified critical free volume, the results can be made quite acceptable.     

Analysis to obtain precise density fluctuation of supercritical fluids by small-angle X-ray scattering

A procedure of analysis for small-angle X-ray scattering (SAXS) data has been established to obtain density fluctuation of supercritical fluids near the critical point. It is indispensable for the certain analysis to utilize both of high-quality SAXS data measured under stable thermodynamic condition and accurate P–ρ–T data in supercritical region. As a standard example, SAXS measurements have been performed for supercritical CO₂, which is a suitable sample satisfying the condition for both experiment and analysis. The measurements were carried out along four isothermal conditions at reduced temperature of T_r = T/T_c = 1.020, 1.022, 1.043 and 1.064. Comparing the experimental density fluctuation with calculated one from the most reliable equation of state, the differences are within 8% at most.     

Transition from patchlike to clusterlike inhomogeneity arising from hydrogen bonding in water

Assembling of water molecules via hydrogen bonding has been studied by molecular dynamics simulations using flexible potential model. The relationship between the number of H-bonds per molecule, n(HB), the size of H-bonded nets, k, and the size of patches of four-bonded molecules, k(4), has been examined for several thermodynamic states of water ranging from ambient to supercritical conditions. Two kinds of structural inhomogeneity have been found: the patchlike associated with the mean n(HB)> 2.0 and the clusterlike for n(HB)< 1.9. In compressed water up to ~473 K patches coexist with less ordered nets, both constituting the gel-like H-bonded network. The size of patches steeply decreases with the increasing temperature and the decreasing density of water. The inhomogeneity resulting from the presence of patches disappears above 473 K. This feature is associated with the rapid increase in the fraction of unbound molecules and with the breakage of the gel-like network into a variety of H-bonded clusters leading to the clusterlike structural inhomogeneity. In contrast to the patchlike inhomogeneity an increase in temperature and a decrease in density make this kind of inhomogeneity more pronounced. A degree of connectivity of H-bonds has been characterized by a parameter P(g) defined as the total fraction of molecules belonging to the H-bonded clusters of size k ≥ 5. The simulation-derived values of P(g) agree well with the predictions of the random bond theory giving the explicit expression for P(g) as a function of the mean n(HB). Going from ambient to supercritical conditions, we have found that the patchlike inhomogeneity is connected with the very slight reduction in P(g), whereas the clusterlike inhomogeneity generates a steep linear decrease of P(g) with the decreasing mean n(HB). The self-diffusion coefficient calculated for the thermodynamic states of water showing the clusterlike inhomogeneity has occurred to be inversely proportional to the density. We have also found that the clusterlike inhomogeneity is associated with the linear correlation between P(g) and the macroscopic properties of water: the static dielectric constant, the viscosity, and the density. The provided relationships allow one to estimate the degree of connectivity of hydrogen bonds from the measured macroscopic quantities.     

Generalized Langevin theory on the dynamics of simple fluids under external fields

A theory on the time development of the density and current fields of simple fluids under an external field is formulated through the generalized Langevin formalism. The theory is applied to the linear solvation dynamics of a fixed solute regarding the solute as the external field on the solvent. The solute-solvent-solvent three-body correlation function is taken into account through the hypernetted-chain integral equation theory, and the time correlation function of the random force is approximated by that in the absence of the solute. The theoretical results are compared with those of molecular-dynamics (MD) simulation and the surrogate theory. As for the transient response of the density field, our theory is shown to be free from the artifact of the surrogate theory that the solvent can penetrate into the repulsive core of the solute during the relaxation. We have also found a large quantitative improvement of the solvation correlation function compared with the surrogate theory. In particular, the short-time part of the solvation correlation function is in almost perfect agreement with that from the MD simulation, reflecting that the short-time expansion of the theoretical solvation correlation function is exact up to t(2) with the exact three-body correlation function. A quantitative improvement is found in the long-time region, too. Our theory is also applied to the force-force time correlation function of a fixed solute, and similar improvement is obtained, which suggests that our present theory can be a basis to improve the mode-coupling theory on the solute diffusion.     

Origin of large Landau-Placzek ratio in a liquid metallic alloy

Dynamic structure factors for Na(c)K(1-c) liquid metallic alloys and pure components are studied by molecular dynamics simulations. Large values of Landau-Placzek ratio for four compositions of the liquid alloy are analyzed by wave-number dependent contributions from relaxation and propagating processes within the generalized collective modes method. The origin of the large Landau-Placzek ratio for liquid alloys is discussed.     

Molecular dynamics simulation of the cybotactic region in gas-expanded methanol-carbon dioxide and acetone-carbon dioxide mixtures

Local solvation and transport effects in gas-expanded liquids (GXLs) are reported based on molecular simulation. GXLs were found to exhibit local density enhancements similar to those seen in supercritical fluids, although less dramatic. This approach was used as an alternative to a multiphase atomistic model for these mixtures by utilizing experimental results to describe the necessary fixed conditions for a locally (quasi-) stable molecular dynamics model of the (single) GXL phase. The local anisotropic pair correlation function, orientational correlation functions, and diffusion rates are reported for two systems: CO2-expanded methanol and CO2-expanded acetone at 298 K and pressures up to 6 MPa.     

Solvation shell dynamics studied by molecular dynamics simulation in relation to the translational and rotational dynamics of supercritical water and benzene

The solvation shell dynamics of supercritical water is analyzed by molecular dynamics simulation with emphasis on its relationship to the translational and rotational dynamics. The relaxation times of the solvation number (tau S), the velocity autocorrelation function (tau D), the angular momentum correlation function (tau J), and the second-order reorientational correlation function (tau 2R) are studied at a supercritical temperature of 400 degrees C over a wide density region of 0.01-1.5 g cm(-3). The relaxation times are decomposed into those conditioned by the solvation number n, and the effect of the short-ranged structure is examined in terms of its probability Pn of occurrence. In the low to medium-density range of 0.01-0.4 g cm(-3), the time scales of water dynamics are in the following sequence: tau D>tau S approximately or > tau J approximately or > tau 2R. This means that the rotation in supercritical water is of the "in-shell" type while the translational diffusion is not. The comparison to supercritical benzene is also performed and the effect of hydrogen bonding is examined. The water diffusion is not of the in-shell type up to the ambient density of 1.0 g cm(-3), which corresponds to the absence of the transition from the collision to the Brownian picture, whereas such transition is present in the case of benzene. The absence of the transition in water comes from the fast reorganization of the hydrogen bonds and the enhanced mobility of the solvation shell in supercritical conditions.     

Light scattering and dielectric relaxation studies of glass forming liquids

Extended Abstract: Glass forming organic liquids and polymers exhibit long range density fluctuations with correlation length ξ in the range of 10–300 nm at temperatures above T_g (1 - 6). This follows from dynamic and static light scattering experiments revealing some unexpected features, which cannot be explained on the basis of conventional liquid state theories: (i) In static light scattering the intensity I(q → 0) is no longer proportional to the isothermal compressibility, (ii) This excess scattering I_exc shows a strong q-dependence (q = (4π/Λ.)sin(θ/2)) corresponding to a correlation length ξ in the above mentioned range, (iii) The Landau-Placzek ratio I_Rayleigh/2I_Brillouin is much too high compared with the results of light scattering theories, (iv) In photon correlation spectroscopy a new ultraslow hydrodynamic mode (Γ ˜ q²) is detected with relaxation rates Γ about 10⁻⁶ to 10⁻⁹ lower than those of the α-process at a given temperature. In order to explain these observations, a two-state fluid model is proposed, which starts from the coexistence of “liquid-like” and “aperiodic solid-like” regions within the liquids. Such ideas have been discussed many times before, so for example A.R. Ubbelohde (7) speculates about “anticrystalline” clusters in liquids. Molecular dynamics simulations of atomic liquids showed that long range orientational fluctuations appear upon supercooling (8). A preferred icosahedral ordering is observed (9) and the number of icosahedral clusters increases with decreasing temperature (10). In connection with the interpretation of the dynamics of supercooled liquids different “two-state” models have been proposed (11 - 15). For the explanation of the light scattering results we propose that the molecules in the different dynamic states (“liquid” or “solid”) aggregate during annealing of the liquid at temperatures above T_g. Experiments showed that the equilibration times can be rather long (3 - 5), but nevertheless the liquids exhibiting long range density fluctuations are in the state of lowest free energy. We claim that our observations are the first experimental proof of the existence of such different dynamic states, which have been discussed many times before. The extended secondary clusters can also be detected by ultra small angle X-ray scattering.     

Solvation dynamics by coherence period resolved transient grating

We report a third-order nonlinear time-domain method, coherence period (τ) resolved transient grating (TRTG), that gives accurate solvation dynamics free from population relaxation in a short data acquisition time. The validity of TRTG is established by theory and experiment. The TRTG signal is shown to follow the transition frequency correlation function by an analytic expression based on the response function theory for delta function pulses and by model numerical calculations including finite pulse durations. TRTG is demonstrated for two cyanine dyes IR144 and IR125 in methanol by using a diffractive-optics based four wave mixing apparatus. Solvation dynamics in methanol obtained from the TRTG are consistent with those reported previously confirming the validity of TRTG.     

Transport properties of Mie(14,7) fluids: molecular dynamics simulation and theory

An extensive computer simulation study is presented for the self-diffusion coefficient, the shear viscosity, and the thermal conductivity of Mie(14,7) fluids. The time-correlation function formalism of Green-Kubo is utilized in conjunction with molecular dynamics (MD) simulations. In addition to molecular simulations, the results of a recent study [A. Eskandari Nasrabad, J. Chem. Phys. 128, 154514 (2008)] for the mean free volume are applied to calculate the self-diffusion coefficients within a free volume theory framework. A detailed comparison between the MD simulation and free volume theory results for the diffusion coefficient is given. The density fluctuation theory of shear viscosity is used to compute the shear viscosity and the results are compared to those from MD simulations. The density and temperature dependences of different time-correlation functions and transport coefficients are studied and discussed.