Source author record

Markus Reiher

Markus Reiher appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

62works
14topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

62 published item(s)

preprint2023arXiv

Universal QM/MM Approaches for General Nanoscale Applications

Quantum mechanics/molecular mechanics (QM/MM) hybrid models allow one to address chemical phenomena in complex molecular environments. Whereas this modeling approach can cope with a large system size at moderate computational costs, the models are often tedious to construct and require manual preprocessing and expertise. As a result, transferability to new application areas can be limited and the many parameters are not easy to adjust to reference data that are typically scarce. Therefore, it is desirable to devise automated procedures of controllable accuracy, which enables such modeling in a standardized and black-box-type manner. Although diverse best-practice protocols have been set up for the construction of individual components of a QM/MM model (e.g., the MM potential, the type of embedding, the choice of the QM region), automated procedures that reconcile all steps of the QM/MM model construction are still rare. Here, we review the state of the art of QM/MM modeling with a focus on automation. We elaborate on MM model parametrization, on atom-economical physically-motivated QM region selection, and on embedding schemes that incorporate mutual polarization as critical components of the QM/MM model. In view of the broad scope of the field, we mostly restrict the discussion to methodologies that build \textit{de novo} models based on first-principles data, on uncertainty quantification, and on error mitigation with a high potential for automation. Ultimately, it is desirable to be able to set up reliable QM/MM models in a fast and efficient automated way without being constrained by specific chemical or technical limitations.

preprint2022arXiv

Autonomous Reaction Network Exploration in Homogeneous and Heterogeneous Catalysis

Autonomous computations that rely on automated reaction network elucidation algorithms may pave the way to make computational catalysis on a par with experimental research in the field. Several advantages of this approach are key to catalysis: (i) Automation allows one to consider orders of magnitude more structures in a systematic and open-ended fashion than what would be accessible by manual inspection. Eventually, full resolution in terms of structural varieties and conformations as well as with respect to the type and number of potentially important elementary reaction steps (including decomposition reactions that determine turnover numbers) may be achieved. (ii) Fast electronic structure methods with uncertainty quantification warrant high efficiency and reliability in order to not only deliver results quickly, but also to allow for predictive work. (iii) A high degree of autonomy reduces the amount of manual human work, processing errors, and human bias. Although being inherently unbiased, it is still steerable with respect to specific regions of an emerging network and with respect to the addition of new reactant species. This allows for a high fidelity of the formalization of some catalytic process and for surprising in silico discoveries. In this work, we first review the state of the art in computational catalysis to embed autonomous explorations into the general field from which it draws its ingredients. We then elaborate on the specific conceptual issues that arise in the context of autonomous computational procedures, some of which we discuss at an example catalytic system.

preprint2022arXiv

Chemoton 2.0: Autonomous Exploration of Chemical Reaction Networks

Fueled by advances in hardware and algorithm design, large-scale automated explorations of chemical reaction space have become possible. Here, we present our approach to an open-source, extensible framework for explorations of chemical reaction mechanisms based on the first principles of quantum mechanics. It is intended to facilitate reaction network explorations for diverse chemical problems with a wide range of goals such as mechanism elucidation, reaction path optimization, retrosynthetic path validation, reagent design, and microkinetic modeling. The stringent first-principles basis of all algorithms in our framework is key for the general applicability that avoids any restrictions to specific chemical systems. Such an agile framework requires multiple specialized software components of which we present three modules in this work. The key module, Chemoton,drives the exploration of reaction networks. For the exploration itself, we introduce two new algorithms for elementary-step searches that are based on Newton trajectories. The performance of these algorithms is assessed for a variety of reactions characterized by a broad chemical diversity in terms of bonding patterns and chemical elements. We reproduce and significantly extend what is known about these reactions and provide the resulting data to be used as a starting point for further explorations and for future reference.

preprint2022arXiv

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. 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.

preprint2022arXiv

Heuristics and Uncertainty Quantification in Rational and Inverse Compound and Catalyst Design

The goal of inverse (quantum) approaches is to devise methods and approaches capable of efficiently searching chemical space in such a way that the design of novel materials and compounds with specific properties is as direct and efficient as possible. Here, we review the current state of the field with a focus on the most recent developments. We discuss the importance of heuristic rules and human intuition for rational compound design. Moreover, we elaborate on options for reliable uncertainty quantification for computational results, which is crucial for a truly predictive application of any in silico method.

preprint2022arXiv

Quantum Proton Effects from Density Matrix Renormalization Group Calculations

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ö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.

preprint2022arXiv

The Transferability Limits of Static Benchmarks

Every practical method to solve the Schrödinger equation for interacting many-particle systems introduces approximations. Such methods are therefore plagued by systematic errors. For computational chemistry, it is decisive to quantify the specific error for some system under consideration. Traditionally, the primary resource for such an error assessment have been benchmarking results, usually taken from the literature. However, their transferability to a specific molecular system, and hence, the reliability of the traditional approach always remains uncertain to some degree. In this communication, we elaborate on the shortcomings of this traditional way of static benchmarking by exploiting statistical analyses at the example of one of the largest quantum chemical benchmark sets available. We demonstrate the uncertainty of error estimates in the light of the choice of reference data selected for a benchmark study. To alleviate the issues with static benchmarks, we advocate to rely instead on a rolling and system-focused approach for rigorously quantifying the uncertainty of a quantum chemical result.

preprint2022arXiv

Ultra-fast Spectroscopy for High-Throughput and Interactive Quantum Chemistry

We present ultra-fast quantum chemical methods for the calculation of infrared and ultraviolet-visible spectra designed to provide fingerprint information during autonomous and interactive explorations of molecular structures. Characteristic spectral signals can serve as diagnostic probes for the identification and characterization of molecular structures. These features often do not require ultimate accuracy with respect to peak position and intensity, which alleviates the accuracy-time dilemma in ultra-fast electronic structure methods. If approximate ultra-fast algorithms are supplemented with an uncertainty quantification scheme for the detection of potentially large prediction errors in signal position and intensity, an offline refinement will always be possible to confirm or discard the predictions of the ultra-fast approach. Here, we present ultra-fast electronic structure methods for such a protocol in order to obtain ground- and excited-state electronic energies, dipole moments, and their derivatives for real-time applications in vibrational spectroscopy and photophysics. As part of this endeavor, we devise an information-inheritance partial Hessian approach for vibrational spectroscopy, a tailored subspace diagonalization approach and a determinant-selection scheme for excited-state calculations.

preprint2021arXiv

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.

preprint2021arXiv

Expansive Quantum Mechanical Exploration of Chemical Reaction Paths

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.

preprint2021arXiv

Molecule-specific Uncertainty Quantification in Quantum Chemical Studies

Solving the electronic Schrödinger equation for changing nuclear coordinates provides access to the Born-Oppenheimer potential energy surface. This surface is the key starting point for almost all theoretical studies of chemical processes in electronic ground and excited states (including molecular structure prediction, reaction mechanism elucidation, molecular property calculations, quantum and molecular dynamics). Electronic structure models aim at a sufficiently accurate approaximation of this surface. They have therefore become a cornerstone of theoretical and computational chemistry, molecular physics, and materials science. In this work, we elaborate on general features of approximate electronic structure models such as accuracy, transferability, and general applicability in order to arrive at a perspective for future developments, of which a vanguard has already arrived. Our quintessential proposition is that meaningful quantum mechanical predictions for chemical phenomena require system-specific uncertainty information for each and every electronic structure calculation, if objective conclusions shall be drawn with confidence.

preprint2021arXiv

Prospects of Quantum Computing for Molecular Sciences

Molecular science is governed by the dynamics of electrons, atomic nuclei, and their interaction with electromagnetic fields. A reliable physicochemical understanding of these processes is crucial for the design and synthesis of chemicals and materials of economic value. Although some problems in this field are adequately addressed by classical mechanics, many require an explicit quantum mechanical description. Such quantum problems represented by exponentially large wave function should naturally benefit from quantum computation on a number of logical qubits that scales only linearly with system size. In this perspective, we focus on the potential of quantum computing for solving relevant problems in the molecular sciences -- molecular physics, chemistry, biochemistry, and materials science.

preprint2021arXiv

Solvation Free Energies in Subsystem Density Functional Theory

For many chemical processes the accurate description of solvent effects are vitally important. Here, we describe a hybrid ansatz for the explicit quantum mechanical description of solute-solvent and solvent-solvent interactions based on subsystem density functional theory and continuum solvation schemes. Since explicit solvent molecules may compromise scalability of the model and transferability of the predicted solvent effect, we aim to retain both, for different solutes as well as for different solvents. The key for the transferability is the consistent subsystem decomposition of solute and solvent. The key for the scalability is the performance of subsystem DFT for increasing numbers of subsystems. We investigate molecular dynamics and stationary point sampling of solvent configurations and compare the resulting (Gibbs) free energies to experiment and theoretical methods. We can show that with our hybrid model reaction barriers and reaction energies are accurately reproduced compared to experimental data.

preprint2020arXiv

Analytically projected rotationally symmetric explicitly correlated Gaussian Functions with one-axis-shifted centers

A new explicitly correlated functional form for expanding the wave function of an N-particle system with arbitrary angular momentum and parity is presented. We develop the projection-based approach, numerically exploited in our previous work [J. Chem. Phys. 149, 184105 (2018)], to explicitly correlated Gausssians with one-axis shifted centers and derive the matrix elements for the Hamiltonian and the angular momentum operators by analytically solving the integral projection operator. Variational few-body calculations without assuming the Born-Oppenheimer approximation are presented for several rotationally excited states of three- and four-particle systems. We show how the new formalism can be used as a unified framework for high-accuracy calculations of properties of small atoms and molecules.

preprint2020arXiv

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 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.

preprint2020arXiv

Immersive Interactive Quantum Mechanics for Teaching and Learning Chemistry

The impossibility of experiencing the molecular world with our senses hampers teaching and understanding chemistry because very abstract concepts (such as atoms, chemical bonds, molecular structure, reactivity) are required for this process. Virtual reality, especially when based on explicit physical modeling (potentially in real time), offers a solution to this dilemma. Chemistry teaching can make use of advanced technologies such as virtual-reality frameworks and haptic devices. We show how an immersive learning setting could be applied to help students understand the core concepts of typical chemical reactions by offering a much more intuitive approach than traditional learning settings. Our setting relies on an interactive exploration and manipulation of a chemical system; this system is simulated in real-time with quantum chemical methods, and therefore, behaves in a physically meaningful way.

preprint2020arXiv

Molassembler: Molecular graph construction, modification and conformer generation for inorganic and organic molecules

We present the graph-based molecule software Molassembler for building organic and inorganic molecules. Molassembler provides algorithms for the construction of molecules built from any set of elements from the periodic table. In particular, poly-nuclear transition metal complexes and clusters can be considered. Structural information is encoded as a graph. Stereocenter configurations are interpretable from Cartesian coordinates into an abstract index of permutation for an extensible set of polyhedral shapes. Substituents are distinguished through a ranking algorithm. Graph and stereocenter representations are freely modifiable and chiral state is propagated where possible through incurred ranking changes. Conformers are generated with full stereoisomer control by four spatial dimension Distance Geometry with a refinement error function including dihedral terms. Molecules are comparable by an extended graph isomorphism and their representation is canonicalizeable. Molassembler is written in C++ and provides Python bindings.

preprint2020arXiv

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 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$

preprint2020arXiv

Self-Parametrizing System-Focused Atomistic Models

Computational studies of chemical reactions in complex environments such as proteins, nanostructures, or on surfaces require accurate and efficient atomistic models applicable to the nanometer scale. In general, an accurate parametrization of the atomistic entities will not be available for arbitrary system classes, but demands a fast automated system-focused parametrization procedure to be quickly applicable, reliable, flexible, and reproducible. Here, we develop and combine an automatically parametrizable quantum chemically derived molecular mechanics model with machine-learned corrections under autonomous uncertainty quantification and refinement. Our approach first generates an accurate, physically motivated model from a minimum energy structure and its corresponding Hessian matrix by a partial Hessian fitting procedure of the force constants. This model is then the starting point to generate a large number of configurations for which additional off-minimum reference data can be evaluated on the fly. A $Δ$-machine learning model is trained on these data to provide a correction to energies and forces including uncertainty estimates. During the procedure, the flexibility of the machine learning model is tailored to the amount of available training data. The parametrization of large systems is enabled by a fragmentation approach. Due to their modular nature, all model construction steps allow for model improvement in a rolling fashion. Our approach may also be employed for the generation of system-focused electrostatic molecular mechanics embedding environments in a quantum-mechanical/molecular-mechanical hybrid model for arbitrary atomistic structures at the nanoscale.

preprint2020arXiv

Semiclassical Dispersion Corrections efficiently improve Multi-Configurational Theory with Short-Range Density-Functional Dynamic Correlation

Multi-configurational wave functions are known to describe electronic structure across a Born-Oppenheimer surface qualitatively correct. However, for quantitative reaction energies, dynamical correlation originating from the many configurations involving excitations out of the restricted orbital space, the active space, must be considered. Standard procedures involve approximations that eventually limit the ultimate accuracy achievable (most prominently, multi-reference perturbation theory). At the same time, the computational cost increase dramatically due to the necessity to obtain higher-order reduced density matrices. It is this disproportion that leads us here to propose a MC-srDFT-D hybrid approach of semiclassical dispersion (D) corrections to cover long-range dynamical correlation in a multi-configurational (MC) wave function theory which includes short-range (sr) dynamical correlation by density functional theory (DFT) without double counting. We demonstrate that the reliability of this approach is very good (at negligible cost), especially when considering that standard second-order multi-reference perturbation theory usually overestimates dispersion interactions.

preprint2020arXiv

Systematic Microsolvation Approach with a Cluster-Continuum Scheme and Conformational Sampling

Solvation is a notoriously difficult and nagging problem for the rigorous theoretical description of chemistry in the liquid phase. Successes and failures of various approaches ranging from implicit solvation modeling through dielectric continuum embedding and microsolvated quantum chemical modeling to explicit molecular dynamics highlight this situation. Here, we focus on quantum chemical microsolvation and discuss an explicit conformational sampling ansatz to make this approach systematic. For this purpose, we introduce an algorithm for the rolling and automated microsolvation of solutes. Our protocol takes conformational sampling and rearrangements in the solvent shell into account. Its reliability is assessed by monitoring the evolution of the spread and average of the observables of interest.

preprint2020arXiv

Tailored Coupled Cluster Theory in Varying Correlation Regimes

The tailored coupled cluster (TCC) approach is a promising ansatz that preserves the simplicity of single-reference coupled cluster theory, while incorporating a multi-reference wave function through amplitudes obtained from a preceding multi-configurational calculation. Here, we present a detailed analysis of the TCC wave function based on model systems, which require an accurate description of both static and dynamic correlation. We investigate the reliability of the TCC approach with respect to the exact wave function. In addition to the error in the electronic energy and standard coupled cluster diagnostics, we exploit the overlap of TCC and full configuration interaction wave functions as a quality measure. We critically review issues such as the required size of the active space, size-consistency, symmetry breaking in the wave function, and the dependence of TCC on the reference wave function. We observe that possible errors caused by symmetry breaking can be mitigated by employing the determinant with the largest weight in the active space as reference for the TCC calculation. We find the TCC model to be promising in calculations with active orbital spaces which include all orbitals with a large single-orbital entropy, even if the active spaces become very large and then may require modern active-space approaches that are not restricted to comparatively small numbers of orbitals. Furthermore, utilizing large active spaces can improve on the TCC wave-function approximation and reduce the size-consistency error, because the presence of highly excited determinants affects the accuracy of the coefficients of low-excited determinants in the active space.

preprint2020arXiv

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 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.

preprint2019arXiv

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 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.

preprint2019arXiv

Large-scale quantum-dynamics with matrix product states

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.

preprint2019arXiv

The Exploration of Chemical Reaction Networks

Modern computational chemistry has reached a stage at which massive exploration into chemical reaction space with unprecedented resolution with respect to the number of potentially relevant molecular structures has become possible. Various algorithmic advances have shown that such structural screenings must and can be automated and routinely carried out. This will replace the standard approach of manually studying a selected and restricted number of molecular structures for a chemical mechanism. The complexity of the task has led to many different approaches. However, all of them address the same general target, namely to produce a complete atomistic picture of the kinetics of a chemical process. It is the purpose of this overview to categorize the problems that are to be targeted and to identify the principle components and challenges of automated exploration machines so that the various existing approaches and future developments can be compared based on well-defined conceptual principles.

preprint2018arXiv

Minimum Energy Paths and Transition States by Curve Optimization

Transition states and minimum energy paths are essential to understand and predict chemical reactivity. Double-ended methods represent a standard approach for their determination. We introduce a new double-ended method that optimizes reaction paths described by curves. Unlike other methods, our approach optimizes the curve parameters rather than distinct structures along the path. With molecular paths represented as continuous curves, the optimization can benefit from the advantages of an integral-based formulation. We call this approach ReaDuct and demonstrate its applicability for molecular paths parametrized by B-spline curves.

preprint2016arXiv

Automated Selection of Active Orbital Spaces

One of the key challenges of quantum-chemical multi-configuration methods is the necessity to manually select orbitals for the active space. This selection requires both expertise and experience and can therefore impose severe limitations on the applicability of this most general class of ab initio methods. A poor choice of the active orbital space may yield even qualitatively wrong results. This is obviously a severe problem, especially for wave function methods that are designed to be systematically improvable. Here, we show how the iterative nature of the density matrix renormalization group combined with its capability to include up to about one hundred orbitals in the active space can be exploited for a systematic assessment and selection of active orbitals. These benefits allow us to implement an automated approach for active orbital space selection, which can turn multi-configuration models into black box approaches.

preprint2016arXiv

Measuring Multi-Configurational Character by Orbital Entanglement

One of the most critical tasks at the very beginning of a quantum chemical investigation is the choice of either a multi- or single-configurational method. Naturally, many proposals exist to define a suitable diagnostic of the multi-configurational character for various types of wave functions in order to assist this crucial decision. Here, we present a new orbital-entanglement based multi-configurational diagnostic termed $Z_{s(1)}$. The correspondence of orbital entanglement and static (or nondynamic) electron correlation permits the definition of such a diagnostic. We chose our diagnostic to meet important requirements such as well-defined limits for pure single-configurational and multi-configurational wave functions. The $Z_{s(1)}$ diagnostic can be evaluated from a partially converged, but qualitatively correct, and therefore inexpensive density matrix renormalization group wave function as in our recently presented automated active orbital selection protocol. Its robustness and the fact that it can be evaluated at low cost make this diagnostic a practical tool for routine applications.

preprint2016arXiv

Molecular propensity as a driver for explorative reactivity studies

Quantum chemical studies of reactivity involve calculations on a large number of molecular structures and comparison of their energies. Already the set-up of these calculations limits the scope of the results that one will obtain, because several system-specific variables such as the charge and spin need to be set prior to the calculation. For a reliable exploration of reaction mechanisms, a considerable number of calculations with varying global parameters must be taken into account, or important facts about the reactivity of the system under consideration can go undetected. For example, one could miss crossings of potential energy surfaces for different spin states or might not note that a molecule is prone to oxidation. Here, we introduce the concept of molecular propensity to account for the predisposition of a molecular system to react across different electronic states in certain nuclear configurations. Within our real-time quantum chemistry framework, we developed an algorithm that allows us to be alerted to such a propensity of a system under consideration.

preprint2016arXiv

Systematic Error Estimation for Chemical Reaction Energies

For the theoretical understanding of the reactivity of complex chemical systems accurate relative energies between intermediates and transition states are required. Despite its popularity, density functional theory (DFT) often fails to provide sufficiently accurate data, especially for molecules containing transition metals. Due to the huge number of intermediates that need to be studied for all but the simplest chemical processes, DFT is to date the only method that is computationally feasible. Here, we present a Bayesian framework for DFT that allows for error estimation of calculated properties. Since the optimal choice of parameters in present-day density functionals is strongly system dependent, we advocate for a system-focused re-parameterization. While, at first sight, this approach conflicts with the first-principles character of DFT that should make it in principle system independent, we deliberately introduce system dependence because we can then assign a stochastically meaningful error to the system-dependent parametrization that makes it non-arbitrary. By re-parameterizing a functional that was derived on a sound physical basis to a chemical system of interest we obtain a functional that yields reliable confidence intervals for reaction energies. We demonstrate our approach at the example of catalytic nitrogen fixation.

preprint2016arXiv

The Delicate Balance of Static and Dynamic Electron Correlation

Multi-configurational approaches yield universal wave function parameterizations that can qualitatively well describe electronic structures along reaction pathways. For quantitative results, multi-reference perturbation theory is required to capture dynamic electron correlation from the otherwise neglected virtual orbitals. Still, the overall accuracy suffers from the finite size and choice of the active orbital space and peculiarities of the perturbation theory. Fortunately, the electronic wave functions at equilibrium structures of reactants and products can often be well described by single-reference methods and hence are accessible to accurate coupled cluster calculations. Here, we calculate the heterolytic double dissociation energy of four 3d-metallocenes with the complete active space self-consistent field method and compare to highly accurate coupled cluster data. Our coupled cluster data are well within the experimental error bars. This accuracy can also be approached by complete active space calculations with an orbital selection based on information entropy measures. The entropy based active space selection is discussed in detail. We find a very subtle balance between static and dynamic electron correlation effects that emphasizes the need for algorithmic active space selection and that differs significantly from restricted active space results for identical active spaces reported in the literature.

preprint2015arXiv

Accelerating Wave Function Convergence in Interactive Quantum Chemical Reactivity Studies

The inherently high computational cost of iterative self-consistent-field (SCF) methods proves to be a critical issue delaying visual and haptic feedback in real-time quantum chemistry. In this work, we introduce two schemes for SCF acceleration. They provide a guess for the initial density matrix of the SCF procedure generated by extrapolation techniques. SCF optimizations then converge in fewer iterations, which decreases the execution time of the SCF optimization procedure. To benchmark the proposed propagation schemes, we developed a test bed for performing quantum chemical calculations on sequences of molecular structures mimicking real-time quantum chemical explorations. Explorations of a set of six model reactions employing the semi-empirical methods PM6 and DFTB3 in this testing environment showed that the proposed propagation schemes achieved speedups of up to thirty percent as a consequence of a reduced number of SCF iterations.

preprint2015arXiv

An Efficient Matrix Product Operator Representation of the Quantum-Chemical Hamiltonian

We describe how to efficiently construct the quantum chemical Hamiltonian operator in matrix product form. We present its implementation as a density matrix renormalization group (DMRG) algorithm for quantum chemical applications in a purely matrix product based framework. Existing implementations of DMRG for quantum chemistry are based on the traditional formulation of the method, which was developed from a viewpoint of Hilbert space decimation and attained a higher performance compared to straightforward implementations of matrix product based DMRG. The latter variationally optimizes a class of ansatz states known as matrix product states (MPS), where operators are correspondingly represented as matrix product operators (MPO). The MPO construction scheme presented here eliminates the previous performance disadvantages while retaining the additional flexibility provided by a matrix product approach; for example, the specification of expectation values becomes an input parameter. In this way, MPOs for different symmetries - abelian and non-abelian - and different relativistic and non-relativistic models may be solved by an otherwise unmodified program.

preprint2015arXiv

Heuristics-Guided Exploration of Reaction Mechanisms

For the investigation of chemical reaction networks, the efficient and accurate determination of all relevant intermediates and elementary reactions is mandatory. The complexity of such a network may grow rapidly, in particular if reactive species are involved that might cause a myriad of side reactions. Without automation, a complete investigation of complex reaction mechanisms is tedious and possibly unfeasible. Therefore, only the expected dominant reaction paths of a chemical reaction network (e.g., a catalytic cycle or an enzymatic cascade) are usually explored in practice. Here, we present a computational protocol that constructs such networks in a parallelized and automated manner. Molecular structures of reactive complexes are generated based on heuristic rules derived from conceptual electronic-structure theory and subsequently optimized by quantum chemical methods to produce stable intermediates of an emerging reaction network. Pairs of intermediates in this network that might be related by an elementary reaction according to some structural similarity measure are then automatically detected and subjected to an automated search for the connecting transition state. The results are visualized as an automatically generated network graph, from which a comprehensive picture of the mechanism of a complex chemical process can be obtained that greatly facilitates the analysis of the whole network. We apply our protocol to the Schrock dinitrogen-fixation catalyst to study alternative pathways of catalytic ammonia production.

preprint2015arXiv

Real-time feedback from iterative electronic structure calculations

Real-time feedback from iterative electronic structure calculations requires to mediate between the inherently unpredictable execution times of the iterative algorithm employed and the necessity to provide data in fixed and short time intervals for real-time rendering. We introduce the concept of a mediator as a component able to deal with infrequent and unpredictable reference data to generate reliable feedback. In the context of real-time quantum chemistry, the mediator takes the form of a surrogate potential that has the same local shape as the first-principles potential and can be evaluated efficiently to deliver atomic forces as real-time feedback. The surrogate potential is updated continuously by electronic structure calculations and guarantees to provide a reliable response to the operator for any molecular structure. To demonstrate the application of iterative electronic structure methods in real-time reactivity exploration, we implement self-consistent semi-empirical methods as the data source and apply the surrogate-potential mediator to deliver reliable real-time feedback.

preprint2015arXiv

Relativistic Kinetic-Balance Condition for Explicitly Correlated Basis Functions

This paper presents the derivation of a kinetic-balance condition for explicitly correlated basis functions employed in semi-classical relativistic calculations. Such a condition is important to ensure variational stability in algorithms based on the first-quantized Dirac theory of 1/2-fermions. We demonstrate that the kinetic-balance condition can be obtained from the row reduction process commonly applied to solve systems of linear equations. The resulting form of kinetic balance establishes a relation for the $4^N$ components of the spinor of an $N$-fermion system to the non-relativistic limit, which is in accordance with recent developments in the field of exact decoupling in relativistic orbital-based many-electron theory.

preprint2014arXiv

Determining Factors for the Accuracy of DMRG in Chemistry

The Density Matrix Renormalization Group (DMRG) algorithm has been a rising star for the accurate ab initio exploration of Born-Oppenheimer potential energy surfaces in theoretical chemistry. However, owing to its iterative numerical nature, pit falls, that can affect the accuracy of DMRG energies, need to be circumvented. Here, after a brief introduction into this quantum chemical method, we discuss criteria that determine the accuracy of DMRG calculations.

preprint2014arXiv

Electric Transition Dipole Moment in pre-Born-Oppenheimer Molecular Structure Theory

This paper presents the calculation of the electric transition dipole moment in a pre-Born-Oppenheimer framework. Electrons and nuclei are treated equally in terms of the parametrization of the non-relativistic total wave function, which is written as a linear combination of basis functions constructed with explicitly correlated Gaussian functions and the global vector representation. The integrals of the electric transition dipole moment are derived corresponding to these basis functions in both the length and the velocity representation. The complete derivation and the calculations are performed in laboratory-fixed Cartesian coordinates without relying on coordinates which separate the center of mass from the translationally invariant degrees of freedom. The effect of the overall motion is eliminated via translationally invariant integral expressions. As a numerical example the electric transition dipole moment is calculated between two rovibronic levels of the H2 molecule assignable to the lowest rovibrational states of the X ^1Sigma^+_g and B ^1Sigma^+_u electronic states in the clamped-nuclei framework. This is the first evaluation of this quantity in a full quantum mechanical treatment without relying on the Born-Oppenheimer approximation.

preprint2014arXiv

Gradient-Driven Molecule Construction: An Inverse Approach Applied to the Design of Small-Molecule Fixating Catalysts

Rational design of molecules and materials usually requires extensive screening of molecular structures for the desired property. The inverse approach to deduce a structure for a predefined property would be highly desirable, but is, unfortunately, not well-defined. However, feasible strategies for such an inverse design process may be successfully developed for specific purposes. We discuss options for calculating 'jacket' potentials that fulfill a predefined target requirement - a concept that we recently introduced [T. Weymuth, M. Reiher, MRS Proceediungs, 2013, 1524, DOI:10.1557/opl.2012.1764]. We consider the case of small-molecule activating transition metal catalysts. As a target requirement we choose the vanishing geometry gradients on all atoms of a subsystem consisting of a metal center binding the small molecule to be activated. The jacket potential can be represented within a full quantum model or by a sequence of approximations of which a field of electrostatic point charges is the simplest. In a second step, the jacket potential needs to be replaced by a chemically viable chelate-ligand structure for which the geometry gradients on all of its atoms are also required to vanish. In order to analyze the feasibility of this approach, we dissect a known dinitrogen-fixating catalyst to study possible design strategies that must eventually produce the known catalyst.

preprint2014arXiv

Hydrogen-Activation Mechanism of [Fe] Hydrogenase Revealed by Multi-Scale Modeling

When investigating the mode of hydrogen activation by [Fe] hydrogenases, not only the chemical reactivity at the active site is of importance but also the large-scale conformational change between the so-called open and closed conformations, which leads to a special spatial arrangement of substrate and iron cofactor. To study H2 activation, a complete model of the solvated and cofactor-bound enzyme in complex with the substrate methenyl-H4MPT+ was constructed. Both the closed and open conformations were simulated with classical molecular dynamics on the 100 ns time scale. Quantum-mechanics/molecular-mechanics calculations on snapshots then revealed the features of the active site that enable the facile H2 cleavage. The hydroxyl group of the pyridinol ligand can easily be deprotonated. With the deprotonated hydroxyl group and the structural arrangement in the closed conformation, H2 coordinated to the Fe center is subject to an ionic and orbital push-pull effect and can be rapidly cleaved with a concerted hydride transfer to methenyl-H4MPT+. An intermediary hydride species is not formed.

preprint2014arXiv

Interactive Chemical Reactivity Exploration

Elucidating chemical reactivity in complex molecular assemblies of a few hundred atoms is, despite the remarkable progress in quantum chemistry, still a major challenge. Black-box search methods to find intermediates and transition-state structures might fail in such situations because of the high-dimensionality of the potential energy surface. Here, we propose the concept of interactive chemical reactivity exploration to effectively introduce the chemist's intuition into the search process. We employ a haptic pointer device with force-feedback to allow the operator the direct manipulation of structures in three dimensions along with simultaneous perception of the quantum mechanical response upon structure modification as forces. We elaborate on the details of how such an interactive exploration should proceed and which technical difficulties need to be overcome. All reactivity-exploration concepts developed for this purpose have been implemented in the Samson programming environment.

preprint2014arXiv

Inverse Quantum Chemistry: Concepts and Strategies for Rational Compound Design

The rational design of molecules and materials is becoming more and more important. With the advent of powerful computer systems and sophisticated algorithms, quantum chemistry plays an important role in rational design. While traditional quantum chemical approaches predict the properties of a predefined molecular structure, the goal of inverse quantum chemistry is to find a structure featuring one or more desired properties. Herein, we review inverse quantum chemical approaches proposed so far and discuss their advantages as well as their weaknesses.

preprint2014arXiv

Mode-Tracking Based Stationary-Point Optimization

In this work, we present a transition-state optimization protocol based on the Mode-Tracking algorithm [J. Chem. Phys. 118 (2003) 1634]. By calculating only the eigenvector of interest instead of diagonalizing the full Hessian matrix and performing an eigenvector following search based on the selectively calculated vector, we can efficiently optimize transition-state structures. The initial guess structures and eigenvectors are either chosen from a linear interpolation between the reactant and product structures, from a nudged-elastic band search, from a constrained-optimization scan, or from the minimum-energy structures. Alternatively, initial guess vectors based on chemical intuition may be defined. We then iteratively refine the selected vectors by the Davidson subspace iteration technique. This procedure accelerates finding transition states for large molecules of a few hundred atoms. It is also beneficial in cases where the starting structure is very different from the transition-state structure or where the desired vector to follow is not the one with lowest eigenvalue. Explorative studies of reaction pathways are feasible by following manually constructed molecular distortions.

preprint2014arXiv

New Benchmark Set of Transition-Metal Coordination Reactions for the Assessment of Density Functionals

We present the WCCR10 data set of ten ligand dissociation energies of large cationic transition metal complexes for the assessment of approximate exchange--correlation functionals. We analyze nine popular functionals, namely BP86, BP86-D3, B3LYP, B3LYP-D3, B97-D-D2, PBE, TPSS, PBE0, and TPSSh by mutual comparison and by comparison to experimental gas-phase data measured with well-known precision. The comparison of all calculated data reveals a large, system-dependent scattering of results with nonnegligible consequences for computational chemistry studies on transition metal compounds. Considering further the comparison with experiment, the non-empirical functionals PBE and TPSS turn out to be among the best functionals for our reference data set. The deviation can be lowered further by including Hartree--Fock exchange. Accordingly, PBE0 and TPSSh are the two most accurate functionals for our test set, but also these functionals exhibit deviations from experiment by up to 50 kJ/mol for individual reactions. As an important result we found no functional to be reliable for all reactions. Furthermore, for some of the ligand dissociation energies studied in this work, dispersion corrections yield results which increase the deviation from experiment. This deviation increases further if structure optimization including dispersion corrections is performed. Finally, we compare our results to other benchmark studies and highlight that the performance assessed for different density functionals depends significantly on the reference molecule set chosen.

preprint2014arXiv

Self-Consistent Embedding of Density-Matrix Renormalization Group Wavefunctions in a Density Functional Environment

We present the first implementation of a density matrix renormalization group algorithm embedded in an environment described by density functional theory. The frozen density embedding scheme is used with a freeze-and-thaw strategy for a self-consistent polarization of the orbital-optimized wavefunction and the environmental densities with respect to each other.

preprint2013arXiv

An efficient implementation of two-component relativistic exact-decoupling methods for large molecules

We present an efficient algorithm for one- and two-component relativistic exact-decoupling calculations. Spin-orbit coupling is thus taken into account for the evaluation of relativistically transformed (one-electron) Hamiltonian. As the relativistic decoupling transformation has to be evaluated with primitive functions, the construction of the relativistic one-electron Hamiltonian becomes the bottleneck of the whole calculation for large molecules. For the established exact-decoupling protocols, a minimal matrix operation count is established and discussed in detail. Furthermore, we apply our recently developed local DLU scheme [J. Chem. Phys. 136 (2012) 244108] to accelerate this step. With our new implementation two-component relativistic density functional calculations can be performed invoking the resolution-of-identity density-fitting approximation and (Abelian as well as non-Abelian) point group symmetry to accelerate both the exact-decoupling and the two-electron part. The capability of our implementation is illustrated at the example of silver clusters with up to 309 atoms, for which the cohesive energy is calculated and extrapolated to the bulk.

preprint2013arXiv

Four-Component Density Matrix Renormalization Group

We present the first implementation of the relativistic quantum chemical two- and four-component density matrix renormalization group algorithm that includes a variational description of scalar-relativistic effects and spin-orbit coupling. Numerical results based on the four-component Dirac-Coulomb Hamiltonian are presented for the standard reference molecule for correlated relativistic benchmarks: thallium hydride.

preprint2013arXiv

Orbital entanglement in bond-formation processes

The accurate calculation of the (differential) correlation energy is central to the quantum chemical description of bond-formation and bond-dissociation processes. In order to estimate the quality of single- and multi-reference approaches for this purpose, various diagnostic tools have been developed. In this work, we elaborate on our previous observation [J. Phys. Chem. Lett. 3, 3129 (2012)] that one- and two-orbital-based entanglement measures provide quantitative means for the assessment and classification of electron correlation effects among molecular orbitals. The dissociation behavior of some prototypical diatomic molecules features all types of correlation effects relevant for chemical bonding. We demonstrate that our entanglement analysis is convenient to dissect these electron correlation effects and to provide a conceptual understanding of bond-forming and bond-breaking processes from the point of view of quantum information theory.

preprint2013arXiv

Unravelling the quantum-entanglement effect of noble gas coordination on the spin ground state of CUO

The accurate description of the complexation of the CUO molecule by Ne and Ar noble gas matrices represents a challenging task for present-day quantum chemistry. Especially, the accurate prediction of the spin ground state of different CUO--noble-gas complexes remains elusive. In this work, the interaction of the CUO unit with the surrounding noble gas matrices is investigated in terms of complexation energies and dissected into its molecular orbital quantum entanglement patterns. Our analysis elucidates the anticipated singlet--triplet ground-state reversal of the CUO molecule diluted in different noble gas matrices and demonstrates that the strongest uranium-noble gas interaction is found for CUOAr4 in its triplet configuration.

preprint2012arXiv

Accurate ab initio spin densities

We present an approach for the calculation of spin density distributions for molecules that require very large active spaces for a qualitatively correct description of their electronic structure. Our approach is based on the density-matrix renormalization group (DMRG) algorithm to calculate the spin density matrix elements as basic quantity for the spatially resolved spin density distribution. The spin density matrix elements are directly determined from the second-quantized elementary operators optimized by the DMRG algorithm. As an analytic convergence criterion for the spin density distribution, we employ our recently developed sampling-reconstruction scheme [J. Chem. Phys. 2011, 134, 224101] to build an accurate complete-active-space configuration-interaction (CASCI) wave function from the optimized matrix product states. The spin density matrix elements can then also be determined as an expectation value employing the reconstructed wave function expansion. Furthermore, the explicit reconstruction of a CASCI-type wave function provides insights into chemically interesting features of the molecule under study such as the distribution of $α$- and $β$-electrons in terms of Slater determinants, CI coefficients, and natural orbitals. The methodology is applied to an iron nitrosyl complex which we have identified as a challenging system for standard approaches [J. Chem. Theory Comput. 2011, 7, 2740].

preprint2012arXiv

Elimination of the Translational Kinetic Energy Contamination in pre-Born-Oppenheimer Calculations

In this paper we present a simple strategy for the elimination of the translational kinetic energy contamination of the total energy in pre-Born--Oppenheimer calculations carried out in laboratory-fixed Cartesian coordinates (LFCCs). The simple expressions for the coordinates and the operators are thus preserved throughout the calculations, while the mathematical form and the parametrisation of the basis functions are chosen so that the translational and rotational invariances are respected. The basis functions are constructed using explicitly correlated Gaussian functions (ECGs) and the global vector representation. First, we observe that it is not possible to parametrise the ECGs so that the system is at rest in LFCCs and at the same time the basis functions are square-integrable with a non-vanishing norm. Then, we work out a practical strategy to circumvent this problem by making use of the properties of the linear transformation between the LFCCs and translationally invariant and center-of-mass Cartesian coordinates as well as the transformation properties of the corresponding basis function parameter matrices. By exploiting these formal mathematical relationships we can identify and separate the translational contamination terms in the matrix representation of the kinetic energy operator in the LFCC formalism. We present numerical examples for the translational contamination and its elimination for the two lowest rotational energy levels of the singlet hydrogen molecule, corresponding to para- and ortho-H2, respectively, treated as four-particle quantum systems.

preprint2012arXiv

Entanglement Measures for Single- and Multi-Reference Correlation Effects

Electron correlation effects are essential for an accurate ab initio description of molecules. A quantitative a priori knowledge of the single- or multi-reference nature of electronic structures as well as of the dominant contributions to the correlation energy can facilitate the decision regarding the optimum quantum chemical method of choice. We propose concepts from quantum information theory as orbital entanglement measures that allow us to evaluate the single- and multi-reference character of any molecular structure in a given orbital basis set. By studying these measures we can detect possible artifacts of small active spaces.

preprint2012arXiv

Local Relativistic Exact Decoupling

We present a systematic hierarchy of approximations for {\it local} exact-decoupling of four-component quantum chemical Hamiltonians based on the Dirac equation. Our ansatz reaches beyond the trivial local approximation that is based on a unitary transformation of only the atomic block-diagonal part of the Hamiltonian. Systematically, off-diagonal Hamiltonian matrix blocks can be subjected to a unitary transformation to yield relativistically corrected matrix elements. The full hierarchy is investigated with respect to the accuracy reached for the electronic energy and molecular properties on a balanced test molecule set that comprises molecules with heavy elements in different bonding situations. Our atomic (local) assembly of the unitary transformation needed for exact decoupling provides an excellent local approximation for any relativistic exact-decoupling approach. Its order-$N^2$ scaling can be further reduced to linear scaling by employing the neighboring-atomic-blocks approximation. Therefore, it is an efficient relativistic method perfectly well suited for relativistic calculations on large molecules. If a large molecule contains many light atoms (typically hydrogen atoms), the computational costs can be further reduced by employing a well-defined non-relativistic approximation for these light atoms without significant loss of accuracy.

preprint2012arXiv

Optimized Unrestricted Kohn-Sham Potentials from Ab Initio Spin Densities

The reconstruction of the exchange-correlation potential from accurate ab initio electron densities can provide insights into the limitations of the currently available approximate functionals and provide guidance for devising improved approximations for density-functional theory (DFT). For open-shell systems, the spin density is introduced as an additional fundamental variable in Spin-DFT. Here, we consider the reconstruction of the corresponding unrestricted Kohn-Sham potentials from accurate ab initio spin densities. In particular, we investigate whether it is possible to reconstruct the spin exchange-correlation potential, which determines the spin density in spin-unrestricted Kohn-Sham-DFT, despite the numerical difficulties inherent to the optimization of potentials with finite orbital basis sets. We find that the recently developed scheme for unambiguously singling out an optimal optimized potential [J. Chem. Phys. 135, 244102 (2011)] can provide such spin potentials accurately. This is demonstrated for two test cases, the lithium atom and the dioxygen molecule, and target (spin) densities from Full-CI and CASSCF calculations, respectively.

preprint2012arXiv

Real-time Quantum Chemistry

Significant progress in the development of efficient and fast algorithms for quantum chemical calculations has been made in the past two decades. The main focus has always been the desire to be able to treat ever larger molecules or molecular assemblies---especially linear and sub-linear scaling techniques are devoted to the accomplishment of this goal. However, as many chemical reactions are rather local, they usually involve only a limited number of atoms so that models of about two hundred (or even less) atoms embedded in a suitable environment are sufficient to study their mechanisms. Thus, the system size does not need to be enlarged, but remains constant for reactions of this type that can be described by less than two hundred atoms. The question then arises how fast one can obtain the quantum chemical results. This question is not directly answered by linear-scaling techniques. In fact, ideas such as haptic quantum chemistry or interactive quantum chemistry require an immediate provision of quantum chemical information which demands the calculation of data in "real time". In this perspective, we aim at a definition of real-time quantum chemistry, explore its realm and eventually discuss applications in the field of haptic quantum chemistry. For the latter we elaborate whether a direct approach is possible by virtue of real-time quantum chemistry.

preprint2012arXiv

Spin in Density-Functional Theory

The accurate description of open-shell molecules, in particular of transition metal complexes and clusters, is still an important challenge for quantum chemistry. While density-functional theory (DFT) is widely applied in this area, the sometimes severe limitations of its currently available approximate realizations often preclude its application as a predictive theory. Here, we review the foundations of DFT applied to open-shell systems, both within the nonrelativistic and the relativistic framework. In particular, we provide an in-depth discussion of the exact theory, with a focus on the role of the spin density and possibilities for targeting specific spin states. It turns out that different options exist for setting up Kohn-Sham DFT schemes for open-shell systems, which imply different definitions of the exchange-correlation energy functional and lead to different exact conditions on this functional. Finally, we suggest some possible directions for future developments.

preprint2011arXiv

Construction of CASCI-type wave functions for very large active spaces

We present a procedure to construct a configuration-interaction expansion containing arbitrary excitations from an underlying full-configuration-interaction-type wave function defined for a very large active space. Our procedure is based on the density-matrix renormalization group (DMRG) algorithm that provides the necessary information in terms of the eigenstates of the reduced density matrices to calculate the coefficient of any basis state in the many-particle Hilbert space. Since the dimension of the Hilbert space scales binomially with the size of the active space, a sophisticated Monte Carlo sampling routine is employed. This sampling algorithm can also construct such configuration-interaction-type wave functions from any other type of tensor network states. The configuration-interaction information obtained serves several purposes. It yields a qualitatively correct description of the molecule's electronic structure, it allows us to analyze DMRG wave functions converged for the same molecular system but with different parameter sets (e.g., different numbers of active-system (block) states), and it can be considered a balanced reference for the application of a subsequent standard multi-reference configuration-interaction method.

preprint2010arXiv

Complete-Graph Tensor Network States: A New Fermionic Wave Function Ansatz for Molecules

We present a new class of tensor network states that are specifically designed to capture the electron correlation of a molecule of arbitrary structure. In this ansatz, the electronic wave function is represented by a Complete-Graph Tensor Network (CGTN) ansatz which implements an efficient reduction of the number of variational parameters by breaking down the complexity of the high-dimensional coefficient tensor of a full-configuration-interaction (FCI) wave function. We demonstrate that CGTN states approximate ground states of molecules accurately by comparison of the CGTN and FCI expansion coefficients. The CGTN parametrization is not biased towards any reference configuration in contrast to many standard quantum chemical methods. This feature allows one to obtain accurate relative energies between CGTN states which is central to molecular physics and chemistry. We discuss the implications for quantum chemistry and focus on the spin-state problem. Our CGTN approach is applied to the energy splitting of states of different spin for methylene and the strongly correlated ozone molecule at a transition state structure. The parameters of the tensor network ansatz are variationally optimized by means of a parallel-tempering Monte Carlo algorithm.