1
|
Ricci E, Minelli M, De Angelis MG. Modelling Sorption and Transport of Gases in Polymeric Membranes across Different Scales: A Review. MEMBRANES 2022; 12:857. [PMID: 36135877 PMCID: PMC9502097 DOI: 10.3390/membranes12090857] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/29/2022] [Revised: 08/24/2022] [Accepted: 08/27/2022] [Indexed: 06/02/2023]
Abstract
Professor Giulio C. Sarti has provided outstanding contributions to the modelling of fluid sorption and transport in polymeric materials, with a special eye on industrial applications such as membrane separation, due to his Chemical Engineering background. He was the co-creator of innovative theories such as the Non-Equilibrium Theory for Glassy Polymers (NET-GP), a flexible tool to estimate the solubility of pure and mixed fluids in a wide range of polymers, and of the Standard Transport Model (STM) for estimating membrane permeability and selectivity. In this review, inspired by his rigorous and original approach to representing membrane fundamentals, we provide an overview of the most significant and up-to-date modeling tools available to estimate the main properties governing polymeric membranes in fluid separation, namely solubility and diffusivity. The paper is not meant to be comprehensive, but it focuses on those contributions that are most relevant or that show the potential to be relevant in the future. We do not restrict our view to the field of macroscopic modelling, which was the main playground of professor Sarti, but also devote our attention to Molecular and Multiscale Hierarchical Modeling. This work proposes a critical evaluation of the different approaches considered, along with their limitations and potentiality.
Collapse
Affiliation(s)
- Eleonora Ricci
- Department of Civil, Chemical, Environmental and Materials Engineering (DICAM), Alma Mater Studiorum—University of Bologna, 40126 Bologna, Italy
| | - Matteo Minelli
- Department of Civil, Chemical, Environmental and Materials Engineering (DICAM), Alma Mater Studiorum—University of Bologna, 40126 Bologna, Italy
| | - Maria Grazia De Angelis
- Institute for Materials and Processes, School of Engineering, University of Edinburgh, Edinburgh EH9 3FB, UK
| |
Collapse
|
2
|
Voyiatzis E, Stroeks A. Atomistic Modeling of Hydrogen and Oxygen Solubility in Semicrystalline PA-6 and HDPE Materials. J Phys Chem B 2022; 126:6102-6111. [PMID: 35921684 DOI: 10.1021/acs.jpcb.2c02854] [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
Hydrogen is a clean and sustainable energy carrier which plays a major role in the transition of the global energy market to a less fossil fuel dependent future. Polymer-based materials are crucial in the production, storage, transportation, and energy extraction of hydrogen. More insights in the hydrogen-polymers interactions are required to guide material design and product development, especially for hydrogen solubility in polymers, which is crucial in many applications. The current study aims at rationalizing the determining factors of hydrogen solubility in two relevant polymers: polyamide-6 (PA-6) and high density polyethylene (HDPE). Based on atomistic molecular dynamics simulations and experimental data, we have reached several conclusions related to hydrogen and oxygen solubility in these two polymers. The crystal phases of PA-6 and HDPE are impenetrable to hydrogen and oxygen at elevated pressures, despite the small molecular size of hydrogen and oxygen. The practical implication for gas barrier applications is that polymer crystals act as impermeable obstacles and gas migration takes place primarily in the amorphous phase. Experimental hydrogen and oxygen solubilities in PA-6 and HDPE at elevated pressures can be predicted in a semiquantitative manner by molecular simulations. The discrepancies between experimental and predicted values could be attributed to neglect of the effect of crystal regions on the amorphous polymer domains. Although hydrogen is smaller than oxygen, it has been experimentally observed that hydrogen has a lower solubility in PA-6 and HDPE than oxygen. This observation has been confirmed by molecular simulations and attributed to the more favorable energetic interactions of oxygen with PA-6 and PE than of hydrogen. These interactions dominate the solubility behavior over the distribution of the accessible volume in the polymers.
Collapse
Affiliation(s)
| | - Alexander Stroeks
- DSM Engineering Materials, Urmonderbaan 22, 6167 RD Geleen, The Netherlands
| |
Collapse
|
3
|
Rahbari A, Hens R, Ramdin M, Moultos OA, Dubbeldam D, Vlugt TJH. Recent advances in the continuous fractional component Monte Carlo methodology. MOLECULAR SIMULATION 2020. [DOI: 10.1080/08927022.2020.1828585] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- A. Rahbari
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, Netherlands
| | - R. Hens
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, Netherlands
| | - M. Ramdin
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, Netherlands
| | - O. A. Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, Netherlands
| | - D. Dubbeldam
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam, Netherlands
| | - T. J. H. Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, Netherlands
| |
Collapse
|
4
|
Vergadou N, Theodorou DN. Molecular Modeling Investigations of Sorption and Diffusion of Small Molecules in Glassy Polymers. MEMBRANES 2019; 9:E98. [PMID: 31398889 PMCID: PMC6723301 DOI: 10.3390/membranes9080098] [Citation(s) in RCA: 38] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/12/2019] [Revised: 07/22/2019] [Accepted: 07/23/2019] [Indexed: 11/16/2022]
Abstract
With a wide range of applications, from energy and environmental engineering, such as in gas separations and water purification, to biomedical engineering and packaging, glassy polymeric materials remain in the core of novel membrane and state-of the art barrier technologies. This review focuses on molecular simulation methodologies implemented for the study of sorption and diffusion of small molecules in dense glassy polymeric systems. Basic concepts are introduced and systematic methods for the generation of realistic polymer configurations are briefly presented. Challenges related to the long length and time scale phenomena that govern the permeation process in the glassy polymer matrix are described and molecular simulation approaches developed to address the multiscale problem at hand are discussed.
Collapse
Affiliation(s)
- Niki Vergadou
- Molecular Thermodynamics and Modelling of Materials Laboratory, Institute of Nanoscience and Nanotechnology, National Center for Scientific Research Demokritos, Aghia Paraskevi Attikis, GR-15310 Athens, Greece.
| | - Doros N Theodorou
- School of Chemical Engineering, National Technical University of Athens, GR 15780 Athens, Greece
| |
Collapse
|
5
|
Molecular-level insight of gas transport in composite poly (4-methyl-2-pentyne) and nanoparticles of titanium dioxide. Eur Polym J 2018. [DOI: 10.1016/j.eurpolymj.2018.03.022] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
6
|
Yang Q, Achenie LEK. Exploration of bulk and interface behavior of gas molecules and 1-butyl-3-methylimidazolium tetrafluoroborate ionic liquid using equilibrium and nonequilibrium molecular dynamics simulation and quantum chemical calculation. Phys Chem Chem Phys 2018; 20:10121-10131. [PMID: 29588998 DOI: 10.1039/c7cp07714a] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Ionic liquids (ILs) show brilliant performance in separating gas impurities, but few researchers have performed an in-depth exploration of the bulk and interface behavior of penetrants and ILs thoroughly. In this research, we have performed a study on both molecular dynamics (MD) simulation and quantum chemical (QC) calculation to explore the transport of acetylene and ethylene in the bulk and interface regions of 1-butyl-3-methylimidazolium tetrafluoroborate ([BMIM]-[BF4]). The diffusivity, solubility and permeability of gas molecules in the bulk were researched with MD simulation first. The subdiffusion behavior of gas molecules is induced by coupling between the motion of gas molecules and the ions, and the relaxation processes of the ions after the disturbance caused by gas molecules. Then, QC calculation was performed to explore the optical geometry of ions, ion pairs and complexes of ions and penetrants, and interaction potential for pairs and complexes. Finally, nonequilibrium MD simulation was performed to explore the interface structure and properties of the IL-gas system and gas molecule behavior in the interface region. The research results may be used in the design of IL separation media.
Collapse
Affiliation(s)
- Quan Yang
- Department of Chemical Engineering, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA.
| | | |
Collapse
|
7
|
Yang Q, Whiting WI. Molecular-level insight of the differences in the diffusion and solubility of penetrants in polypropylene, poly(propylmethylsiloxane) and poly(4-methyl-2-pentyne). J Memb Sci 2018. [DOI: 10.1016/j.memsci.2017.12.011] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
8
|
Anderson LR, Yang Q, Ediger AM. Comparing gas transport in three polymers via molecular dynamics simulation. Phys Chem Chem Phys 2018; 20:22123-22133. [PMID: 30113613 DOI: 10.1039/c8cp02829j] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
Molecular dynamics (MD) simulation was employed to study the transport of methane and n-butane molecules in the bulk and interface region of polyethylene (PE), poly(4-methyl-2-pentyne) (PMP) and polydimethylsiloxane (PDMS).
Collapse
Affiliation(s)
- Luke R. Anderson
- Department of Material Science and Engineering
- Virginia Polytechnic Institute and State University
- Blacksburg
- USA
| | - Quan Yang
- Department of Chemical Engineering
- Virginia Polytechnic Institute and State University
- Blacksburg
- USA
- Sandia National Laboratories
| | - Andrew M. Ediger
- Department of Material Science and Engineering
- Virginia Polytechnic Institute and State University
- Blacksburg
- USA
| |
Collapse
|
9
|
Akkermans RLC. Solvation Free Energy of Regular and Azeotropic Molecular Mixtures. J Phys Chem B 2017; 121:1675-1683. [DOI: 10.1021/acs.jpcb.7b00125] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Reinier L. C. Akkermans
- Dassault Systèmes, BIOVIA Ltd., 334 Cambridge Science Park, Cambridge, CB4 0WN, United Kingdom
| |
Collapse
|
10
|
Donnini S, Ullmann RT, Groenhof G, Grubmüller H. Charge-Neutral Constant pH Molecular Dynamics Simulations Using a Parsimonious Proton Buffer. J Chem Theory Comput 2016; 12:1040-51. [PMID: 26881315 DOI: 10.1021/acs.jctc.5b01160] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In constant pH molecular dynamics simulations, the protonation states of titratable sites can respond to changes of the pH and of their electrostatic environment. Consequently, the number of protons bound to the biomolecule, and therefore the overall charge of the system, fluctuates during the simulation. To avoid artifacts associated with a non-neutral simulation system, we introduce an approach to maintain neutrality of the simulation box in constant pH molecular dynamics simulations, while maintaining an accurate description of all protonation fluctuations. Specifically, we introduce a proton buffer that, like a buffer in experiment, can exchange protons with the biomolecule enabling its charge to fluctuate. To keep the total charge of the system constant, the uptake and release of protons by the buffer are coupled to the titration of the biomolecule with a constraint. We find that, because the fluctuation of the total charge (number of protons) of a typical biomolecule is much smaller than the number of titratable sites of the biomolecule, the number of buffer sites required to maintain overall charge neutrality without compromising the charge fluctuations of the biomolecule, is typically much smaller than the number of titratable sites, implying markedly enhanced simulation and sampling efficiency.
Collapse
Affiliation(s)
- Serena Donnini
- Nanoscience Center and Department of Biological and Environmental Sciences, University of Jyväskylä , P. O. Box 35, 40014 Jyväskylä, Finland
| | - R Thomas Ullmann
- Department of Theoretical and Computational Biophysics, Max Planck Institute for Biophysical Chemistry , Am Faßberg 11, 37077 Göttingen, Germany
| | - Gerrit Groenhof
- Nanoscience Center and Department of Chemistry, University of Jyväskylä , P. O. Box 35, 40014 Jyväskylä, Finland
| | - Helmut Grubmüller
- Department of Theoretical and Computational Biophysics, Max Planck Institute for Biophysical Chemistry , Am Faßberg 11, 37077 Göttingen, Germany
| |
Collapse
|
11
|
Yang Q, Achenie LE, Cai W. Comparing Penetrants Transport in Composite Poly (4-methyl-2-pentyne) and Nanoparticles of Cristobalite Silica and Faujasite Silica through Molecular Dynamics Simulation. Ind Eng Chem Res 2013. [DOI: 10.1021/ie400524k] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Affiliation(s)
- Quan Yang
- Department
of Chemical Engineering, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061, United
States
| | - Luke E. Achenie
- Department
of Chemical Engineering, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061, United
States
| | - Weibin Cai
- Honghao Qingyuan Technology Corp. Ltd., Beijing, 102218, China
| |
Collapse
|
12
|
Donnini S, Tegeler F, Groenhof G, Grubmüller H. Constant pH Molecular Dynamics in Explicit Solvent with λ-Dynamics. J Chem Theory Comput 2011; 7:1962-1978. [PMID: 21687785 PMCID: PMC3114466 DOI: 10.1021/ct200061r] [Citation(s) in RCA: 142] [Impact Index Per Article: 10.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2011] [Indexed: 12/22/2022]
Abstract
pH is an important parameter in condensed-phase systems, because it determines the protonation state of titratable groups and thus influences the structure, dynamics, and function of molecules in solution. In most force field simulation protocols, however, the protonation state of a system (rather than its pH) is kept fixed and cannot adapt to changes of the local environment. Here, we present a method, implemented within the MD package GROMACS, for constant pH molecular dynamics simulations in explicit solvent that is based on the λ-dynamics approach. In the latter, the dynamics of the titration coordinate λ, which interpolates between the protonated and deprotonated states, is driven by generalized forces between the protonated and deprotonated states. The hydration free energy, as a function of pH, is included to facilitate constant pH simulations. The protonation states of titratable groups are allowed to change dynamically during a simulation, thus reproducing average protonation probabilities at a certain pH. The accuracy of the method is tested against titration curves of single amino acids and a dipeptide in explicit solvent.
Collapse
Affiliation(s)
- Serena Donnini
- Department of Theoretical and Computational Biophysics, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| | | | - Gerrit Groenhof
- Department of Theoretical and Computational Biophysics, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| | - Helmut Grubmüller
- Department of Theoretical and Computational Biophysics, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| |
Collapse
|
13
|
Hörstermann H, Hentschke R, Amkreutz M, Hoffmann M, Wirts-Rütters M. Predicting water sorption and volume swelling in dense polymer systems via computer simulation. J Phys Chem B 2010; 114:17013-24. [PMID: 21141921 DOI: 10.1021/jp105210y] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Atomistic model structures of amorphous polyamide 6 (PA-6) and of an adhesive system consisting of the diglycidyl ether of bisphenol A (DGEBA) as epoxy resin and isophorone diamine (IPD) as a curing agent are generated. For the adhesive, we use a new approach for the generation of the cross-linked polymer networks. It takes into account the chemical reaction kinetics of the curing reaction and, therefore, results in more realistic network structures. On the basis of the corresponding model structures, the equilibrium water content and the swelling ratio of amorphous PA-6 and of the DGEBA+IPD networks are calculated via computer simulation for different thermodynamic conditions. A hybrid method is used combining the molecular dynamics technique with an accelerated test particle insertion method. The results are in reasonable agreement with experiments and, in the case of the PA-6 system, with results obtained via other computer simulation methods.
Collapse
Affiliation(s)
- Henning Hörstermann
- Fachbereich Mathematik und Naturwissenschaften, Bergische Universität, D-42097 Wuppertal, Germany
| | | | | | | | | |
Collapse
|
14
|
Gillet G, Vitrac O, Desobry S. Prediction of Partition Coefficients of Plastic Additives between Packaging Materials and Food Simulants. Ind Eng Chem Res 2010. [DOI: 10.1021/ie9010595] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Guillaume Gillet
- Laboratoire National de métrologie et d’Essais, Centre Energie, Matériaux et Emballage, 29 avenue Roger Hennequin, 78197 Trappes CEDEX, France, Institut National de la Recherche Agronomique, Joint Research Unit 1145 “Food Process Engineering” between INRA, Agroparistech and CNAM, 1 avenue des Olympiades, 91300 Massy, France, and Nancy Université, LSGA-ENSAIA-INPL, 2 avenue de la forêt de Haye, BP 172, 54505 Vandoeuvre lès Nancy, France
| | - Olivier Vitrac
- Laboratoire National de métrologie et d’Essais, Centre Energie, Matériaux et Emballage, 29 avenue Roger Hennequin, 78197 Trappes CEDEX, France, Institut National de la Recherche Agronomique, Joint Research Unit 1145 “Food Process Engineering” between INRA, Agroparistech and CNAM, 1 avenue des Olympiades, 91300 Massy, France, and Nancy Université, LSGA-ENSAIA-INPL, 2 avenue de la forêt de Haye, BP 172, 54505 Vandoeuvre lès Nancy, France
| | - Stéphane Desobry
- Laboratoire National de métrologie et d’Essais, Centre Energie, Matériaux et Emballage, 29 avenue Roger Hennequin, 78197 Trappes CEDEX, France, Institut National de la Recherche Agronomique, Joint Research Unit 1145 “Food Process Engineering” between INRA, Agroparistech and CNAM, 1 avenue des Olympiades, 91300 Massy, France, and Nancy Université, LSGA-ENSAIA-INPL, 2 avenue de la forêt de Haye, BP 172, 54505 Vandoeuvre lès Nancy, France
| |
Collapse
|
15
|
|
16
|
Gillet G, Vitrac O, Desobry S. Prediction of Solute Partition Coefficients between Polyolefins and Alcohols Using a Generalized Flory−Huggins Approach. Ind Eng Chem Res 2009. [DOI: 10.1021/ie801141h] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Guillaume Gillet
- Centre Energie, Matériaux et Emballage, Laboratoire National de métrologie et d’Essais, 29 avenue Roger Hennequin, 78197 Trappes Cedex, France, LSGA-ENSAIA-INPL, Nancy Université, 2 avenue de la forêt de Haye, BP 172, 54505 Vandoeuvre lès Nancy, France, and UMR 1145 Génie Industriel Alimentaire, Institut National de la Recherche Agronomique, 1 avenue des Olympiades, 91300 Massy, France
| | - Olivier Vitrac
- Centre Energie, Matériaux et Emballage, Laboratoire National de métrologie et d’Essais, 29 avenue Roger Hennequin, 78197 Trappes Cedex, France, LSGA-ENSAIA-INPL, Nancy Université, 2 avenue de la forêt de Haye, BP 172, 54505 Vandoeuvre lès Nancy, France, and UMR 1145 Génie Industriel Alimentaire, Institut National de la Recherche Agronomique, 1 avenue des Olympiades, 91300 Massy, France
| | - Stéphane Desobry
- Centre Energie, Matériaux et Emballage, Laboratoire National de métrologie et d’Essais, 29 avenue Roger Hennequin, 78197 Trappes Cedex, France, LSGA-ENSAIA-INPL, Nancy Université, 2 avenue de la forêt de Haye, BP 172, 54505 Vandoeuvre lès Nancy, France, and UMR 1145 Génie Industriel Alimentaire, Institut National de la Recherche Agronomique, 1 avenue des Olympiades, 91300 Massy, France
| |
Collapse
|
17
|
Spyriouni T, Boulougouris GC, Theodorou DN. Prediction of Sorption of CO2 in Glassy Atactic Polystyrene at Elevated Pressures Through a New Computational Scheme. Macromolecules 2009. [DOI: 10.1021/ma8015294] [Citation(s) in RCA: 46] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Theodora Spyriouni
- Institute of Physical Chemistry, NCSR “Demokritos”, Aghia Paraskevi Attikis, 153 10 Athens, Greece, Scienomics SARL, 17, Square Eduard VII, 75009 Paris, France, and Department of Materials Science and Engineering, School of Chemical Engineering, National Technical University of Athens, 9 Heroon Polytechniou, Zografou Campus, Athens 15780, Greece
| | - Georgios C. Boulougouris
- Institute of Physical Chemistry, NCSR “Demokritos”, Aghia Paraskevi Attikis, 153 10 Athens, Greece, Scienomics SARL, 17, Square Eduard VII, 75009 Paris, France, and Department of Materials Science and Engineering, School of Chemical Engineering, National Technical University of Athens, 9 Heroon Polytechniou, Zografou Campus, Athens 15780, Greece
| | - Doros N. Theodorou
- Institute of Physical Chemistry, NCSR “Demokritos”, Aghia Paraskevi Attikis, 153 10 Athens, Greece, Scienomics SARL, 17, Square Eduard VII, 75009 Paris, France, and Department of Materials Science and Engineering, School of Chemical Engineering, National Technical University of Athens, 9 Heroon Polytechniou, Zografou Campus, Athens 15780, Greece
| |
Collapse
|
18
|
Özal TA, Peter C, Hess B, van der Vegt NFA. Modeling Solubilities of Additives in Polymer Microstructures: Single-Step Perturbation Method Based on a Soft-Cavity Reference State. Macromolecules 2008. [DOI: 10.1021/ma702329q] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Affiliation(s)
- Tugba A. Özal
- Max Planck Institute for Polymer Research, Ackermannweg 10, D-55128 Mainz, Germany
| | - Christine Peter
- Max Planck Institute for Polymer Research, Ackermannweg 10, D-55128 Mainz, Germany
| | - Berk Hess
- Max Planck Institute for Polymer Research, Ackermannweg 10, D-55128 Mainz, Germany
| | | |
Collapse
|
19
|
Hess B, Peter C, Ozal T, van der Vegt NFA. Fast-Growth Thermodynamic Integration: Calculating Excess Chemical Potentials of Additive Molecules in Polymer Microstructures. Macromolecules 2008. [DOI: 10.1021/ma702070n] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Berk Hess
- Max-Planck Institute for Polymer Research, Ackermannweg 10, D-55128 Mainz, Germany
| | - Christine Peter
- Max-Planck Institute for Polymer Research, Ackermannweg 10, D-55128 Mainz, Germany
| | - Tugba Ozal
- Max-Planck Institute for Polymer Research, Ackermannweg 10, D-55128 Mainz, Germany
| | | |
Collapse
|
20
|
Theodorou DN. A reversible minimum-to-minimum mapping method for the calculation of free-energy differences. J Chem Phys 2007; 124:034109. [PMID: 16438569 DOI: 10.1063/1.2138701] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A general method is introduced for the calculation of the free-energy difference between two systems, 0 and 1, with configuration spaces omega(0), omega(1) of the same dimensionality. The method relies upon establishing a objective mapping between disjoint subsets gamma(i)(0) of omega(0) and corresponding disjoint subsets gamma(i)(1) of omega(1), and averaging a function of the ratio of configurational integrals over gamma(i)(0) and gamma(i)(1) with respect to the probability densities of the two systems. The mapped subsets gamma(i)(0) and gamma(i)(1) need not span the entire configuration spaces omega(0) and omega(1). The method is applied for the calculation of the excess chemical potential mu(ex) in a Lennard-Jones (LJ) fluid. In this case, omega(0) is the configuration space of a (N-1) real molecule plus one ideal-gas molecule system, while omega(1) is the configuration space of a N real molecule system occupying the same volume. Gamma(i)(0) and gamma(i)(1) are constructed from hyperspheres of the same radius centered at minimum-energy configurations of a set of "active" molecules lying within distance a from the ideal-gas molecule and the last real molecule, respectively. An algorithm is described for sampling gamma(i)(0) and gamma(i)(1) given a point in omega(0) or in omega(1). The algorithm encompasses three steps: "quenching" (minimization with respect to the active-molecule degrees of freedom), "mutation" (gradual conversion of the ideal-gas molecule into a real molecule, with simultaneous minimization of the energy with respect to the active-molecule degrees of freedom), and "excitation" (generation of points on a hypersphere centered at the active-molecule energy minimum). These steps are also carried out in reverse, as required by the bijective nature of the mapping. The mutation step, which establishes a reversible mapping between energy minima with respect to the active degrees of freedom of systems 0 and 1, ensures that excluded volume interactions emerging in the process of converting the ideal-gas molecule into a real molecule are relieved through appropriate rearrangement of the surrounding active molecules. Thus, the insertion problem plaguing traditional methods for the calculation of chemical potential at high densities is alleviated. Results are presented at two state points of the LJ system for a variety of radii a of the active domain. It is shown that the estimated values of mu(ex) are correct in all cases and subject to an order of magnitude lower statistical uncertainty than values based on the same number of Widom [J. Chem. Phys. 39, 2808 (1963)] insertions at high fluid densities. Optimal settings for the new algorithm are identified and distributions of the quantities involved in it [number of active molecules, energy at the sampled minima of systems 0 and 1, and free-energy differences between subsets gamma(i)(0) and gamma(i)(1) that are mapped onto each other] are explored.
Collapse
Affiliation(s)
- Doros N Theodorou
- School of Chemical Engineering, National Technical University of Athens, 9 Heroon Polytechniou Street, Zografou Campus, 157 80 Athens, Greece.
| |
Collapse
|
21
|
Abstract
A replica exchange method is presented which requires fewer replicas and is designed to be used for large systems. In this method, dynamically scaled replicas are placed between conventional replicas at broadly spaced temperatures. The potential of the scaled replicas is linearly scaled by a dynamical variable which varies between 0 and 1. When the variable is near either end point the replica can undergo exchanges with one of its neighboring replicas. Two different versions of the method are presented for a model system of a small peptide in water. The scaled replica can replace many replicas and the method can be up to ten times more efficient than conventional replica exchange.
Collapse
Affiliation(s)
- Steven W Rick
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, USA.
| |
Collapse
|
22
|
Siegert MR, Heuchel M, Hofmann D. A generalized direct-particle-deletion scheme for the calculation of chemical potential and solubilities of small- and medium-sized molecules in amorphous polymers. J Comput Chem 2007; 28:877-89. [PMID: 17238170 DOI: 10.1002/jcc.20594] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
The development, validation, and first applications of a generalized version of an inverse Widom method are described. It permits the calculation of solubility coefficients for molecules as large as, e.g., benzene in all polymers for which reasonable forcefield parameters exist. Predicting the solubility is a key to the knowledge-based design of materials utilized to solve permeability related problems. For long time, particle insertion methods, such as the Widom method, were the only way to predict solubilities from molecular models, but they, in most cases, only worked well for rather small penetrants (e.g., H2, O2, N2). Therefore, a few years ago, a new particle deletion algorithm "DPD" was introduced by Boulougouris, Economou, and Theodorou to overcome this problem in principle. The related computer code was, however, only applicable to special, relatively simple model systems. As application examples for the generalized version described here, solubility calculations for nitrogen, oxygen, and benzene in poly(dimethyl siloxane) are presented.
Collapse
Affiliation(s)
- Martin R Siegert
- GKSS Research Center, Institute of Polymer Research, Kantstrasse 55, 14513 Teltow, Germany
| | | | | |
Collapse
|
23
|
Kuwajima S, Kikuchi H, Fukuda M. Molecular-dynamics evaluation of fluid-phase equilibrium properties by a novel free-energy perturbation approach: Application to gas solubility and vapor pressure of liquid hexane. J Chem Phys 2006; 124:124111. [PMID: 16599666 DOI: 10.1063/1.2178321] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A novel free-energy perturbation method is developed for the computation of the free energy of transferring a molecule between fluid phases. The methodology consists in drawing a free-energy profile of the target molecule moving across a binary-phase structure built in the computer. The novelty of the method lies in the difference of the definition of the free-energy profile from the common definition. As an important element of the method, the process of making a correction to the transfer free energy with respect to the cutoff of intermolecular forces is elucidated. In order to examine the performance of the method in the application to fluid-phase equilibrium properties, molecular-dynamics computations are carried out for the evaluation of gas solubility and vapor pressure of liquid n-hexane at 298.15 K. The gas species treated are methane, ethane, propane, and n-butane, with the gas solubility expressed as Henry's constant. It is shown that the method works fine and calculated results are generally in good agreement with experiments. It is found that the cutoff correction is strikingly large, constituting a dominant part of the calculated transfer free energy at the cutoff of 8 A.
Collapse
Affiliation(s)
- Satoru Kuwajima
- NanoSimulation Associates, 825-1 Amado-cho, Villa DE 201, Hanamigawa-ku, Chiba-shi, Chiba 262-0043, Japan.
| | | | | |
Collapse
|
24
|
Kikuchi H, Kuwajima S, Fukuda M. Novel method to estimate solubility of small molecules in cis-polyisoprene by molecular dynamics simulations. J Chem Phys 2001. [DOI: 10.1063/1.1398590] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
25
|
van der Vegt NFA. Temperature Dependence of Gas Transport in Polymer Melts: Molecular Dynamics Simulations of CO2 in Polyethylene. Macromolecules 2000. [DOI: 10.1021/ma991737f] [Citation(s) in RCA: 50] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
26
|
Fukuda M. Solubilities of small molecules in polyethylene evaluated by a test-particle-insertion method. J Chem Phys 2000. [DOI: 10.1063/1.480594] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
27
|
van der Vegt NFA, Briels WJ, Wessling M, Strathmann H. The sorption induced glass transition in amorphous glassy polymers. J Chem Phys 1999. [DOI: 10.1063/1.479042] [Citation(s) in RCA: 54] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|