Tonatiuh Rangel

Structural and excited-state properties of oligoacene crystals from first principles

U.S. Department of Energy, DE-AC02-05CH11231, DE-AC02-05CH11231, Vetenskapsradet, United States - Israel Binational Science Foundation

Band structure of gold from many-body perturbation theory

The band structure of gold is calculated using ab initio many-body perturbation theory. Different approximations within the GW approach are considered. Standard single-shot G(0)W(0) corrections modify the sp-like bands while leaving unchanged the 5d occupied bands. Beyond G(0)W(0), quasiparticle self-consistency on the wave functions lowers the 5d bands. Globally, many-body effects achieve an opening of the 5d-6sp interband gap of similar to 0.4 to similar to 0.8 eV, reducing the discrepancy with the experiment. Finally, the quasiparticle band structure is compared to the one obtained by the widely used HSE (Heyd, Scuseria, and Ernzerhof) hybrid functional.

Evaluating the GW Approximation with CCSD(T) for Charged Excitations Across the Oligoacenes

Charged excitations of the oligoacene family of molecules, relevant for astrophysics and technological applications, are widely studied and therefore provide an excellent system for benchmarking theoretical methods. In this work, we evaluate the performance of many-body perturbation theory within the GW approximation relative to new high-quality CCSD(T) reference data for charged excitations of the acenes. We compare GW calculations with a number of hybrid density functional theory starting points and with eigenvalue self-consistency. Special focus is given to elucidating the trend of GW-predicted excitations with molecule length increasing from benzene to hexacene. We find that GW calculations with starting points based on an optimally tuned range-separated hybrid (OTRSH) density functional and eigenvalue self-consistency can yield quantitative ionization potentials for the acenes. However, for larger acenes, the predicted electron affinities can deviate considerably from reference values. Our work paves the way for predictive and cost-effective GW calculations of charged excitations of molecules and identifies certain limitations of current GW methods used in practice for larger molecules.

Transport properties of molecular junctions from many-body perturbation theory

The conductance of single molecule junctions is calculated using a Landauer approach combined with many-body perturbation theory to account for electron correlation. Contrary to intuition, a mere correction of the density-functional theory eigenvalues is found not to affect noticeably the zero-bias conductance. To improve the agreement with experiments, it is necessary to go beyond this standard procedure by also updating the wave functions. This leads to both the reduction of the molecular character and the increase of the e(g)(d(z)2) character on gold atoms around the Fermi energy.

Origin of Magnetism and Quasiparticles Properties in Cr-Doped TiO2

Combining the local spin density approximation (LSDA)+U and an analysis of superexchange interactions beyond density functional theory, we describe the magnetic ground state of Cr-doped TiO2, an intensively studied and debated dilute magnetic oxide. In parallel, we correct our LSDA+U (+ superexchange) ground state through GW corrections (GW@LSDA+U) that reproduce the position of the impurity states and the band gaps in satisfying agreement with experiments. Because of the different topological coordinations of Cr-Cr bonds in the ground states of rutile and anatase, superexchange interactions induce either ferromagnetic or antiferromagnetic couplings of Cr ions. In Cr-doped anatase, this interaction leads to a new mechanism which stabilizes a (nonrobust) ferromagnetic ground state, in keeping with experimental evidence, without the need to invoke F-center exchange. The interplay between structural defects and vacancies in contributing to the superexchange is also unveiled.

Large Bulk Photovoltaic Effect and Spontaneous Polarization of Single-Layer Monochalcogenides

We use a first-principles density functional theory approach to calculate the shift current and linear absorption of uniformly illuminated single-layer Ge and Sn monochalcogenides. We predict strong absorption in the visible spectrum and a large effective three-dimensional shift current (∼100  μA/V^{2}), larger than has been previously observed in other polar systems. Moreover, we show that the integral of the shift-current tensor is correlated to the large spontaneous effective three-dimensional electric polarization (∼1.9  C/m^{2}). Our calculations indicate that the shift current will be largest in the visible spectrum, suggesting that these monochalcogenides may be promising for polar optoelectronic devices. A Rice-Mele tight-binding model is used to rationalize the shift-current response for these systems, and its dependence on polarization, in general terms with implications for other polar materials.

Structural and vibrational stability of M and Z phases of silicon and germanium from first principles

First-principles calculations were performed to investigate the structural feasibility of M and Z phases (novel monoclinic and orthorhombic structures recently reported for carbon) for silicon and germanium. The lattice parameters, bulk modulus, vibrational properties, and elastic constants are calculated using the local density approximation to describe the exchange-correlation energy, while the optical properties are calculated by using Many-Body Perturbation Theory in the G0W0 approximation. Our results indicate that silicon and germanium with the proposed crystal symmetries are elastically and vibrationally stable and are small band-gap semiconductors. We discuss the possible synthesis of such materials.

Assessment of two hybrid van der Waals density functionals for covalent and non-covalent binding of molecules

Two hybrid van der Waals density functionals (vdW-DFs) are developed using 25% Fock exchange with (i) the consistent-exchange vdW-DF-cx functional [K. Berland and P. Hyldgaard, Phys. Rev. B 89, 035412 (2014)] and (ii) with the vdW-DF2 functional [K. Lee et al., Phys. Rev. B 82, 081101 (2010)]. The ability to describe covalent and non-covalent binding properties of molecules is assessed. For properties related to covalent binding, atomization energies (G2-1 set), molecular reaction energies (G2RC set), and ionization energies (G21IP set) are benchmarked against experimental reference values. We find that hybrid-vdW-DF-cx yields results that are rather similar to those of the standard non-empirical hybrid PBE0 [C. Adamo and V. Barone, J. Chem. Phys. 110, 6158 (1999)], with mean average deviations (MADs) of 4.9 and 5.0 kcal/mol for the G2-1 set, respectively. In this comparison, experimental reference values are used, back corrected by wavefunction-based quantum-chemistry calculations of zero-point energies. Hybrid vdW-DF2 follows somewhat different trends, showing on average significantly larger deviations from the reference energies, with a MAD of 14.5 kcal/mol for the G2-1 set. Non-covalent binding properties of molecules are assessed using the S22 benchmark set of non-covalently bonded dimers and the X40 set of dimers of small halogenated molecules, using wavefunction-based quantum chemistry results as references. For the S22 set, hybrid-vdW-DF-cx performs better than standard vdW-DF-cx for the mostly hydrogen-bonded systems, with MAD dropping from 0.6 to 0.3 kcal/mol, but worse for purely dispersion-bonded systems, with MAD increasing from 0.2 to 0.6 kcal/mol. Hybrid-vdW-DF2 offers a slight improvement over standard vdW-DF2. Similar trends are found for the X40 set, with hybrid-vdW-DF-cx performing particularly well for binding energies involving the strongly polar hydrogen halides, but poorly for systems with tiny binding energies. Our study of the X40 set reveals the potential of mixing Fock exchange with vdW-DF, but also highlights shortcomings of the hybrids constructed here. The solid performance of hybrid-vdW-DF-cx for covalent-bonded systems, as well as the strengths and issues uncovered for non-covalently bonded systems, makes this study a good starting point for developing even more accurate hybrid vdW-DFs.

Can molecular projected density of states (PDOS) be systematically used in electronic conductance analysis

Summary Using benzenediamine and benzenedithiol molecular junctions as benchmarks, we investigate the widespread analysis of the quantum transport conductance in terms of the projected density of states (PDOS) onto molecular orbitals (MOs). We first consider two different methods for identifying the relevant MOs: (1) diagonalization of the Hamiltonian of the isolated molecule and (2) diagonalization of a submatrix of the junction Hamiltonian constructed by considering only basis elements localized on the molecule. We find that these two methods can lead to substantially different MOs and hence PDOS. Furthermore, within Method 1, the PDOS can differ depending on the isolated molecule chosen to represent the molecular junction (e.g., with or without dangling bonds); within Method 2, the PDOS depends on the chosen basis set. We show that these differences can be critical when the PDOS is used to provide a physical interpretation of the conductance (especially when its value is small, as it happens typically at zero bias). In this work, we propose a new approach in an attempt to reconcile the two traditional methods. Although some improvements were achieved, the main problems remain unsolved. Our results raise more general questions and doubts on a PDOS-based analysis of the conductance.

Effects of quantum confinement on excited state properties of SrTiO 3 from ab initio many-body perturbation theory

The Ruddlesden-Popper (RP) homologous series Sr