1
|
Hu YC, Tanaka H. Revealing the role of liquid preordering in crystallisation of supercooled liquids. Nat Commun 2022; 13:4519. [PMID: 35927419 PMCID: PMC9352720 DOI: 10.1038/s41467-022-32241-z] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2022] [Accepted: 07/21/2022] [Indexed: 11/09/2022] Open
Abstract
The recent discovery of non-classical crystal nucleation pathways has revealed the role of fluctuations in the liquid structural order, not considered in classical nucleation theory. On the other hand, classical crystal growth theory states that crystal growth is independent of interfacial energy, but this is questionable. Here we elucidate the role of liquid structural ordering in crystal nucleation and growth using computer simulations of supercooled liquids. We find that suppressing the crystal-like structural order in the supercooled liquid through a new order-killing strategy can reduce the crystallisation rate by several orders of magnitude. This indicates that crystal-like liquid preordering and the associated interfacial energy reduction play an essential role in nucleation and growth processes, forcing critical modifications of the classical crystal growth theory. Furthermore, we evaluate the importance of this additional factor for different types of liquids. These findings shed new light on the fundamental understanding of crystal growth kinetics.
Collapse
Affiliation(s)
- Yuan-Chao Hu
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo, 153-8505, Japan
| | - Hajime Tanaka
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo, 153-8505, Japan. .,Research Center for Advanced Science and Technology, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo, 153-8904, Japan.
| |
Collapse
|
2
|
Banerjee A, Sevilla M, Rudzinski JF, Cortes-Huerto R. Finite-size scaling and thermodynamics of model supercooled liquids: long-range concentration fluctuations and the role of attractive interactions. SOFT MATTER 2022; 18:2373-2382. [PMID: 35258066 DOI: 10.1039/d2sm00089j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
We compute partial structure factors, Kirkwood-Buff integrals (KBIs) and chemical potentials of model supercooled liquids with and without attractive interactions. We aim at investigating whether relatively small differences in the tail of the radial distribution functions result in contrasting thermodynamic properties. Our results suggest that the attractive potential favours the nucleation of long-range structures. Indeed, upon decreasing temperature, Bathia-Thornton structure factors display anomalous behaviour in the k→0 limit. KBIs extrapolated to the thermodynamic limit confirm this picture, and excess coordination numbers identify the anomaly with long-range concentration fluctuations. By contrast, the purely repulsive system remains perfectly miscible for the same temperature interval and only reveals qualitatively similar concentration fluctuations in the crystalline state. Furthermore, differences in both isothermal compressibilities and chemical potentials show that thermodynamics is not entirely governed by the short-range repulsive part of the interaction potential, emphasising the nonperturbative role of attractive interactions. Finally, at higher density, where both systems display nearly identical dynamical properties and repulsive interactions become dominant, the anomaly disappears, and both systems also exhibit similar thermodynamic properties.
Collapse
Affiliation(s)
- Atreyee Banerjee
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany.
| | - Mauricio Sevilla
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany.
| | - Joseph F Rudzinski
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany.
| | | |
Collapse
|
3
|
Machine-Learned Free Energy Surfaces for Capillary Condensation and Evaporation in Mesopores. ENTROPY 2022; 24:e24010097. [PMID: 35052123 PMCID: PMC8774451 DOI: 10.3390/e24010097] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/15/2021] [Revised: 12/29/2021] [Accepted: 01/05/2022] [Indexed: 12/04/2022]
Abstract
Using molecular simulations, we study the processes of capillary condensation and capillary evaporation in model mesopores. To determine the phase transition pathway, as well as the corresponding free energy profile, we carry out enhanced sampling molecular simulations using entropy as a reaction coordinate to map the onset of order during the condensation process and of disorder during the evaporation process. The structural analysis shows the role played by intermediate states, characterized by the onset of capillary liquid bridges and bubbles. We also analyze the dependence of the free energy barrier on the pore width. Furthermore, we propose a method to build a machine learning model for the prediction of the free energy surfaces underlying capillary phase transition processes in mesopores.
Collapse
|
4
|
Desgranges C, Delhommelle J. Entropy scaling close to criticality: From simple to metallic systems. Phys Rev E 2021; 103:052102. [PMID: 34134262 DOI: 10.1103/physreve.103.052102] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2020] [Accepted: 04/09/2021] [Indexed: 11/07/2022]
Abstract
Entropy has recently drawn considerable interest both as a marker to detect the onset of phase transitions and as a reaction coordinate, or collective variable, to span phase transition pathways. We focus here on the behavior of entropy along the vapor-liquid phase coexistence and identify how the difference in entropy between the two coexisting phases vary in ideal and metallic systems along the coexistence curve. Using flat-histogram simulations, we determine the thermodynamic conditions of coexistence, critical parameters, including the critical entropy, and entropies along the binodal. We then apply our analysis to a series of systems that increasingly depart from ideality and adopt a metal-like character, through the gradual onset of the Friedel oscillation in an effective pair potential, and for a series of transition metals modeled with a many-body embedded-atoms force field. Projections of the phase boundary on the entropy-pressure and entropy-temperature planes exhibit two qualitatively different behaviors. While all systems modeled with an effective pair potential lead to an ideal-like behavior, the onset of many-body effects results in a departure from ideality and a markedly greater exponent for the variation of the entropy of vaporization with temperature away from the critical temperature.
Collapse
Affiliation(s)
- Caroline Desgranges
- Department of Chemistry, New York University, New York, New York 10003, USA and Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| | - Jerome Delhommelle
- Department of Chemistry, New York University, New York, New York 10003, USA and Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| |
Collapse
|
5
|
Desgranges C, Delhommelle J. The central role of entropy in adiabatic ensembles and its application to phase transitions in the grand-isobaric adiabatic ensemble. J Chem Phys 2020; 153:094114. [PMID: 32891099 DOI: 10.1063/5.0021488] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Entropy has become increasingly central to characterize, understand, and even guide assembly, self-organization, and phase transition processes. In this work, we build on the analogous role of partition functions (or free energies) in isothermal ensembles and that of entropy in adiabatic ensembles. In particular, we show that the grand-isobaric adiabatic (μ, P, R) ensemble, or Ray ensemble, provides a direct route to determine the entropy. This allows us to follow the variations of entropy with the thermodynamic conditions and thus explore phase transitions. We test this approach by carrying out Monte Carlo simulations on argon and copper in bulk phases and at phase boundaries. We assess the reliability and accuracy of the method through comparisons with the results from flat-histogram simulations in isothermal ensembles and with the experimental data. Advantages of the approach are multifold and include the direct determination of the μ-P relation, without any evaluation of pressure via the virial expression, the precise control of the system size (number of atoms) via the input value of R, and the straightforward computation of enthalpy differences for isentropic processes, which are key quantities to determine the efficiency of thermodynamic cycles. A new insight brought by these simulations is the highly symmetric pattern exhibited by both systems along the transition, as shown by scaled temperature-entropy and pressure-entropy plots.
Collapse
Affiliation(s)
- Caroline Desgranges
- Department of Chemistry, New York University, New York, New York 10003, USA and Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| | - Jerome Delhommelle
- Department of Chemistry, New York University, New York, New York 10003, USA and Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| |
Collapse
|
6
|
Liang Y, Díaz Leines G, Drautz R, Rogal J. Identification of a multi-dimensional reaction coordinate for crystal nucleation in Ni3Al. J Chem Phys 2020; 152:224504. [DOI: 10.1063/5.0010074] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Affiliation(s)
- Yanyan Liang
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Grisell Díaz Leines
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Ralf Drautz
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Jutta Rogal
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
7
|
Desgranges C, Delhommelle J. Nucleation of Capillary Bridges and Bubbles in Nanoconfined CO 2. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2019; 35:15401-15409. [PMID: 31675236 DOI: 10.1021/acs.langmuir.9b01744] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Using molecular simulation, we examine the capillary condensation and the capillary evaporation of CO2 in cylindrical nanopores. More specifically, we employ the recently developed μV T-S method to determine the microscopic mechanism associated with these processes and the corresponding free energy profiles. We calculate the free energy barrier for capillary condensation and identify that the key step consists in the nucleation of a liquid bridge of a critical size. Similarly, the free energy maximum found for the capillary evaporation process is found to correspond to the nucleation of a vapor bubble of a critical size. In addition, we assess the impact of the strength of the wall-fluid on the height of the free energy barrier and on the critical size of liquid bridges (condensation process) and vapor bubbles (evaporation process). We observe that the height of the free energy barrier increases with the strength of the wall-fluid interactions. Finally, we build a theoretical model, based on capillary theory, to rationalize our findings. In particular, the simulation results reveal a linear scaling of the free energy barrier with the critical size, in excellent agreement with the theoretical predictions.
Collapse
Affiliation(s)
- Caroline Desgranges
- Department of Chemistry , New York University , New York , New York 10003 , United States
- Department of Chemistry , University of North Dakota , Grand Forks , North Dakota 58202 , United States
| | - Jerome Delhommelle
- Department of Chemistry , New York University , New York , New York 10003 , United States
- Department of Chemistry , University of North Dakota , Grand Forks , North Dakota 58202 , United States
| |
Collapse
|
8
|
Desgranges C, Delhommelle J. Can Ordered Precursors Promote the Nucleation of Solid Solutions? PHYSICAL REVIEW LETTERS 2019; 123:195701. [PMID: 31765200 DOI: 10.1103/physrevlett.123.195701] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/02/2019] [Indexed: 06/10/2023]
Abstract
Crystallization often proceeds through successive stages that lead to a gradual increase in organization. Using molecular simulation, we determine the nucleation pathway for solid solutions of copper and gold. We identify a new nucleation mechanism (liquid→L1_{2} precursor→solid solution) involving a chemically ordered intermediate that is more organized than the end product. This nucleation pathway arises from the low formation energy of L1_{2} clusters which, in turn, promote crystal nucleation. We also show that this mechanism is composition dependent since the high formation energy of other ordered phases precludes them from acting as precursors.
Collapse
Affiliation(s)
- Caroline Desgranges
- Department of Chemistry, New York University, New York, New York 10003, USA and Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| | - Jerome Delhommelle
- Department of Chemistry, New York University, New York, New York 10003, USA and Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| |
Collapse
|
9
|
Desgranges C, Delhommelle J. Non-monotonic variations of the nucleation free energy in a glass-forming ultra-soft particles fluid. SOFT MATTER 2018; 14:5977-5985. [PMID: 29911716 DOI: 10.1039/c8sm00887f] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Using molecular dynamics simulation, we study the impact of the degree of supercooling on the crystal nucleation of ultra-soft particles, modeled with the Gaussian core potential. Focusing on systems with a high number density, our simulations reveal dramatically different behaviors as the degree of supercooling is varied. In the moderate supercooling regime, crystal nucleation proceeds as expected from classical nucleation theory, with a decrease in the free energy of nucleation, as well as in the size of the critical nucleus, as supercooling is increased. On the other hand, in the large supercooling regime, we observe an unusual reversal of behavior with an increase in the free energy of nucleation and in the critical size, as supercooling is increased. This unexpected result is analyzed in terms of the interplay between the glass transition and the crystal nucleation process. Specifically, medium range order crystal-like domains, with structural features different from that of the crystal nucleus, are found to form throughout the system when the supercooling is very large. These, in turn, play a pivotal role in the increase in the free energy of nucleation, as well as in the critical size, as the temperature gets closer to the glass transition.
Collapse
Affiliation(s)
- Caroline Desgranges
- Department of Chemistry, University of North Dakota, Grand Forks, ND 58202, USA.
| | | |
Collapse
|
10
|
Desgranges C, Delhommelle J. Modeling antigen-antibody nanoparticle bioconjugates and their polymorphs. J Chem Phys 2018; 148:124507. [PMID: 29604830 DOI: 10.1063/1.5018855] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The integration of nanomaterials with biomolecules has recently led to the development of new ways of designing biosensors, and through their assembly, to new hybrid structures for novel and exciting applications. In this work, we develop a coarse-grained model for nanoparticles grafted with antibody molecules and their binding with antigens. In particular, we isolate two possible states for antigen-antibody pairs during the binding process, termed as recognition and anchoring states. Using molecular simulation, we calculate the thermodynamic and structural features of three possible crystal structures or polymorphs, the body-centered cubic, simple cubic, and face-centered cubic phases, and of the melt. This leads us to determine the domain of stability of the three solid phases. In particular, the role played by the switching process between anchoring and recognition states during melting is identified, shedding light on the complex microscopic mechanisms in these systems.
Collapse
Affiliation(s)
- Caroline Desgranges
- Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| | - Jerome Delhommelle
- Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| |
Collapse
|
11
|
Desgranges C, Delhommelle J. Unusual Crystallization Behavior Close to the Glass Transition. PHYSICAL REVIEW LETTERS 2018; 120:115701. [PMID: 29601768 DOI: 10.1103/physrevlett.120.115701] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/05/2017] [Revised: 02/13/2018] [Indexed: 06/08/2023]
Abstract
Using molecular simulations, we shed light on the mechanism underlying crystal nucleation in metal alloys and unravel the interplay between crystal nucleation and glass transition, as the conditions of crystallization lie close to this transition. While decreasing the temperature of crystallization usually results in a lower free energy barrier, we find an unexpected reversal of behavior for glass-forming alloys as the temperature of crystallization approaches the glass transition. For this purpose, we simulate the crystallization process in two glass-forming Copper alloys, Ag_{6}Cu_{4}, which has a positive heat of mixing, and CuZr, characterized by a large negative heat of mixing. Our results allow us to identify this unusual behavior as directly correlated with a nonmonotonic temperature dependence for the formation energy of connected icosahedral structures, which are incompatible with crystalline order and impede the development of the crystal nucleus, leading to an unexpectedly larger free energy barrier at low temperature. This, in turn, promotes the formation of a predominantly closed-packed critical nucleus, with fewer defects, thereby suggesting a new way to control the structure of the crystal nucleus, which is of key importance in catalysis.
Collapse
Affiliation(s)
- Caroline Desgranges
- Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| | - Jerome Delhommelle
- Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| |
Collapse
|
12
|
Bechelli S, Gonzalez B, Piquet V, Essafri I, Desgranges C, Delhommelle J. Free Energy of Nucleation and Interplay between Size and Composition in CuNi Systems. J Phys Chem B 2017; 121:8558-8563. [DOI: 10.1021/acs.jpcb.7b06028] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Solene Bechelli
- Department of Chemistry, University of North Dakota, Grand
Forks, North Dakota 58202, United States
| | - B. Gonzalez
- Department of Chemistry, University of North Dakota, Grand
Forks, North Dakota 58202, United States
| | - Vincent Piquet
- Department of Chemistry, University of North Dakota, Grand
Forks, North Dakota 58202, United States
| | - Ilham Essafri
- Department of Chemistry, University of North Dakota, Grand
Forks, North Dakota 58202, United States
| | - Caroline Desgranges
- Department of Chemistry, University of North Dakota, Grand
Forks, North Dakota 58202, United States
| | - Jerome Delhommelle
- Department of Chemistry, University of North Dakota, Grand
Forks, North Dakota 58202, United States
| |
Collapse
|
13
|
Desgranges C, Delhommelle J. Free energy calculations along entropic pathways. III. Nucleation of capillary bridges and bubbles. J Chem Phys 2017. [DOI: 10.1063/1.4982943] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Caroline Desgranges
- Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| | - Jerome Delhommelle
- Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| |
Collapse
|
14
|
Falk ML, Weihs TP. Suppression of homogeneous crystal nucleation of the NiAl intermetallic by a composition gradient: A molecular dynamics study. J Chem Phys 2017. [DOI: 10.1063/1.4982821] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- Michael L. Falk
- Department of Materials Science and Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA
- Department of Mechanical Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA
- Department of Physics and Astronomy, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Timothy P. Weihs
- Department of Materials Science and Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA
- Department of Mechanical Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA
| |
Collapse
|
15
|
Radu M, Kremer K. Enhanced Crystal Growth in Binary Lennard-Jones Mixtures. PHYSICAL REVIEW LETTERS 2017; 118:055702. [PMID: 28211741 DOI: 10.1103/physrevlett.118.055702] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/21/2016] [Indexed: 06/06/2023]
Abstract
We study the crystal growth in binary Lennard-Jones mixtures by molecular dynamics simulations. Growth dynamics, the structure of the liquid-solid interfaces as well as droplet incorporation into the crystal vary with solution properties. For demixed systems we observe a strongly enhanced crystal growth at the cost of enclosed impurities. Furthermore, we find different interface morphologies depending on solubility. We relate our observations to growth mechanisms based on the Gibbs-Thomson effect as well as to predictions of the Kardar-Parisi-Zhang theory in 2+1 dimensions.
Collapse
Affiliation(s)
- M Radu
- Max-Planck Institut für Polymerforschung, Ackermannweg 10, D-55128 Mainz, Germany
| | - K Kremer
- Max-Planck Institut für Polymerforschung, Ackermannweg 10, D-55128 Mainz, Germany
| |
Collapse
|
16
|
Tane M, Kimizuka H, Ichitsubo T. Two distinct crystallization processes in supercooled liquid. J Chem Phys 2017; 144:194505. [PMID: 27208956 DOI: 10.1063/1.4949329] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
Using molecular dynamics simulations we show that two distinct crystallization processes, depending on the temperature at which crystallization occurs, appear in a supercooled liquid. As a model for glass-forming materials, an Al2O3 model system, in which both the glass transition and crystallization from the supercooled liquid can be well reproduced, is employed. Simulations in the framework of an isothermal-isobaric ensemble indicate that the calculated time-temperature-transformation curve for the crystallization to γ(defect spinel)-Al2O3 exhibited a typical nose shape, as experimentally observed in various glass materials. During annealing above the nose temperature, the structure of the supercooled liquid does not change before the crystallization, because of the high atomic mobility (material transport). Thus, the crystallization is governed by the abrupt crystal nucleation, which results in the formation of a stable crystal structure. In contrast, during annealing below the nose temperature, the structure of the supercooled liquid gradually changes before the crystallization, and the formed crystal structure is less stable than that formed above the nose temperature, because of the restricted material transport.
Collapse
Affiliation(s)
- Masakazu Tane
- The Institute of Scientific and Industrial Research, Osaka University, Ibaraki, Osaka 567-0047, Japan
| | - Hajime Kimizuka
- Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Tetsu Ichitsubo
- Department of Materials Science and Engineering, Kyoto University, Kyoto 606-8501, Japan
| |
Collapse
|
17
|
Desgranges C, Delhommelle J. Free energy calculations along entropic pathways. II. Droplet nucleation in binary mixtures. J Chem Phys 2016; 145:234505. [DOI: 10.1063/1.4972011] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Caroline Desgranges
- Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| | - Jerome Delhommelle
- Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| |
Collapse
|
18
|
Desgranges C, Delhommelle J. Free energy calculations along entropic pathways. I. Homogeneous vapor-liquid nucleation for atomic and molecular systems. J Chem Phys 2016; 145:204112. [DOI: 10.1063/1.4968231] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Caroline Desgranges
- Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| | - Jerome Delhommelle
- Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| |
Collapse
|
19
|
Nandi UK, Banerjee A, Chakrabarty S, Bhattacharyya SM. Composition dependence of the glass forming ability in binary mixtures: The role of demixing entropy. J Chem Phys 2016; 145:034503. [PMID: 27448892 DOI: 10.1063/1.4958630] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We present a comparative study of the glass forming ability of binary systems with varying composition, where the systems have similar global crystalline structure (CsCl+fcc). Biased Monte Carlo simulations using umbrella sampling technique show that the free energy cost to create a CsCl nucleus increases as the composition of the smaller particles is decreased. We find that systems with comparatively lower free energy cost to form CsCl nucleus exhibit more pronounced pre-crystalline demixing near the liquid/crystal interface. The structural frustration between the CsCl and fcc crystal demands this demixing. We show that closer to the equimolar mixture, the entropic penalty for demixing is lower and a glass forming system may crystallize when seeded with a nucleus. This entropic penalty as a function of composition shows a non-monotonic behaviour with a maximum at a composition similar to the well known Kob-Anderson (KA) model. Although the KA model shows the maximum entropic penalty and thus maximum frustration against CsCl formation, it also shows a strong tendency towards crystallization into fcc lattice of the larger "A" particles which can be explained from the study of the energetics. Thus for systems closer to the equimolar mixture although it is the requirement of demixing which provides their stability against crystallization, for KA model it is not demixing but slow dynamics and the presence of the "B" particles make it a good glass former. The locally favoured structure around "B" particles is quite similar to the CsCl structure and the incompatibility of CsCl and fcc hinders the fcc structure growth in the KA model. Although the glass forming binary systems studied here are quite similar, differing only in composition, we find that their glass forming ability cannot be attributed to a single phenomenon.
Collapse
Affiliation(s)
- Ujjwal Kumar Nandi
- Polymer Science and Engineering Division, CSIR-National Chemical Laboratory, Pune-411008, India
| | - Atreyee Banerjee
- Polymer Science and Engineering Division, CSIR-National Chemical Laboratory, Pune-411008, India
| | - Suman Chakrabarty
- Physical and Materials Chemistry Division, CSIR-National Chemical Laboratory, Pune-411008, India
| | | |
Collapse
|