1
|
Midya J, Das SK. Finite-size scaling study of dynamic critical phenomena in a vapor-liquid transition. J Chem Phys 2017; 146:044503. [PMID: 28147549 DOI: 10.1063/1.4974506] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Via a combination of molecular dynamics (MD) simulations and finite-size scaling (FSS) analysis, we study dynamic critical phenomena for the vapor-liquid transition in a three dimensional Lennard-Jones system. The phase behavior of the model has been obtained via the Monte Carlo simulations. The transport properties, viz., the bulk viscosity and the thermal conductivity, are calculated via the Green-Kubo relations, by taking inputs from the MD simulations in the microcanonical ensemble. The critical singularities of these quantities are estimated via the FSS method. The results thus obtained are in nice agreement with the predictions of the dynamic renormalization group and mode-coupling theories.
Collapse
Affiliation(s)
- Jiarul Midya
- Theoretical Sciences Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur P.O., Bangalore 560064, India
| | - Subir K Das
- Theoretical Sciences Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur P.O., Bangalore 560064, India
| |
Collapse
|
2
|
Roy S, Dietrich S, Höfling F. Structure and dynamics of binary liquid mixtures near their continuous demixing transitions. J Chem Phys 2016; 145:134505. [DOI: 10.1063/1.4963771] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
3
|
Roy S, Das SK. Finite-size scaling study of shear viscosity anomaly at liquid-liquid criticality. J Chem Phys 2015; 141:234502. [PMID: 25527943 DOI: 10.1063/1.4903810] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We study the equilibrium dynamics of a symmetrical binary Lennard-Jones fluid mixture near its consolute criticality. Molecular dynamics simulation results for the shear viscosity, η, from a microcanonical ensemble are compared with those from a canonical ensemble with various thermostats. It is observed that the Nosé-Hoover thermostat is a good candidate for this purpose, and is therefore adopted for the quantification of the critical singularity of η, to avoid the temperature fluctuations (or even drifts) that are often encountered in microcanonical simulations. Via a finite-size scaling analysis of our simulation data we have been able to confirm that the shear viscosity exhibits a weak critical singularity in agreement with the theoretical predictions.
Collapse
Affiliation(s)
- Sutapa Roy
- Theoretical Sciences Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur P.O., Bangalore 560064, India
| | - Subir K Das
- Theoretical Sciences Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur P.O., Bangalore 560064, India
| |
Collapse
|
4
|
Roy S, Das SK. Simulation of transport around the coexistence region of a binary fluid. J Chem Phys 2014; 139:064505. [PMID: 23947869 DOI: 10.1063/1.4817777] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We use Monte Carlo and molecular dynamics simulations to study phase behavior and transport properties in a symmetric binary fluid where particles interact via Lennard-Jones potential. Our results for the critical behavior of collective transport properties, with particular emphasis on bulk viscosity, is understood via appropriate application of finite-size scaling technique. It appears that the critical enhancements in these quantities are visible far above the critical point. This result is consistent with an earlier report from computer simulations where, however, the authors do not quantify the critical singularity.
Collapse
Affiliation(s)
- Sutapa Roy
- Theoretical Sciences Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur P.O., Bangalore 560064, India
| | | |
Collapse
|
5
|
Gross M, Varnik F. Critical dynamics of an isothermal compressible nonideal fluid. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 86:061119. [PMID: 23367905 DOI: 10.1103/physreve.86.061119] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/12/2012] [Revised: 08/23/2012] [Indexed: 06/01/2023]
Abstract
A pure fluid at its critical point shows a dramatic slow-down in its dynamics, due to a divergence of the order-parameter susceptibility and the coefficient of heat transport. Under isothermal conditions, however, sound waves provide the only possible relaxation mechanism for order-parameter fluctuations. Here we study the critical dynamics of an isothermal, compressible nonideal fluid via scaling arguments and computer simulations of the corresponding fluctuating hydrodynamics equations. We show that, below a critical dimension of 4, the order-parameter dynamics of an isothermal fluid effectively reduces to "model A," characterized by overdamped sound waves and a divergent bulk viscosity. In contrast, the shear viscosity remains finite above two dimensions. Possible applications of the model are discussed.
Collapse
Affiliation(s)
- Markus Gross
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, Universitätsstr 90a, 44789 Bochum, Germany.
| | | |
Collapse
|
6
|
Gross M, Varnik F. Simulation of static critical phenomena in nonideal fluids with the lattice Boltzmann method. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 85:056707. [PMID: 23004903 DOI: 10.1103/physreve.85.056707] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/19/2011] [Indexed: 06/01/2023]
Abstract
A fluctuating nonideal fluid at its critical point is simulated with the lattice Boltzmann method. It is demonstrated that the method, employing a Ginzburg-Landau free energy functional, correctly reproduces the static critical behavior associated with the Ising universality class. A finite-size scaling analysis is applied to determine the critical exponents related to the order parameter, compressibility and specific heat. A particular focus is put on finite-size effects and issues related to the global conservation of the order parameter.
Collapse
Affiliation(s)
- M Gross
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, Germany.
| | | |
Collapse
|
7
|
|
8
|
Zeng M, Mi J, Zhong C. Density Functional Theory Integrated with Renormalization Group Theory for Criticality of Nanoconfined Fluids. J Phys Chem B 2010; 114:3894-901. [DOI: 10.1021/jp911070a] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Ming Zeng
- Laboratory of Computational Chemistry, Department of Chemical Engineering Beijing University of Chemical Technology, Beijing 100029, China
| | - Jianguo Mi
- Laboratory of Computational Chemistry, Department of Chemical Engineering Beijing University of Chemical Technology, Beijing 100029, China
| | - Chongli Zhong
- Laboratory of Computational Chemistry, Department of Chemical Engineering Beijing University of Chemical Technology, Beijing 100029, China
| |
Collapse
|
9
|
Das SK, Sengers JV, Fisher ME. Simulating critical dynamics in liquid mixtures: Short-range and long-range contributions. J Chem Phys 2007; 127:144506. [DOI: 10.1063/1.2770736] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
10
|
Radhakrishnan A, McConnell H. Composition fluctuations, chemical exchange, and nuclear relaxation in membranes containing cholesterol. J Chem Phys 2007; 126:185101. [PMID: 17508832 DOI: 10.1063/1.2730805] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A thermodynamic model of cholesterol-phospholipid complexes is used as a starting point for calculating fluctuations in membranes containing cholesterol and phospholipids. The calculations describe fluctuations in the concentration of complexes formed between cholesterol and phospholipids with longer saturated fatty acid chains. The fluctuations in complex concentrations arise by two distinct mechanisms. In one, the chemical composition of the sample varies from point to point, and the concentration of the complexes varies according to local chemical equilibrium. In the second, the composition remains fixed, and the complexes form and dissociate according to chemical reaction kinetics. In both cases the nuclear resonance frequency of a deuterium labeled phospholipid undergoes fluctuations and line broadening as a consequence of the formation and dissociation of complexes. For a specific ternary lipid mixture at its critical composition, deuterium nuclear resonance line broadening of chain labeled phospholipids is calculated for temperatures up to 10 degrees -20 degrees above the miscibility critical temperature. This line broadening is associated with fluctuations in the degree of phospholipid chain ordering related to the formation and dissociation of complexes.
Collapse
Affiliation(s)
- Arun Radhakrishnan
- Department of Molecular Genetics, University of Texas Southwestern Medical Center, Dallas, Texas 75390, USA
| | | |
Collapse
|
11
|
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
|
12
|
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
|
13
|
|
14
|
Chen A, Chimowitz EH, De S, Shapir Y. Universal dynamic exponent at the liquid-gas transition from molecular dynamics. PHYSICAL REVIEW LETTERS 2005; 95:255701. [PMID: 16384471 DOI: 10.1103/physrevlett.95.255701] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/08/2003] [Revised: 02/02/2005] [Indexed: 05/05/2023]
Abstract
The liquid-gas system is expected to exhibit distinct dynamic behavior in the fluid's critical region (model H). We present molecular dynamics simulations of a Lennard-Jones fluid model starting from specially designed, near-equilibrium, initial conditions. By following the fluid's relaxation towards equilibrium, we calculate the requisite transport coefficients in the critical region. The results yield the scaling behavior of the thermal diffusion coefficient D(T) approximately xi(-1.023+/-0.018) (xi is the correlation length) and a nonconventional divergent heat conductivity, all of which are in accord with mode-coupling and renormalization group predictions, as well as some experimental data.
Collapse
Affiliation(s)
- A Chen
- Department of Chemical Engineering, University of Rochester, Rochester, New York 14627, USA
| | | | | | | |
Collapse
|
15
|
Santos A, López de Haro M. Demixing can occur in binary hard-sphere mixtures with negative nonadditivity. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2005; 72:010501. [PMID: 16089927 DOI: 10.1103/physreve.72.010501] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/05/2004] [Indexed: 05/03/2023]
Abstract
A binary fluid mixture of nonadditive hard spheres characterized by a size ratio gamma = sigma(2)/sigma(1) < 1 and a nonadditivity parameter Delta = 2 sigma(12)/(sigma(1) + sigma(2)) - 1 is considered in infinitely many dimensions. From the equation of state in the second virial approximation (which is exact in the limit d--> infinity) a demixing transition with a critical consolute point at a packing fraction scaling as eta approximately d2(-d) is found, even for slightly negative nonadditivity, if Delta >-1/8 (ln gamma)(2). Arguments concerning the stability of the demixing with respect to freezing are provided.
Collapse
Affiliation(s)
- A Santos
- Departamento de Física, Universidad de Extremadura, Badajoz, Spain.
| | | |
Collapse
|
16
|
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
|
17
|
Sengers JV, Moldover MR. Comment on "molecular dynamics simulations of a fluid near its critical point". PHYSICAL REVIEW LETTERS 2005; 94:069601. [PMID: 15783786 DOI: 10.1103/physrevlett.94.069601] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/26/2004] [Indexed: 05/24/2023]
Affiliation(s)
- J V Sengers
- University of Maryland College Park, Maryland 20742, USA and National Institute of Standards and Technology Gaithersburg, Maryland 20899-8360, USA
| | | |
Collapse
|
18
|
Santos A, López de Haro M, Yuste SB. Equation of state of nonadditive d-dimensional hard-sphere mixtures. J Chem Phys 2005; 122:024514. [PMID: 15638605 DOI: 10.1063/1.1832591] [Citation(s) in RCA: 37] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An equation of state for a multicomponent mixture of nonadditive hard spheres in d dimensions is proposed. It yields a rather simple density dependence and constitutes a natural extension of the equation of state for additive hard spheres proposed by us [A. Santos, S. B. Yuste, and M. Lopez de Haro, Mol. Phys. 96, 1 (1999)]. The proposal relies on the known exact second and third virial coefficients and requires as input the compressibility factor of the one-component system. A comparison is carried out both with another recent theoretical proposal based on a similar philosophy and with the available exact results and simulation data in d=1, 2, and 3. Good general agreement with the reported values of the virial coefficients and of the compressibility factor of binary mixtures is observed, especially for high asymmetries and/or positive nonadditivities.
Collapse
Affiliation(s)
- A Santos
- Departamento de Física, Universidad de Extremadura, E-06071 Badajoz, Spain.
| | | | | |
Collapse
|
19
|
Abstract
In this paper, we present a cluster algorithm for the numerical simulations of nonadditive hard-core mixtures. This algorithm allows one to simulate and equilibrate systems with a number of particles two orders of magnitude larger than previous simulations. The phase separation for symmetric binary mixtures is studied for different nonadditivities as well as for the Widom-Rowlinson model [B. Widom and J. S. Rowlinson, J. Chem. Phys. 52, 1670 (1970)] in two and three dimensions. The critical densities are determined from finite size scaling. The critical exponents for all the nonadditivities are consistent with the Ising universality class.
Collapse
Affiliation(s)
- Arnaud Buhot
- UMR 5819 (UJF, CNRS, CEA) DRFMC/SI3M, CEA Grenoble, 17 rue des Martyrs, 38054 Grenoble cedex 9, France
| |
Collapse
|