Electric Field Gradient Calculations for Ice VIII and IX using Polarizable Embedding: A Comparative Study on Classical Computers and Quantum Simulators

Dániel Nagy,Peter Reinholdt,Phillip W. K. Jensen,Erik Rosendahl Kjellgren,Karl Michael Ziems,Aaron Fitzpatrick,Stefan Knecht,Jacob Kongsted,Sonia Coriani,Stephan P. A. Sauer
2024-04-23
Abstract:We test the performance of the Polarizable Embedding Variational Quantum Eigensolver Self-Consistent-Field (PE-VQE-SCF) model for computing electric field gradients with comparisons to conventional complete active space self-consistent-field (CASSCF) calculations and experimental results. We compute quadrupole coupling constants for ice VIII and ice IX. We find that the inclusion of the environment is crucial for obtaining results that match the experimental data. The calculations for ice VIII are within the experimental uncertainty for both CASSCF and VQE-SCF for oxygen and lie close to the experimental value for ice IX as well. With the VQE-SCF, which is based on an Adaptive Derivative-Assembled Problem-Tailored (ADAPT) ansatz, we find that the inclusion of the environment and the size of the different basis sets do not directly affect the gate counts. However, by including an explicit environment, the wavefunction and, therefore, the optimization problem becomes more complicated, which usually results in the need to include more operators from the operator pool, thereby increasing the depth of the circuit.
Quantum Physics,Chemical Physics
What problem does this paper attempt to address?
This paper discusses the use of the Polarization Embedding Variational Quantum Eigensolver Self-Consistent Field (PE-VQE-SCF) model to calculate the electric field gradients of ice VIII and ice IX, comparing them with the traditional Complete Active Space Self-Consistent Field (CASSCF) calculations and experimental results. The study found that including the environment is crucial for matching experimental data. For ice VIII, both CASSCF and VQE-SCF calculations are within the experimental error range, while for ice IX, they are close to the experimental values. The VQE-SCF is based on the Adaptive VQE with Specific Problem Hypothesis (ADAPT-VQE), and although the inclusion of the environment does not directly affect gate count, it makes the wave function and optimization problem more complex, possibly requiring more operators and increasing circuit depth. The paper also discusses the theoretical background of electric field gradients, polarization embedding model, and VQE-SCF method, and provides computational details. Quantum simulators were used to calculate ice VIII and ice IX to validate the performance of the PE-VQE-SCF method in determining nuclear quadrupole coupling constants.