Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo method”

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 235 records · Page 13

rmc-discord : reverse Monte Carlo refinement of diffuse scattering and correlated disorder from single crystals

A user-friendly program has been developed to analyze diffuse scattering from single crystals with the reverse Monte Carlo method. The approach allows for refinement of correlated disorder from atomistic supercells with magnetic or structural (occupational and/or displacive) disorder. The program is written in Python and optimized for performance and efficiency. Refinements of two user cases obtained with legacy neutron-scattering data demonstrate the effectiveness of the approach and the developed program. It is shown with bixbyite, a naturally occurring magnetic mineral, that the calculated three-dimensional spin-pair correlations are resolved with finer real-space resolution compared with the pair distribution function calculated directly from the reciprocal-space pattern. With the triangular lattice Ba 3 Co 2 O 6 (CO 3 ) 0.7 , refinements of occupational and displacive disorder are combined to extract the one-dimensional intra-chain correlations of carbonate molecules that move toward neighboring vacant sites to accommodate strain induced by electrostatic interactions. The program is packaged with a graphical user interface and extensible to serve the needs of single-crystal diffractometer instruments that collect diffuse-scattering data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Effects of non-equilibrium ionization and excitation on radiation absorption in plasma plumes induced by ablation of metal targets with nanosecond laser pulses

Ionization and radiation absorption in nanosecond laser-induced plasma plumes are routinely modeled using the Saha–Boltzmann equilibrium ionization model (EQM). However, the equilibrium assumption can be inaccurate during the laser pulse when non-equilibrium effects significantly impact radiation absorption. In the present work, the EQM and non-equilibrium collisional-radiative model (CRM) are compared to reveal the effect of plasma non-equilibrium on radiation absorption in non-homogeneous plumes and degree of plasma shielding. Simulations of plume expansion induced by irradiation of a copper target in 1 atm argon background gas with a 10 ns Gaussian pulse at a fluence from 8 Jcm −2 to 14 Jcm −2 are performed with a hybrid computational model that couples a lumped particle direct simulation Monte Carlo method with either CRM or EQM. The simulations show that the EQM strongly underestimates the effects of ionization and radiation absorption compared to CRM and, contrary to the CRM, predict strong ionization of the background gas. The differences between the models are explained by the qualitatively different coupling between plume expansion and dynamics of ionization and exitation processes in the CRM and EQM under conditions when the characteristic times for most radiation- and electron-induced processes are longer than the pulse duration. The CRM-based predictions are also found to agree much better with available experimental data. In conclusion, these results indicate that the model of Saha–Boltzmann equilibrium cannot be used for reliable prediction of the degree of plasma shielding in plumes induced by nanosecond laser pulses or for processing results of spectroscopic measurements at early stages of expansion of such plumes.

Physics - Plasma physics↗

Toward a Better Understanding of Ni Coarsening in Solid Oxide Cells: NiH on Ni (111) Examined Using a Combined Theoretical Approach

The coarsening of the Ni particles in the hydrogen electrode of solid oxide cells (SOCs) is an important degradation mechanism. Here, in this paper, density-functional theory and kinetic Monte Carlo methods are used to explore our recent hypothesis that the surface diffusion of NiH may cause faster Ni coarsening in electrolysis cell mode under an overpotential. Using both methods, the diffusion constant or diffusivity of NiH on Ni (111) is determined as the product of the surface coverage and single-molecule diffusivity for the first time considering all possible diffusion paths. It is then determined versus overpotential at the triple-phase boundary of the hydrogen electrode assuming a typical operating temperature of the SOC. Under a significant overpotential, the diffusivity of NiH is found to be sufficiently large to support the above hypothesis that NiH may promote Ni coarsening. However, based on the adsorption configurations identified, the dissociation and reformation of NiH on Ni (111) could occur. Thus, more work is needed to develop a model of Ni coarsening considering both molecular and dissociated forms of NiH.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electronic Structure Theory and Novel Materials

This grant supported research on electronic structure and materials theory, with focus on three main issues: (i) novel techniques to deal with correlation in the electronic ground-state, (ii) topological materials, (iii) the phase diagram of lattice spin models. Regarding (i), we applied to the homogeneous electron liquid an approach that we previously developed in the context of molecular systems. In this scheme the electronic occupation probabilities and the natural spin orbitals are used to construct an approximate two-body density matrix for the electronic ground-state. Regarding (ii) we used standard electronic structure methods based on density functional theory to model topological materials and interpret experimental observations. Finally, regarding (iii) we further developed a numerical approach to compute the renormalized couplings within real space renormalization group theory in the context of lattice spin models. The main findings were the following. (i) We found that with our approximate two-body density matrix, which works well for small molecules, is not sufficiently accurate for condensed phase systems. Missing a systematic way of improving on the adopted approximations, we decided not to pursue this approach. (ii) We performed two studies. In one, we investigated the influence of Te defects on the topological properties of a WTe2 monolayer, finding that while Te vacancies, even in modest concentration, destroy the topological character, Te adatoms do not, consistent with a recent experiment. In another study, we predicted Weyl semimetal character and strong anomalous Hall effect in the Heusler compensated ferrimagnet Ti2MnAl. (iii) We developed a new Monte Carlo method to do real space renormalization group calculations for lattice spin models. We subsequently extended the scheme to deal with lattice spin models in presence of quenched disorder, finding that the approach can distinguish systems with finite and strong disorder. In the finite disorder case, the method allows one to find with good approximation the critical coupling distribution and the critical exponents.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Using Best Basis Inventory Data to Direct Strategies for Real Time Monitoring of Hanford High Level Waste – 26226

The potential to accelerate the processing of low- and high-level tank waste by applying real-time monitoring (RTM) of chemical and physical properties has prompted research into the suitability of multiple analytical methods for that purpose. The broad variety of waste stream properties and the large number of analytes of interest (as evidenced by Waste Acceptance Criteria (WAC) and Process Control Limit (PCL) lists) lead to an overwhelming set of possible analytical scenarios. This report describes the use of Best Basis Inventory (BBI) data to find the most relevant analytical targets for the specific case of monitoring the blending of High Level Waste from multiple tanks prior to introduction into a vitrification facility. Campaigns for blending this waste to minimize the risk of exceeding WACs and PCLs have been proposed. However, the predicted compositions of the blended materials do not incorporate any uncertainties that may be associated with the representativeness of the waste layer samples or the laboratory analyses that generated the BBI data. Also, they do not include any uncertainty associated with the precision of collecting highly specific fractions of the layers during a blending campaign or any inhomogeneities that may exist in those layers. Monte Carlo methods are used to apply uncertainties to the compositions of the individual layers specified in the campaign recipes. The resulting variations in the compositions of the blended materials allow estimation of the risks of exceeding WACs and PCLs for each campaign. A critical subset of WACs/PCLs – NOx, NaK, AlFeZr, and S – are especially at risk of being exceeded in multiple campaigns. These analytes should be the focus of instrument development. We also have extracted the expected solid/supernate distribution for these analytes, which establishes important performance criteria for individual analytical methods. The BBI data also permits an understanding of the different chemical forms in which the analytes appear. Thus, the need to establish instrumental sensitivity to these forms can be gainfully addressed. Although concentrating on one specific application – the blending of tank waste - this approach should be generalizable for the analysis of other possible RTM applications for waste processing.

Lascola, Robert [Savannah River National Laborator↗

Controlling Structural, Electronic, and Energy Flow Dynamics of Catalytic Processes through Tailored Nanostructures

MoS 2 (molybdenum disulfide) is a highly-versatile catalyst material for support of numerous reactions from hydrodesulfurization and denitrogenation to the focus of this renewal proposal: hydrogenation of CO/CO 2 towards (higher) alcohols. At the same time, MoS 2 is a non-toxic, environmentally-benign and rather inert material- which under ambient conditions for some time has served as a lubricant and, more recently, as next-generation electronic material. The apparent contrast between inertness and stability in ambient, and catalytic activity under reactive conditions is puzzling and calls for a synergistic theoretical and experimental investigation with the long-term objective of enabling the rational design of MoS 2 -based catalysts for alcohol-formation reactions by providing a microscopic understanding of the environmental factors that determines site activity and selectivity. Our research project seeks answers to the questions (a) what conformation does MoS 2 adopt under reaction conditions (as opposed to that under ultrahigh vacuum and low temperatures)?; (b) what reaction pathways exist on such a material?; (c) how can the local environment of the active sites be manipulated so as to make MoS 2 an efficient catalyst for production of higher alcohol from syngas? In particular, research strategies will explore how the basal plane composed of sulfur atoms can be activated so as to exhibit a reactivity of its own, by addressing three research targets and building on extensive preliminary and enabling work: (1) vacancies and vacancy aggregates on the basal plane; (2) non-local catalyst transformation through alkali doping, hydrogenation and phase transition; (3) fabrication of a metal-nanoparticle-activated MoS 2 system, in which particle anchoring, reactive sites and pathways as well as selectivity are controlled by design, as an example of predictive development of a catalyst material,. All strategies are directed to improve the efficacy of the key reactive sites and selectivity of chemical pathways by design, to replace the inefficient methodology of trial and error in catalyst development. This research project represents a synergistic combination of computational guidance, foundational surface-science-based experiments and validation under reactive conditions that aims at transformative new insights into the working of MoS 2 -based hydrogenation-catalysts, a topic squarely at the center of the interest of DOE BES. Guided and led by Talat Rahman, a computational physicist, this project will apply density functional theory to understand structure, reaction pathways and chemical potentials associated with MoS 2 -based CO/CO 2 hydrogenation, augmented by kinetic Monte Carlo methods for reaction rates and prefactors as well as ab-initio molecular dynamics for evaluation of thermal stability. Complementary experimental input and validation will originate from co-PI Ludwig Bartels, a physical chemist and materials scientist, whose group focuses on local imaging and preparation of MoS 2 materials, from co-PI Peter Dowben, an experimental physicist, whose group is expert in the spectroscopy of occupied and unoccupied electronic states, and from senior collaborator Michael White of Brookhaven National Laboratory, whose group generates high-resolution electronic and activity information on size-selected well-defined metal chalcogenide clusters. This collaborative effort will enable a comprehensive understanding of the correlation of structural integrity and catalytic activity of MoS 2 in forms ranging from extended films to individual particles with known geometries and binding sites. Alcohol formation from syngas is a rapidly emerging application that has great potential through facile, economic, and decentralized biomass gasification. CO 2 activation is one of the most pressing concerns of our time: increasing CO 2 levels in the atmosphere change the climate and expose the globe to environmental transformations with the potential for enormous economic and societal impact. We will investigate CO/CO 2 hydrogenation via an interdisciplinary research collaboration with established synergy – one that, in accordance with the mission of the DOE, involves accredited Hispanic-Serving Institutions and that, through student exchange with international collaborators and National Labs, will directly benefit a broad spectrum of communities and generate human resources in sciences essential to their future.

2D materials↗

Predicting deviations in ex-core detector response in Krsko NPP during an earthquake with Monte-Carlo neutron transport methods

To study the response of ex-core detectors during an earthquake, an approach using Monte Carlo code MCNP was developed. As a first step to predict the response of the ex-core detectors a detailed MCNP model of the Krsko reactor core was created. A detailed MCNP core model was used to calculate 3D power distributions inside the core. Calculated power distributions were verified by comparison to the CORD-2 calculations, which is currently used for core design calculation verification of the Krsko nuclear power plant. Due to the outside the core position of the neutron detectors, the hybrid code ADVANTG is used to generate weight windows to speed up neutron transport outside the reactor core. To be able to use ADVANTG, fixed neutron source had to be reconstructed from the criticality core calculation. It was found that the power range ex-core detector movement has a negligible effect on the value of the thermal neutron flux at the active part of the detector. However, the 5 cm radial movement of the intermediate-range detector leads to 7 % - 8 % change in thermal neutron flux within the active intermediate-range detector region. The analysis continued with the evaluation of the effect of core barrel movement on the ex-core detector response. It was determined that the 2 mm core barrel radial oscillation can lead to ∼ 4% change in thermal neutron flux within the active detector region. It was evaluated that reactor core movement inside the baffle can lead up to ∼ 1.9% and ∼ 1.3% deviation in power-range and intermediate-range detector signal respectively. Analysis showed that the mechanical movement of ex-core neutron detectors could not explain the fluctuations in the ex-core detector signal. However, combined core barrel and reactor core inside baffle oscillations could be a probable reason for the observed fluctuations in the ex-core detector signal during an earthquake. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Explicitly antisymmetrized neural network layers for variational Monte Carlo simulation

The combination of neural networks and quantum Monte Carlo methods has arisen as a promising path forward for highly accurate electronic structure calculations. Previous proposals have combined equivariant neural network layers with a final antisymmetric layer in order to satisfy the antisymmetry requirements of the electronic wavefunction. However, to date it is unclear if one can represent antisymmetric functions of physical interest, and it is difficult to precisely measure the expressiveness of the antisymmetric layer. Here, this work attempts to address this problem by introducing explicitly antisymmetrized universal neural network layers. This approach has a computational cost which increases factorially with respect to the system size, but we are nonetheless able to apply it to small systems to better understand how the structure of the antisymmetric layer affects its performance. We first introduce a generic antisymmetric (GA) neural network layer, which we use to replace the entire antisymmetric layer of the highly accurate ansatz known as the FermiNet. We demonstrate that the resulting FermiNet-GA architecture can yield effectively the exact ground state energy for small atoms and molecules. We then consider a factorized antisymmetric (FA) layer which more directly generalizes the FermiNet by replacing the products of determinants with products of antisymmetrized neural networks. We find, interestingly, that the resulting FermiNet-FA architecture does not significantly outperform the FermiNet. This strongly suggests that the sum of products of antisymmetries is a key limiting aspect of the FermiNet architecture. To explore this further, we investigate a slight modification of the FermiNet, called the full determinant mode, which replaces each product of determinants with a single combined determinant. We find that the full single-determinant FermiNet closes a large part of the gap between the standard single-determinant FermiNet and FermiNet-GA on small atomic and molecular problems. Surprisingly, on the nitrogen molecule at a dissociating bond length of 4.0 Bohr, the full single-determinant FermiNet can outperform the largest standard FermiNet calculation with 64 determinants.

97 MATHEMATICS AND COMPUTING↗

Combining Experimental and Theoretical Techniques to Gain an Atomic Level Understanding of the Defect Binding Mechanism in Hard Carbon Anodes for Sodium Ion Batteries

Sodium ion batteries (NIBs) are an attractive alternative to lithium–ion batteries in applications that require large–scale energy storage due to sodium's high natural abundance and low cost. Hard carbon (HC) is the most promising anode material for NIBs; however, there is a knowledge gap in the understanding of the sodium binding mechanism that prevents a rational design of HC. This study tunes sucrose–derived HC via synthesis temperature then evaluates the structural, physical, and electrochemical properties. Neutron total scattering is used to generate structural models by fitting pair distribution functions (PDF) with a combination of molecular dynamics and reverse Monte Carlo methods. From this model, the number and type of structural features are identified, quantified, and correlated to the galvanostatic charge/discharge. A method of PDF “fingerprinting” binding sites using Na probe atoms is developed and analyzing these PDFs reveals an atomistic view of ion binding sites responsible for “defect” storage mechanisms. Combining these techniques results in an atomic–level study that provides a big picture of the Na–binding mechanism in NIBs, which allows for more precise tuning of the structure–property relationships in the future. Finally, the methodologies developed will also enable new strategies for the analysis of amorphous functional materials.

25 ENERGY STORAGE↗

Structural Properties of Nanometer‐Sized Gold Nanoparticles on a Silicon Substrate

Properties of gold nanoparticles (Au NPs) supported on Si substrates, which have been the focus of recent attention, depend sensitively on the size, structure, and coordination of Au NPs as well as on the processes of interdiffusion and chemical bonding at the Au–Si interface. There is a great need, therefore, for developing structural techniques capable of imaging such nanosystems in atomic detail. Here, it is discussed how the atomic structure of nanometer‐sized Au NPs deposited on a silicon surface can be determined by using high‐energy X‐ray diffraction measurements in combination with an analysis based on atomic‐pair‐distribution functions and advanced reverse Monte Carlo methods. This approach delivers simple and incisive pictures of the structure of Au clusters on Si and reveals direct structural evidence for the emergence of a Au‐silicide interface following room temperature deposition of Au NPs on Si. This study will contribute to designing Au/Si nanostructures with desirable properties.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Temperature-Dependent Configurational Entropy Calculations for Refractory High-Entropy Alloys

The cluster expansion formalism for alloys is used to construct surrogate models for three refractory high-entropy alloys (NbTiVZr, HfNbTaTiZr, and AlHfNbTaTiZr). These cluster expansion models are then used along with Monte Carlo methods and thermodynamic integration to calculate the configurational entropy of these refractory high-entropy alloys as a function of temperature. Many solid solution alloy design guidelines are based on the ideal entropy of mixing, which increases monotonically with N , the number of elements in the alloy. However, our results show that at low temperatures, the configurational entropy of these materials is largely independent of N , and the assumption described above only holds in the high-temperature limit. This suggests that alloy design guidelines based on the ideal entropy of mixing require further examination.

36 MATERIALS SCIENCE↗

High temperature stabilization of ultrafine grain tungsten alloys through synergistic compositional complexities

Thermally-stabilized fine-grained microstructures in tungsten provide a pathway to harnessing enhanced properties such as a reduced ductile-to-brittle transition temperature and improved strength while mitigating the adverse effects of grain growth and recrystallization. Here, in this study, we employ a material design strategy that relies on grain boundary segregation in the nanocrystalline state driven by alloy thermodynamics balanced with impurity scavenging through in situ formation of kinetically-stabilizing metal carbides. A W-Ti-Cr alloy is designed through lattice Monte Carlo methods and subsequently synthesized in a single-phase nanocrystalline state, which upon annealing, evolves into an ultrafine grained microstructure containing chromium grain boundary segregation collectively with a dispersed TiX (X=C,O) phase through the reaction of titanium with carbon and oxygen impurities. In situ synchrotron X-ray diffraction experiments demonstrate that increased Ti promotes stabilization across a larger temperature range but with diminishing returns above 10 at.% Ti. Long-term stability was confirmed through the retention of the ultrafine grained microstructure for a total of 8 days (192 h) at 1300 °C without grain/carbide growth and/or recrystallization. Our results demonstrate that, through strategic tailoring of composition and microstructure, one can harness the benefits of both thermodynamic and kinetic stabilization mechanisms, opening pathways for future alloy formulations that expand the window of stability.

36 MATERIALS SCIENCE↗

SCALE capabilities for high temperature gas-cooled reactor analysis

The SCALE code system’s ability to address stochastic distributions of fuel particles within a graphite matrix has been revisited in both multigroup (MG) features and continuous-energy (CE) Monte Carlo methods. Furthermore, this paper describes current and emergent SCALE capabilities within the CSAS sequence to address double-heterogeneous systems and presents verification and validation studies of these methods and data. Good agreement was obtained for a high temperature gas-cooled reactor (HTGR) fuel pebble model between CSAS MG eigenvalue calculations and corresponding CE reference solutions. Code-to-code comparisons for this HTGR pebble model showed good agreement of CSAS-KENO and CSAS-Shift CE calculations and the Serpent and MCNP codes in terms of eigenvalues and reaction rate ratios. Validation studies based on two HTGR experiments resulted in good agreement between MG and CE results, as well as between experiment and calculation, although the level of agreement was significantly influenced by the applied ENDF/B nuclear data library.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Catalytic resonance theory for the kinetics of photon-promoted catalysis

The illumination of catalytic surfaces with a continuous or pulsed stream of photons dynamically modulates surface chemistry for faster rates, non-equilibrium conversion, or product selectivity control. To establish fundamental principles of dynamic photon-modulated catalysis, the photocatalytic conversion of a generic surface reaction was simulated using the kinetic Monte Carlo method to understand the kinetic implications of an independent stream of photons that promotes surface product desorption. The time-averaged photocatalytic rate at differential conditions for varying photon flux and temperatures indicated three kinetic regimes described by product thermal desorption control, surface reaction control, and an intermediate kinetic regime with a zero slope Arrhenius plot, consistent with a degree of rate control dominated by the photon arrival frequency (i.e., per-site photon flux). Here, the maximum photocatalytic rate occurred orders of magnitude above the Sabatier limit at the resonance frequency, identified as the photon arrival frequency matching the surface reaction rate constant.

Canavan, Jesse R. [University of Minnesota, Minnea↗

Neural network model of neutral beam injection in the EAST tokamak to enable fast transport simulations

The neutral beam injection (NBI) system in EAST produces energetic neutral particles, which collide with electrons and ions in tokamak plasmas and heat the plasmas through Coulomb collisions. Moreover, it drives a non-inductive source of current, due to the charge-exchange collision between neutral particles and ions, and injects toroidal torque, which generates a toroidal rotation of the plasma. The effect caused by the NBI system, such as plasma heating, current drive, total neutron rate, momentum transfer, and shine-through, are modeled by a comprehensive module called NUBEAM. However, NUBEAM is computationally intensive since it relies on Monte Carlo methods. In this work, a neural network model has been developed as a surrogate model for NUBEAM in EAST. The database for neural-network model training, validation and testing is generated by running TRANSP for experimental discharges from recent EAST campaigns (after the latest NBI upgrade) while using the NUBEAM module. Simulation results illustrate that the trained neural network has the capability of replicating the predictions made by NUBEAM while demanding a significantly shorter execution time. Finally, these results indicate that surrogate models like the one proposed in this work could enable fast transport simulations for EAST after integrating them into a control-oriented predictive code such as COTSIM.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First operation and validation of simulations for the divertor cryo-vacuum pump in Wendelstein 7-X

Ten cryo-vacuum pumps (CVPs) were installed in the subdivertor region of each island divertor in the stellarator Wendelstein 7-X (W7-X) and operated for the first time during the recently completed plasma campaign OP2.1. A pumping speed of 70 ± 1 $\frac{m^3}{s}$ was measured during dedicated tests with known hydrogen gas injection. Based on a conductance model, the estimated pumping speed ranges from 86-93 $\frac{m^3}{s}$ for different sticking coefficients between 0.6 and 0.8. After completion of the initial tests the CVPs were operated successfully throughout the campaign, with regeneration performed once a week. Neutral gas pressures in the subdivertor in the range of 10 −4 mbar are well within the molecular flow regime and limit the particle exhaust capabilities of the CVPs. Simulations of the neutral gas pressure in the three-dimensional complex geometry of the subdivertor were performed using the DIVGAS code based on the direct simulation Monte Carlo method and a model implemented in the steady-state thermal package in ANSYS, which are in agreement with the measured values during plasma operation.

Cryo-vacuum pumping↗

Efficient particle control in systems with large density gradients

Simulations of large density gradients present a number of challenges for direct Monte Carlo methods, since they lead to too few particles in dilute regions and prohibitively many in the dense regions. Here, we propose a particle control methodology that gives the user more control of the number of particles per cell by introducing a variable weight for each particle. The proposed scheme is based on the stochastic weighted particle method, requires only small modifications to DSMC, and exactly conserves mass, momentum, and energy. In validation tests of systems with density ratios of order 102-104, we observe 101-102 times less variance in the dilute region compared to a DSMC solution with the same number of system particles, while introducing a moderate additional computational cost.

DSMC↗

A Kaczmarz-inspired approach to accelerate the optimization of neural network wavefunctions

Neural network wavefunctions optimized using the variational Monte Carlo method have been shown to produce highly accurate results for the electronic structure of atoms and small molecules, but the high cost of optimizing such wavefunctions prevents their application to larger systems. We propose the Subsampled Projected-Increment Natural Gradient Descent (SPRING) optimizer to reduce this bottleneck. SPRING combines ideas from the recently introduced minimum-step stochastic reconfiguration optimizer (MinSR) and the classical randomized Kaczmarz method for solving linear least-squares problems. We demonstrate that SPRING outperforms both MinSR and the popular Kronecker-Factored Approximate Curvature method (KFAC) across a number of small atoms and molecules, given that the learning rates of all methods are optimally tuned. For example, on the oxygen atom, SPRING attains chemical accuracy after forty thousand training iterations, whereas both MinSR and KFAC fail to do so even after one hundred thousand iterations.

97 MATHEMATICS AND COMPUTING↗