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.

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.

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.

Design, synthesis, biological evaluation and X-ray structural studies of potent human dihydroorotate dehydrogenase inhibitors based on hydroxylated azole scaffolds.

A new generation of potent hDHODH inhibitors designed by a scaffold-hopping replacement of the quinolinecarboxylate moiety of brequinar, one of the most potent known hDHODH inhibitors, is presented here. Their general structure is characterized by a biphenyl moiety joined through an amide bridge with an acidic hydroxyazole scaffold (hydroxylated thiadiazole, pyrazole and triazole). Molecular modelling suggested that these structures should adopt a brequinar-like binding mode involving interactions with subsites 1, 2 and 4 of the hDHODH binding site. Initially, the inhibitory activity of the compounds was studied on recombinant hDHODH. The most potent compound of the series in the enzymatic assays was the thiadiazole analogue 4 (IC50 16 nM). The activity was found to be dependent on the fluoro substitution pattern at the biphenyl moiety as well as on the choice/substitution of the heterocyclic ring. Structure determination of hDHODH co-crystallized with one representative compound from each series (4, 5 and 6) confirmed the brequinar-like binding mode as suggested by modelling. The specificity of the observed effects of the compound series was tested in cell-based assays for antiproliferation activity using Jurkat cells and PHA-stimulated PBMC. These tests were also verified by addition of exogenous uridine to the culture medium. In particular, the triazole analogue 6 (IC50 against hDHODH: 45 nM) exerted potent in vitro antiproliferative and immunosuppressive activity without affecting cell survival.