Self-diffusion near the liquid–vapor critical point

Transport coefficients and the spinodal of a fluid

Based on the results of computer simulation of transport coefficients in a Lennard-Jones fluid, it has been shown that in (p,T) coordinates, lines of constant values of coefficients of self-diffusion D and excess shear viscosity Δη extended into the region of superheated liquid and supersaturated vapor have an envelope. Each of the ensembles of such lines is divided into two parts by a critical isoline, where the conditions φ = φc = const (φ = D,Δη) are fulfilled. One part of this ensemble has as an envelope the branch of the spinodal that refers to superheated liquid; the other part has the branch of the spinodal of supersaturated vapor. The critical point of liquid-vapor equilibrium is the cusp of the spinodal, at which the critical isolines have a common tangent line with the spinodal and the binodal. Thus, at the approach to the spinodal, the relations (∂D/∂p)(T) → ∞ and (∂Δη/∂p)(T) → ∞ hold.

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.

Static and dynamic critical behavior of a symmetrical binary fluid: A computer simulation

A symmetrical binary, A+B Lennard-Jones mixture is studied by a combination of semi-grand-canonical Monte Carlo (SGMC) and molecular dynamics (MD) methods near a liquid-liquid critical temperature T(c). Choosing equal chemical potentials for the two species, the SGMC switches identities (A-->B-->A) to generate well-equilibrated configurations of the system on the coexistence curve for T<T(c) and at the critical concentration, x(c)=12, for T>T(c). A finite-size scaling analysis of the concentration susceptibility above T(c) and of the order parameter below T(c) is performed, varying the number of particles from N=400 to 12 800. The data are fully compatible with the expected critical exponents of the three-dimensional Ising universality class. The equilibrium configurations from the SGMC runs are used as initial states for microcanonical MD runs, from which transport coefficients are extracted. Self-diffusion coefficients are obtained from the Einstein relation, while the interdiffusion coefficient and the shear viscosity are estimated from Green-Kubo expressions. As expected, the self-diffusion constant does not display a detectable critical anomaly. With appropriate finite-size scaling analysis, we show that the simulation data for the shear viscosity and the mutual diffusion constant are quite consistent both with the theoretically predicted behavior, including the critical exponents and amplitudes, and with the most accurate experimental evidence.

Critical dynamics in a binary fluid: simulations and finite-size scaling

We report comprehensive simulations of the critical dynamics of a symmetric binary Lennard-Jones mixture near its consolute point. The self-diffusion coefficient exhibits no detectable anomaly. The data for the shear viscosity and the mutual-diffusion coefficient are fully consistent with the asymptotic power laws and amplitudes predicted by renormalization-group and mode-coupling theories provided finite-size effects and the background contribution to the relevant Onsager coefficient are suitably accounted for. This resolves a controversy raised by recent molecular simulations.

Friction of -bead macromolecules in solution: effects of the bead-solvent interaction

The role of the bead-solvent interaction has been studied for its influence on the dynamics of an N-bead macromolecule which is immersed into a solution. Using a Fokker-Planck equation for the phase-space distribution function of the macromolecule, we show that all the effects of the solution can be treated entirely in terms of the friction tensors which are assigned to each pair of interacting beads in the chain. For the high-density as well as for the critical solvent, the properties of these tensors are discussed in detail and are calculated by using several (realistic) choices of the bead-solvent potential. From the friction tensors, moreover, an expression for the center-of-mass friction coefficient of a (N-bead) chain macromolecule is derived. Numerical data for this coefficient for "truncated" Lennard-Jones bead-solvent potential are compared with results from molecular dynamic simulations and from the phenomenological theoretical data as found in the literature.

Low-frequency Raman spectra of sub- and supercritical CO2: qualitative analysis of the diffusion coefficient behavior

We report the results of the low-frequency Raman experiments on CO(2) which were carried out in a wide density range, along the liquid-gas coexistence curve in a temperature range of 293-303 K, and on the critical isochore of 94.4 cm(3) mol(-1) in a temperature range of 304-315 K. In our approach, the qualitative behavior of the diffusion coefficient D is predicted, assuming the following: first, that the low-frequency Raman spectra can be interpreted in terms of the translation rotation motions; second, that the random force could be replaced by the total force to calculate the friction coefficient; and finally, that the Einstein frequency is associated with the position of the maximum of the low-frequency Raman spectrum. The results show that the diffusion coefficient increases along the coexistence curve, and its values are almost constant on the critical isochore. The predicted values reproduce qualitatively those obtained by other techniques. The values of D were also calculated by molecular-dynamics simulation and they qualitatively reproduce the behavior of D.

Histogram reweighting method for dynamic properties

The histogram reweighting technique, widely used to analyze Monte Carlo data, is shown to be applicable to dynamic properties obtained from molecular dynamics simulations. The theory presented here is based on the fact that the correlation functions in systems in thermodynamic equilibrium are averages over initial conditions of trajectory functions, the latter depending on the volume of the system, the total number of particles, and the classical Hamiltonian. Thus, the well-known histogram reweighting method can be almost straightforwardly applied to reconstruct the probability distribution of initial states at different thermodynamic conditions, without extra computational effort. Correlation functions and transport coefficients are obtained with this method from few simulation data sets.

Dynamics of fluids near the consolute critical point: A molecular-dynamics study of the Widom–Rowlinson mixture

Molecular-dynamics simulations are presented for the dynamic behavior of the Widom-Rowlinson mixture [B. Widom, and J. S. Rowlinson, J. Chem. Phys. 52, 1670 (1970)] at its critical point. This model consists of two components where like species do not interact and unlike species interact via a hard-core potential. Critical exponents are obtained from a finite-size scaling analysis. The self-diffusion coefficient shows no anomalous behavior near the critical point. The shear viscosity and thermal conductivity show no divergent behavior for the system sizes considered, although there is a significant critical enhancement. The mutual diffusion coefficient, D(AB), vanishes as D(AB) approximately xi(-1.26 +/- 0.08), where xi is the correlation length. This is different from the renormalization-group (D(AB) approximately xi(-1.065)) mode coupling theory (D(AB) approximately xi(-1)) predictions. The theories and simulations can be reconciled if we assume that logarithmic corrections to scaling are important.

Chemical reaction dynamics within anisotropic solvents in time-dependent fields

The dynamics of low-dimensional Brownian particles coupled to time-dependent driven anisotropic heavy particles (mesogens) in a uniform bath (solvent) have been described through the use of a variant of the stochastic Langevin equation. The rotational motion of the mesogens is assumed to follow the motion of an external driving field in the linear response limit. Reaction dynamics have also been probed using a two-state model for the Brownian particles. Analytical expressions for diffusion and reaction rates have been developed and are found to be in good agreement with numerical calculations. When the external field driving the mesogens is held at constant rotational frequency, the model for reaction dynamics predicts that the applied field frequency can be used to control the product composition.

Brownian motion in dynamically disordered media

The motion of Brownian test particles in a model random potential with time dependent correlations is investigated using four methods: renormalized perturbation, perturbation using Martin, Siggia, and Rose functional formalism (MSR), the Edwards variational method on the MSR functional, and renormalization group with the MSR function. The disorder averaged one-particle propagators determined by the renormalized perturbation expansion and MSR perturbation expansion are identical to the second and possibly higher order, and the two-particle propagators determined by these perturbation methods are identical at the first and possibly higher order. The one-particle propagator determined by the Edwards method is identical to the perturbation expansions at the first order, but the second-order analogue of the Edwards method has a more complex expression, which reduces to the second-order perturbation expression with additional higher-order terms. The diffusion constant and two-particle correlations are calculated from these propagators and are used to determine the effects of the random potential on the Brownian particles. Generally, the diffusion rate decreases with the disorder strength and increases with the temporal decay rate. The two competing mechanisms result in an enhancement of the diffusion constant for weak potentials with fast temporal fluctuations. The system exhibits two-particle correlations that are inherently non-Gaussian and indicate clustering behavior. The diffusion constant is also determined from a simple one-loop renormalization group calculation. In the static limit, the diffusion constant calculated by the renormalization group recovers the results of Deem and Chandler [M.W. Deem and D. Chandler, J. Stat. Phys. 76, 911 (1994)].