Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles calculation”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y alloys

This dataset contains input and output files from density functional theory (DFT) simulations used to study the disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y (0 ≤ x + y ≤ 1) alloys and their binary end members MnTe, MnSb, and MnBi. The alloys adopt the hexagonal NiAs-type (nickeline) structure and span ternary (MnTexSb1-x, MnTexBi1-x, MnBixSb1-x), and quaternary compositions across the full MnTe–MnSb–MnBi composition triangle. For each alloy composition, the dataset provides DFT calculations in three magnetic configurations: A-type antiferromagnetic (AFM), C-type AFM, and ferromagnetic (FM). Every magnetic configuration folder contains the fully relaxed crystal structure (CONTCAR), VASP input parameters (INCAR), and the main VASP output file (OUTCAR), from which total electronic energies, Mn magnetic moments, lattice parameters, and percent volume changes between magnetic states are extracted. These data are used to construct compositional phase diagrams, evaluate thermodynamic stability (formability), and map magnetoelastic responses across the alloy space. For A-type AFM and FM configurations, additional data are provided as follows: (i) FORCE_CONSTANTS and thermal_properties.yaml files at the top level of A-type_AFM/ and FM/ folders — present only for compositions marked with an asterisk (*) in Table I of the main text. These are derived from Phonopy finite-displacement calculations on full disordered 128-atom supercells and provide vibrational free energy, entropy (Svib)contribution from explicit disorder calculations. (Table I of the associated main manuscript) (ii) A VCA/ subfolder within A-type_AFM/ and FM/, containing FORCE_CONSTANTS and thermal_properties.yaml from Virtual Crystal Approximation phonon calculations (without spin-orbit coupling). VCA data are available for all compositions and are used to estimate vibrational contributions to the Gibbs free energy across the full composition space. (iii) A SOC/ subfolder containing CONTCAR, INCAR, and OUTCAR from spin-orbit coupling calculations, providing relativistic corrections to electronic energies and lattice parameters (Tables S2–S3 of the SM, and Table I of the main manuscript). (iv) A SOC/VCA/ subfolder containing FORCE_CONSTANTS and thermal_properties.yaml from VCA phonon calculations performed within the SOC framework, combining relativistic and vibrational thermodynamic corrections. The computed properties are used to map the AFM–FM magnetic crossover near MnTe0.75Sb0.25, demonstrate disorder- and spin-induced phonon broadening, identify a semiconductor-to-metal crossover, and quantify the pronounced magnetoelastic volume response near the magnetic phase boundary.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electron-phonon coupling strength from ab initio frozen-phonon approach

In this work, we propose a fast method for high-throughput screening of potential superconducting materials. The method is based on calculating metallic screening of zone-center phonon modes, which provides an accurate estimate for the electron-phonon coupling strength. This method is complementary to the recently proposed Rigid Muffin Tin (RMT) method, which amounts to integrating the electronphonon coupling over the entire Brillouin zone (as opposed to the zone center), but in a relatively inferior approximation. We illustrate the use of this method by applying it to MgB 2 , where the high-temperature superconductivity is known to be driven largely by the zone-center modes, and compare it to a sister compound AlB 2 . We further illustrate the usage of this descriptor by screening a large number of binary hydrides, for which accurate first-principle calculations of electron-phonon coupling have been recently published. Together with the RMT descriptor, this method opens a way to perform initial high-throughput screening in search of conventional superconductors via machine learning or data mining.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Diffusion energy barrier of Au on Bi 2 Se 3 : theory and experiment

The stability and diffusion of ultra-thin thermally deposited Au films on Bi 2 Se 3 was studied using scanning tunneling microscopy and density functional theory calculations. The Au/Bi 2 Se 3 interface is of interest as gold is predicted to provide excellent electrical contact while maintaining the spin-polarized characteristics of the electronic states in Bi 2 Se 3 that make the material attractive for spintronic applications. When deposited at room temperature, Au 10 covers the surface with tightly packed islands of nanometer scale dimension. The surface morphology is stable up to 400K. At this annealing temperature, Au atoms have sufficient energy to diffuse across the surface and aggregate into larger nanostructures. At 550K, the Bi 2 Se 3 surface is only sparsely covered, and the Au has formed clusters with length scales 5-10 times larger than the original islands formed at room temperature. Comparison of the experiment and first principle calculation lead to the conclusion that the diffusion energy barrier for Au on Bi 2 Se 3 is as high as 0.47 eV, 15 which is much larger than diffusion barriers on other van der Waals materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Modeling Short-Range Order in Disordered Rocksalt Cathodes by Pair Distribution Function Analysis

Pair distribution function (PDF) analysis is a powerful technique for the characterization of short-range order (SRO) in disordered materials. Accurate interpretation of experimental PDF data is critically reliant on the development of structural models that can account for local variations in site occupancies and bond lengths. To this end, we outline an approach to model SRO using first-principles calculations based on the cluster-expansion formalism. These methods are validated on neutron scattering data from two disordered rocksalt oxyfluorides, Li 1.3 Mn 0.4 Ti 0.3 O 1.7 F 0.3 and Li 1.3 Mn 0.4 Nb 0.2 Ti 0.1 O 1.7 F 0.3 . For each composition, we demonstrate that an average structure without any SRO fails to reproduce several key features in the experimental PDF. To pinpoint the origin of the suspected SRO in these materials, configurational and displacive effects were separately investigated using two disparate models. Special quasi-random structures were relaxed using density functional theory to account for local changes in bond lengths while maintaining a near-random ionic configuration. This leads to slightly improved accuracy but still misrepresents asymmetry in the first few peaks of the PDF. Monte Carlo simulations were performed to model configurational SRO on a fixed lattice, which by itself is shown to have a minimal influence on the PDF. Instead, we find that it is the bond length relaxations within environments created by SRO which controls the details of the PDF, thereby highlighting the subtle but important coupling between configurational and displacive SRO in disordered materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Perpendicular magnetic anisotropy in half-metallic thin-film Co 2 CrAl

Magnetocrystalline anisotropy (MCA) is one of the key parameters investigated in spin-based electronics (spintronics), e.g. for memory applications. Here, we employ first-principles calculations to study MCA in thin film full Heusler alloy Co 2 CrAl. This material was studied in the past, and has been reported to exhibit half-metallic electronic structure in bulk geometry. In our recent work, we showed that it retains a 100% spin-polarization in thin-film geometry, at CrAl atomic surface termination. Here, we show that the same termination results in a perpendicular magnetic anisotropy, while Co surface termination not only destroys the half-metallicity, but also results in in-plane magnetization orientation. In addition, for films thicker than around 20 nm the contribution from magnetic shape anisotropy may become decisive, resulting in in-plane magnetization orientation. To the best of our knowledge, this is one of the first reports of half-metallic thin-film surfaces with perpendicular magnetic anisotropy. Furthermore, this result may be of interest for potential nano-device applications, and may stimulate a further experimental study of this and similar materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

NEXAFS Spectroscopy of P3HT and PBTTT at the Sulfur K-Edge

The sulfur K-edge near-edge X-ray absorption fine-structure (NEXAFS) spectra of the common conjugated polymers P3HT and PBTTT are studied from both experimental and theoretical perspectives. Experimental angle-resolved spectra are measured to characterize both the dominant peaks and the dichroism of the polymers. First-principles calculations using the density functional theory-based many-body X-ray absorption spectroscopy (MBXAS) method are performed for the two polymers as well as for the thiophene and thienothiophene units that make up the conjugated backbones of these polymers. Through this combined approach, we are able to confidently assign the observed peaks to specific molecular orbitals and identify the orientation of their transition dipole moments (TDMs) with respect to the coordinate frame of the polymer backbone. Here, in particular, we are able to establish the character and orthogonal nature of the three main low-energy peaks at: (i) 2473.5 eV, 1s → (S–C)­π* with TDM along the π-stacking direction; (ii) 2474.1 eV, 1s → (S–C)­σ* with TDM along the backbone; and (iii) 2475.4 eV, 1s → (S–C)­σ* with TDM perpendicular to the first two. By performing both gas-phase and solid-state simulations, and with reference to the NEXAFS spectra of thiophene and thienothiophene building blocks, the influences of polymerization and molecular packing are also explored.

Chantler, Paul Alexander [Monash University, Clayt↗
DOE OSTI2021

Low-frequency Raman signature of Ag-intercalated few-layer MoS 2

Density functional theory based calculations and experimental analysis on a limited number of real samples are performed to study how the presence of silver intercalated in the van der Waals gap of few-layer MoS 2 affects the low-frequency Raman active modes of this material. Silver is found to predominantly affect the breathing-like and shear-like vibrational modes of MoS 2 . These modes correspond to quasi-rigid movements of each individual layer with a restoring force (and, in turn, frequency) that is determined by modulations in the weak interlayer interactions. Noticeable red-shifts with increasing Ag concentration are found for all low-frequency modes. This finding indicates the potential for low-frequency vibrations as useful gauges for practical determination of silver concentration using low-frequency Raman spectroscopy. Finally, this work also describes a semi-classical linear chain model that allows to extrapolate results to a large number of layers. Further, first-principles calculations show how Raman spectroscopy can be used to characterize the quality of the two-dimensional interface between MoS 2 and a silver substrate.

36 MATERIALS SCIENCE↗
DOE OSTI2026

IRMA

IRMA (In)elastic Representation of Materials As S(α,β) evaluations IRMA turns one phonon model into three outputs that usually require three separate tool chains: an evaluated nuclear-data file, predicted neutron-scattering spectra, and scattering kernels for Monte Carlo transport. The three outputs draw on a single, consistent description of the material, so the evaluation, the spectroscopy that can validate it, and the transport that uses it always agree about the physics. Nuclear data. IRMA writes ENDF-6 File 7 thermal scattering evaluations on automatically constructed (α, β) grids. This part reimplements and generalizes NJOY's LEAPR: the classic kernels reproduce freshly generated NJOY2016 tapes digit for digit and published reference tapes to about 1e-4, and the generalized paths add the exact coherent one-phonon term, anisotropic Debye-Waller tensors, coherent elastic for arbitrary crystals, and a per-species partition for polyatomic materials. The tapes feed NJOY, AMPX, FUDGE, and every transport code downstream of them. Neutron spectroscopy. The irma.spectra forward model projects the same physics onto an instrument's kinematics and resolution: INS spectra for VISION and generic indirect geometries, and 2-D S(Q,E) powder maps for direct-geometry spectrometers, from a phonopy model or straight from a phonon DOS. It can be used to predict a proposed measurement before beam time; in analysis, it supplies the calculated single-scattering counterpart of a measured spectrum, from the same material description the evaluation was built from. Monte Carlo transport. The irma.ncrystal exporter writes per-temperature scattering kernels for the companion NCrystal plugin, so McStas, OpenMC, and other NCrystal-aware codes sample the same physics. The exported kernels carry the per-site anisotropic Debye-Waller tensors, keeping directional coherent-elastic physics that NCrystal's standard scalar treatment does not represent. With the same physics inside a transport code, an entire beamline becomes a virtual experiment: IRMA's end-to-end validation ran a custom McStas implementation of the ARCS spectrometer, assembled from the existing McVine and McStas models, against measured data. From a bare crystal structure. The irma mlip front end builds the phonon model itself: a structure file and a choice of potential are enough. Nine pretrained machine-learned interatomic potentials are supported, on a laptop CPU, with no first-principles calculation; an approximate phonon model for a new material costs minutes, not a DFT campaign, and the build emits prefilled inputs for all three outputs. The result is a good starting point rather than a finished evaluation: survey-quality physics with every parameter exposed for review. A converged atomistic calculation enters the same way, as a phonopy model, when higher fidelity is needed.

Ramic, Kemal [Oak Ridge National Laboratory (ORNL)↗
DOE OSTI2020

Finite-Temperature Structures of Supported Subnanometer Catalysts Inferred via Statistical Learning and Genetic Algorithm-Based Optimization

Single-atom catalysts (SACs) minimize noble metal utilization and can alter the activity and selectivity of supported metal nanoparticles. However, the morphology of active centers, including single atoms and subnanometer clusters of a few atoms, remains elusive due to experimental challenges. The computational cost to describe numerous cluster shapes and sizes makes direct first-principles calculations impractical. We present a computational framework to enable structure determination for single-atom and subnanometer cluster catalysts. As a case study, we obtained the low energy structures of Pd n (n = 1-21) clusters supported on CeO 2 (111), which are critical components of automobile three-way catalysts. Trained on density functional theory data, a three-dimensional cluster expansion is established using statistical learning to describe the Hamiltonian and predict energies of supported Pdn clusters of any structure. Low energy stable and metastable structures are identified using a Metropolis Monte Carlo-based genetic algorithm in the canonical ensemble at 300 K. We observe that supported single atoms sinter to form bilayer clusters and large cluster isomers share similarities in both shape and energy, and elucidate the significance of the support and microstructure on cluster stability. We discovered a simple surrogate structure-energy model, where the energy per atom scales with the square root of the average first coordination number, which can be used to estimate energies and compare the stability of clusters. Our framework, applicable to any metal/support system, fills an important methodological gap to predict the stability of supported metal catalysts in the subnanometer regime.

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

Single-Layer Dititanium Oxide Ti 2 O MOene: Multifunctional Promises for Electride, Anode Materials, and Superconductor

Using the first-principles calculations, we report the existence of the single-layer (SL) dititanium oxide Ti 2 O (labeled as MOene) that constructs a novel family of MXene based on transition-metal oxides. This MOene material strongly contrasts the conventional ones consisting of transition-metal carbides and/or nitrides. SL Ti 2 O has high thermal and dynamical stabilities because of the strong Ti–O ionic bonding interactions. Moreover, this material is an intrinsic electride and exhibits extremely low diffusion barriers of ~12.0 and 6.3 meV for Li and Na diffusion, respectively. When applied as anode materials in lithium-ion batteries and sodium-ion batteries, it possesses a high energy storage capacity (960.23 mAhg –1 ), surpassing the traditional MXenes-based anodes. The superb electrochemical performance stems from the existing anionic electron on Ti 2 O surface. Astonishingly, SL Ti 2 O is also determined to be a superconductor with a superconducting transition temperature (T c ) of ~9.8 K, which originates from the soft-mode of the first acoustic phonon branch and enhanced electron–phonon coupling in the low-frequency region. Furthermore, this soft-mode behaves much softer upon applying a compressive strain of 2%, leading to a higher T c of 11.9 K. Furthermore, our finding broadens the family of MXenes and could facilitate more experimental efforts toward future nanodevices.

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

First-principles analysis of the Al-rich corner of Al-Li-Cu phase diagram

The phase diagram of Al-Li-Cu system in the Al-rich region was determined by means of first-principles calculations and statistical mechanics. The mixing enthalpies of many configurations for different lattices in the whole Al-Li-Cu system were determined by density functional theory simulations to find the stable phases in the convex hull. They were fitted with a cluster expansion to calculate the free energy of the configurations with different compositions as a function of temperature in the Al-rich region (Al content > 40 at.%) by means of Monte Carlo simulations. It was found that the ground state phases in the Al-rich part of the Al-Li-Cu phase diagram were α-Al, θ' (Al 2 Cu), $δ$' (Al 3 Li), $δ$ (AlLi) and T 1 (Al 6 Cu 4 Li 3 ), while θ'' (Al 3 Cu), T 1' (Al 2 CuLi) and Al 3 Cu 2 Li were found on the lowest mixing enthalpy surfaces of their lattices and were metastable. α-Al, $δ$ and T 1 are stable phases in the whole temperature range while $δ$' becomes metastable at very low temperature and θ (Al 2 Cu) replaces θ' as the stable phase at approximately 550 K due to the vibrational entropic contribution. In addition, the phase diagram in the Al-rich region was built and it was shown in isothermal sections from 100 to 900 K. They were in good agreement with the limited experimental data in the literature and provided new information regarding the stability, solubility and stoichiometry of the different phases. This information is important to understand the precipitation mechanisms during high temperature aging.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Development of Mechanistic Fission Gas Release and Swelling Models for UN Fuels in BISON

This report describes the work in NEAMS (Nuclear Energy Advanced Modeling and Simulation) to develop a mechanistic fission gas model for uranium nitride fuels in BISON. The existing Sifgrs (Simple integrated fission gas release and swelling) model tracks the average properties of two bubble populations in the bulk and at the grain boundaries. It was recognized that dislocations play a crucial role in the fission gas swelling of UN, and an irradiation-induced dislocation density model was needed, as well as a model describing how fission gas interacts with dislocations creating a third population of bubbles along dislocations. A mechanistic model that tracks dislocation bubbles was implemented in Sifgrs. This model was used to simulate fission gas swelling and release in UN. Lower-length-scale calculations and experimental observations from carbide fuel were leveraged to help populate the model with essential parameters. As a placeholder, an empirical function was formulated for the evolution of the dislocation network. Because of the difference in evolution of this network at different temperatures, the dislocation bubbles are able to capture behavior that the bulk intragranular bubbles cannot. It was found that the bulk bubbles dominate microscopic swelling at low temperatures, and the dislocation bubbles dominate at higher temperatures. Based on this model, the transition between the two bubble types is the main factor behind the breakaway swelling phenomenon in UN. In order to test the model, a couple of assessments were run. For the lower temperature JOYO pins, the model produces reasonable fission gas release and swelling values. For the higher temperature SP1 pin, the model dramatically underestimates the fission gas release. This issue is attributed to the gas being trapped inside dislocation bubbles, unable to escape to grain boundaries to cause fission gas release, and could be remedied by a mechanistic dislocation model that allows the dislocation density to decrease at very high temperatures. A preliminary mechanistic dislocation model was developed supported by first-principle calculations. These calculations provided valuable insight into how interstitial defects cluster in UN in the {110} orientation, and may eventually form dislocation loops, leading to the conclusion that dislocation loops may nucleate from these clusters. An estimate of the dislocation line energy in UN was obtained and will be improved in future work. The free-energy cluster dynamics code Centipede was used to track interstitial and vacancy absorption at dislocation loops and calculate their growth. In addition, improvements to Centipede were made including new convergence criteria for transient simulations and sink driving force updates.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

A New Family of Proton‐Conducting Electrolytes for Reversible Solid Oxide Cells: BaHf x Ce 0.8− x Y 0.1 Yb 0.1 O 3− δ

Abstract Reversible solid oxide cells based on ceramic proton conductors have potential to be the most efficient system for large‐scale energy storage. The performance and long‐term durability of these systems, however, are often limited by the ionic conductivity or stability of the proton‐conducting electrolyte. Here new family of solid oxide electrolytes, BaHf x Ce 0.8− x Y 0.1 Yb 0.1 O 3− δ (BHCYYb), which demonstrate a superior ionic conductivity to stability trade‐off than the state‐of‐the‐art proton conductors, BaZr x Ce 0.8− x Y 0.1 Yb 0.1 O 3− δ (BZCYYb), at similar Zr/Hf concentrations, as confirmed by thermogravimetric analysis, Raman, and X‐ray diffraction analysis of samples over 500 h of testing are reported. The increase in performance is revealed through thermodynamic arguments and first‐principle calculations. In addition, lab scale full cells are fabricated, demonstrating high peak power densities of 1.1, 1.4, and 1.6 W cm −2 at 600, 650, and 700 °C, respectively. Round‐trip efficiencies for steam electrolysis at 1 A cm −2 are 78%, 72%, and 62% at 700, 650, and 600 °C, respectively. Finally, CO 2 H 2 O electrolysis is carried out for over 700 h with no degradation.

Murphy, Ryan↗
DOE OSTI2021

Composition-Dependent Photoluminescence Properties and Anti-Counterfeiting Applications of A 2 AgX 3 (A = Rb, Cs; X = Cl, Br, I)

Copper(I) halides are emerging as attractive alternatives to lead halide perovskites for optical and electronic applications. However, blue-emitting all-inorganic copper(I) halides suffer from poor stability and lack of tunability of their photoluminescence (PL) properties. In this work, the preparation of silver(I) halides A 2 AgX 3 (A = Rb, Cs; X = Cl, Br, I) through solid-state synthesis is reported. In contrast to the Cu(I) analogs, A 2 AgX 3 are broad-band emitters sensitive to A and X site substitutions. First-principle calculations show that defect-bound excitons are responsible for the observed main PL peaks in Rb 2 AgX 3 and that self-trapped excitons (STEs) contribute to a minor PL peak in Rb 2 AgBr 3 . This is in sharp contrast to Rb 2 CuX 3 , in which the PL is dominated by the emission by STEs. Moreover, the replacement of Cu(I) with Ag(I) in A 2 AgX 3 significantly improves photostability and stability in the air under ambient conditions, which enables their consideration for practical applications. Thus, luminescent inks based on A 2 AgX 3 are prepared and successfully used in anti-counterfeiting applications. The excellent light emission properties, significantly improved stability, simple preparation method, and tunable light emission properties demonstrated by A 2 AgX 3 suggest that silver(I) halides may be attractive alternatives to toxic lead halide perovskites and unstable copper(I) halides for optical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Polaronic Conductivity in Cr 2 Ge 2 Te 6 Single Crystals

We report that intrinsic 2D ferromagnetic semiconductors are an important class of materials for spin-charge conversion applications. Cr 2 Ge 2 Te 6 retains long-range magnetic order in the bilayer at cryogenic temperatures and shows complex magnetic interactions with considerable magnetic anisotropy. Here, a series of structural, magnetic, X-ray scattering, electronic, thermal transport and first-principles calculation studies are performed, which reveal that localized electronic charge carriers in Cr 2 Ge 2 Te 6 are dressed by the surrounding lattice and are involved in polaronic transport via hopping that is observed via magnetocrystalline anisotropy. This opens the possibility for manipulation of charge transport in Cr 2 Ge 2 Te 6 —based devices by electron–phonon- and spin–orbit coupling-based tailoring of polaron properties.

2D materials↗
DOE OSTI2021

Accelerated Discovery and Design of Ultralow Lattice Thermal Conductivity Materials Using Chemical Bonding Principles

Semiconductors with very low lattice thermal conductivities are highly desired for applications relevant to thermal energy conversion and management, such as thermoelectrics and thermal barrier coatings. Although the crystal structure and chemical bonding are known to play vital roles in shaping heat transfer behavior, material design approaches of lowering lattice thermal conductivity using chemical bonding principles are uncommon. In this work, an effective strategy of weakening interatomic interactions and therefore suppressing lattice thermal conductivity based on chemical bonding principles is presented and a high-efficiency approach of discovering low κ L materials by screening the local coordination environments of crystalline compounds is developed. The resulting first-principles calculations uncover 30 hitherto unexplored compounds with (ultra)low lattice thermal conductivities from 13 prototype crystal structures contained in the Inorganic Crystal Structure Database. Furthermore, an approach of rationally designing high-performance thermoelectrics is demonstrated by additionally incorporating cations with stereochemically active lone-pair electrons. Here these results not only provide atomic-level insights into the physical origin of the low lattice thermal conductivity in a large family of copper/silver-based compounds but also offer an efficient approach to discover and design materials with targeted thermal transport properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Tuning Spin-Orbit Torques Across the Phase Transition in VO 2 /NiFe Heterostructure

The emergence of spin-orbit torques as a promising approach to energy-efficient magnetic switching has generated large interest in material systems with easily and fully tunable spin-orbit torques. We report current-induced spin-orbit torques in VO 2 /NiFe heterostructures are investigated using spin-torque ferromagnetic resonance, where the VO 2 layer undergoes a prominent insulator-metal transition. A roughly twofold increase in the Gilbert damping parameter, α, with temperature is attributed to the change in the VO 2 /NiFe interface spin absorption across the VO 2 phase transition. More remarkably, a large modulation (±100%) and a sign change of the current-induced spin-orbit torque across the VO 2 phase transition suggest two competing spin-orbit torque generating mechanisms. The bulk spin Hall effect in metallic VO 2 , corroborated by the first-principles calculation of the spin Hall conductivity σ SH ≈ -10 4 ($\frac{\hbar}{e}$) Ω -1 m -1 , is verified as the main source of the spin-orbit torque in the metallic phase. The self-induced/anomalous torque in NiFe, with opposite sign and a similar magnitude to the bulk spin Hall effect in metallic VO 2 , can be the other competing mechanism that dominates as temperature decreases. For applications, the strong tunability of the torque strength and direction opens a new route to tailor spin-orbit torques of materials that undergo phase transitions for new device functionalities.

36 MATERIALS SCIENCE↗