1
|
Dobrautz W, Sokolov IO, Liao K, Ríos PL, Rahm M, Alavi A, Tavernelli I. Toward Real Chemical Accuracy on Current Quantum Hardware Through the Transcorrelated Method. J Chem Theory Comput 2024; 20:4146-4160. [PMID: 38723159 PMCID: PMC11137825 DOI: 10.1021/acs.jctc.4c00070] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2024] [Revised: 04/15/2024] [Accepted: 04/16/2024] [Indexed: 05/29/2024]
Abstract
Quantum computing is emerging as a new computational paradigm with the potential to transform several research fields including quantum chemistry. However, current hardware limitations (including limited coherence times, gate infidelities, and connectivity) hamper the implementation of most quantum algorithms and call for more noise-resilient solutions. We propose an explicitly correlated Ansatz based on the transcorrelated (TC) approach to target these major roadblocks directly. This method transfers, without any approximation, correlations from the wave function directly into the Hamiltonian, thus reducing the resources needed to achieve accurate results with noisy quantum devices. We show that the TC approach allows for shallower circuits and improves the convergence toward the complete basis set limit, providing energies within chemical accuracy to experiment with smaller basis sets and, thus, fewer qubits. We demonstrate our method by computing bond lengths, dissociation energies, and vibrational frequencies close to experimental results for the hydrogen dimer and lithium hydride using two and four qubits, respectively. To demonstrate our approach's current and near-term potential, we perform hardware experiments, where our results confirm that the TC method paves the way toward accurate quantum chemistry calculations already on today's quantum hardware.
Collapse
Affiliation(s)
- Werner Dobrautz
- Department
of Chemistry and Chemical Engineering, Chalmers
University of Technology, 41296 Gothenburg, Sweden
| | - Igor O. Sokolov
- IBM
Quantum, IBM Research Zurich, Säumerstrasse 4, 8803 Rüschlikon, Switzerland
| | - Ke Liao
- Max
Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Pablo López Ríos
- Max
Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Martin Rahm
- Department
of Chemistry and Chemical Engineering, Chalmers
University of Technology, 41296 Gothenburg, Sweden
| | - Ali Alavi
- Max
Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield
Road, Cambridge CB2 1EW, U.K.
| | - Ivano Tavernelli
- IBM
Quantum, IBM Research Zurich, Säumerstrasse 4, 8803 Rüschlikon, Switzerland
| |
Collapse
|
2
|
Giovannini T, Scavino M, Koch H. Time-Dependent Multilevel Density Functional Theory. J Chem Theory Comput 2024; 20:3601-3612. [PMID: 38648031 DOI: 10.1021/acs.jctc.4c00041] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2024]
Abstract
We present a novel three-layer approach based on multilevel density functional theory (MLDFT) and polarizable molecular mechanics to simulate the electronic excitations of chemical systems embedded in an external environment within the time-dependent DFT formalism. In our method, the electronic structure of a target system, the chromophore, is determined in the field of an embedded inactive layer, which is treated as frozen. Long-range interactions are described by employing the polarizable fluctuating charge (FQ) force field. The resulting MLDFT/FQ thus accurately describes both electrostatics (and polarization) and non-electrostatic target-environment interactions. The robustness and reliability of the approach are demonstrated by comparing our results with experimental data reported for various organic molecules in solution.
Collapse
Affiliation(s)
| | - Marco Scavino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Henrik Koch
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| |
Collapse
|
3
|
Ma Y, Li Z, Chen X, Ding B, Li N, Lu T, Zhang B, Suo B, Jin Z. Machine-learning assisted scheduling optimization and its application in quantum chemical calculations. J Comput Chem 2023; 44:1174-1188. [PMID: 36648254 DOI: 10.1002/jcc.27075] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2022] [Revised: 12/16/2022] [Accepted: 12/23/2022] [Indexed: 01/18/2023]
Abstract
Easy and effective usage of computational resources is crucial for scientific calculations, both from the perspectives of timeliness and economic efficiency. This work proposes a bi-level optimization framework to optimize the computational sequences. Machine-learning (ML) assisted static load-balancing, and different dynamic load-balancing algorithms can be integrated. Consequently, the computational and scheduling engine of the ParaEngine is developed to invoke optimized quantum chemical (QC) calculations. Illustrated benchmark calculations include high-throughput drug suit, solvent model, P38 protein, and SARS-CoV-2 systems. The results show that the usage rate of given computational resources for high throughput and large-scale fragmentation QC calculations can primarily profit, and faster accomplishing computational tasks can be expected when employing high-performance computing (HPC) clusters.
Collapse
Affiliation(s)
- Yingjin Ma
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| | - ZhiYing Li
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| | - Xin Chen
- ShenZhen Bay Laboratory, Shenzhen, China
| | - Bowen Ding
- Institute of Chemistry, Chinese Academy of Sciences, Beijing, China
| | - Ning Li
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
- College of Chemistry and Materials Engineering, Wenzhou University, Wen Zhou, China
| | - Teng Lu
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| | - Baohua Zhang
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| | - BingBing Suo
- Department of Physics, Northwest University, Xi'an, China
| | - Zhong Jin
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| |
Collapse
|
4
|
Giovannini T, Marrazzini G, Scavino M, Koch H, Cappelli C. Integrated Multiscale Multilevel Approach to Open Shell Molecular Systems. J Chem Theory Comput 2023; 19:1446-1456. [PMID: 36780359 PMCID: PMC10018740 DOI: 10.1021/acs.jctc.2c00805] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/14/2023]
Abstract
We present a novel multiscale approach to study the electronic structure of open shell molecular systems embedded in an external environment. The method is based on the coupling of multilevel Hartree-Fock (MLHF) and Density Functional Theory (MLDFT), suitably extended to the unrestricted formalism, to Molecular Mechanics (MM) force fields (FF). Within the ML region, the system is divided into active and inactive parts, thus describing the most relevant interactions (electrostatic, polarization, and Pauli repulsion) at the quantum level. The surrounding MM part, which is formulated in terms of nonpolarizable or polarizable FFs, permits a physically consistent treatment of long-range electrostatics and polarization effects. The approach is extended to the calculation of hyperfine coupling constants and applied to selected nitroxyl radicals in an aqueous solution.
Collapse
Affiliation(s)
| | - Gioia Marrazzini
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Marco Scavino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Henrik Koch
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy.,Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
5
|
Abstract
We introduce the concept of fragment localized molecular orbitals (FLMOs), which are Hartree-Fock molecular orbitals localized in specific fragments constituting a molecular system. In physical terms, we minimize the local electronic energies of the different fragments, at the cost of maximizing the repulsion between them. To showcase the approach, we rationalize the main interactions occurring in large biological systems in terms of interactions between the fragments of the system. In particular, we study an anticancer drug intercalated within DNA and retinal in anabaena sensory rhodopsin as prototypes of molecular systems embedded in biological matrixes. Finally, the FLMOs are exploited to rationalize the formation of two oligomers, prototypes of amyloid diseases, such as Parkinson and Alzheimer.
Collapse
Affiliation(s)
| | - Henrik Koch
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy.,Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| |
Collapse
|
6
|
Goletto L, Kjønstad EF, Folkestad SD, Høyvik IM, Koch H. Linear-Scaling Implementation of Multilevel Hartree-Fock Theory. J Chem Theory Comput 2021; 17:7416-7427. [PMID: 34747179 PMCID: PMC8675138 DOI: 10.1021/acs.jctc.1c00299] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We introduce a new algorithm for the construction of the two-electron contributions to the Fock matrix in multilevel Hartree-Fock (MLHF) theory. In MLHF, the density of an active molecular region is optimized, while the density of an inactive region is fixed. The MLHF equations are solved in a reduced molecular orbital (MO) basis localized to the active region. The locality of the MOs can be exploited to reduce the computational cost of the Fock matrix: the cost related to the inactive density becomes linear scaling, while the iterative cost related to the active density is independent of the system size. We demonstrate the performance of this new algorithm on a variety of systems, including amino acid chains, water clusters, and solvated systems.
Collapse
Affiliation(s)
- Linda Goletto
- Department of Chemistry, Norwegian University of Science and Technology, Trondheim 7491, Norway
| | - Eirik F Kjønstad
- Department of Chemistry, Norwegian University of Science and Technology, Trondheim 7491, Norway
| | - Sarai D Folkestad
- Department of Chemistry, Norwegian University of Science and Technology, Trondheim 7491, Norway
| | - Ida-Marie Høyvik
- Department of Chemistry, Norwegian University of Science and Technology, Trondheim 7491, Norway
| | - Henrik Koch
- Department of Chemistry, Norwegian University of Science and Technology, Trondheim 7491, Norway.,Scuola Normale Superiore, Pisa 56126, Italy
| |
Collapse
|