-
Electron energy-loss spectrum and exciton band structure of ${\mathrm{WSe}}_{2}$ monolayer studied by ab initio Bethe-Salpeter equation calculations
Authors:
Yun-Chen Shih,
Fredrik Andreas Nilsson,
Guang-Yu Guo
Abstract:
Bounded excitons in transition metal dichalcogenides monolayers lead to numerous opto-electronic applications, which require a detailed understanding of the exciton dynamics. The dynamical properties of excitons with finite momentum transfer $\textbf{Q}$ can be investigated experimentally using electron energy-loss (EEL) spectroscopy. The EEL spectrum depends on the response function of the materi…
▽ More
Bounded excitons in transition metal dichalcogenides monolayers lead to numerous opto-electronic applications, which require a detailed understanding of the exciton dynamics. The dynamical properties of excitons with finite momentum transfer $\textbf{Q}$ can be investigated experimentally using electron energy-loss (EEL) spectroscopy. The EEL spectrum depends on the response function of the material which in turn is determined by the exciton energies and eigenvectors in the exciton band structure. In this work, we utilize ab initio density-functional theory plus Bethe-Salpeter equation (DFT+BSE) approach to explore the exciton band structure and also $\textbf{Q}$-resolved EEL spectrum in monolayer ${\mathrm{WSe}}_{2}$. In particular, we carefully examine the discrepancies and connections among the existing EEL spectrum formulas for quasi-two-dimensional (2D) systems, and establish a proper definition of the EEL spectrum, which is then used to calculate the EEL spectra of monolayer ${\mathrm{WSe}}_{2}$. We find that remarkably, the dispersion of the calculated lowest-energy EELS peaks for the in-plane momentum transfer follows almost precisely the non-parabolic upper band of the lowest bright A exciton, and also agrees well with the previous experiment. Furthermore, we show that only the bright exciton with its electric dipole being parallel to the direction of the transfered momentum is excited, i.e., EEL spectroscopy selectively probes bright exciton bands. This explains why only the upper band of the A exciton, which is a longitudinal exciton with an in-plane dipole moment, was observed in the previous experiment. Our findings will stimulate further EEL experiments to measure other branches of the exciton band structure, such as the parabolic lower band of the A exciton, and hence will lead to a better understanding of the exciton dynamics in quasi-2D materials.
△ Less
Submitted 20 July, 2024;
originally announced July 2024.
-
GPAW: An open Python package for electronic-structure calculations
Authors:
Jens Jørgen Mortensen,
Ask Hjorth Larsen,
Mikael Kuisma,
Aleksei V. Ivanov,
Alireza Taghizadeh,
Andrew Peterson,
Anubhab Haldar,
Asmus Ougaard Dohn,
Christian Schäfer,
Elvar Örn Jónsson,
Eric D. Hermes,
Fredrik Andreas Nilsson,
Georg Kastlunger,
Gianluca Levi,
Hannes Jónsson,
Hannu Häkkinen,
Jakub Fojt,
Jiban Kangsabanik,
Joachim Sødequist,
Jouko Lehtomäki,
Julian Heske,
Jussi Enkovaara,
Kirsten Trøstrup Winther,
Marcin Dulak,
Marko M. Melander
, et al. (22 additional authors not shown)
Abstract:
We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually indepen…
▽ More
We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE) providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe-Salpeter Equation (BSE), variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn-Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support of GPU acceleration has been achieved with minor modifications of the GPAW code thanks to the CuPy library. We end the review with an outlook describing some future plans for GPAW.
△ Less
Submitted 16 April, 2024; v1 submitted 23 October, 2023;
originally announced October 2023.
-
High-throughput search for triplet point defects with narrow emission lines in 2D materials
Authors:
A. Sajid,
F. Nilsson,
S. Manti,
F. Bertoldo,
J. J. Mortensen,
K. S. Thygesen
Abstract:
We employ a first-principles computational workflow to screen for optically accessible, high-spin point defects in wide band gap two-dimensional (2D) crystals. Starting from an initial set of 5388 point defects, comprising both intrinsic and extrinsic, single and double defects in ten previously synthesised 2D host materials, we identify 596 defects with a triplet ground state. For these defects,…
▽ More
We employ a first-principles computational workflow to screen for optically accessible, high-spin point defects in wide band gap two-dimensional (2D) crystals. Starting from an initial set of 5388 point defects, comprising both intrinsic and extrinsic, single and double defects in ten previously synthesised 2D host materials, we identify 596 defects with a triplet ground state. For these defects, we calculate the defect formation energy, the hyperfine (HF) coupling, and the zero-field splitting (ZFS) tensors. For 39 triplet transitions exhibiting particularly low Huang-Rhys factors, we calculate the full photo-luminescence (PL) spectrum. Our approach reveals many new spin defects with narrow PL line shapes and emission frequencies covering a broad spectral range. Most of the defects are hosted in hexagonal BN, which we ascribe to its high stiffness, but some are also found in MgI2, MoS2, MgBr2 and CaI2. As specific examples, we propose the defects vSMoS0 and NiSMoS0 in MoS2 as interesting candidates with potential applications to magnetic field sensors and quantum information technology. All the data will be made available in the open access database QPOD.
△ Less
Submitted 28 May, 2023;
originally announced June 2023.
-
Exciton Superfluidity in 2D Heterostructures from First Principles: The importance of material specific screening
Authors:
Rune Højlund,
Emil Grovn,
Sahar Pakdel,
Kristian S. Thygesen,
Fredrik Nilsson
Abstract:
Recent theoretical and experimental studies suggest that van der Waals heterostructures with n- and p-doped bilayers of transition metal dichalcogenides are promising facilitators of exciton superfluidity. Exciton superfluidity in such bilayer systems is often modelled by solving a mean-field gap equation defined for only the conduction and valence band of the electron and hole material respective…
▽ More
Recent theoretical and experimental studies suggest that van der Waals heterostructures with n- and p-doped bilayers of transition metal dichalcogenides are promising facilitators of exciton superfluidity. Exciton superfluidity in such bilayer systems is often modelled by solving a mean-field gap equation defined for only the conduction and valence band of the electron and hole material respectively. A key quantity entering the gap equation is the effective Coulomb potential acting as the bare interaction in the subspace of the model. Since the model only includes a few bands around the Fermi energy the effective model interaction is partially screened. Although the screening is a material dependent quantity it has in previous works been accounted for in an ad hoc manner, by assuming a static dielectric constant of 2 for a wide range of different materials. In this work we show that the effective model interaction can be derived from first principles using open source code frameworks. Using this novel ab initio downfolding procedure we show that the material dependent screening is essential to predict both magnitude and trends of exciton binding energies and superfluid properties. Furthermore, we suggest new material platforms of both transition metal oxides and dichalcogenides with superior properties compared to the standard devices with two transition metal dichalcogenide layers.
△ Less
Submitted 9 August, 2023; v1 submitted 16 April, 2023;
originally announced April 2023.
-
Excitonic insulators and superfluidity in 2D bilayers without external fields
Authors:
Fredrik Nilsson,
Mikael Kuisma,
Sahar Pakdel,
Kristian S. Thygesen
Abstract:
We explore a new platform for realising excitonic insulators, namely van der Waals (vdW) bilayers comprising 2D Janus materials. In previous studies, Type-II heterobilayers have been brought to the excitonic insulating regime by tuning the band alignment using external gates. In contrast, the Janus bilayers of the current work represent intrinsic excitonic insulators. We first conduct ab initio ca…
▽ More
We explore a new platform for realising excitonic insulators, namely van der Waals (vdW) bilayers comprising 2D Janus materials. In previous studies, Type-II heterobilayers have been brought to the excitonic insulating regime by tuning the band alignment using external gates. In contrast, the Janus bilayers of the current work represent intrinsic excitonic insulators. We first conduct ab initio calculations to obtain the quasiparticle band structures, screened Coulomb interaction, and interlayer exciton binding energies of the bilayers. These ab initio-derived quantities are then used to construct a BCS-like Hamiltonian of the exciton condensate. By solving the mean-field gap equation, we identify 16 vdW Janus bilayers with insulating ground states and superfluid properties. Our calculations expose a new class of advanced materials that are likely to exhibit novel excitonic phases at low temperatures, and highlight the subtle competition between interlayer hybridization, spin-orbit coupling, and dielectric screening that govern their properties.
△ Less
Submitted 9 August, 2023; v1 submitted 19 December, 2022;
originally announced December 2022.
-
Effects of dynamical screening on the BCS-BEC crossover in double bilayer graphene: Density functional theory for exciton bilayers
Authors:
Fredrik Nilsson,
Ferdi Aryasetiawan
Abstract:
We derive a gap equation for bilayer excitonic systems based on density functional theory and benchmark our results against quantum Monte-Carlo simulations and recent experiments on double bilayer graphene. The gap equation has a mean-field form but includes a consistent treatment of dynamical screening. We show that the gap survives at much higher densities than previously thought from mean-field…
▽ More
We derive a gap equation for bilayer excitonic systems based on density functional theory and benchmark our results against quantum Monte-Carlo simulations and recent experiments on double bilayer graphene. The gap equation has a mean-field form but includes a consistent treatment of dynamical screening. We show that the gap survives at much higher densities than previously thought from mean-field estimates which gives strong indications that the double-bilayer graphene systems at zero magnetic field can be used as model systems to investigate the BCS-BEC crossover. Furthermore, we show that Josephson-like transfer of pairs can be substantial for small band gaps and densities.
△ Less
Submitted 20 January, 2021;
originally announced January 2021.
-
Normal state of Nd$_{1-x}$Sr$_x$NiO$_2$ from self-consistent $GW$+EDMFT
Authors:
Francesco Petocchi,
Viktor Christiansson,
Fredrik Nilsson,
Ferdi Aryasetiawan,
Philipp Werner
Abstract:
The recent discovery of superconductivity in hole-doped NdNiO$_2$ thin films has captivated the condensed matter physics community. Such compounds with a formal Ni$^+$ valence have been theoretically proposed as possible analogues of the cuprates, and the exploration of their electronic structure and pairing mechanism may provide important insights into the phenomenon of unconventional superconduc…
▽ More
The recent discovery of superconductivity in hole-doped NdNiO$_2$ thin films has captivated the condensed matter physics community. Such compounds with a formal Ni$^+$ valence have been theoretically proposed as possible analogues of the cuprates, and the exploration of their electronic structure and pairing mechanism may provide important insights into the phenomenon of unconventional superconductivity. At the modeling level, there are however fundamental issues that need to be resolved. While it is generally agreed that the low-energy properties of cuprates can to a large extent be captured by a single-band model, there has been a controversy in the recent literature about the importance of a multi-band description of the nickelates. The origin of this controversy is that studies based entirely on density functional theory (DFT) calculations miss important correlation and multi-orbital effects induced by Hund coupling, while model calculations or simulations based on the combination of DFT and (extended) dynamical mean field theory ((E)DMFT) involve ad-hoc parameters and double counting corrections that substantially affect the results. Here we use a multi-site extension of the recently developed $GW$+EDMFT method, which is free of adjustable parameters, to self-consistently compute the interaction parameters and electronic structure of hole-doped NdNiO$_2$. This full ab-initio simulation demonstrates the importance of a multi-orbital description, even for the undoped compound, and produces results for the resistivity and Hall conductance in qualitative agreement with experiment.
△ Less
Submitted 30 May, 2020;
originally announced June 2020.
-
Screening from $e_g$ states and antiferromagnetic correlations in $d^{(1,2,3)}$ perovskites: A $GW$+EDMFT investigation
Authors:
Francesco Petocchi,
Fredrik Nilsson,
Ferdi Aryasetiawan,
Philipp Werner
Abstract:
We perform a systematic {\it ab initio} study of the electronic structure of Sr(V,Mo,Mn)O$_3$ perovskites, using the parameter-free $GW$+EDMFT method. This approach self-consistently calculates effective interaction parameters, taking into account screening effects due to nonlocal charge fluctuations. Comparing the results of a 3-band ($t_{2g}$) description to those of a 5-band ($t_{2g}$+$e_g$) mo…
▽ More
We perform a systematic {\it ab initio} study of the electronic structure of Sr(V,Mo,Mn)O$_3$ perovskites, using the parameter-free $GW$+EDMFT method. This approach self-consistently calculates effective interaction parameters, taking into account screening effects due to nonlocal charge fluctuations. Comparing the results of a 3-band ($t_{2g}$) description to those of a 5-band ($t_{2g}$+$e_g$) model, it is shown that the $e_g$ states have little effect on the low-energy properties and the plasmonic features for the first two compounds but play a more active role in SrMnO$_3$. In the case of SrMnO$_3$ paramagnetic $GW$+EDMFT yields a metallic low-temperature solution on the verge of a Mott transition, while antiferromagnetic $GW$+EDMFT produces an insulating solution with the correct gap size. We discuss the possible implications of this result for the nature of the insulating state above the Néel temperature, and the reliability of the $GW$+EDMFT scheme.
△ Less
Submitted 22 November, 2019;
originally announced November 2019.
-
Position Representation of Effective Electron-Electron Interactions in Solids
Authors:
Tor Jonas Sjöstrand,
Fredrik Nilsson,
Christoph Friedrich,
Ferdi Aryasetiawan
Abstract:
An essential ingredient in many model Hamiltonians, such as the Hubbard model, is the effective electron-electron interaction $U$, which enters as matrix elements in some localized basis. These matrix elements provide the necessary information in the model, but the localized basis is incomplete for describing $U$. We present a systematic scheme for computing the manifestly basis-independent dynami…
▽ More
An essential ingredient in many model Hamiltonians, such as the Hubbard model, is the effective electron-electron interaction $U$, which enters as matrix elements in some localized basis. These matrix elements provide the necessary information in the model, but the localized basis is incomplete for describing $U$. We present a systematic scheme for computing the manifestly basis-independent dynamical interaction in position representation, $U({\bf r},{\bf r}';ω)$, and its Fourier transform to time domain, $U({\bf r},{\bf r}';τ)$. These functions can serve as an unbiased tool for the construction of model Hamiltonians. For illustration we apply the scheme within the constrained random-phase approximation to the cuprate parent compounds La$_2$CuO$_4$ and HgBa$_2$CuO$_4$ within the commonly used 1- and 3-band models, and to non-superconducting SrVO$_{3}$ within the $t_{2g}$ model. Our method is used to investigate the shape and strength of screening channels in the compounds. We show that the O 2$p_{x,y}-$Cu 3$d_{x^2-y^2}$ screening gives rise to regions with strong attractive static interaction in the minimal (1-band) model in both cuprates. On the other hand, in the minimal ($t_{2g}$) model of SrVO$_3$ only regions with a minute attractive interaction are found. The temporal interaction exhibits generic damped oscillations in all compounds, and its time-integral is shown to be the potential caused by inserting a frozen point charge at $τ=0$. When studying the latter within the three-band model for the cuprates, short time intervals are found to produce a negative potential.
△ Less
Submitted 4 February, 2019;
originally announced February 2019.
-
Dynamically screened Coulomb interaction in the parent compounds of hole-doped cuprates, trends and exceptions
Authors:
Fredrik Nilsson,
Krister Karlsson,
Ferdi Aryasetiawan
Abstract:
Although the cuprate high-temperature superconductors were discovered already 1986 the origin of the pairing mechanism remains elusive. While the doped compounds are superconducting with high transition temperatures $T_{c}$ the undoped compounds are insulating due to the strong effective Coulomb interaction between the Cu $3d$ holes. We investigate the dependence of the maximum superconducting tra…
▽ More
Although the cuprate high-temperature superconductors were discovered already 1986 the origin of the pairing mechanism remains elusive. While the doped compounds are superconducting with high transition temperatures $T_{c}$ the undoped compounds are insulating due to the strong effective Coulomb interaction between the Cu $3d$ holes. We investigate the dependence of the maximum superconducting transition temperature, $T_{c\text{ max}}$, on the onsite effective Coulomb interaction $U$ using the constrained random-phase approximation. We focus on the commonly used one-band model of the cuprates, including only the antibonding combination of the Cu $d_{x^2-y^2}$ and O $p_x$ and $p_y$ orbitals, and find a clear screening dependent trend between the static value of $U$ and $T_{c\text{ max}}$ for the parent compounds of a large number of hole-doped cuprates. Our results suggest that superconductivity is favored by a large onsite Coulomb repulsion. We analyze both the trend in the static value of $U$ and its frequency dependence in detail and, by comparing to other works, speculate on the mechanisms behind the trend.
△ Less
Submitted 19 December, 2018;
originally announced December 2018.
-
Multitier self-consistent $GW$+EDMFT
Authors:
Fredrik Nilsson,
Lewin Boehnke,
Philipp Werner,
Ferdi Aryasetiawan
Abstract:
We discuss a parameter-free and computationally efficient ab initio simulation approach for moderately and strongly correlated materials, the multitier self-consistent $GW$+EDMFT method. This scheme treats different degrees of freedom, such as high-energy and low-energy bands, or local and nonlocal interactions, within appropriate levels of approximation, and provides a fully self-consistent descr…
▽ More
We discuss a parameter-free and computationally efficient ab initio simulation approach for moderately and strongly correlated materials, the multitier self-consistent $GW$+EDMFT method. This scheme treats different degrees of freedom, such as high-energy and low-energy bands, or local and nonlocal interactions, within appropriate levels of approximation, and provides a fully self-consistent description of correlation and screening effects in the solid. The ab initio input is provided by a one-shot $G^0W^0$ calculation, while the strong-correlation effects originating from narrow bands near the Fermi level are captured by a combined $GW$ plus extended dynamical mean-field (EDMFT) treatment. We present the formalism and technical details of our implementation and discuss some general properties of the effective EDMFT impurity action. In particular, we show that the retarded impurity interactions can have non-causal features, while the physical observables, such as the screened interactions of the lattice system, remain causal. We then turn to stretched sodium as a model system to explore the performance of the multitier self-consistent $GW$+EDMFT method in situations with different degrees of correlation. While the results for the physical lattice spacing $a_0$ show that the scheme is not very accurate for electron-gas like systems, because nonlocal corrections beyond $GW$ are important, it does provide physically correct results in the intermediate correlation regime, and a Mott transition around a lattice spacing of $1.5a_0$. Remarkably, even though the Wannier functions in the stretched compound are less localized, and hence the bare interaction parameters are reduced, the self-consistently computed impurity interactions show the physically expected trend of an increasing interaction strength with increasing lattice spacing.
△ Less
Submitted 29 September, 2017; v1 submitted 21 June, 2017;
originally announced June 2017.
-
When strong correlations become weak: Consistent merging of $GW$ and DMFT
Authors:
Lewin Boehnke,
Fredrik Nilsson,
Ferdi Aryasetiawan,
Philipp Werner
Abstract:
The cubic perovskite SrVO$_3$ is generally considered to be a prototype strongly correlated metal with a characteristic three-peak structure of the $d$-electron spectral function, featuring a renormalized quasiparticle band in between pronounced Hubbard sidebands. Here we show that this interpretation, which has been supported by numerous "ab-initio" simulations, has to be reconsidered. Using a fu…
▽ More
The cubic perovskite SrVO$_3$ is generally considered to be a prototype strongly correlated metal with a characteristic three-peak structure of the $d$-electron spectral function, featuring a renormalized quasiparticle band in between pronounced Hubbard sidebands. Here we show that this interpretation, which has been supported by numerous "ab-initio" simulations, has to be reconsidered. Using a fully self-consistent $GW$+extended dynamical mean-field theory calculation we find that the screening from nonlocal Coulomb interactions substantially reduces the effective local Coulomb repulsion, and at the same time leads to strong plasmonic effects. The resulting effective local interactions are too weak to produce pronounced Hubbard bands in the local spectral function, while prominent plasmon satellites appear at energies which agree with those of the experimentally observed sidebands. Our results demonstrate the important role of nonlocal interactions and dynamical screening in determining the effective interaction strength of correlated compounds.
△ Less
Submitted 14 November, 2016; v1 submitted 7 April, 2016;
originally announced April 2016.
-
Dynamical screening in La2CuO4
Authors:
Philipp Werner,
Rei Sakuma,
Fredrik Nilsson,
Ferdi Aryasetiawan
Abstract:
We show that the dynamical screening of the Coulomb interaction among Cu-d electrons in high-Tc cuprates is very strong and that a proper treatment of this effect is essential for a consistent description of the electronic structure. In particular, we find that ab-initio calculations for undoped La2CuO4 yield an insulator only if the frequency dependence of the Coulomb interaction is taken into ac…
▽ More
We show that the dynamical screening of the Coulomb interaction among Cu-d electrons in high-Tc cuprates is very strong and that a proper treatment of this effect is essential for a consistent description of the electronic structure. In particular, we find that ab-initio calculations for undoped La2CuO4 yield an insulator only if the frequency dependence of the Coulomb interaction is taken into account. We also identify a collective excitation in the screened interaction at 9 eV which is rather localized on the copper site, and which is responsible for a satellite structure at energy -13 eV, located below the p bands.
△ Less
Submitted 14 November, 2014;
originally announced November 2014.