Quantum computing in pharma: A multilayer embedding approach for near future applications

Multiscale Embedding for Quantum Computing

We present a novel multiscale embedding scheme that links conventional QM/MM embedding and bootstrap embedding (BE) to allow simulations of large chemical systems on limited quantum devices. We also propose a mixed-basis BE scheme that facilitates BE calculations on extended systems using classical computers with limited memory resources. Benchmark data suggest the combination of these two strategies as a robust path in attaining the correlation energies of large realistic systems, combining the proven accuracy of BE with chemical and biological systems of interest in a lower computational cost method. Due to the flexible tunability of the resource requirements and systematic fragment construction, future developments in the realization of quantum computers naturally offer improved accuracy for multiscale BE calculations.

Calculating the Energy Profile of an Enzymatic Reaction on a Quantum Computer

Quantum computing (QC) provides a promising avenue for enabling quantum chemistry calculations, which are classically impossible due to computational complexity that increases exponentially with system size. As fully fault-tolerant algorithms and hardware, for which an exponential speedup is predicted, are currently out of reach, recent research efforts have been dedicated to developing and scaling algorithms for Noisy Intermediate-Scale Quantum (NISQ) devices to showcase the practical usefulness of such machines. To demonstrate the usefulness of NISQ devices in the field of chemistry, we apply our recently developed FAST-VQE algorithm and a state-of-the-art quantum gate reduction strategy based on propositional satisfiability together with standard optimization tools for the simulation of the rate-determining proton transfer step for CO2 hydration catalyzed by carbonic anhydrase resulting in the first application of a quantum computing device for the simulation of an enzymatic reaction. To this end, we have combined classical force field simulations with quantum mechanical methods on classical and quantum computers in a hybrid calculation approach. The presented technique significantly enhances the accuracy and capabilities of QC-based molecular modeling and finally pushes it into compelling and realistic applications. The framework is general and can be applied beyond the case of computational enzymology.

Contextual Subspace Auxiliary-Field Quantum Monte Carlo: Improved Bias with Reduced Quantum Resources

Using trial wave functions prepared on quantum devices to reduce the bias of auxiliary-field quantum Monte Carlo (QC-AFQMC) has established itself as a promising hybrid approach to the simulation of strongly correlated many body systems. Here, we further reduce the required quantum resources by decomposing the trial wave function into classical and quantum parts, respectively treated by classical and quantum devices, within the contextual subspace projection formalism. Importantly, we show that our algorithm is compatible with the recently developed matchgate shadow protocol for efficient overlap calculation in QC-AFQMC. Investigating the nitrogen dimer and the reductive decomposition of ethylene carbonate in lithium-based batteries, we observe that our method outperforms a number of established algorithm for ground state energy computations, while reaching chemical precision with less than half of the original number of qubits.

A general framework for active space embedding methods with applications in quantum computing

We developed a general framework for hybrid quantum-classical computing of molecular and periodic embedding approaches based on an orbital space separation of the fragment and environment degrees of freedom. We demonstrate its potential by presenting a specific implementation of periodic range-separated DFT coupled to a quantum circuit ansatz, whereby the variational quantum eigensolver and the quantum equation-of-motion algorithm are used to obtain the low-lying spectrum of the embedded fragment Hamiltonian. The application of this scheme to study localized electronic states in materials is showcased through the accurate prediction of the optical properties of the neutral oxygen vacancy in magnesium oxide (MgO). Despite some discrepancies in the position of the main absorption band, the method demonstrates competitive performance compared to state-of-the-art ab initio approaches, particularly evidenced by the excellent agreement with the experimental photoluminescence emission peak.

Quantum-centric high performance computing for quantum chemistry

High performance computing (HPC) is renowned for its capacity to tackle complex problems. Meanwhile, quantum computing (QC) provides a potential way to accurately and efficiently solve quantum chemistry problems. The emerging field of quantum-centric high performance computing (QCHPC), which merges these two powerful technologies, is anticipated to enhance computational capabilities for solving challenging problems in quantum chemistry. The implementation of QCHPC for quantum chemistry requires interdisciplinary research and collaboration across multiple fields, including quantum chemistry, quantum physics, computer science and so on. This perspective provides an introduction to the quantum algorithms that are suitable for deployment in QCHPC, focusing on conceptual insights rather than technical details. Parallel strategies to implement these algorithms on quantum-centric supercomputers are discussed. We also summarize high performance quantum emulating simulators, which are considered a viable tool to explore QCHPC. We conclude with challenges and outlooks in this field.

A Vision for the Future of Multiscale Modeling

The rise of modern computer science enabled physical chemistry to make enormous progresses in understanding and harnessing natural and artificial phenomena. Nevertheless, despite the advances achieved over past decades, computational resources are still insufficient to thoroughly simulate extended systems from first principles. Indeed, countless biological, catalytic and photophysical processes require ab initio treatments to be properly described, but the breadth of length and time scales involved makes it practically unfeasible. A way to address these issues is to couple theories and algorithms working at different scales by dividing the system into domains treated at different levels of approximation, ranging from quantum mechanics to classical molecular dynamics, even including continuum electrodynamics. This approach is known as multiscale modeling and its use over the past 60 years has led to remarkable results. Considering the rapid advances in theory, algorithm design, and computing power, we believe multiscale modeling will massively grow into a dominant research methodology in the forthcoming years. Hereby we describe the main approaches developed within its realm, highlighting their achievements and current drawbacks, eventually proposing a plausible direction for future developments considering also the emergence of new computational techniques such as machine learning and quantum computing. We then discuss how advanced multiscale modeling methods could be exploited to address critical scientific challenges, focusing on the simulation of complex light-harvesting processes, such as natural photosynthesis. While doing so, we suggest a cutting-edge computational paradigm consisting in performing simultaneous multiscale calculations on a system allowing the various domains, treated with appropriate accuracy, to move and extend while they properly interact with each other. Although this vision is very ambitious, we believe the quick development of computer science will lead to both massive improvements and widespread use of these techniques, resulting in enormous progresses in physical chemistry and, eventually, in our society.

Accelerating reliable multiscale quantum refinement of protein-drug systems enabled by machine learning

Biomacromolecule structures are essential for drug development and biocatalysis. Quantum refinement (QR) methods, which employ reliable quantum mechanics (QM) methods in crystallographic refinement, showed promise in improving the structural quality or even correcting the structure of biomacromolecules. However, vast computational costs and complex quantum mechanics/molecular mechanics (QM/MM) setups limit QR applications. Here we incorporate robust machine learning potentials (MLPs) in multiscale ONIOM(QM:MM) schemes to describe the core parts (e.g., drugs/inhibitors), replacing the expensive QM method. Additionally, two levels of MLPs are combined for the first time to overcome MLP limitations. Our unique MLPs+ONIOM-based QR methods achieve QM-level accuracy with significantly higher efficiency. Furthermore, our refinements provide computational evidence for the existence of bonded and nonbonded forms of the Food and Drug Administration (FDA)-approved drug nirmatrelvir in one SARS-CoV-2 main protease structure. This study highlights that powerful MLPs accelerate QRs for reliable protein-drug complexes, promote broader QR applications and provide more atomistic insights into drug development.

Solvent Distribution Effects on Quantum Chemical Calculations with Quantum Computers

We present a combination of three-dimensional reference interaction site model self-consistent field (3D-RISM-SCF) theory and the variational quantum eigensolver (VQE) to consider the solvent distribution effects within the framework of quantum-classical hybrid computing. The present method, 3D-RISM-VQE, does not include any statistical errors from the solvent configuration sampling owing to the analytical treatment of the statistical solvent distribution. We apply 3D-RISM-VQE to compute the spatial distribution functions of solvent water around a water molecule, the potential and Helmholtz energy curves of NaCl, and to analyze the Helmholtz energy component and related properties of H2O and NH4+. Moreover, we utilize 3D-RISM-VQE to analyze the extent to which solvent effects alter the efficiency of quantum calculations compared with calculations in the gas phase using the L1-norms of molecular electronic Hamiltonians. Our results demonstrate that the efficiency of quantum chemical calculations on a quantum computer in solution is virtually the same as that in the gas phase.

A theoretical spectroscopy study of the photoluminescence properties of narrow band Eu2+-doped phosphors containing multiple candidate doping centers. Prediction of an unprecedented narrow band red phosphor

We have previously presented a computational protocol that is based on an embedded cluster model and operates in the framework of TD-DFT in conjunction with the excited state dynamics (ESD) approach. The protocol is able to predict the experimental absorption and emission spectral shapes of Eu2+-doped phosphors. In this work, the applicability domain of the above protocol is expanded to Eu2+-doped phosphors bearing multiple candidate Eu doping centers. It will be demonstrated that this protocol provides full control of the parameter space that describes the emission process. The stability of Eu doping at various centers is explored through local energy decomposition (LED) analysis of DLPNO-CCSD(T) energies. This enables further development of the understanding of the electronic structure of the targeted phosphors, the diverse interactions between Eu and the local environment, and their impact on Eu doping probability, and control of the emission properties. Hence, it can be employed to systematically improve deficiencies of existing phosphor materials, defined by the presence of various intensity emission bands at undesired frequencies, towards classes of candidate Eu2+-doped phosphors with desired narrow band red emission. For this purpose, the chosen study set consists of three UCr4C4-based narrow-band phosphors, namely the known alkali lithosilicates RbNa[Li3SiO4]2:Eu2+ (RNLSO2), RbNa3[Li3SiO4]4:Eu2+ (RNLSO) and their isotypic nitridolithoaluminate phosphors consisting of CaBa[LiAl3N4]2:Eu2+ (CBLA2) and the proposed Ca3Ba[LiAl3N4]4:Eu2+ (CBLA), respectively. The theoretical analysis presented in this work led us to propose a modification of the CBLA2 phosphor that should have improved and unprecedented narrow band red emission properties. Finally, we believe that the analysis presented here is important for the future rational design of novel Eu2+-doped phosphor materials, with a wide range of applications in science and technology.

Perspective on the Current State-of-the-Art of Quantum Computing for Drug Discovery Applications

Computational chemistry is an essential tool in the pharmaceutical industry. Quantum computing is a fast evolving technology that promises to completely shift the computational capabilities in many areas of chemical research by bringing into reach currently impossible calculations. This perspective illustrates the near-future applicability of quantum computation of molecules to pharmaceutical problems. We briefly summarize and compare the scaling properties of state-of-the-art quantum algorithms and provide novel estimates of the quantum computational cost of simulating progressively larger embedding regions of a pharmaceutically relevant covalent protein-drug complex involving the drug Ibrutinib. Carrying out these calculations requires an error-corrected quantum architecture that we describe. Our estimates showcase that recent developments on quantum phase estimation algorithms have dramatically reduced the quantum resources needed to run fully quantum calculations in active spaces of around 50 orbitals and electrons, from estimated over 1000 years using the Trotterization approach to just a few days with sparse qubitization, painting a picture of fast and exciting progress in this nascent field.