Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density functional theory development”

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

Importance of bulk excitations and coherent electron-photon-phonon scattering in photoemission from PbTe(111): Ab initio theory with experimental comparisons

This paper presents a fully ab initio many-body photoemission framework that includes coherent three-body electron-photon-phonon scattering to predict the transverse momentum distributions and the mean transverse energies (MTEs) of bulk photoelectrons from single-crystal photocathodes. The need to develop such a theory stems from the lack of studies that provide complete understanding of the underlying fundamental processes governing the transverse momentum distribution of photoelectrons emitted from single crystals. For example, initial predictions based on density-functional theory calculations of effective electron masses suggested that the (111) surface of PbTe would produce very small MTEs (≤ 15 meV), whereas our experiments yielded MTEs ten to twenty times larger than these predictions and also exhibited a lower photoemission threshold than predicted. The ab initio framework presented in this paper correctly reproduces both the magnitude of the MTEs from our measurements in PbTe(111) and the observed photoemission below the predicted threshold. Our results show that both photoexcitations into states that propagate in the bulk of the material and coherent many-body electron-photon-phonon scattering processes, which initial predictions ignored, play surprisingly important roles in photoemission from PbTe(111). Lastly, from the lessons learned, we recommend a procedure for rapid computational screening of potential single-crystal photocathodes for applications in next-generation ultrafast electron diffraction and X-ray free-electron lasers, which will enable significant advances in condensed matter research.

36 MATERIALS SCIENCE↗

Global and Local Connectivities Describe Hydrogen Intercalation in Metal Oxides

We introduce global connectivity, Ω, which quantifies the detailed connectivity of a material and accurately captures reactivity behavior. Here, we demonstrate that combining global and local connectivity describes how metal oxides interact with hydrogen. Blending density functional theory, graph theory, and machine learning we built a reactivity model which accurately predicts hydrogen intercalation potentials of different metal oxides experimentally measured in the lab. The use of global connectivity can accelerate materials design through the development of novel structure-property relationships.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational design of high-entropy rare earth aluminum garnets for advanced thermal and environmental barrier coatings

To enhance the protection of Ni-based superalloys in gas turbine engines’ high-temperature environments, it’s crucial to develop advanced thermal/environmental barrier coating (T/EBC) materials with a balanced combination of thermal and mechanical properties. This optimization is essential to safeguard against chemical and thermal challenges. Here, in this study, we harness the power of density functional theory (DFT) in conjunction with combinatorial chemistry methodologies to engineer high-performance high-entropy rare earth disilicates of the RE 3 Al 5 O 12 family (where RE denotes Y, Gd, Er, and Yb). These materials are meticulously designed to exhibit superior phase stability, a targeted coefficient of thermal expansion (CTE), low lattice thermal conductivity, and robust mechanical properties. The determination of CTE values is accomplished through phonon calculations at various volume settings within the quasi-harmonic approximation, while lattice thermal conductivities are rigorously assessed employing the Debye-Callaway model, accounting for three distinct phonon processes. Our findings highlight the remarkable attributes of the solid solution (Y 1/4 Gd 1/4 Er 1/4 Yb 1/4 ) 3 Al 5 O 12 , which displays a reduction in lattice thermal conductivity compared to its individual constituents while maintaining a favorable range of CTE values. The novel T/EBC material, distinguished by their multifaceted functionalities, are poised to use in substantial enhancements in the performance of engines.

36 MATERIALS SCIENCE↗

A computational approach to complex junctions and interfaces

Molecular- and nano- junctions, surfaces, interfaces, and interfacial processes have been the subjects of this DOE supported program for about 20 years. We have investigated electronic and magnetic structure, spin-dependent charge transport, and emerging phenomena at surfaces and interfaces in a number of physical systems ranging from nano-clusters, to nanowires, to two-dimensional (2D) crystals and their interactions, to bulk matter. Eighty-four refereed papers were published [1-84], with a majority effort from the group led by the PI, and over a hundred presentations were made at international conferences. We use first-principles methods in the framework of density functional theory in conjunction with nano-equilibrium Green function techniques, semi-empirical Boltzmann transport theory, tight-binding models, and the Hubbard model and many-body wavefunction approaches. In pursuing our scientific goals, we have also developed a number of computational algorithms and computational packages.

2D materials↗

Towards Elucidating Structure–Spectra Relationships in Rhamnogalacturonan II: Computational Protocols for Accurate 13C and 1H Shifts for Apiose and Its Borate Esters

Apiose is a naturally occurring, uncommon branched-chain pentose found in plant cell walls as part of the complex polysaccharide Rhamnogalacturonan II (RG-II). The structural elucidation of the three-dimensional structure of RG-II by nuclear magnetic resonance (NMR) spectroscopy is significantly complicated by the ability of apiose to cross-link via borate ester linkages to form RG-II dimers. Here, we developed a computational approach to gain insight into the structure–spectra relationships of apio–borate complexes in an effort to complement experimental assignments of NMR signals in RG-II. Our protocol involved structure optimizations using density functional theory (DFT) followed by isotropic magnetic shielding constant calculations using the gauge-invariant atomic orbital (GIAO) approach to predict chemical shifts. We evaluated the accuracy of 23 different functional–basis set (FBS) combinations with and without implicit solvation for predicting the experimental 1 H and 13 C shifts of a methyl apioside and its three borate derivatives. The computed NMR predictions were evaluated on the basis of the overall shift accuracy, relative shift ordering, and the ability to distinguish between dimers and monomers. We demonstrate that the consideration of implicit solvation during geometry optimizations in addition to the magnetic shielding constant calculations greatly increases the accuracy of NMR chemical shift predictions and can correctly reproduce the ordering of the 13 C shifts and yield predictions that are, on average, within 1.50 ppm for 13 C and 0.12 ppm for 1 H shifts for apio–borate compounds.

BASIC BIOLOGICAL SCIENCES,BIOMASS FUELS↗

A Robust Methodology to Elucidate Kinetics of Room Temperature Electrochemical Propane Adsorption on Platinum

Electrocatalytic activation of alkanes can further decarbonize chemical manufacturing by leveraging affordable renewable electricity and readily available shale gas reserves in the United States. Earlier works have identified the unique role of Pt in adsorbing and activating alkanes, like propane, at room temperature in acidic, aqueous electrolytes, revealing spontaneous formation of deeply dehydrogenated propane-derived surface species with an intact C 3 - backbone. Although an adsorption mechanism was hypothesized, it has not been explicitly investigated to date, preventing the quantification of kinetic rate parameters. A robust methodology to investigate and benchmark propane adsorption kinetics on Pt is critical for the rational design of electrocatalysts that exhibit higher selectivity toward desired partially oxidized products. Herein, we analyze an oxidative current transience that appears during the adsorption of propane on Pt in aqueous electrochemical conditions and develop a methodology that elucidates the adsorption mechanism and enables quantification of rate parameters such as order dependences and apparent activation barriers. This method yields an expected first-order dependence with respect to propane concentration at low coverage and reveals a second-order dependence with respect to the concentration of surface active sites. Additionally, the apparent activation barrier for propane adsorption was calculated using an Arrhenius analysis of the current transience under temperature control. The experimentally measured activation barrier of 35 kJ mol –1 is in excellent agreement with the theoretical barrier calculated by density functional theory (DFT). The kinetic analysis was extended, via the use of transition state theory, to extract entropy and enthalpy of activation, yielding consistent results with the proposed two-step adsorption mechanism and DFT calculations. These results demonstrate reliable quantification of kinetic parameters for electrocatalytic activation of C–H bonds in alkanes that can be employed for rational catalyst development for a versatile range of electrocatalytic conditions.

alkane activation↗

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↗

Discovering Dinuclear Dioxygen-Bridged Cobalt(III) Complexes for Selective Binding of O 2 from Air

The design and development of dioxygen activation in porous crystalline materials is a useful avenue for exploring selective adsorption of O 2 that shows significant potential to enable separations of O 2 and N 2 from air. Porous materials featuring redox-active metal centers have received attention regarding selective O 2 adsorption via chemisorption. Drawing inspiration from a dinuclear cobalt material ([(Co(III) 2 (bpbp)O 2 ) 2 bdc](PF 6 ) 4 (CSD code: GAMVIB; bpbp – = 2,6-bis(N,N-bis(2-pyridylmethyl)aminomethyl)-4-tert-butylphenolato; bdc 2– = 1,4-benzenedicarboxylato)) that displays reversible and selective O 2 adsorption, we focus on searching for potential O 2 -selective materials with dinuclear cobalt clusters that have dioxygen-bridged Co(III) complexes. We combine structure screening with a high-level hybrid periodic density functional theory (DFT) workflow to investigate O 2 and N 2 adsorption in materials from validated crystal structure databases (e.g., the CSD database). These calculations identify multiple materials that are predicted to have superior O 2 binding capability relative to GAMVIB. Grand Canonical Monte Carlo (GCMC) simulations based on DFT-developed force fields were performed for selected candidates to estimate the adsorption performance for O 2 /N 2 mixtures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Dual active site tandem catalysis of metal hydroxyl oxides and single atoms for boosting oxygen evolution reaction

We report the high voltage in oxygen evolution reaction (OER) often causes structural change in electrocatalysts and forms multiple active sites. Therefore, exploring the synergy of various active sites is extremely significant to develop catalysts for OER with multiple elementary steps. Herein, we adopt the physically adsorbed metal ions method to successfully synthesize a highly-efficient electrocatalyst containing the dual active sites of Fe-NiOOH and NiC4 single atoms (marked as Ni SAs/Fe-NiOOH). The Ni SAs/Fe-NiOOH displays outstanding OER performance with an overpotential of 269 mV to deliver current density of 10 mA/cm 2 that shows 55 mV superior to commercial IrO 2 /CB at the same condition. Experiments and density functional theory calculations indicate that the excellent OER activity of Ni SAs/Fe-NiOOH catalyst is attributed to the synergy of dual active sites of NiC 4 SAs and Fe-NiOOH. A tandem catalysis mechanism is also proposed to reveal the synergism of two active centers, which makes the potential-determining step more facile and accordingly decreases the OER overpotential. This work offers a new concept of tandem catalysis to develop the electrocatalysts with many elemental steps, like OER and oxygen reduction reaction catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanistic Study of Functional Electrolyte Solvents for High-Voltage Lithium Batteries

The pervasive use of Ni-rich cathode active materials, e.g., LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811), for high-energy-density Li-ion batteries (LIBs) has been hindered by rapid battery capacity decay when cycled with high charge cutoff voltages due to electrolyte decomposition in the conventional carbonate solvent-based electrolytes, oxidative parasitic side reactions at the electrolyte/cathode interface, and irreversible phase changes in the cathode active materials leading to dissolution of transition metals into the electrolytes. Various functional electrolyte solvents have been studied to tackle the above technical challenges, yet the roles of individual solvents in the performance of LIBs remain poorly understood. Here, in this study, we systematically investigate electrochemical performance mechanisms of fluorinated and organosilicon single solvents and cosolvents, for the first time, in high-voltage Li/NMC811 batteries, using electrochemical and analytical characterizations and density functional theory modeling. We observe that some unique combinations of the functional solvents can lead to exceptionally stable high-voltage cycle performance in the Ni-rich cathode-based LIBs. Our mechanistic study reveals that the synergistic effect of solvents plays a vital role in enabling electrochemical stability at both the Ni-rich cathode and the Li metal anode. Understanding the electrochemical performance mechanisms of functional solvents can greatly help in designing and formulating advanced electrolytes that enable the development of high-voltage, high-energy-density, long-cycle-life lithium batteries.

density functional theory modeling↗

[(MeCN)Ni(CF 3 ) 3 ] - and [Ni(CF 3 ) 4 ] 2– : Foundations toward the Development of Trifluoromethylations at Unsupported Nickel

Nickel anions [(MeCN)Ni(CF 3 ) 3 ] - and [Ni(CF 3 ) 4 ] 2– were prepared by the formal addition of 3 and 4 equiv, respectively, of AgCF 3 to [(dme)NiBr 2 ] in the presence of the [PPh 4 ] + counterion. Detailed insights into the electronic properties of these new compounds were obtained through the use of density functional theory (DFT) calculations, spectroscopy-oriented configuration interaction (SORCI) calculations, X-ray absorption spectroscopy, and cyclic voltammetry. The data collectively show that trifluoromethyl complexes of nickel, even in the most common oxidation state of nickel(II), are highly covalent systems whereby a hole is distributed on the trifluoromethyl ligands, surprisingly rendering the metal to a physically more reduced state. In the cases of [(MeCN)Ni(CF 3 ) 3 ] - and [Ni(CF 3 ) 4 ] 2- , these complexes are better physically described as d 9 metal complexes. [(MeCN)Ni(CF 3 ) 3 ] - is electrophilic and reacts with other nucleophiles such as phenoxide to yield the unsupported [(PhO)Ni(CF 3 ) 3 ] 2– salt, revealing the broader potential of [(MeCN)Ni(CF 3 ) 3 ] - in the development of “ligandless” trifluoromethylations at nickel. Proof-in-principle experiments show that the reaction of [(MeCN)Ni(CF 3 ) 3 ] - with an aryl iodonium salt yields trifluoromethylated arene, presumably via a high-valent, unsupported, and formal organonickel(IV) intermediate. Evidence of the feasibility of such intermediates is provided with the structurally characterized [PPh 4 ] 2 [Ni(CF 3 ) 4 (SO 4 )], which was derived through the two-electron oxidation of [Ni(CF 3 ) 4 ] 2– .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Foaming prediction in pure liquids from dimensionless numbers inspired by the theory of fluid behavior for drops

Foaming prediction is critical for selecting materials and designing processes in industries such as bioprocessing and gas processing. Existing models lack the generality needed for a wide range of materials and overlook the foaming behavior in pure liquids. Here, this work presents a novel method for predicting foaming in pure liquids based on their density, surface tension, and viscosity, using Reynolds ( Re ) and Ohnesorge ( Oh ) numbers. A foaming prediction map, leveraging the theory of fluid drop behavior, was developed by plotting these numbers. This map delineates distinct non-foaming and foaming regions, functioning as a binary classifier for foaming predictions. The map was fitted and validated through shake test experiments on 46 liquids, demonstrating reliable predictions, except for a specific region characterized by small Oh and large Re numbers. This region corresponded to relatively low foam stability and high turbulence, making foaming predictions challenging for liquids in this category.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Visualizing Oxidation Mechanisms in Few-Layered Black Phosphorus via In Situ Transmission Electron Microscopy

Layered two-dimensional (2D) black phosphorus (BP) exhibits novel semiconducting properties including a tunable bandgap and high electron mobility. However, the poor stability of BP in ambient environment severely limits potential for application in future electronic and optoelectronic devices. While passivation or encapsulation of BP using inert materials/polymers has emerged as a plausible solution, a detailed fundamental understanding of BP's reaction with oxygen is imperative to rationally advance its use in applications. Here, we use in situ environmental transmission electron microscopy to elucidate atomistic structural changes in mechanically exfoliated few-layered BP during exposure to varying partial pressures of oxygen. An amorphous oxide layer is seen on the actively etching BP edges, and the thickness of this layer increases with increasing oxygen partial pressure, indicating that oxidation proceeds via initial formation of amorphous P x O y species which sublime to result in the etching of the BP crystal. We observe that while few-layered BP is stable under the 80 kV electron beam (e-beam) in vacuum, the lattice oxidizes and degrades at room temperature in the presence of oxygen only in the region under the e-beam. The oxidative etch rate also increases with increasing e-beam dosage, suggesting the presence of an energy barrier for the oxidation reaction. Preferential oxidative etching along the [0 0 1] and [0 0 1] crystallographic directions is observed, in good agreement with density functional theory calculations showing favorable thermodynamic stability of the oxidized BP (0 0 $\bar{1}$) planes compared to the (1 0 0) planes. Finally, we expect the atomistic insights and fundamental understanding obtained here to aid in the development of novel approaches to integrate BP in future applications.

2D materials↗

2D in-Plane Ordered MXene Nanosheets Derived from (Mo 2/3 Er 1/3 ) 2 AlC Rare-Earth i-MAX for Energy Storage Applications

MXenes have become one of the most versatile families of two-dimensional (2D) materials due to their high conductivity, hydrophilicity, and remarkable electrochemical performance. This has stimulated intense efforts to design and synthesize MXenes, including structurally unique in-plane ordered 2D MXenes called i-MXenes. Here, we have synthesized the quaternary rare earth (RE)-based i-MAX phase (Mo 2/3 Er 1/3 ) 2 AlC using an arc melting method, and the corresponding 2D i-MXene was then obtained through a LiF/HCl soft etching process. Literature studies have shown that Al and the RE element are etched out during the etching process, leading to the formation of pure vacancy-ordered Mo1.33C 2D i-MXene. However, our investigation reveals that upon exposure to a fluorine solution, the i-MAX phase forms RE fluoride impurities, which are challenging to remove through HCl−DI water washing and persist in the final product, resulting in impure Mo 1.33 C@Er i-MXene. These results were confirmed by various characterizations such as X-ray diffraction, Raman spectroscopy, X-ray photoelectron spectroscopy, and scanning transmission electron microscopy. Although the Mo 1.33 C@Er electrode showed a 24-fold increase in specific capacitance compared to its parent i-MAX phase, it still exhibited a high charge-transfer resistance arising from the insulating nature of RE fluoride byproducts, which adversely influence the overall capacitance behavior of the synthesized 2D Mo 1.33 C@Er i-MXenes. This study contributes to identifying pathways for the preparation of pure 2D i-MXenes from RE-based i-MAX phases and developing improved synthesis methods. With additional process optimization, the 2D i-MXene holds a strong potential for electrochemical energy storage applications. Additionally, the electronic structures of Mo 1.33 C were theoretically studied using first-principles density functional theory calculations, which revealed that pristine Mo 1.33 C is metallic, and this metallic nature is preserved even with −O, −F, and mixed functionalization.

chemical structure↗

Gold–Thiolate Nanocluster Dynamics and Intercluster Reactions Enabled by a Machine Learned Interatomic Potential

Monolayer protected metal clusters comprise a rich class of molecular systems and are promising candidate materials for a variety of applications. While a growing number of protected nanoclusters have been synthesized and characterized in crystalline forms, their dynamical behavior in solution, including prenucleation cluster formation, is not well understood due to limitations both in characterization and first-principles modeling techniques. Recent advancements in machine-learned interatomic potentials are rapidly enabling the study of complex interactions such as dynamical behavior and reactivity on the nanoscale. Here, we develop an Au-S-C-H atomic cluster expansion (ACE) interatomic potential for efficient and accurate molecular dynamics simulations of thiolate-protected gold nanoclusters (Au n (SCH 3 ) m ). Trained on more than 30,000 density functional theory calculations of gold nanoclusters, the interatomic potential exhibits ab initio level accuracy in energies and forces and replicates nanocluster dynamics including thermal vibration and chiral inversion. Long dynamics simulations (up to 0.1 μs time scale) reveal a mechanism explaining the thermal instability of neutral Au 25 (SR) 18 clusters. Specifically, we observe multiple stages of isomerization of the Au 25 (SR) 18 cluster, including a chiral isomer. Additionally, we simulate coalescence of two Au 25 (SR) 18 clusters and observe series of clusters where the formation mechanisms are critically mediated by ligand exchange in the form of [Au-S] n rings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Growth of Low-Defect WSe 2 Film via High-Purity van der Waals Crystal Precursor

Two-dimensional (2D) semiconducting transition metal dichalcogenides (TMDs) exhibit exceptional electrical and optical properties, empowering their promising prospects for future nanoelectronics. Despite major advances in n-type 2D semiconductors, the field has yet to synthesize high-mobility p-type 2D TMDs, in particular WSe 2 , and systematically query the influence of defects. In this study, we unveil the pivotal role of substitutional impurity defects vis-à-vis the precursor used and growth method employed in defining the quality of 2D p-type WSe 2 . Density functional theory calculations suggest the adverse effect of Fe-, Co-, Ni- and Si-substituted W impurity defects on the mobility of WSe 2 , whereas defects such as O-, S-substituted Se and Mo-substituted W pose negligible impact. Guided by the theory, we pinpoint van der Waals (vdW) crystals, commonly used in mechanical exfoliation, as the optimal precursor, and develop a facile vdW crystal physical vapor deposition (PVD) method to grow high-purity monolayer 2D WSe 2 film (VPVD-WSe 2 ) that is continuous across a centimeter scale. A suite of spectroscopies confirms the markedly reduced defect density of the as-synthesized WSe 2 compared to those by typical chemical vapor deposition methods, and by PVD with commercial or hydrothermal precursors. Scanning tunneling microscopy further evidence the ultralow substitutional impurity defect density of VPVD-WSe 2 , greatly outperforming the control samples and approaching the mechanically exfoliated counterparts. The VPVD-WSe 2 based field-effect transistors exhibit notable electrical performance with record-high field-effect hole mobility up to 112 cm 2 V –1 s –1 at room temperature, exceeding the best-reported monolayer WSe 2 synthesized by chemical vapor deposition and rivaling the mechanically exfoliated 2D WSe 2 flakes.

WSe2↗