1
|
Patra C, Halder S, Maitra R. Projective quantum eigensolver via adiabatically decoupled subsystem evolution: A resource efficient approach to molecular energetics in noisy quantum computers. J Chem Phys 2024; 160:214122. [PMID: 38836451 DOI: 10.1063/5.0210854] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2024] [Accepted: 05/21/2024] [Indexed: 06/06/2024] Open
Abstract
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 framework based on the principles of nonlinear dynamics and synergetics. Such decoupling emulates total parameter optimization in a lower dimensional manifold, while a mutual synergistic relationship among the parameters is exploited to ensure characteristic accuracy via a non-iterative energy correction. 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 establish 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 the efficient utilization of near-term quantum hardware resources.
Collapse
Affiliation(s)
- Chayan Patra
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Sonaldeep Halder
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Rahul Maitra
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
- Centre of Excellence in Quantum Information, Computing, Science and Technology, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
2
|
Patra C, Agarawal V, Halder D, Chakraborty A, Mondal D, Halder S, Maitra R. A Synergistic Approach towards Optimization of Coupled Cluster Amplitudes by Exploiting Dynamical Hierarchy. Chemphyschem 2023; 24:e202200633. [PMID: 36314661 DOI: 10.1002/cphc.202200633] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2022] [Revised: 10/29/2022] [Indexed: 11/06/2022]
Abstract
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 analyzed 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.
Collapse
Affiliation(s)
- Chayan Patra
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai, 400076, India
| | - Valay Agarawal
- Department of Chemistry, University of Chicago, Chicago, IL, USA, 60637
| | - Dipanjali Halder
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai, 400076, India
| | - Anish Chakraborty
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai, 400076, India
| | - Dibyendu Mondal
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai, 400076, India
| | - Sonaldeep Halder
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai, 400076, India
| | - Rahul Maitra
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai, 400076, India
| |
Collapse
|
3
|
Agarawal V, Roy S, Shrawankar KK, Ghogale M, Bharathi S, Yadav A, Maitra R. A hybrid coupled cluster-machine learning algorithm: Development of various regression models and benchmark applications. J Chem Phys 2022; 156:014109. [PMID: 34998340 DOI: 10.1063/5.0072250] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The iterative solution of the coupled cluster equations exhibits a synergistic relationship among the various cluster amplitudes. The iteration scheme is analyzed 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 [Agarawal et al., J. Chem. Phys. 154, 044110 (2021)]. We develop a coupled cluster-machine learning hybrid scheme where various supervised machine learning strategies are introduced to establish the interdependence between the principal and auxiliary amplitudes on-the-fly. While the coupled cluster equations are solved only to determine the principal amplitudes, the auxiliary amplitudes, on the other hand, are determined via regression as unique functionals of the principal amplitudes. This leads to significant reduction in the number of independent degrees of freedom during the iterative optimization, which saves significant computation time. A few different regression techniques have been developed, which have their own advantages and disadvantages. The scheme has been applied to several molecules in their equilibrium and stretched geometries, and our scheme, with all the regression models, shows a significant reduction in computation time over the canonical coupled cluster calculations without unduly sacrificing the accuracy.
Collapse
Affiliation(s)
- Valay Agarawal
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Samrendra Roy
- Department of Energy Science and Engineering, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Kapil K Shrawankar
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | | | - S Bharathi
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Anchal Yadav
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Rahul Maitra
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
4
|
Olivera-Atencio ML, Morillo M, Casado-Pascual J. System size synchronization. Phys Rev E 2021; 104:064204. [PMID: 35030830 DOI: 10.1103/physreve.104.064204] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2021] [Accepted: 11/23/2021] [Indexed: 06/14/2023]
Abstract
In this paper we bring out the existence of a kind of synchronization associated with the size of a complex system. A dichotomic random jump process associated with the dynamics of an externally driven stochastic system with N coupled units is constructed. We define an output frequency and phase diffusion coefficient. System size synchronization occurs when the average output frequency is locked to the external one and the average phase diffusion coefficient shows a very deep minimum for a range of system sizes. Analytical and numerical procedures are introduced to study the phenomenon, and the results describe successfully the existence of system size synchronization.
Collapse
Affiliation(s)
| | - Manuel Morillo
- Física Teórica, Universidad de Sevilla, Apartado de Correos 1065, Sevilla 41080, Spain
| | - Jesús Casado-Pascual
- Física Teórica, Universidad de Sevilla, Apartado de Correos 1065, Sevilla 41080, Spain
| |
Collapse
|
5
|
Agarawal V, Patra C, Maitra R. An approximate coupled cluster theory via nonlinear dynamics and synergetics: The adiabatic decoupling conditions. J Chem Phys 2021; 155:124115. [PMID: 34598582 DOI: 10.1063/5.0060087] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The coupled cluster iteration scheme is analyzed 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 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 is expressed as a unique function 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.
Collapse
Affiliation(s)
- Valay Agarawal
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai, India
| | - Chayan Patra
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai, India
| | - Rahul Maitra
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai, India
| |
Collapse
|
6
|
Agarawal V, Roy S, Chakraborty A, Maitra R. Accelerating coupled cluster calculations with nonlinear dynamics and supervised machine learning. J Chem Phys 2021; 154:044110. [PMID: 33514076 DOI: 10.1063/5.0037090] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
In this paper, the iteration scheme associated with single reference coupled cluster theory has been analyzed using nonlinear dynamics. The phase space analysis indicates the presence of a few significant cluster amplitudes, mostly involving valence excitations, that 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 a polynomial kernel ridge regression model has been employed to express each of the enslaved amplitudes uniquely in terms of the former set of amplitudes. The subsequent coupled cluster iterations are restricted solely to determine those significant excitations, and the enslaved amplitudes are determined through the already established functional mapping. We will show that our hybrid scheme leads to a significant reduction in the computational time without sacrificing the accuracy.
Collapse
Affiliation(s)
- Valay Agarawal
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai, India
| | - Samrendra Roy
- Department of Energy Science and Engineering, Indian Institute of Technology Bombay, Mumbai, India
| | - Anish Chakraborty
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai, India
| | - Rahul Maitra
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai, India
| |
Collapse
|
7
|
Agarawal V, Chakraborty A, Maitra R. Stability analysis of a double similarity transformed coupled cluster theory. J Chem Phys 2020; 153:084113. [PMID: 32872847 DOI: 10.1063/5.0018383] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this paper, we have analyzed 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ödinger equation is cast as a multivariate time-discrete map, and 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.
Collapse
Affiliation(s)
- Valay Agarawal
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Anish Chakraborty
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Rahul Maitra
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
8
|
|
9
|
Inchiosa ME, Bulsara AR. Nonlinear dynamic elements with noisy sinusoidal forcing: Enhancing response via nonlinear coupling. PHYSICAL REVIEW. E, STATISTICAL PHYSICS, PLASMAS, FLUIDS, AND RELATED INTERDISCIPLINARY TOPICS 1995; 52:327-339. [PMID: 9963437 DOI: 10.1103/physreve.52.327] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
|
10
|
Li J, Cao L, Wu D, Xu S. Application of the slaving principle to nonlinear non-Markovian systems. PHYSICAL REVIEW. E, STATISTICAL PHYSICS, PLASMAS, FLUIDS, AND RELATED INTERDISCIPLINARY TOPICS 1994; 50:1862-1873. [PMID: 9962188 DOI: 10.1103/physreve.50.1862] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
|
11
|
Wischert W, Wunderlin A, Pelster A, Olivier M, Groslambert J. Delay-induced instabilities in nonlinear feedback systems. PHYSICAL REVIEW. E, STATISTICAL PHYSICS, PLASMAS, FLUIDS, AND RELATED INTERDISCIPLINARY TOPICS 1994; 49:203-219. [PMID: 9961208 DOI: 10.1103/physreve.49.203] [Citation(s) in RCA: 28] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
|
12
|
Friedrich R, Haken H. Exact stationary probability distribution for a simple model of a nonequilibrium phase transition. ACTA ACUST UNITED AC 1988. [DOI: 10.1007/bf01313938] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
13
|
Sch�ner G, Haken H. A systematic elimination procedure for Ito stochastic differential equations and the adiabatic approximation. ACTA ACUST UNITED AC 1987. [DOI: 10.1007/bf01307868] [Citation(s) in RCA: 18] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
14
|
Schöner G, Haken H. The slaving principle for stratonovich stochastic differential equations. ACTA ACUST UNITED AC 1986. [DOI: 10.1007/bf01726198] [Citation(s) in RCA: 18] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
15
|
|
16
|
|