Search NASA⌕ Search

SEARCH · Search NASA

Results for “particle distribution functions”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Surface enhanced performance of La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3-δ cathodes by infiltration Pr-Ni-Mn-O progress

Herein the present study reports the enhancement of electrochemical oxygen reduction activity of porous La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3–δ (LSCF) cathodes by coating a thin film of Pr-Ni-Mn oxide (PNM5) using a multi-step infiltration process. XRD examination reveals that PNM5 mainly contains a multiphase mixture of Pr 6 O 11 , PrNiO 3 , MnO and NiO. SEM morphology shows a thin PNM5 film and small particles are formed on the surface of LSCF backbone particles. Impedance spectrum analysis indicates that PNM5 infiltrated LSCF exhibits dramatically reduced polarization resistance (R p ), reaching R p of 0.244 Ω cm 2 at 973 K, which is one-half of the baseline LSCF cathode. The activation energy of LSCF cathodes infiltrated with PNM5 is 1.45 eV, slightly lower than the baseline LSCF cathode (1.77 eV). Distribution of relaxation time (DRT) function analysis shows PNM5 infiltration layer significantly promotes the oxygen reduction reaction (ORR) of cathode surface. With the increase in firing temperature, the total resistance increases and R p changes from ion transport to oxygen reduction reaction. Degradation rate of the PNM5-infiltrated LSCF is also lower, 0.02168% vs 0.07093% for the baseline LSCF over a 200-h period. A single cell testing indicated that the peak power density of the PNM5-infiltrated cell was increased by 140.66%. Overall, PNM5 could be a potential catalyst for boosting the performance of a commercial LSCF cathode for solid oxide fuel cells (SOFC).

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nuclear spin features relevant to ab initio nucleon-nucleus elastic scattering

Effective interactions for elastic nucleon-nucleus scattering from first principles require the use of the same nucleon-nucleon interaction in the structure and reaction calculations, as well as a consistent treatment of the relevant operators at each order. Previous work using these interactions has shown good agreement with available data. Here, we study the physical relevance of one of these operators, which involves the spin of the struck nucleon, and examine the interpretation of this quantity in a nuclear structure context. Using the framework of the spectator expansion and the underlying framework of the no-core shell model, we calculate and examine spin-projected, one-body momentum distributions required for effective nucleon-nucleus interactions in $J=0$ nuclear states. The calculated spin-projected, one-body momentum distributions for $^4$He, $^6$He, and $^8$He display characteristic behavior based on the occupation of protons and neutrons in single particle levels, with more nucleons of one type yielding momentum distributions with larger values. Additionally, we find this quantity is strongly correlated to the magnetic moment of the $2^+$ excited state in the ground state rotational band for each nucleus considered. In conclusion, we find that spin-projected, one-body momentum distributions can probe the spin content of a $J=0$ wave function. This feature may allow future ab initio nucleon-nucleus scattering studies to inform spin properties of the underlying nucleon-nucleon interactions. The observed correlation to the magnetic moment of excited states illustrates a previously unknown connection between reaction observables such as the analyzing power and structure observables like the magnetic moment.

6 ≤ A ≤ 19↗
DOE OSTI2026

Microscopic theory of angular momentum distributions across the full range of fission fragments

Modern nuclear theory provides qualitative insights into the fundamental mechanisms of nuclear fission and is increasingly capable of making reliable quantitative predictions. Most quantities of interest pertain to the primary fission fragments, whose subsequent decay is typically modeled using statistical reaction models. Consequently, a key objective of fission theory is to inform these models by predicting the initial conditions of the primary fragments. Here, in this work, we employ a framework that combines joint angular momentum and particle number projection with time-dependent configuration mixing to calculate the angular momentum distributions of primary fragments. Focusing on the benchmark cases of neutron-induced fission of 235 U and 239 Pu , we predict—for the first time—microscopic angular momentum distributions for all fragments observed in experiments. Our results reveal a pronounced sawtooth pattern in the average angular momentum as a function of fragment mass, consistent with recent measurements. Additionally, we observe substantial variations in angular momentum distributions along isobaric chains, indicating that commonly used empirical formulas lack sufficient accuracy. We also quantify a strong correlation between the angular momentum and the deformation of the fragments at scission, and a weak correlation in the magnitude of the angular momentum between fragment partners. The generated data will enable estimation of the impact of microscopic distributions on fission spectra, paving the way toward fission modeling based on microscopic inputs.

Physics - Nuclear physics and radiation physics↗
DOE OSTI2025

HfZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-zirconium (Hf-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

DED Additively Manufactured HEAs Optimized via Parametric Study of Functionally Graded Materials

Complex concentrated alloys (CCAs) are a system of alloys composition of nearly equiatomic elements, creating a state of high configurational entropy in the material. These alloys are of interest for various industrial applications due to their high strength and hardness, wide operational temperature range, creep and diffusion resistance, and radiation resistance. In the current work CoCrFeNi-base CCAs are synthesize by combining three commonly industrial alloys: IN718, SS316L, and 70Co30Cr (commercially known as Stellite 21), in-situ utilizing directed energy deposition (DED) additive manufacturing (AM). One bulk specimen is fabricated with a 1:1:1 ratio of IN718, SS316L, and 70Co30Cr where IN718 is a source of nickel and chrome, SS316 is a source of iron and chrome, and 70Co30Cr is a source of cobalt and supplemental chrome. Additionally, a functionally graded material (FGM) of the identified CCA system from SS316 was fabricated with 10 layers, where the initial five layers are the base material of SS316L and graded by steps of 20% to transition to the CCA system. Overall build quality, microstructure, compositional distribution, phases analysis, and microhardness were studied by laser optical microscopy, scanning election microscopy (SEM), energy dispersive spectroscopy (EDS), X-ray diffraction analysis (XRD), and Vickers microhardness techniques. Overall build quality showed a successful bulk build with acceptable density, however unmelted particles are observed in the as-built alloy. Compositional results show that the system meets the Boltzmann’s hypothesis definition of an HEA, and that an FCC CCA was produced with increased hardness, small grain size, and superior hardness properties to the SS316L substrate. While initial XRD results show predominately FCC structure, EDS analysis of interdendritic regions show the existence of Laves phases in the material. The fabrication of an FGM CCA allowed for potential screening of additional alloy mixes which may be of interest to explore in future research, and the resulting composition is compared to predicted values. This novel method of CCA fabrication results in substantial cost savings of these CCA systems over traditional methods.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Field, frequency, and temperature dependencies of the surface resistance of nitrogen diffused niobium superconducting radio frequency cavities

We investigate the rf performance of several single-cell superconducting radio-frequency cavities subjected to low temperature heat treatment in nitrogen environment. The cavities were treated at temperature 120 – 165 ° C for an extended period of time (24–48 h) either in high vacuum or in a low partial pressure of ultrapure nitrogen. The improvement in Q 0 with a Q rise was observed when nitrogen gas was injected at ∼ 300 ° C during the cavity cooldown from 800 ° C and held at 165 ° C , without any degradation in accelerating gradient over the baseline performance. The treatment was applied to several elliptical cavities with frequency ranging from 0.75 to 3.0 GHz, showing an improved quality factor as a result of low temperature nitrogen treatments. The Q rise feature is similar to that achieved by nitrogen alloying Nb cavities at higher temperature, followed by material removal by electropolishing. The surface modification was confirmed by the change in electronic mean free path and tuned with the temperature and duration of heat treatment. The decrease of the temperature-dependent surface resistance with increasing rf field, resulting in a Q rise, becomes stronger with increasing frequency and decreasing temperature. The data suggest a crossover frequency of ∼ 0.95 GHz above that the Q rise phenomenon occurs at 2 K. Some of these results can be explained qualitatively with an existing model of intrinsic field-dependence of the surface resistance with both equilibrium and nonequilibrium quasiparticle distribution functions. The change in the Q slope below 0.95 GHz may result from masking contribution of trapped magnetic flux to the residual surface resistance. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Structural Explanation of the Dielectric Enhancement of Barium Titanate Nanoparticles Grown under Hydrothermal Conditions

Abstract When synthesized under certain conditions, barium titanate (BaTiO 3 , BTO) nanoparticles are found to have the non‐thermodynamic cubic structure at room temperature. These particles also have a several‐fold enhanced dielectric constant, sometimes exceeding 6000, and are widely used in thin‐layer capacitors. A hydrothermal approach is used to synthesize BTO nanocrystals, which are characterized by a range of methods, including X‐ray Rietveld refinement and the Williamson–Hall approach, revealing the presence of significant inhomogeneous strain associated with the cubic phase. However, X‐ray pair distribution function measurements clearly show the local structure is lower symmetry than cubic. This apparent inconsistency is resolved by examining 3D Bragg coherent diffraction images of selected nanocrystals, which show the existence of ≈50 nm‐sized domains, which are interpreted as tetragonal twins, and yet cause the average crystalline structure to appear cubic. The ability of these twin boundaries to migrate under the influence of electric fields explains the dielectric anomaly for the nanocrystalline phase.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Development and Integration of a Stochastic Clad Damage Propagation Model into PRONGHORN-SC Subchannel Analysis Code

The failure of fuel pins in nuclear reactors is intrinsically stochastic. Typically, a combination of variation in manufacturing that affects the material characteristics and the fuel assembly dimensions, variation in operating conditions, such as local power, coolant flow rate, and irradiation induced changes in material properties lead to a large uncertainty in failure margin of the fuel pins. Failure, therefore, may occur in exceptional pins with adverse combinations of these variations. Upon a metal fuel pin (U-Pu-Zr/HT9) failure, depressurization of the fuel pin takes place by release of fission gas, liquid sodium bond, and potentially solid fuel particles or molten/eutectic fuel droplets through the hole in cladding. The effect of a fission gas jet on neighbor fuel pins and possible propagation of a clad damage during normal operation was studied experimentally in 1970s and it was found that the post-failure fission gas jet insulates the jet impingement area of the target fuel pin surface and could increase the target pin’s surface temperature by as much as 100 – 200 K during the failed pin depressurization. It was concluded that the effect should not lead to fuel pin failure propagation during normal operation. In accident scenarios of sodium and lead fast reactors such as Unprotected Loss-Of-Flow (ULOF) or Unprotected Transient Over Power (UTOP), the fuel pins can be subjected to higher clad temperatures and fuel pin pressures or fuel clad mechanical/chemical interaction where thermal creep margin becomes significantly lower compared to the normal operation conditions. Therefore, possible stochastic failure and the post-failure fission gas/fuel jet impingement could be critical in order to predict fuel pin failure propagation. Pin depressurization due to fission gas release may degrade the heat transfer by formation of a gas blanket on a neighboring pin surface, which is a local phenomenon, and by causing coolant flow deceleration and starvation, which could affect a surrounding region as well. Furthermore, the potential presence of solid fuel particles or molten fuel at the time of clad failure could boost post-failure jet induced degradation even further. The present study models the U-Pu-Zr/HT9 metal fuel pin failure and stochastic clad damage propagation by biased sampling based on a Cumulative Damage Fraction (CDF) type clad failure criterion and the normal distribution of fuel failure probability density as a function of logarithm of Cumulative Damage Fraction. In addition, the effect of post-failure fission gas jet on heat transfer degradation is modeled for the target pins. This model is called stochastic Clad Damage Propagation (CDAP). The CDAP model is now fully integrated into developmental version of PRONGHORN-SC subchannel analysis code, allowing for modeling local failures and its propagation potential. Section 2 describes the components of the CDAP models. Section 3 describes the model implementation to PRONGHORN-SC and input specifications. Section 4 describes the CDAP model validation coupled to PRONGHORN-SC.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

HfTa_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-tantalum (Hf-Ta). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Ta. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect BCC lattice sites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,040 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Structural properties of PbTe quantum dots revealed by high-energy x-ray diffraction

High-energy x-ray diffraction (HE-XRD) experiments combined with an analysis based on atomic-pair-distribution functions can be an effective tool for probing low-dimensional materials. Here, we show how such an analysis can be used to gain insight into structural properties of PbTe nanoparticles (NPs). Here we interpret our HE-XRD data using an orthorhombic Pnma phase of PbTe, which is an orthorhombic distortion of the rocksalt phase. Although local crystal geometry can vary substantially with particle size at scales below 10 nm, and for very small NPs the particle size itself influences x-ray diffraction patterns, our study shows that HE-XRD can provide a unique nano-characterization tool for unraveling structural properties of nanoscale systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Revisiting the Hail Radar Reflectivity–Kinetic Energy Flux Relation by Combining T-Matrix and Discrete Dipole Approximation Calculations to Size Distribution Observations

The retrieval of hail kinetic energy with weather radars or its simulation in numerical models is challenging because of the shape complexity and variable density of hailstones. We combine 3D scans of individual hailstones with measurements of the particle size distributions (PSD) and T-matrix calculations to understand how hail reflectivity Z changes when approximating hailstones as spheroids, as compared to the realistic shapes obtained by 3D scanning technology. Additionally, recent terminal velocity relations are used to compare Z to the hail kinetic energy flux E ˙ . We parameterize the hail backscattering cross sections at L, S, C, and X bands as a function of size between 0.5 and 5.0 cm, matching the range of the observed PSDs. The scattering calculations use the T-matrix method for size parameters below 1.0 and the discrete dipole approximation (DDA) method otherwise. The DDA calculations are done for 48 digital models of realistic hailstones of sizes between 1 and 5 cm. The DDA cross sections are calculated for multiple orientations and averaged assuming a fully random orientation distribution to provide a single value per hailstone. The T-matrix reflectivity assuming solid ice spheres presents negligible differences to DDA results for size parameters below 1.0. Therefore, T matrix was used to fill in the gaps left by the DDA calculations. The results are mapped to the same size bins of the observed PSDs, allowing the calculation of the radar reflectivity. This is then correlated to E ˙ , allowing a potential improvement of past retrieval methods of E ˙ from Z in multiple wavelengths.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Nucleon momentum distributions for local chiral interactions [Slides]

We use quantum Monte Carlo methods to calculate single- and two-nucleon momentum distributions in 4 He, 12 C, and 16 O. We use correlated many-body wave functions optimized for local chiral interactions up to next-to-next-to-leading order (N 2 L0).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Extreme Heating of Minor Ions in Imbalanced Solar-wind Turbulence

Minor ions in the solar corona are heated to extreme temperatures, far in excess of those of the electrons and protons that comprise the bulk of the plasma. These highly nonthermal distributions make minor ions sensitive probes of the collisionless processes that heat the corona and power the solar wind. The recent discovery of the "helicity barrier" offers a mechanism in which imbalanced Alfvénic turbulence in low-β plasmas preferentially heats protons over electrons, generating high-frequency, proton-cyclotron-resonant fluctuations. We use the hybrid-kinetic particle-in-cell code Pegasus++ to drive imbalanced Alfvénic turbulence in a 3D low-β plasma with additional passive ion species, He 2+ and O 5+ . A helicity barrier naturally develops, followed by clear phase-space signatures of oblique proton-cyclotron-wave heating and Landau-resonant heating from the imbalanced Alfvénic fluctuations. The former results in characteristically arced ion velocity distribution functions, whose non-bi-Maxwellian features are shown by linear ALPS calculations to be critical to the heating process. Additional features include a steep transition-range electromagnetic spectrum, proton-cyclotron waves propagating in the direction of the imbalance, significantly enhanced proton-to-electron heating ratios, ion temperatures that are considerably more perpendicular with respect to magnetic field, and extreme heating of heavier species in a manner consistent with mass scalings inferred from spacecraft measurements. None of these features are realized in an otherwise equivalent simulation of balanced turbulence. If seen simultaneously in the fast solar wind, these signatures of the helicity barrier would testify to the necessity of incorporating turbulence imbalance in a complete theory for the evolution of the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Enhanced collisionless laser absorption in strongly magnetized plasmas

We report that strongly magnetizing a plasma adds a range of waves that do not exist in unmagnetized plasmas and enlarges the laser-plasma interaction (LPI) landscape. In this paper, we use particle-in-cell simulations to investigate strongly magnetized LPI in one dimension under conditions relevant for magneto-inertial fusion experiments, focusing on a regime where the electron-cyclotron frequency is greater than the plasma frequency and the magnetic field is at an oblique angle with respect to the wave vectors. We show that when electron-cyclotron-like hybrid wave frequency is about half the laser frequency, the laser light resonantly decays to magnetized plasma waves via primary and secondary instabilities with large growth rates. These distinct magnetic-field-controlled instabilities, which we collectively call two-magnon decays, are analogous to two-plasmon decays in unmagnetized plasmas. Since additional phase mixing mechanisms are introduced by the oblique magnetic field, collisionless damping of large-amplitude magnetized waves substantially broadens the electron distribution function, especially along the direction of the magnetic field. During this process, energy is transferred efficiently from the laser to plasma waves and then to electrons, leading to a large overall absorptivity when strong resonances are present. The enhanced laser energy absorption may explain hotter-than-expected temperatures observed in magnetized laser implosion experiments and may also be exploited to develop more efficient laser-driven x-ray sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Investigation of 3D printed lightweight hybrid composites via theoretical modeling and machine learning

Hybrid composites combine two or more different fillers to achieve multifunctional or advanced material properties, such as lightweight and enhanced mechanical properties. The properties of the composites significantly depend on their microstructures, which can be tailored via advanced 3D printing processes. Understanding the process-structure-property relationships is critical to enable the design and engineering of novel hybrid composites for applications in aerospace, automotive, and protective coatings. Here, for this work, we develop 3D printable and lightweight hybrid composites and leverage the conventional design of experiments, a theoretical hybrid model, and an image-driven machine learning (ML) method to investigate their mechanical behaviors. The hybrid composites are formulated with elastomer matrix, microfillers, and thin-shell particles, enabling a significant degree of design freedom of microstructures with densities and mechanical properties varying up to 70% and 91%, respectively. Our statistical analysis indicates that the 3D printing path direction and the microfibers fraction are dominating process parameters with contribution percentages of 45.3% and 57.7% on the specific stiffness and strength, respectively. A hybrid mechanics model is developed based on a simple Weibull distribution function and classical single-filler models to effectively capture the variations in mechanical properties, however, it overestimates the values due to its statistical constraints and idealization of experimental uncertainty. The image-driven ML model leverages the microscale images directly without losing the structural details, shows more accurate predictions with experimental data, and has 48.6% lower root mean square error than the theoretical model.

3D printing↗
DOE OSTI2026

HfTi_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-titanium (Hf-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Ti. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2026

TiV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys titanium-vanadium (Ti-V). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Ti and V. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗