SEARCH · Engineering Papers
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.
Density functional theory (DFT) calculations in the delta-phase of Pu have long relied on the antiferromagnetic (AFM) spin structure to approximate the material’s unique electronic structure. The resulting tetragonal distortion away from the correct cubic symmetry was believed to be a necessary consequence of allowing an ordered spin structure.
Conditional-probability density functional theory (CP-DFT) is a formally exact method for finding correlation energies from Kohn-Sham DFT without evaluating an explicit energy functional. We present details on how to generate accurate exchange-correlation energies for the ground-state uniform gas. We also use the exchange hole in a CP antiparallel spin calculation to extract the high-density limit. We give a highly accurate analytic solution to the Thomas-Fermi model for this problem, showing its performance relative to Kohn-Sham and may be useful at high temperatures. We explore several approximations to the CP potential. Furthermore, results are compared to accurate parameterizations for both exchange-correlation energies and holes.
The layered material Mn 3 Si 2 Te 6 , with alternating stacking honeycomb and triangular layers, is attracting considerable attention due to its rich physical properties. Here, in this study, using density functional theory and classical Monte Carlo (MC) methods, we systematically study this system with the 3d 5 electronic configuration. Near the Fermi level, the states are mainly contributed by Te 5p orbitals hybridized with Mn 3d orbitals, resembling a charge-transfer system. Furthermore, the spin orientations of the ferrimagnetic (FiM) ground state display different conductive behaviors when along the ab plane or out-of-plane directions: insulating vs metallic states. The energy difference between the FiM [110] insulating and FiM [001] metallic phases is very small (~0.71 meV/Mn) . Changing the angle θ of spin orientation from in-plane to out-of-plane directions, the band gaps of this system are gradually reduced, leading to an insulator-metal transition, resulting in an enhanced electrical conductivity, related to the colossal angular magnetoresistance (MR) effect. Although the three main magnetic couplings were found to be antiferromagnetic, overall the ground state is FiM. In addition, we also constructed the magnetic phase diagram using the classical XY spin model studied with the MC method. Three magnetic phases were obtained, including antiferromagnetic order, noncollinear spin patterns, and FiM order. Moreover, we also investigated the Se and Ge doping into the Mn 3 Si 2 Te 6 system: the FiM state has the lowest energy among the magnetic candidates for both Se- and Ge-doped cases. The magnetic anisotropy energy (MAE) decreases in the Se-doped case because the Mn orbital moment is reduced as the doping x increases. Due to the small spin-orbit-coupling effect of Se, the insulator-metal transition caused by the spin orientation disappears in the Se-doped case, resulting in an insulating phase in the FiM [001] phase. This causes a reduced colossal angular MR. However, both the MAE and the band gap of the Ge-doped case do not change much with increasing doping x . Our results for Mn 3 Si 2 Te 6 could provide guidance to experimentalists and theorists working on this system or related materials.
Here, in this work, we provide a nuanced view of electron correlation in the context of transition metal complexes, reconciling computational characterization via spin and spatial symmetry breaking in single-reference methods with qualitative concepts from ligand-field and molecular orbital theories. These insights provide the tools to reliably diagnose the multi-reference character, and our analysis reveals that while strong (i.e., static) correlation can be found in linear molecules (e.g., diatomics) and weakly bound and antiferromagnetically coupled (monometal-noninnocent ligand or multi-metal) complexes, it is rarely found in the ground-states of mono-transition-metal complexes. This leads to a picture of static correlation that is no more complex for transition metals than it is, e.g., for organic biradicaloids. In contrast, the ability of organometallic species to form more complex interactions, involving both ligand-to-metal σ-donation and metal-to-ligand π-backdonation, places a larger burden on a theory's treatment of dynamic correlation. We hypothesize that chemical bonds in which inter-electron pair correlation is non-negligible cannot be adequately described by theories using MP2 correlation energies and indeed find large errors vs experiment for carbonyl-dissociation energies from double-hybrid density functionals. A theory's description of dynamic correlation (and to a less important extent, delocalization error), which affects relative spin-state energetics and thus spin symmetry breaking, is found to govern the efficacy of its use to diagnose static correlation.
Magic-angle twisted bilayer graphene (MATBG) is notable as a highly tunable platform for investigating strongly correlated phenomena such as unconventional superconductivity and quantum spin liquids, due to easy control of doping level through gating and sensitive dependence of the magic angle on hydrostatic pressure. Experimental observations of correlated insulating states, unconventional superconductivity and ferromagnetism in MATBG indicate that this system exhibits rich exotic phases. In this work, using density functional theory calculations in conjunction with the effective screening medium method, we find the MATBG under pressure at a twisting angle of 2.88 °and simulate how its electronic states evolve when doping level and electric field perpendicular to plane are tuned by gating. Our calculations show that, at doping levels between two electrons and four holes per moiré unit cell, a ferromagnetic (FM) solution with spin density localized at AA stacking sites is lower in energy than the nonmagnetic solution. The magnetic moment of this FM state decreases with both electron and hole doping and vanishes at four electrons/holes doped per moiré unit cell. Hybridization between the flat bands at the Fermi level and the surrounding dispersive bands can take place at finite doping. On increasing the out-of-plane electric field at zero doping, a transition from the FM state to the nonmagnetic one is seen. An investigation of impurity effects shows that both absorption of O 2 molecules and occurrence of Stone–Wales impurities suppress the FM state, and the mechanisms are understood from our calculations. We also analyze the interlayer bonding character due to flat bands via Wannier functions. Finally, we report trivial band topology of the flat bands in the FM state at a certain doping level.
The ability of materials to convert charge current into spin current is fundamental to many spintronics applications. One means of realizing this conversion is via Rashba-Edelstein effect (REE). Using density functional theory simulations, we predict that REE can be induced in the recently discovered family of semiconducting ferroelectrics, CsGe𝑋 3 (𝑋 = I, Br, Cl). The effect is quantified through Rashba-Edelstein coefficients, 𝜒 𝑥𝑦 =−𝜒 𝑦𝑥 , which are nonzero in valence and conduction bands. The largest values, obtained for CsGeI 3 , are 3.45 × 10 10 $\frac{ℏ}{Λ cm}$ and 0.97 × 10 10 $\frac{ℏ}{Λ cm}$ in the conduction and valence bands, respectively. The values are comparable to, and sometimes exceed, those of other inorganic materials, although the maximal values occur away from the band edges. The coefficients' sign couples to the direction of spontaneous polarization, offering opportunities for nonvolatile spin current manipulation via external electric field. Furthermore, these coefficients are highly tunable through strain engineering owing to strain-induced variations in spin textures and energy isosurfaces. Specifically, 𝜒 in the valence band of CsGeI 3 is enhanced to 3.61 × 10 10 $\frac{ℏ}{Λ cm}$ under 5% biaxial strain. Given this potential, we synthesized CsGeI 3 to validate our simulation structure and found excellent agreement between experiment and simulation, thereby allowing extrapolation of our findings to practically significant temperatures. Furthermore, our study identifies promising materials for nonvolatile, multifunctional spintronic applications.
Topotactic reduction of perovskite oxides offers a powerful approach for discovering novel phenomena, such as superconducting infinite-layer nickelates and polar metallicity, and is commonly accompanied by the emergence of multiple valence states and/or complex crystal fields of transition metals. However, understanding the complex interplay between crystal chemistry, electronic structure, and physical properties at the spin- and orbital-resolved levels in these reduced systems remains elusive. Here, in this study, x-ray absorption spectroscopy, resonant inelastic x-ray scattering (RIXS), and density functional theory calculations are used to uncover topotactic metal-insulator transition and orbital-specific crystal field excitations in brownmillerite La 0.67 Ca 0.33 MnO 2.5 thin films. The Mn valence states are found to be Mn 2+ /Mn 3+ , along with their corresponding populations at octahedral and tetrahedral sites, which effectively weaken the Mn-O hybridization compared to the parent perovskite phase. As a result, La 0.67 Ca 0.33 MnO 2.5 films exhibit an antiferromagnetic insulating ground state. Moreover, by combining the RIXS measurements on selected single-valence manganites, specifically MnO, LaMnO 3 , and CaMnO 3 , with orbital- and spin-resolved density-of-states calculations, the study identifies the dd excitations of octahedrally and tetrahedrally coordinated Mn 2+ /Mn 3+ ions, directly linking the microscopic electronic structure to the macroscopic magnetic/electrical properties.
Substrates and inhibitors of Mo-dependent nitrogenase bind and react at Fe ions of the active-site FeMo-cofactor [7Fe–9S–C–Mo–homocitrate] contained within the MoFe protein α-subunit. The cofactor contains a CFe 6 core, a carbon centered within a trigonal prism of six Fe, whose role in catalysis is unknown. Targeted 13 C labeling of the carbon enables electron-nuclear double resonance (ENDOR) spectroscopy to sensitively monitor the electronic properties of the Fe–C bonds and the spin-coupling scheme adopted by the FeMo-cofactor metal ions. Here this report compares 13 CFe 6 ENDOR measurements for (i) the wild-type protein resting state (E 0 ; α-Val 70 ) to those of (ii) α-Ile 70 , (iii) α-Ala 70 -substituted proteins; (iv) crystallographically characterized CO-inhibited “hi-CO” state; (v) E 4 (4H) Janus intermediate, activated for N 2 binding/reduction by accumulation of 4[ e– /H + ]; (vi) E 4 (2H)* state containing a doubly reduced FeMo-cofactor without Fe-bound substrates; and (vii) propargyl alcohol reduction intermediate having allyl alcohol bound as a ferracycle to FeMo-cofactor Fe6. All states examined, both S = 1/2 and 3/2 exhibited near-zero 13 C isotropic hyperfine coupling constants, $^Ca$ = [-1.3 ↔ +2.7] MHz. Density functional theory computations and natural bond orbital analysis of the Fe-C bonds show that this occurs because a (3 spin-up/3 spin-down) spin-exchange configuration of CFe 6 Fe-ion spins produces cancellation of large spin-transfers to carbon in each Fe–C bond. Previous X-ray diffraction and DFT both indicate that trigonal-prismatic geometry around carbon is maintained with high precision in all these states. The persistent structure and Fe–C bonding of the CFe 6 core indicate that it does not provide a functionally dynamic (hemilabile) “beating heart”-instead it acts as “a heart of steel”, stabilizing the structure of the FeMo-cofactor-active site during nitrogenase catalysis.
Models of encapsulated 1/2 nuclear spin H-1 and P-31 atoms in fullerene and diamond nanocrystallite, respectively, are proposed and examined with ab-initio local density functional method for possible applications as single quantum bits (qubits) in solid-state quantum computers. A H-1 atom encapsulated in a fully deuterated fullerene, C(sub 20)D(sub 20), forms the first model system and ab-initio calculation shows that H-1 atom is stable in atomic state at the center of the fullerene with a barrier of about 1 eV to escape. A P-31 atom positioned at the center of a diamond nanocrystallite is the second model system, and 3 1P atom is found to be stable at the substitutional site relative to interstitial sites by 15 eV, Vacancy formation energy is 6 eV in diamond so that substitutional P-31 atom will be stable against diffusion during the formation mechanisms within the nanocrystallite. The coupling between the nuclear spin and weakly bound (valance) donor electron coupling in both systems is found to be suitable for single qubit applications, where as the spatial distributions of (valance) donor electron wave functions are found to be preferentially spread along certain lattice directions facilitating two or more qubit applications. The feasibility of the fabrication pathways for both model solid-state qubit systems within practical quantum computers is discussed with in the context of our proposed solid-state qubits.
A WSe 2 monolayer shows many interesting properties due to its spin-orbit coupling induced spin splitting in bands around the Fermi level and the spin-valley configuration. The orientation of the spin polarization in the relevant bands is crucial for the nature of exciton states and the optical valley selectivity. In this work, we studied the WSe 2 nanoribbons under different mechanical bending curvatures and electron/hole doping with density functional theory and their optical absorption and excitonic states with many-body perturbation GW and BSE (Bethe-Salpeter equation) methods. We found that the WSe 2 nanoribbons can exhibit an enhanced SOC effect and a spatially varying spin polarization in bands around the Fermi level under bending. The spin-polarization can show an anisotropy (or asymmetry) in those nearly degenerate bands, leading to a controllable magnetism via bending and electron/hole doping to the nanoribbons, suggesting a potential application in compact and controllable magnetic nanodevices and spintronics. The optical absorption spectrum of the nanoribbon presents a large tunability with bending within the near infrared region of about 0.4 to 1.5 eV, showing an enhanced absorption at a large bending condition. Finally, the exciton states generally show mixed or various spin configuration in the electron and hole pairs that are controlled by bending, potentially useful for applications in spin-based quantum information processes.
The rotating magnetocaloric effect (RMCE) offers a promising alternative to conventional magnetocaloric configurations by taking advantage of magnetic anisotropy to simplify device architecture and enhance refrigeration efficiency. In this study, the RMCE in a high-quality single crystal of Tb 3 Ni is investigated, a compound previously shown to exhibit significant magnetocaloric behavior along its easy axis of magnetization. By measuring the magnetization and corresponding entropy change along the three principal crystallographic axes using a discontinuous measurement protocol, we verify the anisotropic magnetic properties and derive the RMCE from rotations between hard ( a, b ) and easy ( c ) axes of magnetization. Our results show a maximum value for the rotating entropy change of 19.5 J kg −1 K −1 for µ 0 H = 7 T around 60 K, within the critical temperature window for industrial gas liquefaction applications. Neutron diffraction and magnetic Pair-Distribution Function (mPDF) analysis reveal that the origin of this large anisotropic response lies in a partially ordered incommensurate spin-density wave phase and persistent short-range ferromagnetic (FM) correlations. These complex magnetic states enable the release of a substantial amount of magnetic entropy when the field is applied along the easy c-axis, effectively driving the large RMCE. Comparison with other RMCE materials confirms Tb 3 Ni as one of the most promising candidates in this temperature regime, offering both a large magnetic entropy change and a wide operating window.
The development of quantum gas microscopy for two-dimensional optical lattices has provided an unparalleled tool to study the Fermi-Hubbard model (FHM) with ultracold atoms. Spin-resolved projective measurements, or snapshots, have played a significant role in quantifying correlation functions, theory verification, and thus the uncovering of underlying physical phenomena such as antiferromagnetism at commensurate filling on bipartite lattices and other charge and spin correlations, as well as dynamical properties at various densities. Here we employ a recent concept, the multiscale structural complexity, and show that when computed for the snapshots (of single spin species, local moments, or total density) it can provide a theory-free property, immediately accessible to experiments. Specifically, after benchmarking results for Ising and $XY$ models, we study the structural complexity of snapshots of the repulsive FHM in the two-dimensional square lattice as a function of doping and temperature. We generate projective measurements using determinant quantum Monte Carlo and compare their complexities against those from the experiment. We demonstrate that these complexities are linked to relevant physical observables such as the entropy and double occupancy. Their behaviors capture the development of correlations and relevant length scales in the system. Furthermore, we provide an open-source code in python which can be implemented into data analysis routines in experimental settings for the square lattice.
Abstract Spin defects in wide-bandgap semiconductors provide a promising platform to create qubits for quantum technologies. Their synthesis, however, presents considerable challenges, and the mechanisms responsible for their generation or annihilation are poorly understood. Here, we elucidate spin defect formation processes in a binary crystal for a key qubit candidate—the divacancy complex (VV) in silicon carbide (SiC). Using atomistic models, enhanced sampling simulations, and density functional theory calculations, we find that VV formation is a thermally activated process that competes with the conversion of silicon (V Si ) to carbon monovacancies (V C ), and that VV reorientation can occur without dissociation. We also find that increasing the concentration of V Si relative to V C favors the formation of divacancies. Moreover, we identify pathways to create spin defects consisting of antisite-double vacancy complexes and determine their electronic properties. The detailed view of the mechanisms that underpin the formation and dynamics of spin defects presented here may facilitate the realization of qubits in an industrially relevant material.
We demonstrate for the case of photoexcited [Ru(2,2'-bipyridine) 3 ] 2+ how femtosecond resonant inelastic X-ray scattering (RIXS) at the ligand K-edge allows one to uniquely probe changes in the valence electronic structure following a metal-to-ligand charge-transfer (MLCT) excitation. Metal–ligand hybridization is probed by nitrogen-1s resonances providing information on both the electron-accepting ligand in the MLCT state and the hole density of the metal center. By comparing to spectrum calculations based on density functional theory, we are able to distinguish the electronic structure of the electron-accepting ligand and the other ligands and determine a temporal upper limit of (250 ± 40) fs for electron localization following the charge-transfer excitation. The spin of the localized electron is deduced from the selection rules of the RIXS process establishing new experimental capabilities for probing transient charge and spin densities.
LuFe 4 Ge 2 crystallizes in the ZrFe 4 Si 2 -type structure, hosting chains of Fe tetrahedra giving rise to geometric frustration and low dimensionality. The compound orders antiferromagnetically at around 36 K accompanied by a simultaneous structural transition from a tetragonal phase to an orthorhombic phase. The hydrostatic pressure dependence of the magnetic and structural transitions is investigated using electrical transport, AC magnetic susceptibility, AC calorimetry, Mössbauer, muon-spin relaxation (μSR), and x-ray-diffraction measurements. External pressure suppresses the first-order transition to the antiferromagnetic phase (AFM1) around 1.8 GPa. The structural transition is largely unaffected by pressure and remains between 30 to 35 K for pressures up to 2 GPa. A second antiferromagnetic phase (AFM2) is observed at higher pressures. The transition from the paramagnetic to the AFM2 phase is of second-order nature and appears to be connected to the structural transition. The magnetic volume fraction obtained from μSR and Mössbauer measurements reveal that the entire sample undergoes magnetic ordering in both magnetic phases. In addition, similar low-temperature muon-precession frequencies in AFM1 and AFM2 phases point at similar ordered moments and magnetic structures in both phases. Our results further indicate enhanced magnetic fluctuations in the pressure-induced AFM2 phase. The experimental observations together with density functional theory calculations suggest that the magnetic- and structural-order parameters in LuFe 4 Ge 2 are linked by magnetic frustration, causing the simultaneous magnetostructural transition.
Ternary intermetallic compound \({\text {Pr}}_2 {\text {Co}}_{0.86} {\text {Si}}_{2.88}\) has been synthesized in single phase and characterized by x-ray diffraction, scanning electron microscopy with energy dispersive x-ray spectroscopy (SEM-EDX) analysis, magnetization, heat capacity, neutron diffraction and muon spin rotation/relaxation ( \(\mu\) SR) measurements. The polycrystalline compound was synthesized in single phase by introducing necessary vacancies in Co/Si sites. Magnetic, heat capacity, and zero-field neutron diffraction studies reveal that the system undergoes magnetic transition below \(\sim\) 4 K. Neutron diffraction measurement further reveals that the magnetic ordering is antiferromagnetic in nature with an weak ordered moment. The high temperature magnetic phase has been attributed to glassy in nature consisting of ferromagnetic clusters of itinerant (3 d ) Co moments as evident by the development of internal field in zero-field \(\mu\) SR below 50 K. The density-functional theory (DFT) calculations suggest that the low temperature magnetic transition is associated with antiferromagnetic coupling between Pr 4 f and Co 3 d spins. Pr moments show spin fluctuation along with unconventional orbital moment quenching due to crystal field. The evolution of the symmetry and the crystalline electric field environment of Pr-ions are also studied and compared theoretically between the elemental Pr and when it is coupled with other elements such as Co. The localized moment of Pr 4 f and itinerant moment of Co 3 d compete with each other below \(\sim\) 20 K resulting in an unusual temperature dependence of magnetic coercivity in the system.
The isotropic harmonic oscillator supplemented by a strong spin-orbit interaction has been the cornerstone of nuclear structure since its inception more than seven decades ago. In this paper we introduce—or rather re-introduce—the “Dirac oscillator,” a fully relativistic basis that has all the desired attributes of the ordinary harmonic oscillator while naturally incorporating a strong spin-orbit coupling. To assess the power and flexibility of the Dirac oscillator basis in the solution of nuclear structure problems within the framework of covariant density-functional theory. Here, self-consistent calculations of binding energies and ground-state densities for a selected set of doubly magic nuclei are performed using the Dirac oscillator basis and are then compared against results obtained with the often-used Runge-Kutta method. Results obtained using the Dirac oscillator basis reproduce with high accuracy those derived using the Runge-Kutta method and suggest a clear path for a generalization to systems with axial symmetry. Although the harmonic oscillator with spin-orbit corrections has been the staple of the nuclear shell model since the beginning, the Dirac oscillator is practically unknown among the nuclear physics community. In this paper we illustrate the power and flexibility of the Dirac oscillator and suggest extensions to the study of systems without spherical symmetry, as required in constrained calculations of nuclear excitations.