1
|
Jin J, Pak AJ, Durumeric AEP, Loose TD, Voth GA. Bottom-up Coarse-Graining: Principles and Perspectives. J Chem Theory Comput 2022; 18:5759-5791. [PMID: 36070494 PMCID: PMC9558379 DOI: 10.1021/acs.jctc.2c00643] [Citation(s) in RCA: 72] [Impact Index Per Article: 36.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2022] [Indexed: 01/14/2023]
Abstract
Large-scale computational molecular models provide scientists a means to investigate the effect of microscopic details on emergent mesoscopic behavior. Elucidating the relationship between variations on the molecular scale and macroscopic observable properties facilitates an understanding of the molecular interactions driving the properties of real world materials and complex systems (e.g., those found in biology, chemistry, and materials science). As a result, discovering an explicit, systematic connection between microscopic nature and emergent mesoscopic behavior is a fundamental goal for this type of investigation. The molecular forces critical to driving the behavior of complex heterogeneous systems are often unclear. More problematically, simulations of representative model systems are often prohibitively expensive from both spatial and temporal perspectives, impeding straightforward investigations over possible hypotheses characterizing molecular behavior. While the reduction in resolution of a study, such as moving from an atomistic simulation to that of the resolution of large coarse-grained (CG) groups of atoms, can partially ameliorate the cost of individual simulations, the relationship between the proposed microscopic details and this intermediate resolution is nontrivial and presents new obstacles to study. Small portions of these complex systems can be realistically simulated. Alone, these smaller simulations likely do not provide insight into collectively emergent behavior. However, by proposing that the driving forces in both smaller and larger systems (containing many related copies of the smaller system) have an explicit connection, systematic bottom-up CG techniques can be used to transfer CG hypotheses discovered using a smaller scale system to a larger system of primary interest. The proposed connection between different CG systems is prescribed by (i) the CG representation (mapping) and (ii) the functional form and parameters used to represent the CG energetics, which approximate potentials of mean force (PMFs). As a result, the design of CG methods that facilitate a variety of physically relevant representations, approximations, and force fields is critical to moving the frontier of systematic CG forward. Crucially, the proposed connection between the system used for parametrization and the system of interest is orthogonal to the optimization used to approximate the potential of mean force present in all systematic CG methods. The empirical efficacy of machine learning techniques on a variety of tasks provides strong motivation to consider these approaches for approximating the PMF and analyzing these approximations.
Collapse
Affiliation(s)
- Jaehyeok Jin
- Department of Chemistry,
Chicago Center for Theoretical Chemistry, Institute for Biophysical
Dynamics, and James Franck Institute, The
University of Chicago, Chicago, Illinois 60637, United States
| | - Alexander J. Pak
- Department of Chemistry,
Chicago Center for Theoretical Chemistry, Institute for Biophysical
Dynamics, and James Franck Institute, The
University of Chicago, Chicago, Illinois 60637, United States
| | - Aleksander E. P. Durumeric
- Department of Chemistry,
Chicago Center for Theoretical Chemistry, Institute for Biophysical
Dynamics, and James Franck Institute, The
University of Chicago, Chicago, Illinois 60637, United States
| | - Timothy D. Loose
- Department of Chemistry,
Chicago Center for Theoretical Chemistry, Institute for Biophysical
Dynamics, and James Franck Institute, The
University of Chicago, Chicago, Illinois 60637, United States
| | - Gregory A. Voth
- Department of Chemistry,
Chicago Center for Theoretical Chemistry, Institute for Biophysical
Dynamics, and James Franck Institute, The
University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
2
|
Ge Y, Zhu Q, Li Y, Dong H, Ma J. An electrostatic-variable coarse-grained model for predicting enthalpy of vaporization, surface tension, diffusivity, conductivity, and dielectric constant of aqueous ionic liquid. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2021.118230] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
3
|
Wu C, Li K, Ning X, Zhang L. An Enhanced Scheme for Multiscale Modeling of Thermomechanical Properties of Polymer Bulks. J Phys Chem B 2021; 125:8612-8626. [PMID: 34291641 DOI: 10.1021/acs.jpcb.1c02663] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
While multiscale modeling significantly enhances the capability of molecular simulations of polymer systems, it is well realized that the systematically derived coarse-grained (CG) models generally underestimate the thermomechanical properties. In this work, a charge-based mapping scheme has been adopted to include explicit electrostatic interactions and benchmarked against two typical polymers, atactic poly(methyl methacrylate) (PMMA) and polystyrene (PS). The CG potentials are parameterized against the oligomer bulks of nine monomers per chain to match the essential structural features and the two basic pressure-volume-temperature (PVT) properties, which are obtained from the all-atomistic (AA) molecular dynamics (MD) simulations at a single elevated temperature. The so-parameterized CG potentials are extended with the MD method to simulate the two polymer bulks of one hundred monomers per chain over a wide temperature range. Without any scaling, all the simulated results, including mass densities and bulk moduli at room temperature, thermal expansion coefficients at rubbery and glassy states, and glass transition temperatures (Tg), compare well with the corresponding experimental data. The proposed scheme not only contributes to realistically simulating various thermomechanical properties of both apolar and polar polymers but also allows for directly simulating their electrical properties.
Collapse
Affiliation(s)
- Chaofu Wu
- Hunan Provincial Key Laboratory of Fine Ceramics and Powder Materials, School of Materials and Environmental Engineering, Hunan University of Humanities, Science and Technology, Loudi, Hunan 417000, P. R. China
| | - Kewen Li
- Hunan Provincial Key Laboratory of Fine Ceramics and Powder Materials, School of Materials and Environmental Engineering, Hunan University of Humanities, Science and Technology, Loudi, Hunan 417000, P. R. China
| | - Xutao Ning
- Hunan Provincial Key Laboratory of Fine Ceramics and Powder Materials, School of Materials and Environmental Engineering, Hunan University of Humanities, Science and Technology, Loudi, Hunan 417000, P. R. China
| | - Lei Zhang
- Hunan Provincial Key Laboratory of Fine Ceramics and Powder Materials, School of Materials and Environmental Engineering, Hunan University of Humanities, Science and Technology, Loudi, Hunan 417000, P. R. China
| |
Collapse
|
4
|
Rudzinski JF, Kloth S, Wörner S, Pal T, Kremer K, Bereau T, Vogel M. Dynamical properties across different coarse-grained models for ionic liquids. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2021; 33:224001. [PMID: 33592598 DOI: 10.1088/1361-648x/abe6e1] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/30/2020] [Accepted: 02/16/2021] [Indexed: 06/12/2023]
Abstract
Room-temperature ionic liquids (RTILs) stand out among molecular liquids for their rich physicochemical characteristics, including structural and dynamic heterogeneity. The significance of electrostatic interactions in RTILs results in long characteristic length- and timescales, and has motivated the development of a number of coarse-grained (CG) simulation models. In this study, we aim to better understand the connection between certain CG parameterization strategies and the dynamical properties and transferability of the resulting models. We systematically compare five CG models: a model largely parameterized from experimental thermodynamic observables; a refinement of this model to increase its structural accuracy; and three models that reproduce a given set of structural distribution functions by construction, with varying intramolecular parameterizations and reference temperatures. All five CG models display limited structural transferability over temperature, and also result in various effective dynamical speedup factors, relative to a reference atomistic model. On the other hand, the structure-based CG models tend to result in more consistent cation-anion relative diffusion than the thermodynamic-based models, for a single thermodynamic state point. By linking short- and long-timescale dynamical behaviors, we demonstrate that the varying dynamical properties of the different CG models can be largely collapsed onto a single curve, which provides evidence for a route to constructing dynamically-consistent CG models of RTILs.
Collapse
Affiliation(s)
| | - Sebastian Kloth
- Institute of Condensed Matter Physics, Technische Universität Darmstadt, Hochschulstr. 6, 64289 Darmstadt, Germany
| | - Svenja Wörner
- Max Planck Institute for Polymer Research, 55128 Mainz, Germany
| | - Tamisra Pal
- Institute of Condensed Matter Physics, Technische Universität Darmstadt, Hochschulstr. 6, 64289 Darmstadt, Germany
| | - Kurt Kremer
- Max Planck Institute for Polymer Research, 55128 Mainz, Germany
| | - Tristan Bereau
- Max Planck Institute for Polymer Research, 55128 Mainz, Germany
- Van 't Hoff Institute for Molecular Sciences and Informatics Institute, University of Amsterdam, Amsterdam 1098 XH, The Netherlands
| | - Michael Vogel
- Institute of Condensed Matter Physics, Technische Universität Darmstadt, Hochschulstr. 6, 64289 Darmstadt, Germany
| |
Collapse
|
5
|
Kloth S, Bernhardt MP, van der Vegt NFA, Vogel M. Coarse-grained model of a nanoscale-segregated ionic liquid for simulations of low-temperature structure and dynamics. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2021; 33:204002. [PMID: 33578395 DOI: 10.1088/1361-648x/abe606] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/21/2020] [Accepted: 02/12/2021] [Indexed: 06/12/2023]
Abstract
We perform molecular dynamics simulations to study the structure and dynamics of the ionic liquid [Omim][TFSI] in a broad temperature range. A particular focus is the progressing nanoscale segregation into polar and nonpolar regions upon cooling. As this analysis requires simulations of large systems for long times, we use the iterative Boltzmann inversion method to develop a new coarse-grained (CG) model from a successful all-atom (AA) model. We show that the properties are similar for both levels of description at room temperature, while the CG model shows stronger nanoscale segregation and faster diffusion dynamics than its AA counterpart at low temperatures. Exploiting these features of the CG model, we find that the characteristic length scale of the structural inhomogeneity nearly doubles to ∼3 nm when the temperature is decreased to about 200 K. Moreover, we observe that the nanoscale segregation is characterized by a bicontinuous morphology. In worm-like nonpolar regions, the ends of the octyl rests of the cations preferentially aggregate in the centers, while the other parts of the alkyl chains tend to be aligned parallel on a next-neighbor level and point outward, allowing for an integration of the imidazolium head groups of the cations into polar regions together with the anions, resembling to some degree the molecular arrangement in cylindrical micelles.
Collapse
Affiliation(s)
- Sebastian Kloth
- Institute of Condensed Matter Physics, Technische Universität Darmstadt, Hochschulstr. 6, 64289 Darmstadt, Germany
| | - Marvin P Bernhardt
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, 64287 Darmstadt, Germany
| | - Nico F A van der Vegt
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, 64287 Darmstadt, Germany
| | - Michael Vogel
- Institute of Condensed Matter Physics, Technische Universität Darmstadt, Hochschulstr. 6, 64289 Darmstadt, Germany
| |
Collapse
|
6
|
Wen C, Odle R, Cheng S. Coarse-Grained Molecular Dynamics Modeling of a Branched Polyetherimide. Macromolecules 2021. [DOI: 10.1021/acs.macromol.0c01440] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Affiliation(s)
- Chengyuan Wen
- Department of Physics, Center for Soft Matter and Biological Physics, and Macromolecules Innovation Institute, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061, United States
| | - Roy Odle
- SABIC, 1 Lexan Lane, Mt. Vernon, Indiana 47620, United States
| | - Shengfeng Cheng
- Department of Physics, Center for Soft Matter and Biological Physics, and Macromolecules Innovation Institute, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061, United States
- Department of Mechanical Engineering, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061, United States
| |
Collapse
|