SEARCH · Search NASA

Results for “Observation nudging”

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.

16 records

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Comparing the interactions between particulate matter and cloud properties over two populated cities in Texas using WRF-Chem fine-resolution modeling

Accurate modeling of aerosol-cloud interactions is essential for reliable weather and air quality simulations, given their significant impact on precipitation patterns, cloud dynamics, and aerosol distributions. This study employed the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem) to examine the impact of enhanced meteorological simulations, achieved through advanced microphysics parameterization supported by data assimilation techniques, on air quality across Texas on August 19 and 20, 2022. We tested four distinct configurations: (1) the Morrison two-moment bulk microphysics scheme, (2) Morrison's with observation nudging, (3) the Spectral Bin Microphysics (SBM), and (4) SBM with observation nudging. While the SBM scheme is known for its detailed representation of aerosol-cloud interactions, our focus was on how improvements in meteorological accuracy translate to more precise air quality simulations. Our findings demonstrated a progressive improvement in simulation accuracy, starting with the Morrison's scheme and further enhanced by adopting the SBM scheme, complemented by incorporating observation nudging. Specifically, the combination of the SBM scheme and the nudging substantially enhanced the model's ability to capture convective precipitation events, as shown by better alignment with NEXRAD radar reflectivity, with R increasing from –0.21 to 0.82, IOA from 0.10 to 0.87, and NMB decreasing from 99% to 34% in Houston. The enhanced meteorology translated into more accurate PM 2.5 concentration simulations, particularly through the more accurate representation of aerosol washout during precipitation events. In Houston, the SBM scheme with nudging improved the model's PM 2.5 simulations significantly, with NMB decreasing from –20% to 5% and IOA improving from 0.43 to 0.61. In San Antonio, improvements were also notable, with NMB improved from –27% to –22%, R increased from 0.48 to 0.82, and IOA increased from 0.66 to 0.86. Furthermore, our results underscore the crucial role of accurate meteorological simulations in refining our understanding of aerosol behaviors in relation to precipitation patterns, directly enhancing the reliability and effectiveness of air quality modeling.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Biomass burning aerosol radiative effects in the Southeast Atlantic depend strongly on meteorological forcing method

Biomass burning aerosols (BBAs) from African fires may strongly impact Earth’s radiation budget in the southeast Atlantic (SEA), but the sign and magnitude of the overall radiative effect (RE) remain uncertain. Aerosol–climate models are needed to separately quantify direct, indirect, and semi-direct REs. Here, we evaluate improved simulations with the UK Met Office's Unified Model and explore REs resulting from various methods used to match observed meteorology (nudging or running forecasts reinitialized at different frequencies). REs are calculated as differences in radiative fluxes between simulations with and without smoke emissions and with and without aerosol absorption. All model setups agree on net warming for the SEA dominated by the direct effect. Simulated smoke, clouds, and the direct effect agree better with observations than previous studies using the same model, though biases in aerosol extinction and liquid water path remain. Changes in cloud droplet number concentration due to BBA self-lofting influence how cleanly we can separate cloud effects into semi-direct and indirect effects. Total RE, which remains unaffected, ranges from +3.0 to +7.9 W m −2 . The 4.9 W m −2 spread arises mainly from simulated semi-direct effects. Forecasts three days long or less probably do not allow time for plausible differences in boundary layer properties due to semi-direct effects to accumulate. Free running simulations with and without smoke accumulate differences in meteorology that are likely spurious “butterfly effects”. We recommend future research quantifying BBA REs over weeks to months to use meteorological forcing techniques that allow aerosol absorption to affect the boundary layer.

Giuffrida, Eric [Carnegie Mellon University, Pitts↗

DOE OSTI2025

Prediction of vacancy defect diffusion paths in high entropy alloys via machine learning on molecular dynamics data

Identifying the diffusion path of point defects is a critical step in understanding their evolution and the mechanisms of related phenomena. Defect diffusion occurs at small length and time scales, with impacts on material properties that may continue to evolve over ns to μs, ms, and the continuum scale (s, min, etc., and cm, m, etc.). The time scale accessible to molecular dynamics (MD) simulations is limited by small step sizes, typically in the fs range. Thus, surrogate models of MD simulations through machine learning (ML)-based algorithms are of great interest, especially for complex systems such as high entropy alloys (HEAs). In this work, dynamics governing vacancy migration in HEA were approximated with graph convolutional network (GCN) models as ansatzes for kinetic Monte Carlo (KMC) rate catalogs. Network design considered that diffusion in crystalline solids generally depends on interactions between defects and their immediate neighbor atoms. Graphs represented the vacancy surroundings, MD-generated trajectories provided training and comparison datasets, and unsupervised GCN models approximated interatomic dynamics governing vacancy migration in HEAs as ansatzes for KMC. A proof-of-concept model trained on MD data for the Fe, Ni, Cr, Co, and Cu HEA environment was used with two different neighbor interactions to assess the feasibility of training a GCN to predict vacancy defect transition rates in the HEA environment. The resulting setup rapidly generated MD-formatted synthetic trajectories based on dynamics learned from the MD training set, with a time acceleration of roughly two orders of magnitude and a similar diffusion coefficient to MD observations. Additionally, Nudged Elastic Band (NEB) calculations were performed on randomly generated FeNiCrCoCu HEA structures to determine vacancy migration barriers across nearest-neighbor sites. Transition probabilities for each jump, categorized by atomic type, were extracted from these calculations. NEB-based and GCN-based approaches led to similar outcomes.

Reimer, C↗

DOE OSTI2024

Changing effects of external forcing on Atlantic–Pacific interactions

Recent studies have highlighted the increasingly dominant role of external forcing in driving Atlantic and Pacific Ocean variability during the second half of the 20th century. This paper provides insights into the underlying mechanisms driving interactions between modes of variability over the two basins. We define a set of possible drivers of these interactions and apply causal discovery to reanalysis data, two ensembles of pacemaker simulations where sea surface temperatures in either the tropical Pacific or the North Atlantic are nudged to observations, and a pre-industrial control run. We also utilize large-ensemble means of historical simulations from the Coupled Model Intercomparison Project Phase 6 (CMIP6) to quantify the effect of external forcing and improve the understanding of its impact. A causal analysis of the historical time series between 1950 and 2014 identifies a regime switch in the interactions between major modes of Atlantic and Pacific climate variability in both reanalysis and pacemaker simulations. A sliding window causal analysis reveals a decaying El Niño–Southern Oscillation (ENSO) effect on the Atlantic as the North Atlantic fluctuates towards an anomalously warm state. The causal networks also demonstrate that external forcing contributed to strengthening the Atlantic's negative-sign effect on ENSO since the mid-1980s, where warming tropical Atlantic sea surface temperatures induce a La Niña-like cooling in the equatorial Pacific during the following season through an intensification of the Pacific Walker circulation. The strengthening of this effect is not detected when the historical external forcing signal is removed in the Pacific pacemaker ensemble. The analysis of the pre-industrial control run supports the notion that the Atlantic and Pacific modes of natural climate variability exert contrasting impacts on each other even in the absence of anthropogenic forcing. The interactions are shown to be modulated by the (multi)decadal states of temperature anomalies of both basins with stronger connections when these states are “out of phase”. We show that causal discovery can detect previously documented connections and provides important potential for a deeper understanding of the mechanisms driving changes in regional and global climate variability.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

The Role of Land‐Atmosphere Feedbacks in Midlatitude Wintertime Surface Temperature Variability

Accurately representing synoptic near-surface temperature variability is crucial to predict weather extremes, yet models remain biased. Existing studies primarily attribute wintertime midlatitude near-surface temperature variability to tropospheric large-scale advection, assuming minimal land influence. However, nudging the model's circulation toward observations yields little improvement in wintertime temperature variance over Northern Hemisphere land, suggesting that land-atmosphere interactions also warrant attention. We introduce a new scaling framework for temperature variance that incorporates local land-atmosphere feedbacks. Comparing our framework to the mixing length approach—which links temperature variance to the meridional temperature gradient and air parcel displacement (mixing length)—shows that land-atmosphere feedbacks are inherently embedded in the mixing length, a connection previously overlooked. Roles of land–atmosphere feedbacks are evaluated via model experiments with perturbed circulation, land, or both. We find that longwave radiative damping dominates temperature variance responses over meridional temperature gradient when both land and circulation are perturbed.

atmosheric science↗

DOE OSTI2026

Improving thermodynamic nudging in the E3SM Atmosphere Model version 2 (EAMv2): strategy and hindcast skills on weather systems

Nudging techniques are commonly employed to constrain atmospheric simulations toward observed states, facilitating model evaluation and sensitivity studies. However, if applied improperly – particularly to thermodynamic variables such as temperature and humidity – nudging can distort physical processes and introduce spurious biases, undermining the credibility of the simulations. This study presents an improved nudging implementation that applies vertically modulated tendencies to reduce adverse impacts on model physics. The framework is tested in version 2 of the Energy Exascale Earth System Model (EAMv2) using a suite of hindcast simulations nudged toward ERA5 reanalysis. We systematically evaluate the individual and combined effects of nudging wind, temperature, and humidity fields on the model's ability to represent large-scale atmospheric states and high-impact weather systems. Results show that the revised strategy – particularly when nudging temperature and humidity at selected levels – enhances hindcast skill by improving agreement with ERA5 without degrading the hydrological cycle or precipitation processes. Additional improvements in surface temperature, outgoing longwave radiation, and precipitation biases are achieved through targeted nudging of land-surface variables. The proposed approach strengthens the representation of large-scale conditions relevant to tropical cyclones, atmospheric rivers, and extratropical cyclones in the low-resolution EAMv2. These findings demonstrate that carefully designed thermodynamic nudging, especially of temperature and humidity, improves the realism of constrained simulations and broadens the utility of nudged EAMv2 for atmospheric modeling, machine learning, and high-impact weather research.

Atmospheric river↗

DOE OSTI2025

Hunga Tonga–Hunga Ha′apai Volcano Impact Model Observation Comparison (HTHH-MOC) project: experiment protocol and model descriptions

The 2022 Hunga volcanic eruption injected a significant amount of water vapor and a moderate amount of sulfur dioxide into the stratosphere, causing observable responses in the climate system. We have developed a model–observation comparison project to investigate the evolution of volcanic water and aerosols and their impacts on atmospheric dynamics, chemistry, and climate, using several state-of-the-art chemistry climate models. The project goals are (1) to evaluate the current chemistry–climate models to quantify their performance in comparison to observations and (2) to understand atmospheric responses in the Earth system after this exceptional event and investigate the potential impacts in the projected future. To achieve these goals, we designed specific experiments for direct comparisons to observations, for example from balloons and the Microwave Limb Sounder satellite instrument. Experiment 1 consists of two sets of free-running ensemble experiments from 2022 to 2031: one with fixed sea-surface temperatures and sea ice and one with coupled ocean. These experiments will help to understand the long-term evolution of water vapor and aerosols; quantify HTHH effects on stratospheric and mesospheric temperatures, dynamics, and transport; understand the impact of dynamic changes on ozone chemistry; quantify the net radiative forcings; and evaluate any surface climate impact. Experiment 2 is a nudged-run experiment from 2022 to 2023 using observed meteorology. To allow participation of more climate models with varying complexities of aerosol simulation, we include two sets of simulations in Experiment 2: Experiment 2a is designed for models with internally generated aerosol, while Experiment 2b is designed for models using prescribed aerosol surface area density. This experiment will help to analyze H 2 O and aerosol evolution, quantify the net radiative forcings, understand the impacts on mid-latitude and polar O 3 chemistry, and allow close comparisons with observations.

Zhu, Yunqian [Univ. of Colorado, Boulder, CO (Unit↗

DOE OSTI2025

Water sources and land capacitor effects stimulate observed summer Arctic moistening and warming

The primary sources of recent summer Arctic moistening trends in reanalysis are uncertain, hindering attribution of observed Arctic warming due to radiative effects from water vapor changes. Here, we use a combined online numerical water tracer and circulation nudging approach in the Community Earth System Model to track the sources of water vapor beyond its initial sources. Trends in boreal summer large-scale circulation have driven moistening of the Arctic over recent decades, having a large impact on the Arctic radiative budget, accounting for 94% of the strengthening water vapor radiative feedback. We identify two key regions supplying the Arctic water vapor feedback: Northeast North America and western/central Eurasia. In both regions, anticyclonic circulations over the southwest Atlantic and eastern Europe move moisture from the tropical oceans poleward to high latitude land through precipitation in winter and spring. During summer, evapotranspiration over land releases this water vapor, and it is transported by winds into the Arctic. We refer to this sequence of terrestrial moisture storage and release as the land capacitor effect. Thus, the impacts of circulation changes on poleward moisture transport and land-atmosphere interactions over high latitudes represent the underlying mechanisms of the recent moistening and warming in the Arctic.

Arctic amplification↗

DOE OSTI2025

Regime-based aerosol–cloud interactions from CALIPSO-MODIS and the Energy Exascale Earth System Model version 2 (E3SMv2) over the Eastern North Atlantic

This study investigates aerosol-cloud interactions in marine boundary layer (MBL) clouds using an advanced deep-learning-driven synoptic-regime-based framework, combining satellite data (CALIPSO vertically resolved aerosol extinction and MODIS cloud properties) with 1° nudged Energy Exascale Earth System Model version 2 (E3SMv2) simulation over the Eastern North Atlantic (ENA; ∼10°×10°, 2006–2014). The E3SMv2 captures observed seasonal variations in cloud droplet number concentrations (N d ) and liquid water path (LWP), though it systematically underestimates N d . We then partition ENA meteorology into four synoptic regimes (Pre-Trough, Post-Trough, Ridge, Trough) via a deep-learning clustering of ERA5 reanalysis fields, enabling regime-dependent aerosol-cloud interactions analyses. Both satellite and E3SMv2 exhibit an inverted-V LWP-N d relationship. In Post-Trough and Ridge regimes, the satellite shows stronger negative LWP-N d sensitivities than in Pre-Trough regime. The Trough regime displays a muted satellite LWP response. In comparison, the model predicts more exaggerated LWP responses across regimes, with LWP increasing too quickly at low N d and decreasing more sharply at high N d , especially in Pre-Trough and Trough regimes. These exaggerated model LWP sensitivities may stem from uncertainties in representing drizzle processes, entrainment, and turbulent mixing. As for N d susceptibility to aerosols, N d increases with MBL aerosol extinction in both datasets, but the simulated aerosol-cloud interactions appear oversensitive to meteorological conditions. Overall, E3SMv2 better captures aerosol effects under regimes that favor stratiform clouds (Post-Trough, Ridge), but performance deteriorates for regimes with deeper, dynamically complex clouds (Trough), highlighting the need for improved representations of those cloud processes in climate models.

Environmental sciences↗

DOE OSTI2024

How well are hazards associated with derechos reproduced in regional climate simulations?

Abstract. A 15-member ensemble of convection-permitting regional simulations of the fast-moving and destructive derecho of 29–30 June 2012 that impacted the northeastern urban corridor of the USA is presented. This event generated 1100 reports of damaging winds, generated significant wind gusts over an extensive area of up to 500 000 km2, caused several fatalities, and resulted in widespread loss of electrical power. Extreme events such as this are increasingly being used within pseudo-global-warming experiments to examine the sensitivity of historical, societally important events to global climate non-stationarity and how they may evolve as a result of changing thermodynamic and dynamic contexts. As such it is important to examine the fidelity with which such events are described in hindcast experiments. The regional simulations presented herein are performed using the Weather Research and Forecasting (WRF) model. The resulting ensemble is used to explore simulation fidelity relative to observations for wind gust magnitudes, spatial scales of convection (as is manifest in high composite reflectivity, cREF), and both rainfall and hail production as a function of model configuration (microphysics parameterization, lateral boundary conditions (LBCs), start date, use of nudging, compiler choice, damping, and number of vertical levels). We also examine the degree to which each ensemble member differs with respect to key mesoscale drivers of convective systems (e.g., convective available potential energy and vertical wind shear) and critical manifestations of deep convection, e.g., vertical velocities, cold-pool generation, and how those properties relate to the correct characterization of the associated atmospheric hazards (wind gusts and hail). Use of a double-moment, seven-class scheme with number concentrations for all species (including hail and graupel) results in the greatest fidelity of model-simulated wind gusts and convective structure to the observations of this event. All ensemble members, however, fail to capture the intensity of the event in terms of the spatial extent of convection and the production of high near-surface wind gusts. We further show very high sensitivity to the LBCs employed and specifically that simulation fidelity is higher for simulations nested within ERA-Interim compared to ERA5. Excess convective available potential energy (CAPE) in all ensemble members after the derecho passage leads to excess production of convective cells, wind gusts, cREF > 40 dBZ, and precipitation during a frontal passage on the subsequent day. This event proved very challenging to forecast in real time and to reproduce in the 15-member hindcast simulation ensemble presented here. Future work could examine if simulations with other initial and lateral boundary conditions can achieve greater fidelity.

Shepherd, Tristan (ORCID:0000000186276419)↗

DOE OSTI2026

Leaf Wax Hydrogen Isotopes Reflect Storm Track Position Over Western North America

Proxy records of past climates have yielded powerful insights into regional hydroclimate dynamics. Novel insights may be gained by reconstructing spatiotemporal changes in precipitation isotopologues in past climate states that cannot be gleaned from individual site‐level studies. In this paper, we ask whether latitudinal gradients in the stable hydrogen isotopic composition of precipitation, as inferred from sedimentary leaf wax biomarkers, reflect aspects of large‐scale climate conditions in western North America (WNA). Modern coretop samples from offshore WNA show that leaf wax hydrogen isotopes broadly track the hydrogen isotopic composition of rainfall between 20 and 40N, but with an offset corresponding to a relatively constant “apparent fractionation” value. Poleward of 40N, the leaf wax signal may be complicated by fluvial transport of leaf waxes from the continental interior. Leaf wax‐inferred precipitation hydrogen isotopes show a shift to more negative values between 30 and 40N. We use modern observational data to show that this latitudinal range marks a major transition between the subtropical arid zone and the location of the midlatitude storm tracks over the northeast Pacific. Nudged water isotope‐enabled models capture the location of the arid to mesic climate transition with fidelity. This modern data set suggests that reconstructions of the latitudinal gradient of precipitation hydrogen isotopes can constrain the sensitivity of the storm tracks to shifts in climatic boundary conditions in past climate states. This approach can yield novel insights into past, present, and future hydroclimate variability in this arid region.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Evaluation of Marine Boundary Layer Clouds Over the Northeast Pacific During the CSET Campaign in E3SM Version 2

It is still challenging to reproduce marine boundary layer (MBL) clouds well in large-scale models despite their importance to the Earth's radiation budget and hydrological cycle. This study evaluates representation of the MBL and clouds in the Energy Exascale Earth System Model (E3SM) version 2. This study compares the E3SM simulation results with remote sensing and reanalysis data during the Cloud System Evolution in the Trades (CSET) field campaign to better understand the stratocumulus-to-cumulus cloud transition (SCT) over the northeast Pacific. E3SM results are extracted along the CSET Lagrangian trajectories. The comparison shows that the E3SM simulation applying horizontal wind nudging performs well in reproducing thermodynamic variables of the MBL and evolution trends of cloud variables along the trajectories. However, substantial overestimations of aerosol and cloud drop number (N d ) are observed, which is explained as an issue with version 2 of the model. Cloud fraction (CF) does decrease from the Californian coast to Hawaii in the E3SM simulation, but most CF values indicate either an overcast or almost clear sky, which differs from satellite and reanalysis data. The effect of N d overestimation on CF evolution is assessed via prescribed-N d simulations. Those simulations with N d modifications show negligible CF changes. A comparison of estimated inversion strength (EIS) also shows that the simulated EIS values are similar to those of reanalysis data. Our study suggests that cloud macrophysics and boundary layer processes are more important in improving the simulation to capture the SCT than refining the model's thermodynamics or cloud microphysics.

Choi, Kyoung Ock [Univ. of Washington, Seattle, WA↗

DOE OSTI2025

Data-Driven Closures and Assimilation for Stiff Multiscale Random Dynamics

Here, we introduce a data-driven and physics-informed framework for propagating uncertainty in stiff, multiscale random ordinary differential equations (RODEs) driven by correlated (colored) noise. Unlike systems subjected to Gaussian white noise, a deterministic equation for the joint probability density function (PDF) of RODE state variables does not exist in closed form. Moreover, such an equation would require as many phase-space variables as there are states in the RODE system. To alleviate this curse of dimensionality, we instead derive exact, albeit unclosed, reduced-order PDF (RoPDF) equations for low-dimensional observables/quantities of interest. The unclosed terms take the form of state-dependent conditional expectations, which are directly estimated from data at sparse observation times. However, for systems exhibiting stiff, multiscale dynamics, data sparsity introduces regression discrepancies that compound during RoPDF evolution. This is overcome by introducing a kinetic-like defect term to the RoPDF equation, which is learned by assimilating in sparse, low-fidelity RoPDF estimates. Two assimilation methods are considered, namely nudging and deep neural networks, which are successfully tested against Monte Carlo simulations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Low-energy pathways lead to self-healing defects in CsPbBr 3

Self-regulation of free charge carriers in perovskites via Schottky defect formation has been posited as the origin of the well-known defect-tolerance of metal halide perovskite materials. Understanding the mechanisms of self-regulation promises to lead to the fabrication of better performing solar cell materials with higher efficiencies. We investigate many such mechanisms here for CsPbBr 3 , a popular representative of a more commercially viable all-inorganic metal halide perovskite. We investigate different atomic-level mechanisms and pathways of the diffusion and recombination of neutral and charged interstitials and vacancies (Schottky pairs) in CsPbBr 3 . We use nudged elastic band calculations and ab initio-derived pseudopotentials within quantum ESPRESSO to determine energies of formation and migration and hence the activation energies for these defects. While halide vacancies are known to exhibit low formation energies, the migration of interstitials is less studied. Our calculations uncover interstitial defect pathways capable of producing an activation energy at, or below, the single experimental value of 0.53 eV observed for the slow, temperature-dependent recovery of light-induced conductivity in bulk CsPbBr 3 . Our work reveals the existence of a low-energy diffusion pathway involving a concerted “domino effect” of interstitials, with the net result that interstitials can diffuse more readily over longer distances than expected. This observation suggests that defect self-healing can be promoted if the “domino effect” strategy can be engaged.

14 SOLAR ENERGY↗

DOE OSTI2024

Collective Neutrino Oscillations and Heavy-element Nucleosynthesis in Supernovae: Exploring Potential Effects of Many-body Neutrino Correlations

In high-energy astrophysical processes involving compact objects, such as core-collapse supernovae or binary neutron star mergers, neutrinos play an important role in the synthesis of nuclides. Neutrinos in these environments can experience collective flavor oscillations driven by neutrino–neutrino interactions, including coherent forward scattering and incoherent (collisional) effects. Recently, there has been interest in exploring potential novel behaviors in collective oscillations of neutrinos by going beyond the one-particle effective or "mean-field" treatments. Here, we seek to explore implications of collective neutrino oscillations, in the mean-field treatment and beyond, for the nucleosynthesis yields in supernova environments with different astrophysical conditions and neutrino inputs. We find that collective oscillations can impact the operation of the νp-process and r-process nucleosynthesis in supernovae. The potential impact is particularly strong in high-entropy, proton-rich conditions, where we find that neutrino interactions can nudge an initial νp-process neutron-rich, resulting in a unique combination of proton-rich low-mass nuclei as well as neutron-rich high-mass nuclei. We describe this neutrino-induced neutron-capture process as the "νi-process." In addition, nontrivial quantum correlations among neutrinos, if present significantly, could lead to different nuclide yields compared to the corresponding mean-field oscillation treatments, by virtue of modifying the evolution of the relevant one-body neutrino observables.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Lattice vacancy migration barriers in Fe-Ni alloys, and an indication as to why Ni atoms diffuse slowly: A first-principles study

Lattice vacancy migration barriers in ferromagnetic Fe 𝑥 ⁢Ni 1−𝑥 alloys (0.4 ≤ 𝑥 ≤ 0.6) are accurately quantified within the framework of ab initio electronic structure calculations using the nudged elastic band (NEB) method. Both the atomically disordered (A1) fcc phase, as well as the atomically ordered, tetragonal L⁢1 0 phase—which is under consideration as a material for a rare-earth-free gap magnet for advanced engineering applications—are investigated. Across an ensemble of NEB calculations performed on supercell configurations spanning a range of compositions and containing disordered, partially ordered, and fully ordered structures, we find that Ni-vacancy interchanges encounter significantly higher energetic barriers than do Fe-vacancy interchanges. We contend that this aspect is a key factor in determining the differences in mobility between Fe and Ni atoms in this ferromagnetic alloy. Moreover, we are able to interpret these findings in terms of the ferromagnetic alloy's underlying spin-polarized electronic structure. Specifically, we report a coupling between the size of local lattice distortions and the magnitude of the local electronic spin polarization around vacancies. This causes Fe atoms to relax into lattice vacancies, while Ni atoms remain rigidly fixed to their original lattice positions. These results give atomic-scale insight into the longstanding experimental observation that Ni exhibits remarkably slow atomic diffusion in Fe-Ni alloys.

density functional theory↗