1
|
Lam MA, Khusid B, Kondic L, Meyer WV. Role of diffusion in crystallization of hard-sphere colloids. Phys Rev E 2021; 104:054607. [PMID: 34942784 DOI: 10.1103/physreve.104.054607] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2021] [Accepted: 10/12/2021] [Indexed: 11/07/2022]
Abstract
Vital for a variety of industries, colloids also serve as an excellent model to probe phase transitions at the individual particle level. Despite extensive studies, origins of the glass transition in hard-sphere colloids discovered about 30 y ago remain elusive. Results of our numerical simulations and asymptotic analysis suggest that cessation of long-time particle diffusivity does not suppress crystallization of a metastable liquid-phase hard-sphere colloid. Once a crystallite forms, its growth is then controlled by the particle diffusion in the depletion zone surrounding the crystallite. Using simulations, we evaluate the solid-liquid interface mobility from data on colloidal crystallization in terrestrial and microgravity experiments and demonstrate that there is no drastic difference between the respective mobility values. The insight into the effect of vanishing particle mobility and particle sedimentation on crystallization of colloids will help engineer colloidal materials with controllable structure.
Collapse
Affiliation(s)
- Michael A Lam
- New Jersey Institute of Technology, Newark, New Jersey 07102, USA
| | - Boris Khusid
- New Jersey Institute of Technology, Newark, New Jersey 07102, USA
| | - Lou Kondic
- New Jersey Institute of Technology, Newark, New Jersey 07102, USA
| | - William V Meyer
- Universities Space Research Association at NASA Glenn Research Center, Cleveland, Ohio 44135, USA
| |
Collapse
|
2
|
Mendoza-Méndez P, Lázaro-Lázaro E, Sánchez-Díaz LE, Ramírez-González PE, Pérez-Ángel G, Medina-Noyola M. Crossover from equilibration to aging: Nonequilibrium theory versus simulations. Phys Rev E 2017; 96:022608. [PMID: 28950613 DOI: 10.1103/physreve.96.022608] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2015] [Indexed: 06/07/2023]
Abstract
Understanding glasses and the glass transition requires comprehending the nature of the crossover from the ergodic (or equilibrium) regime, in which the stationary properties of the system have no history dependence, to the mysterious glass transition region, where the measured properties are nonstationary and depend on the protocol of preparation. In this work we use nonequilibrium molecular dynamics simulations to test the main features of the crossover predicted by the molecular version of the recently developed multicomponent nonequilibrium self-consistent generalized Langevin equation theory. According to this theory, the glass transition involves the abrupt passage from the ordinary pattern of full equilibration to the aging scenario characteristic of glass-forming liquids. The same theory explains that this abrupt transition will always be observed as a blurred crossover due to the unavoidable finiteness of the time window of any experimental observation. We find that within their finite waiting-time window, the simulations confirm the general trends predicted by the theory.
Collapse
Affiliation(s)
- P Mendoza-Méndez
- Instituto de Física "Manuel Sandoval Vallarta," Universidad Autónoma de San Luis Potosí, Álvaro Obregón 64, 78000 San Luis Potosí, SLP, México
| | - E Lázaro-Lázaro
- Instituto de Física "Manuel Sandoval Vallarta," Universidad Autónoma de San Luis Potosí, Álvaro Obregón 64, 78000 San Luis Potosí, SLP, México
| | - L E Sánchez-Díaz
- Department of Materials Science and Engineering, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - P E Ramírez-González
- CONACYT-Instituto de Física "Manuel Sandoval Vallarta," Universidad Autónoma de San Luis Potosí, Álvaro Obregón 64, 78000 San Luis Potosí, SLP, México
| | - G Pérez-Ángel
- Departamento de Física Aplicada CINVESTAV-IPN, Unidad Mérida Apartado Postal 73 Cordemex, 97310 Mérida, Yuc., México
| | - M Medina-Noyola
- Instituto de Física "Manuel Sandoval Vallarta," Universidad Autónoma de San Luis Potosí, Álvaro Obregón 64, 78000 San Luis Potosí, SLP, México
| |
Collapse
|
3
|
Callaham J, Machta J. Population annealing simulations of a binary hard-sphere mixture. Phys Rev E 2017; 95:063315. [PMID: 28709207 DOI: 10.1103/physreve.95.063315] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2017] [Indexed: 11/07/2022]
Abstract
Population annealing is a sequential Monte Carlo scheme well suited to simulating equilibrium states of systems with rough free energy landscapes. Here we use population annealing to study a binary mixture of hard spheres. Population annealing is a parallel version of simulated annealing with an extra resampling step that ensures that a population of replicas of the system represents the equilibrium ensemble at every packing fraction in an annealing schedule. The algorithm and its equilibration properties are described, and results are presented for a glass-forming fluid composed of a 50/50 mixture of hard spheres with diameter ratio of 1.4:1. For this system, we obtain precise results for the equation of state in the glassy regime up to packing fractions φ≈0.60 and study deviations from the Boublik-Mansoori-Carnahan-Starling-Leland equation of state. For higher packing fractions, the algorithm falls out of equilibrium and a free volume fit predicts jamming at packing fraction φ≈0.667. We conclude that population annealing is an effective tool for studying equilibrium glassy fluids and the jamming transition.
Collapse
Affiliation(s)
- Jared Callaham
- Department of Physics, University of Massachusetts, Amherst, Massachusetts 01003, USA
| | - Jonathan Machta
- Department of Physics, University of Massachusetts, Amherst, Massachusetts 01003, USA.,Santa Fe Institute, 1399 Hyde Park Road, Santa Fe, New Mexico 87501, USA
| |
Collapse
|
4
|
Ozawa M, Berthier L. Does the configurational entropy of polydisperse particles exist? J Chem Phys 2017; 146:014502. [DOI: 10.1063/1.4972525] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- Misaki Ozawa
- Laboratoire Charles Coulomb, UMR 5221 CNRS-Université de Montpellier, Montpellier, France
| | - Ludovic Berthier
- Laboratoire Charles Coulomb, UMR 5221 CNRS-Université de Montpellier, Montpellier, France
| |
Collapse
|
5
|
Berthier L, Coslovich D, Ninarello A, Ozawa M. Equilibrium Sampling of Hard Spheres up to the Jamming Density and Beyond. PHYSICAL REVIEW LETTERS 2016; 116:238002. [PMID: 27341260 DOI: 10.1103/physrevlett.116.238002] [Citation(s) in RCA: 84] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/08/2015] [Indexed: 06/06/2023]
Abstract
We implement and optimize a particle-swap Monte Carlo algorithm that allows us to thermalize a polydisperse system of hard spheres up to unprecedentedly large volume fractions, where previous algorithms and experiments fail to equilibrate. We show that no glass singularity intervenes before the jamming density, which we independently determine through two distinct nonequilibrium protocols. We demonstrate that equilibrium fluid and nonequilibrium jammed states can have the same density, showing that the jamming transition cannot be the end point of the fluid branch.
Collapse
Affiliation(s)
- Ludovic Berthier
- Laboratoire Charles Coulomb, UMR 5221 CNRS-Université de Montpellier, Montpellier 34095, France
| | - Daniele Coslovich
- Laboratoire Charles Coulomb, UMR 5221 CNRS-Université de Montpellier, Montpellier 34095, France
| | - Andrea Ninarello
- Laboratoire Charles Coulomb, UMR 5221 CNRS-Université de Montpellier, Montpellier 34095, France
| | - Misaki Ozawa
- Laboratoire Charles Coulomb, UMR 5221 CNRS-Université de Montpellier, Montpellier 34095, France
- Department of Physics, Nagoya University, Nagoya 464-8602, Japan
| |
Collapse
|
6
|
Jadrich R, Schweizer KS. Equilibrium theory of the hard sphere fluid and glasses in the metastable regime up to jamming. I. Thermodynamics. J Chem Phys 2013; 139:054501. [DOI: 10.1063/1.4816275] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
7
|
Bautista-Carbajal G, Moncho-Jordá A, Odriozola G. Further details on the phase diagram of hard ellipsoids of revolution. J Chem Phys 2013; 138:064501. [DOI: 10.1063/1.4789957] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
8
|
|
9
|
Pérez-Ángel G, Sánchez-Díaz LE, Ramírez-González PE, Juárez-Maldonado R, Vizcarra-Rendón A, Medina-Noyola M. Equilibration of concentrated hard-sphere fluids. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 83:060501. [PMID: 21797290 DOI: 10.1103/physreve.83.060501] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/02/2010] [Revised: 01/31/2011] [Indexed: 05/31/2023]
Abstract
We report a systematic molecular dynamics study of the isochoric equilibration of hard-sphere fluids in their metastable regime close to the glass transition. The thermalization process starts with the system prepared in a nonequilibrium state with the desired final volume fraction ϕ for which we can obtain a well-defined nonequilibrium static structure factor S(0)(k;ϕ). The evolution of the α-relaxation time τ(α)(k) and long-time self-diffusion coefficient D(L) as a function of the evolution time t(w) is then monitored for an array of volume fractions. For a given waiting time the plot of τ(α)(k;ϕ,t(w)) as a function of ϕ exhibits two regimes corresponding to samples that have fully equilibrated within this waiting time [ϕ≤ϕ(c)(t(w))] and to samples for which equilibration is not yet complete [ϕ≥ϕ(c)(t(w))]. The crossover volume fraction ϕ(c)(t(w)) increases with t(w) but seems to saturate to a value ϕ(a)≡ϕ(c)(t(w)→∞)≈0.582. We also find that the waiting time t(w)(eq)(ϕ) required to equilibrate a system grows faster than the corresponding equilibrium relaxation time, t(w)(eq)(ϕ)≈0.27[τ(α)(eq)(k;ϕ)](1.43), and that both characteristic times increase strongly as ϕ approaches ϕ(a), thus suggesting that the measurement of equilibrium properties at and above ϕ(a) is experimentally impossible.
Collapse
Affiliation(s)
- Gabriel Pérez-Ángel
- Departamento de Física Aplicada, Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional, Unidad Mérida, Mérida, Yucatán, Mexico
| | | | | | | | | | | |
Collapse
|
10
|
Odriozola G, Berthier L. Equilibrium equation of state of a hard sphere binary mixture at very large densities using replica exchange Monte Carlo simulations. J Chem Phys 2011; 134:054504. [DOI: 10.1063/1.3541248] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|