Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2404.16149

$Δ$ADAPT-VQE: Toward Accurate Calculation of Excitation Energies on Quantum Computers for BODIPY Molecules With Application in Photodynamic Therapy

Anton Nykänen, Leander Thiessen, Elsi-Mari Borrelli, Vijay B. Krishna, Stefan Knecht, Fabijan Pavošević

Published Apr 24, 2024• 1 citation
arXivHTMLPDF

Abstract

Quantum chemistry simulations offer a cost-effective way for computational design of BODIPY photosensitizers with potential use in photodynamic therapy (PDT). However, accurate predictions of photophysical properties, such as excitation energies, pose a challenge for the popular time-dependent density functional theory (TDDFT) and equation-of-motion coupled cluster with singles and doubles (EOM-CCSD) methods. By contrast, reliable descriptions can be achieved by multi-reference quantum chemistry methods, though unfortunately, their computational cost grows exponentially with the number of correlated electrons. Alternatively, quantum computing holds a great potential for exact simulation of photophysical properties in a computationally more efficient way. To this end, we introduce the state-specific $Δ$UCCSD-VQE (unitary coupled cluster with singles and doubles variational quantum eigensolver) and $Δ$ADAPT-VQE methods in which the electronically excited state is calculated via a non-Aufbau electronic configuration. The accuracy and capability of the developed methods are assessed against experimentally determined excitation energies for six BODIPY-derivatives. We show that the proposed methods predict accurate vertical excitation energies that are not only in good agreement with experimental reference data but also outperform popular quantum chemistry methods, such as TDDFT and EOM-CCSD. Spurred by its impressive performance and simplicity, we are confident that $Δ$ADAPT will emerge as the method of choice for guiding the rational design of photosensitizers for PDT and photocatalysis in the era of near-term quantum computing.

DOI
10.48550/arxiv.2404.16149

Related papers

Toward Accurate Calculation of Excitation Energies on Quantum Computers with ΔADAPT-VQE: A Case Study of BODIPY Derivatives

Anton Nykänen, Leander Thiessen, Elsi-Mari Borrelli +3 more

Quantum chemistry simulations offer a cost-effective way to computationally design BODIPY photosensitizers. However, accurate predictions of excitation energies pose a challenge for time-dependent density functional theory and equation-of-motion coupled-cluster singles and doubles methods. By contrast, reliable predictions can be achieved by multireference quantum chemistry methods; unfortunately, their computational cost increases exponentially with the number of electrons. Alternatively, quantum computing holds potential for an exact simulation of the photophysical properties in a computationally more efficient way. Herein, we introduce the state-specific ΔUCCSD-VQE (unitary coupled-cluster singles and doubles-variational quantum eigensolver) and ΔADAPT-VQE methods in which the electronically excited state is calculated via a non-Aufbau configuration. We show for six BODIPY derivatives that the proposed methods predict accurate excitation energies that are in good agreement with those from experiments. Due to its performance and simplicity, we believe that ΔADAPT will become a useful approach for the simulation of BODIPY photosensitizers on near-term quantum devices.

The Journal of Physical Chemistry Letters•Jul 2, 2024•9 citations

Quantum Algorithms for Photoreactivity in Cancer-Targeted Photosensitizers

Yanbing Zhou, Pablo A. M. Casares, Diksha Dhawan +5 more

Photodynamic therapy (PDT) is a targeted cancer treatment that uses light-activated photosensitizers to generate reactive oxygen species that selectively destroy tumor cells, generally causing less collateral damage than conventional treatments. However, its clinical success hinges on the availability of photosensitizers with strong optical sensitivity and high efficiency in generating reactive oxygen species. While classical computational methods have provided useful insights into photosensitizer design, they struggle to scale and often lack the accuracy needed for these simulations. In this work, we show how fault-tolerant quantum algorithms can be used to identify promising photosensitizer candidates for PDT. To predict photosensitizer performance, we assess two computational properties. First, we quantify light sensitivity by calculating the cumulative absorption in the therapeutic window with a threshold projection algorithm. Second, we determine the efficiency of reactive oxygen generation by estimating intersystem crossing (ISC) rates using the evolution-proxy approach, complemented by a vibronic dynamic treatment where appropriate. We apply these algorithms to a clinically relevant and actively pursued class of photosensitizers, BODIPY derivatives, including heavy-atom and transition-metal-substituted systems that are challenging for classical methods. Our resource estimates, obtained with PennyLane, suggest that systems with active spaces ranging from 11 to 45 spatial orbitals can be simulated using $180$-$350$ logical qubits and Toffoli gate depths between $10^7$ and $10^9$, placing our algorithms within reach of realistic fault-tolerant quantum devices. This paves the way to an efficient quantum-based workflow for designing photosensitizers that can accelerate the discovery of new PDT agents.

quant-pharXiv:2512.15889•Dec 19, 2025•0 citations

Molecular-excited-state calculations with the qubit-excitation-based adaptive variational quantum eigensolver protocol

Yordan S. Yordanov, C. H. W. Barnes, David R. M. Arvidsson-Shukur

Calculations of molecular spectral properties, like photodissociation rates and absorption bands, rely on knowledge of the excited state energies of the molecule of interest. Protocols based on the variational quantum eigensolver (VQE) are promising candidates to calculate such energies on emerging noisy intermediate-scale quantum (NISQ) computers. The successful implementation of these protocols on NISQ computers, relies on ans\"atze that can accurately approximate the molecular states and that can be implemented by shallow quantum circuits. We introduce the excited qubit-excitation-based adaptive (e-QEB-ADAPT)-VQE protocol to calculate molecular-excited-state energies. The e-QEB-ADAPT-VQE constructs efficient problem-tailored ans\"atze by iteratively appending evolutions of qubit excitation operators. The e-QEB-ADAPT-VQE also improves on previous ADAPT-VQE protocol in that it is designed to be independent on the choice of initial reference state. We perform classical numerical simulations for LiH and ${\mathrm{BeH}}_{2}$ to benchmark the performance of the e-QEB-ADAPT-VQE. We demonstrate that the e-QEB-ADAPT-VQE can construct highly accurate ans\"atze that require at least an order of magnitude fewer cnot gates than standard fixed unitary coupled-cluster ans\"atze, such as the UCCSD and the GUCCSD. We also show that the e-QEB-ADAPT-VQE is more successful in constructing ans\"atze for excited molecular states than other ADAPT-VQE protocols.

Physical Review A•Sep 29, 2022•17 citations

Quantum Simulation of Molecular Electronic States with a Transcorrelated Hamiltonian: Higher Accuracy with Fewer Qubits

Ashutosh Kumar, Ayush Asthana, Conner Masteran +5 more

Simulation of electronic structure is one of the most promising applications on noisy intermediate-scale quantum (NISQ) era devices. However, NISQ devices suffer from a number of challenges like limited qubit connectivity, short coherence times, and sizable gate error rates. Thus, desired quantum algorithms should require shallow circuit depths and low qubit counts to take advantage of these devices. Here, we attempt to reduce quantum resource requirements for molecular simulations on a quantum computer while maintaining the desired accuracy with the help of classical quantum chemical theories of canonical transformation and explicit correlation. In this work, compact ab initio Hamiltonians are generated classically, in the second quantized form, through an approximate similarity transformation of the Hamiltonian with (a) an explicitly correlated two-body unitary operator with generalized pair excitations that remove the Coulombic electron-electron singularities from the Hamiltonian and (b) a unitary one-body operator to efficiently capture the orbital relaxation effects required for accurate description of the excited states. The resulting transcorrelated Hamiltonians are able to describe both the ground and the excited states of molecular systems in a balanced manner. Using the variational quantum eigensolver (VQE) method based on the unitary coupled cluster with singles and doubles (UCCSD) ansatz and only a minimal basis set (ANO-RCC-MB), we demonstrate that the transcorrelated Hamiltonians can produce ground state energies comparable to the reference CCSD energies with the much larger cc-pVTZ basis set. This leads to a reduction in the number of required CNOT gates by more than 3 orders of magnitude for the chemical species studied in this work. Furthermore, using the quantum equation of motion (qEOM) formalism in conjunction with the transcorrelated Hamiltonian, we are able to reduce the deviations in the excitation energies from the reference EOM-CCSD/cc-pVTZ values by an order of magnitude. The transcorrelated Hamiltonians developed here are Hermitian and contain only one- and two-body interaction terms and thus can be easily combined with any quantum algorithm for accurate electronic structure simulations.

Journal of Chemical Theory and Computation•Aug 19, 2022•44 citations

Molecular excited state calculations with the QEB-ADAPT-VQE

Yordan S. Yordanov, C. H. W. Barnes, David R. M. Arvidsson-Shukur

Calculations of molecular spectral properties, like photodissociation rates and absorption bands, rely on knowledge of the excited state energies of the molecule of interest. Protocols based on the variational quantum eigensolver (VQE) are promising candidates to calculate such energies on emerging noisy intermediate scale quantum (NISQ) computers. The successful implementation of these protocols on NISQ computers, relies on ansätze that can accurately approximate the molecular states and that can be implemented by shallow quantum circuits. In this paper, we introduce the excited qubit-excitation-based adaptive (e-QEB-ADAPT)-VQE protocol to calculate molecular excited state energies. The e-QEB-ADAPT-VQE constructs efficient problem-tailored ansätze by iteratively appending evolutions of qubit excitation operators. The e-QEB-ADAPT-VQE is an adaptation of the QEB-ADAPT-VQE protocol, which is designed to be independent on the choice of an initial reference state. We perform classical numerical simulations for LiH and BeH$_2$ to benchmark the performance of the e-QEB-ADAPT-VQE. We demonstrate that the e-QEB-ADAPT-VQE can construct highly accurate ansätze that require at least an order of magnitude fewer $CNOT$s than standard fixed UCC ansätze, such as the UCCSD and the GUCCSD.

arXiv:2106.06296•Jun 11, 2021•1 citation

Spin-Flip Unitary Coupled Cluster Method: Toward Accurate Description of Strong Electron Correlation on Quantum Computers

Fabijan Pavošević, Ivano Tavernelli, Ángel Rubio

Quantum computers have emerged as a promising platform to simulate strong electron correlation that is crucial to catalysis and photochemistry. However, owing to the choice of a trial wave function employed in the variational quantum eigensolver (VQE) algorithm, accurate simulation is restricted to certain classes of correlated phenomena. Herein, we combine the spin-flip (SF) formalism with the unitary coupled cluster with singles and doubles (UCCSD) method via the quantum equation-of-motion (qEOM) approach to allow for an efficient simulation of a large family of strongly correlated problems. We show that the developed qEOM-SF-UCCSD/VQE method outperforms its UCCSD/VQE counterpart for simulation of the cis – trans isomerization of ethylene, and the automerization of cyclobutadiene and the predicted qEOM-SF-UCCSD/VQE barrier heights are in a good agreement with the experimentally determined values. The developments presented herein will further stimulate the investigation of this approach for simulations of other types of correlated/entangled phenomena on quantum computers.

arXiv:2307.07092•The Journal of Physical Chemistry Letters•Aug 28, 2023•18 citations
ΔΔΔADAPT-VQE: Toward Accurate Calculation of Excitation Energies on Quantum Computers for BODIPY Molecules With Application in Photodynamic Therapy
Quantum chemistry simulations offer a cost-effective way for computational design of BODIPY photosensitizers with potential use in photodynamic therapy (PDT). However, accurate predictions of photophysical properties, such as excitation energies, pose a challenge for the popular time-dependent density functional theory (TDDFT) and equation-of-motion coupled cluster with singles and doubles (EOM-CCSD) methods. By contrast, reliable descriptions can be achieved by multi-reference quantum chemistry methods, though unfortunately, their computational cost grows exponentially with the number of correlated electrons. Alternatively, quantum computing holds a great potential for exact simulation of photophysical properties in a computationally more efficient way. To this end, we introduce the state-specific ΔΔΔUCCSD-VQE (unitary coupled cluster with singles and doubles variational quantum eigensolver) and ΔΔΔADAPT-VQE methods in which the electronically excited state is calculated via a non-Aufbau electronic configuration. The accuracy and capability of the developed methods are assessed against experimentally determined excitation energies for six BODIPY-derivatives. We show that the proposed methods predict accurate vertical excitation energies that are not only in good agreement with experimental reference data but also outperform popular quantum chemistry methods, such as TDDFT and EOM-CCSD. Spurred by its impressive performance and simplicity, we are confident that ΔΔΔADAPT will emerge as the method of choice for guiding the rational design of photosensitizers for PDT and photocatalysis in the era of near-term quantum computing.
Quantum chemistry simulations offer a cost-effective way to computationally design BODIPY photosensitizers. However, accurate predictions of excitation energies pose a challenge for time-dependent density functional theory and equation-of-motion coupled-cluster singles and doubles methods. By contrast, reliable predictions can be achieved by multireference quantum chemistry methods; unfortunately, their computational cost increases exponentially with the number of electrons. Alternatively, quantum computing holds potential for an exact simulation of the photophysical properties in a computationally more efficient way. Herein, we introduce the state-specific ΔUCCSD-VQE (unitary coupled-cluster singles and doubles-variational quantum eigensolver) and ΔADAPT-VQE methods in which the electronically excited state is calculated via a non-Aufbau configuration. We show for six BODIPY derivatives that the proposed methods predict accurate excitation energies that are in good agreement with those from experiments. Due to its performance and simplicity, we believe that ΔADAPT will become a useful approach for the simulation of BODIPY photosensitizers on near-term quantum devices.
Photodynamic therapy (PDT) is a targeted cancer treatment that uses light-activated photosensitizers to generate reactive oxygen species that selectively destroy tumor cells, generally causing less collateral damage than conventional treatments. However, its clinical success hinges on the availability of photosensitizers with strong optical sensitivity and high efficiency in generating reactive oxygen species. While classical computational methods have provided useful insights into photosensitizer design, they struggle to scale and often lack the accuracy needed for these simulations. In this work, we show how fault-tolerant quantum algorithms can be used to identify promising photosensitizer candidates for PDT. To predict photosensitizer performance, we assess two computational properties. First, we quantify light sensitivity by calculating the cumulative absorption in the therapeutic window with a threshold projection algorithm. Second, we determine the efficiency of reactive oxygen generation by estimating intersystem crossing (ISC) rates using the evolution-proxy approach, complemented by a vibronic dynamic treatment where appropriate. We apply these algorithms to a clinically relevant and actively pursued class of photosensitizers, BODIPY derivatives, including heavy-atom and transition-metal-substituted systems that are challenging for classical methods. Our resource estimates, obtained with PennyLane, suggest that systems with active spaces ranging from 11 to 45 spatial orbitals can be simulated using 180180180-350350350 logical qubits and Toffoli gate depths between 10710^7107 and 10910^9109, placing our algorithms within reach of realistic fault-tolerant quantum devices. This paves the way to an efficient quantum-based workflow for designing photosensitizers that can accelerate the discovery of new PDT agents.
Calculations of molecular spectral properties, like photodissociation rates and absorption bands, rely on knowledge of the excited state energies of the molecule of interest. Protocols based on the variational quantum eigensolver (VQE) are promising candidates to calculate such energies on emerging noisy intermediate-scale quantum (NISQ) computers. The successful implementation of these protocols on NISQ computers, relies on ans\"atze that can accurately approximate the molecular states and that can be implemented by shallow quantum circuits. We introduce the excited qubit-excitation-based adaptive (e-QEB-ADAPT)-VQE protocol to calculate molecular-excited-state energies. The e-QEB-ADAPT-VQE constructs efficient problem-tailored ans\"atze by iteratively appending evolutions of qubit excitation operators. The e-QEB-ADAPT-VQE also improves on previous ADAPT-VQE protocol in that it is designed to be independent on the choice of initial reference state. We perform classical numerical simulations for LiH and BeH2{\mathrm{BeH}}_{2}BeH2​ to benchmark the performance of the e-QEB-ADAPT-VQE. We demonstrate that the e-QEB-ADAPT-VQE can construct highly accurate ans\"atze that require at least an order of magnitude fewer cnot gates than standard fixed unitary coupled-cluster ans\"atze, such as the UCCSD and the GUCCSD. We also show that the e-QEB-ADAPT-VQE is more successful in constructing ans\"atze for excited molecular states than other ADAPT-VQE protocols.
Simulation of electronic structure is one of the most promising applications on noisy intermediate-scale quantum (NISQ) era devices. However, NISQ devices suffer from a number of challenges like limited qubit connectivity, short coherence times, and sizable gate error rates. Thus, desired quantum algorithms should require shallow circuit depths and low qubit counts to take advantage of these devices. Here, we attempt to reduce quantum resource requirements for molecular simulations on a quantum computer while maintaining the desired accuracy with the help of classical quantum chemical theories of canonical transformation and explicit correlation. In this work, compact ab initio Hamiltonians are generated classically, in the second quantized form, through an approximate similarity transformation of the Hamiltonian with (a) an explicitly correlated two-body unitary operator with generalized pair excitations that remove the Coulombic electron-electron singularities from the Hamiltonian and (b) a unitary one-body operator to efficiently capture the orbital relaxation effects required for accurate description of the excited states. The resulting transcorrelated Hamiltonians are able to describe both the ground and the excited states of molecular systems in a balanced manner. Using the variational quantum eigensolver (VQE) method based on the unitary coupled cluster with singles and doubles (UCCSD) ansatz and only a minimal basis set (ANO-RCC-MB), we demonstrate that the transcorrelated Hamiltonians can produce ground state energies comparable to the reference CCSD energies with the much larger cc-pVTZ basis set. This leads to a reduction in the number of required CNOT gates by more than 3 orders of magnitude for the chemical species studied in this work. Furthermore, using the quantum equation of motion (qEOM) formalism in conjunction with the transcorrelated Hamiltonian, we are able to reduce the deviations in the excitation energies from the reference EOM-CCSD/cc-pVTZ values by an order of magnitude. The transcorrelated Hamiltonians developed here are Hermitian and contain only one- and two-body interaction terms and thus can be easily combined with any quantum algorithm for accurate electronic structure simulations.
Calculations of molecular spectral properties, like photodissociation rates and absorption bands, rely on knowledge of the excited state energies of the molecule of interest. Protocols based on the variational quantum eigensolver (VQE) are promising candidates to calculate such energies on emerging noisy intermediate scale quantum (NISQ) computers. The successful implementation of these protocols on NISQ computers, relies on ansätze that can accurately approximate the molecular states and that can be implemented by shallow quantum circuits. In this paper, we introduce the excited qubit-excitation-based adaptive (e-QEB-ADAPT)-VQE protocol to calculate molecular excited state energies. The e-QEB-ADAPT-VQE constructs efficient problem-tailored ansätze by iteratively appending evolutions of qubit excitation operators. The e-QEB-ADAPT-VQE is an adaptation of the QEB-ADAPT-VQE protocol, which is designed to be independent on the choice of an initial reference state. We perform classical numerical simulations for LiH and BeH2_22​ to benchmark the performance of the e-QEB-ADAPT-VQE. We demonstrate that the e-QEB-ADAPT-VQE can construct highly accurate ansätze that require at least an order of magnitude fewer CNOTCNOTCNOTs than standard fixed UCC ansätze, such as the UCCSD and the GUCCSD.
Quantum computers have emerged as a promising platform to simulate strong electron correlation that is crucial to catalysis and photochemistry. However, owing to the choice of a trial wave function employed in the variational quantum eigensolver (VQE) algorithm, accurate simulation is restricted to certain classes of correlated phenomena. Herein, we combine the spin-flip (SF) formalism with the unitary coupled cluster with singles and doubles (UCCSD) method via the quantum equation-of-motion (qEOM) approach to allow for an efficient simulation of a large family of strongly correlated problems. We show that the developed qEOM-SF-UCCSD/VQE method outperforms its UCCSD/VQE counterpart for simulation of the cis – trans isomerization of ethylene, and the automerization of cyclobutadiene and the predicted qEOM-SF-UCCSD/VQE barrier heights are in a good agreement with the experimentally determined values. The developments presented herein will further stimulate the investigation of this approach for simulations of other types of correlated/entangled phenomena on quantum computers.