Since the 30s the interatomic potential of the beryllium dimer Be$_2$ has been both an experimental and a theoretical challenge. Calculating the ground-state correlation energy of Be$_2$ along its dissociation path is a difficult problem for theory.
We present ab initio many-body perturbation theory calculations of the Be$_2$ interatomic potential using the GW approximation and the Bethe-Salpeter equation (BSE). The ground-state correlation energy is calculated by the trace formula with checks against the adiabatic-connection fluctuation-dissipation theorem formula. We show that inclusion of GW corrections already improves the energy even at the level of the random-phase approximation. At the level of the BSE on top of the GW approximation, our calculation is in surprising agreement with the most accurate theories and with experiment. It even reproduces an experimentally observed flattening of the interatomic potential due to a delicate correlations balance from a competition between covalent and van der Waals bonding.
We present a novel hybrid quantum/classical (QM/MM) approach to the calculation of charged excitations in molecular solids based on the many-body Greens function $GW$ formalism. Molecules described at the $GW$ level are embedded into the crystalline
environment modeled with an accurate classical polarizable scheme. This allows the calculation of electron addition and removal energies in the bulk and at crystal surfaces where charged excitations are probed in photoelectron experiments. By considering the paradigmatic case of pentacene and perfluoropentacene crystals, we discuss the different contributions from intermolecular interactions to electronic energy levels, distinguishing between polarization, which is accounted for combining quantum and classical polarizabilities, and crystal field effects, that can impact energy levels by up to $pm0.6$ eV. After introducing band dispersion, we achieve quantitative agreement (within 0.2 eV) on the ionization potential and electron affinity measured at pentacene and perfluoropentacene crystal surfaces characterized by standing molecules.
Helium atom is the simplest many-body electronic system provided by nature. The exact solution to the Schrodinger equation is known for helium ground and excited states, and represents a workbench for any many-body methodology. Here, we check the ab
initio many-body GW approximation and Bethe-Salpeter equation (BSE) against the exact solution for helium. Starting from Hartree-Fock, we show that GW and BSE yield impressively accurate results on excitation energies and oscillator strength, systematically improving time-dependent Hartree-Fock. These findings suggest that the accuracy of BSE and GW approximations is not significantly limited by self-interaction and self-screening problems even in this few electron limit. We further discuss our results in comparison to those obtained by time-dependent density-functional theory.
Understanding and controlling heat transport in molecular junctions would provide new routes to design nanoscale coupled electronic and phononic devices. Using first principles full quantum calculations, we tune thermal conductance of a molecular jun
ction by mechanically compressing and extending a short alkane chain connected to graphene leads. We find that the thermal conductance of the compressed junction drops by half in comparison to the extended junction, making it possible to turn on and off the heat current. The low conductance of the off state does not vary by further approaching the leads and stems from the suppression of the transmission of the in--plane transverse and longitudinal channels. Furthermore, we show that misalignment of the leads does not reduce the conductance ratio. These results also contribute to the general understanding of thermal transport in molecular junctions.
We report on magneto-transport measurements in InAs nanowires under large magnetic field (up to 55T), providing a direct spectroscopy of the 1D electronic band structure. Large modulations of the magneto-conductance mediated by an accurate control of
the Fermi energy reveal the Landau fragmentation, carrying the fingerprints of the confined InAs material. Our numerical simulations of the magnetic band structure consistently support the experimental results and reveal key parameters of the electronic confinement.
We study within the many-body Greens function GW and Bethe-Salpeter formalisms the excitation energies of a paradigmatic model dipeptide, focusing on the four lowest-lying local and charge-transfer excitations. Our GW calculations are performed at th
e self-consistent level, updating first the quasiparticle energies, and further the single-particle wavefunctions within the static Coulomb-hole plus screened-exchange approximation to the GW self-energy operator. Important level crossings, as compared to the starting Kohn-Sham LDA spectrum, are identified. Our final Bethe-Salpeter singlet excitation energies are found to agree, within 0.07 eV, with CASPT2 reference data, except for one charge-transfer state where the discrepancy can be as large as 0.5 eV. Our results agree best with LC-BLYP and CAM-B3LYP calculations with enhanced long-range exchange, with a 0.1 eV mean absolute error. This has been achieved employing a parameter-free formalism applicable to metallic or insulating extended or finite systems.
The complex scaling method, which consists in continuing spatial coordinates into the complex plane, is a well-established method that allows to compute resonant eigenfunctions of the time-independent Schroedinger operator. Whenever it is desirable t
o apply the complex scaling to investigate resonances in physical systems defined on numerical discrete grids, the most direct approach relies on the application of a similarity transformation to the original, unscaled Hamiltonian. We show that such an approach can be conveniently implemented in the Daubechies wavelet basis set, featuring a very promising level of generality, high accuracy, and no need for artificial convergence parameters. Complex scaling of three dimensional numerical potentials can be efficiently and accurately performed. By carrying out an illustrative resonant state computation in the case of a one-dimensional model potential, we then show that our wavelet-based approach may disclose new exciting opportunities in the field of computational non-Hermitian quantum mechanics.
We study within the many-body Greens function GW and Bethe-Salpeter approaches the neutral singlet excitations of the zinctetraphenylporphyrin and C70 fullerene donor-acceptor complex. The lowest transition is a charge-transfer excitation between the
donor and the acceptor with an energy in excellent agreement with recent constrained density functional theory calculations. Beyond the lowest charge-transfer state, of which the energy can be determined with simple electrostatic models that we validate, the Bethe-Salpeter approach provides the full excitation spectrum. We evidence the existence of hot electron-hole states which are resonant in energy with the lowest donor intramolecular excitation and show an hybrid intramolecular and charge-transfer character, favouring the transition towards charge separation. These findings support the recently proposed scenario for charge separation at donor-acceptor interfaces through delocalized hot charge-transfer states.
We have developed an efficient scalable kernel method for thermal transport in open systems, with which we have computed the thermal conductance of a junction between bulk silicon and silicon nanowires with diameter up to 10 nm. We have devised scali
ng laws for transmission and reflection spectra, which allow us to predict the thermal resistance of bulk-nanowire interfaces with larger cross sections than those achievable with atomistic simulations. Our results indicate the characteristic size beyond which atomistic systems can be treated accurately by mesoscopic theories.
Understanding the behavior of molecular systems under pressure is a fundamental problem in condensed matter physics. In the case of nitrogen, the determination of the phase diagram and in particular of the melting line, are largely open problems. Two
independent experiments have reported the presence of a maximum in the nitrogen melting curve, below 90 GPa, however the position and the interpretation of the origin of such maximum differ. By means of ab initio molecular dynamics simulations based on density functional theory and thermodynamic integration techniques, we have determined the phase diagram of nitrogen in the range between 20 and 100 GPa. We find a maximum in the melting line, related to a transformation in the liquid, from molecular N_2 to polymeric nitrogen accompanied by an insulator-to-metal transition.
