Quantum billiards with branes on product of Einstein spaces

We consider a gravitational model in dimension D with several forms, l scalar fields and a Lambda-term. We study cosmological-type block-diagonal metrics defined on a product of an 1-dimensional interval and n oriented Einstein spaces. As an electromagnetic composite brane ansatz is adopted and certain restrictions on the branes are imposed the conformally covariant Wheeler-DeWitt (WDW) equation for the model is studied. Under certain restrictions, asymptotic solutions to the WDW equation are found in the limit of the formation of the billiard walls. These solutions reduce the problem to the so-called quantum billiard in (n + l - 1)-dimensional hyperbolic space. Several examples of quantum billiards in the model with electric and magnetic branes, e.g. corresponding to hyperbolic Kac-Moody algebras, are considered. In the case n=2 we find a set of basis asymptotic solutions to the WDW equation and derive asymptotic solutions for the metric in the classical case.


Introduction
In this paper we deal with the quantum billiard approach to D -dimensional cosmological-type models defined on a (warped) product manifold (u − , u + )× M 1 ×· · ·×M n , where M i is a smooth oriented Einstein manifold, i = 1, . . . , n .
The billiard approach in classical gravity originally appeared in the dissertation of Chitré [1] for an explanation of the BKL oscillations [2,3] while approaching to a spacelike singularity in the Bianchi-IX model [4]. In this description a simple triangle billiard in the hyperbolic (Lobachevsky) space H 2 was used. The BKL-like behavior near a timelike singularity was studied in [5,6].
In [7], the billiard approach for D = 4 was extended to the quantum case (see also [8]), i.e. the solutions to the Wheeler-DeWitt (WDW) equation [9] were reduced to the problem of finding the spectrum of the Laplace-Beltrami operator on Chitré's triangle billiard.
Chitré's billiard approach was generalized to a multidimensional cosmological model with multicomponent anisotropic ("perfect") fluid [10,11,12] defined on the product of n Einstein factor spaces and 1-dimensional manifold. The search for an oscillating behavior near the singularity was reduced to the problem of proving the finiteness of the billiard volume. At that time the quantum billiard approach to multidimensional cosmology was suggested in [12,13,14].
The classical billiard approach to multidimensional models with fields of forms and scalar fields in the presence of the Λ -term was suggested in [15] along lines suggested earlier in [12]. In ref. [15] rather a general composite electromagnetic ansatz for the fields of forms on a warped product of several Einstein manifolds and an 1-dimensional base manifold M 0 was developed. Reference [15] contained rather a general construction of the billiard approach for the description of the behavior of scale factors and scalar fields near either a spacelike or a timelike singularity, i.e. the metric in [15] ds 2 = wdu 2 +... , contained an arbitrary sign w = ±1 and a coordinate u . Thus, the paper [15] was dealing with cosmological-type solutions, e.g. cosmological, spherically symmetric, and cylindrically symmetric ones. The metric had a block-diagonal form. In [15] the necessary condition for the formation of walls was formulated in terms of inequalities for scalar products of the brane vectors U s : (U s , U s ) > 0 and for the so-called brane sign parameters: ε s > 0 . Inequalities on Kasner parameters, where formulated in terms of linear functions U s (α) , which give either Kasner or oscillatory asymptotic regimes near the singularity. Another advantages of the approach of [15] was in dealing with a wide variety of signatures of Einstein factor space metrics (though restricted by ε s > 0 ). It was shown that the curvatures of the Einstein factor spaces and the Λ -term are irrelevant near the singularity.
Meanwhile the approach of ref. [15] had some restrictive points, since it was dealing with block-diagonal metrics and putting restrictions on brane intersections (for branes corresponding to the same form field) which guaranteed block-diagonal structure of the stress-energy tensor. For some extension of these restrictions see ref. [16].
Some problems of the approach of ref. [15] were overcome in the papers of Damour, Henneaux and Nicolai [17,18,19,20,21] and some other authors. These works were aimed from the very beginning at studying the generic behavior of solutions near a spacelike singularity (a lá BKL) for gravitational and cosmological models with non-diagonal metrics, fields of forms and scalar fields. This approach was based on a wide use of Iwasawa decomposition and hyperbolic Kac-Moody algebras [22,23,24,25]. It was shown in [20] that for certain models (of supergravity) the billiards (or their parts) are related to Weyl chambers of certain hyperbolic Kac-Moody (KM) algebras. This fact has simplified the proof of the finiteness of the billiard volume in certain cases.
In the recent publications [26,27,28] the quantum billiard approach for the multidimensional gravitational model with several forms was considered. The main motivation for the quantum billiard approach in [26,27] was coming from the quantum gravity paradigm; see [29] and references therein. The asymptotic solutions to the WDW equation from [26,27] (in the model without scalar fields) are similar to those obtained earlier in [12] for a multicomponent anisotropic fluid with certain equations of state. In [28] another, conformally covariant form of the WDW equation [30,31,32,33,34] was used. In this case the minisuperspace was enlarged by including the form potentials for electric non-composite branes. In [28] an example of a 9dimensional quantum billiard for D = 11 model with 120 four-forms which mimic spacelike M2 -brane solutions ( SM2 -branes in D = 11 supergravity) was considered. It was shown that the wave function vanishes as y 0 → −∞ (i.e. at the singularity), where y 0 is the "tortoise" timelike coordinate in minisuperspace [28]. In [35] we have generalized the approach of [28] to the case when scalar fields with dilatonic couplings were added into consideration and the composite electromagnetic ansatz for branes was considered instead of the non-composite electric one from [28]. New examples of quantum billiards with electric and magnetic S -branes in D = 11 and D = 10 models were presented. It was found that in the quantum case adding of magnetic branes changes the asymptotic behavior of the wave functions, while it does not change the classical asymptotic oscillating behavior of the scale factors (and scalar field for D = 10 ). It was found that in certain examples the basis wave functions in the "tortoise gauge" vanish as y 0 → −∞ .
In this paper we generalize the approach from [35] to n Einstein factor spaces and a Λ term. We also extend this approach by relaxing the main restriction for brane vectors: (U s , U s ) > 0 . Here we consider examples of billiards in the model with n non-intersecting electric branes, n ≥ 2 . The brane world volumes are volume forms of M i . We show that in the classical case any of these billiards describe the never ending oscillating behavior of the scale factors while approaching a singularity, which may be either spacelike or timelike one. The examples with timelike singularities are supported by the use of either phantom form fields, or extra time-variables. For n = 2 (when 1-dimensional M 1 and M 2 are forbidden) we obtain the basis asymptotic solutions to the WDW equation.
We also generalize the model with n electric branes by adding a form of rank D . This adding does not change the billiard but it changes (e.g. drastically for D ≤ 7 ) the basis asymptotic solutions to the WDW equation for a certain choice of Bessel function.
Here we also consider an example of 4-dimensional quantum billiard in D = 11 model with ten 4-forms [15]. We use the ansatz with ten magnetic non-composite branes with brane worldvolumes of the form We prove the vanishing of the basis wave functions in the "tortoise gauge" as y 0 → −∞ .

The model
Here we study the multidimensional gravitational model governed by the action is an n a -form ( n a ≥ 2 ) on M and λ a is a 1 -form on R l : λ a (ϕ) = λ aα ϕ α , a ∈ ∆ , α = 1, . . . , l . In (2.1) we denote |g| where ∆ is some finite set of (color) indices and S YGH is the standard (York-Gibbons-Hawking) boundary term. In the models with one time and the usual fields of forms all θ obey θ a > 0 when the signature of the metric is (−1, +1, . . . , +1) . For such a choice of signature θ b < 0 corresponds to a "phantom" form field F b . We consider the manifold with the metric where R * = (u − , u + ) , w = ±1 and g i = g i m i n i (y i )dy m i i ⊗ dy n i i is an Einstein metric on M i satisfying the equation . . , n . The functions γ, β i : R * → R are smooth. We denote d i = dimM i ; i = 1, . . . , n and d 0 = 1 ; D = n ν=0 d ν . We put any manifold M i , i = 1, . . . , n , to be oriented and connected. Then the volume d i -form and signature parameter are correctly defined for all i = 1, . . . , n . The cosmological ( S -brane) solutions correspond to w = −1 and positive definite g i for all i , while static configurations (e.g. fluxbranes, wormholes, black branes etc.) may be obtained when w = 1 , g k are Riemannian metrics for all k > 1 and g 1 is the metric of pseudo-Euclidean signature (−, +, ..., +) . Here we may also deal with solutions having several timelike directions.
Here and below the set S consists of elements where a s ∈ ∆ is the color index, v s = e, m is the electromagnetic index, and the set I s ∈ Ω a s ,v s describes the location of the brane. Due to (2.12) and (2.13) we get d(I) = n a − 1, d(J) = D − n a − 1 .
Here we present two restrictions on the sets of branes which guarantee the diagonal form of the energy-momentum tensor [36].
The first restriction for a pair of two (different) branes both electric ( eepair) or magnetic ( mm -pair) with coinciding color index reads for any I, J ∈ Ω a,v , a ∈ ∆ , v = e, m (here d(I) = d(J) ). The second restriction for any pair of two branes with the same color index, which include one electric and one magnetic brane ( em -pair) has the following form: where I ∈ Ω a,e , J ∈ Ω a,m , a ∈ ∆ . These restrictions are satisfied identically in the non-composite case, when there are no two branes corresponding to the same form F a for any a ∈ ∆ . It follows from [36] that the equations of motion for the model (2.1) and the Bianchi identities, dF s = 0 , s ∈ S m , for fields from (2.4), (2.11)-(2.14), when restrictions (R1) and (R2) are imposed, are equivalent to the equations of motion for the σ -model governed by the action where X = (XÂ) = (β i , ϕ α , Φ s ) ∈ R N , N = n + l + m , m = |S| is the number of branes and the minisupermetric G = GÂB(X)dXÂ ⊗ dXB on the minisuperspace M = R N is defined as follows: is the truncated target space metric with 21) and the co-vectors We denote χ e = +1 and χ m = −1 ; is the indicator of i belonging to I : δ iI = 1 for i ∈ I and δ iI = 0 otherwise; and In the electric case (F (a,m,I) = 0) for finite internal space volumes V i the action (2.18) coincides with the action ( In what follows we will use the scalar products of is the matrix inverse to the matrix (2.20). Here (as in [32]) The scalar products (2.26) read [36] ( The potential (2.23) reads as follows: The scalar products of co-vectors U Λ , U (j) , U s are defined by the following relations [36] ( The vector U Λ is a timelike as well as U (i) with d i > 1 (here we deal with U (i) obeying ξ i = 0 ). The vectors U Λ and U (i) with d i > 1 belong to the same light cone (interiour part) due to relations (U (i) , U Λ ) = −1 .

Quantum billiard approach
Here we generalize the quantum billiard approach for asymptotic solutions to the Wheeler-DeWitt (WDW) equation from [35] to the case of a chain of Einstein spaces in the presence of the Λ -term.
Let us denote by S + the subset of all s ∈ S obeying First we put the following additional restrictions on the model: These restrictions are necessary conditions for the formation of infinite "wall" potential in hyperbolic spaces in certain limit (see below). The first restriction excludes phantom scalar fields. We note that in our previous work [28,35,37] we used a more rigid restriction: S = S + .
By fixing the temporal gauge: where f : M → R is a smooth function, we obtain the Lagrange system with the Lagrangian and the energy constraint The set of Lagrange equations with the constraint (3.6) is equivalent to the set of Hamiltonian equations for the Hamiltonian with the constraint H f = 0, (3.8) where PÂ = µe 2f GÂB(X)ẊB are momenta (for fixed gauge) and (GÂB) = (GÂB) −1 .
Here we use the prescriptions of covariant and conformally covariant quantization of the hamiltonian constraint H f = 0 which was suggested initially by Misner [30] and considered afterwards in [31,32,34] and some other papers.
We obtain the Wheeler-DeWitt (WDW) equation, is the wave function corresponding to the f -gauge (3.4) and satisfying the relation In (3.9) we denote by ∆[G f ] and R[G f ] the Laplace-Beltrami operator and the scalar curvature corresponding to the metric respectively. The Wheeler-DeWitt (WDW) equation (3.9)) is conformally covariant. This follows from (3.11) and the relation: where the coefficients a N and b N are well-known in the conformally covariant theory of scalar field. Now we put f = f (σ) and denotē Here we deal with a special class of asymptotic solutions to the WDW equation. Due to restrictions (3.2) and (3.3) the (minisuperspace) metricŝ G , G have pseudo-Euclidean signatures (−, +, . . . , +) . We put and we imposeĜ AB σ A σ B < 0 . With this choice we deal with the so-called "tortoise" time gauge.
Here we use a diagonalization of σ -variables We restrict the WDW equation to the lower light cone V − = {z = (z 0 , z)|z 0 < 0, η ab z a z b < 0} and we introduce the Misner-Chitré-like coordinates where y 0 < 0 and y 2 < 1 .
In these variables we have f 0 = y 0 . In the following we usē where For the wave function we consider the ansatz from [35]: where where Q s = 0 and e iQ s Φ s = exp(i s∈S Q s Φ s ) . Repeating all calculations from [35] we get Here and in what follows U s (σ) = U s A σ A . Now we proceed with the study the asymptotic solutions to WDW equation in the limit y 0 → −∞ . Due to (3.22) and (3.24) this equation reads Here and in what follows we put Q s = 0 for all s ∈ S + . It follows from the analysis of [15] that for a certain choice of diagonalization (3.17), as Here V ∞ is the potential of infinite walls which are produced by branes with (U s , U s ) > 0 : where we denote θ ∞ (x) = +∞ , for x ≥ 0 and θ ∞ (x) = 0 for x < 0 . The vectors v s , s ∈ S + , which belong to where the N 0 -dimensional vectors u s = (u s 0 , u s ) = (u s a ) are obtained from brane U s -vectors using the matrix (S A a ) from (3.17) for all s ∈ S + . In what follows we use a diagonalization (3.17) obeying for all s ∈ S and all i such that ξ i = 0 (and hence . . , n . The diagonalization (3.17) from [15] obeys these conditions and implies (3.28) and (3.29).
The potential V ∞ corresponds to the billiard B in the hyperbolic space (D N 0 −1 , h L ) . This billiard is an open domain in D N 0 −1 obeying the a set of inequalities:

36)
s ∈ S + . The boundary of the billiard ∂B is formed by parts of hyper-spheres with centers in v s and radii r s . The conditions (3.34) are obeyed for the diagonalization (3.35) with where e = (e A ) is a normed timelike vector (e, e) = −1 obeying (e, U Λ ) < 0 and (e, U s ) < 0 for all s ∈ S . Hence (e, U (i) ) < 0 for all i obeying ξ i = 0 . Our choice in [15] was e = U Λ / |(U Λ , U Λ )| . When all factor spaces M i are Ricci-flat, i.e. all ξ i = 0 , brane part of conditions (3.34) may be relaxed, while the curvature part of these conditions should be omitted. In this case we obtain a more general definition of the billiard walls (e.g. for u s0 ≤ 0 ) described in [38].
Thus, as in [35], we are led to the asymptotic relation for the function Ψ 0,L (y 0 , y) By separating the variables, we obtain the following asymptotic relation (for y 0 → −∞ ) We assume that the minus Laplace-Beltrami operator (−∆ L ) with the zero boundary conditions has a spectrum obeying the following inequality: This restriction was proved in [26,27] for a wide class of billiards with finite volumes.
Here we restrict ourselves to the case of negative A -number A < 0 . Solving eq. (3.40) we get for A < 0 the following set of basis solutions: where B iω (z) = I iω (z), K iω (z) are the modified Bessel functions and We denote In the following we impose the restriction on U = (U A ) : (U, U ) < 0 . We have (U, U Λ ) < 0 due to (U s , U Λ ) < 0 for all s .
From now we use a diagonalization with z -variables obeying (3.37) with e = U/ |(U, U )|. (3.46) For U = kU Λ , with k > 0 , such a diagonalization coincides with that of ref. [15]. It was obtained in [35] that as y 0 → −∞ for any fixed y ∈ B and C ± are non-zero constants, "plus" corresponds to B = I and "minus" -to B = K . Here and Now we outline our analysis from [35] of asymptotic behavior of Ψ f 0 as y 0 → −∞ . Here we fix all Φ s ∈ R , s ∈ S .
For B = K , Ψ f 0 → 0 as y 0 → −∞ for fixed y ∈ B . Now let B = I . If 1 2 q > 2|A| , or, equivalently, we also get Ψ f 0 → 0 as y 0 → −∞ for fixed y ∈ B . For 1 2 q = 2|A| , or, equivalently, we get Ψ f 0 → 0 as y 0 → −∞ for fixed y ∈ B \ { 0} . For y = 0 we get |Ψ f 0 | → +∞ as y 0 → −∞ . It may be shown that in this case, when m = N 0 = n + l (i.e. if m is the minimal number of walls which is necessary for the billiard to have a finite volume) we get as y 0 → −∞ , where C 0 = 0 is a constant irrelevant for our consideration. Thus, for m = n + l and for eigenfunction Ψ L ( y) with Ψ L ( 0) = 0 we get a δ -function in the asymptotic of Ψ f 0 . In this case we have an asymptotic localization of Ψ f 0 at the point y = 0 for our choice of gauge ("tortoise" one). When the scalar fields are absent and we use a diagonalization from [32,39] the relation y = 0 implies the isotropization β i = β and we may talk in terms of asymptotic quantum isotropization of the wave function in the temporal gauge under consideration. When 1 2 q < 2|A| , or, equivalently, we get |Ψ f 0 | → +∞ as y 0 → −∞ for y belonging to the open domain Outside the closure of B ∞ we get the zero limit of our wave function and we may talk in terms of the asymptotic localization of Ψ f 0 in B ∞ . With some exceptions we obtain the same results for the asymptotic behavior of the wave function in the harmonic gauge with f = 0 : Ψ = e −by 0 Ψ f 0 in the limit y 0 → −∞ , since the term (−by 0 ) in the exponent is suppressed generically by e −y 0 . The change of gauge (from tortoise to harmonic) may be sensitive for the asymptotic behavior of the Ψ -function in the case (3.51) if y = 0 and in the case (3.53) when y belongs to the border of the domain B ∞ .

Example 1: (n − 1) -dimensional billiards in the models with electric branes
Here we illustrate our approach by considering the model with the Lagrangian Here we deal with the metric g and the forms F s = dA s , s = 1, . . . , n , on the manifold M from (2.3). We use the metric ansatz from (2.4) which deals with a warped product of the interval (u − , u + ) and n Einstein spaces. ∆L is an extra term with fields of forms which will be specified below.

The configuration with n electric branes
Here we put ∆L = 0 and use the following non-composite electric ansatz for the fields of forms: For our configuration of branes (when d 1 > 1 , d 2 > 1 for n = 2 ) the billiard B ⊂ H n has a finite volume. Indeed, let us suppose that B has an infinite volume. Then there exists a set of (real) Kasner-like parameters α = (α 1 , ..., α n ) obeying the relations Let us consider the quasi-Cartan matrix [34] A ss ′ = 2(U s , U s ′ )/(U s ′ , U s ′ ), (4.6) where the scalar products (2.26) read in our case s, s ′ = 1, . . . , n . Thus we are led to the matrix  1, 1, 1) .
According to the classification of hyperbolic KM algebras by Carbone et al. [25] only the ranks n = 2, 3, 4 should be considered here, since for n > 4 there are no Dynkin diagrams where all nodes are connected by lines.
For n = 3 and (d 1 , d 2 , d 3 ) = (1, 1, 1) we obtain the hyperbolic KM algebra by the number 7 in the classification of Saçlioglu [23] (see also [24]), which is number 80 in the table of ref. [25]. In this case A ss ′ = −2 for all s = s ′ . This KM algebra appears for Bianchi-IX cosmology and its billiard coincides with the Chitré one. In the quantum case this billiard was considered in numerous papers; see [7,8,26,27,40] and references therein.
For the case n = 3 , d 1 = d 2 = 1 , d 3 = 2 we get the billiard corresponding to the hyperbolic KM algebra by number 40 in the classification of [25]. Here We have a billiard of finite volume which may describe the never ending oscillating behavior near either spacelike or timelike singularity. For the case of a spacelike singularity we put w = −1 and use all metrics g i of Euclidean signatures and all forms are taken to be ordinary ones. For the case with a timelike singularity we have three (non-equivalent) possibilities with w = 1 : (a) g 1 = −dx 1 ⊗ dx 1 , g 2 = dx 2 ⊗ dx 2 , and g 3 has the signature (+, +) (b) g 1 = dx 1 ⊗ dx 1 , g 2 = dx 2 ⊗ dx 2 , and g 3 is of signature (−, +) ; (c) g 1 = dx 1 ⊗ dx 1 , g 2 = dx 2 ⊗ dx 2 , and g 3 is of signature (−, −) . In the first two cases only one. form should be phantom: F 1 or F 3 in cases (a) or (b), respectively. In the case (c) all three forms are ordinary ones.
For the last example n = 4 , d 1 = d 2 = d 3 = d 4 = 1 we get the hyperbolic KM algebra by number 124 from [25] with A ss ′ = −1 for all s = s ′ . For our model with four ordinary 2-forms we get a diagonal cosmological metric with w = −1 and g i = dx i ⊗ dx i , i = 1, 2, 3, 4 , which describes a never ending oscillating behavior near the spacelike singularity. An analogous behavior will be obtained in approaching a timelike singularity, if w = 1 , g 1 = −dx 1 ⊗dx 1 , g i = dx i ⊗dx i , i = 2, 3, 4 , when the only one 2-form, namely F 1 , is phantom.
For the model under consideration the basis asymptotic solutions for the wave function are given by eq.

26) reads
Since our diagonalization (4.13) uses a timelike co-vector U Λ which coincides with sum of n brane vectors U = U 1 + · · · + U n , the whole of our analysis of the asymptotic behavior from the previous section is relevant.
We get |Ψ f 0 | → 0 as y 0 → −∞ for the basis solutions with modified Bessel function B = K . For the basis solutions with another choice of modified Bessel function B = I we obtain (generically) non-empty "spots" B ∞ ⊂ B (see (3.54)) for some basis functions, where the |Ψ f 0 | → +∞ . These "spots" appear in the model under consideration since the inequality  1, 1) . In this case we may have a point-like "spot" for y = 0 and a delta-function localization of the wave function Ψ f 0 for y 0 → −∞ . This singularity can be eliminated if we change to the harmonic gauge. In all other cases the radius of any "spot" is non-zero and the "spot" cannot be eliminated by a transition to the harmonic gauge. Equations (4.18), (4.21), when substituted into the general formulas of the previous section, will give a solution to the problem in the quantum case.
Here the asymptotic solution to WDW equation are found up to the spectrum of the (minus) Laplace-Beltrami operator (3.41) with the zero boundary conditions imposed. For the n = 2 case this can be done explicitly.
The case n = 2 . Now we consider the case n = 2 , when d 1 ≥ 2 , d 2 ≥ 2 . We use the following diagonalization of variables: (4.14) For 1-dimensional vectors from (3.31) we get . Thus we are led to 1dimensional billiard B = (y 1 , y 2 ) with point-like walls assigned to which obey −1 < y 1 < 0 and 0 < y 2 < 1 . B belongs to the 1d unit "disk" D 1 = (−1, 1) , which is an image of the 1-dimensional hyperbolic space H 1 ⊂ R 1,1 under the stereographic projection from the point (z 0 , z 1 ) = (−1, 0) . The billiard is subcompact, i.e. its completion [y 1 , y 2 ] is compact. We get (y 1 , 2) . In the quantum case the model with two factor spaces is integrable in the asymptotic regime of the formation of billiard walls. Here we have a discrete spectrum of the Laplace-Beltrami operator on B = (y 1 , y 2 ) with the metric h L = 4dy ⊗ dy/(1 − y 2 ) 2 ( y = y 1 ), when the zero boundary conditions at points y 1 , y 2 are imposed. Making the coordinate transformation we reduce the metric to the simple form h L = dx⊗dx and ∆ L = d 2 /dx 2 . We get a discrete spectrum of the Laplace-Beltrami operator on (x 1 = x(y 1 ), x 2 = x(y 2 )) with the zero boundary conditions: ∆ L Ψ L,k = −E k Ψ L,k , Ψ L,k (x i ) = 0 , i = 1, 2 : The calculations give where . (4.20) (For (d 1 , d 2 ) = (2, 2) we get x 2 − x 1 = ln(2 + √ 3) .) Hence the spectrum depends only on the parameter Q and the quantum number k : k = 1, 2, 3, . . . . In the symmetric case d 1 = d 2 = d the eigenfunctions Ψ L,k have a zero at x = 0 , or, equivalently, at y = 0 , only for even k . It may be shown that in the case d 1 = d 2 = d ≥ 4 (i.e. for D ≥ 9 ) the "spot" covers the billiard with the exception of zeros of Ψ L,k . This means that |Ψ f 0 (y)| → +∞ as y 0 → −∞ for B = I and any y , obeying Ψ L,k (y) = 0 . We recall that |Ψ f 0 (y)| → 0 as y 0 → −∞ for B = K and any y . This is valid for Q 1 = 0 , Q 2 = 0 .
In the classical case the model with two factor spaces is also integrable in the asymptotic regime of the formation of billiard walls. This is considered in the Appendix in detail.

The configuration with n + 1 electric branes
We extend the model from the previous subsection by adding an extra term which we put in the Lagrangian (4.1), where F 0 = dA 0 is D -dimensional form, i.e. n 0 = D . We supplement the electric ansatz (4.2) by the following relation We get an additional brane vector U 0 = U Λ and hence U = U 0 + U 1 + . . .+U n = 2U Λ . Since (U 0 , U 0 ) < 0 adding the term (4.22) to the Lagrangian (4.1) does not change the billiard and the asymptotic behavior of the scale factors (near the singularity). Now the basis asymptotic solutions for the wave function from the previous subsection are modified by adding the new variable Φ 0 and using another prefactor, and another A -number In this case the relations d s ≤ D − 3 (following from (U s , U s ) > 0 ) imply Since the A-number from (4.25) differs from (4.10) we are led to a different asymptotic behavior for the wave function Ψ f 0 → 0 as y 0 → −∞ in this case when the electric brane with the brane vector U 0 obeying (U 0 , U 0 ) < 0 is added. Now we consider the relation for "spots" (4.11). We get n s=1 (U s , U s ) + (U Λ , U Λ ) ≥ −4(U Λ , U Λ ), (4.26) or, equivalently, In this case we get a different restriction on "spots" in comparison with the relation (4.12). Indeed, due to (4.27) for D ≤ 7 the "spots" are absent and hence |Ψ f 0 | → 0 as y 0 → −∞ for all basis solutions. Thus the adding of the term (4.22) to the Lagrangian changes (e.g. drastically) the asymptotic behaviour of the wave functions as y 0 → −∞ , while in the classical case this term is irrelevant for the asymptotic behaviour.
In this case we get J (U J , U J ) < −(U, U ) and hence |Ψ f 0 | → 0 as y 0 → −∞ for all basis solutions. In the harmonic time gauge we also get |Ψ f =0 | → 0 as y 0 → −∞ for all basis solutions.
Thus in this example we have a similar asymptotic behavior of the basis wave functions to the case of 9d billiards with a maximal number of SMbranes, either electric [28] or electric plus magnetic [35].

Conclusions
We have generalized the quantum billiard approach from [35] by considering a cosmological-type model with n Einstein factor spaces in the theory with several forms, l scalar fields and a Λ -term. As in [35], after imposing the electromagnetic composite brane ansatz with certain restrictions for brane intersections and parameters of the model we have used the Wheeler-DeWitt (WDW) equation for the model, written in the conformally covariant form.
By imposing restrictions on the parameters of the model, e.g. on brane U s -vectors and using the vanishing of the potential terms coming from curvatures of the Einstein spaces and the Λ -term [15], we have obtained the asymptotic solutions to the WDW equation, in the limit of formation of billiard walls: y 0 → −∞ , which have a form similar to that from [35].
We have studied a subclass of examples of classical and quantum billiards in the model with n non-intersecting electric branes, e.g., for certain hyperbolic KM algebras of ranks n = 2, 3, 4 . In the classical case any of these billiards B has a finite volume and describes a never ending oscillatory behavior of the scale factors while approaching a singularity, which may be either spacelike or timelike.
In the quantum case the asymptotic basis solutions to the WDW equation in the "tortoise" time gauge tends to zero: Ψ f 0 → 0 , as y 0 → −∞ , for the following choice of the Bessel function, B = K , while for another choice of the Bessel function: B = I we have obtained for some basis solutions |Ψ f 0 | → +∞ when y 0 → −∞ in certain domain B ∞ ⊂ B of non-zero measure -"spot" -for all cases but (d 1 , d 2 ) = (2, 2) and (d 1 , d 2 , d 3 ) = (1, 1, 1) . The "spot" does not disappear in the generic cases when the wave function in the harmonic time gauge is considered. For two exceptional cases we have a point-like "spot" at y = 0 for some basis solutions in "tortoise" time gauge, which corresponds to δ -functions in the limit y 0 → −∞ , but this singularity does not take place in the harmonic time gauge. For n = 2 we have found the asymptotic solutions for the metric in the classical case (see the Appendix) as well as its quantum counterpart, i.e. the asymptotic (basis) solutions to the WDW equation.
Here we have considered the branes with general scalar products (U s , U s ) , while in [35] the restriction (U s , U s ) > 0 was used. The presence of branes with (U s , U s ) ≤ 0 has no effect on the billiard B and the asymptotic classical behavior of scale factors and scalar fields (as y 0 → −∞ ) but it changes the asymptotic solutions to WDW equation. We have illustrated this effect by an example with n electric branes on product of n Einstein spaces obeying (U s , U s ) > 0 and one brane obeying (U 0 , U 0 ) < 0 . It is shown that for such configuration of branes the "spots" disappear for D ≤ 7 and hence |Ψ f 0 | → 0 as y 0 → −∞ for all basis solutions to WDW equation. The analogous asymptotic behavior of all basis solutions Ψ f 0 is shown to be valid for the 4 -dimensional quantum billiard in the D = 11 model with ten magnetic branes, which was considered earlier for the classical case in [15]. This result can be extended to the configuration with composite magnetic S-branes in the model with one 4-form, but an open problem here is to include the Chern-Simons term of D = 11 supergravity into the consideration. In the classical case there were some obstacles for doing so [15].
Recently, a certain interest in studying a possible oscillating behavior near a timelike singularity, started by Parnovsky [5,6], appeared after refs. [41,42]. In [42] the authors speculated that such singularities, if occurring in AdS/CFT and being of the chaotic variety, may be interpreted as transient chaotic renormalization group flows which exhibit features reminiscent of chaotic duality cascades. So, the examples of billiards describing an oscillating behavior near a timelike singularity, which were considered in this paper, may be tested for a possible application to the program suggested in ref. [42]. Appendix. Classical assymtotical solution for n = 2 .
Using the continuity conditions for the scale factors at the points τ k we get for B(τ k − 0) = B (k) i , i = 1, 2 : (A.14) We get from (A. 13) ln(B for i = 1, 2 and k = 1, 2, . . . . Thus, we obtain an asymptotic (oscillating) solution for the metric with two scale factors. It does not depend on Λ , the curvatures of the Einstein spaces and the brane charge densities. This approximation works for small enough value of the parameter Cτ 1 = δ .