Search NASA⌕ Search

SEARCH · Search NASA

Results for “VISCOSITY”

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 OSTI2024

Containment strategy for subsurface H 2 storage based on time-dependent soft solids

Here, we propose an innovative containment strategy based on time-dependent soft solids, namely 2 wt% Laponite suspensions, to reinforce natural subsurface seals and engineer flow barriers, with an eye toward making H 2 subsurface storage scalable and geographically agnostic. This suspension can be injected at its initial low viscosity and elasticity into a porous medium, allowing for easy pumping and targeted delivery. Once inside the target zone, it matures into a soft solid with much higher viscosity and elasticity, acting as a potential flow barrier. We discuss the rheological properties of the suspensions and demonstrate that hydrogen does not adversely affect their microstructure, but rather increases the suspensions’ viscosity and elasticity. Moreover, the suspensions enhance the rock samples’ compressive strength, while hydrogen exposure increases their stiffness and ductility. The ability of rock samples saturated with the suspensions to contain higher injected gas pressures is enhanced by aging at higher temperatures.

08 HYDROGEN↗
DOE OSTI2026

Insights into Rotational and Translational Dynamics in Mixtures of Ethylene Glycol and Choline Chloride Using Nuclear Magnetic Resonance Techniques

This work examines molecular dynamics and interactions in ethylene glycol–choline chloride (EG–ChCl) mixtures across 0–33 mol % ChCl, spanning the true eutectic region near 17–20 mol % and the commonly used 1:2 formulation. We combine pulsed-field-gradient (PFG) diffusion, fast-field-cycling (FFC) relaxometry, temperature-dependent 13 C T 1 , and nuclear Overhauser effect spectroscopy (NOESY) to disentangle local from macroscopic dynamics. PFG and FFC show that both translational and average rotational motions largely track the strong increase in viscosity with ChCl content, with ethylene glycol consistently diffusing faster than the choline cation and no global dynamical anomaly at the eutectic composition. More subtle, site-specific composition effects nevertheless emerge. The ratio of the diffusion coefficient of the hydroxyl group of choline to the diffusion coefficient of the methyl group of choline displays a shallow minimum in the 17–25 mol % region, indicating a modest change in how the hydroxyl-bearing end of choline samples the underlying translational motion relative to the methyl groups. 13 C T 1 analysis shows that rotational correlation times at 25 °C generally increase with ChCl, reflecting viscosity-coupled slowing, while the CH 2 –N α site exhibits a small but reproducible deviation from this monotonic trend near the eutectic. NOESY spectra at similar compositions reveal enhanced cross-relaxation between EG and choline protons, consistent with increased headgroup–solvent contact density rather than a wholesale structural rearrangement. Overall, our multitechnique study demonstrates that EG–ChCl dynamics are predominantly viscosity-dominated, with the eutectic region acting as a subtle dynamical crossover where specific choline segments become maximally coupled to the hydrogen-bond network. These insights refine the structure–dynamics picture of choline-chloride DESs and provide practical guidance for tuning composition in electrochemical, separation, and catalytic applications.

diffusion↗
DOE OSTI2025

Sustainable upcycling of polyethylene waste to compatibilizers and valuable chemicals

Controllable functionalization of polyethylene (PE) waste could generate new polymeric materials that are generally difficult to manufacture sustainably while also addressing the growing plastics waste problem. However, these modifications remain challenging due to the inherent stability of the PE backbone. Non-thermal atmospheric plasma enables molecular activation under mild conditions while utilizing renewable energy but is primarily employed for surface modification, as plasmas do not penetrate the bulk of materials. Herein, controllable bulk oxidative functionalization of PE wax (PEW) and low-density PE (LDPE) of varying molecular weights was achieved, with up to 6 mol% oxygen incorporation, by manipulating melt viscosity. This functionalization was accomplished either through temperature adjustment or by introducing a melt viscosity modifier, removable via simple extraction methods, to reduce LDPE viscosity, enhance diffusion and chain mobility, and enable bulk oxidation. The oxidized LDPE induces compatibilization in blends of poly(lactic acid) (PLA) and LDPE with improved interfacial adhesion and mechanical properties, such as a 70% increase in elongation-at-break values vs. the control. These findings pave the way for catalyst-free upcycling of direct plastics waste and plastics waste-derived products, enabling the creation of high-value products across various markets.

Plastic upcycling↗
DOE OSTI2025

Two-dimensional hydrodynamic viscous electron flow in annular Corbino rings

The concept of fluidic viscosity is ubiquitous in condensed-matter systems hosting a continuum where macroscopic properties can emerge. While an important property of liquids and some solids, only recently was the viscosity of an electron shown to play a role in electronic transport experiments. In this Letter, we present nonlocal electronic transport measurements in concentric annular rings formed in high-mobility two-dimensional electron gases, and the resulting data show that viscous hydrodynamic flow can occur far away from the source-drain current region. Our conclusion of viscous electronic transport is further corroborated by simulations of the Navier-Stokes equations that are found to be in agreement with our measurements below T = 1 K . Finally, this work emphasizes the key role played by viscosity via electron-electron ( e − e ) interaction even when the electronic transport is restricted radially, and for which it should have played no major role. Published by the American Physical Society 2025

Vijayakrishnan, Sujatha (ORCID:0009000080933182)↗
DOE OSTI2024

Polymer Macrocycles: A novel topology to control dynamics of rubbery materials (Final Report)

One of the current challenges in polymer physics and engineering is to understand the behavior of circular macromolecules. Though such molecules have been made by dilute solution ring closure methods in the past, such procedures have been limited to relatively small, unentangled rings as well as to small sample quantities. Initial studies were also hampered by issues surrounding purity of the rings, though this problem has been somewhat overcome by the use of liquid chromatography at critical conditions (LCCC) methods that separate rings and linear chains. Furthermore, in spite of the continued progress a full understanding of ring dynamics remains elusive. The present work has addressed multiple aspects of the behavior of ring molecules by using a novel synthetic route: reversible radical recombination redox polymerization (R3P) to produce poly(3,6-dioxa-1,8-octanedithiol) (PolyDODT) molecules. The synthesis conditions were chosen to produce linear (LDODT), cyclic (RDODT) and linear-cyclic mixtures (LRDODT) in order to investigate dilute solution and rheological responses of the materials. Importantly, the R3P synthesis makes larger quantities and higher molecular weights than do the dilute solution ring closure methods and also has the potential to be scaled to industrial quantities, thus becoming commercializable. We emphasize two important outcomes of the work. The first is that the rheological investigation shows that very large molecular weight the PolyDODT rings permitted the investigation of ring dynamics from glassy towards terminal flow for molecular sizes corresponding to approximately 300 entanglements of the linear counterpart. At the same time, these very large sized rings could be diluted in a solvent in order to reduce the molecular entanglement density to less than unity. The results demonstrate that ring viscosity follows a Rouse-like linear dependence on molecular weight or entanglement number to approximately 15 entanglements. This range of "Rouse-like" response is greater than what was achieved in all prior works where synthesis conditions have limited the molecular sizes achievable. After the entanglement number $Z_w$ increases beyond approximately 15, we find that the viscosity-molecular weight scaling undergoes a rapid change from the Rouse-like $η\sim{Z}^1_w$ to a very strong power-law scaling that is approximately $η\sim{Z}^{5.8}_w$. This aspect of the work strongly suggests that prior investigations on highly pure rings have invariably been performed on unentangled systems simply because ring molecules have very high entanglement thresholds. The results are further supported by plateau modulus measurements that show that in the "Rouse-like" viscosity regime, for the RDODT samples there is a corresponding lack of a rubbery entanglement plateau and once entanglement coupling is seen, the rubbery plateau also appears in the dynamic modulus data. Furthermore, upon appropriate scaling, we show that the unentangled polyDODT rings (i.e., the diluted high molecular weight rings) show virtually identical dynamics (upon appropriate scaling to account for the chain dilutions) to the dynamics of similarly unentangled LCCC fractionated polystyrene rings.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The Evolution of Accreting Binaries: From Brown Dwarfs to Supermassive Black Holes

Circumbinary accretion occurs throughout the universe, from the formation of stars and planets to the aftermath of major galactic mergers. We present an extensive investigation of circumbinary accretion disks, studying circular binaries with mass ratios (q ≡ M 2 /M 1 ) from 0.01 to 1 and at each mass ratio probing the effects of disk thickness and viscosity. We study disks with aspect ratios H/r $\in$ {0.1, 0.05, 0.03} and vary both the magnitude and spatial dependence of viscosity. Although thin accretion disks have previously been found to promote rapid inspirals of equal-mass binaries, we find that gravitational torques become weaker at lower mass ratios and most binaries with 0.01 ≤ q ≤ 0.04 outspiral, which may delay the coalescence of black hole binaries formed from minor mergers and cause high-mass exoplanets to migrate outward. However, in a number of cases, the disks accreting onto binaries with mass ratios ~0.07 fail to develop eccentric modes, leading to extremely rapid inspirals. Variability in black hole accretion correlates with disk eccentricity, and we observe variability above the ~10% level even for mass ratios of 0.01. We demonstrate that the spatial dependence of the viscosity (e.g., α vs. constant ν) significantly affects the degree of preferential accretion onto the secondary, resolving discrepancies between previous studies. Colder circumbinary disks remain eccentric even at q ~ 0.01 and sustain deep, asymmetric cavities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Completion of Transport Property Measurements on Multiple Actinide Fluoride Mixtures

One of the missions of the US Department of Energy’s Office of Nuclear Energy (DOE-NE) Molten Salt Reactor (MSR) Campaign under the Advanced Reactor Technology program has been to experimentally measure thermophysical properties of MSR-relevent salt systems, with the intent of supporting the development of the Molten Salt Thermal Properties Database (MSTDB). This database is jointly funded by the DOE-NE Nuclear Energy Advanced Modeling and Simulation Program and the MSR Campaign. Multiple DOE national laboratories, including Oak Ridge National Laboratory (ORNL), have been conducting measurements of thermophysical properties to support MSTDB development and provide MSR developers with access to new data that has been measured using modern methodologies and more advanced sample characterization techniques. These data may either fill gaps in the database or provide updated higher quality data to replace legacy data. Researchers at ORNL have recognized significant gaps in the transport property data of actinide-bearing fluoride salt systems of MSR industry interest. Moreover, for the data present in MSTDB in this category, the uncertainty margins are generally high, leading to questionability in our current understanding of the thermophysical characterization of actinide fluoride mixtures. As such, the focus of this study has been to generate new transport property data of actinide fluoride mixtures that are of immediate interest to MSR developers. Specifically, the mixtures NaF-UF 4 (78 - 22 mol%) and NaF-KF-UF 4 (57-16.04-26.91 mol%) have been studied—NaF-UF 4 for thermal conductivity and viscosity and NaF-KF-UF 4 for viscosity. Thermal conductivity measurements have been conducted with a variable gap apparatus, whereas viscosity has been measured with a rolling ball viscometer. Methodological and calibration details are provided for both measurement processes, along with measurement system updates that have enabled easier manufacturing of components and fewer challenges associated with conducting the measurements themselves. The resultant data collected for NaF-UF 4 (78–22 mol%) and NaF-KF-UF 4 (57-16.04-26.91 mol%) have been compared with relevant mixture data within the thermophysical arm of the MSTDB (MSTDB-TP).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Applying colloidal silica suspensions injection and sequential gelation to block vertical water flow in well annulus: laboratory testing on rheology, gelation, and injection

We evaluated the application of silica suspension injection and sequential gelation to block vertical water flow in the annuli of long-screened wells. First, we studied the viscosity, rheological behavior, and gelation performance of colloidal silica suspensions in batch tests. Then, we tested the injection of silica suspensions and the water flow blocking efficiency of the later formed silica gel in column and bench-scale sandbox experiments. Micron-sized fumed powder silica suspensions and nanosized silica suspensions recovered from geothermal fluids were tested in this work. Fumed silica suspensions showed shear thinning, while nanosized silica suspensions exhibited Newtonian flow behavior. During the gelation process, the nanosized silica suspension changed from a Newtonian fluid to a shear thinning fluid while increasing its overall viscosity. At comparable concentrations, the nanosized silica suspensions have much lower viscosity than that of the fumed silica suspensions. Increases in the Na + concentration and silica particle concentration in these suspensions shortened the gelation time. Silica suspension gelation in sand columns completely blocked the water flow and sustained the injection pressure up to 50 psig (344.7 kPa). A silica suspension was successfully injected into the target zone in the annulus of a bench-scale sandbox mimicking long-screened wells in the field. The silica gel formed in the annulus effectively blocked chemical transport through the gelled zone. Our research reveals that a process using silica suspension injection and sequential gelation technology is promising for blocking the vertical water flow and chemical transport through the filter pack in targeted zones within the annulus of long-screened well systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Ab-initio molecular dynamics study of eutectic chloride salt: MgCl2–NaCl–KCl

Ionic liquid materials are viable candidates as a heat transfer fluid (HTF) in a wide range of applications, notably within concentrated solar power (CSP) technology and molten salt reactors (MSRs). For next-generation CSP and MSR technologies that strive for higher power generation efficiency, a HTF with wide liquid phase range and energy storage capabilities is crucial. Studies have shown that eutectic chloride salts exhibit thermal stability at high temperatures, high heat storage capacity, and are less expensive than nitrate and carbonate salts. However, the experimental data needed to fully evaluate the potential of eutectic chloride salts as a HTF contender are scarce and entail large uncertainties. Considering the high cost and potential hazards associated with the experimental methods used to determine the properties of ionic liquids, molecular modeling can be used as a viable alternative resource. In this study, the eutectic ternary chloride salt MgCl 2 –NaCl–KCl is modeled using ab-initio molecular dynamics simulations (AIMDs) in the liquid phase. Using the simulated data, the thermophysical and transport properties of eutectic chloride salt can be calculated: density, viscosity, heat capacity, diffusion coefficient, and ionic conductivity. For an initial model validation, experimental pair-distribution function data were obtained from X-ray total scattering techniques and compared to the theoretical pair-distribution function. Additionally, theoretical viscosity values are compared to experimental viscosity values for a similar system. The results provide a starting foundation for a MgCl 2 –NaCl–KCl model that can be extended to predict other fundamental properties.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

The Evolution of Accreting Binaries: From Brown Dwarfs to Supermassive Black Holes

Circumbinary accretion occurs throughout the universe, from the formation of stars and planets to the aftermath of major galactic mergers. We present an extensive investigation of circumbinary accretion disks, studying circular binaries with mass ratios (q ≡ M 2 /M 1 ) from 0.01 to 1 and at each mass ratio probing the effects of disk thickness and viscosity. We study disks with aspect ratios H/r ϵ {0.1, 0.05, 0.03} and vary both the magnitude and spatial dependence of viscosity. Although thin accretion disks have previously been found to promote rapid inspirals of equal-mass binaries, we find that gravitational torques become weaker at lower mass ratios and most binaries with 0.01 ≤ q ≤ 0.04 outspiral, which may delay the coalescence of black hole binaries formed from minor mergers and cause high-mass exoplanets to migrate outward. However, in a number of cases, the disks accreting onto binaries with mass ratios ~0.07 fail to develop eccentric modes, leading to extremely rapid inspirals. Variability in black hole accretion correlates with disk eccentricity, and we observe variability above the ~10% level even for mass ratios of 0.01. We demonstrate that the spatial dependence of the viscosity (e.g., α vs. constant ν) significantly affects the degree of preferential accretion onto the secondary, resolving discrepancies between previous studies. Colder circumbinary disks remain eccentric even at q ~ 0.01 and sustain deep, asymmetric cavities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Secondary organic aerosols derived from intermediate-volatility n-alkanes adopt low-viscous phase state

Abstract. Secondary organic aerosol (SOA) derived from n-alkanes, as emitted from vehicles and volatile chemical products, is a major component of anthropogenic particulate matter, yet the chemical composition and phase state are poorly understood and thus poorly constrained in aerosol models. Here we provide a comprehensive analysis of n-alkane SOA by explicit gas-phase chemistry modeling, machine learning, and laboratory experiments to show that n-alkane SOA adopts low-viscous semi-solid or liquid states. Our study underlines the complex interplay of molecular composition and SOA viscosity: n-alkane SOA with a higher carbon number mostly consists of less functionalized first-generation products with lower viscosity, while the SOA with a lower carbon number contains more functionalized multigenerational products with higher viscosity. This study opens up a new avenue for analysis of SOA processes, and the results indicate few kinetic limitations of mass accommodation in SOA formation, supporting the application of equilibrium partitioning for simulating n-alkane SOA formation in large-scale atmospheric models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

High-Fidelity and High-Performance Computational Simulations for Rapid Design Optimization of Sulfur Thermal Energy Storage

Industrial process heating (IPH) accounts for approximately 70% of US manufacturing energy use and is primarily produced by fossil fuel combustion. Approximately 1500 TWht (approximately 60%) of IPH demand is in the temperature range of 100-300. Industrial applications in this temperature range include drying, hydrothermal processing, thermal enhanced oil recovery, food and beverage, bioethanol production, etc. Cost-effective thermal energy storage (TES) that increases the utilization of waste and renewable heat (solar, geothermal, etc.) could provide significant energy savings and reliable heat sources, decrease emissions, and increase US manufacturing competitiveness through reductions in fuel consumption. TES development has historically been dominated by technologies suitable for deployment with concentrating solar power (CSP). State-of-the-art thermal storage deployed commercially with power tower CSP plants uses a 60%/40% NaNO3/KNO3 molten salt and operates between temperatures of approximately 280 degrees Celsius and 570 degrees Celsius using a two-tank configuration. However, these nitrate salts are unsuitable for operation outside of this temperature range due to a high freezing point of approximately 220 degrees Celsius, and limits on high-temperature salt stability and corrosion resistance of containment alloys. Other materials being investigated for TES include those based on: (1) sensible energy storage (various molten salt compositions, inert solid particles, rocks or pebble beds, sulfur, water, concrete, graphite, etc.), (2) latent energy storage in materials that undergo solid-liquid phase change at relevant temperatures (organic materials for low-temperature applications, inorganic salts and/or metals for high-temperature applications), or (3) thermochemical energy storage (hydrides, hydroxides, carbonates, metal oxides, etc.). The application temperature and challenges pertaining to storage material and/or containment cost, energy density, long-term thermal and cyclic stability, and charge/discharge heat transfer effectiveness drive material selection for a given IPH or electricity generation application. Sulfur is a cheap commodity at $80/ton compared to $1100 - 1300/ton for conventional salts. When using a metric of storage cost per kWh, sulfur costs around 2-3 $/kWh. Previous sulfur TES development focused on high temperature (>600 degrees) concentrated solar power applications with sulfur encapsulated in pipes and flow of gaseous HTF (air) in the shell side. However, for lower-temperature IPH applications in the range of approximately 100-300 degrees Celsius Element 16 adopted a compact and scalable TES design with molten sulfur in the shell and HTF pipes submerged in the molten sulfur bath. The low-cost molten sulfur TES for dispatchable IPH has deployment potential for broad applications. The spatial and temporal evolution of the HTF and sulfur temperature is critical to the TES system performance, and thus detailed modeling can improve understanding of the performance and facilitate design improvements. Using high performance computing and computational fluid dynamics (CFD) a low-cost molten sulfur thermal energy storage (TES) system for industrial process heating (IPH) applications was developed. The unique challenges in CFD modeling of sulfur TES are the sharp property changes of sulfur relevant to the working temperatures. Above 159, liquid sulfur undergoes polymerization, and the viscosity of sulfur rapidly increases by several orders of magnitude between 159 degrees Celsius and 188 degrees Celsius, followed by a decrease in viscosity beyond 188 degrees Celsius due to thermal bound dissociation. In addition, various concentrations of H2S impurities can also modify sulfur viscosity. This numerical challenge is especially relevant to transient simulation of the sulfur TES charging and discharging processes as the extreme property variations limit the applicability of traditional heat transfer correlations. Transient CFD simulations including the temperature-dependent sulfur properties and geometric complexity of the TES design were used to predict the effect of natural convection during charging and discharging on the heat transfer process, sulfur temperature uniformity, charge/discharge rates, and performance of the storage devices. The CFD model was validated with experimental results for a full charge and discharge cycle. The work will show 3D and 2D simulation comparisons aimed to facilitate rapid design iterations and a machine learning based design optimization approach.

CFD↗
DOE OSTI2024

Langmuir adsorption model to assess the impact of silane coupling on nano-dispersion of silica in SBR

Surface active agents are often used to improve dispersion of nanoparticles. Quantitative correlation between these surface-active molecules and nanoscale dispersion is absent from the literature partly because a quantitative measure of nanoscale dispersion does not exist. Recently, we have developed the Virial-van der Waals method to quantify dispersion in nanocomposites using virial coefficients. In this paper, the Langmuir adsorption model is used to quantify the influence of surface-active agents on nano-scale dispersion in terms of the effective second virial coefficient B 2 *. The impact of silane coupling agent on the nano-dispersion and silica aggregate structure in precipitated silica/SBR nanocomposites is demonstrated. It is shown that the higher viscosity SBR matrix led to a greater silica aggregate structural breakup, while lower viscosity matrix improved surface silanization. The isomeric content of the SBR, which impacts the dielectric behavior, impacted whether the system could be modeled through a mean-field or specific interactions. We earlier showed that larger aggregates improve dispersion, and this is reaffirmed in these results. After account is made for aggregate size, nano-scale dispersion improves with the addition of silane coupling agent. The behavior is well modeled using Langmuir monolayer adsorption.

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

Versatile Liquid Crystal Elastomer Formulations Using Amine-Acrylate Chemistry and Processing for Advanced Manufacturing

This study aims to tune rheological properties of liquid crystal elastomers (LCEs) through a strategic approach, leveraging the differing reaction rates of thiols and primary and secondary amines with acrylate monomers or the inclusion of processing additives to enable amenability to advanced manufacturing. We found that varying the time of oligomerization at elevated temperatures of amine-functional monomers with acrylate-functional monomers enabled simple control of rheological properties. Additionally, we could also control rheological properties by introducing reinforcing fillers or organic solvents. By varying the oligomerization conditions, we achieved a broad range of viscosities, from 1.1 Pa s to 560 Pa s as measured at 60 °C. When the rheological properties were varied through the incorporation of various solid and liquid inclusions, viscosities could range from less than 1 Pa s to nearly 1000 Pa s when measured at 60 °C. The shape morphing capabilities that are intrinsic to LCEs and the mechanical properties were also characterized. Here, we find that the oligomerization conditions and the addition of fillers influence shape morphing and energy absorption, as characterized by the stress–strain characteristics, further suggesting their potential to be used in energy dissipation applications where performance can be turned. Last, we demonstrate how the range of accessible rheological properties grants formulations amenable to various advanced manufacturing techniques.

Liquid chromatography↗
DOE OSTI2023

Shear and combustion characterization of printable ammonium perchlorate composite propellant formulations

Seven formulations of ammonium perchlorate composite propellant (APCP) are developed and their properties relevant for successful additive manufacture are characterized. Extrusion in a custom-built 3D printing system and spindle viscometry are used to collect viscosity measurements of curing primary and non-curing secondary versions of the formulations. The formulations that behave similarly to a Bingham plastic, with apparent viscosities between 4 and 8 thousand Pa*s at 30 minutes post-mix, are determined to be most suitable for printing applications. Measurements of one material show a yield stress of 1 kPa. Ambient pressure burn rates of 2.0–2.9 mm/s were measured for the compositions tested. All measured burn rates were comparable to the lower end of typical burn rates for APCP. In conclusion, the results of the characterization demonstrate a propellant suitable for use in a 3D printing system.

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

Fragile-to-strong transition in liquid As 2 S 3 under pressure: The effect of melt metallization

The well-known classification of glass-forming melts into fragile and strong liquids has several notable exceptions, including water, silica, and certain phase-change materials (PCMs). These exceptional fluid systems exhibit a fragile-to-strong transition (FST) upon cooling: a transformation from a high-temperature liquid with fast atomic dynamics, low viscosity, and low flow activation energy, to a viscous supercooled melt with high energy barriers near the glass transition temperature T g . This behavior is critically important for non-volatile memories, photonic tensor cores, reconfigurable metamaterials, and other devices, that use PCMs, enabling nanosecond-scale crystallization in the fragile regime and long data retention in the strong regime near or below T g . A significant structural transformation is expected between these two viscosity regimes, along with a semiconductor-metal (SC-M) transition upon heating, driven by high internal pressure and associated density increase. By applying high external pressure to the canonical low-conducting chalcogenide melt As 2 S 3 , we observed both the FST and the SC-M transition, occurring simultaneously within the same domain of the P, T−phase space. These findings suggest that the FST is not limited to a few exceptional liquids but is a common phenomenon, at least in systems that exhibit melt metallization within specific regions of their P, T−phase diagrams.

first-principles molecular dynamics↗
DOE OSTI2024

Neutron spin echo shows pHLIP is capable of retarding membrane thickness fluctuations

In this study, cell membranes are responsible for a range of biological processes that require interactions between lipids and proteins. While the effects of lipids on proteins are becoming better understood, our knowledge of how protein conformational changes influence membrane dynamics remains rudimentary. Here, we performed experiments and computer simulations to study the dynamic response of a lipid membrane to changes in the conformational state of pH-low insertion peptide (pHLIP), which transitions from a surface-associated (SA) state at neutral or basic pH to a transmembrane (TM) α-helix under acidic conditions. Our results show that TM-pHLIP significantly slows down membrane thickness fluctuations due to an increase in effective membrane viscosity. Our findings suggest a possible membrane regulatory mechanism, where the TM helix affects lipid chain conformations, and subsequently alters membrane fluctuations and viscosity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗