Researcher profile

Markus Reiher

Markus Reiher contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
31works
0followers
13topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

31 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

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.

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.