1
|
Lin TW, Sing CE. Effect of penetrant-polymer interactions and shape on the motion of molecular penetrants in dense polymer networks. J Chem Phys 2024; 160:114905. [PMID: 38511661 DOI: 10.1063/5.0197140] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2024] [Accepted: 03/01/2024] [Indexed: 03/22/2024] Open
Abstract
The diffusion of dilute molecular penetrants within polymers plays a crucial role in the advancement of material engineering for applications such as coatings and membrane separations. The potential of highly cross-linked polymer networks in these applications stems from their capacity to adjust the size and shape selectivity through subtle changes in network structures. In this paper, we use molecular dynamics simulation to understand the role of penetrant shape (aspect ratios) and its interaction with polymer networks on its diffusivity. We characterize both local penetrant hopping and the long-time diffusive motion for penetrants and consider different aspect ratios and penetrant-network interaction strengths at a variety of cross-link densities and temperatures. The shape affects the coupling of penetrant motion to the cross-link density- and temperature-dependent structural relaxation of networks and also affects the way a penetrant experiences the confinement from the network meshes. The attractive interaction between the penetrant and network primarily affects the former since only the system of dilute limit is of present interest. These results offer fundamental insights into the intricate interplay between penetrant characteristics and polymer network properties and also suggest future directions for manipulating polymer design to enhance the separation efficiency.
Collapse
Affiliation(s)
- Tsai-Wei Lin
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
- Materials Research Laboratory, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Charles E Sing
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
- Materials Research Laboratory, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
2
|
Mei B, Schweizer KS. Penetrant shape effects on activated dynamics and selectivity in polymer melts and networks based on self-consistent cooperative hopping theory. SOFT MATTER 2023; 19:8744-8763. [PMID: 37937332 DOI: 10.1039/d3sm01139a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/09/2023]
Abstract
We generalize and apply the microscopic self-consistent cooperative hopping theory for activated penetrant dynamics in polymer melts and crosslinked networks to address the role of highly variable non-spherical molecular shape. The focus is on vastly different shaped penetrants that have identical space filling volume in order to isolate how non-spherical shape explicitly modifies dynamics over a wide range of temperature down to the kinetic glass transition temperature. The theory relates intramolecular and intermolecular structure and kinetic constraints, and reveals how different solvation packing of polymer monomers around variable shaped penetrants impact penetrant hopping. A highly shape-dependent penetrant activated relaxation, including alpha time decoupling and trajectory level cooperativity of the hopping process, is predicted in the deeply supercooled regime for relatively larger penetrants which is sensitive to whether the polymer matrix is a melt or heavily crosslinked network. In contrast, for smaller size penetrants or at high/medium temperatures the effect of isochoric penetrant shape is relatively weak. We propose an aspect ratio variable that organizes how penetrant shape influences the activated relaxation times, leading to a (near) collapse or master curve. The relative absolute values of the penetrant relaxation time (inverse hopping rate) in polymer melts versus in crosslinked networks are found to be opposite when compared at a common absolute temperature versus when they are compared at a fixed value of distance from the glass transition based on the variable Tg/T with Tg the glass transition temperature. Quantitative comparison with recent diffusion experiments on chemically complex molecular penetrants of variable shape but fixed volume in crosslinked networks reveals good agreement, and testable new predictions are made. Extension of the theoretical approach to more complex systems of high experimental interest are discussed, including applications to realize selective transport in membrane separation applications.
Collapse
Affiliation(s)
- Baicheng Mei
- Department of Materials Science, University of Illinois, Urbana, IL 61801, USA.
- Materials Research Laboratory, University of Illinois, Urbana, IL 61801, USA
| | - Kenneth S Schweizer
- Department of Materials Science, University of Illinois, Urbana, IL 61801, USA.
- Department of Chemistry, University of Illinois, Urbana, IL 61801, USA
- Department of Chemical and Biomolecular Engineering, University of Illinois, Urbana, IL 61801, USA
- Materials Research Laboratory, University of Illinois, Urbana, IL 61801, USA
| |
Collapse
|
3
|
Ma XJ, Zhang R. Cooperative activated hopping dynamics in binary glass-forming liquids: effects of the size ratio, composition, and interparticle interactions. SOFT MATTER 2023. [PMID: 37317997 DOI: 10.1039/d3sm00312d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
Slow dynamics in supercooled and glassy liquids is an important research topic in soft matter physics. Compared to the traditionally focused one-component systems, glassy dynamics in mixture systems adds in a rich set of new complexities, which are fundamentally interesting and also relevant for many technological applications. In this paper, we apply the recently developed self-consistent cooperative hopping theory (SCCHT) to systematically investigate the effects of the size ratio, composition and interparticle interactions on the cooperative activated hopping dynamics of matrix (in larger size) and penetrant (in smaller size) particles in varied binary sphere mixture model systems, with a specific focus on ultrahigh mixture packing fractions that mimic the deeply supercooled glass transition conditions for molecular/polymeric mixture materials. Analysis shows that in these high activation barrier cases, the long-range elastic distortion associated with a matrix particle hopping over its cage confinement always generates an elastic barrier of a nonnegligible magnitude, although the ratio between the elastic barrier and local barrier contribution is sensitively dependent on all three mixture-specific system factors considered in this work. SCCHT predicts two general scenarios of penetrant-matrix cooperative activated hopping dynamics: matrix/penetrant co-hopping (regime 1) or the penetrant mean barrier hopping time shorter than that of the matrix (regime 2). Increasing the penetrant-to-matrix size ratio or the penetrant-matrix cross-attraction strength is found to universally enlarge the composition window of regime 1. Diverse dynamical properties characterising different aspects of the cooperative activated hopping process, including the penetrant and matrix transient localization lengths, penetrant and matrix hopping jump distances, different types of local and elastic activated barriers, and matrix long-time diffusivity, relaxation time and dynamic fragility are quantitatively studied against a wide range of variations over the three system factors. Of particular interest is the universal "anti-plasticization" phenomenon achievable for sufficiently strong cross-attractive interactions. The prospects this work opens for the exploration of a wide variety of polymer-based mixture materials are briefly discussed at the end.
Collapse
Affiliation(s)
- Xiao-Juan Ma
- South China Advanced Institute for Soft Matter Science and Technology, School of Emergent Soft Matter, South China University of Technology, Guangzhou 510640, China.
- Guangdong Provincial Key Laboratory of Functional and Intelligent Hybrid Materials and Devices, South China University of Technology, Guangzhou 510640, China
| | - Rui Zhang
- South China Advanced Institute for Soft Matter Science and Technology, School of Emergent Soft Matter, South China University of Technology, Guangzhou 510640, China.
- Guangdong Provincial Key Laboratory of Functional and Intelligent Hybrid Materials and Devices, South China University of Technology, Guangzhou 510640, China
| |
Collapse
|
4
|
Mei B, Lin TW, Sheridan GS, Evans CM, Sing CE, Schweizer KS. How Segmental Dynamics and Mesh Confinement Determine the Selective Diffusivity of Molecules in Cross-Linked Dense Polymer Networks. ACS CENTRAL SCIENCE 2023; 9:508-518. [PMID: 36968535 PMCID: PMC10037493 DOI: 10.1021/acscentsci.2c01373] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/20/2022] [Indexed: 06/18/2023]
Abstract
The diffusion of molecules ("penetrants") of variable size, shape, and chemistry through dense cross-linked polymer networks is a fundamental scientific problem broadly relevant in materials, polymer, physical, and biological chemistry. Relevant applications include separation membranes, barrier materials, drug delivery, and nanofiltration. A major open question is the relationship between transport, thermodynamic state, and penetrant and polymer chemical structure. Here we combine experiment, simulation, and theory to unravel these competing effects on penetrant transport in rubbery and supercooled polymer permanent networks over a wide range of cross-link densities, size ratios, and temperatures. The crucial importance of the coupling of local penetrant hopping to polymer structural relaxation and the secondary importance of mesh confinement effects are established. Network cross-links strongly slow down nm-scale polymer relaxation, which greatly retards the activated penetrant diffusion. The demonstrated good agreement between experiment, simulation, and theory provides strong support for the size ratio (penetrant diameter to the polymer Kuhn length) as a key variable and the usefulness of coarse-grained simulation and theoretical models that average over Angstrom scale structure. The developed theory provides an understanding of the physical processes underlying the behaviors observed in experiment and simulation and suggests new strategies for enhancing selective polymer membrane design.
Collapse
Affiliation(s)
- Baicheng Mei
- Department
of Materials Science, University of Illinois, Urbana, Illinois 61801, United States
- Materials
Research Laboratory, University of Illinois, Urbana, Illinois 61801, United States
| | - Tsai-Wei Lin
- Department
of Chemical & Biomolecular Engineering, University of Illinois, Urbana, Illinois 61801, United States
- Materials
Research Laboratory, University of Illinois, Urbana, Illinois 61801, United States
| | - Grant S. Sheridan
- Department
of Materials Science, University of Illinois, Urbana, Illinois 61801, United States
- Materials
Research Laboratory, University of Illinois, Urbana, Illinois 61801, United States
| | - Christopher M. Evans
- Department
of Materials Science, University of Illinois, Urbana, Illinois 61801, United States
- Department
of Chemical & Biomolecular Engineering, University of Illinois, Urbana, Illinois 61801, United States
- Materials
Research Laboratory, University of Illinois, Urbana, Illinois 61801, United States
| | - Charles E. Sing
- Department
of Materials Science, University of Illinois, Urbana, Illinois 61801, United States
- Department
of Chemical & Biomolecular Engineering, University of Illinois, Urbana, Illinois 61801, United States
- Materials
Research Laboratory, University of Illinois, Urbana, Illinois 61801, United States
| | - Kenneth S. Schweizer
- Department
of Materials Science, University of Illinois, Urbana, Illinois 61801, United States
- Department
of Chemistry, University of Illinois, Urbana, Illinois 61801, United States
- Department
of Chemical & Biomolecular Engineering, University of Illinois, Urbana, Illinois 61801, United States
- Materials
Research Laboratory, University of Illinois, Urbana, Illinois 61801, United States
| |
Collapse
|
5
|
Elucidation of the physical factors that control activated transport of penetrants in chemically complex glass-forming liquids. Proc Natl Acad Sci U S A 2022; 119:e2210094119. [PMID: 36194629 PMCID: PMC9565165 DOI: 10.1073/pnas.2210094119] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Understanding the activated transport of penetrant or tracer atoms and molecules in condensed phases is a challenging problem in chemistry, materials science, physics, and biophysics. Many angstrom- and nanometer-scale features enter due to the highly variable shape, size, interaction, and conformational flexibility of the penetrant and matrix species, leading to a dramatic diversity of penetrant dynamics. Based on a minimalist model of a spherical penetrant in equilibrated dense matrices of hard spheres, a recent microscopic theory that relates hopping transport to local structure has predicted a novel correlation between penetrant diffusivity and the matrix thermodynamic dimensionless compressibility, S0(T) (which also quantifies the amplitude of long wavelength density fluctuations), as a consequence of a fundamental statistical mechanical relationship between structure and thermodynamics. Moreover, the penetrant activation barrier is predicted to have a factorized/multiplicative form, scaling as the product of an inverse power law of S0(T) and a linear/logarithmic function of the penetrant-to-matrix size ratio. This implies an enormous reduction in chemical complexity that is verified based solely on experimental data for diverse classes of chemically complex penetrants dissolved in molecular and polymeric liquids over a wide range of temperatures down to the kinetic glass transition. The predicted corollary that the penetrant diffusion constant decreases exponentially with inverse temperature raised to an exponent determined solely by how S0(T) decreases with cooling is also verified experimentally. Our findings are relevant to fundamental questions in glassy dynamics, self-averaging of angstrom-scale chemical features, and applications such as membrane separations, barrier coatings, drug delivery, and self-healing.
Collapse
|
6
|
Mei B, Schweizer KS. Theory of the Effects of Specific Attractions and Chain Connectivity on the Activated Dynamics and Selective Transport of Penetrants in Polymer Melts. Macromolecules 2022. [DOI: 10.1021/acs.macromol.2c01549] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Baicheng Mei
- Department of Materials Science, University of Illinois, Urbana, Illinois61801, United States
- Materials Research Laboratory, University of Illinois, Urbana, Illinois61801, United States
| | - Kenneth S. Schweizer
- Department of Materials Science, University of Illinois, Urbana, Illinois61801, United States
- Department of Chemistry, University of Illinois, Urbana, Illinois61801, United States
- Department of Chemical and Biomolecular Engineering, University of Illinois, Urbana, Illinois61801, United States
- Materials Research Laboratory, University of Illinois, Urbana, Illinois61801, United States
| |
Collapse
|
7
|
Mei B, Lin TW, Sheridan GS, Evans CM, Sing CE, Schweizer KS. Structural Relaxation and Vitrification in Dense Cross-Linked Polymer Networks: Simulation, Theory, and Experiment. Macromolecules 2022. [DOI: 10.1021/acs.macromol.2c00277] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- Baicheng Mei
- Department of Materials Science, University of Illinois, Urbana, Illinois 61801, United States
- Materials Research Laboratory, University of Illinois, Urbana, Illinois 61801, United States
| | - Tsai-Wei Lin
- Department of Chemical & Biomolecular Engineering, University of Illinois, Urbana, Illinois 61801, United States
- Materials Research Laboratory, University of Illinois, Urbana, Illinois 61801, United States
| | - Grant S. Sheridan
- Department of Materials Science, University of Illinois, Urbana, Illinois 61801, United States
- Materials Research Laboratory, University of Illinois, Urbana, Illinois 61801, United States
| | - Christopher M. Evans
- Department of Materials Science, University of Illinois, Urbana, Illinois 61801, United States
- Department of Chemical & Biomolecular Engineering, University of Illinois, Urbana, Illinois 61801, United States
- Materials Research Laboratory, University of Illinois, Urbana, Illinois 61801, United States
| | - Charles E. Sing
- Department of Materials Science, University of Illinois, Urbana, Illinois 61801, United States
- Department of Chemical & Biomolecular Engineering, University of Illinois, Urbana, Illinois 61801, United States
- Materials Research Laboratory, University of Illinois, Urbana, Illinois 61801, United States
| | - Kenneth S. Schweizer
- Department of Materials Science, University of Illinois, Urbana, Illinois 61801, United States
- Department of Chemistry, University of Illinois, Urbana, Illinois 61801, United States
- Department of Chemical & Biomolecular Engineering, University of Illinois, Urbana, Illinois 61801, United States
- Materials Research Laboratory, University of Illinois, Urbana, Illinois 61801, United States
| |
Collapse
|
8
|
Affiliation(s)
- Catalin P. Gainaru
- Chemical Sciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37830, United States
| | - Alexei P. Sokolov
- Chemical Sciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37830, United States
- Department of Chemistry, University of Tennessee, Knoxville, Tennessee 37996, United States
| |
Collapse
|
9
|
Weigl P, Hecksher T, Dyre JC, Walther T, Blochowicz T. Identity of the local and macroscopic dynamic elastic responses in supercooled 1-propanol. Phys Chem Chem Phys 2021; 23:16537-16541. [PMID: 34312639 DOI: 10.1039/d1cp02671b] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Glass-forming liquids are well known to have significant dynamic heterogeneities, leading to spatially grossly varying elastic properties throughout the system. In this paper, we compare the local elastic response of supercooled 1-propanol monitored by triplet state solvation dynamics to the macroscopic dynamic shear modulus measured by a piezo-electric gauge. The time-dependent responses are found to be identical, which means that the dynamic macroscopic shear modulus provides a good measure of the average local elastic properties. Since the macroscopic shear modulus of a dynamically inhomogeneous system in general is not just the average of the local moduli, there was no reason to expect such a result. This surprising finding not only provides constraints for models of dynamical heterogeneities in glass-forming liquids, but also allows for a fairly straightforward check on elastic models for glassy dynamics.
Collapse
Affiliation(s)
- Peter Weigl
- Institute for Condensed Matter Physics, Technical University of Darmstadt, 64289 Darmstadt, Germany.
| | | | | | | | | |
Collapse
|
10
|
Mei B, Schweizer KS. Theory of the effect of external stress on the activated dynamics and transport of dilute penetrants in supercooled liquids and glasses. J Chem Phys 2021; 155:054505. [PMID: 34364324 DOI: 10.1063/5.0056920] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
We generalize the self-consistent cooperative hopping theory for a dilute spherical penetrant or tracer activated dynamics in dense metastable hard sphere fluids and glasses to address the effect of external stress, the consequences of which are systematically established as a function of matrix packing fraction and penetrant-to-matrix size ratio. All relaxation processes speed up under stress, but the difference between the penetrant and matrix hopping (alpha relaxation) times decreases significantly with stress corresponding to less time scale decoupling. A dynamic crossover occurs at a critical "slaving onset" stress beyond which the matrix activated hopping relaxation time controls the penetrant hopping time. This characteristic stress increases (decreases) exponentially with packing fraction (size ratio) and can be well below the absolute yield stress of the matrix. Below the slaving onset, the penetrant hopping time is predicted to vary exponentially with stress, differing from the power law dependence of the pure matrix alpha time due to system-specificity of the stress-induced changes in the penetrant local cage and elastic barriers. An exponential growth of the penetrant alpha relaxation time with size ratio under stress is predicted, and at a fixed matrix packing fraction, the exponential relation between penetrant hopping time and stress for different size ratios can be collapsed onto a master curve. Direct connections between the short- and long-time activated penetrant dynamics and between the penetrant (or matrix) alpha relaxation time and matrix thermodynamic dimensionless compressibility are also predicted. The presented results should be testable in future experiments and simulations.
Collapse
Affiliation(s)
- Baicheng Mei
- Department of Materials Science, University of Illinois, Urbana, Illinois 61801, USA
| | - Kenneth S Schweizer
- Department of Materials Science, University of Illinois, Urbana, Illinois 61801, USA
| |
Collapse
|
11
|
Jhalaria M, Huang Y, Ruzicka E, Tyagi M, Zorn R, Zamponi M, García Sakai V, Benicewicz B, Kumar S. Activated Transport in Polymer Grafted Nanoparticle Melts. Macromolecules 2021. [DOI: 10.1021/acs.macromol.1c00601] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Mayank Jhalaria
- Department of Chemical Engineering, Columbia University, New York, New York 10027, United States
| | - Yucheng Huang
- Department of Chemistry and Biochemistry, University of South Carolina, Columbia, South Carolina 29208, United States
| | - Eric Ruzicka
- Department of Chemistry and Biochemistry, University of South Carolina, Columbia, South Carolina 29208, United States
| | - Madhusudan Tyagi
- NIST Center for Neutron Research, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-6102, United States
- Department of Materials Science and Engineering, University of Maryland, College Park, Maryland 20742, United States
| | - Reiner Zorn
- Forschungszentrum Jülich GmbH, Jülich Centre for Neutron Science (JCNS-1) and Institute for Biological Information Processing (IBI-8), 52425 Jülich, Germany
| | - Michaela Zamponi
- Forschungszentrum Jülich GmbH, Jülich Centre for Neutron Science at MLZ, Lichtenbergstr. 1 85748 Garching, Germany
| | - Victoria García Sakai
- ISIS Neutron and Muon Source, Rutherford Appleton Laboratory, Chilton, Oxfordshire OX11 0QX, UK
| | - Brian Benicewicz
- Department of Chemistry and Biochemistry, University of South Carolina, Columbia, South Carolina 29208, United States
| | - Sanat Kumar
- Department of Chemical Engineering, Columbia University, New York, New York 10027, United States
| |
Collapse
|
12
|
Mei B, Schweizer KS. Activated penetrant dynamics in glass forming liquids: size effects, decoupling, slaving, collective elasticity and correlation with matrix compressibility. SOFT MATTER 2021; 17:2624-2639. [PMID: 33528485 DOI: 10.1039/d0sm02215b] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
We employ the microscopic self-consistent cooperative hopping theory of penetrant activated dynamics in glass forming viscous liquids and colloidal suspensions to address new questions over a wide range of high matrix packing fractions and penetrant-to-matrix particle size ratios. The focus is on the mean activated relaxation time of smaller tracers in a hard sphere fluid of larger particle matrices. This quantity also determines the penetrant diffusion constant and connects directly with the structural relaxation time probed in an incoherent dynamic structure factor measurement. The timescale of the non-activated fast dissipative process is also studied and is predicted to follow power laws with the contact value of the penetrant-matrix pair correlation function and the penetrant-matrix size ratio. For long time penetrant relaxation, in the relatively lower packing fraction metastable regime the local cage barriers are dominant and matrix collective elasticity effects unimportant. As packing fraction and/or penetrant size grows, much higher barriers emerge and the collective elasticity associated with the correlated matrix dynamic displacement that facilitates penetrant hopping becomes important. This results in a non-monotonic variation with packing fraction of the degree of decoupling between the matrix and penetrant alpha relaxation times. The conditions required for penetrant hopping to become slaved to the matrix alpha process are determined, which depend mainly on the penetrant to matrix particle size ratio. By analyzing the absolute and relative importance of the cage and elastic barriers we establish a mechanistic understanding of the origin of the predicted exponential growth of the penetrant hopping time with size ratio predicted at very high packing fractions. A dynamics-thermodynamics power law connection between the penetrant activation barrier and the matrix dimensionless compressibility is established as a prediction of theory, with different scaling exponents depending on whether matrix collective elasticity effects are important. Quantitative comparisons with simulations of the penetrant relaxation time, diffusion constant, and transient localization length of tracers in dense colloidal suspensions and cold viscous liquids reveal good agreements. Multiple new predictions are made that are testable via future experiments and simulations. Extension of the theoretical approach to more complex systems of high experimental interest (nonspherical molecules, semiflexible polymers, crosslinked networks) interacting via variable hard or soft repulsions and/or short range attractions is possible, including under external deformation.
Collapse
Affiliation(s)
- Baicheng Mei
- Department of Materials Science, University of Illinois, Urbana, IL 61801, USA. and Materials Research Laboratory, University of Illinois, Urbana, IL 61801, USA
| | - Kenneth S Schweizer
- Department of Materials Science, University of Illinois, Urbana, IL 61801, USA. and Materials Research Laboratory, University of Illinois, Urbana, IL 61801, USA and Department of Chemistry, University of Illinois, Urbana, IL 61801, USA and Department of Chemical & Biomolecular Engineering, University of Illinois, Urbana, IL 61801, USA
| |
Collapse
|
13
|
Kanduč M, Kim WK, Roa R, Dzubiella J. How the Shape and Chemistry of Molecular Penetrants Control Responsive Hydrogel Permeability. ACS NANO 2021; 15:614-624. [PMID: 33382598 PMCID: PMC7844830 DOI: 10.1021/acsnano.0c06319] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/28/2020] [Accepted: 12/28/2020] [Indexed: 06/12/2023]
Abstract
The permeability of hydrogels for the selective transport of molecular penetrants (drugs, toxins, reactants, etc.) is a central property in the design of soft functional materials, for instance in biomedical, pharmaceutical, and nanocatalysis applications. However, the permeation of dense and hydrated polymer membranes is a complex multifaceted molecular-level phenomenon, and our understanding of the underlying physicochemical principles is still very limited. Here, we uncover the molecular principles of permeability and selectivity in hydrogel permeation. We combine the solution-diffusion model for permeability with comprehensive atomistic simulations of molecules of various shapes and polarities in a responsive hydrogel in different hydration states. We find in particular that dense collapsed states are extremely selective, owing to a delicate balance between the partitioning and diffusivity of the penetrants. These properties are sensitively tuned by the penetrant size, shape, and chemistry, leading to vast cancellation effects, which nontrivially contribute to the permeability. The gained insights enable us to formulate semiempirical rules to quantify and extrapolate the permeability categorized by classes of molecules. They can be used as approximate guiding ("rule-of-thumb") principles to optimize penetrant or membrane physicochemical properties for a desired permeability and membrane functionality.
Collapse
Affiliation(s)
- Matej Kanduč
- Jožef
Stefan Institute, Jamova 39, SI-1000 Ljubljana, Slovenia
| | - Won Kyu Kim
- Korea
Institute for Advanced Study, 85 Hoegiro, Seoul 02455, Republic of Korea
| | - Rafael Roa
- Departamento
de Física Aplicada I, Facultad de Ciencias, Universidad de Málaga, Campus de Teatinos s/n, E-29071 Málaga, Spain
| | - Joachim Dzubiella
- Applied
Theoretical Physics−Computational Physics, Physikalisches Institut, Albert-Ludwigs-Universität Freiburg, Hermann-Herder Strasse 3, D-79104 Freiburg, Germany
- Research
Group for Simulations of Energy Materials, Helmholtz-Zentrum Berlin für Materialien und Energie, Hahn-Meitner-Platz 1, D-14109 Berlin, Germany
| |
Collapse
|
14
|
Abstract
We combine ideas from polymer and glassy liquid physics to construct a new model for the bond-breaking time scale of attractive sticker groups in associating copolymer liquids that form transient networks. The activated event is argued to be a two-step process, involving first the release of the nonsticker dynamic caging constraints that defines the primary alpha relaxation, followed by attractive stickers surmounting an association free-energy barrier subject to a local frictional resistance which can be strongly affected by relaxation-diffusion decoupling. The ideas embedded in the model produce a consistent and good description of the bond-breaking time scale for diverse polymer chemistries and architectures as a function of temperature and fraction of sticky groups. Chemically sensible values for association free energies are deduced. In strong contrast, the existing phenomenological models are shown to incur qualitative failures.
Collapse
|
15
|
Phan AD, Knapik-Kowalczuk J, Paluch M, Hoang TX, Wakabayashi K. Theoretical Model for the Structural Relaxation Time in Coamorphous Drugs. Mol Pharm 2019; 16:2992-2998. [PMID: 31095393 DOI: 10.1021/acs.molpharmaceut.9b00230] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We propose a simple approach to investigate the structural relaxation time and glass transition of amorphous drugs. Amorphous materials are modeled as a set of equal sized hard spheres. The structural relaxation time over many decades in hard-sphere fluids is theoretically calculated using the elastically collective nonlinear Langevin equation theory associated with Kramer's theory. Then, new thermal mapping from a real material to an effective hard-sphere fluid provides temperature-dependent relaxation time, which can be compared to experiments. Numerical results quantitatively agree with previous experiments for pharmaceutical binary mixtures having different weight ratios. We carry out experiments to test our calculations for an ezetimibe-simvastatin-Kollidon VA64 mixture. Our approach would provide a simple but comprehensive description of glassy dynamics in amorphous composites.
Collapse
Affiliation(s)
| | - Justyna Knapik-Kowalczuk
- Institute of Physics , University of Silesia, SMCEBI , 75 Pułku Piechoty 1a , 41-500 Chorzów , Poland
| | - Marian Paluch
- Institute of Physics , University of Silesia, SMCEBI , 75 Pułku Piechoty 1a , 41-500 Chorzów , Poland
| | - Trinh X Hoang
- Institute of Physics , Vietnam Academy of Science and Technology , 10 Dao Tan, Ba Dinh , Hanoi 100000 , Vietnam
| | - Katsunori Wakabayashi
- Department of Nanotechnology for Sustainable Energy, School of Science and Technology , Kwansei Gakuin University , Sanda , Hyogo 669-1337 , Japan
| |
Collapse
|
16
|
Poling-Skutvik R, Roberts RC, Slim AH, Narayanan S, Krishnamoorti R, Palmer JC, Conrad JC. Structure Dominates Localization of Tracers within Aging Nanoparticle Glasses. J Phys Chem Lett 2019; 10:1784-1789. [PMID: 30916569 DOI: 10.1021/acs.jpclett.9b00309] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
We investigate the transport and localization of tracer probes in a glassy matrix as a function of relative size using dynamic X-ray scattering experiments and molecular dynamics simulations. The quiescent relaxations of tracer particles evolve with increasing waiting time, tw. The corresponding relaxation times increase exponentially at small tw and then transition to a power-law behavior at longer tw. As tracer size decreases, the aging behavior weakens and the particles become less localized within the matrix until they delocalize at a critical size ratio δ0 ≈ 0.38. Localization does not vary with sample age even as the relaxations slow by approximately an order of magnitude, suggesting that matrix structure controls tracer localization.
Collapse
Affiliation(s)
- Ryan Poling-Skutvik
- Department of Chemical and Biomolecular Engineering , University of Houston , Houston , Texas 77204-4004 , United States
| | - Ryan C Roberts
- Department of Chemical and Biomolecular Engineering , University of Houston , Houston , Texas 77204-4004 , United States
| | - Ali H Slim
- Department of Chemical and Biomolecular Engineering , University of Houston , Houston , Texas 77204-4004 , United States
| | - Suresh Narayanan
- Advanced Photon Source , Argonne National Laboratory , Argonne , Illinois 60439 , United States
| | - Ramanan Krishnamoorti
- Department of Chemical and Biomolecular Engineering , University of Houston , Houston , Texas 77204-4004 , United States
| | - Jeremy C Palmer
- Department of Chemical and Biomolecular Engineering , University of Houston , Houston , Texas 77204-4004 , United States
| | - Jacinta C Conrad
- Department of Chemical and Biomolecular Engineering , University of Houston , Houston , Texas 77204-4004 , United States
| |
Collapse
|
17
|
Kim WK, Kanduč M, Roa R, Dzubiella J. Tuning the Permeability of Dense Membranes by Shaping Nanoscale Potentials. PHYSICAL REVIEW LETTERS 2019; 122:108001. [PMID: 30932643 DOI: 10.1103/physrevlett.122.108001] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/21/2018] [Revised: 12/17/2018] [Indexed: 06/09/2023]
Abstract
Permeability is one of the most fundamental transport properties in soft matter physics, material engineering, and nanofluidics. Here, we report by means of Langevin simulations of ideal penetrants in a nanoscale membrane made of a fixed lattice of attractive interaction sites, how the permeability can be massively tuned, even minimized or maximized, by tailoring the potential energy landscape for the diffusing penetrants, depending on the membrane attraction, topology, and density. Supported by limiting scaling theories we demonstrate that the observed nonmonotonic behavior and the occurrence of extreme values of the permeability is far from trivial and triggered by a strong anticorrelation and substantial (orders of magnitude) cancellation between penetrant partitioning and diffusivity, especially within dense and highly attractive membranes.
Collapse
Affiliation(s)
- Won Kyu Kim
- Research Group for Simulations of Energy Materials, Helmholtz-Zentrum Berlin, D-14109 Berlin, Germany
| | - Matej Kanduč
- Research Group for Simulations of Energy Materials, Helmholtz-Zentrum Berlin, D-14109 Berlin, Germany
- Jožef Stefan Institute, SI-1000 Ljubljana, Slovenia
| | - Rafael Roa
- Física Aplicada I, Facultad de Ciencias, Universidad de Málaga, 29071 Málaga, Spain
| | - Joachim Dzubiella
- Research Group for Simulations of Energy Materials, Helmholtz-Zentrum Berlin, D-14109 Berlin, Germany
- Applied Theoretical Physics-Computational Physics, Physikalisches Institut, Albert-Ludwigs-Universität Freiburg, D-79104 Freiburg, Germany
| |
Collapse
|
18
|
Phan AD, Schweizer KS. Elastically Collective Nonlinear Langevin Equation Theory of Glass-Forming Liquids: Transient Localization, Thermodynamic Mapping, and Cooperativity. J Phys Chem B 2018; 122:8451-8461. [DOI: 10.1021/acs.jpcb.8b04975] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
19
|
Kanduč M, Kim WK, Roa R, Dzubiella J. Selective Molecular Transport in Thermoresponsive Polymer Membranes: Role of Nanoscale Hydration and Fluctuations. Macromolecules 2018. [DOI: 10.1021/acs.macromol.8b00735] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Matej Kanduč
- Research Group for Simulations of Energy Materials, Helmholtz-Zentrum Berlin für Materialien und Energie, Hahn-Meitner-Platz 1, D-14109 Berlin, Germany
| | - Won Kyu Kim
- Research Group for Simulations of Energy Materials, Helmholtz-Zentrum Berlin für Materialien und Energie, Hahn-Meitner-Platz 1, D-14109 Berlin, Germany
| | - Rafael Roa
- Research Group for Simulations of Energy Materials, Helmholtz-Zentrum Berlin für Materialien und Energie, Hahn-Meitner-Platz 1, D-14109 Berlin, Germany
- Departamento de Física Aplicada I, Facultad de Ciencias, Universidad de Málaga, Campus de Teatinos s/n, E-29071 Málaga, Spain
| | - Joachim Dzubiella
- Research Group for Simulations of Energy Materials, Helmholtz-Zentrum Berlin für Materialien und Energie, Hahn-Meitner-Platz 1, D-14109 Berlin, Germany
- Physikalisches Institut, Albert-Ludwigs-Universität Freiburg, Hermann-Herder Strasse 3, D-79104 Freiburg, Germany
| |
Collapse
|
20
|
Roberts RC, Poling-Skutvik R, Palmer JC, Conrad JC. Tracer Transport Probes Relaxation and Structure of Attractive and Repulsive Glassy Liquids. J Phys Chem Lett 2018; 9:3008-3013. [PMID: 29763547 DOI: 10.1021/acs.jpclett.8b01074] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Dynamic coupling of small penetrants to slow, cooperative relaxations within crowded cells, supercooled liquids, and polymer matrices has broad consequences for applications ranging from drug delivery to nanocomposite processing. Interactions between the constituents of these and other disordered media alter the cooperative relaxations, but their effect on penetrant dynamics remains incompletely understood. We use molecular dynamics simulations to show that the motions of hard-sphere tracer particles probe differences in local structure and cooperative relaxation processes in attractive and repulsive glassy liquid matrices with equal bulk packing fractions and long-time diffusivities. Coupling of the tracer dynamics to collective relaxations in each matrix affects the shape of tracer trajectories, which are fractal within the repulsive matrix and more compact in the attractive. These results reveal that the structure of relaxations controls penetrant transport and dispersion in cooperatively relaxing systems and provide insight into dynamical heterogeneity within glassy liquids.
Collapse
Affiliation(s)
- Ryan C Roberts
- Department of Chemical and Biomolecular Engineering , University of Houston , Houston , Texas 77204-4004 , United States
| | - Ryan Poling-Skutvik
- Department of Chemical and Biomolecular Engineering , University of Houston , Houston , Texas 77204-4004 , United States
| | - Jeremy C Palmer
- Department of Chemical and Biomolecular Engineering , University of Houston , Houston , Texas 77204-4004 , United States
| | - Jacinta C Conrad
- Department of Chemical and Biomolecular Engineering , University of Houston , Houston , Texas 77204-4004 , United States
| |
Collapse
|
21
|
Affiliation(s)
- Yuxing Zhou
- Department of Chemical Engineering, Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Scott T. Milner
- Department of Chemical Engineering, Pennsylvania State University, University Park, Pennsylvania 16802, United States
| |
Collapse
|
22
|
Zhang R, Schweizer KS. Microscopic Theory of Coupled Slow Activated Dynamics in Glass-Forming Binary Mixtures. J Phys Chem B 2018; 122:3465-3479. [DOI: 10.1021/acs.jpcb.7b10568] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
23
|
Zhang K, Meng D, Müller-Plathe F, Kumar SK. Coarse-grained molecular dynamics simulation of activated penetrant transport in glassy polymers. SOFT MATTER 2018; 14:440-447. [PMID: 29261207 DOI: 10.1039/c7sm01941f] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Membrane separations of gas mixtures strive to maximize the permeability of a desired species while keeping out undesired ones. Permeability vs. selectivity data from many polymer membranes for a given gas pair with diameters dA and dB are typically collected in a "Robeson plot"', and are bound from above by a line with a slope λ = (dB/dA)2 - 1. A microscopic understanding of this relationship, especially λ, is still missing. We perform molecular dynamics simulations of penetrant diffusion using three different coarse-grained polymer models over a wide range of penetrant sizes, temperatures, and monomer densities. The empirically relevant λ = (dB/dA)2 - 1 is only found for polymers that are either supercooled liquids with caged segmental dynamics or glasses and when the penetrant size is approximately half the Kuhn length of the chains, for which the penetrant diffusion is an activated process.
Collapse
Affiliation(s)
- Kai Zhang
- Department of Chemical Engineering, Columbia University, New York, New York 10027, USA.
| | | | | | | |
Collapse
|
24
|
Affiliation(s)
- Kai Zhang
- Department of Chemical Engineering, Columbia University, New York, New York 10027, United States
| | - Sanat K. Kumar
- Department of Chemical Engineering, Columbia University, New York, New York 10027, United States
| |
Collapse
|