Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mass diffusivity”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Enhanced accuracy through ensembling of randomly initialized auto-regressive models for dynamical systems

Computational mechanics simulations using traditional finite element methods (FEM) require prohibitively expensive computational resources for real-time engineering applications, design optimization, and digital twin implementations. While machine learning (ML) surrogate models offer significant computational speedups, autoregressive ML models for time-dependent mechanical systems suffer from error accumulation that compromises long-term prediction reliability - a critical concern for engineering applications where accuracy over extended time horizons is essential for safety and performance assessments. Here, we propose a deep ensemble framework specifically designed to address this challenge in computational mechanics applications, where multiple ML surrogate models with random weight initializations are trained in parallel and their predictions aggregated during inference. This approach leverages statistical diversity to maximize information gain from a fixed set of training data and to mitigate error propagation, while maintaining the computational efficiency that makes ML surrogates attractive for engineering practice. We validate the framework on three representative problems spanning critical areas of computational mechanics: stress field evolution in heterogeneous microstructures under complex loading (relevant to advanced materials design and composite analysis), planetary-scale shallow water dynamics (applicable to environmental and geotechnical engineering), and Gray-Scott reaction-diffusion systems (relevant to mass transport and chemical process engineering). Across all test cases, the ensemble approach demonstrates consistent error reduction of 15-33% compared to individual models. The codes for this work are available on GitHub (https://github.com/Graham-Brady-Research-Group/AutoregressiveEnsemble_SpatioTemporal_Evolution).

autoregressive prediction↗
DOE OSTI2021

Exfoliation of surfactant swollen layered MWW zeolites into two-dimensional zeolite nanosheets using telechelic liquid polybutadiene

We report two-dimensional (2D) zeolites have the potential to enhance mass transport by reducing diffusion lengths and thus find applications in catalysis and separation. High aspect ratio 2D zeolite nanosheets are required for fabrication of thin, high-flux zeolite membranes. Preparation of 2D zeolite nanosheets can be achieved by exfoliation of layered zeolite precursors, which usually requires multiple steps including swelling with cationic surfactants, exfoliation by applying a mechanical force such as shearing or sonication, purification and removal of the surfactants. The exfoliation of surfactant swollen layered zeolite precursors has been studied using molten polystyrene and telechelic liquid polybutadiene, but still requires extensive shearing force. Recently, a novel one-step exfoliation method has been developed for nonswollen MWW zeolite precursors using tetrabutylammonium hydroxide (TBAOH) solution. In this work, we study the exfoliation of surfactant swollen MWW layered precursors (with and without aluminum) in commercially available liquid hydroxyl-terminated polybutadiene (HTPB) with the aim to provide fundamental understanding of the interaction between the liquid polybutadiene and the surfactant swollen layered zeolite precursors. It was found that surfactant swollen ITQ-1, (ITQ-1(S)), a pure silica layered precursor with MWW framework after being swollen with a cationic surfactant, can be exfoliated in liquid HTPB at room temperature without applying any additional shearing force. The aluminum containing surfactant swollen MWW precursors can be also exfoliated in HTPB at room temperature, but require additional shearing forces. The exfoliation process was monitored using rheology experiments and studied by tuning the composition of the zeolite precursors. The interaction between the layered zeolite precursor and the hydroxyl group of the liquid polybutadiene allows the polymer to intercalate the precursor and exfoliate it into nanosheets. The presence of framework aluminum and the nature of the organic structure directing agents (OSDAs) are key parameters affecting the exfoliation process. Framework Al is unfavorable for the exfoliation possibly due to the surface charges introduced by the aluminum. Bulky OSDA, trimethyladamantammonium (TMAda + ), in the interlayer spacing may facilitate the exfoliation.

2D nanosheets↗
DOE OSTI2020

Liquid–Solid Boundaries Dominate Activity of CO 2 Reduction on Gas-Diffusion Electrodes

Electrochemical CO 2 electrolysis to produce hydrocarbon fuels or material feedstocks offers a renewable alternative to fossilized carbon sources. Gas-diffusion electrodes (GDEs), composed of solid electrocatalysts on porous supports positioned near the interface of a conducting electrolyte and CO 2 gas, have been able to demonstrate the substantial current densities needed for future commercialization. These higher reaction rates have often been ascribed to the presence of a three-phase interface, where solid, liquid, and gas provide electrons, water, and CO 2 , respectively. Conversely, mechanistic work on electrochemical reactions implicates a fully two-phase reaction interface, where gas molecules reach the electrocatalyst’s surface by dissolution and diffusion through the electrolyte. Because the discrepancy between an atomistic three-phase versus two-phase reaction has substantial implications for the design of catalysts, gas-diffusion layers, and cell architectures, the nuances of nomenclatures and governing phenomena surrounding the three-phase-region require clarification. In this paper we outline the macro, micro, and atomistic phenomena occurring within a gas-diffusion electrode to provide a focused discussion on the architecture of the often-discussed three-phase region for CO 2 electrolysis. From this information, we comment on the outlook for the broader CO 2 electroreduction GDE cell architecture.

CO2 reduction↗
DOE OSTI2021

Small to Full-Height Scale Comparisons of Cesium Ion Exchange Performance with Crystalline Silicotitanate

Crystalline silicotitanate (CST) ion exchanger is planned to be used to remove cesium (137Cs) from the aqueous phases of Hanford tank wastes in preparation for vitrification at the Waste Treatment and Immobilization Plant (WTP) Pretreatment Facility. Column scale up testing was conducted to evaluate performance of Cs exchange onto the CST. Batch contact testing was conducted to assess exchange kinetics of Cs exchange at four different CST PSDs. Column testing was conducted at three column sizes, small (2.5% full height), medium (12% full height), and full height, to assess Cs load performance behavior. Testing at the small scale was compared to actual Hanford tank waste testing to verify the validity of simulant tests to accurately represent full height column performance. A change in CST particle size was essential in scaling the small column dynamics up to full scale. Based on these results, a determination of intraparticle and film diffusion impacts on overall mass transfer coefficients will allow future modeling of the breakthrough performance at a range of process conditions.

Westesen, Amy M.↗
DOE OSTI2020

Preliminary Reporting of Thermophysical Property Measurements for the Ghareb Formation

Accurate knowledge of thermophysical properties of rock is vital to develop meaningful models of high level nuclear waste emplacement scenarios. The Israel Atomic Energy Commission is considering storing high level nuclear waste in the Ghareb formation, a porous kerogen bearing chalk. Sandia is supporting this effort with an evolving lab- based geomechanics testing program. We have completed measurements of thermal properties up to 275C and room temperature hydrostatic compaction measurements. We report thermal conductivity, thermal diffusivity, specific heat, and mass loss from our thermal measurements, and we report bulk moduli and porosity loss from our compaction measurements. These values are crucial for the numerical models to simulate heat transfer and formation compressibility around a heat generating repository.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Improved melt model for power flow

Accelerators that drive z-pinch experiments transport current densities in excess of 1 MA/cm 2 in order to melt or ionize the target and implode it on axis. These high current densities stress the transmission lines upstream from the target, where rapid electrode heating causes plasma formation, melt, and possibly vaporization. These plasmas negatively impact accelerator efficiency by diverting some portion of the current away from the target, referred to as “current loss”. Simulations that are able to reproduce this behavior may be applied to improving the efficiency of existing accelerators and to designing systems operating at ever higher current densities. The relativistic particle-in-cell code CHICAGO ® is the primary code for modeling power flow on Sandia National Laboratories’ Z accelerator. We report here on new algorithms that incorporate vaporization and melt into the standard power-flow simulation framework. Taking a hybrid approach, the CHICAGO® kinetic/multi-fluid treatment has been expanded to include vaporization while the quasi-neutral equation-of-motion has been updated for melt at high current-densities. For vaporization, a new one-dimensional substrate model provides a more accurate calculation of electrode thermal, mass, and magnetic field diffusion as well as a means of emitting absorbed contaminants and vaporized metal ions. A quasi-fluid model has been implemented expressly to mimic the motion of imploding liners for accurate inductance histories. For melt, a multi-ion Hall-MHD option has been implemented and benchmarked against Alegra MHD. This new model is described with sufficient detail to reproduce these algorithms in any hybrid kinetic code. Physics results from the new code are also presented. A CHICAGO ® Hall-MHD simulation of a radial transmission line demonstrates that Hall physics, not included in Alegra, has no significant impact on the diffusion of electrode material. When surface contaminant desorption is mocked in as a hydrogen surface plasma, both the surface and bulk-material plasmas largely compress under the influence of the j × B force. Similar results are seen in Alegra, which also shows magnetic and material diffusion scaling with peak current. Test vaporization simulations using MagLIF and a power-flow experimental geometry show Fe + ions diffuse only a few hundred µm from the electrodes, so present models of Z power flow remain valid.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Comparing the MELCOR Aerosol Deposition Model with Exact Analytical Solutions

Exact analytical solutions are presented for the evolution of the aerosol particle mass density function in a control volume for particle deposition due to gravitational settling, thermophoresis, and diffusion. The solutions are for arbitrary initial mass density functions and are applied for an initial lognormal density function. Integration of these solutions provides the suspended mass in the control volume as a function of time. These solutions serve as an exact benchmark to assess the accuracy of numerical methods. For the numerical algorithm used in MELCOR, excellent agreement is obtained for gravitational settling, diffusive deposition, and thermophoretic deposition for the suspended aerosol mass. In all cases, the default number of discrete particle size bins of 10 is shown to converge, with hardly any advantage to using 20 size bins.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Impact of binary interactions on the diffuse supernova neutrino background

Binary interactions, especially mass transfer and mergers, can strongly influence the evolution of massive stars and change their final properties and the occurrence of supernovae. Here, we investigate how binary interactions affect predictions of the diffuse flux of neutrinos. By performing stellar population syntheses including prescriptions for binary interactions, we show that the resulting detection rates of the diffuse supernova neutrino background is enhanced by 15%–20% compared to estimates without binary considerations. A source of significant uncertainty arises due to the presently sparse knowledge of the evolution of rapidly rotating carbon-oxygen cores, especially those created as a result of mergers near the white dwarf to core collapse boundary. As a result, the enhancement effect may be as small as a few percent if the effects of rotation in postmerger systems are neglected, or as large as 75% if trends are extrapolated. Our estimates serve to highlight that binary effects can be important.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Atomistic processes of surface-diffusion-induced abnormal softening in nanoscale metallic crystals

Abstract Ultrahigh surface-to-volume ratio in nanoscale materials, could dramatically facilitate mass transport, leading to surface-mediated diffusion similar to Coble-type creep in polycrystalline materials. Unfortunately, the Coble creep is just a conceptual model, and the associated physical mechanisms of mass transport have never been revealed at atomic scale. Akin to the ambiguities in Coble creep, atomic surface diffusion in nanoscale crystals remains largely unclear, especially when mediating yielding and plastic flow. Here, by using in situ nanomechanical testing under high-resolution transmission electron microscope, we find that the diffusion-assisted dislocation nucleation induces the transition from a normal to an inverse Hall-Petch-like relation of the strength-size dependence and the surface-creep leads to the abnormal softening in flow stress with the reduction in size of nanoscale silver, contrary to the classical “alternating dislocation starvation” behavior in nanoscale platinum. This work provides insights into the atomic-scale mechanisms of diffusion-mediated deformation in nanoscale materials, and impact on the design for ultrasmall-sized nanomechanical devices.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Single-Molecule Tracking Measurements Reveal the Detailed Mechanisms of Molecular Diffusion in Solvent Mixtures under Nanoconfinement

Understanding mass transport mechanisms in nanopores is important for developing advanced materials for chemical separations, chemical sensing, and energy storage. This paper reports a novel imaging platform that is employed for the first time to investigate the detailed diffusion dynamics of single rhodamine B (RhB) dye molecules confined within solution-filled cylindrical anodic aluminum oxide (AAO) nanopores. The imaging platform relies on illumination of horizontally-oriented AAO nanopores in a highly inclined and laminated optical (HILO) light sheet microscopy geometry. The method was used to investigate the translational and orientational dynamics of single rhodamine B (RhB) molecules within horizontally-oriented 5- and 10-nm diameter AAO nanopores filled with water–ethanol mixtures. The established platform enabled the observation of one-dimensional motion along the pore axis involving occasional short- or long-term immobilization at the single-molecule level. Analysis of cumulative squared-displacement distributions revealed fast (5 – 30 µm²/s), intermediate (1 – 5 µm²/s), and slow (< 1 µm²/s) diffusion components. From the effects of mixture composition and pore size on the contributions of these three components, we inferred that the fast, intermediate, and slow components could be assigned to desorption-mediated hopping, crawling, and wiggling motions, respectively. The platform based on the horizontally-oriented AAO nanopores also permitted single-molecule emission polarization measurements that revealed the negligible steric confinement of individual diffusing RhB molecules within the AAO nanopores. The imaging platform based on AAO membranes and HILO microscopy provided a unique means to investigate how solvation-mediated surface interactions and nanoconfinement govern molecular transport in nanoporous environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Tritium Transport Phenomena in Molten-Salt Reactors: Molten Salt Tritium Transport Experiment Design

Tritium is produced from neutron interactions with both lithium and beryllium. Large quantities of tritium are generated in Molten Salt Reactors (MSRs) which use LiF/BeF2 (FLiBe) as the fuel salt. Tritium is unique among the radionuclide hazards as it readily permeates through metal structural materials at high temperatures. All metal surfaces are potential release paths for tritium. For adequate safety analysis and eventual licensing of new reactors, predictive models for tritium transport and release from MSRs must be developed. These models must account for the multiple transport phenomena involved with tritium: fuel salt phase mass transport, dissociation/recombination reactions on metal surfaces, interstitial diffusion through the metal structure, and salt or gas phase mass transport in the downstream fluid. These models also must also be validated with representative experiments. Our previous report outlined tritium transport phenomena involved in MSRs, made suggestions on gaps in the transport dataset, and proposed an experimental test stand to test combined transport effects – tritium transport through pipe walls in a convective salt flow. In this report, we summarize an updated analysis framework for tritium transport in MSRs, report our results on hydrogen and deuterium permeation through Hastelloy N, and describe the final design of the Molten Salt Tritium Transport Experiment (MSTTE, pronounced “misty”). The MSRE provides the only wholistic experimental data set for tritium transport in MSRs and understanding the transport phenomena involved in the MSRE is crucial for future model development. One set of parameters in our analysis framework was unknown for the MSRE—surface reaction rates for tritium on Hastelloy N. This warranted our hydrogen and deuterium permeation campaign to assess the permeability, diffusivity, and solubility of hydrogen isotopes in clean Hastelloy N. Surface reaction rate constants were probed by low pressure measurements, however, no surface effects were observed in the limits of our permeation apparatus. Permeation experiments on oxidized Hastelloy N were not performed for this report but are planned in future work. The experimental test stand, MSTTE, measures combined transport properties of the salt-metal system. MSTTE is a forced convection FLiBe loop with custom designed test section to measure tritium transport through candidate structural materials. We use MSRE relevant dimensionless numbers to design and scale the test section. Hastelloy N is a candidate loop and test section material due to the relevance for the MSRE and related designs, however, other metals are being considered (e.g. 316H SS) which may better align with current vendor concepts.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Lithium‐Containing Hybrid SEI Layer Enabling High Mass Loading and Anode‐Less Sodium Metal Batteries

The continuous rupturing and rebuilding of unstable solid electrolyte interphase (SEI) layer during cycling would block Na + diffusion and induce Na dendrite formation, ultimately limiting the practical application of high-energy-density sodium metal batteries. Herein, a hybrid SEI layer containing Li-species is dexterously constructed on the surface of sodium metal anode. Li-containing inorganic components (Li 3 N, LiF, and Li 2 CO 3 ) are introduced to stabilize the Na/electrolyte interface and enhance the mechanical and diffusion kinetic properties of the SEI layer, which can reduce the side reactions and gas generation, regulate Na + flux during cycling and promote rapid Na + migration for uniform dendrite-free Na deposition. As a result, the constructed Na symmetric cells achieve low overpotential and long cycle life of 5900, 1800, and 500 h at current densities of 3, 10, and 30 mA cm −2 , respectively. Furthermore, the full cells paired with the Na₃V₂(PO₄)₃ cathode demonstrate high specific capacity and excellent cycle stability, even at an ultra-high cathode loading of 39.3 mg cm −2 and a low N/P ratio (negative/positive electrode capacity ratio of 1.21).

Diffusion kinetic↗
DOE OSTI2022

Electrical Analysis of Pulsed Laser Annealed Poly-Si:Ga/SiOx Passivating Contacts

High-efficiency single crystalline silicon (c-Si) solar cells require precise control of dopant diffusion profiles and highly active doping concentrations through thermal annealing. Conventional furnace annealing has been successfully employed for dopant diffusion in aluminum back-surface field (Al-BSF), passivated emitter and rear contact (PERC), and Topcon cells. However, furnace annealing limits some next-generation polycrystalline silicon (poly-Si) on SiOx passivating contacts. This limitation largely affects p-type passivating contacts by not being able to: 1) control the dopant diffusion accurately to prevent B segregation at the SiOx and c-Si interface; and 2) provide highly activated doping concentrations for low contact resistivity. In this contribution, we use a nanosecond excimer laser to examine the passivation quality and electrical performance of both B- and Ga-doped poly-Si/SiOx passivating contacts. The core of this method is to take advantage of the nonequilibrium nature of the anneal through rapid melting and recrystallizing the poly-Si in a short timescale and to achieve doping concentrations above the solid solubility limit. Simulations were performed on polished surfaces to visualize the doping diffusion profiles under different laser conditions, and secondary ion mass spectrometry (SIMS) was used to verify the experimental diffusion profiles post pulsed laser melting (PLM). The results show the dopant diffusion profiles can be tuned precisely through PLM. The electrical analysis using VdP-Hall measurements reveals that the active doping concentration reached 10^21 cm-3 for B and ~2×10^21 cm-3 for Ga in poly-Si, far exceeding their solid solubility limit in Si, with a nearly 100% dopant activation achieved for B, and 20% for Ga. This highly activated dopant profile results in a low contact resistivity of <15 mO·cm^2 for B-doped contacts and <40 mO·cm^2 for Ga-doped contacts. We demonstrate that the passivation quality of the laser annealed contacts shows a linear dependency with the laser energy density and the number of pulses. PLM was also performed on random pyramid textured samples, and we found that preferential melting exists for the tips of pyramids, which can be pursued further for selective pinhole opening. Additionally, we perform thermal stability tests of these hyperdoped passivating contacts and show that under certain laser conditions, the B-doped samples exhibit good thermal stability with near full activation for temperatures up to 800°C. Lastly, we transfer this technique to a larger area by overlapping laser spots and show uniform passivation quality and electrical performance. We are currently implementing this PLM annealed p-type passivating contacts into double-side back-junction poly-Si/SiOx passivating contact devices. Our results show that PLM can not only enable the next generation of high-efficiency c-Si photovoltaics technology but can also benefit the c-Si integrated circuit industry through hyperdoping other semiconductors with atoms that may otherwise have a low solubility.

doping diffusion profile↗
DOE OSTI2025

Hydrogen Diffusion Through Stainless Steel

SAND2025-14430O Hydrogen Diffusion Through Stainless Steel is a tool that employs a finite difference method algorithm to solve Fick's law and other mass transport models pertinent to the diffusion of hydrogen through stainless steel. Additionally, the software uses particle swarm optimization to determine physical parameters based on experimental data. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Mason, Tyler [Sandia National Lab. (SNL-CA), Liver↗
DOE OSTI2023

Self-consistent hopping theory of activated relaxation and diffusion of dilute penetrants in dense crosslinked polymer networks

We generalize a microscopic statistical mechanical theory of the activated dynamics of dilute spherical penetrants in glass-forming liquids to study the influence of crosslinking in polymer networks on the penetrant relaxation time and diffusivity over a wide range of temperature and crosslink fraction (f n ). Our calculations are relevant to recent experimental studies of a nm-sized molecule diffusing in poly-(n-butyl methacrylate) networks. The theory predicts the penetrant relaxation time increases exponentially with the glass transition temperature, T g (f n ), which grows roughly linearly with the square root of fn due to the coupling of local hopping to longer-range collective elasticity. Moreover, T g is also found to be proportional to a geometric confinement parameter defined as the ratio of the penetrant diameter to the mean network mesh size. The decoupling ratio of the penetrant and Kuhn segment alpha times displays a complex non-monotonic dependence on f n and temperature that is well collapsed based on the variable T g (f n )/T. A model for the penetrant diffusion constant that combines activated relaxation and entropic mesh confinement is proposed, which results in a significantly stronger suppression of mass transport with degree of effective supercooling than predicted for the penetrant alpha time. Furthermore, this behavior corresponds to a new network-based type of “decoupling” of diffusion and relaxation. In contrast to the diffusion of larger nanoparticles in high temperature rubbery networks, our analysis in the supercooled regime suggests that for the penetrants studied the mesh confinement effects are of secondary importance relative to the consequences of crosslink-induced slowing down of activated hopping of glassy physics origin.

Chemistry↗