Non-Hermitian coupled cluster method for non-stationary systems and its interaction-picture reinterpretation

The interaction picture in a non-Hermitian realization is discussed in detail and considered for its practical use in many-body quantum physics. The resulting non-Hermitian interaction-picture (NHIP) description of dynamics, in which both the wave functions and operators belonging to physical observables cease to remain constant in time, is a non-Hermitian generalization of the traditional Dirac picture of standard quantum mechanics, which itself is widely used in quantum ﬁeld theory calculations. Particular attention is paid here to the variational (or, better, bivariational) and dynamical (i.e., non-stationary) aspects that are characteristic of the coupled cluster method (CCM) techniques that nowadays form one of the most versatile and most accurate of all available formulations of quantum many-body theory. In so doing we expose and exploit multiple parallels between the NHIP and the CCM in its time-dependent versions.


Introduction
In a wide variety of branches of applied quantum mechanics, one often needs to know, with a reasonable numerical precision, the value and evolution of the time-dependent wave function ψ(t). Not surprisingly, the construction of this function is relatively straightforward only for a few exceptional, not too complicated (and hence also, generically, not very realistic) selfadjoint Hamiltonians h in the Schrödinger equation (in units whereh = 1 throughout), The task usually requires the use of one of the available sophisticated (e.g., typically, perturbative or variational) numerical methods. In the present paper we recall two of the apparently very different strategies of the variational class, with the express intention of showing that behind the apparent formal differences one can also find multiple parallels that can lead to new ideas and novel constructive approaches. The present study has its roots in our earlier investigation [1] in which we succeeded in comparing some alternative methods of solving Eq. (1) in the stationary regime, in each of which the construction remains reducible to the diagonalization of the Hamiltonian h. In particular we revealed a close structural parallelism between the powerful and versatile coupled cluster method (CCM) approach to the diagonalization of operator h (and see, e.g., Refs. [2][3][4][5][6] for details of the formalism) and the successful interacting boson model (IBM) [7][8][9]. The CCM is, by now, well known and very widely used in such fields as atomic and molecular physics, nuclear physics, quantum chemistry (in which it forms the "gold standard" for accuracy), and many branches of condensed matter physics, while the IBM offers an efficient tool for the calculation of the low-lying spectra of energies E of the heaviest stable atomic nuclei [7][8][9].
In the present paper our aim is to demonstrate that the above-mentioned parallelism finds its further natural extension to the domain of dynamics in which the quantum systems in question are non-stationary. We shall first recall the known CCM ↔ IBM parallels [1], and then show how to generalize them. We do so along lines that are guided by the expectation of a mutual enrichment. Explicitly, for the numerical constructions of energies E and wave functions ψ(t), we pay specific attention here to two specific alternative strategies, viz. the various versions of the CCM approach (see, e.g., Ref. [10] for a compact introductory review) and to standard (or "textbook") quantum theory in its so-called three-Hilbert-space formulation [11,12]. For a comprehensive introduction to some of the most successful phenomenological applications of the three-Hilbert-space formalism the reader is also referred to Refs. [13,14].
We anticipate that the future impact of our considerations might range from a purely mathematical guarantee of the computational feasibility of model-building, up to a deeply physical reinterpretation of relationships between alternative formulations of many-body quantum mechanics. We will emphasize certain emergent features of connections between different means of description and between the simplifying assumptions and dynamics-simulating techniques.

The fundamentals of the CCM approach
None of the existing CCM-based calculations leaves the traditional theoretical framework of quantum mechanics. The correspondence between the IBM and CCM formulations of quantum mechanics may be one-to-one. Indeed, both formulations employ the factorized Ansatz in which the auxiliary (and, generally, non-unitary) operator Ω is intended to serve as the source of the quantum correlations which are unaccounted for in the original approximation ψ 0 . One of the most natural mathematical requirements imposed upon the operator Ω is that it should be invertible. Thus, the element ψ of an initial physical Hilbert space H (initial) is treated as an image of a, presumably, perceivably simpler reference state ψ 0 , which may itself belong, in principle, to another, more user-friendly, auxiliary Hilbert space H (user−friendly) . The original Hamiltonian h of Eq. (1) belonging to H (initial) is similarly replaced by its isospectral (and, in general, non-Hermitian) avatar, that belongs to the auxiliary Hilbert space H (user−friendly) . In the CCM setting, the auxiliary operator Ω in the Ansatz of Eq. (2) takes the very specific exponentiated form, This characteristic feature of the method henceforth ensures that the Goldstone linked cluster theorem [15] is automatically satisfied for the system under study at any level of approximation for the operator S, as we explain in more detail below. In turn, this guarantees the very important property for all many-body systems that at all such levels the system is both size-extensive and size-consistent [10], where size-extensivity is the property that the leading term in the energy of an N -particle self-bound system scales linearly with N as N → ∞, and where size-consistency implies that a many-body wave function dissociates correctly into non-interacting fragments under infinite separation of the fragments. The original version of the CCM, as invented independently by Coester and Kümmel [16,17] andČížek [18,19], is nowadays referred to as the normal CCM (NCCM), in order to distinguish it from the later extended (ECCM) version introduced by Arponen [2]. For present purposes most of what we discuss here is equally relevant to both versions, although we mostly have the NCCM in mind for specific applications.

The normal coupled cluster method (NCCM)
The exact ket and bra wave functions for an arbitrary quantum many-body system are defined in H (initial) to be |ψ(t) and respectively, at an arbitrary time t. We assume that the Hilbert space H (initial) for the system may be described in terms of a normalized, stationary reference (or model) state |ψ 0 (i.e., with ψ 0 |ψ 0 = 1), just as in Eq. (2). Within the CCM we further assume that the reference state acts as a cyclic vector for a corresponding set of mutually commuting multi-configurational creation operators {C + I }, such that |ψ 0 acts as a generalized vacuum state with respect to them, in a notation where we define C − I ≡ (C + I ) † and C + 0 ≡ 1, the unit vector. The reference state |ψ 0 must also be chosen to be non-orthogonal to the actual wave function |ψ(t) of the system, ψ 0 |ψ(t) = 0 ; ∀t .
The index I is a set index and, in general, the multi-configurational creation operator C + I comprises a product of single-particle operators, as we illustrate in Sect. 2.3 for a specific example. The set {I } is complete in the usual sense that the set of states {C + I |ψ 0 } provides 123 a complete basis for the ket Hilbert space. It is also convenient to choose the basis to be orthonormalized, such that where δ I,J is an appropriately defined Kronecker symbol. In this case we have that the identity operator in the Hilbert space may be resolved as The ket and bra states of the many-body system are now formally parametrized independently in the NCCM as follows, in terms of the c-number function k(t), and the NCCM time-dependent correlation operators S(t) andS(t), which are themselves now decomposed as follows, The scale factor k(t) is needed for the temporally evolving intermediate normalization (i.e., ψ 0 |ψ(t) = e k(t) ) but is irrelevant henceforth since it cancels from any expectation value [and the normalization condition ψ (t)|ψ(t) = 1 from Eq. (5) is now preserved for all times t] by the specific parametrizations of Eqs. (11) and (12). Since S(t) andS(t) are henceforth treated as independent operators, the Hermiticity relation between corresponding ket and bra states may be violated when subsequent approximations are made, typically by truncating the complete set of configurations {I } in the sums of Eq. (12). This shortcoming is far outweighed by the fact that the CCM parametrizations of Eqs. (11) and (12) always exactly satisfy the very important Hellmann-Feynman theorem [20,21] at all such levels of approximation [10]. It is precisely this feature that guarantees the robustness and accuracy of numerical results obtained within the CCM framework. The expectation value Q(t) of an arbitrary physical observable described by the Hermitian operator q(t) in H (initial) , which may itself have an intrinsic time dependence, may thus be expressed within the NCCM wholly in terms of the set of c-number correlation coefficients {s I (t),s I (t)} as follows, Such an expectation functional is evaluated by employing the nested commutation expansion for the similarity transform Here, the n-fold nested commutators [q, S] n are defined iteratively, The similarity-transformed Hamiltonian, H (t) ≡ e −S(t) he S(t) , lies at the heart of the NCCM. Its parametrizations of Eqs. (11) and (12) are specifically chosen so that the otherwise infinite sum in Eq. (14b) will actually terminate exactly at a low finite order for any operator q all of whose terms involve only a product of a finite number of single-particle operators, such as is almost always the case for the operators of interest. The reason for this termination lies simply in the twin facts that all components of the operator S(t) in the expansion of Eq. (12) commute with one another from Eq. (6), and the basic single-particle operators form a Lie algebra, together with the model state |ψ 0 being a vacuum state for the operator set {C + I }, as in Eq. (7) (and see, e.g., Refs. [6,10] for further details).
The above same key feature of the NCCM is also responsible for all terms in the expansion of Q(t) from Eq. (13) for any physical operator q being linked . Thus, we may write the NCCM expectation value functional Q(t) in the schematic form, J |q(t)|K 1 , · · ·, K n NCCMsJ s K n · · · s K 1 , (16) where, by comparison with Eqs. (14b) and (15), the NCCM matrix elements may be written as The suffix L in Eq. (17) indicates the NCCM linked structure implied by the nested commutator sum, viz. that from each group of particles characterized by each of the configuration-space indices K j , the amplitude s K j associated with the cluster creation operator C + K j , must have at least one particle line connected to the operator q. After expanding Eq. (13) in powers of S by the use of Eq. (14b), only such terms are retained where at least one link (or contraction, due to a non-vanishing commutator term) exists between each S operator and the operator q. The sum over n in Eq. (16) will extend again to a finite limit depending on the operator q.
It is precisely this linkedness feature (specifically for the energy expectation value functional H (t)), as expressed in Eqs. (16) and (17), that implies immediately that the NCCM automatically satisfies exactly the Goldstone linked-cluster theorem [15] at any level of truncation in the expansions of Eq. (12) for the correlation operators. In turn this also immediately guarantees the size-extensivity of the method at all such levels of approximation. Thus, in the solution of the explicit evolution equations that we will now derive for the basic NCCM amplitudes {s I (t),s ( t)}, the only approximation that is ever made in any CCM calculation is to decide which configurations {I } are to be retained in the expansions of Eq. (12).
The NCCM coefficients {s I (t),s I (t)} are obtained by introducing the action functional A, The (bivariational) stationarity principle, with respect to all independent variations in the bra and ket states, subject to the conditions δ ψ (t i )| = 0 = δ|ψ(t i ) where i = 0, 1, is seen to be equivalent to the time-dependent Schrödinger equations, 123 analogous to Eq. (1). Insertion of the NCCM parametrizations of Eqs. (11) and (12) into Eq. (18) yields the result, The independent (bivariational) stationarity conditions, ∂A/∂s I = 0 = ∂A/∂s I , then yield the evolution equations for the NCCM cluster correlation coefficients, in terms of the energy expectation functional H , given as in Eq. (13) with the replacement q(t) → h. We recognize these equations as (being closely related to) the classical Hamilton's equations of motion for the NCCM cluster correlation coefficients. Thus, the coefficients s I ands I form canonically conjugate pairs in the usual sense of classical mechanics.
We can now also easily derive the equation of motion for the expectation value Q(t) ≡ Q[s I ,s I ; t] of an arbitrary operator q(t) in H (initial) given by Eq. (13). Using Eq. (22), we rewrite the usual chain-rule relation, in the suggestive form, in terms of a generalized classical Poisson bracket , defined between the expectation values A and B of two arbitrary operators in H (user−friendly) as follows, Equation (24) is just the well-known classical equation of motion in the canonical formalism. It is clearly the classical, exactly mapped , counterpart of the usual quantum-mechanical Heisenberg equation of motion for the operator q(t) in the original Hilbert space H (initial) , with the lower-case operator q(t) and its upper-case NCCM similarity-transformed avatar operator Q(t) related by Eq. (14a) The above results can be made even more suggestive by choosing new basic variables to be the generalized c-number fields φ I and their canonically conjugate c-number generalized momentum densities π I , which are defined as follows, The inverse transformations are thus simply given by It is now easy to show that Eq. (25) in the new variables takes the following form, The equations of motion (22) lead to their equivalent counterparts for the new variables, The (completely classical) phase space {φ I , π I } has the canonical symplectic structure, In order to complete this description of an exact classicization of an arbitrary quantum many-body theory, we note that the expectation value of the commutator between an arbitrary pair of operators a and b in the original Hilbert space H (initial) is also exactly mapped into its corresponding Poisson bracket of Eq. (25) or Eq. (28). We may use Eq. (13) to write the expectation value of the product ab of the two operators as and hence, using Eq. (14a), as By inserting a resolution of the identity of the form of Eq. (10) between the operators A(t) and B(t) in the middle expression in the above equation, we find The term with I = 0 in the sum in the first term of the above equation is identically zero, since C + 0 ≡ 1. Next, if we use the explicit form for the expectation value of an arbitrary operator given by Eq. (13), together with the NCCM forms for the operators S(t) andS(t) given by Eq. (12), we finally derive the following relations for the remaining terms, 123 Substitution of Eqs. (34a) and (34b) into the second line of Eq. (33) then yields the following result, where in the third line of Eq. (33) we have inserted a resolution of the identity operator of the form of Eq. (10) between the operators C + I and A(t). The last term in Eq. (35) is symmetric under the interchange A B, as [C + I , C + J ] = 0 from Eq. (6), and hence we arrive at the desired relation, It is very important to note that the consistency of the classical and quantum formulations induced by the NCCM can thus be understood as a manifestation of the correspondence principle (or Ehrenfest theorem) in a suitably generalized form.
In partial summary, we have thus been led to a key result of the NCCM formalism, namely that our (rather arbitrarily specified) quantum many-body problem in a Hilbert space has been exactly mapped onto the classical Hamiltonian mechanics for the complex c-number amplitudes {s I ,s I }, or equivalently the canonical fields {φ I } and their conjugate momentum densities {π I }, in a classical NCCM phase space. These fields are themselves defined as many-body amplitudes in the configuration space labeled by the indices {I } that describe the particular subsets (or clusters) of particles under consideration. In the modern terminology of classical mechanics this phase space is a symplectic differentiable manifold [22]. The differentiability is just a consequence of the fact that the c-number parameters {s I ,s I } or, equivalently, {φ I , π I }, are continuous complex numbers, with respect to which one is allowed to take (possibly functional) derivatives. The symplectic nature of the phase space is a direct consequence of the existence of the generalized Poisson bracket, which is just a skew-symmetric bilinear form that can be used to define a Hamiltonian vector field in the tangent space of the manifold [22]. The set of trajectories defined by the equations of motion will fill the whole of the dynamically allowed region of the phase space.
The exact classicization mapping opens up the possibility of being able to exploit or to extend all of the techniques developed in classical mechanics for use in the quantum manybody problem. For example, one can make easy contact with conservation laws and the associated sum rules by using the Noether currents. It is also intuitively apparent from the existence of this mapping onto classical mechanics that the c-number NCCM amplitudes {φ I , π I } or {s I ,s I }, which completely characterize and decompose our many-body problem, may be viewed as a set of generalized mean fields which describe each subsystem of particles in the interacting many-body system as labeled by the configuration-space indices {I }.
Let us now return briefly to the non-Hermitian nature of the NCCM approach. The requirement that the physical bra and ket states in the original Hilbert space H (initial) should be Hermitian conjugates of one another, implies from Eqs. (5) and (11) that and hence to constraints between the set of amplitudes {s I ,s I }. Thus, the observable physical space is just a submanifold in the full NCCM symplectic phase space. It is, however, an invariant submanifold in the sense that it comprises entire trajectories, such that either none or all of the points of a trajectory belong to the physical submanifold, under the sole assumption that the original Hamiltonian h in H (initial) is Hermitian. Temporal evolution thus leaves this submanifold invariant. Even within this submanifold the respective pairs of amplitudes s I and s I are not in general complex conjugates of one another since the underlying NCCM similarity transformation is not unitary. These pairs of amplitudes are in no simple relationship to one another. Correspondingly, neither the generalized fields φ I nor their canonically conjugate generalized momentum densities π I , as defined in Eq. (27a) will in general be real. It is only by performing a suitable complex canonical transformation (i.e., a symplectomorphism) that one can regain a description in terms of manifestly real coordinates {φ I , π I } instead of complex ones. There are actually (at least) two different approaches to achieving the above aim. The first, and perhaps most direct, has been performed by Arponen [23], who first introduces the set of complex conjugate NCCM amplitudes {s * I ,s * I }, whereby the NCCM phase space is enlarged into a genuine complex manifold, but self-evidently now of too large a dimensionality. Thus, finally, the extra degrees of freedom are explicitly eliminated by using the Dirac bracket method and deriving the induced symplectic structure, in order to construct the physical manifold (otherwise known as the constraint surface or physical shell), in which the remaining independent amplitudes are a minimal set of complex conjugate pairs or, equivalently, a manifestly real set {φ I , π I }. He shows that the NCCM physical shell is actually a Kähler manifold, and the constraint functions themselves are of second class [24], as is precisely the case in a gauge theory where the constraint corresponds precisely to a complete gauge fixing, even though in the NCCM case the constraints do not seem to correspond to any internal (hidden) gauge symmetries of the problem. The resulting reduced phase space of the physical shell is shown to be itself a complex manifold with a symplectic structure, just like the original extended phase space. The symplectic structure is given by an induced symplectic two-form, whose general form is derived without making any approximations or restricting assumptions. It has also been shown in Ref. [23] that the NCCM star product is well defined in the reduced phase space, or the constraint surface, where the star product is defined to give the expectation of a product of operators in terms of the expectations of individual operators, and where it is also shown that it is implemented by a certain tensor field of type (2, 0).
The second approach to "re-Hermitizing" the NCCM formulation by a canonical transformation (i.e., a symplectomorphism) of the set of amplitudes {s I ,s I } (or, equivalently, {φ I , π I }) has been performed more indirectly by diagonalizing the effective Hamiltonian that governs small oscillations around the equilibrium ground state [4]. It has been explicitly shown in Ref. [4] that this diagonalization procedure is equivalent to performing a canonical coordinate transformation into normal coordinates in the symplectic phase space. The normal frequencies are also then precisely the excitation energies. In terms of the new normal coordinates the Hamiltonian functional H thus recovers a manifestly Hermitian form that was lost in the original NCCM form of the functional H [s I ,s I ] due to the nonunitary nature of the similarity transformation that lies at the heart of the NCCM.

The extended coupled cluster method (ECCM)
Let us now briefly conclude our discussion of the CCM fundamentals by considering the extension of the above ideas to the ECCM. As we have seen, any NCCM expectation value functional Q(t) = Q[s I ,s I ; t] comprises only linked (or connected) terms. The exponen-123 tiated form of the NCCM ket-state correlation operator in Eq. (11) also guarantees that all of the coefficients {s I ; ∀I = 0} are themselves fully linked. However, the corresponding canonically conjugate NCCM bra-state coefficients {s I ; ∀I = 0} do contain unlinked terms. If and when needed, this can be remedied by introducing the following ECCM exponentiated form for the NCCM bra-state correlation operator,S(t), It is clear from Eqs. (6), (12) and (13) that the NCCM cluster coefficientss I may be expressed ass I = C + I , which in general are not linked quantities. The linked pieces of these expectation values that remain after the disconnected terms have been discarded are precisely the new ECCM cluster coefficients,σ I = C + I linked . From Eqs. (7) and (38) we deduce thatΣ(t)|ψ 0 = 0 so that we may write the ECCM parametrizations of the ket and bra many-body wave functions as follows, They are the ECCM analogues of the NCCM forms in Eq. (11). The ECCM form for the expectation value Q(t) of an arbitrary physical observable described by the Hermitian operator q(t) to the NCCM form given by Eq. (13) is then given in terms of a double similarity transformation as follows, It is convenient to move from the (linked-cluster) set of c-number coefficients {s I (t),σ I (t)} to a new (and also linked-cluster) ECCM set {σ I (t),σ I (t)} belonging to the pair of operators {Σ(t),Σ(t)}, where the new linked operator Σ(t) is defined as follows, The inverse pair of transformations between the coefficients σ I (t) and s I (t) are thus given ∀I = 0 as follows The ECCM expectation value functional, , of an observable operator q(t) is now considerably more complicated than its NCCM counterpart Q[s I ,s I ; t], due to the presence of the double similarity transformation. A detailed analysis [2,3] has shown, however, that the ECCM counterparts to the NCCM results of Eqs. (16) and (17) are given as follows, The suffix DL in Eq. (44) now indicates the ECCM double linked structure, which is characterized by the following two constraints: (i) from each group of particles characterized by each of the configuration-space indices K j , the amplitude σ K j associated with the creation operator C + K j , must have at least one particle line connected to the operator q (exactly as in the L-linking of the NCCM); and further (ii) for each group of particles characterized by each of the configuration space indices J k , the amplitudeσ J k associated with the destruction operator C − J k either must have at least one particle line connected to the operator q, or there must be connections, by at least one particle line in each case, to at least two separate amplitudes σ K j associated with two different creation operators C + K j . The basic ECCM amplitudes {σ I ,σ I } may thus be viewed as a set of quasilocal classical fields, due to the maximal connectivity feature built into the ECCM. What is meant by quasilocality here is that each of the amplitudes, which now collectively completely characterize the theory and exactly describe the ground ket and bra states, obeys the cluster property in the usual sense of approaching zero in the limit that any one particle or group of particles comprising the many-body cluster becomes far removed from the remainder. In turn this permits applications to, e.g., topological excitations and cases with spontaneous symmetry breaking, that the NCCM would perhaps have difficulties in describing, due to the non-locality of the {s I } amplitudes describing the bra states within the NCCM, which is overcome within the ECCM by the second similarity transformation.
If we now insert the ECCM wave function parametrizations of Eq. (39) into the definition of Eq. (18) of the action functional, and make use of the relationships in Eqs. (41) and (42), it is relatively simple to show that the action now takes the ECCM form, The independent (bivariational) stationarity conditions, ∂A/∂σ I = 0 = ∂A/∂σ I , now yield the evolution equations, for the ECCM cluster correlation coefficients, which are the the precise counterparts of Eq. (22) in the NCCM case. Equation (46) clearly now shows that the ECCM coefficients σ I andσ I are again a canonically conjugate pair in exactly the same sense as are s I ands I in the NCCM case. As promised earlier, we can now see why the operator S was replaced in the ECCM by the operator Σ via Eq. (41).
It is now evident that the NCCM Eqs. (24) and (25) find their completely analogous ECCM counterparts, The corresponding ECCM expression for the expectation value AB of the product of two ECCM (doubly similarity-transformed) operators is now considerably more complicated than its NCCM counterpart in Eq. (35). Nevertheless, it can be shown [3] that the expectation value of the commutator of operators A and B takes the form within the ECCM, which is the exact analogue of its NCCM counterpart in Eq. (36). It should by now be evident that there are many other analogies between the NCCM and ECCM, which we do not now enumerate further however.
In conclusion, we reiterate that both versions of the CCM can, as we have shown, be formulated in terms of a variational (or, more properly, a bivariational) principle for either the stationary (i.e., time-independent) or time-dependent Schrödinger equations [2,5]. For the stationary (S-CCM) cases of both the NCCM and ECCM the bivariational principle is for the ground-state expectation value functional of the Hamiltonian, while for the time-dependent (TD-CCM) cases it is for the action functional.

Illustration of stationary aspects of the CCM: the Rabi model
We will elaborate further on the CCM in Sect. 5 where we discuss the physics that lies behind the introduction of non-Hermitian operators, after presenting more of the mathematical framework. For now, however, let us elucidate the use of the CCM in practice via a simple, yet illustrative, example. Among a plethora of possible applications, we choose the Rabi model, an important and archetypal model in quantum optics where it describes a two-level atom (with energy spacinghω 0 ) coupled via a dipole interaction (of coupling strength proportional to g) to a single mode of a quantized electromagnetic radiation field (with frequency ω) [25]. In units whereh = 1, the model is defined by the Hamiltonian, in terms of pseudo-spins related to the conventional Pauli matrices, The field mode is described in terms of the annihilation and creation operators, b and b † , respectively, which obey the usual bosonic commutation relation, Within the context of nuclear magnetic resonance, the model also mimics a spin interacting with a field of phonons [26]. A close relationship of the model also exists to the static Lee model of nuclear interaction [27]. From our present viewpoint the Hamiltonian of Eq. (50) is suitable not only as a rather transparent model of physical reality but also as a useful example to illustrate the applicability and efficiency of the standard NCCM approach (and see Refs. [28,29] for more details). A convenient choice of reference state for this model is now |ψ 0 = |0 | ↓ , which is just the ground state of the system for g = 0, assuming ω 0 > 0 (viz. an empty field mode and an unexcited atom), in a notation where the first ket refers to zero bosons in the occupation number representation and the second to the lower state of the two-level atom in an obvious pseudo-spin representation with the spin quantized in the z-direction. For example, the generic NCCM correlation operators of Eq. (12) can now be written in the following specific forms, The only approximation that is now made is to truncate the sums in Eqs. (54a) and (54b) at the term with n = N , giving the so-called SUB-N approximation, in which all the coefficients s In the stationary version of the NCCM, one may now readily evaluate the ground-state expectation value, q = Q = Q[s I ,s I ] of any physical observable q for the Rabi model, as described in Sect. 2.1. For example, the ground-state energy, E = h , has been so calculated at various NCCM SUB-N levels of approximation in Ref. [28], to which the interested reader is referred for further details. We note only that with the above choices of reference state and cluster operators the NCCM was found to yield a (probably spurious) phase transition at a value g → g c ≈ 0.665 of the coupling constant, which signaled a breakdown of the (convergence of the) calculation for values g > g c , which is precisely the region where the ground state becomes nearly degenerate in energy with the first excited state of the system [30]. Still, the method provided excellent results for E in the region g < g c . We defer further discussion of applications of the TD-CCM to this and other systems to Sect. 6.

Excited states and non-Hermitian reformulations of conventional quantum mechanics
We shall consider the question of how to construct excited states and the excitation spectrum of a quantum many-body system within the CCM later in Sect. 7. However, we note now that whenever one tries to move beyond the CCM-based precise evaluations of the groundstate characteristics of a generic (i.e., not even necessarily just many-body) quantum system, one may encounter a number of new, and sometimes unexpected, methodical challenges. At an initial level of approach it appears that one must find a way of making any CCMinspired explicit construction of an optimal Dyson map Ω = exp S of Eq. (4) sufficiently less reference-dependent. For example, in the nuclear-physics IBM-related setting of Ref. [31], this appeared to be a decisive technical obstacle. The essence of the difficulty lies in the fact that the descriptive ambitions of the IBM constructions involved not only the single ground state but also whole multiplets of the lowlying excited states. In this setting, a decisive amendment of the variational results (see the thorough review and mathematical clarification and explanation of Scholtz et al. [32]) only appeared after a return of attention to the older Dyson's papers on ferromagnetism [33,34] 123 in which calculational success was essentially based on the use of a non-unitary operator Ω in Eq. (2).
The latter idea has later been identified as yielding a nontrivial form for the auxiliary operator product Ω † Ω which is henceforth called the metric operator of the physical Hilbert space, where 1 is again the identity operator. During the subsequent developments of the field different authors succeeded in finding techniques circumventing the obstacles emerging due to the nontrivial nature of the metric (i.e., that it is not simply the identity operator), which is not encountered in conventional Hermitian quantum mechanics. The computational economy of the Hermitian Schrödinger-picture (SP) formulation of quantum mechanics [35] is very persuasive. It is just this aspect of the formulation that explains its popularity and success in many branches of quantum physics and quantum chemistry. Nevertheless, several existing technical limitations of the conventional SP description of quantum mechanics have almost always forced its users to search for amendments. In this context a short letter by Bender and Boettcher [36] proved to be very influential in returning attention to the less usual but fully admissible possibility of non-Hermitian representations of observables as promoted, a few years earlier, in Ref. [32]. The latter paper itself actually recalled the older idea of mapping certain Hermitian fermion operators onto simpler non-Hermitian boson operators by means of the generalized Dyson mapping [37], with all of these developments being precisely the motivation behind the IBM technique itself [7][8][9].
Such a "crypto-Hermitian" amendment (or perhaps, better, extension) of our understanding of quantum theory has, at present, multiple parallels and continuations [14]. The resulting, truly deep, theoretical reconsideration of the first principles of quantum theory has now become widely accepted by a broad community of physicists. The impact of the idea may nowadays be detected even beyond the domain of quantum theory [38]. Still, its origins may be dated back to the studies of certain truly complicated many-fermion quantum systems. It was, in fact, Dyson [33,34] (and, independently, Maleev [39]) who proposed, more than 60 years ago, and mainly for the purely practical purposes of variational calculations, a key idea lying in the (formally reversible) non-unitary mapping of wave functions as in Eq. (2).
The reader should be aware at this point that up till now we have not precisely defined the term "Hermitian" on the mathematical level. In general, one would need to distinguish with care such alternative concepts related to the term "Hermitian" as "symmetric", "self-adjoint" (on a specified domain), and "essentially self-adjoint" (with a well-defined core). However, in most of the prospective applications that we have in mind we work with a finite-dimensional matrix form of the so-called "Hermitian" operator, in which case there is no longer any need for such mathematical precision. Marginally, let us note here that, even in the infinitedimensional and non-matrix models, the essence of many related apparent paradoxes can be identified as lying in the fact that the physical Hilbert space is exclusively presented via its representation in H (user−friendly) [unphysical] . Thus, the clarification of the Hermiticity/non-Hermiticity misunderstandings becomes simple when one consequently stays in the latter space and characterizes the switch to H (final) [physical] by the mere change of the inner product (and the interested reader can find further useful comments on this trick in Ref. [14]).
Due to the latter trick, phenomenolgical Hamiltonians may be, admittedly, manifestly non-Hermitian in a preselected and, presumably, just unphysical, auxiliary Hilbert space H (user−friendly) [unphysical] . As a concrete example let us recall the most popular, ordinary differential, upper-case Hamiltonian H (imaginary cubic) = − d 2 dx 2 +ix 3 , which is PT -symmetric [13,36] but manifestly non-Hermitian in H (user−friendly) [unphysical] = L 2 (R). After an appropriate amendment of the inner product, according to Bender [13], the underlying quantum system S may be assigned its conventional quantum probabilistic unitary-evolution interpretation via the reconstruction of the physical Hermitian conjugation [c.f., the "physical" conjugation H → H ‡ ≡ Θ −1 H † Θ in place of the conventional, friendlier but unphysical, conjugation H → H † .
There exist many interesting aspects of the model described by H (imaginary cubic) and by its various alternatives, which are explained and discussed in, e.g., the reviews presented in Refs. [12][13][14]. For all of them the underlying probabilistic interpretation of the quantum systems is based on the innovative use of non-Hermitian Hamiltonian avatars H with real energy eigenvalue spectra. Almost without exception, all of these examples are presented, in the current literature, in the framework of what we will call here the generalized (i.e., "non-Hermitian") Schrödinger-picture formalism. In its stationary case we will call it, for the sake of definiteness, the Dyson-Maleev formalism (DMF), for a particularly extensive review of the generalized forms of which we may particularly recommend Ref. [37].
Once one manages to achieve at least a reasonable degree of reference-independence of Ω, a suitable combination of the CCM mathematics with IBM physics could open new construction horizons. Indeed, in the strictly stationary case, i.e., in the case with property Ω(t) = Ω(0) , ∀t, there exists an intimate connection and correspondence between the generality (i.e., non-unitarity) of Ω, the non-Hermiticity of the avatars H of the Hamiltonians given by Eq. (3) and the bivariational nature of the CCM recipes, as we discussed in our previous paper [1].
The DMF approach may be perceived as one of the standard numerical algorithms, which transfers the description of the states from the traditional (i.e., often, fermionic, Fock) Hilbert space H  [31,33,37,[39][40][41][42][43]. For this reason the predictions of the DMF approach are essentially identical to those of the conventional Schrödinger picture, provided only that the third, final space is, by construction, assumed unitarily equivalent to the initial one, H (final) ∼ H (initial) . For this reason the key features of the stationary DMF scheme (typically, in its nuclear-physics IBM implementations) remain transparent. Its structure may be summarized by the following compact illustrative flowchart diagram, One of the most remarkable consequences of the non-unitarity of the mapping Ω, seen in such a diagram lies in the coexistence of the Hermiticity and non-Hermiticity properties of the same upper-case operator H , depending on which Hilbert space one is considering. In the stationary setting the above-mentioned operator H represents an observable quantity, despite being non-Hermitian (in the sense that

Transition to non-stationary quantum systems
Our present study was inspired by the recent progress in the development of what could be called reference-independent IBM-like theories [11,[44][45][46][47]. The message of our present paper should be seen in the transfer of our current understanding of the merits of the mappings of Eq. (2) from the IBM-related context into the bivariational constructive CCM strategies. The NHIP formulation of quantum mechanics will be used, in such an application, in a slightly narrower sense, being reserved to cover only the picture of reality with the non-stationarity property, According to the first consistent introduction of the full-fledged non-stationary NHIP formalism in Ref. [44], its basic idea may be perceived as a time-dependent extension of the old IBM-like variational recipe. While its innovated form is definitely more flexible, it is undoubtedly also much more complicated technically [45]. For this reason it is, therefore, perhaps not too surprising that the TD-CCM/NHIP relationship has not yet been studied. It is precisely this omission that we aim to remedy here. This technical complexity also probably explains why only a relatively few realistic applications of several alternative implementations of the NHIP ideas themselves may yet be found in the current literature [45,[47][48][49][50][51][52][53][54][55][56].
Since the birth of quantum theory in its Heisenberg-picture (HP) formulation [57] and, less than a year later, in its Schrödinger-picture (SP) formulation [35], those attempting to apply the theory have always needed to resolve the conflict between the more intuitive nature of the HP quantization of observables and the maximal economy of the transfer of attention to the wave functions in the SP approach. A partial relief of this tension came, later, with the invention of the more universal interaction-picture (IP) (alias the Dirac-picture) formalism from which the SP and HP descriptions of quantum dynamics could have been deduced as two separate special limiting cases. An enhanced flexibility of the ensuing language then also gave rise to the well-known successes of the manifold IP applications (and predictions) in practice, especially in the context of quantum field theory and quantum many-body theory [58,59].
In Table 1 we present a comparison of the HP, SP and IP "strictly Hermitian" descriptions of unitary (i.e., stable) quantum evolution. As a comment on the table we might emphasize that the same physics is described by the single operator evolution equation in the HP formulation, and by the single ket-vector evolution equation in the SP formulation, as well as by a pair of evolution equations in the IP formulation. For a compensation of the seeming disadvantage of the doubling of the number of evolution equations in the IP formalism, we note again that the latter IP picture contains both the former (HP and SP) ones as special limiting cases.
The ultimate choice between the HP, SP and IP (or, indeed, many other [60]) modelbuilding strategies depends, first of all, on the actual form of our specification of the quantum system in question. For this reason, most standard textbooks usually prefer the SP language, only adding the HP and IP analyses of quantum dynamics at the later stages of explanation. This makes the SP-based specification of quantum dynamics less intuitive but shorter, based on the rather formal introduction of a "physical" Hilbert space H (initial) and of a suitable self-adjoint "Hamiltonian" h = h † defined within that space. Many researchers prefer the use of the SP language in practice, since it requires, in addition, a maximally realistic origin for, and "derivation" of, the latter Hamiltonian operator, which is typically found, e.g., via a "quantization" of its suitable classical-physics counterpart.

Non-Hermitian versions of the Heisenberg and Dirac pictures
The key to the extension of the validity of the IBM-type pattern beyond its stationary DMF version of Eq. (56) has been found in Ref. [44]. First of all, the removal of the existing theoretical obstacles and objections (and see, e.g., Ref. [61]) required a refinement of the terminology. In Ref. [44] it has been emphasized that the initial and observable physical Hamiltonian h(t) (defined in H (initial) [textbook] ; in general it may be time-dependent) merely becomes replaced, as in Eq. (58), by its upper-case isospectral non-stationary avatar, which is non-Hermitian in H (user−friendly) [unphysical] but Hermitian in H  (57) has turned out to be much less straightforward. In the literature, the process of this clarification was both lengthy and tedious [44,[61][62][63][64][65][66]. Fortunately, at the end of this process in 2009 (c.f., Ref. [11] and some later addenda in Refs. [52][53][54]) the eventual outcome has transpired to be both relatively elementary and transparent, as we now briefly explain.
Its brief summary may start from a return to Table 1, which reminds us that even in the conventional Hermitian SP setting the SP → HP transition may be perceived as mediated by a mapping of the form of Eq. (2), but in which the operator Ω would be unitary but manifestly time-dependent. In 2007, the feasibility of the extension of the SP ↔ HP correspondence to non-Hermitian cases was opposed by Mostafazadeh [14,61,63,66]. Fortunately, the apparently insurmountable obstacles and obstructions that were initially raised against the free applicability of the non-stationarity postulate, thereafter appeared in essence to be of a purely terminological nature [44,62,64,65,67,68]. As a consequence, the transition to the generalized "non-Hermitian" HP (viz. the NHHP) was eventually formulated definitively in Ref. [46]. Soon thereafter it was also found to be both feasible and useful in practice in some specific applications [69].
Once we are given a lower-case SP operator q representing an arbitrary observable quantity, the only relevant task for theorists is a prediction of the results of experiments based on the evaluation of its expectation value. This means that we have to evaluate the quantity

ψ(t)|Q(t)|ψ(t) .
The precise meaning of the ket state |ψ(t) and of the bra (or, rather, brabra) state ψ(t)|, are to be specified of course (see Sect. 4.1). The upper-case symbol Q(t) just represents here the operator (isospectral to q), which is defined in the space H (user−friendly) [unphysical] via a similarity transformation in complete analogy to its NCCM counterpart in Eq. (14a). The overall theoretical flowchart then has the following compact form, .

Ξ(t) → equation for any Q(t) G(t) → equation for kets |ψ(t) G † (t) → equation for ketkets |ψ(t)
The precise meaning and definitions of the operators and state vectors in Eq. (59) is given in Sect. 4. On this basis we may expect that the use of the alternative formulations of the dynamical equations will not hide their one-to-one correspondences with those in the Hilbert space H (initial) [textbook] of the topmost-box in the diagram comprising Eq. (59).

Non-Hermitian versions of dynamical evolution equations
Curiously enough, the above-outlined NHHP formalism is "quasi-stationary" because the HP metric itself remains time-independent [45,46], As a consequence, it was necessary to move beyond the constraint (60) (c.f., Ref. [47]). The current state of the art is summarized in Table 2 in which the third "non-Hermitian", fully general formalism of Refs. [44,47] is now assigned the abbreviation NHIP. Clearly, Table 2 is just the general counterpart in our non-Hermitian context of the earlier Table 1 that pertains to the conventional Hermitian-generator description of unitary evolution.

Recommended notation conventions
In the most general NHIP context the use of the non-stationary, Dyson-motivated Ansatz of Eq. (2) seems to open a Pandora's box of thorny problems. First of all, once we accept the fact that Ω = Ω(t) is a map which inter-relates a triplet of Hilbert spaces, we find that the notation is insufficient and/or incomplete. Thus, first of all, it does not inform us that such a map connects the initial-reference ket ψ 0 ∈ H (initial) [textbook] and the final-reference ket ψ ∈ H (user−friendly) [unphysical] , which may then itself be re-read, alternatively, as the correct physical ket ψ ∈ H (final) [physical] . Secondly, one would also like to avoid using the subscript 0 in ψ 0 because the symbols ψ and ψ 0 refer, in fact, to the same quantum state in the NHIP representation.
Both of these inconsistencies of notation were successfully removed in Ref. [11], in which it was shown how all of the unnecessary repetitions of the, otherwise necessary, explanatory comments may be circumvented by using the following three simple amendments of the standard Dirac notation for state vectors, viz. by using the triplet of replacements The defining relation between the ket state |ψ(t) and the ketket state |ψ(t) in Eq. (61) validates our earlier assertion in Eq. (60) that the HP metric operator Θ (HP) is stationary since, by definition (and see Table 2), neither of these states evolves in time in an HP formalism. Continuing with the exposition of our recommended notation convention in the three-Hilbert-space approach, we note that the NHIP version of our fundamental Ansatz of Eq. (2), viz.
in terms of the ket state |ψ(t) , may equivalently be re-written as In the derivation we have employed the definition of Q(t) given in Eq. (58). Using our notations we thus also reveal the correct physical status of the final Hilbert space H (final) [physical] , since Eqs. (64) and (65) clearly display the equality of the measurable quantities in the two respective Hilbert spaces H (initial) The required matrix element in the latter space [i.e., that on the right-hand side of Eq. (66)] is precisely the one that we introduced earlier in the discussion in Sect. 3.2.
We never need to leave the auxiliary Hilbert space. One can say even more: thus, any physical quantum state can be characterized by ket |ψ(t) and metric Θ(t) or, much more economically, by the pair comprising ket |ψ(t) and ketket |ψ(t) . Similarly, one can say that any physical observable can be characterized by its "hiddenly-Hermitian" operator Q(t) and metric Θ(t) or, without an explicit use of the metric, by the pair of operators Q(t) and Q † (t), mutually connected by the metric, Equation (67), which is easily derived from the definitions of Θ(t) and Q(t) in Eqs. (55) and (58), respectively, is just the so-called hidden-Hermiticity (alias quasi-Hermiticity) relation [32,70] for an operator Q(t) in the NHIP formalism that belongs to a physical observable. It is, of course, completely equivalent to the relation which provides simply the obvious generalization of the definition of H ‡ in the flowchart of Eq. (56). It is worth noting that the operator Θ itself is both Hermitian in the ordinary sense and quasi-Hermitian in the sense of Eq. (68). More generally, so is any operator that may be expressed as an arbitrary power series in Θ with real-valued coefficients.
As an aside here, we can now directly compare the present NHIP formalism with the earlier NCCM description. Thus, in the NCCM the reference (or model) state ψ 0 is assumed to be stationary (i.e., time-independent). By comparison of Eqs. (13) and (14a) with Eqs. (64) and (65,), we also see that the NCCM operatorS(t) is (up to a multiplicative normalization constant) precisely equal to the NHIP metric operator, viz.S(t) = Θ(t)/ ψ(t)|Θ(t)|ψ(t) , in our newly recommended NHIP notation.

Evolution equations for states
In the light of the definition, given by Eq. (57), of the isospectral non-Hermitian partner H (t) of the textbook Hamiltonian h, one can now claim that the upper-case operator H (t) is, indeed, quasi-Hermitian, in the sense of Refs. [32,70,71], and in accord with Eq. (67). Such an observable for the Hamiltonian is now to be interpreted as an instantaneous total energy. Its observability status reflects the observability of its (manifestly Hermitian in the usual sense) isospectral partner h = h † . Several authors (ranging from Bíla [50,51] up to his most recent followers [67]) decided to re-assign the status of Hamiltonian to another operator. The explanation of this slightly surprising decision is purely terminological. A detailed disentanglement of the puzzle has been given recently [47], which we now briefly summarize since it is very pertinent to our own further developments here. It starts from the conventional initial-space SP evolution equation given by Eq. (1), now rewritten in our new notation of Eq. (61) as and thence from its replacement by the preconditioned NHIP alternative defined in H (user−friendly) [unphysical] . The latter, equivalent Schrödinger-like evolution equation for ket states |ψ(t) , is easily derived from Eq. (62). It contains, naturally, the generator of evolution which is composed of the energy operator H(t) in combination with the so-called Coriolis operator Ξ(t) [45], as we observe explicitly in Eq. (72). For this reason, the re-assignment of the name of the Hamiltonian from H(t) to the newly constructed operator G(t), as has been suggested by several authors (see, e.g., Refs. [50,51,67]) as mentioned above, is, in our view, both misleading and unfortunate. In our notation any given quantum state is labeled by the same Greek letter (say, ψ). Hence, we conclude that once we are describing the state in terms of a pair of independent vectors |ψ(t) and |ψ(t) , we may utilize Eq. (61), which gives the relationship between them, together with Eqs. (69), (71), and (72), to derive the complementary, second, independent evolution law for ketket states |ψ(t) , In this equation, nobody has, as yet, further re-assigned the traditional name of the Hamiltonian to G † (t). Thus, for our current purposes we will reserve the name of an "evolution generator" for both of the operators G(t) and G † (t).

Evolution equations for observables
Given the NHIP Coriolis operator and the initial (i.e., at time t = 0) values of the metric Θ(0) and of an arbitrary observable operator Q (0), such that we may reconstruct the full time-dependence of Q (t) via the following NHIP version of the Heisenberg equation, We may also work with its conjugate version for complement Q † (t),

123
In the light of Eqs. (67), (75) and (76), the initial condition of Eq. (74) at time t = 0 will guarantee the observability status of Q (t) at all times, as required of course, as we now show explicitly. In particular, if we define the operator, which is not an observable, and hence does not evolve as given by Eq. (75), we may use Eqs. (75) and (76), together with the definition of Θ(t) in Eq. (55), to show directly that its actual evolution equation is instead given by which clearly has the solution Z(t) = 0 for all times t under the initial condition Z(0) = 0, as given by Eq. (74). Indeed, straightforward differentiation of the definition of Q(t) in Eq. (58) (and of its conjugate) with respect to time yields immediately (the two commutatorcontaining) Eqs. (75) and (76). Naturally, we assume in so doing that the inaccessible SP operator q does not vary with time. We do not consider further here the strongly anomalous case of a manifestly time-dependent operator, q = q(t), except to note that it has also been discussed in Ref. [46]. Let us now return to Eq. (71) [skipping an analogous analysis of its adjoint partner in Eq. (73)] and let us notice that it contains the operator G(t) of Eq. (72) which is not only non-Hermitian but even non-quasi-Hermitian in our working Hilbert space H (user−friendly) [unphysical] . It is worth adding that it plays not only the role of the generator of the evolution of the statevector kets, as exhibited in Eq. (71), but also the role of the generator of the evolution of the observable instantaneous energy. It is easy to verify the validity of the evolution equation where the last equality, which follows simply from the definition of G(t) in Eq. (72), is thus completely compatible with Eq. (75). Equation (77)  Even in the case of strictly unitary evolution and even in the HP subcase of the general (and entirely methodical) NHIP scheme, as outlined above, both of the operators Ξ(t) and Ξ † (t) should properly be called Coriolis-force potentials (or, more simply, Coriolis terms) rather than Hamiltonians. The name Hamiltonian should, for reasons thus outlined, remains consistently restricted only to the quasi-Hermitian observable H (t) and/or, whenever needed, to its conjugate operator H † (t). This is only reinforced by the fact that none of the other, rejected candidates (from those introduced above) for possible promotion to Hamiltonian status even need to have a real eigenvalue spectrum in general, as may rather easily be illustrated.
Thus, for example, we may readily prove the relation, remains time-independent.

Physics behind non-stationary non-Hermitian operators
The initial motivation of our interest in the mutual relationship between the CCM approach to quantum many-body theory and the so-called three-Hilbert-space formulation of quantum theory was both formal and physical. We noticed that in both of these approaches one tries to combine the idea of the availability of a straightforward, formally friendly, approximation (or of a sequence of approximations) of ψ by a simpler ψ 0 , with the practical awareness of the weak points (i.e., typically, of a slowness of convergence) of the respective approach. For a formal remedy to the latter weakness one may then immediately turn to the, by now, rather standard mathematical technique of the so-called preconditioning [72][73][74][75], the essence of which lies in the use of the factorized Ansatz of Eq. (2).

Physics of time-dependent correlations
For both the CCM and three-Hilbert-space formalisms, successful examples of the practical implementation of the Hilbert-space-mapping approach abound. In just the specific context of many-fermion quantum physics, for example, we may recall, e.g., the pioneering works of Dyson [33,34], Maleev [39], Coester et al. [16,17],Čížek et al. [18,19,76], and Janssen et al. [37], in all of which the above broad approach was adopted. The subjects and applications of these papers covered a broad spectrum of applied quantum theory, including condensedmatter physics [34,77], nuclear physics [7][8][9]78], and the descriptions of large atoms and molecules in quantum chemistry [79,80].
The essence of all of these techniques of making predictions concerning complicated multiparticle structures may be said to lie in a combination of a systematic mathematics behind approximations (reflecting, typically, the calculation feasibility aspects) with an intuitive insight into the relevance of various competing phenomena (i.e., typically, of the various types of the correlations). Thus, in a way which motivated our present paper, a transfer of the years-long experience from the static to non-stationary systems and calculations may be perceived as a true challenge.
We came to the conclusion that a key component of such a transfer should be seen in a clear abstract formulation of the formalism itself. In this sense, we tried to propose and promote a combination of the NHIP language (emphasizing, in essence, the clear separation of the state-and operator-evolution equations) with the CCM-related constructive efficiency. We believe that such a combination will pave the way towards an increase of efficiency of the calculations, based on the use of creation and annihilation operators which is, simultaneously, linear (i.e., related to the construction of the bases) and non-linear (entering the mimicking of correlations via the judicious forms of the exponential mappings Ω).
As we have already noted, at various levels of approximate implementation of any version of the CCM a loss of manifest Hermiticity between corresponding ket and bra states can arise from one of its basic tenets, viz. that formally the ground ket and bra states are parametrized independently [2,5]. In practice any potential shortcomings that this may give rise to are far outweighed by the fact that in so doing one ensures the exact maintenance of the important Hellmann-Feynman theorem [20,21] at all such levels of approximation (and see Ref. [10] for details). Thus, calculations of the ground-state expectation values of a Hermitian operator q that represents an arbitrary observable quantity are completely compatible with those for the energy expectation value, in the sense that the former can be obtained from the usual (perturbation-theoretical) Goldstone diagrams for the energy by replacing in turn each interaction term in h by the operator q.

The role and interpretation of the creation and annihilation operators
One of the strongest motivations of our present goals of uncovering a reinterpretation and translation of the TD-CCM techniques into the NHIP language (and vice versa) may be sought in the key merits of the CCM recipe in which one relies, heavily, upon the traditional explicit use of creation and annihilation operators, as discussed in detail in Sect. 2. This is a decisive merit which reflects our experience, intuition and pragmatic perception of the underlying phenomenology, especially when one deals with correlations carrying certain characteristic features of the quantum particle clustering. Moreover, the explicit use of the concept of clusters also leads to multiple vital simplifications of the explicit constructive calculations.
In the opposite direction, in the NHIP language one shifts attention to several related abstract concepts, such as the correspondence principle and some of its less obvious consequences [14,81]. The conservation of probabilities, which is guaranteed by the underlying unitarity of the temporal evolution in our non-Hermitian pictures is thence put in a new perspective, especially because in the non-stationary dynamical regime the energy clearly then ceases to be conserved. The related operator will represent just the instantaneous total energy.
The NHIP framework seems more (perhaps, even too) general and abstract in most application-oriented contexts. At the same time we saw previously [1], that the very essential advantage gained by building bridges between the NHIP-related and CCM-related constructions could be viewed as a certain optimal balance between the abstract concepts making use of an explicit reference to the specific features of the systems and operators in question.
We discussed in Sect. 2 how the NCCM and ECCM parametrizations of an arbitrary quantum many-body theory enable it to be mapped exactly onto a classical Hamiltonian mechanics for the many-body, classical (c-number) configuration-space amplitudes that completely and exactly describe the ket and bra ground states independently. As we saw explicitly, these mappings arise fundamentally from an underlying one-to-one correspondence that can be proven to exist between the set of commutators in the original quantum many-body Hilbert space and a set of suitably defined generalized classical Poisson brackets.
It has further been shown [3] how the CCM (particularly in its ECCM form) can be interpreted as an exact generalized mean-field theory (i.e. beyond the Hartree-Fock level) formulation of the given quantum many-body problem. This interpretation is itself closely linked with the additional realization that the ECCM can also be construed as an exact bosonization procedure in which the ECCM states are associated in a one-to-one fashion with a set of generalized coherent states in some suitably defined boson space. This ECCM bosonization procedure differs from other such procedures in the sense that the usual motivation for any bosonization scheme is taken to its logical conclusion, viz. that the resultant generalized coherent boson fields are classical c-number fields with only classical (but highly nonlinear) interactions between them. Being able to reinterpret the ECCM as an exact generalized mean-field theory is then reinforced by being able to show that, within the ECCM bosonization scheme, commutators of operators in the original Hilbert space are mapped only onto the tree-level pieces of the corresponding commutators of the respective mapped operators in their boson image space. The tree level of a commutator is here defined to be a restriction only to such contractions that do not result in closed loops. The subsequent manifest exclusion of all closed-loop diagrams thus acts to further reinforce the fact that the ECCM exactly reformulates the quantum-mechanical many-body system that we start with as a classical generalized mean-field theory.
6 Illustrative examples of the TD-CCM approaches

The Rabi model revisited
In view of our general discussions in Sects. 2 and 4 of the fundamentals of the TD-CCM and NHIP formalisms, respectively, and the overlaps that we have seen exist between them, we can now appreciate how the use of the multi-configurational creation operators {C + I } and their destruction counterparts {C − I } in the NCCM parametrizations of Eqs. (11) and (12) leads to a technically feasible implementation of the NHIP concepts, via the associated use of the c-number cluster correlation coefficients defined in Eq. (12), as we now again illustrate.
Thus, firstly, in Sect. 2.3 we outlined the basic ideas behind a prototypical application of the stationary version of the NCCM to the specific case of the Rabi model, whose Hamiltonian is as specified in Eq. (50). The model has also been studied in the non-stationary regime using the time-dependent version of the NCCM [29], in which the basic coefficients {s Perhaps the most important atomic quantity of interest associated with the Rabi model, particularly since it is amenable to experimental observation [109], is the so-called atomic inversion , viz. σ z (t). It is easy to show that in the NCCM parametrization of Sect. 2.3 it may be calculated as follows, Similarly, the most important field observable is the photon number, whose associated operator is n = b † b. Once again, it is easily shown that its time-dependent expectation value has the simple NCCM parametrization, As was done in Ref. [29], an obvious choice of initial condition at time t = t 0 to illustrate the efficacy of the method is to start the system in the reference state itself, i.e., the state |0 | ↓ , which comprises an empty field mode and an unexcited atom, for which s  As observables, both σ z (t) and n (t) should always be real. However, the SUB-N truncations of the NCCM cluster operators inevitably leads to the exact Hermiticity of the corresponding bra and ket states being broken, which means that observables are not constrained to be manifestly real at any such level of approximation. Calculations show (and 123 see Ref. [29]) that for small values of the coupling g the imaginary parts of the observables are extremely small, and decrease as the truncation index N is increased. Thus, for relatively small couplings the evolutions of the calculated values of σ z (t) and n (t) clearly show [29] a (quasi-)periodic exchange of energy between the atom and the field. Nevertheless, as g is increased, the restriction on the maximum SUB-N level that can be attained in practice for higher couplings, ensures that any spurious complex parts cannot be entirely eliminated. Thus, although the source of the errors is well known, a mathematically robust and wellfounded means to suppress, and eventually eliminate, them remains unknown. One of our hopes now is that a resolution of the error-control puzzle might be provided by establishing a deeper correspondence between the NHIP-related abstract "bookkeeping" of the relevant operator evolution equations, and the various pragmatic CCM solution techniques. A best possible outcome would thereby be to improve upon the standard CCM SUB-N truncation hierarchy. We return to this point below after discussing our second illustrative example of an application of the TD-CCM formalism.
Before doing so, however, we remind the reader that the above simple choice of initial conditions was made purely for ease. Other choices are certainly possible, and for an arbitrary physical state as the starting state, we can always first calculate its NCCM representation to yield its corresponding coefficients {s (i) n ,s (i) n ; i = 1, 2}, which may then be taken as the initial conditions at t = t 0 .

Condensed Bose fluid
We now turn to our second, much more ambitious, example of an application of the TD-CCM, viz. to describe the zero-temperature hydrodynamics of a macroscopic condensed Bose fluid, this time via the ECCM [110]. A convenient choice of reference state for this model is now |ψ 0 = |0 , the bare vacuum. Accordingly, we are thus working in a numbernonconserving formulation, as introduced by Bogolubov [111]. Particle-number conservation is then imposed by working with the grand canonical Hamiltonian k (0) ≡ h (0) − μn rather than with the Hamiltonian, h (0) , where μ is the chemical potential and n is the particle number operator. We consider here a system of N identical bosons, each of mass m and interacting via pairwise potentials, such that (in first quantization), We parametrize the many-body configuration space in terms of real-space coordinates, such that the single-boson creation and destruction operators, b † (x) and b(x), respectively, which act to create or annihilate a particle at the three-space point In the original Hilbert space with fixed particle number N , Eq. (10) now takes the specific form, We can now see immediately that the ECCM formalism of Sect. 2.2 provides a framework in which to consider the general case of a non-stationary and spatially non-uniform condensate, b(x) (t) = ψ 0 |eΣ (t) e −S(t) b(x)e S(t) |ψ 0 . Indeed, by making use of Eqs. (14b) and (15) and the fact that S(t) has a corresponding expansion to that for where the last equality follows from the generic definition of Eq. (42). Thus the condensate wave function b(x) (t) is precisely given by the one-body ECCM amplitude σ 1 (x; t). Thus, it is no surprise to learn that at the lowest SUB-1 (i.e., one-body mean-field) level of truncation of the general SUB-N hierarchy (in which all ECCM amplitudes σ n andσ n with n > N are set to zero) our formalism simply reduces to the well-known Gross-Pitaevskii [112,113] description of the condensate wave function (or one-body order parameter). Hence, our usual SUB-N approximation hierarchy can now be used to improve systematically upon this lowest-order one-body mean-field description.
In order to take full advantage of the local U (1) gauge symmetry that follows from particlenumber conservation it is now very convenient to couple the system to external U (1) gauge fields φ(x, t) and A(x, t) (viz. the scalar and vector potentials, respectively), and hence we now study the new grand canonical Hamiltonian k, where (in first quantization), Equation (84a) may be equivalently rewritten in the second-quantized form, Using this grand canonical Hamiltonian k as the time-evolution operator, one may now explicitly evaluate [110] the equations of motion (46) (with H → K) for the ECCM amplitudes {σ n ,σ n }. The whole treatment may now be developed in a completely gauge-invariant fashion so as to provide a complete (zero-temperature) hydrodynmical description. In particular, one may study the evolution of the expectation values of the off-diagonal one-and two-body density operators. One can thereby derive in a wholly gauge-invariant form the exact balance equations (i.e., the local conservation laws) for the hydrodynamically relevant 123 variables (viz. the particle-number density, current density, and energy density, all considered as functions of spatial position and time), and show how these are in turn related to (and derivable from) the cumulant expansion of the one-body density matrix. It has also been shown [110] that each of the balance equations is exactly obeyed at every SUB-N level of truncation.
We note that the U (1) gauge invariance that we have been at pains to incorporate exactly, can, of course, itself be interpreted as the (differential or) local Galilean invariance. An immediate consequence is that the description correctly separates the average translational (hydrodynamical) kinetic energy from the average kinetic energy in the local rest frame, which has otherwise proven to be very difficult to achieve in a fully microscopic treatment. One can clearly reinterpret this exact feature of the time-dependent ECCM treatment of a condensed Bose fluid with the generic proper treatment within the NHIP of the analogous Coriolis terms Ξ(t) of Eq. (72) in the evolution generator operator G(t). Despite being beyond the scope of the present paper, it would be interesting in this light to pursue this relationship further, since one might expect it to extend our intuitive insight into both methods, as well as to lead possibly to improved computational (i.e., approximation) schemes.
Before closing this discussion, however, we note that the temporal evolution of the system has been described by a trajectory in the ECCM symplectic phase space spanned by the cluster coefficients {σ n ,σ n }. The modern viewpoint of statistical mechanics [22], however, is that a proper qualitative description of the system is obtained, rather, from the entire phase portrait , which is just the totality of all allowed trajectories. From this stance, it is clear that the above single-trajectory ECCM approach might now usefully be extended by focusing more on the geometric properties of the ECCM phase space, as have been discussed in some detail in Ref. [96]. By contrast, we have so far focused more on the algebraic structure of the ECCM phase space (and see Ref. [95] for further details), and in so doing we have been led inexorably to the SUB-N hierarchy of approximations. As we have seen in Sect. 2.3, for example, this scheme sometimes has computational limitations, which stem essentially from the fact that, although the ECCM has left its ultimate origins in perturbation theory far behind, the SUB-N scheme itself does still seem to have some links with perturbation theory. The hope has been expressed [110] that one might be able to exploit the geometric structure of the ECCM phase space to derive wholly new (and potentially more robust and more powerful) approximation schemes based wholly on its geometric properties. Again, we note that, although such an investigation lies beyond the scope of the present paper, any progress in this area will certainly find immediate impact in the NHIP formalism too.

Outlook
In the arena of realistic and predictive CCM calculations, most attention is usually paid to the study of molecules and/or other quantum many-body systems in their ground state. Once our interest shifts to the excited states the strategy has to be modified, as we alluded to in Sect. 3. Within the S-CCM such a modification was first proposed by Emrich [83,84] within the context of the time-independent (i.e., stationary) Schrödinger equation. This was done by a suitable modification of the ground-state parametrization, which did not involve a change of reference state. A later alternative TD-CCM approach to excited states was advocated by Arponen and his co-workers [2,4] within the context of the time-dependent Schrödinger equation, wherein it was shown how excited states could be obtained as the normal modes of a suitably defined dynamical Hamiltonian matrix obtained within linear response theory, in direct analogy to the usual procedure in classical mechanics. The complete equivalence of the two approaches has been formally demonstrated [4,5,10]. In effect, what we have done in the present paper is to build upon, and extend, these results.

Reference-dependence versus reference-independence
In the broader context of constructive quantum many-body theory the key purpose of our present paper may be seen in the description of the close relationship between the referencedependent and reference-independent forms of the Ansatz of Eq. (2). In this sense our present study of parallels between the alternative NHIP and TD-CCM approaches to non-stationary dynamics was preceded by our earlier paper [1] in which the methodical framework was perceivably simplified and reduced to the less sophisticated search for parallels between the stationary DMF and S-CCM methods.
As a concise summary of the message provided by Ref. [1] it can be said that the reconsideration of the S-CCM method from the reference-independent DMF perspective may be well illustrated in diagrammatic form by the flowchart in our present Eq. (56). Therein one assumes an input knowledge of the Hamiltonian operator h and of the Dyson mapping Ω. In this setting the main weakness of the reference-independent DMF recipe lies in the absence of any hint of how we should choose the operator map Ω. Hence, the stationary DMF approach may, somewhat crudely, but simply, be characterized as guesswork, or as an intuition-based implementation of a purely trial-and-error-type strategy.
This does not mean that the DMF constructions were not successful in practice. In fact, the opposite is true [32]. Still, the S-CCM techniques are more systematic because, given a Hamiltonian operator h, the S-CCM specification of the necessary Dyson operator Ω = exp S is not at all arbitrary, with the rationale for the exponentiated form as we have outlined above. Hence, the S-CCM techniques are both "algorithmic" and extremely intuitive. A related major advantage of the S-CCM lies in the use of the very concrete truncated-series expansions of the Dysonian exponent operator S in an operator basis formed by certain multi-configurational creation operators, as alluded to above. This opens the dual possibilities of multiple technical simplifications and extremely powerful approximation hierarchies that are well-defined and physically motivated, both of which are precisely due to the intrinsic reference-dependence of the S-CCM construction.
In the present paper we have pointed out that a similar balance between the merits and shortcomings must also necessarily exist (and should prove useful) when one replaces the S-CCM ↔ DMF stationary correspondence by its TD-CCM ↔ NHIP time-dependent generalization. A priori one may expect that in the generalized setting, any technical simplification (obtained from the mutual enrichment) could play a more decisive role, first of all, because of the enormous overall increase in mathematical complexity that emerges after the respective replacements S-CCM → TD-CCM and/or DMF → NHIP.
Particular attention has to be paid to a deeper theoretical role played by the transfer of mathematical know-how from the DMF and NHIP formalisms to their respective S-CCM and TD-CCM counterparts. In this manner the brute-force S-CCM and TD-CCM constructions appear open to further amendments. A typical example is the possible replacement of a certain ad hoc, auxiliary S-CCM or TD-CCM operator admitted by the principles of bivariationality, and hence needed to characterize the bra states independently from the ket states (viz. the operatorS in the NCCM orΣ in the ECCM) that is thus required for the evaluation of the expectation values of an arbitrary physically observable operator, by alternative forms that maintain the canonical CCM symplectic structure discussed in Sect. 2, which are both inspired and restricted by the implicit mathematical merits (e.g., the Hermiticity reinstallation property) of its metric-operator-provided reinterpretation.

Interaction-picture context
Let us recall the well-known method of Seidewitz [114] for avoiding the no-go consequences of Haag's theorem in quantum field theory (QFT), which states that, under the usual assumptions made in QFT, any field that is unitarily equivalent to a free field must itself be a free field. In particular, Seidewitz shows how Haag's theorem can be avoided when QFT is formulated using an additional invariant path parameter, as well as the usual four-position parameters. His method relies on the removal of the spectral condition, essentially by replacing the usual on-mass-shell operator of relativistic energy and momentum by its off-shell generalization P. Importantly, this is accompanied by the preservation of the traditional IP approach to the Dyson perturbation expansions of scattering matrices, thereby providing a fully consistent basis for performing the usual practical QFT calculations.
For the purposes of our present discussion, what is important in the above Seidewitz construction [114] is that, in order to bypass the limitations imposed by Haag's theorem in the standard formulations of QFT, the frame-dependent zeroth component ofP, which plays the role of the energy operator in conventional QFT, necessarily becomes different from the newly introduced free relativistic Hamiltonian. From a rather abstract point of view the latter idea is reminiscent of and, indeed, completely paralleled by, our present introduction of the distinctions between several alternative candidates to play the role of a non-stationary non-Hermitian Hamiltonian.
As we have seen, only one of these candidates (viz. the IP operator H (t)) is hiddenly Hermitian (i.e., observable). The other ones [viz. the IP operators Ξ(t) and G(t)], together with their conjugate forms are not. In general, the differences [G † (t) − G(t)] and [Ξ † (t) − Ξ(t)] do not vanish. In turn, this opens up a wide space for multiple unconventional dynamical scenarios that provide an interplay and mutual cancellations between the non-unitarily evolving components in states ψ(t) and in observables Q(t), respectively. Along these lines, the ultimate (hidden) unitarity of the whole picture is achieved.
We described the corresponding non-Hermitian generalizations of evolution equations, thereby offering a firm ground for the transfer (NHIP → TD-CCM) of the operator NHIP formalism into its variational TD-CCM parallels. In the present stage of development of the theory we already understand why and how operators G (t) and G † (t) control the timeevolution of states, as well as why and when their Coriolis-operator partners Ξ(t) and Ξ † (t) play the role of the generators of time-evolution of all of the IP-represented observables. As two specific examples of where the opposite transfer (TD-CCM → NHIP) of experience from the TD-CCM to its NHIP counterpart might lead to further insights we recall the specific applications of the former to the Rabi Hamiltonian and the condensed Bose fluid, as discussed in Sect. 6. In the context of much recent research in the area that has been aimed mainly at the SP formalism, we have also clarified here the deep changes in the role of the so-called Hilbert space metric in the generic situations in which it is allowed to vary with time, Θ = Θ(t). Furthermore, having shown here how we can, completely equivalently, work instead solely with the pairs {|ψ(t) , |ψ(t) } of the single-state-representing kets and ketkets, we also thereby demonstrated how the explicit need for the operator Θ(t) itself has certainly been weakened, if not altogether eliminated. We believe that a transfer of this experience to the variational context is potentially extremely productive.

Summary
The formulation of the relationships illustrated in Eq. (59) was motivated, first of all, by the many phenomenological successes of the conventional stationary special case of Eq. (56). In the present paper, the emphasis was aimed rather at an innovative reinterpretation of the existing non-stationary (TD-CCM) extensions of the variational CCM techniques.
The underlying mathematical and technical details required to transfer the non-stationary Dyson-inspired formalism (viz. the NHIP approach) were outlined. It is precisely in this domain where a large part of our present contributions are truly original in their own right. The remaining originality resides in bringing into juxtaposition, within a newly developed conceptual and notational framework that is broad and powerful enough to encompass both, two very powerful methodologies that have hitherto been seen as quite separate. We have been at some pains to draw parallels that might henceforth be exploited to advance both formalisms and their subsequently enhanced arenas of applications. We have described in some detail how the key source of the mathematical inspiration for our study lies in the unconventional non-Hermiticity of Hamiltonians in the formalisms. This gave birth earlier [1] to the explicit, and rather fruitful, description of manifold conceptual parallels between two otherwise seemingly disparate theoretical constructs. These comprise, on one hand, the very successful stationary S-CCM formulations of quantum many-body theory and, on the other, their (at least, in principle) somewhat more ambitious, quasi-Hermitian [32][33][34]70] (otherwise known as pseudo-Hermitian [14] or PT -symmetric [13]) stationary analogues, which, by now, have themselves also become rather widely used in a variety of applications of quantum theory.
A non-stationary, TD-CCM version of the CCM theory has been considered here, therefore. This was done because we believe that, in contrast to the rather universal Dyson-inspired quasi-Hermiticity techniques, a characteristic and specific merit of the less universal CCM theory may be seen in the much larger number of very accurate calculations that it has found in such a wide range of applications in many diverse subfields of quantum physics and chemistry.
Several conclusions from our study also appear to provide a perceptibly deeper insight into abstract quantum theory itself. In particular, we have shown that one can easily remove the representation-framework restrictions as accepted both in Ref. [11] (wherein only wave functions evolved in time, i.e., in the non-Hermitian Schrödinger picture) and in Ref. [46] (where the transition to the non-Hermitian Heisenberg picture, in which wave functions remain constant in time, was analyzed). In other words, our present version of the NHIP formalism (in which both wave functions and operators of observables cease to remain constant in time) may be briefly characterized as an immediate non-Hermitian generalization of the interaction picture (alias the traditional Dirac picture) as described in most standard quantum mechanics textbooks (and see, e.g., Refs. [115][116][117]), and as now almost universally used in quantum field theory calculations (and see, e.g., Refs. [58,59,[118][119][120]).
One of the deep unifying features between the quasi-Hermiticity techniques on the one hand and the CCM on the other is their ability to be formulated in terms of a bivariational principle. For the general case of quasi-Hermitian operators in quantum mechanics, Scholtz et al. [32] showed in particular how the introduction of the metric operator was especially useful for the implementation of a variational principle, which could itself then be used for (approximate) calculational purposes. Exactly the same bivariational formulation of both S-CCM and TD-CCM approaches [2,5] has led to the extremely accurate descriptions of a wide variety of strongly-interacting quantum many-body systems.
Another common feature of the quasi-Hermiticity and CCM (particularly the ECCM) formalisms is their deep relationship to exact bosonization mappings. In the case of the ECCM this has even resulted in its ultimate realization as an exact classicization [3], as we discussed in some detail in Sect. 5.2. The ECCM itself, as we also discussed, introduces two independent preconditioning operators Ω of the form of Eq. (4). We now fully expect that this important feature of the ECCM might also find applications into a further useful, parallel extension of the three-Hilbert-space (alias the NHIP) formalism, which mirrors the NCCM → ECCM extension. However, such a discussion takes us far beyond the aims of the present paper.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.