the interaction between high - energy jets from active galactic nuclei ( agn ) and their surroundings has long been a topic of great astrophysical interest .
it is well known that agn feedback can control the size of a galaxy by influencing star formation , but the mechanism behind this is not well understood .
several recent observations @xcite , as well as numerical studies ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) , have demonstrated that agn feedback can be `` negative '' or `` positive . ''
the exact astrophysical conditions in the jet and cloud are important in determining what direction feedback takes .
jets can be roughly divided according to their fanaroff - riley classification @xcite .
fast , energetic frii jets seem more likely to result in negative feedback . negative feedback curbs or even halts star formation , and is thought to result from the extreme radiative and kinetic energies of the jet , which heat and disperse the star - forming gas . additionally , the kinetic energy of a jet creates turbulence that can prevent ambient gas from cooling and subsequently coalescing ( e.g. * ? ? ?
for example , a study of the system 3c 326 by @xcite found that , despite the strong h@xmath1 line emission and an inferred molecular gas mass of 2 @xmath2 m@xmath3 , the star formation is 20 times lower than predicted by the kennicutt - schmidt law .
they infer that turbulent heating from the jet is inhibiting star formation .
other studies , though , suggest increased star formation may be seen in the cocoon of such jets ( e.g. * ? ? ?
in contrast , fri jets propagate through the ism / igm with energies high enough to create compression in the surrounding gas , but low enough to reduce the chance of significant turbulent heating .
these jets are observed in positive feedback cases , wherein the effect of the jet serves to enhance star formation , including centaurus a ( * ? ? ?
* and references therein ) , 4c 41.17 @xcite , and minkowski s object ( hereafter m.o . )
* hereafter c06 ) .
we give more details on both positive and negative feedback in section [ sec : jetcloud ] .    in this paper , we focus on the case of m.o . , a peculiar star forming object located at a redshift of @xmath4 ( c06 ) that is currently being bombarded by a fr i radio jet from the nearby galaxy ngc 541 .
system is of particular interest due to the lack of evidence for an especially dense ism or igm .
there is also not much evidence for cold gas outside of the jet interaction site as is the case in centaurus a. as a result , it is unlikely that significant star formation would proceed in m.o . without the interaction of the jet .    a strong argument in favor of jet - induced star formation in m.o .
is the morphology of the jet - cloud interaction site . outside of the jet interaction , the gas in m.o .
is warm ( @xmath5 k ) and clumpy . near the jet interaction site
there is a double structure of gas wrapped around the jet and numerous h@xmath1 regions ( c06 ) .
c06 thus determined that it is likely that the jet interaction in m.o .
caused the warm gas to cool into the h@xmath0 regions straddling the jet , in contrast to the pre - existing cold gas regions in centaurus a. also , the star forming regions in m.o . correlate with the jet - cloud morphology ; the region where the star formation is the highest is the center of the jet - cloud interaction , and the star formation rate ( sfr ) decreases laterally from this point ( c06 ) . unless m.o . happens to be a very unique system
, then it may be that jet - induced star formation is a more general phenomenon than previously thought .
the present work can be viewed as an extension of our earlier study of the interactions of radiative shocks with clouds @xcite .
that work focused mostly on the effects of planar shocks overtaking individual ( or small collections of ) warm clumps on the scale of @xmath6 pc . in the current work ,
we explore the much richer problem of a full jet intersecting an inhomogeneous intergalactic cloud on the scale of tens of kpc .
the paper proceeds as follows : section 2 covers the theory behind jet - cloud interactions , section 3 describes the numerical models used to capture the m.o .
system , section 4 details the simulation results , and section 5 concludes the findings .
the basic idea of jet - induced star formation ( i.e. positive feedback ) is that the collision of the jet with the cloud will trigger a series of shocks within the cloud .
the immediate effect of these shocks will be to compress and heat the gas . depending on how the radiative processes scale with density and temperature
, the net result can be to dramatically increase the radiative efficiency within the cloud .
if the temperature dependence is shallower than the density dependence , the cloud can enter a phase of runaway cooling .
this process occurs most quickly in relatively over - dense regions of the original cloud .
these over - dense regions then proceed to collapse at an accelerating pace .
provided some of these clumps start sufficiently close to the jeans limit , this collapse will push them beyond this limit , such that gravitational collapse can take over and the clump will proceed to form stars .
the properties of the jet and cloud are key to controlling this process . for positive feedback to be important
, the initial cloud must be dense enough for some parts to be reasonably close to the jeans limit .
the temperature must also be such that any increase in temperature is met with a dramatic increase in cooling ( the hydrogen cooling edge at @xmath5 k is a good example ) .
it generally helps for the jet to be significantly less dense than the cloud .
finally , the jet velocity needs to be fast enough to trigger shocks in the cloud , yet not so fast that the cloud is disrupted before cooling can have much of an effect .
our numerical simulations are performed using the well - tested _
cosmos++ _ computational astrophysics code @xcite , specifically its newtonian hydrodynamics solvers .
_ cosmos++ _ carries over many of the multi - physics capabilities found in its predecessor code _
cosmos _ @xcite .
the newtonian solvers have previously been utilized to study the bar mode instability in magnetized , rotating neutron stars @xcite and the galactic center g2 event @xcite .
the current work uses the high resolution shock capturing ( hrsc ) scheme , which was described in its relativistic form in @xcite .
as there are few differences between the relativistic and newtonian forms , we do not give a full presentation here , focusing instead on the packages that are most important to this paper : chemistry , cooling , and star formation .
note that , although magnetic fields can play an important role in shock - induced star formation ( cf .
* ) , they are not considered in this work .
most prior numerical studies of stimulated star formation from jet - generated shocks have been hampered by the resolution of the computational mesh ( e.g. * ? ? ?
* ) . in the present work ,
we significantly improve on previous resolution limitations by employing the adaptive mesh refinement ( amr ) capabilities of _ cosmos++_. _ cosmos++ _ employs a _
amr scheme , in which refinement and de - refinement decisions are made on a cell - by - cell basis , using an oct - tree network to traverse the grid hierarchy .
each level of refinement doubles the spatial resolution in each dimension within a given parent cell .
this style of local amr scheme ensures that the refinement and de - refinement conform as closely as possible to the shape of the region of interest , in this case the shocks and unstable cooling fronts triggered inside the cloud by the jet .    in this work
, we present an idealized case of a direct , axially symmetric collision between a jet and a pre - existing spherical cloud . as such , it is a simple , and well controlled , test simulation , though its correspondence with m.o .
is imperfect . in the case of m.o . , it is thought that the collision was between a jet and a stellar bridge connecting the elliptical galaxy , ngc 541 , with the interacting galaxies , ngc 545/547 ( c06 ) .
there is also evidence that the jet is slowly sweeping across m.o .
( c06 ) . despite these differences , the correspondence in parameters between our simulation and m.o .
ensure that our results are applicable and the simulation can be used to better understand the dynamics of this particular object and of jet - induced star formation more generally .
although we performed some two - dimensional simulations to test different code options and explore our parameter space , we mostly focus on reporting the results of our 3d simulations .
the 3d simulations have a base resolution of either @xmath7 or @xmath8 zones to cover a domain that is approximately @xmath9 with reflection boundaries applied in the @xmath10- and @xmath11-directions , so that we only simulate one quadrant of the full problem .
the finest spatial resolution achieved is 19.5 pc per zone , reached by including 2 levels of refinement on top of the @xmath8 base mesh , equivalent to a uniform mesh of @xmath12 zones .
the criterion used for refinement is that any zone with @xmath13 @xmath14 is kept at the maximum refinement , while zones that fall below @xmath15 @xmath14 are allowed to de - refine , provided neighboring zones never differ by more than one level of refinement and no zone is allowed to drop below the base resolution .
zones are checked against the refinement and de - refinement criteria once every ten evolution steps in the numerical code .
we find that a minimum resolution close to our base value is required even in the background in order to get reasonable convergence in the star formation rate ( as will be discussed in section [ sec : sfr ] ) .
the cloud , which represents the parent object of m.o .
, is initialized with a radius of @xmath16 pc .
it is , therefore , somewhat smaller , in terms of projected area , than the real value of 275 kpc@xmath17 ( c06 ) .
the gas within the cloud is distributed using a random , log - normal density , @xmath18 @xmath14 , with mean and median values of @xmath19 @xmath14 and @xmath20 @xmath14 , respectively , where @xmath21 is a randomly drawn variable with a mean of 0 and variance of 1 .
the gas has a mean molecular weight of @xmath22 , appropriate for approximately solar metallicity gas .
the total mass of our cloud ends up being @xmath23 , about a factor of four more than the estimated mass of m.o .
( @xmath24 ; c06 ) .
the cloud is immersed in a background gas with @xmath25 @xmath14 and @xmath26 k. the cloud and the background are initially assumed to be in pressure equilibrium , which is used to set the initial temperature of the cloud .
this gives a value around @xmath27 k. each simulation is run for five sound - crossing - times , @xmath28 , of the cloud , which corresponds to about 73 myr total . to avoid numerical problems caused by the background gas density dropping too far below its initial value
, we set a density floor of @xmath29 @xmath14 .
a jet is introduced into the simulation domain through one end . for this work
, we choose @xmath30 @xmath14 ( with @xmath31 the same as the cloud ) , a jet diameter of @xmath32 pc , and a jet velocity of 10% of the speed of light . for a cylindrical jet ,
this yields a kinetic power , @xmath33 erg s@xmath34 , where @xmath35 is the cross - sectional area of the jet .
our dependence on these model parameters is explored in section [ sec : jet_power ] .      in our simulations , we follow the abundances of 9 atomic and molecular species : , , , , , @xmath36 , h@xmath37 , h@xmath0 , and h@xmath38 .
the evolution of each species is governed by equations of the form @xmath39}}{\partial t } + \nabla \cdot ( \rho^{[m ] } \mathbf{v } )          = { \sum_{i=1}^{n_{s}}}{\sum_{j=1}^{n_{s } } }                               { { k_{ij}(t)}{\rho^{[i]}}{\rho^{[j ] } } }              + \sum_{i=1}^{n_{s } } { i_i ( \nu ) \rho^{[i]}}. \label{eqn : dens_m}\ ] ] these rate equations are solved using a stable , semi - implicit , backward difference scheme that we developed in @xcite , and which has since become a standard method adopted by the general community ( e.g. * ? ? ?
* ) due to its combination of robustness , efficiency , and accuracy .
a total of 27 gas - phase chemical reactions are included in the full network ( @xmath40 ) , including 19 collisional and 8 photoionization / photodissociation processes .
the exact reaction chains are spelled out in @xcite .
the photoionization field is @xmath41 s@xmath34 , appropriate for cosmic uv background radiation at low redshift @xcite .
the photodissociation rate is @xmath42 s@xmath34 , appropriate for the local interstellar medium @xcite .
the energy equation @xmath43          = -\lambda(t , n^{[m]})~ , \label{eqn : dual_en}\ ] ] where @xmath44 is the total energy density , including the internal and kinetic contributions , accounts for the cooling and heating of the gas via a total of seven different mechanisms : collisional - excitation , collisional - ionization , recombination , bremsstrahlung , metal - line cooling ( dominantly carbon , oxygen , neon , and iron ) , molecular - hydrogen cooling , and photoionization heating . the cumulative cooling function is @xmath45 } ) =      \sum_{i=1}^{n_{s } } \sum_{j=1}^{n_{s } }   \dot{e}_{ij}(t ) n^{[i ] } n^{[j ] }      - \sum_{i=1}^{n_{s } } j_i n^{[i ] }      + \dot{e}_m(t ) n^2 ~ , \label{eqn : coolingchem}\ ] ] where @xmath46 are the cooling rates from 2-body interactions between species @xmath47 and @xmath48 and @xmath49 represents the frequency - integrated photoionization and photodissociation heating rates .
@xmath50 is the temperature - dependent cooling rate for metals ( including carbon , oxygen , neon , and iron lines , assuming solar metallicity ) , taken from @xcite .
a cooling floor is set at @xmath51 k , below which only adiabatic cooling is possible . in practice , we rarely see temperatures below @xmath52 k.      in this study , we are particularly interested in tracking the formation of stars within the cloud .
we follow the approach of @xcite in defining a density threshold , @xmath53 , above which star formation is triggered at a rate given by @xmath54 where @xmath55 is the local free - fall time and @xmath56 controls the star formation efficiency .
we use values @xmath57 @xmath14 and @xmath58 for most of this work , though in section [ sec : sfr ] we explore how these choices affect our results .
this star formation appears as a sink term in the continuity equation @xmath59 by tracking how much mass is lost to star formation during each compute cycle , we are able to continuously track the star formation rate throughout the simulation .
we track the most massive ( @xmath60 ) `` stars '' created in this way using tracer particles .
each star particle is given an initial velocity equal to the velocity of the gas in the zone in which the star is created .
each star particle s velocity remains constant for the remainder of the simulation , as these particles are not acted on by any force .
this is not entirely realistic as some momentum may get redistributed during the free - fall process that leads to star formation . on longer timescales , gravitational forces from the rest of the cloud
will also act on the stars , eventually altering their velocities .
these effects are not considered in this work .    at the temperatures typical of star formation in our simulations ( @xmath61 k )
, the jeans length is @xmath62 which is considerably less than the radius of the cloud , yet well above the resolution limits of our simulations .
however , since self - gravity is not included in our current simulations , the clumps can not truly become jeans unstable .
table [ tab : params ] summarizes the 15 simulations presented in this work .
the naming convention refers to the base resolution and how many total grid resolution levels there are .
figure [ fig : volume ] presents a volume visualization from the end time of our highest resolution 3d simulation ( 384x128x128_3level ) .
a number of general features are apparent .
first , the cloud is dense enough and has enough inertia to dramatically slow the propagation of the jet ( represented by its temperature in red ) . at the speed
the jet is traveling , if not for the cloud ( and , to a lesser extent , the background gas ) impeding its progress , it should have traversed 74.6 box lengths ( 298.4 cloud radii ) over the duration of the simulation .
instead , most of the jet material is deflected to the sides of the cloud , though a significant fraction of its energy is deposited within the cloud gas .
the green material represents neutral hydrogen ( ) and shows that the core of the original cloud remains relatively intact .
some cooling ( represented by the blue , h@xmath0 , gas ) is happening on the fringes of the cloud ( at a radius of about @xmath63 in figure [ fig : volume ] ) , where weak shocks are being driven into the cloud by what is effectively now an over - pressure cocoon around the cloud .
more substantial cooling is happening toward the front of the cloud in the post - shock gas ( dense blue filament just downstream of where the jet terminates ) .
cccccccc [ tab : params ] 288x96_2level & @xmath64 & 2 & 4 & 0.1 & 1.5 & 0.1 & 5 + 288x96_3level & @xmath64 & 3 & 4 & 0.1 & 1.5 & 0.1 & 5 + 288x96_4level & @xmath64 & 4 & 4 & 0.1 & 1.5 & 0.1 & 5 + 288x96_5level & @xmath64 & 5 & 4 & 0.1 & 1.5 & 0.1 & 5 + 288x96_3level_n1_e.05 & @xmath64 & 3 & 1 & 0.05 & 1.5 & 0.1 & 5 + 288x96_3level_n4_e.05 & @xmath64 & 3 & 4 & 0.05 & 1.5 & 0.1 & 5 + 288x96_3level_n7_e.05 & @xmath64 & 3 & 7 & 0.05 & 1.5 & 0.1 & 5 + 288x96_3level_rhoj10 & @xmath64 & 3 & 4 & 0.1 & 1.5 & 0.1 & 50 + 288x96_3level_dj2 & @xmath64 & 3 & 4 & 0.1 & 3.0 & 0.1 & 5 + 288x96_3level_vj5 & @xmath64 & 3 & 4 & 0.1 & 1.5 & 0.5 & 5 + 384x128_3level & @xmath65 & 3 & 4 & 0.1 & 1.5 & 0.1 & 5 + 288x96x96_2level & @xmath66 & 2 & 4 & 0.1 & 1.5 & 0.1 & 5 + 288x96x96_3level & @xmath66 & 3 & 4 & 0.1 & 1.5 & 0.1 & 5 + 384x128x128_2level & @xmath8 & 2 & 4 & 0.1 & 1.5 & 0.1 & 5 + 384x128x128_3level & @xmath8 & 3 & 4 & 0.1 & 1.5 & 0.1 & 5 +     myr ) of our highest resolution 3d simulation ( 384x128x128_3level ) .
red represents hot , @xmath67 k , material , mostly the jet , although possibly with some cloud gas mixed in ; green represents regions of the cloud with a neutral hydrogen density , @xmath68 @xmath14 ; and blue represents regions of the cloud with a molecular hydrogen ( our tracer for cold gas ) density , @xmath69 @xmath14 .
data have been reflected across the @xmath70 and @xmath71 planes to create this image .
axes are marked in units of cloud radii , @xmath72 .
[ fig : volume ] ]    figure [ fig:2d ] shows what the evolution looks like in one of our 2d simulations ( 384x128_3level ) .
this evolution shows broad similarities to what is seen in the 3d simulations , although there are notable differences . for one ,
the jet has penetrated much further into the cloud in 2d by the end of the simulation .
we will also discuss in section [ sec : cooling_front ] how the cooling and star formation in 2d take place in a `` wall '' proceeding forward through the cloud , while in 3d , they take place in a ring , centered about the jet axis . in section [ sec : star_particles ] , we also find that the star particles have a broader velocity distribution in 2d than 3d .
such differences illustrate the necessity of 3d simulations .    , i.e. red represents @xmath67 k , green represents @xmath73 @xmath14 , and blue represents @xmath74 @xmath14 .
the panels show the evolution , spaced at intervals of one cloud - crossing time , i.e. @xmath75 .
the dramatic change in the size of the cloud between the first and second panels is due to photoionization .
the overall size and density of the cloud does not change much , but the ionization levels do .
[ fig:2d],title="fig : " ] , i.e. red represents @xmath67 k , green represents @xmath73 @xmath14 , and blue represents @xmath74 @xmath14 .
the panels show the evolution , spaced at intervals of one cloud - crossing time , i.e. @xmath75 . the dramatic change in the size of the cloud between the first and second panels is due to photoionization .
the overall size and density of the cloud does not change much , but the ionization levels do .
[ fig:2d],title="fig : " ] , i.e. red represents @xmath67 k , green represents @xmath73 @xmath14 , and blue represents @xmath74 @xmath14 .
the panels show the evolution , spaced at intervals of one cloud - crossing time , i.e. @xmath75 . the dramatic change in the size of the cloud between the first and second panels is due to photoionization .
the overall size and density of the cloud does not change much , but the ionization levels do .
[ fig:2d],title="fig : " ] + , i.e. red represents @xmath67 k , green represents @xmath73 @xmath14 , and blue represents @xmath74 @xmath14 .
the panels show the evolution , spaced at intervals of one cloud - crossing time , i.e. @xmath75 . the dramatic change in the size of the cloud between the first and second panels is due to photoionization .
the overall size and density of the cloud does not change much , but the ionization levels do .
[ fig:2d],title="fig : " ] , i.e. red represents @xmath67 k , green represents @xmath73 @xmath14 , and blue represents @xmath74 @xmath14 .
the panels show the evolution , spaced at intervals of one cloud - crossing time , i.e. @xmath75 . the dramatic change in the size of the cloud between the first and second panels is due to photoionization .
the overall size and density of the cloud does not change much , but the ionization levels do .
[ fig:2d],title="fig : " ] , i.e. red represents @xmath67 k , green represents @xmath73 @xmath14 , and blue represents @xmath74 @xmath14 .
the panels show the evolution , spaced at intervals of one cloud - crossing time , i.e. @xmath75 . the dramatic change in the size of the cloud between the first and second panels is due to photoionization .
the overall size and density of the cloud does not change much , but the ionization levels do .
[ fig:2d],title="fig : " ]      figure [ fig : shocks ] is a volume visualization showing regions of significant velocity convergence ( @xmath76 ) , as an indicator of where shocks are located . from it , we see that , even after 73 myr , the jet - driven shock still has not reached the center of the cloud . we can compare this plot with the estimated shock velocity in the cloud . for our chosen jet velocity , @xmath77 km s@xmath34 , we can use the shock jump conditions , @xmath78 where @xmath79 and @xmath80 are measured in the rest frame of the pre - shock gas , to estimate the speed of the shocks in the background gas , @xmath81 km s@xmath34 . if the shock in the background is strong , then the postshock pressure is approximately @xmath82 .
if we assume the shocks are also strong inside the cloud and that the postshock background and cloud gas reach pressure equilibrium , then we can estimate the speed of the shock in the cloud @xmath83 for our median cloud density , @xmath84 @xmath14 , this gives us a shock speed in the cloud of @xmath85 km s@xmath34 .
thus , the shock should have traveled approximately 22% of the way through the cloud by the end of our simulation , roughly consistent with figure [ fig : shocks ] .
as expected , material that is ahead of the shock ( downstream ) remains relatively undisturbed , except on the fringes of the cloud where the over - pressure cocoon is driving weaker shocks into the cloud ( blue outline of the cloud in figure [ fig : shocks ] ) .    ) , indicative of shocks , from the final time dump ( @xmath86 myr ) of our highest resolution 3d simulation ( 384x128x128_3level ) .
yellow , orange , and red shades indicate the strong shocks associated with the jet , while blue indicates much weaker shocks .
data have been reflected across the @xmath70 and @xmath71 planes to create this image .
axes are marked in units of cloud radii , @xmath72 .
[ fig : shocks ] ]      as mentioned in section [ sec : jetcloud ] , for the jet feedback to be positive , it is critical for the cooling timescale to be short compared to the shock - crossing time of the cloud .
following @xcite , we estimate the cooling time to be @xmath87 this is _ much
_ shorter than any other relevant timescale in the problem .
interestingly , the greatest amount of cooling ( and associated star formation ) is seen in a ring radiating away from the initial impact point at the head of the cloud . in figure
[ fig : volume ] , you can see evidence of this in the light green ring of just to the left of the cloud s center .
a smaller , denser ring of h@xmath0 follows behind the ring .
figure [ fig : formation ] shows how this cooling front and associated star formation moves over the course of the simulation .
this distribution is significantly different than what one sees in 2d simulations ( figure [ fig:2d ] ) , where cooling and star formation occurs more uniformly across the cloud diameter ( perpendicular to the jet propagation ) .    ) at @xmath88 ( orange ) , 44 ( green ) , 58 ( blue ) , and 73 myr ( purple ) from our highest resolution 3d simulation ( 384x128x128_3level ) .
the left panel shows a side view , while the right panel shows a head - on view .
data have again been reflected across the @xmath70 and @xmath71 planes .
[ fig : formation],title="fig : " ] ) at @xmath88 ( orange ) , 44 ( green ) , 58 ( blue ) , and 73 myr ( purple ) from our highest resolution 3d simulation ( 384x128x128_3level ) .
the left panel shows a side view , while the right panel shows a head - on view .
data have again been reflected across the @xmath70 and @xmath71 planes .
[ fig : formation],title="fig : " ]    figure [ fig : slices ] helps clarify the nature of these rings . note first that , as expected , there is a compression front moving forward through the body of the cloud ( approximately vertical feature extending from @xmath89 to @xmath90 in both panels ) .
there is also a compression front moving radially into the cloud from all sides ( arced feature extending from @xmath91 to @xmath92 in both panels , although it extends around the entire cloud ) , driven by the over - pressured cocoon .
cooling seems to be operating effectively in both of these regions , as the temperature remains @xmath93 k ( blue regions in the right panel ) .
however , in most parts of these fronts , the density does not quite reach the threshold for star formation ( @xmath94 ) .
this does _ not _ appear to be a resolution effect as the fronts are well resolved in most places and our overall s.f.r . seems to have converged ( see figure [ fig : sfr ] ) .
instead , the ring of star formation seen in figures [ fig : formation ] and [ fig : stars ] seems to correspond to regions where the forward propagating compression wave intersects the radial compression front ( most strongly around @xmath95 ) .
the combined effect of the two shocks is enough to push the density over the star formation threshold .    .
the right panel shows temperature in k. note how the densest , coldest gas is found around @xmath96 , where the nearly planar front moving forward through the cloud intersects the nearly spherical front being driven into the cloud .
[ fig : slices],title="fig : " ] .
the right panel shows temperature in k. note how the densest , coldest gas is found around @xmath96 , where the nearly planar front moving forward through the cloud intersects the nearly spherical front being driven into the cloud .
[ fig : slices],title="fig : " ]      as a reminder , whenever more than @xmath97 of gas is converted into stars within a given zone within a single cycle , then a star particle is created to track the properties , such as position , age , and velocity , of that `` star . ''
following this prescription , we created over @xmath98 star particles in our highest resolution simulation ( 384x128x128_3level ) .
the star particles span an age range from 0  72 myr , although the star particles associated with the jet interaction are all younger than about 50 myr ; the oldest stars formed spontaneously at the start of the simulation inside zones whose random density exceeded the star formation threshold .
even the 50 myr age is considerably older than the stellar age range of 7.0 - 8.5 myr estimated for m.o .
much of this discrepancy can be explained by the delay associated with the free - fall time , which is not accounted for in the `` ages '' of our star particles . for our threshold density , @xmath53 ,
the associated free - fall time is @xmath99 myr , bringing the ages more in line with m.o .
furthermore , the stellar ages in c06 may actually be telling us more about the age of the jet .
it may be that m.o .
has not yet evolved for as long as we ran our simulations .
an interesting point about the stellar ages in our simulations is that they show a negative curvature along the direction of jet propagation , that is , the youngest ( currently forming ) stars are found between two populations of slightly older stars , one in the upstream direction and one downstream .
this can be seen for one of our 2d simulations in the right panel of figure [ fig : stars_2d ] .
the upstream population are stars that formed recently , only slightly before the current star formation .
the downstream population are some of the first stars to form from the jet interaction , but they are now actually located ahead of the current star formation front because they received a velocity kick larger than the current shock speed . since the star particles are not coupled to the gas , they can actually pass ahead of the shock as it slows down , giving an apparently older population ahead of the current star formation front .
myr ) of our 384x128_3level 2d simulation . in the left panel ,
the color of each particle indicates its velocity magnitude ( in km s@xmath34 ) , while in the right panel , the color indicates age ( in myr ) .
[ fig : stars_2d],title="fig : " ]   myr ) of our 384x128_3level 2d simulation . in the left panel ,
the color of each particle indicates its velocity magnitude ( in km s@xmath34 ) , while in the right panel , the color indicates age ( in myr ) .
[ fig : stars_2d],title="fig : " ]    by using our results from sec .
[ sec : shock ] and inverting the shock jump condition , we predict a post - shock velocity in the cloud of @xmath100 km s@xmath34 , consistent with the observed 40 km s@xmath34 shear velocity in m.o .
( c06 ) . since
our star particles are assigned their velocity based upon the velocity of the gas from which they form , we expect the simulated star particle velocities to be similar .
figure [ fig : histogram ] shows a histogram of the velocity distribution of all of the star particles .
the distribution peaks at around 50 km s@xmath34 , which is fairly consistent with our crude predictions and the observations of m.o .
figures [ fig : stars_2d ] ( left panel ) and [ fig : stars ] show the spatial distributions of the star particles , colored by the magnitude of their velocities , in 2d and 3d , respectively .
note the significant differences in the distributions , although the magnitudes of the velocities ( and ages of stars ) are comparable . as mentioned above , the fastest moving star particles are on the downstream edges of the distributions .
they are also found further out , i.e. at larger radii , on average in the 3d simulation .
bins are filled with stars that formed in the initial cloud , prior to its interaction with the jet .
the peak of the 3d distribution is roughly consistent with observations of m.o .
[ fig : histogram],title="fig : " ]   bins are filled with stars that formed in the initial cloud , prior to its interaction with the jet .
the peak of the 3d distribution is roughly consistent with observations of m.o .
[ fig : histogram],title="fig : " ]     myr ) of our highest resolution 3d simulation ( 384x128x128_3level ) .
the color of each particle indicates its velocity magnitude ( in km s@xmath34 ) .
the left panel shows a side view , while the right panel shows a head - on view .
data have again been reflected across the @xmath70 and @xmath71 planes .
[ fig : stars],title="fig : " ]   myr ) of our highest resolution 3d simulation ( 384x128x128_3level ) .
the color of each particle indicates its velocity magnitude ( in km s@xmath34 ) .
the left panel shows a side view , while the right panel shows a head - on view .
data have again been reflected across the @xmath70 and @xmath71 planes .
[ fig : stars],title="fig : " ]      figure [ fig : sfr ] shows the measured star formation rate ( s.f.r . ) as a function of time for four different 3d simulations done at different effective resolutions from 57.1 pc per zone ( 288x96x96_2level ) to 19.5 pc per zone ( 384x128x128_3level ) . as might be expected , the simulated s.f.r .
is sensitive to resolution . at higher resolutions ,
more gas is going to be able to reach the density threshold , @xmath53 , of our star formation model .
however , at some point a limit should be reached where the high density filaments are well enough resolved that their size and peak density are no longer functions of resolution .
it appears we may have reached this point in our highest resolution simulation .
note that the peak star formation rate actually goes down compared to the next lower resolution .
also , the time at which star formation begins remains consistent at about 25 myr for the two highest resolutions .
a convergence at a resolution of about 25 pc per zone also seems to be supported by our 2d tests ( right panel of figure [ fig : sfr ] ) , where we were able to consider even higher resolutions .    .
together , the figures suggest that our results are reasonably well converged at resolutions above about 25 pc per zone .
note that the highest resolution 3d simulation produces a steady s.f.r .
of @xmath101 yr@xmath34 ( corresponding to an h@xmath102 luminosity of @xmath103 erg s@xmath34 ) .
this falls short of the observed rate in m.o .
( grey , dashed line ) by about a factor of 6 . the units on the s.f.r . axis in 2d are arbitrary because there is no meaningful way to assign a volume to the cloud .
[ fig : sfr],title="fig : " ] .
together , the figures suggest that our results are reasonably well converged at resolutions above about 25 pc per zone .
note that the highest resolution 3d simulation produces a steady s.f.r .
of @xmath101 yr@xmath34 ( corresponding to an h@xmath102 luminosity of @xmath103 erg s@xmath34 ) .
this falls short of the observed rate in m.o .
( grey , dashed line ) by about a factor of 6 . the units on the s.f.r .
axis in 2d are arbitrary because there is no meaningful way to assign a volume to the cloud .
[ fig : sfr],title="fig : " ]    the s.f.r . , or correspondingly the h@xmath102 luminosity , since @xcite @xmath104 in figure [ fig : sfr ] is about a factor of 6 less than the observed value in m.o .
( @xmath105 yr@xmath34 and @xmath106 erg s@xmath34 ; c06 ) .
this most likely means that one or more of our simulation parameters does not match the situation in m.o . to increase the s.f.r . in the simulation , we could : decrease the star formation density threshold , @xmath53 ; increase the star formation efficiency , @xmath56 ; increase the jet power ; or increase the initial cloud density .
we use a set of 2d simulations to test the impact of changing @xmath53 , @xmath56 , and the jet power . in the rest of this section ,
we focus on @xmath53 and @xmath56 . in the next section , we consider varying the parameters that affect the jet power .
we do not , however , explore the effect of further increasing the cloud density , as it is already high enough to trigger spontaneous star formation .
we do not want such spontaneous star formation to overwhelm the jet - induced star formation that we are studying .
figure [ fig : sfrmodel ] shows the measured s.f.r . for four different 2d simulations using different values for @xmath53 and @xmath56 .
not surprisingly , if we lower the star formation threshold density ( black , dashed curve ) , then star formation occurs much more easily , so much so that the rate of spontaneous star formation ( that not associated with the jet interaction ) is within a factor of 10 of the jet - induced rate .
interestingly , cutting the star formation efficiency in half ( green , dot - dashed curve vs. blue , dotted curve ) does not simply cut the s.f.r . in half .
this illustrates the non - linearity of this star formation process and the need for high fidelity simulations .     and
@xmath56 , for 2d simulations with a base resolution of @xmath64 and 3 total resolution levels .
the dotted , blue curve shows results with our default parameters .
the units on the s.f.r .
axis are arbitrary because in 2d , there is no meaningful way to assign a volume to the cloud .
[ fig : sfrmodel ] ]    ultimately , we find that adjusting the star formation parameters does not really aid in raising our measured s.f.r .
our standard parameters , @xmath57 @xmath14 and @xmath58 , already produce the highest induced s.f.r . of any of these simulations . increasing the star formation efficiency
further might yield some further increase in the s.f.r . , but we are already at what we consider a reasonable upper limit .      using the known correlation between jet power and radio luminosity , @xmath107 @xcite , and the radio luminosity of ngc 541 , @xmath108 erg s@xmath34 @xcite
, we estimate the kinetic power of the jet impacting m.o . to be @xmath109 erg s@xmath34 .
by contrast , our default jet parameters give a power that is more than two orders of magnitude lower . in figure
[ fig : sfrjet ] , we explore how increasing the jet power would impact our results .
increasing the jet density ( black , dashed curve ) increases the s.f.r . and also causes star formation to start earlier .
this is because a denser jet more easily penetrates the background gas to reach the cloud .
increasing the jet diameter ( green , dot - dashed curve ) also increases the s.f.r . , although this happens gradually , not even peaking by the end of the simulation .
finally , increasing the jet velocity ( red , double - dot - dashed curve ) again causes star formation to start earlier , but also makes the shock crossing time comparable to the cooling timescale , such that the cloud is disrupted ( and ultimately blown off the simulation domain ) before a significant number of stars can form ( s.f.r .
goes to zero around @xmath110 myr ) .    , @xmath111 , and @xmath112 , for 2d simulations with a base resolution of @xmath64 and 3 total resolution levels .
the dotted , blue curve shows results with our default parameters .
the units on the s.f.r .
axis are arbitrary because in 2d , there is no meaningful way to assign a volume to the cloud .
[ fig : sfrjet ] ]    from these tests , it appears that the best option for matching the s.f.r . in m.o . would be a simulation with a denser jet . extrapolating from our 2d results , and assuming the same increase is seen in 3d
, it appears that a factor of 25 increase in jet density might be enough to have the peak s.f.r . in our simulation
match the currently observed value in m.o .
this would also bring the kinetic power of our jet more in line .
figure [ fig : history ] tracks the total mass of , h@xmath0 , and stars over the course of our highest resolution simulation . while the simulation produces roughly the correct amount of [ @xmath113 in the simulation vs. @xmath24 for m.o .
( c06 ) ] , it produces far less h@xmath0 [ @xmath114 vs. @xmath115 in m.o .
( lacy et al . , in prep . ) ] , and about a factor of 10 fewer stars [ @xmath116 vs. @xmath117 in m.o .
( c06 ) ] .
we find an electron number density in our cloud of @xmath118 @xmath14 , also somewhat lower than the range of 1 - 10 @xmath14 obtained for m.o .
since the amount of gas is comparable , this suggests that the cooling and star formation in the simulations are not nearly as efficient as they are in m.o . the sizable discrepancy in h@xmath0 could be reduced if we included the effects of self - shielding , which would reduce the amount of photodissociation and increase the yield of h@xmath0 . however , even in simulations where we completely ignore photodissociation , we still find the h@xmath0 mass to be at least two orders of magnitude too low .
another important point about our h@xmath0 values is that we are only forming h@xmath0 through gas phase reactions . in m.o .
, it may be that substantial h@xmath0 formation takes place on dust grains .
however , simulations that include precipitation of h@xmath0 on dust , with standard grain parameters , show only a very modest increase in @xmath119 .
likely , further work will be needed to better understand the h@xmath0 levels in m.o .    , @xmath119 , @xmath120 , and @xmath121 over time for our highest resolution 3d simulation ( 384x128x128_3level ) .
[ fig : history ] ]    figure [ fig : history ] also tracks the `` star formation efficiency , '' @xmath121 , over time . by this measure , our simulation achieves a peak star formation efficiency of 0.6% .
this , too , is lower than the value of 4% measured in m.o .
all of these discrepancies can be traced to the low s.f.r .
seen in our simulations . as mentioned previously ,
the _ peak _ s.f.r . in our simulations ( see figure [ fig : sfr ] ) is about a factor of 6 less than the _ current _ s.f.r in m.o .
, though we showed in the previous section that increasing the jet power in the simulations might be able to remove this discrepancy .
in this work , we use 2d and 3d simulations to model the evolution of a radio jet impacting a single , dwarf - galaxy - scale cloud in a direct , axially - symmetric collision .
our intention is to use these simulations to better understand observations of m.o . , the peculiar starburst galaxy located at the termination point of the radio jet from ngc 541 .    our first , main conclusion is that jet - induced star formation ( i.e. positive feedback ) is possible under this scenario .
figures [ fig : sfr][fig : history ] all show dramatic increases in star formation attributable to the jet interaction .
interestingly , we find that most of the star formation in our 3d simulations takes place in a propagating ring centered about the jet axis ( figure [ fig : formation ] ) .
the ring lies at the intersection of a forward propagating shock driven into the cloud by the jet and a spherical compression shock driven by the over - pressure cocoon ( figure [ fig : slices ] ) .
these rings may potentially explain the torus of star formation observed in m.o .
importantly , this star formation is occurring upstream of an ring ( compare figures [ fig : volume ] and [ fig : formation ] ) , also consistent with observations ( c06 ) .    additionally , the velocity histogram of our star particles peaks at @xmath122 km s@xmath34 ( left panel of figure [ fig : histogram ] ) , which is again consistent with observations ( c06 ) . the spatial distribution of the star particle velocities is also interesting .
the fastest moving star particles are found furthest downstream ( figures [ fig : stars_2d ] and [ fig : stars ] ) .
these are also some of the first star particles to form ( i.e. they are the oldest ) , which is consistent with a slowing of the propagation speed of the star formation front .    despite these successes , the correlation of our results with m.o .
is imperfect .
the main discrepancy is that our measured star formation rate is about a factor of 6 lower than the observed value .
c06 also suggested that the star formation rate in m.o . may have been even higher in the past .
along with not producing the expected mass of stars , our simulations also fall short on the amount of h@xmath0 they produce .    in the paper ,
we have tried to explore possible avenues to rectify these discrepancies .
we first considered our star formation model ; there we found that our results are not strongly sensitive to its parameters .
further , we are already close to what we consider a reasonable upper limit for the star - formation efficiency , so there is little chance to make up the discrepancy that way .
the sensitivity of our results to the jet parameters was another matter .
first , we found that the jet power may have been lower than what is inferred in the case of m.o .
by as much as two orders of magnitude . in our 2d simulations , though , we found a complicated dependence of the s.f.r . on jet power
( figure [ fig : sfrjet ] ) . while increasing the jet diameter by a factor of two or the jet density by a factor of ten roughly doubled the s.f.r .
, increasing the jet velocity by a factor of 5 actually lowered the s.f.r . and destroyed the cloud on a relatively short timescale .
a similar strong sensitivity to the shock velocity was seen in our earlier work @xcite . going the other direction and significantly lowering
the jet velocity would introduce a different problem
star formation would likely proceed too slowly .
once the first generation of massive stars form , there is only a limited amount of time ( of the order a few myr ) for star formation to continue before negative feedback , in the form of heating from the first massive stars , would effectively shut it off @xcite .
this effect might explain why the s.f.r . in m.o .
appears to be lower now that it was in the past  negative feedback may already be kicking in .
as future observations continue to constrain the star formation history of m.o . , we plan to continue to refine our simulations .
future modifications to our setup may include : simulating non - axially - symmetric interactions between the jet and cloud ; having the jet sweep across the cloud ; or adding more feedback mechanisms , such as heating from young stars and supernovae .
this work used the extreme science and engineering discovery environment ( xsede ) , which is supported by national science foundation grant number aci-1053575 .
pcf acknowledges support from national science foundation grants ast-1211230 and ast-1616185 .
jwlw acknowledges support from nrao student observing support grant sospa3 - 020 .
work by pa was performed in part under the auspices of the u.s .
department of energy by lawrence livermore national laboratory under contract de - ac52 - 07na27344 .


inspired by wheeler s quest for the existence of geons @xcite , boson stars were introduced by feinblum and mckinley @xcite , kaup @xcite , and bonazzola and ruffini @xcite . in these
boson stars the electromagnetic vector field was replaced by some tentative scalar field . with the discovery of a higgs - like boson last year at the lhc @xcite
the first fundamental scalar field has been found . but numerous scalar fields have been predicted to exist in high energy physics and cosmology .
boson stars arise as stationary localized solutions of the coupled einstein - klein - gordon equations @xcite .
the physical properties of boson stars depend strongly on the type of scalar field potential employed ( see e.g.  the review articles @xcite ) .
mini boson stars arise , when only a mass term is present but no self - interaction .
they can reach only relatively small masses .
larger boson stars are obtained , when a repulsive quartic self - interaction is included @xcite .
for these two types of boson stars gravity is a necessary ingredient .
in contrast , in the presence of a sextic potential one obtains solitonic boson stars @xcite , which possess a flat space - time limit , where they correspond to non - topological solitons @xcite ( or @xmath0-balls @xcite ) . moreover , these solitonic boson stars can reach even higher masses @xcite .
here we consider boson stars , which are compact in the sense , that the scalar field of these spherically symmetric configurations is finite inside a ball of radius @xmath3 , but vanishes identically outside this radius . in this respect
the compact @xmath0-balls resemble stars @xcite .
obtained from a v - shaped self - interaction potential , these compact boson stars also possess a flat space - time limit , compact @xmath0-balls @xcite .
these represent solutions of the signum - gordon equation .
the study of scalar fields with a v - shaped self - interaction potential has revealed interesting physical phenomena . when coupled to electromagnetism
, the balance of forces allows for shell - like configurations @xcite . in these @xmath0-shells
the scalar field vanishes identically both inside a certain radius @xmath4 and outside a certain radius @xmath3 . the scalar field
thus forms a shell of charged matter , @xmath5 .
when coupling these shells to gravity the resulting boson shells possess an empty minkowski space interior @xmath6 .
however , the shells need not be empty in their interior , they can harbour a black hole in there @xcite .
thus one finds that uniqueness and no - hair theorems for black holes can be avoided in the presence of boson shells @xcite .    whereas these previous studies considered only asymptotically flat solutions , we here include a cosmological constant . on the one hand , a positive cosmological constant is relevant from an observational point of view , since it can model the dark energy of the universe .
since boson stars are very compact objects that can possess very high densities , they have been suggested as alternatives to supermassive black holes , e.g.  in the center of galaxies @xcite .
even if that would be excluded by observations ( for a discussion see e.g.  @xcite ) boson stars could still act as toy models for very compact objects , e.g.  neutron stars .
such a model of a compact star in a space - time with positive cosmological constant should be a more realistic description of compact stars in the universe , since all observations seem to indicate the existence of a form of dark energy .
a negative cosmological constant , on the other hand , leads to solutions that can be interpreted within the ads / cft correspondence @xcite .
recently , the study of boson stars in ads space - time received increasing attention @xcite .
this is related to the fact that within the context of a holographic description of superconductors and superfluids @xcite ( for reviews see @xcite ) the formation of scalar hair on charged solitons in asymptotically ads has been interpreted as an insulator / superconductor phase transition @xcite .
the limit of setting the electric charge of the scalar field @xmath7 to infinity , which due to the scaling symmetries corresponds to setting newton s constant @xmath8 to zero , is called the `` the probe limit '' in this context .
we adapt this nomenclature here and refer to the case , where the matter field equation is solved in a fixed background , as `` the probe limit '' . in the opposite limit @xmath9 ,
the gauge symmetry becomes global and the resulting solutions are uncharged solitons in ads .
these are essentially uncharged boson stars and have been suggested to be the holographic description of glueball condensates @xcite .
boson stars in asymptotic ads are also of interest from another point of view .
it has been suggested that the dynamical formation of a black hole in ads is the dual description of thermalization in a strongly coupled quantum field theory .
as such the stability of ads space - time was studied with respect to perturbations , and it was conjectured that ads is unstable under arbitrarily small scalar perturbations and that eventually a black hole would form due to the reflection of the perturbations on the ads boundary @xcite . however , in @xcite it was shown that boson stars appear to be nonlinearly stable .
if that were true , thermalization in the dual field theory would not occur .
hence , boson stars in ads play an important rle in the context of the nonlinear ( in)stability of ads space - time , and thus of its dual description .    here
we first consider the set of boson star solutions for various space - time dimensions @xmath10 in the probe limit .
interestingly , the presence of a positive cosmological constant allows for the existence of boson shells without an electromagnetic field .
we note , that all objects constructed here are electrically neutral .    by solving the coupled set of the einstein - signum - gordon equations , we subsequently determine the domain of existence of the compact ( a)ds boson stars and shells in @xmath10 space - time dimensions .
we analyze their physical properties and briefly address the stability of the boson stars from a catastrophe theory point of view @xcite .
we also address astrophysical aspects of boson stars and boson shells in four dimensions for the physical value of the cosmological constant .
the paper is organized as follows . in section 2
we present the action , the ansatz , the equations of motion together with the scaling property , the boundary conditions and the global charges .
we present the solutions in the probe limit in section 3 .
the boson stars and shells obtained with the back reaction taken into account are discussed in section 4 .
we end with our conclusions and an outlook in section 5 .
we consider the action of a self - interacting complex scalar field @xmath11 coupled to einstein gravity in @xmath1 dimensions @xmath12 \sqrt{-g } d^dx   , \label{action}\ ] ] with curvature scalar @xmath13 , cosmological constant @xmath14 , newton s constant @xmath8 , and the asterisk denotes complex conjugation . the scalar potential @xmath15 is chosen as @xmath16    variation of the action with respect to the metric and the matter fields leads , respectively , to the einstein equations @xmath17 with stress - energy tensor @xmath18 and the matter field equation , @xmath19 where @xmath20 denotes the covariant derivative .
invariance of the action under the global phase transformation @xmath21 leads to the conserved current @xmath22 and the associated conserved charge @xmath0 .      to construct spherically symmetric boson star solutions
we employ schwarzschild - like coordinates and adopt the spherically symmetric metric @xmath23 where @xmath24 is the metric on the @xmath25 dimensional unit sphere .
the associated ansatz for the boson field takes the form @xmath26 with the frequency @xmath27 .
the conserved scalar charge @xmath0 @xmath28 is then proportional to @xmath27 .
we next introduce dimensionless quantities by @xmath29 thus the coupling strength of gravity is expressed in terms of the coupling constant @xmath30 .
this yields the set of equations @xmath31      let us now specify the boundary conditions for the metric and the boson field . for the metric function @xmath32 we adopt @xmath33 where @xmath34 is the outer radius of the boson star .
since it retains this value to infinity , this fixes the time coordinate .
for the metric function @xmath35 we require at the origin the regularity condition for globally regular ball - like boson star solutions @xmath36 and for globally regular shell - like solutions @xmath37 where @xmath38 is the inner radius of the shell .    for boson stars we require for the boson field function one condition at the origin and two conditions at the outer radius @xmath34 @xmath39 since this is one condition too many , we introduce another auxiliary differential equation , @xmath40 , by treating @xmath34 as a function .
thus @xmath41 is constant , but the value of the constant is adjusted in the numerical scheme such that the boundary conditions eqs .
( [ bc ] ) are satisfied .
this determines the outer radius of the star .    for boson shells , on the other hand
, we require at the inner radius @xmath42 and at the outer radius @xmath43 the conditions @xmath44 we now also make the ratio of inner and outer radius @xmath45 an auxiliary ( constant ) variable .    for the numerical computation we introduce the scaled coordinate @xmath46 , such that the inner radius is at @xmath47 and outer radius is at @xmath48 .
we refer to the solution in the exterior region @xmath49 as the outer solution of the boson stars and boson shells . in the asymptotically flat case ,
the outer solution is given by the schwarzschild solution .
in the presence of a cosmological constant the schwarzschild - de sitter and schwarzschild - anti - de sitter solutions @xmath50 are exact solutions of the odes in the exterior region .
hence the mass parameter of the solutions is given by @xmath51      analogously , we refer to the solution in the interior region @xmath52 as the inner solution of the boson shells . in the asymptotically flat case
, the regular inner solution corresponds to flat minkowski space , whereas in the presence of a cosmological constant the regular inner solutions correspond to either de sitter or anti - de sitter space .
note however , that in general @xmath53 .
thus a rescaling of the time coordinate , @xmath54 , is required to obtain the de sitter or anti - de sitter line element in the standard form .    in principle
, we could replace these regular inner solutions by the appropriate black hole solutions , analogously to the asymptotically flat case @xcite
. then these inner solutions would correspond to schwarzschild - de sitter or schwarzschild - anti - de sitter solutions .
here we present the families of solutions in the so called probe limit .
thus we obtain the solutions for vanishing coupling to gravity , i.e. , @xmath55 , in the respective background .      for vanishing cosmological constant the @xmath0-ball solutions can be found analytically and expressed in terms of bessel functions .
the @xmath56 solution was given in @xcite @xmath57 for @xmath10 dimensions this generalizes according to @xmath58 where @xmath59 .
the constant @xmath60 and the outer radius @xmath41 are determined by the conditions @xmath61 and @xmath62 .
the latter yields @xmath63 .
hence @xmath34 is the smallest non - vanishing zero of @xmath64 . from the first condition it then follows that @xmath65 .
thus we find @xmath66    the properties of the unique solution in @xmath56 dimensions was discussed in @xcite .
when going to higher dimensions , the properties of the respective solutions vary only slowly with @xmath1 , and likewise when going to @xmath67 .
this is seen in fig .
[ fig1 ] , when restricting to a vanishing cosmological constant , @xmath68 .
let us now consider the @xmath0-ball solutions in an ( anti-)de sitter background in @xmath1 dimensions . to obtain these solutions , we have solved the scalar field equation in the respective background numerically , employing a newton - raphson scheme .
when the scaled cosmological constant @xmath69 is varied , the solutions change smoothly from the minkowski background solutions . in fig .
[ fig1 ] we exhibit the dependence of the solutions on @xmath69 . here
we show the value of the scalar field at the origin @xmath70 together with the value of the outer radius @xmath41 of the solutions for @xmath67 , 4 , 5 and 10 dimensions .    as the scaled cosmological constant @xmath69 increases from zero , the value of the scalar field at the origin @xmath70 decreases along with the outer radius @xmath41 .
interestingly , there is a maximal value @xmath71 , which increases with the dimension @xmath1 , for which these solutions exist . at @xmath71 a second branch of solutions
is encountered . moving backwards along this second branch the value of the scalar field at the origin
@xmath70 continues to decrease , until it reaches zero at a critical value @xmath2 .    at @xmath2
the solutions change character and @xmath0-shells arise .
as @xmath69 decreases further , the inner radius @xmath72 increases along with the outer radius @xmath41 . in the limit @xmath73 the size of the shells diverges while the ratio @xmath74 tends to one .
thus there are no @xmath0-shells in a minkowski or anti - de sitter background .    likewise , when the @xmath0-ball solutions are continued to negative values of the cosmological constant , they change smoothly from the minkowski background solutions , as seen in fig .
as the scaled cosmological constant @xmath69 decreases from zero , the value of the scalar field at the origin @xmath70 increases along with the outer radius @xmath41 until a minimal value @xmath75 encountered , beyond which no such solutions exist .
our data indicate , that @xmath76 . for a derivation of this limit in @xmath56 dimensions
see appendix a.
to study the back reaction of the @xmath0-balls and @xmath0-shells on the space - time , we have solved the coupled system of equations for the metric and the scalar field numerically . in the following we first discuss the case of @xmath56 dimensions , and then turn to other dimensions .
let us start by briefly recalling the properties of the single family of asymptotically flat compact boson star solutions found in @xcite . at @xmath55
this family starts from the @xmath0-ball solution in the minkowski background . as seen in fig .
[ fig2 ] , with increasing @xmath30 the value of the scalar field at the origin @xmath70 decreases , until it reaches a finite minimum , and then increases strongly , while @xmath30 undergoes damped oscillations . in other physical quantities these damped oscillations with respect to @xmath30 lead to a spiral - like pattern , as seen for the outer radius @xmath41 or the mass parameter @xmath77 .
such a behaviour is typical for boson stars and neutron stars .    here
we have extended this family of compact boson star solutions to negative values of the coupling constant @xmath30 , as seen in fig .
the physical interpretation of the solutions with negative @xmath30 is that they represent compact solutions made from phantom scalar fields .
thus the negative sign of @xmath30 can be absorbed by the negative lagrangian of the phantom field .
ever since the relevance of dark energy for cosmology became apparent , such phantom fields are found ubiquitously in the literature .
moreover , phantom fields also allow for the formation of various types of wormholes .
let us now turn to compact de sitter boson stars by increasing the value of @xmath69 from zero .
we demonstrate the effect of a positive cosmological constant on the compact solutions in fig .
[ fig2 ] by exhibiting the physical properties of such families of solutions for several values of @xmath69 .
we see that for finite @xmath69 the minimum of @xmath70 reaches zero .
this signals the occurrence of boson shells .
in particular , for a given value of @xmath69 , @xmath70 reaches zero at a critical value @xmath78 . then for @xmath79 boson shells exist .
@xmath78 is exhibited in fig .
for very small positive @xmath69 , the critical value @xmath78 is negative . thus there
exist phantom boson shells in this case , which turn into ordinary boson shells , as @xmath30 increases beyond zero .
extrapolation to @xmath80 indicates , that indeed no shells exist in this limit .    for all values of @xmath69 the compact boson stars exhibit the characteristic spirals . in fig .
[ fig2 ] these are seen for the value of the outer radius @xmath41 , the value of the metric function at the outer radius @xmath81 and the value of the scaled mass @xmath77 .
interestingly , phantom type boson star solutions exist only for small values of @xmath69 .
since the mass @xmath77 has the same sign as @xmath30 , branches with negative mass exist only for small values of @xmath69 .    whereas there is an upper bound @xmath82 for the compact boson stars
, there is no such bound for the boson shells . however , with increasing @xmath30 their outer radius decreases , and tends to a finite limiting value . at the same time ,
the ratio of the inner and outer radii @xmath74 increases , and tends to one .
thus the shells become smaller and thinner , while their scaled mass @xmath77 grows . on the other hand , as @xmath30 is kept fixed while @xmath69 is decreased , the shells grow in size , while the ratio @xmath74 tends to one . here in the limit
@xmath83 the shell size diverges .      in the above subsection
we have constructed the domain of existence of compact boson stars and boson shells in terms of dimensionless quantities .
we can obtain physical solutions with dimensionful quantities by scaling these dimensionless solutions appropriately .
in the case of vanishing cosmological constant we have considered compact stars @xcite . in particular ,
we have shown , that when the mass of these boson stars is on the order of the solar mass then their radius is on the order of ten(s ) of kilometers , thus they can correspond in mass and size to neutron stars . moreover , spirals are also encountered for neutron stars , when they approach the black hole limit .    concerning their stability , we employ arguments from catastrophe theory @xcite . according to catastrophe theory ,
the stability changes only at turning points .
thus when starting from a stable configuration , the stability should change at the maximum of the mass .
therefore solutions inside the spiral should be unstable . for neutron stars or ordinary boson stars
this has been confirmed by a mode analysis .
since the value of the cosmological constant as obtained from cosmology is very small , @xmath84 m@xmath85 in metric units , its presence hardly affects the properties of boson stars that have masses on the order of the mass of the sun .
thus the results for boson stars correspond to those obtained before @xcite .
however , the presence of a positive cosmological constant , no matter how small , does allow for boson shells . but
those boson shells possess cosmological mass and length scales .
it would be interesting to see whether such thin boson shells can be associated with voids , i.e.  with vast regions of empty space surrounded by a shell of matter . to see the effect of @xmath14 on the compact boson stars
let us now consider very big scales . in fig .
[ fig4 ] we show our scaled results for the mass @xmath86 in units of @xmath87 for boson stars and boson shells versus the outer radius @xmath41 in units of gigaparsec . keeping @xmath14 at its physical value ,
we have translated it into a timescale @xmath88 via the relation ( [ scaling ] ) , employing units of gigayears .
shown in addition are the black hole horizon and the cosmological horizon of the corresponding schwarzschild - de sitter space - times .
note , that these two horizons coincide for the extremal configuration with the maximum value of the mass .
these schwarzschild - de sitter values form the boundary , within which all extended objects must remain . as seen in fig .
[ fig4 ] , the boson stars reside well within these bounds .
note , that since we consider only positive values of the gravitational coupling , i.e. , no phantom fields , the boson star curves associated with the smaller values of the parameter @xmath88 consist of two disconnected parts .
the boson shells , in contrast , exist until they reach this cosmological bound .
in particular , all boson shell curves extend precisely to the extremal value of the schwarzschild - de sitter curve , where the two horizons coincide .
when approaching this limiting configuration , the inner radius of the shells approaches the outer radius .
thus the ratio @xmath74 tends to the value one in this limit .
more massive shells can not exist .    converting the values in fig .
[ fig4 ] into numerical values we find that with @xmath89 the mass of the boson stars is on the order of @xmath90 kg and their radius is on the order of gpc .
these are sizes that are beyond those of galaxies and galaxy clusters .
if there were dark matter distributions on such large scales , our solutions would be able to model those .
choosing somewhat smaller values of @xmath88 , on the other hand , we would find sizes relevant for galaxies or galaxy clusters , so that these solutions could be considered to model the dark matter halo of galaxies or the dark matter in galaxy clusters , respectively .
note that in the limit of vanishing mass the radius @xmath91 becomes spurious since the non - scaled boson field vanishes identically .
as expected , compact boson stars exist also for negative values of the cosmological constant .
indeed , the asymptotically minkowski solutions can be smoothly extended to negative values of @xmath69 , thus yielding asymptotically ads boson stars . in fig .
[ fig5 ] we exhibit some of their physical properties versus the coupling constant @xmath30 , for several values of @xmath69 .
we note , that the domain of existence of these compact ads boson stars decreases with decreasing @xmath69 .
this suggests that there is a limiting minimal value for @xmath69 for compact boson stars .
analogously to the ds case , some of the physical properties of ads boson stars exhibit damped oscillations with respect to @xmath30 whereas other properties exhibit spirals .
moreover , as in the ds case , there are phantom boson stars , associated with negative values of @xmath30 .
however , we do not find ads boson shells . in the ads case , the scalar field @xmath92 never reaches the value zero at the origin , necessary for boson shells to arise .
we conclude , that the extra attraction associated with negative @xmath14 inhibits the formation of ads shells even stronger than in the asymptotically flat case , @xmath93 .
the existence of shells needs repulsion , that can be provided either by a positive cosmological constant , as seen in the previous subsection , or by the presence of electric charge @xcite .
here we consider the domain of solutions and their properties for various space - time dimensions .
we have made a complete study for dimensions @xmath67 , 5 and 10 .
since the dependence on @xmath1 is mostly rather smooth , we exhibit only a number of selected cases . as an example of a solution , we show an asymptotically de sitter boson shell solution in fig .
[ fig6 ] .
we start our discussion by considering boson stars and boson shells in @xmath95 dimensions .
some of their properties are exhibited in fig .
[ fig7 ] , and can be compared to those of fig .
[ fig2 ] .
a surprising feature is that in @xmath95 boson shells seem to exist in the asymptotically flat case , @xmath96 .
this follows from the additional ( almost vertical ) line present in fig .
[ fig7 ] for @xmath70 .
however , since this line is reaching zero for a very small negative value of @xmath30 , those shells are phantom shells .
the small branch of phantom shells is seen more clearly in the plot of @xmath97 in the inset . as @xmath83 , the size of these phantom shells diverges .    for finite values of @xmath69 , however , we obtain also ordinary ds boson shells .
the critical value @xmath78 of the transition between the boson stars and shells is seen in fig .
in particular , the resulting phantom ds boson shells continue to exist beyond @xmath55 , where they smoothly turn into ordinary boson shells . as in @xmath56 , at fixed @xmath69 with increasing @xmath30 the outer radius @xmath41 of these boson shells decreases , tending to a finite limiting value , while their ratio @xmath74 of inner and outer radius increases towards one .
the compact ds boson stars , on the other hand , exist only below a maximal value of @xmath30 , which is either given by the onset of the spiral or by the transition to ds boson shells . for a given @xmath69 ,
the domain of existence of ds boson stars in @xmath95 with respect to @xmath30 is smaller than in @xmath56 .
when going to higher dimensions , this trend continues .
in contrast , the domain of existence of phantom ds boson stars increases with increasing @xmath1
.    in @xmath67 dimensions gravity is non - dynamic .
therefore we may expect that boson stars may exhibit a different behavior in @xmath67 dimensions than in higher dimensions .
let us first inspect the dependence of some of the properties of the compact 3-dimensional objects on @xmath30 for several values of @xmath69 .
[ fig8 ] exhibits the scalar field @xmath92 at the origin , @xmath70 , the outer radius @xmath41 , and the scaled mass @xmath77 .
here we observe , that indeed some properties of compact boson stars are very different in @xmath67 dimensions .
first of all , we notice a maximal value of @xmath70 , that decreases with increasing @xmath69 .
this maximal value is encoutered for negative values of @xmath30 , thus these configurations correspond to phantom boson stars .    as this maximal value of @xmath70 two branches of solutions merge .
along one of these branches @xmath70 reaches zero .
thus boson shells emerge at the corresponding critical value @xmath98 .
with increasing @xmath69 the critical value @xmath98 increases ( at least for positive @xmath98 ) , analogously to other dimensions .
the second branch of boson star solutions , however , exhibits a different behaviour from the one observed before
. clearly , the damped oscillations of @xmath70 with @xmath30 are not present
. instead a monotonic decrease of @xmath70 with increasing @xmath30 is observed , and no maximal value of @xmath30 is encountered .
we further observe that for small @xmath69 , @xmath70 and the outer radius @xmath41 depend only weakly on @xmath69 .
the absence of a maximal value of @xmath30 suggests to consider the dependence of the compact boson stars on @xmath69 , choosing fixed large values of @xmath30 . as can be seen from figs .
[ fig8]e and f , for a given @xmath30 boson star solutions exist only up to a maximal value of @xmath69 , that increases with @xmath30 .
finally , we note that the solutions with @xmath68 are not asymptotically flat , since their mass parameter @xmath77 is finite .
this holds for boson stars and boson shells , alike .
let us turn finally to ads boson stars in @xmath1 dimensions .
some properties of asymptotically ads boson stars in @xmath95 are exhibited in fig .
[ fig9 ] . as in four dimensions , there are no ordinary ads boson shells in other than four dimensions .
however , there is a very small region of phantom ads boson shells .
this is seen by inspecting the critical curve @xmath78 for negative values of @xmath30 close to zero . here
the critical curve passes negative values of @xmath69 as well .
like the phantom shells for @xmath99 these phantom ads shells diverge in size as @xmath83 .
otherwise all the basic properties of these ads solutions in @xmath95 are similar to those of the @xmath56 ads solutions discussed above .
moreover , we observe only gradual changes with increasing @xmath1 .
in particular , the domain of existence of ads boson stars with respect to @xmath30 decreases with @xmath1 .
only the domain of phantom ads shells increases slightly , as seen in fig .
[ fig3 ] .
the lower dimensional case @xmath67 is special again , as can be seen in figs .
[ fig8]e and f , where the dependence of @xmath70 and @xmath77 on @xmath69 is shown for several values of @xmath30 .
interestingly , for a given @xmath30 there is no lower bound of @xmath69 encountered .
moreover , the mass becomes practically independent of @xmath30 as @xmath69 becomes sufficiently small .
when @xmath69 is fixed instead , while @xmath30 is varied , no spirals are encountered for the compact 3-dimensional boson stars , whereas spirals are present in all higher dimensions .
this was observed before for ordinary boson stars @xcite . restricting to positive @xmath30
, there is a maximal value of @xmath70 for a given @xmath69 .
however , when allowing for phantom fields , a minimal value of @xmath30 is encountered , beyond which @xmath70 increases without bound , while the mass decreases .
we have studied compact boson stars and shells obtained with a v - shaped interaction potential in @xmath10 dimensions .
the v - shaped potential confines the scalar field to a finite region , which can be ball - like or shell - like .    in the probe limit
, we have given the general analytical solution for @xmath0-balls in a minkowski background . here
no @xmath0-shells exist .
@xmath0-shells arise only beyond a critical value @xmath100 of the cosmological constant , which increases with the number of dimensions @xmath1 .
likewise , @xmath0-balls exist only above a minimal value of the cosmological constant , which seems to correspond to @xmath76 .
subsequently , we have taken the backreaction into account .
the resulting configurations correspond to compact boson stars and boson shells . by solving the coupled set of einstein - scalar field equations , we have obtained the full set of solutions , subject to minkowski , de sitter and anti - de sitter asymptotics for a number of space - time dimensions , ranging from 3 to 10 .    for any dimension @xmath101 there are compact boson stars with all three types of asymptotics .
but concerning their properties , we see a distinct behaviour in three dimensions , that is different from the common behaviour encountered in all higher dimensions . in four and higher dimensions ,
these boson stars exist in a finite intervall @xmath102 . in constrast ,
in three dimensions there is no upper bound on the value of @xmath30 .
also , all boson stars in four and higher dimensions exhibit a spiral - like dependence of the outer radius and the mass on the coupling constant @xmath30 . at the same time , the scalar field value @xmath103 exhibits damped oscillations . in constrast , in three dimensions the respective boson star properties do not exhibit such a spiral - like dependence or damped oscillations .    by exploring the parameter space
, we also find boson stars for negative values of @xmath30 .
these boson stars with negative @xmath30 correspond to phantom boson stars , since the negative sign of @xmath30 can be reinterpreted as a negative sign associated with the scalar field in the lagrangian , and thus with a phantom field .
ample motivation for the consideration of phantom fields is nowadays provided by cosmology . in constrast to boson stars , boson shells do not exist for minkowski asymptotics , if there is no additional force present , balancing the gravitational attraction .
consequently , there are no ordinary ads boson shells .
however , ds boson shells do exist . here
the positive cosmological constant provides the necessary repulsion . on the other hand ,
there exist small regions of phantom ads shells in more than four dimensions .    in four dimensions we have also considered astrophysical aspects of the compact boson stars and boson shells . while we can always adjust the parameters of the solutions to describe compact astrophysical objects with masses and sizes of neutron stars as discussed in @xcite , the influence of the cosmological constant on these objects is negligible , when the physical value of @xmath14 is taken .
the new feature is , however , that in addition to boson stars there exist also boson shells for positive values of @xmath14 .
we have then addressed the question , what the properties of such compact objects would be , if we set the scale by the physical value of @xmath14 .
interestingly , in this case the resulting sets of boson stars reach huge masses and sizes , that are more akin to structures on the largest scales of the universe . the boson shells , on the other hand , can grow in mass until they reach the limit , set by the extremal schwarzschild - de sitter solution , for which the event horizon and the cosmological horizon merge .    for negative @xmath14 one might be tempted to consider the ads / cft correspondence and try to interpret the solutions within this framework .
however , all of our solutions are compact .
the outer solutions are all given in terms of the schwarzschild - ads solutions .
thus the ads boundary does not feel anything of the solutions except for their mass .
consequently , for such compact solutions the concept of holography does not work .
indeed , lots of different compact objects may sit in the bulk , and if they have the same mass , the boundary does not notice a difference .    as our next step
we plan to include rotation @xcite .
rotating boson stars are known for non - compact configurations @xcite .
interestingly , their angular momentum @xmath104 is quantized in terms of their particle number @xmath0 , @xmath105 , where @xmath106 is an integer .
we expect , that the rotation of compact boson stars will lead to interesting new features .
moreover , there may be rotating boson shells in the presence of a cosmological constant .
it should also be interesting to construct interacting compact @xmath0-balls and @xmath0-shells for finite cosmological constant .
these should arise in the presence of several complex scalar fields @xcite .
we would like to thank eugen radu for helpful discussions .
we gratefully acknowledge support by the dfg , in particular , also within the dfg research training group 1620  models of gravity  .
here we consider the ads solutions in the probe limit . we introduce the scaled coordiante @xmath107 , where @xmath108 , and the scaled scalar field @xmath109 .
this yields for the ode of the function @xmath110 @xmath111 where @xmath112 , @xmath113 , and prime denotes the derivative with respect to @xmath114 .
now the general solution of the inhomogeneous ode can be written as @xmath119 where @xmath120 and @xmath121 are constants and @xmath122   \label{inhsol}\ ] ] is a special solution of the inhomogeneous ode .
let us assume that @xmath116 is regular at @xmath123 and @xmath124 . as a consequence the integral in eq .
( [ y2 ] ) diverges and @xmath125 is singular at @xmath123 . however , the special solution of the inhomogeneous ode is regular and vanishes at @xmath123 .
thus to obtain the regular solution with @xmath126 we set @xmath127 and @xmath128 .
next we consider the conditions for compact solutions , i.e. @xmath129 and @xmath130 at some @xmath131 .
@xmath132 = 0 \ , \label{bc1}\\ \psi'(\eta_o ) & = & y_1'(\eta_o ) + \frac{\ell^2}{\hat{\phi}(0 ) } \left[y_2'(\eta_o)\int_0^{\eta_o}\tau^m y_1(\tau)d\tau - y_1'(\eta_0)\int_0^{\eta_o}\tau^m y_2(\tau)d\tau\right ] = 0 \ . \label{bc2}\end{aligned}\ ] ] the linear superpositions @xmath133 and @xmath134 yield @xmath135\int_0^{\eta_o}\tau^m y_1(\tau)d\tau \ , \label{bc1a}\\ 0 & = & \left[(y_1'(\eta_o)y_2(\eta_0)-y_2'(\eta_o)y_1(\eta_o)\right ] \left(1-\frac{\ell^2}{\hat{\phi}(0)}\int_0^{\eta_o}\tau^m y_2(\tau)d\tau\right ) \ ,   \label{bc2a}\end{aligned}\ ] ] respectively , where we set @xmath136 .    since @xmath116 and @xmath125 are solutions of the homogeneous ode it follows that @xmath137 for some constant @xmath138 .
as a consequence @xmath139 \neq 0 $ ] and the conditions ( [ bc1a ] ) and ( [ bc2a ] ) reduce to @xmath140 the first equation determines the point @xmath131 and the second the value of @xmath70 .
now we are left with the problem to determine the range of @xmath141 for which solutions of eq .
( [ bc1b ] ) exist .
although the solutions of the homogeneous ode can be expressed in terms of hypergeometric functions , we did not succeed to determine the minimal values of @xmath141 in the general case , i.e. in all dimensions .
however , in @xmath56 dimensions the solution @xmath116 can be expressed in terms of trigonometric functions , @xmath142 \ , \label{homsol}\ ] ] which simplifies the problem considerably .
here we consider only @xmath143 .
the case @xmath144 needs special treatment .    we will show that eq .
( [ bc1b ] ) has only a solution @xmath131 if @xmath145 .
clearly , the integral in eq .
( [ bc1b ] ) can only vanish if the integrand possesses ( at least ) one zero . to analyze this condition
we introduce a new coordinate @xmath146 , with @xmath147 , and consider the function @xmath148 since @xmath149 does not change sign on the interval @xmath147 it is sufficient to consider the function @xmath150 now we are left with the question for which values of @xmath141 the function @xmath151 does not possess a zero .
let us first restrict to @xmath152 .
expanding the function @xmath151 for small values of @xmath153 we find @xmath154 . on the other hand evaluating @xmath151 at @xmath155 yields @xmath156 for @xmath157 .
thus , @xmath158 possesses at least one zero for @xmath157 .
consequently , we can restrict to @xmath159 .
we rewrite eq .
( [ eqf1 ] ) as @xmath160 \label{eqf2 } \\       & = & -\frac{\ell+1}{2}\left[b\cos a \sin a+ \sin(ba)\right ] \ , \label{eqf3}\end{aligned}\ ] ] where @xmath161 , @xmath162 are restricted to @xmath163 and @xmath164 .
it can be seen from fig .
[ fig11 ] that @xmath165 indeed does not a possess a zero in the region @xmath166 .
consequently there are no compact solutions with @xmath167 .
now we turn to the case @xmath168 . with @xmath169 and @xmath170
the function @xmath151 reads @xmath171 \ . \label{eqf4 } \\       & = &   \frac{1}{2\ell'}\left[(\ell'-1)\sin(z'(\ell'+1))+(\ell'+1 ) \sin(z'(\ell'-1))\right ] \label{eqf5}\end{aligned}\ ] ] note that now @xmath172 .
therefore , @xmath173 and @xmath174 imply @xmath175 and @xmath176 .
this shows that the function @xmath151 does not have a zero for @xmath168 .    to conclude , we have found a lower bound for @xmath141 , i.e. @xmath177 , below which no compact solution can exist .
however , this does not prove that solutions exist for @xmath178 , since the condition eq .
( [ bc1b ] ) is more restrictive than @xmath179 for some @xmath114 . indeed , for @xmath180 eq .
( [ bc1b ] ) yields an equation of the form @xmath181 , which has no solution for @xmath182 .
however , consider the case when the function @xmath158 possesses only one zero on the interval @xmath183 for some @xmath145 .
denote the zero by @xmath184 .
the condition ( [ bc1b ] ) can then be written as @xmath185 the first integral yields a finite negative contribution , since @xmath158 is bounded and negative on @xmath186 .
the second integral on the other hand assumes any positive value between zero and infinity , as @xmath187 ranges between @xmath184 and @xmath188 , since @xmath158 is bounded and positive on @xmath189 .
consequently , there exists a @xmath187 on @xmath189 for which both integrals in eq .
( [ intf ] ) cancel , implying that compact solutions exist for some @xmath190 .    as an example we computed the exact solution for @xmath191 .
we found @xmath192\right.\right .
\nonumber\\ & & \left.\left .             \
\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ -15(\eta^4 + 10\eta^2 - 15)\eta\log(1+\eta^2)\right)\right\ } \ .\end{aligned}\ ] ] for this solution the values @xmath193 and @xmath194 have been computed numerically .
remarkably , they coincide with the corresponding values of our numerically computed solution of the ode up to eight digits .
g.  aad _ et al .
_ [ atlas collaboration ] , phys .
b * 716 * , 1 ( 2012 ) [ arxiv:1207.7214 [ hep - ex ] ] .
s.  chatrchyan _ et al . _
[ cms collaboration ] , phys .
b * 716 * , 30 ( 2012 ) [ arxiv:1207.7235 [ hep - ex ] ] .
e.  w.  mielke and r.  scherzer , phys .
d * 24 * ( 1981 ) 2111 .
r.  friedberg , t.  d.  lee and a.  sirlin , phys .
d * 13 * , 2739 ( 1976 ) .
s.  r.  coleman , nucl .
b * 262 * , 263 ( 1985 ) [ erratum - ibid .
b * 269 * , 744 ( 1986 ) ] .
b.  hartmann , b.  kleihaus , j.  kunz and i.  schaffer , phys .
b * 714 * , 120 ( 2012 ) [ arxiv:1205.0899 [ gr - qc ] ] .
h.  arodz and j.  lis , phys .
d * 77 * , 107702 ( 2008 ) .
[ arxiv:0803.1566 [ hep - th ] ] .
h.  arodz and j.  lis , phys .
d * 79 * , 045002 ( 2009 ) .
[ arxiv:0812.3284 [ hep - th ] ] .
h.  arodz , j.  karkowski and z.  swierczynski , acta phys .
b * 43 * ( 2012 ) 79 [ arxiv:1201.2279 [ hep - th ] ] .
b.  kleihaus , j.  kunz , c.  lammerzahl and m.  list , phys .
b * 675 * , 102 ( 2009 ) . [ arxiv:0902.4799 [ gr - qc ] ] .
b.  kleihaus , j.  kunz , c.  lammerzahl and m.  list , phys .  rev .
d * 82 * , 104050 ( 2010 ) .
[ arxiv:1007.1630 [ gr - qc ] ] .
f.  e.  schunck and a.  r.  liddle , lect .
notes phys .
* 514 * , 285 ( 1998 ) [ arxiv:0811.3764 [ astro - ph ] ] .
j.  m.  maldacena , adv .
* 2 * , 231 ( 1998 ) [ hep - th/9711200 ] .
e.  witten , adv .
* 2 * , 253 ( 1998 ) [ hep - th/9802150 ] .
k.  sakamoto and k.  shiraishi , jhep * 9807 * , 015 ( 1998 ) [ gr - qc/9804067 ] .