Vienna Ab Initio Simulation Package Research Articles

Excavating oxygen vacancies in BaZrO3 to boost photocatalysis via screening vantage crystal planes

Investigations into the design and advancement of oxygen vacancies (OVs) are at the forefront of high-efficiency catalyst research. However, the lack of experimental evidence and a comprehensive mechanistic understanding regarding the promotion of OVs poses challenges. Herein, we investigate the role of crystal planes in elucidating the mechanism of OVs formation and their impact on photocatalytic hydrogen evolution. Initially, the Vienna ab-initio simulation package was utilized to calculate the key characteristics and band parameters of BaZrO3, subsequently, the screening crystal planes of (110) was synthesized. Experimental findings indicated that the (110) crystal plane increased the OVs rate from 32 % to 56 %, consequently improving photocatalytic hydrogen evolution performance from 13.1 to 20.5 mmol·g−1·h−1. Furthermore, theoretical calculations were employed to delve into both OVs formation mechanisms and catalytic activity. This study presents an innovative approach for strategically modifying catalyst OVs for efficient and high-performance photocatalytic applications.

Investigation of Be, Mg, Ti-adsorbed boron-germanene nanoribbons for nano applications.

One-dimensional systems are nanostructures of significant interest in research due to their numerous potential applications. This study focuses on the investigation of one-dimensional boron-germanene nanoribbons (BGeNRs) and BGeNRs doped with Be, Mg, and Ti. Density functional theory (DFT) combined with the Vienna Ab initio Simulation Package (VASP) forms the foundation of this research. The electromagnetic and optical properties of these structures are systematically examined. The findings reveal that all the studied structures exhibit metallic behavior, with differences in their magnetic properties. The magnetic moments of the pristine and Be-doped structures are both zero, whereas the Mg and Ti-doped structures exhibit magnetic moments of 0.012 μB and 2.234 μB, respectively. Partial density of states (PDOS) analyses highlight the contributions of various elements and the complex multi-orbital hybridization among them. The optical properties are investigated through the real and imaginary parts of the dielectric function, along with the absorption coefficient and electron-hole density. This study indicates potential applications in adsorption sensors, the modulation of system magnetism via adsorption, and information transmission technologies.&#xD.

First-principles data for solid solution niobium-tantalum-vanadium alloys with body-centered-cubic structures

We present four open-source datasets that provide results of density functional theory (DFT) calculations of ground-state properties of refractory solid solution binary alloys niobium-tantalum (NbTa), niobium-vanadium (NbV), tantalum-vanadium (TaV), and ternary alloys NbTaV ordered in body-centered-cubic (BCC) structures with 128 Bravais lattice sites. The first-principles code used to run the calculations is the Vienna Ab-Initio Simulation Package. The calculations have been collected by uniformly sampling chemical compositions across the entire compositional range. For each chemical composition, the calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. This sampling methodology resulted in running DFT simulations for a total of 3,100 randomized atomic configurations over 31 chemical compositions for each of the three binary alloys Nb-Ta, Nb-V, Ta-V, and a total of 10,500 randomized atomic structures over 105 chemical compositions for the ternary alloys Nb-Ta-V. For each atomic configuration, geometry optimization has been performed, and the data released contains information about each step of geometry optimization for each atomic configuration.

Mechanistic investigation of methanol-to-olefins conversion catalyzed by H-ZSM-5 zeolite: a DFT study.

The mechanisms for the formation of the first C - C bond and lower olefins on methanol to olefins (MTO) conversion on H-ZSM-5 had been focused in dispute. In this paper, density functional theory has been used to study the reaction mechanisms of methanol to olefins on ZSM-5. The configurations of reactants, intermediates, products and transition state of the numerous reactions involved in such a process have been optimized, as well as the elementary reactions related to these configurations were determined by the calculation of corresponding activation energy barriers and reaction heats. Here, two different kinds of the mechanisms were proposed for the formation of dimethylether (DME), one involving an associative interaction of two methanol molecules with the zeolite Brønsted acid sites and the other occurring via a surface methoxy species and a methanol molecule. A critical intermediate of the methoxy methyl cation was theoretically verified by the reaction of the methoxy species and dimethylether. Besides, it was found that the first intermediates containing a C - C bond were 1,2-dimethoxyethane and 2-methoxy-ethanolare, in which the former was formed from methoxy species with dimethyl ether and the latter was formed from methanol by onium ions((CH3)2O+CH2CH2OCH3), respectively. For the whole reaction mechanism, the results in this paper indicated that the ethene formation is more favorable than propylene formation due to the low activation energy barrier for ethene formation (123.49 vs. 162.09kJ.mol-1). From these calculations, it would be concluded that ethene is the first alkene product that induces the occurrence of the hydrocarbon pool mechanism. All the periodic density function theory (DFT) calculations were performed by the Vienna Ab Initio Simulation package (VASP). The interaction between nucleus and valence electron was described using the pseudopotentials found in the projector augmented wave (PAW) method. PBE-D3 was used in the whole DFT calculations and CI-NEB was used to locate transition state.

Ab Initio Modelling of g-ZnO Deposition on the Si (111) Surface

Recent studies show that zinc oxide (ZnO) nanostructures have promising potential as an absorbing material. In order to improve the optoelectronic properties of the initial system, this paper considers the process of adsorbing multilayer graphene-like ZnO onto a Si (111) surface. The density of electron states for two- and three-layer graphene-like zinc oxide on the Si (111) surface was obtained using the Vienna ab-initio simulation package by the DFT method. A computer model of graphene-like Zinc oxide on a Si (111)-surface was created using the DFT+U approach. One-, two- and three-plane-thick graphene-zinc oxide were deposited on the substrate. An isolated cluster of Zn3O3 was also considered. The compatibility of g-ZnO with the S (100) substrate was tested, and the energetics of deposition were calculated. This study demonstrates that, regardless of the possible configuration of the adsorbing layers, the Si/ZnO structure remains stable at the interface. Calculations indicate that, in combination with lower formation energies, wurtzite-type structures turn out to be more stable and, compared to sphalerite-type structures, wurtzite-type structures form longer interlayers and shorter interplanar distances. It has been shown that during the deposition of the third layer, the growth of a wurtzite-type structure becomes exothermic. Thus, these findings suggest a predictable relationship between the application method and the number of layers, implying that the synthesis process can be modified. Consequently, we believe that such interfaces can be obtained through experimental synthesis.

Aptamer Encapsulated Inside the Array Channel of Ni‐MOF for Bisphenol A Determination in Multi‐interference System

AbstractAs an endocrine disruptor, bisphenol A (BPA) has many adverse effects on environmental safety and human health. Aiming at the existing problems at the current detection methods, herein, a 3D metal‐organic framework (Ni‐MOF) was fabricated as a carrier to achieve BPA aptamers immobilization inside its array channels under Mg2+ regulation. Based on the reversible quenching of 6‐FAM (6‐carboxyfluorescein) fluorescence labeled on the BPA aptamers during the encapsulation process, a fluorescent aptamer sensing platform was constructed for the quantitative detection of BPA in water. Morphology analysis, VASP (Vienna ab‐initio simulation package) and other kinetic simulations were performed to elaborate the mechanism of the influence of pore size, medium ions, etc. on the recognition process. Due to the in‐hole fixation strategy of aptamers, Ni‐MOF played an obvious protection aptamers. Not only the activity of BPA aptamer@MOF composite was maintained more than 50 % in complex environments such as pH 3.0–11.0, 30–70 °C and organic solvents, but also the aptamer was protected from nuclease hydrolysis under physiological conditions. The stability and application range of the sensor are greatly improved, detection limit of 0.34 μM. The result was expected to provide theoretical guidance for the rapid detection of pollutants in complex environments.

Orange-reddish emitting of trivalent samarium doped double perovskite Sr2LuNbO6 with high thermal stability for white-LED application

A series of trivalent samarium ions (Sm3+) doped double perovskite Sr2LuNbO6 (x at% Sm3+:Sr2LuNbO6, x = 0.5–15) phosphors have been prepared via high-temperature solid-state method. The structural characterization of these phosphors was conducted employing X-ray diffraction (XRD), XRD Rietveld refinement, scanning electron microscope (SEM) and transmission electron microscope (TEM), respectively. Density functional theory calculations utilizing the Vienna ab initio simulation package (VASP) were employed to unveil the electronic structures, encompassing the band structure and density of states (DOS). Photoluminescence (PL) spectroscopy was employed to scrutinize the luminescent properties, while the thermal stability of the PL spectra was also evaluated. The most prominent emission peak, exhibiting an orange-reddish hue, was observed at 601 nm upon excitation at 406 nm. The quenching mechanism of Sm3+ ions within the Sr2LuNbO6 host lattice was identified as a dipole-dipole interaction, with the quenching concentration of 3 at%. Moreover, the phosphors demonstrated excellent thermal stability, with an activation energy (Ea) value of 0.126 eV. These findings collectively underscore the promising potential of Sm3+:Sr2LuNbO6 orange-reddish phosphors for applications in white light-emitting diodes (w-LEDs).

A theoretical study of surface lithium effects on the [111] SiC nanowires as anode materials.

Silicon carbide nanowires (SiCNWs) are considered a promising alternative material for application in lithium-ion batteries, with researchers striving to develop new electrode materials that exhibit high capacity and high charge/discharge rate performance. To gain a deeper understanding of the application of SiCNWs in semiconductor material science and energy supply fields, we investigated the effects of nanoscale and surface lithiation on the electrical and mechanical properties of SiCNWs grown along the [111] direction. First-principles calculation was used to study their geometries, electronic structures, and associated electrochemical properties. Herein, we considered SiCNWs with full hydrogen passivation, full lithium passivation, and mixed passivation at different sizes. The formation energy indicates that the stability of SiCNWs increases with the increasing diameter, and the surface-lithiated SiC nanowires (Li-SiCNWs) are found to be energetically stable. The mixed passivated SiCNWs exhibit the properties of indirect band gap with the increase of lithium atoms on the surface, while the fully lithium passivated nanowires exhibit metallic behavior. Charge analysis shows that a portion of the electrons on the lithium atoms are transferred to the surface atoms of the nanowires and electrons prefer to cluster more near the C atoms. Additionally, Li-SiCNWs still have good mechanical resistance during the lithiation process. The stable open-circuit voltage range and theoretical capacity of these SiCNWs indicate their suitability as anode materials. In this study, Materials Studio 8.0 was used to construct the models of the SiCNWs. And all the density functional theory (DFT) calculations were performed by the Vienna ab initio Simulation Package (VASP). The self-consistent field calculations are performed over a Monkhorst-Pack net of 1 × 1 × 6k-points. The energy convergence criteria for the self-consistent field calculation were set to 10-5eV/atom with a cutoff energy of 400eV.

Initial corrosion behavior of 10CrNi3MoV steel and Q235 steel in marine atmosphere

Exploring the corrosion behavior of materials in a marine atmosphere and studying their corrosion resistance mechanisms can provide a basis for developing new weathering steels. This paper primarily introduces the initial corrosion behavior and corrosion resistance mechanism of 10CrNi3MoV steel in a marine atmosphere and explores the effects of Ni on the corrosion resistance. The corrosion resistance of 10CrNi3MoV steel and Q235 steel was compared using the weight loss method. The results show that the corrosion rate of 10CrNi3MoV steel was approximately half that of Q235 steel. The rust layer morphology and elemental distribution of 10CrNi3MoV steel were observed via scanning electron microscopy and energy-dispersive X-ray spectroscopy. The composition of the rust layer was analyzed using X-ray diffraction and X-ray photoelectron spectroscopy, and the protective effect of the rust layer was further tested using an electrochemical test. The reason why Fe2NiO4 is beneficial for improving the corrosion resistance of the rust layer was explained using the Vienna Ab initio Simulation Package. The results reveal that the morphology, composition, element distribution, and electrochemical properties of 10CrNi3MoV steel were considerably stronger than those of Q235 steel because of the formation of Fe2NiO4 in the rust layer. Fe2NiO4 refines the rust layer particles and promotes the transformation of corrosion products to α-FeOOH. In addition, it imparts ion selectivity to the rust layer, effectively blocking the invasion of Cl−.

A DFT study of the physical properties of novel vanadium-based half-Heusler alloys VYAl (Y = Ni,Pd, and Pt) pertinent to their optoelectronics application

We investigated three novel 18-valence electron count (VEC) vanadium-based half-Heusler alloys, VYAl (Y = Ni, Pd, and Pt), for their structural, mechanical, electronic, and optical properties using density functional theory (DFT). The computation was carried out employing Vienna Ab-initio Simulation Package (VASP) and all three alloys have exhibited structural, chemical, mechanical, and thermodynamic stability. We obtained the lattice constants of 5.520, 5.800, and 5.820 Å for VNiAl, VPdAl, and VPtAl, respectively. The elastic parameters imply a ductile and hard nature of the system and exhibit anisotropic behavior. Electronic band characterization demonstrated that all three alloys are semiconductors with direct narrow gaps at the Γ-point which is crucial for photovoltaic applications. In addition to that, the high density of the band states near the Fermi level indicates a promising candidate for thermoelectric conversion. The influence of SOC interaction in the system generates a small shift in the bands. From the optical response, we infer that the studied compounds have exceptional light absorption and reflection quality, making them suitable for optoelectronics purposes.

Role of outer shell electron-nuclear distant of transition metal atoms (TMA) on band gap reduction and optical properties of TiO2 semiconductor

The reduction of the band gap of titanium dioxide (TiO2) via transition metal atoms (TMAs) is the focus of many research groups to increase its absorption to visible region of electromagnetic (EM) radiation. In this modeling it was shown that atoms having near outer electron shells affect strongly on the band gap of TiO2. The TiO2 has exceptional photoactivity, high stability, and cost-effectiveness. The results of the current study emphasized that adjustable TiO2 photocatalyst material can be produced through band gap reduction by Pt, Pd, and Ni dopants. The electronic configurations and optical properties of both undoped and (Ni, Pd and Pt)-doped TiO2 have been studied using first-principles calculations within the framework of Density Functional Theory (DFT) with the aid of the Vienna Ab initio Simulation Package (VASP). The role of electronegativity of TMAs on band gap drop was discussed for the first time. The band gap of clean TiO2 was found to be 3.2 eV using the band structure (BS) and Tauq model. The results indicate that Ni TMA with lowest electronegativity reduces the band gap of TiO2 to the lowest value (0.86 eV). The results of BS and density of state revealed the fact that small outer shell electron-nuclear distant TMA more affected the band gap of TiO2. The band gap determined from the imaginary part of the optical dielectric function closely aligns with those determined from the BS diagram. The increase in the refractive index (n) was observed with the increase of Ni, Pd, and Pt into the TiO2 structure, which indicates an increase in charge density within the TiO2 structure. The reflectance, absorbance, and transmittance results suggest that Pd and Pt-doped TiO2 are responsive to the visible region of electromagnetic radiation, whereas Ni-doped TiO2 exhibits responsiveness in the IR region. In its pure state; TiO2 was found almost transparent in the visible and IR regions of EM radiation.

Tunable electronic and magnetic properties of defective g-ZnO monolayer doping with non-metallic elements.

The electronic and magnetic properties of non-metallic (NM) elements doping defective graphene-like ZnO (g-ZnO) monolayer including O vacancy (VO) and Zn vacancy (VZn) are studied. The results show that VO-g-ZnO is a semiconductor and VZn-g-ZnO is a magnetic semiconductor. B, C, N, Si, P, 2S, and 2Si doping VO-g-ZnO systems present half-metal and magnetic semiconductors, and the magnetism mainly originates from the spin polarization of doping atoms. For single or double NM elements doping VZn-g-ZnO, 2P doping system presents a semiconductor, while other systems present ferromagnetic metal, half-metal, and magnetic semiconductor. The magnetism of single NM elements doping VZn-g-ZnO mainly comes from the spin polarization of O atoms near the defect point. For double NM elements doping VZn-g-ZnO, spin splitting occurs mainly in p orbitals of O atoms, dopant atoms, and d orbitals of Zn atoms. NM elements doping defect g-ZnO can effectively regulate the electronic and magnetic properties of the system. The software package VASP 5.4.1 (Vienna ab initio Simulation Package) is used for calculations in this paper. The local density approximation (LDA) is adopted as an exchange and correlation function to perform the structural optimization and analysis of electronic structure and magnetic properties.

Pressure effects on the structural, electronic, elastic, optical, and vibrational properties of YMg intermetallic compounds: a first-principles study

The properties of YMg in B2 structure have been comprehensively analysed using the first-principles plane-wave pseudopotential method. Specifically, the structural, electronic, elastic, vibrational, and optical properties were investigated using the generalized gradient approximation (GGA) method in the context of density functional theory. The Vienna ab initio simulation package (VASP) was utilized for these calculations. The computed lattice parameter (3.803 Å) and bulk modulus (41.33 GPa) are consistent with the earlier data on ambient pressure. The electronic band structure and energy-dependent density of states reveal the metallic nature of the titled compounds. The Born stability requirements confirmed the mechanical stability. The analysis of Pugh’s and Poisson’s ratios and Cauchy’s pressure reveals that YMg is ductile under the pressures in consideration. According to several anisotropy indices, the compound is noticeably anisotropic both in ambient and under pressure. Our investigation includes an analysis of several fundamental mechanical parameters of the material, including the bulk modulus, the pressure derivative of the Zener anisotropy factor, Poisson’s ratio, isotropic shear modulus and Young’s modulus with a particular focus on their dependence on pressure. We have determined that the elastic constants obtained remain mechanically stable, satisfying the Born Stability conditions even at high pressures of up to 60 GPa. To explore the dynamic stability of YMg, we analysed the material’s phonon dispersion curves. The examined compound displays stability under dynamic conditions from 0 GPa to 30 GPa, as evidenced by its positive vibration frequencies. However, this stability is not sustained under higher pressure, as the compound becomes unstable after 30 GPa to 70 GPa. The electronic band structure and density of states diagrams demonstrate YMg’s metallic properties. At atmospheric pressure (0 GPa), the total density of states (TDOS) near the Fermi level is approximately 1.63 states/eV, with pressure application reducing DOS. The dielectric function, refractive index, and energy loss spectra are examined within the 0–20 eV energy range. YMg exhibits its highest absorption between 4 and 11 eV. The peak optical conductivity is observed around 0.78 eV (equivalent to 1589.5409 nm), while the most significant energy loss occurs at 11.90 eV, roughly corresponding to 2.8 Hz in the ultraviolet spectrum. Moreover, we extensively analyzed the material’s phonon thermodynamic and optical properties, providing insights into its behavior under various conditions. The outcomes acquired at zero pressure are generally coherent with the current theoretical values.

Establishing the NiMo6Se8 Chevrel Phase as a Promising Material Using DFT

AbstractIn this study, the NiMo6Se8 Chevrel phase is analyzed using Density Functional Theory (DFT) and the Vienna Ab‐initio Simulation Package (VASP). The analysis focuses on the phase's structural, electrical, and mechanical characteristics to fill gaps in the current literature. The presence of a rhombohedral crystal structure confirms its thermodynamic stability, as indicated by a negative formation enthalpy, which suggests that it can be synthesized under favorable conditions. The electronic properties of the phase are analyzed, indicating that it exhibits semiconductor characteristics with a bandgap of 1.07 eV. This makes it appropriate for various technological applications. The estimated elastic constants provide an indication of mechanical strength and flexibility, with a noticeable presence of anisotropic elasticity. The confirmation of dynamical stability is achieved by analyzing the phonon dispersion curve, which reveals the absence of any negative frequencies. Furthermore, the material has a low thermal conductivity, increasing its suitability for thermoelectric applications. The analysis emphasizes the versatile capabilities of the NiMo6Se8 Chevrel phase, especially in thermoelectric and energy storage applications, showcasing its promising potential for future technological implementation.

Thermoelectric properties and thermal transport in two-dimensional GaInSe3 and GaInTe3 monolayers: A first-principles study

We here report the electronic structure calculation of GaInSe3 and GaInTe3 monolayers with the P3m1 (no. 156) space group. The electronic structure and thermoelectric properties of the monolayers are calculated through the Vienna Ab initio Simulation Package and BoltzTraP2 codes. The dynamic and thermodynamic stabilities were verified by calculating their phonon spectra and simulating ab initio molecular dynamics. The monolayers were found to have a direct bandgap, with both PBE + SOC and HSE06 + SOC potentials. The lattice thermal conductivity of GaInTe3 monolayer calculated using Phono3py code shows ultra-low values due to enhanced phonon–phonon scattering. Combining electrical and thermal transport, the values have been evaluated. Importantly, the p-type GaInTe3 has excellent thermoelectric properties at 700 K, with a zT value of 2, indicating that the p-type GaInTe3 has potential application in the field of thermoelectricity.

Enhanced response characteristics of NO2 gas sensor based on ultrathin SnS2 nanoplates: Experimental and DFT study

Layered-metal dichalcogenides with extraordinary characteristics of vast surface area, tunable bandgap and superior adsorption capability enable the potential for application in gas sensors. However, the synthesis of effective material for enhanced response performance remains a challenge. Herein, we exploited a fascinating sensitivity and selectivity towards NO2 gas detection using SnS2 nanoflakes prepared via the hydrothermal method. SnS2 nanoflakes with a thickness of 25 nm and an average diameter of approximately 500 nm show the potential for the detection of NO2 gas at low concentrations of ppb levels. The sensing properties of the SnS2 sensors were investigated for different concentrations of NO2 at various operating temperatures. The sensor exhibits the highest gas-sensing response of 161 at 250 οC upon exposure to 5 ppm of NO2 gas with fast response and recovery times. In addition, the sensor shows excellent selectivity with a low detection limit of ppb level. The electronic structure and gas-sensing mechanism are elucidated via finding density of states, charge density, and band structure based on DFT study which is calculated by the Vienna ab-initio simulation package (VASP). The considerable small adsorption energy reveals a physisorption of the NO2 molecules on the SnS2 surface (-0.174 eV), indicating the SnS2 nanoflakes are intriguing candidates for the speedy detection of NO2 gas.

Energy efficiency and performance analysis of a legacy atomic scale materials modeling simulator (VASP)

This work tackles the performance and energy consumption analysis of a legacy scientific application, the VASP (Vienna Ab-initio Simulation Package), an application commonly used by physicists and chemists for modeling materials at the atomic scale. Many of these scientific applications have been implemented in Fortran, where energy metrics instrumentation is not straightforward. We obtained performance figures (execution time and energy consumption) by instrumenting the source code using EML. This energy measurement library has been modified to introduce Fortran interfaces for these metrics. The analysis was carried out using different matrix algebra libraries, parallelization techniques, and hardware platforms, emphasizing on the MPI, OpenMP, and CUDA parallel implementations of the algorithms used in VASP. We employ various material specifications (atomic structures) and molecular sizes of a silicon-based crystal to create a set of benchmarks for these specifications, leading to some recommendations for final users regarding performance improvements. The proposed benchmarking technique assists the user in selecting the right combination of problem size, compilers, and parallelization options available in VASP. For a given system platform, the user will be able to determine not only the architecture to use (GPU or multicore processors), but also the appropriate library and parallelization according to the atomic structure and molecular size.

Precipitation and growth mechanism of TiN on Y2O3 (110) and (100) surfaces of yttrium-containing T4003 stainless steel

The Y2O3/TiN interface of main complex inclusions in T4003 stainless steel with rare earth element Yttrium (Y) content of 0.014 wt% was studied. The chemical composition and morphology of composite inclusions were characterized by scanning electron microscope (SEM) and X-ray energy dispersive spectrometer (EDS), the formation process of Y2O3/TiN composite inclusions was calculated by the first-principles calculations of the Vienna Ab-initio Simulation Package (VASP). Compared with Fe matrix, TiN is more likely to grow in an uneven state, and tends to exist in stainless steel in the form of Y2O3(100)/TiN(110). The main reason for the stability of the interface system is that the p orbitals of the O atoms at the interface forms covalent bonds with the d orbitals of the Ti atoms. It provides a new idea for improving the properties of stainless steel effectively by using heterogeneous nucleation of inclusions in steel.

Atomic Diffusivities of Yttrium, Titanium and Oxygen Calculated by Ab Initio Molecular Dynamics in Molten 316L Oxide-Dispersion-Strengthened Steel Fabricated via Additive Manufacturing.

Oxide-dispersion-strengthened (ODS) steels have long been viewed as a prime solution for harsh environments. However, conventional manufacturing of ODS steels limits the final product geometry, is difficult to scale up to large components, and is expensive due to multiple highly involved, solid-state processing steps required. Additive manufacturing (AM) can directly incorporate dispersion elements (e.g., Y, Ti and O) during component fabrication, thus bypassing the need for an ODS steel supply chain, the scale-up challenges of powder processing routes, the buoyancy challenges associated with casting ODS steels, and the joining issues for net-shape component fabrication. In the AM process, the diffusion of the dispersion elements in the molten steel plays a key role in the precipitation of the oxide particles, thereby influencing the microstructure, thermal stability and high-temperature mechanical properties of the resulting ODS steels. In this work, the atomic diffusivities of Y, Ti, and O in molten 316L stainless steel (SS) as functions of temperature are determined by ab initio molecular dynamics simulations. The latest Vienna Ab initio Simulation Package (VASP) package that incorporates an on-the-fly machine learning force field for accelerated computation is used. At a constant temperature, the time-dependent coordinates of the target atoms in the molten 316L SS were analyzed in the form of mean square displacement in order to obtain diffusivity. The values of the diffusivity at multiple temperatures are then fitted to the Arrhenius form to determine the activation energy and the pre-exponential factor. Given the challenges in experimental measurement of atomic diffusivity at such high temperatures and correspondingly the lack of experimental data, this study provides important physical parameters for future modeling of the oxide precipitation kinetics during AM process.

ElasTool v3.0: Efficient computational and visualization toolkit for elastic and mechanical properties of materials

Efficient computation and visualization of elastic and mechanical properties are crucial in the selection of materials and the design of new materials. The ElasTool v3.0 toolkit marks a significant advancement in the computational analysis and visualization of elastic and mechanical properties of materials, essential in material selection and design. This enhanced version extends beyond standard calculations like elastic tensor, Young's modulus, bulk modulus, and Poisson's ratio. It introduces capabilities for computing minimum thermal conductivity, linear compressibility, rendering the Christoffel equation, and elastic energy density. Notably, it integrates advanced visualization tools, including compatibility with Plotly and Elate web platforms for interactive web-based property exploration. A key feature of ElasTool v3.0 is the implementation of second-order elastic constants (SOECs) for tubular 2D-based nanostructures and nanotubes. Leveraging high-efficiency strain-matrix sets (OHESS), the toolkit now facilitates efficient computation of elastic constants and mechanical properties at both zero and finite temperatures for 1D, 2D, and 3D dimensions. ElasTool is openly accessible on GitHub: https://github.com/elastool. Program summaryProgram Title: ElasTool v3.0CPC Library link to program files:https://doi.org/10.17632/9hpgmspfmk.1Licensing provisions: GNU General Public License, version 3Programming language: Python 3Nature of problem: A computational toolkit that can efficiently compute the elastic and mechanical properties of materials using ab initio methods and provide advanced visualization capabilities is crucial to accurate material analysis.Solution method:ElasTool addresses this need by incorporating high-efficiency strain-matrix sets that require only small strain sets with ab initio methods such as density functional theory and molecular dynamics. This package is versatile and adept at calculating elastic constants for both 2D and 3D structures, including tubular 2D-based nanostructures and nanotubes and 1D nanoribbons at zero and finite temperatures. Beyond basic elastic moduli, ElasTool now boasts advanced visualization capabilities and computes additional critical properties such as minimum thermal conductivity, elastic energy density, and linear compressibility. This enhancement transforms ElasTool into a comprehensive tool for accelerated material design and modeling, and advanced visualization, offering both efficiency and depth in material property analysis.Additional comments including restrictions and unusual features:ElasTool is currently integrated with the Vienna Ab initio Simulation Package (VASP). Adapting to other electronic structures is straightforward.ElasTool offers serial and parallel execution. Depending on the system size, the number of strains, the nature of the calculations (zero or finite temperature), and the computing architecture, the computing time can vary from a minute to as long as 15 hours.