An official website of the United States government
The present work is a novel, systematic study of the effect of density functional theory input parameters on the vacancy formation energy (VFE), migration barrier for diffusion, and electronic structure for each element in the CoCrNi medium-entropy alloy (MEA). In particular, the novelties include: (1) calculating the aforementioned properties of Co, Cr, or Ni, in the CoCrNi MEA using magnetic and non-magnetic states, and two versions of the generalized gradient approximation: Perdew, Burke, and Ernzerhof (PBE) and the PBE version for solids (PBEsol), and (2) a detailed comparison of 0 K activation energy to experimental creep activation energies. First-principles calculations at 0 K are performed using the Vienna ab-initio simulation package. Special quasirandom structures (SQS) and Widom-type substitution are employed. For each element, Co, Cr, or Ni, non-magnetic calculations result in a higher VFE and larger range of calculated values for the configurations studied. The averaged migration barrier is the highest for Co in the CoCrNi for three of four sets of calculation parameters in the configurations studied. Finally, the results indicate that the average 0 K activation energy for diffusion makes up 70–80% of the experimental creep activation energy, depending on the exchange-correlation functional employed.
more »
« less
The binding of atomic hydrogen on graphene from density functional theory and diffusion Monte Carlo calculations
In this work, density functional theory (DFT) and diffusion Monte Carlo (DMC) methods are used to calculate the binding energy of a H atom chemisorbed on the graphene surface. The DMC value of the binding energy is about 16% smaller in magnitude than the Perdew–Burke–Ernzerhof (PBE) result. The inclusion of exact exchange through the use of the Heyd–Scuseria–Ernzerhof functional brings the DFT value of the binding energy closer in line with the DMC result. It is also found that there are significant differences in the charge distributions determined using PBE and DMC approaches.
more »
« less
- Award ID(s):
- 2028826
- PAR ID:
- 10366366
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 156
- Issue:
- 14
- ISSN:
- 0021-9606
- Page Range / eLocation ID:
- Article No. 144702
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The field of 2D materials has grown dramatically in the past two decades. 2D materials can be utilized for a variety of next-generation optoelectronic, spintronic, clean energy, and quantum computing applications. These 2D structures, which are often exfoliated from layered van der Waals materials, possess highly inhomogeneous electron densities and can possess short- and long-range electron correlations. The complexities of 2D materials make them challenging to study with standard mean-field electronic structure methods such as density functional theory (DFT), which relies on approximations for the unknown exchange-correlation functional. To overcome the limitations of DFT, highly accurate many-body electronic structure approaches such as diffusion Monte Carlo (DMC) can be utilized. In the past decade, DMC has been used to calculate accurate magnetic, electronic, excitonic, and topological properties in addition to accurately capturing interlayer interactions and cohesion and adsorption energetics of 2D materials. This approach has been applied to 2D systems of wide interest, including graphene, phosphorene, MoS2, CrI3, VSe2, GaSe, GeSe, borophene, and several others. In this review article, we highlight some successful recent applications of DMC to 2D systems for improved property predictions beyond standard DFT.more » « less
-
The acidity of anatase titania before and after KOH doping was probed by pyridine adsorption in a pulse microreactor and modeled by DFT optimization of the geometry of CO and pyridine adsorption on a periodic slab of (101) and (100) surfaces using a GGA/PBE functional and verified by an example of a single-point calculation of the optimized geometry using an HSE-06 hybrid functional. The anatase (101) surface was slightly more acidic compared to the (100) surface. Both experimental and computational methods show that the acidity of anatase surfaces decreased after KOH doping and increased after the dissociative adsorption of water. Higher acidity of Ti metal centers was indicated by the shortening of the Ti-N, Ti-C, and C-O bond lengths, increasing the IR frequency of CO and pyridine ring vibrations and energy of adsorption. The DFT calculated energy of pyridine adsorption was analyzed in terms of binding energy and the energy of lattice distortion. The latter was used to construct Hammett plots for the adsorption of 4-substituted pyridines with electron-donating and -withdrawing substituents. The Hammett rho constant was obtained and used to characterize the acidity of various metal centers of −1.51 vs. −1.46 on pristine (101) and (100) surfaces, which were lowered to −1.07 and −1.19 values on KOH-doped (101) and (100) surfaces, respectively. The mechanism of lowering surface acidity via KOH doping proceeds through the stabilization of the atomic structure of Lewis acid centers. When an alkaline metal cation binds to several lattice oxygen atoms, the surface structure becomes more rigid. The ability of Ti atoms to move toward the adsorbate is restricted. Consequently, the lattice distortion energy and binding energy are decreased. In contrast, higher flexibility of the outermost layer of Ti atoms as a result of electron density redistribution, for example, in the presence of water on the surface, allows them to move farther outward, make shorter contacts with the adsorbate, and attain higher energies of binding and lattice distortion.more » « less
-
Lattice thermal conductivity (κL) is a crucial characteristic of crystalline solids with significant implications for thermal management, energy conversion, and thermal barrier coating. The advancement of computational tools based on density functional theory (DFT) has enabled the effective utilization of phonon quasi-particle-based approaches to unravel the underlying physics of various crystalline systems. While the higher order of anharmonicity is commonly used for explaining extraordinary heat transfer behaviors in crystals, the impact of exchange-correlation (XC) functionals in DFT on describing anharmonicity has been largely overlooked. The XC functional is essential for determining the accuracy of DFT in describing interactions among electrons/ions in solids and molecules. However, most XC functionals in solid-state physics are primarily focused on computing the properties that only require small atomic displacements from the equilibrium (within the harmonic approximation), such as harmonic phonons and elastic constants, while anharmonicity involves larger atomic displacements. Therefore, it is more challenging for XC functionals to accurately describe atomic interactions at the anharmonicity level. In this study, we systematically investigate the room-temperature κL of 16 binary compounds with rocksalt and zincblende structures using var- ious XC functionals such as local density approximation (LDA), Perdew-Burke-Ernzerhof (PBE), revised PBE for solid and surface (PBEsol), optimized B86b functional (optB86b), revised Tao-Perdew-Staroverov-Scuseria (revTPSS), strongly constrained and appropriately normed functional (SCAN), regularized SCAN (rSCAN) and regularized-restored SCAN (r2SCAN) in combination with different perturbation orders, including phonon within harmonic approximation (HA) plus three- phonon scattering (HA+3ph), phonon calculated using self-consistent phonon theory (SCPH) plus three-phonon scattering (SCPH+3ph), and SCPH phonon plus three- and four-phonon scattering (SCPH+3,4ph). Our results show that the XC functional exhibits strong entanglement with perturbation order and the mean relative absolute error (MRAE) of the computed κL is strongly influenced by both the XC functional and perturbation order, leading to error cancellation or amplification. The minimal (maximal) MRAE is achieved with revTPSS (rSCAN) at the HA+3ph level, SCAN (r2SCAN) at the SCPH+3ph level, and PBEsol (rSCAN) at the SCPH+3,4ph level. Among these functionals, PBEsol exhibits the highest accuracy at the highest perturbation order. The SCAN- related functionals demonstrate moderate accuracy but are suffer from numerical instability and high computational costs. Furthermore, the different impacts of quartic anharmonicity on κL in rocksalt and zincblende structures are identified by all XC functionals, attributed to the distinct lattice anharmonicity in these two structures. These findings serve as a valuable reference for selecting appropriate functionals for describing anharmonic phonons and offer insights into high-order force constant calculations that could facilitate the development of more accurate XC functionals for solid materials.more » « less
-
Abstract Inspired by significant local distortions found near vacancies in a neutron pair distribution function analysis study (G. King et al.,Inorg. Chem.2012, 51, 13060) of Sr2MSbO5.5(M = Ca and Sr), this computational study finds minimum‐energy structures with these and related distortions using density functional theory (DFT) with the Perdew‐Burke‐Ernzerhof (PBE) functional as implemented in the Vienna Ab Initio Simulations Package (VASP) (G. Kresse and J. Furthmüller,Phys. Rev. B, 1996, 54, 11169; G. Kresse and J. Hafner,Phys. Rev. B, 1993, 47, 558; G. Kresse and J. Furthmüller,Comput. Mater. Sci., 1996, 6, 15). All structures were optimized using the conjugate gradient method. The global minima found for both systems featured trigonal bipyramid SbO5structures and edge sharing with M‐centered polyhedra. However, while calcium ions occupied full and partial octahedra, the larger strontium ions were more commonly found in full and partial pentagonal bipyramids. Molecular dynamics with velocity rescaling at1200K revealed movements of the oxygen vacancy via polyhedral rotations. This work highlights the need to consider both square pyramid to trigonal bipyramid rearrangements around small ions and rotational polyhedral movements in simulating oxygen vacancy conduction in oxygen‐deficient double perovskites.more » « less
