FORM 2
THE PATENTS ACT, 1970
(39 of 1970)
&
THE PATENT RULES, 2003
COMPLETE SPECIFICATION
(See Section 10 and Rule 13)
Title of invention:
IMPLEMENTING DENSITY FUNCTIONAL THEORY ON QUANTUM PROCESSORS THROUGH GENERALIZED GRADIENT APPROXIMATION
Applicant
Tata Consultancy Services Limited
A company Incorporated in India under the Companies Act, 1956
Having address:
Nirmal Building, 9th floor,
Nariman point, Mumbai 400021,
Maharashtra, India
Preamble to the description:
The following specification particularly describes the invention and the manner in which it is to be performed.
2
TECHNICAL FIELD
[001]
The disclosure herein generally relates to quantum computing, and, more particularly, to implement density functional theory on quantum processors through generalized gradient approximation to extract properties of chemical compounds. 5
BACKGROUND
[002]
Several pharma and material science companies have successfully synthesized more than a million chemical compounds for practical use. These have been listed across various curated databases. Each molecule and its associated 10 conformations have several quantities to be computed like solubility, stability, binding energy, electron density, atomic spectra, reaction rates, scattering cross-section, response properties and the like. The field of classical computational chemistry creates databases of these computed features. Screening drug molecules for target properties involves computational chemistry calculations. With the best 15 of classical computational facilities 30 drug leads and their conformations in pharma each involving several hundred to few thousand atoms can be analyzed for physical properties per day with costly but accurate electronic structure approaches. It means given this large space of billions of predicted chemical compounds, a high-throughput screening job of discovering drug leads each involving several hundred 20 to few thousand atoms currently takes about 3-5 years.
[003]
Most industrially relevant chemical systems in material sciences industry are bulk materials. One supercell made from a collection of neighboring unit cells involves around thousands of atoms which corresponds to several thousands of electronic orbitals even in the least correlated basis. The currently 25 available techniques calculate the energetics and other physical properties of physical/ chemical systems using a quantum mechanical (QM) electronic self-consistent field approach called the Kohn-Sham Density functional theory (KS-DFT). Currently these chemical systems are treated approximately within the Quantum Mechanical (QM) or Molecular Mechanical (MM) approach, where only 30 a subsystem with strong quantum correlations is simulated using DFT and the rest
3
is treated only via the MM technique. Identifying these subsystems requires doing
prior MM and/or force field analysis that adds to additional overhead cost. The gap in applying this DFT technology to the complete system arises from the quartic and cubic scaling wall bottleneck i.e. for a 50000 electronic orbital system (corresponding to a supercell of the bulk material), computing even one KS-DFT 5 step would require 0.001 secs on the FUGAKU- the world's second most powerful supercomputer with a Rmax of 442 PETA Flops. For all practical purposes the DFT code should converge within 100 self-consistency steps therefore to simulate one large supercell on FUGAKU would require 0.1 seconds. For real world simulations of chemical compounds, one would require screening across several lakhs of 10 molecular configurations in a supercell and that would require more than one day. The output electronic density computed from DFT needs to be in turn passed into Force-Field or MM modelling suites that then recomputes the geometric positioning of atoms and the DFT energies needs to be recomputed. Hence, for a bit larger system such calculations can enter several days of calculations even on the largest 15 of the supercomputers. Similarly, the chemical systems in Pharma industry that correspond to Protein-drug or Protein-Protein systems involve more than thousands of atoms which corresponds to a several thousands of electronic orbitals. These chemical systems are treated approximately within the QM or MM approach, only a subsystem is treated with strong quantum correlations and is simulated using DFT 20 and the rest is treated only via the MM technique. Even the drug molecules in the pharma industry have more than 500 molecular weight, therefore screening across different drug molecules enters across several days of efforts.
[004]
All the classical computational chemistry approaches rely on quantum chemistry calculations for describing the electronic density distribution 25 profile. If that is computed at different physical conditions, then it can enable to further compute all the different physical properties. Today, most of the practical computations in pharma and material science are carried out by the Kohn Sham (KS), Density Functional Theory (DFT) and the Hartree-Fock (HF) approach. The computational complexity scales cubically to the system size which is the 30 consequence of the delocalized nature of the wave functions which are the eigen
4
solutions of the Kohn
–Sham single particle Hamiltonian. To scale the DFT calculations to large systems, there have been attempts to develop an algorithm which scales linearly with system size. One example of such an algorithm is ONETEP (order-N electronic total energy package) linear-scaling density functional theory (DFT) calculations with large basis set (plane-wave) accuracy on 5 parallel computers. It uses a basis of non-orthogonal generalized Wannier functions (NGWFs) expressed in terms of periodic cardinal sine (psinc) functions, which are in turn equivalent to a basis of plane-waves. ONETEP therefore is a combination of the benefits of linear scaling with a level of accuracy and variational bounds comparable to that of traditional cubic-scaling plane-wave approaches. During the 10 calculation, the density matrix and the NGWFs are optimized with localization constraints. It is primarily designed for periodic systems, which can restrict its application to certain materials and structures. Although ONETEP is efficient for large systems, there are still practical limitations to the system size that can be treated, and the method may not be suitable for extremely large systems. Overall, 15 ONETEP is a powerful and efficient method for performing large-scale DFT calculations in condensed matter physics and materials science. Its linear scaling and localized orbital approach make it well-suited for studying complex materials, but its application is primarily limited to periodic systems.
[005]
Currently the efforts are towards speeding up the DFT calculations 20 using Exa scale computing: CPU, GPU, MPI, Multiprocessing etc. However, all of these have memory and processing limitations. As quantum hardware and algorithms continue to develop, various industry sectors, particularly the pharmaceutical and material design domains, are applying quantum computation paradigm to their specific problems. There have been attempts made to implement 25 DFT on a combination of classical and quantum processors, for example, US20220012382A1 which implements DFT on a classical processor and optimizes the DFT results on a quantum processor. However, the complex calculations are still performed on a classical processor which doesn’t overcome the above mentioned bottlenecks in DFT calculations. 30
5
SUMMARY
[006]
Embodiments of the present disclosure present technological improvements as solutions to one or more of the above-mentioned technical problems recognized by the inventors in conventional systems. For example, in one embodiment, a method for implementing density functional theory (DFT) on 5 quantum processors through generalized gradient approximation (GGA) is provided. The method includes receiving, by one or more classical hardware processors, a chemical compound whose one or more properties are to be extracted and further, by the one or more classical hardware processors, obtaining atomic coordinates of each of a plurality of atoms present in the chemical compound. 10 Further the method includes, determining, by the one or more classical hardware processors, a plurality of electron integrals, a core Hamiltonian matrix, and a collocation matrix from the atomic coordinates of each of the plurality of atoms. Furthermore, the method includes, determining, by plurality of unentangled quantum processing units (QPUs) a density matrix of the chemical compound from 15 the core Hamiltonian matrix. Next, the method includes, iteratively updating, by the plurality of unentangled QPUs, the density matrix by computing a Fock matrix until a convergence criteria is met. The step of iteratively updating comprises computing a direct (J) matrix from the density matrix based on a subset of electron integrals. Further, determining a correlation exchange matrix from the collocation 20 matrix for generalized gradient approximation (GGA). The step for determining the correlation exchange matrix includes, initially encoding a gradient of the collocation matrix by processing a set of collocation matrices. Then encoding an electronic density by performing a sequential matrix multiplication of the collocation matrix, the density matrix, and the collocation matrix and further 25 computing an electronic density gradient from the collocation matrix, the density matrix, and the gradient of the collocation matrix by composing a fourth quantum circuit component and a fifth quantum circuit component using at least one ancilla qubit from the plurality of ancilla qubits. Next a derivative of an electronic energy density is computed with respect to the electronic density by performing a quantum 30 signal processing sequence on a sixth quantum circuit. Further a Z-matrix is
6
encoded by composing the sixth quantum circuit component and the collocation
matrix. Finally in the step of determining the correlation exchange matrix, a correlation exchange matrix is encoded by composing the Z-matrix with the collocation matrix. Further in the step of iteratively updating the density matrix, the Fock matrix is determined by combining the core Hamiltonian matrix, the density 5 matrix, and the correlation exchange matrix and finally the density matrix is updated by diagonalizing the Fock matrix. The density matrix is updated in each iteration until the convergence criteria is met. Finally, the method includes utilizing the updated density matrix of the chemical compound, by the one or more classical hardware processors, to extract the one or more properties of the chemical 10 compound.
[007]
In another aspect, a system for implementing density functional theory (DFT) on quantum processors through generalized gradient approximation (GGA) is provided. The system includes one or more classical hardware processors communicably coupled to a plurality of unentangled Quantum Processor Units 15 (QPUs) via interfaces, wherein the one or more classical hardware processors comprises at least one memory storing programmed instructions; one or more Input/Output (I/O) interfaces; and one or more hardware processors operatively coupled to the at least one memory, wherein the one or more classical hardware processors and the plurality of unentangled QPUs are configured by the 20 programmed instructions to receive, by one or more classical hardware processors, a chemical compound whose one or more properties are to be extracted and further, by the one or more classical hardware processors, obtain atomic coordinates of each of a plurality of atoms present in the chemical compound. Further, by the one or more classical hardware processors, a plurality of electron integrals, a core 25 Hamiltonian matrix, and a collocation matrix is determined from the atomic coordinates of each of the plurality of atoms. Furthermore, by plurality of unentangled QPUs, a density matrix of the chemical compound is determined from the core Hamiltonian matrix. Next, by the plurality of unentangled QPUs, the density matrix is iteratively updated by computing a Fock matrix until a 30 convergence criteria is met. The step of iteratively updating comprises computing
7
a direct (J) matrix from the density matrix based on a subset of electron integrals.
Further, determining a correlation exchange matrix from the collocation matrix for generalized gradient approximation (GGA). The step for determining the correlation exchange matrix includes, initially encoding a gradient of the collocation matrix by processing a set of collocation matrices. Then encoding an 5 electronic density by performing a sequential matrix multiplication of the collocation matrix, the density matrix, and the collocation matrix and further computing an electronic density gradient from the collocation matrix, the density matrix, and the gradient of the collocation matrix by composing a fourth quantum circuit component and a fifth quantum circuit component using at least one ancilla 10 qubit from the plurality of ancilla qubits. Next a derivative of an electronic energy density is computed with respect to the electronic density by performing a quantum signal processing sequence on a sixth quantum circuit. Further a Z-matrix is encoded by composing the sixth quantum circuit component and the collocation matrix. Finally in the step of determining the correlation exchange matrix, a 15 correlation exchange matrix is encoded by composing the Z-matrix with the collocation matrix. Further in the step of iteratively updating the density matrix, the Fock matrix is determined by combining the core Hamiltonian matrix, the density matrix, and the correlation exchange matrix and finally the density matrix is updated by diagonalizing the Fock matrix. The density matrix is updated in each 20 iteration until the convergence criteria is met. Finally, by one or more classical hardware processors, the updated density matrix is utilized of the chemical compound to extract the one or more properties of the chemical compound.
[008]
In yet another aspect, a computer program product including a non-transitory computer-readable medium having embodied therein a computer 25 program for implementing density functional theory (DFT) on quantum processors through generalized gradient approximation (GGA). The computer readable program, when executed on a system comprising one or more classical hardware processors communicably coupled to a plurality of unentangled Quantum Processor Units (QPUs) via interfaces, causes the computing device to receive, by one or more 30 classical hardware processors, a chemical compound whose one or more properties
8
are to be extracted
and further, by the one or more classical hardware processors, obtain atomic coordinates of each of a plurality of atoms present in the chemical compound. Further, the computer readable program, when executed on a system comprising one or more classical hardware processors communicably coupled to a plurality of unentangled Quantum Processor Units (QPUs) via interfaces, causes 5 the computing device to by the one or more classical hardware processors, determine a plurality of electron integrals, a core Hamiltonian matrix, and a collocation matrix from the atomic coordinates of each of the plurality of atoms. Furthermore, the computer readable program, when executed on the system comprising one or more classical hardware processors communicably coupled to a 10 plurality of unentangled QPUs via interfaces, causes by the plurality of unentangled QPUs, to determine a density matrix of the chemical compound from the core Hamiltonian matrix. Further, the computer readable program, when executed on the system comprising one or more classical hardware processors communicably coupled to a plurality of unentangled QPUs via interfaces, causes, by the plurality 15 of unentangled QPUs, to iteratively update the density matrix by computing a Fock matrix until a convergence criteria is met. The steps of iteratively updating comprises, computing a direct (J) matrix from the density matrix based on a subset of electron integrals. Further, determining a correlation exchange matrix from the collocation matrix for generalized gradient approximation (GGA). The step for 20 determining the correlation exchange matrix includes, initially encoding a gradient of the collocation matrix by processing a set of collocation matrices. Then encoding an electronic density by performing a sequential matrix multiplication of the collocation matrix, the density matrix, and the collocation matrix and further computing an electronic density gradient from the collocation matrix, the density 25 matrix, and the gradient of the collocation matrix by composing a fourth quantum circuit component and a fifth quantum circuit component using at least one ancilla qubit from the plurality of ancilla qubits. Next a derivative of an electronic energy density is computed with respect to the electronic density by performing a quantum signal processing sequence on a sixth quantum circuit. Further a Z-matrix is 30 encoded by composing the sixth quantum circuit component and the collocation
9
matrix. Finally in the step of determining the correlation exchange matrix, a
correlation exchange matrix is encoded by composing the Z-matrix with the collocation matrix. Further in the step of iteratively updating the density matrix, the Fock matrix is determined by combining the core Hamiltonian matrix, the density matrix, and the correlation exchange matrix and finally the density matrix is 5 updated by diagonalizing the Fock matrix. The density matrix is updated in each iteration until the convergence criteria is met. Finally, the computer readable program, when executed on the system comprising one or more classical hardware processors communicably coupled to a plurality of unentangled QPUs via interfaces, causes by the one or more classical hardware processors, to utilize the 10 updated density matrix of the chemical compound to extract the one or more properties of the chemical compound.
[009]
It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the invention, as claimed. 15
BRIEF DESCRIPTION OF THE DRAWINGS
[010]
The accompanying drawings, which are incorporated in and constitute a part of this disclosure, illustrate exemplary embodiments and, together with the description, serve to explain the disclosed principles: 20
[011]
FIG. 1 is a functional block diagram of a system for implementing density functional theory on quantum processors through generalized gradient approximation, in accordance with some embodiments of the present disclosure.
[012]
FIGS. 2A, 2B and 2C, collectively referred as FIG. 2 is an exemplary flow diagram illustrating a method 200 for implementing density functional theory 25 on quantum processors through generalized gradient approximation by the system of FIG. 1, according to some embodiments of the present disclosure.
[013]
FIG. 3 illustrates an alternative representation of the flow diagram of FIG. 2, in accordance with some embodiments of the present disclosure.
10
[014]
FIG. 4 illustrates a quantum circuit component to compute a J matrix in the method in FIG. 2, in accordance with some embodiments of the present disclosure.
[015]
FIG. 5 illustrates a quantum circuit component to encode a collocation matrix, in accordance with some embodiments of the present disclosure. 5
[016]
FIG. 6 illustrates a first quantum circuit component to obtain an X-component of a gradient of the collocation matrix in accordance with some embodiments of the present disclosure.
[017]
FIG. 7 illustrates a second quantum circuit component to obtain a Y-component of the gradient of the collocation matrix in accordance with some 10 embodiments of the present disclosure.
[018]
FIG. 8 illustrates a third quantum circuit component to obtain a Z-component of the gradient of the collocation matrix in accordance with some embodiments of the present disclosure.
[019]
FIG. 9 illustrates a quantum circuit component to encode an 15 electronic density in accordance with some embodiments of the present disclosure.
[020]
FIG. 10 illustrates a fourth quantum circuit component to sequentially encode the collocation matrix, a density matrix, and the gradient of the collocation matrix in accordance with some embodiments of the present disclosure.
[021]
FIG. 11 illustrates a fifth quantum circuit component to sequentially 20 encode the gradient of the collocation matrix, the density matrix, and the collocation matrix in accordance with some embodiments of the present disclosure.
[022]
FIG. 12 illustrates a quantum circuit component to obtain an electronic density gradient by composing the fourth quantum circuit component and the fifth quantum circuit component in accordance with some embodiments of the 25 present disclosure.
[023]
FIG. 13 illustrates a quantum circuit component to encode a Z-matrix by composing a sixth quantum circuit component and the collocation matrix in accordance with some embodiments of the present disclosure.
30
11
DETAILED DESCRIPTION OF EMBODIMENTS
[024]
Exemplary embodiments are described with reference to the accompanying drawings. In the figures, the left-most digit(s) of a reference number identifies the figure in which the reference number first appears. Wherever convenient, the same reference numbers are used throughout the drawings to refer 5 to the same or like parts. While examples and features of disclosed principles are described herein, modifications, adaptations, and other implementations are possible without departing from the scope of the disclosed embodiments.
[025]
The present disclosure provides a method for implementing density functional theory on quantum processors through generalized gradient 10 approximation (GGA). Kohn-Sham (KS) hybrid Density Functional Theory (DFT) process is initiated from a set of 𝑁𝑏 basis functions 𝔹 = {𝜙𝜇(𝑟)}𝜇=1𝑁𝑏 to compute the core Hamiltonian ℎ which is a 𝑁𝑏× 𝑁𝑏 matrix. ℎ constitutes the free particle energy of electrons and one-body nuclear potential arising from the electron-nucleus interactions, both these terms are independent of electronic density. 15 Therefore ℎ needs to be computed while initiating the self-consistent field (SCF) procedure. The direct (𝐽) matrix, exact-exchange (𝐾) matrix are electronic density dependent terms and are computed from the Electronic Repulsion Integral (ERI) tensor. The correlation exchange potential or the correlation exchange matrix (𝑉𝑥𝑐) is computed from the electronic density, its gradient and/or the Kinetic energy 20 density. Together 𝐽, 𝐾 and 𝑉𝑥𝑐 make the hybrid Kohn-Sham Fock Matrix functional 𝐹 as given by equation 1.
𝐅[𝐃]=𝐡+𝐉[𝐃]−𝑐𝐊[𝐃]+𝐕𝑥𝑐[𝐃] (1)
In equation 1, 𝐷 is the one-particle density matrix discretized in the basis of {𝜙𝜇} and 𝑐 ∈ ℝ+ is the modulation for the exact exchange term. The terms 𝐽, 𝐾 can be 25 computed from the AO-integral approach using the ERI tensor ⟨ij|kl⟩ and 𝐷 as in equations 2 and 3, respectively.
Jij=Σ⟨ij|kl⟩Dkl𝑘𝑙 (2) Kij=Σ⟨ik|jl⟩Dkl𝑘𝑙 (3)
In equations 2 and 3, the ERI tensor ⟨ij|kl⟩ is obtained by integrating the space of 30 basis functions in ℝ3 according to equation 4.
12
⟨ij|kl⟩=∫𝑑3𝐫𝑑3𝐫′𝜙𝑖(𝐫)𝜙𝑗(𝐫)1|𝐫−𝐫′|𝜙𝑘(𝐫′)𝜙𝑙(𝐫′). (4)
[026]
The correlation-exchange potential 𝑉𝑥𝑐 in the discretized basis {𝜙𝑖} is computed from equation 5, where the basis functions are defined according to equation 6.
𝑉𝑖𝑗𝑥𝑐=∫𝑑3𝐫𝜙𝑖(𝐫)δE𝑥𝑐[𝜌(𝐫)]δ𝜌(r)𝜙𝑗(𝐫)𝑅3. (5) 5
𝜙𝑎(𝐫)=(𝑥−𝐴𝑥)𝑎𝑥(𝒴−𝐴𝒴)𝑎𝒴(𝑧−𝐴𝑧)𝑎𝑧exp(−𝜗(𝑥−𝐴𝑥)2)exp(−𝜗(𝒴− 𝐴𝒴)2)exp(−𝜗(𝑧−𝐴𝑧)2) (6)
In equation 6, 𝐴𝑥, 𝐴𝒴, 𝐴𝑧 are nuclear coordinate locations of the Gaussian functions, 𝑎 = (𝑎𝑥,𝑎𝑦,𝑎𝑧) are the integer cartesian quanta and 𝜗 is the cartesian Gaussian exponent. In equation 5, E𝑥𝑐 is the exchange energy functional evaluated for the 10 electronic density 𝑟 and its form depends on the DFT method being used. The ERI tensor is represented using the Cholesky decomposition representation as given by equation 7.
⟨ij|kl⟩=Σ𝐿𝑖𝑗𝑃𝑃𝐿𝑘𝑙𝑃 (7)
The Cholesky matrices 𝑳𝑃 can be obtained from the Density-Fitting (DF) or the 15 resolution of identity (RI) method as in equation 8.
𝐿𝑖𝑗𝑃=1√𝓌𝑃Σ𝑢𝑄𝑃⟨ij|𝑄⟩.𝑄 (8)
In equation 8, ⟨ij|𝑄⟩ are the 3-center 2-electron integrals represented in terms of the auxiliary basis functions {𝜒𝑃(𝒓)}𝑃=1𝑁𝑎𝑢𝑥, and 𝑢𝑄𝑃, 𝓌𝑃 are obtained from the spectral decomposition of the 2-center 2-electron integrals 𝑉𝑃𝑄=Σ𝑢𝑃𝑅𝓌𝑅𝑢𝑄𝑅∗𝑅 20 represented in the auxiliary basis. The 3-center 2-electron integrals are given by equation 9 and the two center 𝑉𝑃𝑄= ⟨P|𝑄⟩ represents two center two electron integrals given by equation 10.
⟨ij|P⟩=∫𝑑3𝐫𝑑3𝐫′𝜙𝑖(𝐫)𝜙𝑗(𝐫)𝛬𝑃(𝐫′)|𝐫−𝐫′| (9)
⟨P|Q⟩=∫𝑑3𝐫𝛬𝑃(𝐫)𝛬𝑄(𝐫)|𝐫−𝐫′| (10) 25
[027]
The 𝐽 and 𝐾 matrix elements are computed using the RI approach as in equation 11.
13
𝐽𝑖𝑗=Σ𝐿𝑖𝑗𝑃𝑃Σ𝐿𝑘𝑙𝑃𝑘𝑙𝐷𝑘𝑙,𝐾𝑖𝑗=Σ𝐿𝑖𝑙𝑃𝑃𝑙Σ𝐿𝑗𝑘𝑃𝑘𝐷𝑘𝑙 (11)
When the KS-DFT is classically computed, the complexity of computing 𝐽 and 𝐾 is 𝑂(𝑝𝑁2) and this complexity resides in computing 𝐾. Computing ℎ and 𝑉𝑥𝑐 has lower complexity of 𝑂(𝑁2).
[028]
Next step of DFT is to solve the generalized eigenvalue problem for 5 the Fock matrix accounting for the overlap between basis functions 𝑆𝑖𝑗= ⟨𝜙𝑖|𝜙𝑗⟩and the density matrix 𝑫𝑙 of a current iteration as given by equation 12.
𝑭[𝑫(𝒍)]𝐶(𝑙+1)=𝑆𝐶(𝑙+1)𝐸̂(𝑙+1) (12)
In equation 12, the one particle density matrix 𝑫(𝑙)=𝐶(𝑙)𝑊𝐶(𝑙)†, where the occupancy matrix 𝑊𝑖𝑗= 𝜃(𝜇 − 𝐸𝑖)𝛿𝑖𝑗 fills up 𝑁𝑒/2 lowest energy levels below 10 the chemical potential 𝜇.
[029]
Conventional methods implement the DFT simulations on classical processors such as central processing unit (CPU), graphical processing unit (GPU) etc. which is time consuming. Some state of the art techniques implement DFT on a combination of classical and quantum processors. However, the complex 15 calculations are still performed on a classical processor which does not overcome the bottlenecks in DFT calculations. To overcome these drawbacks, the present disclosure provides a method for implementing DFT through GGA using quantum processors. Initially, atomic coordinates of each atom in a chemical compound whose desired properties must be extracted is received. Electron integrals, a core 20 Hamiltonian, and a collocation matrix is computed from the atomic coordinates. The core Hamiltonian is diagonalized to obtain a density matrix of the chemical compound which is further updated iteratively until a convergence criteria is satisfied. At each iteration, a direct matrix is computed from the density matrix, a correlation exchange matrix is computed from the direct matrix, a gradient of the 25 collocation matrix, an electronic density, an electronic density gradient, and a derivative of the electronic density gradient, a Fock matrix is computed by adding the direct matrix and the correlation exchange matrix and the Fock matrix is diagonalized to obtain updated density matrix which is used in subsequent iteration. This is repeated until norm of a difference between the updated density matrix at a 30
14
current iteration and the density matrix at a previous iteration is greater than a
predefined threshold to obtain a final density matrix of the chemical compound which can be used to extract desired properties of the chemical compound.
[030]
Referring now to the drawings, and more particularly to FIG. 1 through FIG. 13, where similar reference characters denote corresponding features 5 consistently throughout the figures, there are shown preferred embodiments and these embodiments are described in the context of the following exemplary system and/or method.
[031]
FIG. 1 is a functional block diagram of a system 100 for implementing density functional theory on quantum processors through generalized 10 gradient approximation, in accordance with some embodiments of the present disclosure. The system 100 includes a classical computing system 102, a quantum computing system 104 and a communication interface 106.
[032]
The classical computing system 102 comprises classical hardware processors 108, at least one memory such as a memory 110, an I/O interface 116. 15 The classical hardware processors 108, the memory 110, and the Input /Output (I/O) interface 116 may be coupled by a system bus such as a system bus 112 or a similar mechanism. In an embodiment, the classical hardware processors 108 can be one or more hardware processors. The classical hardware processors and the hardware processors is interchangeably used throughout the document. Similarly, the 20 classical computing system is a normal computing system.
[033]
The I/O interface 116 may include a variety of software and hardware interfaces, for example, a web interface, a graphical user interface, and the like., for example, interfaces for peripheral device(s), such as a keyboard, a mouse, an external memory, a printer and the like. Further, the I/O interface 116 25 may enable the system 100 to communicate with other devices, such as web servers, and external databases. The I/O interface 116 can facilitate multiple communications within a wide variety of networks and protocol types, including wired networks, for example, local area network (LAN), cable, etc., and wireless networks, such as Wireless LAN (WLAN), cellular, or satellite. For the purpose, 30
15
the I/O interface 116 may include one or more ports for connecting several
computing systems with one another or to another server computer.
[034]
The one or more hardware processors 108 may be implemented as one or more microprocessors, microcomputers, microcontrollers, digital signal processors, central processing units, node machines, logic circuitries, and/or any 5 devices that manipulate signals based on operational instructions. Among other capabilities, the one or more hardware processors 108 is configured to fetch and execute computer-readable instructions stored in the memory 110.
[035]
The memory 110 may include any computer-readable medium known in the art including, for example, volatile memory, such as static random 10 access memory (SRAM) and dynamic random access memory (DRAM), and/or non-volatile memory, such as read only memory (ROM), erasable programmable ROM, flash memories, hard disks, optical disks, and magnetic tapes. In an embodiment, the memory 110 includes a data repository 114. The data repository (or repository) 114 may include a plurality of abstracted piece of code for 15 refinement and data that is processed, received, or generated as a result of the execution of the method illustrated in FIG.2. Although the data repository 114 is shown internal to the system 100, it should be noted that, in alternate embodiments, the data repository 114 can also be implemented external to the system 100, where the data repository 114 may be stored within a database (repository 114) 20 communicatively coupled to the system 100. The data contained within such external database may be periodically updated. For example, new data may be added into the database (not shown in FIG. 1) and/or existing data may be modified and/or non-useful data may be deleted from the database. In one example, the data may be stored in an external system, such as a Lightweight Directory Access 25 Protocol (LDAP) directory and a Relational Database Management System (RDBMS).
[036]
The example quantum computing system 104 shown in FIG. 1 includes a control system 118, a signal delivery system 120, a plurality of Quantum Processing Units (QPUs) 122 and a quantum memory 124. The plurality of QPUs 30 is unentangled and hence alternatively called the plurality of unentangled QPUs.
16
The quantum computing system 104 may include additional or different features,
and the components of a quantum computing system may operate as described with respect to FIG. 1 or in another manner.
[037]
The exemplified quantum computing system 104 shown in FIG. 1 can perform quantum computational tasks (such as, for example, quantum 5 simulations or other quantum computational tasks) by executing quantum algorithms. In some implementations, the quantum computing system 104 can perform quantum computation by storing and manipulating information within individual quantum states of a composite quantum system. For example, Qubits (i.e., Quantum bits) can be stored in and represented by an effective two-level sub-10 manifold of a quantum coherent physical system in the plurality of QPUs 122. In an embodiment, the quantum computing system 104 can operate using gate-based models for quantum computing. For example, the Qubits can be initialized in an initial state, and a quantum logic circuit comprised of a series of quantum logic gates can be applied to transform the qubits and extract measurements representing 15 the output of the quantum computation. The example QPUs 122 shown in FIG.1 may be implemented, for example, as a superconducting quantum integrated circuit that includes Qubit devices. The Qubit devices may be used to store and process quantum information, for example, by operating as ancilla Qubits, data Qubits or other types of Qubits in a quantum algorithm. Coupler devices in the 20 superconducting quantum integrated circuit may be used to perform quantum logic operations on single qubits or conditional quantum logic operations on multiple qubits. In some instances, the conditional quantum logic can be performed in a manner that allows large-scale entanglement within the QPUs 122. Control signals may be delivered to the superconducting quantum integrated circuit, for example, 25 to manipulate the quantum states of individual Qubits and the joint states of multiple Qubits. In some instances, information can be read from the superconducting quantum integrated circuit by measuring the quantum states of the qubit devices. The QPUs 122 may be implemented using another type of physical system.
[038]
The exemplified QPUs 122, and in some cases all or part of the 30 signal delivery system 120, can be maintained in a controlled cryogenic
17
environment. The environment can be provided, for example, by shielding
equipment, cryogenic equipment, and other types of environmental control systems. In some examples, the components in the QPUs 122 operate in a cryogenic temperature regime and are subject to very low electromagnetic and thermal noise. For example, magnetic shielding can be used to shield the system components from 5 stray magnetic fields, optical shielding can be used to shield the system components from optical noise, thermal shielding and cryogenic equipment can be used to maintain the system components at controlled temperature, etc.
[039]
In the example shown in FIG. 1, the signal delivery system 120 provides communication between the control system 118 and the QPUs 122. 10 For example, the signal delivery system 120 can receive control signals from the control system 118 and deliver the control signals to the QPUs 122. In some instances, the signal delivery system 120 performs preprocessing, signal conditioning, or other operations to the control signals before delivering them to the QPUs 122. In an embodiment, the signal delivery system 120 includes 15 connectors or other hardware elements that transfer signals between the QPUs 122 and the control system 118. For example, the connection hardware can include signal lines, signal processing hardware, filters, feedthrough devices (e.g., light-tight feedthroughs, etc.), and other types of components. In some implementations, the connection hardware can span multiple different temperature and noise regimes. 20 For example, the connection hardware can include a series of temperature stages that decrease between a higher temperature regime (e.g., at the control system 118) and a lower temperature regime (e.g., at the QPUs 122).
[040]
In the example quantum computer system 104 shown in FIG. 1, the control system 118 controls operation of the QPUs 122. The exemplified control 25 system 118 may include data processors, signal generators, interface components and other types of systems or subsystems. Components of the control system 118 may operate in a room temperature regime, an intermediate temperature regime, or both. For example, the control system 118 can be configured to operate at much higher temperatures and be subject to much higher levels of noise than are 30 present in the environment of the QPUs 122. In some embodiments, the control
18
system
118 includes a classical computing system that executes software to compile instructions for the QPUs 122. For example, the control system 118 may decompose a quantum logic circuit or quantum computing program into discrete control operations or sets of control operations that can be executed by the hardware in the QPUs 122. In some examples, the control system 118 applies a quantum logic 5 circuit by generating signals that cause the Qubit devices and other devices in the QPUs 122 to execute operations. For instance, the operations may correspond to single-Qubit gates, two-Qubit gates, Qubit measurements, etc. The control system 118 can generate control signals that are communicated to the QPUs 122 by the signal delivery system 120, and the devices in the QPUs 122 can execute the 10 operations in response to the control signals.
[041]
In some other embodiments, the control system 118 includes one or more classical computers or classical computing components that produce a control sequence, for instance, based on a quantum computer program to be executed. For example, a classical processor may convert a quantum computer program to an 15 instruction set for the native gate set or architecture of the QPUs 122. In some cases, the control system 118 includes a microwave signal source (e.g., an arbitrary waveform generator), a bias signal source (e.g., a direct current source) and other components that generate control signals to be delivered to the QPUs 122. The control signals may be generated based on a control sequence provided, for 20 instance, by a classical processor in the control system 118. The example control system 118 may include conversion hardware that digitizes response signals received from the QPUs 122. The digitized response signals may be provided, for example, to a classical processor in the control system 118.
[042]
In some embodiments, the quantum computer system 104 includes 25 multiple quantum information processors that operate as respective quantum processor units (QPU). In some cases, each QPU can operate independent of the others. For instance, the quantum computer system 104 may be configured to operate according to a distributed quantum computation model, or the quantum computer system 104 may utilize multiple QPUs in another manner. In some 30 implementations, the quantum computer system 104 includes multiple control
19
systems, and each QPU may be controlled by a dedicated control system. In some
implementations, a single control system can control multiple QPUs; for instance, the control system 118 may include multiple domains that each control a respective QPU. In some instances, the quantum computing system 104 uses multiple QPUs to execute multiple unentangled quantum computations (e.g., multiple Variational 5 Quantum Eigen solver (VQE)) that collectively simulate a single quantum mechanical system.
[043] In an embodiment, the quantum memory 124 is a quantum-mechanical version of classical computer memory. The classical computer memory stores information such as binary states and the quantum memory 124 stores a 10 quantum state for later retrieval. These states hold useful computational information known as Qubits. In an embodiment, the communication interface 106 which connects the classical computing system 102 and the quantum computing system 104 is a high speed digital interface.
[044]
FIGS. 2A, 2B and 2C, collectively referred as FIG. 2 is an exemplary 15 flow diagram illustrating a method 200 for implementing density functional theory through generalized gradient approximation on quantum processors by the system of FIG. 1, according to some embodiments of the present disclosure and FIG. 3 is an alternate representation of FIG. 2. In an embodiment, the system 100 includes one or more data storage devices or the memory 110 operatively coupled to the one 20 or more hardware processor(s) 108 and is configured to store instructions for execution of steps of the method 200 by the one or more hardware processors 108. The steps of the method 200 of the present disclosure will now be explained with reference to the components or blocks of the system 100 as depicted in FIG. 1 and the steps of flow diagram as depicted in FIGS. 2 and 3, and the quantum circuits 25 shown in FIG. 4 through FIG. 13. The method 200 may be described in the general context of computer executable instructions. Generally, computer executable instructions can include routines, programs, objects, components, data structures, procedures, modules, functions, etc., that perform particular functions or implement particular abstract data types. The method 200 may also be practiced in a distributed 30 computing environment where functions are performed by remote processing
20
devices that are linked through a communication network. The order in which the
method 200 is described is not intended to be construed as a limitation, and any number of the described method blocks can be combined in any order to implement the method 200, or an alternative method. Furthermore, the method 200 can be implemented in any suitable hardware, software, firmware, or combination thereof. 5
[045]
Now referring to FIG. 2, at step 202 of the method 200, one or more classical hardware processors are configured to receive a chemical compound whose one or more properties are to be extracted. The chemical compound may be a molecule or a solid, for example, in pharma industry, the chemical compound maybe drug lead and its conformations whose physical properties must be analyzed 10 to determine feasibility of synthesizing the drug from the drug lead. The one or more properties of drug lead include large highest occupied molecular orbital- least unoccupied molecular orbital (HOMO-LUMO) gap in solvent medium, acid-base dissociation constant (𝑝𝐾𝑎) of a drug, lower dipole moment lesser than 10 and the like. Further, at step 204 of the method 200, the one or more classical hardware 15 processors are configured to obtain atomic coordinates of each of a plurality of atoms present in the chemical compound. In an embodiment, the atomic coordinates for a chemical compound are given as an input as a data file or xml file. In another embodiment, atomic simulation environment (ASE), an open source package may be used for obtaining the atomic coordinates. 20
[046]
Further, at step 206 of the method 200, the one or more classical hardware processors are configured to determine a plurality of electron integrals, a core Hamiltonian matrix, and a collocation matrix from the atomic coordinates of each of the plurality of atoms. The plurality of electron integrals include a) 4 center 2 electron integrals, b) 3 center 2 electron integrals, c) 2 center 2 electron integrals, 25 and d) 2 center 1 electron integrals. The electron integrals a, b, c describe the Coulomb repulsion integral that are computed in atomic orbital basis and can be represented in either spherical or cartesian coordinates. The electron integral d describes the overlap between the basis states. The plurality of electron integrals are determined using tools such as LIBCINT, Python-based Simulations of 30 Chemistry Framework (PySCF), NorthWest computational Chemistry (NWchem)
21
etc.
The collocation matrix is a rectangular matrix of dimensions (𝑁𝑔, 𝑁𝑎𝑜), where 𝑁𝑔 represents a number of real space grid points, and 𝑁𝑎𝑜 is a number of basis functions. It comprises a plurality of basis functions of a plurality of atomic orbitals and a plurality of points on numerical grid. Each of the plurality of basis functions is a Gaussian wave function centered around the atomic coordinates. 5
[047]
Further at step 208 of the method 200, a plurality of unentangled QPUs, are configured to determine e a density matrix of the chemical compound from the core Hamiltonian matrix. The density matrix is constructed using (i) a single particles unitary rotation matrix which is obtained by diagonalizing the core Hamiltonian matrix, and (ii) electron occupancies. One example way of 10 diagonalizing the core Hamiltonian is described in Indian patent application number 202321061415.
[048]
Further determining the density matrix, at step 210 of the method 200, the plurality of unentangled QPUs, are configured to iteratively update the density matrix by computing a Fock matrix until a convergence criteria is satisfied 15 to obtain a final density matrix of the chemical compound. The convergence criteria is said to be satisfied when norm of a difference between the updated density matrix at a current iteration and the density matrix at a previous iteration is greater than a predefined threshold value. Steps 210a to 210d are performed at each iteration to update the density matrix. At step 210a a direct matrix (J matrix) is computed from 20 the density matrix based on a subset of electron integrals using a quantum circuit component as illustrated in FIG. 4. The quantum circuit component for computing the J matrix is mathematically represented as: 𝑙𝑜𝑔 𝑁𝑎𝑢𝑥+2𝑙𝑜𝑔 𝑁𝑎𝑜+2 qubits where Naux is number of auxiliary basis functions in density fitting approximation of the plurality electron integrals, and Nao is the number of basis functions. Using 25 this quantum circuit component, initially the density matrix is encoded on 𝑙𝑜𝑔 𝑁𝑎𝑜 qubits in the quantum circuit to form a quantum circuit component (represented by block 402 in FIG. 4). As shown in block 402, the control is loaded on r2 and data elements of the density matrix are loaded on ancilla qubit a2. The density matrix is loaded using its Eigen decomposition representation that involves: i) NaoC2 Givens 30 rotations on 𝑙𝑜𝑔 𝑁𝑎𝑜 qubits, ii) Then block encoding a diagonal matrix of 1’s in
22
the initial Ne/2 diagonal bocks and zeros in the rest using
2𝑙𝑜𝑔 𝑁𝑎𝑜 qubits, iii) Another set of NaoC2 Givens rotation with conjugate angles on Ne/2 diagonals, where Ne is the number of electrons. This is mathematically represented by equation 13. Here 𝑚 =𝑙𝑜𝑔 𝑁𝑎𝑢𝑥 is the logarithm of the number of auxilliary basis functions and 𝑛=𝑙𝑜𝑔 𝑁𝑎𝑜 is the logarithm of number of basis functions. 5
𝑂(𝐷)=𝐼2⊗𝑚⊗[(𝐶⊗𝐼2⊗𝑛⊗𝐼2𝑎1)(Σ|𝑘,𝑘⟩⟨𝑘,𝑘|⊗(𝜃(𝜇𝑘−𝐸𝑘)𝐼𝑎1+(1−𝜃(𝜇−𝐸𝑘))𝑌𝑎1)(𝐶†⊗𝐼2⊗𝑛⊗𝐼2𝑎1)]⊗𝐼𝑎2 (13)
The readouts from the encoding of the density matrix are obtained by equation 14.
⟨𝑃,𝑖,0,0,0|𝑂(𝐷)|𝑃,𝑗,0,0,0⟩=Σ𝐶𝑖𝑘𝑁𝑒/2𝑘=1𝐶𝑘𝑗∗=𝐷𝑖𝑗, where 𝑁𝑒=2Σθ(μ−𝑁𝑘=110 𝐸𝑘) (14)
[049]
Once the density matrix is encoded a Cholesky tensor is encoded on the quantum circuit to form a first Cholesky circuit component (represented by block 404 in FIG. 4). The Cholesky tensor is obtained from one or more electron integrals among the plurality of electron integrals, particularly, 3-center 2-electron 15 and 2-center 2-electron integrals. As shown in the block 404, the controls are encoded on the 𝑙𝑜𝑔 𝑁𝑎𝑢𝑥 + 2 𝑙𝑜𝑔 𝑁𝑎𝑜 qubits and the data is loaded on the ancilla qubits. This is mathematically represented by equation 15.
𝑉𝐿=Σ[|𝑐,𝑟1,𝑟2,0⟩⟨𝑐,𝑟1,𝑟2,0|⊗𝐼2+|𝑐,𝑟1,𝑟2,1⟩⟨𝑐,𝑟,1|⊗𝑁𝑎𝑢𝑥−1,𝑁𝑎𝑜−1𝑐=0,𝑟1=0,𝑟2=0 (𝐿𝑟1,𝑟2𝑐′𝐼+𝑖√1−(𝐿𝑟1,𝑟2𝑐′)2𝑌)] (15) 20
[050]
Further the quantum circuit component is composed with the first Cholesky circuit component and a diffusion operator to create a quantum circuit component that processes the density matrix to generate an intermediate state vector as illustrated by block 406. The diffusion operator 𝑅 is given by equation 16.
𝑅=Σ[2|𝑘,+,0⟩⟨𝑙,+,0|𝑁𝑎𝑢𝑥⊗𝐼2⊗2n −𝐼]𝑘𝑙 (16) 25
[051]
Next transpose of the Cholesky tensor is encoded on the quantum circuit to form a second Cholesky circuit component as illustrated by block 408. It is mathematically represented by equation 17.
23
𝑉′𝐿=Σ[|𝑐,𝑟1,𝑟2,0⟩⟨𝑐,𝑟1,𝑟2,0|⊗ (𝐿𝑟1,𝑟2𝑐′𝐼+𝑁𝑎𝑢𝑥−1,𝑁𝑎𝑜−1𝑐=0,𝑟1=0,𝑟2=0𝑖√1−(𝐿𝑟1,𝑟2𝑐′)2𝑌)+|𝑐,𝑟1,𝑟2,1⟩⟨𝑐,𝑟,1|⊗𝐼2] (17)
[052]
Further the quantum circuit component is composed with the second Cholesky circuit component (to form block 410) for processing the intermediate state vector to obtain a plurality of states at the second set of qubits in the quantum 5 circuit. Sequences of bitstrings are read from the second set of qubits (by block 412) to obtain the direct matrix. This step is mathematically represented by equation 18.
⟨0,𝑖,0,0,0|𝐻⊗𝑚𝑉𝐿′𝑅𝑉𝐿𝑂(𝐷)|0,𝑗,0,1,0⟩=Σ𝐿𝑖𝑗𝑎𝐿𝑘𝑙𝑎𝐷𝑘𝑙=𝐽𝑖𝑗 (18)
[053]
Once the direct matrix is obtained, at step 210b, a correlation exchange matrix is determined from the collocation matrix. Steps 210b1 through 10 210b6 provides an explanation for determining the correlation exchange matrix. The collocation matrix is encoded in a quantum circuit component as illustrated in FIG. 5, where the collocation matrix is encoded in terms of multi-qubit unitary gates in which each entry of the collocation matrix is encoded on a2 ancilla qubits with controls on NaO, Ng sets of qubits. Initially, at step 210b1, a gradient of collocation 15 matrix is encoded by processing a set of collocation matrices. These set of collocation matrices are encoded for a set of positions (x+/-delta,y,z),(x,y+/-delta,z), (x,y,z+/-delta) for the number of basis functions at all (x,y,z) positions. The set of collocation matrices are encoded using a first quantum circuit component as shown in FIG. 6, a second quantum circuit component as shown in FIG. 7, a third 20 quantum circuit component as shown in FIG. 8 on (a) control qubits composed of a first set of qubits and a second set of qubits, and (b) a plurality of ancilla qubits where the first set of qubits is associated with the number of points on the numerical grid, the second set of qubits is associated with the number of basis functions, and the plurality of ancilla qubits are used to store numerical entries for the set of 25 collocation matrices. This is mathematically represented as: log Ng + log Nao + 1+2 qubits, where log Ng represents the first set of qubits, log Nao represents the second set of qubits, and 1+2 qubits represent the plurality of ancilla qubits. In FIGS. 6,7, and 8 a1, a2, and a3 represents the plurality of ancilla qubits. FIG. 6
24
illustrates a first quantum circuit component
, where a first subset of subset of collocation matrices are encoded for a first set of grid points ((+delta x, 0,0) and (-delta x,0,0)) to obtain an X-component of the gradient of the collocation matrix. Equation 19 represents the collocation matrix,
Φ=( 𝜙1(𝑥1,𝑦1,𝑧1),………𝜙1(𝑥𝑁𝑔,𝑦𝑁𝑔,𝑧𝑁𝑔)⋮⋮𝜙𝑁𝑎𝑜(𝑥1,𝑦1,𝑧1),………𝜙𝑁𝑎𝑜(𝑥𝑁𝑔,𝑦𝑁𝑔,𝑧𝑁𝑔)) (19) 5
[054]
The first subset of collocation matrices are represented using equation 20,
Φ±Δ𝑥=( 𝜙1(𝑥1±Δ𝑥,𝑦1,𝑧1),………𝜙1(𝑥𝑁𝑔±Δ𝑥,𝑦𝑁𝑔,𝑧𝑁𝑔)⋮⋮𝜙𝑁𝑎𝑜(𝑥1±Δ𝑥,𝑦1,𝑧1),………𝜙𝑁𝑎𝑜(𝑥𝑁𝑔±Δ𝑥,𝑦𝑁𝑔,𝑧𝑁𝑔)) (20)
Further FIG. 7 illustrates a second quantum circuit component, where a second subset of collocation matrices is encoded for a second set of grid points ((0, +delta 10 y, 0) and (0, -delta y, 0)) to obtain a Y-component of the gradient of the collocation matrix. The second subset of collocation matrices are represented using equation 21,
Φ±Δ𝑦=( 𝜙1(𝑥1,𝑦1±Δ𝑦,𝑧1),………𝜙1(𝑥𝑁𝑔,𝑦𝑁𝑔+±Δ𝑦,𝑧𝑁𝑔)⋮⋮𝜙𝑁𝑎𝑜(𝑥1±Δ𝑥,𝑦1,𝑧1),………𝜙𝑁𝑎𝑜(𝑥𝑁𝑔,𝑦𝑁𝑔±Δ𝑦,𝑧𝑁𝑔)) (21)
15
[055]
Next FIG. 8 illustrates a third quantum circuit component, where a third subset of collocation matrices is encoded for a set of grid points ((0,0, +delta z) and (0,0, -delta z)) to obtain a Z-component of the gradient of the collocation matrix. The third subset of collocation matrices are represented using equation 22,
Φ±Δ𝑧=( 𝜙1(𝑥1,𝑦1±Δ𝑦,𝑧1),………𝜙1(𝑥𝑁𝑔,𝑦𝑁𝑔+±Δ𝑦,𝑧𝑁𝑔)⋮⋮𝜙𝑁𝑎𝑜(𝑥1,𝑦1,𝑧1±Δ𝑧),………𝜙𝑁𝑎𝑜(𝑥𝑁𝑔,𝑦𝑁𝑔±Δ𝑦,𝑧𝑁𝑔±Δ𝑧)) (22) 20
25
[056]
The X-component of the gradient of the collocation matrix is mathematically represented using equations 23.
𝑉∂𝑗Φ=𝐼2⊗(𝑛𝑎𝑜+𝑛𝑔)𝐻𝐼2⊗2𝐻Σ(|𝑐,𝑟,0,0⟩⟨𝑐,𝑟,0,0|⊗[Φ+Δ𝑗,𝑐,𝑟𝐼+𝑟,𝑐𝑖√1−Φ+Δ𝑗,𝑐,𝑟2𝑌] +|𝑐,𝑟,0,1⟩⟨𝑐,𝑟,0,1|⊗𝐼2+|𝑐,𝑟,1,0⟩⟨𝑐,𝑟,1,0|⊗[Φ−Δ𝑗,𝑐,𝑟𝐼+𝑖√1−Φ−Δ𝑗,𝑐,𝑟2𝑌]+5 |𝑐,𝑟,1,1⟩⟨𝑐,𝑟,1,1|⊗𝐼2)𝐼2⊗(𝑛𝑎𝑜+𝑛𝑔)𝐻𝐼2⊗2 (23)
[057]
Similarly, 𝑉∂𝑦Φ and 𝑉∂𝑧Φ (the Y-component of the gradient of the collocation matrix and the Z-component of the gradient of the collocation matrix) are obtained using the second and third quantum circuit component. Using these a quantum circuit for ∇𝜌 is mathematically represented as in equations 24, 25, and 10 26.
∂𝑥𝜌|𝑟𝑖=〈0,𝑖,0,0,0|𝑉Φ𝑅1𝑉𝐷𝑅1𝑉∂𝑥Φ|0,𝑗,0,0,0〉+〈0,𝑖,0,0,0|𝑉∂𝑥Φ𝑅1𝑉𝐷𝑅1𝑉Φ|0,𝑗,0,0,0〉=Σ[Φ𝑖𝑏𝐷𝑎𝑏(ΦΔ𝑋,𝑖𝑎−Φ−Δ𝑋,𝑖𝑎)+(ΦΔ𝑋,𝑖𝑏−Φ−Δ𝑋,𝑖𝑏)𝐷𝑎𝑏Φ𝑖𝑎] (24) 15
∂𝑦𝜌|𝑟𝑖=〈0,𝑖,0,0,0|𝑉Φ𝑅1𝑉𝐷𝑅1𝑉∂𝑦Φ|0,𝑗,0,0,0〉+〈0,𝑖,0,0,0|𝑉∂𝑦Φ𝑅1𝑉𝐷𝑅1𝑉Φ|0,𝑗,0,0,0〉=Σ[Φ𝑖𝑏𝐷𝑎𝑏(ΦΔ𝑌,𝑖𝑎−Φ−Δ𝑌,𝑖𝑎)+(ΦΔ𝑌,𝑖𝑏−Φ−Δ𝑌,𝑖𝑏)𝐷𝑎𝑏Φ𝑖𝑎] (25)
∂𝑧𝜌|𝑟𝑖=〈0,𝑖,0,0,0|𝑉Φ𝑅1𝑉𝐷𝑅1𝑉∂𝑥Φ|0,𝑗,0,0,0〉+〈0,𝑖,0,0,0|𝑉∂𝑥Φ𝑅1𝑉𝐷𝑅1𝑉Φ|0,𝑗,0,0,0〉=ΣΦ𝑖𝑏𝐷𝑎𝑏(ΦΔ𝑍,𝑖𝑎−Φ−Δ𝑍,𝑖𝑎)+20 (ΦΔ𝑍,𝑖𝑏−Φ−Δ𝑍,𝑖𝑏)𝐷𝑎𝑏Φ𝑖𝑎] (26)
Then the first quantum circuit component, the second quantum circuit component, and the third quantum circuit component are composed to obtain the gradient of the collocation matrix.
[058]
Further, at step 210b2, an electronic density is encoded by 25 sequentially performing a matrix multiplication of the collocation matrix, the density matrix, and the collocation matrix which is illustrated in a quantum circuit
26
component represented as
Urho. In FIG. 9, C represents unitary rotation matrix obtained from diagonalizing Fock matrix, M symbolizes the sequence of controlled unitary gates acting on the r2 that encodes the entries of the D matrix on the qubit a2, D symbolizes the occupancy matrix whose first Ne/2 entries are one and the rest is zero, C† is the Hermitian conjugate of C, and Φ =symbolizes the collocation 5 matrix
[059]
Next at step 210b3, an electronic density gradient is computed from the collocation matrix, the density matrix, and the gradient of the collocation matrix by composing a fourth quantum circuit component and a fifth quantum circuit component using at least one ancilla qubit from the plurality of ancilla qubits. The 10 fourth quantum circuit component and a fifth quantum circuit component are shown in FIG. 10 and 11 respectively. The fourth quantum circuit component sequentially encodes the collocation matrix, the density matrix, and the gradient of the collocation matrix. The density matrix is encoded on control qubits composed of the second set of qubits and at least one ancilla qubits from the plurality of ancilla 15 qubits. Equation 27 mathematically represents the encoding done using the fourth quantum circuit component,
𝑉1=|00⟩⟨00|⊗𝑉𝜕𝑥Φ𝑅1𝑉𝐷𝑅1𝑉Φ+|01⟩⟨01|⊗𝑉𝜕𝑦Φ𝑅1𝑉𝐷𝑅1𝑉Φ+|10⟩⟨10|⊗𝑉𝜕𝑧Φ𝑅1𝑉𝐷𝑅1𝑉Φ (27)
The fifth quantum circuit component sequentially encodes the gradient of the 20 collocation matrix, the density matrix, and the collocation matrix. Equation 28 mathematically represents the encoding done using the fifth quantum circuit component,
𝑉2=|00⟩⟨00|⊗𝑉𝜕𝑥Φ𝑅1𝑉𝐷𝑅1𝑉Φ+|01⟩⟨01|⊗𝑉𝜕𝑦Φ𝑅1𝑉𝐷𝑅1𝑉Φ+|10⟩⟨10|⊗𝑉𝜕𝑧Φ𝑅1𝑉𝐷𝑅1𝑉Φ+|11⟩⟨11|⊗𝐼 (28) 25
[060]
Finally, the electronic density gradient is computed by composing the fourth quantum circuit component and the fifth quantum circuit component which is shown in FIG. 12. Equation 29 mathematically represents the electronic density gradient,
𝑉∇𝜌=|0⟩⟨0|⊗𝑉1+|0⟩⟨0|⊗𝑉2 (29) 30
27
[061]
Further at step 210b4, a derivative of an electronic energy density is computed with respect to the electronic density by performing a quantum signal processing sequence on a sixth quantum circuit as shown inside the dashed line box of FIG. 13. The sixth quantum circuit component encodes the electronic density and the electronic density gradient as shown in FIG. 13. In FIG. 13, Urho represents 5 the quantum circuit component for encoding the electronic density as shown in FIG. 9. The quantum signal processing is performed on the sixth quantum circuit component that encodes the electronic density and its gradient, a nonlinear function. The nonlinear function represents the derivative of the electronic energy density with respect to the electronic density. The quantum signal processing sequence is 10 made from a sequence of phase shifted reflectors and the sixth quantum circuit component encoding the electronic density and the electronic density gradient. The phase angles for the phase shifted reflectors are obtained from the Chebyshev expansion of the nonlinear function using a classical computer. These angles can be predetermined and are determined only once. S(phi1),…., S(phik) in FIG. 13 15 represents the phase angles.
[062]
At step 210b5 a Z-matrix is encoded by composing the sixth quantum circuit component and the collocation matrix as illustrated in FIG. 12. Equation 30 below represents the encoding of the Z-matrix,
𝑉𝑍=Π[(𝐻⊗𝑛𝑎𝑜𝑁𝑘=120 ⊗𝐼⊗(𝑛𝑔+2))𝑉Φ𝑅1𝑉𝐷𝑅1𝑉Φ((1+𝑒−𝑖𝜙𝑘)|0,𝑖,0,0⟩⟨0,𝑖,0,0|−1)] Π[(𝐻⊗𝑛𝑎𝑜⊗𝐼⊗(𝑛𝑔+2))𝑉∇Φ𝑅1𝑉𝐷𝑅1𝑉Φ((1+𝑒−𝑖𝜙𝑘)|0,𝑖,0,0⟩⟨0,𝑖,0,0|−1)]𝑁𝑘=1
(𝐻⊗𝑛𝑎𝑜⊗𝐼⊗(𝑛𝑔+2))𝑉Φ𝐻⊗𝑛𝑔 (30)
[063]
Finally at step 210b6, a correlation exchange matrix is encoded by composing the Z-matrix with the collocation matrix. This is represented by equation 25 31 below,
𝑉𝑖𝑗𝑥𝑐=〈𝑖,0,0,0|𝐻⊗𝑛𝑔𝑉Φ𝑉𝑍|𝑗,0,0,0〉=ΣΦ𝑖𝑟𝑟𝑑𝐸(𝜌,∇𝜌)𝑑𝜌(𝑟)Φ𝑗𝑟 (31)
28
[064]
In an embodiment, theorem A is applied on equations 24 through 26 and equations 29 through 31 to implement matrix multiplications and tensor contractions on quantum circuit components with qubit count efficient resources. The theorem A states that: If 𝐴 and 𝐵 are general rectangular matrices of dimensions 𝑑𝑖𝑚(𝐴) = (𝑁,𝑃) and 𝑑𝑖𝑚(𝐵) = (𝑃,𝑀) then there is a unitary 5 operation 𝑈(𝐴,𝐵) of dimension 2𝑛𝑞× 2𝑛𝑞 that operates on a system of 𝑛𝑞=𝑝+max(𝑚,𝑛)+2 qubit registers |⋅⟩𝑝|⋅⟩max(𝑚,𝑛)|⋅⟩𝑎1|⋅⟩𝑎2 (where 𝑛=⌈𝑙𝑜𝑔2𝑁⌉,𝑚=⌈𝑙𝑜𝑔2𝑀⌉,𝑝=⌈𝑙𝑜𝑔2𝑃⌉ )and block encodes the matrix multiplication of 𝐴 and 𝐵 such that:
⟨0|𝑝⟨𝑖|max(𝑚,𝑛)⟨0|𝑎1⟨0|𝑎2𝑈(𝐴,𝐵)|0⟩𝑝|𝑗⟩max(𝑚,𝑛)|1⟩𝑎1|0⟩𝑎2= 10 1𝑚𝑎𝑥(𝑀,𝑁)𝑃|Σ𝐴𝑖𝑘𝑘𝐵𝑘𝑗||𝐴||||𝐵|| (32)
Proof of the theorem A is given below:
Take the normalized matrices 𝐴′=𝐴/(√2||𝐴||), 𝐵′=𝐵/(√2||𝐵||). For these define two unitary operators 𝑉(𝐴),𝑉(𝐵), 15
𝑉𝐴=Σ[|𝑐,𝑟,0⟩⟨𝑐,𝑟,0|⊗ (𝐴𝑟𝑐′𝐼+𝑖√1−𝐴𝑟𝑐′2𝑌)+ |𝑐,𝑟,1⟩⟨𝑐,𝑟,1|⊗2𝑝,2max(𝑚,𝑛)𝑐=0,𝑟=0 𝐼2]
𝑉𝐵=Σ[|𝑐,𝑟,0⟩⟨𝑐,𝑟,0|⊗ 𝐼2+ |𝑐,𝑟,1⟩⟨𝑐,𝑟,1|⊗ (𝐵𝑐𝑟′𝐼+2𝑝,2max(𝑚,𝑛)𝑐=0,𝑟=0𝑖√1−𝐵𝑐𝑟′2𝑌) ]
(33) 20
The classical data of the B matrix is loaded using the state preparation oracle 𝑉𝐵𝐻⊗𝑝 on the initial state |0⟩|𝑗⟩|1⟩|0⟩,
|Φ𝐵⟩= 𝑉𝐵𝐻⊗𝑝|0⟩|𝑗⟩|1⟩|0⟩ = 1√𝑃Σ[𝐵𝑐𝑗′|𝑐,𝑗,1,0⟩+√1−𝐵𝑐𝑗′2|𝑐,𝑗,1,1⟩]𝑐 (34)
The classical data of the A matrix is loaded using the state preparation oracle 𝑉𝐴𝐻⊗𝑝, 25
|Φ𝐴⟩= 𝑉𝐴𝐻⊗𝑝|0⟩|𝑖⟩|0⟩|0⟩ =1√𝑃Σ[𝐴𝑖𝑐′|𝑐,𝑖,0,0⟩ +√1−𝐴𝑖𝑐′2|𝑐,𝑖,0,1⟩]𝑐 (35)
29
The states |Φ𝐴⟩ and |Φ𝐵⟩ are orthogonal,
⟨Φ𝐴|Φ𝐵⟩=0 (36)
Next define diffusion operator $R$ acting on the row registers and the ancillas 𝑎1, 𝑎2,
𝑅=𝐼2⊗𝑝⊗[𝐻⊗max(𝑚,𝑛)⊗𝐻⊗𝐼2)(2|0,0,0⟩⟨0,0,0|−1)𝐻⊗max(𝑚,𝑛)⊗𝐻⊗5 𝐼2)]
=𝐼2⊗𝑝⊗[Σ2|𝑘,+,0⟩⟨𝑙,+,0|𝑚𝑎𝑥(𝑀,𝑁)−𝐼𝑘,𝑙] (37)
Then the overlap between these two states |Φ𝐴⟩ and 𝑅|Φ𝐵⟩ is given by,
⟨Φ𝐴|𝑅|Φ𝐵⟩=⟨0,𝑖,0,0|𝐻𝑐⊗𝑝𝑉𝐴†𝑅𝑉𝐵𝐻𝑐⊗𝑝|0,𝑗,1,0⟩=2Σ𝐴𝑖𝑐′𝑐𝐵𝑐𝑗′𝑚𝑎𝑥(𝑀,𝑁)𝑃=Σ𝐴𝑖𝑐𝑐𝐵𝑐𝑗𝑚𝑎𝑥(𝑀,𝑁)𝑃||𝐴||||𝐵|| 10 (38)
By construction it is proved the existence of $U(A,B)$ that can be defined without any isometry,
𝑈(𝐴,𝐵)=𝐻𝑐⊗𝑝𝑉𝐴†𝑅𝑉𝐵𝐻𝑐⊗𝑝 (39)
[065]
Once the correlation exchange matrix is determined, at step 210c, a 15 Fock matrix is computed by adding the core Hamiltonian matrix, the density matrix, and the correlation exchange matrix. At step 210d, qubitized diagonalization of the Fock matrix is performed to obtain an updated density matrix. The updated density matrix is used in a subsequent iteration. Step 210 is performed iteratively until the convergence criteria is satisfied to get a final density matrix. The convergence 20 criteria is satisfied when norm of difference between the density matrix of a current iteration and the density matrix of a previous iteration is lesser than a pre-defined tolerance value. Finally, at step 212, the one or more classical hardware processors are configured to utilize the final density matrix of the chemical compound to extract the one or more properties of the chemical compound. The final density 25 matrix thus obtained enables computation of one or more properties such as, (i) HOMO-LUMO gap from the energies of the highest occupied Kohn Sham orbital and lowest unoccupied orbital, (ii) the charge distribution from integrating the electronic density in regions around the different atoms, where the electronic
30
density
is in turn obtained from the density matrix, and (iii) dipole moment from the expectation value of the perturbation electric field operator with respect to the density matrix.
EXPERIMENTAL RESULTS
[066]
The J matrix computation was performed on classical processor 5 using conventional methods and quantum processor. Also, the DFT was implemented through GGA using classical processors and quantum processors. Classical complexity is measured in terms of space and time complexity. Similarly, quantum complexity is measured in terms of number of qubits and gate complexity. The results of J-matrix computations and correlation exchange matrix computations 10 are recorded in table 1.
Quantity
Classical Complexity
Quantum Complexity
Space
Time
No.of Qubits
Gate complexity
J(AO)
𝑂(𝑁𝑏4)
𝑂(𝑁𝑏4)
4log𝑁𝑏+2
8log𝑁𝑏
J(DF)
𝑂(𝑁𝑏2𝑁𝑎𝑢𝑥)
𝑂(𝑁𝑏2𝑁𝑎𝑢𝑥)
2log𝑁𝑏+log𝑁𝑎𝑢𝑥
4log𝑁𝑏+2log𝑁𝑎𝑢𝑥
J(sN)
𝑂(𝑁𝑏2𝑁𝑔)
𝑂(𝑁𝑏2𝑁𝑔)
2log𝑁𝑏+log𝑁𝑔
4log𝑁𝑏+2log𝑁𝑔
𝑉𝑥𝑐[𝜌,∇𝜌](𝐺𝐺𝐴)
𝑂(𝑁2𝑁𝑔)
𝑂(𝑁2𝑁𝑔)
log𝑁+𝑙𝑜𝑔𝑁𝑔
4(log𝑁+𝑙𝑜𝑔𝑁𝑔)𝑁𝑏𝑁𝑔
Table 1
[067]
The classical and quantum resources needed to compute the J matrix using the AO integral approach, the Density fitting (DF) approach and the semi-numerical exchange (sN) approach are shown in Table 1. For all these three types 15 of methods to compute J matrix the Quantum circuits of same type as given in the present disclosure is applicable. For the AO integral approach the number of qubits needed to encode the electronic integrals is 4log Nb+2, where Nb=Nao is the number of basis functions. On the other hand for the classical approach the space complexity is 20 𝑂(𝑁𝑏4) to store 𝑁𝑏4 numbers. Log is clearly a sublinear memory scaling on the other hand the polynomial scaling involved is steep in classical approach. The gate complexity needed for computing the J matrix stems from two sources, one is the diffusion quantum circuit block and the multiqubit gates needed to encode the
31
integrals and the density matrix. The first component from the diffusion quantum
circuit block leads to the gate complexity 8log𝑁𝑏. The second component depends on the molecular system, and the gate complexity for that component depends on the floating-point precision with which we are encoding the integrals associated with the molecule. In comparison the classical time complexity scales as 𝑂(𝑁𝑏4). 5
[068]
For the Density fitting approach the number of qubits needed to encode the electronic integrals is 2log𝑁𝑏+log𝑁𝑎𝑢𝑥 , where 𝑁𝑏=𝑁𝑎𝑜 is the number of basis functions and 𝑁𝑎𝑢𝑥 number of auxilliarly basis vectors. On the other hand for the classical approach the space complexity is 𝑂(𝑁𝑏2𝑁𝑎𝑢𝑥) to store that many numbers. Log is clearly a sublinear memory scaling on the other hand 10 the polynomial scaling involved is steep in classical approach. The gate complexity needed for computing the J matrix stems from two sources, one is the diffusion quantum circuit block and the multiqubit gates needed to encode the integrals and the density matrix. The first component from the diffusion quantum circuit block leads to the gate complexity 4log𝑁𝑏+2log𝑁𝑎𝑢𝑥. The second component 15 depends on the molecular system, and the gate complexity for that component depends on the floating-point precision with which we are encoding the integrals associated with the molecule. In comparison the base classical time complexity scales as 𝑂(𝑁𝑏2𝑁𝑎𝑢𝑥).
[069]
For the semi numerical engine approach the number of qubits 20 needed to encode the electronic integrals is 2log𝑁𝑏+log𝑁𝑔, where 𝑁𝑏=𝑁𝑎𝑜 is the number of basis functions and 𝑁𝑔 number of grid points. On the other hand for the classical approach the space complexity is 𝑂(𝑁𝑏2𝑁𝑔)to store that many numbers. Log is clearly a sublinear memory scaling on the other hand the polynomial scaling involved is steep in classical approach. The gate complexity needed for computing 25 the J matrix stems from two sources, one is the diffusion quantum circuit block and the multiqubit gates needed to encode the integrals and the density matrix. The first component from the diffusion quantum circuit block leads to the gate complexity 4log𝑁𝑏+2log𝑁𝑔. The second component depends on the molecular system, and the gate complexity for that component depends on the floating-point precision with 30
32
which we are encoding the integrals associated with the molecule. In comparison
the base classical time complexity scales as 𝑂(𝑁𝑏2𝑁𝑔).
[070]
For the GGA protocol the correlation exchange matrix is a function of the electronic density and its gradient. For computing the electronic density, the number of qubits needed to encode the electronic density is given by 𝑙𝑜𝑔 𝑁 +5 𝑙𝑜𝑔 𝑁𝑔, where 𝑁𝑔 is the number of points on numerical grid and 𝑁 is number of basis functions. Additionally, 2 ancilla qubits are needed to encode the different components of the electronic density and its gradient. The gate complexity is associated with the diffusion operator that enables the multiplication between the collocation matrix, the gradient of the collocation matrix and the density matrix. 10 And it scales logarithmically with the number of basis functions. Classically to store the integrals is 𝑂(𝑁2𝑁𝑔), and the time complexity is 𝑂(𝑁2𝑁𝑔).
[071]
The written description describes the subject matter herein to enable any person skilled in the art to make and use the embodiments. The scope of the subject matter embodiments is defined by the claims and may include other 15 modifications that occur to those skilled in the art. Such other modifications are intended to be within the scope of the claims if they have similar elements that do not differ from the literal language of the claims or if they include equivalent elements with insubstantial differences from the literal language of the claims.
[072]
It is to be understood that the scope of the protection is extended to 20 such a program and in addition to a computer-readable means having a message therein; such computer-readable storage means contain program-code means for implementation of one or more steps of the method, when the program runs on a server or mobile device or any suitable programmable device. The hardware device can be any kind of device which can be programmed including e.g., any kind of 25 computer like a server or a personal computer, or the like, or any combination thereof. The device may also include means which could be e.g., hardware means like e.g., an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA), or a combination of hardware and software means, e.g., an ASIC and an FPGA, or at least one microprocessor and at least one memory with 30 software processing components located therein. Thus, the means can include both
33
hardware means and software means. The method embodiments described herein
could be implemented in hardware and software. The device may also include software means. Alternatively, the embodiments may be implemented on different hardware devices, e.g., using a plurality of CPUs.
[073]
The embodiments herein can comprise hardware and software 5 elements. The embodiments that are implemented in software include but are not limited to, firmware, resident software, microcode, etc. The functions performed by various components described herein may be implemented in other components or combinations of other components. For the purposes of this description, a computer-usable or computer readable medium can be any apparatus that can 10 comprise, store, communicate, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device.
[074]
The illustrated steps are set out to explain the exemplary embodiments shown, and it should be anticipated that ongoing technological development will change the manner in which particular functions are performed. 15 These examples are presented herein for purposes of illustration, and not limitation. Further, the boundaries of the functional building blocks have been arbitrarily defined herein for the convenience of the description. Alternative boundaries can be defined so long as the specified functions and relationships thereof are appropriately performed. Alternatives (including equivalents, extensions, 20 variations, deviations, etc., of those described herein) will be apparent to persons skilled in the relevant art(s) based on the teachings contained herein. Such alternatives fall within the scope of the disclosed embodiments. Also, the words “comprising,” “having,” “containing,” and “including,” and other similar forms are intended to be equivalent in meaning and be open ended in that an item or items 25 following any one of these words is not meant to be an exhaustive listing of such item or items, or meant to be limited to only the listed item or items. It must also be noted that as used herein and in the appended claims, the singular forms “a,” “an,” and “the” include plural references unless the context clearly dictates otherwise.
[075]
Furthermore, one or more computer-readable storage media may be 30 utilized in implementing embodiments consistent with the present disclosure. A
34
computer
-readable storage medium refers to any type of physical memory on which information or data readable by a processor may be stored. Thus, a computer-readable storage medium may store instructions for execution by one or more processors, including instructions for causing the processor(s) to perform steps or stages consistent with the embodiments described herein. The term “computer-5 readable medium” should be understood to include tangible items and exclude carrier waves and transient signals, i.e., be non-transitory. Examples include random access memory (RAM), read-only memory (ROM), volatile memory, nonvolatile memory, hard drives, CD ROMs, DVDs, flash drives, disks, and any other known physical storage media. 10
[076]
It is intended that the disclosure and examples be considered as exemplary only, with a true scope of disclosed embodiments being indicated by the following claims.
We Claim:
1. A quantum simulation method (200), performed by a system comprising one or more classical hardware processors and a plurality of unentangled Quantum Processor Units (QPUs), wherein the one or more classical hardware processors are communicably coupled to the plurality of unentangled QPUs by communication interfaces, wherein the quantum simulation method comprising:
receiving, via the one or more classical hardware processors, a chemical compound whose one or more properties are to be extracted, wherein the chemical compound is one of (i) a molecule, and (ii) a solid (202);
obtaining, via the one or more classical hardware processors, atomic coordinates of each of a plurality of atoms present in the chemical compound (204);
determining, via the one or more classical hardware processors, a plurality of electron integrals, a core Hamiltonian matrix, and a collocation matrix from the atomic coordinates of each of the plurality of atoms, wherein the collocation matrix is a rectangular matrix of dimension Ng and NaO, and wherein Ng represents a number of points on a numerical grid, and NaO represents a number of basis functions (206);
determining, via the plurality of unentangled QPUs, a density matrix of the chemical compound from the core Hamiltonian matrix wherein the density matrix is constructed using (i) a single particles unitary rotation matrix, obtained by diagonalizing the core Hamiltonian matrix, and (ii) electron occupancies (208); and
iteratively updating, via the plurality of unentangled QPUs, the density matrix by computing a Fock matrix until a convergence criteria is met, wherein iteratively updating the density matrix at each iteration comprises (210):
computing a direct (J) matrix from the density matrix based
on a subset of electron integrals (210a);

determining a correlation exchange matrix from the collocation matrix for generalized gradient approximation by (210b),
encoding a gradient of the collocation matrix by processing a set of collocation matrices, wherein the set of collocation matrices are encoded for a set of positions (x+/-delta,y,z),(x,y+/-delta,z), (x,y,z+/-delta) for the number of basis functions at all (x,y,z) positions using a first quantum circuit component, a second quantum circuit component, a third quantum circuit component on (a) control qubits composed of a first set of qubits and a second set of qubits, and (b) a plurality of ancilla qubits wherein the first set of qubits is associated with the number of points on the numerical grid, the second set of qubits is associated with the number of basis functions, and the plurality of ancilla qubits are used to store numerical entries for the set of collocation matrices (210b1);
encoding an electronic density by performing a sequential matrix multiplication of the collocation matrix, the density matrix, and the collocation matrix (210b2);
computing an electronic density gradient from the collocation matrix, the density matrix, and the gradient of the collocation matrix by composing a fourth quantum circuit component and a fifth quantum circuit component using at least one ancilla qubit from the plurality of ancilla qubits (210b3);
computing a derivative of an electronic energy density with respect to the electronic density by performing a quantum signal processing sequence on a sixth quantum circuit, wherein the sixth quantum circuit component

encodes the electronic density and the electronic density gradient (210b4);
encoding a Z-matrix by composing the sixth quantum circuit component and the collocation matrix (210b5); and
encoding a correlation exchange matrix by
composing the Z-matrix with the collocation matrix (210b6);
determining the Fock matrix by combining the core
Hamiltonian matrix, the density matrix, and the correlation
exchange matrix (210c); and
updating the density matrix by diagonalizing the Fock matrix (210d); and utilizing, via the one or more classical hardware processors, the updated density matrix of the chemical compound, to extract the one or more properties of the chemical compound (212).
2. The method as claimed in claim 1, wherein encoding the gradient of the collocation matrix comprises:
encoding, via the plurality of unentangled QPUs, a first subset of collocation matrices for a first set of grid points ((+delta x, 0,0) and (-delta x,0,0)) to obtain an X-component of the gradient of the collocation matrix using a first quantum circuit component;
encoding, via the plurality of unentangled QPUs, a second subset of collocation matrices for a second set of grid points ((0, +delta y, 0) and (0, -delta y, 0)) to obtain a Y-component of the gradient of the collocation matrix using a second quantum circuit component;
encoding, via the plurality of unentangled QPUs, a third subset of collocation matrices for a set of grid points ((0,0, +delta z) and (0,0, -delta z)) to obtain a Z-component of the gradient of the collocation matrix using a third quantum circuit component; and

composing, via the plurality of unentangled QPUs, the first quantum circuit component, the second quantum circuit component, and the third quantum circuit component to obtain the gradient of the collocation matrix.
3. The method as claimed in claim 1, wherein,
the fourth quantum circuit component sequentially encodes the collocation matrix, the density matrix, and the gradient of the collocation matrix, wherein the density matrix is encoded on control qubits composed of the second set of qubits and at least one ancilla qubits from the plurality of ancilla qubits, and
the fifth quantum circuit component sequentially encodes the gradient of the collocation matrix, the density matrix, and the collocation matrix.
4. The method as claimed in claim 1, wherein, the (i) control qubits composed
of the first set of qubits and the second set of qubits, and (ii) the plurality of
ancilla qubits, is mathematically represented as:
log Ng + log Nao + 1+2 qubits wherein log Ng represents the first set of qubits, log Nao represents the second set of qubits, and 1+2 qubits represent the plurality of ancilla qubits.
5. The method as claimed in claim 1, wherein the convergence criteria is satisfied when norm of difference between the density matrix of a current iteration and the density matrix of a previous iteration is lesser than a pre-defined tolerance value.
6. A system (100) comprising:
one or more classical hardware processors (108) and a plurality of unentangled Quantum Processor Units (QPUs) (122), wherein the one or more classical hardware processors (108) are communicably coupled to the

plurality of unentangled QPUs (122) by one or more communication interfaces (106), wherein the one or more classical hardware processors are operatively coupled to at least one memory (110) storing programmed instructions and one or more Input/Output (I/O) interfaces (116); and the plurality of unentangled quantum processors (122) are operatively coupled to the at least one quantum memory (124), wherein the one or more hardware processors (108) and the plurality of unentangled QPUs (122) are configured by the programmed instructions to:
receive a chemical compound whose one or more properties are to be extracted, wherein the chemical compound is one of (i) a molecule, and (ii) a solid;
obtain atomic coordinates of each of a plurality of atoms present in the chemical compound;
determine a plurality of electron integrals, a core Hamiltonian matrix, and a collocation matrix from the atomic coordinates of each of the plurality of atoms, wherein the collocation matrix is a rectangular matrix of dimension Ng and NaO, and wherein Ng represents a number of points on a numerical grid, and NaO represents a number of basis functions;
determine a density matrix of the chemical compound from the core Hamiltonian matrix wherein the density matrix is constructed using (i) a single particles unitary rotation matrix, obtained by diagonalizing the core Hamiltonian matrix, and (ii) electron occupancies; and
iteratively update the density matrix by computing a Fock matrix until a convergence criteria is met, wherein iteratively updating the density matrix at each iteration comprises:
compute a direct (J) matrix from the density matrix based on
a subset of electron integrals;
determine a correlation exchange matrix from the
collocation matrix for generalized gradient approximation by,
encoding a gradient of the collocation matrix by processing a set of collocation matrices, wherein the set of

collocation matrices are encoded for a set of positions (x+/-delta,y,z),(x,y+/-delta,z), (x,y,z+/-delta) for the number of basis functions at all (x,y,z) positions using a first quantum circuit component, a second quantum circuit component, a third quantum circuit component on (a) control qubits composed of a first set of qubits and a second set of qubits, and (b) a plurality of ancilla qubits wherein the first set of qubits is associated with the number of points on the numerical grid, the second set of qubits is associated with the number of basis functions, and the plurality of ancilla qubits are used to store numerical entries for the set of collocation matrices;
encoding an electronic density by performing a sequential matrix multiplication of the collocation matrix, the density matrix, and the collocation matrix;
computing an electronic density gradient from the collocation matrix, the density matrix, and the gradient of the collocation matrix by composing a fourth quantum circuit component and a fifth quantum circuit component using at least one ancilla qubit from the plurality of ancilla qubits;
computing a derivative of an electronic energy density with respect to the electronic density by performing a quantum signal processing sequence on a sixth quantum circuit, wherein the sixth quantum circuit component encodes the electronic density and the electronic density gradient;
encoding a Z-matrix by composing the sixth quantum circuit component and the collocation matrix; and
encoding a correlation exchange matrix by composing the Z-matrix with the collocation matrix;

determining the Fock matrix by combining the core Hamiltonian matrix, the density matrix, and the correlation exchange matrix; and
updating the density matrix by diagonalizing the Fock matrix; and utilizing, via the one or more classical hardware processors, the updated density matrix of the chemical compound, to extract the one or more properties of the chemical compound.
7. The system of claim 6, wherein encoding the gradient of the collocation
matrix comprises:
encoding a first subset of collocation matrices for a first set of grid points ((+delta x, 0,0) and (-delta x,0,0)) to obtain an X-component of the gradient of the collocation matrix using a first quantum circuit component;
encoding a second subset of collocation matrices for a second set of grid points ((0, +delta y, 0) and (0, -delta y, 0)) to obtain a Y-component of the gradient of the collocation matrix using a second quantum circuit component;
encoding a third subset of collocation matrices for a set of grid points ((0,0, +delta z) and (0,0, -delta z)) to obtain a Z-component of the gradient of the collocation matrix using a third quantum circuit component; and
composing the first quantum circuit component, the second quantum circuit component, and the third quantum circuit component to obtain the gradient of the collocation matrix.
8. The system of claim 6, wherein,
the fourth quantum circuit component sequentially encodes the collocation matrix, the density matrix, and the gradient of the collocation matrix, wherein the density matrix is encoded on control qubits composed of the second set of qubits and at least one ancilla qubits from the plurality of ancilla qubits, and

the fifth quantum circuit component sequentially encodes the gradient of the collocation matrix, the density matrix, and the collocation matrix.
9. The system of claim 6, wherein, the (i) control qubits composed of the first
set of qubits and the second set of qubits, and (ii) the plurality of ancilla
qubits, is mathematically represented as:
log Ng + log Nao + 1+2 qubits wherein log Ng represents the first set of qubits, log Nao represents the second set of qubits, and 1+2 qubits represent the plurality of ancilla qubits.
10. The system of claim 6, wherein the convergence criteria is satisfied as norm
of difference between the density matrix of a current iteration and the
density matrix of a previous iteration is lesser than a pre-defined tolerance
value.