Search NASA⌕ Search

SEARCH · Search NASA

Results for “Work function”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

The identification of excited species in arc jet flow

Spectrographic work done at the Atmospheric Reentry Material and Structures Facility (arc jet) located at the Johnson Space Center has led to the identification of several excited molecular and atomic states. The excited molecular states identified are: first positive nitrogen system, second positive nitrogen system, the first negative nitrogen system, the gamma system for nitric oxide, and the 306.4 nm system of OH. Excited atoms identified were nitrogen, oxygen, hydrogen, silicon, copper, sodium, barium, potassium, and calcium. The latter five are considered contaminants. Excited molecular states of oxygen were not seen, suggesting full dissociation of oxygen molecules to oxygen atoms within the arc column and nozzle. Further, evidence exists that O(-) may be present since a background continuum is seen, and because of the existence of positive species (first negative system of N2(+)). Interpretation of spectrographic plates was enhanced by the use of a microdensitometer, and by the application of a second order least squares routine which determined wavelength as a function of plate location. Results of this work will ultimately improve models used in the calculation of heat transfer rates to the space shuttle and the aerobraking orbit transfer vehicles.

Willey, Ronald J.↗
DOE OSTI2024

Metagenomic clustering links specific metabolic functions to globally relevant ecosystems

ABSTRACT Metagenomic sequencing has advanced our understanding of biogeochemical processes by providing an unprecedented view into the microbial composition of different ecosystems. While the amount of metagenomic data has grown rapidly, simple-to-use methods to analyze and compare across studies have lagged behind. Thus, tools expressing the metabolic traits of a community are needed to broaden the utility of existing data. Gene abundance profiles are a relatively low-dimensional embedding of a metagenome’s functional potential and are, thus, tractable for comparison across many samples. Here, we compare the abundance of KEGG Ortholog Groups (KOs) from 6,539 metagenomes from the Joint Genome Institute’s Integrated Microbial Genomes and Metagenomes (JGI IMG/M) database. We find that samples cluster into terrestrial, aquatic, and anaerobic ecosystems with marker KOs reflecting adaptations to these environments. For instance, functional clusters were differentiated by the metabolism of antibiotics, photosynthesis, methanogenesis, and surprisingly GC content. Using this functional gene approach, we reveal the broad-scale patterns shaping microbial communities and demonstrate the utility of ortholog abundance profiles for representing a rapidly expanding body of metagenomic data. IMPORTANCE Metagenomics, or the sequencing of DNA from complex microbiomes, provides a view into the microbial composition of different environments. Metagenome databases were created to compile sequencing data across studies, but it remains challenging to compare and gain insight from these large data sets. Consequently, there is a need to develop accessible approaches to extract knowledge across metagenomes. The abundance of different orthologs (i.e., genes that perform a similar function across species) provides a simplified representation of a metagenome’s metabolic potential that can easily be compared with others. In this study, we cluster the ortholog abundance profiles of thousands of metagenomes from diverse environments and uncover the traits that distinguish them. This work provides a simple to use framework for functional comparison and advances our understanding of how the environment shapes microbial communities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Development of an ab initio learned model of electron deposition range in deuterium-tritium plasmas through time-dependent density functional theory calculations and machine learning

Accurate hydrodynamic modeling for laser-direct-drive (LDD) inertial-confinement-fusion (ICF) relies on precise calculations of the electron thermal conduction in all target materials. The nonlocal stopping range of electrons in ICF plasmas directly influences thermal conduction; yet, no first principles model exists for the electron mean free path in the conduction-zone regime. This work utilized time-dependent stochastic density-functional theory (TD-sDFT) to calculate the electron stopping power in deuterium-tritium (DT) plasmas at (ρ, T) conditions relevant to the conduction zone and the compressed shell in ICF. Using a combination of our TD-sDFT data and already established analytical models, we developed and trained an artificial neural network to create a global model for the nonlocal electron deposition range, λ E . We compared our machine-learning (ML) based model for λ E to the currently-used modified-Lee-More model in LDD radiation-hydrodynamic codes, such as lilac, and saw an overall decrease in the deposition range. To understand the effects of λ E on LDD ICF implosion dynamics, we implemented the ML-based model into lilac; specifically, we looked at designs consistent with a current experiment on the OMEGA laser and for a newly designed LDD-ICF target for the future OMEGA-Next facility. In both cases, we saw an overall drop in predicted ablation pressure, peak areal density, and neutron yield due to the reduced thermal conduction (smaller λ E ) in DT plasmas. Comparisons with the experiment on OMEGA are also made.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Green's Function Wind Turbine Induction Model That Incorporates Complex Inflow Conditions

ABSTRACT In this work, we develop a new analytical turbine induction model that can incorporate complex inflow conditions including cases where the wind velocity and temperature profiles can vary as functions of height. This induction model is derived from the linearized Navier–Stokes and leads to a second‐order ODE that can be solved using a Green's function formulation. The corresponding Green's function for several configurations are found including the infinite domain, semi‐infinite domain with ground plane, and a power law velocity inflow profile. The results of this approach are then compared with simulations of the turbine induction field using the AMR‐Wind CFD solver with a uniformly loaded actuator disk model. These comparisons show that the Green's function approach captures the centerline blockage, three‐dimensional blockage flow field, and streamwise velocity slow down, with very good agreement for lower thrust conditions and at larger distances away from rotor disk. The effects of shear on the turbine blockage were also compared using a power law inflow profile, and we show that this approach matches the CFD predictions for the cases considered.

17 WIND ENERGY↗
DOE OSTI2025

Joint Factorization of QCD and QED Radiation in Lepton-Hadron Scattering

The factorization theorem plays an important role in the analysis of high energy quantum chromodynamic (QCD) processes, separating the nonperturbative hadronic interaction into the universal parton distribution functions (PDFs) and fragmentation functions (FFs) and the process-dependent interactions into short distance perturbative calculations, with any interference power suppressed. With a virtual photon exchange, lepton-hadron deep inelastic scattering (DIS) provides an electromagnetic hard probe for the partonic structure of colliding hadrons and has played an important role in the development of QCD factorization. However, the collision induced QED radiation can change the momentum of the exchanged but unobserved virtual photon, making the photon-hadron frame, where the factorization formalism for DIS and semi-inclusive DIS (SIDIS) was derived, ill defined. A new analogous factorization approach has been introduced to separate the leading power process-independent QED radiative contributions to the single photon exchange by introducing lepton distribution functions (LDFs) and lepton fragmentation functions (LFFs), while process-dependent effects are perturbatively calculated with large logarithms removed [J. High Energ. Phys. 2021, 157 (2021)]. These LDFs and LFFs are considered global, as they appear in many different interactions, such as e+e-, DIS and SIDIS, so data from experiments can be used to fit and describe these functions across a wide range of lepton scattering. In this work, I will apply this new hybrid factorization approach to lepton-hadron DIS and SIDIS. For DIS, I derive the NLO short distance perturbative contribution to the cross section and demonstrate the effects the QED radiation has on the cross section using this approach using the CTEQ parameterization for the QCD functions. As part of the SIDIS analysis, I study the cross-section in two different kinematic regions: (1) the scattered lepton and observed hadron are not near back-to-back, and (2) they are close to back-to-back, where collinear QCD factorization works for (1) and TMD QCD factorization for (2) while collinear QED factorization works for both. As part of this work, I show the effects on the SIDIS cross section using fixed order calculations for the unpolarized structure function by first showing the effect of the radiative corrections on the main kinematic variables, especially how the internal transverse momentum is significantly correlated to the external angular dependence, and then the unpolarized structure function (or cross section) with matching between the descriptions for low and high transverse momentum. This work will impact the calculations for predictions for data from COMPASS and various Jefferson Lab experiments.

Cammarota, Justin [Univ. of Kentucky, Lexington, K↗
DOE OSTI2026

The origin of metallic conductivity in Pt 3 O 4 : a first principles study

The platinum oxide Pt 3 O 4 exhibits metallic conductivity even though it contains square-planar PtO4 units, which in related oxides such as PtO are usually associated with insulating behavior. To identify the electronic origin of this anomalous metallicity, we performed a comprehensive first-principles study using the PBE and r 2 SCAN functionals together with Hubbard U corrections and spin-orbit coupling (SOC). Structural benchmarks show that r 2 SCAN with SOC and a moderate U value (<4 eV) reproduces the experimental lattice constants and formation enthalpy, whereas larger U values (~8 eV) destabilize the cubic structure. Across all functionals and U values considered in this work, Pt 3 O 4 remains metallic. Analyses of the projected density of states, band structures, charge-density isosurfaces, and bonding characteristics demonstrate that the dominant contribution to the metallic character originates from delocalized Pt–O–Pt hybridized antibonding states at the Fermi level. Direct Pt–Pt interactions are present but contribute less strongly to the conductivity. Bader charge analysis reveals only weak Pt charge disproportionation, consistent with mixed Pt II /Pt III character, and a small charge-transfer energy that prevents localization of the Pt 5d electrons even at elevated U. In contrast, PtO develops a Mott or charge-transfer gap under modest U despite having the same PtO 4 coordination environment. These findings demonstrate that persistent Pt–O–Pt covalency is the primary driver of metallicity in Pt 3 O 4 and support the view that this phase can remain conductive under oxygen reduction and oxygen evolution reaction conditions in fuel cell and electrolyzer environments.

36 MATERIALS SCIENCE↗
NASA NTRS2018

Human-Automation Allocations for Current Robotic Space Operations

Within the Human Research Program, one risk delineates the uncertainty surrounding crew working with automation and robotics in spaceflight. The Risk of Inadequate Design of Human and Automation/Robotic Integration (HARI) is concerned with the detrimental effects on crew performance due to ineffective user interfaces, system designs and/or functional task allocation, potentially compromising mission success and safety. Risk arises because we have limited experience with complex automation and robotics. One key gap within HARI, is the gap related to functional allocation. The gap states: We need to evaluate, develop, and validate methods and guidelines for identifying human-automation/robot task information needs, function allocation, and team composition for future long duration, long distance space missions. Allocations determine the human-system performance as it identifies the functions and performance levels required by the automation/robotic system, and in turn, what work the crew is expected to perform and the necessary human performance requirements. Allocations must take into account each of the human, automation, and robotic systems capabilities and limitations. Some functions may be intuitively assigned to the human versus the robot, but to optimize efficiency and effectiveness, purposeful role assignments will be required. The role of automation and robotics will significantly change in future exploration missions, particularly as crew becomes more autonomous from ground controllers. Thus, we must understand the suitability of existing function allocation methods within NASA as well as the existing allocations established by the few robotic systems that are operational in spaceflight. In order to evaluate future methods of robotic allocations, we must first benchmark the allocations and allocation methods that have been used. We will present 1) documentation of human-automation-robotic allocations in existing, operational spaceflight systems; and 2) To gather existing lessons learned and best practices in these role assignments, from spaceflight operational experience of crew and ground teams that may be used to guide development for future systems. NASA and other space agencies have operational spaceflight experience with two key Human-Automation-Robotic (HAR) systems: heavy lift robotic arms and planetary robotic explorers. Additionally, NASA has invested in high-fidelity rover systems that can carry crew, building beyond Apollo's lunar rover. The heavy lift robotic arms reviewed are: Space Station Remote Manipulator System (SSRMS), Japanese Remote Manipulator System (JEMRMS), and the European Robotic Arm (ERA, designed but not deployed in space). The robotic rover systems reviewed are: Mars Exploration Rovers, Mars Science Laboratory rover, and the high-fidelity K10 rovers. Much of the design and operational feedback for these systems have been communicated to flight controllers and robotic design teams. As part of the mitigating the HARI risk for future human spaceflight operations, we must document function allocations between robots and humans that have worked well in practice.

robotic allocation↗
DOE OSTI2025

Convergence analysis for a nonlocal gradient descent method via directional Gaussian smoothing

We analyze the convergence of a nonlocal gradient descent method for minimizing a class of high-dimensional non-convex functions, where a directional Gaussian smoothing (DGS) is proposed to define the nonlocal gradient (also referred to as the DGS gradient). The method was first proposed in [Zhang et al., Enabling long-range exploration in minimization of multimodal functions, UAI 2021], in which multiple numerical experiments showed that replacing the traditional local gradient with the DGS gradient can help the optimizers escape local minima more easily and significantly improve their performance. However, a rigorous theory for the efficiency of the method on nonconvex landscape is lacking. In this work, we investigate the scenario where the objective function is composed of a convex function, perturbed by deterministic oscillating noise. We provide a convergence theory under which the iterates exponentially converge to a tightened neighborhood of the solution, whose size is characterized by the noise wavelength. Here, we also establish a correlation between the optimal values of the Gaussian smoothing radius and the noise wavelength, thus justifying the advantage of using moderate or large smoothing radii with the method. Furthermore, if the noise level decays to zero when approaching the global minimum, we prove that DGS-based optimization converges to the exact global minimum with linear rates, similarly to standard gradient-based methods in optimizing convex functions. Several numerical experiments are provided to confirm our theory and illustrate the superiority of the approach over those based on the local gradient.

Tran, Hoang [Oak Ridge National Laboratory (ORNL),↗
NASA NTRS1999

Electroencephalographic imaging of higher brain function

High temporal resolution is necessary to resolve the rapidly changing patterns of brain activity that underlie mental function. Electroencephalography (EEG) provides temporal resolution in the millisecond range. However, traditional EEG technology and practice provide insufficient spatial detail to identify relationships between brain electrical events and structures and functions visualized by magnetic resonance imaging or positron emission tomography. Recent advances help to overcome this problem by recording EEGs from more electrodes, by registering EEG data with anatomical images, and by correcting the distortion caused by volume conduction of EEG signals through the skull and scalp. In addition, statistical measurements of sub-second interdependences between EEG time-series recorded from different locations can help to generate hypotheses about the instantaneous functional networks that form between different cortical regions during perception, thought and action. Example applications are presented from studies of language, attention and working memory. Along with its unique ability to monitor brain function as people perform everyday activities in the real world, these advances make modern EEG an invaluable complement to other functional neuroimaging modalities.

Non-NASA Center↗
DOE OSTI2026

Engineering an Extremely Hybrid PKS for Adipic Acid Production

Polyketide synthases (PKSs) are modular enzymes with exceptional potential as biocatalysts for producing non-native compounds. Here, we report the first PKS-based pathway for adipic acid (AA), an industrial monomer for nylon production, by engineering one of the most extensively hybridized PKS systems to date. Using a retrobiosynthetic approach, we identified EtnB, a succinyl-CoA-loading module that uniquely retains the terminal carboxyl group, enabling access to dicarboxylic polyketide products, rarely produced by canonical PKSs. EtnB was coupled to a fully reducing extension module through an engineered communication linker, which improved ACP–KS interactions, enhanced titers, and demonstrated selective succinyl-CoA loading in vivo . This construct integrates genes from five organisms─with domains from seven PKS modules joined across six non-natural junctions─and functions in both Escherichia coli and Pseudomonas putida . Additional engineering that included AT domain exchanges, optimization of extender unit supply, and host strain metabolic rewiring further increased AA production. Together, this work demonstrates that highly chimeric PKSs can be rendered functional through rational design, expands the PKS toolkit with a carboxyl-retaining loading module, and establishes a versatile platform for engineering diacids and other noncanonical products through PKS pathways.

biomanufacturing↗
DOE OSTI2025

Machine Learning‐Guided Discovery of High‐Entropy Perovskite Oxide Electrocatalysts via Oxygen Vacancy Engineering

Abstract High‐entropy perovskite oxides (HEPOs) have recently emerged as multifunctional catalysts. However, the HEPOs’ structural and compositional complexity hinders the easy and accurate extrapolation of activity indicators, which are essential for establishing structure‐property correlations. Here, OxiGraphX, is introduced as a novel graph neural network (GNN) model designed to capture the complex relationships among structure, composition, and atomic chemical environments for accurate prediction of oxygen vacancy formation energies (OVFEs) in HEPOs. By integrating machine learning (ML), density functional theory (DFT), and experimental validation, this work demonstrates an efficient framework for rapidly and accurately screening HEPO electrocatalysts for oxygen evolution reaction (OER). The OxiGraphX predicts OVFEs with a precision exceeding existing data, enabling the identification of compositions of higher oxygen vacancy content (OVC) and, thus, higher catalytic activity. Furthermore, the model explores latent spaces that translate effectively into experimental domains, bridging computational predictions with real‐world applications. This approach accelerates the discovery of high‐performance HEPO catalysts while providing deeper insights into their catalytic mechanisms.

Chemistry↗
DOE OSTI2024

Exact enforcement of temporal continuity in sequential physics-informed neural networks

The use of deep learning methods in scientific computing represents a potential paradigm shift in engineering problem solving. One of the most prominent developments is Physics-Informed Neural Networks (PINNs), in which neural networks are trained to satisfy partial differential equations (PDEs). While this method shows promise, the standard version has been shown to struggle in accurately predicting the dynamic behavior of time-dependent problems. To address this challenge, methods have been proposed that decompose the time domain into multiple segments, employing a distinct neural network in each segment and directly incorporating continuity between them in the loss function of the minimization problem. In this work we introduce a method to exactly enforce continuity between successive time segments via a solution ansatz. This hard constrained sequential PINN (HCS-PINN) method is simple to implement and eliminates the need for any loss terms associated with temporal continuity. The method is tested for a number of benchmark problems involving both linear and non-linear PDEs. Examples include various first order time dependent problems in which traditional PINNs struggle, namely advection, Allen–Cahn, and Korteweg–de Vries equations. Furthermore, second and third order time-dependent problems are demonstrated via wave and Jerky dynamics examples, respectively. Notably, the Jerky dynamics problem is chaotic, making the problem especially sensitive to temporal accuracy. Finally, the numerical experiments conducted with the proposed method demonstrated superior convergence and accuracy over both traditional PINNs and the soft-constrained counterparts.

42 ENGINEERING↗
DOE OSTI2025

Resolving Lonsdaleite's decade-long controversy: Atomistic insights into a metastable diamond polymorph

Lonsdaleite, a theoretically proposed hexagonal diamond polymorph, has remained at the center of a five-decade scientific controversy since its 1967 identification. While some studies claim it exhibits superior hardness through compression-induced structural changes, others contend it is merely a stacking-faulted cubic diamond. Meteoritic samples and synthetic preparations have yielded conflicting evidence, with even advanced characterisation techniques like XRD and TEM failing to provide definitive proof. In this work, we employ first-principles density functional theory (DFT) and molecular dynamics (MD) simulations to generate unambiguous theoretical fingerprints through XRD, Raman, and SAED patterns that distinguish true Lonsdaleite from cubic diamond and its defective variants. Our atomistic approach quantifies the thermodynamic metastability of Lonsdaleite under realistic pressure-temperature conditions, reveals distinct spectral signatures through simulated Raman and resolves the structural ambiguity through generalised stacking fault energy analysis. By establishing clear criteria for definitive identification, this study provides long-awaited clarity to the Lonsdaleite debate while offering a robust computational framework for characterising metastable carbon phases in meteoritic, synthetic and industrial materials.

DFT↗
DOE OSTI2025

Data-Efficient Strategies for Probabilistic Voltage Envelopes under Network Contingencies

This work presents an efficient data-driven method to construct probabilistic voltage envelopes (PVE) using power flow learning in grids with network contingencies. First, a network-aware Gaussian process (GP) termed Vertex-Degree Kernel (VDK-GP), developed in prior work, is used to estimate voltage–power functions for a few network configurations. The paper introduces a novel multi-task vertex degree kernel (MT-VDK) that amalgamates the learned VDK-GPs to determine power flows for unseen networks, with a significant reduction in the computational complexity and hyperparameter requirements compared to alternate approaches. Simulations on the IEEE 30-Bus network demonstrate the retention and transfer of power flow knowledge in both N-1 and N-2 contingency scenarios. The MT-VDK-GP approach achieves over 50 % reduction in mean prediction error for novel N-1 contingency network configurations in low training data regimes (50–250 samples) over VDK-GP. Additionally, MT-VDK-GP outperforms a hyper-parameter based transfer learning approach in over 75 % of N-2 contingency network structures, even without historical N-2 outage data. Furthermore, the proposed method demonstrates the ability to achieve PVEs using sixteen times fewer power flow solutions compared to Monte-Carlo sampling-based methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Improved Protein Semi-Synthesis Enables Biophysical Studies of Thioamide Destabilization of β-Sheet Interactions

Abstract Thioamides are natural post-translational modifications of the peptide backbone and can be introduced synthetically to probe protein folding or functionalize peptides for translational applications. In this work, we demonstrate that thioamide-containing peptides with C-terminal thioesters can be efficiently generated using Knorr pyrazole activation and used in subsequent native chemical ligation reactions to generate thioamide-containing proteins. We compare this method to acyl azide activation and find that both routes provide similar yields. We also investigate ultrasound-mediated desulfurization of the ligation site cysteine for potential advantages over chemical radical initiators. Scaling up our syntheses allows us to study thioamide perturbations to the β-sheet region of the B1 domain of protein G (GB1) as well as β-strand interactions in amyloid fibrils of the Parkinson’s disease protein α-synuclein. In both contexts, we observe dramatic destabilization of the β-sheet networks, manifested in decreased GB1 thermal stability and altered folding and slowed aggregation of α-synuclein. These findings illustrate the impact that a single atom substitution can have on cooperative hydrogen-bonding networks and prompt future study of both systems.

Yanagawa, Evan S. K. [University of Pennsylvania ,↗
DOE OSTI2025

Enhancing Dihydrogen Interaction of a Zirconium Metal–Organic Framework by Metal Doping

Advancing efficient hydrogen storage technologies is essential for enabling a sustainable energy future, especially in onboard applications. While hydrogen offers high gravimetric energy density and zero-emission combustion, its low volumetric energy density presents significant storage challenges. Metal-organic frameworks (MOFs), well known for their tunable porosity and high surface areas, have emerged as promising candidates for adsorption-based hydrogen storage. This study investigated a chemically robust zirconium-based MOF, MOF-808, as a representative platform for hydrogen storage enhancement through metal ion doping. Various divalent metal ions were introduced into MOF-808 via one-pot synthesis or post-synthetic modification (PSM) to evaluate their effects on metal doping efficiency, framework stability, and hydrogen adsorption performance. Our findings demonstrate that metal doping enhanced the hydrogen binding affinity of MOF-808 while preserving its structural integrity and excellent stability. A Mg-doped MOF, MOF-808@Mg 2:1, showed a 59% increase in hydrogen uptake, and a Cu-doped MOF, MOF-808-ZrCu, exhibited a 33% increase in isosteric heat of adsorption for H 2 compared to the pristine MOF-808 activated at the same temperature. This work highlights the potential of metal-functionalized stable MOFs for practical hydrogen storage applications. Furthermore, these materials are also being studied for their potential to enhance CO 2 adsorption.

Adsorption↗
DOE OSTI2025

Molecular Modeling and Molecular Dynamics Simulation of a Packed and Intact Bacterial Microcompartment

Bacterial microcompartments (BMCs) are protein-bound organelles found in some bacteria which encapsulate enzymes for enhanced catalytic activity. These compartments spatially sequester enzymes within semipermeable shell proteins and are packed full of enzyme cargoes and metabolites as they fulfill their function. Coupling together recent SAXS and proteomics work, it is possible to develop molecular models for these microcompartments and interrogate enzyme and metabolite dynamics within. Our primary goal of this study is to quantify the permeability of metabolite glyceraldehyde-3-phosphate (G3P) and dihydroxyacetone phosphate (DHAP) across the BMC shell through classical molecular dynamics simulation. The Haliangium ochraceum model of BMC shell (PDB: 6MZX) was used to model an intact BMC of approximately 10 million atoms. Working at this scale presented its own challenges in managing large data sets, with multiple challenges and hardware advances discussed that facilitated this work. Over approximately 750 ns of aggregate simulation, we see multiple permeation events for these metabolites that were added at high concentration through the pores present within BMC shell tiles. When compared to independent permeability estimates for the same metabolites determined through replica exchange umbrella sampling simulations, the permeabilities varied by approximately 3 orders of magnitude. Regardless, the permeability coefficients for both G3P and DHAP are highly similar and very high, such that only very small concentration gradients can be maintained across the BMC shell between the cytosol and BMC interior. The large simulation systems also facilitated comparisons for molecular diffusivity in the crowded environment within the BMC shell. By our estimates, the viscosity within a packed BMC shell is at least 10-fold higher than it would be in neat solution and is the real driver for varying permeability estimates we obtained through simulation. These findings will be used as design inputs for future bioengineering efforts to make products from BMCs, highlighting how permeable BMC shells can be.

Diffusion↗
DOE OSTI2025

Computing Reaction Kinetics with MC-PDFT–OPESf: Combining Multireference Electronic Structure Theory and Enhanced Sampling

Accurate rate constants are crucial for understanding and optimizing catalytic reactions mediated by enzymes, metalloproteins, and heterogeneous catalysts. These systems frequently present a dual computational challenge. Multiconfigurational reaction sites require multireference techniques for the accurate treatment of the electronic structure, and high activation barriers prevent efficient sampling of unbiased reactive transitions. In this work, we combine multiconfiguration pair-density functional theory (MC-PDFT) as an accurate and efficient multireference electronic structure method with on-the-fly probability-enhanced sampling flooding (OPESf) as an enhanced sampling method capable of accelerating reactive transitions. We demonstrate the approach on the Diels–Alder [4+2] cycloaddition between cis-butadiene and ethene as a reaction characterized by a large activation barrier and multireference character. MC-PDFT–OPESf provides reaction rates in agreement with experiments at a fraction of the computational cost required by conventional unbiased ab initio calculations. Here, we propose MC-PDFT–OPESf as an efficient approach for computing kinetics in strongly correlated molecular systems.

Chemical calculations↗