Constructing semilocal approximations for noncollinear spin density functional theory featuring exchange-correlation torques
Not Available
74 ATOMIC AND MOLECULAR PHYSICS↗
SEARCH · Engineering Papers
Results for “function approximation”
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.
High-throughput determination of Hubbard $U$ and Hund $J$ values for transition metal oxides via the linear response formalism
DFT+U provides a convenient, cost-effective correction for the self-interaction error (SIE) that arises when describing correlated electronic states using conventional approximate density functional theory (DFT). The success of a DFT+U(+J) calculation hinges on the accurate determination of its Hubbard U and Hund J parameters, and the linear response (LR) methodology has proven to be computationally effective and accurate for calculating these parameters. This study provides a high-throughput computational analysis of the U and J values for transition metal d-electron states in a representative set of over 1000 magnetic transition metal oxides (TMOs), providing a frame of reference for researchers who use DFT+U to study transition metal oxides. In order to perform this high-throughput study, an ATOMATE workflow is developed for calculating U and J values automatically on massively parallel supercomputing architectures. Here, to demonstrate an application of this workflow, the spin-canting magnetic structure and unit cell parameters of the multiferroic olivine LiNiPO4 are calculated using the computed Hubbard U and Hund J values for Ni-d and O-p states, and are compared with experiment. Both the Ni-d U and J corrections have a strong effect on the Ni-moment canting angle. Additionally, including a O-pU value results in a significantly improved agreement between the computed lattice parameters and experiment
36 MATERIALS SCIENCE↗
Control of Three Degrees-of-Freedom Wave Energy Converters Using Pseudo-Spectral Methods
This study presents a solution to the optimal control problem of a three degrees-of-freedom (3DOF) wave energy converter (WEC). The three modes are the heave, pitch, and surge. The dynamic model is characterized by a coupling between the pitch and surge modes, while the heave is decoupled. The heave, however, excites the pitch motion through nonlinear parametric excitation in the pitch mode. This paper uses Fourier series (FS) as basis functions to approximate the states and the control. A simplified model is first used where the parametric excitation term is neglected and a closed-form solution for the optimal control is developed. For the parametrically excited case, a sequential quadratic programming approach is implemented to solve for the optimal control numerically. Numerical results show that the harvested energy from three modes is greater than three times the harvested energy from the heave mode alone. Moreover, the harvested energy using a control that accounts for the parametric excitation is significantly higher than the energy harvested when neglecting this nonlinear parametric excitation term.
16 TIDAL AND WAVE POWER↗
Field Characteristics of Spin-Torque Diode Sensitivity in the Presence of a Bias Current
The use of the spin-torque diode effect, which is generated upon the current-induced transfer of spin angular momentum in magnetic tunnel junctions, opens the prospect of considerably improving microwave sensitivity in the gigahertz frequency band over that of semiconductor Schottky diodes. In this work, the spin-diode effect of microwave signal rectification in a magnetic tunnel junction upon the resonant excitation of spin waves in the free magnetic layer as a result of the current-induced spin-transfer torque effect is analyzed theoretically. Frequency characteristics of the resonant response of a spin-torque diode to a microwave signal are calculated in the linear macrospin approximation as functions of the direction and value of the applied magnetic field and bias current. It is shown that with zero bias current, the maximum rectified voltage across a junction is obtained in the geometry of mutually perpendicular magnetizations of magnetic layers; as the resonant frequency of oscillations in a magnetic field increases, this voltage drops at the retained equilibrium orientation of spins in the layers. After the bias current is switched on, the resonance amplitude of forced oscillations of the free layer spins upon microwave excitation grows sharply at the critical point of the spin-torque diode’s loss of equilibrium state stability. The linewidth at this point is limited only by nonlinear effects. Improving spin-torque diode sensitivity is important for its use in microwave holographic vision systems.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Improved unconstrained energy functional method for eigensolvers in electronic structure calculations
This paper reports on the performance of a preconditioned conjugate gradient based iterative eigensolver using an unconstrained energy functional minimization scheme. In contrast to standard implementations, this scheme avoids an explicit reorthogonalization of the trial eigenvectors and becomes an attractive alternative for the solution of very large problems. The unconstrained formulation is implemented in the first-principles materials and chemistry CP2K code, which performs electronic structure calculations based on a density functional theory approximation to the solution of the many-body Schrödinger equation. We study the convergence of the unconstrained formulation, as well as its parallel scaling, on a Cray XC40 at the National Energy Research Scientific Computing Center (NERSC). The systems we use in our studies are bulk liquid water, a supramolecular catalyst gold(III)-complex, a bilayer of MoS2-WSe2 and a divacancy point defect in silicon, with the number of atoms ranging from 2,247 to 12,288. We show that the unconstrained formulation with an appropriate preconditioner has good convergence properties and scales well to 230k cores, roughly 38% of the full machine.
Del Ben, M↗
Multi-fidelity Gaussian Process Bandit Optimisation
In many scientific and engineering applications, we are tasked with the maximisation of an expensive to evaluate black box function f. Traditional settings for this problem assume just the availability of this single function. However, in many cases, cheap approximations to f may be obtainable. For example, the expensive real world behaviour of a robot can be approximated by a cheap computer simulation. We can use these approximations to eliminate low function value regions cheaply and use the expensive evaluations of f in a small but promising region and speedily identify the optimum. We formalise this task as a multi-fidelity bandit problem where the target function and its approximations are sampled from a Gaussian process. We develop MF-GP-UCB, a novel method based on upper confidence bound techniques. In our theoretical analysis we demonstrate that it exhibits precisely the above behaviour and achieves better bounds on the regret than strategies which ignore multi-fidelity information. Empirically, MF-GP-UCB outperforms such naive strategies and other multi-fidelity methods on several synthetic and real experiments.
97 MATHEMATICS AND COMPUTING↗
WSBP/NWBP small-world neural network and its performance for wind power forecasting
A small-world neural network has stronger generalization ability with high transfer efficiency than that of the regular neural networks. This paper presents two novel small-world neural networks, the Watts-Strogatz small-world based on a BP neural network (WSBP) and a Newman-Watts small-world neural network based on a BP neural network (NWBP), related to previous research of complex networks. The algorithms are developed separately by adopting WS and NW small-world networks as their topological structures, and their derivation and convergence criterion are progressively discussed. After that, the proposed models are subsequently tested by two typical nonlinear functions which confirm their significant improvement over the regular BP networks and other algorithms. Finally, a wind power prediction system is advanced to verify their generalization abilities, and show that the models are practically feasible and effective with improved accuracy and acceptable forecasting errors caused by wind fluctuation and randomness with a time scale up to 24 h.
17 WIND ENERGY↗
SunSpec Open Solar Data Exchange (oSDX) (Final Report)
Orange Button is a data interchange standard for solar finance and operations practitioners designed to drive down soft costs. To achieve its goals, SunSpec convened 330+ companies to validate Orange Button functional requirements. Approximately 4,000 industry-expert hours were invested to dissect 25 priority use cases and develop 4,500 Orange Button taxonomy entries and accompanying APIs. Orange Button 1.0 is now available on GitHub royalty-free to the general public. Early adopters include highly influential companies such as Aurora Solar, BayWa r.e., Clean Power Research, DNV-GL, LG, Mercatus, Oracle, Salesforce, sPower, SunPower, and Wells Fargo. Finally, the project was completed on time and on budget, including $2M match funding provided by industry.
14 SOLAR ENERGY↗
Checkpoint-Restart and Platform Design Criteria
Long computations often use checkpoint/restart to survive the inevitable failures. Optimal checkpoint/restart strategy can be derived in terms of special functions, and approximated within the regime of interest. As a further step under nominal assumptions, platforms can be optimized for maximum progress rate, computed in terms of non-dimensional parameter JMTTI/Delta, where JMTTI is the Job Mean Time To Interrupt and Delta is the time it takes to write (or read) a checkpoint. Scaling to a larger platform of the same technology requires scaling checkpoint bandwidth quadratically to maintain the same progress rate. This was observed during the design of LANL’s Roadrunner in 2007, led to burst buffers proposed in 2009 and introduced on Trinity in 2015. The need to compute reliably continues to impact the design of future exascale platforms.
97 MATHEMATICS AND COMPUTING↗
Gradient-based Optimization for Regression in the Functional Tensor-Train Format
We consider the task of low-multilinear-rank functional regression, i.e., learning a low-rank parametric representation of functions from scattered real-valued data. Our first contribution is the development and analysis of gradient-based optimization procedures, including stochastic gradient descent and quasi-Newton methods, for learning the parameters of a functional tensor-train (FT) that yields improved accuracy over standard alternating least squares methods. The functional tensor-train uses the tensor-train (TT) representation of low-rank arrays as an ansatz for a class of low-multilinear-rank functions. The FT is represented by a set of matrix-valued functions that contain a set of univariate functions, and the regression task is to learn the parameters of these univariate functions. Our second contribution demonstrates that using nonlinearly parameterized univariate functions, e.g., symmetric kernels with moving centers, within each core can outperform the standard approach of using a linear expansion of basis functions. Our final contributions are new rank adaptation and group-sparsity regularization procedures to minimize overfitting. We use several benchmark problems to demonstrate at least an order of magnitude lower accuracy with gradient-based optimization methods than ALS in the low-sample number regime. We also demonstrate an order of magnitude reduction in accuracy on a test problem resulting from nonlinear parameterizations over linear parameterizations. Finally we compare regression performance with 22 other nonparametric and parametric regression methods on 10 real-world data sets. We achieve top-five accuracy for seven of the data sets and best accuracy for two of the data sets. These rankings are the best amongst parametric models and competetive with the best non-parametric methods.
97 MATHEMATICS AND COMPUTING↗
AIACHNE's contribution for Nuclear Energy Agency Working Party on International Nuclear Data Evaluation Co-operation Subgroup 50
The AIACHNE (AI/ML Informed cAlifornium CHi Nuclear data Experiment) project aims at designing an experiment for the 252 Cf Prompt Fission Neutron Spectrum (PFNS) that explores systematic biases in an experimental database retrieved from the EXFOR databases. To that end, machine learning (ML) methods were applied to pint-point measurement features likely related to bia. From that information, we selected a feature that should be explored by the AIACHNE experiment. Measurement features are metadata encapsulating all pertinent information about the physical measurement and analysis techniques. Examples are, for instance, what neutron and fission detectors were used for the physical metadata, and what background reduction techniques were employed for analysis techniques. Such metadata were retrieved both from EXFOR entries as well as the literature of data sets described in detail in Ref. [2]. The prerequisite for applying machine learning techniques is casting the metadata into a format that can be parsed by the algorithm. This step might seem trivial but requires to find a unique language where metadata that carry the same physics meaning across several experiments must have the same identifier. One example is, for instance, the neutron detector. As seen in Figure 1, the machine learning code identified the use of 6 Li detectors as being related to bias in some datasets of the AIACHNE 252 Cf PFNS experimental database. In fact, here are several experiments that used neutron detectors containing 6Li in the database, for instance for the example below. EXFOR format has a unique keywords describing detectors such as “SCIN” or “GLASD”. One may think that these keywords are already sufficient descriptors for ML to uniquely find an issue. However, “SCIN” (used for [3, 4]) and “GLASD” (used for [5]) fail to inform the algorithm what is the active material in the detector. And, the key common issue leading to bias in 252 Cf related to neutron detectors is not whether it is a glass detector or a scintillator. No, the issue is that 6 Li was within both detector types and that even small mistakes in the detector response functions around approximately 200 keV are amplified by the 6 Li(n,α) resonance there leading to bias in data as highlighted in Fig. 1 and Ref. [1]. Hence, the features describing the neutron detector must call out the active material in the detector, rather than the existing EXFOR detector keyword, that the ML algorithm can find physically meaningful features related to bias. The AIACHNE team used a precursor of the WPEC (Working Party on International Nuclear Data Evaluation Co-operation) SG(Subgroup)-50 format to store the metadata for the ML analysis.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
The Evolution and Mass Dependence of Galaxy Cluster Pressure Profiles at 0.05 ≤ z ≤ 0.60 and 4 × 10 14 M ⊙ ≤ M 500 ≤ 30 × 10 14 M ⊙
Abstract We have combined X-ray observations from Chandra with Sunyaev–Zel’dovich effect data from Planck and Bolocam to measure intracluster medium pressure profiles from 0.03 R 500 ≤ R ≤ 5 R 500 for a sample of 21 low- z galaxy clusters with a median redshift of 〈 z 〉 = 0.08 and a median mass of 〈 M 500 〉 = 6.1 × 10 14 M ⊙ and a sample of 19 mid- z galaxy clusters with 〈 z 〉 = 0.50 and 〈 M 500 〉 = 10.6 × 10 14 M ⊙ . The mean scaled pressure in the low- z sample is lower at small radii and higher at large radii, a trend that is accurately reproduced in similarly selected samples from The Three Hundred simulations. This difference appears to be primarily due to dynamical state at small radii, evolution at intermediate radii, and a combination of evolution and mass dependence at large radii. Furthermore, the overall flattening of the mean scaled pressure profile in the low- z sample compared to the mid- z sample is consistent with expectations due to differences in the mass accretion rate and the fractional impact of feedback mechanisms. In agreement with previous studies, the fractional scatter about the mean scaled pressure profile reaches a minimum of ≃20% near 0.5 R 500 . This scatter is consistent between the low- z and mid- z samples at all radii, suggesting it is not strongly impacted by sample selection, and this general behavior is reproduced in The Three Hundred simulations. Finally, analytic functions that approximately describe the mass and redshift trends in mean pressure profile shape are provided.
79 ASTRONOMY AND ASTROPHYSICS↗
Baryon fraction from the BAO amplitude: a consistent approach to parameterizing perturbation growth
Galaxy clustering constrains the baryon fraction Omega_b/Omega_m through the amplitude of baryon acoustic oscillations and the suppression of perturbations entering the horizon before recombination. This produces a different pre-recombination distribution of baryons and dark matter. After recombination, the gravitational potential responds to both components in proportion to their mass, allowing robust measurement of the baryon fraction. This is independent of new-physics scenarios altering the recombination background (e.g. Early Dark Energy). The accuracy of such measurements does, however, depend on how baryons and CDM are modeled in the power spectrum. Previous template-based splitting relied on approximate transfer functions that neglected part of information. We present a new method that embeds an extra parameter controlling the balance between baryons and dark matter in the growth terms of the perturbation equations in the CAMB Boltzmann solver. This approach captures the baryonic suppression of CDM prior to recombination, avoids inconsistencies, and yields a clean parametrization of the baryon fraction in the linear power spectrum, separating out the simple physics of growth due to the combined matter potential. We implement this framework in an analysis pipeline using Effective Field Theory of Large-Scale Structure with HOD-informed priors and validate it against noiseless LCDM and EDE cosmologies with DESI-like errors. The new scheme achieves comparable precision to previous splitting while reducing systematic biases, providing a more robust way to baryon-fraction measurements. In combination with BBN constraints on the baryon density and Alcock-Paczynski estimates of the matter density, these results strengthen the use of baryon fraction measurements to derive a Hubble constant from energy densities, with future DESI and Euclid data expected to deliver competitive constraints.
Crespi, Andrea [U. Waterloo (main); Waterloo U., I↗
RTN-095: The Vera C. Rubin Observatory Data Preview 1
We present Rubin Data Preview 1 (DP1), the first release of data from the NSF-DOE Vera C. Rubin Observatory, consisting of raw and calibrated single-epoch images, coadds, difference images, detection catalogs, and other derived data products. DP1 is based on 1792 science-grade optical/near-infrared exposures acquired over 48 distinct nights by the Rubin Commissioning Camera, LSSTComCam, on the Simonyi Survey Telescope at the Summit Facility on Cerro Pachón, Chile during the first on-sky commissioning campaign in late 2024. DP1 covers a total of ~15 sq. deg. over seven roughly equally-sized non-contiguous fields, each independently observed in six broad photometric bands, ugrizy, spanning a range of stellar densities and latitudes and overlapping with external reference datasets. The median image quality across all bands, measured by the FWHM of the point-spread function, is approximately 1.13 arcseconds, with the sharpest images reaching about 0.65 arcseconds. DP1 contains approximately 2.3 million distinct astrophysical objects, of which 1.6 million are extended in at least one band, and 431 solar system objects, of which 93 are new discoveries. DP1 is approximately 3.5 TB in size and available to Rubin data rights holders via the Rubin Science Platform, a cloud-based environment for the analysis of petascale astronomical data. While small compared to future LSST releases, its high quality and diversity of data support a broad range of early science investigations across all four LSST themes, providing a valuable opportunity to engage with Rubin data ahead of the start of full operations in late 2025.
79 ASTRONOMY AND ASTROPHYSICS↗
The Eikonal Approximation for the Bilocal Vertex Function and NuW(2) in a Fermion-Neutral Vector Cluon Model
Explore the source record for details and available documents.
Eichten, E.↗
Pseudo-spectral method to control three-degree-of-freedom wave energy converters
The invention provides optimal control of a three-degree-of-freedom wave energy converter using a pseudo-spectral control method. The three modes are the heave, pitch and surge. A dynamic model is characterized by a coupling between the pitch and surge modes, while the heave is decoupled. The heave, however, excites the pitch motion through nonlinear parametric excitation in the pitch mode. The invention can use a Fourier series as basis functions to approximate the states and the control. For the parametric excited case, a sequential quadratic programming approach can be implemented to numerically solve for the optimal control. The numerical results show that the harvested energy from three modes is greater than three times the harvested energy from the heave mode alone. Moreover, the harvested energy using a control that accounts for the parametric excitation is significantly higher than the energy harvested when neglecting this nonlinear parametric excitation term.
Abdelkhalik, Ossama↗
Artificial intelligence “sees” split electrons
Chemical bonds between atoms are stabilized by the exchange-correlation (xc) energy, a quantum-mechanical effect in which “social distancing” by electrons lowers their electrostatic repulsion energy. Kohn-Sham density functional theory (DFT) states that the electron density determines this xc energy, but the density functional must be approximated. Furthermore, this is usually done by satisfying exact constraints of the exact functional (making the approximation predictive), by fitting to data (making it interpolative), or both. Two exact constraints—the ensemble-based piecewise linear variation of the total energy with respect to fractional electron number and fractional electron z-component of spin —require hard-to-control nonlocality. On page 1385 of this issue, Kirkpatrick et al. have taken a big step toward more accurate predictions for chemistry through the machine learning of molecular data plus the fractional charge and spin constraints, expressed as data that a machine can learn.
74 ATOMIC AND MOLECULAR PHYSICS↗
Specific Reaction Parameter Density Functional Based on the Meta-Generalized Gradient Approximation: Application to H 2 + Cu(111) and H 2 + Ag(111)
Specific reaction parameter density functionals (SRP-DFs), which can describe dissociative chemisorption reactions on metals to within chemical accuracy, have so far been based on exchange functionals within the generalized gradient approximation (GGA) and on GGA correlation functionals or van der Waals correlation functionals. These functionals are capable of describing the molecule–metal surface interaction accurately, but they suffer from the general GGA problem that this can be done only at the cost of a rather poor description of the metal. Here, we show that it is possible also to construct SRP-DFs for H 2 dissociation on Cu(111) based on meta-GGA functionals, introducing three new functionals based on the “made-simple” (MS) concept. The exchange parts of the three functionals (MS-PBEl, MS-B86bl, and MS-RPBEl) are based on the expressions for the PBE, B86b, and RPBE exchange functionals. Quasi-classical trajectory (QCT) calculations performed with potential energy surfaces (PESs) obtained with the three MS functionals reproduce molecular beam experiments on H 2 , D 2 + Cu(111) with chemical accuracy. Therefore, these three non-empirical functionals themselves are also capable of describing H2 dissociation on Cu(111) with chemical accuracy. Similarly, QCT calculations performed on the MS-PBEl and MS-B86bl PESs reproduced molecular beam and associative desorption experiments on D 2 , H 2 + Ag(111) more accurately than was possible with the SRP48 density functional for H 2 + Cu(111). Also, the three new MS functionals describe the Cu, Ag, Au, and Pt metals more accurately than the all-purpose Perdew–Burke–Ernzerhof (PBE) functional. The only disadvantage we noted of the new MS functionals is that, as found for the example of H 2 + Cu(111), the reaction barrier height obtained by taking weighted averages of the MS-PBEl and MS-RPBEl functionals is tunable over a smaller range (9 kJ/mol) than possible with the standard GGA PBE and RPBE functionals (33 kJ/mol). As a result of this restricted tunability, it is not possible to construct an SRP-DF for H 2 + Ag(111) on the basis of the three examined MS meta-GGA functionals.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗