We present a novel, relativistic accretion model for accretion onto a Schwarzschild black hole. This consists of a purely hydrodynamical mechanism in which, by breaking spherical symmetry, a radially accreting flow transitions into an inflow-outflow configuration. The spherical symmetry is broken by considering that the accreted material is more concentrated on an equatorial belt, leaving the polar regions relatively under-dense. What we have found is a flux-limited accretion regime in which, for a sufficiently large accretion rate, the incoming material chokes at a gravitational bottleneck and the excess flux is redirected by the density gradient as a bipolar outflow. The threshold value at which the accreting material chokes is of the order of the mass accretion rate found in the spherically symmetric case studied by Bondi and Michel. We describe the choked accretion mechanism first in terms of a general relativistic, analytic toy model based on the assumption of an ultrarelativistic stiff fluid. We then relax this approximation and, by means of numerical simulations, show that this mechanism can operate also for general polytropic fluids. Interestingly, the qualitative inflow-outflow morphology obtained appears as a generic result of the proposed symmetry break, across analytic and numeric results covering both the Newtonian and relativistic regimes. The qualitative change in the resulting steady state flow configuration appears even for a very small equatorial to polar density contrast ($\sim 0.1\,\%$) in the accretion profile. Finally, we discuss the applicability of this model as a jet-launching mechanism in different astrophysical settings.

Using a recent homogeneous sample of 40 high quality velocity dispersion profiles for Galactic globular clusters, we study the low gravitational acceleration regime relevant to the outskirts of these systems. We find that a simple empirical profile having a central Gaussian component and a constant large radius asymptote, $\sigma_{\infty}$, accurately describes the variety of observed velocity dispersion profiles. We use published population synthesis models, carefully tailored to each individual cluster, to estimate mass to light ratios from which total stellar masses, $M$, are inferred. We obtain a clear scaling, reminiscent of the galactic Tully-Fisher relation of $\sigma_{\infty}( km s^{-1})= 0.084^{+0.075}_{-0.040} (M/M_{\odot})^{0.3 \pm 0.051} $, which is interesting to compare to the deep MOND limit of $\sigma_{\infty} (km s^{-1})=0.2(M/M_{\odot})^{0.25}$. Under a Newtonian interpretation, our results constitute a further restriction on models where initial conditions are crafted to yield the outer flattening observed today. Within a modified gravity scheme, as the globular clusters studied are not isolated objects in the deep MOND regime, the results obtained point towards a modified gravity where the external field effect of MOND does not appear, or is much suppressed.

Steady state, spherically symmetric accretion flows are well understood in terms of the Bondi solution. Spherical symmetry however, is necessarily an idealized approximation to reality. Here we explore the consequences of deviations away from spherical symmetry, first through a simple analytic model to motivate the physical processes involved, and then through hydrodynamical, numerical simulations of an ideal fluid accreting onto a Newtonian gravitating object. Specifically, we consider axisymmetric, large-scale, small amplitude deviations in the density field such that the equatorial plane is over dense as compared to the polar regions. We find that the resulting polar density gradient dramatically alters the Bondi result and gives rise to steady state solutions presenting bipolar outflows. As the density contrast increases, more and more material is ejected from the system, attaining speeds larger than the local escape velocities for even modest density contrasts. Interestingly, interior to the outflow region, the flow tends locally towards the Bondi solution, with a resulting total mass accretion rate through the inner boundary $choking$ at a value very close to the corresponding Bondi one. Thus, the numerical experiments performed suggest the appearance of a maximum achievable accretion rate, with any extra material being ejected, even for very small departures from spherical symmetry.

We propose in this letter a relativistic coordinate independent interpretation for Milgrom's acceleration $a_{0}=1.2 \times 10^{-8} \hbox{cm/s}^{2}$ through a geometric constraint obtained from the product of the Kretschmann invariant scalar times the surface area of 2--spheres defined through suitable characteristic length scales for local and cosmic regimes, described by Schwarzschild and Friedman--Lema\^\i tre--Robertson--Walker (FLRW) geometries, respectively. By demanding consistency between these regimes we obtain an appealing expression for the empirical (so far unexplained) relation between the accelerations $a_0$ and $c H_0$. Imposing this covariant geometric criterion upon a FLRW model, yields a dynamical equation for the Hubble scalar whose solution matches, to a very high accuracy, the cosmic expansion rate of the $\Lambda$CDM concordance model fit for cosmic times close to the present epoch. We believe that this geometric interpretation of $a_0$ could provide relevant information for a deeper understanding of gravity

Gamma-ray bursts (GRBs) are spectacularly energetic events, with the potential to inform on the early universe and its evolution, once their redshifts are known. Unfortunately, determining redshifts is a painstaking procedure requiring detailed follow-up multi-wavelength observations often involving various astronomical facilities, which have to be rapidly pointed at these serendipitous events. Here we use Machine Learning algorithms to infer redshifts from a collection of observed temporal and spectral features of GRBs. We obtained a very high correlation coefficient ($0.96$) between the inferred and the observed redshifts, and a small dispersion (with a mean square error of $0.003$) in the test set. The addition of plateau afterglow parameters improves the predictions by $61.4\%$ compared to previous results. The GRB luminosity function and cumulative density rate evolutions, obtained from predicted and observed redshift are in excellent agreement indicating that GRBs are effective distance indicators and a reliable step for the cosmic distance ladder.

Under Newtonian dynamics, the relative motion of the components of a binary star should follow a Keplerian scaling with separation. Once orientation effects and a distribution of ellipticities are accounted for, dynamical evolution can be modelled to include the effects of Galactic tides and stellar mass perturbers, over the lifetime of the solar neighbourhood. This furnishes a prediction for the relative velocity between the components of a binary and their projected separation. Taking a carefully selected small sample of 81 solar neighbourhood wide binaries from the {\it Hipparcos} catalogue, we identify these same stars in the recent Gaia DR2, to test the prediction mentioned using the latest and most accurate astrometry available. The results are consistent with the Newtonian prediction for projected separations below 7000 AU, but inconsistent with it at larger separations, where accelerations are expected to be lower than the critical $a_{0}=1.2 \times 10^{-10} $ { m s$^{-2}$} value of MONDian gravity. This result challenges Newtonian gravity at low accelerations and shows clearly the appearance of gravitational anomalies of the type usually attributed to dark matter at galactic scales, now at much smaller stellar scales.

Within the dark matter paradigm, explaining observed orbital dynamics at galactic level through the inclusion of a dominant dark halo, implies also the necessary appearance of dynamical friction effects. Satellite galaxies, globular clusters and even stars orbiting within these galactic halos, will perturb the equilibrium orbits of dark matter particles encountered, to produce a resulting trailing wake of slightly enhanced dark matter density associated with any perturber in the halo. The principal effect of this gravitational interaction between an orbiting body and the dark matter particles composing it, is the appearance of a frictional drag force slowly removing energy and angular momentum from the perturber. Whilst this effect might be relevant to help bring about the actual merger of the components of interacting forming galaxies, at smaller stellar scales, it becomes negligible. However, the trailing wake will still be present. In this letter I show that the corresponding dark matter wake associated to the Sun, will constitute a small but resonant perturbation on solar system dynamics which can be ruled out, as current laser and radio ranging measurements are now over an order of magnitude more precise than the amplitude of the orbital perturbations which said wake implies. The absence of any such detection implies the nonexistence of the dynamical friction trailing wake on the sun, which in turn strongly disfavours dark matter as an explanation for the observed gravitational anomalies at galactic scales.

Understanding the observations of dynamical tracers and the trajectories of lensed photons at galactic scales within the context of General Relativity (GR), requires the introduction of a hypothetical dark matter dominant component. The onset of these gravitational anomalies, where the Schwarzschild solution no longer describes observations, closely corresponds to regions where accelerations drop below the characteristic $a_{0}$ acceleration of MOND, which occur at a well established mass-dependent radial distance, $R_{c}\propto (GM/a_{0})^{1/2}$. At cosmological scales, inferred dynamics are also inconsistent with GR and the observed distribution of mass. The current accelerated expansion rate requires the introduction of a hypothetical dark energy dominant component. We here show that for a Schwarzschild metric at galactic scales, the scalar curvature, K, multiplied by $(r^{4}/M)$ at the critical MOND transition radius, $r=R_{c}$, has an invariant value of $\kappa_{B}=K(r^{4}/M)=28Ga_{0}/c^{4}$. Further, assuming this condition holds for $r>R_{c}$, is consistent with the full spacetime which under GR corresponds to a dominant isothermal dark matter halo, to within observational precision at galactic level. For a FLRW metric, this same constant bounding curvature condition yields for a spatially flat spacetime a cosmic expansion history which agrees with the $\Lambda$CDM empirical fit for recent epochs, and which similarly tends asymptotically to a de Sitter solution. Thus, a simple covariant purely geometric condition identifies the low acceleration regime of observed gravitational anomalies, and can be used to guide the development of { extended} gravity theories at both galactic and cosmological scales.

The X-ray afterglow plateau emission observed in many Gamma-ray Bursts (GRBs) has been interpreted as either being fueled by fallback onto a newly formed black hole, or by the spin-down luminosity of an ultra-magnetized millisecond neutron star. If the latter model is assumed, GRB X-ray afterglow light curves can be analytically reproduced. We fit a sample of GRB X-ray plateaus, interestingly yielding a distribution in the magnetic field versus spin period (B-P) diagram consistent with $B\propto P^{7/6}$. This is expected from the well-established physics of the spin-up line minimum period for Galactic millisecond pulsars. The normalisation of the relation we obtain perfectly matches spin-up line predictions for the expected masses ($\sim 1 M_{\odot}$) and radii ($\sim 10 {\rm ~km}$) of newly born magnetars, and mass accretion rates consistent with GRB expectations of $10^{-4} M_{\odot}/{\rm s} <\dot{M}< 10^{-1} M_{\odot}/{\rm s}$. Short GRBs with extended emission (SEE) appear towards the high period end of the distribution, while the long GRBs (LGRBs) towards the short period end. This result is consistent with spin-up limit expectations where the total accreted mass determines the position of the neutron star in the B-P diagram. The P-B distribution for LGRBs and SEE are statistically different, further supporting the idea that the fundamental plane relation \citep{dainotti16c,Dainotti2017} is a powerful discriminant among those populations. Our conclusions are robust against suppositions regarding the GRB collimation angle and magnetar breaking index, which shifts the resulting magnetar properties parallel to the spin-up line, and strongly support a magnetar origin for GRBs presenting X-ray plateaus.

Pressure supported systems modeled under MONDian extended gravity are expected to show an outer flattening in their velocity dispersion profiles. A characteristic scaling between the amplitude of the asymptotic velocity dispersion and the radius at which the flattening occurs is also expected. By comprehensively analyzing the dynamical behavior of $\sim$300 extremely low rotating elliptical galaxies from the MaNGA survey, we show this type of pressure supported systems to be consistent with MONDian expectations, for a range of central velocity dispersion values of $60 km/s < \sigma_{central}< 280 km/s $ and asymptotic velocity dispersion values of $28km/s < \sigma_{\infty}<250 km/s$. We find that a universal velocity dispersion profile accurately describes the studied systems; the predicted kinematics of extended gravity are verified for all well observed galaxies.

The following is a short comment on the recent claim made by van Dokkum and collaborators about the existence of a low surface brightness galaxy, NGC1052-DF2, not containing dark matter. A discovery used by the authors to both reject proposals of a failure of Newtonian dynamics in the low acceleration regime (e.g., MOND), and prove the dark matter hypothesis is correct. It is shown here that the claim in untenable.

Long gamma-ray bursts (GRBs) with a plateau phase in their X-ray afterglows obeys a three-dimensional (3D) relation (Dainotti et al. 2016), between the rest-frame time at the end of the plateau, Ta, its corresponding X-ray luminosity, La, and the peak luminosity in the prompt emission, Lpeak, an extension of the two-dimensional Dainotti relation. This 3D relation identifies a GRB fundamental plane whose existence we here confirm. We extend the original analysis with X-ray data untill July 2016 gathering 183 Swift GRBs with afterglow plateaus and known redshifts. We compare several GRB categories, such as shorts with extended emission, SEE, X-ray Flashes, GRBs associated with SNe, a sample of only long-duration GRBs (132), selected from the total sample by excluding GRBs of the previous categories, and the gold sample, composed by GRBs with light curves with good data coverage and relatively flat plateaus. The relation planes for each of these categories are not statistically different from the gold fundamental plane, with the exception of the SEE, which are hence identified as a physically distinct class of objects. The gold fundamental plane has an intrinsic scatter smaller than any plane derived from the other sample categories. Thus, the distance of any particular GRB category from this plane becomes a key parameter. Additionally, we tested this 3D relations by using GRBs observed at high energy, namely the peak luminosity values derived by the Fermi-Gamma Ray Burst Monitor (GBM). The 3D relation is also confirmed for GRBs observed by the GBM, thus showing its independence from the energy range. Furthermore, we computed the several category planes with T*a as a dependent parameter obtaining for each category smaller intrinsic scatters (reaching a reduction of $24\%$ for the long GRBs). The fundamental plane is independent from several prompt and afterglow parameters.

In some instances, e.g. near phase transitions, thermodynamic fluctuations become macroscopically relevant, and relative amplitudes grow far above the standard $N^{-1/2}$ scale, with $N$ the number of particles. Such large fluctuations are characterised by a scale invariant Gaussian power spectrum. In this letter I show that the abrupt drop in the baryonic sound speed across recombination leads to conditions resulting in such large thermodynamic Gaussian fluctuations in the ionisation fraction of the baryons. Under pressure equilibrium, this will result in a mechanism for generating scale invariant density and temperature fluctuations in the baryonic component, inherent to the thermodynamics of the baryons themselves. Within a $\Lambda$CDM framework, this extra random fluctuation source leads to a decoupling of the inflationary relic small wave number spectrum and the amplitude of the Gaussian random fluctuations at frequencies higher than the first acoustic peak, an effect which could explain the mismatch between cosmic microwave background (CMB) inferences and local kinetic determinations of the Hubble constant. Within modified gravity theories in absence of dark matter, the mechanism proposed serves as a source for random Gaussian density fluctuations in the acoustic peak region.

Under Newtonian gravity total masses for dSph galaxies will scale as $M_{T} \propto R_{e} \sigma^{2}$, with $R_{e}$ the effective radius and $\sigma$ their velocity dispersion. When both of the above quantities are available, the resulting masses are compared to observed stellar luminosities to derive Newtonian mass to light ratios, given a physically motivated proportionality constant in the above expression. For local dSphs and the growing sample of ultrafaint such systems, the above results in the largest mass to light ratios of any galactic systems known, with values in the hundreds and even thousands being common. The standard interpretation is for a dominant presence of an as yet undetected dark matter component. If however, reality is closer to a MONDian theory at the extremely low accelerations relevant to such systems, $\sigma$ will scale with { stellar mass} $M_{*}^{1/4}$. This yields an expression for the mass to light ratio which will be obtained under Newtonian assumptions of $(M/L)_{N}=120 R_{e}(\Upsilon_{*}/L)^{1/2}$. Here we compare $(M/L)_{N}$ values from this expression to Newtonian inferences for this ratios for the actual $(R_{e}, \sigma, L)$ observed values for a sample of recently observed ultrafaint dSphs, obtaining good agreement. Then, for systems where no $\sigma$ values have been reported, we give predictions for the $(M/L)_{N}$ values which under a MONDian scheme are expected once kinematical observations become available. For the recently studied Dragonfly 44 { and Crater II systems}, reported $(M/L)_{N}$ values are also in good agreement with MONDian expectations.

Under a MONDian view, astrophysical systems are expected to follow Newtonian dynamics whenever the local acceleration is above the critical $a_{0}=1.2 \times 10^{-10} m s^{-2}$, and enter a modified regime for accelerations below this critical value. Indeed, the dark matter phenomenology on galactic and subgalactic scales appears always, and only, at low accelerations. It is standard to find the $a<a_{0}$ regime towards the low density outskirts of astronomical systems, where under a Newtonian interpretation, dark matter becomes conspicuous. Thus, it is standard to find, and to think, of the dense central regions of observed systems as purely Newtonian. However, under spherical symmetry in the MONDian as in the Newtonian case, the local acceleration will tend to zero as one approaches the very centre of a mass distribution. It is clear that for spherically symmetric systems, an inner $a<a_{0}$ region will necessarily appear interior to a critical radius which will depend on the details of the density profile in question. Here we calculate analytically such a critical radius for a constant density core, and numerically for a cored isothermal profile. Under a Newtonian interpretation, such a central MONDian region will be interpreted as extra mass, analogous to the controversial black holes sometimes inferred to lie at the centres of globular clusters, despite an absence of nuclear activity detected to date. We calculate this effect and give predictions for the "central black hole" mass to be expected under Newtonian interpretations of low density Galactic globular clusters.

For any MONDian extended theory of gravity where the rotation curves of spiral galaxies are explained through a change in physics rather than the hypothesis of dark matter, a generic dynamical behavior is expected for pressure supported systems: an outer flattening of the velocity dispersion profile occurring at a characteristic radius, where both the amplitude of this flat velocity dispersion and the radius at which it appears are predicted to show distinct scalings with the total mass of the system. By carefully analyzing the dynamics of globular clusters and elliptical galaxies, we are able to significantly extend the astronomical diversity of objects in which MONDian gravity has been tested, from spiral galaxies, to the much larger mass range covered by pressure supported systems. We show that a universal projected velocity dispersion profile accurately describes various classes of pressure supported systems, and further, that the expectations of extended gravity are met, across seven orders of magnitude in mass. These observed scalings are not expected under dark matter cosmology, and would require particular explanations tuned at the scales of each distinct astrophysical system.

We discuss the phenomenology of gravitational lensing in the purely metric $f\left(\chi\right)$ gravity, an $f(R)$ gravity where the action of the gravitational field depends on the source mass. We focus on the strong lensing regime in galaxy-galaxy lens systems and in clusters of galaxies. Using an approximate metric solution accurate to second order of the velocity field $v/c$, we show how, in the $f\left(\chi\right)=\chi^{3/2}$ gravity, the same light deflection can be produced by point-like lenses with masses smaller than in General Relativity; this mass difference increases with increasing impact parameter and decreasing lens mass. However, for sufficiently massive point-like lenses and small impact parameters, $f\left(\chi\right)=\chi^{3/2}$ and GR yield indistinguishable light deflection angles: this regime occurs both in observed galaxy-galaxy lens systems and in the central regions of galaxy clusters. In the former systems, the GR and $f\left(\chi\right)$ masses are compatible with the mass of standard stellar populations and little or no dark matter, whereas, on the scales of the core of galaxy clusters, the presence of substantial dark matter is required both in General Relativity, and in our approximate $f\left(\chi\right)=\chi^{3/2}$ point-like lens solution. We thus conclude that our approximate metric solution of $f\left(\chi\right)=\chi^{3/2}$ is unable to describe the observed phenomenology of the strong lensing regime without the aid of dark matter.

NGC 288 is a diffuse Galactic globular cluster, it is remarkable in that its low density results in internal accelerations being below the critical MOND $a_{0}$ acceleration throughout. This makes it an ideal testing ground for MONDian gravity, as the details of the largely unknown transition function between the Newtonian and modified regimes become unimportant. Further, exact analytical solutions exist for isothermal spherical equilibrium structures in MOND, allowing for arbitrary values of the anisotropy parameter, $\beta$. In this paper we use observations of the velocity dispersion profile of NGC 288, which is in fact isothermal, as dynamical constraints on MONDian models for this cluster, where the remaining free parameters are adjusted to fit the observed surface brightness profile. We find the optimal fit requires $\beta =0$, an isotropic solution with a total mass of $3.5 \pm 1.1 \times 10^{4} M_{\odot}$.

Dwarf spheroidal galaxies are the smallest known stellar systems where under Newtonian interpretations, a significant amount of dark matter is required to explain observed kinematics. In fact, they are in this sense the most heavily dark matter dominated objects known. That, plus the increasingly small sizes of the newly discovered ultra faint dwarfs, puts these systems in the regime where dynamical friction on individual stars starts to become relevant. We calculate the dynamical friction timescales for pressure supported isotropic spherical dark matter dominated stellar systems, yielding $\tau_{DF} =0.93 (r_{h}/10 pc)^{2} (\sigma/ kms^{-1}) Gyr$, { where $r_{h}$ is the half-light radius}. For a stellar velocity dispersion value of $3 km/s$, as typical for the smallest of the recently detected ultra faint dwarf spheroidals, dynamical friction timescales becomes smaller than the $10 Gyr$ typical of the stellar ages for these systems, for $r_{h}<19 pc$. Thus, this becomes a theoretical lower limit below which dark matter dominated stellar systems become unstable to dynamical friction. We present a comparison with structural parameters of the smallest ultra faint dwarf spheroidals known, showing that these are already close to the stability limit derived, any future detection of yet smaller such systems would be inconsistent with a particle dark matter hypothesis.

A class of long Gamma-Ray Bursts (GRBs) presenting light curves with an extended plateau phase in their X-ray afterglows obeys a correlation between the rest frame end time of the plateau, $T_a$, and its corresponding X-ray luminosity, $L_{a}$, Dainotti et al. (2008). In this work we perform an analysis of a total sample of 176 {\it Swift} GRBs with known redshifts, exhibiting afterglow plateaus. By adding a third parameter, that is the peak luminosity in the prompt emission, $L_{peak}$, we discover the existence of a new three parameter correlation. The scatter of data about this plane becomes smaller when a class-specific GRB sample is defined. This sample of 122 GRBs is selected from the total sample by excluding GRBs with associated Supernovae (SNe), X-ray flashes and short GRBs with extended emission. {\bf With this sample the three parameter correlation identifies a GRB `fundamental plane'}. Moreover, we further limit our analysis to GRBs with lightcurves having good data coverage and almost flat plateaus, 40 GRBs forming our `gold sample'. The intrinsic scatter, $\sigma_{int}=0.27 \pm 0.04$, for the three-parameter correlation for this last subclass is more than twice smaller than the value for the $L_{a}-T_a$ one, making this the tightest three parameter correlation involving the afterglow plateau phase. Finally, we also show that a slightly less tight correlation is present between $L_{peak}$ and a proxy for the total energy emitted during the plateau phase, $L_a T_a$, confirming the existence of an energy scaling between prompt and afterglow phases.

Removing outbursts from multiwavelength light curves of the blazar Mrk~421, we construct outburstless time series for this system. A model-independent power spectrum light curve analysis in the optical, hard X-ray and gamma-rays for this outburstless state and also the full light-curves, show clear evidence for a periodicity of ~ 310 days across all wavelengths studied. A subsequent full maximum likelihood analysis fitting an eclipse model confirms a periodicity of 310 \pm 1 days. The power spectrum of the signal in the outburstless state of the source does not follow a flicker noise behaviour and so, the system producing it is not self-organised. This and the fact that the periodicity is better defined in the outburstless state, strongly suggests that it is not produced by any internal physical processes associated to the central engine. The simplest physical mechanism to which this periodicity could be ascribed is a dynamical effect produced by an orbiting supermassive black hole companion eclipsing the central engine. Interestingly, the optimal eclipse model infers a brightness enhancement of (136.4 \pm 20 )%, suggesting an eclipse resulting in a gravitational lensing brightening. Consisting with this interpretation, the eclipse occurs for only ( 9.7 \pm 0.2)% of the orbital period.

The following is a comment on the recent letter by Iocco et al. (2015, arXiv:1502.03821) where the authors claim to have found "...convincing proof of the existence of dark matter...". The letter in question presents a compilation of recent rotation curve observations for the Milky Way, together with Newtonian rotation curve estimates based on recent baryonic matter distribution measurements. A mismatch between the former and the latter is then presented as "evidence for dark matter". Here we show that the reported discrepancy is the well known gravitational anomaly which consistently appears when dynamical accelerations approach the critical Milgrom acceleration a_0 = 1.2 \times 10^{-10} m / s^2. Further, using a simple modified gravity force law, the baryonic models presented in Iocco et al. (2015), yield dynamics consistent with the observed rotation values.

Toomre's $Q$ stability parameter has long been shown through various theoretical arguments and numerical simulations, to be the principal determinant of stability against self-gravity in a galactic disk, under classical gravity. Comparison with observations however, has not always confirmed the condition of $Q=Q_{crit}$ to be well correlated with various critical galactic radii. In this paper we derive the analogous critical parameter, $Q_{M}$, under MONDian gravity. The result is a modification by a factor of $(\sigma \Omega)/a_{0}$, $Q_{M}=(\sigma \Omega)^{2}/(a_{0}G \Sigma)$, where $a_{0}$ is the critical acceleration scale of MOND. We then show through a direct comparison to a homogeneous sample of observed disk galaxies with measured brightness profiles, rotation velocity curves and internal velocity dispersion profiles, that the critical radii at which brightness profiles dip below the exponential fit, are much more accurately predicted by $Q_{M}$ than by the $Q$ of classical gravity. This provides a new and completely independent argument supporting the reality of a change in the form of gravity on reaching the low acceleration regime.

We study the dark energy equation of state as a function of redshift in a non-parametric way, without imposing any {\it a priori} $w(z)$ (ratio of pressure over energy density) functional form. As a check of the method, we test our scheme through the use of synthetic data sets produced from different input cosmological models which have the same relative errors and redshift distribution as the real data. Using the luminosity-time $L_{X}-T_{a}$ correlation for GRB X-ray afterglows (the Dainotti et al. correlation), we are able to utilize GRB sample from the {\it Swift} satellite as probes of the expansion history of the Universe out to $z \approx 10$. Within the assumption of a flat FLRW universe and combining SNeIa data with BAO constraints, the resulting maximum likelihood solutions are close to a constant $w=-1$. If one imposes the restriction of a constant $w$, we obtain $w=-0.99 \pm 0.06$ (consistent with a cosmological constant) with the present day Hubble constant as $H_{0}=70.0 \pm 0.6$ ${\rm km} \, {\rm s}^{-1} {\rm Mpc}^{-1}$ and density parameter as $\Omega_{\Lambda 0}=0.723 \pm 0.025$, while non-parametric $w(z)$ solutions give us a probability map which is centred at $H_{0}=70.04 \pm 1$ ${\rm km} \, {\rm s}^{-1} {\rm Mpc}^{-1}$ and $\Omega_{\Lambda 0}=0.724 \pm 0.03$. Our chosen GRB data sample with full correlation matrix allows us to estimate the amount, as well as quality (errors) of data, needed to constrain $w(z)$ in the redshift range extending an order of magnitude in beyond the farthest SNeIa measured.

Recent observations for three types of astrophysical systems severely challenge the GR plus dark matter scenario, showing a phenomenology which is what modified gravity theories predict. Stellar kinematics in the outskirts of globular clusters show the appearance of MOND type dynamics on crossing the $a_{0}$ threshold. Analysis shows a ``Tully-Fisher'' relation in these systems, a scaling of dispersion velocities with the fourth root of their masses. Secondly, an anomaly has been found at the unexpected scales of wide binaries in the solar neighbourhood. Binary orbital velocities cease to fall along Keplerian expectations, and settle at a constant value, exactly on crossing the $a_{0}$ threshold. Finally, the inferred infall velocity of the bullet cluster is inconsistent with the standard cosmological scenario, where much smaller limit encounter velocities appear. This stems from the escape velocity limit present in standard gravity; the ``bullet'' should not hit the ``target'' at more than the escape velocity of the joint system, as it very clearly did. These results are consistent with extended gravity, but would require rather contrived explanations under GR, each. Thus, observations now put us in a situation where modifications to gravity at low acceleration scales cease to be a matter of choice, to now become inevitable.

Whenever in a classical accretion disk the thin disk approximation fails interior to a certain radius, a transition from Keplerian to radial infalling trajectories should occur. We show that this transition is actually expected to occur interior to a certain critical radius, provided surface density profiles are steeper than $\Sigma(R) \propto R^{-1/2}$, and further, that it probably corresponds to the observationally inferred phenomena of thick hot walls internally limiting the extent of many stellar accretion disks. Infalling trajectories will lead to the convergent focusing and concentration of matter towards the very central regions, most of which will simply be swallowed by the central object. We show through a perturbative hydrodynamical analysis, that this will naturally develop a well collimated pair of polar jets. A first analytic treatment of the problem described is given, proving the feasibility of purely hydrodynamical mechanisms of astrophysical jet generation. The purely hydrodynamic jet generation mechanism explored here complements existing ideas focusing on the role of magnetic fields, which probably account for the large-scale collimation and stability of jets.

The inability of primordial baryonic density fluctuations, as observed in the cosmic microwave background (CMB), to grow into the present day astronomical structures is well established, under Newtonian and Einsteinian gravity. It is hence customary to assume the existence of an underlying dark matter component with density fluctuations, $\Delta(M)$, having amplitudes much larger than what CMB observations imply for the baryons. This is in fact one of the recurrent arguments used in support of the dark matter hypothesis. In this letter we prove that the same extended theory of gravity which has been recently shown to accurately reproduce gravitational lensing observations, in absence of any dark matter, and which in the low velocity regime converges to a MONDian force law, implies a sufficiently amplified self-gravity to allow purely baryonic fluctuations with amplitudes in accordance with CMB constraints to naturally grow into the $z=0$ astrophysical structures detected. The linear structure formation scenario which emerges closely resembles the standard concordance cosmology one, as abundantly calibrated over the last decade to match multiple observational constraints at various redshifts. However, in contrast with what occurs in the concordance cosmology, this follows not from a critical dependence on initial conditions and the fine tuning of model parameters, but from the rapid convergence of highly arbitrary initial conditions onto a well defined $\Delta(M,z)$ attractor solution.

In this article we perform a second order perturbation analysis of the gravitational metric theory of gravity $ f(\chi) = \chi^{3/2} $ developed by Bernal et al. (2011). We show that the theory accounts in detail for two observational facts: (1) the phenomenology of flattened rotation curves associated to the Tully-Fisher relation observed in spiral galaxies, and (2) the details of observations of gravitational lensing in galaxies and groups of galaxies, without the need of any dark matter. We show how all dynamical observations on flat rotation curves and gravitational lensing can be synthesised in terms of the empirically required metric coefficients of any metric theory of gravity. We construct the corresponding metric components for the theory presented at second order in perturbation, which are shown to be perfectly compatible with the empirically derived ones. It is also shown that under the theory being presented, in order to obtain a complete full agreement with the observational results, a specific signature of Riemann's tensor has to be chosen. This signature corresponds to the one most widely used nowadays in relativity theory. Also, a computational program, the MEXICAS (Metric EXtended-gravity Incorporated through a Computer Algebraic System) code, developed for its usage in the Computer Algebraic System (CAS) Maxima for working out perturbations on a metric theory of gravity, is presented and made publicly available.

That the stellar halo of the Milky Way has a density profile which to first approximation satisfies $\rho \propto r^{-3}$ has been known for a long time. More recently, it has become clear that M31 also has such an extended stellar halo, which approximately follows the same radial scaling. Studies of distant galaxies have revealed the same phenomenology. Also, we now know that the density profiles of the globular cluster systems of our Galaxy and Andromeda to first approximation follow $\rho \propto r^{-3}$, $\Sigma \propto R^{-2}$ in projection. Recently, diffuse populations of stars have been detected spherically surrounding a number of Galactic globular clusters, extending much beyond the Newtonian tidal radii, often without showing any evidence of tidal features. Within the standard Newtonian and GR scenario, numerous and diverse particular explanations have been suggested, individually tailored to each of the different classes of systems described above. Here we show that in a MONDian gravity scenario, any isothermal tenuous halo of tracer particles forming a small perturbation surrounding a spherically symmetric mass distribution, will have an equilibrium configuration which to first approximation satisfies a $\rho \propto r^{-3}$ scaling.

Elliptical galaxies are systems where dark matter is usually less necessary to explain observed dynamics than in the case of spiral galaxies, however there are some instances where Newtonian gravity and the observable mass are insufficient to explain their observed structure and kinematics. Such is the case of NGC 4649, a massive elliptical galaxy in the Virgo cluster for which recent studies report a high fraction of dark matter, 0.78 at $4R_e$. However this galaxy has been studied within the MOND hypothesis, where a good agreement with the observed values of velocity dispersion is found. Using a MONDian gravity force law, here we model this galaxy as a self-consistent gravitational equilibrium dynamical system. This force law reproduces the MOND phenomenology in the $a<a_{0}$ regime, and reduces to the Newtonian case when $a>a_{0}$. Within the MONDian $a<a_{0}$ scales, centrifugal equilibrium or dispersion velocities become independent of radius, and show a direct proportionality to the fourth root of the total baryonic mass, $V^{4}\propto(M G a_{0})$. We find that the recent detailed observations of the surface brightness profile and the velocity dispersion profile for this galaxy are consistent with the phenomenology expected in MONDian theories of modified gravity, without the need of invoking the presence of any hypothetical dark matter.

Over the past couple of years, a number of observational studies have confirmed the flattening of the radial velocity dispersion profiles for stars in various nearby globular clusters. As the projected radial coordinate is increased, a radius appears beyond which, the measured velocity dispersion ceases to drop and settles at a fixed value, $\sigma_{\infty}$. Under Newtonian gravity, this is explained by invoking tidal heating from the overall Milky Way potential on the outer, more loosely bound stars, of the globular clusters in question. From the point of view of modified gravity theories, such an outer flattening is expected on crossing the critical acceleration threshold $a_{0}$, beyond which, a transition to MONDian dynamics is expected, were equilibrium velocities cease to be a function of distance. In this paper we attempt to sort out between the above competing explanations, by looking at their plausibility in terms of an strictly empirical approach. We determine Newtonian tidal radii using masses accurately calculated through stellar population modelling, and hence independent of any dynamical assumptions, distances, size and orbital determinations for a sample of 16 globular clusters. We show that their Newtonian tidal radii at perigalacticon are generally larger that the radii at which the flattening in the velocity dispersion profiles occurs, by large factors of 4, on average. While this point makes the Newtonian tidal explanation suspect, it is found that the radii at which the flattening is observed on average correlate with the radii where the $a_{0}$ threshold is crossed, and that $\sigma_{\infty}$ values scale with the fourth root of the total masses, all features predicted under modified gravity theories.

We study the evolution of the galactic spin using data of high redshift galaxies in the fields of the Great Observatories Origins Deep Survey (GOODS). Through simple dynamical considerations we estimate the spin for the disc galaxies in our sample and find that its distribution is consistent with that found for nearby galaxies. Defining a dimensionless angular momentum parameter for the disc component of the galaxies ($\lambda_{d}$), we do not find signs of evolution in the redshift range $0.4 \leq z \leq 1.2$. We find that the mass and environmental dependence of the spin of our high redshift galaxies are similar to that of low-$z$ galaxies; showing a strong dependence on mass, in the sense that low-mass systems present higher $\lambda_{d}$ values than high-mass galaxies, with no significant dependence on the environmental density. These results lead us to conclude that, although individual disc galaxies might occasionally suffer from strong evolution, they evolve in such a way that the overall spin distribution of the galactic population remains constant from $z\sim1$ to the present epoch.

Assuming Newton's gravity and GR to be valid at all scales leads to the dark matter hypothesis as a requirement demanded by the observed dynamics and measured baryonic content at galactic and extragalactic scales. Alternatively, modified gravity scenarios where a change of regime appears at acceleration scales $a<a_{0}$ have been proposed. This modified regime at $a<a_{0}$ will generically be characterised by equilibrium velocities which become independent of distance. Here we identify a critical test in this debate and we propose its application to samples of wide binary stars. Since for $1 M_{\odot}$ systems the acceleration drops below $a_{0}$ at scales of around 7000 AU, a statistical survey of wide binaries with relative velocities and separations reaching $10^{4}$ AU and beyond should prove useful to the above debate. We apply the proposed test to the best currently available data. Results show a constant upper limit to the relative velocities in wide binaries which is independent of separation for over three orders of magnitude, in analogy with galactic flat rotation curves in the same $a<a_{0}$ acceleration regime. Our results are suggestive of a breakdown of Kepler's third law beyond $a \approx a_{0}$ scales, in accordance with generic predictions of modified gravity theories designed not to require any dark matter at galactic scales and beyond.

In modified theories of gravity including a critical acceleration scale, $a_{0}$, a critical length scale, $r_{M}=(GM/a_{0})^{1/2}$, will naturally arise, with the transition from the Newtonian to the dark matter mimicking regime occurring for systems larger than $r_{M}$. This adds a second critical scale to gravity, in addition to the one introduced by the criterion $v < c$ of the Shwarzschild radius, $r_{S}=2GM/c^{2}$. The distinct dependencies of the two above length scales give rise to non-trivial phenomenology in the (mass, length) plane for astrophysical structures, which we explore here. Surprisingly, extrapolation to atomic scales suggests gravity should be at the dark matter mimicking regime there.

In the context of theories of gravity modified to account for the observed dynamics of galactic systems without the need to invoke the existence of dark matter, a prediction often appears regarding low acceleration systems: wherever $a$ falls below $a_{0}$ one should expect a transition from the classical to the modified gravity regime.This modified gravity regime will be characterised by equilibrium velocities which become independent of distance, and which scale with the fourth root of the total baryonic mass, $V^{4} \propto M$. The two above conditions are the well known flat rotation curves and Tully-Fisher relations of the galactic regime. Recently however, a similar phenomenology has been hinted at, at the outskirts of Galactic globular clusters, precisely in the region where $a<a_{0}$. Radial profiles of the projected velocity dispersion have been observed to stop decreasing along Keplerian expectations, and to level off at constant values beyond the radii where $a<a_{0}$. We have constructed gravitational equilibrium dynamical models for a number of globular clusters for which the above gravitational anomaly has been reported, using a modified Newtonian force law which yields equilibrium velocities equivalent to MOND. We find models can be easily constructed having an inner Newtonian region and an outer modified gravity regime, which reproduce all observational constraints, surface brightness profiles, total masses and line of sight velocity dispersion profiles. Through the use of detailed single stellar population models tuned individually to each of the globular clusters in question, we derive estimates of the total masses for these systems. Interestingly, we find that the asymptotic values of the velocity dispersion profiles are consistent with scaling with the fourth root of the total masses, as expected under modified gravity scenarios.

Using dimensional analysis techniques we present an extension of Newton's gravitational theory built under the assumption that Milgrom's acceleration constant is a fundamental quantity of nature. The gravitational force converges to Newton's gravity and to a MOND-like description in two different mass and length regimes. It is shown that a modification on the force sector (and not in the dynamical one as MOND does) is more convenient and can reproduce and predict different phenomena usually ascribed to dark matter at the non-relativistic level.

Modified gravity scenarios where a change of regime appears at acceleration scales $a<a_{0}$ have been proposed. Since for $1 M_{\odot}$ systems the acceleration drops below $a_{0}$ at scales of around 7000 AU, a statistical survey of wide binaries with relative velocities and separations reaching $10^{4}$ AU and beyond should prove useful to the above debate. We apply the proposed test to the best currently available data. Results show a constant upper limit to the relative velocities in wide binaries which is independent of separation for over three orders of magnitude, in analogy with galactic flat rotation curves in the same $a<a_{0}$ acceleration regime. Our results are suggestive of a breakdown of Kepler's third law beyond $a \approx a_{0}$ scales, in accordance with generic predictions of modified gravity theories designed not to require any dark matter at galactic scales and beyond.

As the total mass of stellar populations of fixed IMF, age and metallicity tends to infinity, resulting intrinsic $M/L$ values quickly converge to fixed numbers associated to the particulars of the stars in question. When going to very small stellar populations, the stochastic sampling of the IMF will yield an inherent spread for the $M/L$ values, which will become probabilistic quantities. For the recently discovered ultra-faint dwarf spheroidal galaxies, with total luminosities dropping below 10^3 L_V/L_o, it is important to asses the amplitude of the probabilistic spread in inherent $M/L$ values. The total baryonic masses of these systems are usually estimated from their luminosities, assuming fixed, deterministic $M/L$ values, suitable for the infinite population limit. These total baryonic masses are crucial for the testing of structure formation scenarios, as the local ultra-faint dwarf spheroidals represent the most extreme galactic scales known. Also subject to reliable $M/L$ values, is the use of these systems as possible discriminants between dark matter and modified gravity theories. By simulating large collections of stellar populations, each consisting of a particular collection of individual stars, we compute statistical distributions for the resulting $M/L$ values. For total numbers of stars in the range of what is observed for the ultra-faint dwarf spheroidals, inherent $M/L$ values of stellar populations can be expected to vary by factors of upwards of 3, interestingly, systematically skewed towards higher values than what corresponds to the infinite population limit. This can explain part of the spread in reported baryonic masses for these systems, which also appear as shifted systematically towards high dark to baryonic mass ratios at fixed stellar velocity dispersions.

The theoretical understanding of density waves in disk galaxies starts from the classical WKB perturbative analysis of tight-winding perturbations, the key assumption being that the potential due to the density wave is approximately radial. The above has served as a valuable guide in aiding the understanding of both simulated and observed galaxies, in spite of a number of caveats being present. The observed spiral or bar patterns in real galaxies are frequently only marginally consistent with the tight-winding assumption, often in fact, outright inconsistent. Here we derive a complementary formulation to the problem, by treating quasi-radial density waves under simplified assumptions in the linear regime. We assume that the potential due to the density wave is approximately tangential, and derive the corresponding dispersion relation of the problem. We obtain an instability criterion for the onset of quasi-radial density waves, which allows a clear understanding of the increased stability of the higher order modes, which appear at progressively larger radii, as often seen in real galaxies. The theory naturally yields a range of pattern speeds for these arms which appears constrained by the condition $\Omega_{p}<\Omega_{0} \pm \kappa /m$. For the central regions of galaxies where solid body rotation curves might apply, we find weak bars in the oscillatory regime with various pattern speeds, including counter rotating ones, and a prediction for $\Omega_{p}$ to increase towards the centre, as seen in the rapidly rotating bars within bars of some numerical simulations. We complement this study with detailed numerical simulations of galactic disks and careful Fourier analysis of the emergent perturbations, which support the theory presented.

Using an extensive sample of galaxies selected from the Sloan Digital Sky Survey Data Release 5, we compare the angular momentum distribution of active galactic nucleus (AGN) with non-AGN hosting late-type galaxies. To this end we characterize galactic spin through the dimensionless angular momentum parameter $\lambda$, which we estimate approximately through simple dynamical considerations. Using a volume limited sample, we find a considerable difference when comparing the empirical distributions of $\lambda$ for AGNs and non-AGN galaxies, the AGNs showing typically low $\lambda$ values and associated dispersions, while non-AGNs present higher $\lambda$ values and a broader distribution. A more striking difference is found when looking at $\lambda$ distributions in thin $M_{r}$ cuts, while the spin of non-AGN galaxies presents an anti-correlation with $M_{r}$, with bright (massive) galaxies having low spins, AGN host galaxies present uniform values of $\lambda$ at all magnitudes, a behavior probably imposed by the fact that most late-type AGN galaxies present a narrow range in color, with a typical constant $\lambda$ value. We also find that the fraction of AGN hosting galaxies in our sample strongly depends on galactic spin, increasing dramatically for decreasing $\lambda$. For AGN host galaxies we compute the mass of their supermassive black holes and find that this value tends to be higher for low spin galaxies, even at fixed luminosity, a result that could account, to a certain extent, for the spread on the luminosity black-hole mass relation.

We study the effects of the cluster environment on galactic morphology by defining a dimensionless angular momentum parameter $\lambda_{d}$, to obtain a quantitative and objective measure of galaxy type. The use of this physical parameter allows us to take the study of morphological transformations in clusters beyond the measurements of merely qualitative parameters, e.g. S/E ratios, to a more physical footing. To this end, we employ an extensive Sloan Digital Sky Survey sample (Data Release 7), with galaxies associated with Abell galaxy clusters. The sample contains 121 relaxed Abell clusters and over 51,000 individual galaxies, which guarantees a thorough statistical coverage over a wide range of physical parameters. We find that the median $\lambda_{d}$ value tends to decrease as we approach the cluster center, with different dependences according to the mass of the galaxies and the hosting cluster; low and intermediate mass galaxies showing a strong dependence, while massive galaxies seems to show, at all radii, low $\lambda_{d}$ values. By analysing trends in $\lambda_{d}$ as functions of the nearest neighbour environment, clustercentric radius and velocity dispersion of clusters, we can identify clearly the leading physical processes at work. We find that in massive clusters ($\sigma>700$ km/s), the interaction with the cluster central region dominates, whilst in smaller clusters galaxy-galaxy interactions are chiefly responsible for driving galactic morphological transformations.

In the pursuit of a general formulation for a modified gravitational theory at the non-relativistic level and as an alternative to the dark matter hypothesis, we construct a model valid over a wide variety of astrophysical scales. Through the inclusion of Milgrom's acceleration constant into a gravitational theory, we show that very general formulas can be constructed for the acceleration felt by a particle. Dimensional analysis shows that this inclusion naturally leads to the appearance of a mass-length scale in gravity, breaking its scale invariance. A particular form of the modified gravitational force is constructed and tested for consistency with observations over a wide range of astrophysical environments, from solar system to extragalactic scales. We show that over any limited range of physical parameters, which define a specific class of astrophysical objects, the dispersion velocity of a system must be a power law of its mass and size. These powers appear linked together through a natural constraint relation of the theory. This yields a generalised gravitational equilibrium relation valid for all astrophysical systems. A general scheme for treating spherical symmetrical density distributions is presented, which in particular shows that the fundamental plane of elliptical galaxies, the Newtonian virial equilibrium, the Tully-Fisher and the Faber-Jackson relations, as well as the scalings observed in local dwarf spheroidal galaxies, are nothing but particular cases of that relation when applied to the appropriate mass-length scales. We discuss the implications of this approach for a modified theory of gravity and emphasise the advantages of working with the force, instead of altering Newton's second law of motion, in the formulation of a gravitational theory.

We use a specific form of the interpolation function in the MOND formalism, which optimally accounts for the internal structure of dwarf spheroidal (dSph) galaxies, to explore the consequences it has on the scaling relations seen in these systems.}} {{The particular form of the interpolation function we used leads to a law of gravity that does not degrade the good fit of the MOND proposal on galactic scales, and in fact, slightly improves the accordance with observations at dSph scales. This formalism yields a good description of gravitational phenomena without the need of invoking any still undetected and hypothetically dominant dark matter, in the weak field regime probed by local dSph galaxies.}} {{Isothermal equilibrium density profiles then yield projected surface density profiles for the local dSph galaxies in very good agreement with observational determinations, for values of the relevant parameters as inferred from recent observations of these Galactic satellites. The observed scaling relations for these systems are also naturally accounted for within the proposed scheme, including a previously unrecognised correlation of the inferred mass-to-light ratios of local dSph's with the ages of their stellar populations, which is natural in modified gravity schemes in the absence of dark matter.}} {The results shed some light on the form that the MOND interpolating function may have in the most challenging regime, which occurs at moderate accelerations and intermediate mass--weighted lengths.

We study the growth rates of massive black holes in the centres of galaxies from accretion of dark matter from their surrounding haloes. By considering only the accretion due to dark matter particles on orbits unbound to the central black hole, we obtain a firm lower limit to the resulting accretion rate. We find that a runaway accretion regime occurs on a timescale which depends on the three characteristic parameters of the problem: the initial mass of the black hole, and the volume density and velocity dispersion of the dark matter particles in its vicinity. An analytical treatment of the accretion rate yields results implying that for the largest black hole masses inferred from QSO studies ($>10^{9} M_{\odot}$), the runaway regime would be reached on time scales which are shorter than the lifetimes of the haloes in question for central dark matter densities in excess of $250 M_{\odot}$pc$^{-3}$. Since reaching runaway accretion would strongly distort the host dark matter halo, the inferences of QSO black holes in this mass range lead to an upper limit on the central dark matter densities of their host haloes of $\rho_{0} < 250 M_{\odot} $pc$^{-3}$. This limit scales inversely with the assumed central black hole mass. However, thinking of dark matter profiles as universal across galactic populations, as cosmological studies imply, we obtain a firm upper limit for the central density of dark matter in such structures.

We search for correlations between the spin in pairs of spiral galaxies, to study if the angular momentum gain for each galaxy was the result of tidal torques imprint by the same tidal field. To perform our study we made use of a sample of galaxy pairs identified using the Sloan Digital Sky Survey. We find a weak, but statistically significant correlation between the spin magnitude of neighbouring galaxies, but no clear alignment between their orientation. We show that events such as interactions with close neighbours play an important role in the value of the spin for the final configuration, as we find these interactions tend to reduce the value of the $\lambda$ spin parameter of late-type galaxies considerably, with dependence on the morphology of the neighbour. This implies that the original tidal field for each pair could have been similar, but the redistribution of angular momentum at later stages of evolution is important.

To study if the angular momentum gain for each member of a galaxy pair was the result of tidal torques imprinted by the same tidal field, we search for correlations between the spin in pairs of spiral galaxies identified using the Sloan Digital Sky Survey. We find a weak, but statistically significant correlation between the spin magnitude of neighbouring galaxies. We show that events such as close interactions with neighbours play an important role in the value of the spin for the final configuration, as we find these interactions tend to reduce the value of the $\lambda$ spin parameter of late-type galaxies considerably. This implies that the original tidal field for each pair could have been similar, but the redistribution of angular momentum at later stages of evolution is important

We present heuristic arguments suggesting that if EM waves with wavelengths somewhat larger than the Schwarzschild radius of a black hole were fully absorbed by it, the second law of thermodynamics would be violated, under the Bekenstein interpretation of the area of a black hole as a measure of its entropy. Thus, entropy considerations make the well known fact that large wavelengths are only marginally absorbed by black holes, a natural consequence of thermodynamics. We also study numerically the ingoing radial propagation of a scalar field wave in a Schwarzschild metric, relaxing the standard assumption which leads to the eikonal equation, that the wave has zero spatial extent. We find that if these waves have wavelengths larger that the Schwarzschild radius, they are very substantially reflected, fully to numerical accuracy. Interestingly, this critical wavelength approximately coincides with the one derived from entropy considerations of the EM field, and is consistent with well known limit results of scattering in the Schwarzschild metric. The propagation speed is also calculated and seen to differ from the value $c$, for wavelengths larger than $R_{s}$, in the vicinity of $R_{s}$. As in all classical wave phenomena, whenever the wavelength is larger or comparable to the physical size of elements in the system, in this case changes in the metric, the zero extent 'particle' description fails, and the wave nature becomes apparent.

By using a very simple and general model to describe the dynamics of disk galaxies, we estimate the $\lambda$ spin parameter for a sample of observed galaxies and present a study in which we show that several important physical properties are intrinsically related to the dynamics of the systems. Although correlations between average metallicity with magnitude or Hubble type are evident, we obtain equally strong correlations with the spin parameter, where galaxies with low $\lambda $ values present higher abundances and galaxies with high $\lambda$ values are poor in metals. Also, the gas content of the galaxies correlates with $\lambda$, with high $\lambda$ systems showing higher gas mass fractions than low $\lambda$ galaxies, highlighting the important role this parameter plays in the structure of disk galaxies and the proposal of $\lambda$ as a robust and objective physical measure of galactic morphology.

When estimating the abundances which result from a given star formation event, it is customary to treat the IMF as a series of weight factors to be applied to the stellar yields, as a function of mass, implicitly assuming one is dealing with an infinite population. However, when the stellar population is small, the standard procedure would imply the inclusion of fractional numbers of stars at certain masses. We study the effects of small number statistics on the resulting abundances by performing an statistical sampling of the IMF to form a stellar population out of discrete numbers of stars. A chemical evolution code then follows the evolution of the population, and traces the resulting abundances. The process is repeated to obtain an statistical distribution of the resulting abundances and their evolution. We explore the manner in which different elements are affected, and how different abundances converge to the infinite population limit as the total mass increases. We include a discussion of our results in the context of dwarf spheroidal galaxies and show the recently reported internal dispersions in abundance ratios for dSph galaxies might be partly explained through the stochastic effects introduced by a low star formation rate, which can account for dispersions of over 2 dex in [C/O], [N/O], [C/Fe], [N/Fe] and [O/Fe].

We use a sample of galaxies from the Sloan Digital Sky Survey (SDSS) to search for correlations between the $\lambda$ spin parameter and the environment and mass of galaxies. In order to calculate the total value of $\lambda$ for each observed galaxy, we employed a simple model of the dynamical structure of the galaxies, which allows a rough estimate of the value of $\lambda$ using only readily obtainable observables from the luminous galaxies. Use of a large volume-limited sample (upwards of 11,000) allows reliable inferences of mean values and dispersions of $\lambda$ distributions. We find, in agreement with some N-body cosmological simulations, no significant dependence of $\lambda$ on the environmental density of the galaxies. For the case of mass, our results show a marked correlation with $\lambda$, in the sense that low-mass galaxies present both higher mean values of $\lambda$ and associated dispersions, than high-mass galaxies. These results provide interesting constrain on the mechanisms of galaxy formation and acquisition of angular momentum, a valuable test for cosmological models.

We use a sample of galaxies from the Sloan Digital Sky Survey (SDSS) to search for correlations between the $\lambda$ spin parameter and the environment and mass of galaxies. In order to calculate the total value of $\lambda$ for each observed galaxy, we employed a simple model of the dynamical structure of the galaxies which allows a rough estimate of the value of $\lambda$ using only readily obtainable observables from the luminous galaxies. Use of a large volume limited sample (upwards of 11,000) allows reliable inferences of mean values and dispersions of $\lambda$ distributions. We find, in agreement with some N-body cosmological simulations, no significant dependence of $\lambda$ on the environmental density of the galaxies. For the case of mass, our results show a marked correlation with $\lambda$, in the sense that low mass galaxies present both higher mean values of $\lambda$ and associated dispersions, than high mass galaxies. This last direct empirical result, at odds with expectations from N-body cosmological simulations, provides interesting constrain on the mechanisms of galaxy formation and acquisition of angular momentum, a valuable test for cosmological models.

We estimate the timescales for orbital decay of wide binaries embedded within dark matter halos, due to dynamical friction against the dark matter particles. We derive analytical scalings for this decay and calibrate and test them through the extensive use of N-body simulations, which accurately confirm the predicted temporal evolution. For density and velocity dispersion parameters as inferred for the dark matter halos of local dSph galaxies, we show that the decay timescales become shorter than the ages of the dSph stellar populations for binary stars composed of 1 $M_{\odot}$ stars, for initial separations larger than 0.1 pc. Such wide binaries are conspicuous and have been well measured in the solar neighborhood. The prediction of the dark matter hypothesis is that they should now be absent from stellar populations embedded within low velocity dispersion, high density dark mater halos, as currently inferred for the local dSph galaxies, having since evolved into tighter binaries. Relevant empirical determinations of this will become technically feasible in the near future, and could provide evidence to discriminate between dark matter particle halos or modified gravitational theories, to account for the high dispersion velocities measured for stars in local dSph galaxies.

The colour-magnitude diagrams of resolved single stellar populations, such as open and globular clusters, have provided the best natural laboratories to test stellar evolution theory. Whilst a variety of techniques have been used to infer the basic properties of these simple populations, systematic uncertainties arise from the purely geometrical degeneracy produced by the similar shape of isochrones of different ages and metallicities. Here we present an objective and robust statistical technique which lifts this degeneracy to a great extent through the use of a key observable: the number of stars along the isochrone. Through extensive Monte Carlo simulations we show that, for instance, we can infer the four main parameters (age, metallicity, distance and reddening) in an objective way, along with robust confidence intervals and their full covariance matrix. We show that systematic uncertainties due to field contamination, unresolved binaries, initial or present-day stellar mass function are either negligible or well under control. This technique provides, for the first time, a proper way to infer with unprecedented accuracy the fundamental properties of simple stellar populations, in an easy-to-implement algorithm.

We investigate the success and problems of MOdified Newtonian Dynamics (MOND) in explaining the structural parameters and dynamics of remote Galactic globular clusters (GCs) and dwarf spheroidal (dSph) galaxies. Using the MOND value for the mass of the Milky Way as inferred from the Galactic rotation curve, we derive the tidal radii of Galactic GCs, and compare to observed values. Except for Pal 14, the predicted tidal radii of GCs are systematically larger than the observed nominal values. However, the associated uncertainties are so large that tidal radii are consistent on the $1\sigma$ level. We have considered the importance of the Galactic tidal force on the survival of dSphs under MOND. Assuming mass-to-light ratios compatible with a naked stellar population, we found that the present Galactic dSphs preserve their integrity over one Hubble time, except Sextans which may survive the tidal interaction only for several Gyr. Mass-to-light ratios as inferred from the internal kinematics of dSph galaxies can be used, but they appear too large to be accounted for only by the stellar population in Willman 1, Coma Berenice, Ursa Minor, Draco, Ursa Major and possibly Bo\"{o}tes dwarves. Finally, the ability of the Sculptor dwarf to retain the observed population of low-mass X-ray binaries (LMXBs) is examined. Under the MOND paradigm, we find that the retention fraction in Sculptor is likely not larger than a few percent. Compared to the dark matter scenario, it turns out that MOND makes the retention problem worse. We propose that measurements of the radial velocities of the observed LMXBs in Sculptor could provide a way to distinguish between modified gravities or extended and massive dark matter halos.

Using simple dimensional arguments for both spiral and elliptical galaxies, we present formulas to derive an estimate of the halo spin parameter $\lambda$ for any real galaxy, in terms of common observational parameters. This allows a rough estimate of $\lambda$, which we apply to a large volume limited sample of galaxies taken from the SDSS data base. The large numbers involved (11,597) allow the derivation of reliable $\lambda$ distributions, as signal adds up significantly in spite of the errors in the inferences for particular galaxies. We find that if the observed distribution of $\lambda$ is modeled with a log-normal function, as often done for this distribution in dark matter halos that appear in cosmological simulations, we obtain parameters $\lambda_{0}=0.04 \pm 0.005$ and $\sigma_{\lambda}=0.51 \pm 0.05$, interestingly consistent with values derived from simulations. For spirals, we find a good correlation between empirical values of $\lambda$ and visually assigned Hubble types, highlighting the potential of this physical parameter as an objective classification tool.

We show that if the total internal energy of a black hole is constructed as the sum of $N$ photons all having a fixed wavelength chosen to scale with the Schwarzschild radius as $\lambda=\alpha R_{s}$, then $N$ will scale with $R_{s}^{2}$. A statistical mechanical calculation of the configuration proposed yields (\alpha = 4 \pi^2 / \ln(2)) and a total entropy of the system $S=k_{B}N \ln(2)$, in agreement with the Bekenstein entropy of a black hole . It is shown that the critical temperature for Bose-Einstein condensation for relativistic particles of $\lambda=\alpha R_{s}$ is always well below the Hawking temperature of a black hole, in support of the proposed internal configuration. We then examine our results from the point of view of recent loop quantum gravity ideas and find that a natural consistency of both approaches appears. We show that the Jeans criterion for gravitational instability can be generalised to the special and general relativistic regimes and holds for any type of mass--energy distribution.

We investigate the in-spiraling timescales of globular clusters in dwarf spheroidal (dSph) and dwarf elliptical (dE) galaxies, due to dynamical friction. We address the problem of these timescales having been variously estimated in the literature as much shorter than a Hubble time. Using self-consistent two-component (dark matter and stars) models, we explore mechanisms which may yield extended dynamical friction timescales in such systems in order to explain why dwarf galaxies often show globular cluster systems. As a general rule, dark matter and stars both give a comparable contribution to the dynamical drag. By exploring various possibilities for their gravitational make-up, it is shown that these studies help constrain the parameters of the dark matter haloes in these galaxies, as well as to test alternatives to dark matter. Under the assumption of a dark haloes having a constant density core, dynamical friction timescales are naturally extended upwards of a Hubble time. Cuspy dark haloes yield timescales $\lesssim$ 4.5 Gyr, for any dark halo parameters in accordance with observations of stellar line-of-sight velocity dispersion in dwarf spheroidal galaxies. We find that under the hypothesis of MOND dynamics, due to the enhanced dynamical drag of the stars, the dynamical friction timescales would be extremely short. Taking the well-measured structural parameters of the Fornax dSph and its globular cluster system as a case study, we conclude that requiring dynamical friction timescales comparable to the Hubble time strongly favours dark haloes with a central core.

We present a highly simplified model of the dynamical structure of a disk galaxy where only two parameters fully determine the solution, mass and angular momentum. We show through simple physical scalings that once the mass has been fixed, the angular momentum parameter $\lambda$ is expected to regulate such critical galactic disk properties as colour, thickness of the disk and disk to bulge ratio. It is hence expected to be the determinant physical ingredient resulting in a given Hubble type. A simple analytic estimate of $\lambda$ for an observed system is provided. An explicit comparison of the distribution of several galactic parameters against both Hubble type and $\lambda$ is performed using observed galaxies. Both such distributions exhibit highly similar characteristics for all galactic properties studied, suggesting $\lambda$ as a physically motivated classification parameter for disk galaxies.

We present a highly simplified model of the dynamical structure of a disk galaxy where only two parmeters fully determine the solution, mass and angular momentum. We show through simple physical scalings that once the mass has been fixed, the angular momentum parameter $\lambda$ is expected to regulate such critical galactic disk properties as colour, thickness of the disk and disk to bulge ratio. It is hence expected to be the determinant physical ingeredient resulting in a given Hubble type. A simple analytic estimate of $\lambda$ for an observed system is provided. An explicit comparison of the distribution of several galactic parameters against both Hubble type and $\lambda$ is performed using observed galaxies. Both such distributions exhibit highly similar characteristics for all galactic properties studied, suggesting $\lambda$ as a physically motivated classification parameter for disk galaxies.

The idea of a unified model for quasars and micro-quasars has been considered for a long time, despite the different environments and physical conditions where both classes of objects reside. Here we show the existence of a simple scaling law, relating the maximum size of a jet to the properties of the gas medium into which it expands. This appears to be valid for all types of hydrodynamical jets, and can be thought of as a broad unified model. The expansion velocity of the jet and the physical properties of the surrounding gas combine in such a way that a limit to the maximum extent of jets at different scales can be obtained.

Terminal wind velocity and mass loss rate are the most fundamental parameters of stellar winds. Unfortunately, their determination calls for high resolution spectroscopy in a range of wavelengths spanning from the ultraviolet to the infrared. For weak and/or distant objects, this becomes unfeasible. Nevertheless, it is possible to obtain an accurate estimate of these parameters through a simplified study of the formation processes of resonant lines which show P Cyg profiles. In this case, the line profile is a complex function of 6 parameters, with the radiative transport treatment becoming relatively inexpensive. However, preforming a 6 dimensional parameter fit raises a number of problems if one seeks an objective and automatic procedure to yield the optimal values, from which wind velocities and mass loss rates can be estimated. Useing a Likelihood function to construct a well defined statistical estimator of the goodness of fit which corresponds to a given model, we turn to a genetic algorithm through which we find the global maximum of the 6 dimensional Likelihood hyper-surface. We here present the implementation of the method, its successful testing with synthetic line profiles, where the answer is known in advance, together with first results of its application to real data

The gas response to a proposed spiral stellar pattern for our Galaxy is presented here as calculated via 2D hydrodynamic calculations utilizing the ZEUS code in the disk plane. The locus is that found by Drimmel (2000) from emission profiles in the K band and at 240 $\mu m$. The self-consistency of the stellar spiral pattern was studied in previous work (see Martos et al. 2004). It is a sensitive function of the pattern rotation speed, $\Omega_p$, among other parameters which include the mass in the spiral and its pitch angle. Here we further discuss the complex gaseous response found there for plausible values of $\Omega_p$ in our Galaxy, and argue that its value must be close to $20 km s^{-1} kpc^{-1}$ from the strong self-consistency criterion and other recent, independent studies which depend on such parameter. However, other values of $\Omega_p$ that have been used in the literature are explored to study the gas response to the stellar (K band) 2-armed pattern. For our best fit values, the gaseous response to the 2-armed pattern displayed in the K band is a four-armed pattern with complex features in the interarm regions. This response resembles the optical arms observed in the Milky Way and other galaxies with the smooth underlying two-armed pattern of the old stellar disk populations in our interpretation.

We report calculations of the stellar and gaseous response to a Milky Way mass distribution model including a spiral pattern with a locus as traced by K-band observations, over imposed on the axisymmetric components in the plane of the disk. The stellar study extends calculations from previous work concerning the self-consistency of the pattern. The stellar response to the imposed spiral mass is studied via computations of the central family of periodic and nearby orbits as a function of the pattern rotation speed, $\Omega_p$, among other parameters. A fine grid of values of $\Omega_p$ was explored ranging from 12 to 25 $km s^{-1} kpc^{-1}$. Dynamical self-consistency is highly sensitive to $\Omega_p$, with the best fit appearing at 20 $km s^{-1} kpc^{-1}$. We give an account of recent independent pieces of theoretical and observational work that are dependent on the value of $\Omega_p$, all of which are consistent with the value found here; the recent star formation history of the Milky Way, local inferences of cosmic ray flux variations and Galactic abundance patterns. The gaseous response, which is also a function of $\Omega_p$, was calculated via 2D hydrodynamic simulations with the ZEUS code. For $\Omega_p = 20 km s^{-1} kpc^{-1}$, the response to a two-armed pattern is a structured pattern of 4 arms, with bifurcations along the arms and interarm features. The pattern resembles qualitatively the optical arms observed in our Galaxy and other galaxies. The complex gaseous pattern appears to be linked to resonances in stellar orbits. Among these, the 4:1 resonance plays an important role, as it also determines the extent over which the imposed spiral pattern is dynamically self-consistent.

Within the context of the standard structure formation scenario, massive present day elliptical galaxies are sometimes thought of as the result of a major merger of spiral systems. Through extensive SPH simulations of merging spirals, we have explored these processes with the aim of quantifying their relaxation times. This is important, as it sets a minimum time interval between the onset of a merger, and the appearance of an elliptical galaxy. We then compare this constraint with predictions of the hierarchical scenario, computed through Press-Schechter merger trees. We find evidence for elliptical systems which appear not to have been formed by a major merger of spirals.

This paper was withdrawn by the authors due to serious flaws in the physics that render the results meaningless.

We examine isothermal dark matter halos in hydrostatic equilibrium with a cosmological constant Lambda =Omega_\Lambda rho_{crit}c^2, where Omega_\Lambda=0.7, and rho_{crit} is the present value of the critical density with h=0.65. The Newtonian limit of General Relativity yields equilibrium equations that are different from those arising by merely coupling an ``isothermal sphere'' to the Lambda-field within a Newtonian framework. The conditions for the existence and stability of circular geodesic orbits show the existence of (I) an ``isothermal region'' (0<r<r_2), circular orbits are stable and all variables behave almost identically to those of an isothermal sphere; (II) an ``asymptotic region'' (r>r_1) dominated by the Lambda-field, where the Newtonian potential oscillates and circular orbits exist in disconnected patches of the domain of r; (III) a ``transition region'' (r_2<r<r_1), circular orbits exist but are unstable. We also find that no stable configuration exists with central density, rho_c, smaller than 2 Lambda, hence any galactic haloes which virialized at $z< 30$ in must have rho_c >0.008 M_\odot {pc}^3, in agreement with rotation curve studies of dwarf galaxies. Since r_2 marks the largest radius of a stable circular orbit, it provides a ``cut off'' radius. For current estimates of rho_c and velocity dispersion of galactic structures, this is around five times larger than the virialization radius. The effects of the Lambda$--field can hence be ignored in structure formation models, but could be significant in the dynamics of superclusters in the linear regime or of gravitational clustering at large scales.

We calculate chemical evolution models for 4 dwarf spheroidal satellites of the Milky Way (Carina, Ursa Minor, Leo I and Leo II) for which reliable non-parametric star formation histories have been derived. In this way, the independently obtained star formation histories are used to constrain the evolution of the systems we are treating. This allows us to obtain robust inferences on the history of such crucial parameters of galactic evolution as gas infall, gas outflows and global metallicities for these systems. We can then trace the metallicity and abundance ratios of the stars formed, the gas present at any time within the systems and the details of gas ejection, of relevance to enrichment of the galaxies environment. We find that galaxies showing one single burst of star formation (Ursa Minor and Leo II) require a dark halo slightly larger that the current estimates for their tidal radii, or the presence of a metal rich selective wind, which might carry away much of the energy output of their supernovae before this might have interacted and heated the gas content, for the gas to be retained until the observed stellar population have formed. Systems showing extended star formation histories (Carina and Leo I) however, are consistent with the idea that their tidally limited dark haloes provide the necessary gravitational potential wells to retain their gas. The complex time structure of the star formation in these systems, remains difficult to understand. Observations of detailed abundance ratios for Ursa Minor, strongly suggest that the star formation history of this galaxy might in fact resemble the complex picture presented by Carina or Leo I, but localized at a very early epoch.

We use a cosmological galactic evolutionary approach to model the Milky Way. A detailed treatment of the mass aggregation and dynamical history of the growing dark halo is included, together with a self consistent physical treatment for the star formation processes within the growing galactic disc. This allows us to calculate the temporal evolution of star and gas surface densities at all galactic radii, in particular, the star formation history (SFH) at the solar radius. A large range of cosmological mass aggregation histories (MAHs) is capable of producing a galaxy with the present day properties of the Milky Way. The resulting SFHs for the solar neighbourhood bracket the available observational data for this feature, the most probable MAH yielding the optimal comparison with these observations. We also find that the rotation curve for our Galaxy implies the presence of a constant density core in its dark matter halo.

Using a cosmological galactic evolutionary approach to model the Milky Way, we calculate the star formation history (SFH) of the solar neighborhood. The good agreement we obtain with the observational inferences suggests that our physical model describes accurately the long term/large spatial trends of the local and global Milky Way SFH. In this model, star formation is triggered by disk gravitational instabilities and self-regulated by an energy balance in the ISM. The drivers of the SFH are the cosmological gas infall rate and the gas surface density determined by the primordial spin parameter. A LambdaCDM cosmology was used throughout.

We explore the predictions of the standard hierarchical clustering scenario of galaxy formation, regarding the numbers and metallicities of PopIII stars likely to be found within our Galaxy today. By PopIII we shall be referring to stars formed at large redshift ($z>4$), with low metallicities ($[Z/Z_{\odot}]<-2.5$) and in small systems (total mass $\simlt$ $2\times 10^{8} M_{\odot}$) that are extremely sensitive to stellar feedback, and which through a prescribed merging history (Lacey & Cole 1993) end up becoming part of the Milky Way today. An analytic, extended Press-Schechter formalism is used to get the mass functions of halos which will host PopIII stars at a given redshift, and which will end up in Milky Way sized systems today. Each of these is modeled as a mini galaxy, with a detailed treatment of the dark halo structure, angular momentum distribution, final gas temperature and disk instabilities, all of which determine the fraction of the baryons which are subject to star formation. Use of new primordial metallicity stellar evolutionary models allows us to trace the history of the stars formed, give accurate estimates of their expected numbers today, and their location in $L/L_{\odot}$ vs. $T/K$ HR diagrams. A first comparison with observational data suggests that the IMF of the first stars was increasingly high mass weighted towards high redshifts, levelling off at $z\simgt 9$ at a characteristic stellar mass scale $m_s=10-15 M_\odot$.

The unprecedented quality of recent colour-magnitude diagrams of resolved stellar populations in nearby galaxies requires state-of-the-art techniques to infer the star formation histories which gave rise to the observed distributions. We have developed a maximum likelihood technique, which coupled to a variational calculus allows us to make a robust, non-parametric reconstruction of the evolution of the star formation rate. A full Bayesian analysis is also applied to assess whether the best solutions found are also good fits to the data. Applying this new method to WFPC2 observations of five dSph galaxies of the Local Group, we find a wide variety of star formation histories, with no particular epoch being dominant. In the case of the solar neighbourhood observed by Hipparcos we infer, with an unprecedented resolution of 50 Myr, its star formation history over the past 3 Gyr. The surprising regularity of star formation episodes separated by some 0.5 Gyr could be interpreted as the result of interactions with two spiral arms or the Galactic bar. These bursts possibly trigger the formation of massive star clusters which slowly dissolve into the galactic field.

One of the predictions of the standard CDM is that dark haloes have centrally divergent density profiles. An extensive body of rotation curve observations of dwarf and low surface brightness galaxies shows the dark haloes of those systems to be characterized by soft constant density central cores. Several physical processes have been proposed to produce soft cores in dark haloes, each one with different scaling properties. With the aim of discriminating among them we have examined the rotation curves of dark matter dominated dwarf and low surface brightness galaxies and the inner mass profiles of two clusters of galaxies lacking a central cD galaxy and with evidence of soft cores in the centre. The core radii and central densities of these haloes scale in a well defined manner with the depth of their potential wells, as measured through the maximum circular velocity. As a result of our analysis we identify self-interacting CDM as a viable solution to the core problem, where a non-singular isothermal core is formed in the halo center surrounded by a Navarro, Frenk, & White profile in the outer parts. We show that this particular physical situation predicts core radii in agreement with observations. Furthermore, using the observed scalings, we derive an expression for the minimum cross section (\sigma) which has an explicit dependence with the halo dispersion velocity (v). If m_x is the mass of the dark matter particle: \sigma/m_x ~4 10^-25 (v/100 km s^-1)^-1 cm^2/Gev.

Using observational data for a sample of dark matter dominated galaxies and two cluster of galaxies, we have found that the central halo density does not depend on its mass, and the core radius is roughly proportional to the maximum rotation velocity. A good agreement with these scaling laws is obtained for CDM halos whose dense inner parts were expanded by gravothermal instabilities if the particles efficiently self-interact only in these parts. We find that the particle cross-section is inversely proportional to the velocity dispersion.

Within the framework of the cold dark matter (CDM) cosmogony, a central cusp in the density profiles of virialized dark haloes is predicted. This prediction disagrees with the soft inner halo mass distribution inferred from observations of dwarf and low surface brightness galaxies, and some clusters of galaxies. By analysing data for some of these objects, we find that the halo central density is nearly independent of the mass from galactic to galaxy cluster scales with an average value of around 0.02 M_sun/pc^3. We show that soft cores can be produced in the CDM haloes by introducing a lower cut-off in the power spectra of fluctuations and assuming high orbital thermal energies during halo formation. However, the scale invariance of the halo central density is not reproduced in these cases. The introduction of self-interaction in the CDM particles offers the most attractive alternative to the core problem. We propose gravothermal expansion as a possible mechanism to produce soft cores in the CDM haloes with self-interacting particles. A global thermodynamical equilibrium can explain the central density scale invariance. We find a minimum cross section capable of establishing isothermal cores in agreement with the observed shallow cores. If \sigma and m_x are the cross section and mass of the dark matter particle, and v is the halo velocity dispersion, then (\sigma/m_x) \~ 4 10^{-25} (100 km s^{-1}/v) cm^2/GeV.

We use data from the Hipparcos catalogue to construct colour-magnitude diagrams for the solar neighbourhood, which are then treated using advanced Bayesian analysis techniques to derive the star formation history, $SFR(t)$, of this region over the last 3 Gyr. The method we use allows the recovery of the underlying $SFR(t)$ without the need of assuming any {\it a priori} structure or condition on $SFR(t)$, and hence yields a highly objective result. The remarkable accuracy of the data permits the reconstruction of the local $SFR(t)$ with an unprecedented time resolution of $\approx 50 $ Myr. A $SFR(t)$ having an oscillatory component of period $\approx 0.5$ Gyr is found, superimposed on a small level of constant star formation activity. Problems arising from the non-uniform selection function of the Hipparcos satellite are discussed and treated. Detailed statistical tests are then performed on the results, which confirm the inferred $SFR(t)$ to be compatible with the observed distribution of stars.

The evolution of the star formation rate in the Galaxy is one of the key ingredients quantifying the formation and determining the chemical and luminosity evolution of galaxies. Many complementary methods exist to infer the star formation history of the components of the Galaxy, from indirect methods for analysis of low-precision data, to new exact analytic methods for analysis of sufficiently high quality data. We summarise available general constraints on star formation histories, showing that derived star formation rates are in general comparable to those seen today. We then show how colour-magnitude diagrams of volume- and absolute magnitude-limited samples of the solar neighbourhood observed by Hipparcos may be analysed, using variational calculus techniques, to reconstruct the local star formation history. The remarkable accuracy of the data coupled to our maximum-likelihood variational method allows objective quantification of the local star formation history with a time resolution of ~ 50 Myr. Over the past 3Gyr, the solar neighbourhood star formation rate has varied by a factor of ~ 4, with characteristic timescale about 0.5Gyr, possibly triggered by interactions with spiral arms.

We use recent HST colour-magnitude diagrams of the resolved stellar populations of a sample of local dSph galaxies (Carina, LeoI, LeoII, Ursa Minor and Draco) to infer the star formation histories of these systems, $SFR(t)$. Applying a new variational calculus maximum likelihood method which includes a full Bayesian analysis and allows a non-parametric estimate of the function one is solving for, we infer the star formation histories of the systems studied. This method has the advantage of yielding an objective answer, as one need not assume {\it a priori} the form of the function one is trying to recover. The results are checked independently using Saha's $W$ statistic. The total luminosities of the systems are used to normalize the results into physical units and derive SN type II rates. We derive the luminosity weighted mean star formation history of this sample of galaxies.

Late accretion models for formation of the Galactic halo require that many Galactic satellite galaxies have been cannibalised into the halo field. Comparison of the metallicity and age distribution function of stars in the surviving satellites with the apparently exclusively old stars in the field halo can constrain the importance of any such process. We have developed a new objective technique to determine star formation histories in dSph galaxies. We apply this technique to the surviving Galactic satellites, deducing an approximately uniform distribution of ages for the constituents, quite unlike the halo field stars. Thus, late accretion did not play a substantial part in Galactic halo formation.

Objective determination of a star formation history from a colour-magnitude diagram, independently of assumed parametric descriptions, is necessary to determine the evolutionary history of a galaxy. We introduce a new method for solving maximum likelihood problems through variational calculus, and apply it to the case of recovering an unknown star formation history, SFR(t), from a colour-magnitude diagram. This method provides a non-parametric solution with the advantage of requiring no initial assumptions on the SFR(t). As a full maximum likelihood statistical model is used, one may exploit all the information available in the photometric data which can be explained by available isochrones. Extensive tests of the method show it to be reliable under noise conditions comparable to those appropriate for real observations. Implications for future surveys, such as the GAIA mission are significant: accurate determination of the evolutionary history of the Milky Way is limited only by adequate spectroscopy and/or photometry to provide a determination of stellar metallicity and temperature for a sufficient sample of stars of luminosities at and above an old main-sequence turnoff.

We introduce a new method for solving maximum likelihood problems through variational calculus, and apply it to the case of recovering an unknown star formation history, $SFR(t)$, from a resulting HR diagram. This approach allows a totally non-parametric solution which has the advantage of requiring no initial assumptions on the $SFR(t)$. As a full maximum likelihood statistical model is used, we take advantage of all the information available in the HR diagram, rather than concentrating on particular features such as turn off points or luminosity functions. We test the method using a series of synthetic HR diagrams produced from known $SFR(t)$, and find it to be quite successful under noise conditions comparable to those present in current observations. At this point we restrict the analysis to situations where the metallicity of the system is known, as is the case with the resolved populations of the dwarf spheroidal companions to the Milky Way or the solar neighbourhood Hipparcos data. We also include tests to quantify the way uncertainties in the assumed metallicity, binary fraction and IMF affect our inferences.

We present a simplified analytic approach to the problem of the spiraling of a massive body orbiting within the dark halo of a dwarf galaxy. This dark halo is treated as the core region of a King distribution of dark matter particles, in consistency with the observational result of dwarf galaxies having solid body rotation curves. Thus we derive a simple formula which provides a reliable and general first order solution to the problem, totally analogous to the one corresponding to the dynamical friction problem in an isothermal halo. This analytic approach allows a clear handling and a transparent understanding of the physics and the scaling of the problem. A comparison with the isothermal case shows that in the core regions of a King sphere, dynamical friction proceeds at a different rate, and is sensitive to the total core radius. Thus, in principle, observable consequences may result. In order to illustrate the possible effects, we apply this formula to the spiraling of globular cluster orbits in dwarf galaxies, and show how present day globular cluster systems could in principle be used to derive better limits on the structure of dark halos around dwarf galaxies, when the observational situation improves. As a second application, we study the way a massive black hole population forming a fraction of these dark halos would gradually concentrate towards the centre, with the consequent deformation of an originally solid body rotation curve. This effect allows us to set limits on the fraction/mass of any massive black hole minority component of the dark halos of dwarf galaxies. In essence, we take advantage of the way the global matter distribution fixes the local distribution function for the dark matter particles, which in turn determines the dynamical friction problem.

We re-examine the Fall & Efstathiou (1980) scenario for galaxy formation, including the dark halo gravitational reaction to the formation of the baryon disk, as well as continuous variations in the intrinsic halo density profile. The recently published rotation curves of LSB and dwarf galaxies together with previously known scaling relations provide sufficient information on the present day structure of late type disk galaxies to invert the problem. By requiring that the models reproduce all the observational restrictions we can fully constrain the initial conditions of galaxy formation, with a minimum of assumptions, in particular without the need to specify a CDM halo profile. This allows one to solve for all the initial conditions, in terms of the halo density profile, the baryon fraction and the total angular momentum. We find that a unique initial halo shape is sufficient to accurately reproduce the rotation curves of both LSB and normal late type spiral galaxies. This unique halo profile differs substantially from that found in standard CDM models (Navarro, Frenk & White 1996). A galactic baryon fraction of 0.065 is found. The initial value of the dimensionless angular momentum is seen to be the principal discriminator between the galaxy classes we examine. The present day scalings between structural parameters are seen to originate in the initial conditions.