Method and system for quantum computing-enabled molecular first-principles simulations - Patents.com

Abstract

The present disclosure provides methods and systems for using hybrid classical and non-classical (e.g., quantum) computing architectures to calculate the quantum mechanical energy and / or electronic structure of chemical systems, as well as to identify stable conformations of chemical systems (e.g., molecules) and / or to perform first-principles molecular dynamics calculations or simulations on chemical systems.

Description

[Technical field]

[0001] cross reference This application claims the benefit of U.S. Provisional Application No. 62 / 949,263, filed December 17, 2019, and U.S. Provisional Application No. 62 / 847,141, filed May 13, 2019, each of which is incorporated by reference in its entirety for all purposes. [Background technology]

[0002] In chemistry and biology, identifying and predicting the electronic structure and most energetically stable conformer of a molecule is of great importance since molecular function is essentially embedded in the molecular conformation. For example, the reaction rate in a catalytic reaction can vary significantly based on which of several different conformations of the catalyst is used. As another example, a protein is more functional or functional at all when it forms a particular tertiary structure.

[0003] To accurately identify and predict the electronic structure and most stable conformers, highly accurate quantum chemical methods such as Coupled-Cluster (CC) theory or Full Configuration Interaction (Full CI) can be implemented. However, the computational cost of such methods can grow exponentially with the size of the molecule, and they often become intractable when the size of the molecule exceeds about 50 atoms for CC and about 10 atoms for Full CI, even when implemented on some current state-of-the-art classical computers. Therefore, highly efficient and accurate computational frameworks are needed to identify the most stable conformers of industrially relevant chemical compounds and biologically relevant large molecules.

[0004] Quantum computing (QC) techniques may be capable of calculating the quantum mechanical energy and / or electronic structure of molecules with exponentially less computational resources compared to classical computing. Thus, high-precision quantum chemical calculations that are intractable using classical computing may become tractable using QC techniques. However, QC techniques may face challenges such as high cost and scarcity of QC resources. In addition, it is technically difficult to increase the number of qubits in a quantum computer, where the size of the quantum computing device is limited. In addition, qubits are highly sensitive to noise and environmental effects, which can cause qubits to decoherence in a very short time, thus providing a relatively small window to perform meaningful calculations. Summary of the Invention [Problem to be solved by the invention]

[0005] Recognized herein is the need for quantum algorithms and circuits that efficiently leverage current near-term quantum computing systems to solve complex quantum chemical problems. One approach is to decompose an industrial-sized problem into sub-problems, identify the more complex sub-problems, and then use a quantum computer to process a subset of the problem, e.g., the sub-problems that are difficult for classical computers. [Means for solving the problem]

[0006] The systems and methods provided herein utilize problem decomposition (PD) techniques in quantum chemistry to identify and predict the electronic structure and most energetically stable set of conformers of a molecule. Such PD techniques may include Fragment Molecular Orbital (FMO) methods, Divide-and-Conquer (DC) methods, Density Matrix Embedding Theory (DMET) methods, Density Matrix Renormalization Group (DMRG) methods, tensor networks, incremental methods, and the like, as described herein.

[0007] In quantum chemistry, PD techniques have been developed to efficiently calculate molecular energy and / or electronic structure with reasonable accuracy using classical computing. In PD techniques, molecules can be broken down into smaller fragments, so that quantum mechanical energy and / or electronic structure calculations are tractable for each fragment. Quantum mechanical energy and / or electronic structure calculations can then be performed separately for each fragment. The quantum mechanical energy and / or electronic structure calculations resulting from each fragment can be recombined into a solution for the original molecule.

[0008] The systems and methods provided herein for implementing PD techniques on a QC platform may enable quantum mechanical energy and / or electronic structure calculations to be performed with a high level of accuracy for each fragment. Furthermore, the small size of each fragment may enable highly accurate calculations to be performed on QC devices where the scale of calculations is rather limited, thus efficiently and accurately obtaining the energies and / or electronic structures of complex, industrially relevant molecules.

[0009] Identifying the electronic structure and most energetically stable conformer of a molecule is a fundamental process in chemical and biological research and development. Although such a process can be carried out by actually synthesizing the molecule and identifying its electronic structure and conformation using various physicochemical measurements, such an experimental process can require a very large amount of resources, such as human effort and time. Thus, highly efficient and accurate computational methods and systems such as those provided by the present disclosure can significantly reduce the need for such resources and make the general R&D process more efficient. Furthermore, the methods and systems described herein can be applied not only to single chemical system structures (e.g., chemical compounds and biomolecules), but also to molecular assemblies with different associations. For example, the methods and systems disclosed herein can be applied to identifying the most stable binding orientation of a drug candidate to a target protein, determined from a group of potential binding orientations.

[0010] The present disclosure provides methods and systems for efficiently identifying electronic structures and stable conformations of chemical systems (e.g., molecules) using a hybrid architecture of quantum and classical computing processors. The method may include obtaining an index of the molecule, calculating or obtaining a set of conformations of the molecule, and decomposing the chemical system into fragments (subsystems) for each conformation (which may optionally be stored in a list). The method may further include calculating a fermion Hamiltonian (molecular Hamiltonian or electronic Hamiltonian) for each fragment of each conformation of the molecule, converting each fermion Hamiltonian to an equivalent qubit Hamiltonian, converting the qubit Hamiltonian to a quantum circuit, calculating initial states of the qubits involved in the calculation of the total quantum mechanical energy and / or electronic structure, generating (e.g., through computational simulation) molecular quantum mechanical energy and / or electronic structure on quantum hardware or a classical simulator of the quantum circuit, and combining the energies and / or electronic structures of the multiple fragments to obtain an estimate of the total energy of the chemical system. The method may further include repeating these operations for all conformations in the group of conformations and classifying the conformations in the group of conformations based on the estimated total quantum mechanical energy and / or electronic structure. The method may further include providing an index of the classified conformational group (e.g., in a list).

[0011] In one aspect, the disclosure provides a method for performing a quantum mechanical energy or electronic structure calculation on a chemical system, the method being performed by a hybrid computing unit comprising a classical computer and a distributed computing system comprising a plurality of non-classical computers, the method including: (a) decomposing at least one conformation in a group of conformations of the chemical system into a plurality of molecular fragments; (b) using the hybrid computing unit, determining a quantum mechanical energy or electronic structure of each of at least a subset of the plurality of molecular fragments; (c) combining the quantum mechanical energies or electronic structures determined in (b); and (d) electronically outputting a report indicative of the quantum mechanical energies or electronic structures combined in (c).

[0012] In some embodiments, the plurality of non-classical computers includes at least one quantum computer. In some embodiments, the at least one quantum computer includes one or more members selected from the group consisting of a quantum hardware device and a classical simulator of a quantum circuit. In some embodiments, one of the quantum mechanical energies includes nuclear-nuclear repulsion energy.

[0013] In some embodiments, the method further includes providing an input to a hybrid computing unit, the input including a set of atomic coordinates of the chemical system. In some embodiments, the method further includes performing (a)-(c) for two or more conformations in the group of conformations of the chemical system. In some embodiments, the method further includes classifying a combined quantum mechanical energy or electronic structure of at least a subset of the plurality of molecular fragments.

[0014] In some embodiments, (a) includes applying one or more members selected from the group consisting of a fragment molecular orbital (FMO) method, a divide-and-conquer (DC) method, a density matrix embedding theory (DMET) method, a density matrix renormalization group (DMRG) method, a tensor network, and an incremental method.

[0015] In some embodiments, (d) includes determining a fermion Hamiltonian (molecular Hamiltonian or electronic Hamiltonian) for molecular fragments of at least a subset of the plurality of molecular fragments, converting the fermion Hamiltonian to an equivalent qubit Hamiltonian, converting the qubit Hamiltonian to a quantum circuit, and using the quantum circuit to determine the quantum mechanical energy or electronic structure of the molecular fragments. In some embodiments, the method further includes determining the quantum mechanical energy or electronic structure using the molecular Hamiltonian. In some embodiments, the method further includes determining the quantum mechanical energy or electronic structure using the electronic Hamiltonian. In some embodiments, converting the fermion Hamiltonian to an equivalent qubit Hamiltonian includes converting fermion operators of the Hamiltonian to qubit operators.

[0016] In some embodiments, the method further includes performing an ab initio molecular dynamics (AIMD) simulation of the chemical system. In some embodiments, the AIMD simulation includes, prior to (a), acquiring an indication of the chemical system, the indication including coordinates of each particle of a plurality of particles in the chemical system and a velocity of each particle in the chemical system, and, subsequent to (c), (i) determining forces on each particle in the chemical system from the combined energy or electronic structure, (ii) updating the coordinates of each particle in the chemical system and the velocity of each particle in the chemical system, and (iii) electronically outputting a report indicating the coordinates or velocities. In some embodiments, (i) includes applying Jordan's quantum algorithm for numerical gradient estimation to the quantum mechanical energy or electronic structure. In some embodiments, (ii) includes applying one or more members selected from the group consisting of a Vere procedure, a velocity Vere procedure, a symplectic integration, a Runge-Kutta integration, and a Biemann integration.

[0017] In another aspect, a system for performing a quantum mechanical energy or electronic structure calculation on a chemical system may comprise a memory including instructions for performing a quantum mechanical energy or electronic structure calculation on the chemical system, and a hybrid computing unit operably coupled to the memory, the hybrid computing unit comprising a distributed computing system comprising at least one classical computer and a plurality of non-classical computers, the hybrid computing unit configured to execute the instructions to at least (a) decompose at least one conformation in a group of conformations of the chemical system into a plurality of molecular fragments, (b) determine a quantum mechanical energy or electronic structure of at least a subset of the plurality of molecular fragments, (c) combine the quantum mechanical energies or electronic structures determined in (b), and (d) electronically output a report indicative of the combined quantum mechanical energy or electronic structure in (c).

[0018] In another aspect, a non-transitory computer readable medium may comprise machine executable code that, when executed by a hybrid computing unit comprising at least one classical computer and a distributed computing system comprising a plurality of non-classical computers, implements a method for performing quantum mechanical energy or electronic structure calculations on a chemical system, the method including: (a) decomposing at least one conformation in a group of conformations of the chemical system into a plurality of molecular fragments; (b) determining the quantum mechanical energy or electronic structure of at least a subset of the plurality of molecular fragments; (c) combining the quantum mechanical energies or electronic structures determined in (b); and (d) electronically outputting a report indicative of the quantum mechanical energies or electronic structures combined in (c).

[0019] In another aspect, the disclosure provides a method for performing a quantum mechanical energy or electronic structure calculation on a chemical system, the method being performed by a hybrid computing unit comprising a distributed computing system comprising a plurality of classical computers and at least one non-classical computer, the method including: (a) decomposing at least one conformation in a group of conformations of the chemical system into a plurality of molecular fragments; (b) using the hybrid computing unit, determining a quantum mechanical energy or electronic structure of at least a subset of the plurality of molecular fragments; (c) combining the quantum mechanical energies or electronic structures determined in (b); and (d) electronically outputting a report indicative of the quantum mechanical energies or electronic structures combined in (c).

[0020] In some embodiments, the at least one non-classical computer comprises at least one quantum computer. In some embodiments, the at least one quantum computer comprises one or more members selected from the group consisting of a quantum hardware device and a classical simulator of a quantum circuit. In some embodiments, one of the quantum mechanical energies comprises nuclear-nuclear repulsion energy.

[0021] In some embodiments, the method further includes providing an input to a hybrid computing unit, the input including a set of atomic coordinates of the chemical system. In some embodiments, the method further includes performing (a)-(c) for two or more conformations in the group of conformations of the chemical system. In some embodiments, the method further includes classifying a combined quantum mechanical energy or electronic structure of at least a subset of the plurality of molecular fragments.

[0022] In some embodiments, (a) includes applying one or more members selected from the group consisting of a fragment molecular orbital (FMO) method, a divide-and-conquer (DC) method, a density matrix embedding theory (DMET) method, a density matrix renormalization group (DMRG) method, a tensor network, and an incremental method.

[0023] In some embodiments, (b) includes determining a fermion Hamiltonian (molecular Hamiltonian or electronic Hamiltonian) for molecular fragments of at least a subset of the plurality of molecular fragments, converting the fermion Hamiltonian to an equivalent qubit Hamiltonian, converting the qubit Hamiltonian to a quantum circuit, and using the quantum circuit to determine the quantum mechanical energy or electronic structure of the molecular fragments. In some embodiments, the method further includes determining the quantum mechanical energy or electronic structure using the molecular Hamiltonian. In some embodiments, the method further includes determining the quantum mechanical energy or electronic structure using the electronic Hamiltonian. In some embodiments, converting the fermion Hamiltonian to an equivalent qubit Hamiltonian includes converting fermion operators of the Hamiltonian to qubit operators.

[0024] In some embodiments, the method further includes performing an ab initio molecular dynamics (AIMD) simulation of the chemical system. In some embodiments, the AIMD simulation includes, prior to (a), acquiring an indication of the chemical system, the indication including coordinates of each particle of a plurality of particles in the chemical system and a velocity of each particle in the chemical system, and, subsequent to (c), (i) determining forces on each particle in the chemical system from the combined energy or electronic structure, (ii) updating the coordinates of each particle in the chemical system and the velocity of each particle in the chemical system, and (iii) electronically outputting a report indicating the coordinates or velocities. In some embodiments, (i) includes applying Jordan's quantum algorithm for numerical gradient estimation to the quantum mechanical energy or electronic structure. In some embodiments, (ii) includes applying one or more members selected from the group consisting of a Vere procedure, a velocity Vere procedure, a symplectic integration, a Runge-Kutta integration, and a Biemann integration.

[0025] In some embodiments, the method further includes dispatching one or more of the plurality of fragments to one or more remote endpoints and receiving quantum mechanical energy or electronic structure from the one or more remote endpoints. In some embodiments, at least one of the one or more remote endpoints includes a non-classical computer. In some embodiments, the one or more remote endpoints includes a portion of a cloud computing system. In some embodiments, the method further includes, prior to (a), receiving at least one conformation from the client-side library and dispatching the at least one conformation to the first remote endpoint. In some embodiments, at least one of (a) and (c) occurs at the first remote endpoint. In some embodiments, the method further includes dispatching one or more of the plurality of fragments to one or more second remote endpoints and receiving quantum mechanical energy or electronic structure from the one or more second remote endpoints. In some embodiments, the method further includes sending the report to the client-side library. In some embodiments, at least one of the second remote endpoints includes a non-classical computer. In some embodiments, the one or more remote endpoints includes a portion of a cloud computing system.

[0026] In another aspect, a system for performing a quantum mechanical energy or electronic structure calculation on a chemical system may comprise a computer memory including instructions for performing a quantum mechanical energy or electronic structure calculation on the chemical system, and a hybrid computing unit operably coupled to the memory, the hybrid computing unit comprising a distributed computing system comprising a plurality of classical computers and at least one non-classical computer, the hybrid computing unit configured to execute the instructions to at least (a) decompose at least one conformation in a group of conformations of the chemical system into a plurality of molecular fragments, (b) determine a quantum mechanical energy or electronic structure of at least a subset of the plurality of molecular fragments, (c) combine the quantum mechanical energies or electronic structures determined in (b), and (d) electronically output a report indicative of the combined quantum mechanical energy or electronic structure in (c).

[0027] In another aspect, a non-transitory computer readable medium may comprise machine executable code that, when executed by a hybrid computing unit comprising a distributed computing system comprising a plurality of classical computers and at least one non-classical computer, implements a method for performing a quantum mechanical energy or electronic structure calculation on a chemical system, the method including: (a) decomposing at least one conformation in a group of conformations of the chemical system into a plurality of molecular fragments; (b) determining a quantum mechanical energy or electronic structure of at least a subset of the plurality of molecular fragments; (c) combining the quantum mechanical energies or electronic structures determined in (b); and (d) electronically outputting a report indicative of the quantum mechanical energies or electronic structures combined in (c).

[0028] In another aspect, a method for performing quantum mechanical energy or electronic structure calculations on a chemical system is provided. The method may be performed by a hybrid computing unit comprising at least one classical computer and a distributed computing system comprising a plurality of non-classical computers. The method may include (a) decomposing at least one conformation in a group of conformations of the chemical system into a plurality of molecular fragments, (b) using the hybrid computing unit to determine the quantum mechanical energy or electronic structure of at least a subset of the plurality of molecular fragments, (c) combining the quantum mechanical energies or electronic structures determined in (b), and (d) electronically outputting a report indicative of the quantum mechanical energies or electronic structures combined in (c).

[0029] In some embodiments, the plurality of non-classical computers includes at least one quantum computer. In some embodiments, the at least one quantum computer includes one or more members selected from the group consisting of a quantum hardware device and a classical simulator of a quantum circuit. In some embodiments, the plurality of non-classical computers includes different types of non-classical computers. In some embodiments, one of the quantum mechanical energies includes a nuclear-nuclear repulsion energy. In some embodiments, the method further includes providing an input to the hybrid computing unit, the input including a set of atomic coordinates of the chemical system. In some embodiments, the method further includes performing (a)-(c) for two or more conformations within the group of conformations of the chemical system. In some embodiments, the method further includes classifying the combined quantum mechanical energy or electronic structure of the at least the subset of the plurality of molecular fragments.

[0030] In some embodiments, (a) includes applying one or more members selected from the group consisting of a fragment molecular orbital (FMO) method, a divide and conquer (DC) method, a density matrix embedding theory (DMET) method, a density matrix renormalization group (DMRG) method, a tensor network, and an incremental method. In some embodiments, (b) includes determining a fermion Hamiltonian for the molecular fragments of at least the subset of the plurality of molecular fragments, converting the fermion Hamiltonian to an equivalent qubit Hamiltonian, converting the qubit Hamiltonian to a quantum circuit, and determining a quantum mechanical energy or electronic structure of the molecular fragments using the quantum circuit. In some embodiments, the method further includes determining the quantum mechanical energy or electronic structure using a molecular Hamiltonian. In some embodiments, the method further includes determining the quantum mechanical energy or electronic structure using an electronic Hamiltonian. In some embodiments, converting the fermion Hamiltonian to an equivalent qubit Hamiltonian includes converting fermion operators of the Hamiltonian to qubit operators.

[0031] In some embodiments, the method further includes performing an ab initio molecular dynamics (AIMD) simulation of the chemical system. In some embodiments, the AIMD simulation includes, prior to (a), acquiring an indication of a chemical system, the indication including coordinates of each particle of a plurality of particles in the chemical system and a velocity of each particle in the chemical system, and, subsequent to (c), (i) determining forces on each particle in the chemical system from the combined energy or electronic structure, (ii) updating the coordinates of each particle in the chemical system and the velocity of each particle in the chemical system, and (iii) electronically outputting a report indicating the coordinates or the velocities. In some embodiments, (i) includes applying Jordan's quantum algorithm for numerical gradient estimation to the quantum mechanical energy or electronic structure. In some embodiments, (ii) includes applying one or more members selected from the group consisting of a Vere procedure, a velocity Vere procedure, a symplectic integration, a Runge-Kutta integration, and a Biemann integration.

[0032] In some embodiments, the method further includes dispatching one or more of the plurality of fragments to one or more remote endpoints and receiving the quantum mechanical energy or electronic structure from the one or more remote endpoints. In some embodiments, at least one of the one or more remote endpoints comprises a non-classical computer. In some embodiments, the one or more remote endpoints comprise part of a cloud computing system. In some embodiments, the method further includes, prior to (a), receiving the at least one conformation from a client-side library and dispatching the at least one conformation to a first remote endpoint. In some embodiments, at least one of (a) and (c) occurs at the first remote endpoint.

[0033] In some embodiments, the method further includes dispatching one or more of the plurality of fragments to one or more second remote endpoints and receiving the quantum mechanical energy or electronic structure from the one or more second remote endpoints. In some embodiments, the method further includes sending the report to the client-side library. In some embodiments, at least one of the second remote endpoints includes a non-classical computer. In some embodiments, the one or more remote endpoints include portions of a cloud computing system. In some embodiments, the decomposing in (a) is performed using the at least one classical computer. In some embodiments, the determining in (b) is performed using at least one non-classical computer of the plurality of non-classical computers. In some embodiments, the combining in (c) is performed using the at least one classical computer.

[0034] In another aspect, a system for performing quantum mechanical energy or electronic structure calculations on a chemical system is provided, the system may include a hybrid computing unit operably coupled to the memory, the hybrid computing unit including a distributed computing system including at least one classical computer and a plurality of non-classical computers, the hybrid computing unit configured to at least (a) decompose at least one conformation in a group of conformations of the chemical system into a plurality of molecular fragments, (b) determine quantum mechanical energies or electronic structures of at least a subset of the plurality of molecular fragments, (c) combine the quantum mechanical energies or electronic structures determined in (b), and (d) electronically output a report indicative of the quantum mechanical energies or electronic structures combined in (c).

[0035] In some embodiments, the system further comprises a computer memory containing instructions for performing said quantum mechanical energy or electronic structure calculations on said chemical system, and said hybrid computing unit is configured to execute said instructions to perform at least (a)-(d).

[0036] In another aspect, a non-transitory computer readable medium is provided that includes machine executable code that, when executed by a hybrid computing unit comprising at least one classical computer and a distributed computing system comprising a plurality of non-classical computers, implements a method for performing a quantum mechanical energy or electronic structure calculation on a chemical system, the method may include (a) decomposing at least one conformation in a group of conformations of the chemical system into a plurality of molecular fragments, (b) determining a quantum mechanical energy or electronic structure of at least a subset of the plurality of molecular fragments, (c) combining the quantum mechanical energies or electronic structures determined in (b), and (d) electronically outputting a report indicative of the quantum mechanical energies or electronic structures combined in (c).

[0037] In another embodiment, a method for performing quantum mechanical energy or electronic structure calculations on a chemical system is provided. The method may be performed by a hybrid computing unit comprising a distributed computing system comprising a plurality of classical computers and at least one non-classical computer. The method may include (a) decomposing at least one conformation in a group of conformations of the chemical system into a plurality of molecular fragments, (b) using the hybrid computing unit to determine the quantum mechanical energy or electronic structure of at least a subset of the plurality of molecular fragments, (c) combining the quantum mechanical energies or electronic structures determined in (b), and (d) electronically outputting a report indicative of the quantum mechanical energies or electronic structures combined in (c).

[0038] In some embodiments, the at least one non-classical computer comprises at least one quantum computer. In some embodiments, the at least one quantum computer comprises one or more members selected from the group consisting of a quantum hardware device and a classical simulator of a quantum circuit. In some embodiments, the at least one non-classical computer comprises a plurality of different types of non-classical computers. In some embodiments, one of the quantum mechanical energies comprises a nuclear-nuclear repulsion energy. In some embodiments, the method further comprises providing an input to the hybrid computing unit, the input comprising a set of atomic coordinates of the chemical system. In some embodiments, the method further comprises performing (a)-(c) for two or more conformations within the group of conformations of the chemical system. In some embodiments, the method further comprises classifying the combined quantum mechanical energies or electronic structures of the at least the subset of the plurality of molecular fragments.

[0039] In some embodiments, (a) includes applying one or more members selected from the group consisting of a fragment molecular orbital (FMO) method, a divide-and-conquer (DC) method, a density matrix embedding theory (DMET) method, a density matrix renormalization group (DMRG) method, a tensor network, and an incremental method.

[0040] In some embodiments, (b) includes (i) determining a fermion Hamiltonian for molecular fragments of at least the subset of the plurality of molecular fragments, (ii) converting the fermion Hamiltonian to an equivalent qubit Hamiltonian, (iii) converting the qubit Hamiltonian to a quantum circuit, and (iv) determining a quantum mechanical energy or electronic structure of the molecular fragment using the quantum circuit. In some embodiments, the method further includes determining the quantum mechanical energy or electronic structure using a molecular Hamiltonian. In some embodiments, the method further includes determining the quantum mechanical energy or electronic structure using an electronic Hamiltonian. In some embodiments, converting the fermion Hamiltonian to an equivalent qubit Hamiltonian includes converting a fermion operator of a Hamiltonian to a qubit operator.

[0041] In some embodiments, the method further includes performing an ab initio molecular dynamics (AIMD) simulation of the chemical system. In some embodiments, the AIMD simulation includes, prior to (a), acquiring an indication of a chemical system, the indication including coordinates of each particle of a plurality of particles in the chemical system and a velocity of each particle in the chemical system, and, subsequent to (c), (i) determining forces on each particle in the chemical system from the combined energy or electronic structure, (ii) updating the coordinates of each particle in the chemical system and the velocity of each particle in the chemical system, and (iii) electronically outputting a report indicating the coordinates or the velocities. In some embodiments, (i) includes applying Jordan's quantum algorithm for numerical gradient estimation to the quantum mechanical energy or electronic structure. In some embodiments, (ii) includes applying one or more members selected from the group consisting of a Vere procedure, a velocity Vere procedure, a symplectic integration, a Runge-Kutta integration, and a Biemann integration.

[0042] In some embodiments, the method further includes dispatching one or more of the fragments to one or more remote endpoints and receiving the quantum mechanical energy or electronic structure from the one or more remote endpoints. In some embodiments, at least one of the one or more remote endpoints includes a non-classical computer. In some embodiments, the one or more remote endpoints include portions of a cloud computing system.

[0043] In some embodiments, the method further includes, prior to (a), receiving the at least one conformation from a client-side library and dispatching the at least one conformation to a first remote endpoint. In some embodiments, at least one of (a) and (c) occurs at the first remote endpoint. In some embodiments, the method further includes dispatching one or more of the plurality of fragments to one or more second remote endpoints and receiving the quantum mechanical energy or electronic structure from the one or more second remote endpoints. In some embodiments, the method further includes transmitting the report to the client-side library. In some embodiments, at least one of the second remote endpoints comprises a non-classical computer. In some embodiments, the one or more remote endpoints comprise portions of a cloud computing system.

[0044] In some embodiments, the decomposing in (a) is performed using at least one classical computer of the plurality of classical computers. In some embodiments, the determining in (b) is performed using the at least one non-classical computer. In some embodiments, the combining in (c) is performed using at least one classical computer of the plurality of classical computers.

[0045] In another aspect, a system for performing quantum mechanical energy or electronic structure calculations on a chemical system is provided, which may include a hybrid computing unit operably coupled to the memory, the hybrid computing unit comprising a distributed computing system including a plurality of classical computers and at least one non-classical computer, the hybrid computing unit configured to at least (a) resolve at least one conformation in a group of conformations of the chemical system into a plurality of molecular fragments, (b) determine quantum mechanical energies or electronic structures of at least a subset of the plurality of molecular fragments, (c) combine the quantum mechanical energies or electronic structures determined in (b), and (d) electronically output a report indicative of the quantum mechanical energies or electronic structures combined in (c).

[0046] In some embodiments, the system further comprises a computer memory containing instructions for performing said quantum mechanical energy or electronic structure calculations on said chemical system, and said hybrid computing unit is configured to execute said instructions to perform at least (a)-(d).

[0047] In another aspect, a non-transitory computer readable medium is provided that includes machine executable code that, when executed by a hybrid computing unit comprising a distributed computing system comprising a plurality of classical computers and at least one non-classical computer, performs a method for performing a quantum mechanical energy or electronic structure calculation on a chemical system, the method may include (a) decomposing at least one conformation in a group of conformations of the chemical system into a plurality of molecular fragments, (b) determining a quantum mechanical energy or electronic structure of at least a subset of the plurality of molecular fragments, (c) combining the quantum mechanical energies or electronic structures determined in (b), and (d) electronically outputting a report indicative of the quantum mechanical energies or electronic structures combined in (c).

[0048] In another aspect, a method for performing quantum mechanical energy or electronic structure calculations on a chemical system is provided. The method may be performed by a hybrid computing unit comprising at least one classical computer and a distributed computing system comprising a plurality of non-classical computers. The method may include (a) decomposing at least one conformation in a group of conformations of the chemical system into a plurality of molecular fragments, (b) dispatching a subset of the plurality of molecular fragments to a plurality of solvers, (c) using the plurality of solvers to determine quantum mechanical energies or electronic structures of a plurality of molecular fragments of the subset of the plurality of molecular fragments, and (d) electronically outputting a report indicative of the quantum mechanical energies or electronic structures determined in (c).

[0049] In some embodiments, the plurality of non-classical computers includes at least one quantum computer. In some embodiments, the at least one quantum computer includes one or more members selected from the group consisting of a quantum hardware device and a classical simulator of a quantum circuit. In some embodiments, the plurality of non-classical computers includes different types of non-classical computers. In some embodiments, one of the quantum mechanical energies includes a nuclear-nuclear repulsion energy. In some embodiments, the method further includes providing an input to the hybrid computing unit, the input including a set of atomic coordinates of the chemical system. In some embodiments, the method further includes performing (a)-(c) for two or more conformations within the group of conformations of the chemical system. In some embodiments, the method further includes the combined quantum mechanical energies or electronic structures of the at least the subset of the plurality of molecular fragments. In some embodiments, (a) includes applying one or more members selected from the group consisting of a fragment molecular orbital (FMO) method, a divide-and-conquer (DC) method, a density matrix embedding theory (DMET) method, a density matrix renormalization group (DMRG) method, a tensor network, and an incremental method.

[0050] In some embodiments, (c) includes determining a fermion Hamiltonian for a molecular fragment of the plurality of molecular fragments, converting the fermion Hamiltonian to an equivalent qubit Hamiltonian, converting the qubit Hamiltonian to a quantum circuit, and determining a quantum mechanical energy or electronic structure of the molecular fragment using the quantum circuit. In some embodiments, the method further includes determining the quantum mechanical energy or electronic structure using a molecular Hamiltonian. In some embodiments, the method further includes determining the quantum mechanical energy or electronic structure using an electronic Hamiltonian. In some embodiments, converting the fermion Hamiltonian to an equivalent qubit Hamiltonian includes converting a fermion operator of a Hamiltonian to a qubit operator.

[0051] In some embodiments, the method further includes performing an ab initio molecular dynamics (AIMD) simulation of the chemical system. In some embodiments, the AIMD simulation includes, prior to (a), acquiring an indication of a chemical system, the indication including coordinates of each particle of a plurality of particles in the chemical system and a velocity of each particle in the chemical system, and, subsequent to (c), (i) determining forces on each particle in the chemical system from the combined energy or electronic structure, (ii) updating the coordinates of each particle in the chemical system and the velocity of each particle in the chemical system, and (iii) electronically outputting a report indicating the coordinates or the velocities. In some embodiments, (i) includes applying Jordan's quantum algorithm for numerical gradient estimation to the quantum mechanical energy or electronic structure. In some embodiments, (ii) includes applying one or more members selected from the group consisting of a Vere procedure, a velocity Vere procedure, a symplectic integration, a Runge-Kutta integration, and a Biemann integration.

[0052] In some embodiments, the plurality of solvers include one or more remote endpoints. In some embodiments, the method further includes receiving the quantum mechanical energy or electronic structure from the one or more remote endpoints. In some embodiments, at least one of the one or more remote endpoints includes a non-classical computer. In some embodiments, the one or more remote endpoints include portions of a cloud computing system.

[0053] In some embodiments, the method further includes, prior to (a), receiving the at least one conformation from a client-side library and dispatching the at least one conformation to a first remote endpoint. In some embodiments, at least one of (a) and (c) occurs at the first remote endpoint. In some embodiments, the method further includes dispatching one or more of the plurality of fragments to one or more second remote endpoints and receiving the quantum mechanical energy or electronic structure from the one or more second remote endpoints. In some embodiments, the method further includes transmitting the report to the client-side library. In some embodiments, at least one of the second remote endpoints comprises a non-classical computer. In some embodiments, the one or more remote endpoints comprise portions of a cloud computing system.

[0054] In some embodiments, the decomposing in (a) is performed using the at least one classical computer. In some embodiments, the dispatching in (b) is performed using the at least one classical computer. In some embodiments, the dispatching in (b) is performed using a classical computer separate from the at least one classical computer. In some embodiments, the determining in (c) is performed using at least one non-classical computer of the plurality of non-classical computers. In some embodiments, the outputting in (d) is performed using the at least one classical computer.

[0055] In another aspect, a system for performing a quantum mechanical energy or electronic structure calculation on a chemical system is provided. The system may comprise a hybrid computing unit operably coupled to the memory, the hybrid computing unit comprising a distributed computing system comprising a plurality of classical computers and at least one non-classical computer, the hybrid computing unit configured to at least (a) decompose at least one conformation in a group of conformations of the chemical system into a plurality of molecular fragments, (b) dispatch a subset of the plurality of molecular fragments to a plurality of solvers, (c) determine quantum mechanical energies or electronic structures of a plurality of molecular fragments of the subset of the plurality of molecular fragments using the plurality of solvers, and (d) electronically output a report indicative of the quantum mechanical energy or electronic structure combined in (c). In some embodiments, the system further comprises a computer memory including instructions for performing the quantum mechanical energy or electronic structure calculation on the chemical system, the hybrid computing unit configured to execute the instructions to perform at least (a)-(d).

[0056] In another aspect, a non-transitory computer readable medium is provided that includes machine executable code that, when executed by a hybrid computing unit comprising a distributed computing system comprising a plurality of classical computers and at least one non-classical computer, performs a method for performing a quantum mechanical energy or electronic structure calculation on a chemical system, the method may include (a) decomposing at least one conformation in a group of conformations of the chemical system into a plurality of molecular fragments, (b) determining a quantum mechanical energy or electronic structure of at least a subset of the plurality of molecular fragments, (c) combining the quantum mechanical energies or electronic structures determined in (b), and (d) electronically outputting a report indicative of the quantum mechanical energies or electronic structures combined in (c).

[0057] Another aspect of the present disclosure provides a non-transitory computer-readable medium that includes machine-executable code that, when executed by one or more computer processors, performs any of the methods described above or elsewhere herein.

[0058] Another aspect of the present disclosure provides a system comprising one or more computer processors and a computer memory coupled thereto, the computer memory including machine executable code that, when executed by the one or more computer processors, performs any of the methods described above or elsewhere herein.

[0059] Additional aspects and advantages of the present disclosure will become readily apparent to those skilled in the art from the following detailed description, in which only illustrative embodiments of the present disclosure have been shown and described. As will be recognized, the present disclosure is capable of other and different embodiments, and its several details are capable of modification in various obvious respects, all without departing from the present disclosure. Accordingly, the drawings and description are to be regarded as illustrative in nature, and not as restrictive.

[0060] Incorporation by Reference All publications, patents, and patent applications mentioned herein are incorporated by reference to the same extent as if each individual publication, patent, or patent application was specifically and individually indicated to be incorporated by reference. To the extent that the publications and patents or patent applications incorporated by reference conflict with the disclosure contained herein, the present specification is intended to supersede and / or precede any such conflicting material.

[0061] The novel features of the invention are set forth with particularity in the appended claims. A better understanding of the features and advantages of the present invention will be obtained by reference to the following detailed description that sets forth illustrative embodiments, in which the principles of the invention are utilized, and the accompanying drawings (herein referred to as "Figures"), in which: [Brief description of the drawings]

[0062] [Figure 1] 1 is a flow chart for an example method for providing an index to a classified list of conformers of a molecule using problem decomposition techniques on quantum computing hardware, according to some embodiments disclosed herein. [Diagram 2] 1 is a flowchart for an example method for providing indications of quantum mechanical energy and / or electronic structure of a subsystem, as defined by problem decomposition techniques, on quantum computing hardware according to some embodiments disclosed herein. [Diagram 3] 1 is a flowchart for an example method for providing a measure of a Hamiltonian expectation value on quantum computing hardware according to some embodiments disclosed herein. [Figure 4] FIG. 1 is an exemplary diagram of n-heptane, with dotted lines indicating the Bond Detached Atom (BTA) in the Fragment Molecular Orbital (FMO) fragmentation. [Diagram 5]FIG. 1 is an exemplary diagram of n-heptane showing a comparison between the results obtained by canonical CCSD and divide-and-conquer CCSD (DC-CCSD), and between the results obtained by canonical CCSD and fragment molecular orbital CCSD (FMO-CCSD). [Figure 6] An exemplary diagram of n-heptane [left] showing the minimum sphere to accommodate a conformer (dotted circle) and the distance (solid line) between the terminal carbon atoms involved in the dihedral angle (1-4 distance); a plot showing the total quantum mechanical energy (left arrow) versus the diameter of the minimum sphere (right arrow) for each of the conformers classified based on the total quantum mechanical energy [center]; and a plot showing the total quantum mechanical energy (left arrow) versus the minimum 1-4 distance (right arrow) for each conformer [right]. [Figure 7] FIG. 2 is an exemplary diagram of 3-methylheptane, with dotted lines indicating the Bond Detached Atom (BTA) in the Fragment Molecular Orbital (FMO) fragmentation. [Figure 8] FIG. 1 illustrates the quantum mechanical energy distributions of n-heptane (blue) and 3-methylheptane (red). [Figure 9] FIG. 1 is an exemplary diagram of 3-methylheptane showing a comparison between the results obtained by canonical CCSD and divide-and-conquer CCSD (DC-CCSD), and between the results obtained by canonical CCSD and fragment molecular orbital CCSD (FMO-CCSD). [Figure 10] FIG. 1 illustrates a computer control system that is programmed or otherwise configured to carry out the methods provided herein. [Figure 11] 1 is a flow chart for an example of an incremental method for performing problem decomposition. [Figure 12] FIG. 1 illustrates molecular orbitals, atoms, molecular fragments, and molecules used as the basis for the incremental method. [Figure 13]1 is a flow chart for an example method for performing first-principles molecular dynamics (AIMD) for molecules using problem decomposition techniques on quantum computing hardware, according to some embodiments disclosed herein. [Figure 14] 1 is a flowchart for an example of a method for calculating forces on each particle of a system in an ab initio molecular dynamics (AIMD) simulation, according to some embodiments disclosed herein. [Figure 15] FIG. 1 illustrates an example of a system or combination of systems that may be used to solve a problem, such as a quantum chemistry problem or simulation. [Figure 16] 1 is a flow chart for an example of a method for performing quantum mechanical energy or electronic structure calculations on a chemical system using a distributed computing system. [Figure 17] A diagram illustrating an example of the architecture of a distributed computing system comprising non-classical (e.g., quantum computers) and multiple classical computers. [Figure 18] FIG. 1 illustrates a distributed computing system including sequential problem decomposition, according to some embodiments disclosed herein. [Figure 19] FIG. 1 illustrates a distributed computing system including problem dispatch according to some embodiments disclosed herein. [Figure 20] FIG. 1 illustrates an example architecture of a distributed computing system including problem dispatch in a client-side library, according to some embodiments disclosed herein. [Figure 21] FIG. 1 illustrates an example architecture of a distributed computing system including problem dispatch at a remote endpoint, according to certain embodiments disclosed herein. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0063] While various embodiments of the present invention have been shown and described herein, it will be apparent to those skilled in the art that such embodiments are provided by way of example only. Numerous variations, changes, and substitutions may occur to those skilled in the art without departing from the invention. It is understood that various alternatives to the embodiments of the invention described herein may be used.

[0064] Unless otherwise specified, all technical terms used herein have the same meaning as commonly understood by those skilled in the art to which this invention belongs. As used in this specification and the appended claims, the singular forms "a", "an" and "the" include plural references unless the context clearly indicates otherwise. Any reference to "or" in this specification is intended to include "and / or" unless otherwise specified.

[0065] Whenever the terms "at least," "greater than," or "greater than or equal to" precede the first number in a series of two or more numbers, the term "at least," "greater than," or "greater than or equal to" applies to each and every number in the series. For example, 1, 2, or 3 or more equals 1 or more, 2 or more, or 3 or more.

[0066] Whenever the term "no more than," "less than," or "less than or equal to" precedes the first number in a series of two or more numbers, the term "less than," "less than," or "less than or equal to" applies to each and every number in the series. For example, 3, 2, or 1 or less is equal to 3 or less, 2 or less, or 1 or less.

[0067] In the following detailed description, reference is made to the accompanying drawings, which form a part hereof. In the drawings, like numerals typically identify like components unless otherwise noted. The illustrative embodiments described in the detailed description, figures, and claims are not intended to be limiting. Other embodiments may be utilized, and other changes may be made, without departing from the scope of the subject matter presented herein. It will be readily understood that the aspects of the present disclosure as generally described herein and illustrated in the figures can be arranged, substituted, combined, separated, and designed in a wide variety of different configurations, all of which are expressly contemplated herein.

[0068] The present disclosure provides a method for applying problem decomposition (PD) techniques in quantum chemistry to identifying and predicting the quantum mechanical energy and / or electronic structure of a chemical system, or to identifying a set of the most energetically stable conformers of a molecule. The systems and methods provided herein for implementing PD techniques on a QC platform may allow quantum mechanical energy and / or electronic structure calculations to be performed with a high level of accuracy for each fragment. Furthermore, the small size of each fragment may allow highly accurate calculations to be performed on a QC device where the scale of calculations is rather limited, thus efficiently and accurately obtaining the energy and / or electronic structure of complex, industrially relevant molecules. The methods and systems described herein may be applied not only to a single chemical system, but also to molecular assemblies with different association structures. For example, the methods and systems disclosed herein may be applied to identifying the most stable binding orientation of a drug candidate to a target protein from a group of potential binding orientations.

[0069] In some cases, a classical computer may be configured to perform one or more classical algorithms. A classical algorithm (or classical computational task) may include an algorithm (or computational task) that can be performed by one or more classical computers without the use of a quantum computer, a quantum-ready computing service, or a quantum-enabled computing service. A classical algorithm may include a non-quantum algorithm. A classical computer may include a computer that does not include a quantum computer, a quantum-ready computing service, or a quantum-enabled computer. A classical computer may process or store data represented by digital bits (e.g., zeros ("0") and ones ("1")) rather than quantum bits (qubits). Examples of classical computers include, but are not limited to, server computers, desktop computers, laptop computers, notebook computers, sub-notebook computers, netbook computers, netpad computers, set-top computers, media streaming devices, handheld computers, Internet appliances, mobile smartphones, tablet computers, personal digital assistants, home game consoles, and vehicles.

[0070] The hybrid computing unit may comprise a classical computer and a quantum computer. The quantum computer may be configured to implement one or more quantum algorithms for solving a computational problem (e.g., at least a portion of a quantum chemistry simulation). The one or more quantum algorithms may be executed using a quantum computer, a quantum ready computing service, or a quantum enabled computing service. For example, the one or more quantum algorithms may be executed using a system or method described in U.S. Patent Publication No. 2018 / 0107526, entitled "METHODS AND SYSTEMS FOR QUANTUM READY AND QUANTUM ENABLED COMPUTATIONS," which is incorporated herein by reference in its entirety. The classical computer may comprise at least one classical processor and computer memory, and may be configured to implement one or more quantum algorithms for solving a computational problem (e.g., at least a portion of a quantum chemistry simulation). The digital computer may comprise at least one computer processor and computer memory, and the digital computer may include a computer program having instructions executable by the at least one computer processor to render an application. The application may facilitate the use of quantum and / or classical computers by users.

[0071] Some implementations may use quantum computers along with classical computers that operate on bits, such as personal desktops, laptops, supercomputers, distributed computing clusters, cloud-based computing resources, smartphones, or tablets.

[0072] The system may include an interface for a user. In some cases, the interface may include an Application Programming Interface (API). The interface may provide a programmatic model that ignores (e.g., by hiding from the user) the internal details (e.g., architecture and operation) of the quantum computer. In some cases, the interface may minimize the need to update application programs in response to quantum hardware changes. In some cases, the interface may remain unchanged when the quantum computer has changes in its internal structure.

[0073] The present disclosure provides systems and methods that may include non-classical (e.g., quantum) computing or the use of non-classical (e.g., quantum) computing. Quantum computers may be capable of solving certain classes of computational tasks more efficiently than classical computers. However, quantum computing resources may be scarce and expensive and may involve a certain level of expertise to be used efficiently or effectively (e.g., cost-effectively or cost-effectively). Several parameters may be tuned for a quantum computer to deliver its potential computational power.

[0074] A quantum computer (or other type of non-classical computer) may be capable of functioning in parallel with a classical computer as a co-processor. A hybrid architecture (e.g., a computing system) comprising a classical computer and a quantum computer may be highly efficient in addressing complex computational tasks such as quantum chemistry simulations. The systems and methods disclosed herein may be capable of efficiently and accurately decomposing or partitioning a quantum chemistry problem and delegating appropriate components of the quantum chemistry simulation to a quantum computer or a classical computer.

[0075] Although this disclosure refers to quantum computers, the methods and systems of this disclosure may be used for use with other types of computers, which may be non-classical computers. Such non-classical computers may include quantum computers, hybrid quantum computers, quantum-like computers, or other computers that are not classical computers. Examples of non-classical computers may include, but are not limited to, Hitachi Ising solvers, coherent Ising machines based on optical parameters, and other solvers that utilize different physical phenomena to gain additional efficiency in solving certain classes of problems.

[0076] In some cases, the quantum computer may include one or more of adiabatic quantum computers, quantum gate arrays, one-way quantum computers, topological quantum computers, quantum Turing machines, semiconductor-based quantum computers, trapped-ion quantum computers, trapped-atom quantum computers, optical lattices, quantum dot computers, spin-based quantum computers, spatial-based quantum computers, Loss-DiVincenzo quantum computers, Nuclear Magnetic Resonance (NMR)-based quantum computers, solution NMR quantum computers, solid-state NMR quantum computers, solid-state NMR Kane quantum computers, helium surface electron quantum computers, resonator quantum electrodynamics-based quantum computers, molecular magnet quantum computers, fullerene-based quantum computers, linear optical quantum computers, diamond-based quantum computers, Nitrogen Vacancy (NV) diamond-based quantum computers, Bose-Einstein condensate-based quantum computers, transistor-based quantum computers, and rare earth metal ion doped inorganic crystal-based quantum computers. The quantum computer may include one or more of a quantum annealer, an Ising solver, an Optical Parametric Oscillator (OPO), and a gate model of quantum computing.

[0077] In some cases, the non-classical computers of the present disclosure may include noisy intermediate-scale quantum devices. The term Noisy Intermediate-Scale Quantum (NISQ) was proposed by John Preskill in "Quantum Computing in the NISQ era and beyond", arXiv:1801.00862. Here, "noisy" may suggest that there is imperfect control over the qubits, and "intermediate-scale" may refer to the number of qubits, which may range from 50 to hundreds. Several physical systems made from superconducting qubits, artificial atoms, and ion traps have been proposed so far as viable candidates for building NISQ quantum devices and ultimately universal quantum computers.

[0078] In some cases, a classical simulator of the quantum circuit may be used, which may run on a classical computer such as a MacBook Pro laptop, a Windows laptop, or a Linux laptop. In some cases, the classical simulator may run on a cloud computing platform with access to multiple computing nodes in a parallel or distributed fashion. In some cases, all or a portion of the quantum mechanical energy and / or electronic structure calculations may be performed using a classical simulator.

[0079] The methods described herein may be implemented on an analog quantum simulator. The analog quantum simulator may be a quantum mechanical system consisting of multiple manufactured qubits. The analog quantum simulator may be designed to simulate a quantum system by using physically different but mathematically equivalent or nearly equivalent systems. In the analog quantum simulator, each qubit may be realized in an ion of a string of trapped atomic ions in a linear radio frequency trap. A bias source, called a local field bias, may be coupled to each qubit. The local field bias at the qubit may be programmable and controllable. In some cases, a qubit control system comprising a digital processing device is connected to the system of qubits and may program and adjust the local field bias at the qubit.

[0080] Classical Computers In some cases, the systems, media, networks, and methods described herein include a classical computer, or the use of a classical computer. In some cases, the classical computer includes one or more hardware central processing units (CPUs) that perform the functions of the classical computer. In some cases, the classical computer further comprises an operating system (OS) configured to execute executable instructions. In some cases, the classical computer is connected to a computer network. In some cases, the classical computer is connected to the Internet such that it accesses the World Wide Web. In some cases, the classical computer is connected to a cloud computing infrastructure. In some cases, the classical computer is connected to an intranet. In some cases, the classical computer is connected to a data storage device.

[0081] According to the description herein, suitable classical computers may include, by way of non-limiting examples, server computers, desktop computers, laptop computers, notebook computers, sub-notebook computers, netbook computers, netpad computers, set-top computers, media streaming devices, handheld computers, Internet appliances, mobile smartphones, tablet computers, personal digital assistants, home gaming consoles, and vehicles. Smartphones may be suitable for use with the methods and systems described herein. Select televisions, gaming consoles, and digital music players, possibly with computer network connectivity, may be suitable for use in the systems and methods described herein. Suitable tablet computers may include those having booklet, slate, and convertible configurations.

[0082] In some cases, a classical computer includes an operating system configured to execute executable instructions. An operating system may be software, including programs and data, that, for example, manages the device's hardware and provides services for the execution of applications. Suitable server operating systems include, by way of non-limiting examples, FreeBSD, OpenBSD, NetBSD®, Linux®, Apple® Mac OS X Server®, Oracle® Solaris®, Windows Server®, and Novell® NetWare®. Suitable personal computer operating systems may include, by way of non-limiting examples, UNIX-like operating systems, such as Microsoft® Windows®, Apple® Mac OS X®, UNIX®, and GNU / Linux®. In some cases, the operating system is provided by cloud computing. Suitable mobile smartphone operating systems may include, by way of non-limiting example, Nokia® Symbian® OS, Apple® iOS®, Research In Motion® BlackBerry® OS®, Google® Android®, Microsoft® Windows Phone® OS, Microsoft® Windows Mobile® OS, Linux®, and Palm® WebOS®. Suitable media streaming device operating systems may include, by way of non-limiting example, Apple TV®, Roku®, Boxee®, Google TV®, Google Chromecast®, Amazon Fire®, and Samsung® HomeSync®.Suitable home gaming console operating systems may include, by way of non-limiting example, Sony® PS3®, Sony® PS4®, Microsoft® Xbox 360®, Microsoft Xbox One, Nintendo® Wii®, Nintendo® Wii U®, and Ouya®.

[0083] In some cases, the classical computer includes a storage and / or memory device. In some cases, the storage and / or memory device is one or more physical devices used to temporarily or permanently store data or programs. In some cases, the device is volatile memory and requires power to maintain the stored information. In some cases, the device is non-volatile memory and retains the stored information when the classical computer is not powered. In some cases, the non-volatile memory includes flash memory. In some cases, the non-volatile memory includes dynamic random access memory (DRAM). In some cases, the non-volatile memory includes ferroelectric random access memory (FRAM). In some cases, the non-volatile memory includes phase change random access memory (PRAM). In other embodiments, the device is a storage device, including, by way of non-limiting example, CD-ROMs, DVDs, flash memory devices, magnetic disk drives, magnetic tape drives, optical disk drives, and cloud computing based storage. In some cases, the storage and / or memory device is a combination of devices such as those described herein.

[0084] In some cases, classical computers include a display to transmit visual information to a user. In some cases, the display is a cathode ray tube (CRT). In some cases, the display is a liquid crystal display (LCD). In some cases, the display is a thin film transistor liquid crystal display (TFT-LCD). In some cases, the display is an organic light emitting diode (OLED) display. In some cases, the OLED display is a passive matrix OLED (PMOLED) or an active matrix OLED (AMOLED) display. In some cases, the display is a plasma display. In other embodiments, the display is a video projector. In some cases, the display is a combination of devices such as those described herein.

[0085] In some cases, a classical computer includes an input device to receive information from a user. In some cases, the input device is a keyboard. In some cases, the input device is a pointing device, including, by way of non-limiting example, a mouse, a trackball, a track pad, a joystick, a game controller, or a stylus. In some cases, the input device is a touch screen or a multi-touch screen. In some cases, the input device is a microphone for capturing voice or other audio input. In some cases, the input device is a video camera or other sensor for capturing motion or visual input. In some cases, the input device is a Kinect, Leap Motion, or the like. In some cases, the input device is a combination of devices, such as those described herein.

[0086] Non-transitory computer-readable storage medium In some cases, the systems and methods described herein include one or more non-transitory computer-readable storage media encoded with a program including instructions executable by an operating system of an optionally networked digital processing device. In some cases, the computer-readable storage medium is a tangible component of a classical computer. In some cases, the computer-readable storage medium is optionally removable from the classical computer. In some cases, the computer-readable storage medium includes, by way of non-limiting examples, CD-ROMs, DVDs, flash memory devices, solid-state memory, magnetic disk drives, magnetic tape drives, optical disk drives, cloud computing systems and services, and the like. In some cases, the programs and instructions are encoded in the medium permanently, substantially permanently, semi-permanently, or non-transitoryly.

[0087] Embodiments of the disclosed methods for efficiently identifying stable conformations of chemical systems are described below.

[0088] Identification of target chemical systems A hybrid computing unit comprising a classical computer and a quantum computer can be used to perform quantum mechanical energy and / or electronic structure calculations on chemical systems. For example, such a hybrid computing unit can be used to perform a method for efficiently identifying stable conformations of a chemical system (e.g., a molecule).

[0089] FIG. 1 illustrates a flow chart for an example of a method 100 for providing an index to a classified list of conformers of a molecule using problem decomposition techniques on quantum computing hardware.

[0090] Method 100 may include obtaining an index of an input molecule according to operation 102. Method 100 disclosed herein may be applicable to any type of chemical system. Chemical systems may include, for example, organic compounds, inorganic compounds, polymers, peptides, polypeptides, proteins, nucleic acids, carbohydrates, etc. Methods disclosed herein may also be applicable to complexes of molecules, such as one or more protein-drug complexes (including or excluding solvent molecules).

[0091] Generation of conformer families The method 100 may include determining a set of conformations of the chemical system. For example, the method 100 may include generating a set (e.g., a list) of conformers for the input chemical system according to operation 104. A variety of different approaches may be used to enumerate conformers for the chemical system. In some cases, an exhaustive conformational sampler may be used, where conformations of the molecule are sampled by varying all dihedral angles around rotatable bonds in the chemical system. In some cases, Monte Carlo simulations or molecular dynamics simulations may be performed to generate the set of conformers. In another embodiment, the set of conformers of the molecule is as input along with the chemical system information.

[0092] Selection and manipulation of conformers from a group of conformers Method 100 may include selecting a conformer from a group or list (e.g., an ordered list) of conformers according to operation 106 and performing, for a compatible selection, at least any one, two, three, four, five, six, or seven, or up to any seven, six, five, four, three, two, or one, of operations 108, 110, 112, 114, 116, 118, and / or 120. At least any one, two, three, four, five, six, or seven, or up to any seven, six, five, four, three, two, or one, of operations 108, 110, 112, 114, 116, 118, and / or 120 may be performed for each conformer in the group of conformers.

[0093] (a) PD fragmentation of a chemical system Method 100 may include decomposing at least one conformation in the set into a plurality of molecular fragments. For example, method 100 may include decomposing the chemical system into a plurality of smaller fragments or subsystems (e.g., lists) (e.g., performing problem decomposition on the chemical system) according to operation 108. The particular scheme for decomposing the system into subsystems may vary depending on the PD technique used. In general, a suitable PD fragment ("fragment") size may be selected such that the computational resources required to process the fragments do not exceed the capabilities of the quantum-classical hardware to be used.

[0094] A variety of fragmentation techniques for chemical systems may be suitable for use, including, but not limited to, (i) divide-and-conquer (DC), (ii) fragment molecular orbital (FMO), (iii) density matrix embedding theory (DMET), (iv) density matrix renormalization group (DMRG), (v) tensor networks, (vi) incremental methods (as described herein with respect to FIG. 11), and the like.

[0095] For example, the FMO method was first described by Kitaura et al., "Fragment molecular orbital method: an approximate computational method for large molecules", Chemical Physics Letters, 1999, 313, 701, which is incorporated herein by reference in its entirety. The FMO method has been applied to many systems, such as those described by Fedorov et al., "Exploring chemistry with the fragment molecular orbital method", Physical Chemistry Chemical Physics, 2012, 14, 7562, which is incorporated herein by reference in its entirety.

[0096] For example, the DC method was first described by Yang, "Direct calculation of electron density in density-functional theory: Implementation for benzene and a tetrapeptide," Physical Review A, 1991, 44, 7823, which is incorporated herein by reference in its entirety. The DC method has been further developed and described by, for example, Akama et al., "Implementation of divide-and-conquer method including Hartree-Fock exchange interaction," Journal of Computational Chemistry, 2007, 28, 2003; and Kobayashi et al., "Divide-and-conquer approaches to quantum chemistry: Theory and implementation," in Linear-Scaling Techniques in Computational Chemistry and Physics: Methods and Applications, edited by Zalesny et al. (Springer Netherlands, Dordrecht, 2011), 97-127, each of which is incorporated herein by reference in its entirety.

[0097] For example, DMET was first described by Knizia et al., "Density Matrix Embedding: A Simple Alternative to Dynamical Mean-Field Theory," Physical Review Letters, 2012, 109, 186404, which is incorporated herein by reference in its entirety. DMET is further developed and described by Wouters et al., "A Practical Guide to Density Matrix Embedding Theory in Quantum Chemistry," Journal of Chemical Theory and Computation, 2016, 12, 2706, which is incorporated herein by reference in its entirety.

[0098] For example, DMRGs were first described by Steven R. White, "Density Matrix Formulation for Quantum Renormalization Groups," Physical Review Letters, 1992, 69, 2863, which is incorporated herein by reference in its entirety. A review of DMRGs is provided by Ulrich Schollwock, arxiv.org:cond-mat / 0409292 [cond-mat.str-el] or Review of Modern Physics, 2005, 77, 259, which is incorporated herein by reference in its entirety.

[0099] For example, tensor networks can be mathematical representations of quantum many-body states based on their entanglement structure. Different tensor network structures describe different physical situations, such as low-energy states of gapped 1D systems, 2D systems, and scale-invariant systems. Tensor networks can represent quantum states as one or more matrix product states. A review of tensor networks is provided by Roman Orus, "Tensor networks for complex quantum systems," Nature Reviews Physics, 2019, 1, 538, which is incorporated herein by reference in its entirety.

[0100] FIG. 11 illustrates a flow chart for an example of an incremental method 1100 for performing problem decomposition. The method 1100 may be generally referred to as an "incremental method" or may be variously referred to depending on the quantum chemistry method employed. For example, the utilization of a unitary coupled cluster (UCC) method in an incremental method may be referred to as an "incremental unitary coupled cluster" (iUCC) method. The method 1100 may include obtaining an index of the molecule according to operation 1102. The method 1100 disclosed herein may be applicable to any type of molecule. Molecules may include, for example, organic compounds, inorganic compounds, polymers, peptides, polypeptides, proteins, nucleic acids, carbohydrates, and the like. The method disclosed herein may also be applicable to complexes of molecules, such as one or more protein-drug complexes (including or excluding solvent molecules).

[0101] The method 1100 may include performing an incremental evolution of the energy of the molecule, according to operation 1104. The incremental evolution may be performed according to equation (1). E C =Σ i ε i +Σ i>j ε ij +Σ i>j>k ε ijk (1)

[0102] Here, the correlation energy E C(the difference between the total molecular energy and the mean-field Hartree-Fock energy) is expressed as an n-body Bethe-Goldstone expansion. The individual n-body correlation energy contributions are given by equation (2): ε i =E C (i) (2) ε ij =E C (ij)-ε i -ε j ε ijk =E C (ijk)-ε ij -ε ik -ε jk -ε i -ε j -ε k

[0103] Here, E C (i) is the individual one-body correlation energy, E C (ij) is the individual two-body correlation energy, E C (ijk) are the individual three-body correlation energies. The indices i, j, and k may correspond to any number of molecular orbitals, atoms, molecular fragments, or whole molecules. The indices i, j, and k may correspond to any potential combination of molecular orbitals, atoms, molecular fragments, or whole molecules. Thus, the incremental expansion may be expressed in terms of any potential combination of molecular orbitals, atoms, molecular fragments, or whole molecules.

[0104] FIG. 12 depicts the molecular orbitals, atoms, molecular fragments, and molecules used as the basis for the incremental method. When the incremental expansion is expressed in terms of any potential combination of atoms, fragments, and molecules, the resulting framework can be that of the FMO method. When the incremental expansion is expressed in terms of molecular orbitals, the resulting framework can be that of the incremental full configuration interaction (iFCI) method, such as that described by Zimmerman et al., "Strong Correlation in Incremental Full Configuration Interaction," Journal of Chemical Physics, 2017, 146, 224104, which is incorporated herein by reference in its entirety.

[0105] Returning to the description of FIG. 11, the method 1100 may further include, according to operation 1106, solving the Schrödinger equation for each increment described with respect to operation 1104 (e.g., by using a quantum chemistry simulator to solve the quantum chemistry problem according to method 200). In some cases, the solution of the Schrödinger equation may be achieved using a phase estimation procedure. For example, a phase estimation algorithm is described by Aspuru-Guzik et al., "Simulated Quantum Computation of Molecular Energies," Science, 2005, 309, 1704, which is incorporated herein by reference in its entirety. In some cases, the solution of the Schrödinger equation may be achieved using adiabatic quantum simulation. In some cases, the solution of the Schrödinger equation may be achieved by solving a Unitary Coupled Cluster (UCC) problem within a Variational Quantum Eigensolver (VQE). For example, VQE is described by McClean et al., "The theory of variational hybrid quantum-classical algorithms," New Journal of Physics, 2016, 18, 023023, which is incorporated herein by reference in its entirety. In some cases, the UCC ansatz may include every possible excitation for the increment. In such cases, the UCC ansatz may be equivalent to the exact solution of the Schrödinger equation, or to the full configuration interaction (FCI) for each increment. In some cases, truncation of the UCC ansatz to lower excitations may be used to approximate (to any possible approximation) the exact result for each increment. The solution of the Schrödinger equation may be repeated for one or more increments. For example, the solution of the Schrödinger equation may be repeated for any possible subset of all increments, or may be repeated for all increments. In some cases, the solution of the Schrödinger equation for each increment may be parallelized.For example, the solution of the Schrödinger equation for each increment can be parallelized using high performance computing architectures.

[0106] The method 1100 may further include calculating the quantum mechanical molecular electronic energy according to operation 1108. The quantum mechanical molecular electronic energy may be calculated by summing each of the incremental contributions according to equation (1) to obtain the quantum mechanical molecular correlation energy and thus the total quantum mechanical energy of the system under consideration.

[0107] Returning to the discussion of FIG. 1, in some cases, the same fragmentation scheme may be used for all conformers in a group. Such a fragmentation scheme may be appropriate, for example, when hardware capabilities are limited and fragment sizes must be very small. This fragmentation scheme may involve error compensation as described in the examples herein. In some cases, different fragmentation schemes may be used for one or more conformers in a group.

[0108] (b) Calculation of quantum mechanical energy and / or electronic structure for each PD fragment Method 100 may include determining, using a hybrid computing unit, a quantum mechanical energy and / or an electronic structure of each of at least a subset of the plurality of molecular fragments. For example, method 100 may include calculating the quantum mechanical energy and / or electronic structure of one or more subsystems according to acts 110, 112, and 114).

[0109] According to operation 110, the next fragment or subsystem in the list may be selected. Operations 112 and 114 may then be considered for each PD fragment. For example, according to operation 112, the quantum mechanical energy and / or electronic structure of the subsystem may be calculated (e.g., by using a quantum chemistry simulator to solve the quantum chemistry problem according to method 200). According to operation 114, the resulting quantum mechanical energy and / or electronic structure of the subsystem may be stored.

[0110] (i) PD molecular Hamiltonian construction In some cases, using the hybrid computing unit to determine the quantum mechanical energy and / or electronic structure of each of the molecular fragments may include determining the quantum mechanical energy and / or electronic structure of the molecular fragments (e.g., by constructing a molecular Hamiltonian or an electronic Hamiltonian), converting the quantum mechanical energy and / or electronic structure to an equivalent qubit energy and / or electronic structure (e.g., by converting fermion operators of the Hamiltonian to qubit operators), and using a quantum circuit to determine the quantum mechanical energy and / or electronic structure of the molecular fragments.

[0111] FIG. 2 illustrates a flow chart for an example of a method 200 for providing indications of quantum mechanical energy and / or electronic structure of a subsystem, as defined by problem decomposition techniques, on quantum computing hardware.

[0112] Method 200 may include obtaining an index of the subsystem according to operation 202. An approach for solving quantum chemistry problems using classical computing may be to use the Born-Oppenheimer approximation, where electronic and nuclear wave functions are decoupled and the electronic Hamiltonian is solved. However, methods disclosed herein may or may not use the Born-Oppenheimer approximation. Such a choice of using or not using the Born-Oppenheimer approximation may be selected according to operation 204, for example, by input from a user of the system.

[0113] If the Born-Oppenheimer approximation is selected by the user, an electronic Hamiltonian for the fragment may be constructed according to operation 206. If the Born-Oppenheimer approximation is not selected by the user, a molecular Hamiltonian for the fragment may be constructed according to operation 208.

[0114] The qubit Hamiltonian may be constructed fragment-by-fragment using either (a) a first quantization formalism according to operation 212, in which the space within the Hamiltonian is discretized by a lattice of qubits, or (b) a second quantization formalism according to operation 214, in which fermionic operators are converted to qubit operators. Such a choice of using either the first or second quantization formalism to generate the qubit Hamiltonian may be selected according to operation 210, e.g., by input from a user of the system.

[0115] For example, in the case of the second quantization formalism using the Born-Oppenheimer approximation according to act 206, the electron Hamiltonian H el teeth,

[0116]

number

[0117]

number

[0118]

number

[0119] The exact form of the molecular Hamiltonian can vary depending on the PD technique and the framework being used, such as Full Coordination Interaction (Full CI) or Coupled Cluster Theory (CC).

[0120] (ii) PD qubit Hamiltonian construction When the first quantization formalism is selected according to operation 212, the qubit Hamiltonian may be obtained by discretizing the three-dimensional real space into a three-dimensional lattice of qubits. Each lattice point may then be represented by a qubit variable.

[0121] According to operation 214, when the second quantization formalism is selected, the molecular Hamiltonian (which may be based on spin operators) may be transformed into a qubit Hamiltonian. The qubit Hamiltonian is a function of the {σ x ,σ y ,σ z}.

[0122] The transformation from spin to qubit Hamiltonian can be accomplished in a variety of ways, including but not limited to the Jordan-Wigner transformation or the Bravyi-Kitaev transformation.

[0123] For example, the Jordan-Wigner transformation gives the following qubit Hamiltonian:

[0124]

number

[0125] Where:

[0126]

number

[0127]

number

[0128] According to operation 216, time in the Hamiltonian may be discretized in preparation for performing a simulation of the Hamiltonian.

[0129] (iii) Circuit preparation A qubit Hamiltonian may be simulated according to operation 218. The qubit Hamiltonian may be simulated by performing at least one, two, three, or four, or up to four, three, two, or one, of operations 310, 312, 314, and 318 disclosed herein with respect to FIG.

[0130] An obstacle in the process of accurately distinguishing molecular conformations may be performing total quantum mechanical energy and / or electronic structure calculations. To help alleviate this obstacle, PD techniques may be used to decompose the problem into smaller, more manageable elements. In some cases, the total quantum mechanical energy and / or electronic structure calculations for each of the subsets of subproblems may be performed using a quantum computer. In some cases, the quantum computation process for each of the subsets of subproblems may be simulated on a classical computer. The process of calculating the total quantum mechanical energy and / or electronic structure using a quantum computer may include running a quantum algorithm to calculate the lowest eigenvalue of a Hamiltonian that describes the subproblems.

[0131] 3 illustrates a flow chart for an example of a method 218 for providing an indication of the expectation value of a Hamiltonian on quantum computing hardware. The method 218 may include operation 218 of method 200.

[0132] Method 218 may include, in accordance with operation 310, transferring the Hamiltonian to a quantum circuit that matches the characteristics (e.g., the connectivity of the qubits and which gates can be applied) of the computing system (e.g., quantum computing system or hardware, or quantum-classical system or hardware) being used. Techniques for computing the lowest energy eigenvalues ​​of the Hamiltonian may include phase estimation algorithms and variational quantum eigensolvers (VQE). For example, phase estimation algorithms are described by Aspuru-Guzik et al., "Simulated Quantum Computation of Molecular Energies," Science, 2005, 309, 1704, which is incorporated herein by reference in its entirety. For example, VQE is described by McClean et al., "The theory of variational hybrid quantum-classical algorithms," New Journal of Physics, 2016, 18, 023023, which is incorporated herein by reference in its entirety. These algorithms can be implemented to encode molecular or submolecular qubit Hamiltonians into the parameters of quantum circuits.

[0133] (iv) Preparation of the initial state for each PD fragment Method 218 may include preparing an initial state (or initial guess) for the quantum chemistry simulation on the quantum-classical hardware, according to operation 312. The suitable initial state may be a Hartree-Fock wave function. The suitable initial state may be a wave function obtained by a post-Hartree-Fock method. The suitable initial state may be prepared, for example, using any of the systems or methods described in Matsuura et al., “VanQver: The Variational and Adiabatically Navigated Quantum Eigensolver,” arXiv:1810.11511, Oct. 31, 2018, which is incorporated herein by reference in its entirety.

[0134] (v) Simulation of the PD Hamiltonian on quantum-classical hardware Given the quantum circuit (from act 310) and the initial state (from act 312), method 218 may include simulating a qubit Hamiltonian. Method 218 may include compiling and running (e.g., optimizing) the initial state and / or the qubit Hamiltonian on a quantum computer, according to act 314. In some cases, the quantum computer includes a quantum hardware device 316 or a classical simulator of the quantum circuit (e.g., quantum hardware simulator 316). For example, according to act 314, the transformed quantum circuit and the initial qubit state may be sent to quantum hardware device 316 or to quantum hardware simulator 316 to perform a quantum chemical simulation on a fragment-by-fragment basis.

[0135] The transmission of the circuit and initial states to calculate the total quantum mechanical energy and / or electronic structure of the fragments (according to operations 310 and 312, respectively) may occur serially as the transformed fragments are ready, or they may all be calculated and then transmitted in parallel to one or more quantum hardware devices or classical simulators.

[0136] In some cases, the quantum computer may include one or more of adiabatic quantum computers, quantum gate arrays, one-way quantum computers, topological quantum computers, quantum Turing machines, semiconductor-based quantum computers, trapped-ion quantum computers, trapped-atom quantum computers, optical lattices, quantum dot computers, spin-based quantum computers, space-based quantum computers, Loss-DiVincenzo quantum computers, nuclear magnetic resonance (NMR)-based quantum computers, solution NMR quantum computers, solid-state NMR quantum computers, solid-state NMR Kane quantum computers, helium surface electron quantum computers, resonator quantum electrodynamics-based quantum computers, molecular magnet quantum computers, fullerene-based quantum computers, linear optical quantum computers, diamond-based quantum computers, nitrogen-vacancy (NV) diamond-based quantum computers, Bose-Einstein condensate-based quantum computers, transistor-based quantum computers, and rare-earth metal ion-doped inorganic crystal-based quantum computers. The quantum computer may include one or more of a quantum annealer, an Ising solver, an optical parametric oscillator (OPO), and a gate model of quantum computing.

[0137] In some cases, a classical simulator of the quantum circuit may be used, which may run on a classical computer such as a MacBook Pro laptop, a Windows laptop, or a Linux laptop. In some cases, the classical simulator may run on a cloud computing platform with access to multiple computing nodes in a parallel or distributed manner. In some cases, the total quantum mechanical energy and / or electronic structure calculations for a subset of the fragments may be performed using a classical simulator, and the total quantum mechanical energy and / or electronic structure calculations for the remainder of the fragments may be performed using quantum hardware.

[0138] (vi) Measurement of the resulting state Method 218 may include measuring a qubit to provide a classical index of the lowest eigenvalue, according to operation 318. Depending on the algorithm used, the parameters required to provide the electronic structure shape that resulted in that lowest energy eigenvalue may also be provided. The basis of the measurement may be indicated by the Hamiltonian and the quantum algorithm being used. Measurements on the quantum data stored in the qubit may convert that information into classical bits of information. At least a portion of operation 218 may be repeated to provide an accurate estimate of the data being measured. In this case, multiple results obtained from multiple repeated executions of operation 218 may be averaged. Depending on the algorithm used, the parameters required to provide the electronic structure shape that resulted in that lowest energy eigenvalue may also be provided.

[0139] Returning to the discussion of FIG. 2, after operation 218 has been performed one or more times, method 200 may include measuring an indication of the expectation value of the Hamiltonian, according to operation 220.

[0140] Method 200 may include determining whether the Born-Oppenheimer approximation was used (e.g., in act 204), according to act 222. If the Born-Oppenheimer approximation was used, method 200 may include calculating a nucleus-nucleus repulsion energy and then adding the calculated nucleus-nucleus repulsion energy to the measured expectation value, according to act 224. Method 200 may include providing an indication of the quantum mechanical energy and / or electronic structure of the subsystem, thereby terminating the quantum chemical simulation performed by method 200, according to act 226.

[0141] Returning to the discussion of Figure 1, method 100 may include storing the resulting quantum mechanical energies and / or electronic structures of the subsystems, such as in a list of quantum mechanical subsystem energies and / or electronic structures, in accordance with operation 112. Method 100 may include determining whether all subsystems of the conformer have been processed to calculate their quantum mechanical energies and / or electronic structures, in accordance with operation 100, and if not, the next subsystem on the list may be selected (in accordance with operation 110) and operations 112 and 114 may be performed thereon.

[0142] (c) Combining the quantum mechanical energies and / or electronic structures of the PD fragments After one or more molecular fragments of a chemical system have been processed to calculate their quantum mechanical energy and / or electronic structure, a method for performing quantum mechanical energy and / or electronic structure calculations for a chemical system using a hybrid computing unit may include combining the quantum mechanical energies and / or electronic structures determined for the molecular fragments. For example, a method for efficiently identifying stable conformations of a chemical system may include recombining the energies and / or electronic structures obtained for each fragment to obtain a total quantum mechanical energy and / or electronic structure of the conformers of the entire chemical system (e.g., molecule) according to operation 118. The manner in which the recombination of the energies and / or electronic structures of the fragments to obtain a total quantum mechanical energy and / or electronic structure of the conformers of the chemical system depends on, and may be fully described by, the problem decomposition (PD) method used in operation 108. The resulting quantum mechanical energies and / or electronic structures of the conformers may then be stored, such as in a list of quantum mechanical conformer energies and / or electronic structures.

[0143] According to operation 120, it is determined whether all conformers of interest of a chemical system (e.g., a molecule) have been processed to calculate their quantum mechanical energy and / or electronic structure, and if not, the next conformer on the list is selected (according to operation 106) and operations 108, 110, 112, 114, 116, and 118 are performed on it.

[0144] In some cases, operations 108, 110, 112, 114, 116, and 118 are performed until a stopping criterion is satisfied. In one embodiment, the stopping criterion may be convergence of the electronic structure energy of each fragment. In another embodiment, the stopping criterion may be convergence of the molecular properties of each fragment (number of electrons in the fragment, reduced density matrix, etc.). Although the problem decomposition may not change the molecular properties being calculated, e.g., energy, in some cases it may be useful to vary the manner and type of problem decomposition until a stopping criterion is satisfied. For example, the number and / or size of the fragments may be iteratively changed. For example, block decimation may be iteratively changed in density matrix renormalization group techniques. For example, the fragment size may be changed from larger fragments to smaller fragments to increase the speed to convergence of a system with many fragments. For example, in a symmetric system, the number and / or size of the fragments may be changed until sufficient translational invariance between the fragments is found to increase the speed to convergence of the system by exploiting the symmetry.

[0145] Prediction of the most stable conformer After the conformers of interest within the group of conformers have been processed to calculate their quantum mechanical energies and / or electronic structures, the conformers provided within the group of conformers may be sorted in any order, such as sorted in order of least or most stable, according to operation 122, based on the estimate of the total quantum mechanical energy and / or electronic structure of each of the conformers provided by operation 118. According to operation 124, an index of a sorted list of conformers of the chemical system is provided based on the resulting quantum mechanical energies and / or electronic structures, which provides a prediction of the most stable conformer from within the group of conformers of the chemical system.

[0146] PD techniques can generally provide accurate results, as shown by works such as those described by Fedorov et al., "Exploring chemistry with the fragment molecular orbital method," Physical Chemistry Chemical Physics, 2012, 14, 7562; Kobayashi et al., "Divide-and-conquer approaches to quantum chemistry: Theory and implementation," in Linear-Scaling Techniques in Computational Chemistry and Physics: Methods and Applications, edited by Zalesny et al. (Springer Netherlands, Dordrecht, 2011), 97-127; and Wouters et al., "A Practical Guide to Density Matrix Embedding Theory in Quantum Chemistry," Journal of Chemical Theory and Computation, Journal of Chemical Theory and Computation, 2016, 12, 2706, each of which is incorporated herein by reference in its entirety.

[0147] In addition, the following examples illustrate the good correlation between the energies obtained by a particular method (e.g., by CCSD) with and without PD, even when the fragment sizes are very small. Thus, the most stable conformer of a chemical system can be directly identified based on the energies obtained by PD or by using the methods disclosed above to narrow the size of the conformer population for more accurate calculations.

[0148] first principles molecular dynamics The disclosed systems and methods may be used to simulate the evolution of molecular structures over time using first-principles molecular dynamics (AIMD) techniques. In such simulations, quantum-enabled problem decomposition (PD) techniques described herein (e.g., as described herein with respect to FIG. 1, FIG. 2, or FIG. 3) to calculate the quantum mechanical energy and / or electronic structure of the molecule. The quantum mechanical energy and / or electronic structure calculations may serve as the basis for force calculations in the AIMD framework. Forces on particles (such as one or more atoms in the molecule) within the molecule may be determined based on the quantum mechanical energy obtained by the quantum-enabled PD techniques described herein. The positions and velocities of the particles may then be updated using the AIMD techniques.

[0149] FIG. 13 illustrates a flow chart for an example of a method 1300 for performing first-principles molecular dynamics (AIMD) for molecules using problem decomposition techniques on quantum computing hardware.

[0150] Method 1300 may include obtaining an index of an input molecule according to operation 1302. Method 1300 disclosed herein may be applicable to any type of chemical system. Chemical systems may include, for example, organic compounds, inorganic compounds, polymers, peptides, polypeptides, proteins, nucleic acids, carbohydrates, etc. Methods disclosed herein may also be applicable to complexes of molecules, such as one or more protein-drug complexes (including or excluding solvent molecules).

[0151] The method 1300 may include obtaining initial coordinates of the particles in the system according to operation 1304. The initial coordinates of the particles in the system may correspond to, for example, coordinates of atomic nuclei in a molecule. The initial coordinates of the particles in the system may be theoretically derived or experimentally derived. For example, the initial coordinates of the particles in the system may be derived from a predicted molecular structure. The initial coordinates of the particles in the system may be derived from an experimental procedure, such as X-ray crystallography, transmission electron microscopy (TEM), scanning electron microscopy (SEM), scanning tunneling microscopy (STEM), atomic force microscopy (AFM), solution nuclear magnetic resonance (NMR), solid-state NMR, or other experimental procedure. The initial coordinates of the particles in the system may be obtained from a database, such as PubChem, Chemical Entities of Biological Interest (ChEBI), DrugBank, Small Molecule Pathway Database (SMPDB), ChemDB, Protein Data Bank (PDB), or other database.

[0152] The method 1300 may include obtaining an initial velocity of the particles in the system, according to operation 1306. The initial velocity of the particles may be obtained in a variety of manners. For example, the initial velocity of the particles may be obtained by randomly selecting a velocity for each particle from a Maxwell-Boltzmann distribution at a temperature. In some cases, such a procedure may result in a net momentum of the system, resulting in an initial linear motion of the system as a whole. In some cases, the initial linear motion may be eliminated by calculating the net momentum of the system and adjusting the initial velocity of each particle to reduce the net momentum to zero. Similarly, the procedure may result in a net angular momentum of the system, resulting in an initial rotational motion of the system as a whole. In some cases, the initial rotational motion may be eliminated by calculating the net angular momentum of the system and adjusting the initial angular velocity of each particle to reduce the net angular momentum to zero.

[0153] Additional parameters may be set during either operations 1304 or 1306. For example, a target number of molecular dynamics time steps, a time increment, a target temperature, and / or a target pressure may be specified.

[0154] The method 1300 may include calculating a force on each particle in the system, according to operation 1308.

[0155] 14 illustrates a flow chart for an example method 1400 for calculating forces on each particle of a system in an ab initio molecular dynamics (AIMD) simulation. The method may include implementing a method for quantum-enabled PD for quantum mechanical energy and / or electronic structure calculations of the system, such as method 100 described herein.

[0156] The method may further include estimating a force on each particle of the system according to operation 1402. The force on each particle of the system may be calculated from quantum mechanical energy and / or electronic structure calculations of the system. The force on each particle of the system may be calculated by various procedures. For example, the force on each particle of the system may be calculated using Jordan's quantum algorithm for numerical gradient estimation as disclosed in Jordan, "Fast Quantum Algorithm for Numerical Gradient Estimation," Physical Review Letters, 2015, 95, 050501, which is incorporated herein by reference in its entirety. Jordan's quantum algorithm for numerical gradient estimation may be implemented using quantum hardware (such as a quantum computer described herein) or on a quantum simulator (such as a quantum simulator described herein). The force on each particle of the system may be calculated using numerical gradient estimation techniques on classical hardware (such as a classical computer described herein).

[0157] Returning to the discussion of FIG. 13, the method 1300 may further include updating the coordinates and / or velocities of the particles in the system for the next time step, according to operation 1310. The coordinates and / or velocities of the particles in the system may be updated according to various procedures. For example, the coordinates and / or velocities of the particles in the system may be updated using a Bere procedure, where the coordinates and / or velocities are updated using a series expansion of the coordinates and / or velocities based on the most recent or second most recent time step. The coordinates and / or velocities of the particles in the system may be updated using a velocity Bere procedure, where the position is updated based on the most recent velocity, the velocity is partially updated based on the most recent forces, the updated forces are calculated using the updated position, and the velocity is fully updated based on the most recent forces. The coordinates and / or velocities of the particles in the system may be updated by numerical integration of the forces using various integration techniques, such as symplectic integration, Bere-Stormer integration, Runge-Kutta integration, Biemann integration, or other integration techniques. During the updating of the coordinates and / or velocities of the particles in the system, various thermostats and / or barostats can be applied to maintain control of the temperature and / or pressure of the system. For example, a Langevin thermostat and an Anderson barostat can be applied.

[0158] The method 1300 may include storing coordinates and / or velocities of particles in the system, such as in a list of coordinates and / or velocities, according to operation 1312. The list of coordinates and / or velocities may include a trajectory of the system.

[0159] Method 1300 may include exploring the number of molecular dynamics time steps, according to operation 1314. Any one or more of operations 1302, 1304, 1306, 1308, 1310, and 1312 may be repeated until a number of time steps is met, such as a threshold, a predefined number of steps, etc. At such point, method 1300 may be aborted.

[0160] The method 1300 may include providing an indication of a resulting trajectory of the system.

[0161] Distributed Computing In some cases, the quantum mechanical energy and / or electronic structure of each of a subset of the plurality of molecular fragments may be determined or calculated using one or more distributed computing systems, such as one or more cluster or cloud-based computing systems. The distributed computing system may comprise a plurality of non-classical computers (such as any non-classical computers described herein), a plurality of classical computers, or both. For example, the distributed computing system may comprise a plurality of non-classical computers. Each non-classical computer of the plurality of non-classical computers may be assigned (e.g., through a scheduling routine) to determine or calculate the quantum mechanical energy and / or electronic structure of one or more molecular fragments of the plurality of molecular fragments. Each non-classical computer may be configured to determine or calculate the quantum mechanical energy and / or electronic structure of one or more molecular fragments assigned to the non-classical computer in parallel with the determination or calculation of the quantum mechanical energy and / or electronic structure of other molecular fragments assigned to other non-classical computers of the plurality of non-classical computers. In this manner, the determination of quantum mechanical molecular electronic energies and / or molecular electronic structures may be significantly sped up.

[0162] The distributed computing system has at least about 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 200, 300, 400, 500, 600, 700, 800, 900, 1,000, 2,000, 3,000, 4,000, 5,000, 6,000, 7,000, 8,000, 9,000, 10,000, In one embodiment, the present invention may comprise 00, 20,000, 30,000, 40,000, 50,000, 60,000, 70,000, 80,000, 90,000, 100,000, 200,000, 300,000, 400,000, 500,000, 600,000, 700,000, 800,000, 900,000, 1,000,000 or more non-classical computers. Distributed computing systems can reach up to approximately 1,000,000, 900,000, 800,000, 700,000, 600,000, 500,000, 400,000, 300,000, 200,000, 100,000, 90,000, 80,000, 70,000, 60,000, 50,000, 40,000, 30,000, 20,000, The distributed computing system may comprise 10,000, 9,000, 8,000, 7,000, 6,000, 5,000, 4,000, 3,000, 2,000, 1,000, 900, 800, 700, 600, 500, 400, 300, 200, 100, 90, 80, 70, 60, 50, 40, 30, 20, 10, 9, 8, 7, 6, 5, 4, 3, 2, or 1 non-classical computer. The distributed computing system may comprise a number of non-classical computers that are within a range defined by any two of the preceding values.

[0163] Each non-classical computer of the plurality of non-classical computers has at least about 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 200, 300, 400, 500, 600, 700, 800, 900, 1,000, 2,000, 3,000, 4,000, 5,000, 6,000, 7,000, 8,000, 9,000, 10,000, 20,000, 30,000, 40,000, 50,000, 60,000, 70,000, 80,000, 90,000, 100,000, 200,000, 300,000, 400,000, 500,000, 600,000, 700,000, 800,000, 900,000, 1,000,000, 2,000,000, 3,000,000, 4, 000,000, 5,000,000, 6,000,000, 7,000,000, 8,000,000, 9,000,000, 10,000,000, 20,000,000, 30,000,000, 40,000,000, 50,000,000, 60,000,000, 70,000,000, 80,000,000, 90,000,000, 100,000,00 The quantum mechanical energy and / or electronic structure of 0, 200,000,000, 300,000,000, 400,000,000, 500,000,000, 600,000,000, 700,000,000, 800,000,000, 900,000,000, 1,000,000,000, or more fragments may be determined or calculated.Each non-classical computer has a maximum of approximately 1,000,000,000, 900,000,000, 80,000,000, 7000,000,000, 600,000,000, 500,000,000, 400,000,000, 300,000,000, 200,000,000, 100,000,000, 90,000,000, 80,000,000, 70 ,000,000, 60,000,000, 50,000,000, 40,000,000, 30,000,000, 20,000,000, 10,000,000, 9,000,000, 8,000,000, 7,000,000, 6,000,000, 5,000,000, 4,000,000, 3,000,000, 2,000,000, 1,000,0 00, 900,000, 800,000, 700,000, 600,000, 500,000, 400,000, 300,000, 200,000, 100,000, 90,000, 80,000, 70,000, 60,000, 50,000, 40,000, 30,000, 20,000, 10,000, 9,000, 8,000, 7,000, 6,000, The non-classical computers may be configured to determine or calculate the quantum mechanical energies and / or electronic structures of 5,000, 4,000, 3,000, 2,000, 1,000, 900, 800, 700, 600, 500, 400, 300, 200, 100, 90, 80, 70, 60, 50, 40, 30, 20, 10, 9, 8, 7, 6, 5, 4, 3, 2, or 1 fragments. Each non-classical computer may be configured to determine or calculate the quantum mechanical energies and / or electronic structures of several fragments that are within a range defined by any two of the preceding values.

[0164] FIG. 16 illustrates a flow chart for an example of a method 1600 for performing quantum mechanical energy or electronic structure calculations on a chemical system using a distributed computing system.

[0165] Method 1600 may include one or more operations described herein with respect to method 100 of FIG. 1. For example, method 1600 may include obtaining an index for the molecule (operation 102 as described herein with respect to method 100 of FIG. 1, not shown in FIG. 16). Method 1600 may include generating or obtaining a list of conformers for the molecule (operation 104 as described herein with respect to method 100 of FIG. 1, not shown in FIG. 16). Method 1600 may include selecting a next conformer in the list (operation 106 as described herein with respect to method 100 of FIG. 1, not shown in FIG. 16). Method 1600 may include performing a problem decomposition to populate the list of conformer subsystems (operation 108 as described herein with respect to method 100 of FIG. 1).

[0166] 16, the method 1600 can include assigning each subset of a conformer to a job scheduler. Each subset of a conformer can include any molecular fragment described herein. The job scheduler can then assign each subset of a conformer to a distributed computing system described herein.

[0167] Method 1600 may include calculating the energy and / or electronic structure of each subsystem using a distributed computing system (operation 112 as described herein with respect to method 100 of FIG. 1, not shown in FIG. 16 ). Method 1600 may include storing the resulting energy and / or electronic structure of each subsystem (operation 114 as described herein with respect to method 100 of FIG. 1, not shown in FIG. 16 ). Method 1600 may include determining whether all subsystems of a conformer have been evaluated (operation 116 as described herein with respect to method 100 of FIG. 1, not shown in FIG. 16 ). Method 1600 may include combining all subsystem energies and / or electronic structures into a conformer energy and / or electronic structure (operation 118 as described herein with respect to method 100 of FIG. 1, not shown in FIG. 16 ). Method 1600 may include determining whether all conformers of the molecule have been evaluated (operation 120 as described herein with respect to method 100 of FIG. 1, not shown in FIG. 16). Method 1600 may include classifying the conformations based on the resulting energy and / or electronic structure of each conformer (operation 122 as described herein with respect to method 100 of FIG. 1, not shown in FIG. 16). Method 1600 may include providing an indication of a classified list of conformers of the molecule based on the resulting energy and / or electronic structure (operation 124 as described herein with respect to method 100 of FIG. 1, not shown in FIG. 16).

[0168] In some cases, one or more distributed computing systems may be utilized to perform one or more operations of the methods described herein. For example, a distributed computing system may be utilized to perform one or more of operations 102, 104, 106, 108, 110, 112, 114, 116, 118, 120, 122, and 124 of method 100 described herein with respect to FIG. 1, one or more of operations 202, 204, 206, 208, 210, 212, 214, 216, 218, 220, 222, 224, or 226 of method 200 described herein with respect to FIG. 2, one or more of operations 300 described herein with respect to FIG. 11, one or more of operations 1302, 1304, 1306, 1308, 1310, 1312, 1314, or 1316 of method 1300 described herein with respect to FIG. 13, or operation 1402 of method 1400 described herein with respect to FIG. 14. The distributed computing system may comprise a plurality of classical computers (such as any classical computers described herein). Each classical computer of the plurality of classical computers may be assigned (e.g., via a scheduling routine) to perform any one or more of the operations described herein with reference to the distributed computing system for one or more of the plurality of molecular fragments. Each classical computer may be configured to perform any one or more of the operations for one or more molecular fragments assigned to it in parallel with the performance of one or more operations for other molecular fragments assigned to other classical computers of the plurality of classical computers. In this manner, execution of operations may be significantly sped up.

[0169] The distributed computing system has at least about 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 200, 300, 400, 500, 600, 700, 800, 900, 1,000, 2,000, 3,000, 4,000, 5,000, 6,000, 7,000, 8,000, 9,000, 10,000, In one embodiment, the present invention may comprise 00, 20,000, 30,000, 40,000, 50,000, 60,000, 70,000, 80,000, 90,000, 100,000, 200,000, 300,000, 400,000, 500,000, 600,000, 700,000, 800,000, 900,000, 1,000,000 or more classical computers. Distributed computing systems can reach up to approximately 1,000,000, 900,000, 800,000, 700,000, 600,000, 500,000, 400,000, 300,000, 200,000, 100,000, 90,000, 80,000, 70,000, 60,000, 50,000, 40,000, 30,000, 20,000, The distributed computing system may comprise 10,000, 9,000, 8,000, 7,000, 6,000, 5,000, 4,000, 3,000, 2,000, 1,000, 900, 800, 700, 600, 500, 400, 300, 200, 100, 90, 80, 70, 60, 50, 40, 30, 20, 10, 9, 8, 7, 6, 5, 4, 3, 2, or 1 classical computer. The distributed computing system may comprise a number of classical computers that are within a range defined by any two of the preceding values.

[0170] Each classical computer of the plurality of classical computers has at least about 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 200, 300, 400, 500, 600, 700, 800, 900, 1,000, 2,000, 3,000, 4,000, 5,000, 6,000, 7,000, 8,000, 9,000, 10,000, 2 0,000, 30,000, 40,000, 50,000, 60,000, 70,000, 80,000, 90,000, 100,000, 200,000, 300,000, 400,000, 500,000, 600,000, 700,000, 800,000, 900,000, 1,000,000, 2,000,000, 3,000,000, 4,000,000, 5,000, 000, 6,000,000, 7,000,000, 8,000,000, 9,000,000, 10,000,000, 20,000,000, 30,000,000, 40,000,000, 50,000,000, 60,000,000, 70,000,000, 80,000,000, 90,000,000, 100,000,000, 200,000,000, 300,0 In one embodiment, the present invention may be configured to perform any one or more of the operations described herein with reference to a distributed computing system for up to 100,000, 400,000,000, 500,000,000, 600,000,000, 700,000,000, 800,000,000, 900,000,000, 1,000,000,000 or more fragments.Each classical computer can run up to 1,000,000,000, 900,000,000, 80,000,000, 7000,000,000, 600,000,000, 500,000,000, 400,000,000, 300,000,000, 200,000,000, 100,000,000, 90,000,000, 80,000,000, 70,000,000 0, 60,000,000, 50,000,000, 40,000,000, 30,000,000, 20,000,000, 10,000,000, 9,000,000, 8,000,000, 7,000,000, 6,000,000, 5,000,000, 4,000,000, 3,000,000, 2,000,000, 1,000,000, 900,000, 80 0,000, 700,000, 600,000, 500,000, 400,000, 300,000, 200,000, 100,000, 90,000, 80,000, 70,000, 60,000, 50,000, 40,000, 30,000, 20,000, 10,000, 9,000, 8,000, 7,000, 6,000, 5,000, 4,000, 3,000, 2 Each classical computer may be configured to perform any one or more of the operations described herein with reference to a distributed computing system for a number of fragments that are within a range defined by any two of the preceding values.

[0171] In some cases, a distributed computing system may comprise multiple non-classical computers (such as any non-classical computers described herein) and multiple classical computers (such as any classical computers described herein). In some cases, a problem (such as a quantum chemistry problem or simulation) may be solved using a distributed computing system comprising various types or combinations of systems, such as, for example, one or more classical computers, one or more non-classical computers (such as one or more quantum computers), or a combination of one or more classical computers and one or more non-classical computers. For example, FIG. 15 illustrates an example of a system or combination of systems that may be used to solve a problem, such as a quantum chemistry problem or simulation.

[0172] Figure 15 illustrates various types of application layer preprocessing methods and approximations that can be used to solve the problem. On the right, the computation can be performed entirely on a classical computing system. On the left, the computation can be performed entirely on a non-classical computing system (e.g., gate model quantum hardware). In between, various types of quantum simulators and quantum emulators can be used to perform the computation.

[0173] For example, the methods described herein may be implemented on an analog quantum simulator (e.g., a gate-model quantum simulator). The analog quantum simulator may be a quantum mechanical system consisting of a plurality of fabricated qubits. The analog quantum simulator may be designed to simulate a quantum system by using physically different but mathematically equivalent or nearly equivalent systems. For example, each qubit may be realized in an ion of a string of trapped atomic ions in a linear radio frequency trap. A bias source, referred to as a local field bias, may be coupled to each qubit. The local field bias on the qubit may be programmable and controllable. In some cases, a qubit control system comprising a digital processing device is connected to the system of qubits and may program and adjust the local field bias on the qubits.

[0174] An analog quantum simulator may include a set of gates and connections that may be implemented natively on the hardware. An analog quantum simulator may include a set of gates and connections that may not be implemented natively on the hardware (non-native gates). An analog quantum simulator may use a combination of ancillary qubits (e.g., ancilla qubits) and native gates to simulate the behavior of the non-native gates. The problem to be solved may utilize one or more of a qubitized Hamiltonian, a quantum algorithm layer, a circuit compiler and / or optimizer, or an interface to a hardware or simulator backend.

[0175] In some cases, all or a portion of the quantum mechanical energy and / or electronic structure calculations may be performed using a classical simulator (e.g., a classical emulator). The classical simulator may run on a classical computer, such as a MacBook Pro laptop, a Windows laptop, or a Linux laptop. In some cases, the classical simulator may run on a cloud computing platform with access to multiple computing nodes in a parallel or distributed manner, as described herein with respect to Figures 18, 19, 20, and 21.

[0176] In some cases, a classical simulator of a quantum circuit may be used that can simulate a circuit layer of a computation on quantum hardware (e.g., a classical circuit layer emulator). The problem to be solved may utilize one or more of a qubitized Hamiltonian, a quantum algorithm layer, or a circuit compiler, a circuit optimizer, or both. The classical circuit layer emulator may enable testing, prototyping, etc. of quantum machine code without the use of expensive quantum hardware. The classical circuit layer emulator may emulate the gate operation of a gate-model quantum computer, or an analog quantum simulator, or both.

[0177] In some cases, a classical simulator of a quantum algorithm may be used that can simulate the circuit compilation and circuit layers of the computation on quantum hardware (e.g., a classical algorithm emulator). The problem to be solved may utilize one or more of the qubitized Hamiltonians, or the quantum algorithm layers. The classical algorithm emulator may enable testing, prototyping, etc. of the quantum machine code without the use of expensive quantum hardware. The classical algorithm emulator may emulate the gate operations of a gate model quantum computer, or an analog quantum simulator, or both. The classical algorithm emulator may emulate the circuit compilation and optimization of a gate model quantum computer, or an analog quantum simulator, or both.

[0178] In some cases, quantum mechanical energy and / or electronic structure calculations can be performed using classical simulators that can simulate quantum mechanical systems without converting the fermion Hamiltonian to a qubit Hamiltonian (e.g., fermion quantum emulators).

[0179] FIG. 17 illustrates an example of a distributed computing system architecture that includes a non-classical or quantum computer (QC1) and multiple classical computers (CL1 through CLK, where K is an integer). The distributed computing system may include at least one classical computer and at least one non-classical computer. In some cases, the distributed computing system may include all classical computers. In some cases, the distributed computing system may have all quantum computers. The distributed computing system may include (i) at least or multiple (e.g., at least 2, 3, 4, 5, 6, 7, 8, 9, 10, or more) non-classical computers (e.g., quantum computers), and (ii) at least or multiple (e.g., at least 2, 3, 4, 5, 6, 7, 8, 9, 10, or more) classical computers.

[0180] The distributed computing system may be managed by a scheduler. The scheduler may adjust one or more parameters of at least one of the one or more subproblems. The scheduler may identify computing resources in various nodes of the distributed computing network, e.g., one or more non-classical computers, one or more classical computers, one or more virtual machines, one or more cloud-based machines, one or more work stations, one or more supercomputing nodes, one or more servers, etc. The scheduler may organize the subproblems, prioritize the subproblems, distribute the problems to the various computing resources, etc.

[0181] 18, 19, 20, and 21 illustrate various examples of organizations of a distributed computing system of the present disclosure. In some cases, the distributed computing system of any of FIG. 18, 19, 20, or 21 may be used to perform one or more operations of the methods described herein. For example, the distributed computing system of any of FIG. 18, 19, 20, or 21 may perform one or more of operations 102, 104, 106, 108, 110, 112, 114, 116, 118, 120, 122, and 124 of the method 100 described herein with respect to FIG. 1, one or more of operations 202, 204, 206, 208, 210, 212, 214, 216, 218, 220, 222, 224, or 226 of the method 200 described herein with respect to FIG. 2, one or more of operations 202, 204, 206, 208, 210, 212, 214, 216, 218, 220, 222, 224, or 226 of the method 200 described herein with respect to FIG. 13, or operation 1402 of a method 1400 described herein with respect to FIG. 14.

[0182] 18 illustrates a distributed computing system including sequential problem decomposition, according to some embodiments. Quantum-enabled problem decomposition (PD) techniques described herein for computing the quantum mechanical energy and / or electronic structure of a molecule (e.g., as described herein with respect to FIG. 1, FIG. 2, or FIG. 3) may include embodiments, variations, or examples of sequential problem decomposition techniques that may be implemented on the distributed computing system 1800 of FIG.

[0183] The target 1810 may include a chemical system. The chemical system may include, for example, a molecule, a portion of a molecule, a fragment, an aggregate, etc. The target 1810 may be decomposed into one or more fragments. For example, the target 1810 may be decomposed into fragments 1811, 1812, 1813, 1814, and 1815. FIG. 18 illustrates a problem decomposition device 1820. The problem decomposition device may include any of the problem decomposition methods and / or techniques disclosed herein. After decomposition, the fragment 1811 may be encoded as a fragment 1811′ on an electronic structure solver 1830. The electronic structure solver 1830 may comprise a non-classical or quantum computing system as described herein. The electronic structure solver 1830 may calculate the energy E′ of the fragment 1811′. The electronic structure solver 1830 may pass the energy E' to the problem decomposer 1820, and the fragments 1812 may be passed to the electronic structure solver 1830. Each fragment of the target 1810 may be passed sequentially to the electronic structure solver 1830. All or a portion of the multiple fragments may be passed to the electronic structure solver 1830. The multiple fragments may be passed to the electronic structure solver 1830 in any order.

[0184] The problem decomposer 1820 and the electronic structure solver 1830 may comprise portions of a distributed computing system. For example, the problem decomposer 1820 may be local to a user (e.g., on the same machine used by the user, in the same physical location on a separate machine, etc.) and the electronic structure solver 1830 may be remote to the user. For example, the problem decomposer 1820 may comprise portions of a client-side library. For example, the electronic structure solver 1830 may comprise remote endpoints. For example, the problem decomposer 1820 may be remote to a user at a first remote endpoint and the electronic structure solver 1830 may be remote to a user at a second remote endpoint. The first remote endpoint and the second remote endpoint may be the same endpoint. The first remote endpoint and the second remote endpoint may be remote to each other. For example, the first remote endpoint may comprise a remote server and the second remote endpoint may comprise a non-classical computer. The remote endpoint may include a classical computing system, a non-classical computing system, or a hybrid computing unit as disclosed herein.

[0185] 18 shows five fragments and solutions, but the methods of the present disclosure may be used with any number of fragments and solvers. For example, the distributed computing system of FIG. 18 may have at least about 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 200, 300, 400, 500, 600, 700, 800, 900, 1,000, 2,000, 3,000, 4,000, 5,000, 6,000, 7,000, 8,000, 9,000, 1 0,000, 20,000, 30,000, 40,000, 50,000, 60,000, 70,000, 80,000, 90,000, 100,000, 200,000, 300,000, 400,000, 500,000, 600,000, 700,000, 800,000, 900,000, 1,000,000, 2,000,000, 3,000,000, 4,00 0,000, 5,000,000, 6,000,000, 7,000,000, 8,000,000, 9,000,000, 10,000,000, 20,000,000, 30,000,000, 40,000,000, 50,000,000, 60,000,000, 70,000,000, 80,000,000, 90,000,000, 100,000,000, 2 In one embodiment, the present invention may be configured to perform any one or more of the operations described herein for 00,000,000, 300,000,000, 400,000,000, 500,000,000, 600,000,000, 700,000,000, 800,000,000, 900,000,000, 1,000,000,000 or more fragments.For example, the distributed computing system in Figure 18 can achieve up to approximately 1,000,000,000, 900,000,000, 80,000,000, 7000,000,000, 600,000,000, 500,000,000, 400,000,000, 300,000,000, 200,000,000, 100,000,000, 90,000,000, 80, 000,000, 70,000,000, 60,000,000, 50,000,000, 40,000,000, 30,000,000, 20,000,000, 10,000,000, 9,000,000, 8,000,000, 7,000,000, 6,000,000, 5,000,000, 4,000,000, 3,000,000, 2,000,000, 1, 000,000, 900,000, 800,000, 700,000, 600,000, 500,000, 400,000, 300,000, 200,000, 100,000, 90,000, 80,000, 70,000, 60,000, 50,000, 40,000, 30,000, 20,000, 10,000, 9,000, 8,000, 7,000, 6,000 , 5,000, 4,000, 3,000, 2,000, 1,000, 900, 800, 700, 600, 500, 400, 300, 200, 100, 90, 80, 70, 60, 50, 40, 30, 20, 10, 9, 8, 7, 6, 5, 4, 3, 2, or 1 fragment. For example, the distributed computing system of FIG. 18 may be configured to perform any one or more of the operations described herein for a number of fragments that are within a range defined by any two of the preceding values.

[0186] In some cases, a method including sequential problem decomposition may include generating an instance of an electronic structure solver, generating each of a plurality of fragments, and sequentially solving the electronic structure of each fragment with the electronic structure solver to generate an energy for each fragment. Sequential problem decomposition may be improved by parallelization of the problem.

[0187] Parallelization of a problem may be facilitated by problem decomposition and subsequent distribution of problem subsystems across one or more nodes in a high performance computer, which may include one or more non-classical computers. Figures 18, 19, 20, and 21 show examples of distributed computing systems that may implement the methods described herein for performing quantum mechanical energy or electronic structure calculations for chemical systems using a distributed computing system (e.g., as described herein with respect to Figure 16). Method 1600 may include one or more operations of the quantum-enabled problem decomposition (PD) techniques described herein (e.g., as described herein with respect to Figure 1, Figure 2, or Figure 3) for calculating the quantum mechanical energy and / or electronic structure of a molecule. Steps of the methods of Figures 1, 2, 3, 11, 13, or 14 may be implemented on a distributed computing system comprising multiple non-classical computers (such as any non-classical computers described herein) and multiple classical computers (such as any classical computers described herein).

[0188] FIG. 19 illustrates a distributed computing system 1900 with problem dispatch, according to some embodiments. Problem dispatch may be used to distribute problem subsystems to one or more nodes. For example, problem dispatch may provide logic that controls the distribution of the subsystems to the nodes and returns results from the nodes. Problem dispatch may control whether each subsystem of the problem is solved locally (e.g., on the same machine, in the same physical location on a separate machine, etc.), on the cloud, or in a high performance computing cluster.

[0189] The problem decomposition may create subsystems, e.g., fragments, and data structures for each fragment. For example, the data structures may include the type of electronic structure solver and / or parameters to be passed to the solver. In some cases, the data structures may pass input from a user. In some cases, the data structures may pass parameters based on the input problem. The problem dispatch may handle the creation and / or implementation of the electronic structure solver with parameters, such as parameters from a dictionary. The problem dispatch may return output from the solver, e.g., energy. A simple problem dispatch may create and / or implement each solver sequentially. In some examples, the problem dispatch may use a multiprocessing package to create and / or implement each solver using a parallel scheme, e.g., Python's multiprocessing package.

[0190] The target 1910 may include a chemical system. The chemical system may include, for example, a molecule, a portion of a molecule, a fragment, an assembly, etc. The target 1910 may be decomposed into one or more fragments. For example, the target 1910 may be decomposed into fragments 1911, 1912, 1913, 1914, and 1915 in a problem decomposition device. Figure 19 illustrates a problem decomposition device 1920. The problem decomposition device may include instructions for performing any of the problem decomposition methods and / or techniques disclosed herein.

[0191] After decomposition, the fragments may be distributed by problem dispatch 1940 to one or more electronic structure solvers 1931, 1932, 1933, 1934, and 1935. Problem dispatch 1940 may create and / or implement the electronic structure solvers. Problem dispatch 1940 may pass parameters to the solvers, such as input parameters from a user. Problem dispatch may return outputs from the solvers, such as, for example, energies.

[0192] Fragments 1911, 1912, 1913, 1914, and 1915 may be encoded as fragments 1911', 1912', 1913', 1914', and 1915' on electronic structure solvers 1931, 1932, 1933, 1934, and 1935. Electronic structure solvers 1931, 1932, 1933, 1934, and 1935 may comprise one or more non-classical computing systems, one or more quantum computing systems, or one or more hybrid computing units, as described herein. Electronic structure solvers 1931, 1932, 1933, 1934, and 1935 calculate energies E for fragments 1911', 1912', 1913', 1914', and 1915', respectively. 1 , E 2 , E 3 , E 4 , and E 5 The electronic structure solvers 1931, 1932, 1933, 1934, and 1935 can calculate the energy E 1 , E 2 , E 3 , E 4 , and E 5 to problem dispatch 1940. The electronic structure solver may receive and return fragments from problem dispatch.

[0193] All or a portion of the multiple fragments may be handed off to the electronic structure solver. The multiple fragments may be handed off to the electronic structure solver in any order, serially, or in a parallel manner. The problem decomposer 1920, problem dispatcher 1930, and electronic structure solvers 1931, 1932, 1933, 1934, and 1935 may comprise portions of a distributed computing system.

[0194] The distributed computing system with problem dispatch of FIG. 19 may be implemented in a variety of ways. For example, some code may be executed by a “client-side” digital computing device, e.g., a user's digital computing device, and other code may be executed in a portion of a distributed computing system, which may include one or more remote endpoints. The distributed computing system may include one or more cluster or cloud-based computing systems. The distributed computing system may include multiple non-classical computers (such as any non-classical computers described herein), multiple classical computers, or both. The distributed computing system may include one or more endpoints. The one or more endpoints may include a Representational State Transfer (REST) ​​endpoint. The one or more endpoints may execute responses from a client-side library. The REST calls may be executed on one or more classical, hybrid, or non-classical computing devices connected through a network, such as a distributed computing system, e.g., a cloud network.

[0195] 19 shows five fragments and solutions, but the methods of the present disclosure may be used with any number of fragments and solvers. For example, the distributed computing system of FIG. 19 may have at least about 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 200, 300, 400, 500, 600, 700, 800, 900, 1,000, 2,000, 3,000, 4,000, 5,000, 6,000, 7,000, 8,000, 9,000, 1 0,000, 20,000, 30,000, 40,000, 50,000, 60,000, 70,000, 80,000, 90,000, 100,000, 200,000, 300,000, 400,000, 500,000, 600,000, 700,000, 800,000, 900,000, 1,000,000, 2,000,000, 3,000,000, 4,00 0,000, 5,000,000, 6,000,000, 7,000,000, 8,000,000, 9,000,000, 10,000,000, 20,000,000, 30,000,000, 40,000,000, 50,000,000, 60,000,000, 70,000,000, 80,000,000, 90,000,000, 100,000,000, 2 In one embodiment, the present invention may be configured to perform any one or more of the operations described herein for 00,000,000, 300,000,000, 400,000,000, 500,000,000, 600,000,000, 700,000,000, 800,000,000, 900,000,000, 1,000,000,000 or more fragments.For example, the distributed computing system in Figure 19 can achieve up to approximately 1,000,000,000, 900,000,000, 80,000,000, 7000,000,000, 600,000,000, 500,000,000, 400,000,000, 300,000,000, 200,000,000, 100,000,000, 90,000,000, 80, 000,000, 70,000,000, 60,000,000, 50,000,000, 40,000,000, 30,000,000, 20,000,000, 10,000,000, 9,000,000, 8,000,000, 7,000,000, 6,000,000, 5,000,000, 4,000,000, 3,000,000, 2,000,000, 1, 000,000, 900,000, 800,000, 700,000, 600,000, 500,000, 400,000, 300,000, 200,000, 100,000, 90,000, 80,000, 70,000, 60,000, 50,000, 40,000, 30,000, 20,000, 10,000, 9,000, 8,000, 7,000, 6,000 , 5,000, 4,000, 3,000, 2,000, 1,000, 900, 800, 700, 600, 500, 400, 300, 200, 100, 90, 80, 70, 60, 50, 40, 30, 20, 10, 9, 8, 7, 6, 5, 4, 3, 2, or 1 fragment. For example, the distributed computing system of FIG. 19 may be configured to perform any one or more of the operations described herein for a number of fragments that are within a range defined by any two of the preceding values.

[0196] There are many techniques for distributing the problem between a client-side library and one or more endpoints. The distribution of functionality between the client-side library and one or more endpoints can vary based on the problem to be solved and / or the needs of the user.

[0197] In one approach, problem decomposition and problem dispatch may occur within a client-side library, see, e.g., FIG. 20. Problem decomposition and problem dispatch may occur on a user's digital computing device. Problem dispatch may distribute subsystems to electronic structure solvers that may be located at one or more remote endpoints. The remote endpoints may include nodes of one or more high performance computing systems, e.g., one or more non-classical computing systems, one or more classical computing systems, or a combination thereof.

[0198] In another approach, a problem may be sent from a client-side library to a remote endpoint, which may break down the problem and dispatch the problem, see, e.g., FIG. 21. In some cases, a problem may be dispatched from a remote endpoint to one or more second remote endpoints, which may include one or more nodes of a high-performance computing system. In some cases, one or more nodes of the one or more second remote endpoints may be local to the problem dispatch (e.g., on the same machine, in the same physical location on a separate machine, etc.). In some cases, one or more nodes of the one or more second remote endpoints may be remote to the problem dispatch.

[0199] In another approach, the problem may be broken down in a client-side library and sent to a remote endpoint that includes a problem dispatcher and then sent to one or more second remote endpoints that include electronic structure solvers.

[0200] FIG. 20 illustrates an example architecture 2000 of a distributed computing system with problem dispatch in a client-side library, according to some embodiments. In the illustrated embodiment, problem decomposition may be performed by client-side code. The problem decomposition may generate subsystems and / or solver parameters. The client-side system may forward the subsystems and / or solver parameters to the problem dispatch. The problem dispatch may create and execute a REST call, map the return values ​​to the input subsystems, and return the values ​​to the problem decomposition. The architecture of FIG. 20 may be advantageous when a user does not want to disclose the problem to be solved, since a fragment of data may be forwarded to the remote endpoint rather than a more complete implementation of the problem.

[0201] As shown in FIG. 20, the target 2010, the problem decomposer 2020, and the problem dispatcher 2040 may comprise portions of a client-side computing system 2050, such as a client-side library. The target 2010 may comprise a chemical system. The chemical system may comprise, for example, a molecule, a portion of a molecule, a fragment, an assembly, etc. The target 2010 may be decomposed into one or more fragments. For example, the target 2010 may be decomposed into fragments 2011, 2012, 2013, 2014, and 2015 in the problem decomposer. FIG. 20 illustrates the problem decomposer 2020. The problem decomposer may include instructions for performing any of the problem decomposition methods and / or techniques disclosed herein.

[0202] After decomposition, the fragments may be distributed by problem dispatch 2040 to one or more electronic structure solvers 2031, 2032, 2033, 2034, and 2035. Problem dispatch 2040 may create and / or implement the electronic structure solvers. Problem dispatch 2040 may pass parameters, such as input parameters from a user, to one or more solvers. Problem dispatch may return outputs from the solvers, such as, for example, energies.

[0203] Problem dispatch 2040 may be responsible for communicating with one or more solvers. In the illustrated embodiment, the one or more solvers may comprise portions of one or more remote endpoints 2060. The remote endpoints may be local to each other (e.g., on the same machine, the same physical location on separate machines, etc.) or remote to each other. The remote endpoints may include remote servers, cloud networks, portions of distributed computing systems, etc. The solvers invoked by problem dispatch may be specific to a problem type and / or fragment type. Problem dispatch may be responsible, in part, for distributing computational operations to reduce computation time and / or increase computational accuracy.

[0204] Fragments 2011, 2012, 2013, 2014, and 2015 may be encoded as fragments 2011', 2012', 2013', 2014', and 2015' on electronic structure solvers 2031, 2032, 2033, 2034, and 2035. Electronic structure solvers 2031, 2032, 2033, 2034, and 2035 may include one or more non-classical computing systems, one or more quantum computing systems, or one or more hybrid computing units, as described herein. Electronic structure solvers 2031, 2032, 2033, 2034, and 2035 calculate energies E for fragments 2011', 2012', 2013', 2014', and 2015', respectively. 1 , E 2 , E 3 , E 4 , and E 5 The electronic structure solvers 2031, 2032, 2033, 2034, and 2035 can calculate the energy E 1 , E 2 , E 3 , E 4 , and E 5 to problem dispatch 2040. The electronic structure solver may receive and return fragments from problem dispatch.

[0205] All or a portion of the multiple fragments may be presented to the electronic structure solver. The multiple fragments may be presented to the electronic structure solver in any order, serially or in parallel.

[0206] 20 shows five fragments and solutions, but the methods of the present disclosure may be used with any number of fragments and solvers. For example, the distributed computing system of FIG. 20 may have at least about 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 200, 300, 400, 500, 600, 700, 800, 900, 1,000, 2,000, 3,000, 4,000, 5,000, 6,000, 7,000, 8,000, 9,000, 1 0,000, 20,000, 30,000, 40,000, 50,000, 60,000, 70,000, 80,000, 90,000, 100,000, 200,000, 300,000, 400,000, 500,000, 600,000, 700,000, 800,000, 900,000, 1,000,000, 2,000,000, 3,000,000, 4,00 0,000, 5,000,000, 6,000,000, 7,000,000, 8,000,000, 9,000,000, 10,000,000, 20,000,000, 30,000,000, 40,000,000, 50,000,000, 60,000,000, 70,000,000, 80,000,000, 90,000,000, 100,000,000, 2 In one embodiment, the present invention may be configured to perform any one or more of the operations described herein for 00,000,000, 300,000,000, 400,000,000, 500,000,000, 600,000,000, 700,000,000, 800,000,000, 900,000,000, 1,000,000,000 or more fragments.For example, the distributed computing system in Figure 20 can achieve up to approximately 1,000,000,000, 900,000,000, 80,000,000, 7000,000,000, 600,000,000, 500,000,000, 400,000,000, 300,000,000, 200,000,000, 100,000,000, 90,000,000, 80, 000,000, 70,000,000, 60,000,000, 50,000,000, 40,000,000, 30,000,000, 20,000,000, 10,000,000, 9,000,000, 8,000,000, 7,000,000, 6,000,000, 5,000,000, 4,000,000, 3,000,000, 2,000,000, 1, 000,000, 900,000, 800,000, 700,000, 600,000, 500,000, 400,000, 300,000, 200,000, 100,000, 90,000, 80,000, 70,000, 60,000, 50,000, 40,000, 30,000, 20,000, 10,000, 9,000, 8,000, 7,000, 6,000 , 5,000, 4,000, 3,000, 2,000, 1,000, 900, 800, 700, 600, 500, 400, 300, 200, 100, 90, 80, 70, 60, 50, 40, 30, 20, 10, 9, 8, 7, 6, 5, 4, 3, 2, or 1 fragment. For example, the distributed computing system of FIG. 20 may be configured to perform any one or more of the operations described herein for a number of fragments that are within a range defined by any two of the preceding values.

[0207] FIG. 21 illustrates an example architecture 2100 of a distributed computing system with problem dispatch at a remote endpoint, according to some embodiments. In the illustrated embodiment, a problem may be sent from a client-side library to a first remote endpoint. Problem decomposition may be executed by code at the first remote endpoint. The problem decomposition may generate subsystems and / or solver parameters. The first remote endpoint may forward the subsystems and / or solver parameters to the problem dispatch. The problem dispatch may send and receive to one or more second remote endpoints. For example, the problem dispatch may create and execute REST calls to one or more second remote endpoints, map return values ​​from the one or more second remote endpoints to the input subsystems, and return values ​​to the problem decomposition.

[0208] In the illustrated embodiment, the client code may be a thin library containing versions of an electronic structure solver and problem decomposition that makes REST calls to a distributed computing system, e.g., a cloud infrastructure. For example, the REST calls may be made to a problem decomposition endpoint, which performs the problem decomposition and makes calls to a like number of electronic structure solver endpoints using problem dispatch. The problem decomposition remote calls may return to the user in the same or similar manner as local calls.

[0209] The architecture of FIG. 21 may be advantageous for several reasons. For example, restricted client-side libraries may minimize the amount of code located on the client machine. This may allow less powerful client-side machines to execute the code. This may also allow for easy distribution and updating of the software. In a second example, more computationally intensive problem decomposition methods may be implemented by increasing the amount of code on provider-owned hardware. The provider may implement GPU machines or high performance architectures to which users may not have direct access for problem decomposition.

[0210] As shown in Figure 21, a target 2110 may include a portion of a client-side computing system 2150, such as a client-side library. The target 2110 may include a chemical system. The chemical system may include, for example, molecules, portions of molecules, fragments, aggregates, etc. The architecture of Figure 21 may limit the computational load on the client-side computing system.

[0211] The client computing system may be in communication with a first one or more remote endpoints 2160. For example, the client computing system may send one or more of the target, the molecule, the conformation, one or more fragments, and the computation parameters to the first one or more remote endpoints. For example, the client computing system may receive one or more of values ​​corresponding to a solution to the problem (one or more energies, eigenvalues, structures, rates, etc.), information about the state of the computation, parameters regarding the progress of the computation, etc. The first one or more endpoints 2160 may include a problem decomposer 2120 and a problem dispatcher 2140. The target 2110 may be decomposed into one or more fragments in the problem decomposer 2120. For example, the target 2110 may be decomposed into fragments 2111, 2112, 2113, 2114, and 2115 in the problem decomposer. The problem decomposer may include instructions for performing any of the problem decomposition methods and / or techniques disclosed herein.

[0212] After decomposition, the fragments may be distributed by problem dispatch 2140 to one or more electronic structure solvers 2131, 2132, 2133, 2134, and 2135. Problem dispatch 2140 may create and / or implement the electronic structure solvers. Problem dispatch 2040 may pass parameters, such as input parameters from a user, to one or more solvers. Problem dispatch may return outputs from the solvers, such as, for example, energies.

[0213] Problem dispatch 2140 may be responsible for communicating with one or more solvers. In the illustrated embodiment, the one or more solvers may comprise portions of one or more second endpoints 2170. The second endpoints may be local to each other (e.g., on the same machine, the same physical location on separate machines, etc.) or remote to each other. The second endpoints may include remote servers, cloud networks, portions of distributed computing systems, etc. The solvers invoked by problem dispatch may be specific to a problem type and / or fragment type. Problem dispatch may be responsible, in part, for distributing computational operations to reduce computation time and / or increase computational accuracy.

[0214] Fragments 2111, 2112, 2113, 2114, and 2115 may be encoded as fragments 2111', 2112', 2113', 2114', and 2115' on electronic structure solvers 2131, 2132, 2133, 2134, and 2135. Electronic structure solvers 2131, 2132, 2133, 2134, and 2135 may include one or more non-classical computing systems, one or more quantum computing systems, or one or more hybrid computing units, as described herein. Electronic structure solvers 2131, 2132, 2133, 2134, and 2135 calculate energies E for fragments 2111', 2112', 2113', 2114', and 2115', respectively. 1 , E 2 , E 3 , E 4 , and E 5 The electronic structure solvers 2131, 2132, 2133, 2134, and 2135 may calculate the energy E 1 , E 2 , E 3 , E 4 , and E 5 to problem dispatch 2140. The electronic structure solver may receive and return fragments from problem dispatch.

[0215] All or a portion of the multiple fragments may be presented to the electronic structure solver. The multiple fragments may be presented to the electronic structure solver in any order, serially or in parallel.

[0216] 21 shows five fragments and solutions, but the methods of the present disclosure may be used with any number of fragments and solvers. For example, the distributed computing system of FIG. 21 may be configured to generate a distributed computing system that generates at least about 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 200, 300, 400, 500, 600, 700, 800, 900, 1,000, 2,000, 3,000, 4,000, 5,000, 6,000, 7,000, 8,000, 9,000, 1 0,000, 20,000, 30,000, 40,000, 50,000, 60,000, 70,000, 80,000, 90,000, 100,000, 200,000, 300,000, 400,000, 500,000, 600,000, 700,000, 800,000, 900,000, 1,000,000, 2,000,000, 3,000,000, 4,00 0,000, 5,000,000, 6,000,000, 7,000,000, 8,000,000, 9,000,000, 10,000,000, 20,000,000, 30,000,000, 40,000,000, 50,000,000, 60,000,000, 70,000,000, 80,000,000, 90,000,000, 100,000,000, 2 In one embodiment, the present invention may be configured to perform any one or more of the operations described herein for 00,000,000, 300,000,000, 400,000,000, 500,000,000, 600,000,000, 700,000,000, 800,000,000, 900,000,000, 1,000,000,000 or more fragments.For example, the distributed computing system in Figure 21 can achieve up to approximately 1,000,000,000, 900,000,000, 80,000,000, 7000,000,000, 600,000,000, 500,000,000, 400,000,000, 300,000,000, 200,000,000, 100,000,000, 90,000,000, 80, 000,000, 70,000,000, 60,000,000, 50,000,000, 40,000,000, 30,000,000, 20,000,000, 10,000,000, 9,000,000, 8,000,000, 7,000,000, 6,000,000, 5,000,000, 4,000,000, 3,000,000, 2,000,000, 1, 000,000, 900,000, 800,000, 700,000, 600,000, 500,000, 400,000, 300,000, 200,000, 100,000, 90,000, 80,000, 70,000, 60,000, 50,000, 40,000, 30,000, 20,000, 10,000, 9,000, 8,000, 7,000, 6,000 , 5,000, 4,000, 3,000, 2,000, 1,000, 900, 800, 700, 600, 500, 400, 300, 200, 100, 90, 80, 70, 60, 50, 40, 30, 20, 10, 9, 8, 7, 6, 5, 4, 3, 2, or 1 fragment. For example, the distributed computing system of FIG. 21 may be configured to perform any one or more of the operations described herein for a number of fragments that are within a range defined by any two of the preceding values.

[0217] Computer Systems The present disclosure provides computer systems programmed to implement the methods of the present disclosure. Figure 10 illustrates a computer system 1001 programmed or otherwise configured to determine a family of conformations of a chemical system, resolve at least one conformation in the family into a plurality of molecular fragments, determine using a hybrid computing unit a quantum mechanical energy and / or electronic structure for each of at least a subset of the plurality of molecular fragments, combine the determined quantum mechanical energies and / or electronic structures, and electronically output a report indicative of the combined quantum mechanical energies and / or electronic structures.

[0218] The computer system 1001 may define various aspects of the methods and systems of the present disclosure, such as determining a family of conformations of the chemical system, decomposing at least one conformation in the family into a plurality of molecular fragments, using a hybrid computing unit, determining the quantum mechanical energy and / or electronic structure of each of at least a subset of the plurality of molecular fragments, combining the determined quantum mechanical energies and / or electronic structures, and electronically outputting a report indicative of the combined quantum mechanical energies and / or electronic structures.

[0219] The computer system 1001 may be a user's electronic device or a computer system located remotely relative to the electronic device. The electronic device may be a mobile electronic device. The computer system 1001 includes a central processing unit (CPU, also referred to herein as "processor" and "computer processor") 1005, which may be a single-core or multi-core processor, or multiple processors for parallel processing. The computer system 1001 also includes memory or memory locations 1010 (e.g., random access memory, read-only memory, flash memory), electronic storage 1015 (e.g., hard disk), communication interface 1020 (e.g., network adapter) for communicating with one or more other systems, and peripheral devices 1025, such as cache, other memory, data storage, and / or electronic display adapters. The memory 1010, storage 1015, interface 1020, and peripheral devices 1025 are in communication with the CPU 1005 through a communication bus (solid lines), such as a motherboard. The storage device 1015 may be a data storage device (or data repository) for storing data. The computer system 1001 may be operatively coupled to a computer network ("network") 1030 with the aid of a communication interface 1020. The network 1030 may be the Internet, an Internet and / or an extranet, or an intranet and / or an extranet in communication with the Internet.

[0220] The network 1030 is, in some cases, a telecommunications and / or data network. The network 1030 may include one or more computer servers, which may enable distributed computing, such as cloud computing. For example, the one or more computers may enable cloud computing to perform various aspects of the analysis, calculation, and generation of the present disclosure, such as, for example, determining a set of conformations of a chemical system, decomposing at least one conformation in the set into a plurality of molecular fragments, using a hybrid computing unit to determine the quantum mechanical energy and / or electronic structure of each of at least a subset of the plurality of molecular fragments, combining the determined quantum mechanical energies and / or electronic structures, and electronically outputting a report indicating the combined quantum mechanical energy and / or electronic structure, through the network 1030 (the "cloud"). Such cloud computing may be provided, for example, by cloud computing platforms such as Amazon Web Services (AWS), Microsoft Azure, Google Cloud platform, and IBM Cloud. The network 1030 may, in some cases, implement a peer-to-peer network that may enable devices coupled to the computer system 1001 to act as clients or servers with the aid of the computer system 1001. "Cloud" services (including one or more of the cloud platforms mentioned above) may also be used to provide data storage.

[0221] The CPU 1005 can execute sequences of machine-readable instructions, which may be embodied in a program or software. The instructions may be stored in a memory location, such as the memory 1010. The instructions may be directed to the CPU 1005, which may then be programmed or otherwise configured to perform the methods of the present disclosure. Examples of operations performed by the CPU 1005 may include fetch, decode, execute, and writeback.

[0222] The CPU 1005 may be part of a circuit, such as an integrated circuit. One or more other components of the system 1001 may be included in the circuit. In some cases, the circuit is an application specific integrated circuit (ASIC). The CPU 1005 may include one or more general purpose processors, one or more graphics processing units (GPUs), or a combination thereof.

[0223] The storage device 1015 can store files such as drivers, libraries, and saved programs. The storage device 1015 can store (and sometimes exchange data with) user data, such as a family of conformations of a chemical system, a plurality of resolved molecular fragments, quantum mechanical energies and / or electronic structures of molecular fragments, combined quantum mechanical energies and / or electronic structures of conformers, a list of molecular fragments with their quantum mechanical energies and / or electronic structures, a list of conformers of a molecule with their combined quantum mechanical energies and / or electronic structures, and reports showing combined quantum mechanical energies and / or electronic structures. The computer system 1001 may, in some cases, include one or more additional data storage devices that are external to the computer system 1001, such as located on a remote server in communication with the computer system 1001 over an intranet or the Internet.

[0224] The computer system 1001 can communicate with one or more remote computer systems through the network 1030. For example, the computer system 1001 can communicate with a remote computer system of a user. Examples of remote computer systems include a personal computer (e.g., a portable PC), a slate or tablet PC (e.g., an Apple® iPad®, a Samsung® Galaxy Tab), a phone, a smartphone (e.g., an Apple® iPhone®, an Android-enabled device, a Blackberry®), or a personal digital assistant. A user can access the computer system 1001 through the network 1030. A user may control or define various aspects of the disclosed methods and systems, such as determining a family of conformations of a chemical system, decomposing at least one conformation in the family into a plurality of molecular fragments, using a hybrid computing unit to determine the quantum mechanical energy and / or electronic structure of each of at least a subset of the plurality of molecular fragments, combining the determined quantum mechanical energies and / or electronic structures, and electronically outputting a report indicative of the combined quantum mechanical energies and / or electronic structures.

[0225] Methods as described herein may be implemented as machine (e.g., computer processor) executable code stored in electronic storage locations of computer system 1001, such as memory 1010 or electronic storage 1015. Machine executable or machine readable code may be provided in the form of software. In use, the code may be executed by processor 1005. In some cases, the code may be retrieved from storage 1015 and stored in memory 1010 for ready access by processor 1005. In some situations, electronic storage 1015 may be eliminated and machine executable instructions stored in memory 1010.

[0226] The code may be precompiled and configured for use with a machine having a processor adapted to execute the code, or may be compiled during run-time. The code may be supplied in a programming language that may be selected to allow the code to execute in a precompiled or run-time compiled manner.

[0227] Aspects of the systems and methods provided herein, such as the computer system 1001, may be embodied in programming. Various aspects of the present technology may be considered as a "product" or "article of manufacture" that is typically in the form of machine (or processor) executable code and / or associated data executed or embodied on some type of machine-readable medium. The machine executable code may be stored in electronic storage such as memory (e.g., read-only memory, random access memory, flash memory, solid-state memory) or a hard disk. A "storage" type medium may include any or all of the tangible memory of a computer, processor, or the like, or their associated modules, such as various semiconductor memories, tape drives, disk drives, and the like, which may provide non-transitory storage at any time for software programming. All or portions of the software may at times be communicated over the Internet or various other telecommunications networks. Such communication may enable, for example, loading of the software from one computer or processor to another, for example, from an administrative server or host computer to an application server computer platform. Thus, another type of medium that may have a software element includes light waves, radio waves, and electromagnetic waves, such as those used across physical interfaces between local devices, through wired and optical land-line networks, and over various air links. Physical elements that convey such waves, such as wired or wireless links, optical links, or the like, may also be considered media that have software. As used herein, unless limited to non-transitory tangible "storage" media, terms such as computer or machine "readable medium" refer to any medium that participates in providing instructions to a processor for execution.

[0228] Thus, a machine-readable medium, such as a computer executable code, may take many forms, including but not limited to a tangible storage medium, a carrier wave medium, or a physical transmission medium. Non-volatile storage media include, for example, optical or magnetic disks, such as any of the storage devices in any computer, such as may be used to implement the databases, etc., shown in the figures. Volatile storage media include dynamic memory, such as the main memory of a computer platform. Tangible transmission media include coaxial cables, copper wire, and fiber optics, including the wires that comprise a bus within a computer system. Carrier wave transmission media may take the form of electric or electromagnetic signals, or acoustic or light waves, such as those generated during radio frequency (RF) and infrared (IR) data communications. Thus, common forms of computer readable media include, for example, a floppy disk, a flexible disk, a hard disk, magnetic tape, any other magnetic medium, a CD-ROM, a DVD or DVD-ROM, any other optical medium, punch cards paper tape, any other physical storage medium with a pattern of holes, RAM, ROM, PROM and EPROM, FLASH-EPROM, any other memory chip or cartridge, a carrier wave transporting data or instructions, a cable or link transporting such a carrier wave, or any other medium from which a computer may read programming code and / or data. Many of these forms of computer readable media may be involved in carrying one or more sequences of one or more instructions to a processor for execution.

[0229] The computer system 1001 may include or be in communication with an electronic display 1035 that includes a user interface (UI) 1040 for providing, for example, user selection of a group of conformations of a chemical system, conformations within the group for decomposition into a plurality of molecular fragments, at least a subset of the plurality of molecular fragments for determining quantum mechanical energy and / or electronic structure, and use of the Born-Oppenheimer approximation. Examples of UIs include, but are not limited to, graphic user interfaces (GUIs) and web-based user interfaces.

[0230] The computer system 1001 may include or be in communication with a non-classical computer (e.g., a quantum computer) 1045 for performing, for example, quantum algorithms (e.g., quantum mechanical energy and / or electronic structure calculations). The non-classical computer 1045 may be operatively coupled to the central processing unit 1005 and / or a network 1030 (e.g., the cloud).

[0231] The computer system of the present disclosure may be, for example, as described in International Application No. PCT / CA2017 / 050709, U.S. Application No. 15 / 486,960, U.S. Patent No. 9,537,953, and U.S. Patent No. 9,660,859, each of which is incorporated by reference in its entirety herein.

[0232] The methods and systems of the present disclosure may be implemented by one or more algorithms. The algorithm may be implemented by software when executed by the central processing unit 1005. The algorithm may, for example, determine a family of conformations of a chemical system, resolve at least one conformation in the family into a plurality of molecular fragments, determine, using a hybrid computing unit, a quantum mechanical energy and / or electronic structure for each of at least a subset of the plurality of molecular fragments, combine the determined quantum mechanical energies and / or electronic structures, and electronically output a report indicative of the combined quantum mechanical energies and / or electronic structures.

[0233] Although described herein with respect to particular systems, such as hybrid or quantum-classical computing or computing hardware, a problem (such as a quantum chemistry problem or simulation) may be solved using computing systems that comprise various types or combinations of systems, such as, for example, one or more classical computers, one or more non-classical computers (such as one or more quantum computers), or a combination of one or more classical computers and one or more non-classical computers. For example, Figure 15 illustrates an example of a system or combination of systems that may be used to solve a problem, such as a quantum chemistry problem or simulation.

[0234] Further details regarding systems and methods for performing quantum mechanical energy or electronic structure calculations on chemical systems can be found in U.S. Provisional Patent Application No. 62 / 593,060, filed November 30, 2017, and PCT Application No. PCT / CA2018 / 051531, filed November 30, 2018, which are hereby incorporated by reference in their entireties for all purposes. EXAMPLES

[0235] Example 1 (n-heptane) The correlation of the results of the total quantum mechanical energy calculations with and without PD was investigated for different conformations of the compound. Simulation results for fixed conformations with PD may not be within chemical accuracy. However, if this is due to a systematic error, comparing two erroneous results for different conformers of the same molecule can offset this error and provide an accurate relative quantum mechanical energy difference between two conformations of the molecule. Thus, this approach can be used to accurately pick the best conformer (e.g., the most stable conformer) based on their total quantum mechanical energy value, even when one does not have an optimally accurate estimate of the total quantum mechanical energy for each individual conformer. In this approach, a more aggressive PD technique (e.g., DC with a relatively small buffer size) can be used to find the best conformer from the group of available conformers. A more aggressive PD technique may result in smaller submolecules, which may mean that less quantum resources may be required to perform experiments for large molecules. Thus, this approach may enable highly efficient and accurate prediction of the most stable conformer of a chemical system using quantum computing resources.

[0236] In this example, n-heptane is targeted, as shown in Figure 4, where the dotted lines indicate the Bond Detached Atoms (BTAs) in the Fragment Molecular Orbital (FMO) fragmentation. A family of 40 conformations of n-heptane was generated by changing the four dihedral angles by 120 degrees (trans, gauche, gauche') and then symmetrically removing redundant and high energy conformations. To obtain a correlation of the total energy with and without problem decomposition (PD), CCSD was performed as a baseline reference, and two problem decomposition methods, DC-CCSD and FMO-CCSD, were applied to this molecular system. Seven fragments: two terminal CH 3 Group and 5 CH 2groups were considered. For DC, buffer sizes of 3 Å, 4 Å, 5 Å, and 6 Å were investigated. For FMO, two-body and three-body expansions were investigated. All calculations were performed using GAMESS-US with the 6-31G basis set. The GAMESS quantum chemistry package is described in Schmidt et al., "General Atomic and Molecular Electronic Structure System," Journal of Computational Chemistry, 1993, 14, 1347-1363, which is incorporated herein by reference in its entirety. The DC method was tested with buffer sizes smaller than 3 Å, but the calculations for nearly all conformers failed to converge to a solution.

[0237] Figure 5 illustrates a comparison between the canonical CCSD and DC-CCSD results, and between the canonical CCSD and FMO-CCSD results (list of conformer quantum mechanical energy values) for n-heptane. Good correlations are obtained between the results from the canonical CCSD and the CCSD with problem decomposition, with the coefficient of determination R 2 was greater than 0.96, except for FMO with a three-body expansion (FMO_3). Although DC-CCSD provided better results, DC calculations sometimes had difficulty converging on a solution. For n-heptane, FMO provided solutions for all 40 conformers investigated, while DC provided solutions for 35 and 36 conformers using 3 Å and 4 Å buffer sizes, respectively. It should also be noted that the number of spin orbitals required to solve one fragment can vary depending on the conformation in the case of DC calculations, since the buffer region is defined based on the distance from the center of the fragment.

[0238] Referring again to FIG. 5, several clusters of conformers were observed, for example as indicated by the dotted circles in the top right panel. We now briefly discuss why these clusters are observed. First, we investigated the relationship between the canonical CCSD energy and the diameter of the smallest sphere that can accommodate the conformers. This diameter can be considered as a measure of the conformer's structural compactness. As shown in the center panel of FIG. 6, the total quantum mechanical energy generally increases when the conformer becomes structurally compact due to steric repulsion. However, as can be seen, the diameter does not fully explain the conformer clustering in terms of the total quantum mechanical energy. Next, we investigated the relationship between the total quantum mechanical energy and the distance between the two outermost carbon atoms in the dihedral (1-4 distance). As illustrated by the right panel of FIG. 6, which illustrates the relationship between the total quantum mechanical energy and the smallest 1-4 distance for each conformer, the 1-4 distance explains the clustering behavior very well. The 1-4 distance changes depending on whether the dihedral angle is trans, gauche, or gauche'. In the trans case, the 1-4 distance is the longest. The gauche and gauche' dihedral angles have the same 1-4 distance, which is shorter than the trans 1-4 distance, causing a higher (less stable) total quantum mechanical energy due to steric repulsions. This is the main source of discretization of the total energy and the reason for the observed clustering of conformers with respect to the total quantum mechanical energy in this molecular system.

[0239] Example 2 (3-Methylheptane) As observed, both FMO and DC work relatively well for simple polymer systems. Next, diverse energy landscapes were generated for exploration by grafting one methyl group onto the carbon atom at the "3" position of n-heptane to obtain 3-methylheptane, as shown in Figure 7. The introduction of a methyl group at the "3" position puts the molecule into an asymmetric state. As in the case of n-heptane, a family of conformations was generated for 3-methylheptane by changing the four dihedral angles by 120 degrees (trans, gauche, gauche'), and 65 conformations were obtained after removing the high energy conformations.

[0240] FIG. 8 illustrates the quantum mechanical energy distributions (energy relative to lowest) obtained by CCSD to illustrate how one methyl group modulates and diversifies the quantum mechanical energy landscape from that of n-heptane. FIG. 9 illustrates a comparison between the canonical CCSD and DC-CCSD results, and between the canonical CCSD and FMO-CCSD results (list of quantum mechanical conformer energy values) for 3-methylheptane. As shown, the FMO (two-body) approach for 3-methylheptane has an R of 0.94. 2 The total energy obtained by DC with 3 Å buffer is somewhat closer to the canonical CCSD than that obtained by FMO2, and R 2 is lower than that of FMO. The DC approach for 3-methylheptane provides excellent agreement with the canonical CCSD when the buffer is increased to 4 Å. However, it is noted that DC suffers slightly from convergence failure here as well. FMO yielded solutions for all 65 conformers, while DC was able to yield solutions for 38 and 46 conformers with buffer sizes of 3 Å and 4 Å, respectively.

[0241] Example 3 (solver fragments in DMET) Example 3 is an example of DMET problem decomposition with a sequential implementation. Example 3 also shows an implementation in which fragmentation is specified by the user. The user can pass in a list of fragments, a molecule specification, and a mean field specification. The atoms of the molecule can be indexed in the order they return from a call to a quantum chemistry package, for example, a PySCF molar atom call. Each fragment can be specified by the indices of the atoms it contains, and optionally the type of solver the user wants to use to solve the fragment, and / or parameters that the solver can use.

[0242] If the user passes in a solver theory using the fragment atoms, the fragment may be solved with the specified solver, otherwise the fragment may be solved with the instance held by the problem decomposition object.

[0243] An example is shown below. pd = DMETProblemDecomposition() solver = FCISolver() pd.electronic_structure_solver = solver #Specification example: H 4 molecule H4_RING = """ H 0.7071067811865476 0.0 0.0 H 0.0 0.7071067811865476 0.0 H -1.0071067811865476 0.0 0.0 H 0.0 -1.0071067811865476 0.0 """ mol = gto.Mole() mol.atom = H4_RING mol.basis = "3-21g" mol.charge = 0 mol.spin = 0 mol.build() #Fragment specification, example #Example, fragment 1 contains the first two atoms and is solved by the VQE solver fragment1 = ([0,1], {"next_solver": "VQESolver", "solver_params" : {"hardware_backend_type" = "MicrosoftQSharpParametricSolver", "ansatz_type" : "MicrosoftQSharpParametricSolver.Ansatze.UCCSD"}}) #Fragment 2, for example, will be solved by the FCI solver if not specified fragment 2 = ([0,1], None) pd.simulate(mol, [fragment1, fragment2])

[0244] Example 4 (Solver parameters for fragments in DMET) In some cases, the methods for specifying fragments in DMET may not work for incremental methods. For example, in incremental methods, fragments may be generated automatically and they may be defined by the size of the interaction rather than by the atoms in the fragment. In this case, the methods and systems disclosed herein may adapt the fragment specification so that the first member of the tuple is the increment name instead of a list of atoms. Similar techniques may be applied to increasingly higher order perturbation theory methods.

[0245] for example, #1 The field terms can be solved with a VQE solver with parameters. fragment1 = ("1-body", {"next_solver": "VQESolver", "solver_params" : "hardware_backend_type" = "MicrosoftQSharpParametricSolver", "ansatz_type": "MicrosoftQSharpParametricSolver.Ansatze.UCCSD"}}) #The two-body terms can be solved with the default solver (e.g., FCI). fragment2 = ("2-body", None) pd.simulate(mol, [fragment1, fragment2])

[0246] Fragments that do not have a custom solver specified may use the default electronic structure solver held by the problem decomposition object, just as in the DMET case above.

[0247] Example 5 (FNO, nested solvers, and nested QEMIST) Example 5 shows an example where the problem decomposition may not be specified by the user. Example 5 also shows an example implementation of the Frozen Natural Orbital (FNO) method. The frozen natural orbital method may be combined with the coupled cluster (CC) method to increase the speed of the CC calculation. The FNO method may reduce the virtual space of the correlation calculation by at least about half. For example, the FNO method may reduce the computational cost by identifying and removing combinations of virtual orbitals that do not significantly contribute to the CC energy. FNO may be implemented together with the problem decomposition methods disclosed herein. For example, FNO may be implemented before the problem decomposition, between the problem decomposition and the electronic structure solver, or simply before the electronic structure solver. For example, a REST call sent to a distributed computing system, e.g., the cloud, may have a nested structure. In the nested structure, each request before the electronic structure solver may have a “next_solver” parameter that includes a call to the next step in the pipeline. Some examples are as follows:

[0248] FNO before electronic structure solver To use the FNO solver in conjunction with an electronic structure solver, a user may specify the electronic structure solver and its parameters in a call to the FNO solver, which may then make a call to the next step in the pipeline and return the results after the nested solvers have executed.

[0249] for example, next_solver_parameters = {"next_solver": "VQESolver", "solver_params" : {"hardware_backend_type" = "MicrosoftQSharpParametricSolver", "ansatz_type": "MicrosoftQSharpParametricSolver.Ansatze.UCCSD"}} fno = FNOSolver() fno.simulate(molecule, next_solver_parameters)

[0250] FNO before DMET before electronic structure solver To use DMET after FNO, the FNO solver may nest the DMET calls above. Note that the dictionary containing the calls to the two fragments may be inserted into the dictionary of the complete call.

[0251] for example, H4_RING = """ H 0.7071067811865476 0.0 0.0 H 0.0 0.7071067811865476 0.0 H -1.0071067811865476 0.0 0.0 H 0.0 -1.0071067811865476 0.0 """ mol = gto.Mole() mol.atom = H4_RING mol.basis = "3-21g" mol.charge = 0 mol.spin = 0 mol.build() #Generate fragments #The first one is solved with VQE and contains the first two atoms. fragment1 = ([0,1], {"next_solver": "VQESolver", "solver_params" : {"hardware_backend_type" = "MicrosoftQSharpParametricSolver", "ansatz_type": "MicrosoftQSharpParametricSolver.Ansatze.UCCSD"}}) fragment2 = ([2,3], {"next_solver": "FCISolver"}) full_pipeline = {"next_solver": "DMETProblemDecomposition", "solver_parameters": [fragment1, fragment2]}

[0252] DMET before FNO before electronic structure solver To use FNO for each fragment produced by DMET, one can add FNO as the next solver for each DMET fragment, and then add the next solver of the FNO call as the electronic structure solver, reversing the call above. In this example, FNO can be used for some fragments, e.g., fragments that may be more computationally expensive.

[0253] for example, H4_RING = """ H 0.7071067811865476 0.0 0.0 H 0.0 0.7071067811865476 0.0 H -1.0071067811865476 0.0 0.0 H 0.0 -1.0071067811865476 0.0 """ mol = gto.Mole() mol.atom = H4_RING mol.basis = "3-21g" mol.charge = 0 mol.spin = 0 mol.build() es_solver1 = {"next_solver": "VQESolver", "solver_params" : {"hardware_backend_type" ="MicrosoftQSharpParametricSolver", "ansatz_type": "MicrosoftQSharpParametricSolver.Ansatze.UCCSD"} fragment1 =([0,1],{"next_solver": "FNOSolver", "solver_parameters": es_solver1}) es_solver2 = {"next_solver": "FCISolver"} fragment2 =([2,3],{"next_solver": "FNOSolver", "solver_parameters": es_solver2}) pd = DMETProblemDecomposition() pd.simulate(mol, [fragment1, fragment2])

[0254] Although preferred embodiments of the present invention have been shown and described herein, it will be apparent to those skilled in the art that such embodiments are provided by way of example only. It is not intended that the present invention be limited by the specific examples provided herein. Although the present invention has been described with reference to the foregoing specification, the descriptions and illustrations herein are not intended to be construed in a limiting sense. Numerous variations, changes, and substitutions will now occur to those skilled in the art without departing from the present invention. Furthermore, it is to be understood that all aspects of the present invention are not limited to the specific depictions, configurations, or relative proportions set forth herein, which depend upon a variety of conditions and variables. It is to be understood that various alternatives to the embodiments of the invention described herein may be used in practicing the present invention. It is therefore intended that the present invention encompass any such alternatives, modifications, variations, or equivalents. It is intended that the following claims define the scope of the present invention, and that methods and structures within the scope of these claims and their equivalents are thereby covered.

Claims

1. 1. A method for performing quantum mechanical energy or electronic structure calculations on a chemical system, the method being performed by a hybrid computing unit comprising at least one classical computer and a distributed computing system comprising a plurality of electronic structure solvers, the electronic structure solvers comprising a non-classical computer, a quantum computing system, and a classical simulator, the method comprising: (a) decomposing, by a problem decomposition device included in the hybrid computing unit, at least one conformation in the family of conformations of the chemical system into a plurality of molecular fragments; (b) generating, by the problem decomposition device, for each of one or more of the plurality of molecular fragments, a respective data structure for the molecular fragment, the data structure specifying a type of the electronic structure solver to be used to determine the quantum mechanical energy or electronic structure of the molecular fragment and parameters to be passed to the electronic structure solver of the specified type and used by the electronic structure solver of the specified type to determine the quantum mechanical energy or electronic structure of the molecular fragment; (c) providing, by the problem decomposition device, one or more of the plurality of molecular fragments and their respective data structures to a problem dispatch included in the hybrid computing unit; (d) dispatching, by the problem dispatch included in the hybrid computing unit, one or more of the plurality of molecular fragments and their respective data structures to the electronic structure solver included in the distributed computing system of a type specified by each of the data structures, the problem dispatch including logic for creating and implementing the electronic structure solver using each of the data structures; (e) determining quantum mechanical energies or electronic structures of the plurality of molecular fragments according to the parameters included in each of the data structures using the electronic structure solver included in the distributed computing system of the type specified in each of the data structures; (f) combining, by the problem decomposer, the quantum mechanical energies or electronic structures determined in (e); and (g) electronically outputting a report indicative of the quantum mechanical energies or electronic structures combined in (f).

2. The method of claim 1 , wherein the non-classical computer comprises a Hitachi Ising solver or a coherent Ising machine based on optical parameters.

3. The method of claim 2 , wherein the classical simulator comprises a classical simulator of a quantum circuit.

4. The method of claim 1 , wherein the non-classical computers comprise different types of non-classical computers.

5. The method of claim 1 , wherein one of the quantum mechanical energies comprises nuclear-nuclear repulsion energy.

6. The method of claim 1 , further comprising providing an input to the hybrid computing unit, the input comprising a set of atomic coordinates of the chemical system.

7. 2. The method of claim 1, further comprising performing (a)-(g) for two or more conformations in said group of conformations of said chemical system.

8. 8. The method of claim 7, further comprising classifying the combined quantum mechanical energies or electronic structures of the plurality of molecular fragments.

9. 2. The method of claim 1 , wherein (a) comprises applying one or more members selected from the group consisting of a Fragment Molecular Orbital (FMO) method, a Divide-and-Conquer (DC) method, a Density Matrix Embedding Theory (DMET) method, a Density Matrix Renormalization Group (DMRG) method, a tensor network, and an incremental method.

10. (e) is determining a molecular fragment fermion Hamiltonian for the plurality of molecular fragments; Transforming the fermion Hamiltonian into an equivalent qubit Hamiltonian; converting the qubit Hamiltonian into a quantum circuit; and The method of claim 1 , comprising using the quantum circuit to determine the quantum mechanical energy or electronic structure of the molecular fragment.

11. The method of claim 10, further comprising determining the quantum mechanical energy or electronic structure using a molecular Hamiltonian.

12. The method of claim 10 , further comprising determining the quantum mechanical energy or electronic structure using an electronic Hamiltonian.

13. The method of claim 10 , wherein converting the fermion Hamiltonian to an equivalent qubit Hamiltonian comprises converting fermion operators of a Hamiltonian to qubit operators.

14. 10. The method of claim 1, further comprising performing an ab initio Molecular Dynamics (AIMD) simulation of the chemical system based on the quantum mechanical energy or electronic structure combined in (f).

15. Conducting the AIMD simulation includes: prior to (a), acquiring an indicia of a chemical system, the indicia including coordinates of each particle of a plurality of particles in the chemical system and a velocity of each particle in the chemical system; Following (f), (i) determining forces on each particle in the chemical system from the combined quantum mechanical energy or electronic structure and using a numerical gradient estimation technique, the numerical gradient estimation technique being implemented using the quantum computing system or the classical computer included in the distributed computing system; (ii) using the determined forces on particles in the chemical system to update the coordinates of each of the particles in the chemical system and the velocity of each of the particles in the chemical system; and (iii) electronically outputting a report indicating the updated coordinates or velocities; 15. The method of claim 14, wherein updating the coordinates and the velocity of each of the particles comprises updating the coordinates and the velocity of each of the particles by performing a numerical integration of forces on particles in the chemical system or a velocity velocimetry procedure.

16. The method of claim 15 , wherein the numerical gradient estimation technique includes Jordan's quantum algorithm.

17. The method of claim 15 , wherein performing the numerical integration includes performing a symplectic integration, a Runge-Kutta integration, and a Beaman integration.

18. The method of claim 1 , wherein the electronic structure solvers included in the distributed computing system of a type specified in a respective data structure include remote endpoints.

19. 20. The method of claim 18, wherein at least one of the remote endpoints comprises a non-classical computer.

20. 20. The method of claim 18, wherein the remote endpoint comprises a portion of a cloud computing system.

21. The method of claim 1 , further comprising, prior to (a), receiving the at least one conformation from a client-side library and dispatching the at least one conformation to a first remote endpoint.

22. 22. The method of claim 21, wherein at least one of (a)-(d) and (f) occurs at the first remote endpoint.

23. 23. The method of claim 22, wherein the electronic structure solvers included in the distributed computing system of the type specified in the respective data structures include one or more second remote endpoints.

24. The method of claim 23 , further comprising sending the report to the client-side library.

25. 24. The method of claim 23, wherein at least one of the second remote endpoints comprises a non-classical computer.

26. 24. The method of claim 23, wherein the one or more remote endpoints comprise portions of a cloud computing system.

27. 2. The method of claim 1 , wherein the decomposing in (a) is performed using the at least one classical computer.

28. 10. The method of claim 1 , wherein the determining in (e) is performed using at least one non-classical computer.

29. 10. The method of claim 1 , wherein the combining in (f) is performed using at least one classical computer.

30. 1. A system for performing quantum mechanical energy or electronic structure calculations on a chemical system, comprising: A hybrid computing unit operably coupled to a memory, the hybrid computing unit comprising a distributed computing system comprising at least one classical computer and a plurality of electronic structure solvers, the electronic structure solvers comprising a non-classical computer, a quantum computing system, and a classical simulator, the hybrid computing unit comprising at least: (a) decomposing, by a problem decomposition device included in the hybrid computing unit, at least one conformation in the family of conformations of the chemical system into a plurality of molecular fragments; (b) generating, by the problem decomposition device, for each of one or more of the plurality of molecular fragments, a respective data structure for the molecular fragment, the data structure specifying a type of the electronic structure solver to be used to determine the quantum mechanical energy or electronic structure of the molecular fragment and parameters to be passed to the electronic structure solver of the specified type and used by the electronic structure solver of the specified type to determine the quantum mechanical energy or electronic structure of the molecular fragment; (c) providing, by the problem decomposition device, one or more of the plurality of molecular fragments and their respective data structures to a problem dispatch included in the hybrid computing unit; (d) dispatching, by the problem dispatch included in the hybrid computing unit, one or more of the plurality of molecular fragments and their respective data structures to the electronic structure solver included in the distributed computing system of a type specified by each of the data structures, the problem dispatch including logic for creating and implementing the electronic structure solver using each of the data structures; (e) determining quantum mechanical energies or electronic structures of the plurality of molecular fragments according to the parameters included in each of the data structures using the electronic structure solver included in the distributed computing system of the type specified in each of the data structures; (f) combining, by the problem decomposer, the quantum mechanical energies or electronic structures determined in (e); and (g) electronically outputting a report indicative of the quantum mechanical energy or electronic structure combined in (f).

31. 31. The system of claim 30, further comprising a computer memory containing instructions for performing the quantum mechanical energy or electronic structure calculations on the chemical system, the hybrid computing unit being configured to execute the instructions to perform at least (a) through (g).

32. 1. A non-transitory computer readable medium comprising machine executable code that, when executed by a hybrid computing unit comprising at least one classical computer and a distributed computing system comprising a plurality of electronic structure solvers, a non-classical computer, a quantum computing system, and an electronic structure solver comprising a classical simulator, implements a method for performing quantum mechanical energy or electronic structure calculations on a chemical system, the method comprising: (a) decomposing, by a problem decomposition device included in the hybrid computing unit, at least one conformation in the family of conformations of the chemical system into a plurality of molecular fragments; (b) generating, by the problem decomposition device, for each of one or more of the plurality of molecular fragments, a respective data structure for the molecular fragment, the data structure specifying a type of the electronic structure solver to be used to determine the quantum mechanical energy or electronic structure of the molecular fragment and parameters to be passed to the electronic structure solver of the specified type and used by the electronic structure solver of the specified type to determine the quantum mechanical energy or electronic structure of the molecular fragment; (c) providing, by the problem decomposition device, one or more of the plurality of molecular fragments and their respective data structures to a problem dispatch included in the hybrid computing unit; (d) dispatching, by the problem dispatch included in the hybrid computing unit, one or more of the plurality of molecular fragments and their respective data structures to the electronic structure solver included in the distributed computing system of a type specified by each of the data structures, the problem dispatch including logic for creating and implementing the electronic structure solver using each of the data structures; (e) determining quantum mechanical energies or electronic structures of the plurality of molecular fragments according to the parameters included in each of the data structures using the electronic structure solver included in the distributed computing system of the type specified in each of the data structures; (f) combining, by the problem decomposer, the quantum mechanical energies or electronic structures determined in (e); and (g) electronically outputting a report indicative of the quantum mechanical energy or electronic structure combined in (f).

33. 1. A method for performing quantum mechanical energy or electronic structure calculations on a chemical system, the method being performed by a hybrid computing unit comprising a distributed computing system comprising a plurality of classical computers and at least one non-classical computer, the method comprising: (a) decomposing, by a problem decomposition device included in the distributed computing system, at least one conformation in the family of conformations of the chemical system into a plurality of molecular fragments; (b) generating, by the problem decomposition device, for each of one or more of the plurality of molecular fragments, a respective data structure for the molecular fragment, the data structure specifying a type of electronic structure solver to be used to determine the quantum mechanical energy or electronic structure of the molecular fragment and parameters to be passed to the electronic structure solver of the specified type and used by the electronic structure solver of the specified type to determine the quantum mechanical energy or electronic structure of the molecular fragment; (c) providing, by the problem decomposition device, one or more of the plurality of molecular fragments and their respective data structures to a problem dispatch included in the hybrid computing unit; (d) dispatching, by a problem dispatch included in the distributed computing system, one or more of the plurality of molecular fragments and their respective data structures to an electronic structure solver included in the distributed computing system of a type specified by each of the data structures, the problem dispatch including logic for creating and implementing the electronic structure solver using each of the data structures; (e) determining quantum mechanical energies or electronic structures of the plurality of molecular fragments according to the parameters included in each of the data structures using the electronic structure solver included in the distributed computing system of the type specified in each of the data structures; (f) combining, by the problem decomposer, the quantum mechanical energies or electronic structures determined in (e); and (g) electronically outputting a report indicative of the quantum mechanical energies or electronic structures combined in (f).

34. 34. The method of claim 33, wherein the at least one non-classical computer comprises at least one quantum computer.

35. 35. The method of claim 34, wherein the at least one quantum computer comprises one or more members selected from the group consisting of a quantum hardware device and a classical simulator of a quantum circuit.

36. 34. The method of claim 33, wherein the at least one non-classical computer comprises a plurality of different types of non-classical computers.

37. 34. The method of claim 33, wherein one of the quantum mechanical energies comprises nuclear-nuclear repulsion energy.

38. 34. The method of claim 33, further comprising providing an input to the hybrid computing unit, the input comprising a set of atomic coordinates of the chemical system.

39. 34. The method of claim 33, further comprising performing (a)-(f) for two or more conformations in said group of conformations of said chemical system.

40. 40. The method of claim 39, further comprising classifying the combined quantum mechanical energy or electronic structure of the plurality of molecular fragments.

41. 34. The method of claim 33, wherein (a) comprises applying one or more members selected from the group consisting of a fragment molecular orbital (FMO) method, a divide and conquer (DC) method, a density matrix embedding theory (DMET) method, a density matrix renormalization group (DMRG) method, a tensor network, and an incremental method.

42. (e) is (a) determining a molecular fragment fermion Hamiltonian for the plurality of molecular fragments; (b) transforming the fermion Hamiltonian into an equivalent qubit Hamiltonian; (c) converting the qubit Hamiltonian into a quantum circuit; and 34. The method of claim 33, comprising: (d) using the quantum circuit to determine the quantum mechanical energy or electronic structure of the molecular fragment.

43. 43. The method of claim 42, further comprising determining the quantum mechanical energy or electronic structure using a molecular Hamiltonian.

44. 43. The method of claim 42, further comprising determining the quantum mechanical energy or electronic structure using an electronic Hamiltonian.

45. 43. The method of claim 42, wherein converting the fermion Hamiltonian to an equivalent qubit Hamiltonian comprises converting fermion operators of the Hamiltonian to qubit operators.

46. 34. The method of claim 33, further comprising performing an ab initio molecular dynamics (AIMD) simulation of the chemical system based on the quantum mechanical energy or electronic structure combined in (f).

47. Conducting the AIMD simulation includes: prior to (a), acquiring an indicia of a chemical system, the indicia including coordinates of each particle of a plurality of particles in the chemical system and a velocity of each particle in the chemical system; (f) followed by (i) determining forces on each particle in the chemical system from the combined quantum mechanical energy or electronic structure and using a numerical gradient estimation technique, the numerical gradient estimation technique being implemented using a quantum computing system or the classical computer included in the distributed computing system; (ii) using the determined forces on particles in the chemical system to update the coordinates of each of the particles in the chemical system and the velocity of each of the particles in the chemical system; and (iii) electronically outputting a report indicating the updated coordinates or velocities; 47. The method of claim 46, wherein updating the coordinates and the velocity of each of the particles comprises updating the coordinates and the velocity of each of the particles by performing a numerical integration of forces on particles in the chemical system or a velocity vector procedure.

48. 48. The method of claim 47, wherein the numerical gradient estimation technique includes Jordan's quantum algorithm.

49. 48. The method of claim 47, wherein performing the numerical integration includes performing symplectic integration, Runge-Kutta integration, and Beaman integration.

50. 34. The method of claim 33, wherein the electronic structure solver includes a remote end.

51. 51. The method of claim 50, wherein at least one of the remote endpoints comprises a non-classical computer.

52. 51. The method of claim 50, wherein the remote endpoint comprises a portion of a cloud computing system.

53. 34. The method of claim 33, further comprising, prior to (a), receiving the at least one conformation from a client-side library and dispatching the at least one conformation to a first remote endpoint.

54. 54. The method of claim 53, wherein at least one of (a)-(d) and (f) occurs at the first remote endpoint.

55. 55. The method of claim 54, wherein the electronic structure solver includes one or more second remote endpoints.

56. 56. The method of claim 55, further comprising sending the report to the client-side library.

57. 56. The method of claim 55, wherein at least one of the second remote endpoints comprises a non-classical computer.

58. 56. The method of claim 55, wherein the one or more remote endpoints comprise portions of a cloud computing system.

59. 34. The method of claim 33, wherein the decomposing in (a) is performed using at least one classical computer of the plurality of classical computers.

60. 34. The method of claim 33, wherein the determining in (e) is performed using the at least one non-classical computer.

61. 34. The method of claim 33, wherein the combining in (f) is performed using at least one classical computer of the plurality of classical computers.

62. 1. A system for performing quantum mechanical energy or electronic structure calculations on a chemical system, comprising: A hybrid computing unit operably coupled to a memory, the hybrid computing unit comprising a distributed computing system comprising a plurality of classical computers and at least one electronic structure solver, the electronic structure solver comprising a non-classical computer, a quantum computing system, and a classical simulator, the hybrid computing unit comprising at least (a) decomposing, by a problem decomposition device included in a hybrid computing unit, at least one conformation in the family of conformations of the chemical system into a plurality of molecular fragments; (b) generating, by the problem decomposition device, for each of one or more of the plurality of molecular fragments, a respective data structure for the molecular fragment, the data structure specifying a type of the electronic structure solver to be used to determine the quantum mechanical energy or electronic structure of the molecular fragment and parameters to be passed to the electronic structure solver of the specified type and used by the electronic structure solver of the specified type to determine the quantum mechanical energy or electronic structure of the molecular fragment; (c) providing, by the problem decomposition device, one or more of the plurality of molecular fragments and their respective data structures to a problem dispatch included in the hybrid computing unit; (d) dispatching, by the problem dispatch included in the hybrid computing unit, one or more of the plurality of molecular fragments and their respective data structures to the electronic structure solver included in the distributed computing system of a type specified by each of the data structures, the problem dispatch including logic for creating and implementing the electronic structure solver using each of the data structures; (e) determining quantum mechanical energies or electronic structures of the plurality of molecular fragments according to the parameters included in each of the data structures using the electronic structure solver included in the distributed computing system of the type specified in each of the data structures; (f) combining, by the problem decomposer, the quantum mechanical energies or electronic structures determined in (e); and (g) electronically outputting a report indicative of the quantum mechanical energy or electronic structure combined in (f).

63. 63. The system of claim 62, further comprising a computer memory containing instructions for performing the quantum mechanical energy or electronic structure calculations on the chemical system, wherein the hybrid computing unit is configured to execute the instructions to perform at least (a) through (g).

64. 1. A non-transitory computer readable medium comprising machine executable code that, when executed by a hybrid computing unit comprising a plurality of classical computers and at least one electronic structure solver, a distributed computing system comprising an electronic structure solver comprising non-classical computers, a quantum computing system, and a classical simulator, implements a method for performing quantum mechanical energy or electronic structure calculations on a chemical system, the method comprising: (a) decomposing, by a problem decomposition device included in the hybrid computing unit, at least one conformation in the family of conformations of the chemical system into a plurality of molecular fragments; (b) generating, by the problem decomposition device, for each of one or more of the plurality of molecular fragments, a respective data structure for the molecular fragment, the data structure specifying a type of the electronic structure solver to be used to determine the quantum mechanical energy or electronic structure of the molecular fragment and parameters to be passed to the electronic structure solver of the specified type and used by the electronic structure solver of the specified type to determine the quantum mechanical energy or electronic structure of the molecular fragment; (c) providing, by the problem decomposition device, one or more of the plurality of molecular fragments and their respective data structures to a problem dispatch included in the hybrid computing unit; (d) dispatching, by the problem dispatch included in the hybrid computing unit, one or more of the plurality of molecular fragments and their respective data structures to the electronic structure solver included in the distributed computing system of a type specified by each of the data structures, the problem dispatch including logic for creating and implementing the electronic structure solver using each of the data structures; (e) determining quantum mechanical energies or electronic structures of the plurality of molecular fragments according to the parameters included in each of the data structures using the electronic structure solver included in the distributed computing system of the type specified in each of the data structures; (f) combining, by the problem decomposer, the quantum mechanical energies or electronic structures determined in (e); and (g) electronically outputting a report indicative of the quantum mechanical energy or electronic structure combined in (f).

65. The method of claim 1 , further comprising dispatching one or more of the plurality of molecular fragments to one or more of the electronic structure solvers.

66. 34. The method of claim 33, further comprising dispatching one or more of the plurality of molecular fragments to one or more of the electronic structure solvers.

Images