Two-dimensional transition metal dichalcogenides (2D-TMDs) have been proposed as novel optoelectronic materials for space applications due to their relatively light weight. MoS2 has been shown to have excellent semiconducting and photonic properties. Although the strong interaction of ionizing gamma radiation with bulk materials has been demonstrated, understanding its effect on atomically thin materials has scarcely been investigated. Here, we report the effect of gamma irradiation on the structural and electronic properties of a monolayer of MoS2. We perform Raman spectroscopy and X-ray photoelectron spectroscopy (XPS) studies of MoS2, before and after gamma ray irradiation with varying doses and density functional theory (DFT) calculations. The Raman spectra and XPS results demonstrate that point defects dominate after the gamma irradiation of MoS2. DFT calculations elucidate the electronic properties of MoS2 before and after irradiation. Our work makes several contributions to the field of 2D materials research. First, our study of the electronic density of states and the electronic properties of a MoS2 monolayer irradiated by gamma rays sheds light on the properties of a MoS2 monolayer under gamma irradiation. Second, our study confirms that point defects are formed as a result of gamma irradiation. And third, our DFT calculations qualitatively suggest that the conductivity of the MoS2 monolayer may increase after gamma irradiation due to the creation of additional defect states.
more »
« less
This content will become publicly available on September 1, 2026
Toward improved property prediction of 2D materials using many-body quantum Monte Carlo methods
The field of 2D materials has grown dramatically in the past two decades. 2D materials can be utilized for a variety of next-generation optoelectronic, spintronic, clean energy, and quantum computing applications. These 2D structures, which are often exfoliated from layered van der Waals materials, possess highly inhomogeneous electron densities and can possess short- and long-range electron correlations. The complexities of 2D materials make them challenging to study with standard mean-field electronic structure methods such as density functional theory (DFT), which relies on approximations for the unknown exchange-correlation functional. To overcome the limitations of DFT, highly accurate many-body electronic structure approaches such as diffusion Monte Carlo (DMC) can be utilized. In the past decade, DMC has been used to calculate accurate magnetic, electronic, excitonic, and topological properties in addition to accurately capturing interlayer interactions and cohesion and adsorption energetics of 2D materials. This approach has been applied to 2D systems of wide interest, including graphene, phosphorene, MoS2, CrI3, VSe2, GaSe, GeSe, borophene, and several others. In this review article, we highlight some successful recent applications of DMC to 2D systems for improved property predictions beyond standard DFT.
more »
« less
- Award ID(s):
- 2213398
- PAR ID:
- 10642113
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- Applied Physics Reviews
- Volume:
- 12
- Issue:
- 3
- ISSN:
- 1931-9401
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The theorems of density functional theory (DFT) establish bijective maps between the local external potential of a many-body system and its electron density, wavefunction and, therefore, one-particle reduced density matrix. Building on this foundation, we show that machine learning models based on the one-electron reduced density matrix can be used to generate surrogate electronic structure methods. We generate surrogates of local and hybrid DFT, Hartree-Fock and full configuration interaction theories for systems ranging from small molecules such as water to more complex compounds like benzene and propanol. The surrogate models use the one-electron reduced density matrix as the central quantity to be learned. From the predicted density matrices, we show that either standard quantum chemistry or a second machine-learning model can be used to compute molecular observables, energies, and atomic forces. The surrogate models can generate essentially anything that a standard electronic structure method can, ranging from band gaps and Kohn-Sham orbitals to energy-conserving ab-initio molecular dynamics simulations and infrared spectra, which account for anharmonicity and thermal effects, without the need to employ computationally expensive algorithms such as self-consistent field theory. The algorithms are packaged in an efficient and easy to use Python code, QMLearn, accessible on popular platforms.more » « less
-
GW approximation is one of the most popular parameter-free many-body methods that go beyond the limitations of the standard density functional theory (DFT) to determine the excitation spectra for moderately correlated materials and in particular the semiconductors. It is also the first step in developing the diagrammatic Monte Carlo method into an electronic structure tool, which would offer a numerically exact solution to the solid-state problem. While most electronic structure packages offer support for GW calculations for band-insulating materials, the level of support for metallic systems is somewhat limited. This limitation can be partly attributed to the relatively minor differences often observed between GW and DFT results in treating metallic systems, which is not expected to persist to higher orders in perturbation theory. Describing metals within the GW framework presents a challenge, as it requires accurate resolution of Fermi surface singularities, which, in turn, calls for a dense momentum mesh. Here we implement the GW algorithm within the all-electron Linear Augmented Plane Wave framework, where we pay special attention to the metallic systems, the convergence with respect to momentum mesh, and proper treatment of the deep laying core states, as needed for the future variational diagrammatic Monte Carlo implementation. Our improved algorithm for resolving Fermi surface singularities allows us a stable and accurate analytic continuation of imaginary axis data, which is carried out for GW excitation spectra throughout the Brillouin zone in both the metallic and insulating materials and is compared to numerically more stable contour deformation integration technique. We compute band structures for elemental metallic systems Li, Na, and Mg as well as for various narrow and wide bandgap insulators such as Si, BN, SiC, MgO, LiF, ZnS, and CdS and compare our results with previous GW calculations and available experiments data. Our results are in good agreement with the available literature. Thus our software allows users to compute full bandstructures for metals and insulators using all-electron potential without downfolding to Wannier orbital basis.more » « less
-
Discovering new materials with desired properties has been a dominant and crucial topic of interest in the field of materials science in the past few decades. In this work, novel carbon allotropes and ternary B–C–N structures were generated using the state-of-the-art RG 2 code. All structures were fully optimized using density functional theory with first-principles calculations. Several hundred carbon allotropes and ternary B–C–N structures were identified to be superhard materials. The thermodynamic stability of some randomly selected superhard materials was confirmed by evaluating the full phonon dispersions in the Brillouin zone. The new carbon allotropes and ternary B–C–N structures possess a wide range of mechanical properties generally and Vickers hardness specifically. Through 2D Pearson's correlation map, we first reproduced the well-accepted explanation and relationship of the Vickers hardness of the generated structures with other mechanical properties such as shear modulus, bulk modulus, Pugh's ratio, universal anisotropy, and Poisson's ratio. We then propose two fundamentally new descriptors from the electronic level, namely local potential and electron localization function averaged over a unit cell, both of which exhibit a strong correlation with Vickers hardness. More importantly, these descriptors are easy to access from first-principles calculations (at least two orders of magnitude faster than the traditional calculation of elastic constants), and thus can serve as a fast and accurate approach for screening superhard materials. We also combined these new descriptors with known composition and structural descriptors in the machine learning training process. The new descriptors significantly enhance the performance of the trained machine learning model in predicting the Vickers hardness of unknown materials, which provides strong evidence for local potential and electron localization function to be considered in future high-throughput computation. This work unravels more fundamental but previously unexplored knowledge about superhard materials and the newly proposed electronic level descriptors are expected to accelerate the discovery of new superhard materials.more » « less
-
The genomics approach to materials, heralded by increasingly accurate density functional theory (DFT) calculations conducted on thousands of crystalline compounds, has led to accelerated material discovery and property predictions. However, so far, amorphous materials have been largely excluded from this as these systems are notoriously difficult to simulate. Here, we study amorphous Ta2O5 thin films mixed with Al2O3, SiO2, Sc2O3, TiO2, ZnO, ZrO2, Nb2O5, and HfO2 to identify their crystalline structure upon post-deposition annealing in air both experimentally and with simulations. Using the Materials Project open database, phase diagrams based on DFT calculations are constructed for the mixed oxide systems and the annealing process is evaluated via grand potential diagrams with varying oxygen chemical potential. Despite employing calculations based on crystalline bulk materials, the predictions agree well with the experimentally observed crystallized phases of the amorphous thin films. In the absence of ternary phases, the dopant acts as an amorphizer agent increasing the thermal stability of Ta2O5. The least efficient amorphizer agent is found to be Nb2O5, for which the cation has similar chemical properties to those of Ta in Ta2O5. These results show that DFT calculations can be applied for the prediction of crystallized structures of annealed amorphous materials. This could pave the way for accelerated in silico material discovery and property predictions using the powerful genomic approach for amorphous oxide coatings employed in a wide range of applications such as optical coatings, energy storage, and electronic devices.more » « less
An official website of the United States government
