Search NASA⌕ Search

SEARCH · Search NASA

Results for “self-consistent calculations”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Transport barrier and spinning blob dynamics in the tokamak edge

Abstract In this work, we investigate the dynamics of plasma blobs in the edge of magnetic confinement devices using a full-f gyrokinetic particle-in-cell code with X-point geometry. In simulations, the evolution of a seeded blob is followed as it approaches a naturally-forming zonal shear layer near the separatrix, where the blob is stabilized by a large spin induced by the self-consistent adiabatic electron response, and blob bifurcation and trapping are observed during the cross-field propagation of blobs. A new theoretical explanation in both the zonal free and zonal shear layer is constructed, where the dominant E × B spin motion is included. A theoretical condition for a transport barrier induced by the interaction between spinning blobs and the zonal shear layer is obtained, and its scaling is verified with simulations. The new theoretical framework, especially the transport barrier, is applicable to explain and predict various experimental phenomena. In particular, the transport barrier condition calculated with experimental parameters demonstrates that the blob radial transport for H mode is smaller than L mode in experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Quasiparticle and Optical Properties of Carrier-Doped Monolayer MoTe 2 from First Principles

The intrinsic weak and highly nonlocal dielectric screening of two-dimensional materials is well-known to lead to high sensitivity of their optoelectronic properties to environment. Less studied theoretically is the role of free carriers in those properties. In this report we use ab initio GW and Bethe-Salpeter equation calculations, with a rigorous treatment of dynamical screening and local-field effects, to study the doping dependence of the quasiparticle and optical properties of a monolayer transition-metal dichalcogenide, 2H MoTe 2 . Here we predict a quasiparticle band gap renormalization of several hundreds of meV for experimentally attainable carrier densities and a similarly sizable decrease in the exciton binding energy. This results in an almost constant excitation energy for the lowest-energy exciton resonance with an increasing doping density. Using a newly developed and generally applicable plasmon-pole model and a self-consistent solution of the Bethe-Salpeter equation, we reveal the importance of accurately capturing both dynamical and local-field effects to understand detailed photoluminescence measurements.

2D material↗
DOE OSTI2025

Self-consistent modeling of tokamak edge plasma transport with lithium sources

Magnetic confinement fusion devices require effective heat and particle exhaust solutions on the divertor plates to operate sustainably, especially under reactor-relevant conditions. Liquid lithium divertors have been proposed to address two major challenges: control of excessive heat flux to plasma-facing components through vapor shielding and minimization of core plasma contamination from impurities. The National Spherical Torus Experiment-Upgrade (NSTX-U) will explore lithium as a divertor material due to its potential to meet both objectives. We present a self-consistent coupling framework between the plasma boundary transport code UEDGE and the lithium wall transport code Wall–Li to evaluate the feasibility and operational limits of lithium-based divertors. The model aims to optimize lithium sourcing levels to prevent core plasma contamination via fuel dilution while ensuring divertor protection through vapor shielding. This integrated framework, applicable to any tokamak with lithium sources, dynamically adjusts lithium sourcing based on local plasma conditions and surface temperature. The coupled model is tested using NSTX-like geometry and plasma conditions to assess its performance and reliability. Wall–Li calculates lithium fluxes from plasma-facing components, incorporating physical sputtering, thermally enhanced sputtering, and evaporation driven by surface temperature and ion flux. These fluxes are reintroduced into UEDGE as neutral lithium atoms, enabling simulation of their transport and distribution within the plasma. UEDGE computes plasma and neutral transport, surface heat flux, and iteratively feeds this information back to Wall–Li. A small time step is employed to ensure numerical stability and convergence, enabling accurate simulations over typical tokamak discharge durations. This integrated modeling approach provides a robust tool for identifying operational regimes that balance effective lithium sourcing with minimal core plasma contamination, offering critical insights for optimizing lithium-based divertor systems in current and future fusion devices.

Magnetic confinement fusion↗
DOE OSTI2024

Impact of beam–beam effects on absolute luminosity calibrations at the CERN Large Hadron Collider

At the Large Hadron Collider (LHC), absolute luminosity calibrations obtained by the van der Meer (vdM) method are affected by the mutual electromagnetic interaction of the two beams. The colliding bunches experience relative orbit shifts, as well as optical distortions akin to the dynamic-β effect, that both depend on the transverse beam separation and must therefore be corrected for when deriving the absolute luminosity scale. In the vdM regime, the beam–beam parameter is small enough that the orbit shift can be calculated analytically. The dynamic-β corrections to the luminometer calibrations, however, had until the end of Run 2 been estimated in the linear approximation only. In this report, the influence of beam–beam effects on the vdM-based luminosity scale is quantified, together with the associated systematic uncertainties, by means of simulations that fully take into account the non-linearity of the beam–beam force, as well as the resulting non-Gaussian distortions of the transverse beam distributions. Two independent multiparticle simulations, one limited to the weak-strong approximation and one that models strong-strong effects in a self-consistent manner, are found in excellent agreement; both predict a percent-level shift of the absolute pp-luminosity values with respect to those assumed until recently in the physics publications of the LHC experiments. These results also provide guidance regarding further studies aimed at reducing the beam–beam-related systematic uncertainty on beam–beam corrections to absolute luminosity calibrations by the van der Meer method.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Element Formation in Radiation-hydrodynamics Simulations of Kilonovae

Abstract Understanding the details of r -process nucleosynthesis in binary neutron star merger (BNSM) ejecta is key to interpreting kilonova observations and identifying the role of BNSMs in the origin of heavy elements. We present a self-consistent, two-dimensional, ray-by-ray radiation-hydrodynamic evolution of BNSM ejecta with an online nuclear network (NN) up to a timescale of days. For the first time, an initial numerical relativity ejecta profile composed of the dynamical component and spiral-wave and disk winds is evolved including detailed r -process reactions and nuclear heating effects. A simple model for the jet energy deposition is also included. Our simulation highlights that the common approach of relating in postprocessing the final nucleosynthesis yields to the initial thermodynamic profile of the ejecta can lead to inaccurate predictions. Moreover, we find that neglecting the details of the radiation-hydrodynamic evolution of the ejecta in nuclear calculations can introduce deviations of up to 1 order of magnitude in the final abundances of several elements, including very light and second r -process peak elements. The presence of a jet affects element production only in the innermost part of the polar ejecta, and it does not alter the global nucleosynthesis results. Overall, our analysis shows that employing an online NN improves the reliability of nucleosynthesis and kilonova light-curve predictions.

Magistrelli, Fabio (ORCID:0009000509767851)↗
DOE OSTI2025

Electronic structure and exchange interactions in altermagnetic MnGeP 2 in the quasiparticle self-consistent G W approach

The quasiparticle self-consistent 𝐺⁡𝑊 method is used to study the electronic band structure, optical dielectric function, and exchange interactions in chalcopyrite, 𝐼⁢ ‾ 4 ⁢2⁢𝑑, structure MnGeP 2 . The material is found to be an antiferromagnetic semiconductor with a lowest direct gap of 2.44 eV at the Γ point and a lower indirect gap of 1.87 eV from Γ to 𝑀. The material is an altermagnet because the two magnetic atoms of opposite spin are related by a twofold rotation operation perpendicular to the main fourfold rotation inversion axis. The spin splittings along a low symmetry line like 𝑃⁢𝑁 is sizable while at k points on the diagonal mirror planes or on the twofold symmetry axes the spin splitting is zero. The exchange interactions are calculated using a linear response approach. The antiferromagnetic exchange interaction between nearest neighbors in the primitive unit cell is dominating and found to be slightly decreasing upon carrier doping but not sufficiently to change the interaction to become ferromagnetic. The bare (noninteracting) and interacting transverse spin susceptibilities, which provide interatomic site exchange interactions after averaging over the muffin-tin spheres, are calculated from the 𝐺⁡𝑊 band structure and wave functions. From these exchange interactions, the spin wave spectra are obtained along the high symmetry lines and the Néel temperature is calculated using the mean-field and Tyablikov estimations. The dielectric function and the optical absorption spectra are calculated including excitonic effects using the Bethe Salpeter equation. The exchange interactions around Mn Ge defect sites is also studied. While we find it can generate ferromagnetic interactions with neighboring spins, we did not find direct evidence of producing an overall ferromagnetic phase. First, if Mn antisites are introduced by exchanging Mn with a nearby Ge, the interactions stay largely antiferromagnetic. Second, when we add additional Mn, in other words in Mn-rich stoichiometry, the Mn Ge antisites produce a strong ferromagnetic interaction primarily with the Mn in the same basal plane but weaker ferromagnetic interaction with adjacent plane Mn. The interactions between regular lattice Mn stay antiferromagnetic as before and thus favor keeping the antiferromagnetic order along the [001] direction. Adding Mn antisites, however, does lead to a metallic band structure.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH) 4 H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ΔE estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ΔE is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ~5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

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

Towards a self-consistent model of the convective core boundary in upper main sequence stars: I. 2.5D and 3D simulations

There is strong observational evidence that the convective cores of intermediate-mass and massive main sequence stars are substantially larger than those predicted by standard stellar-evolution models. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-timescale process that is likely to be particularly relevant during the slow evolution on the main sequence. We use our low-Mach-number S EVEN -L EAGUE H YDRO code to study this process in 2.5D and 3D geometries. Starting with a chemically homogeneous model of a 15 M ⊙ zero-age main sequence star, we construct a series of simulations with the luminosity increased and opacity decreased by the same factor, ranging from 10 3 to 10 6 . After reaching thermal equilibrium, all of our models show a clear penetration layer; its thickness becomes statistically constant in time and it is shown to converge upon grid refinement. The penetration layer becomes nearly adiabatic with a steep transition to a radiative stratification in simulations at the lower end of our luminosity range. This structure corresponds to the adiabatic ‘step overshoot’ model often employed in stellar-evolution calculations. The simulations with the highest and lowest luminosity differ by less than a factor of two in the penetration distance. The high computational cost of 3D simulations makes our current 3D data set rather sparse. Depending on how we extrapolate the 3D data to the actual luminosity of the initial stellar model, we obtain penetration distances ranging from 0.09 to 0.44 pressure scale heights, which is broadly compatible with observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Mutually magnetically insulated two-species Brillouin flow

In this work, we use the Brillouin flow analytic framework to examine the physics of Magnetically Insulated Transmission Lines (MITL). Here, we derive a model applicable to any particle species, including both positive and negative ions, in planar and cylindrical configurations. We then show how to self-consistently solve for two-species simultaneously, using magnetically insulated electrons and positive ions as an example. We require both layers to be spatially separated and magnetically insulated (mutually magnetically insulated); for a 7.5 cm gap with a 2 MV bias voltage, this condition requires magnetic fields in excess of 2.73 T. We see a close match between mutually insulated MITL performance and “superinsulated” (high degree of magnetic insulation) electron-only theory, as may be expected for these high magnetic fields. However, the presence of ions leads to several novel effects: (1) Opposite to electron-only theory, total electron currents increase rather than decrease as the degree of magnetic insulation becomes stronger. The common assumption of neglecting electrons for superinsulated MITL operation must be revisited when ions are present—we calculate up to 20× current enhancement. (2) The electron flow layer thickness increases up to double, due to ion space-charge enhancement. (3) The contributions from both ions and electrons to the MITL flow impedance are calculated. The flow impedance drops by over 50% when ions fill the gap, which can cause significant reflections at the load if not anticipated and degrade performance. Additional effects and results from the inclusion of the ion layer are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Self-Consistent Theory for Structural Relaxation, Dynamic Bond Exchange Times, and the Glass Transition in Polymeric Vitrimers

Here, we formulate a statistical mechanical theory for how dynamic bond exchange influences the activated hopping-driven relaxation of Kuhn segments in dynamically cross-linked networks or vitrimers over a wide range of temperatures and cross-link densities. The key new methodological aspect is to address in a self-consistent manner the dynamic consequences of bond exchange on the Kuhn segmental alpha relaxation, and vice versa. The predicted temperature dependence of the segmental alpha time of vitrimers at high temperatures remains the same as that of permanent networks, but at lower temperatures, a significant acceleration of relaxation occurs due to bond exchanges. From a mechanistic perspective, the vitrimer local cage barrier is very weakly affected by bond exchange, while the collective elastic barrier contribution decreases significantly in the deeply supercooled regime. The vitrimer glass transition temperature is predicted to grow linearly with the square root of the cross-link density, as previously found for permanent networks. Material-specific chemical effects such as cross-linker size relative to that of the normal Kuhn segment or special attraction of a cross-linker with polymers are crudely considered based on model calculations. The theory is quantitatively applied to recent experiments on dry ethylene vitrimers. Good agreements are found including that the bond exchange time follows an Arrhenius law at high enough temperatures but upward non-Arrhenius deviations emerge in the deeply supercooled regime, and a collapsed master curve of the segment alpha time exists over a wide range of cross-link densities and temperatures. Possible extensions to treat dynamic heterogeneity effects and penetrant transport in vitrimers are briefly discussed.

chemical calculations↗
DOE OSTI2023

Solving the electronic structure problem for over 100000 atoms in real space

Using a real-space high-order finite-difference approach, we investigate the electronic structure of large spherical silicon nanoclusters. Within Kohn-Sham density functional theory and using pseudopotentials, we report the self-consistent field convergence of a system with over 100000 atoms: a Si 107,641 ⁢H 9,084 nanocluster with a diameter of 16 nm. Our approach uses Chebyshev-filtered subspace iteration to speed up the convergence of the eigenspace, and blockwise Hilbert space-filling curves to speed up sparse matrix-vector multiplications, all of which are implemented in the parsec code. For the largest system, we utilized 2048 nodes (114 688 cores) on the Frontera machine in the Texas Advanced Computing Center. Our quantitative analysis of the electronic structure shows how it gradually approaches its bulk counterpart as a function of nanocluster size. The band gap is enlarged due to quantum confinement in nanoclusters, but decreases as the system size increases, as expected. In conclusion, our work serves as a proof of concept for the capacity of the real-space approach in efficiently parallelizing very large calculations using high-performance computer platforms, which can straightforwardly be replicated in other systems with more than 10 5 atoms.

0-dimensional systems↗
DOE OSTI2024

Exploring the effect of ELM and code-coupling frequencies on plasma and material modeling of dynamic recycling in divertors

Abstract Integrated modeling of plasma-surface interactions provides a comprehensive and self-consistent description of the system, moving the field closer to developing predictive and design capabilities for plasma facing components. One such workflow, including descriptions for the scrape-off-layer plasma, ion-surface interactions and the sub-surface evolution, was previously used to address steady-state scenarios and has recently been extended to incorporate time-dependence and two-way information flow. The new model can address dynamic recycling in transient scenarios, such as the application presented in this paper: the evolution of W samples pre-damaged by helium and exposed to ELMy H-mode plasmas in the DIII-D DiMES. A first set of simulations explored the effect of ELM frequency. This study was discussed in detail in this conference’s proceedings and is summarized here. The 2nd set of simulations, which is the focus of this paper, explores the effect of code-coupling frequency. These simulations include initial SOLPS solutions converged to the inter-ELM state, ion impact energy ( E in ) and angles ( A in ) calculated by hPIC2, and an improved heat transfer description in Xolotl. The model predicts increases in particle fluxes and decreases in heat fluxes by 10%–20% with the coupling time-step. Compared with the first set of simulations, the less shallow impact angle leads to smaller reflection rates and significant D implantation. The higher fraction of implanted flux (and deeper), in particular during ELMs, increases the accumulated D content in the W near-surface region. Future expansion of the workflow includes coupling to hPIC2 and GITR to ensure accurate descriptions of E in and A in , and W impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Sensitivities of time-dependent temperature profile predictions for NSTX with the multi-mode model

The Multi-Mode Model (MMM) for turbulent transport was applied to a large set of well-analyzed discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate its sensitivities to a wide range of plasma conditions. MMM calculations were performed for hundreds of milliseconds in each discharge by performing time-dependent predictive simulations with the 1.5D tokamak integrated modeling code TRANSP. A closely related study (Lestz et al 2025 Plasma Phys. Control. Fusion 67 105029) concluded that MMM predicted electron and ion temperature profiles that were in reasonable agreement with NSTX observations, generally outperforming a different reduced transport model, TGLF. This finding motivates the more thorough investigation of the characteristics of the MMM predictions conducted in this work. The simulations with MMM have electron energy transport dominated by electron temperature gradient modes in the examined discharges with relatively low plasma β (ratio of kinetic plasma pressure to magnetic field pressure) and high collisionality, transitioning to a mixture of different modes for higher β and lower collisionality. The thermal ion diffusivity predicted by MMM is much smaller than the neoclassical contribution, in line with previous experimental analysis of NSTX. Nonetheless, the electron and ion temperature profiles are coupled via collisional energy exchange and thus sensitive to which transport channels are predicted. The time-dependent simulations with MMM are robust to the simulation start time, converging to remarkably similar temperature profiles later during the discharge. MMM typically overpredicts confinement relative to NSTX observations, leading to the prediction of overly steep temperature profiles. Plasmas with spatially broader temperature profiles, higher plasma β, and longer energy confinement times tend to be predicted by MMM with better agreement with the experiment. As a result, these findings provide useful context for understanding the regime-dependent tendencies of MMM in anticipation of self-consistent, time-dependent predictive simulations of NSTX-U discharges with these same modeling tools.

MMM↗
DOE OSTI2024

Importance of Electron Correlation on the Geometry and Electronic Structure of [2Fe–2S] Systems: A Benchmark Study of the [Fe 2 S 2 (SCH 3 ) 4 ] 2–,3–,4– , [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S- p -tol) 4 ] 2– , and [Fe 2 S 2 (S- o -xyl) 4 ] 2– Complexes

Iron–sulfur clusters are crucial for biological electron transport and catalysis. Obtaining accurate geometries, energetics, manifolds of their excited electronic states, and reduction energies is important to understand their role in these processes. Here, using a [2Fe–2S] model complex with Fe II and Fe III oxidation states, which leads to different charges, i.e., [Fe 2 S 2 (SMe) 4 ] 2–,3–,4– , we benchmarked a variety of computational methodologies ranging from density functional theory (DFT) to post-Hartree–Fock methods, including complete active space self-consistent field (CASSCF), multireference configuration interaction, the second-order N-electron valence state perturbation theory (NEVPT2), and the linearized integrand approximation of adiabatic connection (AC0) approaches. Additionally, we studied three experimentally well-characterized complexes, [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S-o-tol) 4 ] 2– , and [Fe 2 S 2 (S-o-xyl) 4 ] 2– , via DFT methods. We conclude that the dynamic electron correlation is important for accurately predicting the geometry of these complexes. Broken symmetry (BS) DFT correctly predicts experimental geometries of low-spin multiplicity, while CASSCF does not. However, BS-DFT significantly overestimates the difference between the low- and high-spin electronic states for a given oxidation state. At the same time, CASSCF underestimates it but provides relative energies closer to the reference NEVPT2 results. Finally, AC0 provides energetics of NEVPT2 quality with the additional advantage of being able to use large CASSCF sizes. NEVPT2 gives the best estimates of the Fe III /Fe III → Fe II /Fe III (4.27 eV) and Fe II /F III → Fe II /F II (7.72 eV) reduction energies. The results provide insight into the electronic structure of these complexes and assist in the understanding of their physical properties.

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

Islands in non-minimal dilaton gravity: exploring effective theories for black hole evaporation

We start from (3 + 1)-dimensional Einstein gravity with minimally coupled massless scalar matter, through spherical dimensional reduction, the matter theory is non-minimally coupled with the dilaton in (1 + 1)-dimensions. Despite its simplicity, constructing a self-consistent one-loop effective theory for this model remains a challenge, partially due to a Weyl-invariant ambiguity in the effective action. With a universal splitting property for the one-loop action, the ambiguity can be identified with the state-dependent part of the covariant quantum stress tensor. By introducing on-shell equivalent auxiliary fields to construct minimal candidates of Weyl-invariant terms, we derive a one-parameter family of one-loop actions with unique, regular, and physical stress tensors corresponding to the Boulware, Hartle-Hawking and Unruh states. We further study the back-reacted geometry and the corresponding quantum extremal islands that were inaccessible without a consistent one-loop theory. Along the way, we elaborate on the implications of our construction for the non-minimal dilaton gravity model.

2D Gravity↗
DOE OSTI2024

Strategies to Obtain Reliable Energy Landscapes from Embedded Multireference Correlated Wavefunction Methods for Surface Reactions

Embedded correlated wavefunction (ECW) theory is a powerful tool for studying ground- and excited-state reaction mechanisms and associated energetics in heterogeneous catalysis. Several factors are important to obtaining reliable ECW energies, critically the construction of consistent active spaces (ASs) along reaction pathways when using a multireference correlated wavefunction (CW) method that relies on a subset of orbital spaces in the configuration interaction expansion to account for static electron correlation, e.g., complete AS self-consistent field theory, in addition to the adequate partitioning of the system into a cluster and environment, as well as the choice of a suitable basis set and number of states included in excited-state simulations. Here, in this work, we conducted a series of systematic studies to develop best-practice guidelines for ground- and excited-state ECW theory simulations, utilizing the decomposition of NH 3 on Pd(111) as an example. We determine that ECW theory results are relatively insensitive to cluster size, the aug-cc-pVDZ basis set provides an adequate compromise between computational complexity and accuracy, and that a fixed-clean-surface approximation holds well for the derivation of the embedding potential. Additionally, we demonstrate that a merging approach, which involves generating ASs from the molecular fragments at each configuration, is preferable to a creeping approach, which utilizes ASs from adjacent structures as an initial guess, for the generation of consistent potential energy curves involving open-d-shell metal surfaces, and, finally, we show that it is essential to include bands of excited states in their entirety when simulating excited-state reaction pathways.

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

Use of FLOSIC for understanding anion-solvent interactions

An Achille’s heel of lower-rung density-functional approximations is that the highest-occupied-molecular-orbital energy levels of anions, known to be stable or metastable in nature, are often found to be positive in the worst case or above the lowest-unoccupied-molecular-orbital levels on neighboring complexes that are not expected to accept charge. A trianionic example, [Cr(C2O4)3]3−, is of interest for constraining models linking Cr isotope ratios in rock samples to oxygen levels in Earth’s atmosphere over geological timescales. Here we describe how crowd sourcing can be used to carry out self-consistent Fermi–Löwdin–Orbital-Self-Interaction corrected calculations (FLOSIC) on this trianion in solution. The calculations give a physically correct description of the electronic structure of the trianion and water. In contrast, uncorrected local density approximation (LDA) calculations result in approximately half of the anion charge being transferred to the water bath due to the effects of self-interaction error. Use of group-theory and the intrinsic sparsity of the theory enables calculations roughly 125 times faster than our initial implementation in the large N limit reached here. By integrating charge density densities and Coulomb potentials over regions of space and analyzing core-level shifts of the Cr and O atoms as a function of position and functional, we unambiguously show that FLOSIC, relative to LDA, reverses incorrect solute-solvent charge transfer in the trianion-water complex. In comparison to other functionals investigated herein, including Hartree–Fock and the local density approximation, the FLOSIC Cr 1s eigenvalues provide the best agreement with experimental core ionization energies.

Chemistry↗
DOE OSTI2024

Phonon dynamics in the chiral Kitaev spin liquid

Here, we investigate the effect of a magnetic field on the extended Kitaev spin-liquid state through phonon dynamics. Using a constrained fermionic self-consistent mean field method, we analyze the quantum spin liquid (QSL) ground state for the extended Kitaev model with both the Zeeman term and the perturbative three-spin interaction term 𝜅. Our results demonstrate the dependence of the stability of the Kitaev QSL state on the field direction, consistent with findings in the literature. Additionally, we calculate the phonon dynamics for acoustic phonons coupled to the Majorana fermion excitations of the Kitaev spin-liquid state, discussing the temperature and field evolution of these quantities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗