Nonadiabatic Molecular Dynamics with Fermionic Subspace-Expansion Algorithms on Quantum Computers.

We introduce a novel computational framework for excited-state molecular quantum dynamics simulations driven by quantum-computing-based electronic-structure calculations. This framework leverages the fewest-switches surface-hopping method for simulating the nuclear dynamics and calculates the required excited-state transition properties with different flavors of the quantum subspace expansion and quantum equation-of-motion algorithms. We apply our method to simulate the collision reaction between a hydrogen atom and a hydrogen molecule.

Vibrational Entanglement through the Lens of Quantum Information Measures.

We introduce a quantum information analysis of vibrational wave functions to understand complex vibrational spectra of molecules with strong anharmonic couplings and vibrational resonances. For this purpose, we define one- and two-modal entropies to guide the identification of strongly coupled vibrational modes and to characterize correlations within modal basis sets. We evaluate these descriptors for multiconfigurational vibrational wave functions which we calculate with the -mode vibrational density matrix renormalization group algorithm.

The odd-number cyclo[13]carbon and its dimer, cyclo[26]carbon.

Molecular rings of carbon atoms (cyclo[]carbons, or C) are excellent benchmarking systems for testing quantum chemical theoretical methods and valuable precursors to other carbon-rich materials. Odd- cyclocarbons, which have been elusive to date, are predicted to be even less stable than even- cyclocarbons. We report the on-surface synthesis of cyclo[13]carbon, C, by manipulation of decachlorofluorene with a scanning probe microscope tip.

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 -mode second-quantization formalism. The resulting -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.

Symmetry-Projected Nuclear-Electronic Hartree-Fock: Eliminating Rotational Energy Contamination.

We present a symmetry projection technique for enforcing rotational and parity symmetries in nuclear-electronic Hartree-Fock wave functions, which treat electrons and nuclei on equal footing. The molecular Hamiltonian obeys rotational and parity inversion symmetries, which are, however, broken by expanding in Gaussian basis sets that are fixed in space. We generate a trial wave function with the correct symmetry properties by projecting the wave function onto representations of the three-dimensional rotation group, i.

The OpenMolcas : 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.

Quantum Computing for Molecular Biology.

Molecular biology and biochemistry interpret microscopic processes in the living world in terms of molecular structures and their interactions, which are quantum mechanical by their very nature. Whereas the theoretical foundations of these interactions are well established, the computational solution of the relevant quantum mechanical equations is very hard. However, much of molecular function in biology can be understood in terms of classical mechanics, where the interactions of electrons and nuclei have been mapped onto effective classical surrogate potentials that model the interaction of atoms or even larger entities.

Second-Order Self-Consistent Field Algorithms: From Classical to Quantum Nuclei.

This work presents a general framework for deriving exact and approximate Newton self-consistent field (SCF) orbital optimization algorithms by leveraging concepts borrowed from differential geometry. Within this framework, we extend the augmented Roothaan-Hall (ARH) algorithm to unrestricted electronic and nuclear-electronic calculations. We demonstrate that ARH yields an excellent compromise between stability and computational cost for SCF problems that are hard to converge with conventional first-order optimization strategies.

Explicitly Correlated Electronic Structure Calculations with Transcorrelated Matrix Product Operators.

In this work, we present the first implementation of the transcorrelated electronic Hamiltonian in an optimization procedure for matrix product states by the density matrix renormalization group (DMRG) algorithm. In the transcorrelation ansatz, the electronic Hamiltonian is similarity-transformed with a Jastrow factor to describe the cusp in the wave function at electron-electron coalescence. As a result, the wave function is easier to approximate accurately with the conventional expansion in terms of one-particle basis functions and Slater determinants.

Quantum Proton Effects from Density Matrix Renormalization Group Calculations.

We recently introduced [ 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.

Expansive Quantum Mechanical Exploration of Chemical Reaction Paths.

Quantum mechanical methods have been well-established for the elucidation of reaction paths of chemical processes and for the explicit dynamics of molecular systems. While they are usually deployed in routine manual calculations on reactions for which some insights are already available (typically from experiment), new algorithms and continuously increasing capabilities of modern computer hardware allow for exploratory open-ended computational campaigns that are unbiased and therefore enable unexpected discoveries. Highly efficient and even automated procedures facilitate systematic approaches toward the exploration of uncharted territory in molecular transformations and dynamics.

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.

Simplified State Interaction for Matrix Product State Wave Functions.

We present an approximation to the state-interaction approach for matrix product state (MPS) wave functions (MPSSI) in a nonorthogonal molecular orbital basis, first presented by Knecht et al. [ , 5881], that allows for a significant reduction of the computational cost without significantly compromising its accuracy. The approximation is well-suited if the molecular orbital basis is close to orthogonality, and its reliability may be estimated a priori with a single numerical parameter.

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.

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.

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 .

Interplay of Stereoelectronic and Vibrational Modulation Effects in Tuning the UPS Spectra of Unsaturated Hydrocarbon Cage Compounds.

The UPS spectra of six hydrocarbon cage compounds have been investigated by a Green-function approach in conjunction with a full harmonic treatment of vibrational modulation effects. The remarkable agreement with experimental results points out the reliability of the proposed computational approach and the strong interplay of stereoelectronic and vibrational effects in tuning the overall spectra.

Modern quantum chemistry with [Open]Molcas.

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree-Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.

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.

The density matrix renormalization group in chemistry and molecular physics: Recent developments and new challenges.

In the past two decades, the density matrix renormalization group (DMRG) has emerged as an innovative new method in quantum chemistry relying on a theoretical framework very different from that of traditional electronic structure approaches. The development of the quantum chemical DMRG has been remarkably fast: it has already become one of the reference approaches for large-scale multiconfigurational calculations. This perspective discusses the major features of DMRG, highlighting its strengths and weaknesses also in comparison with other novel approaches.

Approximate Analytical Gradients and Nonadiabatic Couplings for the State-Average Density Matrix Renormalization Group Self-Consistent-Field Method.

We present an approximate scheme for analytical gradients and nonadiabatic couplings for calculating state-average density matrix renormalization group self-consistent-field wave function. Our formalism follows closely the state-average complete active space self-consistent-field (SA-CASSCF) , which employs a Lagrangian, and the corresponding Lagrange multipliers are obtained from a solution of the coupled-perturbed CASSCF (CP-CASSCF) equations. We introduce a definition of the matrix product state (MPS) Lagrange multipliers based on a single-site tensor in a mixed-canonical form of the MPS, such that a sweep procedure is avoided in the solution of the CP-CASSCF equations.

Theory meets experiment for unravelling the C1s X-ray photoelectron spectra of pyridine, 2-fluoropyridine, and 2,6-difluoropyridine.

High resolution X-ray photoelectron spectra of a series of substituted pyridines (pyridine, 2-fluoropyridine, and 2,6-difluoropyridine) have been recorded and rationalized by means of a quantum mechanical approach based on the density functional theory including vibronic effects at the Franck-Condon level. The significant chemical shifts of the C1s binding energies induced by fluorine atoms are reproduced quantitatively by our computational model, as well as the vibrational fine structure and the band shapes. Nonsymmetric normal modes play an important role due to the core-hole localization in the presence of equivalent carbon atoms in pyridine and 2,6-difluoropyridine.

Large-Scale Quantum Dynamics with Matrix Product States.

Dynamical electronic- and vibrational-structure theories have received a growing interest in the past few years due to their ability to simulate spectra recorded with ultrafast experimental techniques. The exact time evolution of a molecular system can, in principle, be obtained from the time-dependent version of full configuration interaction. Such an approach is, however, limited to few-atom systems due to the exponential increase of its cost with the system dimension.

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

An efficient approximation to the full configuration interaction solution can be obtained with the density matrix renormalization group (DMRG) algorithm without a restriction to a predefined excitation level. In a standard DMRG implementation, however, excited states are calculated with a ground-state optimization in the space orthogonal to all lower lying wave function solutions. A trivial parallelization is therefore not possible, and the calculation of highly excited states becomes prohibitively expensive, especially in regions with a high density of states.

Time-Dependent Formulation of Resonance Raman Optical Activity Spectroscopy.

In this work, we extend the theoretical framework recently developed for the simulation of resonance Raman (RR) spectra of medium-to-large sized systems to its chiral counterpart, namely, resonance Raman optical activity (RROA). The theory is based on a time-dependent (TD) formulation, with the transition tensors obtained as half-Fourier transforms of the appropriate cross-correlation functions. The implementation has been kept as general as possible, supporting adiabatic and vertical models for the PES representation, both in Cartesian and internal coordinates, with the possible inclusion of Herzberg-Teller (HT) effects.