Daniele Varsano

yambo: An ab initio tool for excited state calculations ☆

yambo is an 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. 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. 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 110 procedures and is interfaced to several publicly available density functional ground-state codes.

Exact Coulomb cutoff technique for supercell calculations

Received 23 December 2005; revised manuscript received 31 March 2006; published 26 May 2006 We present a reciprocal space analytical method to cut off the long range interactions in supercell calculations for systems that are infinite and periodic in one or two dimensions, generalizing previous work to treat 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 except in a subdomain of the Fourier space that depends on the periodic dimensionality. We show that the divergences that lead to the nonanalytical behavior can be exactly canceled 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 scheme, we apply it to the case of the Coulomb interaction in systems with reduced periodicity as one-dimensional chains and layers. For these test cases, we address the impact of the cutoff on different relevant quantities for ground and excited state properties, namely: the convergence of the ground state properties, the static polarizability of the system, the quasiparticle corrections in the GW scheme, and the binding energy of the excitonic states in the Bethe-Salpeter equation. The results are very promising and easy to implement in all available first-principles codes.

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

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 one-dimensional nanostructures whose absorption and luminescence performance can be controlled by changing both family and edge termination.

First Principles Effective Electronic Couplings for Hole Transfer in Natural and Size-Expanded DNA

Hole transfer processes between base pairs in natural DNA and size-expanded DNA (xDNA) are studied and compared, by means of an accurate first principles evaluation of the effective electronic couplings (also known as transfer integrals), in order to assess the effect of the base augmentation on the efficiency of charge transport through double-stranded DNA. According to our results, the size expansion increases the average electronic coupling, and thus the CT rate, with potential implications in molecular biology and in the implementation of molecular nanoelectronics. Our analysis shows that the effect of the nucleobase expansion on the charge-transfer (CT) rate is sensitive to the sequence of base pairs. Furthermore, we find that conformational variability is an important factor for the modulation of the CT rate. From a theoretical point of view, this work offers a contribution to the CT chemistry in pi-stacked arrays. Indeed, we compare our methodology against other standard computational frameworks that have been adopted to tackle the problem of CT in DNA, and unravel basic principles that should be accounted for in selecting an appropriate theoretical level.

Ab Initio Geometry and Bright Excitation of Carotenoids: Quantum Monte Carlo and Many Body Green’s Function Theory Calculations on Peridinin

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 Greens Function Theory (MBGFT) calculations of the vertical excitation energy of the Bu(+)-like state for the VMC structure (VMC/MBGFT) provide an 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.

Dielectric and thermal effects on the optical properties of natural dyes: a case study on solvated cyanin.

The optical properties of the flavylium state of the cyanin dye are simulated numerically by combining Car-Parrinello molecular dynamics and linear-response time-dependent density functional theory calculations. The spectrum of the dye calculated in the gas phase is characterized by two peaks in the yellow and in the blue (green and violet), using a GGA-PBE (hybrid-B3LYP) DFT functional, which would bring about a greenish (bright orange) color incompatible with the dark purple hue observed in nature. Describing the effect of the water solvent through a polarizable continuum model does not modify qualitatively the resulting picture. An explicit simulation of both solvent and thermal effects using ab initio molecular dynamics results instead in a spectrum that is compatible with the observed coloration. This result is analyzed in terms of the spectroscopic effects of the molecular distortions induced by thermal fluctuations.

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

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 photoemission including multiphoton 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.

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

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.

Photoexcitation of a Light-Harvesting Supramolecular Triad: A Time-Dependent DFT Study

We present the first time-dependent density functional theory (TDDFT) calculation on a light-harvesting triad carotenoid-diaryl-porphyrin-C(60). 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 photoabsorption spectrum of the supramolecular assembly, and we provide an interpretation of the photoexcitation 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 photoinduced charge-transfer mechanism which characterizes the system.

Spin-polarization transition in the two-dimensional electron gas

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.