Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Cambridge University Press eBooks

Cluster-expansion representation of the quantum dynamics

M. Kira, S. W. Koch

Published Nov 17, 2011• 0 citations
Publisher

Abstract

In this book, we encounter the hierarchy problem when we apply the equation-of-motion technique to analyze the quantum dynamics of the coupled light–matter system. To obtain systematic approximations, we use the so-called cluster-expansion method where many-body quantities are systematically grouped into cluster classes based on how important they are to the overall quantum dynamics. With increasing complexity, the clusters contain independent single particles (singlets), interacting pairs (doublets), three-particle terms (triplets), coupled four-particle contributions (quadruplets), as well as higher-order correlations. In this context, the N -particle concept is somewhat formal because it refers to generic N -particle expectation values 〈 N 〉, which may consist of an arbitrary mixture of carrier, photon, and phonon operators, as discussed in Section 14.2. In order to truncate the hierarchy problem at a given level, 〈 N 〉 is approximated through a functional structure that includes all clusters up to the predetermined level while all remaining clusters with a higher rank are omitted. It is natural that the corresponding approximations can be systematically improved by increasing the number of clusters included. To the best of our knowledge, the idea of coupled-clusters approaches was first formulated by Fritz Coester and Hermann Kümmel in the 1950s to describe nuclear many-body phenomena. The approach was then modified for the needs of quantum chemistry by Jiri Cizek 1966 to deal with many-body phenomena in atoms and molecules. Currently, it is one of the most accurate methods to compute molecular eigenstates.

DOI
10.1017/cbo9781139016926.016

Related papers

Coupled Cluster Method

Rajat Kumar Chaudhuri, Sudip Kumar Chattopadhyay

Many-body perturbation theories (MBPT) (both single and multireference), by design, are capable of representing electron correlation very well. However, the practical implementation of these methods for many-electron systems beyond a certain order is extremely difficult, if not impossible. For instance, there are hardly any calculations beyond third order at the multireference and fifth order at the single reference of perturbation theories. Thus, these methods are typically limited to systems where the higher-order electron correlation effect is not significant. A practical solution to this problem is the coupled cluster theory, which was first introduced by Coester [20] and Coester and Kümmel [21] in nuclear physics and was later adopted in the electron structure theory by Čížek [23, 106] and Paldus and Čížek [107]. The coupled-cluster (CC) theory is a non-perturbative all-order theory. This stateof-the-art theory has now been recognized as the basis of modern quantitative quantum chemistry and atomic physics for high-precision calculations. Its superior accuracy is rooted in the exponential ansatz for the wave function, which provides an energy expression that scales correctly with the number of electrons and under certain provisions leads to a faster convergence of the many-body expansion.

CRC Press eBooks•Feb 17, 2017•0 citations

General aspects of semiconductor optics

M. Kira, S. W. Koch

The optically generated excitations in semiconductors constitute a genuine many-body system. To describe its quantum-optical features, we have to expand significantly the theoretical models used so far. However, the important insights of Chapters 16–23 are already presented in a form in which most of them can directly be used and generalized to analyze central properties of the optical excitations in solids. As for atoms, the optical transitions in semiconductors are induced via dipole interaction between photons and electrons. We can thus efficiently construct a systematic quantum-optical theory for semiconductors by following the cluster-expansion approach. One of the main differences from atoms is that the electronic excitations in semiconductors form a strongly interacting many-body system. Thus, we must systematically treat the arising Coulomb-induced hierarchy problem together with the quantum-optical one. Moreover, the coupling of electrons to lattice vibrations, i.e., the phonons, produces yet another hierarchy problem. In addition, in solid-state spectroscopy one often uses multimode light fields such that one cannot rely on the single-mode simplifications to study semiconductor quantum optics. As shown in Chapter 15, the Coulomb-, phonon-, and photon-induced hierarchy problems have formally an identical structure. Thus, we start the analysis by investigating how semiconductor quantum optics emerges from the dynamics of correlated clusters. We first focus on the basic properties of the optical transitions in the classical regime. This means investigating the fundamental optical phenomena resulting from the singlets. The full singlet–doublet approach is presented in Chapters 28–30.

Cambridge University Press eBooks•Nov 17, 2011•0 citations

Introduction

Isaiah Shavitt, Rodney J. Bartlett

Scope The book focuses primarily on many-body (or better, many-electron) methods for electron correlation. These include Rayleigh–Schrödinger perturbation theory (RSPT) , particularly in its diagrammatic representation (referred to as many-body perturbation theory , or MBPT ), and coupled-cluster (CC) theory; their relationship to configuration interaction (CI) is included. Further extensions address properties other than the energy, and also excited states and multireference CC and MBPT methods. The many-body algebraic and diagrammatic methods used in electronic structure theory have their origin in quantum field theory and in the study of nuclear matter and nuclear structure. The second-quantization formalism was first introduced in a treatment of quantized fields by Dirac (1927) and was extended to fermion systems by Jordan and Klein (1927) and by Jordan and Wigner (1928). This formalism is particularly useful in field theory, in scattering problems and in the study of infinite systems because it easily handles problems involving infinite, indefinite or variable numbers of particles. The diagrammatic approach was introduced into field theory by Feynman (1949a,b) and applied to many-body systems by Hugenholtz (1957) and by Goldstone (1957). Many-body perturbation theory and its linked-diagram formalism were first introduced by Brueckner and Levinson (1955) and by Brueckner (1955), and were formalized by Goldstone (1957). Other important contributions to the methodology, first in field theory and then in the theory of nuclear structure, are due to Dyson (1949a,b), Wick (1950), Hubbard (1957, 1958a,b) and Frantz and Mills (1960). Applications to the electronic structure of atoms and molecules began with the work of Kelly (1963, 1964a,b, 1968), and molecular applications using finite analytical basis sets appeared in the work of Bartlett and Silver (1974a, b).

Cambridge University Press eBooks•Aug 6, 2009•26 citations

Method of Correlated Basis Functions

J. W. Clark, Paul A. Westhaus

The method of correlated basis functions provides a means of systematically exploiting our intuitive understanding of real physical systems. Here we seek especially to develop a formalism for describing systems of strongly interacting fermions, for example, nuclei and liquid ${\mathrm{He}}^{3}$, using a complete set of correlated functions ${\ensuremath{\Psi}}_{\mathrm{m}}=\frac{F{\ensuremath{\Phi}}_{\mathrm{m}}}{{({\ensuremath{\Phi}}_{\mathrm{m}}, {F}^{2}{\ensuremath{\Phi}}_{\mathrm{m}})}^{\frac{1}{2}}}$. $F$ is taken as a product of two-body factors, $\ensuremath{\Pi}{i<j}^{}[f({r}_{\mathrm{ij}})]$, the $f({r}_{\mathrm{ij}})$ embodying at least the correlations arising from strong short-range repulsions; and ${{\ensuremath{\Phi}}_{\mathrm{m}}}$, a set of independent-particle-model functions which might serve as a basis in the absence of highly singular interactions. The matrix elements ${\mathfrak{H}}_{\mathrm{mn}}$, ${\mathfrak{N}}_{\mathrm{mn}}$ of the Hamiltonian and identity operators with respect to the correlated basis ${{\ensuremath{\Psi}}_{\mathrm{m}}}$ are evaluated by generalization of cluster-expansion techniques developed originally for the calculation of the partition function of a classical imperfect gas and adapted to the quantum-mechanical many-body problem by Iwamoto and Yamada. Contributions to our expansions for the ${\mathfrak{H}}_{\mathrm{mn}}$ are classified according to order of magnitude in the quantity $\ensuremath{\rho}\ensuremath{\omega}$, where $\ensuremath{\omega}=\ensuremath{\int}({f}^{2}(r)\ensuremath{-}1)d\mathrm{r}$ is a "correlation parameter" and $\ensuremath{\rho}$ the average particle density. To the extent that $|\ensuremath{\rho}\ensuremath{\omega}|$ is small, it provides a true expansion parameter. A detailed technical discussion of the cluster expansions is given, with special emphasis on their symmetry properties upon truncation. Although the many-body problem is soluble by simultaneous diagonalization of the matrices (${\mathfrak{H}}_{\mathrm{mn}}$) and (${\mathfrak{N}}_{\mathrm{mn}}$), it is preferable first to transform to an orthonormal basis of "unperturbed" correlated functions ${{\ensuremath{\Theta}}_{\mathrm{m}}}$ with respect to which the Hamiltonian and identity operators have the matrix representations (${H}_{\mathrm{mn}}$) and $({N}_{\mathrm{mn}})=({\ensuremath{\delta}}_{\mathrm{mn}})$, and then to attempt the separate diagonalization of (${H}_{\mathrm{mn}}$). In general, such a transformation is expedited by the L\"owdin procedure carried out to the desired order in $\ensuremath{\rho}\ensuremath{\omega}$. However, a Schmidt orthogonalization scheme might prove more advantageous in specific cases. Field-theoretic techniques (of both perturbative and nonperturbative nature) may be brought to bear in the diagonalization of (${H}_{\mathrm{mn}}$). The inverse of the unitary transformation which takes the basis ${{\ensuremath{\Phi}}_{\mathrm{m}}}$ into the basis ${{\ensuremath{\Theta}}_{\mathrm{m}}}$, applied to the given singular Hamiltonian, leads---for any $F$ in the class considered---to a non-singular "effective Hamiltonian." Following this formal cue, we set out to construct explicitly a second-quantized effective Hamiltonian which generates the same matrix elements, to prescribed order in $\ensuremath{\rho}\ensuremath{\omega}$, when operating between independent-particle kets of the model occupation-number representation, as does the given Hamiltonian in the configuration-space representation, between the corresponding correlated functions of the L\"owdin basis. The quasiparticles described by this effective Hamiltonian interact via well-behaved two-body, three-body, \ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}, $Q$-body, \ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{} potentials, in contrast to the highly singular two-body interactions of the real particles. For small $|\ensuremath{\rho}\ensuremath{\omega}|$, the effective two-body forces predominate. An avenue to deeper understanding of the applicability of the nuclear-shell model and the dressed-quasiparticle theory of liquid ${\mathrm{He}}^{3}$ has been opened. This article treats only the formalism for a uniform extended medium. Numerical calculations, as well as the extension of the method to finite systems, will follow.

Physical Review•Jan 21, 1966•160 citations

Many-Body Problem for Strongly Interacting Particles. II. Linked Cluster Expansion

K. A. Brueckner

An approximation method developed previously to deal with many particles in strong interaction is examined in further detail. It is shown that the series giving the interaction energy is a development in a sequence of linked or irreducible cluster terms each of which gives a contribution to the energy proportional to the total number of particles. Consequently the convergence of the expansion is independent of the total number of particles. The origin of this simple feature is illustrated by showing that a similar situation exists in the expansion of standard perturbation theory. The numerical convergence of the expansion is quantitatively discussed for the nuclear problem where it is shown that the correction arising from the first cluster term involving three particles is less than the leading term by a factor of about ${10}^{\ensuremath{-}4}$. The smallness of the correction is largely a result of the action of the exclusion principle.

Physical Review•Oct 1, 1955•770 citations

Calculation of properties with the coupled-cluster method

Hendrik J. Monkhorst

The cluster-expansion approach to the correlation problem, pioneered by Cocster, Kümmel, Cizek and Paldus, is extended to calculation of static and dynamic properties of many-fermion systems. Linear, inhomogeneous equations are obtained for properties of any order. A time-dependent formulation gives frequency-dependent properties, yielding excitation energies, transition probabilities, and (possibly) life times reminiscent of Green's function methods.

International Journal of Quantum Chemistry•Jun 18, 2009•349 citations
In this book, we encounter the hierarchy problem when we apply the equation-of-motion technique to analyze the quantum dynamics of the coupled light–matter system. To obtain systematic approximations, we use the so-called cluster-expansion method where many-body quantities are systematically grouped into cluster classes based on how important they are to the overall quantum dynamics. With increasing complexity, the clusters contain independent single particles (singlets), interacting pairs (doublets), three-particle terms (triplets), coupled four-particle contributions (quadruplets), as well as higher-order correlations. In this context, the N -particle concept is somewhat formal because it refers to generic N -particle expectation values 〈 N 〉, which may consist of an arbitrary mixture of carrier, photon, and phonon operators, as discussed in Section 14.2. In order to truncate the hierarchy problem at a given level, 〈 N 〉 is approximated through a functional structure that includes all clusters up to the predetermined level while all remaining clusters with a higher rank are omitted. It is natural that the corresponding approximations can be systematically improved by increasing the number of clusters included. To the best of our knowledge, the idea of coupled-clusters approaches was first formulated by Fritz Coester and Hermann Kümmel in the 1950s to describe nuclear many-body phenomena. The approach was then modified for the needs of quantum chemistry by Jiri Cizek 1966 to deal with many-body phenomena in atoms and molecules. Currently, it is one of the most accurate methods to compute molecular eigenstates.
Many-body perturbation theories (MBPT) (both single and multireference), by design, are capable of representing electron correlation very well. However, the practical implementation of these methods for many-electron systems beyond a certain order is extremely difficult, if not impossible. For instance, there are hardly any calculations beyond third order at the multireference and fifth order at the single reference of perturbation theories. Thus, these methods are typically limited to systems where the higher-order electron correlation effect is not significant. A practical solution to this problem is the coupled cluster theory, which was first introduced by Coester [20] and Coester and Kümmel [21] in nuclear physics and was later adopted in the electron structure theory by Čížek [23, 106] and Paldus and Čížek [107]. The coupled-cluster (CC) theory is a non-perturbative all-order theory. This stateof-the-art theory has now been recognized as the basis of modern quantitative quantum chemistry and atomic physics for high-precision calculations. Its superior accuracy is rooted in the exponential ansatz for the wave function, which provides an energy expression that scales correctly with the number of electrons and under certain provisions leads to a faster convergence of the many-body expansion.
The optically generated excitations in semiconductors constitute a genuine many-body system. To describe its quantum-optical features, we have to expand significantly the theoretical models used so far. However, the important insights of Chapters 16–23 are already presented in a form in which most of them can directly be used and generalized to analyze central properties of the optical excitations in solids. As for atoms, the optical transitions in semiconductors are induced via dipole interaction between photons and electrons. We can thus efficiently construct a systematic quantum-optical theory for semiconductors by following the cluster-expansion approach. One of the main differences from atoms is that the electronic excitations in semiconductors form a strongly interacting many-body system. Thus, we must systematically treat the arising Coulomb-induced hierarchy problem together with the quantum-optical one. Moreover, the coupling of electrons to lattice vibrations, i.e., the phonons, produces yet another hierarchy problem. In addition, in solid-state spectroscopy one often uses multimode light fields such that one cannot rely on the single-mode simplifications to study semiconductor quantum optics. As shown in Chapter 15, the Coulomb-, phonon-, and photon-induced hierarchy problems have formally an identical structure. Thus, we start the analysis by investigating how semiconductor quantum optics emerges from the dynamics of correlated clusters. We first focus on the basic properties of the optical transitions in the classical regime. This means investigating the fundamental optical phenomena resulting from the singlets. The full singlet–doublet approach is presented in Chapters 28–30.
Scope The book focuses primarily on many-body (or better, many-electron) methods for electron correlation. These include Rayleigh–Schrödinger perturbation theory (RSPT) , particularly in its diagrammatic representation (referred to as many-body perturbation theory , or MBPT ), and coupled-cluster (CC) theory; their relationship to configuration interaction (CI) is included. Further extensions address properties other than the energy, and also excited states and multireference CC and MBPT methods. The many-body algebraic and diagrammatic methods used in electronic structure theory have their origin in quantum field theory and in the study of nuclear matter and nuclear structure. The second-quantization formalism was first introduced in a treatment of quantized fields by Dirac (1927) and was extended to fermion systems by Jordan and Klein (1927) and by Jordan and Wigner (1928). This formalism is particularly useful in field theory, in scattering problems and in the study of infinite systems because it easily handles problems involving infinite, indefinite or variable numbers of particles. The diagrammatic approach was introduced into field theory by Feynman (1949a,b) and applied to many-body systems by Hugenholtz (1957) and by Goldstone (1957). Many-body perturbation theory and its linked-diagram formalism were first introduced by Brueckner and Levinson (1955) and by Brueckner (1955), and were formalized by Goldstone (1957). Other important contributions to the methodology, first in field theory and then in the theory of nuclear structure, are due to Dyson (1949a,b), Wick (1950), Hubbard (1957, 1958a,b) and Frantz and Mills (1960). Applications to the electronic structure of atoms and molecules began with the work of Kelly (1963, 1964a,b, 1968), and molecular applications using finite analytical basis sets appeared in the work of Bartlett and Silver (1974a, b).
The method of correlated basis functions provides a means of systematically exploiting our intuitive understanding of real physical systems. Here we seek especially to develop a formalism for describing systems of strongly interacting fermions, for example, nuclei and liquid He3{\mathrm{He}}^{3}He3, using a complete set of correlated functions ${\ensuremath{\Psi}}_{\mathrm{m}}=\frac{F{\ensuremath{\Phi}}_{\mathrm{m}}}{{({\ensuremath{\Phi}}_{\mathrm{m}}, {F}^{2}{\ensuremath{\Phi}}_{\mathrm{m}})}^{\frac{1}{2}}}$. FFF is taken as a product of two-body factors, $\ensuremath{\Pi}{i<j}^{}[f({r}_{\mathrm{ij}})]$, the f(rij)f({r}_{\mathrm{ij}})f(rij​) embodying at least the correlations arising from strong short-range repulsions; and ${{\ensuremath{\Phi}}_{\mathrm{m}}}$, a set of independent-particle-model functions which might serve as a basis in the absence of highly singular interactions. The matrix elements Hmn{\mathfrak{H}}_{\mathrm{mn}}Hmn​, Nmn{\mathfrak{N}}_{\mathrm{mn}}Nmn​ of the Hamiltonian and identity operators with respect to the correlated basis ${{\ensuremath{\Psi}}_{\mathrm{m}}}$ are evaluated by generalization of cluster-expansion techniques developed originally for the calculation of the partition function of a classical imperfect gas and adapted to the quantum-mechanical many-body problem by Iwamoto and Yamada. Contributions to our expansions for the Hmn{\mathfrak{H}}_{\mathrm{mn}}Hmn​ are classified according to order of magnitude in the quantity $\ensuremath{\rho}\ensuremath{\omega}$, where $\ensuremath{\omega}=\ensuremath{\int}({f}^{2}(r)\ensuremath{-}1)d\mathrm{r}$ is a "correlation parameter" and $\ensuremath{\rho}$ the average particle density. To the extent that $|\ensuremath{\rho}\ensuremath{\omega}|$ is small, it provides a true expansion parameter. A detailed technical discussion of the cluster expansions is given, with special emphasis on their symmetry properties upon truncation. Although the many-body problem is soluble by simultaneous diagonalization of the matrices (Hmn{\mathfrak{H}}_{\mathrm{mn}}Hmn​) and (Nmn{\mathfrak{N}}_{\mathrm{mn}}Nmn​), it is preferable first to transform to an orthonormal basis of "unperturbed" correlated functions ${{\ensuremath{\Theta}}_{\mathrm{m}}}$ with respect to which the Hamiltonian and identity operators have the matrix representations (Hmn{H}_{\mathrm{mn}}Hmn​) and $({N}_{\mathrm{mn}})=({\ensuremath{\delta}}_{\mathrm{mn}})$, and then to attempt the separate diagonalization of (Hmn{H}_{\mathrm{mn}}Hmn​). In general, such a transformation is expedited by the L\"owdin procedure carried out to the desired order in $\ensuremath{\rho}\ensuremath{\omega}$. However, a Schmidt orthogonalization scheme might prove more advantageous in specific cases. Field-theoretic techniques (of both perturbative and nonperturbative nature) may be brought to bear in the diagonalization of (Hmn{H}_{\mathrm{mn}}Hmn​). The inverse of the unitary transformation which takes the basis ${{\ensuremath{\Phi}}_{\mathrm{m}}}$ into the basis ${{\ensuremath{\Theta}}_{\mathrm{m}}}$, applied to the given singular Hamiltonian, leads---for any FFF in the class considered---to a non-singular "effective Hamiltonian." Following this formal cue, we set out to construct explicitly a second-quantized effective Hamiltonian which generates the same matrix elements, to prescribed order in $\ensuremath{\rho}\ensuremath{\omega}$, when operating between independent-particle kets of the model occupation-number representation, as does the given Hamiltonian in the configuration-space representation, between the corresponding correlated functions of the L\"owdin basis. The quasiparticles described by this effective Hamiltonian interact via well-behaved two-body, three-body, \ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}, QQQ-body, \ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{} potentials, in contrast to the highly singular two-body interactions of the real particles. For small $|\ensuremath{\rho}\ensuremath{\omega}|$, the effective two-body forces predominate. An avenue to deeper understanding of the applicability of the nuclear-shell model and the dressed-quasiparticle theory of liquid He3{\mathrm{He}}^{3}He3 has been opened. This article treats only the formalism for a uniform extended medium. Numerical calculations, as well as the extension of the method to finite systems, will follow.
An approximation method developed previously to deal with many particles in strong interaction is examined in further detail. It is shown that the series giving the interaction energy is a development in a sequence of linked or irreducible cluster terms each of which gives a contribution to the energy proportional to the total number of particles. Consequently the convergence of the expansion is independent of the total number of particles. The origin of this simple feature is illustrated by showing that a similar situation exists in the expansion of standard perturbation theory. The numerical convergence of the expansion is quantitatively discussed for the nuclear problem where it is shown that the correction arising from the first cluster term involving three particles is less than the leading term by a factor of about ${10}^{\ensuremath{-}4}$. The smallness of the correction is largely a result of the action of the exclusion principle.