Optimization of highly excited matrix product states with an application to vibrational spectroscopy

SCINE-Software for chemical interaction networks

The software for chemical interaction networks (SCINE) project aims at pushing the frontier of quantum chemical calculations on molecular structures to a new level. While calculations on individual structures as well as on simple relations between them have become routine in chemistry, new developments have pushed the frontier in the field to high-throughput calculations. Chemical relations may be created by a search for specific molecular properties in a molecular design attempt, or they can be defined by a set of elementary reaction steps that form a chemical reaction network. The software modules of SCINE have been designed to facilitate such studies. The features of the modules are (i) general applicability of the applied methodologies ranging from electronic structure (no restriction to specific elements of the periodic table) to microkinetic modeling (with little restrictions on molecularity), full modularity so that SCINE modules can also be applied as stand-alone programs or be exchanged for external software packages that fulfill a similar purpose (to increase options for computational campaigns and to provide alternatives in case of tasks that are hard or impossible to accomplish with certain programs), (ii) high stability and autonomous operations so that control and steering by an operator are as easy as possible, and (iii) easy embedding into complex heterogeneous environments for molecular structures taken individually or in the context of a reaction network. A graphical user interface unites all modules and ensures interoperability. All components of the software have been made available as open source and free of charge.

Harmonic and anharmonic vibrational computations for biomolecular building blocks: Benchmarking DFT and basis sets by theoretical and experimental IR spectrum of glycine conformers

Advanced vibrational spectroscopic experiments have reached a level of sophistication that can only be matched by numerical simulations in order to provide an unequivocal analysis, a crucial step to understand the structure-function relationship of biomolecules. While density functional theory (DFT) has become the standard method when targeting medium-size or larger systems, the problem of its reliability and accuracy are well-known and have been abundantly documented. To establish a reliable computational protocol, especially when accuracy is critical, a tailored benchmark is usually required. This is generally done over a short list of known candidates, with the basis set often fixed a priori. In this work, we present a systematic study of the performance of DFT-based hybrid and double-hybrid functionals in the prediction of vibrational energies and infrared intensities at the harmonic level and beyond, considering anharmonic effects through vibrational perturbation theory at the second order. The study is performed for the six-lowest energy glycine conformers, utilizing available "state-of-the-art" accurate theoretical and experimental data as reference. Focusing on the most intense fundamental vibrations in the mid-infrared range of glycine conformers, the role of the basis sets is also investigated considering the balance between computational cost and accuracy. Targeting larger systems, a broad range of hybrid schemes with different computational costs is also tested.

Flexible DMRG-Based Framework for Anharmonic Vibrational Calculations

We present a novel formulation of the vibrational density matrix renormalization group (vDMRG) algorithm tailored to strongly anharmonic molecules described by general, high-dimensional model representations of potential energy surfaces. For this purpose, we extend the vDMRG framework to support vibrational Hamiltonians expressed in the so-called n-mode second-quantization formalism. The resulting n-mode vDMRG method offers full flexibility with respect to both the functional form of the PES and the choice of the single-particle basis set. We leverage this framework to apply, for the first time, vDMRG based on an anharmonic modal basis set optimized with the vibrational self-consistent field algorithm on an on-the-fly constructed PES. We also extend the n-mode vDMRG framework to include excited-state-targeting algorithms in order to efficiently calculate anharmonic transition frequencies. We demonstrate the capabilities of our novel n-mode vDMRG framework for methyloxirane, a challenging molecule with 24 coupled vibrational modes.

Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, Birnoschi L, Blanco-González A, Bokarev SI, Broer R, Cacciari R, Calio PB, Carlson RK, Carvalho Couto R, Cerdán L, Chibotaru LF, Chilton NF, Church JR, Conti I, Coriani S, Cuéllar-Zuquin J, Daoud RE, Dattani N, Decleva P, de Graaf C, Delcey M, De Vico L, Dobrautz W, Dong SS, Feng R, Ferré N, Filatov(Gulak) M, Gagliardi L, Garavelli M, González L, Guan Y, Guo M, Hennefarth MR, Hermes MR, Hoyer CE, Huix-Rotllant M, Jaiswal VK, Kaiser A, Kaliakin DS, Khamesian M, King DS, Kochetov V, Krośnicki M, Kumaar AA, Larsson ED, Lehtola S, Lepetit MB, Lischka H, López Ríos P, Lundberg M, Ma D, Mai S, Marquetand P, Merritt ICD, Montorsi F, Mörchen M, Nenov A, Nguyen VHA, Nishimoto Y, Oakley MS, Olivucci M, Oppel M, Padula D, Pandharkar R, Phung QM, Plasser F, Raggi G, Rebolini E, Reiher M, Rivalta I, Roca-Sanjuán D, Romig T, Safari AA, Sánchez-Mansilla A, Sand AM, Schapiro I, Scott TR, Segarra-Martí J, Segatta F, Sergentu DC, Sharma P, Shepard R, Shu Y, Staab JK, Straatsma TP, Sørensen LK, Tenorio BNC, Truhlar DG, Ungur L, Vacher M, Veryazov V, Voß TA, Weser O, Wu D, Yang X, Yarkony D, Zhou C, Zobel JP, Lindh R. The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

Synthesis of Hidden Subgroup Quantum Algorithms and Quantum Chemical Dynamics

We describe a general formalism for quantum dynamics and show how this formalism subsumes several quantum algorithms, including the Deutsch, Deutsch-Jozsa, Bernstein-Vazirani, Simon, and Shor algorithms as well as the conventional approach to quantum dynamics based on tensor networks. The common framework exposes similarities among quantum algorithms and natural quantum phenomena: we illustrate this connection by showing how the correlated behavior of protons in water wire systems that are common in many biological and materials systems parallels the structure of Shor's algorithm.

Ground and excited states of even-numbered Hubbard ring at half-filling: comparison of the extended Gutzwiller approach with exact diagonalization

It remains a great challenge in condensed matter physics to develop a method to treat strongly correlated many-body systems with balanced accuracy and efficiency. We introduce an extended Gutzwiller (EG) method incorporating a manifold technique, which builds an effective manifold of the many-body Hilbert space, to describe the ground-state (GS) and excited-state (ES) properties of strongly correlated electrons. We systematically apply an EG projector onto the GS and ES of a non-interacting system. Diagonalization of the true Hamiltonian within the manifold formed by the resulting EG wavefunctions gives the approximate GS and ES of the correlated system. To validate this technique, we implement it on even-numbered fermionic Hubbard rings at half-filling with periodic boundary conditions, and compare the results with the exact diagonalization (ED) method. The EG method is capable of generating high-quality GS and low-lying ES wavefunctions, as evidenced by the high overlaps of wavefunctions between the EG and ED methods. Favorable comparisons are also achieved for other quantities including the total energy, the double occupancy, the total spin and the staggered magnetization. With the capability of accessing the ESs, the EG method can capture the essential features of the one-electron removal spectral function that contains contributions from states deep in the excited spectrum. Finally, we provide an outlook on the application of this method on large extended systems.

Ultrafast Charge Relocation Dynamics in Enol-Keto Tautomerization Monitored with a Local Soft-X-ray Probe

Proton-coupled electron transfer (PCET) is the underlying mechanism governing important reactions ranging from water splitting in photosynthesis to oxygen reduction in hydrogen fuel cells. The interplay of proton and electronic charge distribution motions can vary from sequential to concerted schemes, with elementary steps occurring on ultrafast time scales. We demonstrate with a simulation study that femtosecond soft-X-ray spectroscopy provides key insights into the PCET mechanism of a photoinduced intramolecular enol* → keto* tautomerization reaction. A full quantum treatment of the electronic and nuclear dynamics of 2-(2'-hydroxyphenyl)benzothiazole upon electronic excitation reveals how spectral signatures of local excitations from core to frontier orbitals display the distinctly different stages of charge relocation for the H atom, donating, and accepting sites. Our findings indicate that ultraviolet/X-ray pump-probe spectroscopy provides a unique way to probe ultrafast electronic structure rearrangements in photoinduced chemical reactions essential to understanding the mechanism of PCET.

Relativistic Kramers-Unrestricted Exact-Two-Component Density Matrix Renormalization Group

In this work we develop a variational relativistic density matrix renormalization group (DMRG) approach within the exact-two-component (X2C) framework (X2C-DMRG), using spinor orbitals optimized with the two-component relativistic complete active space self-consistent field. We investigate fine-structure splittings of p- (Ga, In, Tl) and d-block (Sc, Y, La) atoms and excitation energies of monohydride molecules (GeH, SnH, and TlH) with X2C-DMRG calculations using an all-electron relativistic Hamiltonian in a Kramers-unrestricted basis. We find that X2C-DMRG yields accurate ²P and ²D splittings compared to multireference configuration interaction with singles and doubles (MRCISD). We also investigated the degree of symmetry breaking in the atomic multiplets and convergence of electron correlation in the total energies. Symmetry breaking can be large in some cases (∼30 meV); however, increasing the number of renormalized block states m for the DMRG optimization recovers the symmetry breaking by several orders of magnitude. Encouragingly, we find the convergence of electron correlation to be close to MRCISDTQ5 quality. Relativistic X2C-DMRG approaches are important for cases where spin-orbit coupling is significant and the underlying reference wave function requires a large determinantal space. We are able to obtain quantitatively correct fine-structure splittings for systems up to 10¹⁹ number of determinants with traditional CI approaches, which are currently unfeasible to converge for the field.

Quantum Proton Effects from Density Matrix Renormalization Group Calculations

We recently introduced [J. Chem. Phys. 2020, 152, 204103] the nuclear-electronic all-particle density matrix renormalization group (NEAP-DMRG) method to solve the molecular Schrödinger equation, based on a stochastically optimized orbital basis, without invoking the Born-Oppenheimer approximation. In this work, we combine the DMRG method with the nuclear-electronic Hartree-Fock (NEHF-DMRG) approach, treating nuclei and electrons on the same footing. Inter- and intraspecies correlations are described within the DMRG method without truncating the excitation degree of the full configuration interaction wave function. We extend the concept of orbital entanglement and mutual information to nuclear-electronic wave functions and demonstrate that they are reliable metrics to detect strong correlation effects. We apply the NEHF-DMRG method to the HeHHe⁺ molecular ion, to obtain accurate proton densities, ground-state total energies, and vibrational transition frequencies by comparison with state-of-the-art data obtained with grid-based approaches and modern configuration interaction methods. For HCN, we improve on the accuracy of the latter approaches with respect to both the ground-state absolute energy and proton density, which is a major challenge for multireference nuclear-electronic state-of-the-art methods.

Multicomponent heat-bath configuration interaction with the perturbative correction for the calculation of protonic excited states

In this study, we extend the multicomponent heat-bath configuration interaction (HCI) method to excited states. Previous multicomponent HCI studies have been performed using only the variational stage of the HCI algorithm as they have largely focused on the calculation of protonic densities. Because this study focuses on energetic quantities, a second-order perturbative correction after the variational stage is essential. Therefore, this study implements the second-order Epstein-Nesbet correction to the variational stage of multicomponent HCI for the first time. Additionally, this study introduces a new procedure for calculating reference excitation energies for multicomponent methods using the Fourier-grid Hamiltonian (FGH) method, which should allow the one-particle electronic basis set errors to be better isolated from errors arising from an incomplete description of electron-proton correlation. The excited-state multicomponent HCI method is benchmarked by computing protonic excitations of the HCN and FHF^- molecules and is shown to be of similar accuracy to previous excited-state multicomponent methods such as the multicomponent time-dependent density-functional theory and equation-of-motion coupled-cluster theory relative to the new FGH reference values.

Computing vibrational energy levels by solving linear equations using a tensor method with an imposed rank

Present day computers do not have enough memory to store the high-dimensional tensors required when using a direct product basis to compute vibrational energy levels of a polyatomic molecule with more than about five atoms. One way to deal with this problem is to represent tensors using a tensor format. In this paper, we use the canonical polyadic (CP) format. Energy levels are computed by building a basis from vectors obtained by solving linear equations. The method can be thought of as a CP realization of a block inverse iteration method with multiple shifts. The CP rank of the tensors is fixed, and the linear equations are solved with an method. There is no need for rank reduction and no need for orthogonalization, and tensors with a rank larger than the fixed rank used to solve the linear equations are never generated. The ideas are tested by computing vibrational energy levels of a 64-D bilinearly coupled model Hamiltonian and of acetonitrile (12-D).

Excited-State DMRG Made Simple with FEAST

We introduce DMRG[FEAST], a new method for optimizing excited-state many-body wave functions with the density matrix renormalization group (DMRG) algorithm. Our approach applies the FEAST algorithm, originally designed for large-scale diagonalization problems, to matrix product state wave functions. We show that DMRG[FEAST] enables the stable optimization of both low- and high-energy eigenstates, therefore overcoming the limitations of state-of-the-art excited-state DMRG algorithms. We demonstrate the reliability of DMRG[FEAST] by calculating anharmonic vibrational excitation energies of molecules with up to 30 fully coupled degrees of freedom.

Electron Dynamics with the Time-Dependent Density Matrix Renormalization Group

In this work, we simulate the electron dynamics in molecular systems with the time-dependent density matrix renormalization group (TD-DMRG) algorithm. We leverage the generality of the so-called tangent-space TD-DMRG formulation and design a computational framework in which the dynamics is driven by the exact nonrelativistic electronic Hamiltonian. We show that by parametrizing the wave function as a matrix product state, we can accurately simulate the dynamics of systems including up to 20 electrons and 32 orbitals. We apply the TD-DMRG algorithm to three problems that are hardly targeted by time-independent methods: the calculation of molecular (hyper)polarizabilities, the simulation of electronic absorption spectra, and the study of ultrafast ionization dynamics.

Transcorrelated density matrix renormalization group

We introduce the transcorrelated Density Matrix Renormalization Group (tcDMRG) theory for the efficient approximation of the energy for strongly correlated systems. tcDMRG encodes the wave function as a product of a fixed Jastrow or Gutzwiller correlator and a matrix product state. The latter is optimized by applying the imaginary-time variant of time-dependent (TD) DMRG to the non-Hermitian transcorrelated Hamiltonian. We demonstrate the efficiency of tcDMRG with the example of the two-dimensional Fermi-Hubbard Hamiltonian, a notoriously difficult target for the DMRG algorithm, for different sizes, occupation numbers, and interaction strengths. We demonstrate fast energy convergence of tcDMRG, which indicates that tcDMRG could increase the efficiency of standard DMRG beyond quasi-monodimensional systems and provides a generally powerful approach toward the dynamic correlation problem of DMRG.

A general automatic method for optimal construction of matrix product operators using bipartite graph theory

Constructing matrix product operators (MPOs) is at the core of the modern density matrix renormalization group (DMRG) and its time dependent formulation. For the DMRG to be conveniently used in different problems described by different Hamiltonians, in this work, we propose a new generic algorithm to construct the MPO of an arbitrary operator with a sum-of-products form based on the bipartite graph theory. We show that the method has the following advantages: (i) it is automatic in that only the definition of the operator is required; (ii) it is symbolic thus free of any numerical error; (iii) the complementary operator technique can be fully employed so that the resulting MPO is globally optimal for any given order of degrees of freedom; and (iv) the symmetry of the system could be fully employed to reduce the dimension of MPO. To demonstrate the effectiveness of the new algorithm, the MPOs of Hamiltonians ranging from the prototypical spin-boson model and the Holstein model to the more complicated ab initio electronic Hamiltonian and the anharmonic vibrational Hamiltonian with the sextic force field are constructed. It is found that for the former three cases, our automatic algorithm can reproduce exactly the same MPOs as the optimally hand-crafted ones already known in the literature.

Hardware efficient quantum algorithms for vibrational structure calculations

We introduce a framework for the calculation of ground and excited state energies of bosonic systems suitable for near-term quantum devices and apply it to molecular vibrational anharmonic Hamiltonians. Our method supports generic reference modal bases and Hamiltonian representations, including the ones that are routinely used in classical vibrational structure calculations. We test different parametrizations of the vibrational wavefunction, which can be encoded in quantum hardware, based either on heuristic circuits or on the bosonic Unitary Coupled Cluster Ansatz. In particular, we define a novel compact heuristic circuit and demonstrate that it provides a good compromise in terms of circuit depth, optimization costs, and accuracy. We evaluate the requirements, number of qubits and circuit depth, for the calculation of vibrational energies on quantum hardware and compare them with state-of-the-art classical vibrational structure algorithms for molecules with up to seven atoms.

Nuclear-electronic all-particle density matrix renormalization group

We introduce the Nuclear-Electronic All-Particle Density Matrix Renormalization Group (NEAP-DMRG) method for solving the time-independent Schrödinger equation simultaneously for electrons and other quantum species. In contrast to the already existing multicomponent approaches, in this work, we construct from the outset a multi-reference trial wave function with stochastically optimized non-orthogonal Gaussian orbitals. By iterative refining of the Gaussians' positions and widths, we obtain a compact multi-reference expansion for the multicomponent wave function. We extend the DMRG algorithm to multicomponent wave functions to take into account inter- and intra-species correlation effects. The efficient parameterization of the total wave function as a matrix product state allows NEAP-DMRG to accurately approximate the full configuration interaction energies of molecular systems with more than three nuclei and 12 particles in total, which is currently a major challenge for other multicomponent approaches. We present the NEAP-DMRG results for two few-body systems, i.e., H₂ and H₃ ⁺, and one larger system, namely, BH₃.

Elucidating the multi-configurational character of the firefly dioxetanone anion and its prototypes in the biradical region using full valence active spaces

We analyzed the near-degenerate states of the firefly dioxetanone anion (FDO-) and its prototypes, especially in the biradical region, using multi-configurational approaches. The importance of utilizing full valence active spaces by means of density-matrix renormalization group self-consistent field (DMRG-SCF) calculations was described. Our results revealed that the neglect of some valence orbitals can affect the quantitative accuracy in later multi-reference calculations or the qualitative conclusion when optimizing conical intersections. Using all of the relevant valence orbitals of FDO-, we confirmed that there were two conical intersections, as reported in previous work, and that the intersecting states were changed when the active space was enlarged. Beyond these, we found that there were strong interactions between states in the biradical regions, in which the changes in entanglements can be used to visualize the interacting state evolution.