1
|
Mei B, Sheridan GS, Evans CM, Schweizer KS. 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] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2022] [Accepted: 09/13/2022] [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
Affiliation(s)
- Baicheng Mei
- Department of Materials Science, University of Illinois at Urbana-Champaign, Urbana, IL 61801
- Material Research Laboratory, University of Illinois at Urbana-Champaign, Urbana, IL 61801
| | - Grant S. Sheridan
- Department of Materials Science, University of Illinois at Urbana-Champaign, Urbana, IL 61801
- Material Research Laboratory, University of Illinois at Urbana-Champaign, Urbana, IL 61801
| | - Christopher M. Evans
- Department of Materials Science, University of Illinois at Urbana-Champaign, Urbana, IL 61801
- Material Research Laboratory, University of Illinois at Urbana-Champaign, Urbana, IL 61801
- Department of Chemical & Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, IL 61801
| | - Kenneth S. Schweizer
- Department of Materials Science, University of Illinois at Urbana-Champaign, Urbana, IL 61801
- Material Research Laboratory, University of Illinois at Urbana-Champaign, Urbana, IL 61801
- Department of Chemical & Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, IL 61801
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, IL 61801
| |
Collapse
|
2
|
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
|
3
|
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
|
4
|
Stubbing JW, McCoustra MRS, Brown WA. A new technique for determining the refractive index of ices at cryogenic temperatures. Phys Chem Chem Phys 2020; 22:25353-25365. [PMID: 33140768 DOI: 10.1039/d0cp02373f] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A reflection-absorption optical (RAO) spectrometer, operating across the ultra-violet/visible (UV/visible) wavelength region, has been developed that allows simultaneous measurements of optical properties and thickness of thin solid films at cryogenic temperatures in ultrahigh vacuum. The RAO spectrometer enables such measurements to be made after ice deposition, as opposed to most current approaches which make measurements during deposition. This allows changes in the optical properties and in the thickness of the film to be determined subsequent to thermal, photon or charged particle processing. This is not possible with current techniques. A data analysis method is presented that allows the wavelength dependent n and k values for ices to be extracted from the reflection-absorption spectra. The validity of this analysis method is shown using model data from the literature. New data are presented for the reflection UV/visible spectra of amorphous and crystalline single component ices of benzene, methyl formate and water adsorbed on a graphite surface. These data show that, for benzene and methyl formate, the crystalline ice has a larger refractive index than amorphous ice, reflecting changes in the electronic environment occurring in the ice during crystallisation. For water, the refractive index does not vary with ice phase.
Collapse
Affiliation(s)
- James W Stubbing
- Department of Chemistry, School of Life Sciences, University of Sussex, Falmer, Brighton, BN1 9QJ, UK.
| | | | | |
Collapse
|
5
|
Du K, Zemerov SD, Hurtado Parra S, Kikkawa JM, Dmochowski IJ. Paramagnetic Organocobalt Capsule Revealing Xenon Host-Guest Chemistry. Inorg Chem 2020; 59:13831-13844. [PMID: 32207611 PMCID: PMC7672707 DOI: 10.1021/acs.inorgchem.9b03634] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
We investigated Xe binding in a previously reported paramagnetic metal-organic tetrahedral capsule, [Co4L6]4-, where L2- = 4,4'-bis[(2-pyridinylmethylene)amino][1,1'-biphenyl]-2,2'-disulfonate. The Xe-inclusion complex, [XeCo4L6]4-, was confirmed by 1H NMR spectroscopy to be the dominant species in aqueous solution saturated with Xe gas. The measured Xe dissociation rate in [XeCo4L6]4-, koff = 4.45(5) × 102 s-1, was at least 40 times greater than that in the analogous [XeFe4L6]4- complex, highlighting the capability of metal-ligand interactions to tune the capsule size and guest permeability. The rapid exchange of 129Xe nuclei in [XeCo4L6]4- produced significant hyperpolarized 129Xe chemical exchange saturation transfer (hyper-CEST) NMR signal at 298 K, detected at a concentration of [XeCo4L6]4- as low as 100 pM, with presaturation at -89 ppm, which was referenced to solvated 129Xe in H2O. The saturation offset was highly temperature-dependent with a slope of -0.41(3) ppm/K, which is attributed to hyperfine interactions between the encapsulated 129Xe nucleus and electron spins on the four CoII centers. As such, [XeCo4L6]4- represents the first example of a paramagnetic hyper-CEST (paraHYPERCEST) sensor. Remarkably, the hyper-CEST 129Xe NMR resonance for [XeCo4L6]4- (δ = -89 ppm) was shifted 105 ppm upfield from the diamagnetic analogue [XeFe4L6]4- (δ = +16 ppm). The Xe inclusion complex was further characterized in the crystal structure of (C(NH2)3)4[Xe0.7Co4L6]·75 H2O (1). Hydrogen bonding between capsule-linker sulfonate groups and exogenous guanidinium cations, (C(NH2)3)+, stabilized capsule-capsule interactions in the solid state and also assisted in trapping a Xe atom (∼42 Å3) in the large (135 Å3) cavity of 1. Magnetic susceptibility measurements confirmed the presence of four noninteracting, magnetically anisotropic high-spin CoII centers in 1. Furthermore, [Co4L6]4- was found to be stable toward aggregation and oxidation, and the CEST performance of [XeCo4L6]4- was unaffected by biological macromolecules in H2O. These results recommend metal-organic capsules for fundamental investigations of Xe host-guest chemistry as well as applications with highly sensitive 129Xe-based sensors.
Collapse
|
6
|
Cubeta US, Sadtchenko V. Glass softening kinetics in the limit of high heating rates. J Chem Phys 2019; 150:094508. [DOI: 10.1063/1.5046304] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ulyana S. Cubeta
- Chemistry Department, The George Washington University, Washington, District of Columbia 20052, USA
| | - Vlad Sadtchenko
- Chemistry Department, The George Washington University, Washington, District of Columbia 20052, USA
| |
Collapse
|
7
|
Salter TL, Stubbing JW, Brigham L, Brown WA. A TPD and RAIRS comparison of the low temperature surface behavior of benzene, toluene, and xylene on graphite. J Chem Phys 2018; 149:164705. [DOI: 10.1063/1.5051134] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Tara L. Salter
- Department of Chemistry, School of Life Sciences, University of Sussex, Falmer, Brighton BN1 9QJ, United Kingdom
| | - James W. Stubbing
- Department of Chemistry, School of Life Sciences, University of Sussex, Falmer, Brighton BN1 9QJ, United Kingdom
| | - Lorna Brigham
- Department of Chemistry, School of Life Sciences, University of Sussex, Falmer, Brighton BN1 9QJ, United Kingdom
| | - Wendy A. Brown
- Department of Chemistry, School of Life Sciences, University of Sussex, Falmer, Brighton BN1 9QJ, United Kingdom
| |
Collapse
|
8
|
Smith RS, Kay BD. Desorption Kinetics of Benzene and Cyclohexane from a Graphene Surface. J Phys Chem B 2017; 122:587-594. [DOI: 10.1021/acs.jpcb.7b05102] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- R. Scott Smith
- Physical and Computational
Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Bruce D. Kay
- Physical and Computational
Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| |
Collapse
|
9
|
Tu W, Chen Z, Li X, Gao Y, Liu R, Wang LM. Revisiting the glass transition and dynamics of supercooled benzene by calorimetric studies. J Chem Phys 2015; 143:164501. [PMID: 26520521 DOI: 10.1063/1.4933425] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The glass transition and dynamics of benzene are studied in binary mixtures of benzene with five glass forming liquids, which can be divided into three groups: (a) o-terphenyl and m-xylene, (b) N-butyl methacrylate, and (c) N,N-dimethylpropionamide and N,N-diethylformamide to represent the weak, moderate, and strong interactions with benzene. The enthalpies of mixing, ΔH(mix), for the benzene mixtures are measured to show positive or negative signs, with which the validity of the extrapolations of the glass transition temperature T(g) to the benzene-rich regions is examined. The extrapolations for the T(g) data in the mixtures are found to converge around the point of 142 K, producing T(g) of pure benzene. The fragility m of benzene is also evaluated by extrapolating the results of the mixtures, and a fragility m ∼ 80 is yielded. The obtained T(g) and m values for benzene allow for the construction of the activation plot in the deeply supercooled region. The poor glass formability of benzene is found to result from the high melting point, which in turn leads to low viscosity in the supercooled liquid.
Collapse
Affiliation(s)
- Wenkang Tu
- State Key Lab of Metastable Materials Science and Technology, and College of Materials Science and Engineering, Yanshan University, Qinhuangdao, Hebei 066004, China
| | - Zeming Chen
- State Key Lab of Metastable Materials Science and Technology, and College of Materials Science and Engineering, Yanshan University, Qinhuangdao, Hebei 066004, China
| | - Xiangqian Li
- State Key Lab of Metastable Materials Science and Technology, and College of Materials Science and Engineering, Yanshan University, Qinhuangdao, Hebei 066004, China
| | - Yanqin Gao
- State Key Lab of Metastable Materials Science and Technology, and College of Materials Science and Engineering, Yanshan University, Qinhuangdao, Hebei 066004, China
| | - Riping Liu
- State Key Lab of Metastable Materials Science and Technology, and College of Materials Science and Engineering, Yanshan University, Qinhuangdao, Hebei 066004, China
| | - Li-Min Wang
- State Key Lab of Metastable Materials Science and Technology, and College of Materials Science and Engineering, Yanshan University, Qinhuangdao, Hebei 066004, China
| |
Collapse
|
10
|
Zhang R, Schweizer KS. Theory of activated penetrant diffusion in viscous fluids and colloidal suspensions. J Chem Phys 2015; 143:144906. [PMID: 26472397 DOI: 10.1063/1.4932679] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We heuristically formulate a microscopic, force level, self-consistent nonlinear Langevin equation theory for activated barrier hopping and non-hydrodynamic diffusion of a hard sphere penetrant in very dense hard sphere fluid matrices. Penetrant dynamics is controlled by a rich competition between force relaxation due to penetrant self-motion and collective matrix structural (alpha) relaxation. In the absence of penetrant-matrix attraction, three activated dynamical regimes are predicted as a function of penetrant-matrix size ratio which are physically distinguished by penetrant jump distance and the nature of matrix motion required to facilitate its hopping. The penetrant diffusion constant decreases the fastest with size ratio for relatively small penetrants where the matrix effectively acts as a vibrating amorphous solid. Increasing penetrant-matrix attraction strength reduces penetrant diffusivity due to physical bonding. For size ratios approaching unity, a distinct dynamical regime emerges associated with strong slaving of penetrant hopping to matrix structural relaxation. A crossover regime at intermediate penetrant-matrix size ratio connects the two limiting behaviors for hard penetrants, but essentially disappears if there are strong attractions with the matrix. Activated penetrant diffusivity decreases strongly with matrix volume fraction in a manner that intensifies as the size ratio increases. We propose and implement a quasi-universal approach for activated diffusion of a rigid atomic/molecular penetrant in a supercooled liquid based on a mapping between the hard sphere system and thermal liquids. Calculations for specific systems agree reasonably well with experiments over a wide range of temperature, covering more than 10 orders of magnitude of variation of the penetrant diffusion constant.
Collapse
Affiliation(s)
- Rui Zhang
- Department of Materials Science and Frederick Seitz Materials Research Laboratory, University of Illinois, 1304 West Green Street, Urbana, Illinois 61801, USA
| | - Kenneth S Schweizer
- Department of Materials Science and Frederick Seitz Materials Research Laboratory, University of Illinois, 1304 West Green Street, Urbana, Illinois 61801, USA
| |
Collapse
|
11
|
Smith RS, May RA, Kay BD. Probing Toluene and Ethylbenzene Stable Glass Formation Using Inert Gas Permeation. J Phys Chem Lett 2015; 6:3639-3644. [PMID: 26722735 DOI: 10.1021/acs.jpclett.5b01611] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Inert gas permeation is used to investigate the formation of stable glasses of toluene and ethylbenzene. The effect of deposition temperature (T(dep)) on the kinetic stability of the vapor deposited glasses is determined using Kr desorption spectra from within sandwich layers of either toluene or ethylbenzene. The results for toluene show that the most stable glass is formed at T(dep) = 0.92 T(g), although glasses with a kinetic stability within 50% of the most stable glass were found with deposition temperatures from 0.85 to 0.95 T(g). Similar results were found for ethylbenzene, which formed its most stable glass at 0.91 T(g) and formed stable glasses from 0.81 to 0.96 T(g). These results are consistent with recent calorimetric studies and demonstrate that the inert gas permeation technique provides a direct method to observe the onset of molecular translation motion that accompanies the glass to supercooled liquid transition.
Collapse
Affiliation(s)
- R Scott Smith
- Fundamental and Computational Sciences Directorate, Pacific Northwest National Laboratory , Richland, Washington 99352, United States
| | - R Alan May
- Fundamental and Computational Sciences Directorate, Pacific Northwest National Laboratory , Richland, Washington 99352, United States
| | - Bruce D Kay
- Fundamental and Computational Sciences Directorate, Pacific Northwest National Laboratory , Richland, Washington 99352, United States
| |
Collapse
|