We present calculations of both the ground-and excited-state energies of spin defects in solids carried out on a quantum computer, using a hybrid classical-quantum protocol.We focus on the negatively charged nitrogen-vacancy center in diamond and on the double vacancy in 4H SiC, which are of interest for the realization of quantum technologies.We employ a recently developed first-principles quantum embedding theory to describe point defects embedded in a periodic crystal and to derive an effective Hamiltonian, which is then transformed to a qubit Hamiltonian by means of a parity transformation.We use the variational quantum eigensolver (VQE) and quantum subspace expansion methods to obtain the ground and excited states of spin qubits, respectively, and we propose a promising strategy for noise mitigation.We show that by combining zero-noise extrapolation techniques and constraints on electron occupation to overcome the unphysical-state problem of the VQE algorithm, one can obtain reasonably accurate results on near-term-noisy architectures for ground-and excited-state properties of spin defects.
Simulating electronic structure on a quantum computer requires encoding of fermionic systems onto qubits. Common encoding methods transform a fermionic system of
N spin-orbitals into an
N-qubit system, but many of the fermionic configurations do not respect the required conditions and symmetries of the system so the qubit Hilbert space in this case may have unphysical states and thus cannot be fully utilized. We propose a generalized qubit-efficient encoding (QEE) scheme that requires the qubit number to be only logarithmic in the number of configurations that satisfy the required conditions and symmetries. For the case of considering only the particle-conserving and singlet configurations, we reduce the qubit count to an upper bound of
O(mlog2N), where
m is the number of particles. This QEE scheme is demonstrated on an
H2 molecule in the 6-31G basis set and a
LiH molecule in the STO-3G basis set using fewer qubits than the common encoding methods. We calculate the ground-state energy surfaces using a variational quantum eigensolver algorithm with a hardware-efficient ansatz circuit. We choose to use a hardware-efficient ansatz since most of the Hilbert space in our scheme is spanned by desired configurations so a heuristic search for an eigenstate is sensible. The simulations are performed on IBM Quantum machines and the Qiskit simulator with a noise model implemented from a IBM Quantum machine. Using the methods of measurement error mitigation and error-free linear extrapolation, we demonstrate that most of the distributions of the extrapolated energies using our QEE scheme agree with the exact results obtained by Hamiltonian diagonalization in the given basis sets within chemical accuracy. Our proposed scheme and results show the feasibility of quantum simulations for larger molecular systems in the noisy intermediate-scale quantum (NISQ) era. The number of terms in the Hamiltonian has an upper bound of $\mathcal{O}(\frac{{N}^{2m+1}}{(m\ensuremath{-}1)!\phantom{\rule{0.16em}{0ex}}m!})$ for the QEE scheme while it scales as
O(N4) for the Jordan-Wigner encoding scheme. Nevertheless, we present several cases where QEE is useful.
One of the first practical applications of quantum computers is expected to be molecular modelling. Performing this task would profoundly affect areas such as chemistry, materials science and drug synthesis. Modelling of molecules, which are classically intractable, can be achieved with just over
30 qubits, whereas state of the art quantum computers already have more than
50 qubits. The Variational Quantum Eigensolver (VQE) algorithm and VQE based protocols, are promising candidates to enable this task on emerging Noisy Intermediate-Scale Quantum (NISQ) computers. These protocols require short quantum circuits and short coherence times, and are particularly resilient to quantum errors. Nevertheless, there is still a significant gap between the accuracy and the coherence times of current NISQ computers, and the hardware requirements of VQE protocols to simulate practically interesting molecules. In this thesis, I present my contribution to narrowing this gap by developing VQE protocols for molecular modelling that are less demanding on quantum hardware. The VQE relies on the Rayleigh-Ritz variational principle to estimate the eigenvalues of a Hamiltonian operator, by minimizing its expectation value with respect to a trial quantum state, prepared by an ansatz. A major challenge for the practical realisation of VQE protocols on NISQ computers is to construct an ansatz that: (1) can accurately approximate the eigenstates of the Hamiltonian; (2) is easy to optimize; and (3) can be implemented by a shallow circuit, within the capabilities of a NISQ computer. The most widely used, unitary coupled cluster (UCC), type of ans\"atze mathematically correspond to a product of unitary evolutions of fermionic excitation operators. Owing to their fermionic structure, UCC ans\"atze preserve the symmetries of electronic wavefunctions, and thus are accurate and easy to optimize. Nevertheless, UCC ans\"atze are implemented by high depth circuits, which severely limit the size of the molecules that can be reliably simulated on NISQ computers. In this thesis, I begin by constructing efficient quantum circuits to perform evolutions of fermionic excitation operators. The circuits are optimized in the number of two-qubit entangling gates, which are the current bottleneck of NISQ computers. Compared to the standard circuits used to implement evolutions of fermionic excitation operators, the circuits derived in this thesis reduce the number of two-qubit entangling gates by more than
70% on average. As an intermediate result, I also derive efficient circuits to perform evolutions of qubit excitation operators (excitation operators that account for qubit, rather than fermionic commutation relations). Even with the fermionic-excitation-evolution circuits derived here, UCC ans\"atze still require very long circuits, with a particularly large number of two-qubit entangling gates. In this thesis, I consider the use of alternative VQE ans\"atze, based on evolutions of qubit excitation operators. Due to not accounting for fermionic anticommutation, evolutions of qubit excitation operators can be performed by circuits that require asymptotically fewer two-qubit entangling gates. Furthermore, qubit excitation operators preserve many of the physical properties of fermionic excitation operators. Performing a number of classical numerical VQE simulations for small molecules, I show that qubit-excitation-based ans\"atze can approximate molecular electronic wavefunctions almost as accurately as fermionic-excitation-based ans\"atze. Hence, I argue that evolutions of qubit excitation operators are more suitable to construct molecular ans\"atze than evolutions of fermionic excitation operators, especially in the era of NISQ computers. Motivated by the advantage of qubit-excitation-based ans\"atze, I introduce the qubit-excitation-based adaptive variational quantum eigensolver (QEB-ADAPT-VQE). The QEB-ADAPT-VQE belongs to a family of ADAPT-VQE protocols for molecular modelling that grow a problem-tailored ansatz by iteratively appending unitary operators sampled from a predefined finite-size pool of operators. The operator at each iteration is sampled based on an ansatz-growing strategy, which aims to achieve the lowest estimate for the Hamiltonian expectation value at each iteration. In this way, ADAPT-VQE protocols construct shallow-circuit, few-parameter ans\"atze tailored specifically to the molecular systems of interest. In the case of the QEB-ADAPT-VQE, the operator pool is defined by a set of evolutions of single and double qubit excitation operators. I benchmark the performance of the QEB-ADAPT-VQE, by performing classical numerical simulations. I demonstrate that it can construct ans\"atze that are several orders of magnitude more accurate, and require significantly shallower circuits, than standard UCC ans\"atze. I also compare the QEB-ADAPT-VQE against the original fermionic-ADAPT-VQE, which utilizes a pool of fermionic excitation evolutions, and the qubit-ADAPT-VQE, which utilizes a pool of Pauli-string evolutions. I demonstrate that, in terms of circuit efficiency and convergence speed, the QEB-ADAPT-VQE systematically outperforms the qubit-ADAPT-VQE, which to my knowledge was the previous most circuit-efficient, scalable VQE protocol for molecular modeling. The QEB-ADAPT-VQE protocol, therefore represents a significant improvement in the field of VQE protocols for molecular modelling and brings us closer to achieving practical quantum advantage. Lastly, I outline a modified version of the QEB-ADAPT-VQE, the excited-QEB-ADAPT-VQE, designed to estimate energies of excited molecular states. The excited-QEB-ADAPT-VQE is more robust to initial simulation conditions, at the expense of increased computational complexity.
Some of the computational limitations in solving the nuclear many-body problem could be overcome by utilizing quantum computers. The nuclear shell-model calculations providing deeper insights into the properties of atomic nuclei, is one such case with high demand for resources as the size of the Hilbert space grows exponentially with the number of particles involved. Quantum algorithms are being developed to overcome these challenges and advance such calculations. To develop quantum circuits for the nuclear shell-model, leveraging the capabilities of noisy intermediate-scale quantum (NISQ) devices. We aim to minimize resource requirements (specifically in terms of qubits and gates) and strive to reduce the impact of noise by employing relevant mitigation techniques. We achieve noise resilience by designing an optimized ansatz for the variational quantum eigensolver (VQE) based on Givens rotations and incorporating qubit-ADAPT-VQE in combination with variational quantum deflation (VQD) to compute ground and excited states incorporating the zero-noise extrapolation mitigation technique. Furthermore, the qubit requirements are significantly reduced by mapping the basis states to qubits using Gray code encoding and generalizing transformations of fermionic operators to efficiently represent manybody states. By employing the noise-resilient protocols, we achieve the ground and excited state energy levels of 38Ar and 6Li with better accuracy. These energy levels are presented for noiseless simulations, noisy conditions, and after applying noise mitigation techniques. Results are compared for Jordan Wigner and Gray code encoding using VQE, qubit-ADAPT-VQE, and VQD. Our work highlights the potential of noise-resilient protocols to leverage the full potential of NISQ devices in scaling the nuclear shell model calculations.
Compact representations of fermionic Hamiltonians are necessary to perform calculations on quantum computers that lack error-correction. A fermionic system is typically defined within a subspace of fixed particle number and spin while unnecessary states are projected out of the Hilbert space. We provide a bijective mapping using combinatoric ranking to bijectively map fermion basis states to qubit basis states and express operators in the standard spin representation. We then evaluate compact mapping using the Variational Quantum Eigensolver (VQE) with the unitary coupled cluster singles and doubles excitations (UCCSD) ansatz in the compact representation. Compactness is beneficial when the orbital filling is well away from half, and we show at 30 spin orbital
H2 calculation with only 8 qubits. We find that the gate depth needed to prepare the compact wavefunction is not much greater than the full configuration space in practice. A notable observation regards the number of calls to the optimizer needed for the compact simulation compared to the full simulation. We find that the compact representation converges faster than the full representation using the ADAM optimizer in all cases. Our analysis demonstrates the effect of compact mapping in practice.
Modified divacancies in the 4H polytype of silicon carbide (SiC) exhibit enhanced charge stability and spin addressability at room temperature, making them attractive for quantum applications. However, their low formation yield and lack of direct structural identification have hindered progress. Here, we demonstrate a controllable method for high-yield engineering and identification of oxygen-related modified divacancy color centers in 4H-SiC via oxygen-ion implantation. Based on their distinct optical and spin-resonance characteristics, we experimentally resolve four types of modified divacancies. Furthermore, by measuring isotope-resolved 17O hyperfine interactions, we identify them as the four crystallographic configurations of oxygen-vacancy (OV) complexes. Remarkably, single OV centers account for over 90% of the total defect population and exhibit superior optical properties and spin coherence compared with defects created by conventional carbon or nitrogen implantation. We characterize the zero-phonon lines of these OV centers and reveal distinct temperature-dependent behavior in spin-readout contrast. By optimizing implantation dose and annealing temperature, we achieve high-density ensembles and observe Rabi-oscillation beating patterns associated with different orientations of basal-type defects. These results establish a high-yield route for scalable engineering of these four oxygen-related modified divacancies in 4H-SiC and clarify their atomic structure, opening new opportunities for solid-state quantum technologies.