Search

Determination of molecular energetics and properties is one of the core challenges in the near-term quantum computing. To this end, hybrid quantum-classical algorithms are preferred for Noisy Intermediate Scale Quantum (NISQ) architectures. The Projective Quantum Eigensolver (PQE) is one such algorithms that optimizes the parameters of the chemistry-inspired unitary coupled cluster (UCC) ansatz using a conventional coupled cluster-like residual minimization. Such a strategy involves the projection of the Schrodinger equation on to linearly independent basis towards the parameter optimization, restricting the ansatz is solely defined in terms of the excitation operators. This warrants the inclusion of high-rank operators for strongly correlated systems, leading to increased utilization of quantum resources. In this manuscript, we develop a methodology for determining the generalized operators in terms of a closed form residual equations in the PQE framework that can be efficiently implemented in a quantum computer with manageable quantum resources. Such a strategy requires the removal of the underlying redundancy in high-rank excited determinants, generated due to the presence of the generalized operators in the ansatz, by projecting them on to an internally contracted lower dimensional manifold. With the application on several molecular systems, we have demonstrated our ansatz achieves similar accuracy to the (disentangled) UCC with singles, doubles and triples (SDT) ansatz, while utilizing an order of magnitude fewer quantum gates. Furthermore, when simulated under stochastic Gaussian noise or depolarizing hardware noise, our method shows significantly improved noise resilience compared to the other members of PQE family and the state-of-the-art variational quantum eigensolver., Comment: 10 pages, 4 figures

Recent quantum algorithms pertaining to electronic structure theory primarily focus on threshold-based dynamic construction of ansatz by selectively including important many-body operators. These methods can be made systematically more accurate by tuning the threshold to include more number of operators into the ansatz. However, such improvements come at the cost of rapid proliferation of the circuit depth, especially for highly correlated molecular systems. In this work, we address this issue by the development of a novel theoretical framework that relies on the segregation of an ansatz into a dynamically selected core principal component, which is, by construction adiabatically decoupled from the remaining operators. This enables us to perform computations involving the principal component using extremely shallow-depth circuits whereas, the effect of the remaining auxiliary component is folded into the energy function via a cost-efficient non-iterative correction, ensuring the requisite accuracy. We propose a formalism that analytically predicts the auxiliary parameters from the principal ones, followed by a suite of non-iterative auxiliary subspace correction techniques with different levels of sophistication. The auxiliary subspace corrections incur no additional quantum resources, yet complement an inadequately expressive core of the ansatz to recover significant amount of electronic correlations. We have numerically validated the resource efficiency and accuracy of our formalism with a number of strongly correlated molecular systems., Comment: 14 pages, 5 figures

Quantum computers hold immense potential in the field of chemistry, ushering new frontiers to solve complex many body problems that are beyond the reach of classical computers. However, noise in the current quantum hardware limits their applicability to large chemical systems. This work encompasses the development of a projective formalism that aims to compute ground-state energies of molecular systems accurately using Noisy Intermediate Scale Quantum (NISQ) hardware in a resource efficient manner. Our approach is reliant upon the formulation of a bipartitely decoupled parameterized ansatz within the disentangled unitary coupled cluster (dUCC) framework based on the principles of synergetics. Such decoupling emulates the total parameter optimization in a lower dimensional manifold, while a mutual synergistic relationship among the parameters is exploited to ensure characteristic accuracy. Without any pre-circuit measurements, our method leads to a highly compact fixed-depth ansatz with shallower circuits and fewer expectation value evaluations. Through analytical and numerical demonstrations, we demonstrate the method's superior performance under noise while concurrently ensuring requisite accuracy in future fault-tolerant systems. This approach enables rapid exploration of emerging chemical spaces by efficient utilization of near-term quantum hardware resources., Comment: 14 pages, 3 figures

Recent advances in quantum information and quantum science have inspired the development of various compact dynamic structured ans\"{a}tze that are expected to be realizable in the Noisy Intermediate-Scale Quantum (NISQ) devices. However, such ans\"{a}tze construction strategies hitherto developed involve considerable measurements, and thus they deviate significantly in NISQ platform from their ideal structures. Therefore, it is imperative that the usage of quantum resources must be minimized while retaining the expressivity and dynamical structure of the ansatz that can adapt itself depending on the degree of correlation. We propose a novel ansatz construction strategy based on the \textit{ab-initio} many-body perturbation theory that requires \textit{no} pre-circuit measurement and thus it remains structurally unaffected by any hardware noise. The accuracy and quantum complexity associated with the ansatz are solely dictated by a pre-defined perturbative order as desired and hence are tunable. Furthermore, the underlying perturbative structure of the ansatz construction pipeline enables us to decompose any high-rank excitation that appears in higher perturbative orders into the product of various low-rank operators, and it thus keeps the execution gate-depth to its minimum. With a number of challenging applications on strongly correlated systems, we demonstrate that our ansatz performs significantly better, both in terms of accuracy, parameter count and circuit depth, in comparison to the allied unitary coupled cluster based ans\"{a}tze., Comment: 12 pages, 3 figures

The development of various dynamic ansatz-constructing techniques has ushered in a new era, rendering the practical exploitation of Noisy Intermediate-Scale Quantum (NISQ) hardware for molecular simulations increasingly viable. However, they exhibit substantial measurement costs during their execution. This work involves the development of a novel protocol that capitalizes on regenerative machine learning methodologies and many-body perturbation theoretic measures to construct a highly expressive and shallow ansatz within the variational quantum eigensolver (VQE) framework. The machine learning methodology is trained with the basis vectors of a low-rank expansion of the N-electron Hilbert space to identify the dominant high-rank excited determinants without requiring a large number of quantum measurements. These selected excited determinants are iteratively incorporated within the ansatz through their low-rank decomposition. The reduction in the number of quantum measurements and ansatz depth manifests in the robustness of our method towards hardware noise, as demonstrated through numerical applications. Furthermore, the proposed method is highly compatible with state-of-the-art neural error mitigation techniques. This approach significantly enhances the feasibility of quantum simulations in molecular systems, paving the way for impactful advancements in quantum computational chemistry., Comment: 10 Pages, 5 Figures

The calculation of molecular excited states is critically important to decipher a plethora of molecular properties. In this manuscript, we develop an equation of motion formalism on top of a bi-exponentially parametrized ground state wavefunction towards the determination of excited states. While the ground state bi-exponential parametrization ensures an accurate description of the wavefunction through the inclusion of high-rank correlation effects, the excited state is parametrized by a novel linear response operator with an effective excitation rank beyond two. To treat the ground and the excited states in the same footings, in addition to the conventional one and two-body response operators, we introduced certain two-body "generalized" response operators with an effective excitation rank of one. We introduce a projective formulation towards the determination of the perturbed amplitudes for the set of "generalized" operators. Our formulation entails a significantly small number of unknown parameters and is shown to be highly accurate compared to allied methods for several difficult chemical systems., Comment: 13 pages

Variational Quantum Eigensolver (VQE) provides a lucrative platform to determine molecular energetics in near-term quantum devices. While the VQE is traditionally tailored to determine the ground state wavefunction with the underlying Rayleigh-Ritz principle, the access to specific symmetry-adapted excited states remains elusive. This often requires high depth circuit or additional ancilla qubits along with prior knowledge of the ground state wavefunction. We propose a unified VQE framework that treats the ground and excited states in the same footings. With the knowledge of the irreducible representations of the spinorbitals, we construct a multi-determinantal reference that is adapted to a given spatial symmetry where additionally, the determinants are entangled through appropriate Clebsch-Gordan coefficients to ensure the desired spin-multiplicity. We introduce the notion of totally symmetric, spin-scalar unitary which maintains the purity of the reference at each step of the optimization. The state-selectivity safeguards the method against any variational collapse while leading to any targeted low-lying eigenroot of arbitrary symmetry. The direct access to the excited states shields our approach from the cumulative error that plagues excited state calculations in a quantum computer and with few parameter count, it is expected to be realized in near-term quantum devices., Comment: 12 pages, 5 figures

The recently developed Projective Quantum Eigensolver (PQE) offers an elegant procedure to evaluate the ground state energies of molecular systems on quantum computers. However, the noise in available quantum hardware can result in significant errors in computed outcomes, limiting the realization of quantum advantage. Although PQE comes equipped with some degree of inherent noise resilience, any practical implementation with apposite accuracy would require additional routines to suppress the errors further. In this work, we propose a way to enhance the efficiency of PQE by developing an optimal framework for introducing Zero Noise Extrapolation (ZNE) in the nonlinear iterative procedure that outlines the PQE; leading to the formulation of ZNE-PQE. For this method, we perform a detailed analysis of how various components involved in it affect the accuracy and efficiency of the reciprocated energy convergence trajectory. Moreover, we investigate the reasons behind the improvements observed in ZNE-PQE over conventional PQE by performing a comparative analysis of their residue norm landscape. This approach is expected to facilitate practical applications of quantum computing in fields related to molecular sciences, where it is essential to determine molecular energies accurately.

The dual exponential coupled cluster (CC) theory proposed by Tribedi et al.[J. Chem. Theory Comput. 2020, 16, 10, 6317-6328] performs significantly better than the coupled cluster theory with singles and doubles excitations (CCSD) due to the implicit inclusion of high-rank excitations. The high-rank excitations are included through the action of a set of vacuum annihilating scattering operators that act non-trivially on certain correlated wavefunctions and are determined via a set of local denominators involving the energy difference between certain excited states. This often leads the theory to be prone to intruders. In this manuscript, we show that restricting the correlated wavefunction, on which the scattering operators act upon, to be spanned by only the singlet paired determinants can avoid the catastrophic breakdown. For the first time, we present two nonequivalent approaches to arrive at the working equations, viz. the projective approach with sufficiency conditions and the amplitude form with many-body expansion. While the effect of the triple excitation is quite small around molecular equilibrium geometry, this scheme leads to a better qualitative description of the energetics in the regions of strong correlation. With a number of pilot numerical applications, we have demonstrated the performance of the dual-exponential scheme with both the proposed solution strategies while restricting the excitation subspaces coupled to the corresponding lowest spin channels.

The recently developed Projective Quantum Eigensolver (PQE) has been demonstrated as an elegant methodology to compute the ground state energy of molecular systems in Noisy Intermdiate Scale Quantum (NISQ) devices. The iterative optimization of the ansatz parameters involves repeated construction of residues on a quantum device. The quintessential pattern of the iteration dynamics, when projected as a time discrete map, suggests a hierarchical structure in the timescale of convergence, effectively partitioning the parameters into two distinct classes. In this work, we have exploited the collective interplay of these two sets of parameters via machine learning techniques to bring out the synergistic inter-relationship among them that triggers a drastic reduction in the number of quantum measurements necessary for the parameter updates while maintaining the characteristic accuracy of PQE. Furthermore the machine learning model may be tuned to capture the noisy data of NISQ devices and thus the predicted energy is shown to be resilient under a given noise model., Comment: 13 pages, 5 figures

Recent advancements in quantum information and quantum technology has stimulated a good deal of interest in the development of quantum algorithms for energetics and properties of many-fermionic systems. While the variational quantum eigensolver is the most optimal algorithm in the Noisy Intermediate Scale Quantum era, it is imperative to develop low depth quantum circuits that are physically realizable in quantum devices. Within the unitary coupled cluster framework, we develop COMPASS, a disentangled ansatz construction protocol that can dynamically tailor an optimal ansatz using the one and two-body cluster operators and a selection of rank-two scatterers. The construction of the ansatz may potentially be performed in parallel quantum architecture through energy sorting and operator commutativity prescreening. With significant reduction in the circuit depth towards the simulation of molecular strong correlation, COMPASS is shown to be highly accurate and resilient to the noisy circumstances of the near-term quantum hardware., Comment: 9 pages, 4 figures

The coupled cluster iteration scheme for determining the cluster amplitudes involves a set of nonlinearly coupled difference equations. In the space spanned by the amplitudes, the set of equations are analysed as a multivariate time-discrete map where the concept of time appears in an implicit manner. With the observation that the cluster amplitudes have difference in their relaxation timescales with respect to the distributions of their magnitudes, the coupled cluster iteration dynamics are considered as a synergistic motion of coexisting slow and fast relaxing modes, manifesting a dynamical hierarchical structure. With the identification of the highly damped auxiliary amplitudes, their time variation can be neglected compared to the principal amplitudes which take much longer time to reach the fixed points. We analytically establish the adiabatic approximation where each of these auxiliary amplitudes are expressed as unique parametric functions of the collective principal amplitudes, allowing us to study the optimization with the latter taken as the independent degrees of freedom. Such decoupling of the amplitudes significantly reduces the computational scaling without sacrificing the accuracy in the ground state energy as demonstrated by a number of challenging molecular applications. A road-map to treat higher order post-adiabatic effects is also discussed., Comment: 10 Pages, 4 figures

In this paper, we have developed a unitary variant of a double exponential coupled cluster theory, which is capable of mimicking the effects of connected excitations of arbitrarily high rank, using only rank-one and rank-two parametrization of the wavefunction ansatz. While its implementation in a classical computer necessitates the construction of an effective Hamiltonian which involves infinite number of terms with arbitrarily high many-body rank, the same can easily be implemented in the hybrid quantum-classical variational quantum eigensolver framework with a reasonably shallow quantum circuit. The method relies upon the nontrivial action of a unitary, containing a set of rank-two scattering operators, on entangled states generated via cluster operators. We have further introduced a number of variants of the ansatz with different degrees of expressibility by judiciously approximating the scattering operators. With a number of applications on strongly correlated molecules, we have shown that all our schemes can perform uniformly well throughout the molecular potential energy surface without significant additional implementation cost and quantum complexity over the unitary coupled cluster approach with single and double excitations.

The iterative quantum phase estimation algorithm, applied to calculating the ground state energies of quantum chemical systems, is theoretically appealing in its wide scope of being able to handle both weakly and strongly correlated regimes. However, the goodness of the initial state that is sent as an input to the algorithm could strongly decide the quality of the results obtained. In this work, we generate the initial state by using the classical-quantum hybrid variational quantum eigensolver algorithm with unitary coupled cluster ansatz. We apply the procedure to obtain the ground state energies of the H4 molecule on a circle, as the system exhibits an interplay of dynamic as well as static correlation effects at different geometries. Furthermore, we argue on the importance of static correlation in construction of the reference determinant, and propose a minimally parametrized unitary coupled cluster ansatz, which drastically reduces number of variational parameters while incorporating the static correlation effects in the wavefunction. We demonstrate that a carefully and appropriately prepared initial state can greatly reduce the effects of noise due to sampling in the estimation of the desired eigenphase., Comment: 12 pages, 5 figures

The iteration dynamics of the coupled cluster equations exhibits a synergistic relationship among the cluster amplitudes. The iteration scheme may be viewed as a multivariate discrete-time propagation of nonlinearly coupled equations, which is dictated by only a few principal cluster amplitudes. These principal amplitudes usually correspond to only a few valence excitations, whereas all other cluster amplitudes are enslaved, and behave as auxiliary variables. Staring with a few trial iterations, we employ a supervised machine learning strategy to establish a mapping of the principal and auxiliary amplitudes. We introduce a machine learning-coupled cluster hybrid scheme where the coupled cluster equations are solved only to determine the principal amplitudes, which saves significant computation time. The auxiliary amplitudes, on the other hand, are determined via regression. Few different regression techniques have been introduced to express the auxiliary amplitudes as functions of the principal amplitudes. The scheme has been applied to several molecules in their equilibrium and stretched geometries, and our scheme, with both the regression models, shows a significant reduction in computation time without unduly sacrificing the accuracy., Comment: 28 pages, 5 figures

The coupled cluster iteration scheme is analysed as a multivariate discrete-time map using nonlinear dynamics and synergetics. The nonlinearly coupled set of equations to determine the cluster amplitudes are driven by a fraction of the entire set of the cluster amplitudes. These driver amplitudes enslave all other amplitudes through a synergistic inter-relationship, where the latter class of amplitudes behave as the auxiliary variables. The driver and the auxiliary variables exhibit vastly different time scales of relaxation during the iteration process to reach the fixed points. The fast varying auxiliary amplitudes are small in magnitude, while the driver amplitudes are large, and they have a much longer time scale of relaxation. Exploiting their difference in relaxation time-scale, we employ an adiabatic decoupling approximation, where each of the fast relaxing auxiliary modes are expressed as unique functional of the principal amplitudes. This results in a tremendous reduction in the independent degrees of freedom. On the other hand, only the driver amplitudes are determined accurately via exact coupled cluster equations. We will demonstrate that the iteration scheme has an order of magnitude reduction in computational scaling than the conventional scheme. With a few pilot numerical examples, we would demonstrate that this scheme can achieve very high accuracy with significant savings in computational time., Comment: 11 Pages, 9 figures, 1 tikz diagram, 1 table

The dynamics associated with the time series of the iteration scheme of coupled cluster theory has been analysed. The phase space analysis indicates the presence of a few significant cluster amplitudes, mostly involving valence excitations, which dictate the dynamics, while all other amplitudes are enslaved. Starting with a few initial iterations to establish the inter-relationship among the cluster amplitudes, a supervised Machine Learning scheme with polynomial Kernel Ridge Regression model has been employed to express each of the enslaved variables uniquely in terms of the master amplitudes. The subsequent coupled cluster iterations are restricted to a reduced dimension only to determine those significant excitations, and the enslaved variables are determined through the already established functional mapping. We will show that our scheme leads to tremendous reduction in computational time without sacrificing the accuracy., Comment: 8 pages, 8 figures (including 1 Tikz Cartoon)

In this paper, we have analysed the time series associated with the iterative scheme of a double similarity transformed Coupled Cluster theory. The coupled iterative scheme to solve the ground state Schr{\"o}dinger equation is cast as a multivariate time-discrete map, the solutions show the universal Feigenbaum dynamics. Using recurrence analysis, it is shown that the dynamics of the iterative process is dictated by a small subgroup of cluster operators, mostly those involving chemically active orbitals, whereas all other cluster operators with smaller amplitudes are enslaved. Using Synergetics, we will indicate how the master-slave dynamics can suitably be exploited to develop a novel coupled-cluster algorithm in a much-reduced dimension., Comment: The main manuscript file contains 8 pages and 6 figures