Table Of ContentFebruary7,2008 4:21 WSPC-ProceedingsTrimSize:9.75inx6.5in main
1
7 Classical and Quantum Aspects of the Inhomogeneous Mixmaster
0
Chaoticity
0
2
R.Benini12†andG.Montani23‡
n
a 1Dipartimento di Fisica - Universit`adi Bologna and INFN
J Sezione di Bologna, viaIrnerio 46, 40126 Bologna, Italy
2ICRA—International CenterforRelativistic Astrophysics c/o Dipartimento di Fisica (G9)
6
Universit`adi Roma “LaSapienza”, Piazza A.Moro 5 00185 Roma, Italy
1 3ENEA C.R. Frascati (U.T.S. Fusione), Via Enrico Fermi 45, 00044 Frascati, Roma, Italy
†[email protected]
1
‡[email protected]
v
4
9 WerefineMisner’sanalysisoftheclassicalandquantumMixmasterinthefullyinhomo-
0 geneouspicture;webothconnectthequantumbehaviortotheensemblerepresentation,
1 bothdescribethepreciseeffectoftheboundaryconditionsonthestructureofthequan-
0 tumstates.
7
Nearthecosmologicalsingularity,thedynamicsofagenericinhomogeneouscos-
0
mological model is reduced by an ADM procedure to the evolution of the physical
/
c degrees of freedom, i.e. the anisotropies of the Universe. In fact asymptotically to
q
- the Big-Bangthe space points dynamically decouple1 because the spatialgradients
gr of the dynamical variables become of higher order2 and we can model the inhomo-
: geneous Mixmaster via the reduced action:
v
i
X
I = d3xdτ(p ∂ u+p ∂ v ǫ), ǫ=v p2 +p2, (1)
ar u τ v τ − u v
Z
p
The dynamics of such a model is equivalent to (the one of) a billiard-ball on a
Lobachevsky plane; this can be shown by the use of the Jacobi metric.1 The man-
ifold described turns out to have a constant negative curvature, where the Ricci
scalaris givenby R= 2/E2:the complex dynamicsofthe genericinhomogeneous
−
model results in a collection of decoupled dynamical systems, one for each point of
the space,andallofthemequivalentto abilliardproblemonaLobachevskyplane.
We wantto investigatethe relationexisting betweenthe classicalandthe semiclas-
sical dynamics, and from this analysis we will derive the correct operator ordering
to be used when quantizing the system.
Let’s write down the Hamilton-Jacobi equation for the system
2 2
2 2 δ 0 δ 0
ǫ =v S + S (2)
δu δv
(cid:18) (cid:19) (cid:18) (cid:19) !
February7,2008 4:21 WSPC-ProceedingsTrimSize:9.75inx6.5in main
2
This can be explicitly solved by separation of constants
ǫ+ ǫ2 k2(ya)v2
S0(u,v)=k(ya)u+ ǫ2−k2(ya)v2−ǫln 2 p ǫ−2v !+c(ya). (3)
p
The semiclassical analysis can be developed furthermore, and the stationary con-
tinuity equation for the distribution function can be worked in order to obtain
informations about the statistical properties of the model; as soon as we restrict
the dynamics to the configuration space, we get the following equation for the dis-
tribution function w˜
∂w˜(u,v;k) E 2 ∂w˜(u,v;k) E2 2k2v2 w˜(u,v;k)
+ 1 + − =0 (4)
∂u s kv − ∂v kv2 E2 (kv)2
(cid:18) (cid:19) −
This can be solved, and the exact distribution function capn be obtained as soon as
we eliminate by integration the constant k
Ev g u+v kE2v22 −1
w˜(u,v)= (cid:18) q (cid:19)dk (5)
−E v√E2 k2v2
Z v −
Itisworthnothinghowinthecaseg =const,themicrocanonicalLiouvillemeasure
w (u,v)= π is recovered.
mc v2
Weexpectthatthedistributionfunctionw˜(u,v)isre-obtainedassoonasthequan-
tum dynamics is investigated to the first order in ~. This can be easily done as a
WKBapproximationtothe quantumdynamicsisconstructed;assoonasweretain
only the lower order in ~, we obtain that: i) the phase S(u,v) coincides with the
Hamilton-Jacobifunction,andii)the probabilitydensityfunction r(u,v) obeysthe
following equation:
∂r(u,v) E 2 ∂r(u,v) a(E2 k2v2) E2
k + k2 + − − r(u,v)=0. (6)
∂u s v − ∂v v2√E2 k2v2
(cid:18) (cid:19) −
that coincides with (4) for a particular choice for the operator-orderingonly, i.e.
∂ ∂
vˆ2pˆ2 ~2 v2 . (7)
v → − ∂v ∂v
(cid:18) (cid:19)
With this result, the problem of the full quantization of the system can be taken
in consideration. The main problem is the presence of the root square function
in the definition of the Hamiltonian (1), but well grounded motivations exist3 to
assume that the real Hamiltonian and the squared one have same eigenfunctions
and squared eigenvalues. This way the solution of the eigenvalue equation can be
obtained
Ψ(u,v)= anKs−1/2(2nπv)sin(2nπu) (8)
n>0
X
The spectrumisobtainedasDiricheletboundaryconditionsonthe domainΓ (see
Q
Fig.1) are taken into account. The condition on the vertical lines can be imposed
February7,2008 4:21 WSPC-ProceedingsTrimSize:9.75inx6.5in main
3
exactly, but the one on the semicircle cannot be solved exactly; so we approximate
it as in Fig.2 with a straight line v =1/π (it preserves the domain area µ=π).
All these imply that s = 1/2+t, with t , and that the spectrum assumes the
∈ ℜ
3
2.5
2
v 1.5
1
0.5
0-2 -1.5 -1 -0.5 0 0.5 1
u
Fig.1. TherealdomainΓQ Fig.2. Theapproximatedomain
following form
1
E2 =( +t2)~2 (9)
t 4
The values of the real parameter t have to be numerically evaluated by solving the
equation K (2n) = 0 for every natural n. This condition implies a discrete but
it
quite complicated shape for the spectrum.
Asymptotic expansions for high occupation numbers can be derived for different
regions of the parameters (t, n)4
References
1. R.Benini, G. Montani, Phys. Rev. D 70, 103527 (2004)
2. A.A. Kirillov, Zh. Eksp. Theor. Fiz. 103, 721-729 (1993)
3. R.Puzio, Class. Quantum Grav. 11, 609-620 (1994)
4. R.Benini, G. Montani, Class. Quantum Grav. in pubblication