Aromaticity is a well-known phenomenon in both physics and chemistry, and is responsible for many unique chemical and physical properties of aromatic molecules. The primary feature contributing to the stability of polycyclic aromatic hydrocarbons is the delocalised π -electron clouds in the 2 p z orbitals of each of the N carbon atoms. While it is known that electrons delocalize among the hybridized sp 2 orbitals, this paper proposes quantum walk as the mechanism by which the delocalization occurs, and also obtains how the functional chemical structures of these molecules arise naturally out of such a construction. We present results of computations performed for some benzoid polycyclic aromatic hydrocarbons in this regard, and show that the quantum walk-based approach does correctly predict the reactive sites and stability order of the molecules considered.
Aromaticity is a well-known phenomenon in both physics and chemistry, and is responsible for many unique chemical and physical properties of aromatic molecules. The primary feature contributing to the stability of polycyclic aromatic hydrocarbons is the delocalised π -electron clouds in the 2 p z orbitals of each of the N carbon atoms. While it is known that electrons delocalize among the hybridized sp 2 orbitals, this paper proposes quantum walk as the mechanism by which the delocalization occurs, and also obtains how the functional chemical structures of these molecules arise naturally out of such a construction. We present results of computations performed for some benzoid polycyclic aromatic hydrocarbons in this regard, and show that the quantum walk-based approach does correctly predict the reactive sites and stability order of the molecules considered.
This letter investigates the applications and extensions of the resource theory of quantum superposition within the realm of quantum chemistry. Specifically, we explore aromaticity, a fundamental concept originally developed to elucidate the structural symmetry, energetic stability, and chemical reactivity of benzene and its derivatives. While aromaticity and its counterpart, antiaromaticity, are associated with the delocalization of electrons between nonorthogonal atomic orbitals, they lack a universally accepted and comprehensive definition. We demonstrate that the genuine quantum superposition exhibited by biorthogonal atomic orbitals effectively captures the aromaticity order of representative monocyclic molecules. These findings reveal that the quantum resource theories hold significant implications, offering fresh insights into our comprehension of chemical bonding phenomena.
Quantum Embedding density functional theory (QmDFT) embedding offers a highly scalable approach to improve treatment for large highly correlated pi conjugated systems. However, estimating advanced electronic structure properties in polycyclic aromatic hydrocarbons (PAHs) needs advanced exchange correlation functionals that frequently trigger convergence instabilities during the embedding cycle.
In this work, we introduce an adaptive damping and direct inversion in the iterative subspace (DIIS) accelerated protocol that stabilizes the embedding procedure, enabling robust integration of hybrid functionals like B3LYP and CAM B3LYP. Using 10 selected PAHs (linear and fused) molecules as a benchmark. We demonstrate a clear functional-dependent ground-state energetics and frontier-orbital gap estimation. While LDA based approaches yield near-quantitative agreement with FCI in DFT reference energies and is further supported by thermochemical isomerization benchmarks, while B3LYP provide significantly improved agreement with experimental E0-0 transition values. This mapping allows us to bypass explicit excited-state calculations for E0 0 values, thereby significantly reducing computational overhead.
Among the hybrid functionals screened, CAM B3LYP offers a balanced overall performance. Our results establish a stable QmDFT framework and provide useful guidance for functional selection for quantum embedding studies of PAHs and related low-dimensional pi-conjugated materials.
Quantum Monte Carlo (QMC) is one of the most promising methods for solving quantum many-body problems. QMC combines known analytical properties of wave functions, results from other methods, and stochastic techniques, into a powerful tool for investigation of the electronic structure of real materials. Furthermore, the method has a favorable scaling ( O[N3] in the number of electrons), perfect scalability on parallel architectures, and can be applied to a wide range of systems. Our work has focused both on expanding the capabilities of QMC for calculation of quantities beyond energies, such as interatomic forces, and on application of the method to challenging problems in nanostructure materials research. For the calculation of interatomic forces we have implemented a finite difference correlated sampling method. The correlated sampling enables us to avoid the problem of statistical noise in evaluating the energy differences and allows us to obtain forces with an average 1% error. We have analyzed the effect of the method's main approximations on the estimate of the forces, both by theoretical arguments and by performing a variety of tests on small systems. Furthermore, we have used this newly developed method for estimating equilibrium geometries in a problem where other methods have failed. QMC enabled us to treat the electron correlation effects with high accuracy and to provide accurate predictions in a number of interesting applications. In particular, for carbon rings of C4N+2 stoichiometry we have resolved a long-standing problem regarding the competition between the aromaticity and second order Jahn-Teller effects. We have found that the dominant mechanism at small and intermediate sizes, which stabilizes the bond angle and bond length alternated geometries, is the secondorder Jahn-Teller effect while the aromatic isomer is always found to be a transition state. These high accuracy calculations which have involved multi-reference wave functions, extensive variational optimizations, and the use of massively parallel platforms demonstrate clearly the unique capabilities of the QMC approach.
Quantum walks serve as novel tools for performing efficient quantum computation and simulation. In a recent experimental demonstration [1] we have realized photonic quantum walks for simulating cyclic quantum systems, such as hexagonal lattices or aromatic molecules like benzene. In that experiment we explored the wave function dynamics and the probability distribution of a quantum particle located on a six-site system (with periodic boundary conditions), alongside with simpler demonstration of three- and four-site systems, under various initial conditions. Localization and revival of the wave function were demonstrated. After revisiting that experiment we will theoretically analyze the case of noisy quantum walks by implementing the bit-phase flip channel. This will allow us to draw conclusions regarding the performance of our photonic quantum simulation in noisy environments. Finally, we will briefly outline some future directions.
The computation of excited electronic states is an important application for quantum computers. In this work, we simulate the excited state spectra of four aromatic heterocycles on IBM superconducting quantum computers, focussing on active spaces of π→π∗ and n→π∗ excitations. We approximate the ground state with the entanglement forging method, a qubit reduction technique that maps a spatial orbital to a single qubit, rather than two qubits. We then determine excited states using the quantum subspace expansion method. We showcase these algorithms on quantum hardware using up to 8 qubits and employing readout and gate error mitigation techniques. Our results demonstrate a successful application of quantum computing in the simulation of active-space electronic wavefunctions of substituted aromatic heterocycles, and outline challenges to be overcome in elucidating the optical properties of organic molecules with hybrid quantum-classical algorithms.
The chemical bond is an important local concept to understand chemical compounds and processes. Unfortunately, like most local concepts, the chemical bond and the bond order do not correspond to any physical observable and thus cannot be determined as an expectation value of a quantum chemical operator. We recently demonstrated [Boguslawski et al., J. Chem. Theory Comput., 2013, 9, 2959-2973] that one- and two-orbital-based entanglement measures can be applied to interpret electronic wave functions in terms of orbital correlation. Orbital entanglement emerged as a powerful tool to provide a qualitative understanding of bond-forming and bond-breaking processes, and allowed for an estimation of bond orders of simple diatomic molecules beyond the classical bonding models. In this article we demonstrate that the orbital entanglement analysis can be extended to polyatomic molecules to understand chemical bonding.