Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complexity theory”

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

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

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

Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces

Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems

DLVO↗
DOE OSTI2026

B -Alkyl-borabicyclo[3.3.1]nonane Reagents Promote Closed-Shell Nickel-Catalyzed Alkylarylation Toward Encoded Cyclooctene Monomers

Access to new tailored monomers is essential to explore unprecedented polymer structures and modulate their properties. In previous work, we disclosed an iteroselective diarylation of 1,5-cyclooctadiene as an attractive platform for preparing diverse cyclooctene monomers; however, it was limited to C(sp 2 ) coupling partners. Established difunctionalization methods with C(sp 3 ) fragments rely on carbon-based radical formation, which is incompatible with 1,5-cyclooctadiene. This work describes a two-electron redox manifold for the nickel-catalyzed alkylarylation of 1,5-cyclooctadiene and discloses B -alkyl-borabicyclo[3.3.1]nonane (alkyl-9-BBN) as an effective transmetalating reagent for maintaining a polar reaction mechanism. The method provides 5,6-alkylarylated cyclooctenes suitable for ring-opening metathesis polymerization to obtain new materials. The properties of these polymers are benchmarked and fine-tuned by variation of coupling partners in the nickel catalysis. Density functional theory calculations revealed that destabilization of the pretransmetalation complex promotes the reactivity of alkyl-9-BBN in transmetalation compared to alkylboronic esters.

alkyls↗
DOE OSTI2025

Computational alchemy clarifies origins of alloy strengthening

Solid solution strengthening (SSS) is widely used to enhance mechanical properties of metals. Originally developed for dilute alloys, classical SSS theories are presently challenged by the rise of complex concentrated alloys (CCA) with nearly equiatomic compositions. Here, we propose and develop a method of “computational alchemy” in which interatomic interactions are modified to systematically vary two key physical parameters defining SSS - atomic size misfit and elastic stiffness misfit - over a maximally wide range of two misfits. The resulting alchemical alloys are subjected to massive (~10 8 atoms) molecular dynamics (MD) simulations reproducing full complexity of plastic strength response. At variance with prevailing views, stiffness misfit is observed to contribute to SSS on par if not more than size misfit. Furthermore, depending on exactly how two misfits are combined, they result in synergistic (amplification) or antagonistic (compensation) effect on alloy strengthening. Unlike real CCAs in which each component element comes with its own specific size and stiffness, our alchemical model alloys span the space of two misfits continuously revealing trends in alloy strengthening unrecognized so far. Our study demonstrates unique value of intentionally unrealistic models for gaining deep physical insights into material behaviors that are difficult to reveal otherwise.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Quantifying uncertainty in machine learning for nuclear binding energy

Techniques from artificial intelligence and machine learning are increasingly employed in nuclear theory; however, the uncertainties that arise from the complex parameter manifold encoded by the neural networks are often overlooked. Epistemic uncertainties arising from training the same network multiple times for an ensemble of initial weight sets offer a first insight into the confidence of machine learning predictions, but they often come with a high computational cost. Instead, we apply a single-model uncertainty quantification method called Δ-UQ that gives epistemic uncertainties with one-time training. Here, we demonstrate our approach on a two-feature model of nuclear binding energies per nucleon with proton and neutron number pairs as inputs. We show that Δ-UQ can produce reliable and self-consistent epistemic uncertainty estimates and can be used to assess the degree of confidence in predictions made with deep neural networks.

Huang, Mengyao [Lawrence Livermore National Labora↗
DOE OSTI2024

Approach for Inferring Full-Scope Human Reliability Data Based on Simplified Simulator Data

This paper proposes a method for inferring full-scope human reliability data based on the Simplified Human Error Experimental Program (SHEEP) data. It mainly focuses on the human errors observed when using simulators with different complexity levels. In the proposed method, the manner in which human error probabilities (HEPs) change as a result of increasing simulator complexity and how simulator complexity levels are quantified represent key information for inferring full-scope data. In the present study, SHEEP error data pertaining to actual professional operators using Rancor Microworld (Rancor) (i.e., a more simplified simulator) and Compact Nuclear Simulator (CNS) (i.e., a less simplified simulator) were compared with the HuREX error data. An approach to quantifying simulator complexity levels was then proposed based on information theory and acquired eye-tracker data.

99 - GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Influence of Ga doping on magnetic properties, magnetocaloric effect, and electronic structure of pseudobinary GdZn 1 - x Ga x ( x = 0 - 0.1 )

Here, we explore the impact of introducing IIIA-group element Ga in place of IIB-group element Zn in binary intermetallic GdZn on its magnetic and magnetocaloric properties, as well as explicate the modified electronic band structure of the compound. The magnetic transition temperature of the compound decreases with the increase of Ga concentration in GdZn 1-x ⁢Ga x ⁢(x = 0–0.1) while the crystal structure (CsCl prototype) and lattice parameters remain unchanged. Our detailed analysis of magnetization and magnetocaloric data conclusively proves that long-ranged magnetic ordering exists in the sample, despite the magnetic interaction considerably weakening with the increase of Ga. The experimental data are rationalized using both a theoretical machine-learning model and first-principle density functional theory. The electronic band structure of GdZn is manifested with some unusual complex features which gradually diminish with Ga doping and the conventional sinusoidal feature of Ruderman-Kittel-Kasuya-Yosida-type interactions also disappears. A mean-field theory model is developed and can successfully describe the overall magnetocaloric behavior of the GdZn 1-x⁢ Ga x series of samples.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Can Dithiolate Ligands Report Electronic Communication in Transuranium Complexes?

The actinide elements, and in particular the transuranic elements, are some of the least studied elements on the periodic table. The degree of covalent bonding between the actinide elements and ligands in is not well understood. Dithiolene/dithiolate ligands provided a revolution in understanding the bonding in transition metals through the geometric differences the ligands exhibit when interacting with electron deficient or abundant metal centers. To test if dithiolene/dithiolate ligands can be utilized to provide an analogous understanding of transuranic element bonding, density functional theory calculations have been performed on U, Np, and Pu complexes of the formula Cp 2 AnS 2 C 2 H 2 . These calculations show that the dithiolate maintains its non-innocent redox active nature when bonding with actinides. In stark contrast with the transition metals, the direction of electron donation is reversed, in the actinide series the direction of electron flow in high fold cases is from the metal center to the ligand. In conclusion, this suggests the communication is between the π* system of the ligand interacting with the f orbitals of the metal center.

Prange, Micah P. [Pacific Northwest National Labor↗
DOE OSTI2025

Advances and Challenges of SCAN and r 2 SCAN Density Functionals in Transition‐Metal Compounds

ABSTRACT Transition‐metal compounds (TMCs) with open‐shell d ‐electrons are characterized by a complex interplay of lattice, charge, orbital, and spin degrees of freedom, giving rise to various fascinating applications. Often exhibiting exotic properties, these compounds are commonly classified as correlated systems due to strong inter‐electronic interactions called Hubbard U . This inherent complexity presents significant challenges to Kohn‐Sham density functional theory (KS‐DFT), the most widely used electronic structure method in condensed matter physics and materials science. While KS‐DFT is, in principle, exact for the ground‐state total energy, its exchange‐correlation energy must be approximated in practice. The mean‐field nature of KS implementations, combined with the limitations of current exchange‐correlation density functional approximations, has led to the perception that DFT is inadequate for correlated systems, particularly TMCs. Consequently, a common workaround involves augmenting DFT with an on‐site Hubbard‐like U correction. In recent years, the strongly constrained and appropriately normed (SCAN) density functional, along with its refined variant r 2 SCAN, has achieved remarkable progress in accurately describing the structural, energetic, electronic, magnetic, and vibrational properties of TMCs, challenging the traditional perception of DFT's limitations. This review explores the design principles of SCAN and r 2 SCAN, highlights their key advancements in studying TMCs, explains the mechanisms driving these improvements, and addresses the remaining challenges in this evolving field.

Zhang, Yubo [Minjiang Collaborative Center for The↗
DOE OSTI2024

Can we mitigate tropical cyclone formation using aerosols? A review of cyclogenesis and aerosol effects as a theoretical basis

Tropical cyclogenesis is a complex phenomenon, involving processes across multiple scales, with various theories proposed to explain its mechanisms. This paper provides an up-to-date summary of the established physical mechanisms and controlling factors of cyclogenesis. Understanding the nature of tropical cyclone formation is necessary to improve tropical cyclone forecasts and inform tropical cyclone mitigation strategies targeting early-stage intervention. Despite the disparities among theories, the critical role of persistent, intensified convection in tropical cyclone initiation is universally acknowledged. Therefore, our hypothesis is that interrupting convective development using aerosol injection is a potential strategy for mitigating tropical cyclone formation. Targetted injection of coarse and fine aerosols at the formation stage are two proposed potential approaches for mitigating tropical cyclones. These methods have completely different microphysical effects but each may ultimately interfere with the formation of the nascent cyclone. In conclusion, these potential interventions open new areas of further research, which are important and necessary to establish a solid scientific foundation for cyclogenesis mitigation.

Tropical cyclone formation↗
DOE OSTI2025

Influence of initial conditions on data-driven model identification and information entropy for ideal mhd problems

Data-driven methods of model identification are able to discern governing dynamics of a system from data. Such methods are well suited to help us learn about systems with unpredictable evolution or systems with ambiguous governing dynamics given our current understanding. Many plasma problems of interest fall into these categories as there are a wide range of models that exist, however each model is only useful in a certain regime and often limited by computational complexity. To ensure data-driven methods align with theory, they must be consistent and predictable when acting on data whose governing dynamics are known. Weak Sparse Identification of Nonlinear Dynamics (WSINDy) is a recently developed data-driven method that has shown promise in learning governing dynamics from data with high noise levels [1]. This work examines how WSINDy acts on ideal MHD test problems as the initial conditions are varied and specifies limiting requirements for successful equation identification. Furthermore, it is hard to recover the governing dynamics from data that emphasize a single dominant behavior. In these low information cases, Shannon information entropy is able to pick up on the redundancies in the data that affect recoverability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Atomic-scale vacancy engineering unlocks basal-plane catalytic activity in metallic WSe 2 for reversible oxygen electrocatalysis

Two-dimensional metallic transition metal dichalcogenides offer high electrical conductivity and large surface areas for electrocatalysis, yet their inherent basal planes are catalytically inert. Here, we present an atomic-scale vacancy engineering strategy to activate the basal surfaces of metallic WSe 2 for reversible oxygen electrocatalysis. This approach, based on intentionally designed substitutional metal doping, promotes the spontaneous formation of selenium vacancies while preserving the metallic 1 T′ phase, thereby creating highly reactive and oxygen-affinitive sites. Density functional theory calculations reveal that these vacancy-mediated metal complexes dramatically lower the energy barriers for initial oxygen adsorption, enabling dissociative oxygen adsorption. Operando and ex-situ spectroscopic analyses confirm that vacancy-mediated metal complexes transform into dynamic Se/W-oxide intermediates under operating conditions. Se/W-oxides on the surface experimentally and theoretically prove electrocatalytic activity and reversibility. Applying this strategy in lithium–oxygen batteries, the basal-plane activated WSe 2 shows high discharge capacities (9868 mA h g −1 , corresponding to 3947 mA h g$^{-1}_{cathode}$), impressive cycle retention over 550 cycles at 1000 mA h g −1 , and outstanding rate–capability over a wide current–density range (100–3000 mA g −1 ) during 256 cycles.

2D materials↗
DOE OSTI2025

Impacts and emerging research opportunities in Vehicle-Grid Integration for transportation: A review

This review provides a comprehensive examination of Vehicle-Grid Integration (VGI) technologies and their impacts on transportation systems, with a particular emphasis on the transportation-energy nexus. It systematically explores how VGI affects key transportation applications such as charging infrastructure planning, electric vehicle (EV) routing, smart charging coordination, shared mobility, and dynamic pricing. By synthesizing recent literature from both transportation and energy systems perspectives, this study highlights how advanced methodologies, such as reinforcement learning, game theory, and optimization techniques, are used to model the complex interactions between EVs, mobility patterns, and distributed energy systems. Furthermore, the review also identifies critical challenges, including behavioral factors, data limitations, and system scalability. Drawing on these insights, the paper outlines emerging research opportunities to support the design of integrated, resilient, and user-centric VGI solutions that advance sustainable mobility and energy system efficiency.

Charging coordination↗
DOE OSTI2024

Density Functional Tight-Binding Models for Band Structures of Transition-Metal Alloys and Surfaces across the d -Block

First-principles electronic structure simulations are an invaluable tool for understanding chemical bonding and reactions. While machine-learning models such as interatomic potentials significantly accelerate the exploration of potential energy surfaces, electronic structure information is generally lost. Particularly in the field of heterogeneous catalysis, simulated electron band structures provide fundamental insights into catalytic reactivity. This ab initio knowledge is preserved in semiempirical methods such as density functional tight binding (DFTB), which extend the accessible computational length and time scales beyond first-principles approaches. In this paper here we present Shell-Optimized Atomic Confinement (SOAC) DFTB electronic-part-only parametrizations for bulk and surface band structures of all d-block transition metals that enable efficient predictions of electronic descriptors for large structures or high-throughput studies on complex systems outside the computational reach of density functional theory.

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

Mechanisms and Energetics of CO 2 Chemisorption in Choline-Based Eutectic Solvents

Understanding the influences of proton sharing and hydrogen transfer in amine and hydroxyl functionalized eutectic solvents is important for tuning CO 2 binding as they relate to the chemisorption capacity and energetics of solvent regeneration in carbon capture. In this study, we examine how the CO 2 binding mechanism and energetics vary in eutectic solvents composed of choline-based hydrogen-bond acceptors (HBAs) with amine and hydroxyl moieties and ethylene glycol (EG), propylene glycol (PG), 1,2-butanediol (BD), and monoethanolamine (MEA)-type hydrogen-bond donors (HBDs), through absorption–desorption experiments and spectral analysis, supported with density functional theory (DFT) calculations. While the HBD-CO 2 complex is the major and the thermodynamically more favored product of absorption in all of the examined eutectic solvents, this complexation is only possible proceeding a proton transfer from HBD to the HBA in diol-based systems (EG, PG, BD). The interplay between HBA-CO 2 and HBD-CO 2 speciation is governed by the HBD concentration, basicity, and steric effects, with increased HBD-CO 2 interactions in the presence of smaller HBDs (size: BD > PG > EG). On the other hand, in MEA-based systems, CO 2 readily binds to the amine followed by a proton transfer from the amine to the HBA as confirmed by DFT calculations. In this way, the MEA-CO 2 chemisorption capacities exceeded that of the conventional amine (up to 0.93 mol of CO 2 per mol of MEA versus 0.5). Furthermore, the increase in CO 2 capacities and enhanced binding in eutectics including MEA or the more basic HBAs also resulted in regeneration temperatures of 70 to 60 °C with calculated reaction energies of −82 to −57 kJ/mol, compared to 50 °C for others with weaker binding, following CO 2 absorption at 25 °C.

absorption↗
DOE OSTI2025

Shell Thickness and Heterogeneity Dependence of Triplet Energy Transfer between Core–Shell Quantum Dots and Adsorbed Molecules

Quantum dot (QD)-sensitized triplet energy transfer (TET) has found promising applications in photon upconversion and photocatalysis. However, the underlying mechanism of TET in the QD-acceptor complex remains unclear despite the well-developed TET theory for the molecular donor–acceptor systems. Herein, the coupling strength of TET from CdSe/CdS core–shell QDs to 9-anthracene carboxylic acid (ACA) was studied by measuring the TET rate as a function of shell thickness with time-resolved photoluminescence. The change of TET-coupling strength with increasing shell thickness was further compared to those of electron and hole transfers from QDs so that we could test whether QD-sensitized TET is mediated by the charge transfer virtual state and can be considered as simultaneous electron and hole transfers as in molecular donor–acceptor systems. The measured coupling strength of TET from the CdSe/CdS QD decreases exponentially with the CdS shell thickness r: |V|(r) = |V|(0)e –βr , with an exponential decay factor β of 0.19 Å –1 , which is smaller than the sum of the measured decay factors for electron transfer to methyl viologen (0.18 Å –1 ) and hole transfer to phenothiazine (0.29 Å –1 ) from the same QD. This inconsistency is explained by the broadening of QD shell thicknesses in the distance dependence study, which significantly modifies the TET-coupling strength and driving force, resulting in a shallower distance dependence of the TET rate constants. This study sheds light on the fundamental mechanisms of QD-sensitized TET reactions.

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

Multistage nucleation pathway in LiF molten salt mirrors the crystal–melt interface structure

Despite over a century of studies, fundamental questions remain about the processes governing crystal nucleation from melts or solutions. Research over the past three decades has presented mounting evidence for kinetic pathways of crystal nucleation that are more complex than envisioned by the simplest forms of classical theory. Such observations have been presented for colloidal and elemental systems with covalent and metallic bonding. Despite the technological and geochemical importance of molten salts, similar studies for these ionically bonded systems are currently lacking. Here we develop a machine learning interatomic potential for a model ionic system: LiF. The potential features quantum-level accuracy for both liquid and multiple solid polymorphs over wide temperature and pressure ranges and accurately reproduces experimentally measured properties. Thanks to the efficiency of the potential, which enables microsecond-scale molecular dynamics simulations, induction times for nucleation of LiF solids from their melts are computed over a range of undercoolings. With the aid of a set of robust local order parameters established here, the simulations reveal that homogeneous crystal nucleation in undercooled melts preferentially initiates from liquid regions showing slow dynamics and high bond orientational order simultaneously, and the second-shell order of both precritical nuclei and the surface of postcritical nuclei is dominated by hexagonal close packing and body-centered cubic local structure, even though the nucleus core is dominated by face-centered cubic structure corresponding to the stable rocksalt crystal structure. Finally, we establish a connection between the crystallization pathway and the equilibrium crystal-melt interface structure.

Applied Physical Sciences↗
DOE OSTI2025

Rapid Initial-State Preparation for the Quantum Simulation of Strongly Correlated Molecules

Studies on quantum algorithms for ground-state energy estimation often assume perfect ground-state preparation; however, in reality the initial state will have imperfect overlap with the true ground state. Here, we address that problem in two ways: by faster preparation of matrix-product-state (MPS) approximations and by more efficient filtering of the prepared state to find the ground-state energy. We show how to achieve unitary synthesis with a Toffoli complexity about 7 × lower than that in prior work and use that to derive a more efficient MPS-preparation method. For filtering, we present two different approaches: sampling and binary search. For both, we use the theory of window functions to avoid large phase errors and minimize the complexity. We find that the binary-search approach provides better scaling with the overlap at the cost of a larger constant factor, such that it will be preferred for overlaps less than about 0.003 . Finally, we estimate the total resources to perform ground-state energy estimation of Fe - S cluster systems, including the Fe Mo cofactor by estimating the overlap of different MPS initial states with potential ground states of the Fe Mo cofactor using an extrapolation procedure. With a modest MPS bond dimension of 4000 , our procedure produces an estimate of approximately 0.9 overlap squared with a candidate ground state of the Fe Mo cofactor, producing a total resource estimate of 7.3 × 10 10 Toffoli gates; neglecting the search over candidates and assuming the accuracy of the extrapolation, this validates prior estimates that have used perfect ground-state overlap. This presents an example of a practical path to prepare states of high overlap in a challenging-to-compute chemical system. Published by the American Physical Society 2025

Berry, Dominic W. (ORCID:0000000334461449)↗
DOE OSTI2025

Probing the Framework Metal Dependent Properties of Actinide-Centered Polyoxoalkoxide Sandwich-Type Complexes

Development of a simple and scalable synthesis of (TBA) 3 [W 5 O 18 MoNO] provides for the formation of the mixed-metal lacunary polyoxoalkoxide, (TBA) 2 [W 4 O 13 (OMe) 4 MoNO][Na(MeOH)]. This complex was used to synthesize a series of polyoxoalkoxide sandwich-type complexes with the general formula (TBA) 2 [M{W 4 O 13 (OMe) 4 MoNO} 2 ], where M = Zr(IV), Hf(IV), Th(IV), U(IV), and Np(IV). Compared to the analogous all-molybdenum complexes, the series have drastically different optical and redox properties. The results indicate that framework metal substitution acts as a tool for “orbital engineering”, with Density Functional Theory (DFT) calculations revealing that the major consequence of incorporation of tungsten into the complexes is localization of LUMO and LUMO+1 on the molybdenum centers remaining in the molecule. The change in the distribution of the frontier orbitals translates to discrepancies in the electronic properties of the series. Given the rarity of polyoxometalate complexes featuring a U(V) ion, one electron oxidation of (TBA) 2 [U(IV){W 4 O 13 (OMe) 4 MoNO} 2 ] was pursued. Isolation of the corresponding U(V) centered sandwich-type complex is reported, only the second example of U(V)-polyoxometalate complex described to date.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗