Ab initio molecular dynamics on quantum computers.

Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, we present an approach to running AIMD simulations on noisy intermediate-scale quantum (NISQ)-era quantum computers. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. Algorithms for computation of analytical gradients entirely on a quantum computer require quantum fault-tolerant hardware, which is beyond NISQ-era. Therefore, we compute the energy gradients numerically using finite differences, the Hellmann-Feynman theorem, and a correlated sampling technique. This method only requires additional classical calculations of electron integrals for each degree of freedom without any additional computations on a quantum computer beyond the initial VQE run. As a proof of concept, AIMD simulations are demonstrated for the H2 molecule on IBM quantum devices. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular systems.

Strongly Coupled Exciton-Surface Lattice Resonances Engineer Long-Range Energy Propagation.

Achieving propagation lengths in hybrid plasmonic systems beyond typical values of tens of micrometers is important for quantum plasmonics applications. We report long-range optical energy propagation due to excitons in semiconductor quantum dots (SQDs) being strongly coupled to surface lattice resonance (SLRs) in silver nanoparticle arrays. Photoluminescence (PL) measurements provide evidence of an exciton-SLR (ESLR) mode extending at least 600 µm from the excitation region. We also observe additional energy propagation with range well beyond the ESLR mode and with dependency on the coupling strength, g, between SQDs and SLR. Cavity quantum electrodynamics calculations capture the nature of the PL spectra for consistent g values, while coupled dipole calculations show a SQD number-dependent electric field decay profile consistent with the experimental spatial PL profile. Our results suggest an exciting direction wherein SLRs mediate long-range interactions between SQDs, having possible applications in optoelectronics, sensing, and quantum information science.

Non-Hermitian approach for quantum plasmonics.

We examine the limits of applicability of a simple non-Hermitian model for exciton/plasmon interactions in the presence of dissipation and dephasing. The model can be used as an alternative to the more complete Lindblad density matrix approach and is computationally and conceptually simpler. We find that optical spectra in the linear regime can be adequately described by this approach. The model can fail, however, under continuous optical driving in some circumstances. In the case of two quantum dots or excitons interacting with a plasmon, the model can also describe coherences and entanglement qualitatively when both dissipation and dephasing are present and quantitatively in the limit with no dephasing. The approach, within a single excitation manifold, is also applied to assess the role of disorder for 50 quantum dots interacting with a plasmon, where we find that, on average, large enough disorder can help stabilize the ensemble average of the open quantum system toward a dark quasi-steady-state much faster than without disorder. While such single excitation manifold calculations in this size limit can readily be done with either the non-Hermitian or Lindblad forms, as one goes to larger Hilbert space sizes, the computational and storage advantages of the non-Hermitian approach can become more useful.

Excited States of Methylene, Polyenes, and Ozone from Heat-Bath Configuration Interaction.

The electronically excited states of methylene (CH2), ethylene (C2H4), butadiene (C4H6), hexatriene (C6H8), and ozone (O3) have long proven challenging due to their complex mixtures of static and dynamic correlations. The semistochastic heat-bath configuration interaction (SHCI) algorithm, which efficiently and systematically approaches the full configuration interaction (FCI) limit, is used to provide close approximations to the FCI energies in these systems. This article presents the largest FCI-level calculation to date on hexatriene, using a polarized double-ζ basis (ANO-L-pVDZ), which gives rise to a Hilbert space containing more than 1038 determinants. These calculations give vertical excitation energies of 5.58 and 5.59 eV, respectively, for the 21Ag and 11Bu states, showing that they are nearly degenerate. The same excitation energies in butadiene/ANO-L-pVDZ were found to be 6.58 and 6.45 eV. In addition to these benchmarks, our calculations strongly support the presence of a previously hypothesized ring-minimum species of ozone that lies 1.3 eV higher than the open-ring-minimum energy structure and is separated from it by a barrier of 1.11 eV.

Fast semistochastic heat-bath configuration interaction.

This paper presents in detail our fast semistochastic heat-bath configuration interaction (SHCI) method for solving the many-body Schrödinger equation. We identify and eliminate computational bottlenecks in both the variational and perturbative steps of the SHCI algorithm. We also describe the parallelization and the key data structures in our implementation, such as the distributed hash table. The improved SHCI algorithm enables us to include in our variational wavefunction two orders of magnitude more determinants than has been reported previously with other selected configuration interaction methods. We use our algorithm to calculate an accurate benchmark energy for the chromium dimer with the X2C relativistic Hamiltonian in the cc-pVDZ-DK basis, correlating 28 electrons in 76 spatial orbitals. Our largest calculation uses two billion Slater determinants in the variational space and semistochastically includes perturbative contributions from at least trillions of additional determinants with better than 10-5 Ha statistical uncertainty.

State Preparation in Quantum Algorithms for Fragment-Based Quantum Chemistry.

State preparation for quantum algorithms is crucial for achieving high accuracy in quantum chemistry and competing with classical algorithms. The localized active space-unitary coupled cluster (LAS-UCC) algorithm iteratively loads a fragment-based multireference wave function onto a quantum computer. In this study, we compare two state preparation methods, quantum phase estimation (QPE) and direct initialization (DI), for each fragment. We test the two state preparation methods on three systems, ranging from a model system, a set of interacting hydrogen molecules, to more realistic chemical problems, like the C-C double bond breaking in transbutadiene and the spin ladder in a bimetallic system. We analyze the impact of QPE parameters, such as the number of ancilla qubits and Trotter steps, on the prepared state. We find a trade-off between the methods, where DI requires fewer resources for smaller fragments, while QPE is more efficient for larger fragments. Our resource estimates highlight the benefits of system fragmentation in state preparation for subsequent quantum chemical calculations. These findings have broad applications for preparing multireference quantum chemical wave functions on quantum circuits that can be used for realistic chemical applications.

The learnability of Pauli noise.

Recently, several quantum benchmarking algorithms have been developed to characterize noisy quantum gates on today's quantum devices. A fundamental issue in benchmarking is that not everything about quantum noise is learnable due to the existence of gauge freedom, leaving open the question what information is learnable and what is not, which is unclear even for a single CNOT gate. Here we give a precise characterization of the learnability of Pauli noise channels attached to Clifford gates using graph theoretical tools. Our results reveal the optimality of cycle benchmarking in the sense that it can extract all learnable information about Pauli noise. We experimentally demonstrate noise characterization of IBM's CNOT gate up to 2 unlearnable degrees of freedom, for which we obtain bounds using physical constraints. In addition, we show that an attempt to extract unlearnable information by ignoring state preparation noise yields unphysical estimates, which is used to lower bound the state preparation noise.

Localized Quantum Chemistry on Quantum Computers.

Quantum chemistry calculations of large, strongly correlated systems are typically limited by the computation cost that scales exponentially with the size of the system. Quantum algorithms, designed specifically for quantum computers, can alleviate this, but the resources required are still too large for today's quantum devices. Here, we present a quantum algorithm that combines a localization of multireference wave functions of chemical systems with quantum phase estimation (QPE) and variational unitary coupled cluster singles and doubles (UCCSD) to compute their ground-state energy. Our algorithm, termed "local active space unitary coupled cluster" (LAS-UCC), scales linearly with the system size for certain geometries, providing a polynomial reduction in the total number of gates compared with QPE, while providing accuracy above that of the variational quantum eigensolver using the UCCSD ansatz and also above that of the classical local active space self-consistent field. The accuracy of LAS-UCC is demonstrated by dissociating (H2)2 into two H2 molecules and by breaking the two double bonds in trans-butadiene, and resource estimates are provided for linear chains of up to 20 H2 molecules.

SNP in Potentially Defunct Tetrahydrocannabinolic Acid Synthase Is a Marker for Cannabigerolic Acid Dominance in Cannabis sativa L.

The regulation of cannabinoid synthesis in Cannabis sativa is of increasing research interest as restrictions around the globe loosen to allow the plant's legal cultivation. Of the major cannabinoids, the regulation of cannabigerolic acid (CBGA) production is the least understood. The purpose of this study was to elucidate the inheritance of CBGA dominance in C. sativa and describe a marker related to this chemotype. We produced two crossing populations, one between a CBGA dominant cultivar and a tetrahydrocannabinolic acid (THCA) dominant cultivar, and one between a CBGA dominant cultivar and a cannabidiolic acid (CBDA) cultivar. Chemical and genotyping analyses confirmed that CBGA dominance is inherited as a single recessive gene, potentially governed by a non-functioning allelic variant of the THCA synthase. The "null" THCAS synthase contains a single nucleotide polymorphism (SNP) that may render the synthase unable to convert CBGA to THCA leading to the accumulation of CBGA. This SNP can be reliably used as a molecular marker for CBGA dominance in the selection and breeding of C. sativa.

Characteristics of the Diploid, Triploid, and Tetraploid Versions of a Cannabigerol-Dominant F1 Hybrid Industrial Hemp Cultivar, Cannabis sativa 'Stem Cell CBG'.

Hemp (Cannabis sativa L.) has recently become an important crop due to the growing market demands for products containing cannabinoids. Unintended cross-pollination of C. sativa crops is one of the most important threats to cannabinoid production and has been shown to reduce cannabinoid yield. Ploidy manipulation has been used in other crops to improve agronomic traits and reduce fertility; however, little is known about the performance of C. sativa polyploids. In this study, colchicine was applied to two proprietary, inbred diploid C. sativa inbred lines, 'TS1-3' and 'P163', to produce the tetraploids 'TS1-3 (4x)' and 'P163 (4x)'. The diploid, triploid, and tetraploid F1 hybrids from 'TS1-3' × 'P163', 'TS1-3 (4x)' × 'P163', and 'TS1-3 (4x)' × 'P163 (4x)' were produced to test their fertilities, crossing compatibilities, and yields. The results indicated a reduction in fertility in the triploids and the tetraploids, relative to their diploid counterparts. When triploids were used as females, seed yields were less than 2% compared to when diploids were used as females; thus, triploids were determined to be female infertile. The triploids resulting from the crosses made herein displayed increases in biomass and inflorescence weight compared to the diploids created from the same parents in a field setting. Statistical increases in cannabinoid concentrations were not observed. Lastly, asymmetric crossing compatibility was observed between the diploids and the tetraploids of the genotypes tested. The results demonstrate the potential benefits of triploid C. sativa cultivars in commercial agriculture.

Uniaxial transition dipole moments in semiconductor quantum rings caused by broken rotational symmetry.

Semiconductor quantum rings are topological structures that support fascinating phenomena such as the Aharonov-Bohm effect and persistent current, which are of high relevance in the research of quantum information devices. The annular shape of quantum rings distinguishes them from other low-dimensional materials, and enables topologically induced properties such as geometry-dependent spin manipulation and emission. While optical transition dipole moments (TDMs) in zero to two-dimensional optical emitters have been well investigated, those in quantum rings remain obscure despite their utmost relevance to the quantum photonic applications of quantum rings. Here, we study the dimensionality and orientation of TDMs in CdSe quantum rings. In contrast to those in other two-dimensional optical emitters, we find that TDMs in CdSe quantum rings show a peculiar in-plane linear distribution. Our theoretical modeling reveals that this uniaxial TDM originates from broken rotational symmetry in the quantum ring geometries.

Pretreatment with intragraft verapamil prior to percutaneous coronary intervention of saphenous vein graft lesions: results of the randomized, controlled vasodilator prevention on no-reflow (VAPOR) trial.

BACKGROUND: Intragraft verapamil is effective in treating no-reflow during saphenous vein graft (SVG) percutaneous coronary intervention (PCI). In this study, we assessed the use of intragraft verapamil given pre-PCI to prevent no-reflow. METHODS: Patients undergoing SVG PCI were randomized to receive intragraft 200 g verapamil or no verapamil immediately prior to PCI. Pre- and post-PCI, vessel flow was assessed using TIMI flow grade and TIMI frame count by blinded angiographic readers. Tissue level perfusion in the graft territory was assessed using the TIMI myocardial perfusion grade (TMPG). CK-MB or troponin I levels were measured 6 12 hours post-PCI. RESULTS: Ten patients were randomized to the verapamil group and 12 were assigned to the placebo group. No-reflow occurred in 33.3% of the placebo group, compared to none of the verapamil patients (p = 0.10). The use of intragraft verapamil prior to SVG PCI increased flow rate in the vessel as assessed by TIMI frame count (53.3 22.4% faster in the verapamil group versus 11.5 38.9% in the placebo group; p = 0.016). There was a trend toward improved myocardial perfusion as assessed by TMPG. There was no difference in the incidence of cardiac biomarker release following PCI. CONCLUSIONS: Intragraft administration of verapamil prior to saphenous vein graft PCI reduces no-reflow and is associated with a trend toward improved myocardial perfusion.