SEARCH · Engineering Papers

Results for “regularity”

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.

SourceResearch connectionFrom yearThrough year

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

At least 811 records · Page 45

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Tunable Conversion of Poly(Acrylic Acid) to Vinyl Methyl Ether Copolymers via Electrochemical Decarboxylation

Recognizing the potential of polymer-to-polymer transformations for plastics valorization, a new electrochemical strategy for the decarboxylative functionalization of polyacrylates has been developed. In contrast to most electrochemical approaches to plastic waste valorization that use small-molecule mediators or only form deconstruction products, this report describes one of the first direct electrochemical redox upconversions of a commodity polymer and involves activation of poly(acrylic acid), a material commonly found in disposable diapers and other personal hygiene products. Via electrochemical decarboxylation, poly(acrylic acid) is converted to methyl ether-containing copolymers. Although methods to prepare the 1:1 regularly alternating poly(vinyl methyl ether-alt-maleic anhydride) are known, routes to generate copolymers with variable levels of vinyl methyl ether and acrylate moieties are unknown, but such polymers may have utility in similar applications, such as dentifrices. Tunable conversions with up to 62 mol% vinyl methyl ether content were achieved via controlled-current electrolyses at >100 mA. This approach enables efficient access to copolymers with broad range of comonomer content from a common waste polymer starting material.

Dauzvardis, Fabian [University of Delaware, Newark↗

DOE OSTI2026

Semi-Dynamic Leach Testing of Densified Silicon-based Iodine Waste Forms

Iodine waste forms (IWF) require a conceptual corrosion release model (CCRM) to provide iodine (I) release rates for performance modeling nuclear waste disposal repositories. To develop a CCRM, an understanding of the corrosion mechanisms of the IWF are required along with data from consistent test methods to parameterize the model. The present study has advanced both areas by providing and assessing the corrosion resistance of IWF types based on their processing history in minor variations of semi-dynamic leach tests. The test included a series of semi-dynamic leach tests using monolithic IWFs in deionized water (leachant). Several experiments were conducted under alternate test conditions with changes to temperature, leachant replacement, leachant pH, leachant volume, masking, and surface finish to elucidate if varying these conditions impacted IWF corrosion behavior. Tests were conducted on two classes of IWFs: (1) I-bearing silver-mordenite (AgZ) materials processed by hot isostatic pressing (HIP) at different temperatures, pressures, sizes, and times; and (2) I-bearing silver-functionalized silica aerogels (SFA) processed by either HIP or spark plasma sintering (SPS). The corrosion susceptibility of AgZ samples was influenced by HIP temperature and pressure. The SPS SFAs retained I far better than HIP SFAs. Additional findings in this study include: (1) The iodine dissolution rate decreased with decreasing temperature, (2) a common ion effect may occur and slow dissolution of the host phase if the leachant is not regularly replaced, (3) pH controls the dissolution rate, and (4) the iodine dissolution rate slows with extended test time (up to 224 days). Based on this work, these parameters should thus be represented when developing a CCRM.

corrosion↗

DOE OSTI2020

Gene editing and CRISPR in the clinic: current and future perspectives

Genome editing technologies, particularly those based on zinc-finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs), and CRISPR (clustered regularly interspaced short palindromic repeat DNA sequences)/Cas9 are rapidly progressing into clinical trials. Most clinical use of CRISPR to date has focused on ex vivo gene editing of cells followed by their re-introduction back into the patient. The ex vivo editing approach is highly effective for many disease states, including cancers and sickle cell disease, but ideally genome editing would also be applied to diseases which require cell modification in vivo. However, in vivo use of CRISPR technologies can be confounded by problems such as off-target editing, inefficient or off-target delivery, and stimulation of counterproductive immune responses. Current research addressing these issues may provide new opportunities for use of CRISPR in the clinical space. In this review, we examine the current status and scientific basis of clinical trials featuring ZFNs, TALENs, and CRISPR-based genome editing, the known limitations of CRISPR use in humans, and the rapidly developing CRISPR engineering space that should lay the groundwork for further translation to clinical application.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Numerical simulation of parametric resonance in point absorbers using a simplified model

Abstract Parametric resonance is a non‐linear phenomenon in which a system can oscillate at a frequency different from its exciting frequency. Some wave energy converters are prone to this phenomenon, which is usually detrimental to their performance. Here, a computationally efficient way of simulating parametric resonance in point absorbers is presented. The model is based on linear potential theory, so the wave forces are evaluated at the mean position of the body. However, the first‐order variation of the body's centres of gravity and buoyancy is taken into account. This gives essentially the same result as a more rigorous approach of keeping terms in the equation of motion up to second order in the body motions. The only difference from a linear model is the presence of non‐zero off‐diagonal elements in the mass matrix. The model is benchmarked against state‐of‐the‐art non‐linear Froude–Krylov and computational fluid dynamics models for free decay, regular wave, and focused wave group cases. It is shown that the simplified model is able to simulate parametric resonance in pitch to a reasonable accuracy even though no non‐linear wave forces are included. The simulation speed on a standard computer is up to two orders of magnitude faster than real time.

16 TIDAL AND WAVE POWER↗

DOE OSTI2022

The throughput calibration of the VERITAS telescopes

Context. The response of imaging atmospheric Cherenkov telescopes to incident γ -ray-initiated showers in the atmosphere changes as the telescopes age due to exposure to light and weather. These aging processes affect the reconstructed energies of the events and γ -ray fluxes. Aims. This work discusses the implementation of signal calibration methods for the Very Energetic Radiation Imaging Telescope Array System (VERITAS) to account for changes in the optical throughput and detector performance over time. Methods. The total throughput of a Cherenkov telescope is the product of camera-dependent factors, such as the photomultiplier tube gains and their quantum efficiencies, and the mirror reflectivity and Winston cone response to incoming radiation. This document summarizes different methods to determine how the camera gains and mirror reflectivity have evolved over time and how we can calibrate this changing throughput in reconstruction pipelines for imaging atmospheric Cherenkov telescopes. The implementation is validated against seven years of observations with the VERITAS telescopes of the Crab Nebula, which is a reference object in very-high-energy astronomy. Results. Regular optical throughput monitoring and the corresponding signal calibrations are found to be critical for the reconstruction of extensive air shower images. The proposed implementation is applied as a correction to the signals of the photomultiplier tubes in the telescope simulation to produce fine-tuned instrument response functions. This method is shown to be effective for calibrating the acquired γ -ray data and for recovering the correct energy of the events and photon fluxes. At the same time, it keeps the computational effort of generating Monte Carlo simulations for instrument response functions affordably low.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

BeyondPlanck X. Bandpass and beam leakage corrections

We discuss the treatment of bandpass and beam leakage corrections in the Bayesian BeyondPlanck CMB analysis pipeline as applied to the Planck LFI measurements. As a preparatory step, we first apply three corrections to the nominal LFI bandpass profiles including removal of a known systematic effect in the ground measuring equipment at 61 GHz; smoothing of standing wave ripples; and edge regularization. The main net impact of these modifications is an overall shift in the 70 GHz bandpass of +0.6 GHz; we argue that any analysis of LFI data products, either from Planck or BeyondPlanck, should use these new bandpasses. In addition, we fit a single free bandpass parameter for each radiometer of the form Δ i = Δ 0 + $δ$ i , where Δ 0 represents an absolute frequency shift per frequency band and $δ$i is a relative shift per detector. The absolute correction is only fitted at 30 GHz with a full $\chi$ 2 -based likelihood, resulting in a correction of Δ 30 = 0.24 ± 0.03 GHz. The relative corrections are fitted using a spurious map approach, fundamentally similar to the method pioneered by the WMAP team, but without introducing many additional degrees of freedom. All bandpass parameters are sampled using a standard Metropolis sampler within the main BeyondPlanck Gibbs chain, and bandpass uncertainties are thus propagated to all other data products in the analysis. In total, we find that our bandpass model significantly reduces leakage effects. For beam leakage corrections, we adopt the official Planck LFI beam estimates without additional degrees of freedom, and only marginalize over the underlying sky model. We note that this is the first-time leakage from beam mismatch has been included for Planck LFI maps.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Galaxy cluster matter profiles - I. Self-similarity, mass calibration, and observable-mass relation validation employing cluster mass posteriors

We present a study of the weak lensing inferred matter profiles ΔΣ(R) of 698 South Pole Telescope (SPT) thermal Sunyaev-Zel’dovich effect (tSZE) selected and MCMF optically confirmed galaxy clusters in the redshift range 0.25 < z < 0.94 that have associated weak gravitational lensing shear profiles from the Dark Energy Survey (DES). Rescaling these profiles to account for the mass dependent size and the redshift dependent density produces average rescaled matter profiles ΔΣ(R/R200c)/(ρcritR200c) with a lower dispersion than the unscaled ΔΣ(R) versions, indicating a significant degree of self-similarity. Galaxy clusters from hydrodynamical simulations also exhibit matter profiles that suggest a high degree of self-similarity, with RMS variation among the average rescaled matter profiles with redshift and mass falling by a factor of approximately six and 23, respectively, compared to the unscaled average matter profiles. We employed this regularity in a new Bayesian method for weak lensing mass calibration that employs the so-called cluster mass posterior P(M200|ζ̂, λ̂, z), which describes the individual cluster masses given their tSZE (ζ̂) and optical (λ̂, z) observables. This method enables simultaneous constraints on richness λ-mass and tSZE detection significance ζ-mass relations using average rescaled cluster matter profiles. We validated the method using realistic mock datasets and present observable-mass relation constraints for the SPT×DES sample, where we constrained the amplitude, mass trend, redshift trend, and intrinsic scatter. Our observable-mass relation results are in agreement with the mass calibration derived from the recent cosmological analysis of the SPT×DES data based on a cluster-by-cluster lensing calibration. Our new mass calibration technique offers a higher efficiency when compared to the single cluster calibration technique. We present new validation tests of the observable-mass relation that indicate the underlying power-law form and scatter are adequate to describe the real cluster sample but that also suggest a redshift variation in the intrinsic scatter of the λ-mass relation may offer a better description. In addition, the average rescaled matter profiles offer high signal-to-noise ratio (S/N) constraints on the shape of real cluster matter profiles, which are in good agreement with available hydrodynamical ΛCDM simulations. This high S/N profile contains information about baryon feedback, the collisional nature of dark matter, and potential deviations from general relativity.Key words: gravitational lensing: weak / galaxies: clusters: general / large-scale structure of Universe

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Quantification of Dynamic Scattering Effects in Molecular Crystals using Large Angle Rocking Beam Electron Diffraction

Electron crystallography provides a pathway to solve structure of small crystals (< 1um) in size, and thus overcomes difficult synthesis constraints involved in growing large crystals. Generally, electron diffraction data is collected in the form of integrated intensity using continuous rotation or by precession. The measured intensities in 3D are utilized for structure solution. Using this approach, structure solution of difficult crystals, such as small crystals of zeolites, metal-organic frameworks, molecular crystals, and proteins , can be solved by electron diffraction. However, electron structure solutions are regularly reported with higher Rvalues than x-ray or neutron diffraction. While similar structures are found despite the high R-values, the differences in the measure intensity and theory calculated intensity limit information that can be extracted by electron diffraction. Previous work demonstrated that including multiple-scattering effects significantly reduce the R-values . Thus, it is critical to be able to quantify the dynamical diffraction effects in molecular crystals.

Crystallography↗

DOE OSTI2020

Automatic log analysis with NLP for the CMS workflow handling

The central Monte-Carlo production of the CMS experiment utilizes the WLCG infrastructure and manages daily thousands of tasks, each up to thousands of jobs. The distributed computing system is bound to sustain a certain rate of failures of various types, which are currently handled by computing operators a posteriori. Within the context of computing operations, and operation intelligence, we propose a Machine Learning technique to learn from the operators with a view to reduce the operational workload and delays. This work is in continuation of CMS work on operation intelligence to try and reach accurate predictions with Machine Learning. We present an approach to consider the log files of the workflows as regular text to leverage modern techniques from Natural Language Processing (NLP). In general, log files contain a substantial amount of text that is not human language. Therefore, different log parsing approaches are studied in order to map the log files’ words to high dimensional vectors. These vectors are then exploited as feature space to train a model that predicts the action that the operator has to take. This approach has the advantage that the information of the log files is extracted automatically and the format of the logs can be arbitrary. In this work the performance of the log file analysis with NLP is presented and compared to previous approaches.

Layer, Lukas↗

DOE OSTI2024

HPC resources for CMS offline computing: An integration and scalability challenge for the Submission Infrastructure

The computing resource needs of LHC experiments are expected to continue growing significantly during the Run 3 and into the HL-LHC era. The landscape of available resources will also evolve, as High Performance Computing (HPC) and Cloud resources will provide a comparable, or even dominant, fraction of the total compute capacity. The future years present a challenge for the experiments’ resource provisioning models, both in terms of scalability and increasing complexity. The CMS Submission Infrastructure (SI) provisions computing resources for CMS workflows. This infrastructure is built on a set of federated HTCondor pools, currently aggregating 400k CPU cores distributed worldwide and supporting the simultaneous execution of over 200k computing tasks. Incorporating HPC resources into CMS computing represents firstly an integration challenge, as HPC centers are much more diverse compared to Grid sites. Secondly, evolving the present SI, dimensioned to harness the current CMS computing capacity, to reach the resource scales required for the HLLHC phase, while maintaining global flexibility and efficiency, will represent an additional challenge for the SI. To preventively address future potential scalability limits, the SI team regularly runs tests to explore the maximum reach of our infrastructure. In this note, the integration of HPC resources into CMS offline computing is summarized, the potential concerns for the SI derived from the increased scale of operations are described, and the most recent results of scalability test on the CMS SI are reported.

Pérez-Calero Yzquierdo, Antonio↗

DOE OSTI2024

Repurposing of the Run 2 CMS High Level Trigger Infrastructure as a Cloud Resource for Offline Computing

The former CMS Run 2 High Level Trigger (HLT) farm is one of the largest contributors to CMS compute resources, providing about 25k job slots for offline computing. This CPU farm was initially employed as an opportunistic resource, exploited during inter-fill periods, in the LHC Run 2. Since then, it has become a nearly transparent extension of the CMS capacity at CERN, being located on-site at the LHC interaction point 5 (P5), where the CMS detector is installed. This resource has been configured to support the execution of critical CMS tasks, such as prompt detector data reconstruction. It can therefore be used in combination with the dedicated Tier 0 capacity at CERN, in order to process and absorb peaks in the stream of data coming from the CMS detector. The initial configuration for this resource, based on statically configured VMs, provided the required level of functionality. However, regular operations of this cluster revealed certain limitations compared to the resource provisioning and use model employed in the case of WLCG sites. A new configuration, based on a vacuum-like model, has been implemented for this resource in order to solve the detected shortcomings. This paper reports about this redeployment work on the permanent cloud for an enhanced support to CMS offline computing, comparing the former and new models’ respective functionalities, along with the commissioning effort for the new setup.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2024

The ATLAS experiment software on ARM

With an increased dataset obtained during the Run 3 of the LHC at CERN and the even larger expected increase of the dataset by more than one order of magnitude for the HL-LHC, the ATLAS experiment is reaching the limits of the current data processing model in terms of traditional CPU resources based on x86_64 architectures and an extensive program for software upgrades towards the HL-LHC has been set up. The ARM architecture is becoming a competitive and energy efficient alternative. Some surveys indicate its increased presence in HPCs and commercial clouds, and some WLCG sites have expressed their interest. Chip makers are also developing their next generation solutions on ARM architectures, sometimes combining ARM and GPU processors in the same chip. Consequently it is important that the ATLAS software embraces the change and is able to successfully exploit this architecture. We report on the successful porting to ARM of the Athena software framework, which is used by ATLAS for both online and offline computing operations. Furthermore we report on the successful validation of simulation workflows running on ARM resources. For this we have set up an ATLAS Grid site using ARM compatible middleware and containers on Amazon Web Services (AWS) ARM resources. The ARM version of Athena is fully integrated in the regular software build system and distributed in the same way as other software releases. In addition, the workflows have been integrated into the HEPscore benchmark suite which is the planned WLCG wide replacement of the HepSpec06 benchmark used for Grid site pledges. In the overall porting process we have used resources on AWS, Google Cloud Platform (GCP) and CERN. A performance comparison of different architectures and resources will be discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

An asymptotically compatible approach for Neumann-type boundary condition on nonlocal problems

In this paper we consider 2D nonlocal diffusion models with a finite nonlocal horizon parameter δ characterizing the range of nonlocal interactions, and consider the treatment of Neumann-like boundary conditions that have proven challenging for discretizations of nonlocal models. We propose a new generalization of classical local Neumann conditions by converting the local flux to a correction term in the nonlocal model, which provides an estimate for the nonlocal interactions of each point with points outside the domain. While existing 2D nonlocal flux boundary conditions have been shown to exhibit at most first order convergence to the local counter part as δ → 0, the proposed Neumann-type boundary formulation recovers the local case as O(δ 2 ) in the L∞(Ω) norm, which is optimal considering the O(δ 2 ) convergence of the nonlocal equation to its local limit away from the boundary. We analyze the application of this new boundary treatment to the nonlocal diffusion problem, and present conditions under which the solution of the nonlocal boundary value problem converges to the solution of the corresponding local Neumann problem as the horizon is reduced. To demonstrate the applicability of this nonlocal flux boundary condition to more complicated scenarios, we extend the approach to less regular domains, numerically verifying that we preserve second-order convergence for non-convex domains with corners. Finally, based on the new formulation for nonlocal boundary condition, we develop an asymptotically compatible meshfree discretization, obtaining a solution to the nonlocal diffusion equation with mixed boundary conditions that converges with O(δ 2 ) convergence.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Automatic Lane-Level Road Network Extraction from Aerial Imagery for Transportation Digital Twins

Accurate road networks are essential for credible traffic microsimulation and transportation digital twins, yet high-definition maps are often difficult to obtain due to limited availability, high cost, or proprietary restrictions. Some build networks from crowdsourced data, such as OpenStreetMap, but these sources often contain geometric and semantic inconsistencies. Others create networks manually, a process that is labor-intensive and difficult to scale. To address these limitations, this work presents an end-to-end pipeline that automatically extracts georeferenced, lane-level road networks from publicly available high-resolution satellite imagery and converts them into simulation-ready assets. The developed end-to-end pipeline has three primary modules: (1) A computer-vision-based module first detects directed lane geometries and intersection layouts. (2) A heuristic-based topology construction module then identifies approach and exit legs and establishes conflict-free lane-to-lane connections. (3) Finally, an automatic simulation-building module converts the extracted network into standard formats, e.g., OpenDRIVE, and generates routable SUMO networks. The framework supports both complete network construction from scratch and local-scale refinement of existing networks through lane-count correction, transition recovery, and geometric regularization. The proposed pipeline provides a practical pathway to generate traffic simulation networks from satellite imagery, significantly reducing manual reconstruction effort and enabling scalable, continuously updated transportation digital twins.

Guo, Hetian [University of Georgia, Athens] (ORCID↗

DOE OSTI2024

Computing in AEC Education: Hindsight, Insight, and Foresight

In the architecture, engineering, and construction (AEC) fields, computing and information technologies play an increasingly prevalent and complex role in day-to-day work. Consequently, educators must adjust and, in many cases, reimagine curricula and teaching methodologies to adapt to the changing landscape. Past research efforts led by the ASCE Computing Division Education Committee, formerly called the Task Committee on Computing Education of the Technical Council on Computing and Information Technology (TCCIT), have regularly surveyed AEC educators to understand computing trends in AEC curricula, with the latest survey taking place nearly a decade ago. This work presents the results of an updated survey that used this prior work as a springboard, providing timely insights into the computing skills and curricular barriers faced by AEC educators today. The results showed that the technical skills used by students have evolved, but the barriers faced in incorporating new skills into curricula have remained largely the same. In addition to comparisons with prior surveys, this work presents the results of an expanded, open-ended portion of the survey that explores educator perspectives on the future of the AEC workforce in a broader lens than used in previous surveys. Thematic analysis of these open-ended responses revealed themes that were common among responses and provided organization to the findings. For example, educators provided their vision of what competencies the future AEC workforce would need, which were thematically organized into a continuum based on the level of interaction between humans and technology. These results suggest an increasingly complex and evolving relationship between the AEC workforce and emerging technology, highlighting the need for educators to encourage the development of technological adaptability and agility. Overall, this work provides a systematic comparison of current educational practices in AEC computing with a decade ago to illustrate educational shifts and adds a prediction of AEC trends from experts in AEC education, providing crucial discussion of curricular transformations that will better position students for success in the workforce.

42 ENGINEERING↗

DOE OSTI2025

X-Ray Microtomography Measurements of Coevolving Particle Size and Shape

This study measures simultaneous changes in particle size and shape during grain crushing and proposes quantitative metrics to track their coevolving statistical distribution. Idealized granular materials with widely different initial particle shapes (i.e., spheres, quasi-ellipsoids, plates, and rods) are compressed oedometrically at vertical stress levels sufficient to induce particle breakage. X-ray microtomography and digital image analysis are then used to track changes in particle-scale characteristics. The observations indicate that more irregular particles tend to exhibit a higher degree of breakage compared to more spherical grains subjected to an equivalent stress level. Despite these differences, the measurements also indicate that all the tested particle sets approach a common attractor (expressed in terms of aspect ratio, flatness, and elongation), regardless of their initial morphology. A new shape evolution index is proposed to quantify such trends. With this new metric, it is shown that in all the tested materials, the particle shape evolves faster than size, especially for the more regular morphologies, thus implying that the shape tends to reach its attractor point before the grading reaches its ultimate particle size distribution. Furthermore, this finding underpins an ultimate stage of crushing-driven compression that occurs with particle size reduction but minimal shape changes (i.e., self-replicating shapes due to fractures).

Oedometric compression↗

DOE OSTI2020

Extension of frozen natural orbital approximation to open-shell references: Theory, implementation, and application to single-molecule magnets

Natural orbitals are often used to achieve a more compact representation of correlated wave-functions. Using natural orbitals computed as eigenstates of the virtual–virtual block of the state density matrix instead of the canonical Hartree–Fock orbitals results in smaller errors when the same fraction of virtual space is frozen. This strategy, termed frozen natural orbital (FNO) approach, is effective in reducing the cost of regular coupled-cluster (CC) calculations and some multistate methods, such as EOM-IP-CC (equation-of-motion CC for ionization potentials). This contribution extends the FNO approach to the EOM-SF-CC ansatz (EOM-CC with spin-flip). In contrast to EOM-IP-CCSD, EOM-SF-CCSD relies on high-spin open-shell references. Using FNOs computed for an open-shell reference leads to an erratic behavior of the EOM-SF-CC energies and properties due to an inconsistent truncation of the α and β orbital spaces. A general solution to problems arising in the EOM-CC calculations utilizing open-shell references, termed OSFNO (open-shell FNO), is proposed. Furthermore, by means of singular value decomposition (SVD) of the overlap matrix of the α and β orbitals, the OSFNO algorithm identifies the corresponding orbitals and determines virtual orbitals corresponding to the singly occupied space. This is followed by SVD of the singlet part of the state density matrix in the remaining virtual orbital subspace. The so-computed FNOs preserve the spin purity of the open-shell orbital subspace to the extent allowed by the original reference, thus facilitating a safe truncation of the virtual space. The performance of OSFNO is benchmarked for selected diradicals and triradicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation

It is shown how the electronic equations of motion in extended Lagrangian Born–Oppenheimer molecular dynamics simulations can be integrated using low-rank approximations of the inverse Jacobian kernel. This kernel determines the metric tensor in the harmonic oscillator extension of the Lagrangian that drives the evolution of the electronic degrees of freedom. The proposed kernel approximation is derived from a pseudoinverse of a low-rank estimate of the Jacobian, which is expressed in terms of a generalized set of directional derivatives with directions that are given from a Krylov subspace approximation. The approach allows a tunable and adaptive approximation that can take advantage of efficient preconditioning techniques. The proposed kernel approximation for the integration of the electronic equations of motion makes it possible to apply extended Lagrangian first-principles molecular dynamics simulations to a broader range of problems, including reactive chemical systems with numerically sensitive and unsteady charge solutions. This can be achieved without requiring exact full calculations of the inverse Jacobian kernel in each time step or relying on iterative non-linear self-consistent field optimization of the electronic ground state prior to the force evaluations as in regular direct Born–Oppenheimer molecular dynamics. We note the low-rank approximation of the Jacobian is directly related to Broyden’s class of quasi-Newton algorithms and Jacobian-free Newton–Krylov methods and provides a complementary formulation for the solution of nonlinear systems of equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗