Accurate background estimation is essential for spectral and temporal analysis in astrophysics. In this work, we construct the in-orbit background model for the High-Energy Telescope (HE) of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). Based on the two-year blank sky observations of Insight-HXMT/HE, we first investigate the basic properties of the background and find that both the background spectral shape and intensity have long-term evolution at different geographical sites. The entire earth globe is then divided into small grids, each with a typical area of 5x5 square degrees in geographical coordinate system. For each grid, an empirical function is used to describe the long-term evolution of each channel of the background spectrum; the intensity of the background can be variable and a modification factor is introduced to account for this variability by measuring the contemporary flux of the blind detector. For a given pointing observation, the background model is accomplished by integrating over the grids that are passed by the track of the satellite in each orbit. Such a background model is tested with both the blank sky observations and campaigns for observations of a series of celestial sources. The results show an average systematic error of 1.5% for the background energy spectrum (26-100 keV) under a typical exposure of 8 ks, and <3% for background light curve estimation (30-150 keV). Therefore, the background model introduced in this paper is included in the Insight-HXMT software as a standard part specialized for both spectral and temporal analyses.

Fast radio bursts (FRBs) are short pulses observed in radio band from cosmological distances, some of which emit repeating bursts. The physical origins of these mysterious events have been subject to wide speculations and heated debates. One class of models invoke soft gamma-ray repeaters (SGRs), or magnetars, as the sources of FRBs. Magnetars are rotating neutron stars with extremely strong magnetic field and can sporadically emit bursts from X-ray ($\sim$keV) to soft gamma-ray ($\sim$sub-MeV) with duration from $10^{-2}$~s to $10^2$~s. However, even though some bright radio bursts have been observed from some magnetars, no FRB-like events had been detected to be associated with any magnetar burst, including one giant flare, and no radio burst has been associated with any X-ray event from any magnetar. Therefore, there is still no observational evidence for magnetar-FRB association up to today. Recently, a pair of FRB-like bursts (FRB~200428 hereafter) separated by 30 milliseconds (ms) were detected from the general direction of the Galactic magnetar SGR~J1935+2154. Here we report the detection of a non-thermal X-ray burst in the 1-250 keV energy band with the Insight-HXMT satellite, which we identify as emitted from SGR~J1935+2154. The burst showed two hard peaks with a separation of $\sim30$ ms, consistent with the separation between the two bursts in FRB~200428. The delay time between the double radio and X-ray peaks is $\sim8.57$ s, fully consistent with the dispersion delay of FRB~200428. We thus identify the non-thermal X-ray burst is associated with FRB~200428 whose high energy counterpart is the two hard peaks in X-ray. Our results suggest that the non-thermal X-ray burst and FRB~200428 share the same physical origin in an explosive event from SGR~J1935+2154.

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.

Distance measurements of supernova remnants (SNRs) are essential and important. Accurate estimates of physical size, dust masses, and some other properties of SNRs depend critically on accurate distance measurements. However, the determination of SNR distances is still a tough task. Red clump stars (RCs) have a long history been used as standard candles. In this work, we take RCs as tracers to determine the distances to a large group of SNRs in the inner disk. We first select RC stars based on the near-infrared (IR) color-magnitude diagram (CMD). Then, the distance to and extinction of RC stars are calculated. To extend the measurable range of distance, we combine near-IR photometric data from the 2MASS survey with the deeper UKIDSS and VVV surveys. With the help of the Gaia parallaxes, we also remove contaminants including dwarfs and giants. Because an SN explosion compresses the surrounding interstellar medium, the SNR region would become denser and exhibit higher extinction than the surroundings. The distance of a SNR is then recognized by the position where the extinction and its gradient is higher than that of the ambient medium. A total of 63 SNRs' distances in the Galactic inner disk are determined and divided into three Levels A, B, and C with decreasing reliability. The distances to 43 SNRs are well determined with reliability A or B. The diameters and dust masses of SNRs are estimated with the obtained distance and extinction.

We use the In data collected during the 2019 outburst from X-ray pulsar 4U 1901+03 to complement the orbital parameters reported by Fermi/GBM. Using the Insight-HXMT, we examine the correlation between the derivative of the intrinsic spin frequency and bolometric flux based on accretion torque models. It was found that the pulse profiles significantly evolve during the outburst. The existence of two types of the profile's pattern discovered in the Insight-HXMT data indicates that this source experienced transition between a super-critical and a sub-critical accretion regime during its 2019 outburst. Based on the evolution of the pulse profiles and the torque model, we derive the distance to 4U 1901+03 as 12.4+-0.2 kpc.

We study the evolution of the temporal properties of MAXI 1820+070 during the 2018 outburst in its hard state from MJD 58190 to 58289 with Insight-HXMT in a broad energy band 1-150 keV. We find different behaviors of the hardness ratio, the fractional rms and time lag before and after MJD 58257, suggesting a transition occurred at around this point. The observed time lags between the soft photons in the 1-5 keV band and the hard photons in higher energy bands, up to 150 keV, are frequency-dependent: the time lags in the low-frequency range, 2-10 mHz, are both soft and hard lags with a timescale of dozens of seconds but without a clear trend along the outburst; the time lags in the high-frequency range, 1-10 Hz, are only hard lags with a timescale of tens of milliseconds; first increase until around MJD 58257 and decrease after this date. The high-frequency time lags are significantly correlated to the photon index derived from the fit to the quasi-simultaneous NICER spectrum in the 1-10 keV band. This result is qualitatively consistent with a model in which the high-frequency time lags are produced by Comptonization in a jet.

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.

The standard $\Lambda$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 $\nu$. 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 $\Phi = - ~ \frac{GM \nu_{ext}}{r}\left(1 + \frac{K_0}{2} \sin^2 \theta \right)$, where $r$ is distance from the mass $M$ which is in an external field that `saturates' the $\nu$ function at the value $\nu_{ext}$, leading to a fixed value of $K_0 \equiv \frac{\partial Ln ~ \nu}{\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 $\theta$ 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.

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 $\alpha$ Centauri. This requires an astrometric accuracy of $\approx 1 \, \mu$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.

With more than 150 blank sky observations at high Galactic latitude, we make a systematic study to the background of the Low Energy Telescope (LE) of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). Both the on-ground simulation and the in-orbit observation indicate that the background spectrum mainly has two components. One is the particle background that dominates above 7 keV and its spectral shape is consistent in every geographical locations. Another is the diffuse X-ray background that dominates below 7 keV and has a stable spectrum less dependent of the sky region. The particle background spectral shape can be obtained from the blind detector data of all the blank sky observations, and the particle background intensity can be measured by the blind detector at 10-12.5 keV. The diffuse X-ray background in the high Galactic latitude can also be obtained from the blank sky spectra after subtracting the particle background. Based on these characteristics, we develop the background model for both the spectrum and the light curve. The systematic error for the background spectrum is investigated with different exposures (T_exp). For the spectrum with T_exp=1 ks, the average systematic errors in 1-7 keV and 1-10 keV are 4.2% and 3.7%, respectively. We also perform the systematic error analyses of the background light curves with different energy bands and time bins. The results show that the systematic errors for the light curves with different time bins are <8% in 1-10 keV.

Glitches correspond to sudden jumps of rotation frequency ($\nu$) and its derivative ($\dot{\nu}$) of pulsars, the origin of which remains not well understood yet, partly because the jump processes of most glitches are not well time-resolved. There are three large glitches of the Crab pulsar, detected in 1989, 1996 and 2017, which were found to have delayed spin-up processes before the normal recovery processes. Here we report two additional glitches of the Crab pulsar occurred in 2004 and 2011 for which we discovered delayed spin up processes, and present refined parameters of the largest glitch occurred in 2017. The initial rising time of the glitch is determined as $<0.48$ hour. We also carried out a statistical study of these five glitches with observed spin-up processes. The two glitches occurred in 2004 and 2011 have delayed spin-up time scales ($\tau_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $\Delta{\nu}$ vs. $|\Delta{\dot\nu}|$ relation of these five glitches is similar to those with no detected delayed spin-up process, indicating that they are similar to the others in nature except that they have larger amplitudes. For these five glitches, the amplitudes of the delayed spin-up process ($|\Delta{\nu}_{\rm d1}|$) and recovery process ($\Delta{\nu}_{\rm d2}$), their time scales ($\tau_{1}$, $\tau_{2}$), and permanent changes in spin frequency ($\Delta{\nu}_{\rm p}$) and total frequency step ($\Delta{\nu}_{\rm g}$) have positive correlations. From these correlations, we suggest that the delayed spin-up processes are common for all glitches, but are too short and thus difficult to be detected for most glitches.

We present the calibration of the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) X-ray satellite, which can be used to perform timing and spectral studies of bright X-ray sources. Insight-HXMT carries three main payloads onboard: the High Energy X-ray telescope (HE), the Medium Energy X-ray telescope (ME) and the Low Energy X-ray telescope (LE). In orbit, the radioactive sources, activated lines, the fluorescence lines and celestial sources are used to calibrate the energy scale and energy resolution of the payloads. The Crab nebular is adopted as the primary effective area calibrator and empirical functions are constructed to modify the simulated effective areas of the three payloads respectively. The systematic errors of HE, compared to the model of the Crab nebular, are less than 2% in 28--120 keV and 2%--10% above 120 keV. The systematic errors of ME are less than 1.5% in 10--35 keV. The systematic errors of LE are less than 1% in 1--7 keV except the Si K--edge (1.839 keV, up to 1.5%) and less than 2% in 7--10 keV.

The Medium Energy X-ray Telescope (ME) is one of the main payloads of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). The background of Insight-HXMT/ME is mainly caused by the environmental charged particles and the background intensity is modulated remarkably by the geomagnetic field, as well as the geographical location. At the same geographical location, the background spectral shape is stable but the intensity varies with the level of the environmental charged particles. In this paper, we develop a model to estimate the ME background based on the ME database that is established with the two-year blank sky observations of the high Galactic latitude. In this model, the entire geographical area covered by Insight-HXMT is divided into grids of $5^{\circ}\times5^{\circ}$ in geographical coordinate system. For each grid, the background spectral shape can be obtained from the background database and the intensity can be corrected by the contemporary count rate of the blind FOV detectors. Thus the background spectrum can be obtained by accumulating the background of all the grids passed by Insight-HXMT during the effective observational time. The model test with the blank sky observations shows that the systematic error of the background estimation in $8.9-44.0$ keV is $\sim1.3\%$ for a pointing observation with an average exposure $\sim5.5$ ks. We also find that the systematic error is anti-correlated with the exposure, which indicates the systematic error is partly contributed by the statistical error of count rate measured by the blind FOV detectors.

We report on our analysis of the 2019 outburst of the X-ray accreting pulsar 4U 1901+03 observed with Insight-HXMT and NICER. Both spectra and pulse profiles evolve significantly in the decaying phase of the outburst. Dozens of flares are observed throughout the outburst. They are more frequent and brighter at the outburst peak. We find that the flares, which have a duration from tens to hundreds of seconds, are generally brighter than the persistent emission by a factor of $\sim$ 1.5. The pulse profile shape during the flares can be significantly different than that of the persistent emission. In particular, a phase shift is clearly observed in many cases. We interpret these findings as direct evidence of changes of the pulsed beam pattern, due to transitions between the sub- and super-critical accretion regimes on a short time scale. We also observe that at comparable luminosities the flares' pulse profiles are rather similar to those of the persistent emission. This indicates that the accretion on the polar cap of the neutron star is mainly determined by the luminosity, i.e., the mass accretion rate.

Theories of gravity that incorporate new scalar degrees of freedom typically require screening mechanisms to ensure consistency with Solar System tests. One widely-studied mechanism -- the chameleon -- can lead to violations of the equivalence principle (EP), as screened and unscreened objects fall differently. If the stars are screened but the surrounding dark matter is not, this leads to asymmetry between leading and trailing streams. We provide analytic estimates of the magnitude of this effect for realistic Galactic mass distributions. Using a restricted N-body code, we simulate 4 satellites with a range of masses and orbits, together with a variety of strengths of the fifth force and screening levels of the Milky Way and satellite. The ratio of the cumulative number function of stars in the leading and trailing stream as a function of longitude from the satellite is computable from simulations, measurable from the stellar data and can provide a direct test. We forecast constraints for streams at large Galactocentric distances, using the specific example case of Hu-Sawicki gravity. Streams with apocentres between 100 and 200 kpc provide attainable constraints at the level of $|f_{R0}| = 10^{-7}$. Still more stringent constraints at the level of $10^{-7.5}$ or even $10^{-8}$ are plausible provided the environmental screening of the satellite is accounted for. These would be among the tightest astrophysical constraints to date. We note further signatures of chameleon gravity: (i) the trailing stellar stream may become detached from the dark matter progenitor if all the stars are lost, (ii) in the extreme fifth force regime, striations in the stellar trailing tail may develop, (iii) if the satellite is fully screened, its orbital frequency is lower than that of the associated dark matter, which is preferentially liberated into the leading tidal tail.

By combining the photometric, spectroscopic, and astrometric information of the stars in the sightline of SNRs, the distances to and the extinctions of 32 Galactic supernova remnants (SNRs) are investigated. The stellar atmospheric parameters are from the SDSS$-$DR14$/$APOGEE and LAMOST$-$DR5$/$LEGUE spectroscopic surveys. The multi-band photometry, from optical to infrared, are collected from the {\it Gaia}, APASS, Pan--STARRS1, 2MASS, and {\it WISE} surveys. With the calibrated {\it Gaia} distances of individual stars, the distances to 15 of 32 SNRs are well determined from their produced extinction and association with molecular clouds. The upper limits of distance are derived for 3 SNRs. The color excess ratios $E(g_{\rm P1}-\lambda) / E(g_{\rm P1}-r_{\rm P1})$ of 32 SNRs are calculated, and their variation with wavebands is fitted by a simple dust model. The inferred dust grain size distribution bifurcates: while the graphite grains have comparable size to the average ISM dust, the silicate grains are generally larger. Along the way, the average extinction law from optical to near-infrared of the Milky Way is derived from the 1.3 million star sample and found to agree with the CCM89 law with $R_{\rm V}=3.15$.

A new black hole X-ray binary (BHXRB) MAXI J1535-571 was discovered by MAXI during its outburst in 2017. Using observations taken by the first Chinese X-ray satellite, the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT), we perform a joint spectral analysis (2-150 keV) in both energy and time domains. The energy spectra provide the essential input for probing the intrinsic Quasi-Periodic Oscillation (QPO) fractional rms spectra (FRS). Our results show that during the intermediate state, the energy spectra are in general consistent with those reported by Swift/XRT and NuSTAR. However, the QPO FRS become harder and the FRS residuals may suggest the presence of either an additional power-law component in the energy spectrum or a turn-over in the intrinsic QPO FRS at high energies.

We study the 2018 outburst of Aql X-1 via the monitor of all sky X-ray image (MAXI) data. We show that the outburst starting in February 2018 is a member of short-low class in the frame of outburst duration and the peak count rate although the outburst morphology is slightly different from the other fast-rise-exponential-decay (FRED) type outbursts with a milder rising stage. We study the partial accretion in the weak propeller stage of Aql X-1 via the MAXI data of the 2018 outburst. We report on the spectral analysis of 3 observations of Aquila X-1 obtained by Insight - hard X-ray modulation telescope (Insight-HXMT) during the late decay stage of the 2018 outburst. We discuss that the data taken by Insight-HXMT is just after the transition to the weak propeller stage. Our analysis shows the necessity of a comptonization component to take into account the existence of an electron cloud resulting photons partly up-scattered.

Impact craters, as "lunar fossils", are the most dominant lunar surface features and occupy most of the Moon's surface. Their formation and evolution record the history of the Solar System. Sixty years of triumphs in the lunar exploration projects accumulated a large amount of lunar data. Currently, there are 9137 existing recognized craters. However, only 1675 of them have been determined age, which is obviously not satisfactory to reveal the evolution of the Moon. Identifying craters is a challenging task due to their enormous difference in size, large variations in shape and vast presence. Furthermore, estimating the age of craters is extraordinarily difficult due to their complex and different morphologies. Here, in order to effectively identify craters and estimate their age, we convert the crater identification problem into a target detection task and crater age estimation into a taxonomy structure. From an initial small number of available craters, we progressively identify craters and estimate their age from Chang'E data by transfer learning (TL) using deep neural networks. For comprehensive identification of multi-scale craters, a two-stage craters detection approach is developed. Thus 117240 unrecognized lunar craters that range in diameter from 532 km to 1 km are identified. Then, a two-stage classification approach is developed to estimate the age of craters by simultaneously extracting their morphological features and stratigraphic information. The age of 79243 craters larger than 3 km in diameter is estimated. These identified and aged craters throughout the mid and low-latitude regions of the Moon are crucial for reconstructing the dynamic evolution process of the Solar System.

As China's first X-ray astronomical satellite, the Hard X-ray Modulation Telescope (HXMT), which was dubbed as Insight-HXMT after the launch on June 15, 2017, is a wide-band (1-250 keV) slat-collimator-based X-ray astronomy satellite with the capability of all-sky monitoring in 0.2-3 MeV. It was designed to perform pointing, scanning and gamma-ray burst (GRB) observations and, based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we give an overview of the mission and its progresses, including payload, core sciences, ground calibration/facility, ground segment, data archive, software, in-orbit performance, calibration, background model, observations and some preliminary results.

We present a detailed timing study of the brightest persistent X-ray source Sco X-1 using the data collected by the Hard X-ray Modulation Telescope ($Insight$-HXMT) from July 2017 to August 2018. A complete $Z$-track hardness-intensity diagram (HID) is obtained. The normal branch oscillations (NBOs) at $\sim$ 6 Hz in the lower part of the normal branch (NB) and the flare branch oscillations (FBOs) at $\sim$ 16 Hz in the beginning part of the flaring branch (FB) are found in observations with the Low Energy X-ray Telescope (LE) and the Medium Energy X-ray Telescope (ME) of $Insight$-HXMT, while the horizontal branch oscillations (HBOs) at $\sim$ 40 Hz and the kilohertz quasi-periodic oscillations (kHz QPOs) at $\sim$ 800 Hz are found simultaneously up to 60 keV for the first time on the horizontal branch (HB) by the High Energy X-ray Telescope (HE) and ME. We find that for all types of the observed QPOs, the centroid frequencies are independent of energy, while the root mean square (rms) increases with energy; the centroid frequencies of both the HBOs and kHz QPOs increase along the $Z$-track from the top to the bottom of the HB; and the NBOs show soft phase lags increasing with energy. A continuous QPO transition from the FB to NB in $\sim$ 200 s are also detected. Our results indicate that the non-thermal emission is the origin of all types of QPOs, the innermost region of the accretion disk is non-thermal in nature, and the corona is nonhomogeneous geometrically.

The evidences for the influence of thermonuclear (type-I) X-ray bursts upon the surrounding environments in neutron star low-mass X-ray binaries (LMXB) were detected previously via spectral and timing analyses. Benefitting from a broad energy coverage of Insight-HXMT, we analyze one photospheric radius expansion (PRE) burst, and find an emission excess at soft X-rays. Our spectral analysis shows that, such an excess is not likely relevant to the disk reflection induced by the burst emission and can be attributed to an enhanced pre-burst/persistent emission. We find that the burst and enhanced persistent emissions sum up to exceed Eddington luminosity by $\sim$ 40 percentages. We speculate that the enhanced emission is from a region beyond the PRE radius, or through the Comptonization of the corona.

Corona cooling was detected previously from stacking a series of short type-I bursts occurred during the low/had state of atoll outburst. Type-I bursts are hence regarded as sharp probe to our better understanding on the basic property of the corona. The launch of the first Chinese X-ray satellite Insight-HXMT has large detection area at hard X-rays which provide almost unique chance to move further in this research field. We report the first detection of the corona cooling by Insight-HXMT from single short type-I burst showing up during {\bf flare} of 4U 1636-536. This type-I X-ray burst has a duration of $\sim$13 seconds and hard X-ray shortage is detected with significance 6.2~$\sigma$ in 40-70 keV. A cross-correlation analysis between the lightcurves of soft and hard X-ray band, shows that the corona shortage lag the burst emission by 1.6 $\pm$1.2~s. These results are consistent with those derived previously from stacking a large amount of bursts detected by RXTE/PCA within a series of {\bf flares} of 4U 1636-536. Moreover, the broad bandwidth of Insight-HXMT allows as well for the first time to infer the burst influence upon the continuum spectrum via performing the spectral fitting of the burst, which ends up with the finding that hard X-ray shortage appears at around 40 keV in the continuum spectrum. These results suggest that the evolution of the corona along with the outburst{\bf /flare} of NS XRB may be traced via looking into a series of embedded type-I bursts by using Insight-HXMT.

The Hard X-ray Modulation Telescope HXMT or also dubbed as Insight HXMT is China s first astronomical satellite . It was launched on 15 th June 2017 in JiuQuan, China and is currently in service smoothly. It was designed to perform point ing , scan ning and gamma ray burst (GRB) observations and , based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we introduce the mission and its progresses in aspects of payload, core sciences, ground calibration/facility , ground segment, data archive, software, in orbit performance, calibration, background model, observations and preliminary results .

The Hard X-ray Modulation Telescope (HXMT) named Insight is China's first X-ray astronomical satellite. The Low Energy X-ray Telescope (LE) is one of its main payloads onboard. The detectors of LE adopt swept charge device CCD236 with L-shaped transfer electrodes. Charges in detection area are read out continuously along specific paths, which leads to a time response distribution of photons readout time. We designed a long exposure readout mode to measure the time response distribution. In this mode, CCD236 firstly performs exposure without readout, then all charges generated in preceding exposure phase are read out completely. Through analysis of the photons readout time in this mode, we obtained the probability distribution of photons readout time.

We present a study of timing properties of the accreting pulsar 2S 1417-624 observed during its 2018 outburst, based on Swift/BAT, Fermi/GBM, Insight-HXMT and NICER observations. We report a dramatic change of the pulse profiles with luminosity. The morphology of the profile in the range 0.2-10.0keV switches from double to triple peaks at $\sim2.5$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$ and from triple to quadruple peaks at $\sim7$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$. The profile at high energies (25-100keV) shows significant evolutions as well. We explain this phenomenon according to existing theoretical models. We argue that the first change is related to the transition from the sub to the super-critical accretion regime, while the second to the transition of the accretion disc from the gas-dominated to the radiation pressure-dominated state. Considering the spin-up as well due to the accretion torque, this interpretation allows to estimate the magnetic field self-consistently at $\sim7\times 10^{12}$G.

The long-term evolution of the centroid energy of the CRSF in Her X-1 is still a mystery. We report a new measurement from a campaign between {\sl Insight}-HXMT and {\sl NuSTAR} performed in February 2018. Generally, the two satellites show well consistent results of timing and spectral properties. The joint spectral analysis confirms that the previously observed long decay phase has ended, and that the line energy instead keeps constant around 37.5 keV after flux correction.

The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

We report on analysis of observations of the bright transient X-ray pulsar \src obtained during its 2017-2018 giant outburst with Insight-HXMT, \emph{NuSTAR}, and \textit{Swift} observatories. We focus on the discovery of a sharp state transition of the timing and spectral properties of the source at super-Eddington accretion rates, which we associate with the transition of the accretion disk to a radiation pressure dominated (RPD) state, the first ever directly observed for magnetized neutron star. This transition occurs at slightly higher luminosity compared to already reported transition of the source from sub- to super-critical accretion regime associate with onset of an accretion column. We argue that this scenario can only be realized for comparatively weakly magnetized neutron star, not dissimilar to other ultra-luminous X-ray pulsars (ULPs), which accrete at similar rates. Further evidence for this conclusion is provided by the non-detection of the transition to the propeller state in quiescence which strongly implies compact magnetosphere and thus rules out magnetar-like fields.

The Be X-ray pulsar (BeXRP) SXP 59 underwent a giant outburst in 2017 with a peak X-ray luminosity of $1.1\times10^{38}$ erg~s$^{-1}$. We report on the X-ray behaviour of SXP 59 with the XMM--Newton and NuSTAR observations collected at the outburst peak, decay, and the low luminosity states. The pulse profiles are energy dependent, the pulse fraction increases with the photon energy and saturates at $\sim$ 65% above 10 keV. It is difficult to constrain the change in the geometry of emitting region with the limited data. Nevertheless, because the pulse shape generally has a double-peaked profile at high luminosity and a single peak profile at low luminosity, we prefer the scenario that the source transited from the super-critical state to the sub-critical regime. This result would further imply that the neutron star (NS) in SXP 59 has a typical magnetic field. We confirm that the soft excess revealed below 2 keV is dominated by a cool thermal component. On the other hand, the NuSTAR spectra can be described as a combination of the non-thermal component from the accretion column, a hot blackbody emission, and an iron emission line. The temperature of the hot thermal component decreases with time, while its size remains constant ($R \sim 0.6$ km). The existence of the hot blackbody at high luminosity cannot be explained with the present accretion theories for BeXRPs. It means that either more sophisticated spectral models are required to describe the X-ray spectra of luminous BeXRPs, or there is non-dipole magnetic field close to the NS surface.

We present that, instead of establishing the equations of motion, one can model-freely reveal the dynamical properties of a black-box system using a learning machine. Trained only by a segment of time series of a state variable recorded at present parameters values, the dynamics of the learning machine at different training stages can be mapped to the dynamics of the target system along a particular path in its parameter space, following an appropriate training strategy that monotonously decreases the cost function. This path is important, because along that, the primary dynamical properties of the target system will subsequently emerge, in the simple-to-complex order, matching closely the evolution law of certain self-evolved systems in nature. Why such a path can be reproduced is attributed to our training strategy. This particular function of the learning machine opens up a novel way to probe the global dynamical properties of a black-box system without artificially establish the equations of motion, and as such it might have countless applications. As an example, this method is applied to infer what dynamical stages a variable star has experienced and how it will evolve in future, by using the light curve observed presently.

In this work, we report the in-orbit demonstration of X-ray pulsar navigation with Insight-Hard X-ray Modulation Telescope (Insight-HXMT), which was launched on Jun. 15th, 2017. The new pulsar navigation method 'Significance Enhancement of Pulse-profile with Orbit-dynamics' (SEPO) is adopted to determine the orbit with observations of only one pulsar. In this test, the Crab pulsar is chosen and observed by Insight-HXMT from Aug. 31th to Sept. 5th in 2017. Using the 5-day-long observation data, the orbit of Insight-HXMT is determined successfully with the three telescopes onboard - High Energy X-ray Telescope (HE), Medium Energy X-ray Telescope (ME) and Low Energy X-ray Telescope (LE) - respectively. Combining all the data, the position and velocity of the Insight-HXMT are pinpointed to within 10 km (3 sigma) and 10 m/s (3 sigma), respectively.

The recently discovered neutron star transient Swift J0243.6+6124 has been monitored by {\it the Hard X-ray Modulation Telescope} ({\it Insight-\rm HXMT). Based on the obtained data, we investigate the broadband spectrum of the source throughout the outburst. We estimate the broadband flux of the source and search for possible cyclotron line in the broadband spectrum. No evidence of line-like features is, however, found up to $\rm 150~keV$. In the absence of any cyclotron line in its energy spectrum, we estimate the magnetic field of the source based on the observed spin evolution of the neutron star by applying two accretion torque models. In both cases, we get consistent results with $B\rm \sim 10^{13}~G$, $D\rm \sim 6~kpc$ and peak luminosity of $\rm >10^{39}~erg~s^{-1}$ which makes the source the first Galactic ultraluminous X-ray source hosting a neutron star.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

The ultra-diffuse dwarf galaxy NGC 1052-DF2 (DF2) has ten (eleven) measured globular clusters (GCs) with a line-of-sight velocity dispersion of $\sigma=7.8^{+5.2}_{-2.2}\,$km/s ($\sigma=10.6^{+3.9}_{-2.3}\,$km/s). Our conventional statistical analysis of the original ten GCs gives $\sigma=8.0^{+4.3}_{-3.0}\,$km/s. The overall distribution of velocities agrees well with a Gaussian of this width. Due to the non-linear Poisson equation in MOND, a dwarf galaxy has weaker self-gravity when in close proximity to a massive host. This external field effect is investigated using a new analytic formulation and fully self-consistent live $N$-body models in MOND. Our formulation agrees well with that of Famaey and McGaugh (2012). These new simulations confirm our analytic results and suggest that DF2 may be in a deep-freeze state unique to MOND. The correctly calculated MOND velocity dispersion agrees with our inferred dispersion and that of van Dokkum et al. (2018b) if DF2 is within 150 kpc of NGC 1052 and both are 20 Mpc away. The GCs of DF2 are however significantly brighter and larger than normal GCs, a problem which disappears if DF2 is significantly closer to us. A distance of 10-13 Mpc makes DF2 a normal dwarf galaxy even more consistent with MOND and the 13 Mpc distance reported by Trujillo et. al. (2019). We discuss the similar dwarf DF4, finding good agreement with MOND. We also discuss possible massive galaxies near DF2 and DF4 along with their distances and peculiar velocities, noting that NGC 1052 may lie at a distance near 10 Mpc.

A great challenge in present-day physics is to understand whether the observed internal dynamics of galaxies is due to dark matter matter or due to a modification of the law of gravity. Recently, van Dokkum et al. reported that the ultra-diffuse dwarf galaxy NGC1052-DF2 lacks dark matter, and they claimed that this would -- paradoxically -- be problematic for modified gravity theories like Milgromian dynamics (MOND). However, NGC1052-DF2 is not isolated, so that a valid prediction of its internal dynamics in MOND cannot be made without properly accounting for the external gravitational fields from neighbouring galaxies. Including this external field effect following Haghi et al. shows that NGC1052-DF2 is consistent with MOND.

We construct a sample of 70 clusters using data from XMM-Newton and Planck to investigate the $Y_{\rm SZ,Planck}-Y_{\rm SZ, XMM}$ scaling relation and the cool-core influences on the relation. $Y_{\rm SZ,XMM}$ is calculated by accurate de-projected temperature and electron number density profiles derived from XMM-Newton. $Y_{\rm SZ,Planck}$ is the latest Planck data restricted to our precise X-ray size $\theta_{\rm 500}$. To study the cool-core influences on$Y_{\rm SZ,Planck}-Y_{\rm SZ, XMM}$ scaling relation, we apply two criteria, limits of central cooling time and classic mass deposition rate, to distinguish cool-core clusters (CCCs) from non-cool-core clusters (NCCCs). We also use $Y_{\rm SZ,Planck}$ from other papers, which are derived from different methods, to confirm our results. The intercept and slope of the$Y_{\rm SZ,Planck}-Y_{\rm SZ, XMM}$ scaling relation are $A=-0.86 \pm 0.30$, $B=0.83 \pm 0.06$. The intrinsic scatter is $\sigma_{\rm ins}= 0.14 \pm 0.03$. The ratio of \mbox{$Y_{\rm SZ,Planck}/Y_{\rm SZ, XMM}$} is $1.03 \pm 0.05$, which is perfectly agreed with unity. Discrepancies of $Y_{\rm SZ,Planck}-Y_{\rm SZ, XMM}$ scaling relation between CCCs and NCCCs are found in observation. They are independent of cool-core classification criteria and $Y_{\rm SZ,Planck}$ calculation methods, although discrepancies are more significant under the classification criteria of classic mass deposition rate. The intrinsic scatter of CCCs (0.04) is quite small compared to that of NCCCs (0.27). The ratio of $Y_{\rm SZ,Planck}/Y_{\rm SZ, XMM}$ for CCCs is $0.89 \pm 0.05$, suggesting that CCCs' $Y_{\rm SZ,XMM}$ may overestimate SZ signal. By contrast, the ratio of $Y_{\rm SZ,Planck}/Y_{\rm SZ, XMM}$ for NCCCs is $1.14 \pm 0.12$, which indicates that NCCCs' $Y_{\rm SZ,XMM}$ may underestimate SZ signal.

We consider disk stability in the quasi-linear formulation of MOND (QUMOND), the basis for some $N$-body integrators. We generalize the Toomre criterion for the stability of disks to tightly wound, axisymmetric perturbations. We apply this to a family of thin exponential disks with different central surface densities. By numerically calculating their QUMOND rotation curves, we obtain the minimum radial velocity dispersion required for stability against local self-gravitating collapse. MOND correctly predicts much higher rotation speeds in low surface brightness galaxies (LSBs) than does Newtonian dynamics without dark matter. Newtonian models thus require putative very massive halos, whose inert nature implies they would strongly stabilize the disk. MOND also increases the stability of galactic disks, but in contradistinction to Newtonian gravity, this extra stability is limited to a factor of 2. MOND is thus rather more conducive to the formation of bars and spiral arms. Therefore, observation of such features in LSBs could be problematic for Newtonian galaxy models. This could constitute a crucial discriminating test. We quantitatively account for these facts in QUMOND. We also compare numerical QUMOND rotation curves of thin exponential disks to those predicted by two algebraic expressions commonly used to calculate MOND rotation curves. For the choice that best approximates QUMOND, we find the circular velocities agree to within 1.5% beyond $\approx 0.5$ disk scale lengths, regardless of the central surface density. The other expression can underestimate the rotational speed by up to 12.5% at one scale length, though rather less so at larger radii.

The possibility that ultra-diffuse galaxies lacking dark matter has recently stimulated interest to check the validity of Modified Newton Dynamics (MOND) predictions on the scale of such galaxies. It has been shown that the External Field Effect (EFE) induced by the close-by galaxy can suppress the velocity dispersion of these systems, so that they appear almost dark matter free in the Newtonian context. Here, following up on this, we are making a priori predictions for the velocity dispersion of 22 ultra-diffuse galaxies in the nearby Universe. This sample can be used to test MOND and the EFE with future follow-up measurements. We construct a catalog of nearby ultra-diffuse galaxies in galaxy group environments, and set upper and lower limits for the possible velocity dispersion allowed in MOND, taking into account possible variations in the mass-to-light ratio of the dwarf and in the distance to the galaxy group. The prediction for the velocity dispersion is made as a function of the three dimensional separation of the dwarf to its host. In 17 out of 22 cases, the EFE plays a crucial role in the prediction.

A central problem of contemporary physics is whether the law of gravity is non-Newtonian on galaxy scales. Rodrigues et al. argue that Milgromian gravitation, which solves the flat rotation curve problem without the need for dark matter particles, is ruled out at > 10{\sigma} significance. To a large extent, this conclusion relies on galaxies with very uncertain distances and/or nearly edge-on orientations, where dust obscuration often becomes significant. Applying appropriate quality cuts to the data leaves only a handful of outliers to the predictions of Milgromian gravitation according to the analysis of Rodrigues et al., but even these outliers can be explained with Milgromian gravitation.

Collisional ring galaxies (CRGs) are formed through off-center collisions between a target galaxy and an intruder dwarf galaxy. We study the mass distribution and kinematics of the CRGs by tuning the bulge-to-disk mass ratio ($B/D$) for the progenitor; i.e., the target galaxy. We find that the lifetime of the ring correlates with the initial impact velocity vertical to the disk plane (i.e., $v_{z0}$). Three orbits for the collisional galaxy pair, on which clear and asymmetric rings form after collisions, are selected to perform the \textit{N}-body simulations at different values of $B/D$ for the progenitor. It is found that the ring structures are the strongest for the CRGs with small values of $B/D$. The S\'{e}rsic index, $n$, of the central remnant in the target galaxy becomes larger after collision. Moreover, the S\'{e}rsic index of a central remnant strongly correlates with the initial value of $B/D$ for the progenitor. A bulge-less progenitor results in a late-type object in the center of the ring galaxy, whereas a bulge-dominated progenitor leads to an early-type central remnant. Progenitors with $B/D\in [0.1,~0.3]$ (i.e., minor bulges) leave central remnants with $n\approx 4$. These results provide a possible explanation for the formation of a recently observed CRG with an early-type central nucleus, SDSS J1634+2049. In addition, we find that the radial and azimuthal velocity profiles for a ring galaxy are more sensitive to the $B/D$ than the initial relative velocity of the progenitor.

We present the X-ray timing results of the new black hole candidate (BHC) MAXI J1535-571 during its 2017 outburst from Hard X-ray Modulation Telescope (\emph{Insight}-HXMT) observations taken from 2017 September 6 to 23. Following the definitions given by \citet{Belloni2010}, we find that the source exhibits state transitions from Low/Hard state (LHS) to Hard Intermediate state (HIMS) and eventually to Soft Intermediate state (SIMS). Quasi-periodic oscillations (QPOs) are found in the intermediate states, which suggest different types of QPOs. With the large effective area of \emph{Insight}-HXMT at high energies, we are able to present the energy dependence of the QPO amplitude and centroid frequency up to 100 keV which is rarely explored by previous satellites. We also find that the phase lag at the type-C QPOs centroid frequency is negative (soft lags) and strongly correlated with the centroid frequency. By assuming a geometrical origin of type-C QPOs, the source is consistent with being a high inclination system.

We report the photometry of six transits of the hot Jupiter HAT-P-29b obtained from 2013 October to 2015 January. We analyze the new light curves, in combination with the published photometric, and Doppler velocimetric, and spectroscopic measurements, finding an updated orbital ephemeris for the HAT-P-29 system, $T_{\rm C}[0]= 2456170.5494(15)\,[\rm{BJD_{TDB}}]$ and $P=5.723390(13) \,{\rm days}$. It is $17.63\,{\rm s}$ ($4.0\,\sigma$) longer than the previously published value, amounting to errors exceeding $2.5\,\rm{hrs}$ at the time of writing (on UTC 2018 June 1). The measured transit mid-times for HAT-P-29b show no compelling evidence of timing anomalies from a linear model, which rules out the presence of a perturbers with masses greater than $0.6$, $0.7$, $0.5$, and $0.4\,{\rm M_\oplus}$ near the $1:2$, $2:3$, $3:2$, and $2:1$ resonances with HAT-P-29b, respectively.

DAMPE (DArk Matter Particle Explorer) is a spaceborne high-energy cosmic ray and gamma-ray detector, successfully launched in December 2015. It is designed to probe astroparticle physics in the broad energy range from few GeV to 100 TeV. The scientific goals of DAMPE include the identification of possible signatures of Dark Matter annihilation or decay, the study of the origin and propagation mechanisms of cosmic-ray particles, and gamma-ray astronomy. DAMPE consists of four sub-detectors: a plastic scintillator strip detector, a Silicon-Tungsten tracKer-converter (STK), a BGO calorimeter and a neutron detector. The STK is composed of six double layers of single-sided silicon micro-strip detectors interleaved with three layers of tungsten for photon conversions into electron-positron pairs. The STK is a crucial component of DAMPE, allowing to determine the direction of incoming photons, to reconstruct tracks of cosmic rays and to estimate their absolute charge (Z). We present the in-flight performance of the STK based on two years of in-flight DAMPE data, which includes the noise behavior, signal response, thermal and mechanical stability, alignment and position resolution.

DArk Matter Particle Explorer (DAMPE) is a general purpose high energy cosmic ray and gamma ray observatory, aiming to detect high energy electrons and gammas in the energy range 5 GeV to 10 TeV and hundreds of TeV for nuclei. This paper provides a method using machine learning to identify electrons and separate them from gammas,protons,helium and heavy nuclei with the DAMPE data from 2016 January 1 to 2017 June 30, in energy range from 10 to 100 GeV.

Interstellar extinction in ultraviolet is the most severe in comparison with optical and infrared wavebands and a precise determination plays an important role in correctly recovering the ultraviolet brightness and colors of objects. By finding the observed bluest colors at given effective temperature and metallicity range of dwarf stars, stellar intrinsic colors, $C^0_{\rm B,V}$, $C^0_{\rm NUV,B}$, $C^0_{\rm FUV,B}$ and $C^0_{\rm FUV,NUV}$, are derived according to the stellar parameters from the LAMOST spectroscopic survey and photometric results from the $GALEX$ and APASS surveys. With the derived intrinsic colors, the ultraviolet color excesses are calculated for about 25,000 A- and F-type dwarf stars. Analysis of the color excess ratios yields the extinction law related to the $GALEX$ UV bands: $E_{{\rm NUV,B}}$/$E_{{\rm B,V}} = 3.77$, $E_{{\rm FUV,B}}$/$E_{{\rm B,V}} = 3.39$, $E_{{\rm FUV,NUV}}$/$E_{{\rm B,V}} = -0.38$. The results agree very well with previous works in the $NUV$ band and in general with the extinction curve derived by Fitzpatrick (1999) for $R_{\rm V}=3.35$.

3XMM J031820.8-663034, first detected by ROSAT in NGC 1313, is one of a few known transient ultraluminous X-ray sources (ULXs). In this paper, we present decades of X-ray data of this source from ROSAT, XMM-Newton, Chandra and the Neil Gehrels Swift Observatory. We find that its X-ray emission experienced four outbursts since 1992, with a typical recurrent time $\sim$ 1800 days, an outburst duration $\sim 240-300$ days, and a nearly constant peak X-ray luminosity $\sim 1.5 \times 10^{39}$ erg/s. The upper limit of X-ray luminosity at the quiescent state is $\sim 5.6 \times 10^{36}$ erg/s, and the total energy radiated during one outburst is $\sim 10^{46}$ erg. The spectra at the high luminosity states can be described with an absorbed disk black-body, and the disk temperature increases with the X-ray luminosity. We compare its outburst properties with other known transient ULXs including ESO 243-49 HLX-1. As its peak luminosity only marginally puts it in the category of ULXs, we also compare it with normal transient black hole binaries. Our results suggest that the source is powered by an accreting massive stellar-mass black hole, and the outbursts are triggered by the thermal-viscous instability.

We attempt to understand the planes of satellite galaxies orbiting the Milky Way (MW) and M31 in the context of Modified Newtonian Dynamics (MOND), which implies a close MW-M31 flyby occurred ${\approx 8}$ Gyr ago. Using the timing argument, we obtain MW-M31 trajectories consistent with cosmological initial conditions and present observations. We adjust the present M31 proper motion within its uncertainty in order to simulate a range of orbital geometries and closest approach distances. Treating the MW and M31 as point masses, we follow the trajectories of surrounding test particle disks, thereby mapping out the tidal debris distribution. Around each galaxy, the resulting tidal debris tends to cluster around a particular orbital pole. We find some models in which these preferred spin vectors align fairly well with those of the corresponding observed satellite planes. The radial distributions of material in the simulated satellite planes are similar to what we observe. Around the MW, our best-fitting model yields a significant fraction (0.22) of counter-rotating material, perhaps explaining why Sculptor counter-rotates within the MW satellite plane. In contrast, our model yields no counter-rotating material around M31. This is testable with proper motions of M31 satellites. In our best model, the MW disk is thickened by the flyby 7.65 Gyr ago to a root mean square height of 0.75 kpc. This is similar to the observed age and thickness of the Galactic thick disk. Thus, the MW thick disk may have formed together with the MW and M31 satellite planes during a past MW-M31 flyby.

Supernova remnants (SNRs) embody the information of the influence on dust properties by the supernova explosion. Based on the color indexes from the 2MASS photometric survey and the stellar parameters from the SDSS$-$DR12$/$APOGEE and LAMOST$-$DR2$/$LEGUE spectroscopic surveys, the near-infrared extinction law and the distance of the Monoceros SNR are derived together with its nearby two nebulas -- the Rosette Nebula and NGC 2264. The distance is found at the position of the sharp increase of interstellar extinction with distance and the nebular extinction is calculated by subtracting the foreground interstellar extinction. The distance of the Monoceros SNR is determined to be $1.98\,$kpc, larger than previous values. Meanwhile, the distance of the Rosette Nebula is $1.55\,$kpc, generally consistent with previous work. The distance between these two nebulas suggests no interaction between them. The distance of NGC 2264, $1.20\,$kpc, exceeds previous values. The color excess ratio, $E_{\rm JH}/E_{\rm JK_S}$, is 0.657 for the Monoceros SNR, consistent with the average value 0.652 for the Milky Way (Xue et al. 2016). The consistency is resulted from that the SNR material is dominated by interstellar dust rather than the supernova ejecta. $E_{\rm JH}/E_{\rm JK_S}$ equals to 0.658 for the Rosette Nebula, further proving the universality of the near-infrared extinction law.

The Be X-ray pulsar SMC X-3 underwent an extra long and ultraluminous giant outburst from 2016 August to 2017 March. The peak X-ray luminosity is up to $\sim 10^{39}$ erg/s, suggesting a mildly super-Eddington accretion onto the strongly magnetized neutron star. It therefore bridges the gap between the Galactic Be/X-ray binaries ($L_{\rm X}^{\rm peak} \leq 10^{38}$ erg/s) and the ultraluminous X-ray pulsars ($L_{\rm X}^{\rm peak} \geq 10^{40}$ erg/s) found in nearby galaxies. A number of observations were carried out to observe the outburst. In this paper, we perform a comprehensive phase-resolved analysis on the high quality data obtained with the Nustar and XMM-Newton, which were observed at a high and intermediate luminosity levels. In order to get a better understanding on the evolution of the whole extreme burst, we take the Swift results at the low luminosity state into account as well. At the early stage of outburst, the source shows a double-peak pulse profile, the second main peak approaches the first one and merges into the single peak at the low luminosity. The second main peak vanishes beyond 20 keV, and its radiation becomes much softer than that of the first main peak. The line widths of fluorescent iron line vary dramatically with phases, indicating a complicated geometry of accretion flows. In contrast to the case at low luminosity, the pulse fraction increases with the photon energy. The significant small pulse fraction detected below 1 keV can be interpreted as the existence of an additional thermal component located at far away from the central neutron star.

The Xuyi Schmidt Telescope Photometric Survey of the Galactic Anti-center (XSTPS-GAC) is a photometric sky survey that covers nearly 6 000 deg^2 towards Galactic anti-center in g r i bands. Half of its survey field locates on the Galactic Anti-center disk, which makes XSTPS-GAC highly suitable for searching new open clusters in the GAC region. In this paper, we report new open cluster candidates discovered in this survey, as well as properties of these open cluster candidates, such as age, distance and reddening, derived by isochrone fitting in the color-magnitude diagram (CMD). These open cluster candidates are stellar density peaks detected in the star density maps by applying the method from Koposov et al. (2008). Each candidate is inspected on its true color image composed from XSTPS-GAC three band images. Then its CMD is checked, in order to identify whether the central region stars have a clear isochrone-like trend differing from the background stars. The parameters derived from isochrone fitting for these candidates are mainly based on three band photometry of XSTPS-GAC. Meanwhile, when these new candidates are able to be seen clearly on 2MASS, their parameters are also derived based on 2MASS (J-H, J) CMD. Finally, there are 320 known open clusters rediscovered and 24 new open cluster candidates discovered in this work. Further more, the parameters of these new candidates, as well as another 11 known recovered open clusters, are properly determined for the first time.

The DArk Matter Particle Explorer (DAMPE) is one of the four satellites within the Strategic Pioneer Research Program in Space Science of the Chinese Academy of Science (CAS). The Silicon-Tungsten Tracker (STK), which is composed of 768 singled-sided silicon microstrip detectors, is one of the four subdetectors in DAMPE, providing track reconstruction and charge identification for relativistic charged particles. The charge response of DAMPE silicon microstrip detectors is complicated, depending on the incident angle and impact position. A new charge reconstruction algorithm for the DAMPE silicon microstrip detector is introduced in this paper. This algorithm can correct the complicated charge response, and was proved applicable by the ion test beam.

The Kepler-9 system harbors three known transiting planets. The system holds significant interest for several reasons. First, the outer two planets exhibit a period ratio that is close to a 2:1 orbital commensurability, with attendant dynamical consequences. Second, both planets lie in the planetary mass "desert" that is generally associated with the rapid gas agglomeration phase of the core accretion process. Third, there exist attractive prospects for accurately measuring both the sky-projected stellar spin-orbit angles as well as the mutual orbital inclination between the planets in the system. Following the original \textit{Kepler} detection announcement in 2010, the initially reported orbital ephemerides for Kepler-9~b and c have degraded significantly, due to the limited time base-line of observations on which the discovery of the system rested. Here, we report new ground-based photometric observations and extensive dynamical modeling of the system. These efforts allow us to photometrically recover the transit of Kepler-9~b, and thereby greatly improve the predictions for upcoming transit mid-times. Accurate ephemerides of this system are important in order to confidently schedule follow-up observations of this system, for both in-transit Doppler measurements as well as for atmospheric transmission spectra taken during transit.

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

We recently showed that several Local Group (LG) galaxies have much higher radial velocities (RVs) than predicted by a 3D dynamical model of the standard cosmological paradigm. Here, we show that 6 of these 7 galaxies define a thin plane with root mean square thickness of only 101 kpc despite a widest extent of nearly 3 Mpc, much larger than the conventional virial radius of the Milky Way (MW) or M31. This plane passes within ${\sim 70}$ kpc of the MW-M31 barycentre and is oriented so the MW-M31 line is inclined by $16^\circ$ to it. We develop a toy model to constrain the scenario whereby a past MW-M31 flyby in Modified Newtonian Dynamics (MOND) forms tidal dwarf galaxies that settle into the recently discovered planes of satellites around the MW and M31. The scenario is viable only for a particular MW-M31 orbital plane. This roughly coincides with the plane of LG dwarfs with anomalously high RVs. Using a restricted $N$-body simulation of the LG in MOND, we show how the once fast-moving MW and M31 gravitationally slingshot test particles outwards at high speeds. The most distant such particles preferentially lie within the MW-M31 orbital plane, probably because the particles ending up with the highest RVs are those flung out almost parallel to the motion of the perturber. This suggests a dynamical reason for our finding of a similar trend in the real LG, something not easily explained as a chance alignment of galaxies with an isotropic or mildly flattened distribution (probability $= {0.0015}$).

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

Finding the electromagnetic (EM) counterpart of binary compact star merger, especially the binary neutron star (BNS) merger, is critically important for gravitational wave (GW) astronomy, cosmology and fundamental physics. On Aug. 17, 2017, Advanced LIGO and \textit{Fermi}/GBM independently triggered the first BNS merger, GW170817, and its high energy EM counterpart, GRB 170817A, respectively, resulting in a global observation campaign covering gamma-ray, X-ray, UV, optical, IR, radio as well as neutrinos. The High Energy X-ray telescope (HE) onboard \textit{Insight}-HXMT (Hard X-ray Modulation Telescope) is the unique high-energy gamma-ray telescope that monitored the entire GW localization area and especially the optical counterpart (SSS17a/AT2017gfo) with very large collection area ($\sim$1000 cm$^2$) and microsecond time resolution in 0.2-5 MeV. In addition, \textit{Insight}-HXMT quickly implemented a Target of Opportunity (ToO) observation to scan the GW localization area for potential X-ray emission from the GW source. Although it did not detect any significant high energy (0.2-5 MeV) radiation from GW170817, its observation helped to confirm the unexpected weak and soft nature of GRB 170817A. Meanwhile, \textit{Insight}-HXMT/HE provides one of the most stringent constraints (~10$^{-7}$ to 10$^{-6}$ erg/cm$^2$/s) for both GRB170817A and any other possible precursor or extended emissions in 0.2-5 MeV, which help us to better understand the properties of EM radiation from this BNS merger. Therefore the observation of \textit{Insight}-HXMT constitutes an important chapter in the full context of multi-wavelength and multi-messenger observation of this historical GW event.

Understanding the effects of dust extinction is important to properly interpret observations. The optical total-to-selective extinction ratio, Rv = Av/E(B-V), is widely used to describe extinction variations in ultraviolet and optical bands. Since the Rv=3.1 extinction curve adequately represents the average extinction law of diffuse regions in the Milky Way, it is commonly used to correct observational measurements along sightlines toward diffuse regions in the interstellar medium. However, the Rv value may vary even along different diffuse interstellar medium sightlines. In this paper, we investigate the optical--mid-infrared (mid-IR) extinction law toward a very diffuse region at l = 165 in the Galactic plane, which was selected based on a CO emission map. Adopting red clump stars as extinction tracers, we determine the optical-to-mid-IR extinction law for our diffuse region in the two APASS bands (B, V), the three XSTPS-GAC bands (g, r, i), the three 2MASS bands (J, H, Ks, and the two WISE bands (W1, W2). Specifically, 18 red clump stars were selected from the APOGEE--RC catalog based on spectroscopic data in order to explore the diversity of the extinction law. We find that the optical extinction curves exhibit appreciable diversity. The corresponding Rv ranges from 1.7 to 3.8, while the mean Rv value of 2.8 is consistent with the widely adopted average value of 3.1 for Galactic diffuse clouds. There is no apparent correlation between Rv value and color excess E(B-V) in the range of interest, from 0.2 to 0.6 mag, or with specific visual extinction per kiloparsec, AV/d.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

It has recently been proposed, by assuming that dark matter is a superfluid, that MOND-like effects can be achieved on small scales whilst preserving the success of $\Lambda$CDM on large scales. Here we aim to provide the first set of spherical models of galaxy clusters in the context of superfluid dark matter. We first outline the theoretical structure of the superfluid core and the surrounding "normal phase" dark halo of quasi-particles in thermal equlibrium. The latter should encompass the largest part of galaxy clusters. Here, we set the SfDM transition at the radius where the density and pressure of the superfluid and normal phase coincides, neglecting the effect of phonons in the suprefluid core. We then apply the theory to a sample of galaxy clusters, and directly compare the SfDM predicted mass profiles to data. We find that the superfluid formulation can reproduce the X-ray dynamical mass profile of clusters. The SfDM fits however display slight under-predictions of the gravity in the central regions which might be partly related to our neglecting of the effect of phonons in these regions. We conclude that this superfluid formulation is successful in describing galaxy clusters, but further work will be needed to determine whether the parameter choice is consistent with galaxies. Our model could be made more realistic by exploring non-sphericity and the SfDM transition condition we impose.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

We present a detection of 89 candidates of ultra-diffuse galaxies (UDGs) in a 4.9 degree$^2$ field centered on the Hickson Compact Group 95 (HCG 95) using deep $g$- and $r$-band images taken with the Chinese Near Object Survey Telescope. This field contains one rich galaxy cluster (Abell 2588 at $z$=0.199) and two poor clusters (Pegasus I at $z$=0.013 and Pegasus II at $z$=0.040). The 89 candidates are likely associated with the two poor clusters, giving about 50 $-$ 60 true UDGs with a half-light radius $r_{\rm e} > 1.5$ kpc and a central surface brightness $\mu(g,0) > 24.0$ mag arcsec$^{-2}$. Deep $z$'-band images are available for 84 of the 89 galaxies from the Dark Energy Camera Legacy Survey (DECaLS), confirming that these galaxies have an extremely low central surface brightness. Moreover, our UDG candidates are spread over a wide range in $g-r$ color, and $\sim$26% are as blue as normal star-forming galaxies, which is suggestive of young UDGs that are still in formation. Interestingly, we find that one UDG linked with HCG 95 is a gas-rich galaxy with H I mass $1.1 \times 10^{9} M_{\odot}$ detected by the Very Large Array, and has a stellar mass of $M_\star \sim 1.8 \times 10^{8}$ $M_{\odot}$. This indicates that UDGs at least partially overlap with the population of nearly dark galaxies found in deep H I surveys. Our results show that the high abundance of blue UDGs in the HCG 95 field is favored by the environment of poor galaxy clusters residing in H I-rich large-scale structures.

Adopting Schwarzschild's orbit-superposition technique, we construct a series of self-consistent galaxy models, embedded in the external field of galaxy clusters in the framework of Milgrom's MOdified Newtonian Dynamics. These models represent relatively massive ellipticals with a Hernquist radial profile at various distances from the cluster centre. Using $N$-body simulations, we perform a first analysis of these models and their evolution. We find that self-gravitating axisymmetric density models, even under a weak external field, lose their symmetry by instability and generally evolve to triaxial configurations. A kinematic analysis suggests that the instability originates from both box and non-classified orbits with low angular momentum. We also consider a self-consistent isolated system which is then placed in a strong external field and allowed to evolve freely. This model, just as the corresponding equilibrium model in the same external field, eventually settles to a triaxial equilibrium as well, but has a higher velocity radial anisotropy and is rounder. The presence of an external field in MOND universe generically predicts some lopsidedness of galaxy shapes.

The Be X-ray pulsar, SMC X-3 underwent a giant outburst from 2016 August to 2017 March, which was monitored with the Swift satellite. During the outburst, its broadband flux increased dramatically, and the unabsorbed X-ray luminosity reached an extreme value of $\sim 10^{39}$ erg/s around August 24. Using the Swift/XRT data, we measure the observed pulse frequency of the neutron star to compute the orbital parameters of the binary system. After applying the orbital corrections to Swift observations, we find that the spin frequency increases steadily from 128.02 mHz on August 10 and approach to spin equilibrium $\sim 128.74$ mHz in 2017 January with an unabsorbed luminosity of $L_{\rm X} \sim 2\times10^{37}$ erg/s, indicating a strong dipolar magnetic field $B \sim 6.8\times10^{12}$ G at the neutron star surface. The spin-up rate is tightly correlated with its X-ray luminosity during the super-Eddington outburst. The pulse profile in the Swift/XRT data is variable, showing double peaks at the early stage of outburst and then merging into a single peak at low luminosity. Additionally, we report that a low temperature ($kT \sim 0.2$ keV) thermal component emerges in the phase-averaged spectra as the flux decays, and it may be produced from the outer truncated disk or the boundary layer between the exterior flow and the magnetosphere.

The DArk Matter Particle Explorer (DAMPE) is one of the four satellites within Strategic Pioneer Research Program in Space Science of the Chinese Academy of Science (CAS). DAMPE can detect electrons, photons in a wide energy range (5 GeV to 10 TeV) and ions up to iron (100GeV to 100 TeV). Silicon-Tungsten Tracker (STK) is one of the four subdetectors in DAMPE, providing photon-electron conversion, track reconstruction and charge identification for ions. Ion beam test was carried out in CERN with 60GeV/u Lead primary beams. Charge reconstruction and charge resolution of STK detectors were investigated.

The Type IIn supernova (SN) 1978K in the nearby galaxy NGC 1313 has remained bright in multiwavelengths for more than $\sim 25$ yr. The archival data of SN 1978K collected with ROSAT, ASCA, XMM-Newton, and Chandra from 1990 to 2006 show no significant variation of the soft X-ray emission but a hint of flux decrease in hard X-rays. In this work, we perform a detailed analysis using more than 15 yr of XMM-Newton observations. Both the 0.5-2 keV and 2-10 keV light-curves decline as $t^{-1}$ from 2000 to 2015. The transition of light curve profiles can be explained in a way that the reverse shock was radiative at an early phase and then became adiabatic at late times. Such a scenario is also supported by the spectral analysis results. We also found a decrease in the absorption column density, which indicates the presence of a cool shell during the radiative phase.

The fraction of binary stars (fb) is one of most valuable tool to probe the star formation and evolution of multiple systems in the Galaxy. We focus on the relationship between fb and stellar metallicity ([Fe/H]) by employing the differential radial velocity (DRV) method and the large sample observed by the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST). Main-sequence stars from A- to K-types in the third data release (DR3) of LAMOST are selected to estimate fb. Contributions to a profile of DRV from radial velocity (RV) error of single stars (\sigma) and orbital motion of binary stars are evaluated from the profile of DRV. Finally, we employ 365,911 stars with randomly repeating spectral observations to present a detailed analysis of fb and \sigma in the two-dimensional (2D) space of Teff and [Fe/H]. The A-type stars are more likely to be companions in binary star systems than other stars. Furthermore, the reverse correlation between fb and [Fe/H] can be shown statistically, which suggests that fb is a joint function of Teff and [Fe/H]. At the same time, \sigma of the sample for different Teff and [Fe/H] are fitted. Metal-rich cold stars in our sample have the best RV measurement.

We model the cluster A1689 in two modified MOND frameworks (EMOND and what we call generalised MOND or GMOND) with the aim of determining whether it is possible to explain the inferred acceleration profile, from gravitational lensing, without the aid of dark matter. We also compare our result to predictions from MOG/STVG and Emergent Gravity. By using a baryonic mass model, we determine the total gravitational acceleration predicted by the modified gravitational equations and compare the result to NFW profiles of dark matter studies on A1689 from the literature. Theory parameters are inferred empirically, with the aid of previous work. We are able to reproduce the desired acceleration profile of A1689 for GMOND, EMOND and MOG, but not Emergent Gravity. There is much more work which needs to be conducted with regards to understanding how the GMOND parameters behave in different environments. Furthermore, we show that the exact baryonic profile becomes very important when undertaking modified gravity modelling rather than {\Lambda}CDM.

A sample of Coma cluster ultra-diffuse galaxies (UDGs) are modelled in the context of Extended Modified Newtonian Dynamics (EMOND) with the aim to explain the large dark matter-like effect observed in these cluster galaxies. We first build a model of the Coma cluster in the context of EMOND using gas and galaxy mass profiles from the literature. Then assuming the dynamical mass of the UDGs satisfies the fundamental manifold of other ellipticals, and that the UDG stellar mass-to-light matches their colour, we can verify the EMOND formulation by comparing two predictions of the baryonic mass of UDGs. We find that EMOND can explain the UDG mass, within the expected modelling errors, if they lie on the fundamental manifold of ellipsoids, however, given that measurements show one UDG lying off the fundamental manifold, observations of more UDGs are needed to confirm this assumption.

he DArk Matter Particle Explorer (DAMPE) is a general purposed satellite-borne high energy $\gamma-$ray and cosmic ray detector, and among the scientific objectives of DAMPE are the searches for the origin of cosmic rays and an understanding of Dark Matter particles. As one of the four detectors in DAMPE, the Plastic Scintillator Detector (PSD) plays an important role in the particle charge measurement and the photons/electrons separation. The PSD has 82 modules, each consists of a long organic plastic scintillator bar and two PMTs at both ends for readout, in two layers and covers an overall active area larger than 82 cm $\times$ 82 cm. It can identify the charge states for relativistic ions from H to Fe, and the detector efficiency for Z=1 particles can reach 0.9999. The PSD has been successfully launched with DAMPE on Dec. 17, 2015. In this paper, the design, the assembly, the qualification tests of the PSD and some of the performance measured on the ground have been described in detail.

We attempt to fit the observed radial velocities (RVs) of ${\,{\sim}\,}30$ Local Group (LG) galaxies using a 3D dynamical model of it and its immediate environment within the context of the standard cosmological paradigm, $\Lambda$CDM. This extends and confirms the basic results of our previous axisymmetric investigation of the LG (MNRAS, 459, 2237). We find that there remains a tendency for observed RVs to exceed those predicted by our best-fitting model. The typical mismatch is slightly higher than in our 2D model, with a root mean square value of ${\,{\sim}\,}50$ km/s. Our main finding is that including the 3D distribution of massive perturbing dark matter halos is unlikely to help greatly with the high velocity galaxy problem. Nonetheless, the 2D and 3D results differ in several other ways such as which galaxies' RVs are most problematic and the preferred values of parameters common to both models. The anomalously high RVs of several LG dwarfs may be better explained if the Milky Way (MW) and Andromeda (M31) were once moving much faster than in our models. This would allow LG dwarfs to gain very high RVs via gravitational slingshot encounters with a massive fast-moving galaxy. Such a scenario is possible in some modified gravity theories, especially those which require the MW and M31 to have previously undergone a close flyby. In a $\Lambda$CDM context, however, this scenario is not feasible as the resulting dynamical friction would cause a rapid merger.

The DArk Matter Particle Explorer (DAMPE) is one of the four satellites within Strategic Pioneer Research Program in Space Science of the Chinese Academy of Science (CAS). DAMPE can detect electrons, photons and ions in a wide energy range (5 GeV to 10 TeV) and ions up to iron (100GeV to 100 TeV). Plastic Scintillator Detector (PSD) is one of the four payloads in DAMPE, providing e/{\gamma} separation and charge identification up to Iron. An ion beam test was carried out for the Qualification Model of PSD in CERN with 40GeV/u Argon primary beams. The Birk's saturation and charge resolution of PSD were investigated.

Modified Newtonian Dynamics (MOND) is a gravitational framework designed to explain the astronomical observations in the Universe without the inclusion of particle dark matter. Modified Newtonian Dynamics, in its current form, cannot explain the missing mass in galaxy clusters without the inclusion of some extra mass, be it in the form of neutrinos or non-luminous baryonic matter. We investigate whether the MOND framework can be generalized to account for the missing mass in galaxy clusters by boosting gravity in high gravitational potential regions. We examine and review Extended MOND (EMOND), which was designed to increase the MOND scale acceleration in high potential regions, thereby boosting the gravity in clusters. We seek to investigate galaxy cluster mass profiles in the context of MOND with the primary aim at explaining the missing mass problem fully without the need for dark matter. Using the assumption that the clusters are in hydrostatic equilibrium, we can compute the dynamical mass of each cluster and compare the result to the predicted mass of the EMOND formalism. We find that EMOND has some success in fitting some clusters, but overall has issues when trying to explain the mass deficit fully. We also investigate an empirical relation to solve the cluster problem, which is found by analysing the cluster data and is based on the MOND paradigm. We discuss the limitations in the text.

The spectral and temporal evolution during X-ray outbursts give important clues on the accretion process and radiation mechanism in black-hole X-ray binaries (BH XRBs). A set of Swift and RXTE observations were executed to monitor the 2008 outburst of the black-hole candidate Swift J1842.5-1124. We investigate these data to explore the accretion physics in BH XRBs. We carry out a comprehensive spectral and timing analysis on all the available pointing observations, including fitting both X-ray spectra and power density spectra, measuring the optical and near-ultraviolet flux density. We also search for correlations among the spectral and timing parameters. The observed properties of Swift J1842.5-1124 are similar to other BH XRBs in many respects, for example the hardness-intensity diagram and hardness-rms diagram. The type-C quasi-periodic oscillations (QPOs) were observed as the source started to transit from the low-hard state to the high-soft state. The frequency of QPOs correlate with intensity and the hard component index, and anti-correlate with the hardness and the total fractional rms. These relations are consistent with the Lense-Thirring precession model. The estimated U-band flux changed with the X-ray flux, while the flux density at the V band remained 0.26 mJy. These results imply that the X-ray reprocessing or the tail of thermal emission from the outer disk contributes a significant fraction of the U-band radiation; alternatively, the companion star or the jet dominates the flux at longer wavelengths.

Intrinsic colors of normal stars are derived in the popularly used infrared bands involving the 2MASS/JHKs, WISE, Spitzer/IRAC and AKARI/S9W filters. Based on three spectroscopic surveys -- LAMOST, RAVE and APOGEE, stars are classified into groups of giants and dwarfs, as well as metal-normal and metal-poor stars. An empirical analytical relation of the intrinsic color is obtained with stellar effective temperature (Teff) for each group of stars after the zero-reddening stars are selected from the blue edge in the $J-\lambda$ versus (Teff) diagram. It is found that metallicity has little effect on the infrared colors. In the near-infrared bands, our results agree with previous work. In addition, the color indexes H-W2 and Ks-W1 that are taken as constant to calculate interstellar extinction are discussed. The intrinsic color of M-type stars are derived separately due to lack of accurate measurement of their effective temperature.

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 $\Lambda 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.

The positions and velocities of galaxies in the Local Group (LG) measure the gravitational field within it. This is mostly due to the Milky Way (MW) and Andromeda (M31). We constrain their masses using distance and radial velocity (RV) measurements of 32 LG galaxies. To do this, we follow the trajectories of many simulated particles starting on a pure Hubble flow at redshift 9. For each observed galaxy, we obtain a trajectory which today is at the same position. Its final velocity is the model prediction for the velocity of that galaxy. Unlike previous simulations based on spherical symmetry, ours are axisymmetric and include gravity from Centaurus A. We find the total LG mass is ${4.33^{+0.37}_{-0.32}\times{10}^{12}M_\odot}$, with $0.14 \pm 0.07$ of this being in the MW. We approximately account for IC 342, M81, the Great Attractor and the Large Magellanic Cloud. No plausible set of initial conditions yields a good match to the RVs of our sample of LG galaxies. Observed RVs systematically exceed those predicted by the best-fitting $\Lambda$CDM model, with a typical disagreement of ${45.1^{+7.0}_{-5.7}}$ km/s and a maximum of ${110 \pm 13}$ km/s for DDO 99. Interactions between LG dwarf galaxies can't easily explain this. One possibility is a past close flyby of the MW and M31. This arises in some modified gravity theories but not in $\Lambda$CDM. Gravitational slingshot encounters of material in the LG with either of these massive fast-moving galaxies could plausibly explain why some non-satellite LG galaxies are moving away from us even faster than a pure Hubble flow.

We investigate the X-ray spectrum, variability and the surrounding ionized bubble of NGC1313 X-2 to explore the physics of super-Eddington accretion. Beyond the Eddington luminosity, the accretion disk of NGC1313 X-2 is truncated at a large radius ($\sim$ 50 times of innermost stable circular orbit), and displays the similar evolution track with both luminous Galactic black-hole and neutron star X-ray binaries. In super-critical accretion, the speed of radiatively driven outflows from the inner disk is mildly relativistic. Such ultra-fast outflows would be over ionized and might produce weak Fe K absorption lines, which may be detected by the coming X-ray mission {\it Astro-H}. If the NGC1313 X-2 is a massive stellar X-ray binary, the high luminosity indicates that an ephemeral feast is held in the source. That is, the source must be accreting at a hyper-Eddington mass rate to give the super-Eddington emission over $\sim 10^{4}-10^{5}$ yr. The expansion of the surrounding bubble nebula with a velocity of $\sim$ 100 km~s$^{-1}$ might indicate that it has existed over $\sim 10^{6}$ yr and is inflated by the radiatively driven outflows from the transient with a duty cycle of activity of $\sim$ a few percent. Alternatively, if the surrounding bubble nebula is produced by line-driven winds, less energy is required than the radiatively driven outflow scenario, and the radius of the Str\"{o}mgren radius agrees with the nebula size. Our results are in favor of the line-driven winds scenario, which can avoid the conflict between the short accretion age and the apparently much longer bubble age inferred from the expansion velocity in the nebula.

The Cherenkov Telescope Array (CTA) is an international next-generation ground-based gamma-ray observatory. CTA will be implemented as southern and northern hemisphere arrays of tens of small, medium and large-sized imaging Cherenkov telescopes with the goal of improving the sensitivity over the current-generation experiments by an order of magnitude. CTA will provide energy coverage from ~20 GeV to more than 300 TeV. The Schwarzschild-Couder (SC) medium size (9.5m) telescopes will feature a novel aplanatic two-mirror optical design capable of accommodating a wide field-of-view with significantly improved angular resolution as compared to the traditional Davies-Cotton optical design. A full-scale prototype SC medium size telescope structure has been designed and will be constructed at the Fred Lawrence Whipple Observatory in southern Arizona during the fall of 2015. concentrate on the novel features of the design.

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

This is a brief rebuttal to arXiv:1502.03821, which claims to provide the first observational proof of dark matter interior to the solar circle. We point out that this result is not new, and can be traced back at least a quarter century.

We use a sample of 62 clusters of galaxies to investigate the discrepancies of gas temperature and total mass within r500 between XMM-Newton and Chandra data. Comparisons of the properties show that: (1) Both the de-projected and projected temperatures determined by Chandra are higher than those of XMM-Newton and there is a good linear relation for the de-projected temperature. (2) The Chandra mass is much higher than XMM-Newton mass with a bias of 0.15. To explore the reasons for the discrepancy in mass, we recalculate the Chandra mass (expressed as M_c) by modifying its temperature with the de-projected temperature relation. The results show that M_c is more close to the XMM-Newton mass with the bias reducing to 0.02. Moreover, M_c are corrected with the r500 measured by XMM-Newton and the intrinsic scatter is significantly improved with the value reducing from 0.20 to 0.12. These mean that the temperature bias may be the main factor causing the mass bias. At last, we find that M_c is consistent with the corresponding XMM-Newton mass derived directly from our mass relation at a given Chandra mass. Thus, the de-projected temperature and mass relations can provide unbiased corrections for galaxy cluster properties derived from Chandra and XMM-Newton.

We present a three dimensional extinction map in $r$ band. The map has a spatial angular resolution, depending on latitude, between 3 -- 9\,arcmin and covers the entire XSTPS-GAC survey area of over 6,000\,$\rm deg^2$ for Galactic longitude $\rm 140 \leq$ $l$ $ \leq 220\deg$ and latitude $\rm -40\leq$ $b$ $ \leq 40\deg$. By cross-matching the photometric catalog of the Xuyi Schmidt Telescope Photometric Survey of the Galactic Anticentre (XSTPS-GAC) with those of 2MASS and WISE, we have built a multi-band photometric stellar sample of about 30 million stars and applied spectral energy distribution (SED) fitting to the sample. By combining photometric data from the optical to the near-infrared, we are able to break the degeneracy between the intrinsic stellar colours and the amounts of extinction by dust grains for stars with high photometric accuracy, and trace the extinction as a function of distance for low Galactic latitude and thus highly extincted regions. This has allowed us to derive the best-fit extinction and distance information of more than 13 million stars, which are used to construct the three dimensional extinction map. We have also applied a Rayleigh-Jeans colour excess (RJCE) method to the data using the 2MASS and WISE colour $(H-W2)$. The resulting RJCE extinction map is consistent with the integrated two dimensional map deduced using the best-fit SED algorithm. However for individual stars, the amounts of extinction yielded by the RJCE method suffer from larger errors than those given by the best-fit SED algorithm.

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.} $\beta < \gamma/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.

We present a catalog of RR Lyrae stars (RRLs) observed by the Xuyi Schmidt Telescope Photometric Survey (XDSS). The area we consider is located in the North Galactic Cap, covering 376.75 sq deg at RA $\approx$ 150 deg and Dec $\approx$ 27 deg down to a magnitude limit of i $\approx$ 19. Using the variability information afforded by the multi-epoch nature of our XDSS data, combined with colors from the Sloan Digital Sky Survey, we are able to identify candidate RRLs. We find 318 candidates, derive distances to them and estimate the detection efficiency. The majority of our candidates have more than 12 observations and for these we are able to calculate periods. These also allows us to estimate our contamination level, which we predict is between 30% to 40%. Finally we use the sample to probe the halo density profile in the 9-49 kpc range and find that it can be well fitted by a double power law. We find good agreement between this model and the models derived for the South Galactic Cap using the Watkins et al. (2009) and Sesar et al. (2010) RRL data-sets, after accounting for possible contamination in our data-set from Sagittarius stream members. We consider non-spherical double power law models of the halo density profile and again find agreement with literature data-sets, although we have limited power to constrain the flattening due to our small survey area. Much tighter constraints will be placed by current and future wide-area surveys, most notably ESA's astrometric Gaia mission. Our analysis demonstrates that surveys with a limited number of epochs can effectively be mined for RRLs. Our complete sample is provided as accompanying online material.

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.

On 13 December 2012, Chang'e-2 conducted a successful flyby of the near-Earth asteroid 4179 Toutatis at a closest distance of 770 $\pm$ 120 meters from the asteroid's surface. The highest-resolution image, with a resolution of better than 3 meters, reveals new discoveries on the asteroid, e.g., a giant basin at the big end, a sharply perpendicular silhouette near the neck region, and direct evidence of boulders and regolith, which suggests that Toutatis may bear a rubble-pile structure. Toutatis' maximum physical length and width are (4.75 $\times$ 1.95 km) $\pm$10$\%$, respectively, and the direction of the +$z$ axis is estimated to be (250$\pm$5$^\circ$, 63$\pm$5$^\circ$) with respect to the J2000 ecliptic coordinate system. The bifurcated configuration is indicative of a contact binary origin for Toutatis, which is composed of two lobes (head and body). Chang'e-2 observations have significantly improved our understanding of the characteristics, formation, and evolution of asteroids in general.

Since it can describe the rotation curves of galaxies without dark matter and can give rise to accelerated expansion, conformal gravity attracts much attention recently. As a theory of modified gravity, it is important to test conformal gravity with astrophysical observations. Here we constrain conformal gravity with SNIa and Hubble parameter data and investigate whether it suffers from an age problem with the age of APM~08279+5255. We find conformal gravity can accommodate the age of APM~08279+5255 at 3 $\sigma$ deviation, unlike most of dark energy models which suffer from an age problem.

We use a sample of 39 galaxy clusters at redshift z < 0.1 observed by XMM-Newton to investigate the relations between X-ray observables and total mass. Based on central cooling time and central temperature drop, the clusters in this sample are divided into two groups: 25 cool core clusters and 14 non-cool core clusters, respectively. We study the scaling relations of Lbol-M500, M500-T, M500-Mg and M500-YX, and also the influences of cool core on these relations. The results show that the M500-YX relation has a slope close to the standard self-similar value, has the smallest scatter and does not vary with the cluster sample. Moreover, the M500-YX relation is not affected by the cool core. Thus, the parameter of YX may be the best mass indicator.

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.

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.

Polar ring galaxies are ideal objects with which to study the three-dimensional shapes of galactic gravitational potentials since two rotation curves can be measured in two perpendicular planes. Observational studies have uncovered systematically larger rotation velocities in the extended polar rings than in the associated host galaxies. In the dark matter context, this can only be explained through dark halos that are systematically flattened along the polar rings. Here, we point out that these objects can also be used as very effective tests of gravity theories, such as those based on Milgromian dynamics (MOND). We run a set of polar ring models using both Milgromian and Newtonian dynamics to predict the expected shapes of the rotation curves in both planes, varying the total mass of the system, the mass of the ring with respect to the host, as well as the size of the hole at the center of the ring. We find that Milgromian dynamics not only naturally leads to rotation velocities being typically higher in the extended polar rings than in the hosts, as would be the case in Newtonian dynamics without dark matter, but that it also gets the shape and amplitude of velocities correct. Milgromian dynamics thus adequately explains this particular property of polar ring galaxies.

We present measurements of the optical broadband atmospheric extinction coefficients and the night sky brightness at the Xuyi Observational Station of Purple Mountain Observatory (PMO). The measurements are based on CCD imaging data taken in the Sloan Digital Sky Survey g, r and i bands with the Xuyi 1.04/1.20m Schmidt Telescope for the Xuyi Schmidt Telescope Photometric Survey of the Galactic Anti-center (XSTPS-GAC), the photometric part of the Digital Sky Survey of the Galactic Anti-center (DSS-GAC). The data were collected in more than 130 winter nights from 2009 to 2011. We find that the atmospheric extinction coefficients for the g, r and i bands are 0.70, 0.55 and 0.38 mag/airmass, respectively, based on observations taken in several photometric nights. The night sky brightness determined from images of good quality has median val- ues of 21.7, 20.8 and 20.0 mag/arcsec2 and reaches 22.1, 21.2 and 20.4 mag/arcsec2 under the best observing conditions for the g, r and i bands, respectively. The relatively large extinction coefficients compared with other good astronomical observing sites are mainly due to the relatively low elevation (i.e. 180 m) and high humidity of the Station.

The research about asteroids attracts more and more attention recently, especially focusing on their physical structures, such as the spin axis, the rotation period and the shape. The long distance between Earth observers and asteroids makes it impossible to get the shape and other parameters of asteroids directly with the exception of the NEAs (Near Earth Asteroids) and others passed by some spacecrafts. Generally photometric measurement is still the main way to obtain the research data for asteroids now, i.e. the lightcurves recording the brightness and positions of asteroids. Supposing that the shape of the asteroid is a triaxial ellipsoid with a stable spinning status, a new method is present in this article to reconstruct the shape models of asteroids from the lightcurves, with the other physical parameters together. By applying a special curvature function, the method calculates the brightness integration on a unit sphere and Lebedev Quadrature is employed for the discretization. At last the method searches the optimal solution by Levenberg-Marquardt algorithm to minimize the residual of the brightness. By adopting this method not only related physical parameters of asteroids can be obtained at a reasonable accuracy, but also a simple shape model of Ellipsoid can be generated for reconstructing more sophisticated shape model further.

We present a unified parameterization of the fitting functions suitable for density profiles of dark matter haloes or elliptical galaxies. A notable feature is that the classical Einasto profile appears naturally as the continuous limiting case of the cored subfamily amongst the double power-law profiles of Zhao (1996). Based on this, we also argue that there is basically no qualitative difference between halo models well-fitted by the Einasto profile and the standard NFW model. This may even be the case quantitatively unless the resolutions of simulations and the precisions of fittings are sufficiently high to make meaningful distinction possible.

We use Schwarzschild's orbit-superposition technique to construct self-consistent models of the Galactic bar. Using $\chi^2$ minimisation, we find that the best-fit Galactic bar model has a pattern speed $\Omega_{\rm p}=60 \rm{km s^{-1} kpc^{-1}}$, disk mass $\rm{M_{\rm d}=1.0\times10^{11}M_{\odot}}$ and bar angle $\theta_{\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.

Using a deprojection technique, we study the X-ray properties of the galaxy cluster Abell 1835 observed with Chandra, including temperature, abundance, electron density, gas mass fraction, and total mass. A comparison with the results without deprojection shows that the properties do not change much. When we compare the results with those of XMM-Newton, the difference between the temperature profiles derived from Chandra and XMM-Newton data still exists, even if the point-spread function effect of XMM-Newton is corrected. To investigate the reasons for the difference, we used the spectra to cross-calibrate the temperatures. They show that the Chandra spectra can be fitted well with XMM-Newton temperatures. Furthermore, we derive the electron density profile from Chandra data with XMM-Newton temperatures and calculate the projected mass, which is consistent with the XMM-Newton mass and a little lower than the weak lensing mass at r_500. Thus, it seems that the temperature derived from XMM-Newton may be more reliable.

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.

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.

(1) We find that comet 107P/WH was active in 1949, 1979, 1992, 2005, and 2009. (2) Its age can be measured. We find T-AGE=4700 comet years, WB-AGE=7800 cy. (3) This is a methuselah comet very near to its dormancy phase, being temporarily rejuvenated due to a diminution of its perihelion distance. (4) We measured the diameter as Deffe=3.67\pm0.06 km, and the rotational period, Prot=6.093\pm0.002 h. (5) We define the region of the graveyard of comets, and three comets belong to the graveyard: 107P/Wilson-Harrington, 133P/Elst-Pizarro and D/1891W1 Blanpain.

The Damocloids are a group of unusual asteroids, recently enrolling a new member of 2010 EJ104. The dynamical evolution for the Damocloids may uncover a connection passage from the Main Belt, the Kuiper Belt and the scattered disk beyond. According to our simulations, two regions may be considered as possible origin of the Damocloids: the scattered disk, or a part of Oort cloud which will be perturbed to a transient region locating between 700 AU and 1000 AU. Based on the potential origin, the Damocloids can be classified into two types, with relation to their semi-major axes, and about 65.5% Damocloids is classified into type I which mainly originate from Oort cloud. Whether the Damocloids is inactive nuclei of Halley Family Comets may rely on their origin.

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.

We study mass models that correspond to MOND (triaxial) potentials for which the Hamilton-Jacobi equation separates in ellipsoidal coordinates. The problem is first discussed in the simpler case of deep-MOND systems, and then generalized to the full MOND regime. We prove that the Kuzmin property for Newtonian gravity still holds, i.e., that the density distribution of separable potentials is fully determined once the density profile along the minor axis is assigned. At variance with the Newtonian case, the fact that a positive density along the minor axis leads to a positive density everywhere remains unproven. We also prove that (i) all regular separable models in MOND have a vanishing density at the origin, so that they would correspond to centrally dark-matter dominated systems in Newtonian gravity; (ii) triaxial separable potentials regular at large radii and associated with finite total mass leads to density distributions that at large radii are not spherical and decline as ln(r)/r^5; (iii) when the triaxial potentials admit a genuine Frobenius expansion with exponent 0<epsilon<1, the density distributions become spherical at large radii, with the profile ln(r)/r^(3+2epsilon). After presenting a suite of positive density distributions associated with MOND separable potentials, we also consider the important family of (non-separable) triaxial potentials V_1 introduced by de Zeeuw and Pfenniger, and we show that, as already known for Newtonian gravity, they obey the Kuzmin property also in MOND. The ordinary differential equation relating their potential and density along the z-axis is an Abel equation of the second kind that, in the oblate case, can be explicitly reduced to canonical form.

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/.

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.

With the unique database from Michelson Doppler Imager aboard the Solar and Heliospheric Observatory in an interval embodying solar cycle 23, the cyclic behavior of solar small-scale magnetic elements is studied. More than 13 million small-scale magnetic elements are selected, and the following results are unclosed. (1) The quiet regions dominated the Sun's magnetic flux for about 8 years in the 12.25 year duration of Cycle 23. They contributed (0.94 - 1.44) $\times 10^{23}$ Mx flux to the Sun from the solar minimum to maximum. The monthly average magnetic flux of the quiet regions is 1.12 times that of active regions in the cycle. (2) The ratio of quiet region flux to that of the total Sun equally characterizes the course of a solar cycle. The 6-month running-average flux ratio of quiet region had been larger than 90.0% for 28 continuous months from July 2007 to October 2009, which characterizes very well the grand solar minima of Cycles 23-24. (3) From the small to large end of the flux spectrum, the variations of numbers and total flux of the network elements show no-correlation, anti-correlation, and correlation with sunspots, respectively. The anti-correlated elements, covering the flux of (2.9 - 32.0)$\times 10^{18}$ Mx, occupies 77.2% of total element number and 37.4% of quiet Sun flux. These results provide insight into reason for anti-correlated variations of small-scale magnetic activity during the solar cycle.

With the unique database from Michelson Doppler Imager aboard the Solar and Heliospheric Observatory in an interval embodying solar cycle 23, the cyclic behavior of solar small-scale magnetic elements is studied. More than 13 million small-scale magnetic elements are selected, and the following results are unclosed. (1) The quiet regions dominated the Sun\textsf{'}s magnetic flux for about 8 years in the 12.25 year duration of Cycle 23. They contributed (0.94 -- 1.44) $\times 10^{23}$ Mx flux to the Sun from the solar minimum to maximum. The monthly average magnetic flux of the quiet regions is 1.12 times that of active regions in the cycle. (2) The ratio of quiet region flux to that of the total Sun equally characterizes the course of a solar cycle. The 6-month running-average flux ratio of quiet region had been larger than 90.0% for 28 continuous months from July 2007 to October 2009, which characterizes very well the grand solar minima of Cycles 23-24. (3) From the small to large end of the flux spectrum, the variations of numbers and total flux of the network elements show no-correlation, anti-correlation, and correlation with sunspots, respectively. The anti-correlated elements, covering the flux of (2.9 - 32.0)$\times 10^{18}$ Mx, occupies 77.2% of total element number and 37.4% of quiet Sun flux. These results provide insight into reason for anti-correlated variations of small-scale magnetic activity during the solar cycle.

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.

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.

We study the interaction of atomic and molecular hydrogen with a surface of tholin, a man-made polymer considered to be an analogue of aerosol particles present in Titan's atmosphere, using thermal programmed desorption at low temperatures below 30 K. The results are fitted and analyzed using a fine-grained rate equation model that describes the diffusion, reaction and desorption processes. We obtain the energy barriers for diffusion and desorption of atomic and molecular hydrogen. These barriers are found to be in the range of 30 to 60 meV, indicating that atom/molecule-surface interactions in this temperature range are dominated by weak adsorption forces. The implications of these results for the understanding of the atmospheric chemistry of Titan are discussed.

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 $\Lambda$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.

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.

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.

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.

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.

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.

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.

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.

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.

We present new results from BRAVA, a large scale radial velocity survey of the Galactic bulge, using M giant stars selected from the Two Micron All Sky Survey catalogue as targets for the Cerro Tololo Inter-American Observatory 4m Hydra multi-object spectrograph. The purpose of this survey is to construct a new generation of self-consistent bar models that conform to these observations. We report the dynamics for fields at the edge of the Galactic bulge at latitudes b=-8 deg. and compare to the dynamics at b=-4 deg. We find that the rotation curve V(r) is the same at b=-8 deg. as at b=-4 deg. That is, the Galactic boxy bulge rotates cylindrically, as do boxy bulges of other galaxies. The summed line of sight velocity distribution at b=-8 deg. is Gaussian, and the binned longitude-velocity plot shows no evidence for either a (disk) population with cold dynamics or for a (classical bulge) population with hot dynamics. The observed kinematics are well modeled by an edge-on N-body bar, in agreement with published structural evidence. Our kinematic observations indicate that the Galactic bulge is a prototypical product of secular evolution in galaxy disks, in contrast with stellar population results that are most easily understood if major mergers were the dominant formation process.

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.

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.

We perform a stability test of triaxial models in MOdified Newtonian Dynamics (MOND) using N-body simulations. The triaxial models considered here have densities that vary with $r^{-1}$ in the center and $r^{-4}$ at large radii. The total mass of the model varies from $10^8\Msun$ to $10^{10}\Msun$, representing the mass scale of dwarfs to medium-mass elliptical galaxies, respectively, from deep MOND to quasi-Newtonian gravity. We build triaxial galaxy models using the Schwarzschild technique, and evolve the systems for 200 Keplerian dynamical times (at the typical length scale of 1.0 kpc). We find that the systems are virial overheating, and in quasi-equilibrium with the relaxation taking approximately 5 Keplerian dynamical times (1.0 kpc). For all systems, the change of the inertial (kinetic) energy is less than 10% (20%) after relaxation. However, the central profile of the model is flattened during the relaxation and the (overall) axis ratios change by roughly 10% within 200 Keplerian dynamical times (at 1.0kpc) in our simulations. We further find that the systems are stable once they reach the equilibrium state.

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.

We look for observational signatures that could discriminate between Newtonian and modified Newtonian (MOND) dynamics in the Milky Way, in view of the advent of large astrometric and spectroscopic surveys. Indeed, a typical signature of MOND is an apparent disk of "phantom" dark matter, which is uniquely correlated with the visible disk-density distribution. Due to this phantom dark disk, Newtonian models with a spherical halo have different signatures from MOND models close to the Galactic plane. The models can thus be differentiated by measuring dynamically (within Newtonian dynamics) the disk surface density at the solar radius, the radial mass gradient within the disk, or the velocity ellipsoid tilt angle above the Galactic plane. Using the most realistic possible baryonic mass model for the Milky Way, we predict that, if MOND applies, the local surface density measured by a Newtonist will be approximately 78 Msun/pc2 within 1.1 kpc of the Galactic plane, the dynamically measured disk scale-length will be enhanced by a factor of 1.25 with respect to the visible disk scale-length, and the local vertical tilt of the velocity ellipsoid at 1 kpc above the plane will be approximately 6 degrees. None of these tests can be conclusive for the present-day accuracy of Milky Way data, but they will be of prime interest with the advent of large surveys such as GAIA.

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.

Whereas considerable effort has been afforded in understanding the properties of galaxies, a full physical picture, connecting their baryonic and dark-matter content, super-massive black holes, and (metric) theories of gravity, is still ill-defined. Strong gravitational lensing furnishes a powerful method to probe gravity in the central regions of galaxies. It can (1) provide a unique detection-channel of dark-matter substructure beyond the local galaxy group, (2) constrain dark-matter physics, complementary to direct-detection experiments, as well as metric theories of gravity, (3) probe central super-massive black holes, and (4) provide crucial insight into galaxy formation processes from the dark matter point of view, independently of the nature and state of dark matter. To seriously address the above questions, a considerable increase in the number of strong gravitational-lens systems is required. In the timeframe 2010-2020, a staged approach with radio (e.g. EVLA, e-MERLIN, LOFAR, SKA phase-I) and optical (e.g. LSST and JDEM) instruments can provide 10^(2-4) new lenses, and up to 10^(4-6) new lens systems from SKA/LSST/JDEM all-sky surveys around ~2020. Follow-up imaging of (radio) lenses is necessary with moderate ground/space-based optical-IR telescopes and with 30-50m telescopes for spectroscopy (e.g. TMT, GMT, ELT). To answer these fundamental questions through strong gravitational lensing, a strong investment in large radio and optical-IR facilities is therefore critical in the coming decade. In particular, only large-scale radio lens surveys (e.g. with SKA) provide the large numbers of high-resolution and high-fidelity images of lenses needed for SMBH and flux-ratio anomaly studies.

The statistics of high speed satellite galaxies, as reported in the recent literature, can be a powerful diagnosis of the depth of the potential well of the host halo, and hence discriminate between competing gravitational theories. Naively one expects that high speed satellites are more common in Modified Newtonian Dynamics (MOND) than in cold dark matter (CDM) since an isolated MONDian system has an infinite potential well, while CDM halos have finite potential wells. In this \textit{Letter} we report on an initial test of this hypothesis in the context of the first generation of cosmological simulations utilising a rigorous MONDian Poisson solver. We find that such high speed encounters are approximately a factor of four more common in MOND than in the concordance $\Lambda$CDM model of cosmic structure formation.

Modern astronomical data on galaxy and cosmological scales have revealed powerfully the existence of certain dark sectors of fundamental physics, i.e., existence of particles and fields outside the standard models and inaccessible by current experiments. Various approaches are taken to modify/extend the standard models. Generic theories introduce multiple de-coupled fields A, B, C, each responsible for the effects of DM (cold supersymmetric particles), DE (Dark Energy) effect, and MG (Modified Gravity) effect respectively. Some theories use adopt vanilla combinations like AB, BC, or CA, and assume A, B, C belong to decoupled sectors of physics. MOND-like MG and Cold DM are often taken as opposite frameworks, e.g. in the debate around the Bullet Cluster. Here we argue that these ad hoc divisions of sectors miss important clues from the data. The data actually suggest that the physics of all dark sectors is likely linked together by a self-interacting oscillating field, which governs a chameleon-like dark fluid, appearing as DM, DE and MG in different settings. It is timely to consider an interdisciplinary approach across all semantic boundaries of dark sectors, treating the dark stress as one identity, hence accounts for several "coincidences" naturally.

We study the amplitude of the weak gravitational lensing signal as a function of stellar mass around a sample of relatively isolated galaxies. This selection of lenses simplifies the interpretation of the observations, which consist of data from the Red- sequence Cluster Survey and the Sloan Digital Sky Survey. We find that the amplitude of the lensing signal as a function of stellar mass is well described by a power law with a best fit slope \alpha= 0.74 \pm 0.08. This result is inconsistent with Modified Newtonian Dynamics, which predicts \alpha = 0.5 (we find \alpha > 0.5 with 99.7% confidence). As a related test, we determine the MOND mass-to-light ratio as a function of luminosity. Our results require dark matter for the most luminous galaxies (L >=10^11 L_sun). We rule out an extended halo of gas or active neutrinos as a way of reconciling our findings with MOND. Although we focus on a single alternative gravity model, we note that our results provide an important test for any alternative theory of gravity.

Strong gravitational lensing by galaxies in MOdified Newtonian Dynamics (MOND) has until now been restricted to spherically symmetric models. These models were able to account for the size of the Einstein ring of observed lenses, but were unable to account for double-imaged systems with collinear images, as well as four-image lenses. Non-spherical models are generally cumbersome to compute numerically in MOND, but we present here a class of analytic non-spherical models that can be applied to fit double-imaged and quadruple-imaged systems. We use them to obtain a reasonable MOND fit to ten double-imaged systems, as well as to the quadruple-imaged system Q2237+030 which is an isolated bulge-disc lens producing an Einstein cross. However, we also find five double-imaged systems and three quadruple-imaged systems for which no reasonable MOND fit can be obtained with our models. We argue that this is mostly due to the intrinsic limitation of the analytic models, even though the presence of small amounts of additional dark mass on galaxy scales in MOND is also plausible.

The diffuse Galactic $\gamma$-rays from EGRET observation shows excesses above 1 GeV in comparison with the expectations from conventional Galactic cosmic ray (CR) propagation model. In the work we try to solve the "GeV excess" problem by dark matter (DM) annihilation in the frame of supersymmetry (SUSY). Compared with previous works, there are three aspects improved in this work: first, the direction-independent "boost factor" for diffuse $\gamma$-rays from dark matter annihilation (DMA) is naturally reproduced by taking the DM substructures into account; second, there is no need for renormalization of the diffuse $\gamma$-ray background produced by CRs; last but not the least, in this work our new propagation model can give consistent results of both diffuse $\gamma$-rays and antiprotons, by directly adding the signals from DMA to the diffuse $\gamma$-ray background. This is a self-consistent model among several possible scenarios at present, and can be tested or optimized by the forthcoming experiments such as GLAST, PAMELA and AMS02.

Intercluster filaments negligibly contribute to the weak lensing signal in general relativity (GR), $\gamma_{N}\sim 10^{-4}-10^{-3}$. In the context of relativistic modified Newtonian dynamics (MOND) introduced by Bekenstein, however, a single filament inclined by $\approx 45^\circ$ from the line of sight can cause substantial distortion of background sources pointing towards the filament's axis ($\kappa=\gamma=(1-A^{-1})/2\sim 0.01$); this is rigorous for infinitely long uniform filaments, but also qualitatively true for short filaments ($\sim 30$Mpc), and even in regions where the projected matter density of the filament is equal to zero. Since galaxies and galaxy clusters are generally embedded in filaments or are projected on such structures, this contribution complicates the interpretation of the weak lensing shear map in the context of MOND. While our analysis is of mainly theoretical interest providing order-of-magnitude estimates only, it seems safe to conclude that when modeling systems with anomalous weak lensing signals, e.g. the "bullet cluster" of Clowe et al., the "cosmic train wreck" of Abell 520 from Mahdavi et al., and the "dark clusters" of Erben et al., filamentary structures might contribute in a significant and likely complex fashion. On the other hand, our predictions of a (conceptual) difference in the weak lensing signal could, in principle, be used to falsify MOND/TeVeS and its variations.

We present a novel solver for an analogue to Poisson's equation in the framework of modified Newtonian dynamics (MOND). This equation is highly non-linear and hence standard codes based upon tree structures and/or FFT's in general are not applicable; one needs to defer to multi-grid relaxation techniques. After a detailed description of the necessary modifications to the cosmological N-body code AMIGA (formerly known as MLAPM) we utilize the new code to revisit the issue of cosmic structure formation under MOND. We find that the proper (numerical) integration of a MONDian Poisson's equation has some noticable effects on the final results when compared against simulations of the same kind but based upon rather ad-hoc assumptions about the properties of the MONDian force field. Namely, we find that the large-scale structure evolution is faster in our revised MOND model leading to an even stronger clustering of galaxies, especially when compared to the standard LCDM paradigm.

Results from the Bulge Radial Velocity Assay (BRAVA) are presented. BRAVA uses M giant stars, selected from the 2MASS catalog to lie within a bound of reddening corrected color and luminosity, as targets for the Cerro Tololo Inter-American Observatory 4-m Hydra multi-object spectrograph. Three years of observations investigate the kinematics of the Galactic bulge major and minor axes with ~3300 radial velocities from 32 bulge fields and one disk field. We construct a longitude-velocity plot for the bulge stars and find that, contrary to previous studies, the bulge does not rotate as a solid body; for |l|<4 degree the rotation curve has a slope of roughly 100 km/s/kpc and flattens considerably at greater l, reaching a maximum rotation of 75 km/s. We compare our rotation curve and velocity dispersion profile to both the self-consistent model of Zhao (1996) and to N-body models; neither satisfactorily fit both measured rotation and dispersion. We place the bulge on the plot of (Vmax/sigma) vs. epsilon and find that the bulge lies near the oblate rotator line, and very close to the parameters of NGC 4565, an edge-on spiral galaxy with a bulge similar to that of the Milky Way. We find that our summed velocity distribution of bulge stars appears to be sampled from a Gaussian distribution, with sigma=116+/-2 km/s for our summed bulge fields.

The tight correlation between galaxy bulges and their central black hole masses likely emerges in a phase of rapid collapse and starburst at high redshift, due to the balance of gravity on gas with the feedback force from starbursts and the wind from the black hole; the average gravity on per unit mass of gas is ~ 2 x 10^-10 m/sec^2 during the star burst phase. This level of gravity could come from the real r^{-1} cusps of Cold Dark Matter (CDM) halos, but the predicted gravity would have a large scatter due to dependence on cosmological parameters and formation histories. Better agreement is found with the gravity from the scalar field in some co-variant versions of MOND, which can create the mirage of a Newtonian effective dark halo of density Pi r^{-1} near the center, where the characteristic surface density Pi=130alpha^{-1} Msun pc^{-2} and alpha is a fundamental constant of order unity fixed by the Lagrangian of the co-variant theory if neglecting environmental effects. We show with a toy analytical model and a hydrodynamical simulation that a constant background gravity due to MOND/TeVeS scalar field implies a critical pressure synchronizing starbursts and the formation of galaxy bulges and ellipticals. A universal threshold for the formation of the brightest regions of galaxies in a MONDian universe suggests that the central BHs, bulges and ellipticals would respect tight correlations like the M_{bulge}-M_{BH}-sigma relations. In general MOND tends to produce tight correlations in galaxy properties because its effective halo has less freedom and scatter than CDM halos.

Clusters of galaxies offer a robust test bed for probing the nature of dark matter that is insensitive to the assumption of the gravity theories. Both Modified Newtonian Dynamics (MOND) and General Relativity (GR) would require similar amounts of non-baryonic matter in clusters as MOND boosts the gravity only mildly on cluster scales. Gravitational lensing allows us to estimate the enclosed mass in clusters on small (20 - 50 kpc) and large (several 100 kpc) scales independent of the assumptions of equilibrium. Here we show for the first time that a combination of strong and weak gravitational lensing effects can set interesting limits on the phase space density of dark matter in the centers of clusters. The phase space densities derived from lensing observations are inconsistent with neutrino masses ranging from 2 - 7 eV, and hence do not support the 2 eV-range particles required by MOND. To survive, plausible modifications for MOND may be either an additional degree of dynamical freedom in a co-variant incarnation or mass-varying theories of neutrinos.

Several rare coincidences of scales in standard particle physics are needed to explain (i) why neutrinos have mass, (ii) why the negative pressure of the cosmological dark energy (DE) coincides with the positive pressure of random motion of dark matter (DM) in bright galaxies, (iii) why Dark Matter in galaxies seems to have a finite phase-space density, and to follow the Tully-Fisher-Milgrom relation of galaxy rotation curves. The old idea of self-interacting DM is given a new spin: we propose that the neutrino spin-gravity coupling could lead to a cosmic neutrino dark fluid with a an internal energy density varying as function of the local acceleration of the neutrino fluid with respect the CMB background. We link the Tully-Fisher-Milgrom relation of spiral galaxies (or MOND) with the relativistic pressure of the neutrino dark fluid without modifying Einsteinian gravity.

We propose a method to study and constrain modified gravity theories for dark matter using CMB temperature anisotropies and polarization. We assume that the theories considered here have already passed the matter power-spectrum test of large-scale structure. With this requirement met, we show that a modified gravity theory can be specified by parametrizing the time evolution of its dark-matter density contrast, which is completely controlled by the dark matter stress history. We calculate how the stress history with a given parametrization affects the CMB observables, and a qualitative discussion of the physical effects involved is supplemented with numerical examples. It is found that, in general, alternative gravity theories can be efficiently constrained by the CMB temperature and polarization spectra. There exist, however, special cases where modified gravity cannot be distinguished from the CDM model even by using both CMB and matter power spectrum observations, nor can they be efficiently restricted by other observables in perturbed cosmologies. Our results show how the stress properties of dark matter, which determine the evolutions of both density perturbations and the gravitational potential, can be effectively investigated using just the general conservation equations and without assuming any specific theoretical gravitational theory within a wide class.

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 $\nu\Lambda$ 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.

We compute the Milky Way potential in different cold dark matter (CDM) based models, and compare these with the modified Newtonian dynamics (MOND) framework. We calculate the axis ratio of the potential in various models, and find that isopotentials are less spherical in MOND than in CDM potentials. As an application of these models, we predict the escape velocity as a function of the position in the Galaxy. This could be useful in comparing with future data from planned or already-underway kinematic surveys (RAVE, SDSS, SEGUE, SIM, GAIA or the hypervelocity stars survey). In addition, the predicted escape velocity is compared with the recently measured high proper motion velocity of the Large Magellanic Cloud (LMC). To bind the LMC to the Galaxy in a MOND model, while still being compatible with the RAVE-measured local escape speed at the Sun's position, we show that an external field modulus of less than $0.03 a_0$ is needed.

Various TeVeS-inspired and f(R)-inspired theories of gravity have added an interesting twist to the search for dark matter and vacuum energy, modifying the landscape of astrophysics day by day. These theories can be together called a {\bf N}on-{\bf u}niform Dark Energy fluid (a Nu-Lambda fluid or a ${\mathbf V\Lambda}$ fluid); a common thread of these theories, according of an up-to-date summary by HZL \cite{Halle}, is a non-uniform vector field, describing an uneven vacuum energy fluid. The so-called "alternative" gravity theories are in fact in the standard GR gravity framework except that the cosmological "constant" is replaced by a non-trivial non-uniform vacuum energy, which couples the effects of Dark Matter and Dark Energy together by a single field. Built initially bottom-up rather than top-down as most gravity theories, TeVeS-inspired theories are healthily rooted on empirical facts. Here I attempt a review of some sanity checks of these fast-developing theories from galaxy rotation curves, gravitational lensing and cosmic acceleration. I will also discuss some theoretical aspects of the vacuum energy, and point out some analogies with electromagnetism and the Casimir effect.

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-$\Lambda$ model) closely resemble the $\Lambda$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.

The Bekenstein-Milgrom gravity theory with a modified Poisson equation is tested here for the existence of triaxial equilibrium solutions. Using the non-negative least square method, we show that self-consistent triaxial galaxies exist for baryonic models with a mild density cusp $\rho \sim {\Sigma \over r}$. Self-consistency is achieved for a wide range of central concentrations, $\Sigma \sim 10-1000\mathrm{M_{\odot}pc^{-2}}$, representing low-to-high surface brightness galaxies. Our results demonstrate for the first time that the orbit superposition technique is fruitful for constructing galaxy models beyond Newtonian gravity, and triaxial cuspy galaxies might exist without the help of Cold dark Matter.

The knowledge of magnetic topology is the key to understand magnetic energy release in astrophysics. Based on observed vector magnetograms, we have determined threedimensional (3D) topology skeleton of the magnetic fields in active region NOAA 10720. The skeleton consists of six 3D magnetic nulls and a network of corresponding spines, fans, and null-null lines. For the first time, we have identified a spiral magnetic null in Sun's corona. The magnetic lines of force twisted around the spine of the null, forming a 'magnetic wreath' with excess of free magnetic energy and resembling observed brightening structures at extraultraviolet (EUV) wavebands. We found clear evidence of topology eruptions which are referred to as the catastrophic changes of topology skeleton associated with a coronal mass ejection (CME) and an explosive X-ray flare. These results shed new lights in exploring the structural complexity and its role in explosive magnetic activity. In solar astrophysics and space science, the concept of flux rope has been widely used in modelling explosive magnetic activity, although their observational identity is obscure or, at least, lacking of necessary details. The current work suggests that the magnetic wreath associated with the 3D spiral null is likely an important class of the physical entity of flux ropes.

We present the analysis of 23 published rotation curves of disk galaxies belonging to the Ursa Major group of galaxies, with kinematics free of irregularities. The rotation curves are analysed in the context of MOND (Modified Newtonian Dynamics). We add an extra component to the rotation curve fits, in addition to the stellar and gaseous disks: a speculative halo of constant density made of, e.g., neutrinos, which would solve the bulk of the problem currently faced by MOND on rich galaxy clusters scales. We find that this additional unseen mass density is poorly constrained (as expected a priori, given that a neutrino halo never dominates the kinematics), but we also find that the best-fit value is non-zero: rho = 3.8 x 10^{-27} g/cm^3, and that a zero-density is marginally excluded with 87% confidence; also, the 95% confidence upper limit for the density is rho = 9.6 x 10^{-27} g/cm^3. These limits are slightly above the expectations from the Tremaine-Gunn phase space constraints on ordinary 2 eV neutrinos, but in accordance with the maximum density expected for one or two species of 5 eV sterile neutrinos.

We summarize the status of Modified Newtonian Dynamics (MOND) in galaxy clusters. The observed acceleration is typically larger than the acceleration threshold of MOND in the central regions, implying that some dark matter is necessary to explain the mass discrepancy there. A plausible resolution of this issue is that the unseen mass in MOND is in the form of ordinary neutrinos with masses just below the experimentally detectable limit. In particular, we show that the lensing mass reconstructions of the clusters 1E0657-56 (the bullet cluster) and Cl0024+17 (the ring) do not pose a new challenge to this scenario. However, the mass discrepancy for cool X-ray emitting groups, in which neutrinos cannot cluster, pose a more serious problem, meaning that dark baryons could present a more satisfactory solution to the problem of unseen mass in MOND clusters.

We present recent results from a Keck study of the composition of the Galactic bulge, as well as results from the bulge Bulge Radial Velocity Assay (BRAVA). Culminating a 10 year investigation, Fulbright, McWilliam, & Rich (2006, 2007) solved the problem of deriving the iron abundance in the Galactic bulge, and find enhanced alpha element abundances, consistent with the earlier work of McWilliam & Rich (1994). We also report on a radial velocity survey of {\sl 2MASS}-selected M giant stars in the Galactic bulge, observed with the CTIO 4m Hydra multi-object spectrograph. This program is to test dynamical models of the bulge and to search for and map any dynamically cold substructure in the Galactic bulge. We show initial results on fields at $-10^{\circ} < l <+10^{\circ}$ and $b=-4^{\circ}$. We construct a longitude-velocity plot for the bulge stars and the model data, and find that contrary to previous studies, the bulge does not rotate as a solid body; from $-5^{\circ}<l<+5^{\circ}$ the rotation curve has a slope of $\approx 100 km s^{-1}$ and flattens considerably at greater $l$ and reaches a maximum rotation of $45 {km s^{-1}}$ (heliocentric) or $\sim 70 {km s^{-1}}$ (Galactocentric). This rotation is slower than that predicted by the dynamical model of Zhao (1996).

We model the electrons/positrons produced by dark matter annihilations in the colliding galaxy cluster system 1E0657-56 (the bullet cluster). These charged particles, confined by the Magnetic field, clearly trace the path of the bullet, which passes through the main cluster with a speed of 3000-5000 km/s. Adding the effect of dark matter substructure (subhalos) in each cluster we find the annihilation rate is enhancedgreatly and the density of positrons in the trail is similar to that within the bullet cluster. This opens the door to unambiguous detection of the annihilation signal through the SZ effect, at significant separation from the thermal electrons.

We investigate how different models that have been proposed for solving the dark matter problem can fit the velocity dispersion observed around elliptical galaxies, on either a small scale (~ 20kpc) with stellar tracers, such as planetary nebulae, or large scale (~ 200kpc) with satellite galaxies as tracers. Predictions of Newtonian gravity, either containing pure baryonic matter, or embedded in massive cold dark matter (CDM) haloes, are compared with predictions of the modified gravity of MOND. The standard CDM model has problems on a small scale, and the Newtonian pure baryonic model has difficulties on a large scale, while a fit with MOND is possible on both scales.

A rare coincidence of scales in standard particle physics is needed to explain why $\Lambda$ or the negative pressure of cosmological dark energy (DE) coincides with the positive pressure $P_0$ of random motion of dark matter (DM) in bright galaxies. Recently Zlosnik et al. (2007) propose to modify the Einsteinian curvature by adding a non-linear pressure from a medium flowing with a four-velocity vector field $U^\mu$. We propose to check whether a smooth extension of GR with a simple kinetic Lagrangian of $U^\mu$ can be constructed, and whether the pressure can bend space-time sufficiently to replace the roles of a $w=-1$ DE, $w=0$ Cold DM and heavy neutrinos in explaining anomalous accelerations at all scales. As a specific proof of concept we find a Vector-for-$\Lambda$ model (${\mathbf V\Lambda}$-model) and its variants. With essentially {\it no free parameters}, these appear broadly consistent with the solar system, gravitational potentials in dwarf spiral galaxies and the bullet cluster of galaxies, early universe with inflation, structure formation and BBN, and late acceleration with a 1:3 ratio of DM:DE.

Gravity is nearly a universal constant in the cusp of an NFW galaxy halo. Inside this external field an isothermal gas sphere will collapse and trigger a starburst if above a critical central pressure. Thus formed spheroidal stellar systems have Sersic-profile and satisfy the Faber-Jackson relation. The process is consistent with observed starbursts. We also recover the M_BH vs. velocity dispersion relation, if the gas collapse is regulated or resisted by the feedback from radiation from the central BH.

The recent SDSS measured velocity distribution of satellite galaxies has been modelled in the context of MOND. We show that even when the extra constraint of adhering to the projected satellite number density profile is added, the two line of sight (los) velocity dispersion profiles presented in Klypin & Prada (2007) can be matched simply with a radially varying anisotropy. Interestingly, the anisotropies required to fit the los velocity dispersions are remarkably similar to the anisotropies generated by dissipationless collapse simulations in MOND. The mass-to-light ratios of the two host galaxies used are sensible and positivity of the distribution function is satisfied.

Within the cold dark matter (CDM) framework tidal dwarf galaxies (TDGs) cannot contain dark matter, so the recent results by Bournaud et al. (2007) that 3 rotating TDGs do show significant evidence for being dark matter dominated is inconsistent with the current concordance cosmological theory unless yet another dark matter component is postulated. We confirm that the TDG rotation curves are consistent with Newtonian dynamics only if either an additional dark matter component is postulated, or if all 3 TDGs happen to be viewed nearly edge-on, which is unlikely given the geometry of the tidal debris. We also find that the observed rotation curves are very naturally explained without any free parameters within the modified Newtonian dynamics (MOND) framework if inclinations are adopted as derived by Bournaud et al. We explore different inclination angles and two different assumptions about the external field effect. The results do not change significantly, and we conclude therefore that Newtonian dynamics has severe problems while MOND does exceedingly well in explaining the observed rotation curves of the 3 TDGs studied by Bournaud et al.

The self-binding energy and stability of a galaxy in MOND-based gravity are curiously decreasing functions of its center of mass acceleration towards neighbouring mass concentrations. A tentative indication of this breaking of the Strong Equivalence Principle in field galaxies is the RAVE-observed escape speed in the Milky Way. Another consequence is that satellites of field galaxies will move on nearly Keplerian orbits at large radii (100 - 500 kpc), with a declining speed below the asymptotically constant naive MOND prediction. But consequences of an environment-sensitive gravity are even more severe in clusters, where member galaxies accelerate fast: no more Dark-Halo-like potential is present to support galaxies, meaning that extended axisymmetric disks of gas and stars are likely unstable. These predicted reappearance of asymptotic Keplerian velocity curves and disappearance of "stereotypic galaxies" in clusters are falsifiable with targeted surveys.

We investigate the possibility of discriminating between Modified Newtonian Dynamics (MOND) and Newtonian gravity with dark matter, by studying the vertical dynamics of disk galaxies. We consider models with the same circular velocity in the equatorial plane (purely baryonic disks in MOND and the same disks in Newtonian gravity embedded in spherical dark matter haloes), and we construct their intrinsic and projected kinematical fields by solving the Jeans equations under the assumption of a two-integral distribution function. We found that the vertical velocity dispersion of deep-MOND disks can be much larger than in the equivalent spherical Newtonian models. However, in the more realistic case of high-surface density disks this effect is significantly reduced, casting doubts on the possibility of discriminating between MOND and Newtonian gravity with dark matter by using current observations.

We present semi-analytical constraint on the amount of dark matter in the merging bullet galaxy cluster using the classical Local Group timing arguments. We consider particle orbits in potential models which fit the lensing data. {\it Marginally consistent} CDM models in Newtonian gravity are found with a total mass M_{CDM} = 1 x 10^{15}Msun of Cold DM: the bullet subhalo can move with V_{DM}=3000km/s, and the "bullet" X-ray gas can move with V_{gas}=4200km/s. These are nearly the {\it maximum speeds} that are accelerable by the gravity of two truncated CDM halos in a Hubble time even without the ram pressure. Consistency breaks down if one adopts higher end of the error bars for the bullet gas speed (5000-5400km/s), and the bullet gas would not be bound by the sub-cluster halo for the Hubble time. Models with V_{DM}~ 4500km/s ~ V_{gas} would invoke unrealistic large amount M_{CDM}=7x 10^{15}Msun of CDM for a cluster containing only ~ 10^{14}Msun of gas. Our results are generalisable beyond General Relativity, e.g., a speed of $4500\kms$ is easily obtained in the relativistic MONDian lensing model of Angus et al. (2007). However, MONDian model with little hot dark matter $M_{HDM} \le 0.6\times 10^{15}\msun$ and CDM model with a small halo mass $\le 1\times 10^{15}\msun$ are barely consistent with lensing and velocity data.

We present a new test of modified Newtonian dynamics (MOND) on galactic scales, based on the escape speed in the solar neighbourhood. This test is independent from other empirical successes of MOND at reproducing the phenomenology of galactic rotation curves. The galactic escape speed in MOND is entirely determined by the baryonic content of the Galaxy and the external field in which it is embedded. We estimate that the external field in which the Milky Way must be embedded to produce the observed local escape speed of 550 km/s is of the order of a_0/100, where a_0 is the dividing acceleration scale below which gravity is boosted in MOND. This is compatible with the external gravitational field actually acting on the Milky Way.

Observations of spiral galaxies strongly support a one-to-one analytical relation between the inferred gravity of dark matter at any radius and the enclosed baryonic mass. It is baffling that baryons manage to settle the dark matter gravitational potential in such a precise way, leaving no "messy" fingerprints of the merging events and "gastrophysical" feedbacks expected in the history of a galaxy in a concordance Universe. This correlation of gravity with baryonic mass can be interpreted from several non-standard angles, especially as a modification of gravity called TeVeS, in which no galactic dark matter is needed. In this theory, the baryon-gravity relation is captured by the dieletric-like function mu of Modified Newtonian Dynamics (MOND), controlling the transition from 1/r^2 attraction in the strong gravity regime to 1/r attraction in the weak regime. Here, we study this mu-function in detail. We investigate the observational constraints upon it from fitting galaxy rotation curves, unveiling the degeneracy between the stellar mass-to-light ratio and the mu-function as well as the importance of the sharpness of transition from the strong to weak gravity regimes. We also numerically address the effects of non-spherical baryon geometry in the framework of non-linear TeVeS, and exhaustively examine how the mu-function connects with the free function of that theory. In that regard, we exhibit the subtle effects and wide implications of renormalizing the gravitational constant. We finally present a discontinuity-free transition between quasi-static galaxies and the evolving Universe for the free function of TeVeS, inevitably leading to a return to 1/r^2 attraction at very low accelerations in isolated galaxies.

We are undertaking a large scale radial velocity survey of the Galactic bulge which uses M giant stars selected from the 2MASS catalog as targets for the CTIO 4m Hydra multi-object spectrograph. The aim of this survey is to test dynamical models of the bulge and to quantify the importance, if any, of cold stellar streams in the bulge and its vicinity. Here we report on the kinematics of a strip of fields at -10 < l <+10 degres and b=-4 degres. We construct a longitude-velocity plot for the bulge stars and the model data, and find that contrary to previous studies, the bulge does not rotate as a solid body. From -5<l<+5 degrees the rotation curve has a slope of roughly 100 km/s/kpc and flattens considerably at greater l and reaches a maximum rotation of 45 km/s. We compare our rotation curve and velocity dispersion profile to both the self-consistent model of Zhao (1996) and to N-body models; neither fits both our observed rotation curve and velocity dispersion profile. The high precision of our radial velocities (3 km/s) yields an unexpected result: hints of cold kinematic features are seen in a number of the line of sight velocity distributions.

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.

Bekenstein's (2004) TeVeS theory has added an interesting twist to the search for dark matter and dark energy, modifying the landscape of gravity-related astronomy day by day. Built bottom-up rather than top-down as most gravity theories, TeVeS-like theories are healthily rooted on empirical facts, hence immediately passing sanity checks on galaxy rotation curves, solar system constraints, even bullet cluster of galaxies and cosmology with the help of 2eV neutrinos. Nonetheless, empirical checks are far from perfect and complete, and groups of different expertises are rapidly increasing the number of falsifiable properties of the theory. The theory has also been made much simpler and more general thanks to the work of Zlosnik, Ferreira, Starkman (gr-qc/0606039, astro-ph/0607411). Here I attempt a tutorial of how to compute lensing convergence, time delays etc in TeVeS-like theories for non-spherical lenses. I gave examples to illustrate a few common caveats of Dark-Matter-guided intuitions.

We develop a new method to predict the density associated with weak lensing maps of (un)relaxed clusters in a range of theories interpolating between GR and MOND (General Relativity and Modified Newtonian Dynamics). We apply it to fit the lensing map of the bullet merging cluster 1E0657-56, in order to constrain more robustly the nature and amount of collisionless matter in clusters {\it beyond} the usual assumption of spherical equilibrium (Pointecouteau & Silk 2005) and the validity of GR on cluster scales (Clowe et al. 2006). Strengthening the proposal of previous authors we show that the bullet cluster is dominated by a collisionless -- most probably non-baryonic -- component in GR as well as in MOND, a result consistent with the dynamics of many X-ray clusters. Our findings add to the number of known pathologies for a purely baryonic MOND, including its inability to fit the latest data from the Wilkinson Microwave Anisotropy Probe. A plausible resolution of all these issues and standard issues of Cold Dark Matter with galaxy rotation curves is the "marriage" of MOND with ordinary hot neutrinos of 2eV. This prediction is just within the GR-independent maximum of neutrino mass from current $\beta$-decay experiments, and is falsifiable by the Karlsruhe Tritium Neutrino (KATRIN) experiment by 2009. Issues of consistency with strong lensing arcs and the large relative velocity of the two clusters comprising the bullet cluster are also addressed.

We calculate the strong lensing probability as a function of the image-separation $\Delta\theta$ in TeVeS (tensor-vector-scalar) cosmology, which is a relativistic version of MOND (MOdified Newtonian Dynamics). The lens, often an elliptical galaxy, is modeled by the Hernquist profile. We assume a flat cosmology with $\Omega_b=1-\Omega_\Lambda=0.04$ and the simplest interpolating function $\mu(x)={\rm min}(1,x)$. For comparison, we recalculated the probabilities for lenses by Singular Isothermal Sphere (SIS) galaxy halos in LCDM with Schechter-fit velocity function. The amplification bias is calculated based on the magnification of the second bright image rather than the total of the two brighter images. Our calculations show that the Hernquist model predicts insufficient but acceptable probabilities in flat TeVeS cosmology compared with the results of the well defined combined sample of Cosmic Lens All-Sky Survey (CLASS) and Jodrell Bank/Very Large Array Astrometric Survey (JVAS); at the same time, it predicts higher probabilities than SIS model in LCDM at small image separations.

The phenomena customly described with the standard $\Lambda$CDM model are broadly reproduced by an extremely simple model in TeVeS, Bekenstein's (2004) modification of General Relativity motivated by galaxy phenomenology. Our model can account for the acceleration of the universe seen in SNeIa distances without a cosmological constant, and the accelerations seen rotation curves of nearby spiral galaxies and gravitational lensing of high-redshift elliptical galaxies without cold dark matter. The model is consistent with BBN and the neutrino mass. The TeVeS scalar field is shown to play the effective dural roles of Dark Matter and Dark Energy with the amplitudes of the effects controled by a $\mu$-function of the scalar field, called the $\mu-$essence here. We also discuss outliers to the theory's predictions on multi-imaged galaxy lenses and outliers on sub-galaxy scale.

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.

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.

A proper test of Modified Newtonian Dynamics (MOND) in systems of non-trivial geometries depends on modelling subtle differences in several versions of its postulated theories. This is especially important for lensing and dynamics of barely virialised galaxy clusters with typical gravity of scale $\sim a_0 \sim 1\AA{\rm s}^{-2}$. The original MOND formula, the classical single field modification of the Poisson equation, and the multi-field general relativistic theory of Bekenstein (TeVeS) all lead to different predictions as we stray from spherical symmetry. In this paper, we study a class of analytical MONDian models for a system with a semi-Hernquist baryonic profile. After presenting the analytical distribution function of the baryons in spherical limits, we develop orbits and gravitational lensing of the models in non-spherical geometries. In particular, we can generate a multi-centred baryonic system with a weak lensing signal resembling that of the merging galaxy cluster 1E 0657-56 with a bullet-like light distribution. We finally present analytical scale-free highly non-spherical models to show the subtle differences between the single field classical MOND theory and the multi-field TeVeS theory.

Dark Matter (DM) theories and mass-tracing-light theories like MOND are by construction nearly degenerate on galactic scales, but not when it comes to the predicted shapes of Roche Lobes of a two-body system (e.g., a globular cluster orbiting a host galaxy). We show that the flattening of the Roche lobe is sensitive to the function mu(g) in modification of the law of gravity. We generalise the analytical results obtained in the deep-MOND limit by Zhao (2005, astro-ph/0511713 and astro-ph/0512425), and consider a binary in the framework of a MOND-like gravity modification function mu(g) or a general non-Keplerian gravity g \propto R^-\zeta. We give analytical expressions for the inner Lagrange point and Robe lobe axis ratios. The Roche lobe volume is proven to scale linearly with the true mass ratio, which applies to any mu(g), hence mass-tracing light models would overpredict the Roche lobe of a DM-poor globular cluster in a DM-rich host galaxy, and underpredict the size of a DM-richer dwarf satellite. The lobes are squashed with the flattening ~ 0.4 in the strong gravity and ~ 0.6 in the weak gravity; a precise measurement of the flattening could be used to verify the anisotropic dilation effect which is generic to MOND-like gravity. We generalise these results for extended mass distribution, and compare predicted Roche radii in different gravity theories with limiting radii of observed globular clusters and dwarf galaxy satellites.

The phenomena customly called Dark Matter or Modified Newtonian Dynamics (MOND) have been argued by Bekenstein (2004) to be the consequences of a covariant scalar field, controlled by a free function (related to the MOND interpolating function) in its Lagrangian density. In the context of this relativistic MOND theory (TeVeS), we examine critically the interpolating function in the transition zone between weak and strong gravity. Bekenstein's toy model produces too gradually varying functions and fits rotation curves less well than the standard MOND interpolating function. However, the latter varies too sharply and implies an implausible external field effect (EFE). These constraints on opposite sides have not yet excluded TeVeS, but made the zone of acceptable interpolating functions narrower. An acceptable "toy" Lagrangian density function with simple analytical properties is singled out for future studies of TeVeS in galaxies. We also suggest how to extend the model to solar system dynamics and cosmology, and compare with strong lensing data (see also astro-ph/0509590).

We propose to use multiple-imaged gravitational lenses to set limits on gravity theories without dark matter, specificly TeVeS (Bekenstein 2004), a theory which is consistent with fundamental relativistic principles and the phenomenology of MOdified Newtonian Dynamics (MOND) theory. After setting the framework for lensing and cosmology, we derive analytically the deflection angle for the point lens and the Hernquist galaxy profile, and fit galaxy-quasar lenses in the CASTLES sample. We do this with three methods, fitting the observed Einstein ring sizes, the image positions, or the flux ratios. In all cases we consistently find that stars in galaxies in MOND/TeVeS provide adequate lensing. Bekenstein's toy $\mu$ function provides more efficient lensing than the standard MOND $\mu$ function. But for a handful of lenses [indicated in Table 2,3, fig 16] a good fit would require a lens mass orders of magnitude larger/smaller than the stellar mass derived from luminosity unless the modification function $\mu$ and modification scale $a_0$ for the universal gravity were allowed to be very different from what spiral galaxy rotation curves normally imply. We discuss the limitation of present data and summarize constraints on the MOND $\mu$ function. We also show that the simplest TeVeS "minimal-matter" cosmology, a baryonic universe with a cosmological constant, can fit the distance-redshift relation from the supernova data, but underpredicts the sound horizon size at the last scattering. We conclude that lensing is a promising approach to differentiate laws of gravity (see also astro-ph/0512425).

MOdified Newtonian Dynamics (MOND) is evolving from an empirical to a decent theory respecting fundamental physics after Bekenstein (2004) showed that lensing and Hubble expansion can be modeled rigourously in a Modified Relativity. The degeneracy of MOND with Dark Matter can be broken if we examine the non-linear MONDian Poisson's equation in detail. Here we study the effect of tides for a binary stellar system or a baryonic satellite-host galaxy system. We show that the Roche lobe is more squashed than the Newtonian case due to the anisotropic dilation effect in deep-MOND. We prove analytically that the Roche lobe volume scales linearly with the "true" baryonic mass ratio in both Newtonian and deep-MOND regimes, insensitive to the modification to the inertia mass. Hence accurate Roche radii of satellites can break the degeneracy of MOND and dark matter theory. Globular clusters and dwarf galaxies of comparable luminosities and distances show a factor of ten scatter in limiting radii; this is difficult to explain in any "mass-tracing-light" universe. The results here are generalizable to the intermediate MOND regime for a wide class of gravity modification function $\mu(g)$ (Zhao and Tian, astro-ph/0511754).

We argue that the Milky Way (MW) contains thousands intermediate mass black holes (IMBHs) and minihalos with a fraction of IMBHs still being enshrouded in extremely dense mini-spikes of dark matter (DM) particles. Each containing $10^6\msun$ of dark particles and no baryons in a sphere of 50 pc radius, the minihalos are dense enough to survive the Milky Way tide with the nearest minihalo within 2 kpc from the Sun. The IMBH is formed off-centre in a minihalo by gas accretion, and its growth adiabatically compresses a finite density of surrounding dark matter into a $r^{-1.5}$ mini-spike. Some IMBHs recentre on their minihalos after dynamical friction, and some IMBHs are ejected by birth kicks. Detectable by GLAST, the mini-spikes and minihalos should stand out the background and dominate the neutralino annihilation in the smooth MW and satellite galaxies. If they are the unidentified EGRET sources, upper limits can be set on the branching ratio of neutralino annihilation. The supermassive BH of the MW, if orginates from an IMBH, is also likely enshrouded with a mini-spike.

Many past attempts to kill MOND have only strengthened the theory. Better data on galaxy velocity curves clearly favor MOND (without fine-tuning) over cold dark matter. The usual critism on the incompleteness of classical MOND has spurred a Modified Relativity (MR) by Bekenstein. After outlining cosmology and lensing in MOND, we review MOND on small scales. We point out some potential problems of MOND in two-body relaxation and tidal truncation. We argue that the tidal field in any MOND-like gravity theory predicts that the Roche lobe sizes of a binary system are simply proportional to the binary baryonic mass ratio to the power 1/3. An immediate application of this result is that the tidal field and tidal truncation radii of million-star globular clusters and million-star dwarf galaxies (e.g., the Milky Way satellites NGC2419 and Carina) would be very similar because of the one-to-one relation between gravity and baryon distribution. This prediction appears, however, inconsistent with the fact that {\it all} globulars are truncated to much smaller sizes than {\it all} dwarf galaxies. Whether tide is uniquely determined by baryons can also be used to falsify any MOND-like gravity theory, whether classical or relativistic.

The promising neutralino dark matter particles generically condense into numerous earth-mass dark halos with smaller substructures suppressed by free streaming\cite{berezinsky}. The recent {\it Letter to Nature}\cite{moore} claims that these 0.01pc-sized dense halos emerged at redshifts 60--26 are rarely destructed inside galaxies, hence the nearest halos at about 0.1pc from Earth are bright in annihilation-powered gamma-rays, but are inaccessible to direct detections. However, most mini-halos reaching solar neighbourhood should experience strong impulses by individual stars in the Galactic disk and bulge, and have been torn into pc-long tidal streams over a Hubble time with only modest overdensity, reducing indirect annihilation signals. Sweeping across the solar system per century, mini-streams leave directional and temporal signatures for direct searches of neutralinos.

How fast a satellite decays its orbit depends on how slowly its mass is lost by tide. Motivated by inner halo satellite remnants like the Sgr and Omega Cen, we develop fully analytical models to study the orbital decay and tidal massloss of satellites. The orbital decay rate is often severely overestimated if applying the ChandraSekhar's formula without correcting for (a) the evaporation and tidal loss of the satellite and (b) the contraction of satellite orbits due to adiabatic growth of the host galaxy potential over the Hubble time. As a satellite migrates inwards, the increasing halo density affects the dynamical friction in two opposite ways:(1) it boosts the number of halo particles swept in the satellite's gravitational "wake", hence increasing the drag on the satellite, and (2) it boosts the tide which "peels off" the satellite, and reduces the amplitude of the wake. These competing processes can be modeled analytically for a satellite with the help of an empirical formula for the massloss history. The analytical model agrees with more numerical simulations of tidal massloss and dynamical friction. Rapid massloss due to increasing tides at ever smaller radius makes it less likely for streams or remnants of infalling satellites to intrude the inner halo (as the Sgr stream and oCen) than to stay in the outer halo (Magellanic stream), hence any intermediate-mass central black holes of the satellites are also likely "hung up" at large distances as well. To avoid "hung-up" the satellites might start with (i) pericenters much smaller than the typical distances to present-day observed satellites and with (ii) central density much higher than in the often seen finite density cores of observed satellites.

A fully analytical formulation is developed to make dynamical friction modeling more realistic. The rate for a satellite to decay its orbit in a host galaxy halo is often severely overestimated when applying ChandraSekhar's formula without correcting for the tidal loss of the satellite and the adiabactic growth of the host galaxy potential over the Hubble time. As a satellite decays to the inner and denser region of the host galaxy, the high ambient density speeds up their exchange of energy and angular momentum, but shrinks the Roche lobe of the satellite by tides. Eventually both processes of orbital decay and tidal stripping hang up altogether once the satellite is light enough. These competing processes can be modeled analytically for a satellite if we parametrize the massloss history by an empirical formula. We also incorporate the adiabatic contraction of orbits due to growth of the host potential well. Observed dwarf galaxies often show a finite density core, which determines how much inwards its remnants could be delivered to the host galaxy. Ghost streams or remnant cores or globular clusters of satellites should populate preferentially the outer halo (e.g., the Magellanic stream and the Ursa Minor dwarf), rather than the inner halo (e.g., the Sgr stream and the Omega-cen cluster). Massloss due to strong tides in the inner galaxy also makes it problematic for any small central black holes in nuclei of satellite galaxies to decay in orbit and merge into the center of the host galaxy.

Although a decade of microlensing searches towards the Large Magellanic Cloud (LMC) has detected 13-25 possible microlensing events, the nature and the location of the lenses, being either halo machos or LMC stars, remains a subject of debate. The star-star lensing models generically predict the existence of a small population (more than about 5 percent of stars with a spatial and kinematic distribution different from the thin, young disc of the LMC. Here we present the results of a large spectroscopic survey of the LMC, consisting of more than 1300 radial velocities measured accurately with the 2dF instrument. In this large sample, no evidence is found for any extraneous population over the expected LMC and Galactic components. Any additional, kinematically distint, population can only be present at less than the 1 percent level. We discuss the significance of this finding for the LMC self-lensing models.

Near- and mid-IR survey data from DENIS and ISOGAL are used to investigate the structure and formation history of the inner 10 degree (1.4 kpc) of the Milky Way galaxy. Synthetic bolometric corrections and extinction coefficients in the near- and mid-IR are derived for stars of different spectral types, to allow the transformation of theoretical isochrones into observable colour-magnitude diagrams. The observed IR colour-magnitude diagrams are used to derive the extinction, metallicity and age for individual stars. The inner galaxy is dominated by an old population (> 7 Gyr). In addition, an intermediate-age population (200 Myr to 7 Gyr) is detected, which is consistent with the presence of a few hundred Asymptotic Giant Branch stars with heavy mass loss. Furthermore, young stars (< 200 Myr) are found across the inner Bulge. The metallicities of these stellar population components are discussed. These results can be interpreted in terms of an early epoch of intense star formation and chemical enrichment which shaped the bulk of the Bulge and nucleus, and a more continuous star formation history which gradually shaped the disk from the accretion of sub-solar metallicity gas from the halo. A possible increase in star formation about 200 Myr ago might have been triggered by a minor merger. Ever since the formation of the first stars, mechanisms have been at play that mix the populations from the nucleus, Bulge and disk. Luminosity functions across the inner galactic plane indicate the presence of an inclined (bar) structure at > 1 kpc from the galactic centre, near the inner Lindblad resonance. The innermost part of the Bulge, within about 1 kpc from the galactic centre, seems azimuthally symmetric.

Isothermal models and other simple smooth models of dark matter halos of gravitational lenses often predict a dimensionless time delay $H_0\Delta t$ much too small to be comfortable with the observed time delays $\Delta t$ and the widely accepted $H_0$ value ($\sim 70$ km/s/Mpc). This conflict or crisis of the CDM has been highlighted by several recent papers of Kochanek, who claims that the standard value of $H_0$ favors a strangely small halo as compact as the stellar light distribution with an overall nearly Keplerian rotation curve. In an earlier paper (Paper I, astro-ph/0209191) we argue that this is not necessarily the case, at least in a perfectly symmetrical Einstein cross system. Here we introduce a {\it new mass degeneracy} of lens systems to give a counter example to Kochanek's claims. We fit the time delay and image positions in the quadruple image system PG1115+080. Equally good fits are found between lens models with flat vs. Keplerian rotation curves. Time delays in both types of models can be fit with the standard value of $H_0$. We demonstrate that it may still be problematic to constrain the size of lens dark halos even if the data image positions are accurately given and the cosmology is precisely specified.

A degeneracy in strong lens model is shown analytically. The observed time delays and quasar image positions might {\it not} uniquely determine the concentration and the extent of the lens galaxy halo mass distribution. Simply hardwiring the Hubble constant ($H_0$) and the cosmology ($\Omega, \Lambda$) to the standard values of $\Lambda$CDM cosmology might {\it not} fully lift this degeneracy, which exists rigourously even with very accurate data. Equally good fits to the images could be found in lens mass models with either a mostly Keplerian or a flat rotation curve. This degeneracy in mass models makes the task of getting reliable $H_0$ and $\Lambda$ from strong lenses even more daunting.

We argue that any violent galactic winds following early epoch of star bursts would significantly weaken the potentials of galaxies, and leave lasting signatures such as a lowered dark halo density and preferentially radial/escaping orbits for halo tracers such as globular clusters. A galaxy is disintegrated if more than half of its dynamical mass is blown off. The presence of dense halos in galaxies and the absence of intergalactic/escaping globulars should imply an upper limit on the amount of baryons lost in galactic winds around 4% of the total mass of the galaxy. This translates to limits on the baryons participating the early star bursts and baryons locked in stellar remnents, such as white dwarfs. The amount of halo white dwarfs claimed in recent proper motion searches and microlensing observations in the Galactic halo are too high to be consistent with our dynamical upper limits. Similar arguments also imply upper limits for the amount of neutron stars and stellar black holes, in galaxy halos. Nevertheless, a milder outflow is desirable, especially in dwarf galaxies, both for lowering their cold dark matter central density and for injecting metals to the intergalactic medium.

We estimate the rate at which stars are captured by supermassive black holes (BHs) in the centres of bulges and elliptical galaxies assuming that these initially had an isothermal cusp (rho ~ r^{-2} with velocity dispersion sigma_*). If radial orbits can be adequately replenished a BH can grow significantly by capturing stars and particles of the dark halo; the latter could be up to 20--40 percent of the total mass captured. A tight relation M_bh 10^8 (sigma_*/200\kms)^5 between BH hole mass and stellar velocity dispersion arises at the high mass end M_bh >= 10^8 M_sun if these giant BHs grow primarily by swallowing stars. For smaller masses a gradual flattening of the relation and an increase of the scatter are expected due to tidal disruption of stars. At moderate redshifts when bright galaxies may still possess isothermal cusps tidal disruption flares may be as frequent as 10^{-4} -- 10^{-2} yr^{-1} per galaxy, a few percent of the supernovae rate. Efficient randomization of the orbits can occur when the dense cores of infalling satellites sink to the sphere of influence of the BH by dynamical friction. Efficient replenishment of the loss cone of radial orbits can also explain why the supermassive binary BHs expected in hierarchically merging galaxies do generally coalesce, as suggested by the observed relation between BH mass and the inferred mass of stars ejected from an isothermal cusp.

Using current observational data and simple dynamical modeling, we demonstrate that Omega Cen is not special among the Galactic globular clusters in its ability to produce and retain the heavy elements dispersed in the AGB phase of stellar evolution. The multiple stellar populations observed in Omega Cen cannot be explained if it had formed as an isolated star cluster. The formation within a progenitor galaxy of the Milky Way is more likely, although the unique properties of Omega Cen still remain a mystery.

The observed rotation curves of dark matter-dominated dwarf galaxies indicate low density cores, contrary to the predictions of CDM models. A possible solution of this problem involves stellar feedback. A strong baryonic wind driven by vigorous star formation can remove a large fraction of the gas, causing the dark matter to expand. We explore the maximum effect of the feedback with an instantaneous removal of the gaseous disk. The energy input depends on the compactness of the disk, hence the specific angular momentum of the disk. For the plausible cosmological parameters and a large range of the disk angular momentum, the feedback is insufficient to destroy the central halo cusp, while the inner density is lowered only by a modest factor of 2 to 6. Any realistic modeling of the feedback would have even lesser impact on dark matter. We find that no star formation effect can resolve the problems of CDM cusps.

We demonstrate that the flux ratios of 4-image lensed quasars provide a powerful means of probing the small scale structure of Dark Matter (DM) halos. A family of smooth lens models can precisely predict certain combinations of flux ratios using only the positions of the images and lens as inputs. Using 5 observed lens systems we show that real galaxies {\it cannot} be described by smooth singular isothermal ellipsoids, nor by the more general elliptical power-law potentials. Large scale distortions from the elliptical models can not yet be ruled out. Nevertheless we find, by comparing with simulations, that the data can be accounted for by a significant ($\sim 5-10%$) amount of dark substructures within a projected distance of several kpc from the center of lenses. This interpretation favors the Cold Dark Matter (CDM) model other the warm or self-interacting DM models.

We simulate the result of orbital decay for a dwarf satellite galaxy in a disk galaxy potential. The peeling-off of the satellite by the strong tidal field within the disk of the host galaxy quickly reduces the satellite's mass, and drastically slows down the decaying of the orbit. The result is generally a cluster on a large halo orbit, like the M54 cluster of the Sagittarius galaxy in the Milky Way and the massive G1 cluster in outer halo of M31. It is generally challenging for this mechanism to deliver globular cluster-mass objects to the observed very close-in orbit of Omega Cen unless the dwarf galaxy was born within a few scale lengths of the Galactic disk.

The failure of standard cosmolocical models in accounting for the statistics of dwarf satellites and the rotation curve of gas-rich dwarf galaxies in detail has led us to examine whether earlier non-equilibrium models of dwarf spheroidal satellites without any dark matter should be reconsidered in more detail. Such models can explain the high dispersion of the dwarf spheroids by the projection of disrupted tidal debris. We show in the case of Milky Way satellites, that these models predict a significant spread in the apparent magnitude of horizontal branch stars which is correlated with sky position and velocity. In particular, the models produce a strong correlation of radial velocity with the long axis of the dwarf. Current data do not set strong enough constraints on the models, but we suggest that photometric and spectroscopic surveys of extra-tidal stars of nearby dwarf spheroids in the Milky Way and Andromeda can falsify these models without dark matter.

We apply standard disk formation theory with adiabatic contraction within cuspy halo models predicted by the standard LCDM cosmology. The resulting models score remarkably well when confronted with the broad range of observational data available for the Milky Way and M31 galaxies, giving a Milky Way virial mass of 1-2x10^12Msun and concentration C=12-17. We consider two types of models, in which: (A) baryons conserve angular momentum and (B) some of the angular momentum of the baryons is transferred to the dark matter. Type-A models produce good agreement with observed rotation curves and obey constraints in the solar neighborhood, but may have too much dark matter in the center to allow a fast rotating bar. The type-B models with angular momentum transport have a slightly more massive disk and less dark matter in the central part, allowing a fast rotating bar to persist. Both classes of models probably have sufficient baryonic mass in the central 3.5kpc to reproduce recent observational values of the optical depth to microlensing events towards the Galactic center. All models require that about 1/2 of all baryons expected inside the virial radius must not be in the disk or bulge. We investigate whether the range of virial masses allowed by our dynamical models is compatible with constraints from the galaxy luminosity function, and find a range of parameter space that is allowed by this constraint. We conclude that rotation curves and dynamical properties of ``normal'' high surface brightness spiral galaxies appear to be consistent with standard LCDM.

The past growth of the central black hole (BH) might have enhanced the density of cold dark matter halo particles at the Galactic center. We compute this effect in realistic growth models of the present (2-3)*10**6 solar mass BH from a low-mass seed BH, with special attention to dynamical modeling in a realistic galaxy environment with merger and orbital decay of a seed BH formed generally outside the exact center of the halo. An intriguing ``very-dense spike'' of dark matter has been claimed in models of Gondolo and Silk with density high enough to contradict with experimental upper bounds of neutralino annihilation radiation. This ``spike'' disappears completely or is greatly weakened when we include important dynamical processes neglected in their idealized/restrictive picture with cold particles surrounding an at-the-center zero-seed adiabaticly-growing BH. For the seed BH to spiral in and settle to the center within a Hubble time by dynamical friction, the seed mass must be at least a significant fraction of the present BH. Any subsequent at-the-center growth of the BH and steepening of the central Keplerian potential well can squeeze the halo density distribution only mildly, whether the squeezing happens adiabatically or instantaneously.

The off-centered "bar" of the Large Magellanic Cloud (LMC) is commonly assumed, for simplicity's sake, to lie in the same plane as the LMC disk. Relaxing this assumption to the extent allowed by the observations, we can drastically change our interpretation of the microlensing events towards the LMC. We propose that the off-centered ``bar'' in the Large Magellanic Cloud is a foreground object slightly misaligned with, and offset from, the plane of the LMC disk. The small displacement and misalignment are consequences of recent tidal interactions with the Small Magellanic Cloud and the Galaxy. This proposal, though radical, is consistent with the kinematics of the LMC and near-infrared star count maps from the DENIS and 2MASS surveys. It does not violate any of the observational limits on the depth and structure of the LMC -- in particular, the reported 25 to 50 degree inclination range of the LMC and the east-west gradient of distance moduli of standard candles. Contributions to LMC microlensing come from the mutual lensing of stars in the disk and the ``bar'', as well as self-lensing of the disk and the ``bar''. The microlensing optical depth of such configurations from self-lensing alone is comparable to the observed microlensing signal reported by the EROS and MACHO groups. Possible observations are suggested to discriminate between our misaligned, offset ``bar'' model and the conventional picture of an off-centered, planar bar.

While the Hubble constant can be derived from observable time delays between images of lensed quasars, the result is often highly sensitive to assumptions and systematic uncertainties in the lensing model. Unlike most previous authors we explore a broad class of models with non-parametric lens radial profile and allow for non-elliptical lens potentials. The axis ratio and position angle of the lens can be determined from the image positions of quadruple gravitational lensed systems by solving the lens equation semi-analytically, independent of the radial profile. We give simple equations for estimating the power-law slope of the lens density directly from the image positions and for estimating the time delay ratios. Our method greatly simplifies the numerics for fitting observations and is fast in exploring the model parameter space. As an illustration we apply the model to PG1115+080. We show that the measured image positions and time delays do not uniquely determine the Hubble constant. The mass-to-light ratio of the lens galaxy increases with radius, higher at image C than at image B.

Evidence for a bar at the center of the Milky Way triggered a renewed enthusiasm for dynamical modelling of the Galactic bar-bulge. Our goal is to compare the kinematics of a sample of tracers, planetary nebulae, widely distributed over the bulge with the corresponding kinematics for a range of models of the inner Galaxy. Three of these models are N-body barred systems arising from the instabilities of a stellar disk (Sellwood, Fux and Kalnajs), and one is a Schwarzschild system constructed to represent the 3D distribution of the COBE/DIRBE near-IR light and then evolved as an N-body system for a few dynamical times (Zhao). For the comparison of our data with the models, we use a new technique developed by Saha (1998). The procedure finds the parameters of each model, i.e. the solar galactocentric distance R_o in model units, the orientation angle phi, the velocity scale (in km/s per model unit), and the solar tangential velocity which best fit the data.

There are lines of evidence suggesting that some of the observed microlensing events in the direction of the Large Magellanic Cloud (LMC) are caused by ordinary star lenses as opposed to dark Machos in the Galactic halo. Efficient lensing by ordinary stars generally requires the presence of one or more additional concentrations of stars along the line of sight to the LMC disk. If such a population behind the LMC disk exists, then the source stars (for lensing by LMC disk objects) will be drawn preferentially from the background population and will show systematic differences from LMC field stars. One such difference is that the (lensed) source stars will be farther away than the average LMC field stars, and this should be reflected in their apparent baseline magnitudes. We focus on red clump stars: these should appear in the color-magnitude diagram at a few tenths of a magnitude fainter than the field red clump. Suggestively, one of the two near-clump confirmed events, MACHO-LMC-1, is a few tenths of magnitude fainter than the clump.

We compute the predicted optical depth and duration distribution of microlensing events towards Baade's window in a model composed of a Galactic disk and a bar. The bar model is a self-consistent dynamical model built out of individual orbits that has been populated to be consistent with the COBE maps of the Galaxy and kinematic observations of the Bulge. We find most of the lenses are at 6.25 kpc in BW. The microlensing optical depth of a 2E10 Msun bar plus a truncated disk is (2.2+-0.3)E-6, consistent with the very large optical depth (3.2+-1.2)E-6 found by Udalski et al. (1994). This prediction doubles that of axisymmetric models by Kiraga and Paczynski (1994) since the bar is elongated along the line-of-sight. The large Einstein radius and small transverse velocity dispersion also predict a longer event duration in the bar model than the KP model. The event rate and duration distribution also depend on the lower mass cutoff of the lens mass function. With a 0.1Msun cutoff, 5-7 events of typical 20 days duration are expected for the OGLE experiment according to our model, while Udalski et al. observed 9 events with durations from 8-62 days. On the other hand, if most of the lenses are brown dwarfs, our model predicts too many short duration events; a KS test finds only 7% probability for the model with 0.01Msun cutoff to be consistent with current data.

We have extended the spectral dynamics formalism introduced by Binney & Spergel, and have implemented a semi-analytic method to represent regular orbits in any potential, making full use of their regularity. We use the spectral analysis code of Carpintero & Aguilar to determine the nature of an orbit (irregular, regular, resonant, periodic) from a short-time numerical integration. If the orbit is regular, we approximate it by a truncated Fourier time series of a few tens of terms per coordinate. Switching to a description in action-angle variables, this corresponds to a reconstruction of the underlying invariant torus. We then relate the uniform distribution of a regular orbit on its torus to the non-uniform distribution in the space of observables by a simple Jacobian transformation between the two sets of coordinates. This allows us to compute, in a cell-independent way, all the physical quantities needed in the study of the orbit, including the density and in the line-of-sight velocity distribution, with much increased accuracy. The resulting flexibility in the determination of the orbital properties, and the drastic reduction of storage space for the orbit library, provide a significant improvement in the practical application of Schwarzschild's orbit superposition method for constructing galaxy models. We test and apply our method to two-dimensional orbits in elongated discs, and to the meridional motion in axisymmetric potentials, and show that for a given accuracy, the spectral dynamics formalism requires an order of magnitude fewer computations than the more traditional approaches.

Deriving the 3-dimensional volume density distribution from a 2-dimensional light distribution of a system yields generally non-unique results. The case for nearby dust-free systems is studied, taking into account the extra constraints from the perspective effect. It is shown analytically that a new form of non-uniqueness exists. We can design a Phantom Spheroid (PS) for a nearby system which preserves the intrinsic mirror symmetry and projected surface brightness of the system while changing the shape and the major-axis orientation of the system. A family of analytical models are given as functions of the distance ($D_0$) to the object and the amount ($\gamma$) of the superimposed PS density. The PS density, different from the well-known konuses of extragalactic systems, makes the luminosity of the system vary slightly with the distance $D_0$. The physical ranges for $\gamma$ and the major axis angles are constrained analytically by requiring a positive volume density everywhere. These models suggest that observations other than surface brightness maps are required to lift the degeneracy in the tilt angles and axis ratio of the central bar of the Milky Way.

We propose an observational test that can break the degeneracy of two main classes of microlensing models to the Magellanic Clouds: (a) the lenses are located in the Galactic halo, and (b) the lenses are located in the LMC disk. The source stars in the latter (self-lensing) models tend to be at the far side or behind the LMC disk, thus experience more reddening and extinction by dust in the LMC disk than ordinary stars in a nearby line of sight. Clearly such bias would not occur in the macho halo lensing models. We show that this reddening effect is at a level readily observable for the present 30 or so microlensing alerts fields, either with multi-band photometry from a good seeing site, or more definitively with ultra-violet (UV) spectroscopy with Space Telescope Imaging Spectrograph (STIS). Stars behind the LMC dust layer should stand out as UV-faint objects (by more than 1 mag than average stars in the LMC). HST can also resolve numerous faint neighbouring stars within a few arcsecs of a lensed source, hence remove blending in these crowded regions and build a reddening map to control the patchiness of dust.

We study the effect of reddening on microlensed stars towards the LMC. If lenses are in the LMC disk, then the source stars should be at the far side or behind the LMC disk. Thus they should experience more reddening and extinction by dust in the LMC disk than typical stars in the immediate neighbouring lines of sight. We simulate this effect in a variety of models for the LMC stars and dust. We stress that the optical depth is not one number, but a function of the reddening of the survey stars. We discuss how these effects could be used to constrain the fraction of machos in the dark halo. The effect of patchiness of dust can be controled by working with faint stars in the smallest patches of the sky around the microlensed stars. This can be done most effectively with the Hubble Space Telescope in the ultra-violet. The non-detection of the reddening bias would be strongly in favor of machos in the Galactic halo.

ESA's Global Astrometric Interferometer for Astrophysics (GAIA) holds the promise of mapping out the detailed phase space structure of the Galactic halo by providing unprecedented annual proper motion and parallax of $1-10\muas$ astrometric accuracy (Gilmore et al. 1998). Unlike NASA's Space Interferometry Mission (SIM), which will achieve similar accuracies but is a pointed instrument, GAIA will be able to construct a global catalogue of the halo. Here we study proper motions of giant branch stars in a tidal debris torn from a small satellite system in the halo. We follow the evolution of a cold stream on a polar orbit between 8-50 kpc in a variety of histories of the Galactic potential, and observe the bright ($V<18$mag) members of the debris tail with GAIA accuracy. We simulate effects due to the growing or flipping of the Galactic disk over the past 4 Gyrs or the perturbation from a massive accreted lump such as the progenitor of the Magellanic Clouds. Our simulations suggest that the results of Johnston, Zhao, Spergel & Hernquist (1999) and Helmi, Zhao & de Zeeuw (1999) for static Galactic potentials are likely to be largely generalizable to realistic time-dependent potentials: a tidal debris remains cold in spite of evolution and non-axial symmetry of the potential. GAIA proper motion measurements of debris stars might be used to probe both Galactic structure and Galactic history. We also study several other factors influencing our ability to identify streams, including accuracy of radial velocity and parallax data from GAIA, and contamination from random field stars. We conclude that nearby, cold streams could be detected with GAIA if these cousins of the Sagittarius stream exist.

The shape of a galaxy is constrained both by mechanisms of formation (dissipational vs. dissipationless) and by the available orbit families (the shape and amount of regular and stochastic orbits). It is shown that, despite the often very flattened shapes of banana and fish orbits, these boxlet orbits generally do not fit a triaxial galaxy in detail because, similar to loop orbits, they spend too little time at the major axis of the model density distribution. This constraint from the shape of fish orbits is relaxed at (large) radii where the density profile of a galaxy is steep.

We propose a radial velocity survey with the aim to resolve the current dispute on the LMC lensing: in the pro-macho hypothesis the lenses are halo white dwarfs or machos in general; in the pro-star hypothesis both the lenses and the sources are stars in various observed or hypothesized structures of the Magellanic Clouds and the Galaxy. Star-star lensing should prefer sources at the backside or behind the LMC disc because lensing is most efficient if the source is located a few kpc behind a dense screen of stars, here the thin disc of the LMC. This signature of self-lensing can be looked for by a radial velocity survey since kinematics of the stars at the back can be markedly different from that of the majority of stars in the cold, rapidly rotating disc of the LMC. Detailed simulations of effect together with optimal strategies of carrying out the proposed survey are reported here. Assuming that the existing 30 or so alerted stars in the LMC are truely microlensed stars, their kinematics can test the two lensing scenarios; the confidence level varies with the still very uncertain structure of the LMC. Spectroscopy of the existing sample and future events requires about two or three good-seeing nights per year at a 4m-8m class southern telescope, either during the amplification phase or long after.

If the Universe has a significant baryonic dark component in the form of compact objects in galaxy halos (machos), then there is a minute chance (about $10^{-7}$) that one of the Galactic machos passes sufficiently close to our line of sight to a star out of some $10^7$ monitored stars in the Magellanic Clouds (MCs) that it brightens by more than 0.3 magnitude due to gravitational focusing. After a brief discussion of the current controversy over the interpretation of the observed events, i.e., whether the lensing is caused by halo white dwarfs or machos in general or by stars in various observed or hypothesized structures of the Clouds and the Galaxy, I propose a few observations to put ideas of the pro-macho camp and the pro-star camp to test. In particular, I propose a radial velocity survey towards the MCs.

Tremendous progress has been made recently in modelling the morphology and kinematics of centers of galaxies. Increasingly realistic models are built for central bar, bulge, nucleus and black hole of galaxies, including our own. The newly revived Schwarzschild method has played a central role in these theoretical modellings. Here I will highlight some recent work at Leiden on extending the Schwarzschild method in a few directions. After a brief discussion of (i) an analytical approach to include stochastic orbits (Zhao 1996), and (ii) the ``pendulum effect'' of loop and boxlet orbits (Zhao, Carollo, de Zeeuw 1999), I will concentrate on the very promising (iii) spectral dynamics method, with which not only can one obtain semi-analytically the actions of individual orbits as previously known, but also many other physical quantities, such as the density in configuration space and the line-of-sight velocity distribution of a superposition of orbits (Copin, Zhao & de Zeeuw 1999). The latter method also represents a drastic reduction of storage space for the orbit library and an increase in accuracy over the grid-based Schwarzschild method.

The shape of a cusped galactic nucleus is constrained by the range of shapes of orbits in its gravitational potential. It is shown for scale-free non-axisymmetric discs that while a plausible elongated density model requires at least some orbits to spend more time near the major axis than anywhere else, both regular boxlets and tube orbits generally cross the major axis too fast for self-consistency. If galaxies host inner nuclear discs or flat bars with a cuspy surface light profile ($\gamma=|{d \log \mu / d \log r}|>0$), their ellipticity $1-{b / a}$ cannot be greater than about ${\gamma / 2}$. Discs or bars with a shallow central profile ($\gamma \le 0.3$) should not be strongly elliptical.

Future astrometric satellites, such as SIM (NASA's Space Interferometric Mission) and GAIA (ESA's Global Astrometric Interferometer for Astrophysics), hold the promise of mapping out the detailed phase space structure of the Galactic halo by providing unprecedented annual proper motion and parallax of $1-10\mu$as astrometric accuracy. Here we show that proper motions of hundred or so giant branch stars in a tidal debris torn from a small satellite (a $10^{5-7}L_\odot$ Galactic dwarf galaxy or globular cluster) in the halo is sensitive to the current Galactic potential and its past evolution. We follow the evolution of a cold (velocity dispersion of 10 km/s) stream on a nearby (between 8-50 kpc) polar orbit in a variety of histories of the potential of the Galaxy, and observe the bright ($V<18$mag) members of the debris tail with GAIA accuracy. We simulate effects due to the growing or flipping of the Galactic disk over the past 4 Gyrs or the perturbation from a massive accreted lump such as the progenitor of the Magellanic Clouds. We study various factors influencing our ability to identify streams, including contamination from field stars, accuracy of radial velocity and distance data and evolution and non-axial symmetry of the potential. Our simulations suggest that nearby, cold streams could be detected with GAIA if these cousins of the Sagittarius stream exist. Results of Johnston, Zhao, Spergel & Hernquist (1999) and Helmi, Zhao & de Zeeuw (1999) for static Galactic potentials are likely to be largely generalizable to moderately time-evolving potentials. SIM and GAIA measurements of debris stars might be used to probe both Galactic structure and Galactic history.

We explore the use of tidal streams from Galactic satellites to recover the potential of the Milky Way. Our study is motivated both by the discovery of the first lengthy stellar stream in the halo (\cite{it98}) and by the prospect of measuring proper motions of stars brighter than 20th magnitude in such a stream with an accuracy of $\sim 4\mu as/$yr, as will be possible with the Space Interferometry Mission (SIM). We assume that the heliocentric radial velocities of these stars can be determined from supporting ground-based spectroscopic surveys, and that the mass and phase-space coordinates of the Galactic satellite with which they are associated will also be known to SIM accuracy. Using results from numerical simulations as trial data sets, we find that, if we assume the correct form for the Galactic potential, we can predict the distances to the stars as a consequence of the narrow distribution of energy expected along the streams. We develop an algorithm to evaluate the accuracy of any adopted potential by requiring that the satellite and stars recombine within a Galactic lifetime when their current phase-space coordinates are integrated backwards. When applied to a four-dimensional grid of triaxial logarithmic potentials, with varying circular velocities, axis ratios and orientation of the major-axis in the disk plane, the algorithm can recover the parameters used for the Milky Way in a simulated data set to within a few percent using only 100 stars in a tidal stream.

We investigate what the proposed ESA astrometric satellite GAIA will reveal by observing the halo of the Milky Way. Specifically, we look for halo streams which are the signatures left by the merging/accretion events experienced by a typical galaxy like the Milky Way. We run numerical simulations of the disruption of satellite galaxies in a Galactic potential to generate artificial GAIA halo catalogs. We recover the streams by searching for peaks in angular momentum space.

How has the ``fluffy'' core of the Sgr dwarf galaxy survived multiple strong shocks from the tidal force of the Galactic halo and disc since the formation of the core a Hubble time ago? A scenario that Sgr was deflected to its current orbit by the Magellanic Clouds after a rendezvous on the north Galactic pole 2-3 Gyrs ago is examined. It is shown that the conditions of the collision fix both the sense of circulation of Sgr and the LMC around the Galaxy and the slope of the Galactic rotation curve. The model argues that the two orthogonal polar circles traced by a dozen or so Galactic halo dwarf galaxies and globular clusters (LMC-SMC-Magellanic Stream-Draco-Ursa Minor along $l \approx 270^o$ and M54-Ter 7-Ter 8-Arp 2-NGC 2419-Pal 15 along $l \approx 0^o$) are streams of tidal relics from two ancient galaxies which was captured on two intersecting polar rosette orbits by the Galaxy. Our results favor the interpretation of microlensing towards the LMC being due to source or lens stars in tidal features of the Magellanic Clouds. We discuss direct and indirect observations to test the collision scenario.

With the steady increase of the sample size of observed microlenses towards the central regions of the Galaxy, the main source of the uncertainty in the lens mass will shift from the simple Poisson noise to the intrinsic non-uniqueness of our dynamical models of the inner Galaxy, particularly, the Galactic bar. We use a set of simple self-consistent bar models to investigate how the microlensing event rate varies as a function of axis ratio, bar angle and velocity distribution. The non-uniqueness of the velocity distribution of the bar model adds a significant uncertainty (by about a factor of 1.5) to any prediction of the lens mass. Kinematic data and self-consistent models are critical to lift the non-uniqueness. We discuss the implications of these results for the interpretation of microlensing observations of the Galactic bulge. In particular we show that Freeman bar models scaled to the mass of the Galactic bulge/bar imply a typical lens mass of around $0.8M_\odot$, a factor of 3-5 times larger than the value from other models.

We summarize recent work on the structure and dynamics of the Galactic bar and inner disk. Current work focusses on constructing a quantitative model which integrates NIR photometry, source count observations, gas kinematics, stellar dynamical observations, and microlensing. Some avenues for future research are discussed.

Increasing evidences suggest that the Galactic halo is lumpy on kpc scales due to the accretion of at least a dozen small galaxies (LMC/SMC, Sgr, Fornax etc.). Faint stars in such lumpy structures can significant microlense a background star with an optical depth $10^{-7}-10^{-6}$, which is comparable to the observed value to the LMC (Alcock et al. 1996c). The observed several microlensing events towards the LMC can be explained by a tidal debris tail lying in the Magellanic plane if the progenitor of the Magellanic Clouds and Stream and other satellite galaxies in the Magellanic plane has a mass about twice that of the disc of the LMC. The LMC stars can either lense stars in the debris tail behind the LMC, or be lensed by stars in the part of debris tail in front of the LMC. The models are consistent with an elementary particle dominated Galactic halo without massive compact halo objects (MACHOs). They also differ from Sahu's (1994) LMC-self-lensing model by predicting a higher optical depth and event rate and less concentration of events to the LMC center.

Galactic nuclei are now generally thought to have density cusps in their centres, and the evidence is mounting that as a consequence they are unlikely to be triaxial. Self-consistent stellar dynamical models of non-axisymmetric cusps would be an interesting counter-argument to this conclusion. We consider 2-d analogues of triaxial cusps: a sequence of non-axisymmetric, cuspy discs first described by Sridhar & Touma (1997). Scale-free models with potential {}$\Phi\propto r^\alpha$ are examined in detail. It is shown analytically for $0<\alpha \la 0.43$ that self-consistent models with positive phase-space density do not exist. Numerical solutions of the combined Vlasov and Poisson equations suggest that the whole sequence of models with $0<\alpha<1$ are also unphysical. Along with existing work on cusps, we conclude on purely theoretical grounds that galactic nuclei are not expected to be triaxial.

Non-uniqueness of deprojecting the integrated light distribution of a nearby or faraway triaxial body is reviewed in the context of deriving the volume density of the Galactic bar from the COBE/DIRBE maps of the Galactic plane. The exact origin of this non-uniqueness is studied. One can write down a sequence of triaxial bar models which appear identical in integrated light from the Sun's perspective, and the whole sequence is mapped out as a function of the Galactocentric distance of the observer from galactic to extragalactic distance scales. While mirror symmetries and perspective effects are compatible with any orientation of the bar in the sequence, weak upper and lower bounds can still be placed on the angles and axis ratios of the bar by positivity and other general requirements. Star count data of bulge giants are ideal for selecting a unique model from the sequence.

The discovery of the Sgr dwarf galaxy (Ibata, Gilmore, Irwin 1994) shows that the Galactic halo contains large clumps which are not full viralized. Stars in such a clump can lense a background extragalactic source. I use parameters of Sgr to demonstrate that microlensing by a disintegrated dwarf galaxy has a probability (optical depth) $10^{-7}-10^{-6}$, which is comparable to that of an isothermal Galactic dark halo of massive compact objects (MACHOs). This may change implications of observed microlensing events towards the LMC. If more than one such clumps are hidden in the Galactic halo, then given their enormous angular size ($\sim 5^o \times 20^o$), there is a significant chance for {\sl any} line of sight, say the LMC direction, to intersect one of them. As an alternative to the half-MACHO dark halo models (Alcock et al. 1996a,c), I discuss a non-MACHO model with lenses of the observed events being faint stars in a hidden Sgr-like disintegrated dwarf in the line-of-sight path to the LMC.

The microlensing probability (optical depth $\tau$) 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) $\tau$ is proportional to the mass of the bar, $M$. (2) $\tau$ 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 $\alpha>45$ degress. (4) $\tau$ is the largest if the angle $\alpha$ and the axis ratio $y_0/x_0$ conspires so that $y_0/x_0=\tan \alpha$. (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\sigma$ 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 $\tau$ argues for a massive ($> 2\times 10^{10}M_\odot$) boxy bar with $y_0/x0\approx \tan\alpha$ and $\alpha<20$ deg and with a flat radial profile up to corotation.

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.

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 $\gamma / \eta$-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.

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\sigma$ 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.

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.

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.

A self-consistent stellar dynamical model for the Galactic bar is constructed from about 500 numerically computed orbits with an extension of the Schwarzschild technique. The model fits the {\sl COBE} found asymmetric boxy light distribution and the observed stellar kinematics of the bulge. The model potential is also consistent with the non-circular motions of the HI and CO velocity maps of the inner Galaxy. We also use the stellar bar model to construct a N-body model to study stability and a microlensing map towards the bulge, which can account for the observed optical depth and the event duration by the MACHO and OGLE collaborations. The technique used here can be applied to interpret light and velocity data of external bulges/bars and galactic nuclei.

We study a sample of 62 Baade's Window (1,-4) K giants that have published proper motions, radial velocity and metallicity. We construct the velocity ellipsoids for the metal rich stars and metal poor stars. We find a vertex deviation characteristic of a non-axisymmetric system. The long axes of the velocity ellipsoids of two samples point to two nearly perpendicular directions (lv,bv)=(-65+-9,+14+-9) and (25+-14,-11+-14) respectively. The vertex deviations of the velocity ellipsoids cannot be consistently explained by ANY oblate model. We populate orbits in a realistic bar potential with a Gaussian velocity distribution, allowing us to simulate and interpret observations. We conclude that the data are consistent with a triaxial bulge pointing towards (l,b) with l<0 and b=0 as suggested by earlier work. We also predict that low latitude abs(b)<4 deg bulge fields should show the vertex deviation more strongly and would therefore be the best locations for future proper motion studies. The metal rich and metal poor stars appear to occupy the direct boxy and retrograde loop orbits. The correlations between the metallicity and the orbit families can develop if the Bulge forms dissipatively on a sufficiently long time scale. However, it is difficult to reconcile the existence of such correlations with a violent history for the potential of the inner Galaxy.