Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solvers”

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 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Numerical simulation of RHIC polarized proton run 17 spin flipper experiments

RHIC nine-magnet spin flipper has been operated successfully during RHIC polarized proton Run 17, with 97% spin flip efficiency achieved. The results show the importance of mirror resonance removal, small spin tune spread, and proper spin flipper driving tune sweep speed. Detailed spin tracking simulations, based on a Lorentz force and Thomas-BMT differential equation numerical solver code for accuracy, have been carried out to understand the experimental results. Agreement within measurement accuracy is obtained at injection energy, 23.8 GeV. It is not as tight at 255 GeV, reasons for that are exposed. These measurements and numerical studies allow to determine the sensitivity of spin-flip efficiency to the dispersion slopes at the two Siberian snakes and to the ac dipole frequency sweep speed. They also provide guidance for future developments at BNL’s electron-ion collider. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Efficiency and Forward Voltage of Blue and Green Lateral LEDs with V-shaped Defects and Random Alloy Fluctuation in Quantum Wells

For nitride-based blue and green light-emitting diodes (LEDs), the forward voltage V for is larger than expected, especially for green LEDs. This is mainly due to the large barriers to vertical carrier transport caused by the total polarization discontinuity at multiple quantum well and quantum barrier interfaces. The natural random alloy fluctuation in quantum wells has proven to be an important factor reducing V for . However, this does not suffice in the case of green LEDs because of their larger polarization-induced barrier. V-shaped defects (V-defects) have been proposed as another key factor in reducing V for to allow lateral injection into multiple quantum wells, thus bypassing the multiple energy barriers incurred by vertical transport. In this paper, to model carrier transport in the whole LED, we consider both random-alloy and V-defect effects. A fully two-dimensional drift-diffusion charge-control solver is used to model both effects. The results indicate that the turn-on voltages for blue and green LEDs are both affected by random alloy fluctuations and the V-defect density. For green LEDs, V for decreases more due to V-defects, where the smaller polarization barrier at the V-defect sidewall is the major path for lateral carrier injection. Then, we discuss how the V-defect density and size affects the results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Emergent moments in a Hund's impurity

Here, motivated by the relevance of Hund's coupling in the context of multiorbital superconductors, we revisit the problem of a multiorbital Kondo impurity with Hund's interaction. Using dynamical large-N techniques, we propose an efficient approach that retains the essential physics at play, while providing a pathway to scalable quantum impurity solvers. We are able to follow the ground state, dynamic, and thermodynamic properties of this system over many decades of temperature. Our approach captures the emergence of large moments and follows the stretched evolution of the physics down to their exponentially suppressed Kondo temperature. We focus our analysis on the intermediate finite temperature phase which presents an alternate paramagnetic state due to the emergent moment and discuss the relevance of this regime to Hund's metals.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Predicting impurity spectral functions using machine learning

The Anderson Impurity Model (AIM) is a canonical model of quantum many-body physics. Here we investigate whether machine learning models, both neural networks (NN) and kernel ridge regression (KRR), can accurately predict the AIM spectral function in all of its regimes, from empty orbital, to mixed valence, to Kondo. To tackle this question, we construct two large spectral databases containing approximately 410 000 and 600 000 spectral functions of the single-channel impurity problem. We show that the NN models can accurately predict the AIM spectral function in all of its regimes, with pointwise mean absolute errors down to 0.003 in normalized units. We find that the trained NN models outperform models based on KRR and enjoy a speedup on the order of 10 5 over traditional AIM solvers. Finally, the required size of the training set of our model can be significantly reduced using farthest point sampling in the AIM parameter space, which is important for generalizing our method to more complicated multichannel impurity problems of relevance to predicting the properties of real materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Differentiating Hund from Mott physics in a three-band Hubbard-Hund model: Temperature dependence of spectral, transport, and thermodynamic properties

We study the interplay between Mott physics, driven by Coulomb repulsion U, and Hund physics, driven by Hund's coupling J, for a minimal model for Hund metals, the orbital-symmetric three-band Hubbard-Hund model (3HHM) for a lattice filling of 1/3. Hund-correlated metals are characterized by spin-orbital separation (SOS), a Hund's-rule-induced two-stage Kondo-type screening process, in which spin screening occurs at much lower energy scales than orbital screening. By contrast, in Mott-correlated metals, lying close to the phase boundary of a metal-insulator transition, the SOS window becomes negligibly small and the Hubbard bands are well separated. Using dynamical mean-field theory and the numerical renormalization group as real-frequency impurity solver, we identify numerous fingerprints distinguishing Hundness from Mottness in the temperature dependence of various physical quantities. Furthermore, These include ARPES-type spectra, the local self-energy, static local orbital and spin susceptibilities, resistivity, thermopower, and lattice and impurity entropies. Our detailed description of the behavior of these quantities within the context of a simple model Hamiltonian will be helpful for distinguishing Hundness from Mottness in experimental and theoretical studies of real materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Real-time evolution of Anderson impurity models via tensor network influence functionals

In this work, we present and analyze two tensor network-based influence functional approaches for simulating the real-time dynamics of quantum impurity models such as the Anderson model. Via comparison with recent numerically exact simulations, we show that such methods accurately capture the long-time nonequilibrium quench dynamics. The two parameters that must be controlled in these tensor network influence functional approaches are a time discretization (Trotter) error and a bond dimension (tensor network truncation) error. We show that the actual numerical uncertainties are controlled by an intricate interplay of these two approximations, which we demonstrate in different regimes. Our work opens the door to using these tensor network influence functional methods as general impurity solvers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Delocalized polaron and Burstein-Moss shift induced by Li in α–V 2 O 5 : A DFT + DMFT study

For this work, we performed density functional theory (DFT)+U and dynamical mean field theory (DMFT) calculations with a continuous-time quantum Monte Carlo impurity solver to investigate the electronic properties of V 2 O 5 and Li x V 2 O 5 (x = 0.125 and 0.25). Pristine V 2 O 5 is a charge-transfer insulator with strong O p–V d hybridization, and it exhibits a large band gap (E gap ) as well as a nonzero conduction-band (CB) gap.We show that the band gap, the number of d electrons of vanadium, N d , and the CB gap for V 2 O 5 obtained from our DMFT calculations are in excellent agreement with the experimental values. While the DFT +U approach replicates the experimental band gap, it overestimates the value of N d and underestimates the CB gap. In the presence of low Li doping, the electronic properties of V 2 O 5 are mainly driven by a polaronic mechanism, and electron spin resonance and electron nuclear double resonance spectroscopies observed the coexistence of free and bound polarons. Notably, our DMFT results identify both polaron types, with the bound polaron being energetically preferred, while the DFT +U method only predicts the free polaron. Our DMFT analysis also reveals that increased Li doping leads to electron filling in the conduction band, shifting the Fermi level. This result is consistent with the observed Burstein-Moss shift upon enhanced Li doping, and we thus demonstrate that the DFT + DMFT approach can be used for an accurate and realistic description of strongly correlated materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Effect of off-diagonal elements in the Wannier Hamiltonian on DFT + DMFT for low-symmetry materials: Study of Li 2 MnO 3

Here, we study the effect of the off-diagonal elements of the Wannier Hamiltonian on the electronic structure of the low-symmetry material Li 2 MnO 3 ( C2/m ), using dynamical mean field theory calculations with a continuous-time quantum Monte Carlo impurity solver. The presence of significant off-diagonal elements leads to a pronounced suppression of the energy gap. The off-diagonal elements are largest when the Wannier projection is used based on the global coordinate, and they remain substantial even with the projection using the local coordinate close to the direction of Mn-O bonds. We show that the energy gap is enhanced by the diagonalization of the Mn d block in the full p-d Hamiltonian with the application of a unitary rotation matrix. Additionally, the inclusion of small double counting energy is crucial for achieving the experimental gap by reducing p-d hybridization. Furthermore, we establish the efficiency of a low-energy (d-only basis) model for studying the electronic structure of Li 2 MnO 3 , as the Wannier basis represents a hybridized state of Mn d and O p orbitals. These findings suggest an appropriate approach for investigating low-symmetry materials using the density functional theory plus dynamical mean field theory (DFT + DMFT) method. We also find that the antiferromagnetic ground state $\Gamma$ 2u is stable with U ≤ 2 eV within density functional theory+$U$ calculations, which is much smaller than the widely used U = 5 eV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Noncollinear ground states of solids with a source-free exchange correlation functional

In this paper, we expand upon the source-free (SF) exchange correlation (XC) functional developed by Sangeeta Sharma and coworkers to plane-wave density functional theory (DFT) based on the projector augmented wave (PAW) method. This constraint is implemented by the current authors within the VASP source code, using a fast Poisson solver that capitalizes on the parallel three-dimensional fast Fourier transforms (FFTs) implemented in VASP. Using this modified XC functional, we explore the improved convergence behavior that results from applying this constraint to the GGA-PBE+U+J functional. In the process, we compare the noncollinear magnetic ground state computed by each functional and their SF counterpart for a select number of magnetic materials in order to provide a metric for comparing with experimentally determined magnetic orderings. We observe significantly improved agreement with experimentally measured magnetic ground-state structures after applying the source-free constraint. Furthermore, we explore the importance of considering probability current densities in spin-polarized systems, even under no applied field. We analyze the XC torque as well, in order to provide theoretical and computational analyses of the net XC magnetic torque induced by the source-free constraint. Along these lines, we highlight the importance of properly considering the real-space integral of the source-free local magnetic XC field. Our analyses on probability currents, net torque, and constant terms draw additional links to the rich body of previous research on spin-current density functional theory (SCDFT), and pave the way for future extensions and corrections to the SF corrected XC functional.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Landscape of pear-shaped even-even nuclei

The phenomenon of reflection-asymmetric nuclear shapes is relevant to nuclear stability, nuclear spectroscopy, nuclear decays and fission, and the search for new physics beyond the standard model. Global surveys of ground-state octupole deformation, performed with a limited number of models, suggest that the number of pear-shaped isotopes is fairly limited across the nuclear landscape. Here, we carry out a global analysis of ground-state octupole deformations for particle-bound even-even nuclei with Z ≤ 110 and N ≤ 210 using nuclear density functional theory (DFT) with several nonrelativistic and covariant energy density functionals. In this way, we can identify the best candidates for reflection-asymmetric shapes. The calculations are performed in the frameworks of axial reflection-asymmetric Hartree-Fock-Bogoliubov theory and relativistic Hartree-Bogoliubov theory using DFT solvers employing harmonic oscillator basis expansion. We consider five Skyrme and four covariant energy density functionals. We predict several regions of ground-state octupole deformation. In addition to the “traditional” regions of neutron-deficient actinide nuclei around 224 Ra and neutron-rich lanthanides around 146 Ba, we identified vast regions of reflection-asymmetric shapes in very neutron-rich nuclei around 200 Gd and 288 Pu, as well as in several nuclei around 112 Ba. Our analysis suggests several promising candidates with stable ground-state octupole deformation, primarily in the neutron-deficient actinide region, that can be reached experimentally. Detailed comparison between Skyrme and covariant models is performed. Octupole shapes predicted in this study are consistent with the current experimental information. This work can serve as the starting point of a systematic search for parity doublets in odd-mass and odd-odd nuclei, which will be of interest in the context of new physics searches.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Operator evolution from the similarity renormalization group and the Magnus expansion

The Magnus expansion is an efficient alternative to solving similarity renormalization group (SRG) flow equations with high-order, memory-intensive ordinary differential equation solvers. The numerical simplifications it offers for operator evolution are particularly valuable for in-medium SRG calculations, though challenges remain for difficult problems involving intruder states. Here we test the Magnus approach in an analogous but more accessible situation, which is the free-space SRG treatment of the spurious bound states arising from a leading-order chiral effective field theory (EFT) potential with very high cutoffs. In this work, we show that the Magnus expansion passes these tests and then use the investigations as a springboard to address various aspects of operator evolution that have renewed relevance in the context of the scale and scheme dependence of nuclear processes. These aspects include SRG operator flow with band- versus block-diagonal generators, universality for chiral EFT Hamiltonians and associated operators with different regularization schemes, and the impact of factorization arising from scale separation. Implications for short-range correlations physics and the possibilities for reconciling high- and low-resolution treatments of nuclear structure and reactions are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Training and projecting: A reduced basis method emulator for many-body physics

Here, we present the reduced basis method as a tool for developing emulators for equations with tun able parameters within the context of the nuclear many-body problem. The method uses a basis expansion informed by a set of solutions for a few values of the model parameters and then projects the equations over a well-chosen low-dimensional subspace. We connect some of the results in the eigenvector continuation literature to the formalism of reduced basis methods and show how these methods can be applied to a broad set of problems. As we illustrate, the possible success of the formalism on such problems can be diagnosed beforehand by a principal component analysis. We apply the reduced basis method to the one-dimensional Gross-Pitaevskii equation with a harmonic trap ping potential and to nuclear density functional theory for 48 Ca, achieving speed-ups of more than x150 and x250, respectively, when compared to traditional solvers. The outstanding performance of the approach, together with its straightforward implementation, show promise for its application to the emulation of computationally demanding calculations, including uncertainty quantification.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

ROSE: A reduced-order scattering emulator for optical models

A new generation of phenomenological optical potentials requires robust calibration and uncertainty quantification, motivating the use of Bayesian statistical methods. These Bayesian methods usually require calculating observables for thousands or even millions of parameter sets, making fast and accurate emulators highly desirable or even essential. Emulating scattering across different energies or with interactions such as optical potentials is challenging because of the nonaffine parameter dependence, meaning the parameters do not all factorize from individual operators. Here we introduce and demonstrate the reduced-order scattering emulator ( ROSE ) framework, a reduced basis emulator that can handle nonaffine problems. ROSE is fully extensible and works within the publicly available band framework software suite for calibration, model mixing, and experimental design. As a demonstration problem, we use ROSE to calibrate a realistic nucleon-target scattering model through the calculation of elastic cross sections. This problem shows the practical value of the ROSE framework for Bayesian uncertainty quantification with controlled trade-offs between emulator speed and accuracy as compared to high-fidelity solvers. Finally, planned extensions of ROSE are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Inference of neutrino flavor evolution through data assimilation and neural differential equations

The evolution of neutrino flavor in dense environments such as core-collapse supernovae and binary compact object mergers constitutes an important and unsolved problem. Its solution has potential implications for the dynamics and heavy-element nucleosynthesis in these environments. In this paper, we build upon recent work to explore inference-based techniques for the estimation of model parameters and neutrino flavor evolution histories. Here, we combine data assimilation, ordinary differential equation solvers, and neural networks to craft an inference approach tailored for nonlinear dynamical systems. Using this architecture, and a simple two-neutrino-beam, two-flavor model, we compare the performances of nine different optimization algorithms and expand upon previous assessments of the efficacy of inference for tackling problems in flavor evolution. We find that employing this new architecture, together with evolutionary optimization algorithms, accurately captures flavor histories in the small-scale model and allows us to quickly explore both model parameters and initial flavor content. In future work we plan to extend these inference techniques to large numbers of neutrinos.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Single classical field description of interacting scalar fields

In this report we test the degree to which interacting Bosonic systems can be approximated by a classical field as total occupation number is increased. This is done with our publicly available code repository, QIBS, a new massively parallel solver for these systems. We use a number of toy models well studied in the literature and track when the classical field description admits quantum corrections, called the quantum breaktime. This allows us to test claims in the literature regarding the rate of convergence of these systems to the classical evolution. We test a number of initial conditions, including coherent states, number eigenstates, and field number states. We find that of these initial conditions, only number eigenstates do not converge to the classical evolution as occupation number is increased. We find that systems most similar to scalar field dark matter exhibit a logarithmic enhancement in the quantum breaktime with total occupation number. Systems with contact interactions or with field number state initial conditions, and linear dispersions, exhibit a power law enhancement. Finally, we find that the breaktime scaling depends on both model interactions and initial conditions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Tidal disruption of solitons in self-interacting ultralight axion dark matter

We report ultralight axions (ULAs) are promising dark matter candidates that can have a distinct impact on the formation and evolution of structure on nonlinear scales relative to the cold, collisionless dark matter (CDM) paradigm. However, most studies of structure formation in ULA models do not include the effects of self-interactions, which are expected to arise generically. Here, we study how the tidal evolution of solitons is affected by ULA self-interaction strength and sign. Specifically, using the pseudospectral solver UltraDark.jl, we simulate the tidal disruption of self-interacting solitonic cores as they orbit a 10 11 M ⊙ Navarro-Frenk-White CDM host halo potential for a range of orbital parameters, assuming a fiducial ULA particle mass of 10 -22 eV . We find that repulsive (attractive) self-interactions significantly accelerate (decelerate) soliton tidal disruption. We also identify a degeneracy between the self-interaction strength and soliton mass that determines the efficiency of tidal disruption, such that disruption timescales are affected at the ~ 50 % level for variations in the dimensionless ULA self-coupling from λ = -10 -92 to λ = 10 -92 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Exploring the growth index γ L : Insights from different CMB dataset combinations and approaches

In this study we investigate the growth index γ L , which characterizes the growth of linear matter perturbations, while analysing different cosmological datasets. We compare the approaches implemented by two different patches of the cosmological solver : and __. In our analysis we uncover a deviation of the growth index from its expected Λ CDM value of 0.55 when utilizing the Planck dataset, both in the case and in the __ case, but in opposite directions. This deviation is accompanied by a change in the direction of correlations with derived cosmological parameters. However, the incorporation of cosmic microwave background lensing data helps reconcile γ L with its Λ -cold dark matter value in both cases. Conversely, the alternative ground-based telescopes Atacama Cosmology Telescope and South Pole Telescope consistently yield growth index values in agreement with γ L = 0.55 . We conclude that the presence of the A lens problem in the Planck dataset contributes to the observed deviations, underscoring the importance of additional datasets in resolving these discrepancies. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

Dynamical friction in self-interacting ultralight dark matter

Here, we explore how dynamical friction in an ultralight dark matter (ULDM) background is affected by dark matter self-interactions. We calculate the force of dynamical friction on a point mass moving through a uniform ULDM background with self-interactions, finding that the force of dynamical friction vanishes for sufficiently strong repulsive self-interactions. Using the pseudospectral solver ultraDark.jl, we show with simulations that reasonable values of the ULDM self-interaction strength and particle mass cause $\mathcal{O}$(1) differences in the acceleration of an object like a supermassive black hole (SMBH) traveling near the center of a soliton, relative to the case with no self-interactions. For example, repulsive self-interactions with λ=10 -90 yield a deceleration due to dynamical friction ≈ 70% smaller than a model with no self-interactions. We discuss the observational implications of our results for SMBHs near soliton centers and for massive satellite galaxies falling into ultralight axion halos and show that outcomes are dependent on whether a self-interaction is present or not.

79 ASTRONOMY AND ASTROPHYSICS↗