Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex systems 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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Nonlinear Topological Photonics: Capturing Nonlinear Dynamics and Optical Thermodynamics

Combining multiple optical resonators or engineering dispersion of complex media has provided an effective method for demonstrating topological physics controlling photons in unprecedented ways such as unidirectional light propagation and spatially localized modes between an interface or on a corner. Further, adding nonlinear responses to those topological photonic systems has enabled achieving diverse phases of photons in both space and time, allowing for more functionalities in photonic devices that provide a new playground for studying dynamic features of nonlinear topological systems. However, most methods for describing nonlinear topological photonic systems rely on linear topological theories, making it challenging to accurately characterize the topology of nonlinear systems. Thus, substantial efforts have focused on rigorously describing nonlinear topological phases and developing effective tools to analyze nonlinear topological effects. Meanwhile, coupled multimode optical waveguides with nonlinear dynamic responses provide an excellent platform for the statistical description of photons, opening a new paradigm called “optical thermodynamics”. This review will introduce the basic concepts of nonlinear topological photonics and the recent development of theoretical approaches focusing on data-driven approaches for creating phase diagrams as well as the spectral localizer framework and the pseudospectrum method for understanding optical nonlinearities in topological systems. In addition, the new concept of optical thermodynamics will be introduced with some recent theoretical works.

Topological photonics↗
DOE OSTI2025

Electronic structure complexity and extremely large magnetoresistance in antiferromagnetic semimetal SmAgSb 2

SmAgSb 2 , a layered magnetic semimetal in the tetragonal 𝑅⁢𝑇⁢ Sb 2 family (𝑅 = Y, Sc, rare earth; 𝑇 = transition metal), is known to exhibit extremely large magnetoresistance (XMR) below its antiferromagnetic (AFM) transition temperature. Here, in this work, we present a comprehensive investigation combining magnetotransport measurements, density functional theory calculations accounting for electron correlation, and angle-resolved photoemission spectroscopy. Our results reveal a complex electronic structure characterized by a multiband Fermi surface and intricate magnetic ground states. We demonstrate that simple two-band models, previously employed in the literature, fail to consistently describe the observed transport phenomena. Notably, we report an XMR of approximately 25200% at 2 K under a 14 T magnetic field, significantly exceeding earlier reports for this material family and rivaling the performance of prominent nonmagnetic XMR systems. This pronounced enhancement below 𝑇 𝑁 suggests that the XMR originates from a combination of multiband electron-hole compensation and enhanced magnetic scattering in this correlated AFM semimetal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Experimental and fuel-surrogates modeling study of the high-pressure pyrolysis of specialty cetane number fuels: implications for fall-off in ethylene unimolecular dissociation

Single pulse shock tube experiments were conducted at 50 atm nominal pressure and 4 ms nominal reaction time over a temperature range of 900–1800 K, to study the pyrolysis speciation of a multi-component jet fuel, F-24, and six cetane number (CN) specialty fuels - CN30, CN35, CN40, CN45, CN50, and CN55. Gas chromatography (GC) was used to qualitatively and quantitatively analyze the post shock gases. The relationship between the formation of key pyrolysis species and the chemically controlled combustion propensity as reflected by the cetane number of each fuel was examined. A surrogate-based mechanism from the CRECK Modelling Group and chemical-functional group based optimized surrogates (CFGO) were used to simulate the pyrolysis speciation results. The model was able to capture the chemistry of most species except two important pyrolysis intermediates – ethylene and acetylene. Chemical kinetic analyses were performed to identify the important reactions which affect the chemistry of these species; however, the rate parameters of critical reactions were found to be unsuitable for simulating the present high-pressure studies. Here, to address this unsuitability, a theory-based fall-off analysis for three reactions representing the decomposition of ethylene and subsequent formation of acetylene was performed, and these are included in an updated version of the CRECK mechanism. This update resolves discrepancies between the experimental results and simulations for ethylene and acetylene. Reaction flux analyses using the updated surrogate model were also performed to identify the important reaction pathways responsible for the formation of crucial species and to provide an analysis of the chemistry of complex multi-component fuel systems. The fundamental reactions responsible for driving pyrolysis chemistry were greatly influenced by the chemical functional groups present in these fuels. In addition to updating the rate parameters of specific reactions to improve modeling, this study also emphasizes the effectiveness of the fuel-surrogate approach, where surrogates representing the chemical functional group composition of the parent fuel serve as a valuable tool for predicting the combustion chemistry of novel fuels.

Chemical Kinetics↗
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 OSTI2021

Time resolved x-ray diffraction in shock compressed systems

The availability of pulsed x rays on short timescales has opened up new avenues of research in the physics and chemistry of shocked materials. The continued installation of shock platforms such as gas guns and high power lasers placed at beamline x-ray facilities has advanced our knowledge of materials shocked to extreme conditions of pressure and temperature. In addition, theoretical advancements have made direct correspondence with high-pressure x-ray experiments more viable, increasing the predictive capability of these models. In this paper, we discuss both recent experimental results and the theory and modeling that has been developed to treat these complex situations. Finally, we discuss the impact that new platforms and increased beam time may have on the future direction of this field.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Application of machine learning interatomic potentials in heterogeneous catalysis

Heterogeneous catalysts are crucial in modern societies as they promote sustainability by enabling lower-energy pathways for various chemical reactions. While Density Functional Theory (DFT) computations can provide critical insights into how heterogeneous catalysts operate at the atomic level, they are limited by computational costs and unfavorable scaling with system size. Recently, machine learning interatomic potentials (MLIPs) have emerged as a promising alternative to DFT, offering near-DFT accuracy at significantly reduced cost. Here, in this perspective, we discuss the application of MLIPs in heterogeneous catalyst modeling as a surrogate for DFT. We detail how MLIPs have been applied in thermal catalysis to probe active sites, enable studying complex metallic and nanoporous catalysts, and investigate the reconstruction of catalytic surfaces. We review the use of MLIPs in electrocatalysis and photocatalysis, emphasizing their capabilities in studying transition metal oxide surfaces and solid–liquid interfaces. We also discuss the current limitations of MLIPs, particularly their challenges with transferability and description of non-local interactions. Finally, we conclude by identifying promising and underexplored domains in which MLIPs can further advance our understanding of heterogeneous catalysts.

Catalytic surfaces↗
DOE OSTI2025

Simulation of 24,000 Electron Dynamics: Real-Time Time-Dependent Density Functional Theory (TDDFT) with the Real-Space Multigrids (RMG)

Here, we present the theory, implementation, and benchmarking of a real-time time-dependent density functional theory (RT-TDDFT) module within the RMG code, designed to simulate the electronic response of molecular systems to external perturbations. Our method offers insights into nonequilibrium dynamics and excited states across a diverse range of systems, from small organic molecules to large metallic nanoparticles. Benchmarking results demonstrate excellent agreement with established TDDFT implementations and showcase the superior stability of our time integration algorithm, enabling long-term simulations with minimal energy drift. The scalability and efficiency of RMG on massively parallel architectures allow for simulations of complex systems, such as plasmonic nanoparticles with thousands of atoms. Future extensions, including nuclear and spin dynamics, will broaden the applicability of this RT-TDDFT implementation, providing a powerful toolset for studies of photoactive materials, nanoscale devices, and other systems where real-time electronic dynamics is essential.

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

SAP-X2C: Optimally-Simple Two-Component Relativistic Hamiltonian with Size-Intensive Picture Change

We present a simple relativistic exact 2-component (X2C) Hamiltonian that models two-electron picture-change effects using Lehtola’s superposition of atomic potentials (SAP) [S. Lehtola, J. Chem. Theory Comput. 15, 1593−1604 (2019)]. The SAP-X2C approach retains the low cost and technical simplicity of the popular 1-electron X2C (1eX2C) predecessor but is significantly more accurate and has a well-defined thermodynamic limit, making it applicable to extended systems (such as large molecules and periodic crystals). The assessment of the SAP-X2C-based Hartree−Fock total and spinor energies, spin−orbit splittings, equilibrium bond distances, and harmonic vibrational frequencies suggests that SAP-X2C is similar to the more complex atomic meanfield (AMF) X2C counterparts in its ability to approximate the 4-component Dirac−Hartree−Fock reference.

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

A comprehensive approach for elucidating the interplay between 4f n +1 and 4f n 5d 1 configurations in Ln 2+ complexes

Lanthanides (Ln) are typically found in the +3 oxidation state. However, in recent decades, their chemistry has been expanded to include the less stable +2 oxidation state across the entire series except promethium (Pm), facilitated by the coordination of ligands such as trimethylsilylcyclopentadienyl, C 5 H 4 SiMe 3 (Cp'). The [LnCp' 3 ] complexes have been the workhorse for the synthesis and theoretical study of the fundamental aspects of divalent lanthanide chemistry, where experimental and computational evidence have suggested the existence of different ground state (GS) configurations, 4f n+1 or 4f n 5d 1 , depending on the specific metal. Standard reduction potentials and 4f n+1 to 4f n 5d 1 promotion energies have been two factors usually considered to rationalize the occurrence of these variable GS configurations, however the driving force behind this phenomenon is still not clear. In this work we present a comprehensive theoretical approach to shed light on this matter using the [LnCp 3 ] - model systems. We begin by calculating 4f n+1 to 4f n 5d 1 promotion energies and successfully correlate them with existing experimental data. Furthermore, we analyze how changes in the GS charge distribution between the Ln ions, LnCp 3 and the reduced [LnCp 3 ] - complexes (Ln = La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu) correlate with experimental trends in redox potentials and the calculated promotion energies. For this purpose, a comprehensive theoretical work that includes relativistic ligand field density functional theory (LFDFT) and relativistic ab initio wavefunction methods was performed. This study will help the rational design of suitable environments to tune the different GS configurations as well as modulating the spectroscopic properties of new Ln 2+ complexes.

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

Element-Specific Study of Magnetic Anisotropy and Hardening in SmCo 5– x Cu x Thin Films

This work investigates the effect of copper substitution on the magnetic properties of SmCo 5 thin films synthesized by molecular beam epitaxy. A series of thin films with varying concentrations of Cu were grown under otherwise identical conditions to disentangle structural and compositional effects on the magnetic behavior. The combined experimental and theoretical studies show that Cu substitution at the Co 3g sites not only stabilizes the formation of the SmCo 5 structure but also enhances magnetic anisotropy and coercivity. Density functional theory calculations indicate that Sm(Co 4 Cu 3g ) 5 possesses a higher single-ion anisotropy as compared to pure SmCo 5 . In addition, X-ray magnetic circular dichroism reveals that Cu substitution causes an increasing decoupling of the Sm 4f and Co 3d moments. Scanning transmission electron microscopy confirms predominantly SmCo 5 phase formation and reveals nanoscale inhomogeneities in the Cu and Co distribution. Our study based on thin film model systems and advanced characterization as well as modeling reveals novel aspects of the complex interplay of intrinsic and extrinsic contributions to magnetic hysteresis in rare-earth-based magnets, i.e., the combination of increased intrinsic anisotropy due to Cu substitution and the extrinsic effect of inhomogeneous elemental distribution of Cu and Co.

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

MBX V1.2: Accelerating Data-Driven Many-Body Molecular Dynamics Simulations

The MBX software provides an advanced platform for molecular dynamics simulations, leveraging state-of-the-art MB-pol and MB-nrg data-driven many-body potential energy functions. Developed over the past decade, these potential energy functions integrate physics-based and machine-learned many-body terms trained on electronic structure data calculated at the "gold standard" coupled-cluster level of theory. Recent advancements in MBX have focused on optimizing its performance, resulting in the release of MBX v1.2. While the inherently many-body nature of MB-pol and MB-nrg ensures high accuracy, it poses computational challenges. MBX v1.2 addresses these challenges with significant performance improvements, including enhanced parallelism that fully harnesses the power of modern multicore CPUs. In conclusion, these advancements enable simulations on nanosecond time scales for condensed-phase systems, significantly expanding the scope of high-accuracy, predictive simulations of complex molecular systems powered by data-driven many-body potential energy functions.

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

The interplay of excitonic delocalization and vibrational localization in optical lineshapes: A variational polaron approach

The dynamics of molecular excitonic systems are complicated by a competition between electronic coupling (which drives delocalization) and vibrational-electronic (vibronic) interactions (which tend to encourage electronic localization). A particular challenge of molecular systems is that they typically possess a large number of independent vibrations, with frequencies often spanning the entire spectrum of relevant electronic energy gaps. Recent spectroscopic observations and numerical simulations on a water-soluble chlorophyll-binding protein (WSCP) reveal a transition between two regimes of vibronic behavior, a Redfield-like regime in which low-frequency vibrations respond to a delocalized excitonic state, and a Förster-like regime where high-frequency vibrations act as incoherent excitations on individual pigments. Although numerical simulations can reproduce these effects, there is a need for a simple, systematic theory that accurately describes the smooth transition between these two regimes in experimental spectra. Here we address this challenge by generalizing the variational polaron transform approach of [Bloemsma et al., Chem. Phys. 481, 250 (2016)] to include arbitrary bath densities for systems with or without symmetry. We benchmark this theory against both numerical matrix-diagonalization methods and experimental 77 K fluorescence spectra for two WSCP variants, obtaining quite satisfactory agreement in both cases. Here, we apply this theory to offer an explanation for the large loss in apparent electronic coupling in the WSCP Q57K mutant and to examine the likely impact of the interplay between excitonic delocalization and vibrational localization on vibrational sideband shapes and apparent coupling strengths in high-resolution optical spectra for chlorophyll-protein complexes such as WSCP.

14 SOLAR ENERGY↗
DOE OSTI2025

Web-Based Tools for Data-Informed Remedy Optimization: Software Theory and User Guide

This report documents the development and application of two web-based decision-support tools for pump-and-treat (P&T) groundwater remediation systems: PTOLEMY (Pump-and-Treat Optimized Location Evaluation to Maximize Yields) and OPTIMA (Optimization for Pump-and-Treat Implementation, Management, & Assessment). These tools enhance remedy design and management by leveraging advanced computational methods – specifically deep learning and multi-objective optimization – within a user-friendly platform. By integrating data-driven models with established hydrogeological knowledge, PTOLEMY and OPTIMA enable more efficient evaluation of well placement and operational strategies, helping site managers balance multiple remediation objectives under complex conditions. Both tools are implemented as modules within the SOCRATES (Suite Of Comprehensive Rapid Analysis Tools for Environmental Sites) web platform, which provides data access, visualization, and analytics to support remedy optimization across sites in the U.S. Department of Energy Office of Environmental Management complex. PTOLEMY is a rapid screening module designed to identify promising locations for new extraction wells. It employs a multi-channel three-dimensional convolutional neural network (MC3D-CNN) trained on high-fidelity simulation data to predict the relative performance (in terms of contaminant mass recovery) of potential well sites. Through an interactive web interface, PTOLEMY visualizes the probability of high performance across a site, highlighting areas where an extraction well is likely to yield above-threshold contaminant removal over a multi-year period. PTOLEMY’s map-based displays and exportable results support transparent communication of screening analyses. By focusing attention on the most favorable candidate locations, the tool augments traditional engineering judgment and physics-based modeling, providing a data informed basis for subsequent detailed evaluations. OPTIMA is a multi objective optimization module designed to find wellfield layouts and operating schedules that meet various cleanup goals. It quickly evaluates thousands of candidate setups – combinations of well locations, timing, and rates – and returns a small set of best trade-off options for comparison. At its core, OPTIMA uses a U-Net-based surrogate model – a deep-learning emulator of a groundwater flow and transport simulator – to dramatically accelerate scenario evaluations. Coupling this fast surrogate with the NSGA-II (Non-dominated Sorting Genetic Algorithm II) evolutionary algorithm, OPTIMA explores a wide decision space of well locations and schedules to identify Pareto-optimal solutions that trade off key objectives (e.g., minimizing cleanup time, maximizing contaminant mass removal, and minimizing plume extent). The tool outputs a family of optimal configurations and visualizes their trade-offs (Pareto frontiers of cleanup metrics and maps of optimized well placements). Site managers can use these results to understand the range of viable strategies and to select candidate designs for more detailed verification. OPTIMA is currently under active development and not yet fully released; this guide provides early documentation to support planning and gather user feedback.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Strategies to Obtain Reliable Energy Landscapes from Embedded Multireference Correlated Wavefunction Methods for Surface Reactions

Embedded correlated wavefunction (ECW) theory is a powerful tool for studying ground- and excited-state reaction mechanisms and associated energetics in heterogeneous catalysis. Several factors are important to obtaining reliable ECW energies, critically the construction of consistent active spaces (ASs) along reaction pathways when using a multireference correlated wavefunction (CW) method that relies on a subset of orbital spaces in the configuration interaction expansion to account for static electron correlation, e.g., complete AS self-consistent field theory, in addition to the adequate partitioning of the system into a cluster and environment, as well as the choice of a suitable basis set and number of states included in excited-state simulations. Here, in this work, we conducted a series of systematic studies to develop best-practice guidelines for ground- and excited-state ECW theory simulations, utilizing the decomposition of NH 3 on Pd(111) as an example. We determine that ECW theory results are relatively insensitive to cluster size, the aug-cc-pVDZ basis set provides an adequate compromise between computational complexity and accuracy, and that a fixed-clean-surface approximation holds well for the derivation of the embedding potential. Additionally, we demonstrate that a merging approach, which involves generating ASs from the molecular fragments at each configuration, is preferable to a creeping approach, which utilizes ASs from adjacent structures as an initial guess, for the generation of consistent potential energy curves involving open-d-shell metal surfaces, and, finally, we show that it is essential to include bands of excited states in their entirety when simulating excited-state reaction pathways.

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

Actinide Molten Salts: A Machine-Learning Potential Molecular Dynamics Study

We know that actinide molten salts represent a class of important materials in nuclear energy. Understanding them at a molecular level is critical to proper and optimal design of relevant technological applications. Yet, owing to the complexity of electronic structure due to the 5f orbitals, computational studies of heavy elements in condensed phases using ab initio potentials to study the structure and dynamics of these elements embedded in molten salts are difficult. This lack of efficient computational protocols makes it difficult to obtain information on properties that require extensive statistical sampling like transport. To tackle this problem, we adopted a machine-learning approach to study ThCl 4 -NaCl and UCl 3 -NaCl binary systems. The machine-learning potential, with the density functional theory accuracy, allows us to obtain long molecular dynamics trajectories (ns) for large systems (10 3 atoms) at a considerably low computing cost, thereby efficiently gaining information about their bonding structures, thermodynamics, and dynamics at a range of temperature. We observed a considerable change in the coordination environments of actinide elements and their characteristic coordination-sphere lifetime. Our study also suggests that actinides in molten salts may not follow well known entropy-scaling laws.

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

Ultrafast Optically Induced Perturbation of Oxygen Octahedral Rotations in Multiferroic BiFeO 3 Thin Films

The functional properties of complex oxides, including magnetism and ferroelectricity, are closely linked to subtle structural distortions. Ultrafast optical excitations provide the means to manipulate structural features and ultimately to affect the functional properties of complex oxides with picosecond-scale precision. Here, we report that the lattice expansion of multiferroic BiFeO 3 following above-bandgap optical excitation leads to distortion of the oxygen octahedral rotation (OOR) pattern. The continuous coupling between OOR and strain was probed using time-resolved X-ray free-electron laser diffraction with femtosecond time resolution. Density functional theory calculations predict a relationship between the OOR and the elastic strain consistent with the experiment, demonstrating a route to employing this approach in a wider range of systems. Ultrafast control of the functional properties of BiFeO 3 thin films is enabled by this approach because the OOR phenomena are related to ferroelectricity, and via the Fe-O-Fe bond angles, the superexchange interaction between Fe atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Spin polarization engineering in 𝑑-wave altermagnets

Altermagnets host unconventional spin-polarized bands despite zero net magnetization, but controlling their spin structure remains challenging. Here, we theoretically propose a multifield approach to engineer spin polarization in 𝑑-wave altermagnets using gating, optical driving, and in-plane electric fields, which enable tunable and switchable polarizations along multiple directions. Optical driving induces out-of-plane (𝑧) polarization, while gating and in-plane fields generate 𝑥 and 𝑦 polarizations via the Edelstein effect, all of which are experimentally detectable. We further find that spin- and band-selective doping induces chiral optical activity, a feature unique to altermagnets. Our approach offers a versatile means to control the direction of spin polarization in altermagnets.

Floquet systems↗
DOE OSTI2024

A computational approach to complex junctions and interfaces

Molecular- and nano- junctions, surfaces, interfaces, and interfacial processes have been the subjects of this DOE supported program for about 20 years. We have investigated electronic and magnetic structure, spin-dependent charge transport, and emerging phenomena at surfaces and interfaces in a number of physical systems ranging from nano-clusters, to nanowires, to two-dimensional (2D) crystals and their interactions, to bulk matter. Eighty-four refereed papers were published [1-84], with a majority effort from the group led by the PI, and over a hundred presentations were made at international conferences. We use first-principles methods in the framework of density functional theory in conjunction with nano-equilibrium Green function techniques, semi-empirical Boltzmann transport theory, tight-binding models, and the Hubbard model and many-body wavefunction approaches. In pursuing our scientific goals, we have also developed a number of computational algorithms and computational packages.

2D materials↗