Aleksei Vladimirovich Ivanov

Quantum simulation of materials is a promising application area of quantum computers. To practically realize this promise, we must reduce quantum resources while maintaining accuracy. In electronic structure calculations on classical computers, resource reduction has been achieved by using the projector augmented-wave method (PAW) and plane wave basis sets. However, the PAW method generalized for many-body states introduces nonorthogonality effects which impede its direct application to quantum computing. In this work, we develop a unitary variant of the PAW (UPAW) that preserves the orthogonality constraints. We provide a linear-combination-of-unitaries decomposition of the UPAW Hamiltonian to enable ground state estimation using qubitized quantum phase estimation. Additionally, we further improve algorithmic efficiency by extending classical down-sampling techniques into the quantum setting. We then estimate quantum resources for crystalline solids to estimate the energy within chemical accuracy with respect to the full basis set limit, and also consider a supercell approach which is more suitable for calculations of defect states. We provide the quantum resources for energy estimation of a nitrogen-vacancy defect center in diamond which is a challenging system for classical algorithms and a quintessential problem in the studies of quantum point defects.

The decomposition of a square matrix into a sum of Pauli strings is a classical pre-processing step required to realize many quantum algorithms. Such a decomposition requires significant computational resources for large matrices. We present an exact and explicit formula for the Pauli string coefficients which inspires an efficient algorithm to compute them. More specifically, we show that up to a permutation of the matrix elements, the decomposition coefficients are related to the original matrix by a multiplication of a generalised Hadamard matrix. This allows one to use the Fast Walsh-Hadamard transform and calculate all Pauli decomposition coefficients in O ( N 2 log ⁡ N ) time and using O ( 1 ) additional memory, for an N × N matrix. A numerical implementation of our equation outperforms currently available solutions.

We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for the implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE), providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe–Salpeter Equation, variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn–Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support for graphics processing unit (GPU) acceleration has been achieved with minor modifications to the GPAW code thanks to the CuPy library. We end the review with an outlook, describing some future plans for GPAW.

Elucidation of the mechanism for optical spin initialization of point defects in solids in the context of quantum applications requires an accurate description of the excited electronic states involved. While variational density functional calculations have been successful in describing the ground state of a great variety of systems, doubts have been expressed in the literature regarding the ability of such calculations to describe electronic excitations of point defects. A direct orbital optimization method is used here to perform time-independent, variational density functional calculations of a prototypical defect, the negatively charged nitrogen-vacancy center in diamond. The calculations include up to 511 atoms subject to periodic boundary conditions and the excited state calculations require similar computational effort as ground state calculations. Contrary to some previous reports, the use of local and semilocal density functionals gives the correct ordering of the low-lying triplet and singlet states, namely {}^{3}A_2 < {}^{1}E < {}^{1}A_1 < {}^{3}E 3 A 2 < 1 E < 1 A 1 < 3 E . Furthermore, the more advanced meta generalized gradient approximation functionals give results that are in remarkably good agreement with high-level, many-body calculations as well as available experimental estimates, even for the excited singlet state which is often referred to as having multireference character. The lowering of the energy in the triplet excited state as the atom coordinates are optimized in accordance with analytical forces is also close to the experimental estimate and the resulting zero-phonon line triplet excitation energy is underestimated by only 0.15 eV. The approach used here is found to be a promising tool for studying electronic excitations of point defects in, for example, systems relevant for quantum technologies.

Theoretical studies of photochemical processes require a description of the energy surfaces of excited electronic states, especially near degeneracies, where transitions between states are most likely. Systems relevant to photochemical applications are typically too large for high-level multireference methods, and while time-dependent density functional theory (TDDFT) is efficient, it can fail to provide the required accuracy. A variational, time-independent density functional approach is applied to the twisting of the double bond and pyramidal distortion in ethylene, the quintessential model for photochemical studies. By allowing for symmetry breaking, the calculated energy surfaces exhibit the correct topology around the twisted-pyramidalized conical intersection even when using a semilocal functional approximation, and by including explicit self-interaction correction, the torsional energy curves are in close agreement with published multireference results. The findings of the present work point to the possibility of using a single determinant time-independent density functional approach to simulate nonadiabatic dynamics, even for large systems where multireference methods are impractical and TDDFT is often not accurate enough.

Fermi–Löwdin orbitals (FLOs) are a special set of localized orbitals, which have become commonly used in combination with the Perdew–Zunger self-interaction correction (SIC) in the FLO-SIC method. The FLOs are obtained for a set of occupied orbitals by specifying a classical position for each electron. These positions are known as Fermi-orbital descriptors (FODs), and they have a clear relation to chemical bonding. In this study, we show how FLOs and FODs can be used to initialize, interpret, and justify SIC solutions in a common chemical picture, both within FLO-SIC and in traditional variational SIC, and to locate distinct local minima in either of these approaches. We demonstrate that FLOs based on Lewis theory lead to symmetry breaking for benzene—the electron density is found to break symmetry already at the symmetric molecular structure—while ones from Linnett’s double-quartet theory reproduce symmetric electron densities and molecular geometries. Introducing a benchmark set of 16 planar cyclic molecules, we show that using Lewis theory as the starting point can lead to artifactual dipole moments of up to 1 D, while Linnett SIC dipole moments are in better agreement with experimental values. We suggest using the dipole moment as a diagnostic of symmetry breaking in SIC and monitoring it in all SIC calculations. We show that Linnett structures can often be seen as superpositions of Lewis structures and propose Linnett structures as a simple way to describe aromatic systems in SIC with reduced symmetry breaking. The role of hovering FODs is also briefly discussed.