1
|
Bhendale M, Indra A, Singh JK. Does freezing induce self-assembly of polymers? A molecular dynamics study. SOFT MATTER 2023; 19:7570-7579. [PMID: 37751160 DOI: 10.1039/d3sm00892d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/27/2023]
Abstract
This work investigates the freezing-induced self-assembly (FISA) of polyvinyl alcohol (PVA) and PVA-like polymers using molecular dynamics simulations. In particular, the effect of the degree of supercooling, degree of polymerization, polymer type, and initial local concentration on the FISA was studied. It was found that the preeminent factor responsible for FISA is not the diffusion of the polymers away from the nucleating ice front, but the increase in the polymer's local concentration upon freezing of the solvent (water). At a higher degree of supercooling, the polymers are engulfed by the growing ice front, impeding their diffusion into the supercooled solution and finally inhibiting their self-assembly. Conversely, at a relatively lower degree of supercooling, the rate of diffusion of the polymers into the supercooled solution is higher, which increases their local concentration and results in FISA. FISA was also observed to depend on the polymer-solvent interactions. Strongly favorable solute-solvent interactions hinder the self-assembly, whereas unfavorable solute-solvent interactions promote the self-assembly. The polymer and aggregate morphology were investigated using the radius of gyration, end-to-end distance, and asphericity analysis. This study brings molecular insights into the quintessential factors governing self-assembly via freezing of the solvent, which is a novel self-assembly technique especially suitable for biomedical applications.
Collapse
Affiliation(s)
- Mangesh Bhendale
- Department of Chemical Engineering, Indian Institute of Technology, Kanpur, Uttar Pradesh 208016, India.
| | - Aindrila Indra
- Department of Chemical Engineering, Indian Institute of Technology, Kanpur, Uttar Pradesh 208016, India.
| | - Jayant K Singh
- Department of Chemical Engineering, Indian Institute of Technology, Kanpur, Uttar Pradesh 208016, India.
- Prescience Insilico Private Limited, 5th floor, Novel MSR Building, Marathalli, Bengaluru, Karnataka 560037, India
| |
Collapse
|
2
|
Srirangam S, Bhendale M, Singh JK. Does supercooled water retain its universal nucleation behavior under shear at high pressure? Phys Chem Chem Phys 2023; 25:21528-21537. [PMID: 37545252 DOI: 10.1039/d3cp01605f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/08/2023]
Abstract
Understanding the nucleation of homogeneous flow systems at high pressures is vital in protein crystallization and cryopreservation, where high pressure prevents the freezing of biological samples. This study examines the behavior of ice nucleation under shear at various pressures and explores the universal nucleation behavior of the sheared systems applied to supercooled water at higher pressures. In this study, the nucleation rates for the TIP4P/Ice model via a seeding method based on extended classical nucleation theory (CNT) are computed at pressures of 1, 100, 500, 700, and 1000 bar and a constant temperature of 240 K. Using extended CNT with explicitly embodying the shear rate, we analyzed the dependence of pressure on the transport and thermodynamic properties. In line with previous studies, we observed that Δμliq-ice and viscosity decrease while diffusivity increases with an increase in pressure. Furthermore, we showed that the dependence of the nucleation rate on shear at higher pressure is non-monotonic, with the maximum at optimal shear rates between 107 and 108 s-1. Our results demonstrate a non-monotonic pressure dependence of the optimal shear rates, which could originate from a violation of the Stokes-Einstein relation.
Collapse
Affiliation(s)
- Snehitha Srirangam
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Uttar Pradesh, 208016, India.
| | - Mangesh Bhendale
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Uttar Pradesh, 208016, India.
| | - Jayant K Singh
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Uttar Pradesh, 208016, India.
- Prescience Insilico Private Limited, 5th Floor, Novel MSR Building, Marathalli, Bengaluru, Karnataka, 560037, India
| |
Collapse
|
3
|
Lodesani F, Menziani MC, Urata S, Pedone A. Biasing Crystallization in Fused Silica: An Assessment of Optimal Metadynamics Parameters. J Chem Phys 2022; 156:194501. [DOI: 10.1063/5.0089183] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Metadynamics is a useful technique to study rare events such as crystallization. It has been only recently applied to study nucleation and crystallization in glass-forming liquids such as silicates but the optimal set of parameters to drive crystallization and obtaining converged Free Energy Surfaces is still unexplored. <p>In this work, we systematically investigated the effects of the simulation conditions to efficiently study the thermodynamics and mechanism of crystallization in highly viscous systems. As a prototype system, we used fused silica, which easily crystallizes to β-cristobalite through MetaD simulations, owing to its simple microstructure. We investigated the influence of the height, width, and bias factor used to define the biasing Gaussian potential, as well as the effects of the temperature and system size on the results. Among these parameters, the bias factor and temperature seem to be most effective to sample the free energy landscape of melt to crystal transition and reach convergence more quickly. We also demonstrate that the temperature rescaling from T > Tm is a reliable approach to recover free energy surfaces below Tm, provided that the temperature gap is below 600 K and the configurational space has been properly sampled. Finally, albeit a complete crystallization is hard to achieve with large simulation boxes, these can be reliably and effectively exploited to study the first stages of nucleation.
Collapse
Affiliation(s)
- Federica Lodesani
- Universita degli Studi di Modena e Reggio Emilia Dipartimento di Scienze Chimiche e Geologiche, Italy
| | - Maria Cristina Menziani
- Universita degli Studi di Modena e Reggio Emilia Dipartimento di Scienze Chimiche e Geologiche, Italy
| | - Shingo Urata
- Innovative Technology Laboratories, AGC Inc., Japan
| | - Alfonso Pedone
- Universita degli Studi di Modena e Reggio Emilia Dipartimento di Scienze Chimiche e Geologiche, Italy
| |
Collapse
|
4
|
Blow KE, Quigley D, Sosso GC. The seven deadly sins: When computing crystal nucleation rates, the devil is in the details. J Chem Phys 2021; 155:040901. [PMID: 34340373 DOI: 10.1063/5.0055248] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
The formation of crystals has proven to be one of the most challenging phase transformations to quantitatively model-let alone to actually understand-be it by means of the latest experimental technique or the full arsenal of enhanced sampling approaches at our disposal. One of the most crucial quantities involved with the crystallization process is the nucleation rate, a single elusive number that is supposed to quantify the average probability for a nucleus of critical size to occur within a certain volume and time span. A substantial amount of effort has been devoted to attempt a connection between the crystal nucleation rates computed by means of atomistic simulations and their experimentally measured counterparts. Sadly, this endeavor almost invariably fails to some extent, with the venerable classical nucleation theory typically blamed as the main culprit. Here, we review some of the recent advances in the field, focusing on a number of perhaps more subtle details that are sometimes overlooked when computing nucleation rates. We believe it is important for the community to be aware of the full impact of aspects, such as finite size effects and slow dynamics, that often introduce inconspicuous and yet non-negligible sources of uncertainty into our simulations. In fact, it is key to obtain robust and reproducible trends to be leveraged so as to shed new light on the kinetics of a process, that of crystal nucleation, which is involved into countless practical applications, from the formulation of pharmaceutical drugs to the manufacturing of nano-electronic devices.
Collapse
Affiliation(s)
- Katarina E Blow
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - David Quigley
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Gabriele C Sosso
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
5
|
Goswami A, Dalal IS, Singh JK. Universal Nucleation Behavior of Sheared Systems. PHYSICAL REVIEW LETTERS 2021; 126:195702. [PMID: 34047572 DOI: 10.1103/physrevlett.126.195702] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/08/2021] [Accepted: 04/15/2021] [Indexed: 06/12/2023]
Abstract
Using molecular simulations and a modified classical nucleation theory, we study the nucleation, under flow, of a variety of liquids: different water models, Lennard-Jones, and hard sphere colloids. Our approach enables us to analyze a wide range of shear rates inaccessible to brute-force simulations. Our results reveal that the variation of the nucleation rate with shear is universal. A simplified version of the theory successfully captures the nonmonotonic temperature dependence of the nucleation behavior, which is shown to originate from the violation of the Stokes-Einstein relation.
Collapse
Affiliation(s)
- Amrita Goswami
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Kanpur, Uttar Pradesh 208016, India
| | - Indranil Saha Dalal
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Kanpur, Uttar Pradesh 208016, India
| | - Jayant K Singh
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Kanpur, Uttar Pradesh 208016, India
| |
Collapse
|
6
|
Goswami A, Singh JK. Homogeneous nucleation of sheared liquids: advances and insights from simulations and theory. Phys Chem Chem Phys 2021; 23:15402-15419. [PMID: 34279013 DOI: 10.1039/d1cp02617h] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
One of the most ubiquitous and technologically important phenomena in nature is the nucleation of homogeneous flowing systems. The microscopic effects of shear on a nucleating system are still imperfectly understood, although in recent years a consistent picture has emerged. The opposing effects of shear can be split into two major contributions for simple atomic and molecular liquids: increase of the energetic cost of nucleation, and enhancement of the kinetics. In this perspective, we describe the latest computational and theoretical techniques which have been developed over the past two decades. We collate and unify the overarching influences of shear, temperature, and supersaturation on the process of homogeneous nucleation. Experimental techniques and capabilities are discussed, against the backdrop of results from simulations and theory. Although we primarily focus on simple systems, we also touch upon the sheared nucleation of more complex systems, including glasses and polymer melts. We speculate on the promising directions and possible advances that could come to fruition in the future.
Collapse
Affiliation(s)
- Amrita Goswami
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Uttar Pradesh 208016, India.
| | - Jayant K Singh
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Uttar Pradesh 208016, India.
| |
Collapse
|