Extension of Exactly-Solvable Hamiltonians Using Symmetries of Lie Algebras.

Exactly solvable Hamiltonians that can be diagonalized by using relatively simple unitary transformations are of great use in quantum computing. They can be employed for the decomposition of interacting Hamiltonians either in Trotter-Suzuki approximations of the evolution operator for the quantum phase estimation algorithm or in the quantum measurement problem for the variational quantum eigensolver. One of the typical forms of exactly solvable Hamiltonians is a linear combination of operators forming a modestly sized Lie algebra. Very frequently, such linear combinations represent noninteracting Hamiltonians and thus are of limited interest for describing interacting cases. Here, we propose an extension in which the coefficients in these combinations are substituted by polynomials of the Lie algebra symmetries. This substitution results in a more general class of solvable Hamiltonians, and for qubit algebras, it is related to the recently proposed noncontextual Pauli Hamiltonians. In fermionic problems, this substitution leads to Hamiltonians with eigenstates that are single Slater determinants but with different sets of single-particle states for different eigenstates. The new class of solvable Hamiltonians can be measured efficiently using quantum circuits with gates that depend on the result of a midcircuit measurement of the symmetries.

Exactly solvable Hamiltonian fragments obtained from a direct sum of Lie algebras.

Exactly solvable Hamiltonians are useful in the study of quantum many-body systems using quantum computers. In the variational quantum eigensolver, a decomposition of the target Hamiltonian into exactly solvable fragments can be used for the evaluation of the energies via repeated quantum measurements. In this work, we apply more general classes of exactly solvable qubit Hamiltonians than previously considered to address the Hamiltonian measurement problem. The most general exactly solvable Hamiltonians we use are defined by the condition that within each simultaneous eigenspace of a set of Pauli symmetries, the Hamiltonian acts effectively as an element of a direct sum of so(N) Lie algebras and can, therefore, be measured using a combination of unitaries in the associated Lie group, Clifford unitaries, and mid-circuit measurements. The application of such Hamiltonians to decomposing molecular electronic Hamiltonians via graph partitioning techniques shows a reduction in the total number of measurements required to estimate the expectation value compared to previously used exactly solvable qubit Hamiltonians.

A quantum computing view on unitary coupled cluster theory.

We present a review of the Unitary Coupled Cluster (UCC) ansatz and related ansätze which are used to variationally solve the electronic structure problem on quantum computers. A brief history of coupled cluster (CC) methods is provided, followed by a broad discussion of the formulation of CC theory. This includes touching on the merits and difficulties of the method and several variants, UCC among them, in the classical context, to motivate their applications on quantum computers. In the core of the text, the UCC ansatz and its implementation on a quantum computer are discussed at length, in addition to a discussion on several derived and related ansätze specific to quantum computing. The review concludes with a unified perspective on the discussed ansätze, attempting to bring them under a common framework, as well as with a reflection upon open problems within the field.

Fidelity Overhead for Nonlocal Measurements in Variational Quantum Algorithms.

Measuring quantum observables by grouping terms that can be rotated to sums of only products of Pauli z operators (Ising form) is proven to be efficient in near term quantum computing algorithms. This approach requires extra unitary transformations to rotate the state of interest so that the measurement of a fragment's Ising form would be equivalent to the measurement of the fragment for the unrotated state. These extra rotations allow one to perform a fewer number of measurements by grouping more terms into the measurable fragments with a lower overall estimator variance. However, previous estimations of the number of measurements did not take into account nonunit fidelity of quantum gates implementing the additional transformations. Through a circuit fidelity reduction, additional transformations introduce extra uncertainty and increase the needed number of measurements. Here we consider a simple model for errors introduced by additional gates needed in schemes involving groupings of commuting Pauli products. For a set of molecular electronic Hamiltonians, we confirm that the numbers of measurements in schemes using nonlocal qubit rotations are still lower than those in their local qubit rotation counterparts, even after accounting for uncertainties introduced by additional gates.

Controlling energy conservation in quantum dynamics with independently moving basis functions: Application to multi-configuration Ehrenfest.

Application of the time-dependent variational principle to a linear combination of frozen-width Gaussians describing the nuclear wavefunction provides a formalism where the total energy is conserved. The computational downside of this formalism is that trajectories of individual Gaussians are solutions of a coupled system of differential equations, limiting implementation to serial propagation algorithms. To allow for parallelization and acceleration of the computation, independent trajectories based on simplified equations of motion were suggested. Unfortunately, within practical realizations involving finite Gaussian bases, this simplification leads to breaking the energy conservation. We offer a solution for this problem by using Lagrange multipliers to ensure the energy and norm conservation regardless of basis function trajectories or basis completeness. We illustrate our approach within the multi-configurational Ehrenfest method considering a linear vibronic coupling model.

Computational approaches to efficient generation of the stationary state for incoherent light excitation.

Light harvesting processes are often computationally studied from a time-dependent viewpoint, in line with ultrafast coherent spectroscopy experiments. Yet, natural processes take place in the presence of incoherent light, which induces a stationary state. Such stationary states can be described using the eigenbasis of the molecular Hamiltonian, but for realistic systems, a full diagonalization is prohibitively expensive. We propose three efficient computational approaches to obtain the stationary state that circumvents system Hamiltonian diagonalization. The connection between the incoherent perturbations, decoherence, and Kraus operators is established.

On the order problem in construction of unitary operators for the variational quantum eigensolver.

One of the main challenges in the variational quantum eigensolver (VQE) framework is construction of the unitary transformation. The dimensionality of the space for unitary rotations of N qubits is 4N- 1, which makes the choice of a polynomial subset of generators an exponentially difficult process. Moreover, due to non-commutativity of generators, the order in which they are used strongly affects results. Choosing the optimal order in a particular subset of generators requires testing the factorial number of combinations. We propose an approach based on the Lie algebra-Lie group connection and corresponding closure relations that systematically eliminates the order problem.

Measurement optimization in the variational quantum eigensolver using a minimum clique cover.

Solving the electronic structure problem using the Variational Quantum Eigensolver (VQE) technique involves the measurement of the Hamiltonian expectation value. The current hardware can perform only projective single-qubit measurements, and thus, the Hamiltonian expectation value is obtained by measuring parts of the Hamiltonian rather than the full Hamiltonian. This restriction makes the measurement process inefficient because the number of terms in the Hamiltonian grows as O(N4) with the size of the system, N. To optimize the VQE measurement, one can try to group as many Hamiltonian terms as possible for their simultaneous measurement. Single-qubit measurements allow one to group only the terms commuting within the corresponding single-qubit subspaces or qubit-wise commuting. We found that the qubit-wise commutativity between the Hamiltonian terms can be expressed as a graph and the problem of the optimal grouping is equivalent to finding a minimum clique cover (MCC) for the Hamiltonian graph. The MCC problem is NP-hard, but there exist several polynomial heuristic algorithms to solve it approximately. Several of these heuristics were tested in this work for a set of molecular electronic Hamiltonians. On average, grouping qubit-wise commuting terms reduced the number of operators to measure three times less compared to the total number of terms in the considered Hamiltonians.

Deuterium isotope effect in fluorescence of gaseous oxazine dyes.

The increased utility of fluorescence-based methods in recent years has highlighted the need for brighter, more efficient fluorophores. In order to design these fluorophores, an improved fundamental understanding is necessary of the structural components that intrinsically effect fluorescence efficiency. Here, we characterize the intrinsic effects of deuteration on fluorescence from gaseous oxazine dyes, without the influence of dye-solvent interactions, by making use of an ion trap mass spectrometer that has been altered to enable optical measurements. Comparison of emission spectra of four oxazine dyes: cresyl violet, oxazine 4, oxazine 170, and darrow red, show little change in profile upon deuteration of amine groups. However, deuteration significantly increases the efficiency of fluorescence with an increase in fluorescence lifetime and brightness by 10-23% for the gaseous dyes. This increase is less than half that of the quantum yield increase observed in deuterated solution. This indicates the large fluorescence efficiency changes for the oxazine dyes in deuterated solution result from a combination of both intrinsic effects as well as substantial contribution from altered fluorophore-solvent interactions. The intrinsic effects behind increased lifetime upon deuteration are explored using time-dependent density functional theory (TD-DFT) calculations of potential energy surfaces (PESs) for ground and low lying excited electronic states. In accord with experimental observations, calculated S1-S0 emission spectra show only minor differences between deuterated and non-deuterated forms indicating that the deuteration does not affect the radiative channel appreciably. Relaxed PES scans along the torsional motions of the amino groups reveal that the increase in lifetimes upon deuteration is likely due to quenching of different radiationless changes channels in different oxazine dyes. Calculations suggest that tunneling to access twisted intramolecular charge transfer states in S1 is critical in several of the oxazines. However, in at least one of the dyes examined, the large isotope effect is more likely due to differences in intersystem crossing rates. Overall, this combined experimental and computational investigation elucidates the photophysics of a well-known fluorescent scaffold and provides insight into how small differences can dramatically affect fluorescence outcomes.

On Construction of Projection Operators.

The problem of construction of projection operators on eigensubspaces of symmetry operators is considered. This problem arises in many approximate methods for solving time-independent and time-dependent quantum problems, and its solution ensures proper physical symmetries in development of approximate methods. The projector form is sought as a function of symmetry operators and their eigenvalues characterizing the eigensubspace of interest. This form is obtained in two steps: (1) identification of algebraic structures within a set of symmetry operators (e.g., groups and Lie algebras) and (2) construction of the projection operators for individual symmetry operators. The first step is crucial for efficient projection operator construction because it allows for use of information on irreducible representations of the present algebraic structure. The discussed approaches have potential to stimulate further developments of variational approaches for electronic structure of strongly correlated systems and in quantum computing.

Exact and approximate symmetry projectors for the electronic structure problem on a quantum computer.

Solving the electronic structure problem on a universal-gate quantum computer within the variational quantum eigensolver (VQE) methodology requires constraining the search procedure to a subspace defined by relevant physical symmetries. Ignoring symmetries results in convergence to the lowest eigenstate of the Fock space for the second quantized electronic Hamiltonian. Moreover, this eigenstate can be symmetry broken due to limitations of the wavefunction ansatz. To address this VQE problem, we introduce and assess methods of exact and approximate projectors to irreducible eigensubspaces of available physical symmetries. Feasibility of symmetry projectors in the VQE framework is discussed, and their efficiency is compared with symmetry constraint optimization procedures. Generally, projectors introduce a higher number of terms for VQE measurement compared to the constraint approach. On the other hand, the projection formalism improves accuracy of the variational wavefunction ansatz without introducing additional unitary transformations, which is beneficial for reducing depths of quantum circuits.

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.

Nonadiabatic Quantum Dynamics with Frozen-Width Gaussians.

We review techniques for simulating fully quantum nonadiabatic dynamics using the frozen-width moving Gaussian basis functions to represent the nuclear wave function. A choice of these basis functions is primarily motivated by the idea of the on-the-fly dynamics that will involve electronic structure calculations done locally in the vicinity of each Gaussian center and thus avoiding the "curse of dimensionality" appearing in large systems. For quantum dynamics involving multiple electronic states there are several aspects that need to be addressed. First, the choice of the electronic-state representation is one of most defining in terms of formulation of resulting equations of motion. We will discuss pros and cons of the standard adiabatic and diabatic representations as well as the relatively new moving crude adiabatic (MCA) representation. Second, if the number of electronic states can be fixed throughout the dynamics, the situation is different for the number of Gaussians needed for an accurate expansion of the total wave function. The latter increases its complexity along the course of the dynamics and a protocol extending the number of Gaussians is needed. We will consider two common approaches for the extension: (1) spawning and (2) cloning. Third, equations of motion for individual Gaussians can be chosen in different ways, implications for the energy conservation related to these ways will be discussed. Finally, to extend the success of moving basis approaches to quantum dynamics of open systems we will consider the Nonstochastic Open System Schrödinger Equation (NOSSE).

On the breakdown of the Ehrenfest method for molecular dynamics on surfaces.

Due to a continuum of electronic states present in periodic systems, the description of molecular dynamics on surfaces poses a serious computational challenge. One of the most used families of approaches in these settings are friction theories, which up to a random fluctuating force term are based on the Ehrenfest approach. Yet, a mean-field treatment of electronic degrees of freedom in the Ehrenfest method makes this approach inaccurate in some cases. Our aim is to clarify when Ehrenfest breaks down for molecular dynamics on surfaces. Answering this question provides limits of applicability for more approximate friction theories derived from Ehrenfest. We assess the Ehrenfest method on one-dimensional, numerically exactly solvable models with a large but finite number of electronic states. Using the Landau-Zener formula and the Massey parameter, an expression that determines when Ehrenfest breaks down is deduced.

Variational nonadiabatic dynamics in the moving crude adiabatic representation: Further merging of nuclear dynamics and electronic structure.

A new methodology of simulating nonadiabatic dynamics using frozen-width Gaussian wavepackets within the moving crude adiabatic representation with the on-the-fly evaluation of electronic structure is presented. The main feature of the new approach is the elimination of any global or local model representation of electronic potential energy surfaces; instead, the electron-nuclear interaction is treated explicitly using the Gaussian integration. As a result, the new scheme does not introduce any uncontrolled approximations. The employed variational principle ensures the energy conservation and leaves the number of electronic and nuclear basis functions as the only parameter determining the accuracy. To assess performance of the approach, a model with two electronic and two nuclear spacial degrees of freedom containing conical intersections between potential energy surfaces has been considered. Dynamical features associated with nonadiabatic transitions and nontrivial geometric (or Berry) phases were successfully reproduced within a limited basis expansion.

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.

Geometric phase effects in excited state dynamics through a conical intersection in large molecules: N-dimensional linear vibronic coupling model study.

We investigate geometric phase (GP) effects in nonadiabatic transitions through a conical intersection (CI) in an N-dimensional linear vibronic coupling (ND-LVC) model. This model allows for the coordinate transformation encompassing all nonadiabatic effects within a two-dimensional (2D) subsystem, while the other N - 2 dimensions form a system of uncoupled harmonic oscillators identical for both electronic states and coupled bi-linearly with the subsystem coordinates. The 2D subsystem governs ultra-fast nonadiabatic dynamics through the CI and provides a convenient model for studying GP effects. Parameters of the original ND-LVC model define the Hamiltonian of the transformed 2D subsystem and thus influence GP effects directly. Our analysis reveals what values of ND-LVC parameters can introduce symmetry breaking in the 2D subsystem that diminishes GP effects.

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.

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.