No Arabic abstract
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.
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.
Recent experiments revealed that monolayer $alpha$-RuCl$_3$ can be obtain by chemical exfoliation method and exfoliation or restacking of nanosheets can manipulate the magnetic properties of the materials. In this present paper, the electronic and magnetic properties of $alpha$-RuCl$_3$ monolayer are investigated by combining first-principles calculations and Monte Carlo simulations. From first-principles calculations, we found that the spin configuration FM corresponds to the ground state for $alpha$-RuCl$_3$, however, the other excited zigzag oriented spin configuration has energy of 5 meV/atom higher than the ground state. Energy band gap has been obtained as $3$ meV using PBE functionals. When spin-orbit coupling effect is taken into account, corresponding energy gap is determined to be as $57$ meV. We also investigate the effect of Hubbard U energy terms on the electronic band structure of $alpha$-RuCl$_3$ monolayer and revealed band gap increases approximately linear with increasing U value. Moreover, spin-spin coupling terms ($J_1$, $J_2$, $J_3$) have been obtained using first principles calculations. By benefiting from these terms, Monte Carlo simulations with single site update Metropolis algorithm have been implemented to elucidate magnetic properties of the considered system. Thermal variations of magnetization, susceptibility and also specific heat curves indicate that monolayer $alpha$-RuCl$_3$ exhibits a phase transition between ordered and disordered phases at the Curie temperature $14.21$ K. We believe that this study can be utilized to improve two-dimensional magnet materials.
We study a three-orbital Su-Schrieffer-Heeger model defined on a two-dimensional Lieb lattice and in the negative charge transfer regime using determinant quantum Monte Carlo. At half-filling (1 hole/unit cell), we observe a bipolaron insulating phase, where the ligand oxygen atoms collapse and expand about alternating cation atoms to produce a bond-disproportionated state. This phase is robust against moderate hole doping but is eventually suppressed at large hole concentrations, leading to a metallic polaron-liquid-like state with fluctuating patches of local distortions. Our results suggest that the polarons are highly disordered in the metallic state and freeze into a periodic array across the metal-to-insulator transition. We also find an $s$-wave superconducting state at finite doping that primarily appears on the oxygen sublattices. Our approach provides an efficient, non-perturbative way to treat bond phonons in higher dimensions and our results have implications for many materials where coupling to bond phonons is the dominant interaction.
Charged defects in 2D materials have emerging applications in quantum technologies such as quantum emitters and quantum computation. Advancement of these technologies requires rational design of ideal defect centers, demanding reliable computation methods for quantitatively accurate prediction of defect properties. We present an accurate, parameter-free and efficient procedure to evaluate quasiparticle defect states and thermodynamic charge transition levels of defects in 2D materials. Importantly, we solve critical issues that stem from the strongly anisotropic screening in 2D materials, that have so far precluded accurate prediction of charge transition levels in these materials. Using this procedure, we investigate various defects in monolayer hexagonal boron nitride (h-BN) for their charge transition levels, stable spin states and optical excitations. We identify $C_BN_V$ (nitrogen vacancy adjacent to carbon substitution of boron) to be the most promising defect candidate for scalable quantum bit and emitter applications.
We report large scale determinant Quantum Monte Carlo calculations of the effective bandwidth, momentum distribution, and magnetic correlations of the square lattice fermion Hubbard Hamiltonian at half-filling. The sharp Fermi surface of the non-interacting limit is significantly broadened by the electronic correlations, but retains signatures of the approach to the edges of the first Brillouin zone as the density increases. Finite size scaling of simulations on large lattices allows us to extract the interaction dependence of the antiferromagnetic order parameter, exhibiting its evolution from weak coupling to the strong coupling Heisenberg limit. Our lattices provide improved resolution of the Greens function in momentum space, allowing a more quantitative comparison with time-of-flight optical lattice experiments.