Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamics of phase separation”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Ion Pairing Mediates Molecular Organization Across Liquid/Liquid Interfaces

Liquid/liquid interfaces play a central role in scientific fields ranging from nanomaterial synthesis and soft matter electronics to nuclear waste remediation and chemical separations. This diversity of functions arises from an interface’s ability to respond to changing conditions in its neighboring bulk phases. Understanding what drives this interfacial flexibility can provide novel avenues for designing new functional interfaces. However, limiting this progress is an inadequate understanding of the subtle intermolecular and interphase interactions taking place at the molecular level. Here, we use surface-specific vibrational sum frequency generation spectroscopy combined with atomistic molecular dynamics simulations to investigate the self-assembly and structure of model ionic oligomers consisting of an oligodimethylsiloxane (ODMS) tail covalently attached to a positively charged methyl imidazolium (MIM + ) head group at buried oil/aqueous interfaces. We show how the presence of seemingly innocuous salts can impart dramatic changes to the ODMS tail conformations in the oil phase via specific ion effects and ion-pairing interactions taking place in the aqueous phase. These specific ion interactions are shown to drive enhanced amphiphile adsorption, induce morphological changes, and disrupt emergent hydrogen-bonding structures at the interface. Tuning these interactions allows for independent control over the oligomer structure in the oil phase versus interfacial population changes and represents key mechanistic insight that is needed to control chemical reactions at liquid/liquid interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Collective behavior of “flexicles”

In recent years the functionality of synthetic active microparticles has edged even closer to that of their biological counterparts. However, we still lack the understanding needed to recreate at the microscale key features of autonomous behavior exhibited by microorganisms or swarms of macroscopic robots. In this study, we propose a model for a three-dimensional deformable cellular composite particle consisting of self-propelled rod-shaped colloids confined within a flexible vesicle—representing a superstructure we call a “flexicle” that couples particle deformation to the internal dynamics of the internal active components. Using molecular dynamics simulations, we investigate the collective behavior of dense systems composed of many flexicles. We show that individual flexicles exhibit shape changes upon collisions with other flexicles that lead to rearrangements of the internal active rods, which slows flexicle motion. This shape deformability gives rise to a diverse set of motility-induced phase separation phenomena and the spontaneous flow of flexicles reminiscent of the migration of cells in dense tissues. Our findings establish a foundation for designing responsive, cell-like active particles and developing strategies for controlling swarm migration and other autonomous swarm behaviors at cellular and colloidal scales.

Science & Technology - Other Topics↗
DOE OSTI2023

Changes in nitrate binding with lanthanides in BLPhen complexes

Preorganized ligands such as bis-lactam-1,10-phenanthroline (BLPhen) offer unique selectivity trends in solvent extraction of rare-earth ions, particularly from aqueous nitrate solutions. However, due to the ligand’s lipophilicity it is experimentally challenging to obtain crystal structures of the corresponding lanthanide (Ln):ligand complexes to describe the local structure of the first coordination shell and the role of the nitrate ions in the organic phase. In this work, we report first principles molecular dynamics (FPMD) simulations of the complexation of trivalent lanthanide ions (La, Nd, and Eu) with the BLPhen ligand in the presence of nitrate ions in the dichloroethane (DCE) solvent. We find that two nitrate anions are present in the first solvation shell and the third nitrate anion is far from the first solvation shell in all three [Ln(BLPhen) 2 ] 3+ complexes examined. Moving along the lanthanide series from La to Nd to Eu, the binding pocket formed by the two BLPhen ligands shrinks in size; as a result, the two nitrate configuration changes from one-monodentate/one-bidentate binding mode for La and Nd to all monodentate for Eu. Meanwhile, the total coordination number drops from 11 for La and Nd to 10 for Eu. More interestingly, the bidentate binding in [La(BLPhen) 2 ] 3+ and [Nd(BLPhen) 2 ] 3+ complexes is highly dynamic, frequently switching to monodentate and back. The FPMD insights into the differing modes and dynamics of nitrates in the first coordination shell will be useful for further atomistic understanding of the complex structure, formation, and stability in the organic phase.

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

Elucidating the speciation of extracted lanthanides by diglycolamides

Many studies over the past few decades have been devoted to addressing the application of diglycolamides (DGAs) for hydrometallurgical-based, f-element separations. Work to date has shown the molecular structure of a DGA derivative can have a significant impact on intra-lanthanide partitioning patterns. More recent studies have pushed towards probing the structure function relationship of the lanthanide-DGA complex to enable the design of more efficient lanthanide separation systems. Spectroscopic techniques, such as UV–Visible, Fourier Transformed Infrared Spectroscopy (FT-IR), Nuclear Magnetic Resonance (NMR), and Extended X-ray Absorption Fine Structure (EXAFS), provide information regarding the inner-sphere coordination of a given lanthanide-DGA complex. Scattering techniques, such as Dynamic Light Scattering (DLS), Small-Angle X-ray Scattering (SAXS), and Small-Angle Neutron Scattering (SANS), address nanoscale structures including aggregate sizes and morphology. Here, this review assesses the current state-of-knowledge regarding lanthanide-DGA hydrometallurgical (i.e., solvent extraction systems) interrogated using spectroscopic and scattering techniques to characterize the extracted Ln 3+ DGA species. Of particular interest to this review is the impact of varied diluents, inclusion and variation of phase modifiers, and DGA derivatization on system characteristics. While there has been extensive literature on the application of DGAs for f-element separations, the literature lacks a collective assessment of the speciation of Ln 3+ in the organic phase. This review provides new insights into the field of DGA separations, explicitly with an application to intra-lanthanide separations. Specifically, this review illustrates the importance of both the co-extraction anion (Cl - , NO 3 - , or ClO 4 - ) as it pertains to both the aggregate size and Ln 3+ distribution coefficient. It is evident the ability of the anion to disrupt the hydrogen bonding network limits both aggregate size and distribution coefficients according to the Hoffmeister series. This suggests the importance of the large, softer anions with a low charge-to-surface area ratio on encouraging hydrogen bond interactions. In addition, the co-extracted cation (H + vs Na + ) is important for mitigating transfer of Ln 3+ from the aqueous to the organic phase through extensive hydrogen bonding networks. These networks are responsible for forming supramolecular aggregates where a change in morphology is observed with increasing concentrations of H + and/or Ln 3+ in the organic phase.

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

Thiocyanate Ions Form Antiparallel Populations at the Concentrated Electrolyte/Charged Surfactant Interface

Anions play significant roles in the separation of lanthanides and actinides. The molecular-scale details of how these anions behave at aqueous interfaces are not well understood, especially at high ionic strengths. Here, we describe the interfacial structure of thiocyanate anions at a soft charged interface up to 5 M bulk concentration with combined classical and phase-sensitive and molecular dynamics (MD) simulations. At low concentrations thiocyanate ions are mostly oriented with their sulfur end pointing toward the charged surfactants. In this work, the VSFG signal reaches a plateau at around 100 mM bulk concentration, followed by significant changes above 1 M. At high concentrations a new thiocyanate population emerges with their sulfur end pointing toward the bulk liquid. The –CN stretch frequency is different for up and down oriented SCN – ions, indicating different coordination environments. These results provide key molecular-level insights for the interfacial behavior of complex anions in highly concentrated solutions.

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

Photoinduced structural recovery dynamics of rare-earth nickelate thin films

Photoinduced structural dynamics of rare-earth nickelate thin films, NdNiO 3 and SmNiO 3 , grown on (001) oriented SrTiO 3 were studied using time-resolved x-ray diffraction. The evolution of the (002) Bragg peak was tracked following laser excitation. The recovery pathways were found to be strongly dependent on laser fluence for NdNiO 3 and distinct for the two rare-earth nickelates. The recovery of the (002) peak shifts was modeled using a one-dimensional thermal diffusion model which showed that the recovery processes are non-thermal at high fluences. For NdNiO 3 , the timescales for the recovery of the (002) peak shift were found to be closely related to Ni magnetism recovery, potentially indicating magneto-structural coupling. Moreover, the evolution of integrated intensity and full-width-half-maximum points towards the presence of a structural phase separation during recovery. In conclusion, our studies highlight structural recovery pathways in nickelates by comparing the distinct responses of the photoinduced metal-insulator transition in NdNiO 3 and SmNiO 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2026

Geometric decoherence time in Lindbladian dynamics

The onset of decoherence in open many-body systems lacks a dynamical timescale grounded in the loss of bipartite entanglement. Here, we introduce the geometric decoherence time, defined as the earliest moment the monotone relation between logarithmic negativity and Rényi-$\frac{1}{2}$ entropy—exactly equal across any bipartition for pure states—breaks down under open-system evolution, signaling entropy growth without accompanying entanglement growth. We establish this criterion in both single-particle Gaussian dynamics and many-body Lindbladian evolution. We show that quantum mutual information provides a complementary long-time diagnostic: Its asymptotic vanishing is equivalent to factorization of the steady state across the bipartition, a condition strictly stronger than separability, and whenever a product steady state is approached exponentially in trace norm, negativity and mutual information share the same decay rate. In the presence of a strong symmetry, this tracking can fail—residual classical correlations can survive after entanglement has vanished. In the Kitaev chain with balanced gain and loss, we derive a closed-form solution and show that the topological phase sustains longer coherence times than the trivial phase at identical dissipation, with a local minimum at the chiral-symmetric point. In the interacting XXZ chain, exact many-body evolution shows that local 𝑍 dephasing preserves residual classical correlations, whereas gain and loss restore the mutual-information tracking of negativity. Furthermore, our results establish the geometric decoherence time as a dynamical scale tracking the onset of decoherence.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Unexpected Kinetic Solvent Effects Enhance Activity and Selectivity in Biphasic Systems

Biphasic dehydration of fructose to 5-hydroxymethylfurfural (HMF) has shown unprecedented increases in productivity, but a mechanistic understanding is lacking. Herein, we couple fast experimental reaction kinetics, multiscale modeling (phase behavior, classical molecular dynamics(MD), and quantum mechanics/molecular mechanics MD), in situ sampling, and IR and 13 C-NMR spectroscopy to elucidate the complex effects of nonpolar extracting organic solvents on the kinetics of fructose dehydration. We show that these organic solvents can reach significant mutual solubility with water at reaction temperatures, enabling the partition of the sugar and catalyst into the extracting phase. In the organic-rich environment, the dehydration of fructose proceeds faster and more selectively than in water due to increased relative abundance of the reactive furanose isomer, enhanced water–catalyst–substrate interactions driven by nanophase separation, and higher product stability stemming from preferential solvation. Furthermore, we demonstrate that these solvent effects impact other critical biphasic reactions in biomass upgrading and provide qualitative principles for solvent selection.

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

Tuning nonequilibrium phase transitions with inertia

In striking contrast to equilibrium systems, inertia can profoundly alter the structure of active systems. Here, we demonstrate that driven systems can exhibit effective equilibrium-like states with increasing particle inertia, despite rigorously violating the fluctuation–dissipation theorem. Increasing inertia progressively eliminates motility-induced phase separation and restores equilibrium crystallization for active Brownian spheres. This effect appears to be general for a wide class of active systems, including those driven by deterministic time-dependent external fields, whose nonequilibrium patterns ultimately disappear with increasing inertia. The path to this effective equilibrium limit can be complex, with finite inertia sometimes acting to accentuate nonequilibrium transitions. The restoration of near equilibrium statistics can be understood through the conversion of active momentum sources to passive-like stresses. Unlike truly equilibrium systems, the effective temperature is now density dependent, the only remnant of the nonequilibrium dynamics. Further, this density-dependent temperature can in principle introduce departures from equilibrium expectations, particularly in response to strong gradients. Our results provide additional insight into the effective temperature ansatz while revealing a mechanism to tune nonequilibrium phase transitions.

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

Twist Angle-Dependent Interlayer Exciton Lifetimes in van der Waals Heterostructures

In van der Waals (vdW) heterostructures formed by stacking two monolayers of transition metal dichalcogenides, multiple exciton resonances with highly tunable properties are formed and subject to both vertical and lateral confinement. We investigate how a unique control knob, the twist angle between the two monolayers, can be used to control the exciton dynamics. We observe that the interlayer exciton lifetimes in MoSe 2 /WSe 2 twisted bilayers (TBLs) change by one order of magnitude when the twist angle is varied from 1° to 3.5°. Additionally, using a low-energy continuum model, we theoretically separate two leading mechanisms that influence interlayer exciton radiative lifetimes. The shift to indirect transitions in the momentum space with an increasing twist angle and the energy modulation from the moiré potential both have a significant impact on interlayer exciton lifetimes. We further predict distinct temperature dependence of interlayer exciton lifetimes in TBLs with different twist angles, which is partially validated by experiments. While many recent studies have highlighted how the twist angle in a vdW TBL can be used to engineer the ground states and quantum phases due to many-body interaction, our studies explore its role in controlling the dynamics of optically excited states, thus, expanding the conceptual applications of “twistronics”.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dielectron production in central Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02TeV

The first measurement of the 𝑒 + ⁢𝑒 − pair production at midrapidity and low invariant mass in central Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV at the Large Hadron Collider is presented. The yield of 𝑒 + ⁢𝑒 − pairs is compared with a cocktail of expected hadronic decay contributions in the invariant mass (𝑚𝑒⁢𝑒) and pair transverse momentum (𝑝 T,𝑒⁢𝑒 ) ranges 𝑚 𝑒⁢𝑒 < 3.5 GeV/𝑐 2 and 𝑝 T,𝑒⁢𝑒 < 8 GeV/𝑐. For 0.18 < 𝑚 𝑒⁢𝑒 < 0.5 GeV/𝑐 2 the ratio of data to the cocktail of hadronic contributions amounts to 1.40 ± 0.11⁢(stat.) ± 0.23 ⁢(syst.) ± 0.16⁢ (cocktail) and 1.42 ± 0.11⁢ (stat.) ± 0.23 (syst.)$^{+0.24}_{−0.29}$ ⁢(cocktail), including or not including medium effects in the estimation of the heavy-flavor background, respectively. It is consistent with predictions from two different models for an additional contribution of thermal 𝑒+⁢𝑒− pairs from the hadronic and partonic phases. In the intermediate-mass range (1.2 < 𝑚 𝑒⁢𝑒 < 2.6 GeV/𝑐 2 ), the pair transverse impact parameter of the 𝑒 + ⁢𝑒 − pairs (DCA𝑒⁢𝑒, where “DCA” denotes “distance of closest approach”) is used for the first time in Pb-Pb collisions to separate displaced dielectrons from heavy-flavor hadron decays from a possible (thermal) contribution produced at the interaction point. The data are consistent with a suppression of 𝑒 + ⁢𝑒 − pairs from $𝑐\bar{𝑐}$ and an additional prompt component. Finally, the first direct-photon measurement in the 10% most central Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV is reported via the study of virtual direct photons in the transverse momentum range 1 < 𝑝 T < 5 GeV/𝑐. A model including prompt photons, as well as photons from the preequilibrium and fluid-dynamic phases, can reproduce the result, while being at the upper edge of the data uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Measurement of the W-boson angular coefficients and transverse momentum in pp collisions at s=13 TeV with the ATLAS detector

The angular distributions of Drell–Yan lepton pairs provide sensitive probes of the underlying dynamics of quantum chromodynamics (QCD) effects in vector-boson production. This paper presents for the first time the measurement of the full set of angular coefficients together with the differential cross-section as a function of the transverse momentum of the W boson, in the full phase space of the decay leptons. The measurements are performed separately for the W-$$W^-$$ and W+$$W^+$$ channels. The analysis uses proton–proton collision data recorded by the ATLAS experiment at the Large Hadron Collider in 2017 and 2018, during special low-luminosity runs with a reduced number of interactions per bunch crossings (pile-up). The data correspond to an integrated luminosity of 338 pb-1$$^{-1}$$ at a centre-of-mass energy of s=13$$\sqrt{s} = 13$$ TeV. The low pile-up environment provides excellent experimental conditions for high-precision measurements of W-boson production. All results agree with theoretical predictions incorporating finite-order QCD corrections up to order αS2$$\alpha _S^2$$.

Aad, G↗
DOE OSTI2021

Phase segregation mechanisms of small molecule-polymer blends unraveled by varying polymer chain architecture

As phase separation between the small-molecule semiconductor and the polymer binder is the key enabler of blend-based organic field-effect transistors (OFETs) fabricated by low-cost solution processing, it is crucial to understand the underlying phase separation mechanisms that determine the phase morphology, which significantly impacts device performance. Beyond the parameter space investigated in previous work, here we investigate the formation of blends by varying the branch architecture of the polymer binder and by shortening the solvent dry time using ultrasonic spray casting. The phase morphologies of the resulting blend films have been thoroughly characterized with a variety of techniques in three dimensions over multiple length scales, including AFM, energy-filtered transmission electron microscope, and neutron reflectivity, and have been correlated with electrical transport performance. From the results, we have inferred that the phase morphology is kinetically determined, limited by the inherent slow movement of polymer macromolecules. The kinetic picture, supported by molecular dynamics modeling, not only consistently explains our observations but also resolves inconsistencies in previous works. The achieved mechanistic understanding will guide further optimization of blend-based organic electronics, such as OFETs and organic photovoltaics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Excitation control of dual-wound machines for coupling mitigation

A dual-wound machine comprises a dual-wound generator supplying power to two separate powered zones. The generator comprises a wound rotor with a field winding and a stator with two sets of phase windings and a field control loop that controls the excitation voltage applied to the field winding and therefore the magnetic field produced by the rotor, in order to maintain a constant field flux in the generator and mitigate dynamic coupling between the two sets of phase windings when supplying power to unbalanced loads.

Rashkin, Lee Joshua↗
DOE OSTI2020

Enhancing building energy performance by effectively using phase change material and dynamic insulation in walls

Deploying phase change materials (PCMs) in building envelopes can be efficacious in reducing space heating/cooling loads and providing load shedding and shifting capacity. However, the full potential of PCM-integrated envelopes can only be harnessed if the PCM undergoes phase change using free ambient heating/cooling, and the stored energy is effectively transferred between the exterior and the interior environments. Traditional thermal insulation (with a fixed thermal resistance) limits PCM utilization, which restrains the energy saving potential of a PCM-integrated envelope to a small percentage. Proposed dynamic insulation material and system (DIMS) provides the option of varying its thermal resistance based on the indoor and outdoor conditions. Although it has been established that employing PCM as well as DIMS in building envelopes separately improves buildings’ energy performance, no prior studies that analyzed the combined influence of both technologies were identified. In this study, we examine a novel wall design, comprising a layer of PCM between two layers of DIMS. We note that the PCM-DIMS-integrated wall provides significantly higher energy saving potential than the DIMS-only integrated wall or the PCM-only integrated wall in all the climates and wall orientations analyzed in this study. Depending on the climate, the PCM-DIMS-integrated wall could provide 15–72% reduction in annual heat gain and 7–38% reduction in annual heat loss. The analysis presented in this study supports the need to develop scalable dynamic insulations combined with thermal energy storage systems for buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Hydrodynamic instabilities and heat transfer characteristics in the duct flow of a fluid in the supercritical thermodynamic regime

The behavior of fluids at supercritical thermodynamic conditions is inherently complex due to large variations in thermodynamic and transport properties. Recent numerical and experimental investigations illustrate ongoing interest for these fluids, especially supercritical CO 2 and supercritical water, for a variety of applications. For example, supercritical water reactors (SCWR) operate in this extreme condition of high-pressure and temperature, resulting in highly dynamic flow fields and unexpected heat transfer regimes. The potential heat transfer benefits in this regime are directly associated with the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. This work characterizes the hydrodynamic instabilities that arise for fluids at supercritical thermodynamic conditions when buoyancy forces are significant. Two specific configurations are considered, a natural convection cavity flow, and a mixed convection, heated, horizontal channel flow. Natural convection flow in a cavity is a classical configuration with expected behavior below the critical point. This configuration aids in characterizing the effect of the variable properties in the supercritical thermodynamic regime. Further, limited studies in the existing literature have been conducted for low-Reynolds and intermediate-Rayleigh numbers, mixed-convection channel flows for supercritical water, which is the focus of the channel flow configuration. To investigate the thermally driven hydrodynamic instabilities in this regime, a high-order fully-implicit numerical method is used. Such strong variations in thermophysical properties (in particular, density) are difficult to simulate and an altogether compressible framework is needed. Therefore, the compressible Navier-Stokes equations are solved without any additional assumptions. The fully implicit, high-order in space and time, reconstructed discontinuous Galerkin method as implemented within the multi-physics code called ALE3D (Arbitrary Lagrangian and Eulerian in 2D and 3D), developed at Lawrence Livermore National Laboratory (LLNL), is used. This fully implicit, L-stable method accurately captures the compressible nature of the ow in the limit of very low Mach number. It has been widely accepted that above the critical point, only one phase is observed. However, recent research has indicated the existence of the distinct gas-like and liquid-like regions separated by the Widom line, the locus of the maxima of the specific heat. Along the Widom line, density decreases 6-fold, viscosity drops by a factor of 2, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film boiling. A pseudo-film at the heated wall of the cavity and the horizontal channel is observed where buoyancy forces induce mixing through the specific configurations. Further the local Rayleigh and Richardson numbers provide maps of the flow field and the buoyancy forces driving the microscopic mixing. In the first chapter, I describe a background of supercritical fluid and the various applications. The second chapter focuses on the mathematical model and numerical method used for simulations, where a description of the equation of state for supercritical water is described. The third chapter focuses on the natural convection cavity with a heated bottom wall. In this cavity a gas-like and a liquid-like flow within the supercritical thermodynamic regime are observed. The fourth chapter focuses on a forced convection, horizontal channel, distinguishing between the gas-like, liquid-like, and mixed flow regimes. Mixed convection flow, with the addition of gravitational forces in the horizontal channel show the influence of variable properties on the hydrodynamic development, heat transfer, and rising instabilities. The last chapter of this research focuses on characterizing the unstable hydrodynamics through time-averaging processes and analysis of the movement of energy through the developing plumes.

42 ENGINEERING↗
DOE OSTI2024

Consistent and reproducible computation of the glass transition temperature from molecular dynamics simulations

In many fields, from semiconductors for opto-electronic applications to ionic liquids (ILs) for separations, the glass transition temperature (Tg) of a material is a useful gauge for its potential use in practical settings. As a result, there is a great deal of interest in predicting Tg using molecular simulations. However, the uncertainty and variation in the trend shift method, a common approach in simulations to predict Tg, can be high. This is due to the need for human intervention in defining a fitting range for linear fits of density with temperature assumed for the liquid and glass phases across the simulated cooling. The definition of such fitting ranges then defines the estimate for the Tg as the intersection of linear fits. We eliminate this need for human intervention by leveraging the Shapiro–Wilk normality test and proposing an algorithm to define the fitting ranges and, consequently, Tg. Through this integration, we incorporate into our automated methodology that residuals must be normally distributed around zero for any fit, a requirement that must be met for any regression problem. Consequently, fitting ranges for realizing linear fits for each phase are statistically defined rather than visually inferred, obtaining an estimate for Tg without any human intervention. The method is also capable of finding multiple linear regimes across density vs temperature curves. We compare the predictions of our proposed method across multiple IL and semiconductor molecular dynamics simulation results from the literature and compare other proposed methods for automatically detecting Tg from density–temperature data. We believe that our proposed method would allow for more consistent predictions of Tg. We make this methodology available and open source through GitHub.

Chemistry↗
DOE OSTI2025

Equilibrium Separation of Siloxanes in Metal–Organic Frameworks

We present an in silico assessment of metal–organic frameworks (MOFs) for the equilibrium separation of linear and cyclic siloxanes. Using a combination of configurational bias/continuous fractional component Monte Carlo (CB/CFCMC) simulations and Ideal Adsorbed Solution Theory (IAST), we investigated the adsorption of both equimolar and nonequimolar mixtures of linear and cyclic siloxanes in a selection of synthesizable MOFs with medium to large pore volumes. We showed that configurational entropy effects drive the preferential adsorption of linear siloxanes over cyclic siloxanes. Based on synthesizability metrics we identified ZIF-70 as a promising adsorbent for the separation of linear and cyclic siloxanes. Self-diffusivities of linear and cyclic siloxanes in ZIF-70 calculated from molecular dynamics simulations show that equilibrium can be reached on reasonable time scales. We explored vacuum-temperature swing adsorption (VTSA) as a potential method for the recovery of adsorbed linear siloxanes from ZIF-70, achieving a 75% reduction in adsorbed phase concentration. Additionally, we demonstrated that supercritical CO 2 offers an alternative desorption strategy by displacing adsorbed linear siloxanes in ZIF-70 at pressures above 200 bar driven by entropy effects.

Adsorption↗