1
|
Yang Y, Wan J, Li J, Zhao G, Shang X. Molecular modeling of interfacial properties of the hydrogen + water + decane mixture in three-phase equilibrium. Phys Chem Chem Phys 2023; 25:29641-29655. [PMID: 37881893 DOI: 10.1039/d3cp04406h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2023]
Abstract
The understanding of interfacial phenomena between H2 and geofluids is of great importance for underground H2 storage, but requires further study. We report the first investigation on the three-phase fluid mixture containing H2, H2O, and n-C10H22. Molecular dynamics simulation and PC-SAFT density gradient theory are employed to estimate the interfacial properties under various conditions (temperature ranges from 298 to 373 K and pressure is up to around 100 MPa). Our results demonstrate that interfacial tensions (IFTs) of the H2-H2O interface in the H2 + H2O + C10H22 three-phase mixture are smaller than IFTs in the H2 + H2O two-phase mixture. This decrement of IFT can be attributed to C10H22 adsorption in the interface. Importantly, H2 accumulates in the H2O-C10H22 interface in the three-phase systems, which leads to weaker increments of IFT with increasing pressure compared to IFTs in the water + C10H22 two-phase mixture. In addition, the IFTs of the H2-C10H22 interface are hardly influenced by H2O due to the limited amount of H2O dissolved in nonaqueous phases. Nevertheless, positive surface excesses of H2O are seen in the H2-C10H22 interfacial region. Furthermore, the values of the spreading coefficient are mostly negative revealing the presence of the three-phase contact for the H2 + H2O + C10H22 mixture under studied conditions.
Collapse
Affiliation(s)
- Yafan Yang
- State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology, Xuzhou, China.
| | - Jingyu Wan
- State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology, Xuzhou, China.
| | - Jingfa Li
- School of Mechanical Engineering and Hydrogen Energy Research Center, Beijing Institute of Petrochemical Technology, Beijing, China
| | - Guangsi Zhao
- State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology, Xuzhou, China.
| | - Xiangyu Shang
- State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology, Xuzhou, China.
| |
Collapse
|
2
|
Shao L, Ma J, Prelesnik JL, Zhou Y, Nguyen M, Zhao M, Jenekhe SA, Kalinin SV, Ferguson AL, Pfaendtner J, Mundy CJ, De Yoreo JJ, Baneyx F, Chen CL. Hierarchical Materials from High Information Content Macromolecular Building Blocks: Construction, Dynamic Interventions, and Prediction. Chem Rev 2022; 122:17397-17478. [PMID: 36260695 DOI: 10.1021/acs.chemrev.2c00220] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Hierarchical materials that exhibit order over multiple length scales are ubiquitous in nature. Because hierarchy gives rise to unique properties and functions, many have sought inspiration from nature when designing and fabricating hierarchical matter. More and more, however, nature's own high-information content building blocks, proteins, peptides, and peptidomimetics, are being coopted to build hierarchy because the information that determines structure, function, and interfacial interactions can be readily encoded in these versatile macromolecules. Here, we take stock of recent progress in the rational design and characterization of hierarchical materials produced from high-information content blocks with a focus on stimuli-responsive and "smart" architectures. We also review advances in the use of computational simulations and data-driven predictions to shed light on how the side chain chemistry and conformational flexibility of macromolecular blocks drive the emergence of order and the acquisition of hierarchy and also on how ionic, solvent, and surface effects influence the outcomes of assembly. Continued progress in the above areas will ultimately usher in an era where an understanding of designed interactions, surface effects, and solution conditions can be harnessed to achieve predictive materials synthesis across scale and drive emergent phenomena in the self-assembly and reconfiguration of high-information content building blocks.
Collapse
Affiliation(s)
- Li Shao
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Jinrong Ma
- Molecular Engineering and Sciences Institute, University of Washington, Seattle, Washington 98195, United States
| | - Jesse L Prelesnik
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Yicheng Zhou
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Mary Nguyen
- Department of Chemical Engineering, University of Washington, Seattle, Washington 98195, United States.,Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Mingfei Zhao
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Samson A Jenekhe
- Department of Chemical Engineering, University of Washington, Seattle, Washington 98195, United States.,Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Sergei V Kalinin
- Department of Materials Science and Engineering, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Andrew L Ferguson
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Jim Pfaendtner
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States.,Materials Science and Engineering, University of Washington, Seattle, Washington 98195, United States
| | - Christopher J Mundy
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States.,Department of Chemical Engineering, University of Washington, Seattle, Washington 98195, United States
| | - James J De Yoreo
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States.,Materials Science and Engineering, University of Washington, Seattle, Washington 98195, United States
| | - François Baneyx
- Molecular Engineering and Sciences Institute, University of Washington, Seattle, Washington 98195, United States.,Department of Chemical Engineering, University of Washington, Seattle, Washington 98195, United States
| | - Chun-Long Chen
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States.,Department of Chemical Engineering, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
3
|
Fayaz-Torshizi M, Graham EJ, Adjiman CS, Galindo A, Jackson G, Müller EA. SAFT- γ Force Field for the Simulation of Molecular Fluids 9: Coarse-Grained Models for Polyaromatic Hydrocarbons Describing Thermodynamic, Interfacial, Structural, and Transport Properties. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120827] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
4
|
Nkepsu Mbitou RL, Goujon F, Dequidt A, Latour B, Devémy J, Blaak R, Martzel N, Emeriau-Viard C, Tchoufag J, Garruchet S, Munch E, Hauret P, Malfreyt P. Consistent and Transferable Force Fields for Statistical Copolymer Systems at the Mesoscale. J Chem Theory Comput 2022; 18:6940-6951. [PMID: 36205431 DOI: 10.1021/acs.jctc.2c00945] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The statistical trajectory matching (STM) method was applied successfully to derive coarse grain (CG) models for bulk properties of homopolymers. The extension of the methodology for building CG models for statistical copolymer systems is much more challenging. We present here the strategy for developing CG models for styrene-butadiene-rubber, and we compare the quality of the resulting CG force fields on the structure and thermodynamics at different chemical compositions. The CG models are used through the use of a genuine mesoscopic method called the dissipative particle dynamics method and compared to high-resolution molecular dynamics simulations. We conclude that the STM method is able to produce coarse-grained potentials that are transferable in composition by using only a few reference systems. Additionally, this methodology can be applied on any copolymer system.
Collapse
Affiliation(s)
- R L Nkepsu Mbitou
- Université Clermont Auvergne, CNRS, Clermont Auvergne INP, Institut de Chimie de Clermont-Ferrand, F-63000Clermont-Ferrand, France.,Manufacture Française des Pneumatiques Michelin, 23, Place des Carmes, 63040Clermont-Ferrand, France
| | - F Goujon
- Université Clermont Auvergne, CNRS, Clermont Auvergne INP, Institut de Chimie de Clermont-Ferrand, F-63000Clermont-Ferrand, France
| | - A Dequidt
- Université Clermont Auvergne, CNRS, Clermont Auvergne INP, Institut de Chimie de Clermont-Ferrand, F-63000Clermont-Ferrand, France
| | - B Latour
- Manufacture Française des Pneumatiques Michelin, 23, Place des Carmes, 63040Clermont-Ferrand, France
| | - J Devémy
- Université Clermont Auvergne, CNRS, Clermont Auvergne INP, Institut de Chimie de Clermont-Ferrand, F-63000Clermont-Ferrand, France
| | - R Blaak
- Université Clermont Auvergne, CNRS, Clermont Auvergne INP, Institut de Chimie de Clermont-Ferrand, F-63000Clermont-Ferrand, France
| | - N Martzel
- Manufacture Française des Pneumatiques Michelin, 23, Place des Carmes, 63040Clermont-Ferrand, France
| | - C Emeriau-Viard
- Manufacture Française des Pneumatiques Michelin, 23, Place des Carmes, 63040Clermont-Ferrand, France
| | - J Tchoufag
- Manufacture Française des Pneumatiques Michelin, 23, Place des Carmes, 63040Clermont-Ferrand, France
| | - S Garruchet
- Manufacture Française des Pneumatiques Michelin, 23, Place des Carmes, 63040Clermont-Ferrand, France
| | - E Munch
- Manufacture Française des Pneumatiques Michelin, 23, Place des Carmes, 63040Clermont-Ferrand, France
| | - P Hauret
- Manufacture Française des Pneumatiques Michelin, 23, Place des Carmes, 63040Clermont-Ferrand, France
| | - P Malfreyt
- Université Clermont Auvergne, CNRS, Clermont Auvergne INP, Institut de Chimie de Clermont-Ferrand, F-63000Clermont-Ferrand, France
| |
Collapse
|
5
|
Wu Z, Müller-Plathe F. Slip-Spring Hybrid Particle-Field Molecular Dynamics for Coarse-Graining Branched Polymer Melts: Polystyrene Melts as an Example. J Chem Theory Comput 2022; 18:3814-3828. [PMID: 35617016 DOI: 10.1021/acs.jctc.2c00107] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The topology of chains significantly modifies the dynamical properties of polymer melts. Here, we extend a recently developed efficient simulation method, namely the slip-spring hybrid particle-field (SS-hPF) model, to study the structural and dynamical properties of branched polymer melts over large spatial-temporal scales. In the coarse-grained SS-hPF simulation of polymers, the bonded potentials are derived by iterative Boltzmann inversion from the underlying fine-grained model. The nonbonded potentials are computed from a density functional field instead of pairwise interactions used in standard molecular dynamics simulations, which increases the computational efficiency by a factor of 10-20. The entangled dynamics is lost due to the soft-core nature of density functional field interactions. It is recovered by a multichain slip-spring model that is rigorously parametrized from existing experimental or simulation data. To quantitatively predict the relaxation and diffusion of branched polymers, which are dominated by arm retraction rather than chain reptation, the slip-spring algorithm is augmented to improve the polymer dynamics near the branch point. Multiple dynamical observables, e.g., diffusion coefficients, arm relaxations, and tube survival probabilities, are characterized in an example coarse-grained model of symmetric and asymmetric star-shaped polystyrene melts. Consistent dynamical behaviors are identified and compared with theoretical predictions. With a single rescaling factor, the prediction of diffusion coefficients agrees well with the available experimental measurements. In this work, an efficient approach is provided to build chemistry-specific coarse-grained models for predicting the dynamics of branched polymers.
Collapse
Affiliation(s)
- Zhenghao Wu
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Alarich-Weiss-Strasse 8, 64287 Darmstadt, Germany
| | - Florian Müller-Plathe
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Alarich-Weiss-Strasse 8, 64287 Darmstadt, Germany
| |
Collapse
|
6
|
Boothroyd S, Madin OC, Mobley DL, Wang LP, Chodera JD, Shirts MR. Improving Force Field Accuracy by Training against Condensed-Phase Mixture Properties. J Chem Theory Comput 2022; 18:3577-3592. [PMID: 35533269 DOI: 10.1021/acs.jctc.1c01268] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Developing a sufficiently accurate classical force field representation of molecules is key to realizing the full potential of molecular simulations as a route to gaining a fundamental insight into a broad spectrum of chemical and biological phenomena. This is only possible, however, if the many complex interactions between molecules of different species in the system are accurately captured by the model. Historically, the intermolecular van der Waals (vdW) interactions have primarily been trained against densities and enthalpies of vaporization of pure (single-component) systems, with occasional usage of hydration free energies. In this study, we demonstrate how including physical property data of binary mixtures can better inform these parameters, encoding more information about the underlying physics of the system in complex chemical mixtures. To demonstrate this, we retrain a select number of Lennard-Jones parameters describing the vdW interactions of the OpenFF 1.0.0 (Parsley) fixed charge force field against training sets composed of densities and enthalpies of mixing for binary liquid mixtures as well as densities and enthalpies of vaporization of pure liquid systems and assess the performance of each of these combinations. We show that retraining against the mixture data improves the force field's ability to reproduce mixture properties, including solvation free energies, correcting some systematic errors that exist when training vdW interactions against properties of pure systems only.
Collapse
Affiliation(s)
- Simon Boothroyd
- Boothroyd Scientific Consulting Ltd., 71-75 Shelton Street, London WC2H 9JQ, Greater London, U.K
| | - Owen C Madin
- Department of Chemical and Biological Engineering, University of Colorado Boulder, Boulder, Colorado 80309, United States
| | - David L Mobley
- Department of Chemistry, University of California, Irvine, California 92617, United States.,Department of Pharmaceutical Sciences, University of California, Irvine, California 92617, United States
| | - Lee-Ping Wang
- Department of Chemistry, University of California, Davis, California 95616, United States
| | - John D Chodera
- Computational & Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, New York 10065, United States
| | - Michael R Shirts
- Department of Chemical and Biological Engineering, University of Colorado Boulder, Boulder, Colorado 80309, United States
| |
Collapse
|
7
|
Dong Y, Warsahartana HG, Hammad F, Masters A. SAFT
‐γ Mie model for ionic liquids. AIChE J 2022. [DOI: 10.1002/aic.17697] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Yichun Dong
- School of Chemistry and Chemical Engineering Hefei University of Technology Hefei China
| | | | - Faisal Hammad
- Department of Chemical Engineering The University of Manchester Manchester UK
| | - Andrew Masters
- Department of Chemical Engineering The University of Manchester Manchester UK
| |
Collapse
|
8
|
Fayaz-Torshizi M, Xu W, Vella JR, Marshall BD, Ravikovitch PI, Müller EA. Use of Boundary-Driven Nonequilibrium Molecular Dynamics for Determining Transport Diffusivities of Multicomponent Mixtures in Nanoporous Materials. J Phys Chem B 2022; 126:1085-1100. [PMID: 35104134 PMCID: PMC9007456 DOI: 10.1021/acs.jpcb.1c09159] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
The boundary-driven molecular modeling
strategy to evaluate mass
transport coefficients of fluids in nanoconfined media is revisited
and expanded to multicomponent mixtures. The method requires setting
up a simulation with bulk fluid reservoirs upstream and downstream
of a porous media. A fluid flow is induced by applying an external
force at the periodic boundary between the upstream and downstream
reservoirs. The relationship between the resulting flow and the density
gradient of the adsorbed fluid at the entrance/exit of the porous
media provides for a direct path for the calculation of the transport
diffusivities. It is shown how the transport diffusivities found this
way relate to the collective, Onsager, and self-diffusion coefficients,
typically used in other contexts to describe fluid transport in porous
media. Examples are provided by calculating the diffusion coefficients
of a Lennard-Jones (LJ) fluid and mixtures of differently sized LJ
particles in slit pores, a realistic model of methane in carbon-based
slit pores, and binary mixtures of methane with hypothetical counterparts
having different attractions to the solid. The method is seen to be
robust and particularly suited for the study of study of transport
of dense fluids and liquids in nanoconfined media.
Collapse
Affiliation(s)
- Maziar Fayaz-Torshizi
- Department of Chemical Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| | - Weilun Xu
- Department of Chemical Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| | - Joseph R Vella
- ExxonMobil Research and Engineering Company, Irving, Texas 75039-2298, United States
| | - Bennett D Marshall
- ExxonMobil Research and Engineering Company, Annandale, New Jersey 08801, United States
| | - Peter I Ravikovitch
- ExxonMobil Research and Engineering Company, Annandale, New Jersey 08801, United States
| | - Erich A Müller
- Department of Chemical Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| |
Collapse
|
9
|
Snider N. A Simple Equation for the Free Energy of Liquids. J Phys Chem B 2021; 125:12476-12485. [PMID: 34752086 DOI: 10.1021/acs.jpcb.1c05913] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A simple, three-parameter equation is proposed for the free energy of classical fluids composed of molecules which are not hydrogen bonded. Agreement with experiment is much the same as has been obtained from more complicated three-parameter equations. For spherical and nearly spherical molecules, semiquantitative relations between intermolecular pair potential parameters and the parameters of the proposed equation were found. For alkanes, relations between the parameters and carbon number were noted. These relationships contribute to the understanding of fluid state properties in terms of the properties of the molecules of which they are composed.
Collapse
Affiliation(s)
- Neil Snider
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208, United States
| |
Collapse
|
10
|
|
11
|
Fayaz‐Torshizi M, Müller EA. Coarse‐Grained Molecular Simulation of Polymers Supported by the Use of the SAFT‐γ$\gamma$ Mie Equation of State. MACROMOL THEOR SIMUL 2021. [DOI: 10.1002/mats.202100031] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
| | - Erich A. Müller
- Department of Chemical Engineering Imperial College London London SW7 2AZ UK
| |
Collapse
|
12
|
Walker CC, Genzer J, Santiso EE. Effect of Poly(vinyl butyral) Comonomer Sequence on Adhesion to Amorphous Silica: A Coarse-Grained Molecular Dynamics Study. ACS APPLIED MATERIALS & INTERFACES 2020; 12:47879-47890. [PMID: 32921047 DOI: 10.1021/acsami.0c10747] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Modulating a comonomer sequence, in addition to the overall chemical composition, is the key to unlocking the true potential of many existing commercial copolymers. We employ coarse-grained molecular dynamics (MD) simulations to study the behavior of random-blocky poly(vinyl butyral-co-vinyl alcohol) (PVB) melts in contact with an amorphous silica surface, representing the interface found in laminated safety glass. Our two-pronged coarse-graining approach utilizes both macroscopic thermophysical data and all-atom MD simulation data. Polymer-polymer nonbonded interactions are described by the fused-sphere SAFT-γ Mie equation of state, while bonded interactions are derived using Boltzmann inversion to match the bond and angle distributions from all-atom PVB chains. Spatially dependent polymer-surface interactions are mapped from a hydroxylated all-atom amorphous silica slab model and all-atom monomers to an external potential acting on the coarse-grained sites. We discovered an unexpected complex relationship between the blockiness parameter and the adhesion energy. The adhesion strength between PVB copolymers with intermediate VA content and silica was found to be maximal for random-blocky copolymers with a moderately high degree of blockiness rather than for diblock copolymers. We attribute this to two main factors: (1) changes in morphology, which dramatically alter the number of VA beads interacting with the surface and (2) a non-negligible contribution of vinyl butyral (VB) monomers to adhesion energy because of their preference to adsorb to zones with low hydroxyl density on the silica surface.
Collapse
Affiliation(s)
- Christopher C Walker
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - Jan Genzer
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - Erik E Santiso
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27695, United States
| |
Collapse
|
13
|
Zhu K, Müller EA. Generating a Machine-Learned Equation of State for Fluid Properties. J Phys Chem B 2020; 124:8628-8639. [DOI: 10.1021/acs.jpcb.0c05806] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Kezheng Zhu
- Department of Chemical Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| | - Erich A. Müller
- Department of Chemical Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| |
Collapse
|
14
|
Surrogate Models for Studying the Wettability of Nanoscale Natural Rough Surfaces Using Molecular Dynamics. ENERGIES 2020. [DOI: 10.3390/en13112770] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
A molecular modeling methodology is presented to analyze the wetting behavior of natural surfaces exhibiting roughness at the nanoscale. Using atomic force microscopy, the surface topology of a Ketton carbonate is measured with a nanometer resolution, and a mapped model is constructed with the aid of coarse-grained beads. A surrogate model is presented in which surfaces are represented by two-dimensional sinusoidal functions defined by both an amplitude and a wavelength. The wetting of the reconstructed surface by a fluid, obtained through equilibrium molecular dynamics simulations, is compared to that observed by the different realizations of the surrogate model. A least-squares fitting method is implemented to identify the apparent static contact angle, and the droplet curvature, relative to the effective plane of the solid surface. The apparent contact angle and curvature of the droplet are then used as wetting metrics. The nanoscale contact angle is seen to vary significantly with the surface roughness. In the particular case studied, a variation of over 65° is observed between the contact angle on a flat surface and on a highly spiked (Cassie–Baxter) limit. This work proposes a strategy for systematically studying the influence of nanoscale topography and, eventually, chemical heterogeneity on the wettability of surfaces.
Collapse
|
15
|
Hamani AWS, Bazile JP, Hoang H, Luc HT, Daridon JL, Galliero G. Thermophysical properties of simple molecular liquid mixtures: On the limitations of some force fields. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.112663] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
16
|
Walker CC, Genzer J, Santiso EE. Extending the fused-sphere SAFT-γ Mie force field parameterization approach to poly(vinyl butyral) copolymers. J Chem Phys 2020; 152:044903. [DOI: 10.1063/1.5126213] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
Affiliation(s)
- Christopher C. Walker
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27695, USA
| | - Jan Genzer
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27695, USA
| | - Erik E. Santiso
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27695, USA
| |
Collapse
|
17
|
Theodorakis PE, Smith ER, Müller EA. Spreading of aqueous droplets with common and superspreading surfactants. A molecular dynamics study. Colloids Surf A Physicochem Eng Asp 2019. [DOI: 10.1016/j.colsurfa.2019.123810] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
18
|
Molecular Dynamics Simulation of the Superspreading of Surfactant-Laden Droplets. A Review. FLUIDS 2019. [DOI: 10.3390/fluids4040176] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Superspreading is the rapid and complete spreading of surfactant-laden droplets on hydrophobic substrates. This phenomenon has been studied for many decades by experiment, theory, and simulation, but it has been only recently that molecular-level simulation has provided significant insights into the underlying mechanisms of superspreading thanks to the development of accurate force-fields and the increase of computational capabilities. Here, we review the main advances in this area that have surfaced from Molecular Dynamics simulation of all-atom and coarse-grained models highlighting and contrasting the main results and discussing various elements of the proposed mechanisms for superspreading. We anticipate that this review will stimulate further research on the interpretation of experimental results and the design of surfactants for applications requiring efficient spreading, such as coating technology.
Collapse
|
19
|
Jiménez-Serratos G, Cárdenas H, Müller EA. Extension of the effective solid-fluid Steele potential for Mie force fields. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1669836] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
| | - Harry Cárdenas
- Department of Chemical Engineering, Imperial College London, London, UK
| | - Erich A. Müller
- Department of Chemical Engineering, Imperial College London, London, UK
| |
Collapse
|
20
|
Kondratyuk N. Contributions of force field interaction forms to Green-Kubo viscosity integral in n-alkane case. J Chem Phys 2019; 151:074502. [PMID: 31438709 DOI: 10.1063/1.5103265] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Decades of molecular simulation history proved that the Green-Kubo method for shear viscosity converges without any problems in atomic and simple molecular liquids, unlike liquids with high values of viscosity. In the case of highly viscous liquids, the time decomposition method was developed in 2015 by Maginn and co-authors [Y. Zhang, A. Otani, and E. J. Maginn, J. Chem. Theory Comput. 11, 3537-3546 (2015)] which allows us to improve the convergence of the Green-Kubo integral. In this paper, the contributions of intramolecular and intermolecular force field parts to the viscosity integral are discovered to gain the understanding of the Green-Kubo method. The n-alkanes from n-ethane to n-pentane at 330 K in the optimized potentials for liquid simulations-all atom force field are used as reference models. The dependencies of these contributions and decay times of the corresponding correlation functions on the chain length are observed. The nonequilibrium simulations are carried out to verify the Green-Kubo results. The obtained values of viscosity are compared with experimental data.
Collapse
Affiliation(s)
- N Kondratyuk
- Moscow Institute of Physics and Technology (National Research University), 141700 Dolgoprudnyi, Russia and Joint Institute for High Temperatures of the Russian Academy of Sciences, 125412 Moscow, Russia
| |
Collapse
|
21
|
Pervaje AK, Walker CC, Santiso EE. Molecular simulation of polymers with a SAFT-γ Mie approach. MOLECULAR SIMULATION 2019. [DOI: 10.1080/08927022.2019.1645331] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- Amulya K. Pervaje
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC, USA
| | - Christopher C. Walker
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC, USA
| | - Erik E. Santiso
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC, USA
| |
Collapse
|
22
|
Aggregation Behavior of Model Asphaltenes Revealed from Large-Scale Coarse-Grained Molecular Simulations. J Phys Chem B 2019; 123:2380-2396. [PMID: 30735393 DOI: 10.1021/acs.jpcb.8b12295] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Fully atomistic simulations of models of asphaltenes in simple solvents have allowed the study of trends in aggregation phenomena to understand the underlying role played by molecular structure. The detail included at this scale of molecular modeling is, however, at odds with the required spatial and temporal resolution needed to fully understand asphaltene aggregation. The computational cost required to explore the relevant scales can be reduced by employing coarse-grained (CG) models, which consist of lumping a few atoms into a single segment that is characterized by effective interactions. In this work, CG force fields developed via the statistical associating fluid theory (SAFT-γ) [ Müller , E. A. ; Jackson , G. Annu. Rev. Chem. Biomol. Eng. 5 , 2014 , 405 - 427 ] equation of state (EoS) provide a reliable pathway to link the molecular description with macroscopic thermophysical data. A recent modification of the SAFT-VR EoS [ Müller , E. A. ; Mejía , A. Langmuir 33 , 2017 , 11518 - 11529 ], which allows for the parameterization of homonuclear rings, is selected as the starting point to develop CG models for polycyclic aromatic hydrocarbons. The new aromatic-core models, along with others published for simpler organic molecules, are adopted for the construction of asphaltene models by combining different chemical moieties in a group-contribution fashion. We apply the procedure to two previously reported asphaltene models and perform molecular dynamics simulations to validate the coarse-grained representation against benchmark systems of 27 asphaltenes in a pure solvent (toluene or heptane) described in a fully atomistic fashion. An excellent match between both levels of description is observed for the cluster size, radii of gyration, and relative-shape-anisotropy-factor distributions. We exploit the advantages of the CG representation by simulating systems containing up to 2000 asphaltene molecules in an explicit solvent investigating the effect of asphaltene concentration, solvent composition, and temperature on aggregation. By studying large systems facilitated by the use of CG models, we observe stable continuous distributions of molecular aggregates at conditions away from the two-phase precipitation point. As a further example application, a widely accepted interpretation of cluster-size distributions in asphaltenic systems is challenged by performing system-size tests, reversibility checks, and a time-dependence analysis. The proposed coarse-graining procedure is seen to be general and predictive and, hence, can be applied to other asphaltenic molecular structures.
Collapse
|
23
|
Shahruddin S, Jiménez-Serratos G, Britovsek GJP, Matar OK, Müller EA. Fluid-solid phase transition of n-alkane mixtures: Coarse-grained molecular dynamics simulations and diffusion-ordered spectroscopy nuclear magnetic resonance. Sci Rep 2019; 9:1002. [PMID: 30700804 PMCID: PMC6353884 DOI: 10.1038/s41598-018-37799-7] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2018] [Accepted: 12/12/2018] [Indexed: 11/10/2022] Open
Abstract
Wax appearance temperature (WAT), defined as the temperature at which the first solid paraffin crystal appears in a crude oil, is one of the key flow assurance indicators in the oil industry. Although there are several commonly-used experimental techniques to determine WAT, none provides unambiguous molecular-level information to characterize the phase transition between the homogeneous fluid and the underlying solid phase. Molecular Dynamics (MD) simulations employing the statistical associating fluid theory (SAFT) force field are used to interrogate the incipient solidification states of models for long-chain alkanes cooled from a melt to an arrested state. We monitor the phase change of pure long chain n-alkanes: tetracosane (C24H50) and triacontane (C30H62), and an 8-component surrogate n-alkane mixture (C12-C33) built upon the compositional information of a waxy crude. Comparison to Diffusion Ordered Spectroscopy Nuclear Magnetic Resonance (DOSY NMR) results allows the assessment of the limitations of the coarse-grained models proposed. We show that upon approach to freezing, the heavier components restrict their motion first while the lighter ones retain their mobility and help fluidize the mixture. We further demonstrate that upon sub-cooling of long n-alkane fluids and mixtures, a discontinuity arises in the slope of the self-diffusion coefficient with decreasing temperature, which can be employed as a marker for the appearance of an arrested state commensurate with conventional WAT measurements.
Collapse
Affiliation(s)
- S Shahruddin
- PETRONAS Research Sdn. Bhd, Lot 3288 & 3289 Off Jalan Ayer Itam, Kawasan Institusi Bangi, 43000, Kajang, Malaysia.,Department of Chemical Engineering, Imperial College London, South Kensington Campus, London, SW7 2AZ, UK
| | - G Jiménez-Serratos
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London, SW7 2AZ, UK
| | - G J P Britovsek
- Department of Chemistry, Imperial College London, South Kensington Campus, London, SW7 2AZ, UK
| | - O K Matar
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London, SW7 2AZ, UK
| | - E A Müller
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London, SW7 2AZ, UK.
| |
Collapse
|
24
|
An Y, Bejagam KK, Deshmukh SA. Development of Transferable Nonbonded Interactions between Coarse-Grained Hydrocarbon and Water Models. J Phys Chem B 2019; 123:909-921. [DOI: 10.1021/acs.jpcb.8b07990] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
|
25
|
Jiménez-Serratos G, Herdes C, Haslam AJ, Jackson G, Müller EA. Group Contribution Coarse-Grained Molecular Simulations of Polystyrene Melts and Polystyrene Solutions in Alkanes Using the SAFT-γ Force Field. Macromolecules 2017. [DOI: 10.1021/acs.macromol.6b02072] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
| | - Carmelo Herdes
- Department
of Chemical Engineering, University of Bath, Bath BA2 7AY, U.K
| | - Andrew J. Haslam
- Department
of Chemical Engineering, Imperial College London, London SW7 2AZ, U.K
| | - George Jackson
- Department
of Chemical Engineering, Imperial College London, London SW7 2AZ, U.K
| | - Erich A. Müller
- Department
of Chemical Engineering, Imperial College London, London SW7 2AZ, U.K
| |
Collapse
|