Source author record

Hongsheng Zhao

Hongsheng Zhao appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

44works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

44 published item(s)

preprint2022arXiv

From galactic bars to the Hubble tension: weighing up the astrophysical evidence for Milgromian gravity

Astronomical observations reveal a major deficiency in our understanding of physics $-$ the detectable mass is insufficient to explain the observed motions in a huge variety of systems given our current understanding of gravity, Einstein's General theory of Relativity (GR). This missing gravity problem may indicate a breakdown of GR at low accelerations, as postulated by Milgromian dynamics (MOND). We review the MOND theory and its consequences, including in a cosmological context where we advocate a hybrid approach involving light sterile neutrinos to address MOND's cluster-scale issues. We then test the novel predictions of MOND using evidence from galaxies, galaxy groups, galaxy clusters, and the large-scale structure of the Universe. We also consider whether the standard cosmological paradigm ($Λ$CDM) can explain the observations and review several previously published highly significant falsifications of it. Our overall assessment considers both the extent to which the data agree with each theory and how much flexibility each has when accommodating the data, with the gold standard being a clear $a~priori$ prediction not informed by the data in question. Our conclusion is that MOND is favoured by a wealth of data across a huge range of astrophysical scales, ranging from the kpc scales of galactic bars to the Gpc scale of the local supervoid and the Hubble tension, which is alleviated in MOND through enhanced cosmic variance. We also consider several future tests, mostly on scales much smaller than galaxies.

preprint2022arXiv

Overestimated inclinations of Milgromian disc galaxies: the case of the ultradiffuse galaxy AGC 114905

We present hydrodynamical star-forming simulations in the Milgromian dynamics (MOND) framework of a gas-rich disc galaxy with properties similar to AGC 114905, which has recently been argued to have a rotation curve (RC) that is inconsistent with the MOND prediction. Our first model considers the galaxy in isolation, while our second model includes an external field of $0.05 \, a_{_0}$, the estimated gravitational field from large-scale structure. We show that isophotes in the face-on view can differ from circular at the 50% level. This could mislead observers into overestimating the inclination $i$ between disc and sky planes. Because RCs require a correction factor of $1/\sin i$, the actual RC could be much higher than reported by observers. This plausibly reconciles AGC 114905 with MOND expectations.

preprint2020arXiv

Testing Gravity with wide binary stars like $α$ Centauri

We consider the feasibility of testing Newtonian gravity at low accelerations using wide binary (WB) stars separated by $\ge 3$ kAU. These systems probe the accelerations at which galaxy rotation curves unexpectedly flatline, possibly due to Modified Newtonian Dynamics (MOND). We conduct Newtonian and MOND simulations of WBs covering a grid of model parameters in the system mass, semi-major axis, eccentricity and orbital plane. We self-consistently include the external field (EF) from the rest of the Galaxy on the Solar neighbourhood using an axisymmetric algorithm. For a given projected separation, WB relative velocities reach larger values in MOND. The excess is ${\approx 20\%}$ adopting its simple interpolating function, as works best with a range of Galactic and extragalactic observations. This causes noticeable MOND effects in accurate observations of ${\approx 500}$ WBs, even without radial velocity measurements. We show that the proposed Theia mission may be able to directly measure the orbital acceleration of Proxima Centauri towards the 13 kAU-distant $α$ Centauri. This requires an astrometric accuracy of $\approx 1 \, μ$as over 5 years. We also consider the long-term orbital stability of WBs with different orbital planes. As each system rotates around the Galaxy, it experiences a time-varying EF because this is directed towards the Galactic Centre. We demonstrate approximate conservation of the angular momentum component along this direction, a consequence of the WB orbit adiabatically adjusting to the much slower Galactic orbit. WBs with very little angular momentum in this direction are less stable over Gyr periods. This novel direction-dependent effect might allow for further tests of MOND.

preprint2020arXiv

Testing the equivalence principle via the shadow of black holes

We study the equivalence principle, regarded as the cornerstone of general relativity, by analyzing the deformation observable of black hole shadows. Such deformation can arise from new physics and may be expressed as a phenomenological violation of the equivalence principle. Specifically, we assume that there is an additional background vector field that couples to the photons around the supermassive black hole. This type of coupling yields impact on the way the system depends on initial conditions, and affects the black hole shadow at different wavelengths by a different amount, and therefore observations of the shadow in different wavelengths could constrain such couplings. This can be tested by future multi-band observations. Adopting a specific form of the vector field, we obtain constraints on model parameters from Event Horizon Telescope observations and measurements of gas/stellar orbits.

preprint2020arXiv

The Escape Velocity Curve of the Milky Way in Modified Newtonian Dynamics

We determine the escape velocity from the Milky Way (MW) at a range of Galactocentric radii in the context of Modified Newtonian Dynamics (MOND). Due to its non-linear nature, escape is possible if the MW is considered embedded in a constant external gravitational field (EF) from distant objects. We model this situation using a fully self-consistent method based on a direct solution of the governing equations out to several thousand disk scale lengths. We try out a range of EF strengths and mass models for the MW in an attempt to match the escape velocity measurements of Williams et al. (2017). A reasonable match is found if the EF on the MW is ${\sim 0.03 a_{_0}}$, towards the higher end of the range considered. Our models include a hot gas corona surrounding the MW, but our results suggest that this should have a very low mass of ${\sim 2 \times 10^{10} M_\odot}$ to avoid pushing the escape velocity too high. Our analysis favours a slightly lower baryonic disk mass than the ${\sim 7 \times 10^{10} M_\odot}$ required to explain its rotation curve in MOND. However, given the uncertainties, MOND is consistent with both the locally measured amplitude of the MW rotation curve and its escape velocity over Galactocentric distances of 8$-$50 kpc.

preprint2020arXiv

The External Field Dominated Solution In QUMOND & AQUAL: Application To Tidal Streams

The standard $Λ$CDM paradigm seems to describe cosmology and large scale structure formation very well. However, a number of puzzling observations remain on galactic scales. An example is the anisotropic distribution of satellite galaxies in the Local Group. This has led to suggestions that a modified gravity theory might provide a better explanation than Newtonian gravity supplemented by dark matter. One of the leading modified gravity theories is Modified Newtonian Dynamics (MOND). For an isolated point mass, it boosts gravity by an acceleration-dependent factor of $ν$. Recently, a much more computer-friendly quasi-linear formulation of MOND (QUMOND) has become available. We investigate analytically the solution for a point mass embedded in a constant external field of $\mathbf{g}_{ext}$. We find that the potential is $Φ= - ~ \frac{GM ν_{ext}}{r}\left(1 + \frac{K_0}{2} \sin^2 θ\right)$, where $r$ is distance from the mass $M$ which is in an external field that `saturates' the $ν$ function at the value $ν_{ext}$, leading to a fixed value of $K_0 \equiv \frac{\partial Ln ~ ν}{\partial Ln ~ g_{ext}}$. In a very weak gravitational field $\left(\left| \mathbf{g}_{ext} \right| \ll a_0 \right)$, $K_0 = -\frac{1}{2}$. The angle $θ$ is that between the external field direction and the direction towards the mass. Our results are quite close to the more traditional aquadratic Lagrangian (AQUAL) formulation of MOND. We apply both theories to a simple model of the Sagittarius tidal stream. We find that they give very similar results, with the tidal stream seeming to spread slightly further in AQUAL.

preprint2016arXiv

Effects of Lens Motion and Uneven Magnification on Image Spectra

Counter to intuition, the images of an extended galaxy lensed by a moving galaxy cluster should have slightly different spectra in any metric gravity theory. This is mainly for two reasons. One relies on the gravitational potential of a moving lens being time-dependent (the $\text{Moving}$ $\text{Cluster}$ $\text{Effect}$, $\text{MCE}$). The other is due to uneven magnification across the extended, rotating source (the $\text{Differential}$ $\text{Magnification}$ $\text{Effect}$, $\text{DME}$). The time delay between the images can also cause their redshifts to differ because of cosmological expansion. This Differential Expansion Effect is likely to be small. Using a simple model, we derive these effects from first principles. One application would be to the Bullet Cluster, whose large tangential velocity may be inconsistent with the $ΛCDM$ paradigm. This velocity can be estimated with complicated hydrodynamic models. Uncertainties with such models can be avoided using the MCE. We argue that the MCE should be observable with ALMA. However, such measurements can be corrupted by the DME if typical spiral galaxies are used as sources. Fortunately, we find that if detailed spectral line profiles were available, then the DME and MCE could be distinguished. It might also be feasible to calculate how much the DME should affect the mean redshift of each image. Resolved observations of the source would be required to do this accurately. The DME is of order the source angular size divided by the Einstein radius times the redshift variation across the source. Thus, it mostly affects nearly edge-on spiral galaxies in certain orientations. This suggests that observers should reduce the DME by careful choice of target, a possibility we discuss in some detail.

preprint2014arXiv

On radial anisotropy limits in stellar systems

Following earlier authors, we re-examine constraints on the radial velocity anisotropy of generic stellar systems using arguments for phase space density positivity, stability, and separability. It is known that although the majority of commonly used systems have an maximum anisotropy of less than half of the logarithmic density slope \emph{i.e.} $β< γ/2$, there are exceptions for separable models with large central anisotropy. Here we present a new exceptional case with above-threshold anisotropy locally but with an isotropic center nevertheless. These models are non-separable and we maintain positivity. Our analysis suggests that regions of above-threshold anisotropy are more related to regions of possible secular instability, which might be observed in self-consistent galaxies in a short-lived phase. s an equilibrium solution for a collisionless system.

preprint2014arXiv

Stirring N-body systems II: Necessary conditions for the dark matter attractor

We study the evolution of the phase-space of collisionless N-body systems under repeated stirrings or perturbations, which has been shown to lead to a convergence towards a limited group of end states. This so-called attractor was previously shown to be independent of the initial system and environmental conditions. However the fundamental reason for its appearance is still unclear. It has been suggested that the origin of the attractor may be either radial infall (RI), the radial orbit instability (ROI), or energy exchange which, for instance, happens during violent relaxation. Here we examine the effects of a set of controlled perturbations, referred to as `kicks', which act in addition to the standard collisionless dynamics by allowing pre-specified instantaneous perturbations in phase-space. We first demonstrate that the attractor persists in the absence of RI and ROI by forcing the system to expand. We then consider radial velocity kicks in a rigid potential and isotropic velocity kicks, since there are no energy exchanges in these two recipes of kicks. We find that these kicks do not lead to the attractor, indicating that the energy exchange in a dynamic potential is indeed the physical mechanism responsible for the attractor.

preprint2013arXiv

Gravity theories, Transverse Doppler and Gravitational Redshifts in Galaxy Clusters

There is growing interest in testing alternative gravity theories using the subtle gravitational redshifts in clusters of galaxies. However, current models all neglect a transverse Doppler redshift of similar magnitude, and some models are not self-consistent. An equilibrium model would fix the gravitational and transverse Doppler velocity shifts to be about 6sigma^2/c and 3sigma^2/2c in order to fit the observed velocity dispersion sigma self-consistently. This result comes from the Virial Theorem for a spherical isotropic cluster, and is insensitive to the theory of gravity. A gravitational redshift signal also does not directly distinguish between the Einsteinian and f(R) gravity theories, because each theory requires different dark halo mass function to keep the clusters in equilibrium. When this constraint is imposed, the gravitational redshift has no sensitivity to theory. Indeed our N-body simulations show that the halo mass function differs in f(R), and that the transverse Doppler effect is stronger than analytically predicted due to non-equilibrium.

preprint2013arXiv

Local Group timing in Milgromian dynamics. A past Milky Way-Andromeda encounter at z>0.8

The Local Group timing has been one of the first historical probes of the missing mass problem. Whilst modern cosmological probes indicate that pure baryonic dynamics is not sufficient on the largest scales, nearby galaxies and small galaxy groups persistently obey Milgrom's MOND law, which implies that dynamics at small scales is possibly entirely predicted by the baryons. Here, we investigate the Local Group timing in this context of Milgromian dynamics. Making use of the latest measured proper motions and radial velocities for Andromeda and the Magellanic clouds, we integrate backwards their orbits by making use of the Milgromian two-body equation of motion. We find that, with the currently measured proper motions and radial velocity of M31, MOND would imply that the Milky Way and M31 first moved apart via Hubble expansion after birth, but then necessarily got attracted again by the Milgromian gravitational attraction, and had a past fly-by encounter before coming to their present positions. This encounter would most probably have happened 7 to 11 Gyr ago (0.8<z<3). The absence of a dark matter halo and its associated dynamical friction is necessary for such a close encounter not to trigger a merger. Observational arguments which could exclude or favour such a past encounter would thus be very important in view of falsifying or vindicating Milgromian dynamics on the scale of the Local Group. Interestingly, the closest approach of the encounter is small enough (<55 kpc) to have severe consequences on the disk dynamics, including perhaps thick disk formation, and on the satellite systems of both galaxies. The ages of the satellite galaxies and of the young halo globular clusters, all of which form the vast polar structure around the Milky Way, are consistent with these objects having been born in this encounter.

preprint2012arXiv

A New Model for the Milky Way Bar

We use Schwarzschild's orbit-superposition technique to construct self-consistent models of the Galactic bar. Using $χ^2$ minimisation, we find that the best-fit Galactic bar model has a pattern speed $Ω_{\rm p}=60 \rm{km s^{-1} kpc^{-1}}$, disk mass $\rm{M_{\rm d}=1.0\times10^{11}M_{\odot}}$ and bar angle $θ_{\rm bar}=20^{\circ}$ for an adopted bar mass $\rm{M_{\rm bar}=2\times10^{10}M_{\odot}}$. The model can reproduce not only the three-dimensional and projected density distributions but also velocity and velocity dispersion data from the BRAVA survey. We also predict the proper motions in the range $l=[-12^{\circ},12^{\circ}]$, $b=[-10^{\circ},10^{\circ}]$, which appear to be higher than observations in the longitudinal direction. The model is stable within a timescale of 0.5 Gyr, but appears to deviate from steady-state on longer timescales. Our model can be further tested by future observations such as those from GAIA.

preprint2012arXiv

Stirring N-body systems: Universality of end states

We study the evolution of the phase-space of collisionless N-body systems under repeated stirrings or perturbations. We find convergence towards a limited solution group, in accordance with Hansen 2010, that is independent of the initial system and environmental conditions, paying particular attention to the assumed gravitational paradigm (Newtonian and MOND). We examine the effects of changes to the perturbation scheme and in doing so identify a large group of perturbations featuring radial orbit instability (ROI) which always lead to convergence. The attractor is thus found to be a robust and reproducible effect under a variety of circumstances.

preprint2012arXiv

Testing quasilinear modified Newtonian dynamics in the Solar System

A unique signature of the modified Newtonian dynamics (MOND) paradigm is its peculiar behavior in the vicinity of the points where the total Newtonian acceleration exactly cancels. In the Solar System, these are the saddle points of the gravitational potential near the planets. Typically, such points are embedded into low-acceleration bubbles where modified gravity theories a` la MOND predict significant deviations from Newton's laws. As has been pointed out recently, the Earth-Sun bubble may be visited by the LISA Pathfinder spacecraft in the near future, providing a unique occasion to put these theories to a direct test. In this work, we present a high-precision model of the Solar System's gravitational potential to determine accurate positions and motions of these saddle points and study the predicted dynamical anomalies within the framework of quasi-linear MOND. Considering the expected sensitivity of the LISA Pathfinder probe, we argue that interpolation functions which exhibit a "faster" transition between the two dynamical regimes have a good chance of surviving a null result. An example of such a function is the QMOND analog of the so-called simple interpolating function which agrees well with much of the extragalactic phenomenology. We have also discovered that several of Saturn's outermost satellites periodically intersect the Saturn-Sun bubble, providing the first example of Solar System objects that regularly undergo the MOND regime.

preprint2012arXiv

Unifying all mass discrepancies with one effective gravity law?

A remarkably tight relation is observed between the Newtonian gravity sourced by the baryons and the actual gravity in galaxies of all sizes. This can be interpreted as the effect of a single effective force law depending on acceleration. This is however not the case in larger systems with much deeper potential wells, such as galaxy clusters. Here we explore the possibility of an effective force law reproducing mass discrepancies in all extragalactic systems when depending on both acceleration and the depth of the potential well. We exhibit, at least at a phenomenological level, one such possible construction in the classical gravitational potential theory. Interestingly, the framework, dubbed EMOND, is able to reproduce the observed mass discrepancies in both galaxies and galaxy clusters, and to produce multi-center systems with offsets between the peaks of gravity and the peaks of the baryonic distribution.

preprint2011arXiv

The Bulge Radial Velocity Assay (BRAVA): II. Complete Sample and Data Release

We present new radial velocity measurements from the Bulge Radial Velocity Assay (BRAVA), a large scale spectroscopic survey of M-type giants in the Galactic bulge/bar region. The sample of ~4500 new radial velocities, mostly in the region -10 deg < l < +10 deg and b ~ -6 deg more than doubles the existent published data set. Our new data extend our rotation curve and velocity dispersion profile to +20 deg, which is ~2.8 kpc from the Galactic Center. The new data confirm the cylindrical rotation observed at -6 deg and -8 deg, and are an excellent fit to the Shen et al. (2010) N-body bar model. We measure the strength of the TiO molecular band as a first step towards a metallicity ranking of the stellar sample, from which we confirm the presence of a vertical abundance gradient. Our survey finds no strong evidence of previously unknown kinematic streams. We also publish our complete catalog of radial velocities, photometry, TiO band strengths, and spectra, which is available at the IRSA archive: http://irsa.ipac.caltech.edu/ as well as at UCLA: http://brava.astro.ucla.edu/.

preprint2011arXiv

The life cycle of star cluster in a tidal field

The evolution of globular clusters due to 2-body relaxation results in an outward flow of energy and at some stage all clusters need a central energy source to sustain their evolution. Henon provided the insight that we do not need to know the details of the energy production in order to understand the relaxation-driven evolution of the cluster, at least outside the core. He provided two self-similar solutions for the evolution of clusters based on the view that the cluster as a whole determines the amount of energy that is produced in the core: steady expansion for isolated clusters, and homologous contraction for clusters evaporating in a tidal field. We combine these models: the half-mass radius increases during the first half of the evolution, and decreases in the second half; while the escape rate approaches a constant value set by the tidal field. We refer to these phases as `expansion dominated' and `evaporation dominated'. These simple analytical solutions immediately allow us to construct evolutionary tracks and isochrones in terms of cluster half-mass density, cluster mass and galacto-centric radius. From a comparison to the Milky Way globular clusters we find that roughly 1/3 of them are in the second, evaporation-dominated phase and for these clusters the density inside the half-mass radius varies with the galactocentric distance R as rho_h ~ 1/R^2. The remaining 2/3 are still in the first, expansion-dominated phase and their isochrones follow the environment-independent scaling rho_h ~ M^2; that is, a constant relaxation time-scale. We find substantial agreement between Milky Way globular cluster parameters and the isochrones, which suggests that there is, as Henon suggested, a balance between the flow of energy and the central energy production for almost all globular clusters.

preprint2010arXiv

A Semi-analytic Ray-tracing Algorithm for Weak Lensing

We propose a new ray-tracing algorithm to measure the weak lensing shear and convergence fields directly from N-body simulations. We calculate the deflection of the light rays lensed by the 3-D mass density field or gravitational potential along the line of sight on a grid-by-grid basis, rather than using the projected 2-D lens planes. Our algorithm uses simple analytic formulae instead of numerical integrations in the computation of the projected density field along the line of sight, and so is computationally efficient, accurate and straightforward to implement. This will prove valuable in the interpretation of data from the next generation of surveys that will image many thousands of square degrees of sky.

preprint2010arXiv

Comparing different realizations of modified Newtonian dynamics: virial theorem and elliptical shells

There exists several modified gravity theories designed to reproduce the empirical Milgrom's formula (MOND). Here we derive analytical results in the context of the static weak-field limit of two of them (BIMOND, leading for a given set of parameters to QUMOND, and TeVeS). In this limit, these theories are constructed to give the same force field for spherical symmetry, but their predictions generally differ out of it. However, for certain realizations of these theories (characterized by specific choices for their free functions), the binding potential-energy of a system is increased, compared to its Newtonian counterpart, by a constant amount independent of the shape and size of the system. In that case, the virial theorem is exactly the same in these two theories, for the whole gravity regime and even outside of spherical symmetry, although the exact force fields are different. We explicitly show this for the force field generated by the two theories inside an elliptical shell. For more general free functions, the virial theorems are however not identical in these two theories. We finally explore the consequences of these analytical results for the two-body force.

preprint2010arXiv

Lopsidedness of cluster galaxies in modified gravity

We point out an interesting theoretical prediction for elliptical galaxies residing inside galaxy clusters in the framework of modified Newtonian dynamics (MOND), that could be used to test this paradigm. Apart from the central brightest cluster galaxy, other galaxies close enough to the centre experience a strong gravitational influence from the other galaxies of the cluster. This influence manifests itself only as tides in standard Newtonian gravity, meaning that the systematic acceleration of the centre of mass of the galaxy has no consequence. However, in the context of MOND, a consequence of the breaking of the strong equivalence principle is that the systematic acceleration changes the own self-gravity of the galaxy. We show here that, in this framework, initially axisymmetric elliptical galaxies become lopsided along the external field's direction, and that the centroid of the galaxy, defined by the outer density contours, is shifted by a few hundreds parsecs with respect to the densest point.

preprint2010arXiv

Mass discrepancy in galaxy clusters as a result of the offset between dark matter and baryon distributions

Recent studies of lensing clusters reveal that it might be fairly common for a galaxy cluster that the X-ray center has an obvious offset from its gravitational center which is measured by strong lensing. We argue that if these offsets exist, then X-rays and lensing are indeed measuring different regions of a cluster, and may thus naturally result in a discrepancy in the measured gravitational masses by the two different methods. Here we investigate theoretically the dynamical effects of such lensing-X-ray offsets, and compare with observational data. We find that for typical values, the offset alone can give rise to a factor of two difference between the lensing and X-ray determined masses for the core regions of a cluster, suggesting that such "offset effect" may play an important role and should not be ignored in our dynamical measurements of clusters.

preprint2010arXiv

Mass of Galactic Lenses in Modified Newtonian Dynamics

Using strong lensing data Milgrom's MOdified Newtonian Dynamics (MOND) or its covariant TeVeS (Tensor-Vector-Scalar Theory) is being examined here as an alternative to the conventional $Λ$CDM paradigm. We examine 10 double-image gravitational lensing systems, in which the lens masses have been estimated by stellar population synthesis models. While mild deviations exist, we do not find out that strong cases for outliers to the TeVeS theory.

preprint2010arXiv

Modified Kepler's Law, Escape Speed and Two-body Problem in MOND-like Theories

We derive a simple analytical expression for the two-body force in a sub-class of MOND-like theories and make testable predictions in the modification to the two-body orbital period, shape, and precession rate, and escape speed etc. We demonstrate the applications of the modified Kepler's law in the timing of satellite orbits around the Milky Way, and checking the feasibility of MOND in the orbit of Large Magellanic Cloud, the M31 galaxy, and the merging Bullet Clusters. MOND appears to be consistent with satellite orbits although with a tight margin. Our results on two-bodies are also generalized to restricted three-body, many-body problems, rings and shells.

preprint2010arXiv

Offset between dark matter and ordinary matter: evidence from a sample of 38 lensing clusters of galaxies

We compile a sample of 38 galaxy clusters which have both X-ray and strong lensing observations, and study for each cluster the projected offset between the dominant component of baryonic matter center (measured by X-rays) and the gravitational center (measured by strong lensing). Among the total sample, 45% clusters have offsets >10". The >10" separations are significant, considering the arcsecond precision in the measurement of the lensing/X-ray centers. This suggests that it might be a common phenomenon in unrelaxed galaxy clusters that gravitational field is separated spatially from the dominant component of baryonic matter. It also has consequences for lensing models of unrelaxed clusters since the gas mass distribution may differ from the dark matter distribution and give perturbations to the modeling. Such offsets can be used as a statistical tool for comparison with the results of Lambda-CDM simulations and to test the modified dynamics.

preprint2010arXiv

Structure Formation by Fifth Force: Power Spectrum from N-Body Simulations

We lay out the framework to numerically study nonlinear structure formation in the context of scalar-field-coupled cold dark matter models (phiCDM models) where the scalar field phi serves as dynamical dark energy. Adopting parameters for the scalar field which yield a realistic CMB spectrum, we generate the initial conditions for our Nbody simulations, which follow the spatial distributions of dark matter and the scalar field, by solving their equations of motion using the multilevel adaptive grid technique. We show that the spatial configuration of the scalar field tracks both the voids and clusters of dark matter. The phiCDM model differs from standard LCDM at small scales with observable modifications of, e.g., the mass function of halos as well as the matter power spectrum. Nevertheless, the predictions of both models for the Hubble expansion and the CMB spectrum are virtually indistinguishable. Hence, galaxy cluster counts and weak lensing observations, which probe structure formation at small scales, are needed to falsify this class of models.

preprint2010arXiv

Structure Formation by the Fifth Force III: Segregation of Baryons and Dark Matter

In this paper we present the results of $N$-body simulations with a scalar field coupled differently to cold dark matter (CDM) and baryons. The scalar field potential and coupling function are chosen such that the scalar field acquires a heavy mass in regions with high CDM density and thus behaves like a chameleon. We focus on how the existence of the scalar field affects the formation of nonlinear large-scale structure, and how the different couplings of the scalar field to baryons and CDM particles lead to different distributions and evolutions for these two matter species, both on large scales and inside virialized halos. As expected, the baryon-CDM segregation increases in regions where the fifth force is strong, and little segregation in dense regions. We also introduce an approximation method to identify the virialized halos in coupled scalar field models which takes into account the scalar field coupling and which is easy to implement numerically. It is find that the chameleon nature of the scalar field makes the internal density profiles of halos dependent on the environment in a very nontrivial way.

preprint2010arXiv

Substructure lensing in galaxy clusters as a constraint on low-mass sterile neutrinos in tensor-vector-scalar theory: The straight arc of Abell 2390

Certain covariant theories of the modified Newtonian dynamics paradigm seem to require an additional hot dark matter (HDM) component - in the form of either heavy ordinary neutrinos or more recently light sterile neutrinos (SNs) with a mass around 11eV - to be relieved of problems ranging from cosmological scales down to intermediate ones relevant for galaxy clusters. Here we suggest using gravitational lensing by galaxy clusters to test such a marriage of neutrino HDM and modified gravity, adopting the framework of tensor-vector-scalar theory (TeVeS). Unlike conventional cold dark matter (CDM), such HDM is subject to strong phase-space constraints, which allows one to check cluster lens models inferred within the modified framework for consistency. Since the considered HDM particles cannot collapse into arbitrarily dense clumps and only form structures well above the galactic scale, systems which indicate the need for dark substructure are of particular interest. As a first example, we study the cluster lens Abell 2390 and its impressive straight arc with the help of numerical simulations. Based on our results, we outline a general and systematic approach to model cluster lenses in TeVeS which significantly reduces the calculation complexity. We further consider a simple bimodal lens configuration, capable of producing the straight arc, to demonstrate our approach. We find that such a model is marginally consistent with the hypothesis of 11eV SNs. Future work including more detailed and realistic lens models may further constrain the necessary SN distribution and help to conclusively assess this point. Cluster lenses could therefore provide an interesting discriminator between CDM and such modified gravity scenarios supplemented by SNs or other choices of HDM.

preprint2009arXiv

A Realistic Cosmology Without Cold Dark Matter

We propose a new framework unifying cold dark matter (CDM) and Modified Newtonian Dynamics (MOND) to solve their respective problems on galactic scales and large scale structure formation. In our framework the dark matter clusters on large scales but not on galactic scales. This environment dependence of the dark matter behaviors is controlled by a vector field, which also produces the MOND effects in galaxies. We find that in this framework only a single mass scale needs to be introduced to produce the phenomena of CDM, MOND and also dark energy.

preprint2009arXiv

On the separation between baryonic and dark matter: evidence for phantom dark matter?

The recent years have seen combined measurements of X-ray and (weak) lensing contours for colliding galaxy clusters such as, for instance, the famous "Bullet" cluster. These observations have revealed offsets in the peaks of the baryonic and (dominant) gravitational matter component of order ~(100-200) kpc. Such discrepancies are difficult to explain using modified theories for gravity other than dark matter. Or are they not? Here we use the concept of "phantom dark matter" that is based upon a Newtonian interpretation of the MONDian gravitational potential. We show that this idea is in fact capable of producing substantial offsets in idealistic density configurations, involving a uniform external field. However, when analysed in a MONDian cosmological framework we deduce that the size (and probablity) of the effect is too small to explain the observed offsets found in the most recent observations, at least in the simplest incarnation of phantom dark matter as discussed here. The lensing centers in merging galaxy clusters are likely very close to the centers of true mass even in a MONDian cosmology. This gives the support to the idea that neutrino-like non-collisional matter might be responsible for the observed offsets of lensing and X-ray peaks.

preprint2009arXiv

Structure Formation by Fifth Force I: N-Body vs. Linear Simulations

We lay out the frameworks to numerically study the structure formation in both linear and nonlinear regimes in general dark-matter-coupled scalar field models, and give an explicit example where the scalar field serves as a dynamical dark energy. Adopting parameters of the scalar field which yield a realistic CMB spectrum, we generate the initial conditions for our N-body simulations, which follow the spatial distributions of the dark matter and the scalar field by solving their equations of motion using the multilevel adaptive grid technique. We show that the spatial configuration of the scalar field tracks well the voids and clusters of dark matter. Indeed, the propagation of scalar degree of freedom effectively acts a fifth force on dark matter particles, whose range and magnitude are determined by the two model parameters (mu, gamma), local dark matter density as well as the background value for the scalar field. The model behaves like the LCDM paradigm on scales relevant to the CMB spectrum, which are well beyond the probe of the local fifth force and thus not significantly affected by the matter-scalar coupling. On scales comparable or shorter than the range of the local fifth force, the fifth force is perfectly parallel to gravity and their strengths have a fixed ratio 2gamma^2 determined by the matter-scalar coupling, provided that the chameleon effect is weak; if on the other hand there is a strong chameleon effect (i.e., the scalar field almost resides at its effective potential minimum everywhere in the space), the fifth force indeed has suppressed effects in high density regions and shows no obvious correlation with gravity, which means that the dark-matter-scalar-field coupling is not simply equivalent to a rescaling of the gravitational constant or the mass of the dark matter particles.

preprint2009arXiv

The Galactic potential and the asymmetric distribution of hypervelocity stars

In recent years several hypervelocity stars (HVSs) have been observed in the halo of our Galaxy. Such HVSs have possibly been ejected from the Galactic center and then propagated in the Galactic potential up to their current position. The recent survey for candidate HVSs show an asymmetry in the kinematics of candidate HVSs (position and velocity vectors), where more outgoing stars than ingoing stars (i.e. positive Galactocentric velocities vs. negative ones) are observed. We show that such kinematic asymmetry, which is likely due to the finite lifetime of the stars and Galactic potential structure, could be used in a novel method to probe and constrain the Galactic potential, identify the stellar type of the stars in the survey and estimate the number of HVSs. Kinematics-independent identification of the stellar types of the stars in such surveys (e.g. spectroscopic identification) could further improve these results. We find that the observed asymmetry between ingoing and outgoing stars favors specific Galactic potential models. It also implies a lower limit of ~54+-8 main sequence HVSs in the survey sample (>=648+-96 in the Galaxy), assuming that all of the main sequence stars in the survey originate from the Galactic center. The other stars in the survey are likely to be hot blue horizontal branch stars born in the halo rather than stars ejected from the Galactic center.

preprint2009arXiv

Universality of galactic surface densities within one dark halo scale-length

It was recently discovered that the mean dark matter surface density within one dark halo scale length - the radius within which the volume density profile of dark matter remains approximately flat - is constant across a wide range of galaxies. This scaling relation holds for galaxies spanning a luminosity range of 14 magnitudes and the whole Hubble sequence. Here we report that the luminous matter surface density is also constant within one scale length of the dark halo. This means that the gravitational acceleration generated by the luminous component in galaxies is always the same at this radius. Although the total luminous-to-dark matter ratio is not constant, within one halo scale length it is constant. Our finding can be interpreted as a close correlation between the enclosed surface densities of luminous and dark matter in galaxies.

preprint2008arXiv

Dark Fluid: Towards a unification of empirical theories of galaxy rotation, Inflation and Dark Energy

Empirical theories of Dark Matter like MOND gravity and of Dark Energy like f(R) gravity were motivated by astronomical data. But could these theories be branches rooted from a more general hence natural framework? Here we propose the natural Lagrangian of such a framework based on simple dimensional analysis and co-variant symmetry requirements, and explore various outcomes in a top-down fashion. Our framework preserves the co-variant formulation of GR, but allows the expanding physical metric be bent by a single new species of Dark Fluid flowing in space-time. Its non-uniform stress tensor and current vector are simply functions of a vector field of variable norm, resembling the 4-vector electromagnetic potential description for the photon fluid, but is dark (e.g., by very early decoupling from the baryon-radiation fluid). The Dark Fluid framework naturally branches into a continuous spectrum of theories with Dark Energy and Dark Matter effects, including the $f(R)$ gravity, TeVeS-like theories, Einstein-Aether and $νΛ$ theories as limiting cases. When the vector field degenerates into a pure Higgs-like scalar field, we obtain the physics for inflaton and quintessence. In this broad setting we emphasize the non-constant dynamical field behind the cosmological constant effect, and highlight plausible corrections beyond the classical MOND predictions. Choices of parameters can be made to pass BBN, PPN, and causality constraints. The Dark Fluid is inspired to unify/simplify the astronomically successful ingredients of previous constructions: the desired effects of inflaton plus quintessence plus Cold DM particle fields or MOND-like scalar field(s) are shown largely achievable by one vector field only.

preprint2008arXiv

Perturbations In A Non-Uniform Dark Energy Fluid: Equations Reveal Effects of Modified Gravity and Dark Matter

We propose a unified single-field description of the galactic Dark Matter and various uniform scalar fields for the inflation and cosmological constant. The two types of effects could originate from a fluid of both spatially and temporally varying Vacuum Energy if the vacuum has an uneven pressure caused by a photon-like vector field (of perhaps an unstable massive boson). We propose a most general Lagrangian with a {\bf N}on-{\bf u}niform Cosmological Constant for this vacuum fluid (dubbed as a Nu-Lambda fluid), working within the framework of Einsteinian gravity. This theory includes a continuous spectrum of plausible dark energy theories and gravity theories, e.g., inflation, quintessence, k-essence, f(R), Generalized Einstein-Aether f(K), MOND, TeVeS, BSTV etc. theories. It also suggests new models such as a certain f(K+R) model, which suggests intriguing corrections to MOND depending of redshift and density. Some specific constructions of the Nu-Lambda fluid (e.g., Zhao's V-$Λ$ model) closely resemble the $Λ$CDM cosmology on large scale, but fit galaxy rotation curves as good as MOND. Perturbed Einstein Equations in a simple $f(K_4)$ model are solvable and show effects of a DM coupled to DE. Incorporating the perturbation equations here into standard simulations for cosmological structure growth offers a chance to falsify examples of the Nu-Lambda theories.

preprint2007arXiv

Galaxy Bulge Formation: Interplay with Dark Matter Halo and Central Supermassive Black Hole

We present a simple scenario where the formation of galactic bulges was regulated by the dark halo gravity and regulated the growth of the central supermassive black hole. Assuming the angular momentum is low, we suggest that bulges form in a runaway collapse due to the "gravothermal" instability once the central gas density or pressure exceeds certain threshold (Xu & Zhao 2007). We emphasize that the threshold is nearly universal, set by the background NFW dark matter gravity $g_{DM} \sim 1.2 \times 10^{-8}{\rm cm} {\rm sec}^{-2}$ in the central cusps of halos. Unlike known thresholds for gradual formation of galaxy disks, we show that the universal "halo-regulated" star formation threshold for spheroids matches the very high star formation rate and star formation efficiency shown in high-redshift observations of central starburst regions. The starburst feedback also builds up a pressure shortly after the collapse. This large pressure could both act outward to halt further infall of gas from larger scale, and act inward to counter the Compton-thick wind launched from the central black hole in an Eddington accretion. Assuming the feedback balancing inward and outward forces, our scenario naturally gives rise to the black hole-bulge relationships observed in the local universe.

preprint2006arXiv

Analysis of galactic tides and stars on CDM microhalos

A special purpose N-body simulation has been built to understand the tidal heating of the smallest dark matter substructures ($10^{-6}\msun$ and 0.01pc) from the grainy potential of the Milky Way due to individual stars in the disk and the bulge. To test the method we first run simulations of single encounters of microhalos with an isolated star, and compare with analytical predictions of the dark particle bound fraction as a function of impact parameter. We then follow the orbits of a set of microhalos in a realistic flattened Milky Way potential. We concentrate on (detectable) microhalos passing near the Sun with a range of pericenter and apocenter. Stellar perturbers near the orbital path of a microhalo would exert stochachstic impulses, which we apply in a Monte Carlo fashion according to the Besancon model for the distribution of stars of different masses and ages in our Galaxy. Also incorporated are the usual pericenter tidal heating and disk-shocking heating. We give a detailed diagnosis of typical microhalos and find microhalos with internal tangential anisotropy are slightly more robust than the ones with radial anisotropy. In addition, the dark particles generally go through of a random walk in velocity space and diffuse out of the microhalos. We show that the typical destruction time scales are strongly correlated with the stellar density averaged along a microhalo's orbit over the age of the stellar disk. We also present the morphology of a microhalo at several epochs which may hold the key to dark matter detections.

preprint2006arXiv

Tidal Destruction of The First Dark Microhalos

We point out that the usual self-similarity in cold dark matter models is broken by encounters with individual normal galactic stars on sub-pc scale. Tidal heating and stripping must have redefined the density and velocity structures of the population of the Earth-mass dark matter halos, which are likely to have been the first bound structures to form in the Universe. The disruption rate depends strongly on {\it galaxy types} and the orbital distribution of the microhalos; in the Milky Way, stochastic radial orbits are destroyed first by stars in the triaxial bulge, microhalos on non-planar retrograde orbits with large pericenters and/or apocenters survive the longest. The final microhalo distribution in the {\it solar neighborhood} is better described as a superposition of filamentry microstreams rather than as a set of discrete spherical clumps in an otherwise homogeneous medium. We discuss its important consequences to our detections of microhalos by direct recoil signal and indirect annihilation signal.

preprint1996arXiv

Analytical Dynamical Models for Double-Power-Law Galactic Nuclei

Motivated by the finding that the observed surface brightness profile of many galactic nuclei are well fit by double-power-laws, we explore a range of spherical self-consistent dynamical models with this light profile. We find that the corresponding deprojected volume density profile, phase space density and line-of-sight velocity distribution of these models are well fit by simple analytic approximations. We illustrate the application of these models by fitting a sample of about 25 galactic nuclei observed by Hubble Space Telescope. We give the derived volume density, phase space density, velocity dispersion and line profile parameters in tables. The models are useful for predicting kinematic properties of these galaxies for comparison with future observations. They can also be easily applied to seed N-body simulations of galactic nuclei with realistic density profiles for studying the evolution of these systems.

preprint1996arXiv

On the Microlensing Optical Depth of the Galactic Bar

The microlensing probability (optical depth $τ$) toward the Galactic center carries information about the mass distribution of the Galactic bulge/bar, so can be used to constrain the very uncertain shape parameters of the bar. We find $tau$ depends on the bar mass, radial profile, angle, axis scale lengths and boxyness by a few simple analytical formulae, which shows: (1) $τ$ is proportional to the mass of the bar, $M$. (2) $τ$ falls along the minor axis with a strong gradient. (3) An oblate bulge can have more optical depth than a triaxial bar if the bar angle $α>45$ degress. (4) $τ$ is the largest if the angle $α$ and the axis ratio $y_0/x_0$ conspires so that $y_0/x_0=\tan α$. (5) At a fixed field on the minor axis but away from the center, boxy bars with a flat density profile tend to give a larger optical depth than ellipsoidal bars with a steep profile. (6) Main sequence sources should have a significantly lower (20-50\% lower) optical depth than red clump giants if main sequence stars are not observed as deep as the bright clump giants. An application to four COBE-constrained models (Dwek et al. 1994) shows most models produce optical depth $2σ$ lower than MACHO and OGLE observed values even with both a massive bar $2.8\times 10^{10} M_\odot$ and a full disk. The high $τ$ argues for a massive ($> 2\times 10^{10}M_\odot$) boxy bar with $y_0/x0\approx \tanα$ and $α<20$ deg and with a flat radial profile up to corotation.

preprint1995arXiv

A Consistent Microlensing Model for the Galactic Bar

We compute a microlensing map for the Galactic bar. The predicted event rate and event duration distribution are consistent with the $55$ events recently reported by the MACHO and OGLE collaborations. Most of the events are due to lensing by stars in the near end of the bar. Lens mass functions with about 30-60\% of lens mass as brown dwarfs are rejected at $2-6σ$ levels. To make our model useful for other workers, we tabulate the bar's optical depth and average event duration (scaled to $1M_\odot$ lenses) on a grid of Galactic coordinates. The distance and the proper motions of the lens and the source are derived from a consistent dynamical model of the stellar bar, which has originally been built to fit data on the stellar light and stellar/gas kinematics of the bar. We explore several alternative models and we find that our standard model best matches observations.

preprint1995arXiv

A Self-Consistent Dynamical Model for the COBE Detected Galactic Bar

A 3D steady state stellar dynamical model for the Galactic bar is constructed with 485 orbit building blocks using an extension of Schwarzschild technique. The weights of the orbits are assigned using non-negative least square method. The model fits the density profile of the COBE light distribution, the observed solid body stellar rotation curve, the fall-off of minor axis velocity dispersion and the velocity ellipsoid at Baade's window. We show that the model is stable. Maps and tables of observable velocity moments are made for easy comparisons with observation. The model can also be used to set up equilibrium initial conditions for N-body simulations to study stability. The technique used here can be applied to interpret high quality velocity data of external bulges/bars and galactic nuclei.

preprint1995arXiv

Analytical Models For Galactic Nuclei

I present a general family of dynamical models with simple analytical potential-density pairs suited to model galactic bulges and nuclei with double power-law radial density profiles and an optional central black hole. Analytical expressions for the potential and velocity dispersion of these models are given in the spherical case. Many previously known analytical spherical models, including also the recent $γ/ η$-models by Dehnen (1993) and Tremaine et al. (1994), are special cases of this family. This family also forms a complete set for constructing general galaxy models or solving Poisson's equation in the non-spherical case. In particular, a generalized Clutton-Brock (1973) and Hernquist-Ostriker (1992) orthogonal basis set is given.

preprint1995arXiv

New Parametrizations of Non-Gaussian Line-of-sight Velocity Distribution

A five-parameter fitting formula for the line-of-sight stellar velocity distributions of steady state systems is proposed. It can faithfully reproduce velocity distributions of theoretical models ranging from nearly Gaussian profiles to strongly skewed or mildly double-peaked profiles. In contrast to van der Marel and Franx (1993) and Kuijken and Merrifield (1993), the line profiles are required to have neither multi-peaks nor negative velocity wings. Information of the profile is mostly specified by five physically meaningful and nearly orthogonal fitting parameters.