Localized Quantum Chemistry on Quantum Computers

Quantum chemistry calculations of large, strongly correlated systems are typically limited by the computation cost that scales exponentially with the size of the system. Quantum algorithms, designed specifically for quantum computers, can alleviate this, but the resources required are still too large for today's quantum devices. Here, we present a quantum algorithm that combines a localization of multireference wave functions of chemical systems with quantum phase estimation (QPE) and variational unitary coupled cluster singles and doubles (UCCSD) to compute their ground-state energy. Our algorithm, termed "local active space unitary coupled cluster" (LAS-UCC), scales linearly with the system size for certain geometries, providing a polynomial reduction in the total number of gates compared with QPE, while providing accuracy above that of the variational quantum eigensolver using the UCCSD ansatz and also above that of the classical local active space self-consistent field. The accuracy of LAS-UCC is demonstrated by dissociating (H₂)₂ into two H₂ molecules and by breaking the two double bonds in trans-butadiene, and resource estimates are provided for linear chains of up to 20 H₂ molecules.

The evolution of paramagnetic NMR as a tool in structural biology

Paramagnetic NMR data contain extremely accurate long-range information on metalloprotein structures and, when used in the frame of integrative structural biology approaches, they allow for the retrieval of structural details to a resolution that is not achievable using other techniques. Paramagnetic data thus represent an extremely powerful tool to refine protein models in solution, especially when coupled to X-ray or cryoelectron microscopy data, to monitor the formation of complexes and determine the relative arrangements of their components, and to highlight the presence of conformational heterogeneity. More recently, theoretical and computational advancements in quantum chemical calculations of paramagnetic NMR observables are progressively opening new routes in structural biology, because they allow for the determination of the structure within the coordination sphere of the metal center, thus acting as a loupe on sites that are difficult to observe but very important for protein function.

A combined first- and second-order optimization method for improving convergence of Hartree–Fock and Kohn–Sham calculations

In this work, we investigate the optimization of Hartree–Fock (HF) orbitals with our recently proposed combined first- and second-order (SO-SCI) method, which was originally developed for multi-configuration self-consistent field (MCSCF) and complete active space SCF (CASSCF) calculations. In MCSCF/CASSCF, it unites a second-order optimization of the active orbitals with a Fock-based first-order treatment of the remaining closed-virtual orbital rotations. In the case of the single-determinant wavefunctions, the active space is replaced by a preselected “second-order domain,” and all rotations involving orbitals in this subspace are treated at second-order. The method has been implemented for spin-restricted and spin-unrestricted Hartree–Fock (RHF, UHF), configuration-averaged Hartree–Fock (CAHF), as well as Kohn–Sham (KS) density functional theory (RKS, UKS). For each of these cases, various choices of the second-order domain have been tested, and appropriate defaults are proposed. The performance of the method is demonstrated for several transition metal complexes. It is shown that the SO-SCI optimization provides faster and more robust convergence than the standard SCF procedure but requires, in many cases, even less computation time. In difficult cases, the SO-SCI method not only speeds up convergence but also avoids convergence to saddle-points. Furthermore, it helps to find spin-symmetry broken solutions in the cases of UHF or UKS. In the case of CAHF, convergence can also be significantly improved as compared to a previous SCF implementation. This is particularly important for multi-center cases with two or more equal heavy atoms. The performance is demonstrated for various two-center complexes with different lanthanide atoms.

Are lanthanide-transition metal direct bonds a route to achieving new generation {3d-4f} SMMs?

Lanthanide based single-molecule magnets are gaining wide attention due to their potential applications in emerging technologies. One of the main challenges in this area is quenching quantum tunnelling of magnetisation (QTM), which often undercuts the magnetisation reversal barrier. Among the several strategies employed, enhancing exchange coupling has been studied in detail, with large exchanges resulting in stronger quenching of QTM effects. Lanthanides, however, suffer from weak exchanges offered by the deeply buried 4f orbitals and the numerous attempts to enhance the exchange coupling in the {3d-4f} pairs have not exceeded values larger than 30 cm^-1. In this work, using a combination of DFT and the ab initio CASSCF/RASSI-SO method, we have explored lanthanide-transition metal direct bonds as a tool to quench QTM effects. In this direction, we have modelled [PyCp₂LnMCp(CO)₂] (Ln = Gd(III), Dy(III), and Er(III) and M = V(0), Mn(0), Co(0) and Fe(I) and here PyCp₂ = [2,6-(CH₂C₅H₃)₂C₅H₃N]₂^- using [PyCp₂DyFeCp(CO)₂] as an example as reported by Nippe et al. (C. P. Burns, X. Yang, J. D. Wofford, N. S. Bhuvanesh, M. B. Hall and M. Nippe, Angew. Chem., Int. Ed. 2018, 57, 8144). Bonding analysis reveals a dative Ln-TM bond with a donation of π(V/Mnd_xy-π*CO) to 5d_z² (Gd) in the case of Gd-V and Gd-Mn and 4s(Co) to 5d_xy/5d_yz (Gd) for Gd-Co with the transition metal ion being found in the low-spin S = ½ configurations in all the cases. B3LYP/TZV (Gd;CSDZ) calculations on [PyCp₂GdMCp(CO)₂] yield J_Gd-V = -46.1 cm^-1, J_Gd-Mn = -57.1 cm^-1, J_Gd-Co = +55.3 cm^-1, J_Gd-Fe⁺ = +13.9 cm^-1, J_Gd-Vhs = -162.1 cm^-1 and J_Gd-Mnhs = -343.9 cm^-1 and unveiling record-high J values for {3d-4f} complexes. The mechanism of magnetic coupling is developed, which discloses the dominating presence of strong 3d-4f orbital overlaps in most of the cases studied, leading to antiferromagnetic exchange. When these overlaps are weaker and 3d to Gd(5d_z²), charge transfer dominates, yielding a ferromagnetic coupling for the Gd-Co/Gd-Fe⁺ complexes. Calculations performed on the anisotropic Dy(III) and Er(III) complexes reveal that the ground state g_zz axis lies along the Cp-Ln-Cp axis and the Ln-TM bonds, respectively. Thus the Ln-TM bond hinders the single-ion anisotropy of Dy(III) by offering equatorial ligation and lowering the m_J = ±½ state energy, and at the same time, helping in enhancing the axiality of Er(III). When strong {3d-4f} exchange couplings are introduced, record-high barrier heights as high as 229 cm^-1 were accomplished. Furthermore, the exchange coupling annihilates the QTM effects and suggests the lanthanide-transition metal direct bond as a viable alternative to enhance exchange coupling to bring {3d-4f} complexes back in the race for high-blocking SMMs.

NMR for Single Ion Magnets

Nuclear Magnetic Resonance is particularly sensitive to the electronic structure of matter and is thus a powerful tool to characterize in-depth the magnetic properties of a system. NMR is indeed increasingly recognized as an ideal tool to add precious structural information for the development of Single Ion Magnets, small complexes that are recently gaining much popularity due to their quantum computing and spintronics applications. In this review, we recall the theoretical principles of paramagnetic NMR, with particular attention to lanthanoids, and we give an overview of the recent advances in this field.

Spin dynamics in single-molecule magnets and molecular qubits

Over recent decades, much effort has been made to lengthen spin relaxation/decoherence times of single-molecule magnets and molecular qubits by following different chemical design rules such as maximizing the total spin value, controlling symmetry, enhancing the ligand field or inhibiting key vibrational modes. Simultaneously, electronic structure calculations have been employed to provide an understanding of the processes involved in the spin dynamics of molecular systems and served to refine or introduce new design rules. This review focuses on contemporary theoretical approaches focused on the calculation of spin relaxation/decoherence times, highlighting their main features and scope. Fundamental aspects of experimental techniques for the determination of key Single Molecule Magnet/Spin Qubit properties are also reviewed.

The Molpro quantum chemistry package

Molpro is a general purpose quantum chemistry software package with a long development history. It was originally focused on accurate wavefunction calculations for small molecules but now has many additional distinctive capabilities that include, inter alia, local correlation approximations combined with explicit correlation, highly efficient implementations of single-reference correlation methods, robust and efficient multireference methods for large molecules, projection embedding, and anharmonic vibrational spectra. In addition to conventional input-file specification of calculations, Molpro calculations can now be specified and analyzed via a new graphical user interface and through a Python framework.

Magnetic exchange interactions in symmetric lanthanide dimetallics

Multi-frequency EPR spectra and CASSCF-SO calculations on two symmetric homo-dimetallic lanthanide complexes are used to determine the magnetic exchange coupling in the low-lying states.