Variational Quantum Eigensolver Boosted by Adiabatic Connection.

In this work, we integrate the variational quantum eigensolver (VQE) with the adiabatic connection (AC) method for efficient simulations of chemical problems on near-term quantum computers. Orbital-optimized VQE methods are employed to capture the strong correlation within an active space, and classical AC corrections recover the dynamical correlation effects comprising electrons outside of the active space. On two challenging strongly correlated problems, namely, the dissociation of N2 and the electronic structure of the tetramethyleneethane biradical, we show that the combined VQE-AC approach enhances the performance of VQE dramatically. Moreover, since the AC corrections do not bring any additional requirements on quantum resources or measurements, they can actually boost the VQE algorithms. Our work paves the way toward quantum simulations of real-life problems on near-term quantum computers.

On-Surface Synthesis and Determination of the Open-Shell Singlet Ground State of Tridecacene.

The character of the electronic structure of acenes has been the subject of longstanding discussion. However, convincing experimental evidence of their open-shell character has so far been missing. Here, we present the on-surface synthesis of tridecacene molecules by thermal annealing of octahydrotridecacene on a Au(111) surface. We characterized the electronic structure of the tridecacene by scanning probe microscopy, which reveals the presence of an inelastic signal at 126 meV. We attribute the inelastic signal to spin excitation from the singlet diradical ground state to the triplet excited state. To rationalize the experimental findings, we carried out many-body ab initio calculations as well as model Hamiltonians to take into account the effect of the metallic substrate. Moreover, we provide a detailed analysis of how the dynamic electron correlation and virtual charge fluctuation between the molecule and metallic surface reduces the singlet-triplet band gap. Thus, this work provides the first experimental confirmation of the magnetic character of tridecacene.

Toward more accurate adiabatic connection approach for multireference wavefunctions.

A multiconfigurational adiabatic connection (AC) formalism is an attractive approach to compute the dynamic correlation within the complete active space self-consistent field and density matrix renormalization group (DMRG) models. Practical realizations of AC have been based on two approximations: (i) fixing one- and two-electron reduced density matrices (1- and 2-RDMs) at the zero-coupling constant limit and (ii) extended random phase approximation (ERPA). This work investigates the effect of removing the "fixed-RDM" approximation in AC. The analysis is carried out for two electronic Hamiltonian partitionings: the group product function- and the Dyall Hamiltonians. Exact reference AC integrands are generated from the DMRG full configuration interaction solver. Two AC models are investigated, employing either exact 1- and 2-RDMs or their second-order expansions in the coupling constant in the ERPA equations. Calculations for model molecules indicate that lifting the fixed-RDM approximation is a viable way toward improving the accuracy of existing AC approximations.

Hilbert space multireference coupled cluster tailored by matrix product states.

In the past decade, the quantum chemical version of the density matrix renormalization group method has established itself as the method of choice for strongly correlated molecular systems. However, despite its favorable scaling, in practice, it is not suitable for computations of dynamic correlation. Several approaches to include that in post-DMRG methods exist; in our group, we focused on the tailored coupled cluster (TCC) approach. This method works well in many situations; however, in exactly degenerate cases (with two or more determinants of equal weight), it exhibits a bias toward the reference determinant representing the Fermi vacuum. Although sometimes it is possible to use a compensation scheme to avoid this bias for energy differences, it is certainly a drawback. In order to overcome this bias of the TCC method, we have developed a Hilbert-space multireference version of tailored CC, which can treat several determinants on an equal footing. We have implemented and compared the performance of three Hilbert-space multireference coupled cluster (MRCC) variants-the state universal one and the Brillouin-Wigner and Mukherjee's state specific ones. We have assessed these approaches on the cyclobutadiene and tetramethyleneethane molecules, which are both diradicals with exactly degenerate determinants at a certain geometry. We have also investigated the sensitivity of the results on the orbital rotation of the highest occupied and lowest unoccupied molecular orbital (HOMO-LUMO) pair, as it is well known that Hilbert-space MRCC methods are not invariant to such transformations.

Interplay between π-Conjugation and Exchange Magnetism in One-Dimensional Porphyrinoid Polymers.

The synthesis of novel polymeric materials with porphyrinoid compounds as key components of the repeating units attracts widespread interest from several scientific fields in view of their extraordinary variety of functional properties with potential applications in a wide range of highly significant technologies. The vast majority of such polymers present a closed-shell ground state, and, only recently, as the result of improved synthetic strategies, the engineering of open-shell porphyrinoid polymers with spin delocalization along the conjugation length has been achieved. Here, we present a combined strategy toward the fabrication of one-dimensional porphyrinoid-based polymers homocoupled via surface-catalyzed [3 + 3] cycloaromatization of isopropyl substituents on Au(111). Scanning tunneling microscopy and noncontact atomic force microscopy describe the thermal-activated intra- and intermolecular oxidative ring closure reactions as well as the controlled tip-induced hydrogen dissociation from the porphyrinoid units. In addition, scanning tunneling spectroscopy measurements, complemented by computational investigations, reveal the open-shell character, that is, the antiferromagnetic singlet ground state (S = 0) of the formed polymers, characterized by singlet-triplet inelastic excitations observed between spins of adjacent porphyrinoid units. Our approach sheds light on the crucial relevance of the π-conjugation in the correlations between spins, while expanding the on-surface synthesis toolbox and opening avenues toward the synthesis of innovative functional nanomaterials with prospects in carbon-based spintronics.

On-Surface Synthesis and Characterization of [7]Triangulene Quantum Ring.

The ability to engineer geometrically well-defined antidots in large triangulene homologues allows for creating an entire family of triangulene quantum rings (TQRs) with tunable high-spin ground state, crucial for next-generation molecular spintronic devices. Herein, we report the synthesis of an open-shell [7]triangulene quantum ring ([7]TQR) molecule on Au(111) through the surface-assisted cyclodehydrogenation of a rationally designed kekulene derivative. Bond-resolved scanning tunneling microscopy (BR-STM) unambiguously imaged the molecular backbone of a single [7]TQR with a triangular zigzag edge topology, which can be viewed as [7]triangulene decorated with a coronene-like antidot in the center. Additionally, dI/dV mapping reveals that both inner and outer zigzag edges contribute to the edge-localized and spin-polarized electronic states of [7]TQR. Both experimental results and spin-polarized density functional theory calculations indicate that [7]TQR retains its open-shell septuple ground state (S = 3) on Au(111). This work demonstrates a new route for the design of high-spin graphene quantum rings for future quantum devices.

On-Surface Strain-Driven Synthesis of Nonalternant Non-Benzenoid Aromatic Compounds Containing Four- to Eight-Membered Rings.

The synthesis of polycyclic aromatic hydrocarbons containing various non-benzenoid rings remains a big challenge facing contemporary organic chemistry despite a considerable effort made over the last decades. Herein, we present a novel route, employing on-surface chemistry, to synthesize nonalternant polycyclic aromatic hydrocarbons containing up to four distinct kinds of non-benzenoid rings. We show that the surface-induced mechanical constraints imposed on strained helical reactants play a decisive role leading to the formation of products, energetically unfavorable in solution, with a peculiar ring current stabilizing the aromatic character of the π-conjugated system. Determination of the chemical and electronic structures of the most frequent product reveals its closed-shell character and low band gap. The present study renders a new route for the synthesis of novel nonalternant polycyclic aromatic hydrocarbons or other hydrocarbons driven by internal stress imposed by the surface not available by traditional approaches of organic chemistry in solution.

Massively parallel quantum chemical density matrix renormalization group method.

We present, to the best of our knowledge, the first attempt to exploit the super-computer platform for quantum chemical density matrix renormalization group (QC-DMRG) calculations. We have developed the parallel scheme based on the in-house MPI global memory library, which combines operator and symmetry sector parallelisms, and tested its performance on three different molecules, all typical candidates for QC-DMRG calculations. In case of the largest calculation, which is the nitrogenase FeMo cofactor cluster with the active space comprising 113 electrons in 76 orbitals and bond dimension equal to 6000, our parallel approach scales up to approximately 2000 CPU cores.

Quantum Chemistry in the Age of Quantum Computing.

Practical challenges in simulating quantum systems on classical computers have been widely recognized in the quantum physics and quantum chemistry communities over the past century. Although many approximation methods have been introduced, the complexity of quantum mechanics remains hard to appease. The advent of quantum computation brings new pathways to navigate this challenging and complex landscape. By manipulating quantum states of matter and taking advantage of their unique features such as superposition and entanglement, quantum computers promise to efficiently deliver accurate results for many important problems in quantum chemistry, such as the electronic structure of molecules. In the past two decades, significant advances have been made in developing algorithms and physical hardware for quantum computing, heralding a revolution in simulation of quantum systems. This Review provides an overview of the algorithms and results that are relevant for quantum chemistry. The intended audience is both quantum chemists who seek to learn more about quantum computing and quantum computing researchers who would like to explore applications in quantum chemistry.

Ground state of the Fe(ii)-porphyrin model system corresponds to quintet: a DFT and DMRG-based tailored CC study.

Fe(ii)-porphyrins play an important role in many reactions relevant to material science and biological processes, due to their closely lying spin states. Although the prevalent opinion is that these systems posses the triplet ground state, the recent experiment on Fe(ii)-phthalocyanine under conditions matching those of an isolated molecule points toward the quintet ground state. We present a thorough DFT and DMRG-based tailored CC study of Fe(ii)-porphyrin model, in which we address all previously discussed correlation effects. We examine the importance of geometrical parameters, the Fe-N distances in particular, and conclude that the system possesses the quintet ground state.

Toward DMRG-tailored coupled cluster method in the 4c-relativistic domain.

There are three essential problems in computational relativistic chemistry: Electrons moving at relativistic speeds, close lying states, and dynamical correlation. Currently available quantum-chemical methods are capable of solving systems with one or two of these issues. However, there is a significant class of molecules in which all the three effects are present. These are the heavier transition metal compounds, lanthanides, and actinides with open d or f shells. For such systems, sufficiently accurate numerical methods are not available, which hinders the application of theoretical chemistry in this field. In this paper, we combine two numerical methods in order to address this challenging class of molecules. These are the relativistic versions of coupled cluster methods and the density matrix renormalization group (DMRG) method. To the best of our knowledge, this is the first relativistic implementation of the coupled cluster method externally corrected by DMRG. The method brings a significant reduction of computational costs as we demonstrate on the system of TlH, AsH, and SbH.

Toward the efficient local tailored coupled cluster approximation and the peculiar case of oxo-Mn(Salen).

We introduce a new implementation of the coupled cluster method with single and double excitations tailored by the matrix product state wave functions (DMRG-TCCSD), which employs the local pair natural orbital (LPNO) approach. By exploiting locality in the coupled cluster stage of the calculation, we were able to remove some of the limitations that hindered the application of the canonical version of the method to larger systems and/or with larger basis sets. We assessed the accuracy of the approximation using two systems: tetramethyleneethane (TME) and oxo-Mn(Salen). Using the default cut-off parameters, we were able to recover over 99.7% and 99.8% of the canonical correlation energy for the triplet and singlet state of TME, respectively. In the case of oxo-Mn(Salen), we found that the amount of retrieved canonical correlation energy depends on the size of the complete active space (CAS)-we retrieved over 99.6% for the larger 27 orbital CAS and over 99.8% for the smaller 22 orbital CAS. The use of LPNO-TCCSD allowed us to perform these calculations up to quadruple-ζ basis set, amounting to 1178 basis functions. Moreover, we examined dependence of the ground state of oxo-Mn(Salen) on the CAS composition. We found that the inclusion of 4dxy orbital plays an important role in stabilizing the singlet state at the DMRG-CASSCF level via double-shell effect. However, by including dynamic correlation, the ground state was found to be triplet regardless of the size of the basis set or the composition of CAS, which is in agreement with previous findings by canonical DMRG-TCCSD in smaller basis.

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.

An Isolated Molecule of Iron(II) Phthalocyanin Exhibits Quintet Ground-State: A Nexus between Theory and Experiment.

Iron(II) phthalocyanine (FePc) is an important member of the phthalocyanines family with potential applications in the fields of electrocatalysis, magnetic switching, electrochemical sensing, and phototheranostics. Despite the importance of electronic properties of FePc in these applications, a reliable determination of its ground-state is still challenging. Here we present combined state of the art computational methods and experimental approaches, that is, Mössbauer spectroscopy and Superconducting Quantum Interference Device (SQUID) magnetic measurements to identify the ground state of FePc. While the nature of the ground state obtained with density functional theory (DFT) depends on the functional, giving mostly the triplet state, multi-reference complete active space second-order perturbation theory (CASPT2) and density matrix renormalization group (DMRG) methods assign quintet as the FePc ground-state in gas-phase. This has been confirmed by the hyperfine parameters obtained from 57 Fe Mössbauer spectroscopy performed in frozen monochlorobenzene. The use of monochlorobenzene guarantees an isolated nature of the FePc as indicated by a zero Weiss temperature. The results open doors for exploring the ground state of other metal porphyrin molecules and their controlled spin transitions via external stimuli.

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.

Adiabatic state preparation study of methylene.

Quantum computers attract much attention as they promise to outperform their classical counterparts in solving certain type of problems. One of them with practical applications in quantum chemistry is simulation of complex quantum systems. An essential ingredient of efficient quantum simulation algorithms are initial guesses of the exact wave functions with high enough fidelity. As was proposed in Aspuru-Guzik et al. [Science 309, 1704 (2005)], the exact ground states can in principle be prepared by the adiabatic state preparation method. Here, we apply this approach to preparation of the lowest lying multireference singlet electronic state of methylene and numerically investigate preparation of this state at different molecular geometries. We then propose modifications that lead to speeding up the preparation process. Finally, we decompose the minimal adiabatic state preparation employing the direct mapping in terms of two-qubit interactions.

Highly entangled polyradical nanographene with coexisting strong correlation and topological frustration.

Open-shell nanographenes exhibit unconventional π-magnetism arising from topological frustration or strong electron-electron interaction. However, conventional design approaches are typically limited to a single magnetic origin, which can restrict the number of correlated spins or the type of magnetic ordering in open-shell nanographenes. Here we present a design strategy that combines topological frustration and electron-electron interactions to fabricate a large fully fused 'butterfly'-shaped tetraradical nanographene on Au(111). We employ bond-resolved scanning tunnelling microscopy and spin-excitation spectroscopy to resolve the molecular backbone and reveal the strongly correlated open-shell character, respectively. This nanographene contains four unpaired electrons with both ferromagnetic and anti-ferromagnetic interactions, harbouring a many-body singlet ground state and strong multi-spin entanglement, which is well described by many-body calculations. Furthermore, we study the magnetic properties and spin states in the nanographene using a nickelocene magnetic probe. The ability to imprint and characterize many-body strongly correlated spins in polyradical nanographenes paves the way for future advancements in quantum information technologies.

Projection-Based Density Matrix Renormalization Group in Density Functional Theory Embedding.

The density matrix renormalization group (DMRG) method has already proved itself as a very efficient and accurate computational method, which can treat large active spaces and capture the major part of strong correlation. Its application on larger molecules is, however, limited by its own computational scaling as well as demands of methods for treatment of the missing dynamical electron correlation. In this work, we present the first step in the direction of combining DMRG with density functional theory (DFT), one of the most employed quantum chemical methods with favorable scaling, by means of the projection-based wave function (WF)-in-DFT embedding. On two proof-of-concept but important molecular examples, we demonstrate that the developed DMRG-in-DFT approach provides a very accurate description of molecules with a strongly correlated fragment.

Role of Spin Polarization and Dynamic Correlation in Singlet-Triplet Gap Inversion of Heptazine Derivatives.

The new generation of proposed light-emitting molecules for organic light-emitting diodes (OLEDs) has raised considerable research interest due to its exceptional feature─a negative singlet-triplet (ST) gap violating Hund's multiplicity rule in the excited S1 and T1 states. We investigate the role of spin polarization in the mechanism of ST gap inversion. Spin polarization is associated with doubly excited determinants of certain types, whose presence in the wave function expansion favors the energy of the singlet state more than that of the triplet. Using a perturbation theory-based model for spin polarization, we propose a simple descriptor for prescreening of candidate molecules with negative ST gaps and prove its usefulness for heptazine-type molecules. Numerical results show that the quantitative effect of spin polarization decreases linearly with the increasing highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) exchange integral. Comparison of single- and multireference coupled-cluster predictions of ST gaps shows that the former methods provide good accuracy by correctly balancing the effects of doubly excited determinants and dynamic correlation. We also show that accurate ST gaps may be obtained using a complete active space model supplemented with dynamic correlation from multireference adiabatic connection theory.