SEARCH · Engineering Papers

Results for “Orbit determination”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2024

A static quantum embedding scheme based on coupled cluster theory

Here, we develop a static quantum embedding scheme that utilizes different levels of approximations to coupled cluster (CC) theory for an active fragment region and its environment. To reduce the computational cost, we solve the local fragment problem using a high-level CC method and address the environment problem with a lower-level Møller–Plesset (MP) perturbative method. This embedding approach inherits many conceptual developments from the hybrid second-order Møller–Plesset (MP2) and CC works by Nooijen [J. Chem. Phys. 111, 10815 (1999)] and Bochevarov and Sherrill [J. Chem. Phys. 122, 234110 (2005)]. We go beyond those works here by primarily targeting a specific localized fragment of a molecule and also introducing an alternative mechanism to relax the environment within this framework. We will call this approach MP-CC. We demonstrate the effectiveness of MP-CC on several potential energy curves and a set of thermochemical reaction energies, using CC with singles and doubles as the fragment solver, and MP2-like treatments of the environment. The results are substantially improved by the inclusion of orbital relaxation in the environment. Using localized bonds as the active fragment, we also report results for N=N bond breaking in azomethane and for the central C–C bond torsion in butadiene. We find that when the fragment Hilbert space size remains fixed (e.g., when determined by an intrinsic atomic orbital approach), the method achieves comparable accuracy with both a small and a large basis set. Additionally, our results indicate that increasing the fragment Hilbert space size systematically enhances the accuracy of observables, approaching the precision of the full CC solver.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Vibrational Spin-Orbit Coupling Contributions to Excited State Decay of Ligand-to-Ligand Charge Transfer States

Controlling excited state relaxation processes is important in a variety of photochemical and photophysical processes, including the generation of ground and excited state spin polarization for quantum information science applications. Here, we address how specific static distortions – based on vibrational spin-orbit active modes at C2v symmetry determined by group theory – in a series of low-symmetry ligand-to-ligand charge transfer complexes enable direct spin-orbit coupling contributions to T1 → S0 excited state decay. These results are used to address spin-vibronic coupling contributions to T1 → S0 decay in a high-symmetry (tBu2bpy)Pt(S,S) (tBu2bpy = 4,4’-di-tert-butyl-2,2’-bipyridine and S,S = benzene-1,2-dithiolate) ligand-to-ligand charge transfer complex with effective C2v symmetry, where T1 relaxation is both spin- and orbitally forbidden due to the direct spin-orbit coupling matrix element being zero by symmetry. Low-frequency vibrations that involve a pyridine-pyridine twisting motion within the bpy ligand generate large ∂/∂Qi values that will contribute significantly to T1 → S0 relaxation. The work advances ligand design strategies for the generation of tailored T1 → S0 relaxation rates, which can be utilized to optimize the generation of electron spin polarization in radical-elaborated ligand-to-ligand charge transfer complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Single-neutron adding on $^{34}$S

For this work, single-neutron adding data was collected in order to determine the distribution of the single-neutron strength of the 0ƒ 7/2 , 1p 3/2 , 1p 1/2 and 0ƒ 5/2 orbitals outside of Z = 16, N = 18, 34 S. The 34 S(d,p) 35 S reaction has been measured at 8 MeV/u to investigate cross sections to excited states in 35 S. Outgoing proton yields and momenta were analyzed by the Super-Enge Split-Pole Spectrograph in conjunction with the CeBrA demonstrator located at the John D. Fox Laboratory at Florida State University. Angular distributions were compared with Distorted Wave Born Approximation calculations in order to extract single-neutron spectroscopic overlaps. Spectroscopic overlaps and strengths were determined for states in S up through 6 MeV in excitation energy. Each orbital was observed to have fragmented strength where a single level carried the majority. The single-neutron centroids of the 0ƒ 7/2 , 1p 3/2 , 1p 1/2 and 0ƒ 5/2 orbitals were determined to be 2360$^{+90}_{-40}$ keV, 3280$^{+80}_{-50}$ keV, 4780$^{+60}_{-40}$ keV, and ≳7500 keV, respectively. A previous discrepancy in the literature with respect to the distribution of the neutron 1p 1/2 strength was resolved. The integration of the normalized spectroscopic strengths, up to 5.1 MeV in excitation energy, revealed fully-vacant occupancies for the 0ƒ 7/2 , 1p 3/2 , and 1p 1/2 orbitals, as expected. The spacing in the single-neutron energies highlighted a reduction in the traditional N = 28 shell-gap, relative to both the 1p spin-orbit energy difference (N = 32) and the lower limit on the N = 34 shell spacing.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

The Juno mission as a probe of long-range new physics

Orbits of celestial objects, especially the geocentric and heliocentric ones, have been well explored to constrain new long-range forces beyond the Standard Model (SM), often referred to as fifth forces. In this paper, for the first time, we apply the motion of a spacecraft around Jupiter to probe fifth forces that don’t violate the equivalence principle. The spacecraft is the Juno orbiter, and ten of its early orbits already allow a precise determination of the Jovian gravitational field. We use the shift in the precession angle as a proxy to test non-gravitational interactions between Juno and Jupiter. Requiring that the contribution from the fifth force does not exceed the uncertainty of the precession shift inferred from data, we find that a new parameter space with the mass of the fifth-force mediator around 10 −14 eV is excluded at 95% C.L.

new light particles↗

DOE OSTI2024

ELECTRONIC STRUCTURE METHODS AND PROTOCOLS WITH APPLICATION TO DYNAMICS, KINETICS AND THERMOCHEMISTRY

Hydrocarbon combustion involves the reaction dynamics of a tremendous number of species beginning with many-component fuel mixtures and proceeding via a complex system of intermediates to form primary and secondary products. Combustion conditions corresponding to new advanced engines and/or alternative fuels rely increasingly on autoignition and low-temperature-combustion chemistry. In these regimes various transient radical species such as HO2, ROO·, ·QOOH, HCO, NO2, HOCO, and Criegee intermediates play important roles in determining the detailed as well as more general dynamics. A clear understanding and accurate representation of these processes is needed for effective modeling. Given the difficulties associated with making reliable experimental measurements of these systems, computation can play an important role in developing these energy technologies. Accurate calculations have their own challenges since even within the simplest dynamical approximations such as transition state theory, the rates depend exponentially on critical barrier heights and these may be sensitive to the level of quantum chemistry. Moreover, it is well-known that in many cases it is necessary to go beyond statistical theories and consider the dynamics. Quantum tunneling, resonances, radiative transitions, and non-adiabatic effects governed by spin-orbit or derivative coupling can be determining factors in those dynamics. Building upon progress made during a period of prior support through the DOE Early Career Program, this project combines developments in the areas of potential energy surface (PES) fitting and multistate multireference quantum chemistry to allow spectroscopically and dynamically/kinetically accurate investigations of key molecular systems (such as those mentioned above), many of which are radicals with strong multireference character and have the possibility of multiple electronic states contributing to the observed dynamics. An ongoing area of investigation is to develop general strategies for robustly convergent electronic structure theory for global multichannel reactive surfaces including diabatization of energy and other relevant surfaces such as dipole transition. Combining advances in ab initio methods with automated interpolative PES fitting allows the construction of high-quality PESs (incorporating thousands of high-level data) to be done rapidly through parallel processing on high-performance computing (HPC) clusters. In addition, new methods and approaches to electronic structure theory will be developed and tested through applications. This project will explore limitations in traditional multireference calculations (e.g., MRCI) such as those imposed by internal contraction, lack of high-order correlation treatment and poor scaling. Methods such as DMRG-based extended active-space CASSCF and various Quantum Monte Carlo (QMC) methods will be applied (including VMC/DMC and FCIQMC). Insight into the relative significance of different orbital spaces and the robustness of application of these approaches on leadership class computing architectures will be gained. Synergy with other components of this research program such as automated PES fitting and multireference quantum chemistry will be used to address challenges encountered by the standard approaches to computational thermochemistry (those being single-reference quantum chemistry and perturbative treatments of the anharmonic vibrational energy, which break down for some cases of electronic structure or floppy strongly coupled vibrational modes).

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

Anion photoelectron spectroscopy and chemical bonding of ThS 2 – and ThSO –

Anion photoelectron spectra of ThSO – and ThS 2 – were recorded using the third (355 nm) harmonic of an Nd-YAG laser; these provided the measured vertical detachment energies of each anion. The experiments are supported by extensive coupled cluster calculations on ThSO, ThSO – , ThS 2 , and ThS 2 – , as well as the oxygen congeners ThO 2 and ThO 2 – . The ab initio calculations, which included complete basis set extrapolations, spin–orbit effects using four-component coupled cluster, and higher-order correlation contributions through CCSDT(Q), yielded an adiabatic electron affinity for ThO 2 – that was within 0.02 eV of the previously determined experimental value. The singly occupied molecular orbital (SOMO) in all three anions corresponds primarily to the 7s orbital on Th. Successive substitution of S for each O in ThO 2 leads to larger electron affinities and smaller bond angles in the neutral molecules, but larger angles in the anions. As demonstrated by Franck–Condon simulations of the spectra using the CCSD(T) spectroscopic constants, substitution of O by S significantly complicates the resulting detachment spectra due to the lower vibrational frequencies in the sulfur species. Altogether the calculated vertical detachment energies are in very good agreement with the experiment. In addition to the adiabatic electron affinities of each species, atomization energies and heats of formation have also been determined via the FPD approach with expected uncertainties of 1–2 kcal/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Nonadiabatic dynamics of photoexcited thiopyridone isomers: An interplay between El-Sayed’s conditions and energy gap law

We investigated the non-adiabatic dynamics of photoexcited thiopyridone systems across their ortho-, meta-, and para-isomeric forms. The relaxation pathways of the three isomers in both gas phase and solvent environments are mapped using surface hopping dynamics based on time-dependent density functional theory. Our analysis highlights the influence of isomeric structures on photophysical behavior, offering insights into design principles to control photochemical phenomena. The simulations suggest a systematic reduction in the rate of intersystem crossing (ISC) from ortho- to meta- to para-isomer. Comparisons with multiconfigurational wave function methods in the gas phase further demonstrate how electronic structure influences the predicted dynamical pathways. The simulated dynamics demonstrates that the spin–orbit coupling strength alone does not determine the rate of ISC, as both state energetics and underlying electronic and structural features play decisive roles. These aspects explain the much slower ISC in the para-isomer, as well as the non-negligible role of the El-Sayed forbidden pathway in the computed ISC dynamics.

Complete-active space self-consistent field↗

DOE OSTI2024

Theoretical Studies of Polar Systems near Ferroelectrics Quantum Critical Points (Final Report)

The PI has studied emergent quantum phases near polar quantum critical points. A first challenge was to develop a theoretical description of the observation that many polar materials undergo classical first-order transitions while displaying quantum criticality. The possibility of novel metallic states near polar quantum critical points was then studied; identification and characterization of a non-Fermi liquid phase was made in a multiband system. In dilute quantum critical polar metals, electronic coupling to polar energy fluctuations was shown to result in attractive electron-electron interactions and to superconductivity. A spin-phonon resonance measurement in applied magnetic field was proposed to determine the magnitude of spin-orbit mediated electron-phonon coupling; furthermore in the polar phase new field-induced phonon collective modes were identified with specific signatures for experiment. In order to keep current, the PI has also worked on other types of strongly correlated problems. With experimental groups, she has developed phenomenologies to identify and explain observed dynamical critical behaviors. The PI has studied the interplay of topology, fractionalization and deconfinement in a minimalist spin model. An emergent phase transition in a 1 + 1 frustrated spin nanotube has also been identified and characterized. The PI has also begun exploring criticality out of equilibrium, characterizing rich dynamical phases associated with photoinduced polar transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Low-lying states and total internal partition sums of CH

The electronic structure and spin-orbit states of the CH radical have been systematically investigated using multi-reference configuration interaction (MRCI) and single-reference coupled-cluster (CC) methods. These calculations were performed in conjunction with large correlation-consistent basis sets of quadruple-, quintuple-, and sextuple-ζ quality. To achieve high accuracy, electronic energies for all states were extrapolated to the complete basis set (CBS) limit, enabling the detailed construction of potential energy curves and determination of reliable spectroscopic constants. Spin-orbit coupling effects were explicitly incorporated, and vibrational energy levels were computed via Numerov analysis. Furthermore, the resulting values exhibit good to excellent agreement with available experimental data. Dipole moment and transition dipole moment curves were evaluated to assess the opacity characteristics of CH, revealing that transitions such as Χ 2 Π (u′′ = 0) → Α 2 Δ (u′ = 0), Χ 2 Π (u′′ = 0) → Β 2 Σ − (u′ = 0), Χ 2 Π (u′′ = 0) → C 2 Σ + (u′ = 0), and Χ 2 Π (u′′ = 0) → D 2 Σ + (u′ = 3) are particularly probable. Finally, the total internal partition function sum (TIPS) of CH was computed over a broad temperature range (10–30,000 K) based on our high-accuracy ab initio results.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

Birefringence tests of gravity with multimessenger binaries

Extensions to General Relativity (GR) allow the polarization of gravitational waves (GW) from astrophysical sources to suffer from amplitude and velocity birefringence, which respectively induce changes in the ellipticity and orientation of the polarization tensor. Here, we introduce a multi-messenger approach to test this polarization behavior of GWs during their cosmological propagation using binary sources, for which the initial polarization is determined by the inclination and orientation angles of the orbital angular momentum vector with respect to the line of sight. In particular, we use spatially resolved radio imaging of the jet from a binary neutron star (BNS) merger to constrain the orientation angle and hence the emitted polarization orientation of the GW signal at the site of the merger, and compare to that observed on Earth by GW detectors. For GW170817, using past measurements of the inclination angle, we constrain the deviation from GR due to amplitude birefringence to κ A = $-0.12^{+0.60}_ {-0.61}$, while the velocity birefringence parameter κ V remains unconstrained. The inability to constrain κ V is due to the low amplitude of GW170817 in the Virgo detector, and measurements of the polarization orientation require information from a combination of multiple detectors with different alignments. For this reason, we also mock future BNS mergers with resolved afterglow proper motion and project that κ V could be constrained to a precision of 5 rad (corresponding to an angular shift of the GW polarization of δ φV ≈ 0.2 rad for a BNS at 100 Mpc) by a future network of third-generation ground-based GW detectors such as Cosmic Explorer and the radio High Sensitivity Array. Crucially, this velocity birefringence effect cannot be constrained with dark binary mergers as it requires polarization information at the emission time, which can be provided only by electromagnetic emission.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Influence of temperature, doping, and amorphization on the electronic structure and magnetic damping of iron

Hybrid magnonic quantum systems have drawn increased attention in recent years for coherent quantum information processing, but too large magnetic damping is a persistent concern when metallic magnets are used. Their intrinsic damping is largely determined by electron-magnon scattering induced by spin-orbit interactions. In the low scattering limit, damping is dominated by intra-band electronic transitions, which has been theoretically shown to be proportional to the electronic density of states at the Fermi level. In this work, we focus on body-centered-cubic iron as a paradigmatic ferromagnetic material. We comprehensively study its electronic structure using first-principles density functional theory simulations and account for finite lattice temperature, boron (B) doping, and structure amorphization. Our results indicate that temperature induced atomic disorder and amorphous atomic geometries only have a minor influence. Instead, boron doping noticeably decreases the density of states near the Fermi level with an optimal doping level of 6.25%. In addition, we show that this reduction varies significantly for different atomic geometries and report that the highest reduction correlates with a large magnetization of the material. Furthermore, this may suggest materials growth under external magnetic fields as a route to explore in experiment.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Energetic Electron Transport in Magnetized Plasma with Magnetic Islands (Final Technical Report)

The main scientific goal of this project was to investigate the interactions between energetic electrons and magnetic islands in magnetized plasmas from lab to space. The proposed plan of addressed the following specific objectives: (1) Determine if energetic electrons are trapped by magnetic islands. Specifically, test the hypothesis that energetic electrons are trapped near island O-points due to stable orbits and de-confined near island X-points due to X-point tangles. (2) Determine how electron transport changes with island width and location. Specifically, test the hypothesis that electrons are accelerated by contracting magnetic islands through a Fermi acceleration process. (3) Assess how electron transport is affected by island dynamics, including island rotation, bifurcation, and/or island overlap. Specifically, test the hypothesis that electrons are deconfined and possibly accelerated during island rotations, bifurcations, and overlap due to stochastization of the magnetic field lines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Time-domain study of coupled collective excitations in quantum materials

Quantum materials hold immense promises for future applications due to their intriguing electronic, magnetic, thermal, and mechanical properties that often arise from a complex interplay between microscopic degrees of freedom. Important insights of such interactions come from studying the collective excitations of electrons, spins, orbitals, and lattice, whose cooperative motions play a crucial role in determining the novel behavior of these systems and offer us a key tuning knob to modify material properties on-demand through external perturbations. In this regard, ultrafast light-matter interaction has shown great potential in controlling the couplings of collective excitations, and rapid progress in a plethora of time-resolved techniques down to the attosecond regime has significantly advanced our understanding of the coupling mechanisms and guided us in manipulating the dynamical properties of quantum materials. This review aims to highlight recent experiments on visualizing collective excitations in the time domain, focusing on the coupling mechanisms between different collective modes such as phonon-phonon, phonon-magnon, phonon-exciton, magnon-magnon, magnon-exciton, and various polaritons. We introduce how these collective modes are excited by an ultrashort laser pulse and probed by different ultrafast techniques, and we explain how the coupling between collective excitations governs the ensuing nonequilibrium dynamics. We also provide some perspectives on future studies that can lead to discoveries of the emergent properties of quantum materials both in and out of equilibrium.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Multireference diffusion Monte Carlo reaches 2D materials

Abstract Quantum confinement in 2D materials strongly enhances electronic correlation effects. Therefore, predicting the properties of these unique materials, with both a high level of accuracy and computational efficiency, without relying on adjustable parameters or functionals, remains an outstanding theoretical challenge. The majority of theoretical studies are based on the approximations of density functional theory (DFT). The reliability of DFT predictions are heavily dependent on the choice of an approximated exchange-correlation functional. Here, we estimate the magnitude of impact of correlation on the total energy for the quintessential 2D material, graphene, by performing and comparing state-of-the-art selected CI and quantum Monte Carlo extrapolated calculations for a single unit cell at the$$\Gamma$$point. We demonstrate that Self-Healing Diffusion Monte Carlo (SHDMC) obtains a very compact, but high-quality wavefunction for this system that lacks the strong basis set dependence displayed by state of the art quantum chemistry methods. The SHDMC wavefunction is of higher quality compared to that obtained from sCI, in the same orbital basis, while being$$\sim$$ 1000 times smaller in terms of determinant count compared to sCI. We also demonstrate that extrapolating SHDMC results to the infinite determinant limit compares extremely well with complete basis set extrapolated sCI. Our work paves the way for future validation of SHDMC applied to challenging 2D materials.

Science & Technology - Other Topics↗

DOE OSTI2024

The Santa Barbara Binary-disk Code Comparison

We have performed numerical calculations of a binary interacting with a gas disk, using 11 different numerical methods and a standard binary-disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed state of the binary-disk system. We find that all codes can agree on a converged solution (depending on the diagnostic being measured). The zone spacing required for most codes to reach a converged measurement of the torques applied to the binary by the disk is roughly 1% of the binary separation in the vicinity of the binary components. For our disk model to reach a relaxed state, codes must be run for at least 200 binary orbits, corresponding to about a viscous time for our parameters, 0.2(a 2 Ω B /v) is the kinematic viscosity. The largest discrepancies between codes resulted from the dimensionality of the setup (3D vs. 2D disks). We find good agreement in the total torque on the binary between codes, although the partition of this torque between the gravitational torque, orbital accretion torque, and spin accretion torque depends sensitively on the sink prescriptions employed. In agreement with previous studies, we find a modest difference in torques and accretion variability between 2D and 3D disk models. We find cavity precession rates to be appreciably faster in 3D than in 2D.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Tuning the magnetic properties of the spin-split antiferromagnet MnTe through pressure

The hexagonal antiferromagnet MnTe has attracted enormous interest as a prototypical example of a spin-compensated magnet in which the combination of crystal and spin symmetries lifts the spin degeneracy of the electron bands without the need for spin-orbit coupling, a phenomenon called nonrelativistic spin splitting (NRSS). Subgroups of NRSS are determined by the specific spin-interconverting symmetry that connects the two opposite-spin sublattices. In MnTe, this symmetry is rotation, leading to the subgroup with spin splitting away from the Brillouin zone center, often called altermagnetism. MnTe also has the largest spontaneous magnetovolume effect of any known antiferromagnet, implying strong coupling between the magnetic moment and volume. This magnetostructural coupling offers a potential knob for tuning the spin-splitting properties of MnTe. Here, we use neutron diffraction with in situ applied pressure to determine the effects of pressure on the magnetic properties of MnTe and further explore this magnetostructural coupling. We find that applying pressure significantly increases the Néel temperature, but decreases the ordered magnetic moment. We explain this as a consequence of strengthened magnetic exchange interactions under pressure, resulting in higher 𝑇 N , with a simultaneous reduction of the local moment of individual Mn atoms, described here via density functional theory. This reflects the increased orbital hybridization and electron delocalization with pressure. In conclusion, these results shed light on the competition between magnetic exchange interactions and the strength of individual magnetic moments and show that the magnetic properties of MnTe can be controlled by pressure, opening the door to improved properties for spintronic applications through tuning via physical or chemical pressure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

The Curious Case of Pu + : Insight on 5f Orbital Activity from Inductively Coupled Plasma Tandem Mass Spectrometry (ICP-MS/MS) Reactions

A comprehensive understanding of when and how 5f orbitals participate in complex chemical bonding is important for a variety of applications. The actinides are unique in that they possess 5f orbitals and can access high oxidation states, which make them attractive for use in catalysis. Fundamental studies of actinide–ligand interactions offer a mechanism to examine the activation of the 5f orbitals so that the selectivity of 5f orbitals can be assessed. A previous study examined the reaction of Pu + + CO 2 and determined that the reaction efficiency is restricted by a barrier, namely, promotion from the Pu + ground-state configuration, 5f 6 7s, to a reactive-state configuration, 5f 5 6d 2 . Here, the present study illustrates the benefit of activation of Pu’s 5f orbitals when studying the reaction of Pu + + NO. In this reaction, PuO + forms in an exothermic, barrierless process. The 5f orbitals can and do participate in forming a linear intermediate, [N-Pu-O] + , and this drives the exothermic reaction. Understanding the conditions under which 5f orbitals are active in chemical bonding is the key to exploiting the actinides’ selective catalytic capabilities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Microscopic origin of temperature-dependent magnetism in spin-orbit-coupled transition metal compounds

A few 4 d and 5 d transition metal compounds with various electron fillings were recently found to exhibit magnetic susceptibilities χ and magnetic moments that deviate from the well-established Kotani model. This model has been considered for decades to be the canonical expression for descriing the temperature dependence of magnetism in systems with nonnegligible spin-orbit coupling effects. In this paper, we uncover the origin of such discrepancies and determine the applicability and limitations of the Kotani model by calculating the temperature dependence of the magnetic moments of a series of 4 d (Ru-based) and 5 d (W-based) systems at different electron fillings. For this purpose, we perform exact diagonalization of -derived relativistic multiorbital Hubbard models on finite clusters and compute their magnetic susceptibilities. Comparison with experimentally measured magnetic properties indicates that contributions such as a temperature-independent χ 0 background, crystal field effects, Coulomb and Hund's couplings, and intersite interactions—not included in the Kotani model—are especially crucial for correctly describing the temperature dependence of χ and magnetic moments at various electron fillings in these systems. Based on our results, we propose a generalized approach beyond the Kotani model to accurately describe their magnetism. Published by the American Physical Society 2025

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗