1
|
|
2
|
Mutoru JW, Smith W, O’Hern CS, Firoozabadi A. Molecular dynamics simulations of diffusion and clustering along critical isotherms of medium-chain n-alkanes. J Chem Phys 2013; 138:024317. [DOI: 10.1063/1.4773282] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
3
|
Baidakov VG. Transport coefficients and the spinodal of a fluid. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 86:021201. [PMID: 23005752 DOI: 10.1103/physreve.86.021201] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/29/2012] [Indexed: 06/01/2023]
Abstract
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.
Collapse
Affiliation(s)
- Vladimir G Baidakov
- Institute of Thermal Physics, Ural Branch of the Russian Academy of Sciences, 106, Amundsen Street 106, 620016 Ekaterinburg, Russia.
| |
Collapse
|
4
|
Baidakov V, Kozlova Z. The self-diffusion coefficient in metastable states of a Lennard–Jones fluid. Chem Phys Lett 2010. [DOI: 10.1016/j.cplett.2010.09.066] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
5
|
|
6
|
Yoshida K, Matubayasi N, Nakahara M. Solvation shell dynamics studied by molecular dynamics simulation in relation to the translational and rotational dynamics of supercritical water and benzene. J Chem Phys 2008; 127:174509. [PMID: 17994829 DOI: 10.1063/1.2780871] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Ken Yoshida
- Institute for Chemical Research, Kyoto University, Uji, Kyoto 611-0011, Japan
| | | | | |
Collapse
|
7
|
Das SK, Horbach J, Binder K, Fisher ME, Sengers JV. Static and dynamic critical behavior of a symmetrical binary fluid: A computer simulation. J Chem Phys 2006; 125:24506. [PMID: 16848591 DOI: 10.1063/1.2215613] [Citation(s) in RCA: 75] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Subir K Das
- Institute for Physical Science and Technology, University of Maryland, College Park, MD 20742, USA
| | | | | | | | | |
Collapse
|
8
|
Das SK, Fisher ME, Sengers JV, Horbach J, Binder K. Critical dynamics in a binary fluid: simulations and finite-size scaling. PHYSICAL REVIEW LETTERS 2006; 97:025702. [PMID: 16907461 DOI: 10.1103/physrevlett.97.025702] [Citation(s) in RCA: 52] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/23/2006] [Indexed: 05/11/2023]
Abstract
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.
Collapse
Affiliation(s)
- Subir K Das
- Institute for Physical Science and Technology, University of Maryland, College Park, 20742, USA
| | | | | | | | | |
Collapse
|
9
|
Uvarov A, Fritzsche S. Friction of -bead macromolecules in solution: effects of the bead-solvent interaction. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2006; 73:011111. [PMID: 16486126 DOI: 10.1103/physreve.73.011111] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/11/2005] [Indexed: 05/06/2023]
Abstract
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.
Collapse
Affiliation(s)
- Alexander Uvarov
- Fachbereich Physik, Universität Kassel, Heinrich-Plett-Strasse 40, D-34132 Kassel, Germany.
| | | |
Collapse
|
10
|
Idrissi A, Longelin S, Damay P, Leclercq F. Low-frequency Raman spectra of sub- and supercritical CO2: qualitative analysis of the diffusion coefficient behavior. J Chem Phys 2005; 123:94501. [PMID: 16164347 DOI: 10.1063/1.1992476] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- A Idrissi
- Centre d'Etudes et de Recherches Lasers et Applications, Université des Sciences et Technologies de Lille, 59655 Villeneuve d'Ascq Cedex, France.
| | | | | | | |
Collapse
|
11
|
Nieto-Draghi C, Pérez-Pellitero J, Avalos JB. Histogram reweighting method for dynamic properties. PHYSICAL REVIEW LETTERS 2005; 95:040603. [PMID: 16090793 DOI: 10.1103/physrevlett.95.040603] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/19/2005] [Indexed: 05/03/2023]
Abstract
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.
Collapse
Affiliation(s)
- Carlos Nieto-Draghi
- Departament d'Enginyeria Química, ETSEQ, Universitat Rovira i Virgili, Tarragona, Spain
| | | | | |
Collapse
|
12
|
Jagannathan K, Yethiraj A. Dynamics of fluids near the consolute critical point: A molecular-dynamics study of the Widom–Rowlinson mixture. J Chem Phys 2005; 122:244506. [PMID: 16035781 DOI: 10.1063/1.1940053] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
|
13
|
A semi-phenomenological approach to the transport and diffusion of small spherical macromolecules in solution. Chem Phys Lett 2005. [DOI: 10.1016/j.cplett.2004.11.062] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
14
|
Hershkovits E, Hernandez R. Chemical reaction dynamics within anisotropic solvents in time-dependent fields. J Chem Phys 2005; 122:14509. [PMID: 15638676 DOI: 10.1063/1.1829252] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Eli Hershkovits
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, GA 30332-0400, USA
| | | |
Collapse
|
15
|
|
16
|
Masia M, Rey R. Reaction Rate Theory Approach to Thermodynamic State Dependence of Hydration Shell Exchange for Li+(Aq). J Phys Chem B 2003. [DOI: 10.1021/jp022271x] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Affiliation(s)
- Marco Masia
- Departament de Fısica i Enginyeria Nuclear, Universitat Politècnica de Catalunya, Campus Nord B4-B5, Barcelona 08034, Spain
| | - Rossend Rey
- Departament de Fısica i Enginyeria Nuclear, Universitat Politècnica de Catalunya, Campus Nord B4-B5, Barcelona 08034, Spain
| |
Collapse
|
17
|
Witkoskie JB, Yang S, Cao J. Brownian motion in dynamically disordered media. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2002; 66:051111. [PMID: 12513471 DOI: 10.1103/physreve.66.051111] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/22/2002] [Revised: 08/12/2002] [Indexed: 05/24/2023]
Abstract
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)].
Collapse
Affiliation(s)
- James B Witkoskie
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge 02139, USA
| | | | | |
Collapse
|
18
|
Patel N, Biswas R, Maroncelli M. Solvation and Friction in Supercritical Fluids: Simulation−Experiment Comparisons in Diphenyl Polyene/CO2 Systems. J Phys Chem B 2002. [DOI: 10.1021/jp0206727] [Citation(s) in RCA: 34] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- N. Patel
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802
| | - R. Biswas
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802
| | - M. Maroncelli
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802
| |
Collapse
|
19
|
Drozdov AN, Tucker SC. Response to “Comment on ‘Self-diffusion near the liquid–vapor critical point’ ” [J. Chem. Phys. 116, 6379 (2002)]. J Chem Phys 2002. [DOI: 10.1063/1.1458929] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
20
|
Harris KR. Comment on “Self-diffusion near the liquid–vapor critical point” [J. Chem. Phys. 114, 4912 (2001)]. J Chem Phys 2002. [DOI: 10.1063/1.1458928] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
21
|
Cherayil BJ. Anomalies in the self-diffusion coefficient near the critical point. J Chem Phys 2002. [DOI: 10.1063/1.1471549] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
22
|
|
23
|
Affiliation(s)
| | - Susan C. Tucker
- Department of Chemistry, University of California, Davis, California 95616
| |
Collapse
|