Showing 1–44 of 44 results for author: Varsano, D

Distinguishing different stackings in layered materials via luminescence spectroscopy

Authors: Matteo Zanfrognini, Alexandre Plaud, Ingrid Stenger, Frédéric Fossard, Lorenzo Sponza, Léonard Schué, Fulvio Paleari, Elisa Molinari, Daniele Varsano, Ludger Wirtz, François Ducastelle, Annick Loiseau, Julien Barjon

Abstract: Despite its simple crystal structure, layered boron nitride features a surprisingly complex variety of phonon-assisted luminescence peaks. We present a combined experimental and theoretical study on ultraviolet-light emission in hexagonal and rhombohedral bulk boron nitride crystals. Emission spectra of high-quality samples are measured via cathodoluminescence spectroscopy, displaying characterist… ▽ More Despite its simple crystal structure, layered boron nitride features a surprisingly complex variety of phonon-assisted luminescence peaks. We present a combined experimental and theoretical study on ultraviolet-light emission in hexagonal and rhombohedral bulk boron nitride crystals. Emission spectra of high-quality samples are measured via cathodoluminescence spectroscopy, displaying characteristic differences between the two polytypes. These differences are explained using a fully first-principles computational technique that takes into account radiative emission from ``indirect'', finite-momentum, excitons via coupling to finite-momentum phonons. We show that the differences in peak positions, number of peaks and relative intensities can be qualitatively and quantitatively explained, once a full integration over all relevant momenta of excitons and phonons is performed. △ Less

Submitted 27 May, 2023; originally announced May 2023.

Comments: Main: 6 pages and 4 figures, Supplementary: 6 pages and 7 figures

Efficient GW calculations via the interpolation of the screened interaction in momentum and frequency space: The case of graphene

Authors: Alberto Guandalini, Dario A. Leon, Pino D'Amico, Claudia Cardoso, Andrea Ferretti, Massimo Rontani, Daniele Varsano

Abstract: The GW self-energy may become computationally challenging to evaluate because of frequency and momentum convolutions. These difficulties were recently addressed by the development of the multipole approximation (MPA) and the W-av methods: MPA accurately approximates full-frequency response functions using a small number of poles, while W-av improves the convergence with respect to the k-point samp… ▽ More The GW self-energy may become computationally challenging to evaluate because of frequency and momentum convolutions. These difficulties were recently addressed by the development of the multipole approximation (MPA) and the W-av methods: MPA accurately approximates full-frequency response functions using a small number of poles, while W-av improves the convergence with respect to the k-point sampling in 2D materials. In this work we (i) present a theoretical scheme to combine them, and (ii) apply the newly developed approach to the paradigmatic case of graphene. Our findings show an excellent agreement of the calculated QP band structure with angle resolved photoemission spectroscopy (ARPES) data. Furthermore, the computational efficiency of MPA and W-av allows us to explore the logarithmic renormalization of the Dirac cone. To this aim, we develop an analytical model, derived from a Dirac Hamiltonian, that we parameterize using ab-initio data. The comparison of the models obtained with PPA and MPA results highlights an important role of the dynamical screening in the cone renormalization. △ Less

Submitted 9 October, 2023; v1 submitted 21 April, 2023; originally announced April 2023.

Comments: A. Guandalini and D. A. Leon contributed equally to this work

Quenching of low-energy optical absorption in bilayer C$_3$N polytypes

Authors: Matteo Zanfrognini, Miki Bonacci, Fulvio Paleari, Elisa Molinari, Alice Ruini, Andrea Ferretti, Marilia J. Caldas, Daniele Varsano

Abstract: In this work we provide a first principles description of the electronic and optical properties of bilayers C$_3$N, with different stacking motifs AB, AB$'$ and AA$'$. Starting from quasi-particle electronic band-structures, we solve the Bethe Salpeter Equation (BSE) to access the excitonic properties of these bilayers. For all stacking sequences, we see strong optical absorption at energies lower… ▽ More In this work we provide a first principles description of the electronic and optical properties of bilayers C$_3$N, with different stacking motifs AB, AB$'$ and AA$'$. Starting from quasi-particle electronic band-structures, we solve the Bethe Salpeter Equation (BSE) to access the excitonic properties of these bilayers. For all stacking sequences, we see strong optical absorption at energies lower than but close to that of the monolayer. Most relevant, we predict a strong quenching of the low-energy optical absorption, with negligible oscillator strength of low-lying bound excitons. This is a unique phenomenology that does not arise in the monolayer case, nor in other common homo-bilayers. We explain these findings in terms of the small interband dipole matrix elements associated to the valence-conduction transitions involved in these excitons, and discuss them in view of the different stacking motifs. △ Less

Submitted 18 February, 2023; originally announced February 2023.

Comments: 13 pages, 8 figures

Excitonic effects in energy loss spectra of freestanding graphene

Authors: Alberto Guandalini, Ryosuke Senga, Yung-Chang Lin, Kazu Suenaga, Andrea Ferretti, Daniele Varsano, Andrea Recchia, Paolo Barone, Francesco Mauri, Thomas Pichler, Christian Kramberger

Abstract: In this work we perform electron energy-loss spectroscopy (EELS) of freestanding graphene with high energy and momentum resolution to disentangle the quasielastic scattering from the excitation gap of Dirac electrons close to the optical limit. We show the importance of many-body effects on electronic excitations at finite transferred momentum by comparing measured EELS with ab initio calculations… ▽ More In this work we perform electron energy-loss spectroscopy (EELS) of freestanding graphene with high energy and momentum resolution to disentangle the quasielastic scattering from the excitation gap of Dirac electrons close to the optical limit. We show the importance of many-body effects on electronic excitations at finite transferred momentum by comparing measured EELS with ab initio calculations at increasing levels of theory. Quasi-particle corrections and excitonic effects are addressed within the GW approximation and Bethe-Salpeter equation, respectively. Both effects are essential in the description of the EEL spectra to obtain a quantitative agreement with experiments, with the position, dispersion, and shape of both the excitation gap and the $π$ plasmon being significantly affected by excitonic effects. △ Less

Submitted 10 October, 2023; v1 submitted 13 February, 2023; originally announced February 2023.

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

Authors: Miki Bonacci, Junfeng Qiao, Nicola Spallanzani, Antimo Marrazzo, Giovanni Pizzi, Elisa Molinari, Daniele Varsano, Andrea Ferretti, Deborah Prezzi

Abstract: 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) calc… ▽ More 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, a novel 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. △ Less

Submitted 18 January, 2023; v1 submitted 16 January, 2023; originally announced January 2023.

Comments: 14 pages, 6 figures

Journal ref: npj Computational Materials volume 9, Article number: 74 (2023)

Efficient full frequency GW for metals using a multipole approach for the dielectric screening

Authors: Dario A. Leon, Andrea Ferretti, Daniele Varsano, Elisa Molinari, Claudia Cardoso

Abstract: The properties of metallic systems with important and structured excitations at low energies, such as Cu, are challenging to describe with simple models like the plasmon pole approximation (PPA), and more accurate and sometimes prohibitive full frequency approaches are usually required. In this paper we propose a numerical approach to $GW$ calculations on metals that takes into account the frequen… ▽ More The properties of metallic systems with important and structured excitations at low energies, such as Cu, are challenging to describe with simple models like the plasmon pole approximation (PPA), and more accurate and sometimes prohibitive full frequency approaches are usually required. In this paper we propose a numerical approach to $GW$ calculations on metals that takes into account the frequency dependence of the screening via the multipole approximation (MPA), an accurate and efficient alternative to current full-frequency methods that was recently developed and validated for semiconductors and overcomes several limitations of PPA. We now demonstrate that MPA can be successfully extended to metallic systems by optimizing the frequency sampling for this class of materials and introducing a simple method to include the $\mathbf{q}\to 0$ limit of the intra-band contributions. The good agreement between MPA and full frequency results for the calculations of quasi-particle energies, polarizability, self-energy and spectral functions in different metallic systems confirms the accuracy and computational efficiency of the method. Finally, we discuss the physical interpretation of the MPA poles through a comparison with experimental electron energy loss spectra for Cu. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comments: 21 pages

Graphene decoupling through oxygen intercalation on Gr/Co and Gr/Co/Ir interfaces

Authors: Dario A. Leon, Andrea Ferretti, Daniele Varsano, Elisa Molinari, Claudia Cardoso

Abstract: We perform a density functional theory study of the effects of oxygen adsorption on the structural and electronic properties of Gr/Co(0001) and Gr/Co/Ir(111) interfaces. In both interfaces, the graphene-Co distance increases with increasing O concentration. The oxygen intercalation effectively decreases the electronic interaction, preventing the hybridization of graphene states with Co $d$-orbital… ▽ More We perform a density functional theory study of the effects of oxygen adsorption on the structural and electronic properties of Gr/Co(0001) and Gr/Co/Ir(111) interfaces. In both interfaces, the graphene-Co distance increases with increasing O concentration. The oxygen intercalation effectively decreases the electronic interaction, preventing the hybridization of graphene states with Co $d$-orbitals, hence (partly) restoring the typical Dirac cone of pristine graphene. In the case of graphene/Co 1ML/Ir(111), which presents a moiré pattern, the interplay between the O distribution and the continuous change of the graphene-Co registry can be used to tune graphene corrugation and electronic properties. The computed electronic properties are in very good agreement with previously reported angle resolved photoemission spectroscopy and photoemission electron microscopy measurements for Gr/Co(0001). △ Less

Submitted 4 August, 2022; originally announced August 2022.

Comments: 15 pages, 15 figures

Journal ref: Phys. Rev. Materials 6, 064004 (2022)

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

Authors: Alberto Guandalini, Pino D'Amico, Andrea Ferretti, Daniele Varsano

Abstract: Many-body perturbation theory methods, such as the $G_0W_0$ 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 (… ▽ More Many-body perturbation theory methods, such as the $G_0W_0$ 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 ($\mathbf{q} \to 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 (MoS$_2$), 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. △ Less

Submitted 22 June, 2022; v1 submitted 24 May, 2022; originally announced May 2022.

Excitonic effects in graphene-like C$_3$N

Authors: Miki Bonacci, Matteo Zanfrognini, Elisa Molinari, Alice Ruini, Marilia J. Caldas, Andrea Ferretti, Daniele Varsano

Abstract: Monolayer C$_3$N is an emerging two-dimensional indirect band gap semiconductor with interesting mechanical, thermal, and electronic properties. In this work we present a description of C$_3$N electronic and dielectric properties, focusing on the so-called momentum-resolved exciton band structure. Excitation energies and oscillator strengths are computed in order to characterize bright and dark st… ▽ More Monolayer C$_3$N is an emerging two-dimensional indirect band gap semiconductor with interesting mechanical, thermal, and electronic properties. In this work we present a description of C$_3$N electronic and dielectric properties, focusing on the so-called momentum-resolved exciton band structure. Excitation energies and oscillator strengths are computed in order to characterize bright and dark states, and discussed also with respect to the crystal symmetry. Activation of excitonic states is observed for finite transferred momenta: Indeed, we find an active indirect exciton at $\sim$ 0.9 eV, significantly lower than the direct optical gap of 1.96 eV, with excitonic binding energies in the range 0.6-0.9 eV for the lowest states. As for other 2D materials, we find a quasi-linear excitonic dispersion close to $Γ$, which however shows a downward convexity related to the indirect band gap of C$_3$N as well as to the dark nature of the involved excitons. △ Less

Submitted 23 April, 2022; originally announced April 2022.

Comments: 8 pages; 7 figures

Journal ref: Phys. Rev. Materials, 6, 3, 034009 (2022)

Anomalous screening in narrow-gap carbon nanotubes

Authors: Giacomo Sesti, Daniele Varsano, Elisa Molinari, Massimo Rontani

Abstract: The screening of Coulomb interaction controls many-body physics in carbon nanotubes, as it tunes the range and strength of the force that acts on charge carriers and binds electron-hole pairs into excitons. In doped tubes, the effective Coulomb interaction drives the competition between Luttinger liquid and Wigner crystal, whereas in undoped narrow-gap tubes it dictates the Mott or excitonic natur… ▽ More The screening of Coulomb interaction controls many-body physics in carbon nanotubes, as it tunes the range and strength of the force that acts on charge carriers and binds electron-hole pairs into excitons. In doped tubes, the effective Coulomb interaction drives the competition between Luttinger liquid and Wigner crystal, whereas in undoped narrow-gap tubes it dictates the Mott or excitonic nature of the correlated insulator observed at low temperature. Here, by computing the dielectric function of selected narrow- and zero-gap tubes from first principles, we show that the standard effective-mass model of screening systematically underestimates the interaction strength at long wavelength, hence missing the binding of low-energy excitons. The reason is that the model critically lacks the full three-dimensional topology of the tube, being adapted from graphene theory. As ab inito calculations are limited to small tubes, we develop a two-band model dielectric function based on the plane-wave expansion of Bloch states and the exact truncated Coulomb cutoff technique. We demonstrate that our -- computationally cheap -- approach provides the correct screening for narrow-gap tubes of any size and chirality. A striking result is that the screened interaction remains long-ranged even in gapless tubes, as an effect of the microscopic local fields generated by the electrons moving on the curved tube surface. As an application, we show that the effective electron-electron force that is felt at distances relevant to quantum transport experiments is super Coulombic. △ Less

Submitted 14 April, 2022; v1 submitted 25 January, 2022; originally announced January 2022.

Non-adiabatic suppression of 3D excitonic screening in black phosphorus by mid-infrared pulses

Authors: Angela Montanaro, Francesca Giusti, Matteo Zanfrognini, Paola Di Pietro, Filippo Glerean, Giacomo Jarc, Enrico Maria Rigoni, Shahla Y. Mathengattil, Daniele Varsano, Massimo Rontani, Andrea Perucchi, Elisa Molinari, Daniele Fausti

Abstract: The competition between the electron-hole Coulomb attraction and the three-dimensional dielectric screening dictates the optical properties of layered semiconductors. In low-dimensional materials, the equilibrium dielectric environment can be significantly altered by the ultrafast excitation of photo-carriers, leading to renormalized band gap and exciton binding energies. Recently, black phosphoru… ▽ More The competition between the electron-hole Coulomb attraction and the three-dimensional dielectric screening dictates the optical properties of layered semiconductors. In low-dimensional materials, the equilibrium dielectric environment can be significantly altered by the ultrafast excitation of photo-carriers, leading to renormalized band gap and exciton binding energies. Recently, black phosphorus emerged as a 2D material with strongly layer-dependent electronic properties. Here, we resolve the coherent response of screening to sub-gap photo-excitation in bulk black phosphorus and find that mid-infrared pulses tuned across the band gap drive a transient non-thermal suppression of the dielectric screening, which is revealed by the emergence of the single-layer exciton resonance. Our work exposes the role of interlayer interactions in determining the electronic properties of 2D materials and discloses the possibility of optically manipulate them, which is of great relevance for the engineering of versatile van der Waals low-dimensional materials. △ Less

Submitted 24 November, 2021; originally announced November 2021.

Evidence for equilibrium excitons and exciton condensation in monolayer WTe2

Authors: Bosong Sun, Wenjin Zhao, Tauno Palomaki, Zaiyao Fei, Elliott Runburg, Paul Malinowski, Xiong Huang, John Cenker, Yong-Tao Cui, Jiun-Haw Chu, Xiaodong Xu, S. Samaneh Ataei, Daniele Varsano, Maurizia Palummo, Elisa Molinari, Massimo Rontani, David H. Cobden

Abstract: A single monolayer of the layered semimetal WTe2 behaves as a two-dimensional topological insulator, with helical conducting edge modes surrounding a bulk state that becomes insulating at low temperatures. Here we present evidence that the bulk state has a very unusual nature, containing electrons and holes bound by Coulomb attraction (excitons) that spontaneously form in thermal equilibrium. On c… ▽ More A single monolayer of the layered semimetal WTe2 behaves as a two-dimensional topological insulator, with helical conducting edge modes surrounding a bulk state that becomes insulating at low temperatures. Here we present evidence that the bulk state has a very unusual nature, containing electrons and holes bound by Coulomb attraction (excitons) that spontaneously form in thermal equilibrium. On cooling from room temperature to 100 K the conductivity develops a V-shaped dependence on electrostatic doping, while the chemical potential develops a ~43 meV step at the neutral point. These features are much sharper than is possible in an independent-electron picture, but they can be largely accounted for by positing that some of the electrons and holes are paired in equilibrium. Our calculations from first principles show that the exciton binding energy is larger than 100 meV and the radius as small as 4 nm, explaining their formation at high temperature and doping levels. Below 100 K more strongly insulating behavior is seen, suggesting that a charge-ordered state forms. The observed absence of charge density waves in this state appears surprising within an excitonic insulator picture, but we show that it can be explained by the symmetries of the exciton wave function. Monolayer WTe2 therefore presents an exceptional combination of topological properties and strong correlations over a wide temperature range. △ Less

Submitted 26 October, 2021; originally announced October 2021.

Comments: Initial submission. Includes supplementary information. Revised version to appear in Nature Physics

Frequency dependence in GW made simple using a multi-pole approximation

Authors: Dario A. Leon, Claudia Cardoso, Tommaso Chiarotti, Daniele Varsano, Elisa Molinari, Andrea Ferretti

Abstract: In the $GW$ approximation, the screened interaction $W$ is a non-local and dynamical potential that usually has a complex frequency dependence. A full description of such dependence is possible but often computationally demanding. For this reason, it is still common practice to approximate $W(ω)$ using a plasmon pole (PP) model. Such approach, however, may deliver an accuracy limited by its simpli… ▽ More In the $GW$ approximation, the screened interaction $W$ is a non-local and dynamical potential that usually has a complex frequency dependence. A full description of such dependence is possible but often computationally demanding. For this reason, it is still common practice to approximate $W(ω)$ using a plasmon pole (PP) model. Such approach, however, may deliver an accuracy limited by its simplistic description of the frequency dependence of the polarizability, i.e. of $W$. In this work we explore a multi-pole approach (MPA) and develop an effective representation of the frequency dependence of $W$. We show that an appropriate sampling of the polarizability in the frequency complex plane and a multi-pole interpolation can lead to a level of accuracy comparable with full-frequency methods at much lower computational cost. Moreover, both accuracy and cost are controllable by the number of poles used in MPA. Eventually we validate the MPA approach in selected prototype systems, showing that full-frequency quality results can be obtained with a limited number of poles. △ Less

Submitted 3 September, 2021; originally announced September 2021.

Comments: 7 figures, 18 pages + 2 supplementary pages

Journal ref: Phys. Rev. B 104, 115157 (2021)

Narrowing of d bands of FeCo layers intercalated under graphene

Authors: Daniela Pacilè, Claudia Cardoso, Giulia Avvisati, Ivana Vobornik, Carlo Mariani, Dario A. Leon, Pietro Bonfà, Daniele Varsano, Andrea Ferretti, Maria Grazia Betti

Abstract: We report on the electronic properties of an artificial system obtained by the intercalation of equiatomic FeCo layers under graphene grown on Ir(111). Upon intercalation, the FeCo film grows epitaxially on Ir(111), resulting in a lattice-mismatched system. By performing Density Functional Theory calculations, we show that the intercalated FeCo layer leads to a pronounced corrugation of the graphe… ▽ More We report on the electronic properties of an artificial system obtained by the intercalation of equiatomic FeCo layers under graphene grown on Ir(111). Upon intercalation, the FeCo film grows epitaxially on Ir(111), resulting in a lattice-mismatched system. By performing Density Functional Theory calculations, we show that the intercalated FeCo layer leads to a pronounced corrugation of the graphene film. At the same time, the FeCo intercalated layers induce a clear transition from a nearly undisturbed to a strongly hybridized graphene π-band, as measured by angle-resolved photoemission spectroscopy. A comparison of experimental results with the computed band structure and the projected density of states unveils a spin-selective hybridization between the π band of graphene and FeCo-3d states. Our results demonstrate that the reduced dimensionality, as well as the hybridization within the FeCo layers, induce a narrowing and a clear splitting of Fe 3d-up and Fe 3d-down spin bands of the confined FeCo layers with respect to bulk Fe and Co. △ Less

Submitted 29 March, 2021; originally announced March 2021.

Comments: 5 pages, 4 figures

Journal ref: Appl. Phys. Lett. 118, 121602, 2021

Magnetic response and electronic states of well defined Graphene/Fe/Ir(111) heterostructure

Authors: Claudia Cardoso, Giulia Avvisati, Pierluigi Gargiani, Marco Sbroscia, Madan S. Jagadeesh, Carlo Mariani, Dario A. Leon, Daniele Varsano, Andrea Ferretti, Maria Grazia Betti

Abstract: We investigate a well defined heterostructure constituted by magnetic Fe layers sandwiched between graphene (Gr) and Ir(111). The challenging task to avoid Fe-C solubility and Fe-Ir intermixing has been achieved with atomic controlled Fe intercalation at moderate temperature below 500 K. Upon intercalation of a single ordered Fe layer in registry with the Ir substrate, an intermixing of the Gr ban… ▽ More We investigate a well defined heterostructure constituted by magnetic Fe layers sandwiched between graphene (Gr) and Ir(111). The challenging task to avoid Fe-C solubility and Fe-Ir intermixing has been achieved with atomic controlled Fe intercalation at moderate temperature below 500 K. Upon intercalation of a single ordered Fe layer in registry with the Ir substrate, an intermixing of the Gr bands and Fe d states breaks the symmetry of the Dirac cone, with a downshift in energy of the apex by about 3 eV, and well-localized Fe intermixed states induced in the energy region just below the Fermi level. First principles electronic structure calculations show a large spin splitting of the Fe states, resulting in a majority spin channel almost fully occupied and strongly hybridized with Gr π states. X-ray magnetic circular dichroism on the Gr/Fe/Ir heterostructure reveals an ordered spin configuration with a ferromagnetic response of Fe layer(s), with enhanced spin and orbital configurations with respect to the bcc-Fe bulk values. The magnetization switches from a perpendicular easy magnetization axis when the Fe single layer is lattice matched with the Ir(111) surface to a parallel one when the Fe thin film is almost commensurate with graphene. △ Less

Submitted 14 January, 2021; originally announced January 2021.

Journal ref: Physical Review Materials 5, 014405 (2021)

Evidence of ideal excitonic insulator in bulk MoS2 under pressure

Authors: S. Samaneh Ataei, Daniele Varsano, Elisa Molinari, Massimo Rontani

Abstract: Spontaneous condensation of excitons is a long sought phenomenon analogous to the condensation of Cooper pairs in a superconductor. It is expected to occur in a semiconductor at thermodynamic equilibrium if the binding energy of the excitons---electron (e) and hole (h) pairs interacting by Coulomb force---overcomes the band gap, giving rise to a new phase: the 'excitonic insulator' (EI). Transitio… ▽ More Spontaneous condensation of excitons is a long sought phenomenon analogous to the condensation of Cooper pairs in a superconductor. It is expected to occur in a semiconductor at thermodynamic equilibrium if the binding energy of the excitons---electron (e) and hole (h) pairs interacting by Coulomb force---overcomes the band gap, giving rise to a new phase: the 'excitonic insulator' (EI). Transition metal dichalcogenides are excellent candidates for the EI realization because of reduced Coulomb screening, and indeed a structural phase transition was observed in few-layer systems. However, previous work could not disentangle to which extent the origin of the transition was in the formation of bound excitons or in the softening of a phonon. Here we focus on bulk MoS2 and demonstrate theoretically that at high pressure it is prone to the condensation of genuine excitons of finite momentum, whereas the phonon dispersion remains regular. Starting from first-principles many-body perturbation theory, we also predict that the self-consistent electronic charge density of the EI sustains an out-of-plane permanent electric dipole moment with an antiferroelectric texture in the layer plane: At the onset of the EI phase, those optical phonons that share the exciton momentum provide a unique Raman fingerprint for the EI formation. Finally, we identify such fingerprint in a Raman feature that was previously observed experimentally, thus providing direct spectroscopic confirmation of an ideal excitonic insulator phase in bulk MoS2 above 30 GPa. △ Less

Submitted 4 November, 2020; originally announced November 2020.

Journal ref: PNAS vol. 118 no. 13 e2010110118 (2021)

Intrinsic edge excitons in two-dimensional MoS$_2$

Authors: Pino D'Amico, Marco Gibertini, Deborah Prezzi, Daniele Varsano, Andrea Ferretti, Nicola Marzari, Elisa Molinari

Abstract: Using accurate first-principles calculations based on many-body perturbation theory we predict that two-dimensional MoS$_2$ hosts edge excitons with universal character, intrinsic to the existence of edges and lying well below the onset of bulk features. These excitons are largely insensitive to edge terminations or orientation, persisting even in the presence of metallic screening at zigzag edges… ▽ More Using accurate first-principles calculations based on many-body perturbation theory we predict that two-dimensional MoS$_2$ hosts edge excitons with universal character, intrinsic to the existence of edges and lying well below the onset of bulk features. These excitons are largely insensitive to edge terminations or orientation, persisting even in the presence of metallic screening at zigzag edges, with large binding energies of $\sim$0.4 eV. Additional excitons can also emerge in ultranarrow ribbons, or as a function of the chemical nature of the termination. The chemical, structural, and electronic similarities with Se- or W-based transition-metal dichalcogenides suggest that these optical features could be common in this class of materials. △ Less

Submitted 4 September, 2019; originally announced September 2019.

Comments: 7 pages main text, 5 pages supplemental material, 12 figures

Journal ref: Phys. Rev. B 101, 161410 (2020)

A monolayer transition metal dichalcogenide as a topological excitonic insulator

Authors: Daniele Varsano, Maurizia Palummo, Elisa Molinari, Massimo Rontani

Abstract: Monolayer transition metal dichalcogenides in the T' phase promise to realize the quantum spin Hall (QSH) effect at room temperature, because they exhibit a prominent spin-orbit gap between inverted bands in the bulk. Here we show that the binding energy of electron-hole pairs excited through this gap is larger than the gap itself in MoS2, a paradigmatic material that we investigate from first pri… ▽ More Monolayer transition metal dichalcogenides in the T' phase promise to realize the quantum spin Hall (QSH) effect at room temperature, because they exhibit a prominent spin-orbit gap between inverted bands in the bulk. Here we show that the binding energy of electron-hole pairs excited through this gap is larger than the gap itself in MoS2, a paradigmatic material that we investigate from first principles by many-body perturbation theory (MBPT). This paradoxical result hints at the instability of the T' phase against the spontaneous generation of excitons, and indeed we find that it gives rise to a recostructed `excitonic insulator' ground state. Importantly, we show that in this system topological and excitonic order cooperatively enhance the bulk gap by breaking the crystal inversion symmetry, in contrast to the case of bilayers where the frustration between the two orders is relieved by breaking time reversal symmetry. The excitonic topological insulator departs distinctively from the bare topological phase as it lifts the band spin degeneracy, which results in circular dichroism. A moderate biaxial strain applied to the system leads to two additional excitonic phases, different in their topological character but both ferroelectric as an effect of electron-electron interactions. △ Less

Submitted 19 June, 2019; originally announced June 2019.

Comments: five colour figures

Journal ref: This is a pre-print of an article published in Nature Nanotechnology (2020)

Quantifying the Plasmonic Character of Optical Excitations in a Molecular J-Aggregate

Authors: Michele Guerrini, Arrigo Calzolari, Daniele Varsano, Stefano Corni

Abstract: The definition of plasmon at the microscopic scale is far from being understood. Yet, it is very important to recognize plasmonic features in optical excitations, as they can inspire new applications and trigger new discoveries by analogy with the rich phenomenology of metal nanoparticle plasmons. Recently, the concepts of plasmonicity index and the generalized plasmonicity index (GPI) have been d… ▽ More The definition of plasmon at the microscopic scale is far from being understood. Yet, it is very important to recognize plasmonic features in optical excitations, as they can inspire new applications and trigger new discoveries by analogy with the rich phenomenology of metal nanoparticle plasmons. Recently, the concepts of plasmonicity index and the generalized plasmonicity index (GPI) have been devised as computational tools to quantify the plasmonic nature of optical excitations. The question may arise whether any strong absorption band, possibly with some sort of collective character in its microscopic origin, shares the status of plasmon. Here we demonstrate that this is not always the case, by considering a well-known class of systems represented by J-aggregates molecular crystals, characterized by the intense J band of absorption. By means of first-principles simulations, based on a many-body perturbation theory formalism, we investigate the optical properties of a J-aggregate made of push-pull organic dyes. We show that the effect of aggregation is to lower the GPI associated with the J-band with respect to the isolated dye one, which corresponds to a nonplasmonic character of the electronic excitations. In order to rationalize our finding, we then propose a simplified one-dimensional theoretical model of the J-aggregate. A useful microscopic picture of what discriminates a collective molecular crystal excitation from a plasmon is eventually obtained. △ Less

Submitted 16 May, 2019; originally announced May 2019.

Comments: Published by ACS under ACS AuthorChoice license

Journal ref: J. Chem. Theory Comput., 2019, 15 (5), pp 3197-3203

Interplay between Intra- and Intermolecular Charge Transfer in the Optical Excitations of J-Aggregates

Authors: Michele Guerrini, Caterina Cocchi, Arrigo Calzolari, Daniele Varsano, Stefano Corni

Abstract: In a first-principles study based on density functional theory and many-body perturbation theory, we address the interplay between intra- and intermolecular interactions in a J-aggregate formed by push-pull organic dyes by investigating its electronic and optical properties. We find that the most intense excitation dominating the spectral onset of the aggregate, i.e., the J-band, exhibits a combin… ▽ More In a first-principles study based on density functional theory and many-body perturbation theory, we address the interplay between intra- and intermolecular interactions in a J-aggregate formed by push-pull organic dyes by investigating its electronic and optical properties. We find that the most intense excitation dominating the spectral onset of the aggregate, i.e., the J-band, exhibits a combination of intramolecular charge transfer, coming from the push-pull character of the constituting dyes, and intermolecular charge transfer, due to the dense molecular packing. We also show the presence of a pure intermolecular charge-transfer excitation within the J-band, which is expected to play a relevant role in the emission properties of the J-aggregate. Our results shed light on the microscopic character of optical excitations of J-aggregates and offer new perspectives to further understand the nature of collective excitations in organic semiconductors. △ Less

Submitted 24 April, 2019; originally announced April 2019.

Comments: published under ACS Authorchoice license

Journal ref: J. Phys. Chem. C, 2019, 123 (11), pp 6831-6838

Reproducibility in $G_0W_0$ Calculations for Solids

Authors: Tonatiuh Rangel, Mauro Del Ben, Daniele Varsano, Gabriel Antonius, Fabien Bruneval, Felipe H. da Jornada, Michiel J. van Setten, Okan K. Orhan, David D. O'Regan, Andrew Canning, Andrea Ferretti, Andrea Marini, Gian-Marco Rignanese, Jack Deslippe, Steven G. Louie, Jeffrey B. Neaton

Abstract: Ab initio many-body perturbation theory within the $GW$ approximation is a Green's function formalism widely used in the calculation of quasiparticle excitation energies of solids. In what has become an increasingly standard approach, Kohn-Sham eigenenergies, generated from a DFT calculation with a strategically-chosen exchange correlation functional ``starting point'', are used to construct $G$ a… ▽ More Ab initio many-body perturbation theory within the $GW$ approximation is a Green's function formalism widely used in the calculation of quasiparticle excitation energies of solids. In what has become an increasingly standard approach, Kohn-Sham eigenenergies, generated from a DFT calculation with a strategically-chosen exchange correlation functional ``starting point'', are used to construct $G$ and $W$, and then perturbatively corrected by the resultant $GW$ self-energy. In practice, there are several ways to construct the $GW$ self-energy, and these can lead to variations in predicted quasiparticle energies. For example, for ZnO and TiO$_2$, reported $GW$ fundamental gaps can vary by more than 1 eV. In this work, we address the convergence and key approximations in contemporary $G_0W_0$ calculations, including frequency-integration schemes and the treatment of the Coulomb divergence in the exact-exchange term. We study several systems,and compare three different $GW$ codes: BerkeleyGW, Abinit and Yambo. We demonstrate, for the first time, that the same quasiparticle energies for systems in the condensed phase can be obtained with different codes, and we provide a comprehensive assessment of implementations of the $GW$ approximation. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Journal ref: Comput. Phys. Commun. 255, 107242 (2020)

Many-body perturbation theory calculations using the yambo code

Authors: D. Sangalli, A. Ferretti, H. Miranda, C. Attaccalite, I. Marri, E. Cannuccia, P. Melo, M. Marsili, F. Paleari, A. Marrazzo, G. Prandini, P. Bonfà, M. O. Atambo, F. Affinito, M. Palummo, A. Molina-Sánchez, C. Hogan, M. Grüning, D. Varsano, A. Marini

Abstract: yambo is an open source project aimed at studying excited state properties of condensed matter systems from first principles using many-body methods. As input, yambo requires ground state electronic structure data as computed by density functional theory codes such as quantum-espresso and abinit. yambo's capabilities include the calculation of linear response quantities (both independent-particle… ▽ More yambo is an open source project aimed at studying excited state properties of condensed matter systems from first principles using many-body methods. As input, yambo requires ground state electronic structure data as computed by density functional theory codes such as quantum-espresso and abinit. yambo's capabilities include the calculation of linear response quantities (both independent-particle and including electron-hole interactions), quasi-particle corrections based on the GW formalism, optical absorption, and other spectroscopic quantities. Here we describe recent developments ranging from the inclusion of important but oft-neglected physical effects such as electron-phonon interactions to the implementation of a real-time propagation scheme for simulating linear and non-linear optical properties. Improvements to numerical algorithms and the user interface are outlined. Particular emphasis is given to the new and efficient parallel structure that makes it possible to exploit modern high performance computing architectures. Finally, we demonstrate the possibility to automate workflows by interfacing with the yambopy and AiiDA software tools. △ Less

Submitted 7 June, 2019; v1 submitted 11 February, 2019; originally announced February 2019.

Journal ref: Journal of Physics: Condensed Matter, Volume 31, Number 32 (2019)

Electronic and optical properties of doped TiO2 by many-body perturbation theory

Authors: Michael O. Atambo, Daniele Varsano, Andrea Ferretti, S. Samaneh Ataei, Marilia J. Caldas, Elisa Molinari, Annabella Selloni

Abstract: Doping is one of the most common strategies for improving the photocatalytic and solar energy conversion properties of TiO$_2$, hence an accurate theoretical description of the electronic and optical properties of doped TiO$_2$ is of both scientific and practical interest. In this work we use many-body perturbation theory techniques to investigate two typical n-type dopants, Niobium and Hydrogen,… ▽ More Doping is one of the most common strategies for improving the photocatalytic and solar energy conversion properties of TiO$_2$, hence an accurate theoretical description of the electronic and optical properties of doped TiO$_2$ is of both scientific and practical interest. In this work we use many-body perturbation theory techniques to investigate two typical n-type dopants, Niobium and Hydrogen, in TiO$_2$ rutile. Using the GW approximation to determine band edges and defect energy levels, and the Bethe Salpeter equation for the calculation of the absorption spectra, we find that the defect energy levels form non-dispersive bands %associated with localized states lying $\simeq 2.2 eV$ above the top of the corresponding valence bands ($\simeq 0.9 eV$ below the conduction bands of the {\it pristine} material). The defect states are also responsible for the appearance of low energy absorption peaks that enhance the solar spectrum absorption of rutile. The spatial distributions of the excitonic wavefunctions associated with these low energy excitations are very different for the two dopants, suggesting a larger mobility of photoexcited electrons in Nb-TiO$_2$. △ Less

Submitted 16 October, 2019; v1 submitted 10 November, 2018; originally announced November 2018.

Comments: 8 pages, 6 figures

Journal ref: Phys. Rev. Materials 3, 045401 (2019)

Ferromagnetic and Antiferromagnetic Coupling of Spin Molecular Interfaces with High Thermal Stability

Authors: Giulia Avvisati, Claudia Cardoso, Daniele Varsano, Andrea Ferretti, Pierluigi Gargiani, Maria Grazia Betti

Abstract: We report an advanced organic spin-interface architecture with magnetic remanence at room temperature, constituted by metal phthalocyanine molecules magnetically coupled with Co layer(s), mediated by graphene. Fe- and Cu-phthalocyanines assembled on graphene/Co have identical structural configurations, but FePc couples antiferromagnetically with Co up to room temperature, while CuPc couples ferrom… ▽ More We report an advanced organic spin-interface architecture with magnetic remanence at room temperature, constituted by metal phthalocyanine molecules magnetically coupled with Co layer(s), mediated by graphene. Fe- and Cu-phthalocyanines assembled on graphene/Co have identical structural configurations, but FePc couples antiferromagnetically with Co up to room temperature, while CuPc couples ferromagnetically with weaker coupling and thermal stability, as deduced by element-selective X-ray magnetic circular dichroic signals. The robust antiferromagnetic coupling is stabilized by a superexchange interaction, driven by the out-of-plane molecular orbitals responsible of the magnetic ground state and electronically decoupled from the underlying metal via the graphene layer, as confirmed by ab initio theoretical predictions. These archetypal spin interfaces can be prototypes to demonstrate how antiferromagnetic and/or ferromagnetic coupling can be optimized by selecting the molecular orbital symmetry. △ Less

Submitted 15 October, 2018; originally announced October 2018.

Comments: This document is the Accepted Manuscript version of a Published Work that appeared in final form in Nano. Lett. copyright American Chemical Society after peer review and technical editing by the publisher. To access the final edited and published work see http://dx.doi.org/10.1021/acs.nanolett.7b04836

Journal ref: Nano Lett., 2018, 18 (4), pp 2268

Interaction-driven giant orbital magnetic moments in carbon nanotubes

Authors: Joshua O. Island, Marvin Ostermann, Lee Aspitarte, Ethan D. Minot, Daniele Varsano, Elisa Molinari, Massimo Rontani, Gary A. Steele

Abstract: Carbon nanotubes continue to be model systems for studies of confinement and interactions. This is particularly true in the case of so-called "ultra-clean" carbon nanotube devices offering the study of quantum dots with extremely low disorder. The quality of such systems, however, has increasingly revealed glaring discrepancies between experiment and theory. Here we address the outstanding anomaly… ▽ More Carbon nanotubes continue to be model systems for studies of confinement and interactions. This is particularly true in the case of so-called "ultra-clean" carbon nanotube devices offering the study of quantum dots with extremely low disorder. The quality of such systems, however, has increasingly revealed glaring discrepancies between experiment and theory. Here we address the outstanding anomaly of exceptionally large orbital magnetic moments in carbon nanotube quantum dots. We perform low temperature magneto-transport measurements of the orbital magnetic moment and find it is up to seven times larger than expected from the conventional semiclassical model. Moreover, the magnitude of the magnetic moment monotonically drops with the addition of each electron to the quantum dot directly contradicting the widely accepted shell filling picture of single-particle levels. We carry out quasiparticle calculations, both from first principles and within the effective-mass approximation, and find the giant magnetic moments can only be captured by considering a self-energy correction to the electronic band structure due to electron-electron interactions. △ Less

Submitted 21 September, 2018; originally announced September 2018.

Comments: 15 pages and 10 figures including supplemental materials

Journal ref: Phys. Rev. Lett. 121, 127704 (2018)

Role of Quantum-Confinement in Anatase Nanosheets

Authors: Daniele Varsano, Giacomo Giorgi, Koichi Yamashita, Maurizia Palummo

Abstract: Despite most of the applications of anatase nanostructures rely on photoexcited charge processes, yet profound theoretical understanding of fundamental related properties is lacking. Here, by means of ab initio ground and excited-state calculations, we reveal, in an unambiguous way, the role of quantum confinement effect and of the surface orientation, on the electronic and optical properties of a… ▽ More Despite most of the applications of anatase nanostructures rely on photoexcited charge processes, yet profound theoretical understanding of fundamental related properties is lacking. Here, by means of ab initio ground and excited-state calculations, we reveal, in an unambiguous way, the role of quantum confinement effect and of the surface orientation, on the electronic and optical properties of anatase nanosheets (NSs). The presence of bound excitons extremely localized along the (001) direction, whose existence has been recently proven also in anatase bulk, explains the different optical behavior found for the two orientations (001) and (101) when the NS thickness increases. We suggest also that the almost two-dimensional nature of these excitons can be related to the improved photoconversion efficiency observed when a high percentage of (001) facet is present in anatase nanocrystals. △ Less

Submitted 21 August, 2018; originally announced August 2018.

Comments: This document is the Accepted Manuscript version of a Published Work that appeared in final form in J. Phys. Chem. Lett. copyright American Chemical Society after peer review and technical editing by the publisher. To access the final edited and published work see http://dx.doi.org/10.1021/acs.jpclett.7b01717

Journal ref: J. Phys. Chem. Lett., 2017, 8 3867

Theoretical S1 <- S0 absorption energies of the anionic forms of oxyluciferin by Variational Monte Carlo and Many Body Green's Function Theory

Authors: Emanuele Coccia, Daniele Varsano, Leonardo Guidoni

Abstract: The structures of three negatively charged forms (anionic keto-1 and enol-1, dianonic enol-2) of oxyluciferin (OxyLuc), which are the most probable emitters responsible for the firefly bioluminescence, have been fully relaxed at the variational Monte Carlo (VMC) level. Absorption energies of the S1<-S0 vertical transition have been computed using different levels of theory, such as TDDFT, CC2 and… ▽ More The structures of three negatively charged forms (anionic keto-1 and enol-1, dianonic enol-2) of oxyluciferin (OxyLuc), which are the most probable emitters responsible for the firefly bioluminescence, have been fully relaxed at the variational Monte Carlo (VMC) level. Absorption energies of the S1<-S0 vertical transition have been computed using different levels of theory, such as TDDFT, CC2 and many body Green's function Theory (MBGFT). The use of MBGFT, by means of the Bethe-Salpeter (BS) formalism, on VMC structures provides results in excellent agreement with the value (2.26(8) eV) obtained by action spectroscopy experiments for the keto-1 form (2.32 eV). To unravel the role of the quality of the optimized ground state geometry, BS excitation energies have also been computed on CASSCF geometries, inducing a non negligible blue shift (0.08 and 0.07 eV for keto-1 and enol-1 forms, respectively) with respect to the VMC ones. Structural effects have been analyzed in terms of over- or under-correlation along the conjugated bonds of OxyLuc by using different methods for the ground-state optimization. The relative stability of the S$_1$ state for the keto-1 and enol-1 forms depends on the method chosen for the excited state calculation, thus representing a fundamental caveat for any theoretical study on these systems. Finally, Kohn-Sham HOMO and LUMO orbitals of enol-2 are (nearly) bound only when the dianion is embedded into a solvent (water and toluene in the present work); excited state calculations are therefore meaningful only in presence of a dielectric medium which localizes the electronic density. The combination of VMC for the ground state geometry and BS formalism for the absorption spectra clearly outperforms standard TDDFT and quantum chemistry approaches. △ Less

Submitted 24 July, 2018; originally announced July 2018.

Comments: This document is the Accepted Manuscript version of a Published Work that appeared in final form in J. Chem. Theory Comput. copyright American Chemical Society after peer review and technical editing by the publisher. To access the final edited and published work see http://dx.doi.org/10.1021/acs.jctc.7b00505

Journal ref: J. Chem. Theory Comput. 13, 4357 (2017)

FePc Adsorption on the Moiré Superstructure of Graphene Intercalated with a Co Layer

Authors: Giulia Avvisati, Simone Lisi, Pierluigi Gargiani, Ada Della Pia, Oreste De Luca, Daniela Pacilé, Claudia Cardoso, Daniele Varsano, Deborah Prezzi, Andrea Ferretti, Maria Grazia Betti

Abstract: The moiré superstructure of graphene grown on metals can drive the assembly of molecular architectures, as iron-phthalocyanine (FePc) molecules, allowing for the production of artificial molecular configurations. A detailed analysis of the Gr/Co interaction upon intercalation (including a modelling of the resulting moiré pattern) is performed here by density functional theory, which provides an ac… ▽ More The moiré superstructure of graphene grown on metals can drive the assembly of molecular architectures, as iron-phthalocyanine (FePc) molecules, allowing for the production of artificial molecular configurations. A detailed analysis of the Gr/Co interaction upon intercalation (including a modelling of the resulting moiré pattern) is performed here by density functional theory, which provides an accurate description of the template as a function of the corrugation parameters. The theoretical results are a preliminary step to describe the interaction process of the FePc molecules adsorption on the Gr/Co system. Core level photoemission and absorption spectroscopies have been employed to control the preferential adsorption regions of the FePc on the graphene moiré superstructure and the interaction of the central Fe ion with the underlying Co. Our results show that upon molecular adsorption the distance of C atoms from the Co template mainly drives the strength of the molecules-substrate interaction, thereby allowing for locally different electronic properties within the corrugated interface. △ Less

Submitted 29 May, 2018; originally announced May 2018.

Comments: This document is the Accepted Manuscript version of a Published Work that appeared in final form in J. Phys. Chem. C , copyright \c{opyright} American Chemical Society after peer review and technical editing by the publisher. To access the final edited and published work see http://dx.doi.org/10.1021/acs.jpcc.6b09875

Journal ref: J. Phys. Chem. C, 2017, 121, 1639

Electronic Structure Evolution during the Growth of Graphene Nanoribbons on Au(110)

Authors: Ada Della Pia, Giulia Avvisati, Oualid Ourdjini, Claudia Cardoso, Daniele Varsano, Deborah Prezzi, Andrea Ferretti, Carlo Mariani, Maria Grazia Betti

Abstract: Surface-assisted polymerization of molecular monomers into extended chains can be used as the seed of graphene nanoribbon (GNR) formation, resulting from a subsequent cyclo-dehydrogenation process. By means of valence-band photoemission and ab-initio density-functional theory (DFT) calculations, we investigate the evolution of molecular states from monomer 10,10'-dibromo-9,9'bianthracene (DBBA) pr… ▽ More Surface-assisted polymerization of molecular monomers into extended chains can be used as the seed of graphene nanoribbon (GNR) formation, resulting from a subsequent cyclo-dehydrogenation process. By means of valence-band photoemission and ab-initio density-functional theory (DFT) calculations, we investigate the evolution of molecular states from monomer 10,10'-dibromo-9,9'bianthracene (DBBA) precursors to polyanthryl polymers, and eventually to GNRs, as driven by the Au(110) surface. The molecular orbitals and the energy level alignment at the metal-organic interface are studied in depth for the DBBA precursors deposited at room temperature. On this basis, we can identify a spectral fingerprint of C-Au interaction in both DBBA single-layer and polymerized chains obtained upon heating. Furthermore, DFT calculations help us evidencing that GNRs interact more strongly than DBBA and polyanthryl with the Au(110) substrate, as a result of their flatter conformation. △ Less

Submitted 20 April, 2017; originally announced April 2017.

Comments: This document is the Accepted Manuscript version of a Published Work that appeared in final form in J. Phys. Chem C after peer review and technical editing by the publisher. To access the final edited and published work see: http://pubsdc3.acs.org/articlesonrequest/AOR-IB49uBxIS3X6JdF3ieQ9

Journal ref: J. Phys. Chem. C 120, 7323 (2016)

Bonds, lone pairs, and shells probed by means of on-top dynamical correlations

Authors: Stefano Pittalis, Daniele Varsano, Alain Delgado, Carlo Andrea Rozzi

Abstract: The Electron Localization Function (ELF) by Becke and Edgecombe [J. Chem. Phys. {\bf 92}, 5397 (1990)] is routinely adopted as a descriptor of atomic shells and covalent bonds. Since the ELF and its related quantities find useful exploitation also in the construction of modern density functionals, the interest in complementing the ELF is linked to both the quests of improving electronic structure… ▽ More The Electron Localization Function (ELF) by Becke and Edgecombe [J. Chem. Phys. {\bf 92}, 5397 (1990)] is routinely adopted as a descriptor of atomic shells and covalent bonds. Since the ELF and its related quantities find useful exploitation also in the construction of modern density functionals, the interest in complementing the ELF is linked to both the quests of improving electronic structure descriptors and density functional approximations. The ELF uses information which is available by considering parallel-spin electron pairs in single-reference many-body states. In this work, we complement this construction with information obtained by considering antiparallel-spin pairs whose short-range correlations are modeled by a density functional approximation. As a result, the approach requires only a contained computational effort. Applications to a variety of systems show that, in this way, we gain a spatial description of the bond in H$_2$ (which is not available with the ELF) together with some trends not optimally captured by the ELF in other prototypical situations. △ Less

Submitted 7 June, 2018; v1 submitted 17 April, 2017; originally announced April 2017.

Carbon nanotubes as excitonic insulators

Authors: Daniele Varsano, Sandro Sorella, Davide Sangalli, Matteo Barborini, Stefano Corni, Elisa Molinari, Massimo Rontani

Abstract: Fifty years ago Walter Kohn speculated that a zero-gap semiconductor might be unstable against the spontaneous generation of excitons---electron-hole pairs bound together by Coulomb attraction. The reconstructed ground state would then open a gap breaking the symmetry of the underlying lattice, a genuine consequence of electronic correlations. Here we show that this excitonic insulator is realized… ▽ More Fifty years ago Walter Kohn speculated that a zero-gap semiconductor might be unstable against the spontaneous generation of excitons---electron-hole pairs bound together by Coulomb attraction. The reconstructed ground state would then open a gap breaking the symmetry of the underlying lattice, a genuine consequence of electronic correlations. Here we show that this excitonic insulator is realized in zero-gap carbon nanotubes by performing first-principles calculations through many-body perturbation theory as well as quantum Monte Carlo. The excitonic order modulates the charge between the two carbon sublattices opening an experimentally observable gap, which scales as the inverse of the tube radius and weakly depends on the axial magnetic field. Our findings call into question the Luttinger liquid paradigm for nanotubes and provide tests to experimentally discriminate between excitonic and Mott insulator. △ Less

Submitted 15 November, 2017; v1 submitted 27 March, 2017; originally announced March 2017.

Comments: This document includes the Supplementary Information. We have corrected typos in Eq. (56) and sublattice vector definitions that were spotted by Dr. Giacomo Sesti after the paper was published

Journal ref: Nature Communications 8, 1461 (2017)

Breathing bands due to molecular order in CH3NH3PbI3

Authors: Malgorzata Wierzbowska, Juan Jose Melendez, Daniele Varsano

Abstract: CH3NH3PbI3 perovskite is nowadays amongst the most promising photovoltaic materials for energy conversion. We have studied by ab-initio calculations, using several levels of approximation - namely density functional theory including spin-orbit coupling and quasi-particle corrections by means of the GW method, as well as pseudopotential self-interaction corrections -, the role of the methylammonium… ▽ More CH3NH3PbI3 perovskite is nowadays amongst the most promising photovoltaic materials for energy conversion. We have studied by ab-initio calculations, using several levels of approximation - namely density functional theory including spin-orbit coupling and quasi-particle corrections by means of the GW method, as well as pseudopotential self-interaction corrections -, the role of the methylammonium orientation on the electronic structure of this perovskite. We have considered many molecular arrangements within 2x2x2 supercells, showing that the relative orientation of the organic molecules is responsible for a huge band gap variation up to 2 eV. The band gap sizes are related to distortions of the PbI3 cage, which are in turn due to electrostatic interactions between this inorganic frame and the molecules. The strong dependence of the band gap on the mutual molecular orientation is confirmed at all levels of approximations. Our results suggest then that the coupling between the molecular motion and the interactions of the molecules with the inorganic cage could help to explain the widening of the absorption spectrum of CH3NH3PbI3 perovskite, consistent with the observed white spectrum. △ Less

Submitted 9 November, 2017; v1 submitted 8 December, 2016; originally announced December 2016.

Journal ref: Computational Materials Science 142 (2018) 361-371

No-go theorem for the description of Mott phenomena with conventional Density Functional Theory methods

Authors: Zu-Jian Ying, Valentina Brosco, Giorgia Maria Lopez, Daniele Varsano, Paola Gori-Giorgi, José Lorenzana

Abstract: Density functional theory provides the most widespread framework for the realistic description of the electronic structure of solids, but the description of strongly-correlated systems has remained so far elusive. Here we consider a particular limit of electrons in a periodic ionic potential in which a one-band description becomes exact all the way from the weakly-correlated metallic regime to the… ▽ More Density functional theory provides the most widespread framework for the realistic description of the electronic structure of solids, but the description of strongly-correlated systems has remained so far elusive. Here we consider a particular limit of electrons in a periodic ionic potential in which a one-band description becomes exact all the way from the weakly-correlated metallic regime to the strongly-correlated Mott-Hubbard regime. We provide a necessary condition a density functional should fulfill to describe Mott-Hubbard behavior and show that it is not satisfied by standard and widely used local, semilocal and hybrid functionals. We illustrate the condition in the case of a few-atom system and provide an analytic approximation to the exact exchange-correlation potential based on a variational wave function which shows explicitly the correct behavior providing a robust scheme to combine lattice and continuum methods. △ Less

Submitted 12 October, 2015; originally announced October 2015.

Comments: 29 pages, 8 figures

Journal ref: Phys. Rev. B 94, 075154 (2016)

Ab initio geometry and bright excitation of carotenoids: Quantum Monte Carlo and Many Body Green's Function Theory calculations on peridinin

Authors: Emanuele Coccia, Daniele Varsano, Leonardo Guidoni

Abstract: In this letter we report the singlet ground state structure of the full carotenoid peridinin by means of variational Monte Carlo (VMC) calculations. The VMC relaxed geometry has an average bond length alternation of 0.1165(10) Å, larger than the values obtained by DFT (PBE, B3LYP and CAM-B3LYP) and shorter than that calculated at the Hartree-Fock (HF) level. TDDFT and EOM-CCSD calculations on a re… ▽ More In this letter we report the singlet ground state structure of the full carotenoid peridinin by means of variational Monte Carlo (VMC) calculations. The VMC relaxed geometry has an average bond length alternation of 0.1165(10) Å, larger than the values obtained by DFT (PBE, B3LYP and CAM-B3LYP) and shorter than that calculated at the Hartree-Fock (HF) level. TDDFT and EOM-CCSD calculations on a reduced peridinin model confirm the HOMO-LUMO major contribution of the Bu+-like (S2) bright excited state. Many Body Green's Function Theory (MBGFT) calculations of the vertical excitation energy of the Bu+-like state for the VMC structure (VMC/MBGFT) provide excitation energy of 2.62 eV, in agreement with experimental results in n-hexane (2.72 eV). The dependence of the excitation energy on the bond length alternation in the MBGFT and TDDFT calculations with different functionals is discussed. △ Less

Submitted 15 April, 2014; originally announced April 2014.

Journal ref: J. Chem. Theory Comput., vol 10, 501-506 (2014)

Ab-initio angle and energy resolved photoelectron spectroscopy with time-dependent density-functional theory

Authors: U. De Giovannini, D. Varsano, M. A. L. Marques, H. Appel, E. K. U. Gross, A. Rubio

Abstract: We present a time-dependent density-functional method able to describe the photoelectron spectrum of atoms and molecules when excited by laser pulses. This computationally feasible scheme is based on a geometrical partitioning that efficiently gives access to photoelectron spectroscopy in time-dependent density-functional calculations. By using a geometrical approach, we provide a simple descripti… ▽ More We present a time-dependent density-functional method able to describe the photoelectron spectrum of atoms and molecules when excited by laser pulses. This computationally feasible scheme is based on a geometrical partitioning that efficiently gives access to photoelectron spectroscopy in time-dependent density-functional calculations. By using a geometrical approach, we provide a simple description of momentum-resolved photoe- mission including multi-photon effects. The approach is validated by comparison with results in the literature and exact calculations. Furthermore, we present numerical photoelectron angular distributions for randomly oriented nitrogen molecules in a short near infrared intense laser pulse and helium-(I) angular spectra for aligned carbon monoxide and benzene. △ Less

Submitted 25 June, 2012; originally announced June 2012.

Comments: Accepted for publication on Phys. Rev. A

Journal ref: Phys. Rev. A 85, 062515 (2012)

Quantum-dot states and optical excitations in edge-modulated graphene nanoribbons

Authors: Deborah Prezzi, Daniele Varsano, Alice Ruini, Elisa Molinari

Abstract: We investigate from first principles the electronic and optical properties of edge-modulated armchair graphene nanoribbons, including both quasi-particle corrections and excitonic effects. Exploiting the oscillating behavior of the ribbon energy gap, we show that minimal width modulations are sufficient to obtain confinement of both electrons and holes, thus forming optically active quantum dots w… ▽ More We investigate from first principles the electronic and optical properties of edge-modulated armchair graphene nanoribbons, including both quasi-particle corrections and excitonic effects. Exploiting the oscillating behavior of the ribbon energy gap, we show that minimal width modulations are sufficient to obtain confinement of both electrons and holes, thus forming optically active quantum dots with unique properties, such as cohexistence of dot-like and extended excitations and fine tunability of optical spectra, with great potential for optoelectronic applications. △ Less

Submitted 18 April, 2011; originally announced April 2011.

Photo-excitation of a light-harvesting supra-molecular triad: a Time-Dependent DFT study

Authors: N. Spallanzani, C. A. Rozzi, D. Varsano, T. Baruah, M. R. Pederson, F. Manghi, A. Rubio

Abstract: We present the first time-dependent density-functional theory (TDDFT) calculation on a light harvesting triad carotenoid-diaryl-porphyrin-C60. Besides the numerical challenge that the ab initio study of the electronic structure of such a large system presents, we show that TDDFT is able to provide an accurate description of the excited state properties of the system. In particular we calculate t… ▽ More We present the first time-dependent density-functional theory (TDDFT) calculation on a light harvesting triad carotenoid-diaryl-porphyrin-C60. Besides the numerical challenge that the ab initio study of the electronic structure of such a large system presents, we show that TDDFT is able to provide an accurate description of the excited state properties of the system. In particular we calculate the photo-absorption spectrum of the supra-molecular assembly, and we provide an interpretation of the photo-excitation mechanism in terms of the properties of the component moieties. The spectrum is in good agreement with experimental data, and provides useful insight on the photo-induced charge transfer mechanism which characterizes the system. △ Less

Submitted 12 March, 2009; v1 submitted 20 January, 2009; originally announced January 2009.

Comments: Accepted for publication on JPC, March 09th 2009

Yambo: an \textit{ab initio} tool for excited state calculations

Authors: Andrea Marini, Conor Hogan, Myrta Grüning, Daniele Varsano

Abstract: {\tt yambo} is an {\it ab initio} code for calculating quasiparticle energies and optical properties of electronic systems within the framework of many-body perturbation theory and time-dependent density functional theory. Quasiparticle energies are calculated within the $GW$ approximation for the self-energy. Optical properties are evaluated either by solving the Bethe--Salpeter equation or by… ▽ More {\tt yambo} is an {\it ab initio} code for calculating quasiparticle energies and optical properties of electronic systems within the framework of many-body perturbation theory and time-dependent density functional theory. Quasiparticle energies are calculated within the $GW$ approximation for the self-energy. Optical properties are evaluated either by solving the Bethe--Salpeter equation or by using the adiabatic local density approximation. {\tt yambo} is a plane-wave code that, although particularly suited for calculations of periodic bulk systems, has been applied to a large variety of physical systems. {\tt yambo} relies on efficient numerical techniques devised to treat systems with reduced dimensionality, or with a large number of degrees of freedom. The code has a user-friendly command-line based interface, flexible I/O procedures and is interfaced to several publicly available density functional ground-state codes. △ Less

Submitted 3 February, 2009; v1 submitted 17 October, 2008; originally announced October 2008.

Comments: This paper describes the features of the Yambo code, whose source is available under the GPL license at www.yambo-code.org

Journal ref: Comp. Phys. Comm. 180, 1392 (2009)

Optical saturation driven by exciton confinement in molecular-chains: a TDDFT study

Authors: Daniele Varsano, Andrea Marini, Angel Rubio

Abstract: We have identified excitonic confinement in one-dimensional molecular chains (i.e. polyacetylene and H$_2$) as the main driving force for the saturation of the chain polarizability as a function of the number of molecular units. This conclusion is based on first principles time--dependent density functional theory calculations performed with a new derived exchange--correlation kernel. The failur… ▽ More We have identified excitonic confinement in one-dimensional molecular chains (i.e. polyacetylene and H$_2$) as the main driving force for the saturation of the chain polarizability as a function of the number of molecular units. This conclusion is based on first principles time--dependent density functional theory calculations performed with a new derived exchange--correlation kernel. The failure of simple local and semi--local functionals is shown to be related to the lack of memory effects, spatial ultranonlocality, and self--interaction corrections. These effects get smaller as the gap of the system reduces, in which case such simple approximations do perform better. △ Less

Submitted 10 October, 2007; originally announced October 2007.

Optical properties of graphene nanoribbons: the role of many-body effects

Authors: D. Prezzi, D. Varsano, A. Ruini, A. Marini, E. Molinari

Abstract: We investigate from first principles the optoelectronic properties of nanometer-sized armchair graphene nanoribbons (GNRs). We show that many-body effects are essential to correctly describe both energy gaps and optical response. As a signature of the confined geometry, we observe strongly bound excitons dominating the optical spectra, with a clear family dependent binding energy. Our results de… ▽ More We investigate from first principles the optoelectronic properties of nanometer-sized armchair graphene nanoribbons (GNRs). We show that many-body effects are essential to correctly describe both energy gaps and optical response. As a signature of the confined geometry, we observe strongly bound excitons dominating the optical spectra, with a clear family dependent binding energy. Our results demonstrate that GNRs constitute 1D nanostructures whose absorption and luminescence performance can be controlled by changing both family and edge termination. △ Less

Submitted 6 June, 2007; originally announced June 2007.

Comments: 5 pages, 3 figures

A TDDFT study of the excited states of DNA bases and their assemblies

Authors: D. Varsano, R. Di Felice, M. A. L. Marques, A. Rubio

Abstract: We present a detailed study of the optical absorption spectra of DNA bases and base pairs, carried out by means of time dependent density functional theory. The spectra for the isolated bases are compared to available theoretical and experimental data and used to assess the accuracy of the method and the quality of the exchange-correlation functional: Our approach turns out to be a reliable tool… ▽ More We present a detailed study of the optical absorption spectra of DNA bases and base pairs, carried out by means of time dependent density functional theory. The spectra for the isolated bases are compared to available theoretical and experimental data and used to assess the accuracy of the method and the quality of the exchange-correlation functional: Our approach turns out to be a reliable tool to describe the response of the nucleobases. Furthermore, we analyze in detail the impact of hydrogen bonding and $π$-stacking in the calculated spectra for both Watson-Crick base pairs and Watson-Crick stacked assemblies. We show that the reduction of the UV absorption intensity (hypochromicity) for light polarized along the base-pair plane depends strongly on the type of interaction. For light polarized perpendicular to the basal plane, the hypochromicity effect is reduced, but another characteristic is found, namely a blue shift of the optical spectrum of the base-assembly compared to that of the isolated bases. The use of optical tools as fingerprints for the characterization of the structure (and type of interaction) is extensively discussed. △ Less

Submitted 3 March, 2006; originally announced March 2006.

Comments: 31 pages, 8 figures

An exact Coulomb cutoff technique for supercell calculations

Authors: Carlo A. Rozzi, Daniele Varsano, Andrea Marini, Eberhard K. U. Gross, Angel Rubio

Abstract: We present a new reciprocal space analytical method to cutoff the long range interactions in supercell calculations for systems that are infinite and periodic in 1 or 2 dimensions, extending previous works for finite systems. The proposed cutoffs are functions in Fourier space, that are used as a multiplicative factor to screen the bare Coulomb interaction. The functions are analytic everywhere… ▽ More We present a new reciprocal space analytical method to cutoff the long range interactions in supercell calculations for systems that are infinite and periodic in 1 or 2 dimensions, extending previous works for finite systems. The proposed cutoffs are functions in Fourier space, that are used as a multiplicative factor to screen the bare Coulomb interaction. The functions are analytic everywhere but in a sub-domain of the Fourier space that depends on the periodic dimensionality. We show that the divergences that lead to the non-analytical behaviour can be exactly cancelled when both the ionic and the Hartree potential are properly screened. This technique is exact, fast, and very easy to implement in already existing supercell codes. To illustrate the performance of the new scheme, we apply it to the case of the Coulomb interaction in systems with reduced periodicity (as one-dimensional chains and layers). For those test cases we address the impact of the cutoff in different relevant quantities for ground and excited state properties, namely: the convergence of the ground state properties, the static polarisability of the system, the quasiparticle corrections in the GW scheme and in the binding energy of the excitonic states in the Bethe-Salpeter equation. The results are very promising. △ Less

Submitted 2 January, 2006; originally announced January 2006.

Comments: Submitted to Physical Review B on Dec 23rd 2005

Time and energy-resolved two photon-photoemission of the Cu(100) and Cu(111) metal surfaces

Authors: Daniele Varsano, M. A. L. Marques, Angel Rubio

Abstract: We present calculations on energy- and time-resolved two-photon photoemission spectra of images states in Cu(100) and Cu(111) surfaces. The surface is modeled by a 1D effective potential and the states are propagated within a real-space, real-time method. To obtain the energy resolved spectra we employ a geometrical approach based on a subdivision of space into two regions. We treat electronic i… ▽ More We present calculations on energy- and time-resolved two-photon photoemission spectra of images states in Cu(100) and Cu(111) surfaces. The surface is modeled by a 1D effective potential and the states are propagated within a real-space, real-time method. To obtain the energy resolved spectra we employ a geometrical approach based on a subdivision of space into two regions. We treat electronic inelastic effects by taking into account the scattering rates calculated within a GW scheme. To get further insight into the decaying mechanism we have also studied the effect of the variation of the classical Hartree potential during the excitation. This effect turns out to be small. △ Less

Submitted 31 March, 2003; originally announced March 2003.

Comments: 11 pages, 7 figures

Spin-Polarization transition in the two dimensional electron gas

Authors: D. Varsano, S. Moroni, G. Senatore

Abstract: We present a numerical study of magnetic phases of the 2D electron gas near freezing. The calculations are performed by diffusion Monte Carlo in the fixed node approximation. At variance with the 3D case we find no evidence for the stability of a partially polarized phase. With plane wave nodes in the trial function, the polarization transition takes place at Rs=20, whereas the best available es… ▽ More We present a numerical study of magnetic phases of the 2D electron gas near freezing. The calculations are performed by diffusion Monte Carlo in the fixed node approximation. At variance with the 3D case we find no evidence for the stability of a partially polarized phase. With plane wave nodes in the trial function, the polarization transition takes place at Rs=20, whereas the best available estimates locate Wigner crystallization around Rs=35. Using an improved nodal structure, featuring optimized backflow correlations, we confirm the existence of a stability range for the polarized phase, although somewhat shrunk, at densities achievable nowadays in 2 dimensional hole gases in semiconductor heterostructures . The spin susceptibility of the unpolarized phase at the magnetic transition is approximately 30 times the Pauli susceptibility. △ Less

Submitted 26 June, 2000; originally announced June 2000.

Comments: 7 pages, 4 figures