First-principles theory of the nitrogen interstitial in hBN: a plausible model for the blue emitter.

Color centers in hexagonal boron nitride (hBN) have attracted considerable attention due to their remarkable optical properties enabling robust room temperature photonics and quantum optics applications in the visible spectral range. On the other hand, identification of the microscopic origin of color centers in hBN has turned out to be a great challenge that hinders the in-depth theoretical characterization, on-demand fabrication, and development of integrated photonic devices. This is also true for the blue emitter, which is a result of irradiation damage in hBN, emitting at 436 nm wavelength with desirable properties. Here, we propose the negatively charged nitrogen split interstitial defect in hBN as a plausible microscopic model for the blue emitter. To this end, we carried out a comprehensive first-principles theoretical study of the nitrogen interstitial. We carefully analyzed the accuracy of first-principles methods and showed that the commonly used HSE hybrid exchange-correlation functional fails to describe the electronic structure of this defect. Using the generalized Koopman's theorem, we fine-tuned the functional and obtained a zero-phonon photoluminescence (ZPL) energy in the blue spectral range. We showed that the defect exhibits a high emission rate in the ZPL line and features a characteristic phonon side band that resembles the blue emitter's spectrum. Furthermore, we studied the electric field dependence of the ZPL and numerically showed that the defect exhibits a quadratic Stark shift that is perpendicular to plane electric fields, making the emitter insensitive to electric field fluctuations in the first order. Our work emphasizes the need for assessing the accuracy of common first-principles methods in hBN and exemplifies a workaround methodology. Furthermore, our work is a step towards understanding the structure of the blue emitter and utilizing it in photonics applications.

Predicting the FCI Energy of Large Systems to Chemical Accuracy from Restricted Active Space Density Matrix Renormalization Group Calculations.

We theoretically derive and validate with large scale simulations a remarkably accurate power law scaling of errors for the restricted active space density matrix renormalization group (DMRG-RAS) method [J. Phys. Chem. A 126, 9709] in electronic structure calculations. This yields a new extrapolation method, DMRG-RAS-X, which reaches chemical accuracy for strongly correlated systems such as the chromium dimer, dicarbon up to a large cc-pVQZ basis and even a large chemical complex such as the FeMoco with significantly lower computational demands than those of previous methods. The method is free of empirical parameters, performed robustly and reliably in all examples we tested, and has the potential to become a vital alternative method for electronic structure calculations in quantum chemistry and more generally for the computation of strong correlations in nuclear and condensed matter physics.

Toward Large-Scale Restricted Active Space Calculations Inspired by the Schmidt Decomposition.

We present an alternative, memory-efficient, Schmidt decomposition-based description of the inherently bipartite restricted active space (RAS) scheme, which can be implemented effortlessly within the density matrix renormalization group (DMRG) method via the dynamically extended active space procedure. Benchmark calculations are compared against state-of-the-art results of C2 and Cr2, which are notorious for their multireference character. Our results for ground and excited states together with spectroscopic constants demonstrate that the proposed novel approach, dubbed as DMRG-RAS, which is variational and free of uncontrolled method errors, has the potential to outperfom conventional methods for strongly correlated molecules.

Sensitivity of coupled cluster electronic properties on the reference determinant: Can Kohn-Sham orbitals be more beneficial than Hartree-Fock orbitals?

Coupled cluster calculations are traditionally performed over Hartree-Fock reference orbitals (HF-CC methodology). However, in the literature it has been repeatedly raised whether the use of a Kohn-Sham reference (KS-CC methodology) might result in improved performance relative to HF-CC. In the present study, we re-examine the relation of HF-CC and KS-CC methods by comparing the results of widely applied truncated CC calculations (CCSD, CCSD(T), CCSDT) to the limit of full configuration interaction (FCI), which serves as an undebatable reference point of accuracy. Based on a series of CC calculations on diatoms and transition metal complexes, we demonstrate that no systematic improvement of coupled cluster electronic energies, densities and chemical reaction energies is expected when changing from HF to a KS reference. Nevertheless, fortuitous error cancellations might occasionally result in illusory improvement compared to HF-CC. Altogether, the application of KS-CC is not advantageous over HF-CC, but it is also not unreasonable as the choice of reference has negligible influence on the results at sufficiently high CC levels. In addition, KS-CC can be a particularly useful alternative if difficulties are encountered in HF or HF-CC convergence. It is also notable that KS-CC results are found to be practically independent of the chosen density functional, which implies that almost any KS-CC method can be used in place of HF-CC.

DMRG on Top of Plane-Wave Kohn-Sham Orbitals: A Case Study of Defected Boron Nitride.

In this paper, we analyze the numerical aspects of the inherent multireference density matrix renormalization group (DMRG) calculations on top of the periodic Kohn-Sham density functional theory using the complete active space approach. The potential of the framework is illustrated by studying hexagonal boron nitride nanoflakes embedding a charged single boron vacancy point defect by revealing a vertical energy spectrum with a prominent multireference character. We investigate the consistency of the DMRG energy spectrum from the perspective of sample size, basis size, and active space selection protocol. Results obtained from standard quantum chemical atom-centered basis calculations and plane-wave based counterparts show excellent agreement. Furthermore, we also discuss the spectrum of the periodic sheet which is in good agreement with extrapolated data of finite clusters. These results pave the way toward applying the DMRG method in extended correlated solid-state systems, such as point defect qubit in wide band gap semiconductors.

Quantum information-based analysis of electron-deficient bonds.

Recently, the correlation theory of the chemical bond was developed, which applies concepts of quantum information theory for the characterization of chemical bonds, based on the multiorbital correlations within the molecule. Here, for the first time, we extend the use of this mathematical toolbox for the description of electron-deficient bonds. We start by verifying the theory on the textbook example of a molecule with three-center two-electron bonds, namely, diborane(6). We then show that the correlation theory of the chemical bond is able to properly describe the bonding situation in more exotic molecules which have been synthesized and characterized only recently, in particular, the diborane molecule with four hydrogen atoms [diborane(4)] and a neutral zerovalent s-block beryllium complex, whose surprising stability was attributed to a strong three-center two-electron π bond stretching across the C-Be-C core. Our approach is of high importance especially in the light of a constant chase after novel compounds with extraordinary properties where the bonding is expected to be unusual.

The correlation theory of the chemical bond.

The quantum mechanical description of the chemical bond is generally given in terms of delocalized bonding orbitals, or, alternatively, in terms of correlations of occupations of localised orbitals. However, in the latter case, multiorbital correlations were treated only in terms of two-orbital correlations, although the structure of multiorbital correlations is far richer; and, in the case of bonds established by more than two electrons, multiorbital correlations represent a more natural point of view. Here, for the first time, we introduce the true multiorbital correlation theory, consisting of a framework for handling the structure of multiorbital correlations, a toolbox of true multiorbital correlation measures, and the formulation of the multiorbital correlation clustering, together with an algorithm for obtaining that. These make it possible to characterise quantitatively, how well a bonding picture describes the chemical system. As proof of concept, we apply the theory for the investigation of the bond structures of several molecules. We show that the non-existence of well-defined multiorbital correlation clustering provides a reason for debated bonding picture.

Hückel-Hubbard-Ohno modeling of π-bonds in ethene and ethyne with application to trans-polyacetylene.

Quantum chemistry calculations provide the potential energy between two carbon atoms in ethane (H3C-CH3), ethene (H2C[double bond, length as m-dash]CH2), and ethyne (HC[triple bond, length as m-dash]CH) as a function of the atomic distance. Based on the energy function for the σ-bond in ethane, Vσ(r), we use the Hückel model with Hubbard-Ohno interaction for the π electrons to describe the energies Vσπ(r) and Vσππ(r) for the σπ double bond in ethene and the σππ triple bond in ethyne, respectively. The fit of the force functions shows that the electron transfer matrix element and the Peierls coupling can be estimated with some precision whereas the Hubbard-Ohno parameters are insignificant at the distances under consideration. We apply the Hückel-Hubbard-Ohno model to describe the bond lengths and the energies of elementary electronic excitations of trans-polyacetylene, (CH)n, whereby we adjust the σ-bond potential for conjugated polymers.

Orbital Entanglement in Bond-Formation Processes.

The accurate calculation of the (differential) correlation energy is central to the quantum chemical description of bond-formation and bond-dissociation processes. In order to estimate the quality of single- and multireference approaches for this purpose, various diagnostic tools have been developed. In this work, we elaborate on our previous observation [J. Phys. Chem. Lett.2012, 3, 3129] that one- and two-orbital-based entanglement measures provide quantitative means for the assessment and classification of electron correlation effects among molecular orbitals. The dissociation behavior of some prototypical diatomic molecules features all types of correlation effects relevant for chemical bonding. We demonstrate that our entanglement analysis is convenient to dissect these electron correlation effects and to provide a conceptual understanding of bond-forming and bond-breaking processes from the point of view of quantum information theory.