Search

The Bethe-Salpeter Equation (BSE) can be applied to compute from first-principles optical spectra that include the effects of screened electron-hole interactions. As input, BSE calculations require single-particle states, quasiparticle energy levels and the screened Coulomb interaction, which are typically obtained with many-body perturbation theory, whose cost limits the scope of possible applications. This work tries to address this practical limitation, instead deriving spectral energies from Koopmans-compliant functionals and introducing a new methodology for handling the screened Coulomb interaction. The explicit calculation of the $W$ matrix is bypassed via a direct minimization scheme applied on top of a maximally localised Wannier function basis. We validate and benchmark this approach by computing the low-lying excited states of the molecules in Thiel's set, and the optical absorption spectrum of a $\text{C}_{60}$ fullerene. The results show the same trends as quantum chemical methods and are in excellent agreement with previous simulations carried out at the TD-DFT or $G_{0}W_{0}$-\text{BSE} level. Conveniently, the new framework reduces the parameter space controlling the accuracy of the calculation, thereby simplifying the simulation of charge-neutral excitations, offering the potential to expand the applicability of first-principles spectroscopies to larger systems of applied interest.

We present SIMPLE, a code developed to calculate optical properties of metallic and insulating extended systems using the optimal basis method originally proposed by E. L. Shirley in 1996. Two different approaches for the evaluation of the complex dielectric function are implemented: (i) the independent-particle approximation considering both interband and intraband contributions for metals and (ii) the Bethe-Salpeter equation for insulators. Since, notably, the optimal basis set is systematically improvable, accurate results can be obtained at a strongly reduced computational cost. Moreover, the simplicity of the method allows for a straightforward use of the code; improvement of the optimal basis and thus the overall precision of the simulations is simply controlled by one (for metals) or two (for insulators) input thresholds. The code is extensively tested, in terms of verification and performance, on bulk silver for metals and on bulk silicon for insulators.

Excitonic effects due to electron-hole coupling play a fundamental role in renormalising energy levels in dye sensitised and organic solar cells determining the driving force for electron extraction. We show that first-principles calculations based on many-body perturbation theory within the GW-BSE approach provide a quantitative picture of interfacial excited state energetics in organic dye-sensitized TiO2 , delivering a general rule for evaluating relevant energy levels.To perform GW-BSE calculations in such large systems we introduce a new scheme based on maximally localized Wannier functions. With this method the overall scaling of GW-BSE calculations is reduced from O(N4) to O(N3).

Thermal transport coefficients are independent of the specific microscopic expression for the energy density and current from which they can be derived through the Green-Kubo formula. We discuss this independence in terms of a kind of gauge invariance resulting from energy conservation and extensivity, and demonstrate it numerically for a Lennard-Jones fluid, where different forms of the microscopic energy density lead to different time correlation functions for the heat flux, all of them, however, resulting in the same value for the thermal conductivity., Comment: 7 pages, 3 figures

Quantum simulation methods based on density-functional theory are currently deemed unfit to cope with atomic heat transport within the Green-Kubo formalism, because quantum-mechanical energy densities and currents are inherently ill-defined at the atomic scale. We show that, while this difficulty would also affect classical simulations, thermal conductivity is indeed insensitive to such ill-definedness by virtue of a sort of gauge invariance resulting from energy extensivity and conservation. Based on these findings, we derive an expression for the adiabatic energy flux from density-functional theory, which allows heat transport to be simulated using ab-initio equilibrium molecular dynamics. Our methodology is demonstrated by comparing its predictions with those of classical equilibrium and ab-initio non-equilibrium (M\"uller-Plathe) simulations for a liquid-Argon model, and finally applied to heavy water at ambient conditions., Comment: 9 two-column pages, 4 figures

Raman spectroscopy is a useful experimental tool to investigate local deformation and structural changes in SiO$_2$-based glasses. Using a semi-classical modelling of Raman spectra in large samples of silica glasses, we show in this paper that shear plastic flow affects the Raman measurement in the upper part of the spectrum. We relate these changes to structural modifications, as well as a detailed analysis of the vibration modes computed in the same frequency range. These results opens the door to in situ monitoring of plastic damage in silica-based structures.

We present a first-principles investigation of the excited-state properties of electron acceptors in organic photovoltaics including C$_{60}$, C$_{70}$, [6,6]-phenyl-C$_{61}$-butyric-acid-methyl-ester ([C$_{60}$]PCBM), and bis-[C$_{60}$]PCBM using many-body perturbation theory within the Hedin's $G_0W_0$ approximation and an efficient Lanczos approach. Calculated vertical ionization potentials (VIP) and vertical electron affinities (VEA) of C$_{60}$ and C$_{70}$ agree very well with experimental values measured in gas phase. The density of states of all three molecules is also compared to photoemission and inverse photoemission spectra measured on thin-films, exhibiting a close agreement - a rigid energy-gap renormalization owing to intermolecular interactions in the thin-films. In addition, it is shown that the low-lying unoccupied states of [C$_{60}$]PCBM are all derived from the highest-occupied molecular orbitals and the lowest-unoccupied molecular orbitals of fullerene C$_{60}$. The functional side group in [C$_{60}$]PCBM introduces a slight electron transfer to the fullerene cage, resulting in small decreases of both VIP and VEA. This small change of VEA provides a solid justification for the increase of open-circuit voltage when replacing fullerene C$_{60}$ with [C$_{60}$]PCBM as the electron acceptor in bulk heterojunction polymer solar cells., Comment: 9 pages, 4 figures, and 7 tables

Hybrid AMX3 perovskites (A=Cs, CH3NH3; M=Sn, Pb; X=halide) have revolutionized the scenario of emerging photovoltaic technologies. Introduced in 2009 by Kojima et al., a rapid evolution very recently led to 15% efficient solar cells. CH3NH3PbI3 has so far dominated the field, while the similar CH3NH3SnI3 has not been explored for photovoltaic applications, despite the reduced band-gap. Replacement of Pb by the more environment-friendly Sn would facilitate the large uptake of perovskite-based photovoltaics. Despite the extremely fast progress, the materials electronic properties which are key to the photovoltaic performance are relatively little understood. Here we develop an effective GW method incorporating spin-orbit coupling which allows us to accurately model the electronic, optical and transport properties of CH3NH3SnI3 and CH3NH3PbI3, opening the way to new materials design. The different CH3NH3SnI3 and CH3NH3PbI3 properties are discussed in light of their exploitation for solar cells, and found to be entirely due to relativistic effects., Comment: 16 pages, 4 figures

The photoelectron properties of DNA and RNA bases are studied using many-body perturbation theory within the GW approximation, together with a recently developed Lanczos-chain approach. Calculated vertical ionization potentials, electron affinities, and total density of states are in good agreement with experimental values and photoemission spectra. The convergence benchmark demonstrates the importance of using an optimal polarizability basis in the GW calculations. A detailed analysis of the role of exchange and correlation in both many-body and density-functional theory calculations shows that while self-energy corrections are strongly orbital-dependent, they nevertheless remain almost constant for states that share the same bonding character. Finally, we report on the inverse lifetimes of DNA and RNA bases, that are found to depend linearly on quasi-particle energies for all deep valence states. In general, our G0W0-Lanczos approach provides an efficient yet accurate and fully converged description of quasiparticle properties of five DNA and RNA bases.

We calculate the linear and non-linear susceptibilities of periodic longitudinal chains of hydrogen dimers with different bond-length alternations using a diffusion quantum Monte Carlo approach. These quantities are derived from the changes in electronic polarization as a function of applied finite electric field - an approach we recently introduced and made possible by the use of a Berry-phase, many-body electric-enthalpy functional. Calculated susceptibilities and hyper-susceptibilities are found to be in excellent agreement with the best estimates available from quantum chemistry - usually extrapolations to the infinite-chain limit of calculations for chains of finite length. It is found that while exchange effects dominate the proper description of the susceptibilities, second hyper-susceptibilities are greatly affected by electronic correlations. We also assess how different approximations to the nodal surface of the many-body wavefunction affect the accuracy of the calculated susceptibilities., Comment: 9 pages; accepted on J. Chem. Phys