Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT Calculations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Davis Computational Spectroscopy Workflow—From Structure to Spectra

Here, we describe an automated workflow that connects a series of atomic simulation tools to investigate the relationship between atomic structure, lattice dynamics, materials properties, and inelastic neutron scattering (INS) spectra. Starting from the atomic simulation environment (ASE) as an interface, we demonstrate the use of a selection of calculators, including density functional theory (DFT) and density functional tight binding (DFTB), to optimize the structures and calculate interatomic force constants. We present the use of our workflow to compute the phonon frequencies and eigenvectors, which are required to accurately simulate the INS spectra in crystalline solids like diamond and graphite as well as molecular solids like rubrene. We have also implemented a machine-learning force field based on Chebyshev polynomials called the Chebyshev interaction model for efficient simulation (ChIMES) to improve the accuracy of the DFTB simulations. We then explore the transferability of our DFTB/ChIMES models by comparing simulations derived from different training sets. We show that DFTB/ChIMES demonstrates ~100× reduction in computational expense while retaining most of the accuracy of DFT as well as yielding high accuracy for different materials outside of our training sets. The DFTB/ChIMES method within the workflow expands the possibilities to use simulations to accurately predict materials properties of increasingly complex structures that would be unfeasible with ab initio methods.

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

Understanding Lithium Local Environments in LiMn 0.5 Ni 0.5 O 2 Cathodes: A DFT-Supported 6 Li Solid-State NMR Study

LiMn 0.5 Ni 0.5 O 2 has been plagued by the structural failure driven by the irreversible phase change upon subsequent cycling and the structural degradation caused by the seemingly unavoidable Li/Ni exchange. This begs for a better understanding of the correlation between lithium local structures, synthesis conditions, and overall electrochemical properties of materials. In this study, we use density functional theory (DFT) to assist experimental 6 Li solid-state magic angle spinning (MAS) nuclear magnetic resonance (NMR) spectroscopy to understand the nature and change in lithium local structures with extended annealing times (168 vs 15 h) synthesized at 900 °C. To correlate the NMR peak changes after different synthesis conditions, we calculated the NMR spectra of representative low-energy model configurations of LiMn 0.5 Ni 0.5 O 2 composition–zigzag, row, and flower structures–as well as the influence of Ni/Mn ordering and Li/Ni exchange on each spectrum using a supercell that is large enough to allow us to consider the effect of different relative configurations of adjacent transition metal (TM) layers. The analysis based on a combination of the calculated NMR and the experimental spectra suggests that most configurations contributing to the experimental spectra are in fact composed of a blend of configurations higher in energy than the classical well-ordered ground-state structures. Further, the extended annealing of LiMn 0.5 Ni 0.5 O 2 slightly enhances Ni/Mn and Li/Mn orderings while reducing the Li/Ni mixing ratio. This study shows how DFT calculations are crucial in providing a better understanding of the experimental characterization data and therefore local environments and domain structures in oxide materials, such as Li-ion battery cathodes.

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

Electronic structure of LaNiO 2 and CaCuO 2 from a self-consistent vertex-corrected GW approach

In this work, the electronic structure of one of the nickelates (LaNiO 2 ) and one of the cuprates (CaCuO 2 ) is studied with three self-consistent GW-based methods: scGW, sc(GW+vertex), and quasiparticle self-consistent GW. Low-energy features obtained in our study are in many respects similar to the features reported in previous density functional theory plus dynamical mean-field theory (DFT+DMFT) studies. Consistent with the DFT+DMFT conclusion, we find LaNiO 2 to be more correlated than CaCuO 2 . However, correlation effects included in our study change the DFT Fermi surface near the Γ point differently from that reported in DMFT studies. Features that are a few electronvolts away from the Fermi level are broader in our calculations than in the DFT+DMFT, which reflects the differences between the DFT and the GW methods. Our results are in qualitative agreement with previous G0W0 results, but the self-consistency results in quantitative differences. Generally, correlation effects are found to be sufficiently weak in both materials, which allows one to use totally ab initio diagrammatic approaches such as sc(GW+Vertex) and to avoid the methods with adjustable parameters (DFT+U or DFT+DMFT). However, the possibility of some strong correlations at low energy that cannot be captured by perturbative methods cannot be completely excluded. For instance, differences in the Fermi surface should be resolved, thus experimental studies are necessary.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

High temperature elastic properties of sub-stoichiometric yttrium dihydrides

Yttrium hydrides are considered as candidate materials for neutron moderation applied in microreactors (akin transportable tiny nuclear reactors) owing to their superior thermal stability and hydrogen retention. The evolution of elastic properties of these materials at elevated temperatures, needed for predicting the thermomechanical response and performance of the moderator during in-service reactor conditions, however, is lacking. Here, we report the Young’s and shear elastic moduli of three stoichiometries of bulk yttrium hydride (YH x , x = 1.61, 1.82, and 1.84) from room temperature to 1000°C. In situ temperature-dependent measurements of the longitudinal and shear wave velocities were performed using a laser ultrasonic technique while heating the sample in a vacuum-pumped heating stage. The elastic moduli increased linearly with increasing hydrogen content and decreased by ~10% during heating from room temperature to 1000°C in the three YH x compositions. The linear relationship between the elastic moduli and the hydrogen content in yttrium hydride was verified by atomistic calculations based on density functional theory (DFT). The absence of abrupt changes in the temperature-dependent measurements of elastic modulus of the YH x samples suggested negligible loss of hydrogen at elevated temperatures. Excellent agreement was found between the measured and calculated dependence of the elastic moduli on the stoichiometry, thereby providing a new approach for investigating the effects of fabrication-induced parameters (such as porosity) on the elastic moduli. Furthermore, this study demonstrates the utility of the combined approach involving DFT-based atomistic calculations and measurements of the elastic moduli for the informative development of metal hydrides and can be used as a metric for novel moderator materials investigations for emerging microreactors and beyond.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Effects of spin-orbit coupling and very large supercells on the description of acceptors in CdTe

Predicting accurate shallow donor and acceptor levels in semiconductors has been quite challenging using periodic boundary conditions as implemented in current density functional theory codes. The reason is that the wave functions associated with the shallow centers are quite extended and are not fully contained in the supercell (typically with a few hundred atoms) so impurities in the periodically repeated image cells interact with each other. Errors of ~0.1 eV are expected, and these are of the same order of magnitude as the ionization energies themselves. In the case of acceptors in CdTe, this problem is exacerbated by the strong spin-orbit coupling that split the Te-related states at the top of the valence band, making the calculations at least 8 times more expensive. CdTe is an important solar-cell material with record high efficiency of 22%. One of the main limiting factors to increasing the efficiency towards the theoretical limit of ~30% is the often reported very low hole concentration. Pushing to the limit of computational capability by using very large supercells with spin-orbit coupling we report the results of hybrid functional calculations of group-V acceptors in CdTe. We show that extrapolation to the dilute limit leads to an interpretation of the experimental data that is qualitatively different from previous DFT and hybrid functional calculation reports. We find that the group-V impurities indeed behave as shallow acceptors and that the corresponding compensating AX-centers are unstable and do not limit p-type doping. We address the differences between our results and previous theoretical predictions and show that our calculated ionization energies predict hole concentrations that are in excellent agreement with recent temperature-dependent Hall measurements on high-quality single-crystal samples.

14 SOLAR ENERGY↗
DOE OSTI2022

A flexible and scalable scheme for mixing computed formation energies from different levels of theory

Abstract Computational materials discovery efforts are enabled by large databases of properties derived from high-throughput density functional theory (DFT), which now contain millions of calculations at the generalized gradient approximation (GGA) level of theory. It is now feasible to carry out high-throughput calculations using more accurate methods, such as meta-GGA DFT; however recomputing an entire database with a higher-fidelity method would not effectively leverage the enormous investment of computational resources embodied in existing (GGA) calculations. Instead, we propose here a general procedure by which higher-fidelity, low-coverage calculations (e.g., meta-GGA calculations for selected chemical systems) can be combined with lower-fidelity, high-coverage calculations (e.g., an existing database of GGA calculations) in a robust and scalable manner. We then use legacy PBE(+ U ) GGA calculations and new r 2 SCAN meta-GGA calculations from the Materials Project database to demonstrate that our scheme improves solid and aqueous phase stability predictions, and discuss practical considerations for its implementation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

First-principles investigation of the thermophysical properties of NaCl, PuCl 3 , and NaCl-PuCl 3 Molten salts

Molten salts have a variety of applications that span the nuclear and solar industries, and which involve thermal storage and heat transfer. There is a present knowledge gap in the thermophysical properties of molten salts, which limits the readiness level of molten salt applications. This is especially pertinent for molten salt reactors where the fissile material is dissolved within the molten salt. Ab initio molecular dynamics is a common method employed to investigate the structural and thermophysical properties at elevated temperatures. NaCl-PuCl 3 is a candidate fuel salt in molten salt reactors. The scope of this study is to investigate the seven unique compositions and a range of temperatures of NaCl-PuCl 3 and calculate the density, heat capacity, compressibility, enthalpy of mixing, and coefficient of thermal expansion. Within this work, the van der Waals (vdW) interactions are primarily handled with the vdW-DF2 functional but spot-checked with both the dDsC and DFT-D3 methods. The calculated densities are in good agreement with the NaCl and the eutectic compositions, and no experimental values exist for PuCl 3 at these temperatures. The densities are fit to a first-order Redlich-Kister expansion as a function of both composition and temperature. The heat capacity in the literature is scattered but the calculated values agree well with one of the experimental results. Here, the heat capacity increases as a linear function with respect to concentration at a rate of about 7.5 J/mol-K per 10 mol% PuCl 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2020

Electronic energy gap closure and metal-insulator transition in dense liquid hydrogen

Here, using quantum Monte Carlo (QMC) calculations, we investigate the insulator-metal transition observed in liquid hydrogen at high pressure. Below the critical temperature of the transition from the molecular to the atomic liquid, the fundamental electronic gap closure occurs abruptly, with a small discontinuity reflecting the weak first-order transition in the thermodynamic equation of state. Above the critical temperature, molecular dissociation sets in while the gap is still open. When the gap closes, the decay of the off-diagonal reduced density matrix shows that the liquid enters a gapless, but localized, phase: there is a crossover between the insulating and the metallic liquids. Compared to different density functional theory (DFT) functionals, our QMC calculations provide larger values for the fundamental gap and the electronic density of states close to the band edges, indicating that optical properties from DFT potentially benefit from error cancellations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Developing reliable machine learning interatomic potential for Fe–Cr–Ni austenitic alloys

Gaining atomistic understanding of mechanical behavior of heat-resistant structural materials such as Fe–Cr–Ni-based alloys requires an approach with an accuracy close to density functional theory (DFT) that considers the intrinsic properties of the bulk lattice and important defects such as stacking faults, grain boundaries, and surfaces. This work aims to develop reliable machine learning interatomic potential (MLIAP) at cross-scale for Fe–Cr–Ni ternary alloys with a focus on the face-centered-cubic (fcc) solid solution structure. Leveraging the advantages of moment tensor potentials, which typically necessitate a relatively small training dataset and enable rapid calculations using the large-scale atomic/molecular massively parallel simulator package, we ensure the stability and accuracy of the trained potentials. Important defects such as stacking faults, grain boundaries, and surfaces for wide-range compositions are investigated. Structural, thermal, elastic, and defect properties are determined from molecular dynamics simulations comprising several thousand atoms, generated via canonical Monte Carlo simulations guided by the trained potential. The trained potential allows efficient atomic simulations of structural, thermal, and mechanical properties of fcc Fe–Cr–Ni solid solution alloys as a function of composition and temperature. Therefore, the MLIAP approach represents a major advancement from DFT calculations that are limited to small simulation sizes and traditional molecular dynamics simulations using relatively low accuracy potentials. Furthermore, this work outlines a practical foundation for further investigating the structural evolution and mechanical behavior of austenitic stainless steel and nickel-based alloys in a wide array of applications in extreme environments.

Crystal structure↗
DOE OSTI2020

PyCDFT: A Python package for constrained density functional theory

In this paper, we present PyCDFT, a Python package to compute diabatic states using constrained density functional theory (CDFT). PyCDFT provides an object-oriented, customizable implementation of CDFT, and allows for both single-point self-consistent-field calculations and geometry optimizations. PyCDFT is designed to interface with existing density functional theory (DFT) codes to perform CDFT calculations where constraint potentials are added to the Kohn–Sham Hamiltonian. Here, we demonstrate the use of PyCDFT by performing calculations with a massively parallel first-principles molecular dynamics code, Qbox, and we benchmark its accuracy by computing the electronic coupling between diabatic states for a set of organic molecules. We show that PyCDFT yields results in agreement with existing implementations and is a robust and flexible package for performing CDFT calculations. The program is available at https://dx.doi.org/10.5281/zenodo.3821097.

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

Møller–Plesset and Density-Fixed Adiabatic Connections for a Model Diatomic System at Different Correlation Regimes

In recent years, adiabatic connection (AC) interpolations developed within density functional theory (DFT) have been found to provide good performances in the calculation of interaction energies when used with Hartree–Fock (HF) ingredients. The physical and mathematical reasons for such unanticipated performance have been clarified, to some extent, by studying the strong-interaction limit of the Møller–Plesset (MP) AC. In this work, we calculate both the MP and the DFT AC integrand for the asymmetric Hubbard dimer, which allows for a systematic investigation of different correlation regimes by varying two simple parameters in the Hamiltonian: the external potential, Δv, and the interaction strength, U. Notably, we find that, while the DFT AC integrand appears to be convex in the full parameter space, the MP integrand may change curvature twice. Furthermore, we discuss different aspects of the second-order expansion of the correlation energy in each AC, and we demonstrate why the derivative of the λ-dependent density in the MP AC at λ = 0 (i.e., at the HF density) is zero in the model. Concerning the strong-interaction limit of both ACs in the Hubbard dimer setting, we show that the asymptotic value of the MP AC, W ∞ HF , is lower than (or equal to) its DFT analogue, W ∞ KS , if the two are compared at a given density, just like in real space. Furthermore, we also show that this is not always the case if the two quantities are compared at a given external potential.

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

Energy and mechanical properties predictions in Fe-Ni binary system by ab initio calculations

Here, the fcc random solid solutions in Fe-Ni binary system were investigated using ab initio calculations. Density functional theory (DFT) was applied along with supercells generated by the Special Quasi-Random Structure (SQS) method to predict the enthalpy of mixing and magnetic and elastic properties. The results were compared with relevant experimental and computational data produced by other groups. The effect of SQS supercell size and relaxation procedure on the predicted properties of the alloys were also studied. Our calculations show that 16-atoms SQS model provides adequate approach to the mixing energy of the alloy. However, calculation of mechanical properties requires larger supercell. Excellent agreement between experimental and calculated elastic properties was found in fcc region of Fe-Ni system.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Prediction of the structures and heats of formation of MO 2 , MO 3 , and M 2 O 5 for M = V, Nb, Ta, Pa

Structures for the mono-, di-, and tri-bridge isomers of M 2 O 5 as well as those for the MO 2 and MO 3 fragments for M = V, Nb, Ta, and Pa were optimized at the density functional theory (DFT) level. Single point CCSD(T) calculations extrapolated to the complete basis set (CBS) limit at the DFT geometries were used to predict the energetics. Here, the lowest energy dimer isomer was the di-bridge for M = V and Nb and the tri-bridge for M = Ta and Pa. The di-bridge isomers were predicted to be composed of MO 2 + and MO 3 - fragments, whereas the mono- and tri-bridge are two MO 2 + fragments linked by an O 2 - . The heats of formation of M 2 O 5 dimers, as well as MO 2 and MO 3 neutral and ionic species were predicted using the Feller–Peterson–Dixon (FPD) approach. The heats of formation of the MF 5 species were calculated to provide additional benchmarks. Dimerization energies to form the M 2 O 5 dimers are predicted to become more negative going down group 5 and range from -29 to -45 kcal mol -1 . The ionization energies (IEs) for VO 2 and TaO 2 are essentially the same at 8.75 eV whereas the IEs for NbO 2 and PaO 2 are 8.10 and 6.25 eV, respectively. The predicted adiabatic electron affinities (AEAs) range from 3.75 eV to 4.45 eV for the MO 3 species and vertical detachment energies from 4.21 to 4.59 eV for MO 3 - . The calculated M[double bond, length as m-dash]O bond dissociation energies increase from 143 kcal mol -1 for M = V to ~170 kcal mol -1 for M = Nb and Ta to ~200 kcal mol -1 for M = Pa. The M–O bond dissociation energies are all similar ranging from 97 to 107 kcal mol -1 . Natural bond analysis provided insights into the types of chemical bonds in terms of their ionic character. Pa 2 O 5 is predicted to behave like an actinyl species dominated by the interactions of approximately linear PaO 2 + groups.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2006

Electrostatic Properties of PE and PTFE Subjected to Atmospheric Pressure Plasma Treatment; Correlation of Experimental Results with Atomistic Modeling

The use of an atmospheric pressure glow discharge (APGD) plasma was used at KSC to increase the hydrophilicity of spaceport materials to enhance their surface charge dissipation and prevent possible ESD in spaceport operations. Significant decreases in charge decay times were observed after tribocharging the materials using the standard KSC tribocharging test. The polarity and amount of charge transferred was dependent upon the effective work function differences between the respective materials. In this study, polyethylene (PE) and polytetrafluoroethylene (PTFE) were exposed to a He+O2 APGD. The pre and post treatment surface chemistry was analyzed by X-ray photoelectron spectroscopy and contact angle measurements. Semi-empirical and ab initio calculations were performed to correlate the experimental results with some plausible molecular and electronic structure features of the oxidation process. For the PE, significant surface oxidation was observed, as indicated by XPS showing C-O, C=O, and O-C=O bonding, and a decrease in the surface contact angle from 98.9 deg to 61.2 deg. For the PTFE, no C-O bonding appeared and the surface contact angle increased indicating the APGD only succeeded in cleaning the PTFE surface without affecting the surface structure. The calculations using the PM3 and DFT methods were performed on single and multiple oligomers to simulate a wide variety of oxidation scenarios. Calculated work function results suggest that regardless of oxidation mechanism, e.g. -OH, =0 or a combination thereof, the experimentally observed levels of surface oxidation are unlikely to lead to a significant change in the electronic structure of PE and that its increased hydrophilic properties are the primary reason for the observed changes in its electrostatic behavior. The calculations for PTFE argue strongly against significant oxidation of that material, as confirmed by the XPS results.

Trigwell, Steve↗
DOE OSTI2021

Regularized machine learning on molecular graph model explains systematic error in DFT enthalpies

Abstract A major goal of materials research is the discovery of novel and efficient heterogeneous catalysts for various chemical processes. In such studies, the candidate catalyst material is modeled using tens to thousands of chemical species and elementary reactions. Density Functional Theory (DFT) is widely used to calculate the thermochemistry of these species which might be surface species or gas-phase molecules. The use of an approximate exchange correlation functional in the DFT framework introduces an important source of error in such models. This is especially true in the calculation of gas phase molecules whose thermochemistry is calculated using the same planewave basis set as the rest of the surface mechanism. Unfortunately, the nature and magnitude of these errors is unknown for most practical molecules. Here, we investigate the error in the enthalpy of formation for 1676 gaseous species using two different DFT levels of theory and the ‘ground truth values’ obtained from the NIST database. We featurize molecules using graph theory. We use a regularized algorithm to discover a sparse model of the error and identify important molecular fragments that drive this error. The model is robust to rigorous statistical tests and is used to correct DFT thermochemistry, achieving more than an order of magnitude improvement.

36 MATERIALS SCIENCE↗
DOE OSTI2022

How the Bioinspired Fe 2 Mo 6 S 8 Chevrel Breaks Electrocatalytic Nitrogen Reduction Scaling Relations

The nitrogen reduction reaction (NRR) is a renewable alternative to the energy- and CO 2 -intensive Haber–Bosch NH 3 synthesis process but is severely limited by the low activity and selectivity of studied electrocatalysts. The Chevrel phase Fe 2 Mo 6 S 8 has a surface Fe–S–Mo coordination environment that mimics the nitrogenase FeMo-cofactor and was recently shown to provide state-of-the-art activity and selectivity for NRR. Here, we elucidate the previously unknown NRR mechanism on Fe 2 Mo 6 S 8 via grand-canonical density functional theory (GC-DFT) that realistically models solvated and biased surfaces. Fe sites of Fe 2 Mo 6 S 8 selectively stabilize the key *NNH intermediate via a narrow band of free-atom-like surface d-states that selectively hybridize with p-states of *NNH, which results in Fe sites breaking NRR scaling relationships. These sharp d-states arise from an Fe–S bond dissociation during N 2 adsorption that mimics the mechanism of the nitrogenase FeMo-cofactor. Furthermore, we developed a new GC-DFT-based approach for calculating transition states as a function of bias (GC-NEB) and applied it to produce a microkinetic model for NRR at Fe 2 Mo 6 S 8 that predicts high activity and selectivity, in close agreement with experiments. Furthermore, our results suggest new design principles that may identify effective NRR electrocatalysts that minimize the barriers for *N 2 protonation and *NH 3 desorption and that may be broadly applied to the rational discovery of stable, multinary electrocatalysts for other reactions where narrow bands of surface d-states can be tuned to selectively stabilize key reaction intermediates and guide selectivity toward a target product. Furthermore, our results highlight the importance of using GC-DFT and GC-NEB to accurately model electrocatalytic reactions.

10 SYNTHETIC FUELS↗
DOE OSTI2021

Experimental and calculated infrared spectra of disubstituted naphthoquinones

In recent years there has been interest in incorporating substituted 1,4-naphthoquinones (NQs) into the A 1 binding site in photosystem I (PSI) photosynthetic protein complexes. This interest in part stems from the considerably altered bioenergetics of electron transfer that occur in PSI with such substitutions. Time resolved FTIR studies of PSI complexes with disubstituted NQs incorporated have and currently are being undertaken, and with this in mind it is worth considering FTIR absorption spectra of these disubstituted NQs in solution. Here we present FTIR absorbance spectra for 2-bromo-3-methyl-1,4-naphthoquinone (BrMeNQ), 2-chloromethyl-3-methyl-1,4-naphthoquinone (CMMeNQ) and 2-ethylthio-3-methyl-1,4-naphthoquinone (ETMeNQ) in tetrahydrofuran (THF). The FTIR spectra of these di-substituted naphthoquinones (NQs) were compared to FTIR spectra of 2-methyl-3-phytyl-1,4-naphthoquinone [phylloquinone (PhQ)], 2,3-dimethyl-1,4-naphthoquinone (DMNQ), and 2-methyl-1,4-naphthoquinone (2MNQ). To aid in the assignment of bands in the experimental spectra, density functional theory (DFT) based vibrational frequency calculations for all the substituted NQs in solution were undertaken. The calculated and experimental spectra agree well. By calculating normal mode potential energy distributions, unambiguous quantitative band assignments were made. The calculated and experimental spectra together make predictions about what may be observable in time resolved FTIR difference spectra obtained using PSI with the different NQs incorporated. As a result, time resolved FTIR difference spectra are presented that support these predictions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Robustness of magnetocrystalline anisotropy and coercivity in Fe–Co–B

As a potential candidate for critical chemical element free permanent magnet, (Fe 0.7 Co 0.3 ) 2 B, shows promising intrinsic magnetic properties, including large magnetization, high Curie temperature and moderate magnetocrystalline anisotropy (MCA). Still, the coercivity that has been achieved is relatively low. To understand and improve coercivity, here we study the robustness of MCA related to the fluctuation of chemical composition and the structural distortion and its effect on coercivity in (Fe 0.7 Co 0.3 ) 2 B from first principles DFT (density functional theory) calculation and micromagnetic simulation. A fluctuation of cobalt content of 10%, or 1% isotropic lattice strain, or 2% tetragonal distortion reduces MCA energy up to 30%. By considering randomly distributed local MCA energy reduction up to 30%, the calculated coercivity of anisotropic bulk magnets reduced from 3.6 kOe to 2.2 kOe. Rational microstructure design, such as grain size refinement and development of strong crystal texture of MCA easy-axis (i.e., degree of alignment for (Fe 0.7 Co 0.3 ) 2 B grains), can enhance coercivity and mitigate the effect of local MCA reduction. In addition to acceptable intrinsic magnetic properties, the promising hard magnetic phase needs good MCA robustness to facilitate the development of practical magnets.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗