ABINIT: Overview and focus on selected capabilities.

abinit is probably the first electronic-structure package to have been released under an open-source license about 20 years ago. It implements density functional theory, density-functional perturbation theory (DFPT), many-body perturbation theory (GW approximation and Bethe-Salpeter equation), and more specific or advanced formalisms, such as dynamical mean-field theory (DMFT) and the "temperature-dependent effective potential" approach for anharmonic effects. Relying on planewaves for the representation of wavefunctions, density, and other space-dependent quantities, with pseudopotentials or projector-augmented waves (PAWs), it is well suited for the study of periodic materials, although nanostructures and molecules can be treated with the supercell technique. The present article starts with a brief description of the project, a summary of the theories upon which abinit relies, and a list of the associated capabilities. It then focuses on selected capabilities that might not be present in the majority of electronic structure packages either among planewave codes or, in general, treatment of strongly correlated materials using DMFT; materials under finite electric fields; properties at nuclei (electric field gradient, Mössbauer shifts, and orbital magnetization); positron annihilation; Raman intensities and electro-optic effect; and DFPT calculations of response to strain perturbation (elastic constants and piezoelectricity), spatial dispersion (flexoelectricity), electronic mobility, temperature dependence of the gap, and spin-magnetic-field perturbation. The abinit DFPT implementation is very general, including systems with van der Waals interaction or with noncollinear magnetism. Community projects are also described: generation of pseudopotential and PAW datasets, high-throughput calculations (databases of phonon band structure, second-harmonic generation, and GW computations of bandgaps), and the library libpaw. abinit has strong links with many other software projects that are briefly mentioned.

Correction: Enhancement of CO2 binding and mechanical properties upon diamine functionalization of M2(dobpdc) metal-organic frameworks.

[This corrects the article DOI: 10.1039/C7SC05217K.].

Resonance Raman Characterization of Tetracene Monomer and Nanocrystals: Excited State Lattice Distortions With Implications For Efficient Singlet Fission.

The characterization of specific phonon modes and exciton states that lead to efficient singlet fission (SF) may be instrumental in the design of the next generation of high-efficiency photovoltaic devices. To this end, we analyze the absolute resonance Raman (RR) cross sections for tetracene (Tc) both as a monomer in solution and as a crystalline solid in an aqueous suspension of nanocrystals. For both systems, a time-dependent wavepacket model is developed that is consistent with the absolute RR cross sections, the magnitude of the absorption cross sections, and the vibronic line shapes of the fluorescence. In the monomer, the intramolecular reorganization energy is between 1500 and 1800 cm-1 and the solvent reorganization energy is 70 cm-1. In nanocrystals, the total reorganization is diminished to less than 600 cm-1. The lowest energy exciton has an estimated intramolecular reorganization energy between 300 and 500 cm-1 while intermolecular librational phonons have a reorganization energy of about 130 cm-1. The diminished reorganization energy of the nanocrystal is interpreted in the context of the delocalization of the band-edge exciton onto about ∼7 molecules. When electron and electron-hole correlations are included within many-body perturbation theory, the polarized absorption spectra of crystalline Tc are calculated and found to be in agreement with experiment. The low-lying exciton states and optically active phonons that contribute to the polarized crystal absorption are identified. The likely role of coherent exciton phonon evolution in the SF process is discussed.

Enhancement of CO2 binding and mechanical properties upon diamine functionalization of M2(dobpdc) metal-organic frameworks.

The family of diamine-appended metal-organic frameworks exemplified by compounds of the type mmen-M2(dobpdc) (mmen = N,N'-dimethylethylenediamine; M = Mg, Mn, Fe, Co, Zn; dobpdc4- = 4,4'-dioxidobiphenyl-3,3'-dicarboxylate) are adsorbents with significant potential for carbon capture, due to their high working capacities and strong selectivity for CO2 that stem from a cooperative adsorption mechanism. Herein, we use first-principles density functional theory (DFT) calculations to quantitatively investigate the role of mmen ligands in dictating the framework properties. Our van der Waals-corrected DFT calculations indicate that electrostatic interactions between ammonium carbamate units significantly enhance the CO2 binding strength relative to the unfunctionalized frameworks. Additionally, our computed energetics show that mmen-M2(dobpdc) materials can selectively adsorb CO2 under humid conditions, in agreement with experimental observations. The calculations further predict an increase of 112% and 124% in the orientationally-averaged Young's modulus E and shear modulus G, respectively, for mmen-Zn2(dobpdc) compared to Zn2(dobpdc), revealing a dramatic enhancement of mechanical properties associated with diamine functionalization. Taken together, our calculations demonstrate how functionalization with mmen ligands can enhance framework gas adsorption and mechanical properties.

Low-lying excited states in crystalline perylene.

Organic materials are promising candidates for advanced optoelectronics and are used in light-emitting diodes and photovoltaics. However, the underlying mechanisms allowing the formation of excited states responsible for device functionality, such as exciton generation and charge separation, are insufficiently understood. This is partly due to the wide range of existing crystalline polymorphs depending on sample preparation conditions. Here, we determine the linear optical response of thin-film single-crystal perylene samples of distinct polymorphs in transmission and reflection geometries. The sample quality allows for unprecedented high-resolution spectroscopy, which offers an ideal opportunity for judicious comparison between theory and experiment. Excellent agreement with first-principles calculations for the absorption based on the GW plus Bethe-Salpeter equation (GW-BSE) approach of many-body perturbation theory (MBPT) is obtained, from which a clear picture of the low-lying excitations in perylene emerges, including evidence of an exciton-polariton stopband, as well as an assessment of the commonly used Tamm-Dancoff approximation to the GW-BSE approach. Our findings on this well-controlled system can guide understanding and development of advanced molecular solids and functionalization for applications.

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.

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.

An assessment of low-lying excitation energies and triplet instabilities of organic molecules with an ab initio Bethe-Salpeter equation approach and the Tamm-Dancoff approximation.

The accurate prediction of singlet and triplet excitation energies is an area of intense research of significant fundamental interest and critical for many applications. Most calculations of singlet and triplet energies use time-dependent density functional theory (TDDFT) in conjunction with an approximate exchange-correlation functional. In this work, we examine and critically assess an alternative method for predicting low-lying neutral excitations with similar computational cost, the ab initio Bethe-Salpeter equation (BSE) approach, and compare results against high-accuracy wavefunction-based methods. We consider singlet and triplet excitations of 27 prototypical organic molecules, including members of Thiel's set, the acene series, and several aromatic hydrocarbons exhibiting charge-transfer-like excitations. Analogous to its impact in TDDFT, we find that the Tamm-Dancoff approximation (TDA) overcomes triplet instabilities in the BSE approach, improving both triplet and singlet energetics relative to higher level theories. Finally, we find that BSE-TDA calculations built on effective DFT starting points, such as those utilizing optimally tuned range-separated hybrid functionals, can yield accurate singlet and triplet excitation energies for gas-phase organic molecules.

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.

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

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.