Quantum Chemical Applications Research Articles

Tuned layered double hydroxide-based catalysts inducing singlet oxygen evolution: Reactive oxygen species evolution mechanism exploration, norfloxacin degradation and catalysts screen based on machine learning

Heterogeneous Electro-Fenton (Hetero-EF) is pursued as a booming technique for real effluent treatments. Here, a novel Simulation-Experiment-Prediction framework has been conceived to reveal atomic-level ROS evolution mechanisms and screen robust catalysts for Hetero-EF reaction. Dynamics-static calculations unveil that the spontaneous conversion pathway of 1O2 generation, which elucidate the transformation from *O2, *HO2, *O2− to *1O2, and confirm optimal Co3Fe2-LDH catalysts with lowest rate-determining step and energy barrier (0.82 eV) of 1O2 generation providing a reasonable and feasible strategy fabricating and tuning Hetero-EF catalysts for realistic contaminants. Subsequently, the Co3Fe2-LDH was verified with best catalytic performance and 1O2 yield, which increased by exceeding 20% on NFXN decay and the treatment time was shortened ∼2 h. Machine learning establishes a Gradient Boosting Regressor model for exploring appreciable catalysts for environmental remediation. Our results inspired the exploration and application of quantum chemistry and artificial intelligence for environmental catalysis, remediation and real wastewater purification.

Long-range electron-electron interactions in quantum dot systems and applications in quantum chemistry

Long-range interactions play a key role in several phenomena of quantum physics and chemistry. To study these phenomena, analog quantum simulators provide an appealing alternative to classical numerical methods. Gate-defined quantum dots have been established as a platform for quantum simulation, but for those experiments the effect of long-range interactions between the electrons did not play a crucial role. Here we present the first detailed experimental characterization of long-range electron-electron interactions in an array of gate-defined semiconductor quantum dots. We demonstrate significant interaction strength among electrons that are separated by up to four sites, and show that our theoretical prediction of the screening effects matches well the experimental results. Based on these findings, we investigate how long-range interactions in quantum-dot arrays may be utilized for analog simulations of artificial quantum matter. We numerically show that about ten quantum dots are sufficient to observe binding for a one-dimensional $H_2$-like molecule. These combined experimental and theoretical results pave the way for future quantum simulations with quantum dot arrays and benchmarks of numerical methods in quantum chemistry.

The Tutte polynomial of phenylene systems with given number of branching hexagons

AbstractPolynomial graph invariants have been confirmed to have important applications in quantum chemistry and biological information. One of the famous polynomial graph invariants is the Tutte polynomial which gives abundant graph‐theoretical information of the underlying graph. In this paper, we first give a simpler and more efficient method to get the Tutte polynomials of alternating polycyclic chains. Then we obtain the explicit expressions for the Tutte polynomials and the number of spanning trees of phenylene systems with given number of branching hexagons. Moreover, we determine the extremal values of the number of spanning trees among the phenylene systems with given one branching hexagon and two branching hexagons. The corresponding extremal phenylene systems are characterized, respectively.

Quantum chemistry calculations using energy derivatives on quantum computers

Quantum chemistry calculations such as the prediction of molecular properties and modeling of chemical reactions are a few of the critical areas where near-term quantum computers can showcase quantum advantage. We present a method to calculate energy derivatives for both ground state and excited state energies with respect to the parameters of a chemical system based on the framework of the variational quantum eigensolver (VQE). A low-depth implementation of quantum circuits within the hybrid variational paradigm is designed, and their computational costs are analyzed. We showcase the effectiveness of our method by incorporating it in some key quantum chemistry applications of energy derivatives, such as to perform minimum energy configuration search and estimate molecular response properties estimation of H$_2$ molecule, and also to find the transition state of H$_2$ + H $\leftrightarrow$ H + H$_2$ reaction. The obtained results are shown to be in complete agreement with their respective full configuration interaction (FCI) values.

On the fly swapping algorithm for ordering of degrees of freedom in density matrix renormalization group

Density matrix renormalization group (DMRG) and its time-dependent variants have found widespread applications in quantum chemistry, including ab initio electronic structure of complex bio-molecules, spectroscopy for molecular aggregates, and charge transport in bulk organic semiconductors. The underlying wavefunction ansatz for DMRG, matrix product state (MPS), requires mapping degrees of freedom (DOF) into a one-dimensional topology. DOF ordering becomes a crucial factor for DMRG accuracy. In this work, we propose swapping neighboring DOFs during the DMRG sweeps for DOF ordering, which we term ‘on the fly swapping’ (OFS) algorithm. We show that OFS is universal for both static and time-dependent DMRG with minimum computational overhead. Examples are given for one dimensional antiferromagnetic Heisenberg model, ab initio electronic structure of N2 molecule, and the S1/S2 internal conversion dynamics of pyrazine molecule. It is found that OFS can indeed improve accuracy by finding better DOF ordering in all cases.

EPR Spectroscopy of Cu(II) Complexes: Prediction of g-Tensors Using Double-Hybrid Density Functional Theory

Computational electron paramagnetic resonance (EPR) spectroscopy is an important field of applied quantum chemistry that contributes greatly to connecting spectroscopic observations with the fundamental description of electronic structure for open-shell molecules. However, not all EPR parameters can be predicted accurately and reliably for all chemical systems. Among transition metal ions, Cu(II) centers in inorganic chemistry and biology, and their associated EPR properties such as hyperfine coupling and g-tensors, pose exceptional difficulties for all levels of quantum chemistry. In the present work, we approach the problem of Cu(II) g-tensor calculations using double-hybrid density functional theory (DHDFT). Using a reference set of 18 structurally and spectroscopically characterized Cu(II) complexes, we evaluate a wide range of modern double-hybrid density functionals (DHDFs) that have not been applied previously to this problem. Our results suggest that the current generation of DHDFs consistently and systematically outperform other computational approaches. The B2GP-PLYP and PBE0-DH functionals are singled out as the best DHDFs on average for the prediction of Cu(II) g-tensors. The performance of the different functionals is discussed and suggestions are made for practical applications and future methodological developments.

Optimization of the variational quantum eigensolver for quantum chemistry applications

This work studies the variational quantum eigensolver (VQE) algorithm, which is designed to determine the ground state of a quantum mechanical system by combining classical and quantum hardware. Two methods of reducing the number of required qubit manipulations, prone to induce errors, for the variational quantum eigensolver are studied. First, we formally justify the multiple ℤ2 symmetry qubit reduction scheme first sketched by Bravyi et al. [arXiv:1701.08213 (2017)]. Second, we show that even in small, but non-trivial systems such as H2, LiH, and H2O, the choice of entangling methods (gate based or native) gives rise to varying rates of convergence to the ground state of the system. Through both the removal of qubits and the choice of entangler, the demands on the quantum hardware can be reduced. We find that in general, analyzing the VQE problem is complex, where the number of qubits, the method of entangling, and the depth of the search space all interact. In specific cases however, concrete results can be shown, and an entangling method can be recommended over others as it outperforms in terms of difference from the ground state energy.

Ground states of Heisenberg spin clusters from projected Hartree-Fock theory

We apply the projected Hartree-Fock theory (PHF) for approximating ground states of Heisenberg spin clusters. Spin-rotational, point-group, and complex-conjugation symmetry are variationally restored from a broken-symmetry mean-field reference, where the latter corresponds to a product of local spin states. A fermionic formulation of the Heisenberg model furnishes a conceptual connection to PHF applications in quantum chemistry and detailed equations for a self-consistent field optimization of the reference state are provided. Different PHF variants are benchmarked for ground-state energies and spin-pair correlation functions of antiferromagnetic spin rings and three different polyhedra, with various values of the local spin-quantum number $s$. Although PHF is not suitable to study the thermodynamic limit (where it reduces to the conventional HF results), the low computational cost and the compact wave-function representation make PHF a promising complement to existing approaches for ground states of finite spin clusters, particularly for large local spin $s$ and a moderately large number of sites $N$. The present work may also motivate future explorations of more accurate post-PHF methods for Heisenberg spin clusters.

Considerations for evaluating thermodynamic properties with hybrid quantum-classical computing work flows

Quantum chemistry applications on quantum computers currently rely heavily on the variational quantum eigensolver (VQE) algorithm. This hybrid quantum-classical algorithm aims at finding ground-state solutions of molecular systems based on the variational principle. VQE calculations can be systematically implemented for perturbations to each molecular degree of freedom, generating a Born-Oppenheimer potential-energy surface (PES) for the molecule. The PES can then be used to derive thermodynamic properties, which are often desirable for applications in chemical engineering and materials design. It is clear from this process that quantum chemistry applications contain a substantial classical computing component in addition to steps that can be performed using a quantum computer. In order to design efficient work flows that take full advantage of each hardware type, it is critical to consider the entire process so that the high-accuracy electronic energies possible from quantum computing are not squandered in the process of calculating thermodynamic properties. We present a summary of the hybrid quantum-classical work flow to compute thermodynamic properties. This work flow contains many options that can significantly affect the efficiency and the accuracy of the results, including classical optimizer attributes, number of ansatz repetitions, and how the vibrational Schr\"odinger equation is solved to determine vibrational modes. We also analyze the effects of these options by employing robust statistics along with simulations and experiments on actual quantum hardware. We show that, through careful selection of work flow options, nearly order-of-magnitude increases in accuracy are possible at equivalent computing time.

Factorized Structure of the Long-Range Two-Electron Integrals Tensor and its Application in Quantum Chemistry

Factorized Structure of the Long-Range Two-Electron Integrals Tensor and its Application in Quantum Chemistry

Electronic Structure Theory Calculations Using Modern Architectures: KNL vs Haswell

The time to solution and parallel efficiency of several commonly used electronic structure methods (Hartree-Fock, density functional theory, second order perturbation theory, resolution of the identity second order perturbation theory, coupled cluster) are evaluated on both the Intel Xeon Haswell and the Intel Xeon Phi Knights Landing (KNL) architectures. The Haswell completes the benchmark calculations with a faster time to solution than the KNL for all molecules and methods tested. While the Haswell exhibits an average speedup of at least 3.5 relative to the KNL for all nonthreaded computations, the KNL has a better parallel efficiency than the Haswell with increasing core counts. The architectures are further tested using a more computationally costly coupled cluster method on a transition state reaction. The Haswell appears to be the best choice to minimize the time to solution, though for very large systems and high levels of theory that require memory intensive processes the superior memory hierarchy and larger on node memory of the KNL can make it a better choice. These results are used to showcase aspects of novel architectures that will increase efficiency for quantum chemistry applications.

Solving the Schrödinger equation using program synthesis.

We demonstrate that a program synthesis approach based on a linear code representation can be used to generate algorithms that approximate the ground-state solutions of one-dimensional time-independent Schrödinger equations constructed with bound polynomial potential energy surfaces (PESs). Here, an algorithm is constructed as a linear series of instructions operating on a set of input vectors, matrices, and constants that define the problem characteristics, such as the PES. Discrete optimization is performed using simulated annealing in order to identify sequences of code-lines, operating on the program inputs that can reproduce the expected ground-state wavefunctions ψ(x) for a set of target PESs. The outcome of this optimization is not simply a mathematical function approximating ψ(x) but is, instead, a complete algorithm that converts the input vectors describing the system into a ground-state solution of the Schrödinger equation. These initial results pointthe way toward an alternative route for developing novel algorithms for quantum chemistry applications.

Concentration for Random Product Formulas

Quantum simulation has wide applications in quantum chemistry and physics. Recently, scientists have begun exploring the use of randomized methods for accelerating quantum simulation. Among them, a simple and powerful technique, called qDRIFT, is known to generate random product formulas for which the average quantum channel approximates the ideal evolution. qDRIFT achieves a gate count that does not explicitly depend on the number of terms in the Hamiltonian, which contrasts with Suzuki formulas. This work aims to understand the origin of this speed-up by comprehensively analyzing a single realization of the random product formula produced by qDRIFT. The main results prove that a typical realization of the randomized product formula approximates the ideal unitary evolution up to a small diamond-norm error. The gate complexity is already independent of the number of terms in the Hamiltonian, but it depends on the system size and the sum of the interaction strengths in the Hamiltonian. Remarkably, the same random evolution starting from an arbitrary, but fixed, input state yields a much shorter circuit suitable for that input state. In contrast, in deterministic settings, such an improvement usually requires initial state knowledge. The proofs depend on concentration inequalities for vector and matrix martingales, and the framework is applicable to other randomized product formulas. Our bounds are saturated by certain commuting Hamiltonians.

Effect of the deformation parameter on the nonrelativistic energy spectra of the q-deformed Hulthen-quadratic exponential-type potential

In this study, an approximate solution of the Schr�dinger equation for the q-deformed Hulthen-quadratic exponential-type potential model within the framework of the Nikiforov�Uvarov method was obtained. The bound state energy equation and the corresponding eigenfunction was obtained. The energy spectrum is applied to study H2, HCl, CO and LiH diatomic molecules. The effect of the deformation parameters and other potential parameters on the energy spectra of the system were graphically and numerically analyzed in detail. Special cases were considered when the potential parameters were altered, resulting in deformed Hulthen potential, Hulthen potential, deformed quadratic exponential-type potential and quadratic exponential-type potential. The energy eigenvalues expressions agreed with what obtained in literature. Finally, the results can find many applications in quantum chemistry, atomic and molecular physics.

Lower Bounds for Nonrelativistic Atomic Energies.

A recently developed lower bound theory for Coulombic problems (E. Pollak, R. Martinazzo, J. Chem. Theory Comput.2021, 17, 1535) is further developed and applied to the highly accurate calculation of the ground-state energy of two- (He, Li+, and H–) and three- (Li) electron atoms. The method has been implemented with explicitly correlated many-particle basis sets of Gaussian type, on the basis of the highly accurate (Ritz) upper bounds they can provide with relatively small numbers of functions. The use of explicitly correlated Gaussians is developed further for computing the variances, and the necessary modifications are here discussed. The computed lower bounds are of submilli-Hartree (parts per million relative) precision and for Li represent the best lower bounds ever obtained. Although not yet as accurate as the corresponding (Ritz) upper bounds, the computed bounds are orders of magnitude tighter than those obtained with other lower bound methods, thereby demonstrating that the proposed method is viable for lower bound calculations in quantum chemistry applications. Among several aspects, the optimization of the wave function is shown to play a key role for both the optimal solution of the lower bound problem and the internal check of the theory.

Magnetic exchange coupling in Cu dimers studied with modern multireference methods and broken-symmetry coupled cluster theory

Spin-state energetics of exchange-coupled copper complexes pose a persistent challenge for applied quantum chemistry. Here, we provide a comprehensive comparison of all available theoretical approaches to the problem of exchange coupling in two antiferromagnetically coupled bis-μ-hydroxo Cu(II) dimers. The evaluated methods include multireference methods based on the density matrix renormalization group (DMRG), multireference methods that incorporate dynamic electron correlation either perturbatively, such as the N-electron valence state perturbation theory, or variationally, such as the difference-dedicated configuration interaction. In addition, we contrast the multireference results with those obtained using broken-symmetry approaches that utilize either density functional theory or, as demonstrated here for the first time in such systems, a local implementation of coupled cluster theory. The results show that the spin-state energetics of these copper dimers are dominated by dynamic electron correlation and represent an impossible challenge for multireference methods that rely on brute-force expansion of the active space to recover correlation energy. Therefore, DMRG-based methods even at the limit of their applicability cannot describe quantitatively the antiferromagnetic exchange coupling in these dimers, in contrast to dinuclear complexes of earlier transition metal ions. The convergence of the broken-symmetry coupled cluster approach is studied and shown to be a limiting factor for the practical application of the method. The advantages and disadvantages of all approaches are discussed, and recommendations are made for future developments.

Symmetric and antisymmetric kernels for machine learning problems in quantum physics and chemistry

We derive symmetric and antisymmetric kernels by symmetrizing and antisymmetrizing conventional kernels and analyze their properties. In particular, we compute the feature space dimensions of the resulting polynomial kernels, prove that the reproducing kernel Hilbert spaces induced by symmetric and antisymmetric Gaussian kernels are dense in the space of symmetric and antisymmetric functions, and propose a Slater determinant representation of the antisymmetric Gaussian kernel, which allows for an efficient evaluation even if the state space is high-dimensional. Furthermore, we show that by exploiting symmetries or antisymmetries the size of the training data set can be significantly reduced. The results are illustrated with guiding examples and simple quantum physics and chemistry applications.

Spectroscopic and Quantum Chemical Investigations of Hypothetical m-Diacetylaminoazopyrimidine and Its Photovoltaic Properties

Abstract: Here, an attempt is made to theoretically study and predict the electronic and spectroscopic (UV-Vis and IR) and structural properties, quantum chemical descriptors and subsequent application of diacetylaminoazopyrimidine in dye-sensitized solar cells (DSSCs). Ground- and excited-state time-dependent density functional theory (TD-DFT) calculations were carried out using material studio and ORCA software, respectively. The computed ground-state energy gap, chemical hardness, chemical softness, chemical potential, electronegativity and electrophilicity index are: 3.60 eV, 1.80 eV, 0.56 eV, 4.49 eV, -4.49 eV and 5.68, respectively. Conversely, the DFT-predicted excited-state quantum chemical descriptors are: 1.67 eV, 0.83 eV, 1.20 eV, 4.71 eV and -4.71 eV, corresponding to the energy gap, chemical hardness, chemical softness, chemical potential and electronegativity, respectively. Furthermore, vibrational frequency calculations confirm the presence of some key functional groups (N=N, C=O, C-H) present in the dye molecules. The computed optoelectronic parameters, such as light-harvesting efficiency, electron injection and open-circuit voltage are 0.06 eV, -8.59 eV and -5.75 eV, respectively. Overall, the dye possesses a relatively good current conversion efficiency as compared to other dyes studied in the literature; hence, it could be used as a novel material for photovoltaic technological applications. Keywords: Diacetylaminoazopyrimidine, DFT, Excited state, Spectroscopy, DSSCs.

SEQCROW: A ChimeraX bundle to facilitate quantum chemical applications to complex molecular systems.

We describe a bundle for UCSF ChimeraX called SEQCROW that provides advanced structure editing capabilities and quantum chemistry utilities designed for complex organic and organometallic compounds. SEQCROW includes graphical presets and bond editing tools that facilitate the generation of publication-quality molecular structure figures while also allowing users to build molecular structures quickly and efficiently by mapping new ligands onto existing organometallic complexes as well as adding rings and substituents. Other capabilities include the ability to visualize vibrational modes and simulated IR spectra, to compute and visualize molecular descriptors including percent buried volume, ligand cone angles, and Sterimol parameters, to process thermochemical corrections from quantum mechanical computations, to generate input files for ORCA, Psi4, and Gaussian, and to run and manage computational jobs.

Application of topological quantum chemistry in electrides

The recently developed theory of topological quantum chemistry (TQC) has built a close connection between band representations in momentum space and orbital characters in real space. It provides an effective way to diagnose topological materials, leading to the discovery of lots of topological materials after the screening of all known nonmagnetic compounds. On the other hand, it can also efficiently reveal spacial orbital characters, including average charge centers and site-symmetry characters. By using TQC theory with the computed irreducible representations in the first-principles calculations, we demonstrate that the electrides with excess electrons serving as anions at vacancies can be well identified by analyzing band representations (BRs), which cannot be expressed as a sum of atomic-orbital-induced band representations (aBRs). In fact, the floating bands (formed by the excess electrons) belong to the BRs induced from the "pseudo-orbitals" centered at vacancies. In other words, the electrides are proved to be unconventional ionic crystals, where a set of occupied bands is not a sum of aBRs but necessarily contains a BR from vacancies. The TQC theory provides a promising avenue to pursue more electride candidates in ionic crystals.