-
The first-order structural phase transition at low-temperature in GaPt$_{5}$P and its rapid enhancement with pressure
Authors:
A. Sapkota,
T. J. Slade,
S. Huyan,
N. K. Nepal,
J. M. Wilde,
N. Furukawa,
S. H. Laupidus,
L. -L. Wang,
S. L. Bud'ko,
P. C. Canfield
Abstract:
Single crystals of XPt$_{5}$P (X = Al, Ga, and In) were grown from a Pt-P solution at high temperatures, and ambient-pressure measurements of temperature-dependent magnetization, resistivity, and X-ray diffraction were made. Also, the ambient-pressure Hall resistivity and temperature-dependent resistance under pressure were measured on GaPt$_{5}$P. All three compounds have tetragonal $P4/mmm$ crys…
▽ More
Single crystals of XPt$_{5}$P (X = Al, Ga, and In) were grown from a Pt-P solution at high temperatures, and ambient-pressure measurements of temperature-dependent magnetization, resistivity, and X-ray diffraction were made. Also, the ambient-pressure Hall resistivity and temperature-dependent resistance under pressure were measured on GaPt$_{5}$P. All three compounds have tetragonal $P4/mmm$ crystal structure at room-temperature with metallic transport and weak diamagnetism over the $2-300$~K temperature range. Surprisingly, at ambient pressure, both the transport and magnetization measurements on GaPt$_{5}$P show a step-like feature in $70-90$~K region suggesting a possible structural phase transition, and no such features were observed in (Al/In)Pt$_{5}$P. Both the hysteretic nature and sharpness of the feature suggest the first-order transition, and single-crystal X-ray diffraction measurements provided further details of the structural transition with a crystal symmetry likely different than $P4/mmm$ below transition. The transition is characterized by anisotropic changes in the lattice parameters, a volume collapse, and satellite peaks at two distinct wave-vectors. Density functional theory calculations present phonon softening as a possible driving mechanism. Additionally, the structural transition temperature increases rapidly with increasing pressure, reaching room temperature by $\sim 2.2$~GPa, highlighting the high degree of pressure sensitivity and fragile nature of GaPt$_{5}$P room-temperature structure. Although the volume collapse and extreme pressure sensitivity suggest chemical pressure should drive a similar structural change in AlPt$_{5}$P, with smaller unit cell dimensions and volume, its structure is found to be $P4/mmm$ as well. Overall, GaPt$_{5}$P stands out as a sole member of the 1-5-1 family of compounds with a temperature-driven structural change.
△ Less
Submitted 10 June, 2024;
originally announced June 2024.
-
Imaginary phonon modes and phonon-mediated superconductivity in Y2C3
Authors:
Niraj K. Nepal,
Paul C. Canfield,
Lin-Lin Wang
Abstract:
For Y$_2$C$_3$ with a superconducting critical temperature (T$_c$) $\sim$18 K, zone-center imaginary optical phonon modes have been found for the high-symmetry $I$-$43d$ structure due to C dimer wobbling motion and electronic instability from a flat band near Fermi energy. After lattice distortion to the more stable lowest symmetry $P1$ structure, these stabilized low-energy phonon modes with a mi…
▽ More
For Y$_2$C$_3$ with a superconducting critical temperature (T$_c$) $\sim$18 K, zone-center imaginary optical phonon modes have been found for the high-symmetry $I$-$43d$ structure due to C dimer wobbling motion and electronic instability from a flat band near Fermi energy. After lattice distortion to the more stable lowest symmetry $P1$ structure, these stabilized low-energy phonon modes with a mixed C and Y character carry a strong electron-phonon coupling to give arise to the observed sizable T$_c$. Our work shows that compounds with the calculated dynamical instability should not be simply excluded in high-throughput search for new phonon-mediated superconductors. Moreover, we have studied the phase stability of the $I$-$43d$ structure by calculating the enthalpy of different structural motifs of binary compounds containing group IV elements at the 2:3 composition and also exploring the energy landscapes via $ab$ $initio$ molecular dynamics near and out of the $I$-$43d$ structure. Our results show that the $I$-$43d$ type structures with C dimers are preferred in the low to medium pressure range. Because of the wobbling motion of the C dimers, there are many local energy minimums with degenerated energies. Thus, the ensemble average of many $I$-$43d$-distorted structures with C dimer wobbling motion at finite temperature still gives an overall $I$-$43d$ structure.
△ Less
Submitted 22 February, 2024; v1 submitted 31 July, 2023;
originally announced August 2023.
-
Origin of Charge Density Wave in Topological Semimetals SrAl4 and EuAl4
Authors:
Lin-Lin Wang,
Niraj K. Nepal,
Paul C. Canfield
Abstract:
Topological semimetals in BaAl4-type structure show many interesting behaviors, such as charge density wave (CDW) in SrAl4 and EuAl4, but not the isostructural and isovalent BaAl4, SrGa4 and BaGa4. Here using Wannier functions based on density functional theory, we calculate the susceptibility functions with millions of k-points to reach the small q-vector and study the origin and driving force be…
▽ More
Topological semimetals in BaAl4-type structure show many interesting behaviors, such as charge density wave (CDW) in SrAl4 and EuAl4, but not the isostructural and isovalent BaAl4, SrGa4 and BaGa4. Here using Wannier functions based on density functional theory, we calculate the susceptibility functions with millions of k-points to reach the small q-vector and study the origin and driving force behind the CDW. Our comparative study reveals that the origin of the CDW in SrAl4 and EuAl4 is the strong electron-phonon coupling interaction for the transverse acoustic mode at small q-vector along the \{Gamma}-Z direction besides the maximum of the real part of the susceptibility function from the nested Fermi surfaces of the Dirac-like bands, which explains well the absence of CDW in the other closely related compounds in a good agreement with experiment. We also connect the different CDW behaviors in the Al compounds to the macroscopic elastic properties.
△ Less
Submitted 29 March, 2024; v1 submitted 26 June, 2023;
originally announced June 2023.
-
Bending as a control knob for the electronic and optical properties of phosphorene nanoribbons
Authors:
Bimal Neupane,
Hong Tang,
Niraj K. Nepal,
Adrienn Ruzsinszky
Abstract:
We have assessed mechanical bending as a powerful controlling tool for the electronic structure and optical properties of phosphorene nanoribbons. We use state-of-the-art density functional approximations in our work. The overall performance of the recently developed meta-generalized gradient approximation (meta-GGA) mTASK [Phys. Rev. Mater. \textbf{5}, 063803 (2021)] functional establishes the me…
▽ More
We have assessed mechanical bending as a powerful controlling tool for the electronic structure and optical properties of phosphorene nanoribbons. We use state-of-the-art density functional approximations in our work. The overall performance of the recently developed meta-generalized gradient approximation (meta-GGA) mTASK [Phys. Rev. Mater. \textbf{5}, 063803 (2021)] functional establishes the method as a useful alternative to the screened hybrid HSE06 for better band gaps of phosphorene nanoribbons. We present a detailed and novel analysis to interpret the optical absorption of bent phosphorene nanoribbons using the GW-Bethe-Salpeter approximation (GW-BSE). We demonstrate the important role of the unoccupied in-gap state and conclude that this in-gap state in armchair nanoribbons introduced by bending can significantly affect the properties of low-energy excitons, and add some useful opportunities for applications as optoelectronic devices.
△ Less
Submitted 7 December, 2021;
originally announced December 2021.
-
Tunable band gaps and optical absorption properties of bent MoS$_2$ nanoribbons
Authors:
Hong Tang,
Bimal Neupane,
Santosh Neupane,
Shiqi Ruan,
Niraj K. Nepal,
Adrienn Ruzsinszky
Abstract:
The large tunability of band gaps and optical absorptions of armchair MoS$_2$ nanoribbons of different widths under bending is studied using density functional theory and many-body perturbation GW and Bethe-Salpeter equation approaches. We find that there are two critical bending curvatures, and the non-edge and edge band gaps generally show a non-monotonic trend with bending. The non-degenerate e…
▽ More
The large tunability of band gaps and optical absorptions of armchair MoS$_2$ nanoribbons of different widths under bending is studied using density functional theory and many-body perturbation GW and Bethe-Salpeter equation approaches. We find that there are two critical bending curvatures, and the non-edge and edge band gaps generally show a non-monotonic trend with bending. The non-degenerate edge gap splits show an oscillating feature with ribbon width n, with a period delta_n=3, due to quantum confinement effects. The complex strain patterns on the bent nanoribbons control the varying features of band structures and band gaps that result in varying exciton formations and optical properties. The binding energy and the spin singlet-triplet split of the exciton forming the lowest absorption peak generally decrease with bending curvatures. The large tunability of optical properties of bent MoS$_2$ nanoribbons is promising and will find applications in tunable optoelectronic nanodevices.
△ Less
Submitted 31 October, 2022; v1 submitted 9 September, 2021;
originally announced September 2021.
-
First-principles wavevector- and frequency-dependent exchange-correlation kernel for jellium at all densities
Authors:
Aaron D. Kaplan,
Niraj K. Nepal,
Adrienn Ruzsinszky,
Pietro Ballone,
John P. Perdew
Abstract:
We propose a spatially and temporally nonlocal exchange-correlation (xc) kernel for the spin-unpolarized fluid phase of ground-state jellium, for use in time-dependent density functional and linear response calculations. The kernel is constructed to satisfy known properties of the exact xc kernel, to accurately describe the correlation energies of bulk jellium, and to satisfy frequency-moment sum…
▽ More
We propose a spatially and temporally nonlocal exchange-correlation (xc) kernel for the spin-unpolarized fluid phase of ground-state jellium, for use in time-dependent density functional and linear response calculations. The kernel is constructed to satisfy known properties of the exact xc kernel, to accurately describe the correlation energies of bulk jellium, and to satisfy frequency-moment sum rules at a wide range of bulk jellium densities, including those low densities that display strong correlation and symmetry breaking. These effects are easier to understand in the simple jellium model than in real systems. All exact constraints satisfied by the recent MCP07 kernel [A. Ruzsinszky, et al., Phys. Rev. B 101, 245135 (2020)] are maintained in the new revised MCP07 (rMCP07) kernel, while others are added. The revision $f_\mathrm{xc}^\mathrm{rMCP07}(q,ω)$ differs from MCP07 only for non-zero frequencies $ω$. Only at densities much lower than those of real bulk metals is the frequency dependence of the kernel important for the correlation energy of jellium. As the wavevector $q$ tends to zero, the kernel has a $-4πα(ω)/q^2$ divergence whose frequency-dependent ultranonlocality coefficient $α(ω)$ vanishes in jellium, and is predicted by rMCP07 to be extremely small for the real metals Al and Na.}
△ Less
Submitted 22 November, 2021; v1 submitted 5 July, 2021;
originally announced July 2021.
-
Progress towards understanding ultranonlocality through the wavevector and frequency dependence of approximate exchange-correlation kernels
Authors:
Niraj K. Nepal,
Aaron D. Kaplan,
J. M. Pitarke,
Adrienn Ruzsinszky
Abstract:
In the framework of time-dependent density functional theory (TDDFT), the exact exchange-correlation (xc) kernel $f_{xc}(n,q,ω)$ determines the ground-state energy, excited-state energies, lifetimes, and the time-dependent linear density response of any many-electron system. The recently developed MCP07 xc kernel $f_{xc}(n,q,ω)$ of A. Ruzsinszky et al. [Phys. Rev. B 101, 245135 (2020)] yields exce…
▽ More
In the framework of time-dependent density functional theory (TDDFT), the exact exchange-correlation (xc) kernel $f_{xc}(n,q,ω)$ determines the ground-state energy, excited-state energies, lifetimes, and the time-dependent linear density response of any many-electron system. The recently developed MCP07 xc kernel $f_{xc}(n,q,ω)$ of A. Ruzsinszky et al. [Phys. Rev. B 101, 245135 (2020)] yields excellent uniform electron gas (UEG) ground-state energies and plausible plasmon lifetimes. As MCP07 is constructed to describe $f_{xc}$ of the UEG, it cannot capture optical properties of real materials. To verify this claim, we follow Nazarov et al. [Phys. Rev. Lett. 102, 113001 (2009)] to construct the long-range, dynamic xc kernel, $\lim_{q\to 0}f_{xc}(n,q,ω) = -α(ω)e^2/q^2$, of a weakly inhomogeneous electron gas, using MCP07 and other common xc kernels. The strong wavevector and frequency dependence of the "ultranonlocality" coefficient $α(ω)$ is demonstrated for a variety of simple metals and semiconductors. We examine how imposing exact constraints on an approximate kernel shapes $α(ω)$. Comparisons to kernels derived from correlated-wavefunction calculations are drawn.
△ Less
Submitted 30 July, 2021; v1 submitted 24 May, 2021;
originally announced May 2021.
-
Opening band gaps of low-dimensional materials at the meta-GGA level of density functional approximations
Authors:
Bimal Neupane,
Hong Tang,
Niraj K. Nepal,
Santosh Adhikari,
Adrienn Ruzsinszky
Abstract:
The recent TASK meta-GGA density functional [Phys. Rev. Research, 1, 033082 (2019)] is constructed with an enhanced nonlocality in the generalized Kohn-Sham scheme, and therefore harbors great opportunities for band gap prediction. Although this approximation was found to yield excellent band gaps of bulk solids, this accuracy cannot be straightforwardly transferred to low-dimensional materials. T…
▽ More
The recent TASK meta-GGA density functional [Phys. Rev. Research, 1, 033082 (2019)] is constructed with an enhanced nonlocality in the generalized Kohn-Sham scheme, and therefore harbors great opportunities for band gap prediction. Although this approximation was found to yield excellent band gaps of bulk solids, this accuracy cannot be straightforwardly transferred to low-dimensional materials. The reduced screening of these materials results in larger band gaps compared to their bulk counterparts, as an additional barrier to overcome. In this work, we demonstrate how the alteration of exact physical constraints in this functional affects the band gaps of monolayers and nanoribbons, and present accurate band gaps competing with the HSE06 approximation. In order to achieve this goal, we have modified the TASK functional (a) by changing the tight upper-bound for one or two-electron systems ($h_X^0$) from 1.174 to 1.29 (b) by changing the limit of interpolation function $f_X (α\rightarrow \infty$) of the TASK functional that interpolates the exchange enhancement factor $F_X (s,α)$ from $α=$ 0 to 1. The resulting modified TASK (mTASK) was tested for various materials from 3D to 2D to 1D (nanoribbons), and was compared with the results of the higher-level hybrid functional HSE06 or with the G$_0$W$_0$ approximation within many-body perturbation theory. We find that mTASK greatly improves the band gaps and band structures of 2D and 1D systems, without significantly affecting the accuracy of the original TASK for the bulk 3D materials, when compared to the PBE-GGA and SCAN meta-GGA. We further demonstrate the applicability of mTASK by assessing the band structures of TMD nanoribbons with respect to various bending curvatures.
△ Less
Submitted 8 April, 2021;
originally announced April 2021.
-
Resolving the structure-energy dilemma at organic-inorganic interfaces: Adsorption of benzene, thiophene, and xenon over coinage metal surfaces
Authors:
Santosh Adhikari,
Niraj K. Nepal,
Hong Tang,
Adrienn Ruzsinszky
Abstract:
Semilocal (SL) density functional approximations (DFAs) are widely applied but have limitations due to their inability to incorporate long-range van der Waals (vdW) interaction. Non-local functionals (vdW-DF, VV10, rVV10) or empirical methods (DFT+D, DFT+vdW, DFT+MBD) are used with SL-DFAs to account for such missing interaction. The physisorption of a molecule on the surface of the coinage metals…
▽ More
Semilocal (SL) density functional approximations (DFAs) are widely applied but have limitations due to their inability to incorporate long-range van der Waals (vdW) interaction. Non-local functionals (vdW-DF, VV10, rVV10) or empirical methods (DFT+D, DFT+vdW, DFT+MBD) are used with SL-DFAs to account for such missing interaction. The physisorption of a molecule on the surface of the coinage metals (Cu, Ag, and Au) is a typical example of systems where vdW interaction is significant. However, it is difficult to find a general method that reasonably describes both adsorption energy and geometry of even the simple prototypes of cyclic and heterocyclic aromatic molecules like benzene (C6H6) and thiophene (C4H4S) respectively, with reasonable accuracy. In this work, we present an alternative scheme based on Zaremba-Kohn's theory, called DFT+vdW-dZK. We show that, unlike other popular methods, DFT+vdW-dZK and particularly SCAN+vdW-dZK gives an accurate description of the physisorption of a rare-gas atom (Xe) and two small albeit diverse prototype organic molecules on the (111) surfaces of the coinage metals.
△ Less
Submitted 3 January, 2021;
originally announced January 2021.
-
The Fermi-Löwdin self-interaction correction for ionization energies of organic molecules
Authors:
Santosh Adhikari,
Biswajit Santra,
Shiqi Ruan,
Puskar Bhattarai,
Niraj K. Nepal,
Koblar A. Jackson,
Adrienn Ruzsinszky
Abstract:
(Semi)-local density functional approximations (DFAs) suffer from self-interaction error (SIE). When the first ionization energy (IE) is computed as the negative of the highest-occupied orbital (HO) eigenvalue, DFAs notoriously underestimate them compared to quasi-particle calculations. The inaccuracy for the HO is attributed to SIE inherent in DFAs. We assessed the IE based on Perdew-Zunger self-…
▽ More
(Semi)-local density functional approximations (DFAs) suffer from self-interaction error (SIE). When the first ionization energy (IE) is computed as the negative of the highest-occupied orbital (HO) eigenvalue, DFAs notoriously underestimate them compared to quasi-particle calculations. The inaccuracy for the HO is attributed to SIE inherent in DFAs. We assessed the IE based on Perdew-Zunger self-interaction corrections on 14 small to moderate-sized organic molecules relevant in organic electronics and polymer donor materials. Though self-interaction corrected DFAs were found to significantly improve the IE relative to the uncorrected DFAs, they overestimate. However, when the self-interaction correction is interiorly scaled using a function of the iso-orbital indicator zσ, only the regions where SIE is significant get a correction. We discuss these approaches and show how these methods significantly improve the description of the HO eigenvalue for the organic molecules.
△ Less
Submitted 11 August, 2020;
originally announced August 2020.
-
Interpretations of ground-state symmetry breaking and strong correlation in wavefunction and density functional theories
Authors:
John P. Perdew,
Adrienn Ruzsinszky,
Jianwei Sun,
Niraj K. Nepal,
Aaron D. Kaplan
Abstract:
Strong correlations within a symmetry-unbroken ground-state wavefunction can show up in approximate density functional theory as symmetry-broken spin-densities or total densities, which are sometimes observable. They can arise from soft modes of fluctuations (sometimes collective excitations) such as spin-density or charge-density waves at non-zero wavevector. In this sense, an approximate density…
▽ More
Strong correlations within a symmetry-unbroken ground-state wavefunction can show up in approximate density functional theory as symmetry-broken spin-densities or total densities, which are sometimes observable. They can arise from soft modes of fluctuations (sometimes collective excitations) such as spin-density or charge-density waves at non-zero wavevector. In this sense, an approximate density functional for exchange and correlation that breaks symmetry can be more revealing (albeit less accurate) than an exact functional that does not. The examples discussed here include the stretched H$_2$ molecule, antiferromagnetic solids, and the static charge-density wave/Wigner crystal phase of a low-density jellium. It is shown that (and in what sense) the static charge density wave is a soft plasmon.
△ Less
Submitted 23 July, 2020;
originally announced July 2020.
-
Formation energy puzzle in intermetallic alloys: Random phase approximation fails to predict accurate formation energies
Authors:
Niraj K. Nepal,
Santosh Adhikari,
Bimal Neupane,
Adrienn Ruzsinszky
Abstract:
We performed density functional calculations to estimate the formation energies of intermetallic alloys. We used two semilocal approximations, the generalized gradient approximation (GGA) by Perdew-Burke-Ernzerhof (PBE) and the strongly constrained and appropriately normed (SCAN) meta-GGA. In addition, we utilized two nonlocal DFT functionals, the hybrid HSE06, and the state-of-the-art random phas…
▽ More
We performed density functional calculations to estimate the formation energies of intermetallic alloys. We used two semilocal approximations, the generalized gradient approximation (GGA) by Perdew-Burke-Ernzerhof (PBE) and the strongly constrained and appropriately normed (SCAN) meta-GGA. In addition, we utilized two nonlocal DFT functionals, the hybrid HSE06, and the state-of-the-art random phase approximation (RPA). The nonlocal functionals such as HSE06 and RPA yield accurate formation energies of binary alloys with completely-filled d-band metals, where semilocal functionals underperform. The accuracy at the nonlocal functionals is greatly reduced when a partially-filled d-band metal is present in an alloy, while PBE-GGA outperforms in these cases. We show that the accurate prediction of formation energies by any DFT method depends on its ability to predict the accurate electronic properties, e.g., valence d-band contribution to the density of states (DOS). The SCAN meta-GGA often corrects the PBE-DOS, however, it does not provide accurate formation energies compared to PBE. This is assumed to be due to the lack of proper error cancellation that should be expected due to the similar bulk nature of both alloys and their constituents, which may improve with the modification of meta-GGA ingredients. RPA yields too negative formation energies of alloys with partially-filled d-band metals. RPA results can be corrected by restoring the exchange-correlation kernel, thereby improving the short-range electron-electron correlation in metallic densities.
△ Less
Submitted 1 October, 2020; v1 submitted 21 July, 2020;
originally announced July 2020.
-
Understanding plasmon dispersion in nearly-free-electron metals: the relevance of exact constraints for novel exchange-correlation kernels within time-dependent density functional theory
Authors:
Niraj K. Nepal,
Santosh Adhikari,
Bimal Neupane,
Shiqi Ruan,
Santosh Neupane,
Adrienn Ruzsinszky
Abstract:
Small-wavevector excitations in Coulomb-interacting systems can be decomposed into the high-energy collective longitudinal plasmon and the low-energy single-electron excitations. At the critical wavevector and corresponding frequency where the plasmon branch merges with the single-electron excitation region, the collective energy of the plasmon dissipates into single electron-hole excitations. The…
▽ More
Small-wavevector excitations in Coulomb-interacting systems can be decomposed into the high-energy collective longitudinal plasmon and the low-energy single-electron excitations. At the critical wavevector and corresponding frequency where the plasmon branch merges with the single-electron excitation region, the collective energy of the plasmon dissipates into single electron-hole excitations. The jellium model provides a reasonable description of the electron-energy-loss spectrum (EELS) of metals close to the free-electron limit. The random phase approximation (RPA) is exact in the high-density limit but can capture the plasmonic dispersion reasonably even for densities with rs > 1. RPA and all beyond-RPA methods investigated here, result in a wrong infinite plasmon lifetime for a wavevector smaller than the critical one where the plasmon dispersion curve runs into particle-hole excitations. Exchange-correlation kernel corrections to RPA modify the plasmon dispersion curve. There is however a large difference in the construction and form of the kernels investigated earlier. Our current work introduces recent model exchange-only and exchange-correlation kernels and discusses the relevance of some exact constraints in the construction of the kernel. We show that, because the plasmon dispersion samples a range of wavevectors smaller than the range sampled by the correlation energy, different kernels can make a strong difference for the correlation energy and a weak difference for the plasmon dispersion. This work completes our understanding about the plasmon dispersion in realistic metals, such as Cs, where a negative plasmon dispersion has been observed. We find only positive plasmon dispersion in jellium at the density for Cs.
△ Less
Submitted 15 March, 2020;
originally announced March 2020.
-
Treating different bonding situations: Revisiting Au-Cu alloys using the random phase approximation
Authors:
Niraj K. Nepal,
Santosh Adhikari,
Jefferson E. Bates,
Adrienn Ruzsinszky
Abstract:
The ground state equilibrium properties of copper-gold alloys have been explored with the state of art random phase approximation (RPA). Our estimated lattice constants agree with the experiment within a mean absolute percentage error (MAPE) of 1.4 percent. Semi-local functionals such as the generalized gradient approximation (GGA) of Perdew, Burke, and Ernzerhof (PBE) and strongly constrained and…
▽ More
The ground state equilibrium properties of copper-gold alloys have been explored with the state of art random phase approximation (RPA). Our estimated lattice constants agree with the experiment within a mean absolute percentage error (MAPE) of 1.4 percent. Semi-local functionals such as the generalized gradient approximation (GGA) of Perdew, Burke, and Ernzerhof (PBE) and strongly constrained and appropriately normed (SCAN) fail to provide accurate bulk moduli, which indicate their inability to describe the system in a stretched or compressed state with respect to the equilibrium geometry. We find that the non-locality present in RPA is able to describe the transition between two delocalized electron densities (bulk elemental constituents to crystallized alloys), as required to provide accurate formation energies. Based on our results, we conclude that it is difficult to find a universal density functional which can give accurate results for a wide range of properties of intermetallic alloys. However, RPA can capture different bonding situations, often better than other density functionals. It gives accurate results for a wide range of ground state properties for the alloys, generated from metals with completely filled d-shells.
△ Less
Submitted 23 May, 2019; v1 submitted 22 May, 2019;
originally announced May 2019.
-
First-principles study of mechanical and electronic properties of bent monolayer transition metal dichalcogenides
Authors:
Niraj K. Nepal,
Liping Yu,
Qimin Yan,
Adrienn Ruzsinszky
Abstract:
The mechanical and electronic properties of transition metal dichalcogenide (TMD) monolayers corresponding to transition groups IV, VI, and X are explored under mechanical bending from first principles calculations using the strongly constrained and appropriately normed (SCAN) meta-GGA (MGGA). SCAN provides an accurate description of the phase stability of the TMD monolayers. Our calculated lattic…
▽ More
The mechanical and electronic properties of transition metal dichalcogenide (TMD) monolayers corresponding to transition groups IV, VI, and X are explored under mechanical bending from first principles calculations using the strongly constrained and appropriately normed (SCAN) meta-GGA (MGGA). SCAN provides an accurate description of the phase stability of the TMD monolayers. Our calculated lattice parameters and other structural parameters agree well with experiment. We find that bending stiffness (or flexural rigidity) increases as the transition metal group goes from IV to X to VI, with the exception of PdTe$_2$. Variation in mechanical properties (local strain, physical thickness) and electronic properties (local charge density, band structure) with bending curvature is discussed. The local strain profile of these TMD monolayers under mechanical bending is highly non-uniform. The mechanical bending tunes not only the thickness of the TMD monolayers but also the local charge distribution as well as the band structures, adding more functionalization options to these materials.
△ Less
Submitted 11 June, 2019; v1 submitted 10 April, 2019;
originally announced April 2019.
-
Stretched or noded orbital densities and self-interaction correction in density functional theory
Authors:
Chandra Shahi,
Puskar Bhattarai,
Kamal Wagle,
Biswajit Santra,
Sebastian Schwalbe,
Torsten Hahn,
Jens Kortus,
Koblar A. Jackson,
Juan E. Peralta,
Kai Trepte,
Susi Lehtola,
Niraj K. Nepal,
Hemanadhan Myneni,
Bimal Neupane,
Santosh Adhikari,
Adrienn Ruzsinszky,
Yoh Yamamoto,
Tunna Baruah,
Rajendra R. Zope,
John P. Perdew
Abstract:
Semi-local approximations to the density functional for the exchange-correlation energy of a many-electron system necessarily fail for lobed one-electron densities, including not only the familiar stretched densities but also the less familiar but closely-related noded ones. The Perdew-Zunger (PZ) self-interaction correction (SIC) to a semi-local approximation makes that approximation exact for al…
▽ More
Semi-local approximations to the density functional for the exchange-correlation energy of a many-electron system necessarily fail for lobed one-electron densities, including not only the familiar stretched densities but also the less familiar but closely-related noded ones. The Perdew-Zunger (PZ) self-interaction correction (SIC) to a semi-local approximation makes that approximation exact for all one-electron ground- or excited-state densities and accurate for stretched bonds. When the minimization of the PZ total energy is made over real localized orbitals, the orbital densities can be noded, leading to energy errors in many-electron systems. Minimization over complex localized orbitals yields nodeless orbital densities, which reduce but typically do not eliminate the SIC errors of atomization energies. Other errors of PZ SIC remain, attributable to the loss of the exact constraints and appropriate norms that the semi-local approximations satisfy, and suggesting the need for a generalized SIC. These conclusions are supported by calculations for one-electron densities, and for many-electron molecules. While PZ SIC raises and improves the energy barriers of standard generalized gradient approximations (GGA's) and meta-GGA's, it reduces and often worsens the atomization energies of molecules. Thus PZ SIC raises the energy more as the nodality of the valence localized orbitals increases from atoms to molecules to transition states. PZ SIC is applied here in particular to the SCAN meta-GGA, for which the correlation part is already self-interaction-free. That property makes SCAN a natural first candidate for a generalized SIC.
△ Less
Submitted 2 April, 2019; v1 submitted 1 March, 2019;
originally announced March 2019.