Unified Dark Energy and Dust Dark Matter Dual to Quadratic Purely Kinetic K-Essence

We consider a modified gravity plus single-scalar-field model, where the scalar Lagrangian couples symmetrically both to the standard Riemannian volume-form (spacetime integration measure density) given by the square-root of the determinant of the Riemannian metric, as well as to another non-Riemannian volume-form in terms of an auxiliary maximal-rank antisymmetric tensor gauge field. As shown in a previous paper, the pertinent scalar field dynamics provides an exact unified description of both dark energy via dynamical generation of a cosmological constant, and dark matter as a"dust"fluid with geodesic flow as a result of a hidden Noether symmetry. Here we extend the discussion by considering a non-trivial modification of the purely gravitational action in the form of f(R) = R - \alpha R^2 generalized gravity. Upon deriving the corresponding"Einstein-frame"effective action of the latter modified gravity-scalar-field theory we find explicit duality (in the sense of weak versus strong coupling) between the original model of unified dynamical dark energy and dust fluid dark matter, on one hand, and a specific quadratic purely kinetic"k-essence"gravity-matter model with special dependence of its coupling constants on only two independent parameters, on the other hand. The canonical Hamiltonian treatment and Wheeler-DeWitt quantization of the dual purely kinetic"k-essence"gravity-matter model is also briefly discussed.


Introduction
Unified description of dark energy and dark matter as a manifestation of a single entity of matter has been an important challenge in cosmology in the last decade or so (for extensive reviews of dark energy see [1]- [3], and for reviews of dark matter see [4]- [6]).
Originally, unified treatment of dark energy and dark matter was proposed in the "Chaplygin gas" models [7]- [10]. Another trend aimed at unifying dark energy and dark matter is based on the class of "k-essence" models [11]- [14], in particular, on the so called "purely kinetic k-essence" models [15] (for further developments, see [16]- [20]), which successfully avoid difficulties inherent in the generalized Chaplygin gas models related to nonnegligible sound speed.
Also, recently a lot of interest has been attracted by the so called "mimetic" dark matter model proposed in [21,22]. The latter employs a special covariant isolation of the conformal degree of freedom in Einstein gravity, whose dynamics mimics cold dark matter as a pressureless "dust". Further generalizations and extensions of "mimetic" gravity are studied in Refs. [23,24].
Models of explicitly coupled dark matter and dark energy described in terms of two different scalar fields were proposed in Ref. [25].
In this paper we study a class of generalized models of gravity interacting with a single scalar field employing the method of non-Riemannian volume-forms on the pertinent spacetime manifold, i.e., generally covariant integration measure densities independent of the standard Riemannian one given in terms of the square root of the determinant of the metric [26]- [29]. (for further developments, see Ref. [30]). In this general class of models, also called "two-measure gravity theories", the non-Riemannian volume-forms are defined in terms of auxiliary maximal-rank antisymmetric tensor gauge fields ("measure gauge fields").
The introduction of the two integration measures (one standard Riemannian and another non-Riemannian one) opens the possibility to obtain both dark energy and dark matter from a single scalar field dynamics, as it was already observed in Ref. [31]. Subsequently, in Ref. [32] we have gone further and have discovered the fundamental reason that a class of models generalizing those studied in [31] describes a unification of dark matter and dark energy as an exact sum of two separate contributions in the pertinent energy-momentum tensor (see also the review in Section 2 below). This is because of: (i) The appearance of an arbitrary integration constant from a dynamical constraint on the scalar Lagrangian as a result of the equations of motion for the auxiliary "measure" gauge field. This integration constant is identified as a dynamically generated cosmological constant which provides the dark energy component.
(ii) The existence of a hidden Noether symmetry of the non-Riemannian-measure-modified scalar Lagrangian implying a conserved Noether current, which produces a dark matter component as a "dust" fluid flowing along geodesics.
This behavior is totally independent of specific form of the scalar field Lagrangian, as long as the scalar field couples in a symmetric way to both of the measures. The latter also ensures that the hidden "dust" Noether symmetry holds. In addition, the fact that the dynamically generated cosmological constant arises as an arbitrary integration constant makes the observed vacuum energy density totally decoupled from the parameters of the initial scalar Lagrangian.
In the present paper our main object of study is a non-trivial extension of the modified gravity plus single scalar field model considered in [32], describing unification of dark matter and dark energy. Namely, we now modify the purely gravitational part of the action by introducing an extended f (R)-gravity action with f (R) = R − αR 2 within the first-order Palatini formalism 1 .
These new non-Riemannian-measure-modified gravity--scalar-field theory has some very interesting features. On one hand, the hidden "dust" Noether symmetry of 1 Let us recall that R + R 2 gravity within the second order formalism (which was also the first inflationary model) was originally proposed in Ref. [33]. the scalar field Lagrangian remains intact and so does the picture of unified description through the pertinent energy-momentum tensor of both the dynamical dark energy (via the dynamical generation of the cosmological constant) as well as the "dust" dark matter fluid, i.e., they appear as an exact sum of two separate contributions to the corresponding energy density.
On the other hand, upon performing transition to the effective "Einstein frame" and upon appropriate scalar field redefinition we arrive at a dual theory of a specific quadratic purely kinetic "k-essence" form. In the dual purely kinetic "k-essence" formulation the original "dust" Noether symmetry is replaced by a simple shift symmetry of the transformed scalar field.
It is essential to stress on the following important properties of the dual theory: (a) All three constant coefficients in the quadratic purely kinetic "k-essence" action are given in terms of only two independent parameters: α (the R 2 coupling constant) and an arbitrary integration constant M produced by a dynamical constraint resulting from the equation of motion of the "measure" gauge field in the original theory.
(b) Let us emphasize that in the present context applying the notion of "duality" is justified in the usual sense of weak coupling in the original theory versus strong coupling in the dual theory. Indeed, as we will see below (Eqs. (39)) all coupling constants in the dual quadratic purely kinetic "k-essence" theory are functions of 1/α. (c) In the limit α → 0 (strong coupling limit in the dual "k-essence" theory) the physical quantitiesenergy density, pressure etc. have well-defined smooth limiting values in spite of the singularities in the coefficients of the "k-essence" action. The latter limiting values coincide with the corresponding values in the original gravity-matter theory with a standard Einstein-Hilbert gravity action (where α = 0). Thus, the established duality reveals an important feature of the approximate description of unified dark energy and dark matter via the quadratic purely kinetic "k-essence" model (Eq.(40) below) -it becomes an exact unification of both dark species in the strong coupling limit α → 0.
(d) In the limit α → 0 the dual theory parameter M precisely coincides with the dynamically generated effective cosmological constant in the original theory [32]. Thus, in this sense we can say that the dual purely kinetic "k-essence" theory is (partially) dynamically generated.
The above list of properties epitomizes the new features in our treatment of the purely quadratic "k-essence" theory w.r.t. to the previous treatment in [15]. Now, an important remark is in order. In a number of papers dealing with "dust-like" dark matter one uses the so called "Lagrangian multiplier gravity" formalism [34,35]. We would like to point out that this formalism is in fact a special particular case of the above mentioned more general and powerful approach based on non-Riemannian spacetime volume-forms (being used also in the present paper), which appeared earlier and which has profound impact in any (field theory) models with general coordinate reparametrization invariance, such as general relativity and its extensions [26]- [29], [30]- [32], [36]- [37]; strings and (higherdimensional) membranes [38] and supergravity [39] Indeed, dynamical constraints like the one on the scalar field Lagrangian in Eq.(8) below, which routinely appear in all instances of applying the non-Riemannian volume-form method in gravity-matter theories, resembles at first sight analogous constraints on scalar field Lagrangians in the "Lagrangian multiplier gravity" [34,35]. However, employing the concept of non-Riemannian volume-forms in the form of dual field strengths of auxiliary maximal-rank tensor gauge fields ("measure gauge fields") as in Eq.(3) below, instead of bare Lagrange multiplier fields, has certain essential advantages: (i) Dynamical constraints in "two-measure" gravitymatter theories result from the equations of motion of the auxiliary "measure" gauge fields and, thus, they always involve an arbitrary integration constant like M in Eq.(8) below, as opposed to picking some a priori fixed constant within the "Lagrange multiplier gravity" formalism. Depending on the specific gravity-matter theories with non-Riemannian volume-forms under consideration, the pertinent arbitrary integration constant acquires the meaning of a dynamically generated cosmological constant like the integration constant M below (cf. Refs. [26,39]) (ii) Employing the canonical Hamiltonian formalism for Dirac-constrained systems we find that the auxiliary "measure" gauge fields are in fact almost pure gauge degrees of freedom except for the above mentioned arbitrary integration constants which are identified with the conserved Dirac-constrained canonical momenta conjugated to the "magnetic" components of the "measure" gauge fields (see Appendix A in Ref. [37]).
(iii) Upon applying the non-Riemannian volumeform formalism to minimal N = 1 supergravity the dynamically generated cosmological constant triggers spontaneous supersymmetry breaking and mass generation for the gravitino (supersymmetric Brout-Englert-Higgs effect) [39]. Applying the same formalism to antide Sitter supergravity allows to produce simultaneously a very large physical gravitino mass and a very small positive observable cosmological constant [39] in accordance with modern cosmological scenarios for slowly expanding universe of the present epoch [1,2,3].
(iv) Employing two independent non-Riemannian volume-forms in generalized gravity-gauge+scalar-field models [36], thanks to the appearance of several arbitrary integration constants through the equations of motion w.r.t. the "measure" gauge fields, we obtain a remarkable effective scalar potential with two infinitely large flat regions (one for large negative and another one for large positive values of the scalar field ϕ) with vastly different scales appropriate for a unified description of both the early and late universe' evolution. An interesting feature is the existence of a stable initial phase of non-singular universe creation preceding the inflationary phase -stable "emergent universe" without "Big-Bang" [36].
The plan of the paper is as follows. In Section 2, for self-consistency of the exposition, we briefly review [32], namely, the basic aspects of the non-Riemannian volume-form approach leading to dynamical generation of a cosmological constant (dynamical dark energy) and revealing a hidden Noether symmetry allowing for a "dust" fluid representation of dark matter, such that both dark species appear as a sum of two separate contributions to the energy-momentum tensor. Sections 3 contains the main result -establishing duality (in the standard sense of weak versus strong coupling) between the original quadratic f (R)-gravity plus modified-measure scalar-field model, whose matter part delivers an exact unified description of dynamical dark energy and dust fluid dark matter, on one hand, and a specific quadratic purely kinetic "k-essence" gravity-matter model, on the other hand. In Section 4 we analyze the above models within the canonical Dirac-constrained Hamiltonian formalism and derive the corresponding quantum Wheeler-DeWitt equations in the cosmological Friedmann-Lemaitre-Robertson-Walker framework.

Non-Riemannian Volume-Form Formalism
We start by considering the following non-conventional gravity-scalar-field action -a particular case of the general class of the so called "two-measure" gravity-matter theories [27]- [29] (for simplicity we use units with the Newton constant G N = 1/16π): The notations used here and below are as follows: -The first term in (1) is the purely gravitational action in the first order (Palatini) formalism, where the Riemannian metric g µν and the affine connection Γ λ µν are a priori independent variables. In the previous paper [32] we have studied (1) with the simplest choice of S grav : where R denotes the scalar curvature -in this case the first order (Palatini) formalism is equivalent to the ordinary second-order (metric) formalism. -The second term in (1) -the scalar field action is constructed in terms of two mutually independent volume-forms: (a) √ −g ≡ − det g µν is the standard Riemannian integration measure density (spacetime volumeform); (b) Φ(B) denotes an alternative non-Riemannian generally covariant integration measure density defining an alternative non-Riemannian volume-form: where B µνλ is an auxiliary maximal rank antisymmetric tensor gauge field independent of the Riemannian metric. B µνλ (3) is called "measure gauge field". -L(ϕ, X) is general-coordinate invariant Lagrangian of a single scalar field ϕ(x), which can be of an arbitrary generic "k-essence" form [11,12,13,14]: i.e., a nonlinear (in general) function of the scalar kinetic term X.
In this Section we will concentrate on the scalar field action -the second term in (1). First, due to generalcoordinate invariance we have covariant conservation of the pertinent energy-momentum tensor: It follows from the equation of motion w.r.t. ϕ: Further, variation of the action (1) w.r.t. "measure" gauge field B µνλ reads: i.e., the B µνλ -equations of motion yield the following dynamical constraint on the scalar field Lagrangian: where M is arbitrary integration constant. The factor 2 in front of M is for later convenience, moreover, we will take M > 0 in view of its interpretation as a dynamically generated cosmological constant (see Eqs. (10) and (16), (19) below).
It is important to stress that the scalar field dynamics is determined entirely by the first-order differential equation -the dynamical constraint Eq.(8) which, in the simplest case of (4) (L(ϕ, X) = X − V (ϕ)) to be considered henceforth for simplicity, implies: The standard second order differential equation (6) is in fact a consequence of (8) together with the energymomentum conservation ∇ µ T µν = 0 with: The physical meaning of the "measure" gauge field B µνλ (3) as well as the meaning of the integration constant M are most straightforwardly seen within the canonical Hamiltonian treatment of (the scalar field part of) (1) -this is systematically derived in Sec.2 of Ref. [32]. Namely, using short-hand notations for the components of B µνλ (3): we obtain for the canonically conjugated momenta π B , π B i the set of Dirac first class constraints: where we have π B = L(ϕ, X) straightforwardly obtainable from the explicit form of the action (1) taking into account the representation of Φ(B) in (11). The above Dirac constraints imply that all components of 5 the "measure" gauge field B µνλ (12) are pure-gauge (non-propagating) degrees of freedom. The last relation in (13) is the canonical Hamiltonian analog of Eq. (8) within the Lagrangian formalism -in other words, the integration constant M in the B µνλ -equations of motion is (modulo a trivial numerical factor) a conserved Dirac-constrained canonical momentum π B conjugated to the "magnetic" component B of the "measure"-gauge field B µνλ (12).
For more details about the canonical Hamiltonian treatment of general gravity-matter theories with (several independent) non-Riemannian volume-forms we refer to [37].

"Dust" Fluid Conservation Laws
In Ref. [32] we have shown that the scalar field action in (1) possesses a hidden Noether symmetry, namely (1) is invariant (up to a total derivative) under the following nonlinear symmetry transformations: Then, the standard Noether procedure yields the conserved current: Let us stress at this point, that the existence of the hidden symmetry (14) of the action (1) does not depend on the specific form of the scalar field Lagrangian (4). The only requirement is that the kinetic term X must be positive.
T µν (5) and J µ (15) can be cast into the a relativistic hydrodynamical form (taking into account (8)): where: For the pressure p and energy density ρ we have accordingly: where the integration constant M appears as dynamically generated cosmological constant.
Let us note that constancy (19) of the pressure p = −2M in (16) together with the covariant conservation of T µν : upon projecting (21) along the "velocity" vector u µ and orthogonally w.r.t. the latter by Π µλ = g µλ + u µ u λ ), immediately implies both the covariant conservation of J µ (15): as well as the geodesic flow equation: As discussed in [32] the energy-momentum tensor (16) consists of two parts with the following interpretation according to the standard Λ-CDM model [40,41,42] (using notations p = p DM + p DE and ρ = ρ DM + ρ DE in (19)-(20)): -Dark energy part given by the first cosmological constant term in T µν (16), which arises due to the dynamical constraint on the scalar field Lagrangian (8) with p DE = −2M , ρ DE = 2M ; -Dark matter part given by the second term in (16) with p DM = 0 , ρ DM = ρ 0 (ρ 0 as in (17)), which in fact describes a dust fluid. According to the general definitions, see i.e. [43], ρ DM = ρ 0 (17) and ρ DE = 2M are the rest-mass and internal fluid energy densities, so that the Noether conservation law (15) describes dust dark matter "particle number" conservation.
3 Quadratic Gravity Interacting with a Dark Energy -Dark Matter Unifying Scalar field -Dual to Purely Kinetic "K-Essence"

Derivation of the Dual Kinetic Pure "K-Essence" Theory
Let us now consider a modification of the gravitational part of the gravity-scalar-field action (2) as follows: where we have introduced f (R) = R − αR 2 extended gravity action in the first-order Palatini formalism: i.e., with a priori independent metric g µν and affine connection Γ µ νλ . Clearly, since the scalar field action -the second term in (24) -remains the same as in the original action (2), all results in subsection 2.2 remain valid 2 . In other words, the modified gravity-scalar-field action (24) possesses "hidden" Noether symmetry (14) producing "dust" fluid energy density conserved current (22) and interpretation of ϕ as describing simultaneously dark energy and dust dark matter with geodesic dust fluid flow (23).
The gravitational equations of motion resulting from (24) are, however, not of the standard Einstein form. The equations of motion w.r.t. g µν read: with T µν the same as in (10). Trace of Eqs.(26) yields: Variation of (24) w.r.t. Γ µ νλ yields: which shows, following the analogous derivation in the Ref. [27], that Γ µ νλ becomes a Levi-Civita connection: w.r.t. to the Weyl-rescaled metric g µν : Before going over to the physical "Einstein-frame" it is useful to perform the following ϕ-field redefinition: where the last relation follows from the Lagrangian dynamical constraint (9) together with (32).  (29), which together with (34) imply: as well as Eqs.(31)- (32) and (34) we can rewrite all equations of motion resulting from (24), in particular the quadratic f (R)-gravity Eqs. (26), in terms of the new metric g µν (32) and the new scalar field ϕ (33) in the standard form of Einstein gravity equations: with the following notations: -Here R µν and R are the standard Ricci tensor and scalar curvature of the Einstein-frame metric (32). -The Einstein-frame energy-momentum tensor: is given in terms of the following effective ϕ-scalar field Lagrangian of a specific quadratic purely kinetic "k-essence" form: Let us stress that the three constant coefficients in (38) depend only on two independent parameters (α, M ), the second one being a dynamically generated integration constant in the original theory (24). Thus, we have established a duality between the modified-measure gravity-scalar-field theory (24) within the original g µν -frame and the special quadratic purely kinetic "k-essence" theory within the conformally-rescaled g µν -frame (Einstein-frame): with a matter Lagrangian (38)- (39). The Einstein-frame effective energy-momentum-tensor (37) in the perfect fluid representation reads (taking into account (38)-(39)): Because of the obvious Noether symmetry of (40) under constant shift of ϕ: ϕ → ϕ + const (44) the corresponding Noether conservation law is identical to the ϕ-equations of motion: where ∇ µ indicates covariant derivative w.r.t. Levi-Civita connection in the g µν -(Einstein) frame. Introducing now the standard thermodynamical notions of enthalpy per unit particle h and particle number density n within the above Einstein-frame system, where ρ + p = h n, and taking into account the expressions (42)-(43), the Noether current conservation law (45) and the covariant conservation ∇ ν T µν = 0 of (41) can be written, respectively, as: where: Comparing relations (46) in the Einstein frame with the corresponding relations in the original g µν -frame (22)-(23) we conclude that: -Conservation of the dust dark matter energy density current (22) in the g µν -frame is dual to the conservation of the particle number density current within the Einstein frame -first Eq. (46), which in fact is the standard ϕ-equation of motion resulting from the action (40). That is, the "hidden" nonlinear Noether symmetry of (24) is dual to the shiftsymmetry (44) of (40). -In the original g µν -frame the dust dark matter flows along geodesics (23), whereas in the dual Einstein frame the dual purely kinetic "k-essence" fluid does not any more flow along geodesics (second Eq. (46)).
The purely kinetic "k-essence" theory (40) apart from the trivial vacuums ϕ = const possesses in addition a non-trivial "kinetic vacuum" solution X vac of the equations of motion (45): implying the dark energy property (cf. (41)): with effective cosmological constant: The explicit form of (48) reads: It has the form of the standard Hamilton-Jacobi equation for the action of a massive relativistic point-particle moving in a g µν -background with mass squared: In other words ϕ vac (x) = m 0 T , where T is the propertime for the particle to reach the spacetime point x from some reference point x (0) .

FLRW Reduction of the Dual Kinetic Pure "K-Essence" Theory
Let us now consider a reduction of the dual quadratic purely kinetic "k-essence" gravity-scalar-field model (40) for the Friedman-Lemaitre-Robertson-Walker (FLRW) class of metrics: Then the action (40) acquires the form (using again short-hand notations (39); in what follows we will take the spatial FLRW curvature K = 0 for simplicity): with . φ≡ d ϕ dt . The φ ≡ ϕ-equation of motion from (54) yields (using henceforth the gauge N = 1): where p φ is the constant conserved canonically conjugated momentum of φ ≡ ϕ. The equation of motion w.r.t. a resulting from (54) -the Friedmann equation -reads: with ρ(p φ /a 3 ) denoting the energy density as function of p φ a 3 : 8 where in the second line of (57) relation (55) was used. Here . φ(p φ /a 3 ) ≡ y is one of the roots of the cubic equation (55): which explicitly read: where: Similarly, for the pressure we have: where again in the second line the cubic equation (58) has been used. The Friedmann equation (56) can be solved approximately for small and large values of a using expressions (59)-(60): Accordingly, we have for the energy density (57): i.e. radiation domination for small a, and (using notations (50), (52)): i.e., dark energy domination plus a subleading "dust" dark matter contribution for large a. An important property of the above derivation of the duality between the original g µν -frame quadratic f (R)-gravity plus non-Riemannian-modified-measure scalar field action (24), on one hand, and the Einstein-frame special quadratic purely kinetic "k-essence" theory (40), on the other hand, is that there exists a smooth limit α → 0 of the energy density (57) and pressure (61) in the latter theory in spite of the singularity at the strong coupling limit α → 0 in all kinetic "k-essence" coefficients (38)- (39). The corresponding limiting values at α = 0 of the energy density and pressure of the dual purely kinetic "k-essence" theory (40) are those of the original theory (24) with the standard Einstein gravity action (α = 0), i.e., the theory (2) [32]. In particular, in the limit α → 0 we get precisely the expression for the energy density being an exact sum of dark energy and dust dark matter contributions produced by (2) [32].
The cosmological implications of general purely kinetic "k-essence" models have been previously studied extensively in Refs. [15]- [20]. In particular, an important strong inequality (Eqs.(26)- (27) in [15]) involving the parameters of a generic quadratic purely kinetic "kessence" theory was derived from the requirement that the onset of dark matter behavior must occur before the epoch of equal matter and radiation, and also that in the present stage of evolution the dark energy component must exceed twice the dark matter component in the pertinent energy density, The counterpart of the above Scherrer's inequality in the present settings becomes the requirement: which completely conforms to the result (68). From Eqs.(57)-(61) we obtain for the squared sound speed: with ρ as in (57)-(58). Obviously c 2 s → 0 for α → 0, which conforms to the limiting values (68)-(69). Also, c 2 s (71) is well-defined (positive) and reasonably small for ρ > 2M (1 − 8αM ) −1 , which is obviously fulfilled for small α according to (68).
The new aspects in our current treatment w.r.t. Refs. [15]- [20] are as follows: 9 -We have derived an explicit duality between the original g µν -frame quadratic f (R)-gravity plus non-Riemannian-modified-measure scalar field action (24), on one hand, and the Einstein-frame special quadratic purely kinetic "k-essence" theory (40). -Unlike the general purely kinetic "k-essence" treatment [15]- [20] all three coefficients in our kinetic "kessence" action (38)-(39) explicitly depend on only two independent parameters (α, M ), the latter being a dynamically generated integration constant. -In spite of the singularity of all "k-essence" coefficients A 1 , A 2 , A 3 (39) for small α (the coupling constant of the R 2 in the original theory (24)), the energy density and pressure in the dual purely kinetic "k-essence" theory have smooth limit for α → 0, whereby their limiting values (68) and (69) exactly coincide with the corresponding values in the original theory (24) with α = 0, i.e., (2) [32]. -The established duality between (24) and (40) elucidates the origin of the unified description of dark energy and dark matter within the approach based on "purely kinetic k-essence" models [15].

Canonical Hamiltonian Formalism and Wheeler-DeWitt Equation
For the canonical Hamiltonian formalism applied to generic generalized gravity-matter models with one or more non-Riemannian spacetime volume-forms we refer to [37]. In particular, for the theory (2) the canonical Hamiltonian analysis was discussed in detail in Sec.2 of Ref. [32].
Here we will consider specifically the Hamiltonian treatment and quantization of the action (54) -reduction of the purely kinetic "k-essence" model (40) for the Friedman-Lemaitre-Robertson-Walker (FLRW) class of metrics (53).
From the explicit form of the FLRW action (54) we deduce the canonically conjugated momenta p a , π N and p φ w.r.t. a, N and φ ≡ ϕ: where the second relation for π N is a primary Dirac first-class constraint. The total canonical Hamiltonian becomes: with ρ(p φ /a 3 ) as in (57), thus H total by itself is a secondary Dirac first-class constraint with N playing the role ot its Lagrange multiplier. Since φ ≡ ϕ is a cyclic variable its canonically conjugated momentum p φ is conserved. Quantization according to Dirac constrained Hamiltonian formalism proceeds by imposing the quantized operator version of the Hamiltonian constraint -the expression in the square brackets in (73) -on the quantum wave function Ψ (a, p φ ) -the Wheeler-DeWitt equation. The ordering ambiguity in the quantized version of the first term there is resolved by changing variables: and taking the special operator ordering: Therefore, the Wheeler-DeWitt equation acquires the form of Schrödinger equation for zero energy eigenvalue: with effective potential: (ρ(p φ /a 3 ) as in (57), and a and a related as in (74)). Explicitly: where f ( a, p φ ) ≡ y is a root of the cubic equation (cf. Eqs.(55), (58)-(59): Using the asymptotics (64)-(65) we obtain for the small a and large a of V eff ( a, p φ ): for small a, and: for large a. Using (79) we find that V eff ( a, p φ ) (78) is negative for all a ∈ (0, ∞) and diverges to −∞ at both ends of the interval according to (80)-(81).
According to Scherrer's inequality (70) α must be very small. Thus, substituting (68) into (77) and accounting for the change of variables (74) brings the Wheeler-DeWitt effective potential to the following form: valid for any a ∈ (0, ∞). In other words, for small α (82) is a sum of inverted harmonic oscillator potential with negative frequency squared: where Λ eff is the effective cosmological constant (50) for small α, plus an inverse square potential and with an "energy" eigenvalue E = √ 2p φ (1 + 4αM ). Thus, in the limit α → 0 the Wheeler-DeWitt equation (76) becomes Schrödinger equation for an inverted harmonic oscillator with negative frequency squared (83) and with "energy" eigenvalue E = √ 2p φ : The inverted harmonic oscillator was extensively studied in Ref. [44] (for a more recent account and further references, see [45]). In particular, the inverted oscillator was applied in [46] to study the quantum mechanical dynamics of the scalar field in the so called "new inflationary" scenario. Since the energy eigenvalue spectrum of the inverted harmonic oscillator is continuous (E ∈ (−∞, +∞)) and the corresponding energy eigenfunctions are not square-integrable, its application in the context of cosmology [46] required employment of wave-packets w.r.t. E instead of energy eigenfunctions.
Similarly, in the present case the value of p φ -the conserved canonical momentum of the kinetic "k-essence" scalar field φ ≡ ϕ in (54) -is the analog of energy eigenvalue E in the Wheeler-DeWitt Schrödinger-like equation (84) (modulo the factor √ 2). Therefore, now in the strong coupling limit α → 0 the "k-essence" field φ ≡ ϕ plays the role of a Wheeler-DeWitt time τ ≡ 1 √ 2 φ and Eq.(84) acquires the form of a "time-dependent" Schrödinger-like equation for the inverted harmonic oscillator upon Fourier transforming the "energy"-eigenvalue Eq.(84): Following [46] the appropriate normalized to unity (on the semiaxis a ∈ (0, ∞)) wave-packet solution of the "time-dependent" Wheeler-DeWitt equation (85) is of the form (using notation ω (83) for α = 0): where b is an integration constant describing the width of the wave packet. Accordingly, the average value of the FLRW scale factor a = √ 3 4 a 2/3 (cf. (74)-(75)) is given by: exhibiting no singularity ( a → 0) at any "time" τ .

Conclusions
In the present paper we have discussed in some detail the main properties of a generalized model of gravity interacting with a single scalar field, where we have employed the method of non-Riemannian spacetime volumeforms (alternative generally-covariant integration measure densities) constructed in terms of auxiliary maximal rank tensor gauge fields ("measure" gauge fields).
In the preceding paper [32] we have shown that the non-Riemannian-measure-modified scalar field action (2) yields a simple unified description of dark energy and dust dark matter. Namely, the corresponding energy density arises as an exact sum of a dark energy component in the form of a dynamically generated cosmological constant appearing as an arbitrary integration constant in the solution of the "measure" gauge field equations, and a dark matter component produced by a hidden Noether symmetry (not affecting the gravity part) giving rise to a Noether conserved current, which identifies the scalar field dynamics as a dust fluid motion along geodesics.
Here we extended the above treatment by coupling the non-Riemannian-measure-modified scalar field dynamics to quadratic f (R) gravity. We have found an explicit duality in the usual sense of "weak versus strong coupling" between the original non-standard gravityscalar-field model providing exact unified description of dynamical dark energy and dust fluid dark matter in the matter sector, on one hand, and a quadratic purely kinetic "k-essence" gravity-matter model, on the other hand. The latter dual theory arises as the "Einsteinframe" theory of its original counterpart. It is special in a sense that the couplings in the dual quadratic kinetic "k-essence" action are given in terms of only two parameters (α, M ) -the R 2 -coupling constant in the original action (1) and a dynamically generated integration constant M upon solving the equations of motion for the auxiliary "measure" gauge field in (1). Moreover, in spite of the divergence of the "k-essence" coupling constants when α → 0 (strong coupling limit), both the "k-essence" energy density and "k-essence" pressure have smooth limits at α = 0 with the limiting values coinciding with their respective values in (2)the weak coupling limit of the original theory (1), delivering an explicit unified description of dark energy and dust dark matter as an exact sum of two separate contributions to the total energy density.
The established duality in the present paper explains the ability of the purely kinetic "k-essence" models [15] to provide approximately a unified description of dark energy and dark matter and reveals that this unified description becomes exact in the strong coupling limit of a special type of quadratic purely kinetic "kessence" theory.
Finally, we have used the standard Dirac approach to constrained Hamiltonian systems for a canonical Hamiltonian treatment of the non-Riemannian-measure-modified gravity-scalar-field theory, specifically for the reduction of the latter in the case of FLRW class of cosmological spacetime metrics. In the limit of vanishing R 2 coupling the associated Wheeler-DeWitt equation acquires the form of a Schrödinger-like equation with an effective potential of inverted harmonic oscillator. The quantum average value of the FLRW scale factor does not exhibit any singularities in its time evolution.