1
|
Lauricella M, Ghaani MR, Nandi PK, Meloni S, Kvamme B, English NJ. Double Life of Methanol: Experimental Studies and Nonequilibrium Molecular-Dynamics Simulation of Methanol Effects on Methane-Hydrate Nucleation. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2022; 126:6075-6081. [PMID: 35422892 PMCID: PMC8996238 DOI: 10.1021/acs.jpcc.2c00329] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/27/2022] [Revised: 03/09/2022] [Indexed: 06/14/2023]
Abstract
We have investigated systematically and statistically methanol-concentration effects on methane-hydrate nucleation using both experiment and restrained molecular-dynamics simulation, employing simple observables to achieve an initially homogeneous methane-supersaturated solution particularly favorable for nucleation realization in reasonable simulation times. We observe the pronounced "bifurcated" character of the nucleation rate upon methanol concentration in both experiments and simulation, with promotion at low concentrations and switching to industrially familiar inhibition at higher concentrations. Higher methanol concentrations suppress hydrate growth by in-lattice methanol incorporation, resulting in the formation of "defects", increasing the energy of the nucleus. At low concentrations, on the contrary, the detrimental effect of defects is more than compensated for by the beneficial contribution of CH3 in easing methane incorporation in the cages or replacing it altogether.
Collapse
Affiliation(s)
- Marco Lauricella
- School
of Physics, University College Dublin, Belfield, Dublin 4 D04
V1W8, Ireland
- Istituto
per le Applicazioni del Calcolo, Consiglio
Nazionale delle Ricerche, 00185 Rome, Italy
| | - Mohammad Reza Ghaani
- School
of Chemical and Bioprocess Engineering, University College Dublin, Belfield, Dublin 4, Ireland
| | - Prithwish K. Nandi
- School
of Chemical and Bioprocess Engineering, University College Dublin, Belfield, Dublin 4, Ireland
| | - Simone Meloni
- School
of Physics, University College Dublin, Belfield, Dublin 4 D04
V1W8, Ireland
- Dipartimento
di Scienze Chimiche, Farmaceutiche e Agrarie (DOCPAS), University of Ferrara, 44121 Ferrara, Italy
| | - Bjorn Kvamme
- Hyzen
Energy, Laguna Hills, California 92656, United States
| | - Niall J. English
- School
of Chemical and Bioprocess Engineering, University College Dublin, Belfield, Dublin 4, Ireland
| |
Collapse
|
2
|
Yong X, Burnham CJ, English NJ, Tse JS. Classical and path-integral molecular-dynamics study on liquid water and ice melting using non-empirical TTM2.1-F model. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1652774] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Xue Yong
- Department of Physics and Engineering Physics, University of Saskatchewan, Saskatoon, Canada
| | - Christian J. Burnham
- School of Chemical and Bioprocess Engineering, University College Dublin, Dublin, Ireland
| | - Niall J. English
- School of Chemical and Bioprocess Engineering, University College Dublin, Dublin, Ireland
| | - John S. Tse
- Department of Physics and Engineering Physics, University of Saskatchewan, Saskatoon, Canada
| |
Collapse
|
3
|
Guo Q, Ghaani MR, Nandi PK, English NJ. Pressure-Induced Densification of Ice I h under Triaxial Mechanical Compression: Dissociation versus Retention of Crystallinity for Intermediate States in Atomistic and Coarse-Grained Water Models. J Phys Chem Lett 2018; 9:5267-5274. [PMID: 30145899 DOI: 10.1021/acs.jpclett.8b02270] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Molecular-dynamics (MD) simulation of triaxially pressurized ice Ih up to 30 kbar at 240 K (with sudden mechanical pressurization from its ambient-pressure structure) has been carried out with both the single-particle mW and atomistic TIP4P-Ice water potentials on systems of up to ∼1 million molecules, for times of the order of 100 ns. It was found that the TIP4P-Ice systems adopted a high-density liquid state above ∼7 kbar, while densification of the mW systems retained essentially crystalline order, owing to a failure for the tetrahedral network to break down appreciably from its ice Ih lattice structure. Both are intermediate states adopted along the path toward respective thermodynamically stable states (and with pressure removal show reversion to Ih for mW and to supercooled liquid for TIP4P-Ice), similar to recent ice electro-freezing simulations in "No Man's Land". Densification kinetics showed faster mW-system adaptation.
Collapse
Affiliation(s)
- Qiang Guo
- Collaborative Innovation Center of Chemical Science and Engineering (Tianjin), Tianjin Key Laboratory of Applied Catalysis Science and Technology, College of Chemical Engineering and Technology , Tianjin University , Tianjin 300072 , P.R. China
- School of Chemical and Bioprocess Engineering , University College Dublin , Belfield , Dublin 4 , Ireland
| | - Mohammad Reza Ghaani
- School of Chemical and Bioprocess Engineering , University College Dublin , Belfield , Dublin 4 , Ireland
| | - Prithwish K Nandi
- School of Chemical and Bioprocess Engineering , University College Dublin , Belfield , Dublin 4 , Ireland
- Irish Centre for High-End Computing , Grand Canal Quay , Dublin 2 , Ireland
| | - Niall J English
- School of Chemical and Bioprocess Engineering , University College Dublin , Belfield , Dublin 4 , Ireland
| |
Collapse
|
4
|
Egorov AV, Brodskaya EN, Laaksonen A. The Effect of Single-Atomic Ions on the Melting of Microscopic Ice Particles According to Molecular Dynamics Data. COLLOID JOURNAL 2018. [DOI: 10.1134/s1061933x1805006x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
5
|
Nada H. Anisotropy in geometrically rough structure of ice prismatic plane interface during growth: Development of a modified six-site model of H 2O and a molecular dynamics simulation. J Chem Phys 2017; 145:244706. [PMID: 28049310 DOI: 10.1063/1.4973000] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
This paper presents a modified version of the six-site model of H2O [H. Nada and J. P. J. M. van der Eerden, J. Chem. Phys. 118, 7401 (2003)]. Although the original six-site model was optimized by assuming the cut-off of the Coulomb interaction at an intermolecular distance of 10 Å, the modified model is optimized by using the Ewald method for estimating the Coulomb interaction. Molecular dynamics (MD) simulations of an ice-water interface suggest that the melting point of ice at 1 atm in the modified model is approximately 274.5 K, in good agreement with the real melting point of 273.15 K. MD simulations of bulk ice and water suggest that the modified model reproduces not only the structures and density curves of ice and water, but also the diffusion coefficient of water molecules in water near the melting point at 1 atm. Using the modified model, a large-scale MD simulation of the growth at an ice-water interface of the prismatic plane is performed to elucidate the anisotropy in the interface structure during growth. Simulation results indicate that the geometrical roughness of the ice growth front at the interface is greater in the c-axis direction than in the direction normal to the c-axis when it is analyzed along the axes parallel to the prismatic plane. In addition, during the growth at the interface, the transient appearance of specific crystallographic planes, such as a {202¯1} pyramidal plane, occurs preferentially at the ice growth front. The effect of different ensembles with different simulation systems on the anisotropy in the interface structure is also investigated.
Collapse
Affiliation(s)
- Hiroki Nada
- National Institute of Advanced Industrial Science and Technology (AIST), 16-1 Onogawa, Tsukuba, Ibaraki 305-8569, Japan
| |
Collapse
|
6
|
Ambler M, Vorselaars B, Allen MP, Quigley D. Solid–liquid interfacial free energy of ice Ih, ice Ic, and ice 0 within a mono-atomic model of water via the capillary wave method. J Chem Phys 2017; 146:074701. [DOI: 10.1063/1.4975776] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Affiliation(s)
- Michael Ambler
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Bart Vorselaars
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Michael P. Allen
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
- H. H. Wills Physics Laboratory,
Royal Fort, Tyndall Avenue, Bristol BS8 1TL, United Kingdom
| | - David Quigley
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
7
|
Xu Y, Petrik NG, Smith RS, Kay BD, Kimmel GA. Growth rate of crystalline ice and the diffusivity of supercooled water from 126 to 262 K. Proc Natl Acad Sci U S A 2016; 113:14921-14925. [PMID: 27956609 PMCID: PMC5206540 DOI: 10.1073/pnas.1611395114] [Citation(s) in RCA: 85] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Understanding deeply supercooled water is key to unraveling many of water's anomalous properties. However, developing this understanding has proven difficult due to rapid and uncontrolled crystallization. Using a pulsed-laser-heating technique, we measure the growth rate of crystalline ice, G(T), for 180 K < T < 262 K, that is, deep within water's "no man's land" in ultrahigh-vacuum conditions. Isothermal measurements of G(T) are also made for 126 K ≤ T ≤ 151 K. The self-diffusion of supercooled liquid water, D(T), is obtained from G(T) using the Wilson-Frenkel model of crystal growth. For T > 237 K and P ∼ 10-8 Pa, G(T) and D(T) have super-Arrhenius ("fragile") temperature dependences, but both cross over to Arrhenius ("strong") behavior with a large activation energy in no man's land. The fact that G(T) and D(T) are smoothly varying rules out the hypothesis that liquid water's properties have a singularity at or near 228 K at ambient pressures. However, the results are consistent with a previous prediction for D(T) that assumed no thermodynamic transitions occur in no man's land.
Collapse
Affiliation(s)
- Yuntao Xu
- Chemical Physics & Analysis, Physical Sciences Division, Physical & Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, WA 99352
| | - Nikolay G Petrik
- Chemical Physics & Analysis, Physical Sciences Division, Physical & Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, WA 99352
| | - R Scott Smith
- Chemical Physics & Analysis, Physical Sciences Division, Physical & Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, WA 99352
| | - Bruce D Kay
- Chemical Physics & Analysis, Physical Sciences Division, Physical & Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, WA 99352
| | - Greg A Kimmel
- Chemical Physics & Analysis, Physical Sciences Division, Physical & Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, WA 99352
| |
Collapse
|
8
|
|
9
|
English NJ, Waldron CJ. Perspectives on external electric fields in molecular simulation: progress, prospects and challenges. Phys Chem Chem Phys 2016; 17:12407-40. [PMID: 25903011 DOI: 10.1039/c5cp00629e] [Citation(s) in RCA: 142] [Impact Index Per Article: 17.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
In this review, the application of a wide variety of external electric fields in molecular simulation shall be discussed, including time-varying and electromagnetic, as well as the utility and potential impact and prospects for exploitation of such simulations for real-world and industrial end use. In particular, non-equilibrium molecular dynamics will be discussed, as well as challenges in addressing adequate thermostatting and scaling field amplitudes to more experimentally relevant levels. Attention shall be devoted to recent progress and advances in external fields in ab initio molecular simulation and dynamics, as well as elusive challenges thereof (and, to some extent, for molecular dynamics from empirical potentials), such as timescales required to observe low-frequency and intensity field effects. The challenge of deterministic molecular dynamics in external fields in sampling phase space shall be discussed, along with prospects for application of fields in enhanced-sampling simulations. Finally, the application of external electric fields to a wide variety of aqueous, nanoscale and biological systems will be discussed, often motivated by the possibility of exploitation in real-world applications, which serve to underpin our molecular-level understanding of field effects in terms of microscopic mechanisms, and possibly with a view to control thereof.
Collapse
Affiliation(s)
- Niall J English
- School of Chemical and Bioprocess Engineering, University College Dublin, Belfield, Dublin 4, Ireland.
| | | |
Collapse
|
10
|
Małolepsza E, Keyes T. Pathways through Equilibrated States with Coexisting Phases for Gas Hydrate Formation. J Phys Chem B 2015; 119:15857-65. [DOI: 10.1021/acs.jpcb.5b06832] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Edyta Małolepsza
- Department
of Chemistry, Boston University, Boston, Massachusetts 02215-2521, United States
| | - Tom Keyes
- Department
of Chemistry, Boston University, Boston, Massachusetts 02215-2521, United States
| |
Collapse
|
11
|
English NJ, Tse JS. Massively parallel molecular dynamics simulation of formation of ice-crystallite precursors in supercooled water: incipient-nucleation behavior and role of system size. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:032132. [PMID: 26465451 DOI: 10.1103/physreve.92.032132] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/09/2015] [Indexed: 06/05/2023]
Abstract
Ice-crystallite precursor formation in supercooled water was studied via molecular dynamics for systems ranging from ∼10^{6} to 8.6×10^{6} molecules, using a tetrahedrally biased single-site "mW" model. This has established system-size effects in the early onset of nucleation, so as to study often-transient precursors' beguiling propensity to "flicker" into instantaneous locally ordered molecular arrangements redolent of ice. In addition, the adoption of solidlike and liquidlike bimodal local configurational-energy distributions was observed, characteristic of early nucleation. Larger systems favored a higher probability of precursor formation, although such ones were not usually longer lived relative to those in smaller systems (which themselves are rather transient). It was concluded tentatively that subtle effects of differences in systemwide density fluctuations and accessible lower-frequency modes tend to favor precursor formation in larger systems, although not necessarily the precursor's kinetic stability.
Collapse
Affiliation(s)
- Niall J English
- School of Chemical and Bioprocess Engineering, University College Dublin, Belfield, Dublin 4, Ireland
| | - John S Tse
- Department of Physics and Engineering Physics, University of Saskatchewan, Saskatoon, Saskatchewan, Canada S7N 5E2
| |
Collapse
|
12
|
Lauricella M, Meloni S, Liang S, English NJ, Kusalik PG, Ciccotti G. Clathrate structure-type recognition: Application to hydrate nucleation and crystallisation. J Chem Phys 2015; 142:244503. [DOI: 10.1063/1.4922696] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Affiliation(s)
- Marco Lauricella
- School of Physics, University College Dublin, Belfield, Dublin 4, Ireland
- Instituto per le Applicazioni del Calcolo, Consiglio Nazionale delle Ricerche, Rome, Italy
| | - Simone Meloni
- School of Physics, University College Dublin, Belfield, Dublin 4, Ireland
- Laboratory of Computational Chemistry and Biochemistry, École Polytechnique Fedérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Shuai Liang
- Deptartment of Chemistry, University of Calgary, Calgary, Alberta T2N 1N4, Canada
- Key Laboratory of Gas Hydrate, Guangzhou Institute of Energy Conversion, Chinese Academy of Sciences, Guangzhou, China
| | - Niall J. English
- School of Chemical and Bioprocess Engineering, University College Dublin, Belfield, Dublin 4, Ireland
| | - Peter G. Kusalik
- Deptartment of Chemistry, University of Calgary, Calgary, Alberta T2N 1N4, Canada
| | - Giovanni Ciccotti
- School of Physics, University College Dublin, Belfield, Dublin 4, Ireland
- Dipartimento di Fisica and CNISM, Università La Sapienza, P. le A. Moro 5, 00185 Rome, Italy
| |
Collapse
|