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.

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 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.