1
|
Cui X, Yang C, Sun Q, Zhang W, Wang X. Investigating Shear Stress of Ice Accumulated on Surfaces with Various Roughnesses: Effects of a Quasi-Water Layer. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2024. [PMID: 38954504 DOI: 10.1021/acs.langmuir.4c00617] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/04/2024]
Abstract
The investigation of the anti-icing/deicing is essential because the icing phenomenon deteriorates the natural environment and various projects. By conducting molecular dynamics simulation, this work analyzes the effect of the quasi-water layer on the ice shear stress over smooth and rough surfaces, along with the underlying physics of the quasi-water layer. The results indicate that the thickness of the quasi-water layer monotonically increases with temperature, resulting in a monotonic decrease in the ice shear stress on the smooth surface. Due to the joint effects of the smooth surface wettability and the quasi-water layer, the ice shear stress increases and then decreases to almost a constant value when the surface changes from a hydrophobic to a hydrophilic one. For rough surfaces with stripe nanostructures, when the width of the bump for one case equals the depression for the other case, the variations of shear stress with height for these two cases are almost the same. The rough surface is effective in reducing the ice shear stress compared to the smooth surface due to the thickening of the quasi-water layer. Each molecule in the quasi-water layer and its four nearest neighboring molecules gradually form a tetrahedral ice-like structure along the direction away from the surface. The radial distribution function also shows that the quasi-water layer resembles the liquid water rather than the ice structure. These findings shed light on developing anti-icing and deicing techniques.
Collapse
Affiliation(s)
- Xinjiao Cui
- Institute of Thermal Science and Technology, Shandong University, Jinan 250061, China
- Institute for Advanced Technology, Shandong University, Jinan 250061, China
| | - Chao Yang
- Institute of Thermal Science and Technology, Shandong University, Jinan 250061, China
| | - Qiangqiang Sun
- Faculty of Engineering, University of Nottingham, Nottingham NG7 2RD, U.K
| | - Wenqiang Zhang
- School of Mechatronical Engineering, Beijing Institute of Technology, Beijing 100081, China
| | - Xinyu Wang
- Institute of Thermal Science and Technology, Shandong University, Jinan 250061, China
- Shenzhen Research Institute of Shandong University, Shenzhen 518057, China
| |
Collapse
|
2
|
Lupi L, Gallo P. Mode coupling behavior and fragile to strong transition of trehalose in a binary mixture with water upon supercooling. J Chem Phys 2024; 160:244501. [PMID: 38912627 DOI: 10.1063/5.0218369] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2024] [Accepted: 06/06/2024] [Indexed: 06/25/2024] Open
Abstract
We perform molecular dynamics simulations of a binary mixture of water and trehalose with the TIP4P/Ice water model. We analyze the slow dynamics of trehalose molecules in the mildly supercooled region for concentrations of 3.66 and 18.57 wt. %. We previously studied the dynamics of water in the same mixtures. Supercooled TIP4P/Ice water solvating trehalose molecules was found to follow the Mode Coupling Theory (MCT) and to undergo a transition from a fragile to a strong behavior for both concentrations. Here, we show that also the dynamics of trehalose molecules follows the MCT and displays a fragile to strong crossover (FSC). The results show that trehalose in binary mixtures with water shares with it the dynamical behavior typical of glass forming liquids. Moreover, the FSC for trehalose structural relaxation times is found to occur at temperatures close to those previously obtained for water in the same solutions, showing that the dynamics of the solute is strongly coupled to that of the solvent. We also perform a MCT test showing that the trehalose dynamics obeys the MCT time-temperature superposition principle and that the exponents derived from the theory and the ones obtained from fitting procedure of the relaxation times are comparable, confirming that trehalose molecules in supercooled water solutions follow the MCT of glassy dynamics. Moreover, as predicted by the theory, trehalose particles have MCT parameters comparable to those of water in the same mixtures. This is an important result, given that MCT was originally formulated for monoatomic particles.
Collapse
Affiliation(s)
- Laura Lupi
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Rome, Italy
| | - Paola Gallo
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Rome, Italy
| |
Collapse
|
3
|
Tavagnacco L, Zanatta M, Buratti E, Bertoldo M, Chiessi E, Appel M, Natali F, Orecchini A, Zaccarelli E. Water slowing down drives the occurrence of the low temperature dynamical transition in microgels. Chem Sci 2024; 15:9249-9257. [PMID: 38903230 PMCID: PMC11186305 DOI: 10.1039/d4sc02650k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2024] [Accepted: 05/05/2024] [Indexed: 06/22/2024] Open
Abstract
The protein dynamical transition marks an increase in atomic mobility and the onset of anharmonic motions at a critical temperature (T d), which is considered relevant for protein functionality. This phenomenon is ubiquitous, regardless of protein composition, structure and biological function and typically occurs at large protein content, to avoid water crystallization. Recently, a dynamical transition has also been reported in non-biological macromolecules, such as poly(N-isopropyl acrylamide) (PNIPAM) microgels, bearing many similarities to proteins. While the generality of this phenomenon is well-established, the role of water in the transition remains a subject of debate. In this study, we use atomistic molecular dynamics (MD) simulations and elastic incoherent neutron scattering (EINS) experiments with selective deuteration to investigate the microscopic origin of the dynamical transition and distinguish water and PNIPAM roles. While a standard analysis of EINS experiments would suggest that the dynamical transition occurs in PNIPAM and water at a similar temperature, simulations reveal a different perspective, also qualitatively supported by experiments. From room temperature down to about 180 K, PNIPAM exhibits only modest changes of dynamics, while water, being mainly hydration water under the probed extreme confinement, significantly slows down and undergoes a mode-coupling transition from diffusive to activated. Our findings therefore challenge the traditional view of the dynamical transition, demonstrating that it occurs in proximity of the water mode-coupling transition, shedding light on the intricate interplay between polymer and water dynamics.
Collapse
Affiliation(s)
- Letizia Tavagnacco
- CNR Institute of Complex Systems, Uos Sapienza Piazzale Aldo Moro 2 00185 Rome Italy
- Department of Physics, Sapienza University of Rome Piazzale Aldo Moro 2 00185 Rome Italy
| | - Marco Zanatta
- Department of Physics, University of Trento Via Sommarive 14 38123 Trento Italy
| | - Elena Buratti
- Department of Chemical, Pharmaceutical and Agricultural Sciences, University of Ferrara Via L. Borsari 46 44121 Ferrara Italy
| | - Monica Bertoldo
- Department of Chemical, Pharmaceutical and Agricultural Sciences, University of Ferrara Via L. Borsari 46 44121 Ferrara Italy
| | - Ester Chiessi
- Department of Chemical Science and Technologies, University of Rome Tor Vergata Via della Ricerca Scientifica I 00133 Rome Italy
| | - Markus Appel
- Institut Laue-Langevin 71 avenue des Martyrs, CS 20156 38042 Grenoble Cedex 9 France
| | - Francesca Natali
- CNR-IOM, Operative Group Grenoble (OGG), Institut Laue Langevin F-38042 Grenoble France
| | - Andrea Orecchini
- Dipartimento di Fisica e Geologia, Università di Perugia Via Alessandro Pascoli 06123 Perugia Italy
- CNR-IOM c/o Dipartimento di Fisica e Geologia, Università di Perugia Via Alessandro Pascoli 06123 Perugia Italy
| | - Emanuela Zaccarelli
- CNR Institute of Complex Systems, Uos Sapienza Piazzale Aldo Moro 2 00185 Rome Italy
- Department of Physics, Sapienza University of Rome Piazzale Aldo Moro 2 00185 Rome Italy
| |
Collapse
|
4
|
Liu L, Guan D, Lu Y, Sun M, Liu Y, Zhao J, Yang L. A Molecular Dynamics Study on Xe/Kr Separation Mechanisms Using Crystal Growth Method. ACS OMEGA 2024; 9:25822-25831. [PMID: 38911791 PMCID: PMC11191100 DOI: 10.1021/acsomega.4c00108] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/04/2024] [Revised: 03/31/2024] [Accepted: 05/23/2024] [Indexed: 06/25/2024]
Abstract
The separation of xenon/krypton gas mixtures is a valuable but challenging endeavor in the gas industry due to their similar physical characteristics and closely sized molecules. To address this, we investigated the effectiveness of the hydrate-based gas separation method for mixed Xe-Kr gas via molecular dynamics (MD) simulations. The formation process of hydrates facilitates the encapsulation of guest molecules within hydrate cages, offering a potential strategy for gas separation. Higher temperatures and pressures are advantageous for accelerating the hydrate growth rate. The final occupancy of guest molecules and empty cages within 512, 51264, and all hydrate cages were thoroughly examined. An increase in the pressure and temperature enhanced the occupancy rates of Xe in both 512 and 51264 cages, whereas elevated pressure alone improved the occupancy of Kr in 51264 cages. However, the impact of temperature and pressure on Kr occupancy within 512 cages was found to be minimal. Elevated temperature and pressure resulted in a reduced occupancy of empty cages. Predominantly, 51264 cages were occupied by Xe, whereas Kr showed a propensity to occupy the 512 cages. With increasing simulated pressure, the final occupancy of Xe molecules in all cages rose from 0.37 to 0.41 for simulations at 260 K, while the final occupancy of empty cages decreased from 0.24 to 0.2.
Collapse
Affiliation(s)
- Liangliang Liu
- Shenyang
Aircraft Design Institute Shenyang 110042, China
| | - Dawei Guan
- Key
Laboratory of Ocean Energy Utilization and Energy Conservation of
Ministry of Education, Dalian University
of Technology, Dalian 116024, China
| | - Yi Lu
- Shenyang
Aircraft Design Institute Shenyang 110042, China
| | - Mingrui Sun
- Key
Laboratory of Ocean Energy Utilization and Energy Conservation of
Ministry of Education, Dalian University
of Technology, Dalian 116024, China
| | - Yu Liu
- Key
Laboratory of Ocean Energy Utilization and Energy Conservation of
Ministry of Education, Dalian University
of Technology, Dalian 116024, China
| | - Jiafei Zhao
- Key
Laboratory of Ocean Energy Utilization and Energy Conservation of
Ministry of Education, Dalian University
of Technology, Dalian 116024, China
| | - Lei Yang
- Key
Laboratory of Ocean Energy Utilization and Energy Conservation of
Ministry of Education, Dalian University
of Technology, Dalian 116024, China
| |
Collapse
|
5
|
Miles CM, Cullen S, Kenaan H, Gu W, Andrews GP, Sosso GC, Tian Y. Unravelling the interactions between small molecules and liposomal bilayers via molecular dynamics and thermodynamic modelling. Int J Pharm 2024; 660:124367. [PMID: 38901537 DOI: 10.1016/j.ijpharm.2024.124367] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2024] [Revised: 06/06/2024] [Accepted: 06/17/2024] [Indexed: 06/22/2024]
Abstract
Lipid-based drug delivery systems hold immense promise in addressing critical medical needs, from cancer and neurodegenerative diseases to infectious diseases. By encapsulating active pharmaceutical ingredients - ranging from small molecule drugs to proteins and nucleic acids - these nanocarriers enhance treatment efficacy and safety. However, their commercial success faces hurdles, such as the lack of a systematic design approach and the issues related to scalability and reproducibility. This work aims to provide insights into the drug-phospholipid interaction by combining molecular dynamic simulations and thermodynamic modelling techniques. In particular, we have made a connection between the structural properties of the drug-phospholipid system and the physicochemical performance of the drug-loaded liposomal nanoformulations. We have considered two prototypical drugs, felodipine (FEL) and naproxen (NPX), and one model hydrogenated soy phosphatidylcholine (HSPC) bilayer membrane. Molecular dynamic simulations revealed which regions within the phospholipid bilayers are most and least favoured by the drug molecules. NPX tends to reside at the water-phospholipid interface and is characterized by a lower free energy barrier for bilayer membrane permeation. Meanwhile, FEL prefers to sit within the hydrophobic tails of the phospholipids and is characterized by a higher free energy barrier for membrane permeation. Flory-Huggins thermodynamic modelling, small angle X-ray scattering, dynamic light scattering, TEM, and drug release studies of these liposomal nanoformulations confirmed this drug-phospholipid structural difference. The naproxen-phospholipid system has a lower free energy barrier for permeation, higher drug miscibility with the bilayer, larger liposomal nanoparticle size, and faster drug release in the aqueous medium than felodipine. We suggest that this combination of molecular dynamics and thermodynamics approach may offer a new tool for designing and developing lipid-based nanocarriers for unmet medical applications.
Collapse
Affiliation(s)
- Christopher M Miles
- Department of Chemistry, University of Warwick, Coventry, CV4 7AL, United Kingdom
| | - Shane Cullen
- School of Pharmacy, Queen's University Belfast, 97 Lisburn Road, Belfast BT9 7BL, United Kingdom
| | - Hussein Kenaan
- School of Pharmacy, Queen's University Belfast, 97 Lisburn Road, Belfast BT9 7BL, United Kingdom
| | - Wenjie Gu
- School of Pharmacy, Queen's University Belfast, 97 Lisburn Road, Belfast BT9 7BL, United Kingdom
| | - Gavin P Andrews
- School of Pharmacy, Queen's University Belfast, 97 Lisburn Road, Belfast BT9 7BL, United Kingdom
| | - Gabriele C Sosso
- Department of Chemistry, University of Warwick, Coventry, CV4 7AL, United Kingdom.
| | - Yiwei Tian
- School of Pharmacy, Queen's University Belfast, 97 Lisburn Road, Belfast BT9 7BL, United Kingdom.
| |
Collapse
|
6
|
Wu D, Zhao Z, Lin B, Song Y, Qi J, Jiang J, Yuan Z, Cheng B, Zhao M, Tian Y, Wang Z, Wu M, Bian K, Liu KH, Xu LM, Zeng XC, Wang EG, Jiang Y. Probing structural superlubricity of two-dimensional water transport with atomic resolution. Science 2024; 384:1254-1259. [PMID: 38870285 DOI: 10.1126/science.ado1544] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2024] [Accepted: 05/01/2024] [Indexed: 06/15/2024]
Abstract
Low-dimensional water transport can be drastically enhanced under atomic-scale confinement. However, its microscopic origin is still under debate. In this work, we directly imaged the atomic structure and transport of two-dimensional water islands on graphene and hexagonal boron nitride surfaces using qPlus-based atomic force microscopy. The lattice of the water island was incommensurate with the graphene surface but commensurate with the boron nitride surface owing to different surface electrostatics. The area-normalized static friction on the graphene diminished as the island area was increased by a power of ~-0.58, suggesting superlubricity behavior. By contrast, the friction on the boron nitride appeared insensitive to the area. Molecular dynamic simulations further showed that the friction coefficient of the water islands on the graphene could reduce to <0.01.
Collapse
Affiliation(s)
- Da Wu
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, China
| | - Zhengpu Zhao
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, China
| | - Bo Lin
- Department of Materials Science and Engineering, City University of Hong Kong, Hong Kong 999077, China
| | - Yizhi Song
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, China
| | - Jiajie Qi
- State Key Laboratory for Mesoscopic Physics, Frontiers Science Center for Nano-optoelectronics, School of Physics, Peking University, Beijing 100871, China
| | - Jian Jiang
- Department of Materials Science and Engineering, City University of Hong Kong, Hong Kong 999077, China
| | - Zifeng Yuan
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, China
| | - Bowei Cheng
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, China
| | - Mengze Zhao
- State Key Laboratory for Mesoscopic Physics, Frontiers Science Center for Nano-optoelectronics, School of Physics, Peking University, Beijing 100871, China
| | - Ye Tian
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, China
| | - Zhichang Wang
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, China
| | - Muhong Wu
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, China
- Interdisciplinary Institute of Light-Element Quantum Materials and Research Centre for Light-Element Advanced Materials, Peking University, Beijing 100871, China
- Songshan Lake Materials Laboratory, Institute of Physics, CAS and School of Physics, Liaoning University, Shenyang 110036, China
| | - Ke Bian
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, China
| | - Kai-Hui Liu
- State Key Laboratory for Mesoscopic Physics, Frontiers Science Center for Nano-optoelectronics, School of Physics, Peking University, Beijing 100871, China
- Interdisciplinary Institute of Light-Element Quantum Materials and Research Centre for Light-Element Advanced Materials, Peking University, Beijing 100871, China
- Collaborative Innovation Center of Quantum Matter, Beijing 100871, China
| | - Li-Mei Xu
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, China
- Interdisciplinary Institute of Light-Element Quantum Materials and Research Centre for Light-Element Advanced Materials, Peking University, Beijing 100871, China
- Collaborative Innovation Center of Quantum Matter, Beijing 100871, China
| | - Xiao Cheng Zeng
- Department of Materials Science and Engineering, City University of Hong Kong, Hong Kong 999077, China
- Hong Kong Institute for Clean Energy, City University of Hong Kong, Hong Kong 999077, China
| | - En-Ge Wang
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, China
- Interdisciplinary Institute of Light-Element Quantum Materials and Research Centre for Light-Element Advanced Materials, Peking University, Beijing 100871, China
- Songshan Lake Materials Laboratory, Institute of Physics, CAS and School of Physics, Liaoning University, Shenyang 110036, China
- Tsientang Institute for Advanced Study, Zhejiang 310024, China
| | - Ying Jiang
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, China
- Interdisciplinary Institute of Light-Element Quantum Materials and Research Centre for Light-Element Advanced Materials, Peking University, Beijing 100871, China
- Collaborative Innovation Center of Quantum Matter, Beijing 100871, China
- New Cornerstone Science Laboratory, Peking University, Beijing 100871, China
| |
Collapse
|
7
|
Thosar AU, Cai Y, Marks SM, Vicars Z, Choi J, Pallath A, Patel AJ. On the engulfment of antifreeze proteins by ice. Proc Natl Acad Sci U S A 2024; 121:e2320205121. [PMID: 38833468 PMCID: PMC11181090 DOI: 10.1073/pnas.2320205121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2023] [Accepted: 04/16/2024] [Indexed: 06/06/2024] Open
Abstract
Antifreeze proteins (AFPs) are remarkable biomolecules that suppress ice formation at trace concentrations. To inhibit ice growth, AFPs must not only bind to ice crystals, but also resist engulfment by ice. The highest supercooling, [Formula: see text], for which AFPs are able to resist engulfment is widely believed to scale as the inverse of the separation, [Formula: see text], between bound AFPs, whereas its dependence on the molecular characteristics of the AFP remains poorly understood. By using specialized molecular simulations and interfacial thermodynamics, here, we show that in contrast with conventional wisdom, [Formula: see text] scales as [Formula: see text] and not as [Formula: see text]. We further show that [Formula: see text] is proportional to AFP size and that diverse naturally occurring AFPs are optimal at resisting engulfment by ice. By facilitating the development of AFP structure-function relationships, we hope that our findings will pave the way for the rational design of AFPs.
Collapse
Affiliation(s)
- Aniket U. Thosar
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, PA19104
| | - Yusheng Cai
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, PA19104
| | - Sean M. Marks
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, PA19104
| | - Zachariah Vicars
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, PA19104
| | - Jeongmoon Choi
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, PA19104
| | - Akash Pallath
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, PA19104
| | - Amish J. Patel
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, PA19104
| |
Collapse
|
8
|
Hong J, Tian Y, Liang T, Liu X, Song Y, Guan D, Yan Z, Guo J, Tang B, Cao D, Guo J, Chen J, Pan D, Xu LM, Wang EG, Jiang Y. Imaging surface structure and premelting of ice Ih with atomic resolution. Nature 2024; 630:375-380. [PMID: 38778112 DOI: 10.1038/s41586-024-07427-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2023] [Accepted: 04/16/2024] [Indexed: 05/25/2024]
Abstract
Ice surfaces are closely relevant to many physical and chemical properties, such as melting, freezing, friction, gas uptake and atmospheric reaction1-8. Despite extensive experimental and theoretical investigations9-17, the exact atomic structures of ice interfaces remain elusive owing to the vulnerable hydrogen-bonding network and the complicated premelting process. Here we realize atomic-resolution imaging of the basal (0001) surface structure of hexagonal water ice (ice Ih) by using qPlus-based cryogenic atomic force microscopy with a carbon monoxide-functionalized tip. We find that the crystalline ice-Ih surface consists of mixed Ih- and cubic (Ic)-stacking nanodomains, forming 19 × 19 periodic superstructures. Density functional theory reveals that this reconstructed surface is stabilized over the ideal ice surface mainly by minimizing the electrostatic repulsion between dangling OH bonds. Moreover, we observe that the ice surface gradually becomes disordered with increasing temperature (above 120 Kelvin), indicating the onset of the premelting process. The surface premelting occurs from the defective boundaries between the Ih and Ic domains and can be promoted by the formation of a planar local structure. These results put an end to the longstanding debate on ice surface structures and shed light on the molecular origin of ice premelting, which may lead to a paradigm shift in the understanding of ice physics and chemistry.
Collapse
Affiliation(s)
- Jiani Hong
- International Center for Quantum Materials, School of Physics, Peking University, Beijing, People's Republic of China
| | - Ye Tian
- International Center for Quantum Materials, School of Physics, Peking University, Beijing, People's Republic of China.
| | - Tiancheng Liang
- International Center for Quantum Materials, School of Physics, Peking University, Beijing, People's Republic of China
| | - Xinmeng Liu
- International Center for Quantum Materials, School of Physics, Peking University, Beijing, People's Republic of China
| | - Yizhi Song
- International Center for Quantum Materials, School of Physics, Peking University, Beijing, People's Republic of China
| | - Dong Guan
- International Center for Quantum Materials, School of Physics, Peking University, Beijing, People's Republic of China
| | - Zixiang Yan
- International Center for Quantum Materials, School of Physics, Peking University, Beijing, People's Republic of China
| | - Jiadong Guo
- International Center for Quantum Materials, School of Physics, Peking University, Beijing, People's Republic of China
| | - Binze Tang
- International Center for Quantum Materials, School of Physics, Peking University, Beijing, People's Republic of China
| | - Duanyun Cao
- Beijing Key Laboratory of Environmental Science and Engineering, School of Materials Science and Engineering, Beijing Institute of Technology, Beijing, People's Republic of China
- Beijing Institute of Technology Chongqing Innovation Center, Chongqing, People's Republic of China
| | - Jing Guo
- College of Chemistry, Beijing Normal University, Beijing, People's Republic of China
| | - Ji Chen
- School of Physics, Peking University, Beijing, People's Republic of China
- Interdisciplinary Institute of Light-Element Quantum Materials and Research Center for Light-Element Advanced Materials, Peking University, Beijing, People's Republic of China
| | - Ding Pan
- Department of Physics and Department of Chemistry, The Hong Kong University of Science and Technology, Hong Kong, People's Republic of China
| | - Li-Mei Xu
- International Center for Quantum Materials, School of Physics, Peking University, Beijing, People's Republic of China.
- Interdisciplinary Institute of Light-Element Quantum Materials and Research Center for Light-Element Advanced Materials, Peking University, Beijing, People's Republic of China.
- Collaborative Innovation Center of Quantum Matter, Beijing, People's Republic of China.
| | - En-Ge Wang
- International Center for Quantum Materials, School of Physics, Peking University, Beijing, People's Republic of China.
- Interdisciplinary Institute of Light-Element Quantum Materials and Research Center for Light-Element Advanced Materials, Peking University, Beijing, People's Republic of China.
- Collaborative Innovation Center of Quantum Matter, Beijing, People's Republic of China.
- Tsientang Institute for Advanced Study, Zhejiang, People's Republic of China.
| | - Ying Jiang
- International Center for Quantum Materials, School of Physics, Peking University, Beijing, People's Republic of China.
- Interdisciplinary Institute of Light-Element Quantum Materials and Research Center for Light-Element Advanced Materials, Peking University, Beijing, People's Republic of China.
- Collaborative Innovation Center of Quantum Matter, Beijing, People's Republic of China.
- New Cornerstone Science Laboratory, Peking University, Beijing, People's Republic of China.
| |
Collapse
|
9
|
Yasuda I, Endo K, Arai N, Yasuoka K. In-layer inhomogeneity of molecular dynamics in quasi-liquid layers of ice. Commun Chem 2024; 7:117. [PMID: 38811834 PMCID: PMC11136980 DOI: 10.1038/s42004-024-01197-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2023] [Accepted: 05/02/2024] [Indexed: 05/31/2024] Open
Abstract
Quasi-liquid layers (QLLs) are present on the surface of ice and play a significant role in its distinctive chemical and physical properties. These layers exhibit considerable heterogeneity across different scales ranging from nanometers to millimeters. Although the formation of partially ice-like structures has been proposed, the molecular-level understanding of this heterogeneity remains unclear. Here, we examined the heterogeneity of molecular dynamics on QLLs based on molecular dynamics simulations and machine learning analysis of the simulation data. We demonstrated that the molecular dynamics of QLLs do not comprise a mixture of solid- and liquid water molecules. Rather, molecules having similar behaviors form dynamical domains that are associated with the dynamical heterogeneity of supercooled water. Nonetheless, molecules in the domains frequently switch their dynamical state. Furthermore, while there is no observable characteristic domain size, the long-range ordering strongly depends on the temperature and crystal face. Instead of a mixture of static solid- and liquid-like regions, our results indicate the presence of heterogeneous molecular dynamics in QLLs, which offers molecular-level insights into the surface properties of ice.
Collapse
Affiliation(s)
- Ikki Yasuda
- Department of Mechanical Engineering, Keio University, Yokohama, Japan
| | - Katsuhiro Endo
- Department of Mechanical Engineering, Keio University, Yokohama, Japan
| | - Noriyoshi Arai
- Department of Mechanical Engineering, Keio University, Yokohama, Japan
| | - Kenji Yasuoka
- Department of Mechanical Engineering, Keio University, Yokohama, Japan.
| |
Collapse
|
10
|
Elías-Domínguez A, Alvarado JFJ, Pérez-Villaseñor F, Ortíz-Arroyo A, Castro-Agüero Á, López-Medina F, Medina-Velázquez DY. Computer Simulation of Three-Phase Equilibria for Some Water/ n-Alkane Binary Systems. J Phys Chem B 2024; 128:5044-5054. [PMID: 38727627 DOI: 10.1021/acs.jpcb.4c00992] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/24/2024]
Abstract
In this work, the vapor-liquid-liquid equilibrium (VLLE) of the water/n-pentane, water/n-hexane, water/n-octane, and water/n-decane binary systems is calculated by computer simulation using the NVT-Gibbs ensemble (in the version of three simulation boxes) combined with the configurational bias Monte Carlo method. The combination of both methods, the molecular potential models used, and the simulation details allowed us to calculate the triphasic equilibrium properties of the systems studied: the densities of the three phases in equilibrium, their compositions, and potential energies. In previous works, these simulations were not carried out at a temperature range nor water/n-alkanes systems simulated in this work, probably because they are highly nonideal systems; so, to the best of our knowledge, this is the first time that this phenomenon is studied in detail. The results from VLLE simulations of the water/n-pentane system for temperatures from 343.2 to 435 K, the water/n-hexane system for temperatures from 373.11 to 473.15 K, the water/n-octane system for temperatures from 310.9 to 500 K, and for the water/n-decane system for temperatures from 374.15 to 525 K are reported here. The temperature range was selected in concordance with the experimental data available for an adequate study of the VLLE simulation results. The subcritical densities (vapor and liquid rich in n-alkane phases) at various temperatures fit well with the scaling law and the law of rectilinear diameters, allowing the estimation of upper critical end point temperature and density of the VLLE. The simulation results show a good prediction with experimental data reports in the literature.
Collapse
Affiliation(s)
- Arturo Elías-Domínguez
- Facultad de Ciencias Básicas, Ingeniería y Tecnología, Universidad Autónoma de Tlaxcala, Ángel Solana S/N, San Luis Apizaquito, Apizaco, Tlaxcala CP 90341, México
| | - Juan F J Alvarado
- Departamento de Ingeniería Química, Instituto Tecnológico de Celaya, Av. Tecnológico y A. García Cubas S/N, Celaya, Guanajuato CP 38010, México
| | - Fernando Pérez-Villaseñor
- Facultad de Ciencias Básicas, Ingeniería y Tecnología, Universidad Autónoma de Tlaxcala, Ángel Solana S/N, San Luis Apizaquito, Apizaco, Tlaxcala CP 90341, México
| | - Arturo Ortíz-Arroyo
- Facultad de Ciencias Básicas, Ingeniería y Tecnología, Universidad Autónoma de Tlaxcala, Ángel Solana S/N, San Luis Apizaquito, Apizaco, Tlaxcala CP 90341, México
| | - Ángel Castro-Agüero
- Facultad de Ciencias Básicas, Ingeniería y Tecnología, Universidad Autónoma de Tlaxcala, Ángel Solana S/N, San Luis Apizaquito, Apizaco, Tlaxcala CP 90341, México
| | - Friné López-Medina
- Facultad de Ciencias Básicas, Ingeniería y Tecnología, Universidad Autónoma de Tlaxcala, Ángel Solana S/N, San Luis Apizaquito, Apizaco, Tlaxcala CP 90341, México
| | - Dulce Y Medina-Velázquez
- Facultad de Ciencias Básicas, Ingeniería y Tecnología, Universidad Autónoma de Tlaxcala, Ángel Solana S/N, San Luis Apizaquito, Apizaco, Tlaxcala CP 90341, México
| |
Collapse
|
11
|
Bassani CL, Engel M, Sum AK. Mesomorphology of clathrate hydrates from molecular ordering. J Chem Phys 2024; 160:190901. [PMID: 38767264 DOI: 10.1063/5.0200516] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2024] [Accepted: 03/13/2024] [Indexed: 05/22/2024] Open
Abstract
Clathrate hydrates are crystals formed by guest molecules that stabilize cages of hydrogen-bonded water molecules. Whereas thermodynamic equilibrium is well described via the van der Waals and Platteeuw approach, the increasing concerns with global warming and energy transition require extending the knowledge to non-equilibrium conditions in multiphase, sheared systems, in a multiscale framework. Potential macro-applications concern the storage of carbon dioxide in the form of clathrates, and the reduction of hydrate inhibition additives currently required in hydrocarbon production. We evidence porous mesomorphologies as key to bridging the molecular scales to macro-applications of low solubility guests. We discuss the coupling of molecular ordering with the mesoscales, including (i) the emergence of porous patterns as a combined factor from the walk over the free energy landscape and 3D competitive nucleation and growth and (ii) the role of molecular attachment rates in crystallization-diffusion models that allow predicting the timescale of pore sealing. This is a perspective study that discusses the use of discrete models (molecular dynamics) to build continuum models (phase field models, crystallization laws, and transport phenomena) to predict multiscale manifestations at a feasible computational cost. Several advances in correlated fields (ice, polymers, alloys, and nanoparticles) are discussed in the scenario of clathrate hydrates, as well as the challenges and necessary developments to push the field forward.
Collapse
Affiliation(s)
- Carlos L Bassani
- Institute for Multiscale Simulation, Friedrich-Alexander-Universität Erlangen-Nürnberg, 91058 Erlangen, Germany
| | - Michael Engel
- Institute for Multiscale Simulation, Friedrich-Alexander-Universität Erlangen-Nürnberg, 91058 Erlangen, Germany
| | - Amadeu K Sum
- Phases to Flow Laboratory, Chemical and Biological Engineering Department, Colorado School of Mines, Golden, Colorado 80401, USA
| |
Collapse
|
12
|
Jia H, Fan F, Wang Q, Shen Z, Wang Y, Sun H, Pei P, Li C, Lv K, Huang P. Molecular Insights into the Dual Promotion-Inhibition Effects of NaCl at Various Concentrations on the CO 2 Hydrate Growth: A Molecular Simulation Study. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2024; 40:9012-9019. [PMID: 38625688 DOI: 10.1021/acs.langmuir.4c00315] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/17/2024]
Abstract
Hydrate-based CO2 storage in the ocean is considered a potential method for mitigating the greenhouse effect. Numerous studies demonstrated that NaCl exhibited the dual effects of promotion and inhibition in the nucleation and growth processes of CO2 hydrate, whose mechanisms remain unclear. In this study, the effects of NaCl at various concentrations on the CO2 hydrate growth and crystal are investigated. The independent gradient model based on Hirshfeld partition, electrostatic potential, and binding energy is conducted to study the interaction between ions and water molecules. The motion trajectories of ions are observed at the molecular level to reflect the impact of ion motion on hydrate growth. The results show that the influence of NaCl on hydrate growth depends on a delicate balance of dual promotion-inhibition effects. NaCl can combine more water molecules and provide a transport channel of CO2 to promote hydrate growth at low concentrations. Meanwhile, the promoting effects shift toward inhibition with increasing NaCl concentrations. In a word, this paper proposes a novel mechanism for the dual promotion-inhibition effects of NaCl on hydrate growth, which is significant for further research on hydrate-based CO2 storage in the ocean.
Collapse
Affiliation(s)
- Han Jia
- Shandong Key Laboratory of Oilfield Chemistry, School of Petroleum Engineering, China University of Petroleum (East China), Qingdao 266580, China
- Key Laboratory of Unconventional Oil & Gas Development, China University of Petroleum (East China), Ministry of Education, Qingdao 266580, China
| | - Fangning Fan
- Shandong Key Laboratory of Oilfield Chemistry, School of Petroleum Engineering, China University of Petroleum (East China), Qingdao 266580, China
- Key Laboratory of Unconventional Oil & Gas Development, China University of Petroleum (East China), Ministry of Education, Qingdao 266580, China
| | - Qiuxia Wang
- Bohai Oilfield Research Institute, Tianjin Branch, CNOOC China Limited, Tianjin 300459, China
| | - Zhihao Shen
- Shandong Key Laboratory of Oilfield Chemistry, School of Petroleum Engineering, China University of Petroleum (East China), Qingdao 266580, China
- Key Laboratory of Unconventional Oil & Gas Development, China University of Petroleum (East China), Ministry of Education, Qingdao 266580, China
| | - Yuanbo Wang
- Shandong Key Laboratory of Oilfield Chemistry, School of Petroleum Engineering, China University of Petroleum (East China), Qingdao 266580, China
- Key Laboratory of Unconventional Oil & Gas Development, China University of Petroleum (East China), Ministry of Education, Qingdao 266580, China
| | - Han Sun
- Shandong Key Laboratory of Oilfield Chemistry, School of Petroleum Engineering, China University of Petroleum (East China), Qingdao 266580, China
- Key Laboratory of Unconventional Oil & Gas Development, China University of Petroleum (East China), Ministry of Education, Qingdao 266580, China
| | - Pingan Pei
- Shandong Key Laboratory of Oilfield Chemistry, School of Petroleum Engineering, China University of Petroleum (East China), Qingdao 266580, China
- Key Laboratory of Unconventional Oil & Gas Development, China University of Petroleum (East China), Ministry of Education, Qingdao 266580, China
| | - Chuanqi Li
- Shandong Key Laboratory of Oilfield Chemistry, School of Petroleum Engineering, China University of Petroleum (East China), Qingdao 266580, China
- Key Laboratory of Unconventional Oil & Gas Development, China University of Petroleum (East China), Ministry of Education, Qingdao 266580, China
| | - Kaihe Lv
- Shandong Key Laboratory of Oilfield Chemistry, School of Petroleum Engineering, China University of Petroleum (East China), Qingdao 266580, China
- Key Laboratory of Unconventional Oil & Gas Development, China University of Petroleum (East China), Ministry of Education, Qingdao 266580, China
| | - Pan Huang
- Department of Chemical and Materials Engineering, University of Alberta, Edmonton, AB T6G 1H9, Canada
| |
Collapse
|
13
|
Algaba J, Blazquez S, Feria E, Míguez JM, Conde MM, Blas FJ. Three-phase equilibria of hydrates from computer simulation. II. Finite-size effects in the carbon dioxide hydrate. J Chem Phys 2024; 160:164722. [PMID: 38687000 DOI: 10.1063/5.0201306] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2024] [Accepted: 04/03/2024] [Indexed: 05/02/2024] Open
Abstract
In this work, the effects of finite size on the determination of the three-phase coexistence temperature (T3) of the carbon dioxide (CO2) hydrate have been studied by molecular dynamic simulations and using the direct coexistence technique. According to this technique, the three phases involved (hydrate-aqueous solution-liquid CO2) are placed together in the same simulation box. By varying the number of molecules of each phase, it is possible to analyze the effect of simulation size and stoichiometry on the T3 determination. In this work, we have determined the T3 value at 8 different pressures (from 100 to 6000 bar) and using 6 different simulation boxes with different numbers of molecules and sizes. In two of these configurations, the ratio of the number of water and CO2 molecules in the aqueous solution and the liquid CO2 phase is the same as in the hydrate (stoichiometric configuration). In both stoichiometric configurations, the formation of a liquid drop of CO2 in the aqueous phase is observed. This drop, which has a cylindrical geometry, increases the amount of CO2 available in the aqueous solution and can in some cases lead to the crystallization of the hydrate at temperatures above T3, overestimating the T3 value obtained from direct coexistence simulations. The simulation results obtained for the CO2 hydrate confirm the sensitivity of T3 depending on the size and composition of the system, explaining the discrepancies observed in the original work by Míguez et al. [J. Chem Phys. 142, 124505 (2015)]. Non-stoichiometric configurations with larger unit cells show a convergence of T3 values, suggesting that finite-size effects for these system sizes, regardless of drop formation, can be safely neglected. The results obtained in this work highlight that the choice of a correct initial configuration is essential to accurately estimate the three-phase coexistence temperature of hydrates by direct coexistence simulations.
Collapse
Affiliation(s)
- J Algaba
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21007 Huelva, Spain
| | - S Blazquez
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - E Feria
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21007 Huelva, Spain
| | - J M Míguez
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21007 Huelva, Spain
| | - M M Conde
- Departamento de Ingeniería Química Industrial y del Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006 Madrid, Spain
| | - F J Blas
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21007 Huelva, Spain
| |
Collapse
|
14
|
Blazquez S, Algaba J, Míguez JM, Vega C, Blas FJ, Conde MM. Three-phase equilibria of hydrates from computer simulation. I. Finite-size effects in the methane hydrate. J Chem Phys 2024; 160:164721. [PMID: 38686998 DOI: 10.1063/5.0201295] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2024] [Accepted: 04/01/2024] [Indexed: 05/02/2024] Open
Abstract
Clathrate hydrates are vital in energy research and environmental applications. Understanding their stability is crucial for harnessing their potential. In this work, we employ direct coexistence simulations to study finite-size effects in the determination of the three-phase equilibrium temperature (T3) for methane hydrates. Two popular water models, TIP4P/Ice and TIP4P/2005, are employed, exploring various system sizes by varying the number of molecules in the hydrate, liquid, and gas phases. The results reveal that finite-size effects play a crucial role in determining T3. The study includes nine configurations with varying system sizes, demonstrating that smaller systems, particularly those leading to stoichiometric conditions and bubble formation, may yield inaccurate T3 values. The emergence of methane bubbles within the liquid phase, observed in smaller configurations, significantly influences the behavior of the system and can lead to erroneous temperature estimations. Our findings reveal finite-size effects on the calculation of T3 by direct coexistence simulations and clarify the system size convergence for both models, shedding light on discrepancies found in the literature. The results contribute to a deeper understanding of the phase equilibrium of gas hydrates and offer valuable information for future research in this field.
Collapse
Affiliation(s)
- S Blazquez
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J Algaba
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - J M Míguez
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - C Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - F J Blas
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - M M Conde
- Departamento de Ingeniería Química Industrial y del Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006 Madrid, Spain
| |
Collapse
|
15
|
Algaba J, Blazquez S, Míguez JM, Conde MM, Blas FJ. Three-phase equilibria of hydrates from computer simulation. III. Effect of dispersive interactions in the methane and carbon dioxide hydrates. J Chem Phys 2024; 160:164723. [PMID: 38686999 DOI: 10.1063/5.0201309] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2024] [Accepted: 04/02/2024] [Indexed: 05/02/2024] Open
Abstract
In this work, the effect of the range of dispersive interactions in determining the three-phase coexistence line of the CO2 and CH4 hydrates has been studied. In particular, the temperature (T3) at which solid hydrate, water, and liquid CO2/gas CH4 coexist has been determined through molecular dynamics simulations using different cutoff values (from 0.9 to 1.6 nm) for dispersive interactions. The T3 of both hydrates has been determined using the direct coexistence simulation technique. Following this method, the three phases in equilibrium are put together in the same simulation box, the pressure is fixed, and simulations are performed at different temperatures T. If the hydrate melts, then T > T3. Conversely, if the hydrate grows, then T < T3. The effect of the cutoff distance on the dissociation temperature has been analyzed at three different pressures for CO2 hydrate: 100, 400, and 1000 bar. Then, we have changed the guest and studied the effect of the cutoff distance on the dissociation temperature of the CH4 hydrate at 400 bar. Moreover, the effect of long-range corrections for dispersive interactions has been analyzed by running simulations with homo- and inhomogeneous corrections and a cutoff value of 0.9 nm. The results obtained in this work highlight that the cutoff distance for the dispersive interactions affects the stability conditions of these hydrates. This effect is enhanced when the pressure is decreased, displacing the T3 about 2-4 K depending on the system and the pressure.
Collapse
Affiliation(s)
- J Algaba
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - S Blazquez
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J M Míguez
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - M M Conde
- Departamento de Ingeniería Química Industrial y del Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006 Madrid, Spain
| | - F J Blas
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| |
Collapse
|
16
|
Jing X, Zhou L, Ma Y, Fu Z, Huang Q, Zhang Z. Physical breakdown of CH 4 hydrate under stress: a molecular dynamics simulation study. BMC Chem 2024; 18:86. [PMID: 38678235 PMCID: PMC11056054 DOI: 10.1186/s13065-024-01191-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2023] [Accepted: 04/17/2024] [Indexed: 04/29/2024] Open
Abstract
As a solid energy source, CH4 hydrate will inevitably break down physically as the result of geological movement or exploitation. Here, the molecular dynamics method was employed to simulate the uniaxial-deformation behavior of structure I (sI type) CH4 hydrate under stress. The stress increases regardless of whether the hydrate is stretched or squeezed, and other physical parameters also changed, such as hydrate cage numbers, order parameters, and the number of water molecules. A noticeable difference is observed between the two systems. Upon stretching, the stress immediately recovers to 0 GPa once the hydrate is completely stretched apart. During the squeeze process, the stress is ultimately not zero since solid and liquid are always in contact. When the hydrate is stretched apart, about 5% of water molecules change from solid to liquid, about 7.8% of CH4 molecules lose their shelter and become free due to the disintegration of water cages. While in the squeezing process, large cages (51262) are crushed more easily than small cages (512); in the end, about 93.5% of large cages and 73% of small cages are crushed, and approximately 87.5% CH4 is released from the cages. In mining CH4 hydrates, caution must be exercised, as if the hydrates break as a result of stress, a large release of CH4 may pose a security risk.
Collapse
Affiliation(s)
- Xianwu Jing
- Research Institute of Natural Gas Technology, PetroChina Southwest Oil and Gasfield Company, Chengdu, 610213, Sichuan, People's Republic of China.
- Shale Gas Evaluation and Exploitation Key Laboratory of Sichuan Province, Sichuan Provincial Department of Science and Technology, Chengdu, 610051, Sichuan, People's Republic of China.
| | - Li Zhou
- Research Institute of Natural Gas Technology, PetroChina Southwest Oil and Gasfield Company, Chengdu, 610213, Sichuan, People's Republic of China
| | - Yong Ma
- Engineering Technology Department, PetroChina Southwest Oil and Gasfield Company, Chengdu, 610081, Sichuan, People's Republic of China
| | - Ziyi Fu
- Research Institute of Natural Gas Technology, PetroChina Southwest Oil and Gasfield Company, Chengdu, 610213, Sichuan, People's Republic of China
| | - Qian Huang
- PetroChina Planning and Engineering Institute, Beijing, 100083, China
| | - Zhe Zhang
- PetroChina Planning and Engineering Institute, Beijing, 100083, China
| |
Collapse
|
17
|
Fernández-Fernández ÁM, Bárcena Á, Conde MM, Pérez-Sánchez G, Pérez-Rodríguez M, Piñeiro MM. Modeling oceanic sedimentary methane hydrate growth through molecular dynamics simulation. J Chem Phys 2024; 160:144107. [PMID: 38591679 DOI: 10.1063/5.0203116] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2024] [Accepted: 03/23/2024] [Indexed: 04/10/2024] Open
Abstract
The crystallization process of methane hydrates in a confined geometry resembling seabed porous silica sedimentary conditions has been studied using molecular dynamics simulations. With this objective in mind, a fully atomistic quartz silica slit pore has been designed, and the temperature stability of a methane hydrate crystalline seed in the presence of water and guest molecule methane has been analyzed. NaCl ion pairs have been added in different concentrations, simulating salinity conditions up to values higher than average oceanic conditions. The structure obtained when the hydrate crystallizes inside the pore is discussed, paying special attention to the presence of ionic doping inside the hydrate and the subsequent induced structural distortion. The shift in the hydrate stability conditions due to the increasing water salinity is discussed and compared with the case of unconfined hydrate, concluding that the influence of the confinement geometry and pore hydrophilicity produces a larger deviation in the confined hydrate phase equilibria.
Collapse
Affiliation(s)
| | - Álvaro Bárcena
- Dpto. de Física Aplicada, Univ. de Vigo, Vigo 36310, Spain
| | - María M Conde
- Dpto. de Ingeniería Química Industrial y Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, Madrid 28006, Spain
| | | | | | | |
Collapse
|
18
|
Weldon R, Wang F. Water Potential from Adaptive Force Matching for Ice and Liquid with Revised Dispersion Predicts Supercooled Liquid Anomalies in Good Agreement with Two Independent Experimental Fits. J Phys Chem B 2024; 128:3398-3407. [PMID: 38536126 PMCID: PMC11017247 DOI: 10.1021/acs.jpcb.3c06495] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2023] [Revised: 02/15/2024] [Accepted: 03/11/2024] [Indexed: 04/12/2024]
Abstract
A revised version of the Water potential from Adaptive force matching for Ice and Liquid (WAIL) was developed by using the previous data set for fitting the WAIL model but with a dispersion term calculated using symmetry adapted perturbation theory (SAPT). The model has no adjustable parameters and relies solely on fitting first-principles information. The new model, named revised WAIL (rWAIL), shows improved predictions of most properties of water when compared to the previously published WAIL model. The rWAIL model also compares favorably to other first-principles-derived water models, such as MB-Pol, at only a fraction of the computational cost. The rWAIL model is used to study the properties of supercooled water. The model shows evidence of a liquid-liquid phase transition (LLPT) in the supercooled regimes with the liquid-liquid critical point (LLCP) at 203 K and 90 MPa. This estimate is in good agreement with a recent polynomial fit to the experimental density of water. Also, the fit to the surface tension of supercooled water based on the rWAIL model shows excellent agreement with the corresponding fit to the experimental data. Consistent with previously published molecular dynamics and experimental data, the surface tension of water exhibits exponential growth in the supercooled regime, which is likely a result of the emergence of a low-density liquid form of water. The simulation thus unites two separate experimental fits with one first-principles-based model, lending strong evidence of an LLPT in real water.
Collapse
Affiliation(s)
- Raymond Weldon
- Department of Chemistry and
Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, United States
| | - Feng Wang
- Department of Chemistry and
Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, United States
| |
Collapse
|
19
|
Montero de Hijes P, Dellago C, Jinnouchi R, Schmiedmayer B, Kresse G. Comparing machine learning potentials for water: Kernel-based regression and Behler-Parrinello neural networks. J Chem Phys 2024; 160:114107. [PMID: 38506284 DOI: 10.1063/5.0197105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2024] [Accepted: 03/03/2024] [Indexed: 03/21/2024] Open
Abstract
In this paper, we investigate the performance of different machine learning potentials (MLPs) in predicting key thermodynamic properties of water using RPBE + D3. Specifically, we scrutinize kernel-based regression and high-dimensional neural networks trained on a highly accurate dataset consisting of about 1500 structures, as well as a smaller dataset, about half the size, obtained using only on-the-fly learning. This study reveals that despite minor differences between the MLPs, their agreement on observables such as the diffusion constant and pair-correlation functions is excellent, especially for the large training dataset. Variations in the predicted density isobars, albeit somewhat larger, are also acceptable, particularly given the errors inherent to approximate density functional theory. Overall, this study emphasizes the relevance of the database over the fitting method. Finally, this study underscores the limitations of root mean square errors and the need for comprehensive testing, advocating the use of multiple MLPs for enhanced certainty, particularly when simulating complex thermodynamic properties that may not be fully captured by simpler tests.
Collapse
Affiliation(s)
- Pablo Montero de Hijes
- University of Vienna, Faculty of Physics, Kolingasse 14, A-1090 Vienna, Austria
- University of Vienna, Faculty of Earth Sciences, Geography and Astronomy, Josef-Holaubuek-Platz 2, 1090 Vienna, Austria
| | - Christoph Dellago
- University of Vienna, Faculty of Physics, Kolingasse 14, A-1090 Vienna, Austria
| | - Ryosuke Jinnouchi
- Toyota Central R&D Labs., Inc., 41-1 Yokomichi, Nagakute, Aichi 480-1192, Japan
| | | | - Georg Kresse
- University of Vienna, Faculty of Physics, Kolingasse 14, A-1090 Vienna, Austria
- VASP Software GmbH, Berggasse 21, A-1090 Vienna, Austria
| |
Collapse
|
20
|
Cao P, Wu J, Ning F. Mechanical properties of amorphous CO 2 hydrates: insights from molecular simulations. Phys Chem Chem Phys 2024; 26:9388-9398. [PMID: 38444360 DOI: 10.1039/d4cp00203b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2024]
Abstract
Understanding physicochemical properties of amorphous gas hydrate systems is of great significance to reveal structural stabilities of polycrystalline gas hydrate systems. Furthermore, amorphous gas hydrates can occur ordinarily in the nucleation events of gas hydrate systems. Herein, the mechanical properties of amorphous carbon dioxide hydrates are examined by means of all-atom classical molecular dynamic simulations. Our molecular simulation results reveal that mechanical strengths of amorphous carbon dioxide hydrates are evidently governed by temperatures, confining pressures, and ratios of water to carbon dioxide molecules. Notably, under compressive loads, amorphous carbon dioxide hydrates firstly exhibit monotonic strain hardening, followed by an interesting distinct phenomenon characterized by a steady flow stress at further large deformation strains. Furthermore, structural evolutions of amorphous carbon dioxide hydrates are analyzed on the basis of the N-Hbond DOP order parameter. These important findings can not only contribute to our understanding of the structural stabilities of amorphous gas hydrate systems, but also help to develop fundamental understandings about grain boundaries of gas hydrate systems.
Collapse
Affiliation(s)
- Pinqiang Cao
- School of Resource and Environmental Engineering, Wuhan University of Science and Technology, Wuhan, Hubei 430081, China.
| | - Jianyang Wu
- Department of Physics, Research Institute for Biomimetics and Soft Matter, Fujian Provincial Key Laboratory for Soft Functional Materials Research, Xiamen University, Xiamen 361005, China.
| | - Fulong Ning
- Faculty of Engineering, China University of Geosciences, Wuhan, Hubei 430074, China.
| |
Collapse
|
21
|
Song H, Bei Z, Voronin AS, Umaiya Kunjaram UP, Truscott TT, Schwingenschlögl U, Vrouwenvelder JS, Gan Q. A robust thin-film droplet-induced electricity generator. iScience 2024; 27:109291. [PMID: 38450151 PMCID: PMC10915600 DOI: 10.1016/j.isci.2024.109291] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2023] [Revised: 02/08/2024] [Accepted: 02/16/2024] [Indexed: 03/08/2024] Open
Abstract
The pursuit of cost-effective, high-voltage electricity generators activated by droplets represents a new frontier in hydropower technology. This study presents an economical method for crafting droplet generators using common materials such as solid polytetrafluoroethylene (PTFE) films and readily available tapes, eliminating the need for specialized cleanroom facilities. A thorough investigation into voltage-limiting factors, encompassing device capacitance and induced electrode charges, reveals specific areas with potential for optimization. A substantial enhancement in the open-circuit voltage (Voc) was achieved, reaching approximately 282.2 ± 27.9 V-an impressive increase of around 60 V compared to earlier benchmarks. One device showcased its capability to power 100 LEDs concurrently, underscoring its efficacy. Ten such devices created diverse luminous patterns with uniform light intensity for each LED, showcasing the practical potential of the approach. The methodology's cost-effectiveness results in a remarkable cost reduction compared to solution-based materials, paving the way for the widespread adoption of large-scale water droplet energy harvesting.
Collapse
Affiliation(s)
- Haomin Song
- Material Science Engineering, Physical Science Engineering Division, King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia
| | - Zongmin Bei
- Shared Instrumentation Laboratories, School of Engineering & Applied Sciences, The State University of New York at Buffalo, Buffalo, NY 14260, USA
- Department of Electrical Engineering, The State University of New York at Buffalo, Buffalo, NY 14260, USA
| | - Aleksandr S. Voronin
- Applied Physics, Physical Science Engineering Division, King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia
| | | | - Tadd T. Truscott
- Mechanical Engineering, Physical Science Engineering Division, King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia
| | - Udo Schwingenschlögl
- Applied Physics, Physical Science Engineering Division, King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia
| | - Johannes S. Vrouwenvelder
- Water Desalination and Reuse Center, Division of Biological Sciences and Engineering, King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia
| | - Qiaoqiang Gan
- Material Science Engineering, Physical Science Engineering Division, King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia
- Department of Electrical Engineering, The State University of New York at Buffalo, Buffalo, NY 14260, USA
| |
Collapse
|
22
|
Quoika PK, Zacharias M. Liquid-Vapor Coexistence and Spontaneous Evaporation at Atmospheric Pressure of Common Rigid Three-Point Water Models in Molecular Simulations. J Phys Chem B 2024; 128:2457-2468. [PMID: 38427971 PMCID: PMC10945489 DOI: 10.1021/acs.jpcb.3c08183] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2023] [Revised: 02/21/2024] [Accepted: 02/22/2024] [Indexed: 03/03/2024]
Abstract
Molecular dynamics (MD) simulations are widely used to investigate molecular systems at atomic resolution including biomolecular structures, drug-receptor interactions, and novel materials. Frequently, MD simulations are performed in an aqueous solution with explicit models of water molecules. Commonly, such models are parameterized to reproduce the liquid phase of water under ambient conditions. However, often, simulations at significantly higher temperatures are also of interest. Hence, it is important to investigate the equilibrium of the liquid and vapor phases of molecular models of water at elevated temperatures. Here, we evaluate the behavior of 11 common rigid three-point water models over a wide range of temperatures. From liquid-vapor coexistence simulations, we estimated the critical points and studied the spontaneous evaporation of these water models. Moreover, we investigated the influence of the system size, choice of the pressure-coupling algorithm, and rate of heating on the process and compared them with the experimental data. We found that modern rigid three-point water models reproduce the critical point surprisingly well. Furthermore, we discovered that the critical temperature correlates with the quadrupole moment of the respective water model. This indicates that the spatial arrangement of the partial charges is important for reproducing the liquid-vapor phase transition. Our findings may guide the selection of water models for simulations conducted at high temperatures.
Collapse
Affiliation(s)
- Patrick K. Quoika
- Center for Functional Protein
Assemblies, Technical University of Munich, Ernst-Otto-Fischer-Str. 8, Garching 85748, Germany
| | - Martin Zacharias
- Center for Functional Protein
Assemblies, Technical University of Munich, Ernst-Otto-Fischer-Str. 8, Garching 85748, Germany
| |
Collapse
|
23
|
Grabowska J, Kuffel A, Zielkiewicz J. Long-range, water-mediated interaction between a moderately active antifreeze protein molecule and the surface of ice. J Chem Phys 2024; 160:095101. [PMID: 38445741 DOI: 10.1063/5.0187663] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2023] [Accepted: 02/16/2024] [Indexed: 03/07/2024] Open
Abstract
Using molecular dynamics simulations, we show that a molecule of moderately active antifreeze protein (type III AFP, QAE HPLC-12 isoform) is able to interact with ice in an indirect manner. This interaction occurs between the ice binding site (IBS) of the AFP III molecule and the surface of ice, and it is mediated by liquid water, which separates these surfaces. As a result, the AFP III molecule positions itself at a specific orientation and distance relative to the surface of ice, which enables the effective binding (via hydrogen bonds) of the molecule with the nascent ice surface. Our results show that the final adsorption of the AFP III molecule on the surface of ice is not achieved by chaotic diffusion movements, but it is preceded by a remote, water-mediated interaction between the IBS and the surface of ice. The key factor that determines the existence of this interaction is the ability of water molecules to spontaneously form large, high-volume aggregates that can be anchored to both the IBS of the AFP molecule and the surface of ice. The results presented in this work for AFP III are in full agreement with the ones obtained by us previously for hyperactive CfAFP, which indicates that the mechanism of the remote interaction of these molecules with ice remains unchanged despite significant differences in the molecular structure of their ice binding sites. For that reason, we can expect that also other types of AFPs interact with the ice surface according to an analogous mechanism.
Collapse
Affiliation(s)
- Joanna Grabowska
- Faculty of Chemistry, Department of Physical Chemistry, Gdańsk University of Technology, Narutowicza 11/12, 80-233 Gdańsk, Poland
| | - Anna Kuffel
- Faculty of Chemistry, Department of Physical Chemistry, Gdańsk University of Technology, Narutowicza 11/12, 80-233 Gdańsk, Poland
| | - Jan Zielkiewicz
- Faculty of Chemistry, Department of Physical Chemistry, Gdańsk University of Technology, Narutowicza 11/12, 80-233 Gdańsk, Poland
| |
Collapse
|
24
|
Foffi R, Sciortino F. Identification of local structures in water from supercooled to ambient conditions. J Chem Phys 2024; 160:094504. [PMID: 38436442 DOI: 10.1063/5.0188764] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2023] [Accepted: 02/09/2024] [Indexed: 03/05/2024] Open
Abstract
Studies of water thermodynamics have long been tied to the identification of two distinct families of local structures, whose competition could explain the origin of the many thermodynamic anomalies and the hypothesized liquid-liquid critical point in water. Despite the many successes and insights gained, the structural indicators proposed throughout the years were not able to unequivocally identify these two families over a wide range of conditions. We show that a recently introduced indicator, Ψ, which exploits information on the hydrogen bond network connectivity, can reliably identify these two distinct local environments over a wide range of thermodynamic conditions (188-300 K and 0-13 kbar) and that close to the liquid-liquid critical point, the spatial correlations of density fluctuations are identical to those of the Ψ indicator. Our results strongly support the idea that water thermodynamic properties arise from the competition between two distinct and identifiable local environments.
Collapse
Affiliation(s)
- Riccardo Foffi
- Department of Civil, Environmental and Geomatic Engineering, Institute for Environmental Engineering, ETH Zürich, Laura-Hezner-Weg 7, 8093 Zürich, Switzerland
| | - Francesco Sciortino
- Dipartimento di Fisica, Sapienza Università di Roma, Piazzale Aldo Moro 5, I-00185 Rome, Italy
| |
Collapse
|
25
|
Roudsari G, Lbadaoui-Darvas M, Welti A, Nenes A, Laaksonen A. The molecular scale mechanism of deposition ice nucleation on silver iodide. ENVIRONMENTAL SCIENCE: ATMOSPHERES 2024; 4:243-251. [PMID: 38371604 PMCID: PMC10867811 DOI: 10.1039/d3ea00140g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/21/2023] [Accepted: 12/02/2023] [Indexed: 02/20/2024]
Abstract
Heterogeneous ice nucleation is a ubiquitous process in the natural and built environment. Deposition ice nucleation, i.e. heterogeneous ice nucleation that - according to the traditional view - occurs in a subsaturated water vapor environment and in the absence of supercooled water on the solid, ice-forming surface, is among the most important ice formation processes in high-altitude cirrus and mixed-phase clouds. Despite its importance, very little is known about the mechanism of deposition ice nucleation at the microscopic level. This study puts forward an adsorption-based mechanism for deposition ice nucleation through results from a combination of atomistic simulations, experiments and theoretical modelling. One of the most potent laboratory surrogates of ice nucleating particles, silver iodide, is used as a substrate for the simulations. We find that water initially adsorbs in clusters which merge and grow over time to form layers of supercooled water. Ice nucleation on silver iodide requires at minimum the adsorption of 4 molecular layers of water. Guided by the simulations we propose the following fundamental freezing steps: (1) Water molecules adsorb on the surface, forming nanodroplets. (2) The supercooled water nanodroplets merge into a continuous multilayer when they grow to about 3 molecular layers thick. (3) The layer continues to grow until the critical thickness for freezing is reached. (4) The critical ice cluster continues to grow.
Collapse
Affiliation(s)
| | - Mária Lbadaoui-Darvas
- Laboratory of Atmospheric Processes and their Impacts, ENAC, Ecole Polytechnique Fédérale de Lausanne Lausanne Switzerland
- Institute of Chemical Engineering Sciences, Foundation for Research and Technology Hellas (FORTH/ICE-HT) 26504 Patras Greece
| | - André Welti
- Finnish Meteorological Institute FI-00101 Helsinki Finland
| | - Athanasios Nenes
- Laboratory of Atmospheric Processes and their Impacts, ENAC, Ecole Polytechnique Fédérale de Lausanne Lausanne Switzerland
- Institute of Chemical Engineering Sciences, Foundation for Research and Technology Hellas (FORTH/ICE-HT) 26504 Patras Greece
| | - Ari Laaksonen
- Finnish Meteorological Institute FI-00101 Helsinki Finland
- Department of Applied Physics, University of Eastern Finland Kuopio 70211 Finland
| |
Collapse
|
26
|
Midya US, Bandyopadhyay S. Ice Recrystallization Unveils the Binding Mechanism Operating at a Diffused Interface. J Phys Chem B 2024; 128:1170-1178. [PMID: 38287221 DOI: 10.1021/acs.jpcb.3c05934] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2024]
Abstract
Recrystallization of ice is a natural phenomenon that causes adverse effects in cryopreservation, agriculture, and in frozen food industry. It has long been recognized that ice recrystallization occurs through the Ostwald ripening and accretion processes. However, neither of these processes has been explored in microscopic detail by state-of-the-art experimental techniques. We carried out atomistic molecular dynamics (MD) simulations to explore ice recrystallization through the accretion process. Attempts have been made to elucidate the binding mechanism that is operating at the diffused ice-water interface. It is demonstrated that two ice crystals spontaneously recognize each other and bind together to form a large crystal in liquid water, resulting in ice recrystallization by accretion. Interestingly, the study reveals that the binding occurs due to the freezing of the interfacial water layer present between the two ice planes, even at a temperature above the melting point of the ice crystal. The synergistically enhanced ordering effect of two ice surfaces on the interfacial water leads to such freezing occurring during the binding process. However, proper crystallographic alignment is not necessarily required for the binding of the two crystals. Simulations have also been carried out to study the binding between an ice crystal and the model ice-binding surface (IBS) of an antifreeze protein above the melting point of the ice crystal. It is found that such binding at the IBS is accompanied by freezing of the interfacial water. This establishes that the synergetic ordering-driven freezing of interfacial water is a common binding mechanism at the diffused surfaces of ice crystals. We believe that this mechanism will provide a microscopic understanding of the process of recrystallization inhibition and thus help in designing suitable materials for potent applications in recrystallization inhibition.
Collapse
Affiliation(s)
| | - Sanjoy Bandyopadhyay
- Molecular Modeling Laboratory, Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| |
Collapse
|
27
|
Baran Ł, MacDowell LG. Confinement enhanced viscosity vs shear thinning in lubricated ice friction. J Chem Phys 2024; 160:056101. [PMID: 38341698 DOI: 10.1063/5.0180337] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2023] [Accepted: 01/15/2024] [Indexed: 02/13/2024] Open
Abstract
The ice surface is known for presenting a very small kinetic friction coefficient, but the origin of this property remains highly controversial to date. In this work, we revisit recent computer simulations of ice sliding on atomically smooth substrates, using newly calculated bulk viscosities for the TIP4P/ice water model. The results show that spontaneously formed premelting films in static conditions exhibit an effective viscosity that is about twice the bulk viscosity. However, upon approaching sliding speeds in the order of m/s, the shear rate becomes very large, and the viscosities decrease by several orders of magnitude. This shows that premelting films can act as an efficient lubrication layer despite their small thickness and illustrates an interesting interplay between confinement enhanced viscosities and shear thinning. Our results suggest that the strongly thinned viscosities that operate under the high speed skating regime could largely reduce the amount of frictional heating.
Collapse
Affiliation(s)
- Łukasz Baran
- Department of Theoretical Chemistry, Institute of Chemical Sciences, Faculty of Chemistry, Maria Curie-Sklodowska University in Lublin, Lublin, Poland
| | - Luis G MacDowell
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, Madrid 28040, Spain
| |
Collapse
|
28
|
Blow KE, Whale TF, Quigley D, Sosso GC. Understanding the impact of ammonium ion substitutions on heterogeneous ice nucleation. Faraday Discuss 2024; 249:114-132. [PMID: 37782066 DOI: 10.1039/d3fd00097d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/03/2023]
Abstract
Understanding the mechanisms underpinning heterogeneous ice nucleation in the presence of ionic inclusions is important for fields such as cryopreservation and for improved models of climate and weather prediction. Feldspar and ammonium are both present in significant quantities in the atmosphere, and experimental evidence has shown that feldspar can nucleate ice from ammonium-containing solutions at temperatures warmer than water alone. In recent work, Whale hypothesised that this increase in nucleation temperature is due to an increase in configurational entropy when an ammonium ion is included in the ice hydrogen bond network (T. F. Whale, J. Chem. Phys., 2022, 156, 144503). In this work, we investigate the impact of the inclusion of an ammonium ion on the hydrogen bond network by direct enumeration of the number of structures found using Rick's algorithm. We also determine the energy of these systems and thus compare the effects of enthalpy and entropy to test Whale's hypothesis. We find that the inclusion of an ammonium ion increases the total number of configurations under conditions consistent with a realistic surface charge. We also find that the enthalpic contribution is dominant in determining the location of the ammonium ion within the structure, although we note that this neglects other practicalities of ice nucleation.
Collapse
Affiliation(s)
- Katarina E Blow
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, UK.
| | - Thomas F Whale
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, UK.
| | - David Quigley
- Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, UK.
| | - Gabriele C Sosso
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, UK.
| |
Collapse
|
29
|
Hayton JA, Davies MB, Whale TF, Michaelides A, Cox SJ. The limit of macroscopic homogeneous ice nucleation at the nanoscale. Faraday Discuss 2024; 249:210-228. [PMID: 37791990 DOI: 10.1039/d3fd00099k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
Nucleation in small volumes of water has garnered renewed interest due to the relevance of pore condensation and freezing under conditions of low partial pressures of water, such as in the upper troposphere. Molecular simulations can in principle provide insight on this process at the molecular scale that is challenging to achieve experimentally. However, there are discrepancies in the literature as to whether the rate in confined systems is enhanced or suppressed relative to bulk water at the same temperature and pressure. In this study, we investigate the extent to which the size of the critical nucleus and the rate at which it grows in thin films of water are affected by the thickness of the film. Our results suggest that nucleation remains bulk-like in films that are barely large enough accommodate a critical nucleus. This conclusion seems robust to the presence of physical confining boundaries. We also discuss the difficulties in unambiguously determining homogeneous nucleation rates in nanoscale systems, owing to the challenges in defining the volume. Our results suggest any impact on a film's thickness on the rate is largely inconsequential for present day experiments.
Collapse
Affiliation(s)
- John A Hayton
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK.
| | - Michael B Davies
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK.
- Department of Physics and Astronomy, University College London, London WC1E 6BT, UK
| | - Thomas F Whale
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, UK
- School of Earth and Environment, University of Leeds, Leeds, UK
| | - Angelos Michaelides
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK.
| | - Stephen J Cox
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK.
| |
Collapse
|
30
|
Gao Y, Wang S, Jiang J, Liu Y, Francisco JS, Zeng XC. Evidence of Spontaneous Formation of Two-Dimensional Amorphous Clathrates on Superhydrophilic Surfaces. J Am Chem Soc 2024; 146:2503-2513. [PMID: 38237042 DOI: 10.1021/jacs.3c10701] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2024]
Abstract
Clathrate hydrates reserved in the seabed are often dispersed in the pores of coarse-grained sediments; hence, their formation typically occurs under nanoconfinement. Herein, we show the first molecular dynamics (MD) simulation evidence of the spontaneous formation of two-dimensional (2D) clathrate hydrates on crystal surfaces without conventional nanoconfinement. The kinetic process of 2D clathrate formation is illustrated via simulated single-molecule deposition. 2D amorphous patterns are observed on various superhydrophilic face-centered cubic surfaces. Notably, the formation of 2D amorphous clathrate can occur over a wide range of temperatures, even at room temperature. The strong water-surface interaction, the characteristic properties of guest-gas molecules, and the underlying surface structure dictate the formation of 2D amorphous clathrates. Semiquantitative phase diagrams of 2D clathrates are constructed where representative patterns of 2D clathrates for characteristic gas molecules on prototypical Pd(111) and Pt(111) surfaces are confirmed by independent MD simulations. A tunable pattern of 2D amorphous clathrates is demonstrated by changing the lattice strain of the underlying substrate. Moreover, ab initio MD simulations confirm the stability of 2D amorphous clathrate. The underlining physical mechanism for 2D clathrate formation on superhydrophilic surfaces is elucidated, which offers deeper insight into the crucial role of water-surface interaction.
Collapse
Affiliation(s)
- Yurui Gao
- Laboratory of Theoretical and Computational Nanoscience, National Center for Nanoscience and Technology, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Shixian Wang
- Laboratory of Theoretical and Computational Nanoscience, National Center for Nanoscience and Technology, Chinese Academy of Sciences, Beijing 100190, China
| | - Jian Jiang
- Department of Materials Science & Engineering, City University of Hong Kong, Kowloon 999077, Hong Kong
| | - Yuan Liu
- School of Chemical Engineering and Technology, Sun Yat-Sen University, Zhuhai 519082, China
| | - Joseph S Francisco
- Department of Earth & Environmental Science, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Xiao Cheng Zeng
- Department of Materials Science & Engineering, City University of Hong Kong, Kowloon 999077, Hong Kong
| |
Collapse
|
31
|
Lee M, Lee SY, Kang MH, Won TK, Kang S, Kim J, Park J, Ahn DJ. Observing growth and interfacial dynamics of nanocrystalline ice in thin amorphous ice films. Nat Commun 2024; 15:908. [PMID: 38291035 PMCID: PMC10827800 DOI: 10.1038/s41467-024-45234-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2023] [Accepted: 01/16/2024] [Indexed: 02/01/2024] Open
Abstract
Ice crystals at low temperatures exhibit structural polymorphs including hexagonal ice, cubic ice, or a hetero-crystalline mixture of the two phases. Despite the significant implications of structure-dependent roles of ice, mechanisms behind the growths of each polymorph have been difficult to access quantitatively. Using in-situ cryo-electron microscopy and computational ice-dynamics simulations, we directly observe crystalline ice growth in an amorphous ice film of nanoscale thickness, which exhibits three-dimensional ice nucleation and subsequent two-dimensional ice growth. We reveal that nanoscale ice crystals exhibit polymorph-dependent growth kinetics, while hetero-crystalline ice exhibits anisotropic growth, with accelerated growth occurring at the prismatic planes. Fast-growing facets are associated with low-density interfaces that possess higher surface energy, driving tetrahedral ordering of interfacial H2O molecules and accelerating ice growth. These findings, based on nanoscale observations, improve our understanding on early stages of ice formation and mechanistic roles of the ice interface.
Collapse
Affiliation(s)
- Minyoung Lee
- School of Chemical and Biological Engineering, and Institute of Chemical Processes, Seoul National University, Seoul, 08826, Republic of Korea
- Center for Nanoparticle Research, Institute of Basic Science (IBS), Seoul, 08826, Republic of Korea
| | - Sang Yup Lee
- Department of Chemical and Biological Engineering, Korea University, Seoul, 02841, Republic of Korea
- KU-KIST Graduate school of Converging Science and Technology, Korea University, Seoul, 02841, Republic of Korea
- The w:i Interface Augmentation Center, Korea University, Seoul, 02841, Republic of Korea
| | - Min-Ho Kang
- Department of Biomedical-Chemical Engineering, The Catholic University of Korea, Bucheon-si, 14662, Republic of Korea
- Department of Biotechnology, The Catholic University of Korea, Bucheon-si, 14662, Republic of Korea
| | - Tae Kyung Won
- Department of Chemical and Biological Engineering, Korea University, Seoul, 02841, Republic of Korea
- The w:i Interface Augmentation Center, Korea University, Seoul, 02841, Republic of Korea
| | - Sungsu Kang
- School of Chemical and Biological Engineering, and Institute of Chemical Processes, Seoul National University, Seoul, 08826, Republic of Korea
- Center for Nanoparticle Research, Institute of Basic Science (IBS), Seoul, 08826, Republic of Korea
| | - Joodeok Kim
- School of Chemical and Biological Engineering, and Institute of Chemical Processes, Seoul National University, Seoul, 08826, Republic of Korea
- Center for Nanoparticle Research, Institute of Basic Science (IBS), Seoul, 08826, Republic of Korea
| | - Jungwon Park
- School of Chemical and Biological Engineering, and Institute of Chemical Processes, Seoul National University, Seoul, 08826, Republic of Korea.
- Center for Nanoparticle Research, Institute of Basic Science (IBS), Seoul, 08826, Republic of Korea.
- Institute of Engineering Research, College of Engineering, Seoul National University, Seoul, 08826, Republic of Korea.
- Advanced Institutes of Convergence Technology, Seoul National University, Suwon-si, 16229, Republic of Korea.
| | - Dong June Ahn
- Department of Chemical and Biological Engineering, Korea University, Seoul, 02841, Republic of Korea.
- KU-KIST Graduate school of Converging Science and Technology, Korea University, Seoul, 02841, Republic of Korea.
- The w:i Interface Augmentation Center, Korea University, Seoul, 02841, Republic of Korea.
| |
Collapse
|
32
|
Ishiai S, Yasuda I, Endo K, Yasuoka K. Graph-Neural-Network-Based Unsupervised Learning of the Temporal Similarity of Structural Features Observed in Molecular Dynamics Simulations. J Chem Theory Comput 2024; 20:819-831. [PMID: 38190503 DOI: 10.1021/acs.jctc.3c00995] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2024]
Abstract
Classification of molecular structures is a crucial step in molecular dynamics (MD) simulations to detect various structures and phases within systems. Molecular structures, which are commonly identified using order parameters, were recently identified using machine learning (ML), that is, the ML models acquire structural features using labeled crystals or phases via supervised learning. However, these approaches may not identify unlabeled or unknown structures, such as the imperfect crystal structures observed in nonequilibrium systems and interfaces. In this study, we proposed the use of a novel unsupervised learning framework, denoted temporal self-supervised learning (TSSL), to learn structural features and design their parameters. In TSSL, the ML models learn that the structural similarity is learned via contrastive learning based on minor short-term variations caused by perturbations in MD simulations. This learning framework is applied to a sophisticated architecture of graph neural network models that use bond angle and length data of the neighboring atoms. TSSL successfully classifies water and ice crystals based on high local ordering, and furthermore, it detects imperfect structures typical of interfaces such as the water-ice and ice-vapor interfaces.
Collapse
Affiliation(s)
- Satoki Ishiai
- Department of Mechanical Engineering, Keio University, Yokohama 223-8522, Japan
| | - Ikki Yasuda
- Department of Mechanical Engineering, Keio University, Yokohama 223-8522, Japan
| | - Katsuhiro Endo
- Department of Mechanical Engineering, Keio University, Yokohama 223-8522, Japan
- National Institute of Advanced Industrial Science and Technology (AIST), Ibaraki 305-8568, Japan
| | - Kenji Yasuoka
- Department of Mechanical Engineering, Keio University, Yokohama 223-8522, Japan
| |
Collapse
|
33
|
Luo S, Misra RP, Blankschtein D. Water Electric Field Induced Modulation of the Wetting of Hexagonal Boron Nitride: Insights from Multiscale Modeling of Many-Body Polarization. ACS NANO 2024; 18:1629-1646. [PMID: 38169482 DOI: 10.1021/acsnano.3c09811] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
Understanding the behavior of water contacting two-dimensional materials, such as hexagonal boron nitride (hBN), is important in practical applications, including seawater desalination and energy harvesting. Water, being a polar solvent, can strongly polarize the hBN surface via the electric fields that it generates. However, there is a lack of molecular-level understanding about the role of polarization effects at the hBN/water interface, including its effect on the wetting properties of water. In this study, we develop a theoretical framework that introduces an all-atomistic polarizable force field to accurately model the interactions of water molecules with hBN surfaces. The force field is then utilized to self-consistently describe the water-induced polarization of hBN using the classical Drude oscillator model, including predicting the hBN-water binding energies which are found to be in excellent agreement with diffusion Monte Carlo (DMC) predictions. By carrying out molecular dynamics (MD) simulations, we demonstrate that the polarizable force field yields a water contact angle on multilayered hBN which is in close agreement with the recent experimentally reported values. Conversely, an implicit modeling of the hBN-water polarization energy utilizing a Lennard-Jones (LJ) potential, a commonly utilized approximation in previous MD simulation studies, leads to a considerably lower water contact angle. This difference in the predicted contact angles is attributed to the significant energy-entropy compensation resulting from the incorporation of polarization effects at the hBN-water interface. Our work highlights the importance of self-consistently modeling the hBN-water polarization energy and offers insights into the wetting-related interfacial phenomena of water on polarizable materials.
Collapse
Affiliation(s)
- Shuang Luo
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Rahul Prasanna Misra
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Daniel Blankschtein
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
34
|
Khrapak SA. System Size Dependence of the Diffusion Coefficients in MD Simulations: A Simple Correction Formula for Pure Dense Fluids. J Phys Chem B 2024; 128:287-290. [PMID: 38166419 DOI: 10.1021/acs.jpcb.3c07184] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2024]
Abstract
A practical correction formula relating the self-diffusion coefficient of dense liquids from molecular dynamics (MD) simulations with periodic boundary conditions to the self-diffusion coefficient in the thermodynamic limit is discussed. This formula applies to pure dense fluids and has a very simple form D = D0(1 - γN-1/3), where D0 is the self-diffusion coefficient in the thermodynamic limit and N is the number of particles in the simulation. The numerical factor γ is dependent on the geometry of the simulation cell. Remarkably, γ ≃ 1.0 is the most popular cubic geometry. The success of this formula is supported by results from MD simulations, including very recent simulations with a "magic" simulation geometry.
Collapse
Affiliation(s)
- Sergey A Khrapak
- Joint Institute for High Temperatures, Russian Academy of Sciences, Moscow 125412, Russia
| |
Collapse
|
35
|
Hudait A. Multiscale Molecular Dynamics Simulations of Ice-Binding Proteins. Methods Mol Biol 2024; 2730:185-202. [PMID: 37943459 DOI: 10.1007/978-1-0716-3503-2_13] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2023]
Abstract
Ice-binding proteins (IBPs) are a diverse class of proteins that are essential for the survival of organisms in cold conditions. IBPs are diverse in their function and can prevent or promote ice growth and selectively bind to specific crystallographic planes of the growing ice lattice. Moreover, IBPs are widely utilized to modulate ice crystal growth and recrystallization in the food industry and as cryoprotectants to preserve biological matter. A key unresolved aspect of the mode of action is how the ice-binding sites of these proteins distinguish between ice and water and interact with multiple crystal facets of the ice. The use of molecular dynamics (MD) simulation allows us to thoroughly investigate the binding mechanism and energetics of ice-binding proteins, to complement and expand on the mechanistic understandings gained from experiments. In this chapter, we describe a series of molecular dynamics simulation methodologies to investigate the mechanism of action of ice-binding proteins. Specifically, we provide detailed instructions to set up MD simulations to study the binding and interaction of ice-binding proteins using atomistic and coarse-grained simulations.
Collapse
Affiliation(s)
- Arpa Hudait
- Department of Chemistry, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, IL, USA.
| |
Collapse
|
36
|
Ji S, Yang Z, Lei L, Galindo Torres SA, Li L. Estimation of the ice melting point in molecular dynamics simulations based on the finite-size effects. Phys Rev E 2024; 109:014108. [PMID: 38366460 DOI: 10.1103/physreve.109.014108] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2023] [Accepted: 11/22/2023] [Indexed: 02/18/2024]
Abstract
Predicting the ice melting point using molecular dynamics (MD) simulations is nontrivial due to uncertainty associated with the stochastic nature of the simulation and effect of finite domain sizes on the simulated ice-water phase transition. We developed a method based on the percolation theory to make use of the finite size effects to allow determination of a unique critical phase transition temperature as the melting point. The method involves construction of melting/freezing probability curves from multiple simulations with varying temperatures for different domain sizes. While the domain sizes affect the apparent melting/freezing probability and hence generate different curves with a wider probability distribution for a smaller size, the intersection of these curves is unique and locates the melting point. Based on MD simulations using the Tip4p/Ice water model, we tested and demonstrated the effectiveness of this method in locating the critical ice-water phase transition at a melting temperature of 268.78 K. Our analysis also showed that the apparent melting probability at this critical point is ∼0.69, not 0.5 assumed in the ad hoc method used previously. Our method, making no assumption about the system size, may provide a generic framework for analyzing phase transitions influenced by the finite size effects in MD simulations.
Collapse
Affiliation(s)
- Shuting Ji
- Department of Environmental Science and Engineering, Fudan University, Shanghai 200438, China
- School of Engineering, Westlake University, Hangzhou 310030, Zhejiang Province, China
| | - Zhenlei Yang
- School of Engineering, Westlake University, Hangzhou 310030, Zhejiang Province, China
- Key Laboratory of Coastal Environment and Resources of Zhejiang Province, Westlake University, Hangzhou 310030, Zhejiang Province, China
| | - Liang Lei
- School of Engineering, Westlake University, Hangzhou 310030, Zhejiang Province, China
- Key Laboratory of Coastal Environment and Resources of Zhejiang Province, Westlake University, Hangzhou 310030, Zhejiang Province, China
| | - Sergio Andres Galindo Torres
- School of Engineering, Westlake University, Hangzhou 310030, Zhejiang Province, China
- Key Laboratory of Coastal Environment and Resources of Zhejiang Province, Westlake University, Hangzhou 310030, Zhejiang Province, China
| | - Ling Li
- School of Engineering, Westlake University, Hangzhou 310030, Zhejiang Province, China
- Key Laboratory of Coastal Environment and Resources of Zhejiang Province, Westlake University, Hangzhou 310030, Zhejiang Province, China
| |
Collapse
|
37
|
Perrone M, Capelli R, Empereur-mot C, Hassanali A, Pavan GM. Lessons Learned from Multiobjective Automatic Optimizations of Classical Three-Site Rigid Water Models Using Microscopic and Macroscopic Target Experimental Observables. JOURNAL OF CHEMICAL AND ENGINEERING DATA 2023; 68:3228-3241. [PMID: 38115916 PMCID: PMC10726314 DOI: 10.1021/acs.jced.3c00538] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/08/2023] [Revised: 11/17/2023] [Accepted: 11/21/2023] [Indexed: 12/21/2023]
Abstract
The development of accurate water models is of primary importance for molecular simulations. Despite their intrinsic approximations, three-site rigid water models are still ubiquitously used to simulate a variety of molecular systems. Automatic optimization approaches have been recently used to iteratively refine three-site water models to fit macroscopic (average) thermodynamic properties, providing state-of-the-art three-site models that still present some deviations from the liquid water properties. Here, we show the results obtained by automatically optimizing three-site rigid water models to fit a combination of microscopic and macroscopic experimental observables. We use Swarm-CG, a multiobjective particle-swarm-optimization algorithm, for training the models to reproduce the experimental radial distribution functions of liquid water at various temperatures (rich in microscopic-level information on, e.g., the local orientation and interactions of the water molecules). We systematically analyze the agreement of these models with experimental observables and the effect of adding macroscopic information to the training set. Our results demonstrate how adding microscopic-rich information in the training of water models allows one to achieve state-of-the-art accuracy in an efficient way. Limitations in the approach and in the approximated description of water in these three-site models are also discussed, providing a demonstrative case useful for the optimization of approximated molecular models, in general.
Collapse
Affiliation(s)
- Mattia Perrone
- Department
of Applied Science and Technology, Politecnico
di Torino, Corso Duca degli Abruzzi 24, Torino I-10129, Italy
| | - Riccardo Capelli
- Department
of Biosciences, Università degli
Studi di Milano, Via Celoria 26, Milano I-20133, Italy
| | - Charly Empereur-mot
- Department
of Innovative Technologies, University of Applied Sciences and Arts
of Southern Switzerland, Polo Universitario
Lugano, Campus Est, Via
la Santa 1, Lugano-Viganello CH-6962, Switzerland
| | - Ali Hassanali
- The
Abdus Salam International Center for Theoretical Physics, Strada Costiera 11, Trieste 34151, Italy
| | - Giovanni M. Pavan
- Department
of Applied Science and Technology, Politecnico
di Torino, Corso Duca degli Abruzzi 24, Torino I-10129, Italy
- Department
of Innovative Technologies, University of Applied Sciences and Arts
of Southern Switzerland, Polo Universitario
Lugano, Campus Est, Via
la Santa 1, Lugano-Viganello CH-6962, Switzerland
| |
Collapse
|
38
|
Joliat J, Picaud S, Jedlovszky P. Competitive Adsorption of Trace Gases on Ice at Tropospheric Temperatures: A Grand Canonical Monte Carlo Simulation Study. J Phys Chem A 2023; 127:10223-10232. [PMID: 38000079 DOI: 10.1021/acs.jpca.3c04789] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2023]
Abstract
The coadsorption of two atmospheric trace gases on ice is characterized by using, for the first time, grand canonical Monte Carlo (GCMC) simulations performed in conditions similar to those of the corresponding experiments. Adsorption isotherms are simulated at tropospheric temperatures by considering two different gas mixtures of 1-butanol and acetic acid molecules, and selectivity of the ice surface with respect to these species is interpreted at the molecular scale as resulting from a competition process between these molecules for being adsorbed at the ice surface. It is thus shown that the trapping of acetic acid molecules on ice is always favored with respect to that of 1-butanol at low pressures, corresponding to low coverage of the surface, whereas the adsorption of the acid species is significantly modified by the presence of the alcohol molecules in the saturated portion of the adsorption isotherm, in accordance with the experimental observations. The present GCMC simulations thus confirm that competitive adsorption effects have to be taken into consideration in real situations when gas mixtures present in the troposphere interact with the surface of ice particles.
Collapse
Affiliation(s)
- Julien Joliat
- Institut UTINAM─UMR 6213, CNRS/Université de Franche-Comté, 25000 Besançon, France
| | - Sylvain Picaud
- Institut UTINAM─UMR 6213, CNRS/Université de Franche-Comté, 25000 Besançon, France
| | - Pál Jedlovszky
- Department of Chemistry, Eszterházy Károly Catholic University, Leányka U. 6, H-3300 Eger, Hungary
| |
Collapse
|
39
|
Li K, Chen B, Li M, Jiang L, Song Y, Yang M. Facilitation of Hydrate Dissociation and Structural Evolution by Major Marine Anions under Static Electric Fields. J Phys Chem B 2023; 127:10447-10457. [PMID: 37991934 DOI: 10.1021/acs.jpcb.3c06012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2023]
Abstract
Electric fields have been proven to be capable of significantly affecting the equilibrium state of hydrates. In this study, the thermodynamic properties and structural changes of methane hydrate (MH) in various anion solutions in an electric field at 0.7 V/nm were investigated by molecular dynamics simulations. The presence of anions significantly enhances the instability of methane hydrates under electric fields, leading to a staged dissociation process. First, the anions coexist with MH to form a temporary metastable structure under the action of an electric field. Then, the migration of anions causes the dissociation of nearby hydrates and the formation of flow channels in the hydrate layer, which leads to the complete dissociation of MH after a period. The promotive effects of F-, Br-, I-, and Cl- ions were close, while SO42- was relatively weak. The anions are still in hydration shells in the MH phase, but the structure of the hydration shells differs slightly from that in solution (the coordination numbers of I- and SO42- ions increased). The migration resistances of multiple anions to cross the surface of the hydrate layer are similar. However, inside the hydrate phase, the anions with a larger radius have a higher migration resistance. It is difficult for SO42- ions to migrate inside the hydrate phase, and they tend to form a metastable structure on the hydrate surface. Combining our previous studies, SrCl2 solution has the best hydrate promotion under an electric field environment.
Collapse
Affiliation(s)
- Kehan Li
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian 116024, China
| | - Bingbing Chen
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian 116024, China
| | - Mingjun Li
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian 116024, China
| | - Lanlan Jiang
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian 116024, China
| | - Yongchen Song
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian 116024, China
| | - Mingjun Yang
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian 116024, China
| |
Collapse
|
40
|
Fu Z, Jing X, Zhou L, Luo Q, Zhang P. Molecular simulation of imperfect structure I CO 2 hydrate growth in brine. Acta Crystallogr C Struct Chem 2023; 79:513-519. [PMID: 38019214 DOI: 10.1107/s2053229623010148] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2023] [Accepted: 11/22/2023] [Indexed: 11/30/2023] Open
Abstract
In order to investigate the viability of carbon dioxide (CO2) storage in seawater, molecular dynamics techniques were employed to study the dynamic evolution of CO2 hydrate in saline water. The simulation was conducted under specific conditions: a temperature of 275 K, a pressure of 10 MPa and a simulated marine environment achieved using a 3.4 wt% sodium chloride (NaCl) solution. The total simulation time was 1000 ns. The results of the simulation indicate that the pre-existence of CO2 hydrate crystals as seeds leads to rapid growth of CO2 hydrate. However, analysis of the F3 and F4 order parameters reveals that the hydrate does not meet the standard values of the perfect structure I (sI) type, confirming the existence of an imperfect structure during the simulation. Additionally, the changes in the number of different phase states of water molecules during the hydrate growth process shows that there are always some liquid water molecules, which means some water molecules fail to form solid water cages. Further investigation suggests that the presence of Na+ and Cl- hampers the hydrogen bonds between water molecules, resulting in incomplete cage structures. By analyzing the density variations in the system, it is observed that CO2 hydrate, with a density of around 1.133 g cm-3, forms rapidly, surpassing the average density of seawater. This density increase facilitates the efficient and swift containment of CO2 on the seabed, thereby supporting the feasibility of the CO2 storage theory.
Collapse
Affiliation(s)
- Ziyi Fu
- Research Institute of Natural Gas Technology, PetroChina Southwest Oil and Gasfield Company, Chengdu, Sichuan 610213, People's Republic of China
| | - Xianwu Jing
- Research Institute of Natural Gas Technology, PetroChina Southwest Oil and Gasfield Company, Chengdu, Sichuan 610213, People's Republic of China
| | - Li Zhou
- Research Institute of Natural Gas Technology, PetroChina Southwest Oil and Gasfield Company, Chengdu, Sichuan 610213, People's Republic of China
| | - Qin Luo
- Research Institute of Natural Gas Technology, PetroChina Southwest Oil and Gasfield Company, Chengdu, Sichuan 610213, People's Republic of China
| | - Pengfei Zhang
- Research Institute of Natural Gas Technology, PetroChina Southwest Oil and Gasfield Company, Chengdu, Sichuan 610213, People's Republic of China
| |
Collapse
|
41
|
Gattan HS, Mahmoud Alawi M, Bajrai LH, Alandijany TA, Alsaady IM, El-Daly MM, Dwivedi VD, Azhar EI. A Multifaceted Computational Approach to Understanding the MERS-CoV Main Protease and Brown Algae Compounds' Interaction. Mar Drugs 2023; 21:626. [PMID: 38132947 PMCID: PMC10744363 DOI: 10.3390/md21120626] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2023] [Revised: 11/22/2023] [Accepted: 11/28/2023] [Indexed: 12/23/2023] Open
Abstract
Middle East Respiratory Syndrome (MERS) is a viral respiratory disease caused b a special type of coronavirus called MERS-CoV. In the search for effective substances against the MERS-CoV main protease, we looked into compounds from brown algae, known for their medicinal benefits. From a set of 1212 such compounds, our computer-based screening highlighted four-CMNPD27819, CMNPD1843, CMNPD4184, and CMNPD3156. These showed good potential in how they might attach to the MERS-CoV protease, comparable to a known inhibitor. We confirmed these results with multiple computer tests. Studies on the dynamics and steadiness of these compounds with the MERS-CoV protease were performed using molecular dynamics (MD) simulations. Metrics like RMSD and RMSF showed their stability. We also studied how these compounds and the protease interact in detail. An analysis technique, PCA, showed changes in atomic positions over time. Overall, our computer studies suggest brown algae compounds could be valuable in fighting MERS. However, experimental validation is needed to prove their real-world effectiveness.
Collapse
Affiliation(s)
- Hattan S. Gattan
- Special Infectious Agents Unit-BSL3, King Fahd Medical Research Center, King Abdulaziz University, Jeddah 21362, Saudi Arabia; (H.S.G.); (M.M.A.); (L.H.B.); (T.A.A.); (M.M.E.-D.)
- Department of Medical Laboratory Sciences, Faculty of Applied Medical Sciences, King Abdulaziz University, Jeddah 21362, Saudi Arabia
| | - Maha Mahmoud Alawi
- Special Infectious Agents Unit-BSL3, King Fahd Medical Research Center, King Abdulaziz University, Jeddah 21362, Saudi Arabia; (H.S.G.); (M.M.A.); (L.H.B.); (T.A.A.); (M.M.E.-D.)
- Department of Clinical Microbiology and Immunology, Faculty of Medicine, King Abdulaziz University, Jeddah 21589, Saudi Arabia
- Infection Control & Environmental Health Unit, King Abdulaziz University Hospital, King Abdulaziz University, Jeddah 21589, Saudi Arabia
| | - Leena H. Bajrai
- Special Infectious Agents Unit-BSL3, King Fahd Medical Research Center, King Abdulaziz University, Jeddah 21362, Saudi Arabia; (H.S.G.); (M.M.A.); (L.H.B.); (T.A.A.); (M.M.E.-D.)
- Biochemistry Department, Faculty of Sciences, King Abdulaziz University, Jeddah 21362, Saudi Arabia
| | - Thamir A. Alandijany
- Special Infectious Agents Unit-BSL3, King Fahd Medical Research Center, King Abdulaziz University, Jeddah 21362, Saudi Arabia; (H.S.G.); (M.M.A.); (L.H.B.); (T.A.A.); (M.M.E.-D.)
- Department of Medical Laboratory Sciences, Faculty of Applied Medical Sciences, King Abdulaziz University, Jeddah 21362, Saudi Arabia
| | - Isra M. Alsaady
- Special Infectious Agents Unit-BSL3, King Fahd Medical Research Center, King Abdulaziz University, Jeddah 21362, Saudi Arabia; (H.S.G.); (M.M.A.); (L.H.B.); (T.A.A.); (M.M.E.-D.)
- Department of Medical Laboratory Sciences, Faculty of Applied Medical Sciences, King Abdulaziz University, Jeddah 21362, Saudi Arabia
| | - Mai M. El-Daly
- Special Infectious Agents Unit-BSL3, King Fahd Medical Research Center, King Abdulaziz University, Jeddah 21362, Saudi Arabia; (H.S.G.); (M.M.A.); (L.H.B.); (T.A.A.); (M.M.E.-D.)
- Department of Medical Laboratory Sciences, Faculty of Applied Medical Sciences, King Abdulaziz University, Jeddah 21362, Saudi Arabia
| | - Vivek Dhar Dwivedi
- Center for Global Health Research, Saveetha Institute of Medical and Technical Sciences, Saveetha Medical College and Hospitals, Saveetha University, Chennai 605102, India
- Bioinformatics Research Division, Quanta Calculus, Greater Noida 201310, India
| | - Esam I. Azhar
- Special Infectious Agents Unit-BSL3, King Fahd Medical Research Center, King Abdulaziz University, Jeddah 21362, Saudi Arabia; (H.S.G.); (M.M.A.); (L.H.B.); (T.A.A.); (M.M.E.-D.)
- Department of Medical Laboratory Sciences, Faculty of Applied Medical Sciences, King Abdulaziz University, Jeddah 21362, Saudi Arabia
| |
Collapse
|
42
|
Wang L, Kusalik PG. Understanding why constant energy or constant temperature may affect nucleation behavior in MD simulations: A study of gas hydrate nucleation. J Chem Phys 2023; 159:184501. [PMID: 37947514 DOI: 10.1063/5.0169669] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2023] [Accepted: 10/18/2023] [Indexed: 11/12/2023] Open
Abstract
Molecular dynamics simulations have been widely used in exploring the nucleation behavior of many systems, including gas hydrates. Gas hydrates are ice-like solids in which gas molecules are trapped in water cages. During hydrate formation, a considerable amount of heat is released, and previous work has reported that the choice of temperature control scheme may affect the behavior of hydrate formation. The origins of this effect have remained an open question. To address this question, extensive NVE simulations and thermostatted (NPT and NVT) simulations with different temperature coupling strengths have been performed and compared for systems where a water nanodroplet is immersed in a H2S liquid. Detailed analysis of the hydrate structures and their mechanisms of formation has been carried out. Slower nucleation rates in NVE simulations in comparison to NPT simulations have been observed in agreement with previous studies. Probability distributions for various temperature measures along with their spatial distributions have been examined. Interestingly, a comparison of these temperature distributions reveals a small yet noticeable difference in the widths of the distributions for water. The somewhat reduced fluctuations in the temperature for the water species in the NVE simulations appear to be responsible for reducing the hydrate nucleation rate. We further conjecture that the NVE-impeded nucleation rate may be the result of the finite size of the surroundings (here the liquid H2S portion of the system). Additionally, a local spatial temperature gradient arising from the heat released during hydrate formation could not be detected.
Collapse
Affiliation(s)
- Lei Wang
- Department of Chemistry, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Peter G Kusalik
- Department of Chemistry, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| |
Collapse
|
43
|
Song H, Dong H, Dong W, Luo Y. Atomic-Level Insights into Hollow Silica-Based Materials for Drug Delivery: Effects of Wettability and Porosity. ACS Biomater Sci Eng 2023; 9:6156-6164. [PMID: 37831542 DOI: 10.1021/acsbiomaterials.3c01063] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2023]
Abstract
Experimental evidence has demonstrated that the drug carrier capacity can be significantly enhanced through the use of hollow silica particles. Nevertheless, the effects of varying functional drug carrier surfaces and porous structures remain ambiguous. This study employs molecular dynamics simulations to examine the effects of varying the surface wettability, pore size, and flow velocity on the transfer process. The different levels of wettability of the silica surface with the coarse-grained water model is illustrated by adjusted interaction parameters. The effect of wettability is investigated. With weak interactions, the flow molecules form a nanodroplet to transfer through the porous structure. A strong interaction will lead to molecules flowing as a liquid film to transfer through the structure. Interestingly, the "contradiction effect" is observed when the flow molecules fail to penetrate the porous structure with weak interactions, during which surface tension dominates their flow behavior. Moreover, different porous structures are considered. The flow behaviors are divided into three processes: (1) fast flowing, (2) transient point, and (3) penetration flowing. Furthermore, the concept of surface molecules is defined to quantitatively measure the effect of porosity. A recommended contact angle is proposed. The results will pave the way for more carrier structures in medical engineering.
Collapse
Affiliation(s)
- Haoxin Song
- Department of Pharmacy, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an 710061, China
| | - Haiyan Dong
- Department of Pharmacy, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an 710061, China
| | - Weihua Dong
- Department of Pharmacy, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an 710061, China
| | - Yu Luo
- Department of Pharmacy, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an 710061, China
| |
Collapse
|
44
|
Li L, Wang X, Yan Y, Francisco JS, Zhang J, Zeng XC, Zhong J. Resolving Temperature-Dependent Hydrate Nucleation Pathway: The Role of "Transition Layer". J Am Chem Soc 2023; 145:24166-24174. [PMID: 37874937 DOI: 10.1021/jacs.3c08246] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2023]
Abstract
Understanding the nucleation of natural gas hydrate (NGH) at different conditions has important implications to NGH recovery and other industrial applications, such as gas storage and separation. Herein, vast numbers of hydrate nucleation events are traced via molecular dynamics (MD) simulations at different degrees of supercooling (or driving forces). Specifically, to precisely characterize a hydrate nucleus from an aqueous system during the MD simulation, we develop an evolutionary order parameter (OP) to recognize the nucleus size and shape. Subsequently, the free energy landscapes of hydrate during nucleation are explored by using the newly developed OP. The results suggest that at 270 K (or 0.92 Tm supercooling, where Tm is the melting point), the near-rounded nucleus prevails during the nucleation, as described from the classical nucleation theory. In contrast, at relatively strong driving forces of 0.85 and 0.88 Tm, nonclassical nucleation events arise. Specifically, the pathway toward an elongated nucleus becomes as important as the pathway toward a near-rounded nucleus. To explain the distinct nucleation phenomena at different supercoolings, a notion of a "transition layer" (or liquid-blob-like layer) is proposed. Here, the transition layer is to describe the interfacial region between the nucleus and aqueous solution, and this layer entails two functionalities: (1) it tends to retain CH4 depending on the degrees of supercooling and (2) it facilitates collision among CH4, which thus promote the incorporation of CH4 into nucleus. Our simulation indicates that compared to the near-rounded nucleus, the transition layer surrounding the elongated nucleus is more evident with the higher collision rate among CH4 molecules. As such, the transition layer tends to promote the elongated nucleus pathway, while offsetting the cost of larger surface free energy associated with the elongated nucleus. At 0.92 Tm, however, the transition layer gradually disappears, and classical nucleation events dominate. Overall, the notion of "transition layer" offers deeper insight into the NGH nucleation at different degrees of supercooling and could be extended to describe other types of hydrate nucleation.
Collapse
Affiliation(s)
- Liwen Li
- School of Petroleum Engineering and School of Materials Science and Engineering, China University of Petroleum (East China), Qingdao, Shandong 266580, China
- Department of Materials Science & Engineering, City University of Hong Kong, Hong Kong 999077, China
| | - Xiao Wang
- School of Petroleum Engineering and School of Materials Science and Engineering, China University of Petroleum (East China), Qingdao, Shandong 266580, China
| | - Youguo Yan
- School of Petroleum Engineering and School of Materials Science and Engineering, China University of Petroleum (East China), Qingdao, Shandong 266580, China
| | - Joseph S Francisco
- Department of Earth and Environmental Science, University of Pennsylvania, Philadelphia, Pennsylvania 19104-6316, United States
| | - Jun Zhang
- School of Petroleum Engineering and School of Materials Science and Engineering, China University of Petroleum (East China), Qingdao, Shandong 266580, China
| | - Xiao Cheng Zeng
- Department of Materials Science & Engineering, City University of Hong Kong, Hong Kong 999077, China
- Department of Chemistry, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, United States
| | - Jie Zhong
- School of Petroleum Engineering and School of Materials Science and Engineering, China University of Petroleum (East China), Qingdao, Shandong 266580, China
| |
Collapse
|
45
|
Yamaguchi S, Takayama T, Otosu T. Appraisal of TIP4P-type models at water surface. J Chem Phys 2023; 159:171101. [PMID: 37909448 DOI: 10.1063/5.0171999] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2023] [Accepted: 10/03/2023] [Indexed: 11/03/2023] Open
Abstract
In view of the current situation in which non-polarizable rigid water models have been scarcely examined against surface-specific properties, we appraise TIP4P-type models at the liquid water surface on the basis of heterodyne-detected sum frequency generation (HD-SFG) spectroscopy. We find in the HD-SFG spectrum of the water surface that the peak frequency of the hydrogen-bonded OH band, the half width at half maximum of the hydrogen-bonded OH band, and the full width at half maximum of the free OH band are best reproduced by TIP4P, TIP4P/Ew, and TIP4P/Ice, respectively, whereas it is already well known that TIP4P/2005 best reproduces the surface tension. These TIP4P-type models perform better at the water surface in terms of the present appraisal items than some polarizable models in the literature.
Collapse
Affiliation(s)
- Shoichi Yamaguchi
- Department of Applied Chemistry, Graduate School of Science and Engineering, Saitama University, 255 Shimo-Okubo, Sakura, Saitama 338-8570, Japan
| | - Tetsuyuki Takayama
- Department of Applied Chemistry, Graduate School of Science and Engineering, Saitama University, 255 Shimo-Okubo, Sakura, Saitama 338-8570, Japan
| | - Takuhiro Otosu
- Department of Applied Chemistry, Graduate School of Science and Engineering, Saitama University, 255 Shimo-Okubo, Sakura, Saitama 338-8570, Japan
| |
Collapse
|
46
|
Muniz MC, Car R, Panagiotopoulos AZ. Neural Network Water Model Based on the MB-Pol Many-Body Potential. J Phys Chem B 2023; 127:9165-9171. [PMID: 37824703 DOI: 10.1021/acs.jpcb.3c04629] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/14/2023]
Abstract
The MB-pol many-body potential accurately predicts many properties of water, including cluster, liquid phase, and vapor-liquid equilibrium properties, but its high computational cost can make applying it in large-scale simulations quite challenging. In order to address this limitation, we developed a "deep potential" neural network (DPMD) model based on the MB-pol potential for water. We find that a DPMD model trained on mostly liquid configurations yields a good description of the bulk liquid phase but severely underpredicts vapor-liquid coexistence densities. By contrast, adding cluster configurations to the neural network training set leads to a good agreement for the vapor coexistence densities. Liquid phase densities under supercooled conditions are also represented well, even though they were not included in the training set. These results confirm that neural network models can combine accuracy and transferability if sufficient attention is given to the construction of a representative training set for the target system.
Collapse
Affiliation(s)
- Maria Carolina Muniz
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Roberto Car
- Department of Chemistry, Department of Physics, Program in Applied and Computational Mathematics, and Princeton Materials Institute, Princeton University, Princeton, New Jersey 08544, United States
| | | |
Collapse
|
47
|
Kim YD, Jung WH, Ahn DJ, Lim DK. Self-Assembled Nanostructures of Homo-Oligopeptide as a Potent Ice Growth Inhibitor. NANO LETTERS 2023; 23:9500-9507. [PMID: 37843112 DOI: 10.1021/acs.nanolett.3c03059] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/17/2023]
Abstract
This study reports the formation of self-assembled nanostructures with homo-oligopeptides consisting of amino acids (i.e., alanine, threonine, valine, and tyrosine), the resulting morphologies (i.e., spherical shape, layered structure, and wire structure) in aqueous solution, and their potential as ice growth inhibitors. Among the homo-oligopeptides investigated, an alanine homo-oligopeptide (n = 5) with a spherical nanostructure showed the highest ice recrystallization inhibition (IRI) activity without showing a burst ice growth property and with low ice nucleation activity. The presence of nanoscale self-assembled structures in the solution showed superior IRI activity compared to an amino acid monomer because of the higher binding affinity of structures on the growing ice crystal plane. Simulation results revealed that the presence of nanostructures induced a significant inhibition of ice growth and increased lifetime of hydrogen bonding compared with unassembled homo-oligopeptide. These results envision extraordinary performance for self-assembled nanostructures as a desirable and potent ice growth inhibitor.
Collapse
Affiliation(s)
- Yong Duk Kim
- KU-KIST Graduate School of Converging Science and Technology, Korea University, 145 Anam-ro, Seongbuk-gu, Seoul 02841, Republic of Korea
| | - Woo Hyuk Jung
- Department of Chemical and Biological Engineering, Korea University, 145 Anam-ro, Seongbuk-gu, Seoul 02841, Republic of Korea
| | - Dong June Ahn
- KU-KIST Graduate School of Converging Science and Technology, Korea University, 145 Anam-ro, Seongbuk-gu, Seoul 02841, Republic of Korea
- Department of Chemical and Biological Engineering, Korea University, 145 Anam-ro, Seongbuk-gu, Seoul 02841, Republic of Korea
| | - Dong-Kwon Lim
- KU-KIST Graduate School of Converging Science and Technology, Korea University, 145 Anam-ro, Seongbuk-gu, Seoul 02841, Republic of Korea
- Department of Integrative Energy Engineering, Korea University, 145 Anam-ro, Seongbuk-gu, Seoul 02841, Republic of Korea
- Brain Science Institute, Korea Institute of Science and Technology (KIST), 5, Hwarang-ro 14-gil, Seongbuk-gu, Seoul 02792, Republic of Korea
| |
Collapse
|
48
|
Hao X, Li C, Meng Q, Sun J, Huang L, Bu Q, Li C. Molecular Dynamics Simulation of the Three-Phase Equilibrium Line of CO 2 Hydrate with OPC Water Model. ACS OMEGA 2023; 8:39847-39854. [PMID: 37901483 PMCID: PMC10601413 DOI: 10.1021/acsomega.3c05673] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/16/2023] [Accepted: 09/26/2023] [Indexed: 10/31/2023]
Abstract
The three-phase coexistence line of the CO2 hydrate was determined using molecular dynamics (MD) simulations. By using the classical and modified Lorentz-Berthelot (LB) parameters, the simulations were carried out at 10 different pressures from 3 to 500 MPa. For the OPC water model, simulations with the classic and the modified LB parameters both showed negative deviations from the experimental values. For the TIP4P/Ice water model, good agreement with experimental equilibrium data can be achieved when the LB parameter is adjusted based on the solubility of CO2 in water. Our results also show that the influence of the water model on the equilibrium prediction is much larger than the CO2 model. Current simulations indicated that the H2O-H2O and H2O-CO2 cross-interactions' parameters might contribute equally to the accurate prediction of T3. According to our simulations, the prediction of T3 values showed relatively higher accuracy while using the combination of TIP4P/Ice water and EPM2 CO2 with modified LB parameter. Furthermore, varied χ values are recommended for accurate T3 estimation over a wide pressure range. The knowledge obtained in this study will be helpful for further accurate MD simulation of the process of CO2/CH4 replacement.
Collapse
Affiliation(s)
- Xiluo Hao
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Chengfeng Li
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Qingguo Meng
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Jianye Sun
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Li Huang
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Qingtao Bu
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Congying Li
- Center
of Deep Sea Research, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| |
Collapse
|
49
|
Zhu X, Vandamme M, Jiang Z, Brochard L. Molecular simulation of the confined crystallization of ice in cement nanopore. J Chem Phys 2023; 159:154704. [PMID: 37850696 DOI: 10.1063/5.0169783] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2023] [Accepted: 09/26/2023] [Indexed: 10/19/2023] Open
Abstract
Freezing of water under nanoconfinement exhibits physical peculiarities with respect to the bulk water. However, experimental observations are extremely challenging at this scale, which limits our understanding of the effect of confinement on water properties upon freezing. In this study, we use molecular dynamic simulations to investigate how confinement affects the kinetics of growth of ice and the thermodynamic equilibrium of ice-liquid coexistence. TIP4P/Ice water model and CSH-FF model were applied to simulate ice crystallization in a confined cement system at temperatures down to 220 K. We adapted an interface detection algorithm and reparameterized the CHILL/CHILL+ algorithm to capture ice growth. The confinement leads to a shift of the maximum growth rate of ice to a higher temperature than for bulk water. Both the confinement and surface impurities contribute to slowing down the ice growth. For the ice-liquid coexistence at equilibrium, we derive a formulation of Thomson's equation adapted to statistical physics quantities accessible by molecular simulation, and we show that this adapted equation predicts accurately the melting line of bulk and confined ice Ih as a function of pressure. The confinement decreases systematically the melting temperature of ice of about 5 K compared with bulk ice Ih. A premelted water film about 1 nm thick is observed between the solid wall and ice, and its thickness is found to decrease continuously as temperature is lowered. We note that the surface impurities are key to the formation of the premelted water nanofilm when the temperature is lower than 250 K.
Collapse
Affiliation(s)
- Xinping Zhu
- Key Laboratory of Advanced Civil Engineering Materials of Ministry of Education, School of Materials Science and Engineering, Tongji University, Shanghai 201804, China
- Navier, Ecole des Ponts, Univ. Gustave Eiffel, CNRS, Marne-la-Vallée, France
| | - Matthieu Vandamme
- Navier, Ecole des Ponts, Univ. Gustave Eiffel, CNRS, Marne-la-Vallée, France
| | - Zhengwu Jiang
- Key Laboratory of Advanced Civil Engineering Materials of Ministry of Education, School of Materials Science and Engineering, Tongji University, Shanghai 201804, China
| | - Laurent Brochard
- Navier, Ecole des Ponts, Univ. Gustave Eiffel, CNRS, Marne-la-Vallée, France
| |
Collapse
|
50
|
Lupi L, Gallo P. Glassy dynamics of water in TIP4P/Ice aqueous solutions of trehalose in comparison with the bulk phase. J Chem Phys 2023; 159:154504. [PMID: 37850697 DOI: 10.1063/5.0168933] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2023] [Accepted: 09/25/2023] [Indexed: 10/19/2023] Open
Abstract
We perform molecular dynamics simulations of TIP4P/Ice water in solution with trehalose for 3.65 and 18.57 wt. % concentrations and of bulk TIP4P/Ice water at ambient pressure, to characterize the structure and dynamics of water in a sugar aqueous solution in the supercooled region. We find here that TIP4P/Ice water in solution with trehalose molecules follows the Mode Coupling Theory and undergoes a fragile to strong transition up to the highest concentration investigated, similar to the bulk. Moreover, we perform a Mode Coupling Theory test, showing that the Time Temperature Superposition principle holds for both bulk TIP4P/Ice water and for TIP4P/Ice water in the solutions and we calculate the exponents of the theory. The direct comparison of the dynamical results for bulk water and water in the solutions shows upon cooling along the isobar a fastening of water dynamics for lower temperatures, T < 240 K. We found that the counter-intuitive behavior for the low temperature solutions can be explained with the diffusion anomaly of water leading us to the conclusion that the fastening observed below T = 240 K in water dynamics is only fictitious, due to the fact that the density of water molecules in the solutions is higher than the density of the bulk at the same temperature and pressure. This result should be taken into account in experimental investigations which are often carried out at constant pressure.
Collapse
Affiliation(s)
- Laura Lupi
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Rome, Italy
| | - Paola Gallo
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Rome, Italy
| |
Collapse
|