Search

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., Comment: Main: 6 pages and 4 figures, Supplementary: 6 pages and 7 figures

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., Comment: A. Guandalini and D. A. Leon contributed equally to this work

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., Comment: 13 pages, 8 figures

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 $\pi$ plasmon being significantly affected by excitonic effects.

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., Comment: 14 pages, 6 figures

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., Comment: 21 pages

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\'e 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)., Comment: 15 pages, 15 figures

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.

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 $\Gamma$, 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., Comment: 8 pages; 7 figures

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.

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.

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., Comment: Initial submission. Includes supplementary information. Revised version to appear in Nature Physics

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(\omega)$ 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., Comment: 7 figures, 18 pages + 2 supplementary pages

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 {\pi}-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 {\pi} 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., Comment: 5 pages, 4 figures

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 {\pi} 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.

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.

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., Comment: 7 pages main text, 5 pages supplemental material, 12 figures

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., Comment: five colour figures

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., Comment: Published by ACS under ACS AuthorChoice license

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., Comment: published under ACS Authorchoice license

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.

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$., Comment: 8 pages, 6 figures

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., Comment: 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

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., Comment: 15 pages and 10 figures including supplemental materials

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., Comment: 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

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., Comment: 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

The moir\'e 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\'e 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\'e 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., Comment: 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

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., Comment: 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

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.

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., Comment: 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

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.

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., Comment: 29 pages, 8 figures

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) {\AA}, 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.

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.

{\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., Comment: This paper describes the features of the Yambo code, whose source is available under the GPL license at www.yambo-code.org

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.

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., Comment: Submitted to Physical Review B on Dec 23rd 2005

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., Comment: 11 pages, 7 figures