Search

We present an implementation of a perturbative triples correction for coupled cluster including single and double excitations based on the transcorrelated Hamiltonian. Transcorrelation introduces explicit electron correlation in the electronic Hamiltonian through similarity transformation with a correlation factor. Due to this transformation, the transcorrelated Hamiltonian includes up to three-body couplings and becomes non-hermitian. Since conventional coupled cluster is solved by projection, it is well suited to harbor non-hermitian Hamiltonians. The arising three-body operator, however, creates a huge memory bottleneck and increases the runtime scaling of the coupled cluster equations. As it has been shown that the three-body operator can be approximated, by expressing the Hamiltonian in the normal-ordered form, we investigate this approximation for the perturbative triples correction. Results are compared with code-generation based transcorrelated coupled cluster up to quadruple excitations., Comment: 31 pages, 4 figures, 3 tables

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 multi-configurational vibrational wave functions which we calculate with the n-mode vibrational density matrix renormalization group algorithm. Based on the quantum information measures, we present a vibrational entanglement analysis of the vibrational ground and excited states of CO2, which display strong anharmonic effects due to the symmetry-induced and accidental (near-) degeneracies. We investigate the entanglement signature of the Fermi resonance and discuss the maximally entangled state arising from the two degenerate bending modes., Comment: 23 pages, 3 figures, 1 ToC graphic

We provide a taxonomy of concurrency models for BDI frameworks, elicited by analysing state-of-the-art technologies, and aimed at helping both BDI designers and developers in making informed decisions. Comparison among BDI technologies w.r.t. concurrency models reveals heterogeneous support, and low customisability.

The execution of Belief-Desire-Intention (BDI) agents in a Multi-Agent System (MAS) can be practically implemented on top of low-level concurrency mechanisms that impact on efficiency, determinism, and reproducibility. We argue that developers should specify the MAS behaviour independently of the execution model, and choose or configure the concurrency model later on, according to the specific needs of their target domain, leaving the MAS specification unaffected. We identify patterns for mapping the agent execution over the underlying concurrency abstractions, and investigate which concurrency models are supported by some of the most commonly used BDI platforms. Although most frameworks support multiple concurrency models, we find that they mostly hide them under the hood, making them opaque to the developer, and actually limiting the possibility of fine-tuning the MAS.

Hybrid Tensor Networks (hTN) offer a promising solution for encoding variational quantum states beyond the capabilities of efficient classical methods or noisy quantum computers alone. However, their practical usefulness and many operational aspects of hTN-based algorithms, like the optimization of hTNs, the generalization of standard contraction rules to an hybrid setting, and the design of application-oriented architectures have not been thoroughly investigated yet. In this work, we introduce a novel algorithm to perform ground state optimizations with hybrid Tree Tensor Networks (hTTNs), discussing its advantages and roadblocks, and identifying a set of promising applications. We benchmark our approach on two paradigmatic models, namely the Ising model at the critical point and the Toric code Hamiltonian. In both cases, we successfully demonstrate that hTTNs can improve upon classical equivalents with equal bond dimension in the classical part., Comment: 19 pages, 16 figures

We introduce a novel computational framework for excited-states 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. For this system, we critically compare the accuracy and efficiency of different quantum subspace expansion and equation-of-motion algorithms and show that only methods that can capture both weak and strong electron correlation effects can properly describe the non-adiabatic effects that tune the reactive event.

Quantum computation is one of the most promising new paradigms for the simulation of physical systems composed of electrons and atomic nuclei, with applications in chemistry, solid-state physics, materials science, and molecular biology. This requires a truncated representation of the electronic structure Hamiltonian using a finite number of orbitals. While it is, in principle, obvious how to improve on the representation by including more orbitals, this is usually unfeasible in practice (e.g., because of the limited number of qubits available) and severely compromises the accuracy of the obtained results. Here, we propose a quantum algorithm that improves on the representation of the physical problem by virtue of second-order perturbation theory. In particular, our quantum algorithm evaluates the second-order energy correction through a series of time-evolution steps under the unperturbed Hamiltonian. An important application is to go beyond the active-space approximation, allowing to include corrections of virtual orbitals, known as multireference perturbation theory. Here, we exploit that the unperturbed Hamiltonian is diagonal for virtual orbitals and show that the number of qubits is independent of the number of virtual orbitals. This gives rise to more accurate energy estimates without increasing the number of qubits. Moreover, we demonstrate numerically for realistic chemical systems that the total runtime has highly favorable scaling in the number of virtual orbitals compared to previous work. Numerical calculations confirm the necessity of the multireference perturbation theory energy corrections to reach accurate ground state energy estimates. Our perturbation theory quantum algorithm can also be applied to symmetry-adapted perturbation theory. As such, we demonstrate that perturbation theory can help to reduce the quantum hardware requirements for quantum chemistry., Comment: 25 pages, 6 figures

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., Comment: 48 pages, 6 figures, 2 tables

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.e., the special orthogonal group in three dimensions SO(3). As a consequence, the wave function becomes an eigenfunction of the angular momentum operator which (i) eliminates the contamination of the ground state wave function by highly excited rotational states arising from the broken rotational symmetry, and (ii) enables the targeting of specific rotational states of the molecule. We demonstrate the efficiency of the symmetry projection technique by calculating energies of the low-lying rotational states of the H$_2$ and H$_3^+$ molecules., Comment: 37 pages, 3 tables

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 very 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. The simple mathematical structure of these potentials offers huge computational advantages; however, this comes at the cost that all quantum correlations and the rigorous many-particle nature of the interactions are omitted. In this work, we discuss how quantum computation may advance the practical usefulness of the quantum foundations of molecular biology by offering computational advantages for simulations of biomolecules. We not only discuss typical quantum mechanical problems of the electronic structure of biomolecules in this context, but also consider the dominating classical problems (such as protein folding and drug design) as well as data-driven approaches of bioinformatics and the degree to which they might become amenable to quantum simulation and quantum computation., Comment: 76 pages, 7 figures

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. In the electronic case, we show that ARH overcomes the slow convergence of orbitals in strongly-correlated molecules with the example of several iron-sulfur clusters. For nuclear-electronic calculations, ARH significantly enhances the convergence already for small molecules, as demonstrated for a series of protonated water clusters., Comment: 40 pages, 4 figures, 3 tables

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. The transcorrelated Hamiltonian in first quantization comprises up to three-body interactions, which we deal with in the standard way by applying robust density fitting to two- and three-body integrals entering the second-quantized representation of this Hamiltonian. The lack of hermiticity of the transcorrelated Hamiltonian is taken care of along the lines of the first work on transcorrelated DMRG [J. Chem. Phys. 153, 164115 (2020)] by encoding it as a matrix product operator and optimizing the corresponding ground state wave function with imaginary-time time-dependent DMRG. We demonstrate our quantum chemical transcorrelated DMRG approach at the example of several atoms and first-row diatomic molecules. We show that transcorrelation improves the convergence rate to the complete basis set limit in comparison to conventional DMRG. Moreover, we study extensions of our approach that aim at reducing the cost of handling the matrix product operator representation of the transcorrelated Hamiltonian., Comment: 52 pages, 6 figures, 4 tables

Quantum mechanical methods have been devised for the elucidation and clarification of reaction paths of chemical processes over decades. While they are typically deployed in routine calculations on systems for which some insights have already been gained, new algorithms and ever increasing computer hardware capabilities allow now for exploratory open-ended computational campaigns that bear the option for unexpected discoveries and the systematic approach toward uncharted territory of molecular transformations and dynamics. In this work, we discuss the current state of the art of such exploratory techniques., Comment: 24 pages, 6 figures

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., Comment: 44 pages, 9 figures, 4 tables

This review elaborates on the foundation, the advantages, and the prospects of tensor network representations for quantum states in vibrational spectroscopy. The focus is on the recently introduced matrix product state decomposition of nuclear quantum states and its optimization by the density matrix renormalization group algorithm., Comment: 65 pages, 20 figures

We recently introduced [J. Chem. Phys. 152 2020, 204103] the nuclear-electronic all-particle density matrix renormalization group method (NEAP-DMRG) to solve the molecular Schr\"{o}dinger equation, based on a stochastically optimized orbital basis, without invoking the Born-Oppenheimer approximation. In this work, we combine the DMRG with nuclear-electronic Hartree-Fock (NEHF-DMRG), treating nuclei and electrons on the same footing. Inter- and intra-species correlations are described within the DMRG 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 NEHF-DMRG 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 ground-state absolute energy and proton density which is a major challenge for multi-reference nuclear-electronic state-of-the-art methods., Comment: 55 pages, 8 figures, 5 tables

A new and rapidly growing econometric literature is making advances in the problem of using machine learning methods for causal inference questions. Yet, the empirical economics literature has not started to fully exploit the strengths of these modern methods. We revisit influential empirical studies with causal machine learning methods and identify several advantages of using these techniques. We show that these advantages and their implications are empirically relevant and that the use of these methods can improve the credibility of causal analysis.

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 non-relativistic 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., Comment: 47 pages, 11 figures

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 at 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., Comment: 31 pages, 4 figures, 4 tables

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 wave function, 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 the best 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.

We introduce the Nuclear Electronic All-Particle Density Matrix Renormalization Group (NEAP-DMRG) method for solving the time-independent Schr\"odinger equation simultaneously for electrons and other quantum species. In contrast to 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 parametrization of the total wave function as a matrix product state allows NEAP-DMRG to accurately approximate full configuration interaction energies of molecular systems with more than three nuclei and twelve particles in total, which is currently a major challenge for other multicomponent approaches. We present NEAP-DMRG results for two few-body systems, i.e., H$_2$ and H$_3^+$, and one larger system, namely BH$_3$, Comment: 37 pages, 5 figures, 2 tables

Decentralized Online Social Networks (DOSNs) have been proposed as an alternative solution to the current centralized Online Social Networks (OSNs). Online Social Networks are based on centralized architecture (e.g., Facebook, Twitter, or Google+), while DOSNs do not have a service provider that acts as central authority and users have more control over their information. Several DOSNs have been proposed during the last years. However, the decentralization of the OSN requires efficient solutions for protecting the privacy of users, and to evaluate the trust between users. Blockchain represents a disruptive technology which has been applied to several fields, among these also to Social Networks. In this paper, we propose a manageable, user-driven and auditable access control framework for DOSNs using blockchain technology. In the proposed approach, the blockchain is used as a support for the definition of privacy policies. The resource owner uses the public key of the subject to define flexible role-based access control policies, while the private key associated with the subject's Ethereum account is used to decrypt the private data once access permission is validated on the blockchain. We evaluate our solution by exploiting the Rinkeby Ethereum testnet to deploy the smart contract, and to evaluate its performance. Experimental results show the feasibility of the proposed scheme in achieving auditable and user-driven access control via smart contract deployed on the Blockchain.

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 to other novel approaches. The method is presented following its historical development, starting from its original formulation up to its most recent applications. Possible routes to recover dynamical correlation are discussed in detail. Emerging new fields of applications of DMRG are explored, in particular its time-dependent formulation and the application to vibrational spectroscopy., Comment: 25 pages, 2 figures

We present an approximate scheme for analytical gradients and nonadiabatic couplings for calculating state-average density matrix renormalization group self-consistent-field wavefunction. Our formalism follows closely the state-average complete active space self-consistent-field (SA-CASSCF) \emph{ansatz}, 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. We apply our implementation to the optimization of a conical intersection in 1,2-dioxetanone, where we are able to fully reproduce the SA-CASSCF result up to arbitrary accuracy., Comment: 17 pages, 3 figures, 2 tables

Dynamical electronic- and vibrational-structure theories have received a growing interest in the last 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. In the present work, we overcome this unfavorable scaling by employing the time-dependent density matrix renormalization group (TD-DMRG) which parametrizes the time-dependent wavefunction as a matrix product state. The time-dependent Schroedinger equation is then integrated with a sweep-based algorithm, as in standard time-independent DMRG. Unlike other TD-DMRG approaches, the one presented here leads to a set of coupled equations that can be integrated exactly. The resulting theory enables us to study real- and imaginary-time evolutions of Hamiltonians comprising more than 20 degrees of freedom that are challenging for current state-of-the-art quantum dynamics algorithms. We apply our algorithm to the simulation of quantum dynamics of models of increasing complexity, ranging from simple excitonic Hamiltonians to more complex ab-initio vibronic ones., Comment: 60 pages, 14 figures