SEARCH · Search NASA

Results for “molecular dynamics (MD)”

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.

SourceResearch connectionFrom yearThrough year

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

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Natural Charge-Transfer Analysis: Eliminating Spurious Charge-Transfer States in Time-Dependent Density Functional Theory via Diabatization, with Application to Projection-Based Embedding

For many types of vertical excitation energies, linear-response time-dependent density functional theory (LR-TDDFT) offers a useful degree of accuracy combined with unrivaled computational efficiency, although charge-transfer excitation energies are often systematically and dramatically underestimated, especially for large systems and those that contain explicit solvent. As a result, low energy electronic spectra of solution-phase chromophores often contain tens to hundreds of spurious charge-transfer states, making LR-TDDFT needlessly expensive in bulk solution. More nefariously, intensity borrowing by the low-energy charge-transfer states can affect intensities of the valence excitations even if those excitation energies are accurate. At higher excitation energies, it is difficult to distinguish spurious CT states from genuine charge-transfer-to-solvent (CTTS) excitations. In this work, we introduce an automated diabatization scheme that enables fast and effective screening of the CTTS acceptor space in bulk solution. Our procedure introduces the concept of “natural charge-transfer orbitals”, which provide a means to isolate characteristic pairs of orbitals that are most likely to participate in a CTTS excitation. The projection of these orbitals onto solvent-centered virtual orbitals provides a criterion for defining the most important solvent molecules in a given excitation. We apply this method to analyze an ab initio molecular dynamics (MD) trajectory of I-(aq) and report the lowest-energy CTTS band in the absorption spectrum. Our results are in excellent agreement with experimental measurements for bulk I-(aq), and only one-third of the water molecules in the I-(H2O)96 simulation cell need to be described with LR-TDDFT in order to obtain excitation energies that are converged to < 0.1 eV. The tools introduced herein will improve the accuracy, efficiency, and usability of LR-TDDFT in solution-phase environments

Carter-Fenk, Kevin D.↗

DOE OSTI2024

Effect of CO 2 on the water slip flow at silica surfaces for nanometer slit pores of talc

In our previous study, Molecular Dynamics (MD) simulated water flow with no dissolved gas in nanometer slit pores of the hydrophobic talc surfaces revealed a slip length of about 0.5 nm, which is close to the size of “water exclusion zone” but much less than experimental slip length. Atomic Force Microscopy (AFM) imaging confirmed the presence of pancake shape nanobubbles at the hydrophobic talc (001) surface (size from tens to hundreds of nanometers). The effect of CO 2 on water slip flow in talc nanopores is reported to reveal the importance of dissolved gas for confined fluid flow. MD simulated 7 nm CO 2 nanobubbles were found to attach and spread at the talc (001) surface, consistent with the AFM observation. Simulated CO 2 saturated water flow in the slit pores of talc (001) surfaces predicted the generation of CO 2 nanobubbles at the end of the slit pore, due to the stabilization of a CO 2 film at the talc (001) surfaces and a water barrier at the edge surfaces of talc. A “critical thickness” of about 1.5 nm was found for CO 2 nanobubbles to be stabilized at the end of the talc (001) surface. With 1.5 nm CO 2 nanobubbles at the talc (001) surface of a 6 nm slit pore, a water slip length of 1.8 nm was determined. Furthermore, the nanometer slit pore simulations with dissolved CO 2 in water revealed an increased slip length. It is expected that an increased slip length value would be found with an increased CO 2 film thickness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Structural analyses of the PKA RIIβ holoenzyme containing the oncogenic DnaJB1-PKAc fusion protein reveal protomer asymmetry and fusion-induced allosteric perturbations in fibrolamellar hepatocellular carcinoma

When the J-domain of the heat shock protein DnaJB1 is fused to the catalytic (C) subunit of cAMP-dependent protein kinase (PKA), replacing exon 1, this fusion protein, J-C subunit (J-C), becomes the driver of fibrolamellar hepatocellular carcinoma (FL-HCC). Here, we use cryo-electron microscopy (cryo-EM) to characterize J-C bound to RIIβ, the major PKA regulatory (R) subunit in liver, thus reporting the first cryo-EM structure of any PKA holoenzyme. We report several differences in both structure and dynamics that could not be captured by the conventional crystallography approaches used to obtain prior structures. Most striking is the asymmetry caused by the absence of the second cyclic nucleotide binding (CNB) domain and the J-domain in one of the RIIβ:J-C protomers. Using molecular dynamics (MD) simulations, we discovered that this asymmetry is already present in the wild-type (WT) RIIβ 2 C 2 but had been masked in the previous crystal structure. This asymmetry may link to the intrinsic allosteric regulation of all PKA holoenzymes and could also explain why most disease mutations in PKA regulatory subunits are dominant negative. The cryo-EM structure, combined with small-angle X-ray scattering (SAXS), also allowed us to predict the general position of the Dimerization/Docking (D/D) domain, which is essential for localization and interacting with membrane-anchored A-Kinase-Anchoring Proteins (AKAPs). This position provides a multivalent mechanism for interaction of the RIIβ holoenzyme with membranes and would be perturbed in the oncogenic fusion protein. The J-domain also alters several biochemical properties of the RIIβ holoenzyme: It is easier to activate with cAMP, and the cooperativity is reduced. These results provide new insights into how the finely tuned allosteric PKA signaling network is disrupted by the oncogenic J-C subunit, ultimately leading to the development of FL-HCC.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Nanocompression of 20 nm Silver Nanoparticles: In situ Aberration‐Corrected TEM and Atomistic Simulations

Single‐crystalline nanoparticles play an increasingly important role in a wide variety of fields including pharmaceuticals, advanced materials, catalysts for fuel cells, energy materials, as well as environmental detection and monitoring. Yet, the deformation mechanisms of very small nanoparticles are still poorly understood, in particular the role played by single dislocations and their interaction with surfaces. Here, in this work, silver nanoparticles with particularly small dimensions (≈20 nanometers in diameter) are compressed in situ in an aberration‐corrected transmission electron microscopy (TEM) and molecular dynamics (MD) simulations. During compression, the emergence of both dislocations and nanotwins are observed. However, these defects prove to be unstable and disappear upon removal of the indenter. Atomistic simulations confirm the role played by image stresses associated with the nearby surfaces and the reduction in dislocation line length as it approaches the free surface, thereby supporting the experimental observations. These results provide justification for the frequent observation of the absence of dislocations in nanoparticles of a few nanometers in size during in situ experiments, even after significant deformation. This phenomenon contributes to the self‐healing of samples through dislocation ejection toward the surfaces.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Disentangling the evolution of electrons and holes in photoexcited ZnO nanoparticles

The evolution of charge carriers in photoexcited room temperature ZnO nanoparticles in solution is investigated using ultrafast ultraviolet photoluminescence spectroscopy, ultrafast Zn K-edge absorption spectroscopy, and ab initio molecular dynamics (MD) simulations. The photoluminescence is excited at 4.66 eV, well above the band edge, and shows that electron cooling in the conduction band and exciton formation occur in <500 fs, in excellent agreement with theoretical predictions. The x-ray absorption measurements, obtained upon excitation close to the band edge at 3.49 eV, are sensitive to the migration and trapping of holes. They reveal that the 2 ps transient largely reproduces the previously reported transient obtained at 100 ps time delay in synchrotron studies. In addition, the x-ray absorption signal is found to rise in ∼1.4 ps, which we attribute to the diffusion of holes through the lattice prior to their trapping at singly charged oxygen vacancies. Indeed, the MD simulations show that impulsive trapping of holes induces an ultrafast expansion of the cage of Zn atoms in <200 fs, followed by an oscillatory response at a frequency of ∼100 cm−1, which corresponds to a phonon mode of the system involving the Zn sub-lattice.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

An atomistic study connecting underlying dislocation behavior with superior mechanical properties of NiCoCr medium entropy alloy

NiCoCr-based medium-entropy alloy (MEA) with a simple face-centered cubic crystal phase exhibits excellent mechanical properties, often attributed to the synergy of multiple deformation mechanisms. However, the atomistic origin of their outstanding mechanical response, including microstructural evolution and dislocation behavior under varying strain-rates and orientation, remains unclear. In this work, we employ large-scale molecular dynamics (MD) simulations to investigate the changes in deformation mechanisms along three distinct orientations ([110], [111], [100]) under varying strain rates (1 ×10 8 /sec, 1 ×10 10 /sec, 1 ×10 12 /sec) in the NiCoCr MEA. The presence of the stair-rod and the Shockley partial dislocations under uniaxial tensile strain are found to play a key role in the formation of deformation twinning and ε-martensite, which positively correlates with strain-rate dependent dislocation analysis. These findings further establish the role of the dislocations in controlling the superior mechanical response and excellent fracture toughness of the NiCoCr MEA. Systematic transmission-electron microscopy tests performed on the [111]-oriented crystals, deformed at different strain levels, at room temperature provide clear evidence of both the extended stacking-fault and the stair rods, confirming the predicted microstructural features. Finally, this study offers key insights into the complex nucleation mechanisms of deformation twinning and ε-martensite, such as twinning – and transformation–induced plasticity (TWIP-TRIP), providing valuable guidelines for studying similar material classes.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Formation of large loop-hydrogen complexes and related effects on mechanical properties of Zirconium investigated with molecular dynamics method

Here, using molecular dynamics (MD) method, interaction between hydrogen (H) atoms and interstitial or vacancy dislocation loops in Zr is extensively studied. The results show that the binding of hydrogen atoms to dislocation loops is anisotropic (attractive or repulsive), depending on the position of the hydrogen atoms with respect to the loop’s core. The anisotropy is governed by excess volume and stress field induced by the loop. For a 1/3[$2\bar1$$\bar1$$0$] interstitial dislocation loop, hydrogen atoms preferentially segregate near the outside of the loop and close to the decomposition knots of the loop. While for a 1/3[$2\bar1$$\bar1$$0$] or 1/6[$20$$\bar2$$3$] vacancy loop, hydrogen atoms prefer to segregate near the inside of the loop. The present work indicates large loop-hydrogen complexes could form in Zr. The binding of hydrogen to loops is stable up to 600 K, while dissociation is observed at 900 K. Formation of loop-hydrogen complexes affect the activation of slip system, resulting in higher yield stress of Zr. The effect increases with the increasing ratio of H to vacancies or interstitials in the vacancy or interstitial loops. The results imply the effects of large loop-hydrogen complexes should be considered to fully understand radiation damage in Zr.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Nanotwinned alloys under high pressure

Nanotwinned alloys are of interest due to their high strength and ductility, but twin boundaries may not be stable under shear. Computational studies indicate that high hydrostatic pressure may suppress detwinning mechanisms. Here, in this study, we investigate the microstructural changes of nanotwinned-nanocrystalline copper-nickel and Inconel 725 alloys under quasi-hydrostatic pressures up to 50 gigapascals (GPa). The alloys are compressed in a diamond anvil cell. In-situ x-ray diffraction (XRD) and ex-situ transmission electron microscopy (TEM) were employed to monitor microstructural changes. Twin boundary deformation and grain growth occur at 11.4 GPa quasi-hydrostatic pressure in the copper-nickel alloy. Molecular dynamics (MD) simulations reveal that hydrostatic pressure causes elevated local shear stress at grain boundaries, which leads to atomic rearrangements. A superposition of hydrostatic and deviatoric pressures lead to partial dislocation mediated twin boundary migration. In contrast, the Inconel 725 alloy showed stable twin and grain boundaries up to a quasi-hydrostatic pressure of 12.7 GPa. Texture, high solid solution strengthening, and low stacking fault energy are hypothesized to the enhanced microstructural stability in Inconel 725.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Molecular-dynamics study of diffusional creep in uranium mononitride

Uranium mononitride (UN) is a promising advanced nuclear fuel due to its high thermal conductivity and high fissile density. Yet, many aspects of its mechanical behavior and microstructural features are currently unknown. In this paper, molecular dynamics (MD) simulations are used to study UN's diffusional creep. Nanometer-sized polycrystals are used to simulate diffusional creep and to calculate an effective GB width. It is found that Nabarro-Herring creep is not dominant in the temperature range of 1700–2000 K and that the dominant diffusional creep mechanism is Coble creep with an activation energy of 2.28 ± 0.09 eV. A method is proposed to calculate the diffusional GB width and its temperature dependence in polycrystals. The effective GB width of UN is calculated as 2.69 ± 0.08 nm. This value fits very well with the prefactor of the phenomenological Coble creep formula. It is demonstrated that the most comprehensive thermal creep model for UN can be represented as the combination of our Coble creep model and the dislocation creep model proposed by Hayes et al.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

Solvent oligomerization pathways facilitated by electrolyte additives during solid-electrolyte interphase formation

The solid–electrolyte interphase (SEI) layer formation is known to play an important role in determining the lifetime of lithium-ion batteries. A thin, stable SEI layer is linked to overall improved battery performance and longevity, however, the factors and mechanisms that lead to optimal SEI morphology and composition are not well understood. Inclusion of electrolyte additives (fluoroethylene carbonate, FEC; and vinylene carbonate, VC) is often necessary for improving SEI characteristics. To understand how these electrolyte additives impact SEI formation, we employed molecular dynamics (MD) and density functional theory (DFT) simulations to study the reaction networks and oligomerization pathways, respectively, for three systems containing ethylene carbonate (EC), a lithium ion, and FEC or VC. MD simulations suggest radical oligomerization pathways analogous to traditional oligomerization with nucleophilic alkoxide species via SN1 reaction mechanisms. Both S N 1 and S N 2 mechanisms were studied for all three systems using DFT. Oligomerization reactions were studied with both a standard alkoxide species and a ring-opened EC radical as the nucleophiles and EC, FEC, and VC as the electrophiles. For all cases, FEC and VC exhibited lower free energy barriers and more stable adducts when compared with EC. We conclude that one of the role of additives is to modify the oligomerization process of EC by introducing branching points (FEC) or termination points (VC).

25 ENERGY STORAGE↗

DOE OSTI2024

Effect of stoichiometry and hydration level on water domain size and transport in poly(aryl piperidinium) alkaline anion-exchange membranes

Alkaline water electrolysis holds promise in decarbonizing the global economy by enabling renewable hydrogen production with non-precious group metal catalysts. Anion exchange membranes are an important component of alkaline water electrolyzers and would ideally be durable while allowing for high hydroxide conductivity. The poly(aryl piperidinium) (PAP) class of polymers has attracted recent interest due to their good mechanical robustness and high ionic conductivity. Here, in this work, we perform atomistic molecular dynamics (MD) simulations of several PAP polymers at experimentally relevant hydration levels and polymer ion exchange capacities (IECs) to gain nanoscale insight into their properties and to help elucidate the trade-offs that result from tuning the IECs through the polymer stoichiometry. Our MD-predicted macroscopic polymer properties were found to be in good agreement with experimentally available polymer swelling ratios, water-occupied volumes, X-ray scattering, and ionic conductivities. The models show that for hydration levels greater than 8H 2 O per cation a single water cluster will form that percolates through the system. The growth in water cluster size results in large polymer swelling, the creation of larger channels with widths of 7 Å or larger, and nanophase separation between the hydrophilic domains and the polymer with characteristic length scales of approximately 20–30 Å. The experimentally observed lack of a strong X-ray scattering peak at low wavevectors can be explained by a cancellation between the polymer-polymer/water-water and polymer-water correlations and not a loss in nanophase separation. The overlap in coordination environments of the hydroxide oxygen and polymer nitrogen atoms implies that vehicular diffusion between cationic groups could play a role in hydroxide transport. The polymers' hydroxide and water diffusion constants increase by approximately an order of magnitude between hydration levels of 8 and 20H 2 O per cation. However, there are diminishing returns in hydroxide diffusion constant once the IEC exceeds 2.4 meq/g.

08 HYDROGEN↗

DOE OSTI2024

Influence of the Dielectric Constant on the Ionic Current Rectification of Bipolar Nanopores

In this paper, we investigate how the dielectric constant, ϵ, of an electrolyte solvent influences the current rectification characteristics of bipolar nanopores. It is well recognized that bipolar nanopores with two oppositely charged regions rectify current when exposed to an alternating electric potential difference. Here, we consider dilute electrolytes with NaCl only and with a mixture of NaCl and charged nanoparticles. These systems are studied using two levels of description, all-atom explicit water molecular dynamics (MD) simulations and coarse-grained implicit solvent MD simulations. The charge density and electric potential profiles and current-voltage relationship predicted by the implicit solvent simulations with ϵ = 11.3 show good agreement with the predictions from the explicit water simulations. Under nonequilibrium conditions, the predictions of the implicit solvent simulations with a dielectric constant closer to the one of bulk water are significantly different from the predictions obtained with the explicit water model. Further, these findings are closely aligned with experimental data on the dielectric constant of water when confined to nanometric spaces, which suggests that ϵ decreases significantly compared to its value in the bulk. Moreover, the largest electric current rectification is observed in systems containing nanoparticles when ϵ = 78.8. Using enhanced sampling, we have shown that this larger rectification arises from the presence of a significantly deeper minimum in the free energy of the system with a larger ϵ, and when a negative voltage bias is applied. Since implicit solvent models and mean-field continuum theories are often used to design Janus membranes based on bipolar nanopores, this work highlights the importance of properly accounting for the effects of confinement on the dielectric constant of the electrolyte solvent. The results presented here indicate that the dielectric constant in implicit solvent simulations may be used as an adjustable parameter to approximately account for the effects of nanometric confinement on aqueous electrolyte solvents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Coarse-grained molecular dynamics simulation of solvent-dependent cellulose nanofiber interactions

Associations between cellulose are important both in biofuel production and in the use of cellulose for biomaterials. Cellulose nanofibers (CNFs) are sustainable, strong, light-weight alternatives to traditional materials in manufacturing, but are challenging to obtain due to irreversible aggregation in solution during preparative fibrillation. Therefore, it is imperative to understand the underlying factors driving aggregation with a view to designing solvents that can effectively compete with interfiber interactions, hence reducing aggregation. Molecular dynamics (MD) simulation at atomic detail can provide useful information on local interactions. However, the length and timescales accessible are too short to fully capture association processes. Here, we provide a method for accessing the longer length and timescales required using coarse-grained (CG) MD simulations with a MARTINI force field to calculate the interaction behavior of CNFs in three selected solvents: NaOH-urea-water, acetone, and neat water. The CG results are consistent with our prior all-atom MD and with previous experimental results. While acetone is found not to be an effective solvent, urea and ionic moieties in NaOH-urea-water not only solvate the fibrils but also improve the confinement of water molecules around them as shown by the solvent residence times and mean-square displacements. Overall, the presence of urea and ions reduces the likelihood of aggregation in multi-CNF systems relative to neat water irrespective of whether the hydrophobic or hydrophilic CNF surfaces are interacting. In conclusion, the CG method shows clear promise for selecting potential high-performance solvents for experimental prioritization in bioenergy and biomaterials research in a relatively fast manner as well as for understanding the aggregation and rheological behavior of CNF-solvent systems.

aggregation↗

DOE OSTI2022

Molecular dynamics study of silicon atomic layer etching by chlorine gas and argon ions

Classical molecular dynamics (MD) is used to simulate atomic layer etching processes of silicon by alternating exposure to chlorine gas and argon ions. In order to validate our model, a rigorous comparison is done with ion beam experiments found in the literature [Park et al., Jpn. J. Appl. Phys. 44, 389 (2005)]. Here it is shown that the etch per cycle (EPC) as a function of argon ion energy from simulations is in quantitative agreement with experimental results if the correct argon ion fluence is used in the simulations. The EPC as a function of ion irradiation time and amount of chlorine exposure also show good agreement with the experiment. The MD simulations also show the formation of an amorphous silicon region with chlorine atoms mixed uniformly throughout following ion bombardment. Finally, the etch products during the ion irradiation step are analyzed and discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2022

Molecular simulations and NMR reveal how lipid fluctuations affect membrane mechanics

Lipid bilayers form the main matrix of functional cell membranes, and their dynamics underlie a host of physical and biological processes. Here we show that elastic membrane properties and collective molecular dynamics (MD) are related by the mean-square amplitudes (order parameters) and relaxation rates (correlation times) of lipid acyl chain motions. We performed all-atom MD simulations of liquid-crystalline bilayers that allow direct comparison with carbon-hydrogen (CH) bond relaxations measured with NMR spectroscopy. Previous computational and theoretical approaches have assumed isotropic relaxation, which yields inaccurate description of lipid chain dynamics and incorrect data interpretation. Instead, the new framework includes a fixed bilayer normal (director axis) and restricted anisotropic motion of the CH bonds in accord with their segmental order parameters, enabling robust validation of lipid force fields. Simulated spectral densities of thermally excited CH bond fluctuations exhibited well-defined spin-lattice (Zeeman) relaxations analogous to those in NMR measurements. Their frequency signature could be fit to a simple power-law function, indicative of nematic-like collective dynamics. Moreover, calculated relaxation rates scaled as the squared order parameters yielding an apparent K C modulus for bilayer bending. Our results show a strong correlation with K C values obtained from solid-state NMR studies of bilayers without and with cholesterol as validated by neutron spin-echo measurements of membrane elasticity. The simulations uncover a critical role of interleaflet coupling in membrane mechanics and thus provide important insights into molecular sites of emerging elastic properties within lipid bilayers.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Predicting the structural basis of targeted protein degradation by integrating molecular dynamics simulations with structural mass spectrometry

Targeted protein degradation (TPD) is a promising approach in drug discovery for degrading proteins implicated in diseases. A key step in this process is the formation of a ternary complex where a heterobifunctional molecule induces proximity of an E3 ligase to a protein of interest (POI), thus facilitating ubiquitin transfer to the POI. In this work, we characterize 3 steps in the TPD process. (1) We simulate the ternary complex formation of SMARCA2 bromodomain and VHL E3 ligase by combining hydrogen-deuterium exchange mass spectrometry with weighted ensemble molecular dynamics (MD). (2) We characterize the conformational heterogeneity of the ternary complex using Hamiltonian replica exchange simulations and small-angle X-ray scattering. (3) We assess the ubiquitination of the POI in the context of the full Cullin-RING Ligase, confirming experimental ubiquitinomics results. Differences in degradation efficiency can be explained by the proximity of lysine residues on the POI relative to ubiquitin.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Unraveling the role of temperature on the onset of ejecta formation at atomic scales

Shock wave interactions with perturbations on a free surface can lead to the inversion and growth of the perturbation and eventual ejection of a jet of material, referred to as an ejecta microjet, from shocked surfaces. Here, this study carries out large-scale molecular dynamics (MD) simulations using single-crystal Cu and Sn systems with a pre-existing groove to characterize the localized gradients in temperatures and pressures generated that render the microjet formation. MD simulations are carried out for three loading orientations (along the [001], [011], and [111] directions) and shock pressures ranging from 16 GPa to 100 GPa to understand the role of temperatures generated at the groove vertex on the formation of ejecta. The simulations suggest that the interaction of the shock wave with the groove results in a localized increase in temperatures, leading to localized softening at the groove vertex and the generation of an ejecta microjet. For Cu systems, the simulations suggest that the orientation effects on shock wave structures, velocities, and localized softening affect the ejecta formation at low pressures. Jetting is only observed when temperatures at the groove vertex are high enough to induce localized softening (close to melting temperature), and the jet velocity increases with shock pressure. In contrast, the loading orientation is rendered inconsequential for Sn systems due to the melting of the material during shock compression at the pressures chosen. The jet velocities are similar, regardless of crystal orientation in Sn systems.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Control of solvent adsorption in porous liquids through solvent-solvent interactions

Type 3 Porous Liquids (PLs) are a new class of materials that offer transformative potential for gas capture and utilization. These PLs are comprised of a bulky solvent with a suspended empty nanoporous material that enables selective and high-capacity capture from dilute and complex gas steams. The relationship between the nanoporous material, the structure of the solvent molecules, and the time dependence of solvent adsorption into the nanoporous material governs long term adsorption properties. Herein, molecular dynamics (MD) calculations evaluated the solvent dynamics of nine neat solvents with a broad range of chemical identities and experimental density measurements probed the time-dependent adsorption of these solvents into the ZIF-8 metal-organic frameworks (MOFs) to form PL dispersions over three weeks. Two of the nine dispersion, using glyceryl triacetate and 2’-hydroxyacetopheneone as solvents, presented sufficiently low solvent infiltration—characterized by solvent adsorption less than 40% of the ZIF-8 pore volume—to be viable as PLs. The experimental solvent adsorption data for ZIF-8 dispersed in the four aromatic solvents (acetophenone, methylbenzoate, 2-isopropylphenol, and 2’-hydroxyacetophenone) was fit to a kinetic model to quantify rate of solvent adsorption into ZIF-8. The rates of adsorption of these four similar-sized solvents demonstrated that solvent adsorption is slowed and limited by solvent clustering dynamics that are, in turn, driven by intermolecular hydrogen bonding. Ultimately, experimental solvent sorption measurements indicated that glyceryl triacetate and 2’-hydroxyacetophenone with ZIF-8 formed PLs with the most stable porosity and the greatest potential for gas capture. In conclusion, the combined experimental-computational materials exploration framework used here reveals a novel relationship between molecular scale solvent dynamics and macroscale solvent adsorption kinetics critical for the discovery and rapid evaluation of Type 3 PL compositions.

Hydrogen bonding↗