Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Functional coupling”

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

Nonadiabatic couplings from a variational excited state method based on constrained DFT

Excited Costrained Density Functional Theory (XCDFT) is a variational excited state method that extends ground state DFT to the computation of low-lying excited states. It borrows much of the machinery of Constrained DFT (CDFT) with a crucial difference: the constraint imposes a population of one electron in the Hilbert space spanned by the virtuals of a reference ground state. Here, we present theory and implementation for evaluating nonadiabatic coupling vectors (NACVs) between the first excited state computed with XCDFT and the ground state. Our NACVs are computed analytically using density functional perturbation theory with a formalism that is general enough that could be applied to CDFT diabatic states. We showcase the new method with pilot NACV calculations for the conical intersection in H 3 , the avoided crossing in selenoacrolein, and the NACV magnitudes in azobenzene. Despite complications from the nonorthogonality of the wavefunctions, XCDFT’s energy surfaces and NACVs reproduce benchmark values and respect known sum rules within a reasonable degree. This shows that XCDFT is a viable method for nonadiabatic dynamics simulations.

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

Attaining regularization length insensitivity in phase-field models of ductile failure

In this work, a cohesive phase-field model of ductile fracture in a finite-deformation setting is presented. The model is based on a free-energy function in which both elastic and plastic work contributions are coupled to damage. Using a strictly variational framework, the field evolution equations, damage kinetics, and flow rule are jointly derived from a scalar least-action principle. Particular emphasis is placed on the use of a rational function for the stress degradation that maintains a fixed effective strength with decreasing regularization length. The model is employed to examine crack growth in pure mode-I problems through the generation of crack growth resistance (J-R) curves. In contrast to alternative models, the current formulation gives rise to J-R curves that are insensitive to the regularization length. Numerical evidence suggests convergence of local fields with respect to diminishing regularization length as well.

42 ENGINEERING↗
DOE OSTI2026

Single-Phase to Split-Phase Inverters with Advanced Grid Support Functions for Grid-Interactive Applications

This work presents a cost-effective single-phase to split-phase inverter with a reduced switch count, achieving grid interactive performance while maintaining operational efficiency. The proposed system integrates an Andronov-Hopf oscillator based secondary controller, which inherently embeds a nonlinear resistive droop architecture, ensuring rapid dynamic response. A Lyapunov energy function-based primary control enhances transient stability and regulation, while an internal model-based point of common coupling voltage estimation enables cost optimization without additional sensors. Equipped with advanced grid support functionalities, the inverter facilitates seamless distribution system operation with enhanced robustness. The effectiveness of the proposed architecture and control strategy is validated through MATLAB/Simulink and PLECS simulations, demonstrating its feasibility for high-performance grid-supportive applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Theoretical framework for new magnetic materials for quantum computing and information storage. Final report for the Award No. DE-SC0018910

The focus of this grant was on molecular magnetic materials for information storage and quantum computing. We have been developing robust, first-principle methods for computing relevant electronic and magnetic properties of molecular building blocks (SMMs) of novel magnetic materials and quantum computers. These tools enable theoretical modeling of SMMs’ behavior, facilitating the interpretation of experimental studies and aiding the design of novel magnetic materials. Our strategy is based on the spin-flip (SF) approach, which extends the hierarchy of black-box single-reference methods to strongly correlated systems. Specifically, we developed general scalable algorithms and computer codes for calculating molecular properties, with an emphasis on spin-related properties, such as zero-field splittings, hyperfine couplings, and g-tensors. While our primary focus was on SF wave functions and SF-TDDFT, the underlying theory and computer codes were formulated using reduced density matrices, such that these tools are applicable to a broader class of methods. To extend the scope of applicability of wave-function-based SF methods to larger systems, we developed reduced-scaling approaches for the equation-of-motion coupled-cluster (EOM-CC) methods and continue developing libtensor (our open-source general tensor contraction library for many-body methods). We carried out extensive benchmarks and also carried out several applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Coupling Noah-Multiparameterization land-surface Model with Energy Research and Forecasting Model

The Energy Research and Forecasting (ERF) model is a high-performance atmospheric model built on the AMReX adaptive mesh refinement (AMR) framework, enabling efficient simulations on heterogeneous computing platforms that combine multicore processors with hardware accelerators. To support land–atmosphere interactions within ERF’s AMR-based environment, a land-surface model must be capable of operating directly on hierarchically refined meshes. In this work, we present a methodology for coupling the Fortran-based Noah-Multiparameterization (Noah-MP) land-surface model with ERF’s C++ codebase. Rather than rewriting Noah-MP, we construct a Fortran–C interoperability layer using CodeScribe, a tool that leverages large language models (LLMs) to automate the generation of interface code. CodeScribe applies structured prompting techniques to generate bindings that support efficient data exchange and function calls between ERF and Noah-MP. The coupling framework also incorporates AMR-aware data handling strategies, allowing NoahMP to operate seamlessly within ERF’s hierarchical mesh structure. This work provides a structured approach for integrating legacy Fortran models into modern C++-based modeling systems using LLM-assisted code generation.

54 ENVIRONMENTAL SCIENCES↗
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↗
DOE OSTI2025

16 O Electroweak Response Functions from First Principles

We present calculations of various electroweak response functions for the 16 O nucleus obtained using coupled-cluster theory in conjunction with the Lorentz integral transform method. We employ nuclear forces derived at next-to-leading order and next-to-next-to-leading order in chiral effective field theory and perform a Bayesian analysis to assess uncertainties. Our results are in good agreement with available electron-scattering data at |𝐪|≈326 MeV/c. Additionally, we provide several predictions for the weak response functions in the quasielastic peak region at |𝐪| =300 and 400 MeV/c, which are critical for long-baseline neutrino experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Computational methods based on density functional theory for reactions and processes involving electronic spin (Final Technical Report)

This award supports one post-doctoral researcher for 1.5 years. Publications that acknowledge this grant: Refs. 1–14. Refs. 3,10,12 assess current methodology for the evaluation of magnetic exchange couplings in transition metal complexes. In particular, Ref. 10 validates the use of an approximate (non-iterative) Green’s function approach for the calculation of magnetic exchange couplings and will be the foundation for Thrust 2 in this proposal. Refs. 3 and 12 focus on widely used density functional approaches based on the standard energy differences methodology for the particular case of oxo-bridged Fe(III) complexes. Refs. 2,4–7,11 apply current methodologies to problems of practical interest in molecular magnetism. Ref. 13 presents a methodology to explicitly simulate the dynamics of open quantum systems within density functional theory (DFT) calculations based on the Liouville-von Neumann equation of motion for quantum systems driven out-of-equilibrium. Ref. 8 uses non-collinear spin DFT to explain the mechanical behavior of magnetic mono-atomic Pt wires produced in break-junction experiments in the presence of a magnetic field.

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

Constructing tensor network influence functionals for general quantum dynamics

Here, we describe an iterative formalism to compute influence functionals that describe the general quantum dynamics of a subsystem beyond the assumption of linear coupling to a quadratic bath. We use a space-time tensor network representation of the influence functional and investigate its approximability in terms of its bond dimension and time-like entanglement in the tensor network description. We study two numerical models, the spin-boson model and a model of interacting hard-core bosons in a 1D harmonic trap. We find that the influence functional and the intermediates involved in its construction can be efficiently approximated by low bond dimension tensor networks in certain dynamical regimes, which allows the quantum dynamics to be accurately computed for longer times than with direct time evolution methods. However, as one iteratively integrates out the bath, the correlations in the influence functional can first increase before decreasing, indicating that the final compressibility of the influence functional is achieved via non-trivial cancellation.

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

Local Structural Coherence and Interfacial Charge Transfer in Cu 2 ⁢S/Mo⁢S 2 Heterostructure

Precise control over electronic coupling at nanoscale interfaces is critical for designing materials with tunable charge-transfer behavior and catalytic function. Heterostructures with locally coherent interfaces provide a platform for interrogating interfacial charge redistribution in coupled material systems. Here, we report Cu 2 ⁢S/Mo⁢S 2 heterostructures exhibiting nanoscale crystallographic alignment, which are synthesized through a rapid thermal transformation pathway. We employed electrochemical reduction reactions to probe interfacial charge transfer, revealing shifts in product distribution attributable to modified interfacial energetics, even in the absence of optimized catalytic performance. The observed formate Faradaic efficiency suggests that interfacial electronic modulation in the heterostructure shifts product selectivity towards formate, highlighting how interface-driven electronic modulation can direct reaction pathways and influence product selectivity. Optimizing catalyst loading, architecture, and reactor configuration will be critical for future improvement. Structural and compositional integrity were confirmed through powder x-ray diffraction, x-ray photoelectron spectroscopy, and high-resolution transmission electron microscopy. Electron transfer between Cu 2 ⁢S and Mo⁢S 2 domains was further evaluated by electrochemical impedance spectroscopy, while selected-area electron diffraction revealed local crystallographic alignment consistent with a local epitaxial relationship at the Cu 2 S/Mo⁢S 2 heterointerface. To illustrate the broader applicability of this approach, a Zn⁢S/Mo⁢S 2 heterostructure was also synthesized using the same microwave strategy, confirming the generalizability of interfacial engineering principles across metal sulfide-Mo⁢S 2 systems. Collectively, these findings demonstrate that controlled local epitaxial alignment serves as an effective design principle for tuning interfacial energetics and catalytic reactivity in complex heterostructure materials.

carbon capture & utilization↗
DOE OSTI2023

Cosmic variance and the inhomogeneous UV luminosity function of galaxies during reionization

When the first galaxies formed and starlight escaped into the intergalactic medium to reionize it, galaxy formation and reionization were both highly inhomogeneous in time and space, and fully coupled by mutual feedback. To show how this imprinted the UV luminosity function (UVLF) of reionization-era galaxies, we use our large-scale, radiation-hydrodynamics simulation CoDa II to derive the time- and space-varying halo mass function and UVLF, from z ≃ 6–15. That UVLF correlates strongly with local reionization redshift: earlier-reionizing regions have UVLFs that are higher, more extended to brighter magnitudes, and flatter at the faint end than later-reionizing regions observed at the same z. In general, as a region reionizes, the faint-end slope of its local UVLF flattens, and, by z = 6 (when reionization ended), the global UVLF, too, exhibits a flattened faint-end slope, ‘rolling-over’ at M UV ≳ -17. CoDa II’s UVLF is broadly consistent with cluster-lensed galaxy observations of the Hubble Frontier Fields at z = 6–8, including the faint end, except for the faintest data point at z = 6, based on one galaxy at M UV = -12.5. According to CoDa II, the probability of observing the latter is $\sim 5~{{\ \rm per\ cent}}$. However, the effective volume searched at this magnitude is very small, and is thus subject to significant cosmic variance. Here, we find that previous methods adopted to calculate the uncertainty due to cosmic variance underestimated it on such small scales by a factor of 2–4, primarily by underestimating the variance in halo abundance when the sample volume is small.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Controlled n –Doping of Naphthalene Diimide–Based Two–Dimensional Polymers

Two-dimensional polymers (2DPs) are promising as structurally well-defined, permanently porous, organic semiconductors. However, 2DPs are nearly always isolated as closed shell organic species with limited charge carriers, which leads to low bulk conductivities. Here, we enhance the bulk conductivity of two naphthalene diimide (NDI)-containing 2DP semiconductors by controllably n-doping the NDI units using cobaltocene (CoCp 2 ). Optical and transient microwave spectroscopy reveals that both as-prepared NDI-containing 2DPs are semiconducting with sub-2 eV optical bandgaps and photoexcited charge-carrier lifetimes of tens of nanoseconds. Following reduction with CoCp 2 , both 2DPs largely retain their periodic structures and exhibit optical and electron-spin resonance spectroscopic features consistent with the presence of NDI-radical anions. Here, while the native NDI-based 2DPs are electronically insulating, maximum bulk conductivities of >10 –4 S cm –1 are achieved by substoichiometric levels of n-doping. Density functional theory calculations show that the strongest electronic couplings in these 2DPs exist in the out-of-plane (π-stacking) crystallographic directions, which indicates that cross-plane electronic transport through NDI stacks is primarily responsible for the observed electronic conductivity. Taken together, this study underlines that controlled molecular doping is a useful approach to access structurally well-defined, paramagnetic, 2DP n-type semiconductors with measurable bulk electronic conductivities of interest for electronic or spintronic devices.

2D polymers↗
DOE OSTI2021

Approximate Bacon-Shor code and holography

We explicitly construct a class of holographic quantum error correction codes with non-trivial centers in the code subalgebra. Specifically, we use the Bacon-Shor codes and perfect tensors to construct a gauge code (or a stabilizer code with gauge-fixing), which we call the holographic hybrid code. This code admits a local log-depth encoding/decoding circuit, and can be represented as a holographic tensor network which satisfies an analog of the Ryu-Takayanagi formula and reproduces features of the sub-region duality. We then construct approximate versions of the holographic hybrid codes by “skewing” the code subspace, where the size of skewing is analogous to the size of the gravitational constant in holography. These approximate hybrid codes are not necessarily stabilizer codes, but they can be expressed as the superposition of holographic tensor networks that are stabilizer codes. For such constructions, different logical states, representing different bulk matter content, can “back-react” on the emergent geometry, resembling a key feature of gravity. The locality of the bulk degrees of freedom becomes subspace-dependent and approximate. Such subspace-dependence is manifest from the point of view of the “entanglement wedge” and bulk operator reconstruction from the boundary. Exact complementary error correction breaks down for certain bipartition of the boundary degrees of freedom; however, a limited, state-dependent form is preserved for particular subspaces. We also construct an example where the connected two-point correlation functions can have a power-law decay. Coupled with known constraints from holography, a weakly back-reacting bulk also forces these skewed tensor network models to the “large N limit” where they are built by concatenating a large N number of copies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Bending-induced isostructural transitions in ultrathin layers of van der Waals ferrielectrics

Using Landau-Ginzburg-Devonshire (LGD) phenomenological approach we analyze the bending-induced re-distribution of electric polarization and field, elastic stresses and strains inside ultrathin layers of van der Waals ferrielectrics. We consider a CuInP2S6 (CIPS) thin layer with fixed edges and suspended central part, the bending of which is induced by external forces. The unique aspect of CIPS is the existence of two ferrielectric states, FI1 and FI2, corresponding to big and small polarization values, which arise due to the specific four-well potential of the eighth-order LGD functional. When the CIPS layer is flat, the single-domain FI1 state is stable in the central part of the layer, and the FI2 states are stable near the fixed edges. With an increase of the layer bending below the critical value, the sizes of the FI2 states near the fixed edges decreases, and the size of the FI1 region increases. When the bending exceeds the critical value, the edge FI2 states disappear being substituted by the FI1 state, but they appear abruptly near the inflection regions and expand as the bending increases. The bending-induced isostructural FI1-FI2 transition is specific for the bended van der Waals ferrielectrics described by the eighth (or higher) order LGD functional with consideration of linear and nonlinear electrostriction couplings. Furthermore, the isostructural transition, which is revealed in the vicinity of room temperature, can significantly reduce the coercive voltage of ferroelectric polarization reversal in CIPS nanoflakes, allowing for the curvature-engineering control of various flexible nanodevices.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Beyond Component Optimization: Systems Level Biodesign for Lanthanide Recovery

Global demand for lanthanides (Ln) is projected to rise sharply over the next decade, while geographically concentrated supply chains and the low concentrations and matrix complexity of secondary feedstocks limit the reach of conventional hydro- and pyrometallurgical separation. Engineered biological systems offer a selective, low-energy alternative, and component-level advances in Ln-binding proteins, AI-designed selective scaffolds, and cell-surface display platforms now rival synthetic chelators in affinity and selectivity. These components, however, remain functionally isolated. Currently, there are no engineered chassis coupling recognition, intracellular trafficking, accumulation, and controlled release into an end-to-end pipeline. Here, we outline how new biodesign strategies and chassis selection must move beyond bioleaching to encompass the full recovery pathway. Achieving this requires integrating AI/ML-guided design, genome-scale build tools, high-throughput phenotyping, and biophysical transport modeling within a Design–Build–Test–Learn cycle tuned to recognition, trafficking, accumulation, and release.

Biodesign↗
DOE OSTI2020

Heat flow enhancement in a nanoscale plasmonic junction induced by Kondo resonances and electron-phonon coupling

Here recently, we showed that plasmon-exciton coupling can increase entropy current through a bridge coupled to plasmonic metal nanoparticles. Here we show that electron-phonon coupling can also be used to control the entropy current in similar systems. Entropy current tends to decrease due to electron-phonon coupling and to exhibit a monotonous decrease upon temperature ramping. However, an anomaly affecting the current where it is enhanced by electron-phonon coupling is indicated at around 42 times the system's Kondo temperature. We therefore report means to control heat flow by tuning the Kondo resonance through the electron-phonon coupling. We analyze the conditions that bring about these trends due to electron-phonon coupling by employing non-equilibrium Green's function formulation addressing the entropy current and the derived heat flow.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Defect Stability in CdTe Based on Formation Energies and Migration Barriers

Native point defects are thought to play a key role in CdTe, either as compensation centers in intentionally doped material, as a source of conductivity in nominally undoped material, or as electron–hole recombination centers. Here, the discussion of their concentration and impact has often centered only on formation energies and transition levels. Using hybrid density functional calculations, including the effects of spin–orbit coupling (SOC), we discuss the stability of native point defects in CdTe based on their formation energies and migration barriers. We show that although Cd interstitials are the lowest energy donor defects, they are unstable at room temperature due to a low migration barrier. They are important for maintaining charge neutrality during growth or annealing at high temperatures, but once the material is brought to room temperature, they are not frozen in as often assumed and are expected to anneal out, leaving the other more stable defects to determine the conductivity. Taking this into account in the solution of the charge neutrality equation, we are able to predict the conductivity type and carrier concentrations that are in good agreement with experimental observations.

14 SOLAR ENERGY↗