Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA pathway towards high throughput Quantum Monte Carlo simulations for alloys: A case study of two-dimensional (2D) $GaS_xSe_{1-x}$
115
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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 high-entropy alloys Al$_{x}$CrFeCoNi exist over a broad range of Al concentrations ($0 < x < 2$). With increasing Al content their structure is changed from the fcc to bcc phase. We investigate the effect of such structural changes on transport properties including the residual resistivity and the anomalous Hall resistivity. We have performed a detailed comparison of the first-principles simulations with available experimental data. We show that the calculated residual resistivities for all studied alloy compositions are in a fair agreement with available experimental data as concerns both the resistivity values and concentration trends. We emphasize that a good agreement with experiment was obtained also for the anomalous Hall resistivity. We have completed study by estimation of the anisotropic magnetoresistance, spin-disorder resistivity, and Gilbert damping. The obtained results prove that the main scattering mechanism is due to the intrinsic chemical disorder whereas the effect of spin polarization on the residual resistivity is appreciably weaker.
We use the diffusion quantum Monte Carlo (DMC) method to calculate the ground state phase diagram of solid molecular hydrogen and examine the stability of the most important insulating phases relative to metallic crystalline molecular hydrogen. We develop a new method to account for finite-size errors by combining the use of twist-averaged boundary conditions with corrections obtained using the Kwee-Zhang-Krakauer (KZK) functional in density functional theory. To study band-gap closure and find the metallization pressure, we perform accurate quasi-particle many-body calculations using the $GW$ method. In the static approximation, our DMC simulations indicate a transition from the insulating Cmca-12 structure to the metallic Cmca structure at around 375 GPa. The $GW$ band gap of Cmca-12 closes at roughly the same pressure. In the dynamic DMC phase diagram, which includes the effects of zero-point energy, the Cmca-12 structure remains stable up to 430 GPa, well above the pressure at which the $GW$ band gap closes. Our results predict that the semimetallic state observed experimentally at around 360 GPa [Phys. Rev. Lett. {bf 108}, 146402 (2012)] may correspond to the Cmca-12 structure near the pressure at which the band gap closes. The dynamic DMC phase diagram indicates that the hexagonal close packed $P6_3/m$ structure, which has the largest band gap of the insulating structures considered, is stable up to 220 GPa. This is consistent with recent X-ray data taken at pressures up to 183 GPa [Phys. Rev. B {bf 82}, 060101(R) (2010)], which also reported a hexagonal close packed arrangement of hydrogen molecules.
We introduce an algorithm for treating growth on surfaces which combines important features of continuum methods (such as the level-set method) and Kinetic Monte Carlo (KMC) simulations. We treat the motion of adatoms in continuum theory, but attach them to islands one atom at a time. The technique is borrowed from the Dielectric Breakdown Model. Our method allows us to give a realistic account of fluctuations in island shape, which is lacking in deterministic continuum treatments and which is an important physical effect. Our method should be most important for problems close to equilibrium where KMC becomes impractically slow.
While the ongoing search to discover new high-entropy systems is slowly expanding beyond metals, a rational and effective method for predicting in silico the solid solution forming ability of multi-component systems remains yet to be developed. In this article, we propose a novel high-throughput approach, called LTVC, for estimating the transition temperature of a solid solution: ab-initio energies are incorporated into a mean field statistical mechanical model where an order parameter follows the evolution of disorder. The LTVC method is corroborated by Monte Carlo simulations and the results from the current most reliable data for binary, ternary, quaternary and quinary systems (96.6%; 90.7%; 100% and 100%, of correct solid solution predictions, respectively). By scanning through the many thousands of systems available in the AFLOW consortium repository, it is possible to predict a plethora of previously unknown potential quaternary and quinary solid solutions for future experimental validation.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
