Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “path integrals”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Electric current paths in a Si:P delta-doped device imaged by nitrogen-vacancy diamond magnetic microscopy

The recently-developed ability to control phosphorous-doping of silicon at an atomic level using scanning tunneling microscopy, a technique known as atomic precision advanced manufacturing (APAM), has allowed us to tailor electronic devices with atomic precision, and thus has emerged as a way to explore new possibilities in Si electronics. In these applications, critical questions include where current flow is actually occurring in or near APAM structures as well as whether leakage currents are present. In general, detection and mapping of current flow in APAM structures are valuable diagnostic tools to obtain reliable devices in digital-enhanced applications. In this report, we used nitrogen-vacancy (NV) centers in diamond for wide-field magnetic imaging (with a few-mm field of view and micron-scale resolution) of magnetic fields from surface currents flowing in an APAM test device made of a P delta-doped layer on a Si substrate, a standard APAM witness material. We integrated a diamond having a surface NV ensemble with the device (patterned in two parallel mm-sized ribbons), then mapped the magnetic field from the DC current injected in the APAM device in a home-built NV wide-field microscope. The 2D magnetic field maps were used to reconstruct the surface current densities, allowing us to obtain information on current paths, device failures such as choke points where current flow is impeded, and current leakages outside the APAM-defined P-doped regions. Analysis on the current density reconstructed map showed a projected sensitivity of ~0.03 A m -1 , corresponding to a smallest-detectable current in the 200 μm wide APAM ribbon of ~6 μA. These results demonstrate the failure analysis capability of NV wide-field magnetometry for APAM materials, opening the possibility to investigate other cutting-edge microelectronic devices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Report of the 2025 Workshop on Next-Generation Ecosystems for Scientific Computing: Harnessing Community, Software, and AI for Cross-Disciplinary Team Science

This report summarizes insights from the 2025 Workshop on Next-Generation Ecosystems for Scientific Computing: Harnessing Community, Software, and AI for Cross-Disciplinary Team Science, which convened more than 40 experts from national laboratories, academia, industry, and community organizations to chart a path toward more powerful, sustainable, and collaborative scientific software ecosystems. To address urgent challenges at the intersection of high-performance computing (HPC), AI, and scientific software, participants envisioned agile, robust ecosystems built through socio-technical co-design—the intentional integration of social and technical components as interdependent parts of a unified strategy. This approach combines advances in AI, HPC, and software with new models for cross-disciplinary collaboration, training, and workforce development. Key recommendations include building modular, trustworthy AI-enabled scientific software systems; enabling scientific teams to integrate AI systems into their workflows while preserving human creativity, trust, and scientific rigor; and creating innovative training pipelines that keep pace with rapid technological change. Pilot projects were identified as near-term catalysts, with initial priorities focused on hybrid AI/HPC infrastructure, cross-disciplinary collaboration and pedagogy, responsible AI guidelines, and prototyping of public-private partnerships. This report presents a vision of next-generation ecosystems for scientific computing where AI, software, hardware, and human expertise are interwoven to drive discovery, expand access, strengthen the workforce, and accelerate scientific progress.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

High-precision phase control of an optical lattice with up to 50 dB noise suppression

An optical lattice is a periodic light crystal constructed from the standing-wave interference patterns of laser beams. It can be used to store and manipulate quantum degenerate atoms and is an ideal platform for the quantum simulation of many-body physics. A principal feature is that optical lattices are flexible and possess a variety of multidimensional geometries with modifiable band structure. An even richer landscape emerges when control functions can be applied to the lattice by modulation of the position or amplitude with Floquet driving. However, the desire of realizing high-modulation bandwidths while preserving extreme lattice stability has been difficult to achieve. In this paper, we demonstrate an effective solution that consists of overlapping two counterpropagating lattice beams and controlling the phase and intensity of each with independent acousto-optic modulators. Our phase controller mixes sampled light from both lattice beams with a common optical reference. This dual heterodyne locking method allows exquisite determination of the lattice position, while also removing parasitic phase noise accrued as the beams travel along separate paths. Here, we report up to 50 dB suppression in lattice phase noise in the 0.1–1 Hz band, along with significant suppression spanning more than 4 decades of frequency. The absolute phase diffusion of the lattice position is only 10 Å when integrated over 10 s. This method permits precise, high-bandwidth modulation (greater than 50 kHz) of the optical lattice intensity and phase. We demonstrate the efficacy of this approach by executing intricate time-varying phase profiles for atom interferometry.

Atom interferometry↗
DOE OSTI2023

Control of power converters having integrated capacitor blocked transistor cells

Various examples of power converters including Integrated Capacitor Blocked Transistor (ICBT) cells and methods of control of power converters having ICBT cells are described. In one example, a power converter includes an upper arm including a plurality of upper ICBT cells connected in series to form a series connection path and a lower arm including a plurality of lower ICBT cells connected in series in the series connection path. A controller can be configured to provide a control signal pair to each of the upper ICBT cells and a complementary control signal pair to each of the lower ICBT cells to control the converter output. A capacitor voltage controller can be configured to balance a voltage potential among ICBT capacitors in at least one of the upper arm and the lower arm.

Yu, Jianghui↗
DOE OSTI2025

UltraLiM: In-Memory Boolean Logic Architecture Using UltraRAM

Conventional computing architectures encounter ‘von Neumann’ and ‘memory wall’ bottlenecks which arise due to the back-and-forth data movement between the physically separate memory and processing units and the speed mismatch between them, respectively. These bottlenecks hurt both energy efficiency and the throughput of computing systems. To address these challenges, in-memory computing architectures have emerged as a promising alternative. They reduce the need for frequent data movement by executing different computing tasks inside the memory system. Here, we present UltraLiM, a logic-in-memory architecture using the UltraRAM-based memory system. UltraRAM holds the promise of developing a ‘universal memory’, overcoming the limitations of charge-based memories thanks to their non-volatile behavior with lower operating voltage. This work presents an in-memory computing architecture that integrates an UltraRAM-based memory array with a custom-designed peripheral circuitry. With this architecture, we can perform various in-memory Boolean logic operations (such as NOT, NAND, NOR, and XOR) in a single cycle. Leveraging the separate read-write paths in the UltraRAM-based memory array, we optimize read operations without encountering design conflicts. This optimization enhances the sense margin, enabling the use of simpler peripheral circuitry for in-memory logic operations.

Alam, Shamiul [University of Tennessee, Knoxville ↗
DOE OSTI2025

Metal-Insulation REBCO Pancake Coil With Solder Surface Shunt: Testing and Modeling

Rare-earth barium copper oxide (REBCO) high-temperature superconducting (HTS) magnets can generate very high magnetic fields with superior engineering current density and robust mechanical integrity. The no-insulation (NI) winding technique, including its metal-insulation (MI) variant, enables “self-protection” against localized overheating through inter-turn current sharing, improving defect tolerance and allowing temporary current overloads. In REBCO pancake coils, the current-sharing path is predominantly confined to the edges of the REBCO tapes due to insulating buffer layers. Applying a surface shunt using solder can effectively enhance current sharing in NI/MI coils via the REBCO tape edges. Here, in this paper, we present both experimental and numerical research on the current distribution in surface-shunt metal-insulation (SSMI) REBCO pancake coils. We tested a stainless steel co-wound MI REBCO double-pancake coil, before and after applying a soldered surface shunt, in a saturated liquid nitrogen bath. The SSMI coil was modeled using the H-formulation and rotated anisotropic resistivity to study its dynamic behavior, including inter-turn current sharing and intra-turn screening currents. Benchmarked by experimental data, the numerical model is able to study the defect tolerance of the SSMI coil during energization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Investigating the interconnectedness of active transportation and public transit usage as a primer for Mobility-as-a-Service adoption and deployment

With the advent of Mobility-as-a-Service packages to reduce car usage and (by extension) greenhouse gas emissions, it is crucial that researchers and practitioners consider mutual determinants and outcomes that link the adoption of multiple alternative modes. This study therefore investigates the joint usage of active travel (walking, cycling, bike-sharing) and public transit (bus, rail) modes with respect to personal and situational contexts. Online survey data (n = 826) were collected across six Midwestern U.S. states using Amazon MTurk. Respondents indicated their average weekly usage of eight travel modes across three trip purposes so that multimodality could be assessed. Several psychological constructs were extracted via (a) the stages of change framework, used to indicate willingness to adopt new behavior, and (b) confirmatory factor analysis conducted on Likert scales rooted in theories of community, identity, norms, personality, and well-being. A multiple-indicators multiple-causes structural equation model is then employed to investigate the process of adopting a modality style that incorporates active and transit modes. The model confirms that compatible physical and social contexts, as well as navigational skills and openness to learning, are key primers of multimodalism. However, a path juncture stemming from neighborhood support for mobility innovation illustrates a potential polarity in outcomes between individuals and communities. In addition, the stage of active mobility adoption is linked to identity and norm activation, offering further guidance on what influences readiness for change. The seamless integration of mobility services is critical to matching the convenience and comfort of the private vehicle; understanding potential pathways to sustainable mobility, though, requires analyses and interventions that are driven by well-being outcomes and grounded in rigorous psychological frameworks. Here, the research findings offer practical guidance for identifying intervention opportunities to be linked with MaaS enrollment while demonstrating the need to illuminate how mobility might relate to social cohesion, identity expression, and various sources of satisfaction.

42 ENGINEERING↗
DOE OSTI2025

Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM SciDAC) (Technical Final Report)

Runaway electrons can severely damage the plasma facing components on ITER during a major disruption and pose a major risk for tokamak fusion. It has been recognized that an adequate disruption mitigation system (DMS) is essential for the safe operation of ITER. The United States is responsible for the design and implementation of the disruption mitigation system on ITER, and in July 2016 the Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) was launched by DOE, in a joint Fusion Energy Sciences (FES) and Advanced Scientific Computing Research (ASCR) collaboration. SCREAM was a comprehensive theory and simulation SciDAC center that provided physics guidance in the avoidance and mitigation of runaway electrons, and in tandem with domestic and international experiments, helped establish the qualitative and quantitative bases for safe operational scenarios and viable mitigation techniques. The SCREAM center assembled a national team of experts in runaway electron physics, tokamak disruptions, magnetohydrodynamic (MHD) simulation, and advanced algorithms and computing. The team combined advanced simulation and analysis capability facilitated by direct participation of ASCR SciDAC institutes with theoretical models and code development by FES scientists to focus on the runaway risk for ITER and tokamaks in general. The research scope was focussed on integrated simulations of kinetic runaway electrons, including MHD and fluid models of impurity transport, within a research plan guided by theory. The specific research tasks were (1) establish the fundamental physics of runaway generation, saturation, and dynamical evolution in a tokamak; (2) examine the critical path toward runaway avoidance; and (3) investigate the viability and effectiveness of the leading candidate schemes for runaway mitigation. In all three areas, members of the team carried out scoping studies that established the readiness for rapid and critical advances, especially in the deployment and further development of large-to extreme-scale simulation tools. Our multi-pronged computational approach included (1) relativistic Fokker-Planck solvers with discretization in phase space, (2) self-consistent particle-in-cell techniques, (3) particle-based Monte-Carlo, and (4) MHD-particle hybrid simulations. Cross-check between these different methods provided an additional means for verification and further bolstered the fidelity of our physics prediction. Validation against experimental results brings confidence to the predictive capability for ITER and frequently leads to new ideas for understanding and mitigating the thermal quench driven runaway electron phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The physics basis to integrate an MHD stable, high-power hybrid scenario to a cool divertor for steady-state reactor operation

Abstract Coupling a high-performance core to a low heat flux divertor is a crucial step for ITER and a Fusion Pilot Plant or DEMO. Experiments in DIII-D recently expanded the steady-state hybrid scenario to high density and divertor impurity injection to study the feasibility of a radiating mantle solution. This work presents the physics basis for trade-offs between density, current drive and stability to tearing modes (TMs) at high β. EC power is crucial to tailor the plasma profiles into a passively stable state, and to eject impurities from the core. Off-axis EC depositions decrease the heating efficiency, but calculated electron heat transport coefficients show that this effect is partially mitigated by improved confinement inside the EC deposition. Additionally, the reduction in pressure is recovered by increasing the density. This favourable scaling of confinement with density was observed in high power plasmas for years, and this work provides a comprehensive explanation. ELITE predictions indicate that a path in peeling-ballooning stability opens up for certain conditions of density, power, q 95 and shaping, allowing the edge pressure to continue increasing without encountering a limit. In the core, calculated anomalous fast-ion diffusion coefficients are consistent with density fluctuation measurements in the toroidicity-induced Alfvén eigenmode range, showing that smaller fast-ion losses contribute to the enhanced confinement at high density. The edge integration study shows that divertor heat loads can be reduced with Ne and Ar injection, but this eventually triggers a cascade of n = 1, 2, 3 core TMs. We can now show that impurity radiation in the core is small and it is not the cause for the drop in confinement at high Ar and Ne injection rates. The overlap between the core TMs is consistent with the loss of pressure as estimated by the Belt model for the coupled rational surfaces. Optimization of these trade-offs has achieved plasmas with sustained H 98y2 = 1.7, f GW = 0.7 and ∼85% mantle radiation. The scenario and its variations at higher density and on- vs off-axis EC heating has been studied as a candidate for an integrated solution for several reactor designs, such as ITER, ARC, and the ARIES-ACT1 case, showing promising results in terms of fusion power and gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

SST-GPU: A Scalable SST GPU Component for Performance Modeling and Profiling

Programmable accelerators have become commonplace in modern computing systems. Advances in programming models and the availability of unprecedented amounts of data have created a space for massively parallel accelerators capable of maintaining context for thousands of concurrent threads resident on-chip. These threads are grouped and interleaved on a cycle-by-cycle basis among several massively parallel computing cores. One path for the design of future supercomputers relies on an ability to model the performance of these massively parallel cores at scale. The SST framework has been proven to scale up to run simulations containing tens of thousands of nodes. A previous report described the initial integration of the open-source, execution-driven GPU simulator, GPGPU-Sim, into the SST framework. This report discusses the results of the integration and how to use the new GPU component in SST. It also provides examples of what it can be used to analyze and a correlation study showing how closely the execution matches that of a Nvidia V100 GPU when running kernels and mini-apps.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Interconnect architecture for three-dimensional processing systems

A processing system includes a plurality of processor cores formed in a first layer of an integrated circuit device and a plurality of partitions of memory formed in one or more second layers of the integrated circuit device. The one or more second layers are deployed in a stacked configuration with the first layer. Each of the partitions is associated with a subset of the processor cores that have overlapping footprints with the partitions. The processing system also includes first memory paths between the processor cores and their corresponding subsets of partitions. The processing system further includes second memory paths between the processor cores and the partitions.

Jayasena, Nuwan S.↗
DOE OSTI2020

Panel Session 110B: US DOE High-Level Radioactive Waste (HLW) Interpretation

This panel focused on the US DOE High-Level radioactive Waste (HLW) definition interpretation. Discussion integrated topics such as the issuance of the HLW interpretation Federal Registration Notice, the status of the draft Environmental Assessment that analyzes the treatment and commercial disposal of up to 10,000-gallons of SRS Defense Waste Processing Facility recycle wastewater in South Carolina, and the path forward. Panelists with presentations: Status on the Environmental Assessment for the Commercial Disposal of Defense Waste Processing Facility Recycle Wastewater from SRS (Theresa Kliczewski); Enhancing Stakeholder Engagement on the HLW interpretation (Kara Colton, Rick McLeod); Technical Aspects of Potential Disposal of SRS DWPF Recycle Wastewater (Kent Rosenberger)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Cybersecurity Considerations and Research Pathways for Grid-Interactive Efficient Buildings

Federal facilities serve critical missions and functions that require safe, reliable, and efficient operations. Digitization of several facility operations has increased the cost-effectiveness of energy usage and optimization of energy system performance. As the building controls landscape shifts to become more connected and smarter, building operators now face unique opportunities and challenges to adopt smart enabled devices that can lower energy usage while also optimizing building system performance. The grid-interactive efficient buildings (GEB) initiative aims to make buildings cleaner and more flexible through these smart devices. Smart enabled devices allow greater connectivity and control through remote operations and provide crucial data for analytics and increased efficiency. GEBs enable demand flexibility that has the potential to reduce electrical costs and transform the grid edge where buildings connect to power grids. This operation of interconnected systems, if not designed with cybersecurity practices, causes security gaps and introduces potential attack paths by adversarial and non-adversarial entities leading to disruption of operations.

building controls↗
DOE OSTI2024

A collision-based hybrid method for the BGK equation

In this article, we apply the collision-based hybrid method introduced by Hauck and McClarren to the Boltzmann equation with the BGK operator and a hyperbolic scaling. An implicit treatment of the source term is used to handle stiffness associated with the BGK operator. Although it helps the numerical scheme become stable with a large time step size, it is still not obvious to achieve the desired order of accuracy due to the relationship between the size of the spatial cell and the mean free path. Without asymptotic preserving property, a very restricted grid size is required to resolve the mean free path, which is not practical. Our approaches are based on the noncollision-collision decomposition of the BGK equation. We introduce the arbitrary order of nodal discontinuous Galerkin (DG) discretization in space with a semi-implicit time-stepping method; we employ the backward Euler time integration for the uncollided equation and the 2nd order predictor-corrector scheme for the collided equation, i.e., both source terms in uncollided and collided equations are treated implicitly and only streaming term in the collided equation is solved explicitly. This improves the computational efficiency without the complexity of the numerical implementation. Numerical results are presented for various Knudsen numbers to present the effectiveness and accuracy of our hybrid method. Also, we compare the solutions of the hybrid and non-hybrid schemes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Alfalfa Virtual Building Service: Software Engineering Best Practices Applied to Runtime Interaction with Building Energy Models

Buildings are active participants in increasingly complex energy systems. Building Energy Modeling (BEM) has a key role to play in planning and de-risking an equitable energy transition, with BEM-backed "virtual buildings" critical path for diverse applications that include workforce training tools, Hardware-in-the-Loop (HIL) experimentation to study equipment performance under a range of conditions, Control-Hardware-in-the-Loop (CHIL) experimentation to de-risk commercial control implementations at equipment through grid orchestration levels, and integration of dynamic load profiles into grid modeling tools for energy system experimentation at the urban scale. Modeling requirements vary across these applications, but many software engineering tasks do not. The Alfalfa Virtual Building Service (AVBS, see https://github.com/NREL/alfalfa/wiki) is an open-source web service that solves these common tasks robustly in one place, providing a foundational platform for power users to bootstrap their own applications. AVBS abstracts the specifics of runtime interaction with OpenStudio, Modelica, and Spawn of EnergyPlus models behind a unified REST API. Additionally, AVBS provides resources for cloud deployment and scaling to 100s of parallel simulations, a growing library of modular Operational Technology (OT) integrations for emulation of real-world interfaces, and scripts to automate the population of communities of virtual buildings from URBANopt, ResStock and ComStock.

building automation↗
DOE OSTI2022

Modeling Needs to Support the Reconstruction Strategy of Ukraine

Ukraine has been facing difficult times. Russia started a full-scale invasion on February 24, 2022, and has caused huge human and infrastructure losses. As a result of the invasion, more than 8 million people left their homes. The Russian army has targeted Ukrainian energy sources, destroying or damaging over 50% of thermal power generation capacity (coal and gas), 30% of solar generation, and 90% of wind generation. The Russians are also threatening Ukrainian energy security, as they currently control the Zaporizhzhia nuclear power plant – Europe’s largest. In addition, Russia has destroyed 121,000 buildings, including 5,300 multistory buildings - more buildings than Ukraine built during the past five years. The occupiers destroyed or damaged 343 district heat boiler stations (destroyed 12 / damaged 331) and 8 CHPs (destroyed 4/ damaged 4). There is no doubt that Ukraine will be able to protect its sovereignty and restore its territorial integrity; however, this will be a long process. The Government of Ukraine needs to simultaneously focus on many issues, from waging war to the restoration of electricity, water, and heat supply in the territories close to the front line. In addition to these time-sensitive issues, the government also needs to plan a long term multi-sector strategy for the reconstruction of Ukraine’s economy. Ukraine initiated a broad-based political process for recovery at the International Ukraine Recovery Conference (URC 2022) in Lugano, Switzerland. Prior to URC 2022, Ukraine applied for EU membership on February 28 and was granted the status of candidate country by the European Council on June 23, 2022. More than 40 nations are fully committed to supporting Ukraine throughout its path from early to long-term recovery. In addition, these countries link a successful recovery to Ukraine’s European future (URC, 2022a). The reconstruction process will be long. The government also recognizes that it is not enough to build energy infrastructure back – it should be built back better. The government is committed to building a more sustainable and cleaner energy sector. Energy modeling can provide important insights into paths to decarbonize Ukraine’s economy while rebuilding it after the war. Pacific Northwest National Laboratory (PNNL) and the Institute for Economics and Forecasting of the National Academy of Sciences of Ukraine (IEF) work on modeling decarbonization scenarios in the energy sector of Ukraine. The teams focus on the heat sector, which includes district heat, autonomous and individual heating, and energy consumption by final energy consumption sectors, including buildings. Since the teams use world-class integrated assessment models, other sectors will be analyzed as well. The purpose of the document is to describe the needs of key Ukrainian ministries in modeling the reconstruction process. The scope of this document is limited to modeling energy consumption by the power and heat sectors, final energy consumption by buildings, and GHG emissions from these two sectors. The rest of this policy memo is organized as follows. After this introduction, Chapter 2 discusses key provisions of the reconstruction strategy with a focus on energy, buildings, heat supply, environment, and European integration aspirations. Section 3 provides an overview of recent developments and key decisions of the Government of Ukraine in the area of energy efficiency and decarbonization. Finally, Section 4 provides key suggestions on how PNNL and IEF can help the Government of Ukraine in developing strategic documents.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Charge sharing in pixelated semiconductor sensors

The charge sharing between neighboring pixels in pixelated sensors can be used to measure particle or x-ray coordinates with accuracy better than the pixel pitch. The accurate model of the charge distribution shape is essential to achieve ultimate coordinate accuracy. The charge sharing is caused by charge carriers' diffusion on the path from the generation point to pixels. This paper is focused on the diffusion of the initially compact charge cloud in the field free region. The diffusion equation solutions are obtained using separation of variable and Fourier synthesis method for different initial conditions and resulting charge distributions are integrated over pixel areas. In conclusion, the look up table containing pre-calculated values for pixel charge fractions is proposed to speed up numerical calculations.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Advancing density functional tight-binding method for large organic molecules through equivariant neural networks

Semi-empirical quantum-mechanical (QM) methods have become valuable tools for studying complex (bio)molecular systems due to their balance between computational efficiency and accuracy. A key aspect of these methods is their parameterization, which not only governs the reliability of the results but also provides an opportunity to enhance their overall performance. In our previous work [J. Phys. Chem. Lett., 2021, 11, 16], we advanced the third-order semi-empirical density functional tight-binding (DFTB3) method for computing multiple properties of small molecules by developing the machine learning (ML) potential NN rep to bridge the gap between DFTB3 electronic components and those of the hybrid DFT-PBE0 functional. To overcome the limitations of NN rep , we introduce the EquiDTB framework, which leverages physics-inspired equivariant neural networks (NN) to parameterize scalable and transferable many-body Δ TB potentials, replacing the standard pairwise DFTB repulsive potential. This advancement extends the applicability of our ML-corrected DFTB approach to larger molecules and non-covalent systems (including only C, N, O, and H atoms), going beyond the chemical space represented in the training QM datasets. The enhanced performance of EquiDTB over the standard TB methods is demonstrated by the accurate computation of the atomic forces of S66x8 molecular dimers, as well as their interaction energies. Moreover, EquiDTB can be effectively employed to explore the potential energy surfaces of large and flexible drug-like molecules—for example, to determine the minimum energy path between isomers, analyze structural transitions during dynamical simulations, compute vibrational modes, and investigate energetic rankings. The performance for single molecules slightly decreases when the DFTB electronic energy is reduced to first-order but remains superior to standard TB methods. Our work thus demonstrates that an optimal integration of an equivariant NN with QM datasets can advance the DFTB method while maintaining high efficiency, paving the way for reliable (bio)molecular simulations.

Medrano Sandonas, Leonardo [Technische Universität↗