Showing 1–17 of 17 results for author: Faulstich, F M

Exploring Ground and Excited States via Single Reference Coupled-Cluster Theory and Algebraic Geometry

Authors: Svala Sverrisdóttir, Fabian M. Faulstich

Abstract: The exploration of the root structure of coupled cluster equations holds both foundational and practical significance for computational quantum chemistry. This study provides insight into the intricate root structures of these non-linear equations at both the CCD and CCSD level of theory. We utilize computational techniques from algebraic geometry, specifically the monodromy and parametric homotop… ▽ More The exploration of the root structure of coupled cluster equations holds both foundational and practical significance for computational quantum chemistry. This study provides insight into the intricate root structures of these non-linear equations at both the CCD and CCSD level of theory. We utilize computational techniques from algebraic geometry, specifically the monodromy and parametric homotopy continuation methods, to calculate the full solution set. We compare the computed CC roots against various established theoretical upper bounds, shedding light on the accuracy and efficiency of these bounds. We hereby focus on the dissociation processes of four-electron systems such as (H$_2$)$_2$ in both D$_{2{\rm h}}$ and D$_{\infty {\rm h}}$ configurations, H$_4$ symmetrically distorted on a circle, and lithium hydride. We moreover investigate the ability of single-reference coupled cluster solutions to approximate excited state energies. We find that multiple CC roots describe energies of excited states with high accuracy. Remarkably, our investigations reveal that for systems like lithium hydride, CC not only provides high-accuracy approximations to several excited state energies but also to the states themselves. △ Less

Submitted 16 May, 2024; originally announced May 2024.

A static quantum embedding scheme based on coupled cluster theory

Authors: Avijit Shee, Fabian M. Faulstich, Birgitta Whaley, Lin Lin, Martin Head-Gordon

Abstract: We develop a static quantum embedding scheme, utilizing projection equations to solve coupled cluster (CC) amplitudes. To reduce the computational cost (for example, of a large basis set calculation), we solve the local fragment problem using a high-level coupled cluster method and address the environment problem with a lower-level Møller-Plesset (MP) perturbative method. This embedding approach i… ▽ More We develop a static quantum embedding scheme, utilizing projection equations to solve coupled cluster (CC) amplitudes. To reduce the computational cost (for example, of a large basis set calculation), we solve the local fragment problem using a high-level coupled cluster method and address the environment problem with a lower-level Møller-Plesset (MP) perturbative method. This embedding approach is consistently formulated within the coupled cluster framework and will be called MP-CC. We demonstrate the effectiveness of our method through several prototypical molecular examples by analyzing a global quantity, that is, the total correlation energy of the system in the study of potential energy curves (PEC) and thermochemical reaction energies. We have shown that our method can achieve comparable accuracy both with a small and large basis set when the fragment Hilbert space size remains the same. Additionally, our results indicate that increasing the fragment size can systematically enhance the accuracy of observables, approaching the precision of the full coupled cluster solver. △ Less

Submitted 13 April, 2024; originally announced April 2024.

Augmented Lagrangian method for coupled-cluster

Authors: Fabian M. Faulstich, Yuehaw Khoo, Kangbo Li

Abstract: We propose to improve the convergence properties of the single-reference coupled cluster (CC) method through an augmented Lagrangian formalism. The conventional CC method changes a linear high-dimensional eigenvalue problem with exponential size into a problem of determining the roots of a nonlinear system of equations that has a manageable size. However, current numerical procedures for solving t… ▽ More We propose to improve the convergence properties of the single-reference coupled cluster (CC) method through an augmented Lagrangian formalism. The conventional CC method changes a linear high-dimensional eigenvalue problem with exponential size into a problem of determining the roots of a nonlinear system of equations that has a manageable size. However, current numerical procedures for solving this system of equations to get the lowest eigenvalue suffer from two practical issues: First, solving the CC equations may not converge, and second, when converging, they may converge to other -- potentially unphysical -- states, which are stationary points of the CC energy expression. We show that both issues can be dealt with when a suitably defined energy is minimized in addition to solving the original CC equations. We further propose an augmented Lagrangian method for coupled cluster (alm-CC) to solve the resulting constrained optimization problem. We numerically investigate the proposed augmented Lagrangian formulation showing that the convergence towards the ground state is significantly more stable and that the optimization procedure is less susceptible to local minima. Furthermore, the computational cost of alm-CC is comparable to the conventional CC method. △ Less

Submitted 24 March, 2024; originally announced March 2024.

Comments: 17 pages, 8 figures

Recent mathematical advances in coupled cluster theory

Authors: Fabian M. Faulstich

Abstract: This article presents an in-depth educational overview of the latest mathematical developments in coupled cluster (CC) theory, beginning with Schneider's seminal work from 2009 that introduced the first local analysis of CC theory. We offer a tutorial review of second quantization and the CC ansatz, laying the groundwork for understanding the mathematical basis of the theory. This is followed by a… ▽ More This article presents an in-depth educational overview of the latest mathematical developments in coupled cluster (CC) theory, beginning with Schneider's seminal work from 2009 that introduced the first local analysis of CC theory. We offer a tutorial review of second quantization and the CC ansatz, laying the groundwork for understanding the mathematical basis of the theory. This is followed by a detailed exploration of the most recent mathematical advancements in CC theory.Our review starts with an in-depth look at the local analysis pioneered by Schneider which has since been applied to analyze various CC methods. We then move on to discuss the graph-based framework for CC methods developed by Csirik and Laestadius. This framework provides a comprehensive platform for comparing different CC methods, including multireference approaches. Next, we delve into the latest numerical analysis results analyzing the single reference CC method developed by Hassan, Maday, and Wang. This very general approach is based on the invertibility of the CC function's Fréchet derivative. We conclude the article with a discussion on the recent incorporation of algebraic geometry into CC theory, highlighting how this novel and fundamentally different mathematical perspective has furthered our understanding and provides exciting pathways to new computational approaches. △ Less

Submitted 14 January, 2024; originally announced January 2024.

Homotopy continuation methods for coupled-cluster theory in quantum chemistry

Authors: Fabian M. Faulstich, Andre Laestadius

Abstract: Homotopy methods have proven to be a powerful tool for understanding the multitude of solutions provided by the coupled-cluster polynomial equations. This endeavor has been pioneered by quantum chemists that have undertaken both elaborate numerical as well as mathematical investigations. Recently, from the perspective of applied mathematics, new interest in these approaches has emerged using both… ▽ More Homotopy methods have proven to be a powerful tool for understanding the multitude of solutions provided by the coupled-cluster polynomial equations. This endeavor has been pioneered by quantum chemists that have undertaken both elaborate numerical as well as mathematical investigations. Recently, from the perspective of applied mathematics, new interest in these approaches has emerged using both topological degree theory and algebraically oriented tools. This article provides an overview of describing the latter development. △ Less

Submitted 23 June, 2023; originally announced June 2023.

Some mathematical insights on Density Matrix Embedding Theory

Authors: Eric Cancès, Fabian M. Faulstich, Alfred Kirsch, Eloïse Letournel, Antoine Levitt

Abstract: This article provides the first mathematical analysis of the Density Matrix Embedding Theory (DMET) method. We prove that, under certain assumptions, (i) the exact ground-state density matrix is a fixed-point of the DMET map for non-interacting systems, (ii) there exists a unique physical solution in the weakly-interacting regime, and (iii) DMET is exact at first order in the coupling parameter. W… ▽ More This article provides the first mathematical analysis of the Density Matrix Embedding Theory (DMET) method. We prove that, under certain assumptions, (i) the exact ground-state density matrix is a fixed-point of the DMET map for non-interacting systems, (ii) there exists a unique physical solution in the weakly-interacting regime, and (iii) DMET is exact at first order in the coupling parameter. We provide numerical simulations to support our results and comment on the physical meaning of the assumptions under which they hold true. We show that the violation of these assumptions may yield multiple solutions of the DMET equations. We moreover introduce and discuss a specific N-representability problem inherent to DMET. △ Less

Submitted 2 October, 2023; v1 submitted 25 May, 2023; originally announced May 2023.

Comments: 40 pages, 11 figures

MSC Class: 35Q40; 39B42; 49M27

The $S$-diagnostic -- an a posteriori error assessment for single-reference coupled-cluster methods

Authors: Fabian M. Faulstich, Håkon E. Kristiansen, Mihaly A. Csirik, Simen Kvaal, Thomas Bondo Pedersen, Andre Laestadius

Abstract: We propose a novel a posteriori error assessment for the single-reference coupled-cluster (SRCC) method called the $S$-diagnostic. We provide a derivation of the $S$-diagnostic that is rooted in the mathematical analysis of different SRCC variants. We numerically scrutinized the $S$-diagnostic, testing its performance for (1) geometry optimizations, (2) electronic correlation simulations of system… ▽ More We propose a novel a posteriori error assessment for the single-reference coupled-cluster (SRCC) method called the $S$-diagnostic. We provide a derivation of the $S$-diagnostic that is rooted in the mathematical analysis of different SRCC variants. We numerically scrutinized the $S$-diagnostic, testing its performance for (1) geometry optimizations, (2) electronic correlation simulations of systems with varying numerical difficulty, and (3) the square-planar copper complexes [CuCl$_4$]$^{2-}$, [Cu(NH$_3$)$_4$]$^{2+}$, and [Cu(H$_2$O)$_4$]$^{2+}$. Throughout the numerical investigations, the $S$-diagnostic is compared to other SRCC diagnostic procedures, that is, the $T_1$, $D_1$, and $D_2$ diagnostics as well as different indices of multi-determinantal and multi-reference character in coupled-cluster theory. Our numerical investigations show that the $S$-diagnostic outperforms the $T_1$, $D_1$, and $D_2$ diagnostics and is comparable to the indices of multi-determinantal and multi-reference character in coupled-cluster theory in their individual fields of applicability. The experiments investigating the performance of the $S$-diagnostic for geometry optimizations using SRCC reveal that the $S$-diagnostic correlates well with different error measures at a high level of statistical relevance. The experiments investigating the performance of the $S$-diagnostic for electronic correlation simulations show that the $S$-diagnostic correctly predicts strong multi-reference regimes. The $S$-diagnostic moreover correctly detects the successful SRCC computations for [CuCl$_4$]$^{2-}$, [Cu(NH$_3$)$_4$]$^{2+}$, and [Cu(H$_2$O)$_4$]$^{2+}$, which have been known to be misdiagnosed by $T_1$ and $D_1$ diagnostics in the past. This shows that the $S$-diagnostic is a promising candidate for an a posteriori diagnostic for SRCC calculations. △ Less

Submitted 26 January, 2023; originally announced January 2023.

Coupled cluster theory: Towards an algebraic geometry formulation

Authors: Fabian M. Faulstich, Mathias Oster

Abstract: Coupled cluster theory produced arguably the most widely used high-accuracy computational quantum chemistry methods. Despite the approach's overall great computational success, its mathematical understanding is so far limited to results within the realm of functional analysis. The coupled cluster amplitudes, which are the targeted objects in coupled cluster theory, correspond to solutions to the c… ▽ More Coupled cluster theory produced arguably the most widely used high-accuracy computational quantum chemistry methods. Despite the approach's overall great computational success, its mathematical understanding is so far limited to results within the realm of functional analysis. The coupled cluster amplitudes, which are the targeted objects in coupled cluster theory, correspond to solutions to the coupled cluster equations, which is a system of polynomial equations of at most degree four. The high dimensionality of the electronic Schrödinger equation and the non-linearity of the coupled cluster ansatz have so far stalled a formal analysis of this polynomial system. In this article, we present algebraic investigations that shed light on the coupled cluster equations and the root structure of this ansatz. This is of importance for the a posteriori evaluation of coupled cluster calculations. To that end, we investigate the root structure by means of Newton polytopes. We derive a general v-description, which is subsequently turned into an h-description for explicit examples. This perspective reveals an apparent connection between Pauli's exclusion principle and the geometrical structure of the Newton polytopes. We also propose an alternative characterization of the coupled cluster equations projected onto singles and doubles as cubic polynomials on an algebraic variety with certain sparsity patterns. Moreover, we provide numerical simulations of two computationally tractable systems, namely, the two electrons in four spin-orbitals system and the three electrons in six spin-orbitals system. These simulations provide novel insight into the root structure of the coupled cluster solutions when the coupled cluster ansatz is truncated. △ Less

Submitted 28 March, 2024; v1 submitted 18 November, 2022; originally announced November 2022.

MSC Class: 12D10; 14Q20; 90C53; 81-08; 81-10

Journal ref: SIAM Journal on Applied Algebra and Geometry 8.1 (2024): 138-188

Interacting models for twisted bilayer graphene: a quantum chemistry approach

Authors: Fabian M. Faulstich, Kevin D. Stubbs, Qinyi Zhu, Tomohiro Soejima, Rohit Dilip, Huanchen Zhai, Raehyun Kim, Michael P. Zaletel, Garnet Kin-Lic Chan, Lin Lin

Abstract: The nature of correlated states in twisted bilayer graphene (TBG) at the magic angle has received intense attention in recent years. We present a numerical study of an interacting Bistritzer-MacDonald (IBM) model of TBG using a suite of methods in quantum chemistry, including Hartree-Fock, coupled cluster singles, doubles (CCSD), and perturbative triples (CCSD(T)), as well as a quantum chemistry f… ▽ More The nature of correlated states in twisted bilayer graphene (TBG) at the magic angle has received intense attention in recent years. We present a numerical study of an interacting Bistritzer-MacDonald (IBM) model of TBG using a suite of methods in quantum chemistry, including Hartree-Fock, coupled cluster singles, doubles (CCSD), and perturbative triples (CCSD(T)), as well as a quantum chemistry formulation of the density matrix renormalization group method (DMRG). Our treatment of TBG is agnostic to gauge choices, and hence we present a new gauge-invariant formulation to detect the spontaneous symmetry breaking in interacting models. To benchmark our approach, we focus on a simplified spinless, valleyless IBM model. At integer filling ($ν=0$), all numerical methods agree in terms of energy and $C_{2z} \mathcal{T}$ symmetry breaking. Additionally, as part of our benchmarking, we explore the impact of different schemes for removing ``double-counting'' in the IBM model. Our results at integer filling suggest that cross-validation of different IBM models may be needed for future studies of the TBG system. After benchmarking our approach at integer filling, we perform the first systematic study of the IBM model near integer filling (for $|ν|< 0.2$). In this regime, we find that the ground state can be in a metallic and $C_{2z} \mathcal{T}$ symmetry breaking phase. The ground state appears to have low entropy, and therefore can be relatively well approximated by a single Slater determinant. Furthermore, we observe many low entropy states with energies very close to the ground state energy in the near integer filling regime. △ Less

Submitted 16 November, 2022; originally announced November 2022.

On the pure state $v$-representability of density matrix embedding theory

Authors: Fabian M. Faulstich, Raehyun Kim, Zhi-Hao Cui, Zaiwen Wen, Garnet Kin-Lic Chan, Lin Lin

Abstract: Density matrix embedding theory (DMET) formally requires the matching of density matrix blocks obtained from high-level and low-level theories, but this is sometimes not achievable in practical calculations. In such a case, the global band gap of the low-level theory vanishes, and this can require additional numerical considerations. We find that both the violation of the exact matching condition… ▽ More Density matrix embedding theory (DMET) formally requires the matching of density matrix blocks obtained from high-level and low-level theories, but this is sometimes not achievable in practical calculations. In such a case, the global band gap of the low-level theory vanishes, and this can require additional numerical considerations. We find that both the violation of the exact matching condition and the vanishing low-level gap are related to the assumption that the high-level density matrix blocks are non-interacting pure-state $v$-representable (NI-PS-V), which assumes that the low-level density matrix is constructed following the Aufbau principle. In order to relax the NI-PS-V condition, we develop an augmented Lagrangian method to match the density matrix blocks without referring to the Aufbau principle. Numerical results for 2D Hubbard and hydrogen model systems indicate that in some challenging scenarios, the relaxation of the Aufbau principle directly leads to exact matching of the density matrix blocks, which also yields improved accuracy. △ Less

Submitted 13 January, 2022; v1 submitted 18 October, 2021; originally announced October 2021.

Discontinuous Galerkin method with Voronoi partitioning for Quantum Simulation of Chemistry

Authors: Fabian M. Faulstich, Xiaojie Wu, Lin Lin

Abstract: Molecular orbitals based on the linear combination of Gaussian type orbitals are arguably the most employed discretization in quantum chemistry simulations, both on quantum and classical devices. To circumvent a potentially dense two-body interaction tensor and obtain lower asymptotic costs for quantum simulations of chemistry, the discontinuous Galerkin (DG) procedure using a rectangular partitio… ▽ More Molecular orbitals based on the linear combination of Gaussian type orbitals are arguably the most employed discretization in quantum chemistry simulations, both on quantum and classical devices. To circumvent a potentially dense two-body interaction tensor and obtain lower asymptotic costs for quantum simulations of chemistry, the discontinuous Galerkin (DG) procedure using a rectangular partitioning strategy was recently piloted [McClean et al, New J. Phys. 22, 093015, 2020]. The DG approach interpolates in a controllable way between a compact description of the two-body interaction tensor through molecular orbitals and a diagonal characterization through primitive basis sets, such as a planewave dual basis set. The DG procedure gives rise to a block-diagonal representation of the two-body interaction with reduced number of two-electron repulsion integrals, which in turn reduces the cost of quantum simulations. In the present work we extend this approach to be applicable to molecular and crystalline systems of arbitrary geometry. We take advantage of the flexibility of the planewave dual basis set, and combine the discontinuous Galerkin procedure with a general partitioning strategy based on the Voronoi decomposition. We numerically investigate the performance, at the mean-field and correlated levels, with quasi-1D, 2D and 3D partitions using hydrogen chains, H$_4$, CH$_4$ as examples, respectively. We also apply the method to graphene as a prototypical example of crystalline systems. △ Less

Submitted 31 October, 2020; originally announced November 2020.

One-Dimensional Lieb-Oxford Bounds

Authors: Andre Laestadius, Fabian M Faulstich

Abstract: We investigate and prove Lieb-Oxford bounds in one dimension by studying convex potentials that approximate the ill-defined Coulomb potential. A Lieb-Oxford inequality establishes a bound of the indirect interaction energy for electrons in terms of the one-body particle density $ρ_ψ$ of a wave function $ψ$. Our results include modified soft Coulomb potential and regularized Coulomb potential. For… ▽ More We investigate and prove Lieb-Oxford bounds in one dimension by studying convex potentials that approximate the ill-defined Coulomb potential. A Lieb-Oxford inequality establishes a bound of the indirect interaction energy for electrons in terms of the one-body particle density $ρ_ψ$ of a wave function $ψ$. Our results include modified soft Coulomb potential and regularized Coulomb potential. For these potentials, we establish Lieb-Oxford-type bounds utilizing logarithmic expressions of the particle density. Furthermore, a previous conjectured form $I_\mathrm{xc}(ψ)\geq - C_1 \int_{\mathbb R} ρ_ψ(x)^{2} \mathrm{d}x$ is discussed for different convex potentials. △ Less

Submitted 4 June, 2020; v1 submitted 4 October, 2019; originally announced October 2019.

Discontinuous Galerkin discretization for quantum simulation of chemistry

Authors: Jarrod R. McClean, Fabian M. Faulstich, Qinyi Zhu, Bryan O'Gorman, Yiheng Qiu, Steven R. White, Ryan Babbush, Lin Lin

Abstract: Methods for electronic structure based on Gaussian and molecular orbital discretizations offer a well established, compact representation that forms much of the foundation of correlated quantum chemistry calculations on both classical and quantum computers. Despite their ability to describe essential physics with relatively few basis functions, these representations can suffer from a quartic growt… ▽ More Methods for electronic structure based on Gaussian and molecular orbital discretizations offer a well established, compact representation that forms much of the foundation of correlated quantum chemistry calculations on both classical and quantum computers. Despite their ability to describe essential physics with relatively few basis functions, these representations can suffer from a quartic growth of the number of integrals. Recent results have shown that, for some quantum and classical algorithms, moving to representations with diagonal two-body operators can result in dramatically lower asymptotic costs, even if the number of functions required increases significantly. We introduce a way to interpolate between the two regimes in a systematic and controllable manner, such that the number of functions is minimized while maintaining a block diagonal structure of the two-body operator and desirable properties of an original, primitive basis. Techniques are analyzed for leveraging the structure of this new representation on quantum computers. Empirical results for hydrogen chains suggest a scaling improvement from $O(N^{4.5})$ in molecular orbital representations to $O(N^{2.6})$ in our representation for quantum evolution in a fault-tolerant setting, and exhibit a constant factor crossover at 15 to 20 atoms. Moreover, we test these methods using modern density matrix renormalization group methods classically, and achieve excellent accuracy with respect to the complete basis set limit with a speedup of 1-2 orders of magnitude with respect to using the primitive or Gaussian basis sets alone. These results suggest our representation provides significant cost reductions while maintaining accuracy relative to molecular orbital or strictly diagonal approaches for modest-sized systems in both classical and quantum computation for correlated systems. △ Less

Submitted 30 August, 2019; originally announced September 2019.

Numerical and Theoretical Aspects of the DMRG-TCC Method Exemplified by the Nitrogen Dimer

Authors: Fabian M. Faulstich, Mihály Máté, Andre Laestadius, Mihály András Csirik, Libor Veis, Andrej Antalik, Jiří Brabec, Reinhold Schneider, Jiří Pittner, Simen Kvaal, Örs Legeza

Abstract: In this article, we investigate the numerical and theoretical aspects of the coupled-cluster method tailored by matrix-product states. We investigate chemical properties of the used method, such as energy size extensivity and the equivalence of linked and unlinked formulation. The existing mathematical analysis is here elaborated in a quantum chemical framework. In particular, we highlight the use… ▽ More In this article, we investigate the numerical and theoretical aspects of the coupled-cluster method tailored by matrix-product states. We investigate chemical properties of the used method, such as energy size extensivity and the equivalence of linked and unlinked formulation. The existing mathematical analysis is here elaborated in a quantum chemical framework. In particular, we highlight the use of a so-called CAS-ext gap describing the basis splitting between the complete active space and the external part. Moreover, the behavior of the energy error as a function of the optimal basis splitting is discussed. We show numerical investigations on the robustness with respect to the bond dimensions of the single orbital entropy and the mutual information, which are quantities that are used to choose the complete active space. Furthermore, we extend the mathematical analysis with a numerical study on the complete active space dependence of the error. △ Less

Submitted 20 September, 2018; originally announced September 2018.

The Coupled-Cluster Formalism - a Mathematical Perspective

Authors: Andre Laestadius, Fabian M. Faulstich

Abstract: The Coupled-Cluster theory is one of the most successful high precision methods used to solve the stationary Schrödinger equation. In this article, we address the mathematical foundation of this theory with focus on the advances made in the past decade. Rather than solely relying on spectral gap assumptions (non-degeneracy of the ground state), we highlight the importance of coercivity assumptions… ▽ More The Coupled-Cluster theory is one of the most successful high precision methods used to solve the stationary Schrödinger equation. In this article, we address the mathematical foundation of this theory with focus on the advances made in the past decade. Rather than solely relying on spectral gap assumptions (non-degeneracy of the ground state), we highlight the importance of coercivity assumptions - Gårding type inequalities - for the local uniqueness of the Coupled-Cluster solution. Based on local strong monotonicity, different sufficient conditions for a local unique solution are suggested. One of the criteria assumes the relative smallness of the total cluster amplitudes (after possibly removing the single amplitudes) compared to the Gårding constants. In the extended Coupled-Cluster theory the Lagrange multipliers are wave function parameters and, by means of the bivariational principle, we here derive a connection between the exact cluster amplitudes and the Lagrange multipliers. This relation might prove useful when determining the quality of a Coupled-Cluster solution. Furthermore, the use of an Aubin-Nitsche duality type method in different Coupled-Cluster approaches is discussed and contrasted with the bivariational principle. △ Less

Submitted 7 December, 2018; v1 submitted 22 April, 2018; originally announced April 2018.

Analysis of The Tailored Coupled-Cluster Method in Quantum Chemistry

Authors: Fabian M. Faulstich, Andre Laestadius, Örs Legeza, Reinhold Schneider, Simen Kvaal

Abstract: In quantum chemistry, one of the most important challenges is the static correlation problem when solving the electronic Schrödinger equation for molecules in the Born--Oppenheimer approximation. In this article, we analyze the tailored coupled-cluster method (TCC), one particular and promising method for treating molecular electronic-structure problems with static correlation. The TCC method comb… ▽ More In quantum chemistry, one of the most important challenges is the static correlation problem when solving the electronic Schrödinger equation for molecules in the Born--Oppenheimer approximation. In this article, we analyze the tailored coupled-cluster method (TCC), one particular and promising method for treating molecular electronic-structure problems with static correlation. The TCC method combines the single-reference coupled-cluster (CC) approach with an approximate reference calculation in a subspace [complete active space (CAS)] of the considered Hilbert space that covers the static correlation. A one-particle spectral gap assumption is introduced, separating the CAS from the remaining Hilbert space. This replaces the nonexisting or nearly nonexisting gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital usually encountered in standard single-reference quantum chemistry. The analysis covers, in particular, CC methods tailored by tensor-network states (TNS-TCC methods). The problem is formulated in a nonlinear functional analysis framework, and, under certain conditions such as the aforementioned gap, local uniqueness and existence are proved using Zarantonello's lemma. From the Aubin--Nitsche-duality method, a quadratic error bound valid for TNS-TCC methods is derived, e.g., for linear-tensor-network TCC schemes using the density matrix renormalization group method. △ Less

Submitted 13 October, 2019; v1 submitted 15 February, 2018; originally announced February 2018.

Unraveling Mirror Properties in Time-Delayed Quantum Feedback Scenarios

Authors: Fabian M. Faulstich, Manuel Kraft, Alexander Carmele

Abstract: We derive in the Heisenberg picture a widely used phenomenological coupling element to treat feedback effects in quantum optical platforms. Our derivation is based on a microscopic Hamiltonian, which describes the mirror-emitter dynamics based on a dielectric, a mediating fully quantized electromagnetic field, and a single two-level system in front of the dielectric. The dielectric is modeled as a… ▽ More We derive in the Heisenberg picture a widely used phenomenological coupling element to treat feedback effects in quantum optical platforms. Our derivation is based on a microscopic Hamiltonian, which describes the mirror-emitter dynamics based on a dielectric, a mediating fully quantized electromagnetic field, and a single two-level system in front of the dielectric. The dielectric is modeled as a a system of identical two-state atoms. The Heisenberg equation yields a system of describing differential operator equations, which we solve in the Weisskopf-Wigner limit. Due to a finite round-trip time between emitter and dielectric, we yield delay differential operator equations. Our derivation motivates and justifies the typical phenomenological assumed coupling element and allows, furthermore, a generalization to a variety of mirrors, such as dissipative mirrors or mirrors with gain dynamics. △ Less

Submitted 7 September, 2017; v1 submitted 17 March, 2017; originally announced March 2017.