Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RANDOM PHASE 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.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Effects of quasiparticle-vibration coupling on Gamow-Teller strength and 𝛽 decay with the Skyrme proton-neutron finite-amplitude method

We adapt the proton-neutron finite-amplitude method, which in its original form is an efficient implementation of the Skyrme quasiparticle random phase approximation, to include the coupling of quasiparticles to like-particle phonons. The approach allows us to add beyond-quasiparticle random-phase approximation correlations to computations of Gamow-Teller strength and 𝛽-decay rates in deformed nuclei for the first time. We test the approach in several deformed isotopes for which measured strength distributions are available. Furthermore, the additional correlations dramatically improve agreement with the data, and will lead to improved global 𝛽-decay rates.

Beta decay↗
DOE OSTI2022

Coupled Hubbard ladders at weak coupling: Pairing and spin excitations

The Hubbard model provides a simple framework in which one can study how certain aspects of the electronic structure of strongly interacting systems can be tuned to optimize the superconducting pairing correlations and how these changes affect the mechanisms giving rise to them. Here we use a weak-coupling random phase approximation to study a two-dimensional Hubbard model with a unidirectional modulation of the hopping amplitudes as the system evolves from the uniform square lattice to an array of weakly coupled two-leg ladders. We find that the pairing correlations retain their dominant d x 2 - y 2- wavelike structure and that they are significantly enhanced for a slightly modulated lattice. This enhancement is traced backed to an increase in the strength of the spin-fluctuation pairing interaction due to favorable Fermi surface nesting in the modulated system. We then use a random-phase approximation BCS framework to examine the evolution of the neutron resonance in the superconducting state. We find that it changes only weakly for moderate modulations, but breaks up into two distinct resonances at incommensurate wave vectors in the limit of weakly coupled ladders.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Magnetic correlations and pairing tendencies of the hybrid stacking nickelate superlattice La 7 Ni 5 O 17 (La 3 ⁢Ni 2 ⁢O 7 /La 4 Ni 3 ⁢O 10 ) under pressure

Motivated by the recent rapid progress in high-𝑇 𝑐 nickelate superconductors, we comprehensively study the physical properties of the alternating bilayer trilayer stacking nickelate La 7 ⁢Ni 5 ⁢O 17 . The high-symmetry phase of this material, without the tilting of oxygen octahedra, is not stable at ambient conditions but becomes stable under high pressure, where a small hole pocket 𝛾 0 , composed of the 𝑑 3⁢𝑧 2 −𝑟 2 states in the trilayer sublattice, appears. Here, this pocket was identified in our previous work for trilayer La 4 ⁢Ni 3 ⁢O 10 as important to develop superconductivity. Moreover, using random-phase approximation calculations, we find a leading 𝑠 ± pairing state for the high-symmetry phase under pressure with similar pairing strength as that obtained previously for the bilayer La 3⁢ Ni 2 ⁢O 7 compound, suggesting a similar or higher superconducting transition temperature 𝑇 𝑐 , at the random-phase approximation level. In addition, we find that the dominant magnetic fluctuations in the system driving this pairing state have antiferromagnetic structure both in-plane and between the planes of the top and bottom trilayer and bilayer sublattices, while the middle trilayer is magnetically decoupled.

Zhang, Yang [Univ. of Tennessee, Knoxville, TN (Un↗
DOE OSTI2024

Magnetic phase diagram of a two-orbital model for bilayer nickelates with varying doping

Motivated by the recently discovered high-T c bilayer nickelate superconductor La 3⁢ Ni 2 ⁢O 7 , we comprehensively research a bilayer 2×2×2 cluster for different electronic densities n by using the Lanczos method. We also employ the random-phase approximation to quantify the first magnetic instability with increasing Hubbard coupling strength, also varying n. Based on the spin structure factor S(q), we have obtained a rich magnetic phase diagram in the plane defined by n and U/W, at fixed Hund coupling, where U is the Hubbard strength and W the bandwidth. We have observed numerous states, such as A-AFM, Stripes, G-AFM, and C-AFM. At half-filling, n=2 (two electrons per Ni site, corresponding to N=16 electrons), the canonical superexchange interaction leads to a robust G-AFM state (π,π,π) with antiferromagnetic couplings both in-plane and between layers. By increasing or decreasing electronic densities, ferromagnetic tendencies emerge from the “half-empty” and “half-full” mechanisms, leading to many other interesting magnetic tendencies. In addition, the spin-spin correlations become weaker both in the hole or electron doping regions compared with half-filling. At n=1.5 (or N=12), density corresponding to La 3 ⁢Ni 2 ⁢O 7 , we obtained the “Stripe 2” ground state (antiferromagnetic coupling in one in-plane direction, ferromagnetic coupling in the other, and antiferromagnetic coupling along the z axis) in the 2×2×2 cluster. In addition, we obtained a much stronger AFM coupling along the z axis than the magnetic coupling in the xy plane. The random-phase approximation calculations with varying n give very similar results as Lanczos, even though both techniques are based on quite different procedures. Additionally, a state with q/π=(0.6,0.6,1) close to the E-phase wavevector is found in our RPA calculations by slightly reducing the filling to n=1.25, possibly responsible for the E-phase SDW recently observed in experiments. In conclusion, our predictions can be tested by chemically doping La 3 ⁢Ni 2 ⁢O 7 .

36 MATERIALS SCIENCE↗
DOE OSTI2024

Magnetic dipole γ-ray strength functions in the crossover from spherical to deformed neodymium isotopes

We calculate the magnetic dipole $\gamma$-ray strength functions in a chain of even-mass neodymium isotopes $^{144-152}$Nd in the framework of the configuration-interaction (CI) shell model. We infer the strength function by applying the maximum entropy method (MEM) to the exact imaginary-time response function calculated with the shell-model Monte Carlo (SMMC) method. The success of the MEM depends on the choice of a good strength function as a prior distribution. We investigate two choices for the prior strength function: the static path approximation (SPA) and the quasiparticle random-phase approximation (QRPA). We find that the QRPA is a better approximation at low temperatures (i.e., near the ground state), while the SPA is a better choice at finite temperatures. We identify a low-energy enhancement (LEE) in the MEM deexcitation $M1$ strength functions of the even-mass neodymium isotopes and compare with recent experimental results for the total deexcitation $\gamma$-ray strength functions. The LEE is already seen in the SPA strength function but not in the QRPA strength function, indicating the importance of large-amplitude static fluctuations around the mean field in reproducing the LEE. Our method is currently the only one which can reproduce LEE in heavy open-shell nuclei where conventional CI shell model calculations are prohibited. With the onset of deformation as number of neutrons increases along the chain of neodymium isotopes, we observe that some of the LEE strength transfers to a low-energy excitation, which we interpret as a finite-temperature ``scissors'' mode. Here, we also observe a finite-temperature spin-flip mode.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Spin-crossover complexes: Self-interaction correction vs density correction

Complexes containing a transition metal atom with a 3d 4 –3d 7 electron configuration typically have two low-lying, high-spin (HS) and low-spin (LS) states. The adiabatic energy difference between these states, known as the spin-crossover energy, is small enough to pose a challenge even for electronic structure methods that are well known for their accuracy and reliability. In this work, we analyze the quality of electronic structure approximations for spin-crossover energies of iron complexes with four different ligands by comparing energies from self-consistent and post-self-consistent calculations for methods based on the random phase approximation and the Fermi–Löwdin self-interaction correction. Considering that Hartree–Fock densities were found by Song et al., J. Chem. Theory Comput. 14, 2304 (2018), to eliminate the density error to a large extent, and that the Hartree–Fock method and the Perdew–Zunger-type self-interaction correction share some physics, we compare the densities obtained with these methods to learn their resemblance. Here, we find that evaluating non-empirical exchange-correlation energy functionals on the corresponding self-interaction-corrected densities can mitigate the strong density errors and improves the accuracy of the adiabatic energy differences between HS and LS states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Scattering evidence of positional charge correlations in polyelectrolyte complexes

Polyelectrolyte complexation plays an important role in materials science and biology. The internal structure of the resultant polyelectrolyte complex (PEC) phase dictates properties such as physical state, response to external stimuli, and dynamics. Small-angle scattering experiments with X-rays and neutrons have revealed structural similarities between PECs and semidilute solutions of neutral polymers, where the total scattering function exhibits an Ornstein–Zernike form. In spite of consensus among different theoretical predictions, the existence of positional correlations between polyanion and polycation charges has not been confirmed experimentally. Here, we present small-angle neutron scattering profiles where the polycation scattering length density is matched to that of the solvent to extract positional correlations among anionic monomers. The polyanion scattering functions exhibit a peak at the inverse polymer screening radius of Coulomb interactions, q* ≈ 0.2 Å –1 . This peak, attributed to Coulomb repulsions between the fragments of polyanions and their attractions to polycations, is even more pronounced in the calculated charge scattering function that quantifies positional correlations of all polymer charges within the PEC. Screening of electrostatic interactions by adding salt leads to the gradual disappearance of this correlation peak, and the scattering functions regain an Ornstein–Zernike form. Experimental scattering results are consistent with those calculated from the random phase approximation, a scaling analysis, and molecular simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Daejeon16 interaction with contact-term corrections for heavy nuclear systems

The Daejeon16 two-nucleon interaction is employed in many-body approaches based on the mean-field approximation. The perturbative character of Daejeon16 is verified by comparing results for 16 O from the Hartree–Fock (HF) approximation and from the no-core shell model and by examining the magnitude of perturbative corrections to the HF energy in light and heavy nuclei. Here, in order to approximately describe energies and radii across the nuclear chart, a phenomenological correction in the form of a two-plus-three-nucleon contact interaction is introduced. With fitted parameters we achieve a very good description of medium-mass nuclei in terms of energy and size and also in terms of the centroid energy of the giant monopole resonance and the dipole polarizability calculated within the random-phase approximation. Our results provide further justification for the use of Daejeon16 augmented with phenomenological corrections as an effective interaction of perturbative character in a variety of applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Surrogate Hessian accelerated structural optimization for stochastic electronic structure theories

In this work, we present an efficient energy-based method for structural optimization with stochastic electronic structure theories, such as diffusion quantum Monte Carlo (DMC). This method is based on robust line-search energy minimization in reduced parameter space, exploiting approximate but accurate Hessian information from a surrogate theory, such as density functional theory. The surrogate theory is also used to characterize the potential energy surface, allowing for simple but reliable ways to maximize statistical efficiency while retaining controllable accuracy. We demonstrate the method by finding the minimum DMC energy structures of the selected flake-like aromatic molecules, such as benzene, coronene, and ovalene, represented by 2, 6, and 19 structural parameters, respectively. In each case, the energy minimum is found within two parallel line-search iterations. The method is near-optimal for a line-search technique and suitable for a broad range of applications. It is easily generalized to any electronic structure method where forces and stresses are still under active development and implementation, such as diffusion Monte Carlo, auxiliary-field Monte Carlo, and stochastic configuration interaction, as well as deterministic approaches such as the random-phase approximation. Accurate and efficient means of geometry optimization could shed light on a broad class of materials and molecules, showing high sensitivity of induced properties to structural variables.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Weak-coupling theory of neutron scattering as a probe of altermagnetism

Inelastic neutron scattering provides a powerful probe of the magnetic excitations of quantum magnets. Altermagnets have recently emerged as a new class of magnets with vanishing net magnetization characteristic of antiferromagnets and with a spin-split electronic structure typical of ferromagnets. Here we introduce a minimal Hubbard model with two-sublattice orthorhombic anisotropy as a framework to study altermagnetism. Using unrestricted Hartree-Fock calculations, we find an altermagnetic state for this model that evolves from a metallic state to an insulating state with increasing Hubbard-U Coulomb repulsion. We then examine the inelastic neutron scattering response in these states using random-phase approximation calculations of the dynamic spin susceptibility χ"(q,ω). We find that the magnetic excitation spectrum depends on its chirality for q along certain directions in reciprocal space, an observation that may be used in inelastic neutron scattering experiments as a probe of altermagnetism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Magnetic ground state and perturbations of the distorted kagome Ising metal TmAgGe

Here, we present the magnetic orders and excitations of the distorted kagome intermetallic magnet TmAgGe. Using neutron single crystal diffraction we identify the propagation vectors k = ($\frac{1}{2}$ 0 0) and k = (0 0 0) and determine the magnetic structures of the zero-field and magnetic field-induced phases for H along the $\mathcal{a}$ and [–110] crystal directions. We determine the experimental magnetic field- temperature (H, T)-phase diagram and reproduce it by Monte Carlo simulations of an effective spin exchange Hamiltonian for one distorted kagome layer. Our model includes a strong axial single-ion anisotropy and significantly smaller exchange couplings, which span up to the third-nearest neighbors within the layer. Single crystal inelastic neutron scattering (INS) measurements reveal an almost flat, only weakly dispersive mode around 7 meV that we use alongside bulk magnetization data to deduce the crystal-electric field (CEF) scheme for the Tm 3+ ions. Random phase approximation (RPA) calculations based on the determined CEF wave functions of the two lowest quasidoublets enable an estimation of the interlayer coupling that is compatible with the experimental INS spectra. No evidence for low-energy spin waves associated to the magnetic order was found, which is consistent with the strongly Ising nature of the ground state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Weak decays in superheavy nuclei

Superheavy nuclei represent the extreme atoms and nuclides known at the limit of mass and charge. The observed superheavy nuclei are all proton-rich; they decay primarily by emitting 𝛼 particles and by fission with a possible small electron capture (EC) branch. Here, due to the huge atomic numbers and associated relativistic effects, EC decays of superheavy systems are expected to differ from what is known in lighter nuclei. In this letter, using the quantified relativistic nuclear density functional theory and the quasiparticle random-phase approximation with the interaction optimized to experimental EC/𝛽 ± -decay half-lives, and Gamow-Teller resonance energies, we study the EC/𝛽 ± -decays in 𝑍=101–118 nuclei. Both allowed (1 + ) and first-forbidden (0 − ,1 − and 2 − ) transitions are considered. We show that the first-forbidden 1 − transitions dominate the decay rates in almost all studied nuclei. For proton-rich nuclei, EC dominates over 𝛽 + decay. Based on calculations with two relativistic energy density functionals, we identify 45 candidate nuclei in which a competition between weak decays and 𝛼 decay and spontaneous fission is expected.

A ≥ 220↗
DOE OSTI2021

Isoscalar and isovector giant resonances in 44 Ca, 54 Fe, 64,68 Zn and 56,58,60,68 Ni

We have studied the uncharacteristic behavior of the measured values of the isoscalar and isovector centroid energies, E CEN , of nuclear giant resonances of multipolarity from L = 0 to L = 3 in 44 Ca, 54 Fe, 64,68 Zn and 56,58,60,68 Ni. For this purpose, we carried out calculations of E CEN within the spherical Hartree-Fock (HF)-based random phase approximation (RPA) theory with 33 distinct Skyrme-like effective nucleon-nucleon interactions. We have also determined the Pearson linear correlation coefficients between centroid energies, obtained from the HF-RPA, and the various properties of nuclear matter (NM) of each interaction and determined the sensitivity of E CEN to NM properties. We compared the theoretical values of E CEN obtained from the HF-RPA calculations with experimental data and discuss the results, pointing out significant disagreements between theoretical and experimental values. We note in particular, that we obtain good agreement for the theoretical E CEN of the isovector giant dipole resonance and the available experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Electronic structure, self-doping, and superconducting instability in the alternating single-layer trilayer stacking nickelates La 3 ⁢Ni 2 ⁢O 7

Motivated by the recently proposed alternating single-layer trilayer stacking structure for the nickelate La 3 ⁢Ni 2 ⁢O 7 , we comprehensively study this system using ab initio and random-phase approximation techniques. Here, our analysis unveils similarities between this novel La 3⁢ Ni 2 ⁢O 7 structure and other Ruddlesden-Popper nickelate superconductors, such as a similar charge-transfer gap value and orbital-selective behavior of the eg orbitals. Pressure primarily increases the bandwidths of the Ni e g bands, suggesting an enhancement of the itinerant properties of those e g states. By changing the cell volume ratio V/V 0 from 0.9 to 1.10, we found that the bilayer structure in La 3 ⁢Ni 2 ⁢O 7 always has lower energy than the single-layer trilayer stacking La 3 ⁢Ni 2 ⁢O 7 . In addition, we observe a “self-doping” effect (compared to the average 1.5 electrons per eg orbital per site of the entire structure) from the trilayer to the single-layer sublattices and this effect will be enhanced by overall electron doping. Moreover, we find a leading d x 2 -y 2 -wave pairing state that is restricted to the single layer. Because the effective coupling between the single layers is very weak, due to the nonsuperconducting trilayer in-between, this suggests that the superconducting transition temperature T c in this structure should be much lower than in the bilayer structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Surrogate models for linear response

Linear response theory is a well-established method in physics and chemistry for exploring excitations of many-body systems. In particular, the quasiparticle random-phase approximation (QRPA) provides a powerful microscopic framework by building excitations on top of the mean-field vacuum; however, its high computational cost limits model calibration and uncertainty quantification studies. Here, we present two complementary QRPA surrogate models and apply them to study response functions of finite nuclei. One is a reduced-order model that exploits the underlying QRPA structure, while the other utilizes the recently developed parametric matrix model algorithm to construct a map between the system’s Hamiltonian and observables. Our benchmark applications, the calculation of the electric dipole polarizability of 180 Yb and the 𝛽-decay half-life of 80 Ni, show that both emulators can achieve 0.1%–1% accuracy while offering a 6–7 orders of magnitude speedup compared to state-of-the-art QRPA solvers. These results demonstrate that the developed QRPA emulators are well positioned to enable Bayesian calibration and large-scale studies of computationally expensive physics models describing the properties of many-body systems.

Beta decay↗
DOE OSTI2021

Plasmons in two-dimensional topological insulators

Here, we analyze collective excitations in models of two-dimensional topological insulators using the random phase approximation. In a two-dimensional extension of the Su-Schrieffer-Heeger model, edge plasmonic excitations with induced charge-density distributions localized at the boundaries of the system are found in the topologically nontrivial phase, dispersing similarly as one-dimensional bulk plasmons in the conventional Su-Schrieffer-Heeger chain. For two-dimensional bulk collective modes, we reveal regimes of enhanced interband wave function correlations, leading to characteristic hardening and softening of inter- and intraband bulk plasmonic branches, respectively. In the two-dimensional Haldane Chern insulator model, chiral, unidirectional edge plasmons in nano-ribbon architectures are observed, which can be characterized by an effective Coulomb interaction cross section. Bulk collective excitations in the two-dimensional Haldane model are shown to be originated by single-particle band structure details in different topological phases.

2-dimensional systems↗
DOE OSTI2024

Prediction of s ± -Wave Superconductivity Enhanced by Electronic Doping in Trilayer Nickelates La 4 Ni 3 O 10 under Pressure

Motivated by the recently reported signatures of superconductivity in trilayer La 4 ⁢Ni 3 ⁢O 10 under pressure, here we comprehensively study this system using ab initio and random-phase approximation techniques. Without electronic interactions, the Ni d 3z 2 –r 2 orbitals show a bonding-antibonding and nonbonding splitting behavior via the O p z orbitals inducing a “trimer” lattice in La 4⁢ Ni 3 ⁢O 10 , analogous to the dimers of La 3 ⁢Ni 2 ⁢O 7 . The Fermi surface consists of three electron sheets with mixed e g orbitals, and a hole and an electron pocket made up of the d 3⁢z 2 –r 2 orbital, suggesting a Ni two-orbital minimum model. In addition, we find that superconducting pairing is induced in the s ± -wave channel due to partial nesting between the M = (π,π) centered pockets and portions of the Fermi surface centered at the Γ = (0,0) point. With changing electronic density n, the s ± instability remains leading and its pairing strength shows a domelike behavior with a maximum around n = 4.2 ( ~6.7% electron doping). The superconducting instability disappears at the same electronic density as that in the new 1313 stacking La 3 ⁢Ni 2 ⁢O 7 , correlated with the vanishing of the hole pocket that arises from the trilayer sublattice, suggesting that the high-T c superconductivity of La 3⁢ Ni 2 ⁢O 7 does not originate from a trilayer and monolayer structure. Furthermore, we confirm the experimentally proposed spin state in La 4 ⁢Ni 3 ⁢O 10 with an in-plane (π, π) order and antiferromagnetic coupling between the top and bottom Ni layers, and spin zero in the middle layer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Plasmons in Z 2 topological insulators

Here, we study plasmonic excitations in the Kane-Mele model, a two-dimensional Z 2 topological insulator on the honeycomb lattice, using the random phase approximation (RPA). In the topologically nontrivial phase, the model has conducting edge states that traverse the bulk energy gap and display spin-momentum locking. Such a state of matter is called the quantum spin hall (QSH) phase, which is robust against time-reversal (TR) invariant perturbations. We find that in the QSH phase, gapless spin-polarized plasmons can be excited on the edges of the system. The propagation of these plasmons is chiral for each individual spin component and shows spin-momentum locking for both spin components on the same edge. Moreover, we study the effect of external magnetic fields on the gapless edge plasmons. Specifically, out-of-plane magnetic fields delocalize edge plasmons propagating in one direction without affecting the other one, while an in-plane magnetic field can be applied to selectively excite a specific spin-plasmon branch with proper doping or gating to the system. Our findings may have potential applications in novel plasmonic and spintronic devices. We also investigate plasmons in the Kane-Mele model on a finite-sized diamond-shaped nanoflake and observe low-energy plasmons circulating the boundary of the material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗