Efficient Construction of Involutory Linear Combinations of Anticommuting Pauli Generators for Large-Scale Iterative Qubit Coupled Cluster Calculations.

We present an efficient method for construction of a fully anticommutative set of Pauli generators (elements of the Pauli group) from a commutative set of operators that are composed exclusively from Pauli x̂ operators (purely X generators) and sorted by an associated numerical measure, such as absolute energy gradients. Our approach uses the Gauss-Jordan elimination applied to a binary matrix that encodes the set of X generators to bring it to the reduced row-echelon form, followed by the construction of an anticommutative system in a standard basis by means of a modified Jordan-Wigner transformation and returning to the original basis. The algorithm complexity is linear in the size of the X set and quadratic in the number of qubits. The resulting anticommutative sets are used to construct the qubit coupled cluster Ansatz with involutory linear combinations of anticommuting Paulis (QCC-ILCAP) proposed in J. Chem. Theory Comput. 2021, 17 (1), 66-78. We applied the iterative qubit coupled cluster method with the QCC-ILCAP Ansatz to calculations of ground-state potential energy curves for symmetric stretching of the water molecule (36 qubits) and dissociation of N2 (56 qubits).

On the Importance of Well-Defined Thermal Correlation Functions in Simulating Vibronic Spectra.

Two difficulties associated with the computations of thermal vibrational correlation functions are discussed. The first one is the lack of a well-behaved expression that is valid at both high-temperature and T → 0 K limits. Specifically, if the partition function and the propagator are considered separately, then thermal vibrational correlation functions may have an indeterminate form 0/0 in the limit T → 0 K. This difficulty is resolved when the partition function and the propagator are jointly considered in the harmonic approximation, which allows a problematic term that emanates from the zero-point energy to be canceled out, thereby producing a thermal correlation function with a determinate form in the T → 0 K limit. The second difficulty is related to the multivaluedness of the vibrational correlation function. We show numerically that an improper selection of branch leads to discontinuities in the computed correlation function and incorrect vibronic spectra. We propose a phase tracking procedure that ensures continuity of both real and imaginary parts of the correlation function to recover the correct spectra. We support our findings by simulating the UV-vis absorption spectra of pentacene at 4 K and benzene at 298 K. Both are found to be in good agreement with their experimental counterparts.

Estimating Phosphorescent Emission Energies in IrIII Complexes Using Large-Scale Quantum Computing Simulations.

Here we calculate T1 →S0 transition energies in nine phosphorescent iridium complexes using the iterative qubit coupled cluster (iQCC) method to determine if quantum simulations have any advantages over classical methods. These simulations would require a gate-based quantum computer with at least 72 fully-connected logical qubits. Since such devices do not yet exist, we demonstrate the iQCC method using a purpose-built quantum simulator on classical hardware. The results are compared to a selection of common DFT functionals, ab initio methods, and empirical data. iQCC is found to match the accuracy of the best DFT functionals, but with a better correlation coefficient, demonstrating that it is better at predicting the structure-property relationship. Results indicate that the iQCC method has the required accuracy to design organometallic complexes when deployed on emerging quantum hardware and sets an industrially relevant target for demonstrating quantum advantage.

Unitary Transformation of the Electronic Hamiltonian with an Exact Quadratic Truncation of the Baker-Campbell-Hausdorff Expansion.

The application of current and near-term quantum hardware to the electronic structure problem is highly limited by qubit counts, coherence times, and gate fidelities. To address these restrictions within the variational quantum eigensolver (VQE) framework, many recent contributions have suggested dressing the electronic Hamiltonian to include a part of electron correlation, leaving the rest to VQE state preparation. We present a new dressing scheme that combines the preservation of the Hamiltonian hermiticity and an exact quadratic truncation of the Baker-Campbell-Hausdorff expansion. The new transformation is constructed as the exponent of an involutory linear combination (ILC) of anti-commuting Pauli products. It incorporates important strong correlation effects in the dressed Hamiltonian and can be viewed as a classical preprocessing step to alleviate the resource requirements of the subsequent VQE application. The assessment of the new computational scheme for the electronic structure of the LiH, H2O, and N2 molecules shows a significant increase in efficiency compared to the conventional qubit coupled cluster dressings.

Iterative Qubit Coupled Cluster Approach with Efficient Screening of Generators.

An iterative version of the qubit coupled cluster (QCC) method [I. G. Ryabinkin et al., J. Chem. Theory Comput. 2019, 14, 6317] is proposed. The new method seeks to find ground electronic energies of molecules on noisy intermediate-scale quantum devices. Each iteration involves a canonical transformation of the Hamiltonian and employs constant-size quantum circuits at the expense of increasing the Hamiltonian size. We numerically studied the convergence of the method on ground-state calculations for LiH, H2O, and N2 molecules and found that the exact ground-state energies can be systematically approached only if the generators of the QCC ansatz are sampled from a specific set of operators. We report an algorithm for constructing this set that scales linearly with the size of the Hamiltonian.

Revising the measurement process in the variational quantum eigensolver: is it possible to reduce the number of separately measured operators?

Current implementations of the Variational Quantum Eigensolver (VQE) technique for solving the electronic structure problem involve splitting the system qubit Hamiltonian into parts whose elements commute within their single qubit subspaces. The number of such parts rapidly grows with the size of the molecule. This increases the computational cost and can increase uncertainty in the measurement of the energy expectation value because elements from different parts need to be measured independently. To address this problem we introduce a more efficient partitioning of the qubit Hamiltonian using fewer parts that need to be measured separately. The new partitioning scheme is based on two ideas: (1) grouping terms into parts whose eigenstates have a single-qubit product structure, and (2) devising multi-qubit unitary transformations for the Hamiltonian or its parts to produce less entangled operators. The first condition allows the new parts to be measured in the number of involved qubit consequential one-particle measurements. Advantages of the new partitioning scheme resulting in severalfold reduction of separately measured terms are illustrated with regard to the H2 and LiH problems.

Constrained Variational Quantum Eigensolver: Quantum Computer Search Engine in the Fock Space.

Variational quantum eigensolver (VQE) is an efficient computational method promising chemical accuracy in electronic structure calculations on a universal-gate quantum computer. However, such a simple task as computing the electronic energy of a hydrogen molecular cation, H2+, is not possible for a general VQE protocol because the calculation will invariably collapse to a lower energy of the corresponding neutral form, H2. The origin of the problem is that VQE effectively performs an unconstrained energy optimization in the Fock space of the original electronic problem. We show how this can be avoided by introducing necessary constraints directing VQE toward the electronic state of interest. The proposed constrained VQE can find an electronic state with a certain number of electrons, a certain spin, or any other property. Moreover, the new algorithm naturally removes unphysical kinks in potential energy surfaces (PESs), which frequently appeared in the regular VQE and required significant additional quantum resources for their removal. We demonstrate the performance of the constrained VQE by simulating PESs of various states of H2 and H2O on Rigetti Computing Inc.'s 19Q-Acorn quantum processor.

Relation between fermionic and qubit mean fields in the electronic structure problem.

For quantum computing applications, the electronic Hamiltonian for the electronic structure problem needs to be unitarily transformed into a qubit form. We found that mean-field procedures on the original electronic Hamiltonian and on its transformed qubit counterpart can give different results. We establish conditions of when fermionic and qubit mean fields provide the same or different energies. In cases when the fermionic mean-field (Hartree-Fock) approach provides an accurate description (electronic correlation effects are small), the choice of molecular orbitals for the electron Hamiltonian representation becomes the determining factor in whether the qubit mean-field energy will be equal to or higher than that of the fermionic counterpart. In strongly correlated cases, the qubit mean-field approach has a higher chance to undergo symmetry breaking and lower its energy below the fermionic counterpart.

Qubit Coupled Cluster Method: A Systematic Approach to Quantum Chemistry on a Quantum Computer.

A unitary coupled cluster (UCC) form for the wave function in the variational quantum eigensolver has been suggested as a systematic way to go beyond the mean-field approximation and include electron correlation in solving quantum chemistry problems on a quantum computer. Although being exact in the limit of including all possible coupled cluster excitations, practically, the accuracy of this approach depends on the number and type of terms are included in the wave function parametrization. Another difficulty of UCC is a growth of the number of simultaneously entangled qubits even at the fixed Fermionic excitation rank. Not all quantum computing architectures can cope with this growth. To address both problems, we introduce a qubit coupled cluster (QCC) method that starts directly in the qubit space and uses energy response estimates for ranking the importance of individual entanglers for the variational energy minimization. Also, we provide an exact factorization of a unitary rotation of more than two qubits to a product of two-qubit unitary rotations. Thus, the QCC method with the factorization technique can be limited to only two-qubit entanglement gates and allows for very efficient use of quantum resources in terms of the number of coupled cluster operators. The method performance is illustrated by calculating ground-state potential energy curves of H2 and LiH molecules with chemical accuracy, ≤1 kcal/mol, and a symmetric water dissociation curve.

Exact exchange-correlation potentials of singlet two-electron systems.

We suggest a non-iterative analytic method for constructing the exchange-correlation potential, vXC(r), of any singlet ground-state two-electron system. The method is based on a convenient formula for vXC(r) in terms of quantities determined only by the system's electronic wave function, exact or approximate, and is essentially different from the Kohn-Sham inversion technique. When applied to Gaussian-basis-set wave functions, the method yields finite-basis-set approximations to the corresponding basis-set-limit vXC(r), whereas the Kohn-Sham inversion produces physically inappropriate (oscillatory and divergent) potentials. The effectiveness of the procedure is demonstrated by computing accurate exchange-correlation potentials of several two-electron systems (helium isoelectronic series, H2, H3+) using common ab initio methods and Gaussian basis sets.

Geometric Phase Effects in Nonadiabatic Dynamics near Conical Intersections.

Dynamical consideration that goes beyond the common Born-Oppenheimer approximation (BOA) becomes necessary when energy differences between electronic potential energy surfaces become small or vanish. One of the typical scenarios of the BOA breakdown in molecules beyond diatomics is a conical intersection (CI) of electronic potential energy surfaces. CIs provide an efficient mechanism for radiationless electronic transitions: acting as "funnels" for the nuclear wave function, they enable rapid conversion of the excessive electronic energy into the nuclear motion. In addition, CIs introduce nontrivial geometric phases (GPs) for both electronic and nuclear wave functions. These phases manifest themselves in change of the wave function signs if one considers an evolution of the system around the CI. This sign change is independent of the shape of the encircling contour and thus has a topological character. How these extra phases affect nonadiabatic dynamics is the main question that is addressed in this Account. We start by considering the simplest model providing the CI topology: two-dimensional two-state linear vibronic coupling model. Selecting this model instead of a real molecule has the advantage that various dynamical regimes can be easily modeled in the model by varying parameters, whereas any fixed molecule provides the system specific behavior that may not be very illustrative. After demonstrating when GP effects are important and how they modify the dynamics for two sets of initial conditions (starting from the ground and excited electronic states), we give examples of molecular systems where the described GP effects are crucial for adequate description of nonadiabatic dynamics. Interestingly, although the GP has a topological character, the extent to which accounting for GPs affect nuclear dynamics profoundly depends on topography of potential energy surfaces. Understanding an extent of changes introduced by the GP in chemical dynamics poses a problem of capturing GP effects by approximate methods of simulating nonadiabatic dynamics that can go beyond simple models. We assess the performance of both fully quantum (wave packet dynamics) and quantum-classical (surface-hopping, Ehrenfest, and quantum-classical Liouville equation) approaches in various cases where GP effects are important. It has been identified that the key to success in approximate methods is a method organization that prevents the quantum nuclear kinetic energy operator to act directly on adiabatic electronic wave functions.

Improved method for generating exchange-correlation potentials from electronic wave functions.

Ryabinkin, Kohut, and Staroverov (RKS) [Phys. Rev. Lett. 115, 083001 (2015)] devised an iterative method for reducing many-electron wave functions to Kohn-Sham exchange-correlation potentials, vXC(𝐫). For a given type of wave function, the RKS method is exact (Kohn-Sham-compliant) in the basis-set limit; in a finite basis set, it produces an approximation to the corresponding basis-set-limit vXC(𝐫). The original RKS procedure works very well for large basis sets but sometimes fails for commonly used (small and medium) sets. We derive a modification of the method's working equation that makes the RKS procedure robust for all Gaussian basis sets and increases the accuracy of the resulting exchange-correlation potentials with respect to the basis-set limit.

Topologically Correct Quantum Nonadiabatic Formalism for On-the-Fly Dynamics.

On-the-fly quantum nonadiabatic dynamics for large systems greatly benefits from the adiabatic representation readily available from electronic structure programs. However, conical intersections frequently occurring in this representation introduce nontrivial geometric or Berry phases which require a special treatment for adequate modeling of the nuclear dynamics. We analyze two approaches for nonadiabatic dynamics using the time-dependent variational principle and the adiabatic representation. The first approach employs adiabatic electronic functions with global parametric dependence on the nuclear coordinates. The second approach uses adiabatic electronic functions obtained only at the centers of moving localized nuclear basis functions (e.g., frozen-width Gaussians). Unless a gauge transformation is used to enforce single-valued boundary conditions, the first approach fails to capture the geometric phase. In contrast, the second approach accounts for the geometric phase naturally because of the absence of the global nuclear coordinate dependence in the electronic functions.

Mixed Quantum-Classical Dynamics Using Collective Electronic Variables: A Better Alternative to Electronic Friction Theories.

An accurate description of nonadiabatic dynamics of molecular species on metallic surfaces poses a serious computational challenge associated with a multitude of closely spaced electronic states. We propose a mixed quantum-classical scheme that addresses this challenge by introducing collective electronic variables. These variables are defined through analytic block-diagonalization applied to the time-dependent Hamiltonian matrix governing the electronic dynamics. We compare our scheme with a simplified Ehrenfest approach and with a full-memory electronic friction model on a 1D "adatom + atomic chain" model. Our simulations demonstrate that collective-mode dynamics with only a few (two to three) electronic variables is robust and can describe a variety of situations: from a chemisorbed atom on an insulator to an atom on a metallic surface. Our molecular model also reveals that the friction approach is prone to unpredictable and catastrophic failures.

On the inclusion of the diagonal Born-Oppenheimer correction in surface hopping methods.

The diagonal Born-Oppenheimer correction (DBOC) stems from the diagonal second derivative coupling term in the adiabatic representation, and it can have an arbitrary large magnitude when a gap between neighbouring Born-Oppenheimer (BO) potential energy surfaces (PESs) is closing. Nevertheless, DBOC is typically neglected in mixed quantum-classical methods of simulating nonadiabaticdynamics (e.g., fewest-switch surface hopping (FSSH) method). A straightforward addition of DBOC to BO PESs in the FSSH method, FSSH+D, has been shown to lead to numerically much inferior results for models containing conical intersections. More sophisticated variation of the DBOC inclusion, phase-space surface-hopping (PSSH) was more successful than FSSH+D but on model problems without conical intersections. This work comprehensively assesses the role of DBOC in nonadiabaticdynamics of two electronic state problems and the performance of FSSH, FSSH+D, and PSSH methods in variety of one- and two-dimensional models. Our results show that the inclusion of DBOC can enhance the accuracy of surface hopping simulations when two conditions are simultaneously satisfied: (1) nuclei have kinetic energy lower than DBOC and (2) PESs are not strongly nonadiabatically coupled. The inclusion of DBOC is detrimental in situations where its energy scale becomes very high or even diverges, because in these regions PESs are also very strongly coupled. In this case, the true quantum formalism heavily relies on an interplay between diagonal and off-diagonal nonadiabatic couplings while surface hopping approaches treat diagonal terms as PESs and off-diagonal ones stochastically.

Shape-Dependent Interactions of Palladium Nanocrystals with Hydrogen.

Elucidation of the nature of hydrogen interactions with palladium nanoparticles is expected to play an important role in the development of new catalysts and hydrogen-storage nanomaterials. A facile scaled-up synthesis of uniformly sized single-crystalline palladium nanoparticles with various shapes, including regular nanocubes, nanocubes with protruded edges, rhombic dodecahedra, and branched nanoparticles, all stabilized with a mesoporous silica shell is developed. Interaction of hydrogen with these nanoparticles is studied by using temperature-programmed desorption technique and by performing density functional theory modeling. It is found that due to favorable arrangement of Pd atoms on their surface, rhombic dodecahedral palladium nanoparticles enclosed by {110} planes release a larger volume of hydrogen and have a lower desorption energy than palladium nanocubes and branched nanoparticles. These results underline the important role of {110} surfaces in palladium nanoparticles in their interaction with hydrogen. This work provides insight into the mechanism of catalysis of hydrogenation/dehydrogenation reactions by palladium nanoparticles with different shapes.

Why Do Mixed Quantum-Classical Methods Describe Short-Time Dynamics through Conical Intersections So Well? Analysis of Geometric Phase Effects.

Adequate simulation of nonadiabatic dynamics through conical intersection requires accounting for a nontrivial geometric phase (GP) emerging in electronic and nuclear wave functions in the adiabatic representation. Popular mixed quantum-classical (MQC) methods, surface hopping and Ehrenfest, do not carry a nuclear wave function to be able to incorporate the GP into nuclear dynamics. Surprisingly, the MQC methods reproduce ultrafast interstate crossing dynamics generated with the exact quantum propagation so well as if they contained information about the GP. Using two-dimensional linear vibronic coupling models we unravel how the MQC methods can effectively mimic the most significant dynamical GP effects: (1) compensation for repulsive diagonal second-order nonadiabatic couplings and (2) transfer enhancement for a fully cylindrically symmetric component of a nuclear distribution.

Fast Numerical Evaluation of Time-Derivative Nonadiabatic Couplings for Mixed Quantum-Classical Methods.

We have developed a numerical differentiation scheme that eliminates evaluation of overlap determinants in calculating the time-derivative nonadiabatic couplings (TDNACs). Evaluation of these determinants was the bottleneck in previous implementations of mixed quantum-classical methods using numerical differentiation of electronic wave functions in the Slater determinant representation. The central idea of our approach is, first, to reduce the analytic time derivatives of Slater determinants to time derivatives of molecular orbitals and then to apply a finite-difference formula. Benchmark calculations prove the efficiency of the proposed scheme showing impressive several-order-of-magnitude speedups of the TDNAC calculation step for midsize molecules.