Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy density functional”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Analytic Nuclear Gradients for Complete Active Space Linearized Pair-Density Functional Theory

Accurately modeling photochemical reactions is difficult due to the presence of conical intersections and locally avoided crossings, as well as the inherently multiconfigurational character of excited states. As such, one needs a multistate method that incorporates state interaction in order to accurately model the potential energy surface at all nuclear coordinates. The recently developed linearized pair-density functional theory (L-PDFT) is a multistate extension of multiconfiguration PDFT, and it has been shown to be a cost-effective post-MCSCF method (as compared to more traditional and expensive multireference many-body perturbation methods or multireference configuration interaction methods) that can accurately model potential energy surfaces in regions of strong nuclear–electronic coupling in addition to accurately predicting Franck–Condon vertical excitations. Here, in this paper, we report the derivation of analytic gradients for L-PDFT and their implementation in the PySCF-forge software, and we illustrate the utility of these gradients for predicting ground- and excited-state equilibrium geometries and adiabatic excitation energies for formaldehyde, s-trans-butadiene, phenol, and cytosine.

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

Surface Charge and Electrostatic Spin Crossover Effects in CoN 4 Electrocatalysts

Carbon materials doped with nitrogen and 3d transition metals have attracted a great deal of interest for catalyzing electrochemical reactions such as water splitting, oxygen reduction, and carbon dioxide reduction. Here, we employed density functional theory to study Co–N-doped carbon as electrocatalysts for the oxygen reduction and oxygen evolution reactions. Specifically, we investigated the interplay among adsorption energies, the spin state of the CoN 4 active center, and the applied potential. We found that adsorption energies strongly depend on both the applied potential and the spin state of the Co center. Furthermore, spin state transitions induced by the applied potential also play an important role in determining the adsorption energies. Here, this effect originates from a different potential of zero charge and capacitance of each spin state.

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

Computational Approaches for Clean Energy Materials

Currently, 80% of the global final energy consumption occurs in form of fuels and only 20% as electricity. On the other hand, renewable energy additions come almost exclusively in the form of electricity (dominantly photovoltaics and wind). Thus, a successful energy transition will require enormous growth in renewables, sufficient to convert excess electricity into fuels, as well as the development of non-electricity based solar fuel technologies. As much as photovoltaic capacities have grown over the past 20 years, it is far from clear that current technologies and materials are up to the task to grow from here by yet another factor 100 until 2050. Therefore, sustained research efforts on emerging inorganic semiconductors for solar electricity and fuels are essential for facing the double challenge of climate change and energy security. Computational materials science can make important contributions, guiding and supporting research activities through both materials search and discovery and through detailed studies that help to develop a mechanistic understanding of materials performance and bottlenecks. This presentation will highlight three recent computational projects with relevance for photovoltaics and solar fuels (1) Defect graph neural networks (dGNN) for materials discovery in solar thermochemical hydrogen (STCH) [1]. The dGNN approach facilitates broad and fast materials screening for defect properties. (2) Modeling highly off-stoichiometric systems by evaluating the free energy of defect interaction [2]. This approach allows quantitative prediction of H2 production in complex STCH oxides. (3) First-principles atomic structure prediction for interfaces [3]. This work showed how an atomically thin CdCl2 interlayer phase enables in principle ideal electron transport across the incommensurate SnO2/CdTe interface. [1] M.D. Witman, A. Goyal, T. Ogitsu, A.H. McDaniel, S. Lany, Nat. Comput. Sci. (2023). https://doi.org/10.1038/s43588-023-00495-2. [2] A. Goyal, M.D. Sanders, R.P. O'Hayre, S. Lany, PRX Energy 3, 013008 (2024). https://doi.org/10.1103/PRXEnergy.3.013008. [3] A. Sharan, M. Nardone, D. Krasikov, N. Singh, S. Lany, Appl. Phys. Rev. 9, 041411 (2022). https://doi.org/10.1063/5.0104008.

density functional theory↗
DOE OSTI2021

Using distributions to understand neutron and x-ray production in ICF ignition capsules and other high energy density plasmas

Here this paper describes how x-ray and neutron distribution functions can be useful tools to visualize the conditions measured in many types of plasma physics experiments. In particular, we model a standard inertial confinement fusion ignition capsule that consists of a Si doped plastic ablator surrounding a layer of deuterium–tritium (DT) ice as the yield varies from 18 kJ to 16.7 MJ and use the distribution functions to show that neutrons and high energy x rays (15 keV) are produced under similar conditions when the yield is low. However, as the capsule starts to support a propagating burn due to alpha heating, the x rays and neutrons are produced under somewhat different conditions in different parts of the plasma. In particular, the x-ray production takes place mainly in the hot plastic ablator for the full yield ignition capsule under quite different plasma conditions from the DT region producing the 14 MeV neutrons, which results in x-ray images with larger radii than the corresponding neutron images. These same distribution functions can be applied to many other plasma physics experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Using data-science approaches to unravel insights for enhanced transport of lithium ions in single-ion conducting polymer electrolyte

Solid polymer electrolytes have yet to achieve the an ionic conductivity > 1 mS/cm at room temperature for realistic applications. This target implies the need to reduce the effective energy barriers of ion transport in polymer electrolytes to around 20 kJ/mol. In this work, we combine information extracted from existing experimental results with theoretical calculations to provide insights into ion transport in single-ion conductors (SICs) with a focus on lithium ion SICs. Through the analysis of temperature-dependent ionic conductivity data obtained from the literature, we evaluate different methods of extracting energy barriers for lithium transport. The traditional Arrhenius fit to the temperature-dependent ionic conductivity data indicates that the Meyer-Neldel rule holds for SICs. However, the values of the fitting parameters remain unphysical. Our modified approach based on recent work (Macromolecules, 56, 15, 6051(2023)), which incorporates a fixed pre-exponential factor, reveals that the energy barriers exhibit temperature dependence over a wide range of temperatures. Using this approach, we identify a series of anions leading to the energy barriers less than 30 kJ/mol, which include trifluoromethane sulfonimide (TFSI), fluoromethane sulfonimide (FSI), and boron-based organic anions. In our efforts to design the next generation of anions, which can exhibit the energy barriers less than 20 kJ/mol, we focused on boron-containing SICs, and performed density functional theory (DFT) based calculations to connect the chemical structures via the binding energy of cation (lithium)-anion pairs with the experimentally derived effective energy barriers for ion transport. Not only have we identified a correlation between the binding energy and the energy barriers, but we also propose a strategy to design new boron-based anions by using the correlation. This combined approach involving experiments and theoretical calculations is capable of facilitating the identification of promising new anions, which can exhibit ionic conductivity $> 1$ mS/cm near room temperature, thereby expediting the development of novel superionic single-ion conducting polymer electrolytes. The published datasets include all the temperature-dependent ionic conductivity collected from the literature with literature DOIs, DFT calculated binding energies, and python scripts to analyze data, construct statistical models, and generate plots.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Design Strategies of Spinel Oxide Frameworks Enabling Reversible Mg-Ion Intercalation

In this study, reversible Mg 2+ intercalation in metal oxides frameworks is a key enabler for an operational Mg ion battery with high energy density needed for the next generation of energy storage technologies. While functional Mg-ion batteries have been achieved in structures with soft anions (e.g., S 2- and Se 2- ), they do not meet energy density requirements to compete with the current rechargeable lithium-ion batteries due to their low insertion potentials, emphasizing the necessity of finding an oxide-based cathode that operates at high potentials. A leading hypothesis to explain the limited availability of oxide Mg-ion cathodes is the belief that Mg 2+ has sluggish diffusion kinetics in oxides due to strong electrostatic interactions between the Mg 2+ ions and oxide anions in the lattice. From this assessment, it can be hypothesized that such rate limiting kinetic shortcomings can be mitigated by tailoring an oxide framework through creating less stable Mg 2+ - O 2- coordination.

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

Assessment of exchange-correlation functionals on oxygen vacancy formation energies of metal oxides

Oxygen vacancies are a common defect in oxides and play vital roles in many technological applications such as oxygen separation, catalytic reactors, solid oxide fuel cells, and solar thermochemical water splitting. The oxygen vacancy formation energy is directly related to the oxide reduction enthalpy and the vacancy concentration and, therefore, is a key quantity underlying these applications. The complexity of measuring oxygen vacancy formation energies experimentally suggests the utility of calculations based on density functional theory (DFT). However, the calculated results are strongly dependent on the exchange-correlation functionals and other parameters such as the Hubbard U . In this work, we compare the performance of the meta-GGA strongly constrained and appropriately normed (SCAN) functional and the commonly used semilocal generalized gradient approximation (GGA) functionals with experimental values of structural parameters, band gaps, magnetic structures, and oxygen vacancy formation energies for six representative oxides, i.e., CaMnO 3 , SrMnO 3 , LaMnO 3 , YMnO 3 , LaFeO 3 , and CeO 2 . Further, our results show that SCAN usually has better agreement with the experimental lattice constants and band gaps and larger magnetic moments than the commonly used GGA functionals. Although SCAN overestimates the oxygen vacancy formation energies of transition metal oxides and therefore requires unusually large Hubbard U values to reproduce the experimental reduction enthalpies, it does predict the correct oxygen vacancy formation energy of CeO 2 , which challenges the commonly used GGA and hybrid functionals. Our results underscore the challenges that exist in describing these complex oxides by DFT and may shed light on developing more accurate exchange-correlation functionals for oxygen vacancy formation energy calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Asymmetric equilibrium core structures of pyramidal-II < c + a > dislocations in ten hexagonal-close-packed metals

The structures of pyramidal-II < c + a > dislocations, one of the most important defects in structural hexagonal-close-packed (HCP) metals, have not been fully characterized for many of the HCP metals in use today. Here, we employ ab initio informed phase-field dislocation dynamics to determine the minimum energy structure of pyramidal {1¯1¯22} < 11¯23 > dislocations in ten HCP metals, including Be, Co, Mg, Re, Ti, Zn, Cd, Hf, Y, and Zr. As input for the simulations, we calculate, using first-principles density functional theory, the {1¯1¯22} generalized stacking fault energy (GSFE) curves for all ten metals. From these calculations, it is found that magnetism in Co is necessary for achieving a local minimum in the GSFE curve. We observe in simulations that edge and screw character dislocations split into two partials separated by a low-energy intrinsic stacking fault. The splitting distance is shown to scale inversely with the local minimum energy normalized by the product of its shear modulus and Burgers vector. Interestingly, some HCP metals exhibit an asymmetric structure, with either unequal partial Burgers vectors or widths, in contrast to the symmetric configuration expected from linear elastic dislocation theory. We explain these structures by properties of the local maxima in their GSFE curves. Metals with larger degrees of elastic anisotropy result in dislocations with larger splitting distances than would be expected under the commonly used assumption of elastic isotropy. Furthermore, these findings on the sizes and asymmetry in the structures of pyramidal-II < c + a > dislocations are fundamental to understanding how these dislocations glide and interact or react with other defects when these metals are mechanically strained.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Cosmological decoherence from thermal gravitons

A bstract We study the effects of gravitationally-driven decoherence on tunneling processes associated with false vacuum decays, such as the Coleman-De Luccia instanton. We compute the thermal graviton-induced decoherence rate for a wave function describing a perfect fluid of nonzero energy density in a finite region. When the effective cosmological constant is positive, the thermal graviton background sourced by a de Sitter horizon provides an unavoidable decoherence effect, which may have important consequences for tunneling processes in cosmological history. We discuss generalizations and consequences of this effect and comment on its observability and applications to black hole physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Application of machine learning interatomic potentials in heterogeneous catalysis

Heterogeneous catalysts are crucial in modern societies as they promote sustainability by enabling lower-energy pathways for various chemical reactions. While Density Functional Theory (DFT) computations can provide critical insights into how heterogeneous catalysts operate at the atomic level, they are limited by computational costs and unfavorable scaling with system size. Recently, machine learning interatomic potentials (MLIPs) have emerged as a promising alternative to DFT, offering near-DFT accuracy at significantly reduced cost. Here, in this perspective, we discuss the application of MLIPs in heterogeneous catalyst modeling as a surrogate for DFT. We detail how MLIPs have been applied in thermal catalysis to probe active sites, enable studying complex metallic and nanoporous catalysts, and investigate the reconstruction of catalytic surfaces. We review the use of MLIPs in electrocatalysis and photocatalysis, emphasizing their capabilities in studying transition metal oxide surfaces and solid–liquid interfaces. We also discuss the current limitations of MLIPs, particularly their challenges with transferability and description of non-local interactions. Finally, we conclude by identifying promising and underexplored domains in which MLIPs can further advance our understanding of heterogeneous catalysts.

Catalytic surfaces↗
DOE OSTI2021

Effect of static local distortions vs. dynamic motions on the stability and band gaps of cubic oxide and halide perovskites

Ternary ABX 3 perovskites made of corner-sharing BX 6 octahedra have long featured prominently in solid-state chemistry and condensed matter physics. Still, the joint understanding of their two main subgroups—halides and oxides—has not been fully developed. Indeed, unlike the case in simpler compounds having a single, robust repeated motif (“monomorphous”), certain cubic perovskites can manifest a non-thermal (= intrinsic) distribution of local motifs (“polymorphous networks”). Such static deformations can include positional degrees of freedom (e.g., atomic displacements and octahedral tilting) or magnetic moment degrees of freedom in paramagnets. Unlike thermal motion, such static distortions do not time-average to zero, being an expression of the intrinsic symmetry breaking preference of the chemical bonding. Here, the present study compares electronic structure features of oxide and halide perovskites starting from the static polymorphous distribution of motifs described by Density Functional Theory (DFT) minimization of the internal energy, continuing to finite temperature thermal disorder modeled via finite temperature DFT molecular dynamics. We find that (i) different oxide vs. halide ABX 3 compounds adopt different energy-lowering symmetry-breaking modes. The calculated pair distribution function (PDF) of SrTiO 3 from the first-principles agrees with recently measured PDF. (ii) In both oxides and halides, such static distortions lead to band gap blueshifts with respect to undistorted cubic Pm-3m structure. (iii) For oxide perovskites, high-temperature molecular dynamics simulations initiated from the statically distorted polymorphous structures reveal that the thermally-induced distortions can lead to a band gap redshift. (iv) In contrast, for cubic halide perovskite CsPbI 3 , both the intrinsic distortions and the thermal distortions contribute in tandem to band gap blueshift, the former, intrinsic effect being dominant. (v) In the oxide SrTiO 3 and CaTiO 3 (but not in halide) perovskites, octahedral tilting leads to the emergence of a distinct Γ–Γ direct band gap component as a secondary valley minimum to the well-known indirect R–Γ gap. Understanding such intrinsic vs. thermal effects on oxide vs. halide perovskites holds the potential for designing target electronic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

MnRhBi 3 : A Cleavable Antiferromagnetic Metal

Cleavable metallic antiferromagnets may be of use for low-dissipation spintronic devices; however, few are currently known. Here, in this study, we present orthorhombic MnRhBi 3 as one such compound and present a thorough study of its physical properties. Exfoliation is demonstrated experimentally, and the cleavage energy and electronic structure are examined by density functional theory calculations. It is concluded that MnRhBi 3 is a van der Waals-layered material that cleaves easily between neighboring Bi layers and that the Bi atoms have lone pairs extending into the van der Waals gaps. A series of four phase transitions are observed below room temperature, and neutron diffraction shows that at least two of the transitions involve the formation of antiferromagnetic order. Anomalous thermal expansion points to a crystallographic phase transition and/or strong magnetoelastic coupling. This work reveals a complex phase evolution in MnRhBi 3 and establishes this cleavable antiferromagnetic metal as an interesting material for studying the interplay of structure, magnetism, and transport in the bulk and ultrathin limits, as well as the role of lone pair electrons in interface chemistry and proximity effects in van der Waals heterostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Vibrational Spectroscopy of Hexahalo Complexes

Halogenated inorganic complexes A x [MHal y ] (A = alkali metal or alkaline earth, M = transition or main group metal, x = 1–3, and y = 2–9) are an archetypal class of compounds that provide entry points to large areas of inorganic and physical chemistry. All of the hexahalo complexes adopt an octahedral, O h , symmetry (or nearly so). Consequently, one of the bending modes is forbidden in both the infrared and Raman spectra. In the solid state, many of the complexes crystallize in the cubic space group Fm3̅m, which preserves the octahedral symmetry. Even for those that are not cubic, the octahedral symmetry of the [MHal 6 ] n- ion is largely retained and, to a good approximation, so are the selection rules. In the present work, we show that by using the additional information provided by neutron vibrational spectroscopy, in combination with conventional optical spectroscopies, we can generate complete and unambiguous assignments for all the modes. Comparison of the experimental and calculated transition energies for the systems where periodic-density functional theory was possible (i.e., those for which the crystal structure is known) shows that the agreement is almost quantitative. We also provide a linear relationship that enables the prediction of the forbidden mode.

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

Gas-Phase Complexes of Americium and Lanthanides with a Bis-triazinyl Pyridine: Reactivity and Bonding of Archetypes for F-Element Separations

Bis-triazinyl pyridines (BTPs) exhibit solution selectivity for trivalent americium over lanthanides (Ln), the origins of which remain uncertain. Here in this paper, electrospray ionization was used to generate gas-phase complexes [ML 3 ] 3+ , where M = La, Lu, or Am and L is EtBTP 2,6-bis(5,6-diethyl-1,2,4-triazin-3-yl)-pyridine. Collision-induced dissociation (CID) of [ML 3 ] 3+ in the presence of H 2 O yielded a protonated ligand [L(H)] + and hydroxide [ML 2 (OH)] 2+ or hydrate [ML(L–H)(H 2 O)] 2+ , where (L–H) - is a deprotonated ligand. Although solution affinities indicate stronger binding of BTPs toward Am 3+ versus Ln 3+ , the observed CID process is contrastingly more facile for M = Am versus Ln. To understand the disparity, density functional theory was employed to compute potential energy surfaces for two possible CID processes, for M = La and Am. In accordance with the CID results, both the rate determining transition state barrier and the net energy are lower for [AmL 3 ] 3+ versus [LaL 3 ] 3+ and for both product isomers, [ML 2 (OH)] 2+ and [ML(L–H)(H 2 O)] 2+ . More facile removal of a ligand from [AmL 3 ] 3+ by CID does not necessarily contradict stronger Am 3+ –L binding, as inferred from solution behavior. In particular, the formation of new bonds in the products can distort kinetics and thermodynamics expected for simple bond cleavage reactions. In addition to correctly predicting the seemingly anomalous CID behavior, the computational results indicate greater participation of Am 5f versus La 4f orbitals in metal–ligand bonding.

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

Small Molecule Inhibitor-Modulated Al 2 O 3 Atomic Layer Deposition on Monolayer MoS 2 for Controlled Nucleation

The integration of ultrathin dielectrics on two-dimensional (2D) semiconductors is essential for advancing beyond-Si electronics. However, the intrinsic inertness of van der Waals 2D basal planes remains a primary bottleneck to achieving uniform dielectric nucleation and growth. Here, in this study, we introduce a small molecule inhibitor (SMI)-modulated thermal atomic layer deposition (ALD) strategy, exemplified by aluminum oxide (Al 2 O 3 ) ALD on monolayer molybdenum disulfide (1L MoS 2 ) with acetic acid (HAc) SMI. The ABC-type sequence comprises HAc inhibitor (A), trimethylaluminum (TMA) precursor (B), and deionized H 2 O coreactant (C). In situ quartz crystal microbalance (QCM) studies reveal robust HAc adsorption on Al 2 O 3 and suppression of subsequent oxide growth on HAc-passivated surfaces. When applied to 1L MoS 2 , this inhibitory pathway enables HAc to selectively passivate nascent Al 2 O 3 nuclei formed on the MoS 2 surface, limiting their three-dimensional (3D) island coarsening and redirecting precursor adsorption toward the uncovered basal plane. Consequently, nearly continuous ultrathin (∼1.5 nm) Al 2 O 3 films are achieved on 1L MoS 2 with markedly improved uniformity compared to standard Al 2 O 3 ALD using TMA and H 2 O, as validated by atomic force microscopy (AFM), cross-sectional scanning transmission electron microscopy (STEM), and energy-dispersive X-ray spectroscopy (EDS). Density functional theory (DFT) calculations further provide atomistic insight into HAc-modulated Al 2 O 3 nucleation, corroborating the energetic preference of HAc for Al 2 O 3 over MoS 2 and attenuated TMA adsorption on HAc-passivated surfaces. Spatially resolved Raman spectroscopy also confirms that the HAc-modulated process preserves the structural integrity of 1L MoS 2 , with only minimal strain and doping perturbations observed after dielectric deposition. This SMI-modulated approach offers a broadly applicable framework for controlling ALD nucleation across various inhibitors, ALD chemistries, and 2D materials, opening opportunities for reliable dielectric integration in next-generation nanoelectronics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Off-Equilibrium Reactivity of Boron-Enriched Metal Diboride Surfaces in Electroreduction Conditions

Boron-based materials, featuring B-dependent reactivity and diverse phases, are emerging as promising catalyst systems. However, the catalytic mechanism on many borides remains poorly understood due to complex surface reconstructions under reaction conditions. Here, we investigate the MoB 2 surface in conditions of hydrogen evolution reaction in acidic media, using grand canonical global optimization, grand canonical density functional theory, ab initio molecular dynamics, free energy surface sampling, and an analytical model for electrochemical barrier evaluation. We propose a boron-enrichment strategy to tune the surface reactivity of the hexagonal face of MoB 2 . We reveal the dynamic nature of the B-enriched surface under H coverage and kinetic trapping of the system in the metastable regime with an extensive examination of the deactivation pathways. The metastable center B site on B-enriched surfaces, featuring buckled-up configuration and a usual relaxation effect, is found to be highly active toward HER via the Volmer–Heyrovsky mechanism. In conclusion, this work demonstrates how off-equilibrium behaviors can arise from the interplay between adsorbate coverage and surface reconstruction on a seemingly simple surface, and we present a theoretical framework and computational workflows to address these behaviors, along with other realistic complexities, in kinetics simulations.

Adsorption↗
DOE OSTI2024

Molecular-scale Insights into Cooperativity Switching of x TAB Adsorption on Gold Nanoparticles

Quantifying adsorption behaviors is crucial for various applications such as catalysis, separation, and sensing, yet it is generally challenging to access in solution. Here, we report a combined experimental and computational study of the adsorption behaviors of alkyl-trimethylammonium bromides (xTAB), a class of ligands important for colloidal nanoparticle stabilization and shape control, with various alkyl chain lengths x on Au nanoparticles. We use density functional theory (DFT) to calculate xTAB binding energies on Au{111} and Au{110} surfaces with standing-up and lying-down configurations, which provides insights into the adsorption affinity and cooperativity differences of xTAB on these two facets. We demonstrate the key role of van der Waals interactions in determining the xTAB adsorption behavior. These computational results predict and explain the experimental discovery of xTAB’s adsorption behavior switch from stronger affinity, negative cooperativity to weaker affinity, positive cooperativity when the concentration of xTAB increases in solution. We also show that in the standing-up configuration, bilayer adsorption may occur on both facets, which can lead to different differential binding energies and consequently adsorption crossover between the two facets when the ligand concentration increases. Our combined experimental and computational approaches demonstrate a paradigm for gaining molecular-scale insights into adsorbate–surface interactions.

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

Stabilized Synthesis of 2D Verbeekite: Monoclinic PdSe 2 Crystals with High Mobility and In-Plane Optical and Electrical Anisotropy

PdSe 2 has a layered structure with an unusual, puckered Cairo pentagonal tiling. Its atomic bond configuration features planar 4-fold-coordinated Pd atoms and intralayer Se–Se bonds that enable polymorphic phases with distinct electronic and quantum properties, especially when atomically thin. PdSe 2 is conventionally orthorhombic, and direct synthesis of its metastable polymorphic phases is still a challenge. Here, we report an ambient-pressure chemical vapor deposition approach to synthesize metastable monoclinic PdSe 2 . Monoclinic PdSe 2 is shown to be synthesized selectively under Se-deficient conditions that induce Se vacancies. These defects are shown by first-principles density functional theory calculations to reduce the free energy of the metastable monoclinic phase, thereby stabilizing it during synthesis. Further, the structure and composition of the monoclinic PdSe 2 crystals are identified and characterized by scanning transmission electron microscopy imaging, convergent beam electron diffraction, and electron energy loss spectroscopy. Polarized Raman spectroscopy of the monoclinic PdSe 2 flakes reveals their strong in-plane optical anisotropy. Electrical transport measurements show that the monoclinic PdSe 2 exhibits n-type charge carrier conduction with electron mobilities up to ~298 cm 2 V –1 s –1 and a strong in-plane electron mobility anisotropy of ~1.9. The defect-mediated growth pathway identified in this work is promising for phase-selective direct synthesis of other 2D transition metal dichalcogenides.

2D material↗