Square-torsion gravity, dark matter halos and the baryonic Tully-Fisher relation

Square-torsion gravity is applied to the long standing dark matter problem. In this context the theory reduces to General Relativity complemented by a dark stress-energy tensor due to the torsion of spacetime and is studied under the simplifying assumption of spherical symmetry. The dark stress-energy tensor is found to satisfy an anisotropic structure equation. In vacuum this is shown to be equivalent to a wave equation with sources. A natural class of exact solutions is found which explicitly perturbs any seed spacetime metric by a conformal factor satisfying a (1+1)-dimensional wave equation. This leads to the concept of dark coating. The static solutions are then used to construct structures that model dark matter halos surrounding baryonic bodies. In the Newtonian r\'egime the baryonic mass $m_b$ and the flat rotation curve velocity $v_f$ are found to be related by the baryonic Tully-Fisher relation $m_b\propto v_f^4$. The present work proposes thus a possible theoretical motivation of this hitherto purely empirical result. The example of a dark halo on the Schwarzschild geometry is made as a toy model for a galaxy. All qualitative an quantitative features of galactic rotation curves are recovered. A dark halo surrounding a Schwarzschild black hole is found to possess a boundary of staticity called torsion sphere placed between the photon sphere and the event horizon. The phenomenon of dark radiation is briefly exposed. The way for cosmological applications is then opened by showing how Hubble expansion is a natural feature of the theory.


INTRODUCTION
The present work investigates an alternative theory of gravity closely related to General Relativity. Both theories have a common mathematical and aesthetic foundation: gravitational phenomena are described by the geometry of spacetime. The main difference between these sibling theories is that while Einstein's theory derives from an action functional containing only the Riemann curvature tensor, Square-torsion Gravity completes such functional by a term quadratic in the torsion tensor. Our main motivation for the investigation of alternative theories of gravitation is the hope of incorporating into the same framework phenomena that have a long standing problematic interpretations in terms of General Rel-ativity. This article focusses on the so called dark matter problem.

The dark matter problem
The phenomenon of dark mater has been know for almost a century now, its scientific history is thus a subtle and fascinating matter in itself. We hereby give a very concise historical perspective and invite the reader to refer to the very complete historical review [1] and its rich references.

Early history
The first signs of a new open problem arose as early as 1933, when Zwicky estimated the mass of the Coma galaxy cluster using Virial's theorem and the relation between redshift and distance. He obtained a very high mass-to-light ratio. This suggested the presence of dark matter in some form. This fact could as well be explained by the possible non equilibrium of the cluster. For this reason during the 1950s much effort was focussed on determining mass-to-light ratios of various clusters. It became however more and more difficult to reject the presence of dark matter using the hypothesis of non equilibrium. The main perplexity of the researchers came from the great age of the galaxies forming the clusters. The idea thus formed that 90 to 99% of the mass of these clusters could be in the form of extragalactic material. Indeed by the 1970s the relaxation process of clusters seemed to support this idea.
Further enquiries of x-ray emissions put however very stringent upper bounds to intergalactic hot gas, orders of magnitude less than the gas required to explain gravitational observations. This opened the way to more exotic possibilities.
Then another major observational difficulty arose, namely that of galactic rotation curves.
The velocity of the stars orbiting around the galactic centre as a function of the distance from the centre itself can in principle be used to infer the mass distribution in the galaxy. The problem was that, assuming the mass distribution closely follows the luminosity distribution, one expects rotation curves that decrease in a Keplerian way after some critical value of the distance from the galactic centre. What one observes is typically very different, namely flat rotation curves. By the mid 1970s enough evidence were collected from radio astronomy to convince many that the outskirts of galaxies should contain very large amounts of invisible mass.
Naturally, many hypotheses began to be formulated regarding the exact nature of this dark matter. Some suggested it could comprise planets, brown dwarfs, red dwarfs, neutron stars, and black holes. Another paradigm was that of weakly interacting fundamental particles yet to be identified. However, gravitational micro-lensing observations of galactic halos seemed to rule out the first possibility.
The term dark matter has nowadays become synonymous of particles that may account for the missing matter in our Universe.
Taking this point of view, one has then two possibilities to model such particles, as a hot or as a cold gas. These two models lead to very different structure formation on large scales.
The most widely accepted model is that of cold dark matter, ruling therefore out candidate particles of the standard model such as relativistic neutrinos. This idea gained consensus as early as the end of the 1980s.
The picture handed to the 21 st century is hence that of a pervasive dark dust, of yet unknown nature. In spite of all the refined techniques which have been devised during the last decades and in spite of the dedicated search undertaken by many prized experimentalists of these new particles, no successful candidate was found.

Recent developments
The last two decades shed a new light onto two additional elements worthy of meditation: • The baryonic mass m b of a galaxy is strongly correlated to the magnitude of the flat rotation velocity v f . This relation, known as baryonic Tully-Fisher relation, takes the f . As reported in [2] there does not seem to be any dependence on other properties like galactic size or surface brightness. This points out a possible deeper connection between the baryons and the physics which determines v f .
• Data from hundreds of very diverse spiral galaxies seem to point out a fundamental correlation between the baryonic matter distribution and that of dark matter, irrespective of the fact that the galaxy be baryon or dark matter dominated [3].
A theoretically oriented mind could hence be tempted to embrace the possibility that the phenomena hitherto analysed using the dark matter paradigm might be explained by funda-mental laws of gravitation yet to be understood. We hereby explore this very possibility and attempt a purely geometrical description rather than one derived from fundamental particles. It is however important to note that there are fundamental observational differences between galaxies and clusters of galaxies. To name one, the baryonic Tully-Fisher relation can only be observed in galaxies and not in galactic clusters. In spite of this fact, recent observational research indicates that there might be a nearly fundamental constant which parametrises the dark-baryonic matter interplay. We invite the interested reader to refer to [4] and its rich references.

Square-torsion gravity?
All this being said, a specific theory has to be chosen. A promising candidate is the so called square-torsion theory of gravity. It was recently shown in [5] that in this theory some components of the torsion tensor remain to be fixed a priori as though they were external sources. This analogy with matter becomes mathematically exact for a spacetime with vanishing spin density. The free degrees of freedom then appear in the Einstein equations as a stress-energy tensor due to torsion, hereafter christened "dark stress-energy tensor".
Given that such stress-energy can take very diverse forms, e.g. that of dust, this theory is a natural choice when trying to geometrise dark matter.

SQUARE-TORSION GRAVITY IN BRIEF
We hereby give a very brief overview of the main features of Square-torsion Gravity, we refer again to [5] and references therein for a more detailed account.

Main ingredients and notation
Our mathematical model for spacetime is that of an orientable four-dimensional differentiable manifold M endowed with a Lorentzian metric g. We use a local orthonormal basis e := {e I } I∈{0,1,2,3} of the tangent space T M. Let us call θ a local basis of the cotangent space T * M defined by θ I (e J ) = δ I J , or in abstract index notation θ I µ e J µ = δ I J , where Greek indices run from 0 to 3 and denote ordinary indices of T M and T * M respectively. With this structure, one point-wise recovers all the machinery of special relativity. In particular, each affine space is endowed with a Minkowski metric g IJ = η IJ = diag(−1; +1; +1; +1) constant throughout the computations. The inverse frame field θ I µ pulls back this metric metric g IJ to define a metric on spacetime. The inertial structure of spacetime is encoded in an affine connection, equivalent to the definition of a covariant derivative of our reference frame as ∇ µ e I := e J ω J Iµ . (1) The torsion 2-form is defined as while the Riemann curvature 2-form is defined as where we use the common bracket notation to denote a skew-symmetrisation of indices.
Relation between the Levi-Civita and a general torsionful connection: It is well known that there is a unique metric-compatible connection which is torsion-free. This is called Levi-Civita connection and we denote it by Γ I Jµ . From definition (2) of the torsion 2-form one can easily demonstrate that Notice how we lowered indices for practical purposes using the metric g IJ and that the brackets are Lie brackets. In the same way for a generic torsionful connection ω I Jµ we can relate it to Γ I Jµ as follows where is called the contorsion tensor.

Action functional and equations of motion
The action we consider is where L s is the Lagrangian density for all external sources coupled to geometry and d 4 xdet(θ) the volume form. This is the Einstein-Hilbert action complemented by a quadratic dynamical term for the torsion of spacetime. It is conceived as a first order formalism action principle, in which the connection ω IJ µ and the tetrad θ I µ are independent variables. The equations of motion derived from this action by its variation with respect to the connection yield equations that have intrinsic angular momentum as source. In most macroscopic applications, when dealing with perfect fluids or more general fluids that do not include covariant derivatives in their structure equation, we can assume said spin density to vanish. This in turn means that T S IS = 0 The only non vanishing irreducible component of the torsion tensor is thus that which is traceless and has vanishing antisymmetric part. This irreducible component shall be denoted by T .
Let us now come to the equations of motion having the stress-energy tensor of matter P IJ as source. As said, we are interested in non-nematic phases of macroscopic matter and cosmological fluids. In this setting the equations of motion derived from the action functional by its variation with respect to the tetrad are simply where G IJ is the usual Einstein tensor. This is very interesting because we still have 16 free torsion components (up to gauge choice) and, if one interprets these equations as the Einstein equations, the tensor plays the role of a dark content. Indeed, taking the covariant Levi-Civita divergence of (9) we find If we reason in terms of the Levi-Civita connection we are hence forced to consider D IJ as a stress-energy tensor due to the torsional degrees of freedom of our geometry because it is included in the conservation of energy and momentum. Notice that isolated point particles follow geodesics, because semi-classical matter fields are governed by their own equations of motion, which in the eikonal approximation imply geodesic motion.
It is useful to work with a local decomposition of the torsion tensor in terms of its "electric" and "magnetic" parts. To this end we choose a preferred time direction θ 0 µ and denote spacelike indices with lowercase Latin indices running from 1 to 3. We can thus uniquely define In turn, it is useful to decompose the 3-dimensional matrices E ij and B ij into their traceless symmetric and skew-symmetric parts: Using the equations of motion one can finally write: (Equation (15) corrects a misprint in [5].)

THE ASSUMPTION OF SPHERICAL SYMMETRY
From what we have seen in the previous section, the structure of D remains possibly quite complex. For this reason we specialise our discussion to configurations with spherical symmetry, which, in addition to being a strongly simplifying assumption, is of great physical interest.
What do we mean by spherical symmetry in this context? We adopt the common definition according to which a Lorentzian manifold (M, g) in spherically symmetric if it has the group of three dimensional rotations SO(3) as isometry group. We moreover require that SO(3) should act in such a way that its orbits be two dimensional spacelike surfaces.
It can be demonstrated, see e.g. [6], that the requirement of spherical symmetry formulated in this fashion implies that (M, g) be locally a warped product where (M,g) is a two dimensional Lorentzian manifold and (S 2 , g S 2 ) denotes the two-sphere.
R :M → R is the so called warping function, which can be though as "radius function" given that the warped product metric by definition takes the form Let us now turn to the torsion tensor. Spherical symmetry implies the impossibility to single out a preferred element belonging to the distribution of vectors tangent to S 2 . If e.g. we introduce spherical angular coordinates so that g S 2 = d ϑ ⊗ d ϑ + sin 2 ϑ d ϕ ⊗ d ϕ and use an adapted orthonormal frame in which e 2 := R −1 ∂ ϑ and e 3 : (recall that e ij and b ij are traceless). The dark stress energy thus becomes Notice how both the energy density component and pressure can have a priori any sign. Let us then define ρ := D 00 , P := D 11 and P ⊥ := D 22 . The dark stress-energy tensor then satisfies the structure equation Remark on staticity In view of the following applications it is important to dwell a moment on the expression we obtained for the energy current density component J := D 01 . At first sight it might seem completely innocuous. It has however important implications e.g. for the study of vacuum configurations. In vacuum the stress-energy tensor due to matter vanishes. It then becomes the only source for the Einstein tensor. This implies that the static character of vacuum solutions is governed by J. We hence state the following Compatibility condition on J • if a dark stress-energy tensor has J = 0 then ρ and P have opposite sign (or vanish) or equivalently • if a dark stress-energy tensor has ρ and P with the same sign (and non vanishing) The two statements being logically equivalent, let us demonstrate the second point. We start by assuming that ρ > 0 and P > 0. Expression (19) then implies |β| > |ε| and |b| > |e|.
We can thus write the little chain of inequalities |eβ| = |e||β| > |e||ε| > |b||ε| = |bε|. From this we deduce that J = 0. The case where both ρ and P are negative can be demonstrated in the very same way by swapping the role of the electric and magnetic components.
In Square-torsion gravity one can hence have spherically symmetric vacuum solutions that are not static. This fact should be clear if one considers that the present theory is equivalent to General Relativity in the presence of a fluid. Birkhoff's theorem does hence not apply in this setting.

Remark on gauge freedom
Due to the fact that D is quadratic in the torsion tensor, the relation between the components of torsion and those of the dark stress-energy tensor is not invertible. The choice of torsion tensor remains thus not unique. This freedom can be parametrised e.g. by two hyperbolic angles (η, τ ) as follows: • for ρ > 0 and P > 0 • the case of different signs is obvious.
These hyperbolic angles are however not independent. We still have to solve for the energy current density J. It is easy to see e.g. that if J = 0 then τ = ±η.

THE EQUATIONS OF MOTION IN VACUUM
In order to compute the equation of motion for our spacetime geometry in vacuum we now proceed with the computation of the Einstein tensor for our warped product. To this end, let us introduce a pair of coordinates (t, x) onM and an adapted local orthonormal frameθ := (θĨ)Ĩ =0,1 , so that a tilde singles out directions tangent toM. Let then a hat single out directions tangent to S 2 in an analogous wayθ := (θÎ)Î =2,3 . Let us also change variables from the warping function R to its logarithm ψ := log(R). It is a standard task to compute the Riemann and Ricci tensors of the Levi-Civita connection associated to such warped product, we refer again to [6] for more details. LetRĨJ be the components of the Ricci tensor ofg, with this notation one finds the following expression for the components R I J of the Ricci tensor of g: where we used the fact that the Ricci tensor of S 2 takes the simple form 1 2 δÎĴ . Let us now translate the structure equation (20) to the components of the Ricci tensor.
In terms of the Einstein tensor we clearly have and thus This is just the right balance to eliminate the term quadratic in ψ ,S . We find The last expression can be formulated in an even tamer way by remembering that the Lorentzian metricg is two-dimensional and is thus locally equivalent to a conformally flat . The curvature scalar RTT associated to such metric can thus be expressed as We can thus rephrase (25) as and changing variables to Σ := 2φ + ψ and ∆ := φ − ψ we finally come to the result Notice that this equation, together with the compatibility condition on J, is the only we need to solve. This is due to the fact that, apart from the structure equation, ρ and P can take any desired form descending from a solution of (28). This is exactly what we expected, in spherically symmetric Square-torsion Gravity one has one free parameter in the dark stress-energy tensor which can be used to model the observed dark matter density.
An important question then arises. Is there any such parametrisation of dark matter which can be theoretically preferred?

CONFORMAL COATINGS
Our evolution equation (28) is a (1+1)-dimensional wave equation with source. The wave function Σ is independent from the source function ∆. This means that, given any seed solution (Σ 0 , ∆ 0 ) of (28) the couple (Σ 0 + 3 ς, ∆ 0 ) is still a solution if which, rewritten in terms of our new variables (Σ, ∆) reads It is now clear that the above mentioned symmetry in the space of solutions has a direct geometrical interpretations as harmonic conformal rescalings of solutions. This leads naturally to the concept of conformal coating.

Conformal coatings in vacuum
If we have a seed spherically symmetric solution g 0 satisfying the vacuum equation (28) any metric g = e 2 ς g 0 is also a solution to (28), provided that the compatibility condition on J be verified.
In (32) we have emphasised the fact that ς satisfies the wave equation for the twodimensional Lorentzian metricg 0 . In particular, it is immediate to write local expressions of solutions to (32) by means of D'Alembert's formula which shows that sufficiently regular and bounded initial perturbations ς(0, x) remain such in any chart. This class of exact metric perturbations is hence of great physical interest.
This is arguably the most natural class of solutions to (28) because it descends from a natural symmetry of said equation. If we have a physically preferred solution g, equation (32) provides its "dark coating", so to speak. Moreover, we are about to see that this coating method extends directly to non empty Einsteinian solutions.
Example regarding more general vacuum scenarios: Let us make an explicit example in order to show how dangerous more general perturbations of (28) can be. Let us choose g 0 very simply to be the Minkowsky metric. If we perturb also ∆ 0 with a perturbing function δ our equation becomes Let us now assume 2 d := e 2δ − 1 to depend linearly on ς so as to have for some α ∈ R. L α is evidently a Schrödinger operator for the potential α/x 2 , generally known as Calogero potential. It suffices then to study the time independent Schrödinger equation associated to L α . Loosely speaking, from (35) one deduces that a negative part of the spectrum of L α will result in unstable solutions. This problem is studied in great detail in [7]. In particular, for α < −1/4 such negative spectrum does not only exist, but is unbounded from below! Notice that we made an exact consideration, no linearisation was needed, we just chose a simple ad hoc perturbing direction in an appropriately well behaved subspace of solutions to (28).

General conformal coatings
For ease of discussion and future use, let us preliminarily express the Ricci and Einstein tensors for a conformally coated metric in function of their uncoated relatives. To this end it suffices to express a coated orthonormal frame of our choice in therms of its uncoated relative as ϑÏ = e ς θ I , thus denoting coated indices with a diaeresis. Using the vanishing torsion condition one can write from which we deduce that and this enables us to express the Riemann curvature two-form as This expression for the Ricci tensor is of great relevance because it shows how an extended version of our result on conformal coatings holds: If we have a spherically symmetric seed square-torsion gravitational solution g 0 , possibly describing an Einsteinian spacetime filled with matter, then g = e 2ς g 0 is also a square-torsion gravitational solution, provided that the compatibility condition on J be verified.
Let us demonstrate this fact by checking whether the structure equation holds for ς R I J . Notice first of all that because of spherical symmetry ς :M → R and one clearly has ς ,Î = 0. This in turn implies that g 0 ς = 0 is equivalent to ς ;SS = 0. We hence have ς ; T T = ς ;ŜŜ , while ς ,T ς ,T = ς ,S ς ,S . We can thus write the following chain of equalities: From this expression of the Ricci tensor it is now easy to derive that of the Einstein tensor: more so thanks to the fact that the contracted Bianchi identity holds separately for the two components. Let us express this fact more eloquently. Let P IJ be the total stress-energy tensor. As done before, let us write it as a sum of a baryonic contribution P b IJ and a dark one D IJ . We then write our field equations as and write the two separate equations This last step is meaningful because of the restricted identity To prove this, let us first express the coated covariant divergence of a generic symmetric tensor PÏJ in terms of uncoated objects. By means of (37) it is easy to see that If we now insert the definition (42) of Notice now that we are using the uncoated Levi-Civita connection, which we explicitly expressed in (4) as Lie brackets of frame fields. We use this explicitly and write and (45) is thus proved.
Thanks to the coated contracted Bianchi identity we can finally express the coated baryonic conservation equation It is hence physically meaningful to separate the two Einstein tensor components as above.
Notice how this is a coated expression, it hence depends on the dark stress-energy tensor.
This is what in the present model gives rise to dark matter-related effects such as the departure from Keplerianity of rotation curves.

APPLICATIONS
Now that we know how to construct spherically symmetric dark coatings on any base solution of our choice, we turn to physically relevant examples. The choice of a base is quite solid, it is sufficient to take any relevant Einsteinian setting and use it as a plinth.

Static setting
We shall now focus on static solutions. We will hence use a chart ofM with coordinates (t, r) in which we can express staticity by the existence of a timelike killing filed K = ∂ t .

Spherically symmetric cluster
Due to its conceptual and practical importance, we now turn to a Spherically symmetric cluster surrounded by vacuum. Such an idealisation could conceivably be used as a toy model for both spherically symmetric galactic clusters and galaxies. We then imagine to have a continuum model of dust. Let ρ 0 be the rest frame dust mass density and u I the four-velocity field components of our dust continuum. The general form of the stress-energy tensor associated to such fluid is well known to be Let us furthermore assume that this mesoscopic model arises from a microscopic model of point-like particles bound in approximately circular orbits. This fact can be minimally implemented by considering u I to be a random variable on each p ∈ M with a spherically symmetric probability measure giving rise to the following moments where v is the angular velocity of the dust particles with respect to our reference frame.
This expression simply states that at each point one is equally likely to find a dust particle swarming in any angular direction, that the second moment of the radial component of the velocity is negligible with respect to that of the angular component (approximately circular orbits) and that the mesoscopic magnitude v of the orbiting velocity and the rest frame mass density ρ 0 are no random variables, this to stick to the essential features only. Notice that the macroscopic stress energy tensor < P I J > retains its normalisation P S S = −ρ 0 as it should. If we introduce the same notation used for the dark stress-energy tensor, the macroscopic stress-energy tensor components due to baryons are then Let us then introduce Schwarzschild-like coordinates (t, r) and an uncoated orthonormal frame such that θ 0 = e Φ d t and θ 1 = e Λ d r. The expression of the connection and curvature for a fluid-filled spherically symmetric spacetime can be found in any General Relativity textbook, we refer e.g. to [8] for a very complete account. Applied to our case, one finds in Uncoated case: Let us preliminarily discuss the uncoated case. The only missing element in our discussion is that of solving the Einstein equations for P b ⊥ . To this end, let us write the uncoated conservation of stress-energy and focus on the only non trivial component More explicitly one has which leads to the anisotropic Tolman-Oppenheimer-Volkoff equation By specializing (57) to our case, one finally finds We then recover not only Keplerian rotation curves in the Newtonian limit (v, Gm/r << 1), but the full geodesic motion of a test particle in the geometry generated by the surrounding dust swarm.
Coated case: In order to study the coated case we first proceed with the solution of ς ;SS = 0 in the static case, namely for some scale factor R > 0.
Komar mass: The question now arises as to which sign to choose in expression (60). To this end, let us recall that we have a timelike Killing vector filed ∂ t . Thanks to the Ricci identity for Killing fields and introducing the 1-form kÏ canonically associated to our Killing vector field we can write see e.g. [9] for reference. This enables one to define the notion of Komar mass, which for stationary spacetimes is a conserved quantity associated to the time translations generated by ∂ t . In our case one can thus introduce the Komar mass density expressed as Thanks to the clear cut separation of the dark component of the Ricci tensor from the rest one can thus introduce the notion of dark Komar mass density expressed as This can be rephrased more explicitly thanks to (39) and noticing that k0 = −e ς+Φ 4πG e −ς ς ρ k = ± Γ 01 0 + 2Γ 21 If we require this quantity to be positive over the whole chart without any condition on the function m(r) we must then choose the positive solution for ς ,1 .
We can finally apply the conservation law of baryonic stress-energy (50) using its explicit form (46) and write the coated anisotropic Tolman-Oppenheimer-Volkoff equation which, applied to our fluid model, yields Let us now focus on the Newtonian régime. We define this as the limit in which the quantities Gm/r, Φ, r/ R << 1 and all their combined higher integer powers are negligible with respect to the first. In this case, (66) can be expanded at first order as In this régime we thus have a correction to the Keplerian behaviour of rotation curves.
Interpreted from a Newtonian point of view, this is equivalent to a Newtonian dark matter density ρ D N taking the form and corresponds to what one would get from the expression of the Komar dark matter density (64) in this limit. It is interesting to note that this essentially agrees with the central density profile found in computer simulations of cosmological halo formation in the cold dark matter paradigm, this is known as the Navarro-Frenk-White profile see e.g. [10] and references therein.
and inserting this into (67) we finally end up with We hence found that if a rotation curve is flat at the outskirts of a spherically symmetric astrophysical object coated by dark matter in the Newtonian régime, then there is a universal relation between the total baryonic mass and the flat rotational velocity. This is a relation which is known to hold in a surprisingly wide range of astrophysical situations, we refer again to [2] for a full account of the observational peculiarities of this fact and references to further literature on the subject. In its wide survey of literature on the subject, [4] reports the value for galaxies This result was originally stated in terms of solar masses M . Here we used 2GM = 2.9 10 3 m. To put things into perspective, this roughly corresponds to twice the diameter of the observable universe.
Notice that in our toy model we are only able to predict a universal proportionality law, the proportionality constant, being an integration constant, has still to be determined empirically and depends in principle on the particular celestial object of application. The fact that a seemingly universal constant is experimentally found might indicate (within the present scheme of thought) that the dark coating paradigm might have a range of application beyond spherical symmetry, thus predicting the same coating constant R for the region of the universe we live in. The scale of R could be seen as a further tantalising evidence of this.
As a didactic idealised scenario, let us now determine the baryonic mass density which in the Newtonian régime corresponds to a completely flat rotation curve. It suffices to consider the flatness condition as a first order differential equation for the function m(r), namely which is easily solved by for which (70) holds for the total mass enclosed in the ball of positive mass density. In this example we see how the centre of our spherical object is dominated by baryonic matter, while the outskirts are dominated by dark matter.

Schwarzschildean limit
If in our toy model we take the ρ 0 −→ 0 limit for the swarm of dust, while keeping the total enclosed mass m(r) = m b constant we clearly obtain the dark coating of the exterior Schwarzschild spacetime. This is to say, the Schwarzschild metric gets conformally rescaled by a conformal factor e 2ς , where which is nothing but the so called tortoise coordinate. It is then clear from (42) that the conformal factor e 2ς in general generates a curvature singularity at the (now shattered) Schwarzschild horizon. This problem is only apparent though. We must not forget that the compatibility condition on J still has to be satisfied. While this was no issue in the Newtonian régime, it is one if we insist in approaching the Schwarzschild radius. To examine this more closely, let us explicitly write ρ D , J D and P D in this context. From (42) we deduce and using we obtain The scalar curvature gains hence a singularity at r = 2Gm b . However the compatibility condition on J D demands that ρ D < 0 given that P D > 0 on the whole r > 2Gm b region. Let us then rephrase the ρ D < 0 condition in terms of the dimensionless quantities µ := Gm b / R and x = r/ R . One finds This means that our solution is not acceptable within a spherical region of coordinate radius For a typical ultra compact object of astrophysical relevance we have µ << 1, so that for all practical purposes one has Let us now turn to the rotation curves obtained in this limit setting. It suffices to specialise (66) to the case where m(r) = m b . For later convenience, let us use the dimensionless quantities µ and x as above. One has It is then instructive to visualise this for a point-like analogue of the Milky Way, for which the greatest part of the baryonic matter is confined in a very compact region around the origin.
In what follows we are only concerned by the qualitative and quantitative correspondence of out toy model with reality only as far as orders of magnitude are concerned due to the idealised nature of our description. We refer to [11] for more details regarding astrophysical data. For our purposes it is sufficient to note that the ordinary baryonic content of our galaxy is estimated to be around 10 11 M (compared to 10 12 M obtained including dark matter-related phenomena), which corresponds to µ = 6 · 10 −14 . The rim of the galactic It is already evident that x r is of order √ µ and that we thus are in the flat rotation curve region. In FIG. 2. we drew our rotation curve against its Keplerian relative (dashed line), highlighting the position that our sun would have. We thus find the orbital velocity in the outer region to be 8 · 10 −4 , in S.I. units this is equivalent to 2 · 10 5 ms −1 , exactly of the same order as the typical stellar velocity in the disc region of our galaxy. The picture turns then out to be consistent, with no need for exotic forms of matter.

Dynamical setting
Let us now turn to a dynamical setting, i.e. devoid of timelike Killing fields. A general solution to (40) can be formally written in terms of local lightcone coordinates (U, V ) for the two dimensional Lorentzian metricg 0 as follows:

Dark radiation
In general, the f -component can be seen as an ingoing spherically symmetric conformal wave, while the g-component as an outgoing spherically symmetric conformal wave. If one constructs a dark coating of the solution corresponding to a spherically symmetric pulsating star surrounded by vacuum one expects the boundary with the vacuum region to provide boundary conditions compatible with outgoing waves. However, the problem of the details of e.g. stellar pulsations and that of the appropriate boundary conditions go beyond the scope of the present work. A useful exercise in that direction is to compute from (42) the stress-energy tensor P g of a g wave coating of the Schwarzschild exterior solution. From an uncoated perspective it is easy to see that such wave carries the following stress-energy tensor components for which it is easy to see that the compatibility condition on J always holds true. Notice how this stress-energy tensor has a traceless far field limit which can be interpreted as the stress-energy tensor of a swarm of lightlike particles. This analogy becomes more evident if we consider a periodic g with mean values over a period T taking the values < g > T = < g > T = 0 and < g 2 > T = γ (e.g. a sine wave). Let us call k I the normalised future directed null vector pointing in the radial direction (k 0 = k 1 = 1). Thanks to the time translation invariance of the uncoated base, it is meaningful to compute the period-averaged stress-energy tensor < P g I J > at every event. We find This has thus the exact form of the stress energy tensor of a swarm of lightlike particles, whose associated wavenumber gets redshifted by the correct factor 1 − 2Gm b r . This form of dark radiation would be a promising observational signal and test of the present theory given its stark contrast with the Einsteinian setting, where no spherically symmetric gravitational radiation is present.

Hubble expansion
As a final remark, it is interesting to note that a linear Hubble expansion with Hubble constant H can always be implemented by a dark coating by choosing linear functions g and f so that f = −g = H. If we take again the Schwarzschild spacetime and take an uncoated point of view one finds ρ H = 1 8πG This satisfies the compatibility condition on J. Notice that it is an isotropic stress-energy tensor, with equation of state P = −ρ/3 and negative pressure. This is just the critical value of the perfect fluid constant w = −1/3 which underlies many cosmological discussions and makes room for constant expansion. Such possibilities are hence expected to be very promising on the cosmological arena and provide natural toy models of point-like perturbation of the FLRW cosmological paradigm [12].

CONCLUSIONS AND OPEN QUESTIONS
The present theory seems to offer a promising framework for the interpretation of dark matter related phenomena in terms of geometry. In particular, the concept of conformal dark coating, pointed out by the symmetries of the equations of motion, seems of wide applicability, from galactic dynamics to cosmological scales. Further work has to be done in order to determine the robustness of this paradigm outside the spherically symmetric setting. Axially symmetric models will hopefully shed more light on galactic dynamics, thus enabling to produce actual fits of real rotation curves. The question of the non uniqueness of solutions, especially in the non static setting, remains wide open though. The question of possible realistic physical sources of dark radiation has also to be addressed through a thorough implementation of e.g. a stellar model involving pulsations able to generate such signals. A further development of paramount importance will also be the application of the present ideas to cosmological dynamics, given that the theory so naturally includes isotropic Hubble expansion within its possibilities.