Real-time approach to the optical properties of solids and nanostructures: Time-dependent Bethe-Salpeter equation

Many-body effects are known to play a crucial role in the electronic and optical properties of solids and nanostructures. Nevertheless, the majority of theoretical and numerical approaches able to capture the influence of Coulomb correlations are restricted to the linear response regime. In this work, we introduce an approach based on a real-time solution of the electronic dynamics. The proposed approach reduces to the well-known Bethe-Salpeter equation in the linear limit regime and it makes it possible, at the same time, to investigate correlation effects in nonlinear phenomena. We show the flexibility and numerical stability of the proposed approach by calculating the dielectric constants and the effect of a strong pulse excitation in bulk h-BN.

Electromagnetic beam profile dynamics in collisional plasmas

The compression of a finite extent Gaussian laser pulse in collisional plasma is investigated. An analytical model is employed to describe the spatiotemporal evolution of a laser pulse propagating through the plasma medium. The pulse geometry is modeled via an appropriate ansatz which takes into account both beam radius (in space) and pulse width (in time). Compression and self-focusing are taken into account via appropriated group velocity dispersion and nonlinearity terms. The competition among the collisional nonlinearity in the plasma and the effect of divergence due to diffraction is pointed out and investigated numerically. Our results suggest that laser pulse compression and intensity localization is enhanced by plasma collisionality. In specific, a pulse width compression by an order of magnitude approximately is observed, for typical collisional laser plasma parameters, along with a significant increase in the intensity.

Local-field and excitonic effects in the optical response of alpha-alumina

This paper introduces key ingredients of the dielectric response of a-alumina that go beyond an independent-particle (IP) treatment of the valence-electron excitations. The optical-response functions were calculated from first-principles both at the Bethe-Salpeter and the random-phase approximation (RPA) levels. Excitonic effects obtained within the Bethe-Salpeter framework were found essential for reproducing the low-energy part of the experimental spectra (below 15 eV) and the bound exciton in particular. For higher energies, local-field effects introduced through the RPA modified considerably the IP results and provided a satisfactory account of the reflectivity spectra and of the position and shape of the dominant bulk plasmon resonance in the electron energy-loss spectra.

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.

Implementation and testing of Lanczos-based algorithms for Random-Phase Approximation eigenproblems

The treatment of the Random-Phase Approximation Hamiltonians, encountered in different frameworks, like time-dependent density functional theory or Bethe-Salpeter equation, is complicated by their non-Hermicity. Compared to their Hermitian Hamiltonian counterparts, computational methods for the treatment of non-Hermitian Hamiltonians are often less efficient and less stable, sometimes leading to the breakdown of the method. Recently [Gruning et al. Nano Lett. 8 (2009) 28201, we have identified that such Hamiltonians are usually pseudo-Hermitian. Exploiting this property, we have implemented an algorithm of the Lanczos type for Random-Phase Approximation Hamiltonians that benefits from the same stability and computational load as its Hermitian counterpart, and applied it to the study of the optical response of carbon nanotubes. We present here the related theoretical grounds and technical details, and study the performance of the algorithm for the calculation of the optical absorption of a molecule within the Bethe-Salpeter equation framework. (C) 2011 Elsevier B.V. All rights reserved.

Exchange-correlation energy and potential as approximate functionals of occupied and virtual Kohn-Sham orbitals: Application to dissociating H-2

The standard local density approximation and generalized gradient approximations fail to properly describe the dissociation of an electron pair bond, yielding large errors (on the order of 50 kcal/mol) at long bond distances. To remedy this failure, a self-consistent Kohn-Sham (KS) method is proposed with the exchange-correlation (xc) energy and potential depending on both occupied and virtual KS orbitals. The xc energy functional of Buijse and Baerends [Mol. Phys. 100, 401 (2002); Phys. Rev. Lett. 87, 133004 (2001)] is employed, which, based on an ansatz for the xc-hole amplitude, is able to reproduce the important dynamical and nondynamical effects of Coulomb correlation through the efficient use of virtual orbitals. Self-consistent calculations require the corresponding xc potential to be obtained, to which end the optimized effective potential (OEP) method is used within the common energy denominator approximation for the static orbital Green's function. The problem of the asymptotic divergence of the xc potential of the OEP when a finite number of virtual orbitals is used is addressed. The self-consistent calculations reproduce very well the entire H-2 potential curve, describing correctly the gradual buildup of strong left-right correlation in stretched H-2. (C) 2003 American Institute of Physics.

Absolute cross sections for dissociative electron attachment to HCCCN

New absolute cross sections for dissociative electron attachment to HCCCN (cyanoacetylene or propiolonitrile) in the range of 0-10 eV electron energy are presented here, which have been determined from a new analysis of previously reported data (Graupner et al 2006 New J. Phys. 8 117). The highest cross sections are observed for the formation of CN^- at 5.3 eV and CCCN^- at 5.1 eV; approximately 0.06 Ų and 0.05 Ų respectively. As part of the re-analysis, it was necessary to determine absolute cross sections for electron-impact ionization of HCCCN with the binary-encounter Bethe method. These electron-impact ionization absolute cross sections for HCCCN are also presented here; the maximum value was found to be ∼6.6 Ų at ∼80 eV.