Search

High-throughput density functional theory (DFT) calculations allow for a systematic search for conventional superconductors. With the recent interest in two-dimensional (2D) superconductors, we used a high-throughput workflow to screen over 1,000 2D materials in the JARVIS-DFT database and performed electron-phonon coupling calculations, using the McMillan-Allen-Dynes formula to calculate the superconducting transition temperature ($T_c$) for 165 of them. Of these 165 materials, we identify 34 dynamically stable structures with transition temperatures above 5 K, including materials such as W$_2$N$_3$, NbO$_2$, ZrBrO, TiClO, NaSn$_2$S$_4$, Mg$_2$B$_4$C$_2$ and the previously unreported Mg$_2$B$_4$N$_2$ ($T_c$ = 21.8 K). Finally, we performed experiments to determine the $T_c$ of selected layered superconductors (2H-NbSe$_2$, 2H-NbS$_2$, ZrSiS, FeSe) and discuss the measured results within the context of our DFT results. We aim that the outcome of this workflow can guide future computational and experimental studies of new and emerging 2D superconductors by providing a roadmap of high-throughput DFT data.

The study of alloys using computational methods has been a difficult task due to the usually unknown stoichiometry and local atomic ordering of the different structures experimentally. In order to combat this, first-principles methods have been coupled with statistical methods such as the Cluster Expansion formalism in order to construct the energy hull diagram, which helps to determine if an alloyed structure can exist in nature. Traditionally, density functional theory (DFT) has been used in such workflows. In this work we propose to use chemically accurate many-body variational Monte Carlo (VMC) and diffusion Monte Carlo (DMC) methods to construct the energy hull diagram of an alloy system, due to the fact that such methods have a weaker dependence on the starting wavefunction and density functional, scale similarly to DFT with the number of electrons, and have had demonstrated success for a variety of materials. To carry out these simulations in a high-throughput manner, we propose a method called Jastrow sharing, which involves recycling the optimized Jastrow parameters between alloys with different stoichiometries. We show that this eliminates the need for extra VMC Jastrow optimization calculations and results in a significant computational cost savings (on average 1/4 savings of total computational time). Since it is a novel post-transition metal chalcogenide alloy series that has been synthesized in its few-layer form, we used monolayer $GaS_xSe_{1-x}$ as a case study for our workflow. By extensively testing our Jastrow sharing procedure for monolayer $GaS_xSe_{1-x}$ and quantifying the cost savings, we demonstrate how a pathway towards chemically accurate high-throughput simulations of alloys can be achieved using many-body VMC and DMC methods.

Two-dimensional (2D) post-transition metal chalcogenides (PTMC) have attracted attention due to their suitable band gaps and lower exciton binding energies, making them more appropriate for electronic, optical and water-splitting devices than graphene and monolayer transition metal dichalcogenides (TMDs). Of the predicted 2D PTMCs, GaSe has been reliably synthesized and experimentally characterized. Despite this fact, quantities such as lattice parameters and band character vary significantly depending on which density functional theory (DFT) functional is used. Although many-body perturbation theory (GW approximation) has been used to correct the electronic structure and obtain the excited state properties of 2D GaSe, and solving the Bethe-Salpeter equation (BSE) has been used to find the optical gap, we find that the results depend strongly on the starting wavefunction. In attempt to correct these discrepancies, we employed the many-body Diffusion Monte Carlo (DMC) method to calculate the ground and excited state properties of GaSe because DMC has a weaker dependence on the trial wavefunction. We benchmark these results with available experimental data, DFT [local-density approximation, Perdew-Burke-Ernzerhof (PBE), strongly constrained and appropriately normed (SCAN) meta-GGA, and hybrid (HSE06) functionals] and GW-BSE (using PBE and SCAN wavefunctions) results. Our findings confirm monolayer GaSe is an indirect gap semiconductor (Gamma-M) with a quasiparticle electronic gap in close agreement with experiment and low exciton binding energy. We also benchmark the optimal lattice parameter, cohesive energy and ground state charge density with DMC and various DFT methods. We aim to present a terminal theoretical benchmark for pristine monolayer GaSe, which will aide in the further study of 2D PTMCs using DMC methods.