Many-body Perturbation Theory Research Articles

Exploring the structural properties and the optoelectronic features of RbPbX3 (X [formula omitted] Cl, F) perovskite crystals for solar cells solicitations: Showcasing the DFT predictions.

A comprehensive study of the structural, optical, and electronic properties of RbPbX3 (where X is Cl or F) perovskite crystals is reported herein. The calculations of the different physical properties of these two distinguishable RbPbCl3 and RbPbF3 perovskite crystals have been carried within the framework of the Density Functional Theory (DFT) approach for the cubic temperature-dependent phase. The study is motivated by the exploration of innovative perovskites materials for solar cells application due to the marvelous augmentation in power conversion efficiency (PCE) obtained by perovskite crystals (≈ 25% up to date). The investigation were undertaken for determining the perceivable role of the X element in their band gap energy features as well as the band structure proprieties. The simulations were commenced by determining the optimal lattice constants and the lowest total energy of both perovskite crystals. The exploration of the electronic properties of RbPbCl3 and RbPbF3 perovskite crystals have been studied by estimating their variation upon employing spin–orbit-coupling (SOC) and without taking into consideration the SOC interaction by means of LDA-pz and GGA-PBE exchange potentials within the framework of projector augmented wave method (PAW) and the Perdew–Burke–Ernzerh of LDA and (PBE) variant of the generalized gradient approximation (GGA). The integrated suite quantum-ESPRESSO is used for materials modeling, structure calculations, and electronic properties of the studied RbPbCl3 and RbPbF3 perovskite crystals. The thermo_pw driver is employed for the determination of the elastic properties whereas the optical properties were elucidated by aid of Yambo code through solving the many-body perturbation theory (MBPT) approach as introduced in the open-source code. The optical band gap of cubic RbPbCl3 and RbPbCF3 within LDA and PBE, were observed at (2.18, 2.2) eV and (2.4, 2.5) eV respectively. Compared and contrasted with those revealed by the electronic band gap energy of the systems, these obtained values are very consistent with the experimental. The two maximum peaks for the absorption coefficients are found to be within the visible range rending the two perovskite crystals in promising candidate for solar cell applications. Despite some quantitative differences, the material’s enhancement potential is well demonstrated by the absorption coefficient spectrum, refractive index, and conductivity. The study holds promise for solar cell applications and may pave the way for consideration of the RbPbCl3 and RbPbF3 as well as their doped analogs for perovskites-based solar cells as the next generation in the photovoltaic industry.

Correlation balance for describing carbenes: An NOF study

Carbenes are a remarkable type of molecules because of the two electrons that might be paired or unpaired, providing a singlet or a triplet state, respectively. The preference for one of these states can be measured by the singlet-triplet energy gap; however, theoretical prediction is challenging when static correlation arises, and an adequate balance with dynamic correlation is required to achieve correct predictions. Piris Natural Orbital Functionals (PNOFs) have been used before to deal with other static correlation problems, but they have suffered from lack of dynamic correlation, which has been solved by coupling PNOF7 with many-body perturbation theory. Recently, the development of a Global Natural Orbital Functional (GNOF) has been proposed with the aim of including dynamic correlation without the need to couple with perturbation theory, thus becoming a promising alternative to study challenging chemical problems. In this work, we applied the PNOF family of functionals to compute the adiabatic singlet-triplet energy gaps of a set of simple carbenes and compared them with those computed with coupled-cluster methods and experimental values when available. We have found that the GNOF achieves promising results due to an intrinsic balance of static and dynamic correlation. In this regard, the GNOF presents errors that are lower than those of PNOF7, comparable to those of NOF-c-MP2, and is capable of predicting the tendency of substituent effects, proving to be suitable for further application to predict general singlet-triplet energy gaps.

Charge transfer and electronic relaxation effects in the photoemission of EMIM-DCA ionic liquid vapor

Ultraviolet photoelectron spectroscopy (UPS) measurements of EMIM-DCA vapors have been conducted. Density functional theory (DFT), Møller–Plesset perturbation theory (MP2) and many-body perturbation theory calculations (GW) have been used to interpret the experimentally measured spectra. The vapor of ionic liquids normally consists of isolated ion-pairs and the corresponding photoemission spectra should reflect the electronic binding energies of the molecular orbitals. The UPS spectra have revealed that many photolines are shifted when compared to ab initio calculation results. The hybrid DFT functionals, which performed very good in the description of the UPS spectra of many other ionic liquid vapors have shortcomings in the EMIM-DCA case. The MP2 calculation and the quasiparticle GW approach are shown to offer a better description of the experimental UPS spectra. The charge distribution analysis showed that there seems to be a favorable charge transfer from the cation to the anion upon ionization. The new data provides strong support to the argument that in some cases DFT still has difficulties in the modeling of IL ion-pairs.

Spin–orbit effects on the electronic and optical properties of lead iodide

Lead iodide (PbI2) has gained much interest due to its direct electronic gap in the visible range and layered crystal structure. It has thereby been considered as a promising material for applications in atomically thin optoelectronic devices. In this work, we present a detailed investigation of the effect of spin–orbit coupling (SOC) that arises from the presence of heavy atoms on the electronic and optical properties of PbI2 using first-principles calculations based on density-functional theory and many-body perturbation theory. We find that SOC not only alters the bandgap but also induces the mixing of orbital characters, resulting in a significant change in the overall band structure and charge carrier effective masses. Moreover, the band orbital mixing caused by SOC results in the dramatic change in optical transition matrix elements and, correspondingly, the absorption spectrum. Our experimentally measured absorption spectra validate the calculation results and demonstrate the importance of SOC in the optical processes of PbI2. Our findings provide insights that are important for the potential use of PbI2 as a material platform for visible optoelectronic devices.

Towards high-throughput many-body perturbation theory: efficient algorithms and automated workflows

The automation of ab initio simulations is essential in view of performing high-throughput (HT) computational screenings oriented to the discovery of novel materials with desired physical properties. In this work, we propose algorithms and implementations that are relevant to extend this approach beyond density functional theory (DFT), in order to automate many-body perturbation theory (MBPT) calculations. Notably, an algorithm pursuing the goal of an efficient and robust convergence procedure for GW and BSE simulations is provided, together with its implementation in a fully automated framework. This is accompanied by an automatic GW band interpolation scheme based on maximally localized Wannier functions, aiming at a reduction of the computational burden of quasiparticle band structures while preserving high accuracy. The proposed developments are validated on a set of representative semiconductor and metallic systems.

Wigner Gaussian dynamics: Simulating the anharmonic and quantum ionic motion

The atomic motion controls important features of materials, such as thermal transport, phase transitions, and vibrational spectra. However, the simulation of ionic dynamics is exceptionally challenging when quantum fluctuations are relevant (e.g., at low temperatures or with light atoms) and the energy landscape is anharmonic. In this work, we formulate the Time-Dependent Self-Consistent Harmonic Approximation (TDSCHA) in the Wigner framework, paving the way for the efficient computation of the nuclear motion in systems with sizable quantum and thermal anharmonic fluctuations. Besides the improved numerical efficiency, the Wigner formalism unveils the classical limit of TDSCHA and provides a link with the many-body perturbation theory of Feynman diagrams. We further extend the method to account for the non-linear couplings between phonons and photons, responsible, e.g., for a nonvanishing Raman signal in high-symmetry Raman inactive crystals, firstly discussed by Rasetti and Fermi. We benchmark the method in phase III of high-pressure hydrogen ab initio. The nonlinear photon-phonon coupling reshapes the IR spectra and explains the high-frequency shoulder of the H2 vibron observed in experiments. The Wigner TDSCHA is computationally cheap and derived from first principles: it is unbiased by assumptions on the phonon-phonon and phonon-photon scattering and does not depend on empirical parameters. Therefore, the method can be adopted in unsupervised high-throughput calculations.

Dielectric screening at TMD:hBN interfaces: Monolayer-to-bulk transition, local-field effect, and spatial dependence

The dielectric effects of a substrate have been shown to be important for modulating the electronic properties of an adsorbate, especially in van der Waals heterostructures. Here, using the first-principles dielectric embedding $GW$ approach within the framework of many-body perturbation theory, we perform a case study on the dielectric screening effects of hexagonal boron nitride (hBN) on various transition-metal dichalcogenides (TMDs). We consider three systems, monolayer ${\mathrm{MoS}}_{2}$, bilayer ${\mathrm{MoS}}_{2}$, and a mixed ${\mathrm{WS}}_{2}/{\mathrm{MoS}}_{2}$ bilayer adsorbed on hBN, and examine three aspects of the substrate dielectric screening: (i) the thickness dependence and the monolayer-to-bulk transition, for which we consider the effects of one-, two-, three-, and four-layer hBN; (ii) the local-field effect, for which we numerically assess a common approximation of neglecting the in-plane local-field components of the substrate polarizability; and (iii) the spatial dependence, for which we consider the mixed ${\mathrm{WS}}_{2}/{\mathrm{MoS}}_{2}$ bilayer adsorbed on hBN with either side facing the substrate. Our results provide quantitative insight into how the substrate screening effects can be leveraged for band structure engineering.

Modewise Johnson–Lindenstrauss embeddings for nuclear many-body theory

In this work, we explore modewise Johnson-Lindenstrauss embeddings (JLEs) as a tool to reduce the computational cost and memory requirements of nuclear many-body methods. JLEs are randomized projections of high-dimensional data tensors onto low-dimensional subspaces that preserve key structural features. Such embeddings allow for the oblivious and incremental compression of large tensors, e.g., the nuclear Hamiltonian, into significantly smaller random sketches that still allow for the accurate calculation of ground-state energies and other observables. Their oblivious character makes it possible to compress a tensor without knowing in advance exactly what observables one might want to approximate at a later time. This opens the door for the use of tensors that are much too large to store in memory, e.g., complete two-plus three-nucleon Hamiltonians in large, symmetry-unrestricted bases. Such compressed Hamiltonians can be stored and used later on with relative ease. As a first step, we analyze the JLE's impact on the second-order Many-Body Perturbation Theory (MBPT) corrections for nuclear ground-state observables. Numerical experiments for a wide range of closed-shell nuclei, model spaces and state-of-the-art nuclear interactions demonstrate the validity and potential of the proposed approach: We can compress nuclear Hamiltonians hundred- to thousand-fold while only incurring mean relative errors of 1\% or less in ground-state observables. Importantly, we show that JLEs capture the relevant physical information contained in the highly structured Hamiltonian tensor despite their random characteristics. In addition to the significant storage savings, the achieved compressions imply multiple order-of magnitude reductions in computational effort when the compressed Hamiltonians are used in higher-order MBPT or nonperturbative many-body methods.

Optical Tunable Moiré Excitons in Twisted Hexagonal GaTe Bilayers

Optical fine-tunable layer-hybridized Moiré excitons are highly in demand for emerging many-body states in two-dimensional semiconductors. We report naturally confined layer-hybridized bright Moiré excitons with long lifetimes in twisted hexagonal GaTe bilayers, using ab initio many-body perturbation theory and the Bethe–Salpeter equation. Due to the hybridization of electrons and holes between layers, which enhances the brightness of excitons, the twisted bilayer system becomes attractive for optical applications. We find that in both R and H-type stacking Moiré superlattices, more than 200 meV lateral quantum confinements occur on exciton energies, which results in two scenarios: (1) The ground state bright excitons XA are found to be trapped at two high-symmetry points, with opposite electric dipoles in the R-stacking Moiré supercell, forming a honeycomb superlattice of nearest-neighbor dipolar attraction. (2) For H-stacking case, the XA is found to be trapped at only one high-symmetry point exhibiting a triangular superlattice. Our results suggest that twisted h-GaTe bilayer is one of the promising systems for optical fine-tunable excitonic devices and provide an ideal platform for realizing strong correlated Bose–Hubbard physics.

Stacking up Electron-Rich and Electron-Deficient Monolayers to Achieve Extraordinary Mid- to Far-Infrared Excitonic Absorption: Interlayer Excitons in the C3B/C3N Bilayer

Our ability to efficiently detect and generate far-infrared (i.e., terahertz) radiation is vital in areas spanning from biomedical imaging to interstellar spectroscopy. Despite decades of intense research, bridging the terahertz gap between electronics and optics remains a major challenge due to the lack of robust materials that can efficiently operate in this frequency range, and two-dimensional (2D) type-II heterostructures may be ideal candidates to fill this gap. Herein, using highly accurate many-body perturbation theory within the GW plus Bethe-Salpeter equation approach, we predict that a type-II heterostructure consisting of an electron-rich ${\mathrm{C}}_{3}\mathrm{N}$ and an electron deficient ${\mathrm{C}}_{3}\mathrm{B}$ monolayers can give rise to extraordinary optical activities in the mid- to far-infrared range. ${\mathrm{C}}_{3}\mathrm{N}$ and ${\mathrm{C}}_{3}\mathrm{B}$ are two graphene-derived 2D materials that have attracted increasing research attention. Although both ${\mathrm{C}}_{3}\mathrm{N}$ and ${\mathrm{C}}_{3}\mathrm{B}$ monolayers are moderate-gap 2D materials, and they couple only through the rather weak van der Waals interactions, the bilayer heterostructure surprisingly supports extremely bright, low-energy interlayer excitons with large binding energies of 0.2--0.4 eV, offering an ideal material with interlayer excitonic states for mid- to far-infrared applications at room temperature. We also investigate in detail the properties and formation mechanism of the inter- and intralayer excitons.

Nonlocal correlation effects due to virtual spin-flip processes in itinerant electron ferromagnets

We present an ab initio method for electronic structure calculations, which accounts for the interaction of electrons and magnons in ferromagnets. While it is based on a many-body perturbation theory we approximate numerically complex quantities with quantities from time-dependent density functional theory. This results in a simple and affordable algorithm which allows us to consider more complex materials than those usually studied in this context ($3d$ ferromagnets) while still being able to account for the nonlocality of the self-energy. Furthermore, our approach allows for a relatively simple way to incorporate self-consistency. Our results are in a good agreement with experimental and theoretical findings for iron and nickel. Especially the experimental exchange splitting of nickel is predicted accurately within our theory. Additionally, we study the half-metallic ferromagnet NiMnSb concerning its nonquasiparticle states appearing in the band gap due to spin-flip excitations.

First-Principle Studies of the Structural, Electronic, and Optical Properties of Double-Walled Carbon Boron Nitride Nanostructures Heterosystem under Various Interwall Distances

Structural, electronic, and optical properties of a new combined system of carbon and boron nitride nanotubes are studied using the DFT first principles as implemented in Quantum ESPRESSO codes. The corrections to the quasi-particle energies were studied via GW hybrid functional implemented in the YAMBO code within the many-body perturbation theory. The studies were performed under different interwall distances of 3.0 nm, 2.5 nm, and 1.5 nm between CNTs and BNNTs. The results showed that the structural properties demonstrated high stability of the double-walled carbon boron nitride nanotube (DWCBNNT) systems under interwall distance (IWD) of 3.00 nm, 2.50 nm, and 1.50 nm. Results also demonstrated an inverse variation between the IWD and the diameter of the DWCBNNT system. In terms of the electronic properties, all three configurations of the DWCBNNTs reveal semiconducting behavior under KS-DFT showing a direct band gap of 3.30 eV, 1.79 eV, and 0.81 eV under IWD of 3.0 nm, 2.5 nm, and 1.5 nm, respectively. Furthermore, the band gap of the DWCBNNT increases with an increase in IWD (decrease in inner tube diameter) and decreases with a decrease in IWD (increase in inner tube diameter). In all three cases, the bands are formed by the molecular orbitals of the armchair CBNNT which are transformed to a series of continuous energy levels; the behaviors of electrons that formed the heterostructure are related to the behavior of electrons in B, C, and N atoms. From the optical properties perspective, the studies were conducted in parallel and perpendicular directions to the nanotubes’ axes. The presence of static dielectric functions in parallel direction at 3.3, 3.4, and 4.5 for nanotubes under 3.0 nm, 2.5 nm, and 1.5 nm demonstrated optical refraction. Refractions were also observed in directions perpendicular to the nanotubes. Furthermore, optical reflections occur when there is a higher absorption. The ability of these CBNNT hybrid systems to refract in all directions revealed the most exciting properties of the armchair CBNNT suitable to be used in magnifying glass materials. The findings further imply that the optical absorption coefficient is inversely related to the diameter of the nanotubes and is directly correlated to the band gap.

Efficient GW calculations in two dimensional materials through a stochastic integration of the screened potential

Many-body perturbation theory methods, such as the G0W0 approximation, are able to accurately predict quasiparticle (QP) properties of several classes of materials. However, the calculation of the QP band structure of two-dimensional (2D) semiconductors is known to require a very dense BZ sampling, due to the sharp q-dependence of the dielectric matrix in the long-wavelength limit (q → 0). In this work, we show how the convergence of the QP corrections of 2D semiconductors with respect to the BZ sampling can be drastically improved, by combining a Monte Carlo integration with an interpolation scheme able to represent the screened potential between the calculated grid points. The method has been validated by computing the band gap of three different prototype monolayer materials: a transition metal dichalcogenide (MoS2), a wide band gap insulator (hBN) and an anisotropic semiconductor (phosphorene). The proposed scheme shows that the convergence of the gap for these three materials up to 50meV is achieved by using k-point grids comparable to those needed by DFT calculations, while keeping the grid uniform.

Assessing the Role of the Kohn-Sham Density in the Calculation of the Low-Lying Bethe-Salpeter Excitation Energies.

We adopt the GW many-body perturbation theory in conjunction with the Bethe-Salpeter equation (BSE) to compute 57 excitation energies of a set of 37 molecules. By using the PBEh global hybrid functional and a self-consistent scheme on the eigenvalues in GW, we show a strong dependence of the BSE energy on the starting Kohn-Sham (KS) density functional. This arises from both the quasiparticle energies and the spatial localization of the frozen KS orbitals employed to compute the BSE. In order to address the arbitrariness in the mean field choice, we adopt an orbital-tuning scheme where the amount of Fock exchange, α, is tuned to impose that the KS HOMO matches the GW quasiparticle eigenvalue, thus fulfilling the ionization potential theorem in DFT. The performance of the proposed scheme yields excellent results and it is similar to M06-2X and PBEh with α = 75%, consistent with tuned values of α ranging between 60% and 80%.

Effects of an external electric field on the electronic properties and optical excitations of germanane and silicane monolayers

Using density functional theory in conjunction with many-body perturbation theory, we theoretically investigated the electronic structures of monolayers germanane and silicane in an applied out-of-plane uniform electric field. Our results show that although the band structures of both monolayers are affected by the electric field, the band gap width cannot be reduced to zero even for high field-strengths. Moreover, excitons are shown to be robust under electric fields, so that Stark shifts for the fundamental exciton peak is only of the order of a few meV for fields of 1 V Å−1. The electric field has also no significant effect on electron probability distribution, as the exciton dissociation into free electron–hole pairs is not observed even at high electric field strengths. Franz-Keldysh effect is also studied in monolayers germanane and silicane. We found that, due to the shielding effect, the external field is prevented to induce absorption in the spectral region below the gap and only above-gap oscillatory spectral features are allowed. One can benefit from such a characteristic where the absorption near the band edge is not altered by the presence of an electric field, especially since these materials have excitonic peaks in the visible range.

Dynamical renormalization of electron-phonon coupling in conventional superconductors

The adiabatic Born-Oppenheimer approximation is considered to be a robust approach that very rarely breaks down. Consequently, it is predominantly utilized to address various electron-phonon properties in condensed matter physics. By combining many-body perturbation and density functional theories we demonstrate the importance of dynamical (nonadiabatic) effects in estimating superconducting properties in various bulk and two-dimensional materials. Apart from the expected long-wavelength nonadiabatic effects, we found sizable nonadiabatic Kohn anomalies away from the Brillouin zone center for materials with strong intervalley electron-phonon scatterings. Compared to the adiabatic result, these dynamical phonon anomalies can significantly modify electron-phonon coupling strength $\ensuremath{\lambda}$ and superconducting transition temperature ${T}_{c}$. Further, the dynamically induced modifications of $\ensuremath{\lambda}$ have a strong impact on transport properties, where probably the most interesting is the rescaling of the low-temperature and low-frequency regime of the scattering time $1/\ensuremath{\tau}$ from about ${T}^{3}$ to about ${T}^{2}$, resembling the Fermi liquid result for electron-electron scattering. Our goal is to point out the potential implications of these nonadiabatic effects and reestablish their pivotal role in computational estimations of electron-phonon properties.

Optical, electronic, and structural properties of ScAlMgO4

Magnesium aluminate scandium oxide (ScAlMgO4) is a promising lattice-matched substrate material for GaN- and ZnO-based optoelectronic devices. Yet, despite its clear advantages over substrates commonly used in heteroepitaxial growth, several fundamental properties of ScAlMgO4 remain unsettled. Here, we provide a comprehensive picture of its optical, electronic and structural properties by studying ScAlMgO4 single crystals grown by the Czochralski method. We use variable angle spectroscopic ellipsometry to determine complex in-plane and out-of-plane refractive indices in the range from 193 to 1690 nm. An oscillator-based model provides a phenomenological description of the ellipsometric spectra with excellent agreement over the entire range of wavelengths. For convenience, we supply the reader also with Cauchy formulas describing the real part of the anisotropic refractive index for wavelengths above 400 nm. Ab initio many-body perturbation theory modeling provides information about the electronic structure of ScAlMgO4, and successfully validated experimentally obtained refractive index values. Simulations also show exciton binding energy as large as a few hundred of meV, indicating ScAlMgO4 as a promising material for implementation in low-threshold, deep-UV lasing devices operating at room temperature. X-ray diffraction measurements confirm lattice constants of ScAlMgO4 previously reported, but in addition, reveal that dominant crystallographic planes (001) are mutually inclined by about 0.009{\deg}. In view of our work, ScAlMgO4 is a highly transparent, low refractive index, birefringent material similar to a sapphire, but with a much more favorable lattice constant and simpler processing.

Anharmonic lattice dynamics and the origin of intrinsic ultralow thermal conductivity in AgI materials

Ionic conductors such as AgI with ultralow thermal conductivities $({\ensuremath{\kappa}}_{l})$ are of increasing interest because of their excellent thermoelectric properties. However, the origin of their intrinsic low ${\ensuremath{\kappa}}_{l}$ values remain elusive. In this study, comprehensive theoretical calculations of the lattice dynamics and the thermal transport properties of $\ensuremath{\gamma}\text{\ensuremath{-}}\mathrm{AgI}$ (zinc-blende structure) and $\ensuremath{\beta}\text{\ensuremath{-}}\mathrm{AgI}$ (wurtzite structure) as functions of temperature were carried out based on many-body perturbation theory and phonon Boltzmann transport theory. First, the mean-squared displacements (MSDs) of ${\mathrm{Ag}}^{+}$ were significantly larger than those of ${\mathrm{I}}^{\ensuremath{-}}$ in both $\ensuremath{\gamma}\text{\ensuremath{-}}$ and \ensuremath{\beta}-phases below the order-disorder phase transition temperature $({T}_{\mathrm{c}})$, which led to a characteristic ``rattling'' feature and low-frequency, nearly flat local phonon vibrations. According to our previous work [Xie et al., Phys. Rev. Lett. 125, 245901 (2020)], such nondispersive flat phonon band structures are expected to give rise to four-phonon resonance and result in a dramatic increase in the four-phonon scattering over the conventional three-phonon scattering. For $\ensuremath{\gamma}\text{\ensuremath{-}}\mathrm{AgI}$, similar four-phonon resonance behavior was also discovered for the low-lying transverse acoustic phonon branches, and it was found that their four-phonon scattering rates were an order of magnitude larger than the corresponding three-phonon scattering rates. Considering the four-phonon scattering, the theoretical ${\ensuremath{\kappa}}_{l}$ of $\ensuremath{\gamma}\text{\ensuremath{-}}\mathrm{AgI}$ was predicted to be $\ensuremath{\sim}0.32$ W/m K at 300 K, which was in good agreement with the value deduced from our experiments ($\ensuremath{\sim}0.36$ W/m K at 300 K). Compared to $\ensuremath{\gamma}\text{\ensuremath{-}}\mathrm{AgI}$, the acoustic phonons in $\ensuremath{\beta}\text{\ensuremath{-}}\mathrm{AgI}$ were more dispersive, and they intertwined with low-energy optical phonons at the zone boundaries. It was found that three-phonon resonance became as important as four-phonon resonance for the nearly flat longitudinal phonon band. The theoretical ${\ensuremath{\kappa}}_{l}$ for $\ensuremath{\beta}\text{\ensuremath{-}}\mathrm{AgI}$ was determined to be around $\ensuremath{\sim}0.32$ W/m K at room temperature, closely reproducing our measurement value $\ensuremath{\sim}0.29$ W/m K. Our results for AgI demonstrate the strong quartic anharmonicity in materials characterized by the rattling of weak bonding atoms as well as dispersionless phonon band structures. It is believed that this intimate relationship between the low-${\ensuremath{\kappa}}_{l}$ and flat phonon dispersion can be employed as a good indicator when searching for material systems with ultralow ${\ensuremath{\kappa}}_{l}$ values, e.g., cagelike rattling structures, quasi-two-dimensional structures, and chainlike structures.

Electronic and Optical Properties of Alkaline Earth Metal Fluoride Crystals with the Inclusion of Many-Body Effects: A Comparative Study on Rutile MgF2 and Cubic SrF2

We conducted a systematic investigation using state-of-the-art techniques on the electronic and optical properties of two crystals of alkaline earth metal fluorides, namely rutile MgF2 and cubic SrF2. For these two crystals of different symmetry, we present density functional theory (DFT), many-body perturbation theory (MBPT), and Bethe–Salpeter equation (BSE) calculations. We calculated a variety of properties, namely ground-state energies, band-energy gaps, and optical absorption spectra with the inclusion of excitonic effects. The quantities were obtained with a high degree of convergence regarding all bulk electronic and optical properties. Bulk rutile MgF2 has distinguished ground-state and excited-state properties with respect to the other cubic fluoride SrF2 and the other members of the alkaline earth metal fluoride family. The nature of the fundamental gaps and estimates of the self-energy and excitonic effects for the two compounds are presented and discussed in detail. Our results are in good accordance with the measurements and other theoretical–computational data. A comparison is made between the excitation and optical properties of bulk rutile MgF2, cubic SrF2, and the corresponding clusters, for which calculations have recently been published, confirming strong excitonic effects in finite-sized systems.

On the Chemical Potential of Many-Body Perturbation Theory in Extended Systems.

Finite-temperature many-body perturbation theory in the grand-canonical ensemble is fundamental to numerous methods for computing electronic properties at nonzero temperature, such as finite-temperature coupled-cluster. In most applications it is the average number of electrons that is known rather than the chemical potential. Expensive correlation calculations must be repeated iteratively in search for the interacting chemical potential that yields the desired average number of electrons. In extended systems with mobile charges the situation is particular, however. Long-ranged electrostatic forces drive the charges such that the average ratio of negative and positive charges is one for any finite chemical potential. All properties per electron are expected to be virtually independent of the chemical potential, as they are in an electric wire at different voltage potentials. This work shows that per electron, the exchange-correlation free energy and the exchange-correlation grand potential indeed agree in the infinite-size limit. Thus, only one expensive correlation calculation suffices for each system size, sparing the search for the interacting chemical potential. This work also demonstrates the importance of regularizing the Coulomb interaction such that each electron on average interacts only with as many electrons as there are electrons in the simulation, avoiding interactions with periodic images. Numerical calculations of the warm uniform electron gas have been conducted with the Spencer-Alavi regularization employing the finite-temperature Hartree approximation for the self-consistent field and linearized finite-temperature direct-ring coupled-cluster doubles for treating correlation.