1
|
Yagasaki T, Matubayasi N. High Antifouling Performance of Weakly Hydrophilic Polymer Brushes: A Molecular Dynamics Study. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2024; 40:15046-15058. [PMID: 39004900 DOI: 10.1021/acs.langmuir.4c01365] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/16/2024]
Abstract
The antifouling performance of polymer brushes usually improves with increasing hydrophilicity of the grafted polymer. However, in some cases, less hydrophilic polymers show comparable or better antifouling performance than do more hydrophilic polymers. We investigate the mechanism of this anomalous behavior using molecular dynamics (MD) simulations of coarse-grained (CG) models of weakly and strongly hydrophilic polymers. The antifouling performance is evaluated from the potential of mean force of a model protein. The strongly hydrophilic polymer exhibits a better antifouling performance than the weakly hydrophilic polymer when the substrate of the polymer brush is repulsive. However, when the substrate is sufficiently attractive, the weakly hydrophilic polymer brush becomes more effective than the strongly hydrophilic brush in a certain range of grafting density. This is because the weakly hydrophilic polymer chains form a tightly packed layer that prevents the adsorbate molecule from contacting the substrate. We also perform all-atom (AA) MD simulations for several standard polymers to examine the correspondence with the CG polymer models. The weakly hydrophilic CG polymer is found to be similar to poly[N-(2-hydroxypropyl)methacrylamide] and poly(2-hydroxyethyl methacrylate), both of which have a hydroxyl group in a monomer unit. The strongly hydrophilic CG polymer resembles zwitterionic poly(carboxybetaine methacrylate). A discussion referring to the adsorption free energies of proteins on surfaces calculated in previous AA MD studies suggests that the higher antifouling performance of less hydrophilic polymer brushes can be realized for various combinations of protein and surface.
Collapse
Affiliation(s)
- Takuma Yagasaki
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka 560-8531, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka 560-8531, Japan
| |
Collapse
|
2
|
Matsumoto M, Yagasaki T, Tanaka H. GenIce-core: Efficient algorithm for generation of hydrogen-disordered ice structures. J Chem Phys 2024; 160:094101. [PMID: 38426513 DOI: 10.1063/5.0198056] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2024] [Accepted: 02/09/2024] [Indexed: 03/02/2024] Open
Abstract
Ice is different from ordinary crystals because it contains randomness, which means that statistical treatment based on ensemble averaging is essential. Ice structures are constrained by topological rules known as the ice rules, which give them unique anomalous properties. These properties become more apparent when the system size is large. For this reason, there is a need to produce a large number of sufficiently large crystals that are homogeneously random and satisfy the ice rules. We have developed an algorithm to quickly generate ice structures containing ions and defects. This algorithm is provided as an independent software module that can be incorporated into crystal structure generation software. By doing so, it becomes possible to simulate ice crystals on a previously impossible scale.
Collapse
Affiliation(s)
- Masakazu Matsumoto
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Takuma Yagasaki
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Osaka 560-8531, Japan
| | - Hideki Tanaka
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
- Toyota Physical and Chemical Research Institute, Nagakute 480-1192, Japan
| |
Collapse
|
3
|
Yagasaki T, Matubayasi N. Molecular Dynamics Study of the Antifouling Mechanism of Hydrophilic Polymer Brushes. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2023; 39:13158-13168. [PMID: 37672759 DOI: 10.1021/acs.langmuir.3c01552] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/08/2023]
Abstract
We perform all-atom molecular dynamics simulations of the adsorption of amino acid side-chain analogues on polymer brushes. The analogues examined are nonpolar isobutane, polar propionamide, negatively charged propionate ion, and positively charged butylammonium ion. The polymer brushes consist of a sheet of graphene and strongly hydrophilic poly(carboxybetaine methacrylate) (PCBMA) or weakly hydrophilic poly(2-hydroxyethyl methacrylate) (PHEMA). The effective interactions between isobutane and polymer chains are repulsive for PCBMA and attractive for PHEMA. Gibbs energy decomposition analysis shows that this is due to the abundance of water in the PCBMA brush, which increases the steric repulsion and decreases the Lennard-Jones attraction. The affinity of the hydrophilic analogues is low for both PCBMA and PHEMA chains, but the balance between the components of the Gibbs energy is different for the two polymers. The simulations are performed at several θ, where θ is the degree of overlap of polymer chains. The antifouling performance against the neutral analogues is better for PCBMA than for PHEMA in the low and high θ regimes. However, in the middle θ regime, the antifouling performance of PHEMA is close to or better than that of PCBMA. This is attributed to the formation of a dense layer of PHEMA on the graphene surface that inhibits direct adsorption of analogue molecules on graphene. The charged analogues do not bind to either the PHEMA or PCBMA brush irrespective of θ.
Collapse
Affiliation(s)
- Takuma Yagasaki
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka 560-8531, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka 560-8531, Japan
| |
Collapse
|
4
|
Phan A, Stamatakis M, Koh CA, Striolo A. Microscopic insights on clathrate hydrate growth from non-equilibrium molecular dynamics simulations. J Colloid Interface Sci 2023; 649:185-193. [PMID: 37348338 DOI: 10.1016/j.jcis.2023.06.032] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2023] [Revised: 05/03/2023] [Accepted: 06/06/2023] [Indexed: 06/24/2023]
Abstract
Clathrate hydrates form and grow at interfaces. Understanding the relevant molecular processes is crucial for developing hydrate-based technologies. Many computational studies focus on hydrate growth within the aqueous phase using the 'direct coexistence method', which is limited in its ability to investigate hydrate film growth at hydrocarbon-water interfaces. To overcome this shortcoming, a new simulation setup is presented here, which allows us to study the growth of a methane hydrate nucleus in a system where oil-water, hydrate-water, and hydrate-oil interfaces are all simultaneously present, thereby mimicking experimental setups. Using this setup, hydrate growth is studied here under the influence of two additives, a polyvinylcaprolactam oligomer and sodium dodecyl sulfate, at varying concentrations. Our results confirm that hydrate films grow along the oil-water interface, in general agreement with visual experimental observations; growth, albeit slower, also occurs at the hydrate-water interface, the interface most often interrogated via simulations. The results obtained demonstrate that the additives present within curved interfaces control the solubility of methane in the aqueous phase, which correlates with hydrate growth rate. Building on our simulation insights, we suggest that by combining data for the potential of mean force profile for methane transport across the oil-water interface and for the average free energy required to perturb a flat interface, it is possible to predict the performance of additives used to control hydrate growth. These insights could be helpful to achieve optimal methane storage in hydrates, one of many applications which are attracting significant fundamental and applied interests.
Collapse
Affiliation(s)
- Anh Phan
- School of Chemistry and Chemical Engineering, Faculty of Engineering and Physical Sciences, University of Surrey, Guildford, Surrey GU2 7XH, UK.
| | - Michail Stamatakis
- Department of Chemical Engineering, University College London, London WC1E 7JE, UK
| | - Carolyn A Koh
- Center for Hydrate Research, Chemical & Biological Engineering Department, Colorado School of Mines, Golden, CO 80401, United States
| | - Alberto Striolo
- Department of Chemical Engineering, University College London, London WC1E 7JE, UK; School of Chemical, Biological and Materials Engineering, University of Oklahoma, Norman, OK 73019, United States.
| |
Collapse
|
5
|
Alanine rich amphiphilic peptides as green substitutes for hydrate inhibitors: a molecular simulation study. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.121008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
|
6
|
Yagasaki T, Matubayasi N. Molecular dynamics study of the interactions between a hydrophilic polymer brush on graphene and amino acid side chain analogues in water. Phys Chem Chem Phys 2022; 24:22877-22888. [PMID: 36124732 DOI: 10.1039/d2cp03112d] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We perform all-atom molecular dynamics simulations of poly(2-hydroxyethyl methacrylate) (PHEMA) brushes in aqueous solutions of isobutane, propionamide, and sodium propionate. These solutes are side chain analogues to leucine, glutamine, and glutamic acid, respectively. We compute the Gibbs energy profile of the solute's adsorption to the polymer brush and decompose it into the contributions from the steric repulsion, van der Waals interaction, and Coulomb interaction to reveal the energetic origin of repulsion or attraction of the solute by the polymer brush. The Henry adsorption constant is the amount of adsorption normalized by the concentration in aqueous solution. We examine the dependence of this quantity on the grafting density and chain length. Our results suggest that the concurrent primary and ternary adsorption mechanism may be more important than previously expected when the solute is hydrophobic.
Collapse
Affiliation(s)
- Takuma Yagasaki
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka 560-8531, Japan.
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka 560-8531, Japan.
| |
Collapse
|
7
|
Wang D, Li D, Kelland MA, Cai H, Wang J, Xu Y, Lu P, Dong J. Unraveling Amphiphilic Poly( N-vinylcaprolactam)/Water Interface by Nuclear Magnetic Resonance Relaxometry: Control of Clathrate Hydrate Formation Kinetics. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2022; 38:4774-4784. [PMID: 35380846 DOI: 10.1021/acs.langmuir.2c00472] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Water-soluble amphiphilic polymers are vital chemicals in the oil and gas industry to retard crystal growth of hydrocarbon hydrate via surface adsorption and suppress nucleation of a pristine hydrate nucleus, thereby preventing formation of hydrate blockages in flow lines during oil and natural gas production. Apart from a few theoretical modeling studies, an experimental method to study the polymer/water interface in the crystal growth is critically needed. Here, water motions in the hydration shells of an exemplary kinetic inhibitor, poly(N-vinylcaprolactam), during hydrate formation from the tetrahydrofuran/water system are revealed via nuclear magnetic resonance relaxometry. Unequivocal experiments show that the pivotal interfacial water in the tightly bound state gradually freezes at rates depending on the polymer molecular weight (MW). This is supported by nonfreezable water analysis, which is correlated to the inhibition time. The polymers tune the kinetics of the hydration process via interaction with and perturbation of the water molecules. The free water component in the polymer solution crystallizes at a very slow rate when in partially restricted mobility, whereas the bound water component increases in the reaction, with the polymer/water interface serving as the reaction sites. The appropriate MW (including average MW and polydispersity values) of the inhibitive polymers can give rise to maximal retardation of the hydrate crystal growth. This work will help control other multiphase crystallization kinetic processes through the design of inhibitors or promoters functioning in the interface.
Collapse
Affiliation(s)
- Dong Wang
- College of Chemistry and Chemical Engineering, Shaoxing University, Shaoxing, Zhejiang Province 312000, China
| | - Dongfang Li
- College of Chemistry and Chemical Engineering, Shaoxing University, Shaoxing, Zhejiang Province 312000, China
| | - Malcolm A Kelland
- Department of Chemistry, Bioscience and Environmental Engineering, Faculty of Science and Technology, University of Stavanger, Stavanger N-4036, Norway
| | - Haokun Cai
- Ningbo Academy of Product and Food Quality Inspection (Ningbo Fiber Inspection Institute), Ningbo, Zhejiang Province 315048, China
| | - Jie Wang
- College of Chemistry and Chemical Engineering, Shaoxing University, Shaoxing, Zhejiang Province 312000, China
| | - Ying Xu
- College of Chemistry and Chemical Engineering, Shaoxing University, Shaoxing, Zhejiang Province 312000, China
| | - Ping Lu
- College of Chemistry and Chemical Engineering, Shaoxing University, Shaoxing, Zhejiang Province 312000, China
| | - Jian Dong
- College of Chemistry and Chemical Engineering, Shaoxing University, Shaoxing, Zhejiang Province 312000, China
| |
Collapse
|
8
|
Xu K, Lin Y, Shi Q, Li T, Zhang Z, Wu J. Role of mechanical deformation in the thermal transport of sI-type methane hydrate. Phys Chem Chem Phys 2022; 24:5479-5488. [PMID: 35171155 DOI: 10.1039/d1cp04189d] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Natural gas hydrates (NGHs) are rising as an unconventional energy resource. The fundamental thermal characteristics of NGHs are of importance for natural gas exploitation from permafrost and oceanic sediments that are geomechanically deformed. Here, utilizing classic molecular dynamics simulations with all-atom (AA) and coarse-grained (CG) models of the methane guest molecule, the effects of mechanical strain on the thermal conductivity of sI-type methane hydrate are for the first time examined. Upon triaxial tension and compression, methane hydrate exhibits strong asymmetry in the stress responses. As the triaxial loads go from compression to tension, a reduction trend in the thermal conductivity is revealed for methane hydrate with both AA and CG models of methane, within a maximum reduction of over 44%. This reduction is because triaxial strain from compression to tension softens the phonon modes. Interestingly, there is a sudden rise in thermal conductivity at critical triaxial strain of 0.06, originating from that, at which, the phonon modes are hardened and the peaks of radial distribution functions are shifted back. This study provides important information on the thermal conductivity of methane hydrate, which is helpful for the practical production of natural gas from geo-deformed NGH-bearing sediments via a heating technique as well as evaluating their stability.
Collapse
Affiliation(s)
- Ke Xu
- Department of Physics, Research Institute for Biomimetics and Soft Matter, Jiujiang Research Institute and Fujian Provincial Key Laboratory for Soft Functional Materials Research, Xiamen University, Xiamen 361005, P. R. China.
| | - Yanwen Lin
- Department of Physics, Research Institute for Biomimetics and Soft Matter, Jiujiang Research Institute and Fujian Provincial Key Laboratory for Soft Functional Materials Research, Xiamen University, Xiamen 361005, P. R. China.
| | - Qiao Shi
- Department of Physics, Research Institute for Biomimetics and Soft Matter, Jiujiang Research Institute and Fujian Provincial Key Laboratory for Soft Functional Materials Research, Xiamen University, Xiamen 361005, P. R. China.
| | - Tong Li
- Key Laboratory of Marine Materials and Related Technologies, Zhejiang Key Laboratory of Marine Materials and Protective Technologies, Ningbo Institute of Materials Technology and Engineering, Chinese Academy of Sciences, Ningbo 315201, P. R. China
| | - Zhisen Zhang
- Department of Physics, Research Institute for Biomimetics and Soft Matter, Jiujiang Research Institute and Fujian Provincial Key Laboratory for Soft Functional Materials Research, Xiamen University, Xiamen 361005, P. R. China.
| | - Jianyang Wu
- Department of Physics, Research Institute for Biomimetics and Soft Matter, Jiujiang Research Institute and Fujian Provincial Key Laboratory for Soft Functional Materials Research, Xiamen University, Xiamen 361005, P. R. China. .,NTNU Nanomechanical Lab, Norwegian University of Science and Technology (NTNU), Trondheim 7491, Norway
| |
Collapse
|
9
|
Liu J, Yan Y, Chen G, Zhang J, Liu S. Kinetic hydrate inhibitor performance and adsorption characteristics of poly(N-alkyl-N-vinyl acetamide)s: A first-principles study. Colloids Surf A Physicochem Eng Asp 2022. [DOI: 10.1016/j.colsurfa.2021.128097] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
|
10
|
Laroui A, Kelland MA, Wang D, Xu S, Xu Y, Lu P, Dong J. Kinetic Inhibition of Clathrate Hydrate by Copolymers Based on N-Vinylcaprolactam and N-Acryloylpyrrolidine: Optimization Effect of Interfacial Nonfreezable Water of Polymers. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2022; 38:1522-1532. [PMID: 35067060 DOI: 10.1021/acs.langmuir.1c02903] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Amphiphilic polymers have now been designed to achieve an icephobic performance and have been used for ice adhesion prevention. They may function by forming a strongly bonded but nonfreezable water shell which serves as a self-lubricating interfacial layer that weakens the adhesion strength between ice and the surface. Here, an analogous concept is built to prevent the formation of clathrate hydrate compounds during oil and natural gas production, in which amphiphilic water-soluble polymers act as efficient kinetic hydrate inhibitors (KHIs). A novel group of copolymers with N-vinylcaprolactam and N-acryloylpyrrolidine structural units are investigated in this study. The relationships among the amphiphilicity, lower critical solution temperature, nonfreezable bound water, and kinetic hydrate inhibition time are analyzed in terms of the copolymer compositions. Low-field NMR relaxometry revealed the crucial interfacial water in tightly bound dynamic states which led to crystal growth rates changing with the copolymer compositions, in accord with the rotational rheometric analysis results. The nonfreezable bound water layer confirmed by a calorimetry analysis also changes with the polymer amphiphilicity. Therefore, in the interface between the KHI polymers and hydrate, water surrounding the polymers plays a critical role by helping to delay the nucleation and growth of embryonic ice/hydrates. Appropriate amphiphilicity of the copolymers can achieve the optimal interfacial properties for slowing down hydrate crystal growth.
Collapse
Affiliation(s)
- Abdelatif Laroui
- College of Chemistry and Chemical Engineering, Shaoxing University, Shaoxing, Zhejiang Province 312000, China
| | - Malcolm A Kelland
- Department of Chemistry, Bioscience and Environmental Engineering, Faculty of Science and Technology, University of Stavanger, N-4036 Stavanger, Norway
| | - Dong Wang
- College of Chemistry and Chemical Engineering, Shaoxing University, Shaoxing, Zhejiang Province 312000, China
| | - Siyuan Xu
- College of Chemistry and Chemical Engineering, Shaoxing University, Shaoxing, Zhejiang Province 312000, China
| | - Ying Xu
- College of Chemistry and Chemical Engineering, Shaoxing University, Shaoxing, Zhejiang Province 312000, China
| | - Ping Lu
- College of Chemistry and Chemical Engineering, Shaoxing University, Shaoxing, Zhejiang Province 312000, China
| | - Jian Dong
- College of Chemistry and Chemical Engineering, Shaoxing University, Shaoxing, Zhejiang Province 312000, China
| |
Collapse
|
11
|
Liu J, Yan Y, Chen G, Zhang J, Liu S. Adsorption behavior of kinetic inhibitors on hydrate surfaces and its relation to the inhibition performance. Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2021.139108] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
|
12
|
Molecular dynamics simulation of the formation of methane hydrates in the presence of KHIs. Chem Eng Sci 2021. [DOI: 10.1016/j.ces.2021.116508] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
13
|
Bak IG, Heo CH, Kelland MA, Lee E, Kang BG, Lee JS. Clathrate Hydrate Inhibition by Polyisocyanate with Diethylammonium Group. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2021; 37:4147-4153. [PMID: 33794088 DOI: 10.1021/acs.langmuir.0c03663] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Polymers containing amide groups have been used as kinetic hydrate inhibitors (KHIs). The amide group has good performance for hydrate nucleus adsorption, resulting in inhibition of hydrate growth. Polyisocyanates composed of an amide backbone can be KHI candidates; however, the use of polyisocyanates as KHIs has not yet been reported. Herein, we prepared water-soluble poly[3-[[2-(diethylamino)ethyl]thio]-1-propyl isocyanate-ran-hexyl isocyanate] (P(DETPIC-ran-HIC)) to investigate the ability of polyisocyanates to inhibit hydrate formation. In the tetrahydrofuran (THF) clathrate hydrate crystal growth inhibition tests, P(DETPIC-ran-HIC) showed better performance than the polyamide, poly(N-vinylpyrrolidone) (PVP).
Collapse
Affiliation(s)
- In Gyu Bak
- School of Materials Science and Engineering, Gwangju Institute of Science and Technology (GIST), 123 Cheomdangwagi-ro, Buk-gu, Gwangju 61005, Korea
| | - Chi-Ho Heo
- School of Materials Science and Engineering, Gwangju Institute of Science and Technology (GIST), 123 Cheomdangwagi-ro, Buk-gu, Gwangju 61005, Korea
| | - Malcolm A Kelland
- Department of Chemistry, Bioscience and Environmental Engineering, University of Stavanger, N-4036 Stavanger, Norway
| | - Eunji Lee
- School of Materials Science and Engineering, Gwangju Institute of Science and Technology (GIST), 123 Cheomdangwagi-ro, Buk-gu, Gwangju 61005, Korea
| | - Beom-Goo Kang
- Department of Chemical Engineering, Soongsil University, Seoul 06978, Korea
| | - Jae-Suk Lee
- School of Materials Science and Engineering, Gwangju Institute of Science and Technology (GIST), 123 Cheomdangwagi-ro, Buk-gu, Gwangju 61005, Korea
| |
Collapse
|
14
|
Matsumoto M, Yagasaki T, Tanaka H. Formation of hot ice caused by carbon nanobrushes. II. Dependency on the radius of nanotubes. J Chem Phys 2021; 154:094502. [PMID: 33685157 DOI: 10.1063/5.0044300] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Stable crystalline structures of confined water can be different from bulk ice. In Paper I [T. Yagasaki et al., J. Chem. Phys. 151, 064702 (2019)] of this study, it was shown, using molecular dynamics (MD) simulations, that a zeolite-like ice structure forms in nanobrushes consisting of (6,6) carbon nanotubes (CNTs) when the CNTs are located in a triangle arrangement. The melting temperature of the zeolite-like ice structure is much higher than the melting temperature of ice Ih when the distance between the surfaces of CNTs is ∼0.94 nm, which is the best spacing for the bilayer structure of water. In this paper, we perform MD simulations of nanobrushes of CNTs that are different from (6,6) CNTs in radius. Several new porous ice structures form spontaneously in the MD simulations. A stable porous ice forms when the radius of its cavities matches the radius of the CNTs well. All cylindrical porous ice structures found in this study can be decomposed into a small number of structural blocks. We provide a new protocol to classify cylindrical porous ice crystals on the basis of this decomposition.
Collapse
Affiliation(s)
- Masakazu Matsumoto
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Takuma Yagasaki
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Osaka 560-8531, Japan
| | - Hideki Tanaka
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| |
Collapse
|
15
|
Mohr S, Hoevelmann F, Wylde J, Schelero N, Sarria J, Purkayastha N, Ward Z, Navarro Acero P, Michalis VK. Ranking the Efficiency of Gas Hydrate Anti-agglomerants through Molecular Dynamic Simulations. J Phys Chem B 2021; 125:1487-1502. [PMID: 33529037 DOI: 10.1021/acs.jpcb.0c08969] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Using both computational and experimental methods, the capacity of four different surfactant molecules to inhibit the agglomeration of sII hydrate particles was assessed. The computational simulations were carried out using both steered and non-steered molecular dynamics (MD), simulating the coalescence process of a hydrate slab and a water droplet, both covered with surfactant molecules. The surfactants were ranked according to free energy calculations (steered MD) and the number of agglomeration events (non-steered MD). The experimental work was based on rocking cell measurements, determining the minimum effective dose necessary to inhibit agglomeration. Overall, good agreement was obtained between the performance predicted by the simulations and the experimental measurements. Moreover, the simulations allowed us to gain additional insights that are not directly accessible via experiments, such as an analysis of the mass density profiles, the diffusion coefficients, or the orientations of the long tails.
Collapse
Affiliation(s)
- Stephan Mohr
- Nextmol (Bytelab Solutions SL), Barcelona 08018, Spain.,Barcelona Supercomputing Center (BSC), Barcelona 08034, Spain
| | | | - Jonathan Wylde
- Clariant Oil Services, Clariant Corporation, Houston, Texas 77258, United States.,Heriot-Watt University, Edinburgh EH14 4AS, Scotland, U.K
| | | | - Juan Sarria
- Clariant Produkte (Deutschland) GmbH, Frankfurt 65933, Germany
| | | | - Zachary Ward
- Clariant Oil Services, Clariant Corporation, Houston, Texas 77258, United States
| | | | | |
Collapse
|
16
|
Mochizuki K. On-Off of Co-non-solvency for Poly( N-vinylcaprolactam) in Alcohol-Water Mixtures: A Molecular Dynamics Study. J Phys Chem B 2020; 124:9951-9957. [PMID: 33086006 DOI: 10.1021/acs.jpcb.0c07188] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Poly(N-vinylcaprolactam) (PVCL) exhibits co-non-solvency in aqueous solutions of 2-propanol but not in methanol. What distinguishes the impact of these two cosolvents on the polymer conformational stability? We report a molecular dynamics simulation study on PVCL 50-mer and monomers dissolved in methanol-water and 2-propanol-water mixtures. We show that the alcohol-concentration dependence of the effective attraction between a pair of PVCL monomers closely resembles the conformational changes in a single PVCL 50-mer as well as the experimentally observed behavior for PVCL chains. We also found that, at the co-non-solvency maximum, the monomer-monomer attraction works over a long-range beyond the solvent-separated distance. Then, we correlate the long-range attraction to the appearance of a dense alcohol concentration accumulated between the monomers. Furthermore, we distinctly demonstrate that the co-non-solvency of PVCL monomers can be switched on/off by artificially tuning the alcohol size while keeping the energetic parameters. Thus, we conclude that the magnitude of the excluded volume effect in alcohol accompanying the gain in translational entropy of the solvent is crucial to the occurrence of PVCL polymer co-non-solvency.
Collapse
Affiliation(s)
- Kenji Mochizuki
- Department of Chemistry, Zhejiang University, 148 Tianmushan Road, Hangzhou, Zhejiang 310028, P. R. China
| |
Collapse
|
17
|
Yagasaki T, Matsumoto M, Tanaka H. Molecular dynamics study of grain boundaries and triple junctions in ice. J Chem Phys 2020; 153:124502. [DOI: 10.1063/5.0021635] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Takuma Yagasaki
- Research Institute for Interdisciplinary Science and Department of Chemistry, Faculty of Science, Okayama University, Okayama 700-8530, Japan
| | - Masakazu Matsumoto
- Research Institute for Interdisciplinary Science and Department of Chemistry, Faculty of Science, Okayama University, Okayama 700-8530, Japan
| | - Hideki Tanaka
- Research Institute for Interdisciplinary Science and Department of Chemistry, Faculty of Science, Okayama University, Okayama 700-8530, Japan
| |
Collapse
|
18
|
Mohsenzade H, Foroutan S, Dashti A, Ramezanian N, Roosta H. Vinyl lactam-based copolymers and terpolymers as high cloud point kinetic hydrate inhibitors in methane-THF-water system. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.113068] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
19
|
Mochizuki K. Reduction of water-mediated repulsion drives poly(N-vinylcaprolactam) collapse upon heating. Phys Chem Chem Phys 2020; 22:1053-1060. [PMID: 31867584 DOI: 10.1039/c9cp05491j] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Thermo-sensitive aqueous polymers undergo a coil-to-globule transition on heating, with drastic chemical and structural changes. We performed molecular dynamics simulations for PVCL in water to study the driving forces for the polymer's collapse.
Collapse
Affiliation(s)
- Kenji Mochizuki
- Department of Chemistry and Materials
- Faculty of Textile Science and Technology
- Shinshu University
- Nagano 386-8567
- Japan
| |
Collapse
|
20
|
Fang B, Ning F, Hu S, Guo D, Ou W, Wang C, Wen J, Sun J, Liu Z, Koh CA. The effect of surfactants on hydrate particle agglomeration in liquid hydrocarbon continuous systems: a molecular dynamics simulation study. RSC Adv 2020; 10:31027-31038. [PMID: 35520650 PMCID: PMC9056346 DOI: 10.1039/d0ra04088f] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2020] [Accepted: 07/14/2020] [Indexed: 12/02/2022] Open
Abstract
Anti-agglomerants (AAs), both natural and commercial, are currently being considered for gas hydrate risk management of petroleum pipelines in offshore operations. However, the molecular mechanisms of the interaction between the AAs and gas hydrate surfaces and the prevention of hydrate agglomeration remain critical and complex questions that need to be addressed to advance this technology. Here, we use molecular dynamics (MD) simulations to investigate the effect of model surfactant molecules (polynuclear aromatic carboxylic acids) on the agglomeration behaviour of gas hydrate particles and disruption of the capillary liquid bridge between hydrate particles. The results show that the anti-agglomeration pathway can be divided into two processes: the spontaneous adsorption effect of surfactant molecules onto the hydrate surface and the weakening effect of the intensity of the liquid bridge between attracted hydrate particles. The MD simulation results also indicate that the anti-agglomeration effectiveness of surfactants is determined by the intrinsic nature of their molecular functional groups. Additionally, we find that surfactant molecules can affect hydrate growth, which decreases hydrate particle size and correspondingly lower the risk of hydrate agglomeration. This study provides molecular-level insights into the anti-agglomeration mechanism of surfactant molecules, which can aid in the ultimate application of natural or commercial AAs with optimal anti-agglomeration properties. Schematic of anti-agglomeration effect of surfactants promoting gas hydrate particle dispersion.![]()
Collapse
|
21
|
Inhibition Effect of Kinetic Hydrate Inhibitors on the Growth of Methane Hydrate in Gas–Liquid Phase Separation State. ENERGIES 2019. [DOI: 10.3390/en12234482] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
The effect of kinetic hydrate inhibitors (KHIs) on the growth of methane hydrate in the gas–liquid phase separation state is studied at the molecular level. The simulation results show that the kinetic inhibitors, named PVP and PVP-A, show good inhibitory effects on the growth of methane hydrate under the gas–liquid phase separation state, and the initial position of the kinetic hydrate inhibitors has a major effect on the growth of methane hydrates. In addition, inhibitors at different locations exhibit different inhibition performances. When the inhibitor molecules are located at the gas–liquid phase interface, increasing the contact area between the groups of the inhibitor molecules and methane is beneficial to enhance the inhibitory performance of the inhibitors. When inhibitor molecules are located at the solid–liquid phase interface, the inhibitor molecules adsorbed on the surface of the hydrate nucleus and decreased the direct contact of hydrate nucleus with the surrounding water and methane molecules, which would delay the growth of hydrate nucleus. Moreover, the increase of hydrate surface curvature and the Gibbs–Thomson effect caused by inhibitors can also reduce the growth rate of methane hydrate.
Collapse
|
22
|
Wang Y, Fan S, Lang X. Reviews of gas hydrate inhibitors in gas-dominant pipelines and application of kinetic hydrate inhibitors in China. Chin J Chem Eng 2019. [DOI: 10.1016/j.cjche.2019.02.023] [Citation(s) in RCA: 49] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
|
23
|
Yagasaki T, Matsumoto M, Tanaka H. Liquid-liquid separation of aqueous solutions: A molecular dynamics study. J Chem Phys 2019; 150:214506. [DOI: 10.1063/1.5096429] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Takuma Yagasaki
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Masakazu Matsumoto
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Hideki Tanaka
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| |
Collapse
|
24
|
Naullage P, Bertolazzo AA, Molinero V. How Do Surfactants Control the Agglomeration of Clathrate Hydrates? ACS CENTRAL SCIENCE 2019; 5:428-439. [PMID: 30937370 PMCID: PMC6439454 DOI: 10.1021/acscentsci.8b00755] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/16/2018] [Indexed: 05/14/2023]
Abstract
Clathrate hydrates can spontaneously form under typical conditions found in oil and gas pipelines. The agglomeration of clathrates into large solid masses plugs the pipelines, posing adverse safety, economic, and environmental threats. Surfactants are customarily used to prevent the aggregation of clathrate particles and their coalescence with water droplets. It is generally assumed that a large contact angle between the surfactant-covered clathrate and water is a key predictor of the antiagglomerant performance of the surfactant. Here we use molecular dynamic simulations to investigate the structure and dynamics of surfactant films at the clathrate-oil interface, and their impact on the contact angle and coalescence between water droplets and hydrate particles. In agreement with the experiments, the simulations predict that surfactant-covered clathrate-oil interfaces are oil wet but super-hydrophobic to water. Although the water contact angle determines the driving force for coalescence, we find that a large contact angle is not sufficient to predict good antiagglomerant performance of a surfactant. We conclude that the length of the surfactant molecules, the density of the interfacial film, and the strength of binding of its molecules to the clathrate surface are the main factors in preventing the coalescence and agglomeration of clathrate particles with water droplets in oil. Our analysis provides a molecular foundation to guide the molecular design of effective clathrate antiagglomerants.
Collapse
Affiliation(s)
- Pavithra
M. Naullage
- Department
of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
| | - Andressa A. Bertolazzo
- Department
of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
- Departamento
de Ciências Exatas e Educação, Universidade Federal de Santa Catarina, Blumenau, Santa Catarina, Brazil
| | - Valeria Molinero
- Department
of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
- E-mail:
| |
Collapse
|
25
|
Veesam SK, Ravipati S, Punnathanam SN. Recent advances in thermodynamics and nucleation of gas hydrates using molecular modeling. Curr Opin Chem Eng 2019. [DOI: 10.1016/j.coche.2019.01.003] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
26
|
Ma S, Sun L, Kelland MA, Wang Q, Li D, Zhang Y, Dong J. Hydrophobic hydration affects growth of clathrate hydrate: insight from an NMR relaxometric and calorimetric study. Chem Commun (Camb) 2019; 55:2936-2939. [DOI: 10.1039/c8cc09587f] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Water tightly bound to the kinetic inhibitors of tetrahydrofuran hydrate is related to the hydrophobic hydration effect of the inhibitors.
Collapse
Affiliation(s)
- Shang Ma
- College of Chemistry and Chemical Engineering
- Shaoxing University
- Shaoxing
- China
| | - Li Sun
- College of Chemistry and Chemical Engineering
- Shaoxing University
- Shaoxing
- China
- School of Materials Science and Chemical Engineering
| | - Malcolm A. Kelland
- Department of Chemistry, Bioscience and Environmental Engineering
- Faculty of Science and Technology
- University of Stavanger
- N-4036 Stavanger
- Norway
| | - Qingyu Wang
- College of Chemistry and Chemical Engineering
- Shaoxing University
- Shaoxing
- China
| | - Dongfang Li
- College of Chemistry and Chemical Engineering
- Shaoxing University
- Shaoxing
- China
| | - Yixian Zhang
- College of Chemistry and Chemical Engineering
- Shaoxing University
- Shaoxing
- China
| | - Jian Dong
- College of Chemistry and Chemical Engineering
- Shaoxing University
- Shaoxing
- China
| |
Collapse
|
27
|
Kondo H, Mochizuki K, Bayer-Giraldi M. Multiple binding modes of a moderate ice-binding protein from a polar microalga. Phys Chem Chem Phys 2018; 20:25295-25303. [PMID: 30255887 DOI: 10.1039/c8cp04727h] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Ice-binding proteins (IBPs) produced by cold-tolerant organisms interact with ice and strongly control crystal growth. The molecular basis for the different magnitudes of activity displayed by various IBPs (moderate and hyperactive) has not yet been clarified. Previous studies questioned whether the moderate activity of some IBPs relies on their weaker binding modus to the ice surface, compared to hyperactive IBPs, rather than relying on binding only to selected faces of the ice crystal. We present the structure of one moderate IBP from the sea-ice diatom Fragilariopsis cylindrus (fcIBP) as determined by X-ray crystallography and investigate the protein's binding modes to the growing ice-water interface using molecular dynamics simulations. The structure of fcIBP is the IBP-1 fold, defined by a discontinuous β-solenoid delimitated by three faces (A, B and C-faces) and braced by an α-helix. The fcIBP structure shows capping loops on both N- and C-terminal parts of the solenoid. We show that the protein adsorbs on both the prism and the basal faces of ice crystals, confirming experimental results. The fcIBP binds irreversibly to the prism face using the loop between the B and the C-faces, involving also the B-face in water immobilization despite its irregular structure. The α-helix attaches the protein to the basal face with a partly reversible modus. Our results suggest that fcIBP has a looser attachment to ice and that this weaker binding modus is the basis to explain the moderate activity of fcIBP.
Collapse
Affiliation(s)
- Hidemasa Kondo
- Bioproduction Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), Sapporo 062-8517, Japan
| | | | | |
Collapse
|
28
|
Yagasaki T, Matsumoto M, Tanaka H. Phase Diagrams of TIP4P/2005, SPC/E, and TIP5P Water at High Pressure. J Phys Chem B 2018; 122:7718-7725. [DOI: 10.1021/acs.jpcb.8b04441] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Takuma Yagasaki
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Masakazu Matsumoto
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Hideki Tanaka
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| |
Collapse
|
29
|
Bertolazzo AA, Naullage PM, Peters B, Molinero V. The Clathrate-Water Interface Is Oleophilic. J Phys Chem Lett 2018; 9:3224-3231. [PMID: 29812945 DOI: 10.1021/acs.jpclett.8b01210] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/23/2023]
Abstract
The slow nucleation of clathrate hydrates is a central challenge for their use in the storage and transportation of natural gas. Molecules that strongly adsorb to the clathrate-water interface decrease the crystal-water surface tension, lowering the barrier for clathrate nucleation. Surfactants are widely used to promote the nucleation and growth of clathrate hydrates. It has been proposed that these amphiphilic molecules bind to the clathrate surface via hydrogen bonding. However, recent studies reveal that PVCap, an amphiphilic polymer, binds to clathrates through hydrophobic moieties. Here we use molecular dynamic simulations and theory to investigate the mode and strength of binding of surfactants to the clathrate-water interface and their effect on the nucleation rate. We find that the surfactants bind to the clathrate-water interface exclusively through their hydrophobic tails. The binding is strong, driven by the entropy of dehydration of the alkyl chain, as it penetrates empty cavities at the hydrate surface. The hydrophobic attraction of alkyl groups to the clathrate surface also results in strong adsorption of alkanes. We identify two regimes for the binding of surfactants as a function of their density at the hydrate surface, which we interpret to correspond to the two steps of the Langmuir adsorption isotherm observed in experiments. Our results indicate that hydrophobic attraction to the clathrate-water interface is key for the design of soluble additives that promote the nucleation of hydrates. We use the calculated adsorption coefficients to estimate the concentration of sodium dodecyl sulfate (SDS) required to reach nucleation rates for methane hydrate consistent with those measured in experiments. To our knowledge, this study is the first to quantify the effect of surfactant concentration in the nucleation rate of clathrate hydrates.
Collapse
Affiliation(s)
- Andressa A Bertolazzo
- Department of Chemistry , The University of Utah , Salt Lake City , Utah 84112-0580 , United States
- Departamento de Ciências Exatas e Educação , Universidade Federal de Santa Catarina , Blumenau , State of Santa Catarina 88040-900 , Brazil
| | - Pavithra M Naullage
- Department of Chemistry , The University of Utah , Salt Lake City , Utah 84112-0580 , United States
| | - Baron Peters
- Department of Chemical Engineering , University of California , Santa Barbara , California 93106 , United States
| | - Valeria Molinero
- Department of Chemistry , The University of Utah , Salt Lake City , Utah 84112-0580 , United States
| |
Collapse
|