The odd-number cyclo[13]carbon and its dimer, cyclo[26]carbon.

Molecular rings of N carbon atoms (cyclo[N]carbons, or CN) are excellent benchmarking systems for testing quantum chemical theoretical methods and valuable precursors to other carbon-rich materials. Odd-N cyclocarbons, which have been elusive to date, are predicted to be even less stable than even-N cyclocarbons. We report the on-surface synthesis of cyclo[13]carbon, C13, by manipulation of decachlorofluorene with a scanning probe microscope tip. We elucidated the properties of C13 by experiment and theoretical modeling. C13 adopts an open-shell configuration with a triplet ground state and a kinked geometry, which shows different extents of distortion and carbene localization depending on the molecular environment. Moreover, we prepared and characterized the C13 dimer, cyclo[26]carbon, demonstrating the potential of cyclocarbons and their precursors as building blocks for carbon allotropes.

Toward Real Chemical Accuracy on Current Quantum Hardware Through the Transcorrelated Method.

Quantum computing is emerging as a new computational paradigm with the potential to transform several research fields including quantum chemistry. However, current hardware limitations (including limited coherence times, gate infidelities, and connectivity) hamper the implementation of most quantum algorithms and call for more noise-resilient solutions. We propose an explicitly correlated Ansatz based on the transcorrelated (TC) approach to target these major roadblocks directly. This method transfers, without any approximation, correlations from the wave function directly into the Hamiltonian, thus reducing the resources needed to achieve accurate results with noisy quantum devices. We show that the TC approach allows for shallower circuits and improves the convergence toward the complete basis set limit, providing energies within chemical accuracy to experiment with smaller basis sets and, thus, fewer qubits. We demonstrate our method by computing bond lengths, dissociation energies, and vibrational frequencies close to experimental results for the hydrogen dimer and lithium hydride using two and four qubits, respectively. To demonstrate our approach's current and near-term potential, we perform hardware experiments, where our results confirm that the TC method paves the way toward accurate quantum chemistry calculations already on today's quantum hardware.

Ultrafast fragmentation of highly-excited doubly-ionized deoxyribose: role of the liquid water environment.

Ab initio molecular dynamics simulations are used to investigate the fragmentation dynamics following the double ionization of 2-deoxy-D-ribose (DR), a major component in the DNA chain. Different ionization scenarios are considered to provide a complete picture. First focusing on isolated DR2+, fragmentation patterns are determined for the ground electronic state, adding randomly distributed excitation energy to the nuclei. These patterns differ for the two isomers studied. To compare thermal and electronic excitation effects, Ehrenfest dynamics are also performed, allowing to remove the two electrons from selected molecular orbitals. Two intermediate-energy orbitals, localized on the carbon chain, were selected. The dissociation pattern corresponds to the most frequent pattern obtained when adding thermal excitation. On the contrary, targeting the four deepest orbitals, localized on the oxygen atoms, leads to selective ultrafast C-O and/or O-H bond dissociation. To probe the role of environment, a system consisting of a DR molecule embedded in liquid water is then studied. The two electrons are removed from either the DR or the water molecules directly linked to the sugar through hydrogen bonds. Although the dynamics onset is similar to that of isolated DR when removing the same deep orbitals localized on the sugar oxygen atoms, the subsequent fragmentation patterns differ. Sugar damage also occurs following the Coulomb explosion of neighboring H2O2+ molecules due to interaction with the emitted O or H atoms.

Toward Accurate Post-Born-Oppenheimer Molecular Simulations on Quantum Computers: An Adaptive Variational Eigensolver with Nuclear-Electronic Frozen Natural Orbitals.

Nuclear quantum effects such as zero-point energy and hydrogen tunneling play a central role in many biological and chemical processes. The nuclear-electronic orbital (NEO) approach captures these effects by treating selected nuclei quantum mechanically on the same footing as electrons. On classical computers, the resources required for an exact solution of NEO-based models grow exponentially with system size. By contrast, quantum computers offer a means of solving this problem with polynomial scaling. However, due to the limitations of current quantum devices, NEO simulations are confined to the smallest systems described by minimal basis sets, whereas realistic simulations beyond the Born-Oppenheimer approximation require more sophisticated basis sets. For this purpose, we herein extend a hardware-efficient ADAPT-VQE method to the NEO framework in the frozen natural orbital (FNO) basis. We demonstrate on H2 and D2 molecules that the NEO-FNO-ADAPT-VQE method reduces the CNOT count by several orders of magnitude relative to the NEO unitary coupled cluster method with singles and doubles while maintaining the desired accuracy. This extreme reduction in the CNOT gate count is sufficient to permit practical computations employing the NEO method─an important step toward accurate simulations involving nonclassical nuclei and non-Born-Oppenheimer effects on near-term quantum devices. We further show that the method can capture isotope effects, and we demonstrate that inclusion of correlation energy systematically improves the prediction of difference in the zero-point energy (ΔZPE) between isotopes.

Aromaticity Reversal Induced by Vibrations in Cyclo[16]carbon.

Aromaticity is typically regarded as an intrinsic property of a molecule, correlated with electron delocalization, stability, and other properties. Small variations in the molecular geometry usually result in small changes in aromaticity, in line with Hammond's postulate. For example, introducing bond-length alternation in benzene and square cyclobutadiene by modulating the geometry along the Kekulé vibration gradually decreases the magnitude of their ring currents, making them less aromatic and less antiaromatic, respectively. A sign change in the ring current, corresponding to a reversal of aromaticity, typically requires a gross perturbation such as electronic excitation, addition or removal of two electrons, or a dramatic change in the molecular geometry. Here, we use multireference calculations to show how movement along the Kekulé vibration, which controls bond-length alternation, induces a sudden reversal in the ring current of cyclo[16]carbon, C16. This reversal occurs when the two orthogonal π systems of C16 sustain opposing currents. These results are rationalized by a Hückel model which includes bond-length alternation, and which is combined with a minimal model accounting for orbital contributions to the ring current. Finally, we successfully describe the electronic structure of C16 with a "divide-and-conquer" approach suitable for execution on a quantum computer.

On-surface synthesis of a doubly anti-aromatic carbon allotrope.

Synthetic carbon allotropes such as graphene1, carbon nanotubes2 and fullerenes3 have revolutionized materials science and led to new technologies. Many hypothetical carbon allotropes have been discussed4, but few have been studied experimentally. Recently, unconventional synthetic strategies such as dynamic covalent chemistry5 and on-surface synthesis6 have been used to create new forms of carbon, including γ-graphyne7, fullerene polymers8, biphenylene networks9 and cyclocarbons10,11. Cyclo[N]carbons are molecular rings consisting of N carbon atoms12,13; the three that have been reported to date (N = 10, 14 and 18)10,11 are doubly aromatic, which prompts the question: is it possible to prepare doubly anti-aromatic versions? Here we report the synthesis and characterization of an anti-aromatic carbon allotrope, cyclo[16]carbon, by using tip-induced on-surface chemistry6. In addition to structural information from atomic force microscopy, we probed its electronic structure by recording orbital density maps14 with scanning tunnelling microscopy. The observation of bond-length alternation in cyclo[16]carbon confirms its double anti-aromaticity, in concordance with theory. The simple structure of C16 renders it an interesting model system for studying the limits of aromaticity, and its high reactivity makes it a promising precursor to novel carbon allotropes15.

Nonadiabatic Nuclear-Electron Dynamics: A Quantum Computing Approach.

Coupled quantum electron-nuclear dynamics is often associated with the Born-Huang expansion of the molecular wave function and the appearance of nonadiabatic effects as a perturbation. On the other hand, native multicomponent representations of electrons and nuclei also exist, which do not rely on any a priori approximation. However, their implementation is hampered by prohibitive scaling. Consequently, quantum computers offer a unique opportunity for extending their use to larger systems. Here, we propose a quantum algorithm for simulating the time-evolution of molecular systems and apply it to proton transfer dynamics in malonaldehyde, described as a rigid scaffold. The proposed quantum algorithm can be easily generalized to include the explicit dynamics of the classically described molecular scaffold. We show how entanglement between electronic and nuclear degrees of freedom can persist over long times if electrons do not follow the nuclear displacement adiabatically. The proposed quantum algorithm may become a valid candidate for the study of such phenomena when sufficiently powerful quantum computers become available.

Spin-Flip Unitary Coupled Cluster Method: Toward Accurate Description of Strong Electron Correlation on Quantum Computers.

Quantum computers have emerged as a promising platform to simulate strong electron correlation that is crucial to catalysis and photochemistry. However, owing to the choice of a trial wave function employed in the variational quantum eigensolver (VQE) algorithm, accurate simulation is restricted to certain classes of correlated phenomena. Herein, we combine the spin-flip (SF) formalism with the unitary coupled cluster with singles and doubles (UCCSD) method via the quantum equation-of-motion (qEOM) approach to allow for an efficient simulation of a large family of strongly correlated problems. We show that the developed qEOM-SF-UCCSD/VQE method outperforms its UCCSD/VQE counterpart for simulation of the cis-trans isomerization of ethylene, and the automerization of cyclobutadiene and the predicted qEOM-SF-UCCSD/VQE barrier heights are in a good agreement with the experimentally determined values. The developments presented herein will further stimulate the investigation of this approach for simulations of other types of correlated/entangled phenomena on quantum computers.

A quantum computing implementation of nuclearelectronic orbital (NEO) theory: Toward an exact pre-Born-Oppenheimer formulation of molecular quantum systems.

Nuclear quantum phenomena beyond the Born-Oppenheimer approximation are known to play an important role in a growing number of chemical and biological processes. While there exists no unique consensus on a rigorous and efficient implementation of coupled electron-nuclear quantum dynamics, it is recognized that these problems scale exponentially with system size on classical processors and, therefore, may benefit from quantum computing implementations. Here, we introduce a methodology for the efficient quantum treatment of the electron-nuclear problem on near-term quantum computers, based upon the Nuclear-Electronic Orbital (NEO) approach. We generalize the electronic two-qubit tapering scheme to include nuclei by exploiting symmetries inherent in the NEO framework, thereby reducing the Hamiltonian dimension, number of qubits, gates, and measurements needed for calculations. We also develop parameter transfer and initialization techniques, which improve convergence behavior relative to conventional initialization. These techniques are applied to H2 and malonaldehyde for which results agree with NEO full configuration interaction and NEO complete active space configuration interaction benchmarks for ground state energy to within 10-6 hartree and entanglement entropy to within 10-4. These implementations therefore significantly reduce resource requirements for full quantum simulations of molecules on near-term quantum devices while maintaining high accuracy.

Optimizing Quantum Classification Algorithms on Classical Benchmark Datasets.

The discovery of quantum algorithms offering provable advantages over the best known classical alternatives, together with the parallel ongoing revolution brought about by classical artificial intelligence, motivates a search for applications of quantum information processing methods to machine learning. Among several proposals in this domain, quantum kernel methods have emerged as particularly promising candidates. However, while some rigorous speedups on certain highly specific problems have been formally proven, only empirical proof-of-principle results have been reported so far for real-world datasets. Moreover, no systematic procedure is known, in general, to fine tune and optimize the performances of kernel-based quantum classification algorithms. At the same time, certain limitations such as kernel concentration effects-hindering the trainability of quantum classifiers-have also been recently pointed out. In this work, we propose several general-purpose optimization methods and best practices designed to enhance the practical usefulness of fidelity-based quantum classification algorithms. Specifically, we first describe a data pre-processing strategy that, by preserving the relevant relationships between data points when processed through quantum feature maps, substantially alleviates the effect of kernel concentration on structured datasets. We also introduce a classical post-processing method that, based on standard fidelity measures estimated on a quantum processor, yields non-linear decision boundaries in the feature Hilbert space, thus achieving the quantum counterpart of the radial basis functions technique that is widely employed in classical kernel methods. Finally, we apply the so-called quantum metric learning protocol to engineer and adjust trainable quantum embeddings, demonstrating substantial performance improvements on several paradigmatic real-world classification tasks.

Quantum Embedding Method for the Simulation of Strongly Correlated Systems on Quantum Computers.

Quantum computing has emerged as a promising platform for simulating strongly correlated systems in chemistry, for which the standard quantum chemistry methods are either qualitatively inaccurate or too expensive. However, due to the hardware limitations of the available noisy near-term quantum devices, their application is currently limited only to small chemical systems. One way for extending the range of applicability can be achieved within the quantum embedding approach. Herein, we employ the projection-based embedding method for combining the variational quantum eigensolver (VQE) algorithm, although not limited to, with density functional theory (DFT). The developed VQE-in-DFT method is then implemented efficiently on a real quantum device and employed for simulating the triple bond breaking process in butyronitrile. The results presented herein show that the developed method is a promising approach for simulating systems with a strongly correlated fragment on a quantum computer.

Quantum algorithms for quantum dynamics.

Among the many computational challenges faced across different disciplines, quantum-mechanical systems pose some of the hardest ones and offer a natural playground for the growing field of quantum technologies. In this Perspective, we discuss quantum algorithmic solutions for quantum dynamics, reporting on the latest developments and offering a viewpoint on their potential and current limitations. We present some of the most promising areas of application and identify possible research directions for the coming years.

Selectivity in single-molecule reactions by tip-induced redox chemistry.

Controlling selectivity of reactions is an ongoing quest in chemistry. In this work, we demonstrate reversible and selective bond formation and dissociation promoted by tip-induced reduction-oxidation reactions on a surface. Molecular rearrangements leading to different constitutional isomers are selected by the polarity and magnitude of applied voltage pulses from the tip of a combined scanning tunneling and atomic force microscope. Characterization of voltage dependence of the reactions and determination of reaction rates demonstrate selectivity in constitutional isomerization reactions and provide insight into the underlying mechanisms. With support of density functional theory calculations, we find that the energy landscape of the isomers in different charge states is important to rationalize the selectivity. Tip-induced selective single-molecule reactions increase our understanding of redox chemistry and could lead to novel molecular machines.

Molecular Quantum Dynamics: A Quantum Computing Perspective.

ConspectusSimulating molecular dynamics (MD) within a comprehensive quantum framework has been a long-standing challenge in computational chemistry. An exponential scaling of computational cost renders solving the time dependent Schrödinger equation (TDSE) of a molecular Hamiltonian, including both electronic and nuclear degrees of freedom (DOFs), as well as their couplings, infeasible for more than a few DOFs. In the Born-Oppenheimer (BO), or adiabatic, picture, electronic and nuclear parts of the wave function are decoupled and treated separately. Within this framework, the nuclear wave function evolves along potential energy surfaces (PESs) computed as solutions to the electronic Schrödinger equation parametrized in the nuclear DOFs. This approximation, together with increasingly elaborate numerical approaches to solve the nuclear time dependent Schrödinger equation (TDSE), enabled the treatment of up to a few dozens of degrees of freedom (DOFs). However, for particular applications, such as photochemistry, the BO approximation breaks down. In this regime of non-adiabatic dynamics, solving the full molecular problem including electron-nuclear couplings becomes essential, further increasing the complexity of the numerical solution. Although valuable methods such as multiconfigurational time-dependent Hartree (MCTDH) have been proposed for the solution of the coupled electron-nuclear dynamics, they remain hampered by an exponential scaling in the number of nuclear DOFs and by the difficulty of finding universal variational forms.In this Account, we present a perspective on novel quantum computational algorithms, aiming to alleviate the exponential scaling inherent to the simulation of many-body quantum dynamics. In particular, we focus on the derivation and application of quantum algorithms for adiabatic and non-adiabatic quantum dynamics, which include efficient approaches for the calculation of the BO potential energy surfaces (PESs). Thereafter, we study the time-evolution of a model system consisting of two coupled PESs in first and second quantization. In a first application, we discuss a recently introduced quantum algorithm for the evolution of a wavepacket in first quantization and exploit the potential quantum advantage of mapping its spatial grid representation to logarithmically many qubits. For the second demonstration, we move to the second quantization framework and review the scaling properties of two alternative time-evolution algorithms, namely, a variational quantum algorithm (VQA) (based on the McLachlan variational principle) and conventional Trotter-type evolution (based on a Lie-Trotter-Suzuki formula). Both methods clearly demonstrate the potential of quantum algorithms and their favorable scaling compared to the available classical approaches. However, a clear demonstration of quantum advantage in the context of molecular quantum dynamics may require the implementation of these algorithms in fault-tolerant quantum computers, while their application in near-term, noisy quantum devices is still unclear and deserves further investigation.

Quantum algorithm for alchemical optimization in material design.

The development of tailored materials for specific applications is an active field of research in chemistry, material science and drug discovery. The number of possible molecules obtainable from a set of atomic species grow exponentially with the size of the system, limiting the efficiency of classical sampling algorithms. On the other hand, quantum computers can provide an efficient solution to the sampling of the chemical compound space for the optimization of a given molecular property. In this work, we propose a quantum algorithm for addressing the material design problem with a favourable scaling. The core of this approach is the representation of the space of candidate structures as a linear superposition of all possible atomic compositions. The corresponding 'alchemical' Hamiltonian drives the optimization in both the atomic and electronic spaces leading to the selection of the best fitting molecule, which optimizes a given property of the system, e.g., the interaction with an external potential as in drug design. The quantum advantage resides in the efficient calculation of the electronic structure properties together with the sampling of the exponentially large chemical compound space. We demonstrate both in simulations and with IBM Quantum hardware the efficiency of our scheme and highlight the results in a few test cases. This preliminary study can serve as a basis for the development of further material design quantum algorithms for near-term quantum computers.

Improved Accuracy on Noisy Devices by Nonunitary Variational Quantum Eigensolver for Chemistry Applications.

We propose a modification of the Variational Quantum Eigensolver algorithm for electronic structure optimization using quantum computers, named nonunitary Variational Quantum Eigensolver (nu-VQE), in which a nonunitary operator is combined with the original system Hamiltonian leading to a new variational problem with a simplified wave function ansatz. In the present work, as nonunitary operator, we use the Jastrow factor, inspired from classical Quantum Monte Carlo techniques for simulation of strongly correlated electrons. The method is applied to prototypical molecular Hamiltonians for which we obtain accurate ground-state energies with shallower circuits, at the cost of an increased number of measurements. Finally, we also show that this method achieves an important error mitigation effect that drastically improves the quality of the results for VQE optimizations on today's noisy quantum computers. The absolute error in the calculated energy within our scheme is 1 order of magnitude smaller than the corresponding result using traditional VQE methods, with the same circuit depth.

Probing Molecular Excited States by Atomic Force Microscopy.

By employing single charge injections with an atomic force microscope, we investigated redox reactions of a molecule on a multilayer insulating film. First, we charged the molecule positively by attaching a single hole. Then we neutralized it by attaching an electron and observed three channels for the neutralization. We rationalize that the three channels correspond to transitions to the neutral ground state, to the lowest energy triplet excited states and to the lowest energy singlet excited states. By single-electron tunneling spectroscopy we measured the energy differences between the transitions obtaining triplet and singlet excited state energies. The experimental values are compared with density functional theory calculations of the excited state energies. Our results show that molecules in excited states can be prepared and that energies of optical gaps can be quantified by controlled single-charge injections. Our work demonstrates the access to, and provides insight into, ubiquitous electron-attachment processes related to excited-state transitions important in electron transfer and molecular optoelectronics phenomena on surfaces.

Quantum HF/DFT-embedding algorithms for electronic structure calculations: Scaling up to complex molecular systems.

In the near future, material and drug design may be aided by quantum computer assisted simulations. These have the potential to target chemical systems intractable by the most powerful classical computers. However, the resources offered by contemporary quantum computers are still limited, restricting the simulations to very simple molecules. In order to rapidly scale up to more interesting molecular systems, we propose the embedding of the quantum electronic structure calculation into a classically computed environment obtained at the Hartree-Fock (HF) or density functional theory (DFT) level of theory. This result is achieved by constructing an effective Hamiltonian that incorporates a mean field potential describing the action of the inactive electrons on a selected Active Space (AS). The ground state of the AS Hamiltonian is then determined by means of the variational quantum eigensolver algorithm. We show that with the proposed HF and DFT embedding schemes, we can obtain significant energy corrections to the reference HF and DFT calculations for a number of simple molecules in their strongly correlated limit (the dissociation regime) as well as for systems of the size of the oxirane molecule.

Quantum equilibration of the double-proton transfer in a model system porphine.

There is a renewed interest in the derivation of statistical mechanics from the dynamics of closed quantum systems. A central part of this program is to understand how closed quantum systems, i.e., in the absence of a thermal bath, initialized far-from-equilibrium can share a dynamics that is typical to the relaxation towards thermal equilibrium. Equilibration dynamics has been traditionally studied with a focus on the so-called quenches of large-scale many-body systems. We consider here the equilibration of a two-dimensional molecular model system describing the double proton transfer reaction in porphine. Using numerical simulations, we show that equilibration indeed takes place very rapidly (∼200 fs) for initial states induced by pump-dump laser pulse control with energies well above the synchronous barrier. The resulting equilibration state is characterized by a strong delocalization of the probability density of the protons that can be explained, mechanistically, as the result of (i) an initial state consisting of a large superposition of vibrational states, and (ii) the presence of a very effective dephasing mechanism.

Redox Properties of Native and Damaged DNA from Mixed Quantum Mechanical/Molecular Mechanics Molecular Dynamics Simulations.

The redox properties of two large DNA fragments composed of 39 base pairs, differing only by an 8-oxoguanine (8oxoG) defect replacing a guanine (G), were investigated in physiological conditions using mixed quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations. The quantum region of the native fragment comprised 3 G-C base pairs, while one G was replaced by an 8oxoG in the defect fragment. The calculated values for the redox free energy are 6.55 ± 0.28 eV and 5.62 ± 0.30 eV for the native and the 8oxoG-containing fragment, respectively. The respective estimates for the reorganization free energy are 1.25 ± 0.18 eV and 1.00 ± 0.18 eV. Both reactions follow the Marcus theory for electron transfer. The large difference in redox potential between the two fragments shows that replacement of a G by an 8oxoG renders the DNA more easily oxidizable. This finding is in agreement with the suggestion that DNA fragments containing an 8oxoG defect can act as sinks of oxidative damage that protect the rest of the genome from assault. In addition, the difference in redox potential between the native and the defect DNA fragment indicates that a charge transfer-based mechanism for the recognition of DNA defects might be feasible, in line with recent suggestions based on experimental observations.