The observed line-of-sight velocity dispersion $\sigma_{los}$ of the ultra diffuse galaxy Dragonfly 44 (DF44) requires a Newtonian dynamical mass-to-light ratio of $M_{dyn}/L_I=26^{+7}_{-6}$ Solar units. This is well outside the acceptable limits of our stellar population synthesis (SPS) models, which we construct using the integrated galactic initial mass function (IGIMF) theory. Assuming DF44 is in isolation and using Jeans analysis, we calculate $\sigma_{los}$ profiles of DF44 in Milgromian dynamics (MOND) and modified gravity (MOG) theories without invoking dark matter. Comparing with the observed kinematics, the best-fitting MOND model has $M_{dyn}/L_I = 3.6_{-1.2}^{+1.6}$ and a constant orbital anisotropy of $\beta=-0.5_{-1.6}^{+0.4}$. In MOG, we first fix its two theoretical parameters $\alpha$ and $\mu$ based on previous fits to the observed rotation curve data of The HI Nearby Galaxy Survey (THINGS). The DF44 $\sigma_{los}$ profile is best fit with $M_{dyn}/L_I=7.4_{-1.4}^{+1.5}$, larger than plausible SPS values. MOG produces a $\sigma_{los}$ profile for DF44 with acceptable $M_{dyn}/L_I$ and isotropic orbits if $\alpha$ and $\mu$ are allowed to vary. MOND with the canonical $a_0$ can explain DF44 at the $2.40 \sigma$ confidence level (1.66\%) if considering both its observed kinematics and typical star formation histories in an IGIMF context. However, MOG is ruled out at $5.49 \sigma$ ($P$-value of $4.07 \times 10^{-8}$) if its free parameters are fixed at the highest values consistent with THINGS data.

Recent observations of rotationally supported galaxies show a tight correlation between the observed radial acceleration at every radius and the Newtonian acceleration generated by the baryonic mass distribution, the so-called radial acceleration relation (RAR). The rotation curves (RCs) of the SPARC sample of disk galaxies with different morphologies, masses, sizes and gas fractions are investigated in the context of modified Newtonian dynamics (MOND). We include the effect of cold dark baryons by scaling the measured mass in the atomic form by a factor of $c$ in the mass budget of galaxies. In addition to the standard interpolating function, we also fit the RCs and the RAR with the empirical RAR-inspired interpolating function. Slightly better fits for about $47\%$ of galaxies in our sample are achieved in the presence of dark baryons ($c>1$) with the mean value of $c = 2.4\pm 1.3$. Although the MOND fits are not significantly improved by including dark baryons, it results in a decrease in the characteristic acceleration $g_\dag$ by $40\%$. We find no correlation between the MOND critical acceleration $a_0$ and the central surface brightness of the stellar disk, $\mu_{3.6}$. This supports $a_0$ being a universal constant for all galaxies.

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.

We study the dynamical evolution of star clusters on eccentric orbits using a semi-analytical approach. In particular we adapt and extend the equations of EMACSS code, introduced by Gieles et al. (2014), to work with eccentric orbits. We follow the evolution of star clusters in terms of mass, half-mass radius, core radius, Jacobi radius and the total energy over their dissolution time. Moreover, we compare the results of our semi-analytical models against $N$-body computations of clusters with various initial half-mass radius, number of stars and orbital eccentricity to cover both tidally filling and under-filling systems. The evolution profiles of clusters obtained by our semi-analytical approach closely follow those of $N$-body simulations in different evolutionary phases of star clusters, from pre-collapse to post-collapse. Given that the average runtime of our semi-analytical models is significantly less than that of $N$-body models, our approach makes it feasible to study the evolution of large samples of globular clusters on eccentric orbits.

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.

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.

The observed present-day stellar mass function (PDMF) of the solar neighborhood is a mixture of stellar populations born in star-forming events that occurred over the life-time of the thin disk of the Galaxy. Assuming stars form in embedded clusters which have stellar initial mass functions (IMFs) which depend on the metallicity and density of the star-forming gas clumps, the integrated galaxy-wide IMF (IGIMF) can be calculated. The shape of the IGIMF thus depends on the SFR and metallicity. Here, the shape of the PDMF for stars more massive than $1\,M_\odot$ in combination with the mass density in low-mass stars is used to constrain the current star-formation rate (SFR), the star formation history (SFH) and the current stellar plus remnant mass ($M_*$) in the Galactic thin disk. This yields the current SFR, $\dot{M}_*= 4.1^{+3.1}_{-2.8}~M_\odot$yr$^{-1}$, a declining SFH and $M_*=2.1^{+3.0}_{-1.5}\times 10^{11}M_\odot$, respectively, with a V-band stellar mass-to-light ratio of $M_*/L_V=2.79^{+0.48}_{-0.38}$.These values are consistent with independent measurements. We also quantify the surface density of black holes and neutron stars in the Galactic thin disk. The invariant canonical IMF can reproduce the PDMF of the Galaxy as well as the IGIMF, but in the universal IMF framework it is not possible to constrain any of the above Galactic properties. Assuming the IGMF theory is the correct framework and in combination with the vertical velocity dispersion data of stars, it follows that the Milky Way would have appeared as a chain galaxy at high redshift.

The stellar initial mass function (IMF) is commonly assumed to be an invariant probability density distribution function of initial stellar masses being represented by the canonical IMF. As a consequence the galaxy-wide IMF (gwIMF), defined as the sum of the IMFs of all star forming regions, should also be invariant. Recent observational and theoretical results challenge the hypothesis that the gwIMF is invariant. In order to study the possible reasons for this variation we use the IMF determined in resolved star clusters and apply the IGIMF-theory to calculate a grid of gwIMF models for metallicities, -3<[Fe/H]<1, and galaxy-wide star formation rates, $10^{-5}$<SFR<$10^{5}\,\mathrm{M_{\odot}/yr}$. For a galaxy with metallicy [Fe/H]$<0$ and SFR$\,> 1\,M_\odot$/yr, which is a common condition in the early Universe, we find that the gwIMF is top-heavy (more massive stars), when compared to the canonical IMF. For a SFR $< 1\,\mathrm{M_{\odot}/yr}$ the gwIMF becomes top-light regardless of the metallicity. For metallicities $\mathrm{[Fe/H]} > 0$ the gwIMF can become bottom-heavy regardless of the SFR. The IGIMF models predict that massive elliptical galaxies should have formed with a gwIMF that is top-heavy within the first few hundred Myr of the galaxy's life and that it evolves into a bottom-heavy gwIMF in the metal-enriched galactic center. We study the SFR$-$H$\alpha$ relation, its dependency on metallicity and the SFR, the correction factors to the Kennicutt SFR$_{\rm K}-$H$\alpha$ relation, and provide new fitting functions Late-type dwarf galaxies show significantly higher SFRs with respect to Kennicutt SFRs, while star forming massive galaxies have significantly lower SFRs than hitherto thought. This has implications for the gas-consumption time scales and for the main sequence of galaxies. The Leo P and ultra-faint dwarf galaxies are discussed explicitly. [abridged]

Mass models of a sample of 171 low- and high-surface brightness galaxies are presented in the context of the cold dark matter (CDM) theory using the NFW dark matter halo density distribution to extract a new concentration-viral mass relation ($c-M_{vir}$). The rotation curves (RCs) are calculated from the total baryonic matter based on the 3.6 $\mu m $-band surface photometry, the observed distribution of neutral hydrogen, and the dark halo, in which the three adjustable parameters are the stellar mass-to-light ratio, halo concentration and virial mass. Although accounting for a NFW dark halo profile can explain rotation curve observations, the implied $c-M_{vir}$ relation from RC analysis strongly disagrees with that resulting from different cosmological simulations. Also, the $M/L-$color correlation of the studied galaxies is inconsistent with that expected from stellar population synthesis models with different stellar initial mass functions. Moreover, we show that the best-fitting stellar $M/L-$ ratios of 51 galaxies (30\% of our sample) have unphysically negative values in the framework of the $\Lambda $CDM theory. This can be interpreted as a serious crisis for this theory. This suggests either that the commonly used NFW halo profile, which is a natural result of $\Lambda $CDM cosmological structure formation, is not an appropriate profile for the dark halos of galaxies, or, new dark matter physics or alternative gravity models are needed to explain the rotational velocities of disk galaxies.

This is the second of a series of papers in which we present a new solution to reconcile the prediction of single stellar population (SSP) models with the observed stellar mass-to-light ($M/L$) ratios of globular clusters (GCs) in M31 and its trend with respect to $\mathrm{[Fe/H]}$. In the present work our focus is on the empirical relation between age and metallicity for GCs and its effect on the $M/L$ ratio. Assuming that there is an anti-correlation between the age of M31 GCs and their metallicity, we evolve dynamical SSP models of GCs to establish a relation between the $M/L$ ratio (in the $V$ and $K$ band) and metallicity. We then demonstrate that the established $M/L-\mathrm{[Fe/H]}$ relation is in perfect agreement with that of M31 GCs. In our models we consider both the canonical initial mass function (IMF) and the top-heavy IMF depending on cluster birth density and metallicity as derived independently from Galactic GCs and ultra-compact dwarf galaxies by Marks et al. Our results signify that the combination of the density- and metallicity-dependent top-heavy IMF, the anti-correlation between age and metallicity, stellar evolution and standard dynamical evolution yields the best possible agreement with the observed trend of $M/L-\mathrm{[Fe/H]}$ for M31 GCs.

Palomar 4 is a low-density globular cluster with a current mass $\approx30000 M_{\odot}$ in the outer halo of the Milky Way with a two-body relaxation time of the order of a Hubble time. Yet, it is strongly mass segregated and contains a stellar mass function depleted of low-mass stars. Pal 4 was either born this way or it is a result of extraordinary dynamical evolution. Since two-body relaxation cannot explain these signatures alone, enhanced mass loss through tidal shocking may have had a strong influence on Pal 4. Here, we compute a grid of direct N-body simulations to model Pal 4 on various eccentric orbits within the Milky Way potential to find likely initial conditions that reproduce its observed mass, half-light radius, stellar MF-slope and line-of-sight velocity dispersion. We find that Pal 4 is most likely orbiting on an eccentric orbit with an eccentricity of $e\approx 0.9$ and pericentric distance of $R_p\approx5$ kpc. In this scenario, the required 3D half-mass radius at birth is similar to the average sizes of typical GCs ($R_h\approx4-5$ pc), while its birth mass is about $M_0\approx10^5 M_{\odot}$. We also find a high degree of primordial mass segregation among the cluster stars, which seems to be necessary in every scenario we considered. Thus, using the tidal effect to constrain the perigalactic distance of the orbit of Pal 4, we predict that the proper motion of Pal 4 should be in the range $-0.52\leq\mu_\delta\leq-0.38$ mas\,yr$^{-1}$ and $-0.30\leq\mu_{\alpha\cos\delta}\leq-0.15$ mas\,yr$^{-1}$.

The observed mass-to-light ($M/L$) ratios of a large sample of GCs in M31 show an inverse trend with metallicity compared to what is expected from Simple Stellar Population (SSP) models with an invariant canonical stellar IMF, in the sense that the observed $M/L$ ratios decrease with increasing metallicity. We show that incorporating the effect of dynamical evolution the SSP models with a canonical IMF can not explain the decreasing $M/L$ ratios with increasing metallicity for the M31 GCs. The recently derived top-heavy IMF as a function of metallicity and embedded cluster density is proposed to explain the lower than expected $M/L$ ratios of metal-rich GCs. We find that the SSP models with a top-heavy IMF, retaining a metallicity- and cluster mass- dependent fraction of the remnants within the clusters, and taking standard dynamical evolution into account can successfully explain the observed $M/L-[Fe/H]$ relation of M31 GCs. Thus we propose that the kinematical data of GCs can be used to constrain the top-heaviness of the IMF in GCs.

Unlike Newtonian dynamics which is linear and obeys the strong equivalence principle, in any nonlinear gravitation such as Milgromian dynamics (MOND), the strong version of the equivalence principle is violated and the gravitational dynamics of a system is influenced by the external gravitational field in which it is embedded. This so called External Field Effect (EFE) is one of the important implications of MOND and provides a special context to test Milgromian dynamics. Here, we study the rotation curves (RCs) of 18 spiral galaxies and find that their shapes constrain the EFE. We show that the EFE can successfully remedy the overestimation of rotation velocities in 80\% of the sample galaxies in Milgromian dynamics fits by decreasing the velocity in the outer part of the RCs. We compare the implied external field with the gravitational field for non-negligible nearby sources of each individual galaxy and find that in many cases it is compatible with the EFE within the uncertainties. We therefore argue that in the framework of Milgromian dynamics, one can constrain the gravitational field induced from the environment of galaxies using their RCs. We finally show that taking into account the EFE yields more realistic values for the stellar mass-to-light ratio in terms of stellar population synthesis than the ones implied without the EFE.

We have obtained membership probabilities of stars within a field of radius $\sim3^\circ$ around the centre of the open cluster Alpha Persei using proper motions and photometry from the PPMXL and WISE catalogues. We have identified 810 possible stellar members of Alpha Persei. We derived the global and radial present-day mass function (MF) of the cluster and found that they are well matched by two-stage power-law relations with different slopes at different radii. The global MF of Alpha Persei shows a turnover at $m=0.62\,\mathrm{M}_{\odot}$ with low and high-mass slopes of $\alpha_\mathrm{low}=0.50\pm0.09$ ($0.1<m/\mathrm{M}_{\odot}<0.62$) and $\alpha_\mathrm{high}=2.32\pm0.14$ ($0.62\leq m/\mathrm{M}_{\odot}<4.68$) respectively. The high-mass slope of the cluster increases from $2.01$ inside $1\hbox{$.\!\!^\circ$}10$ to $2.63$ outside $2\hbox{$.\!\!^\circ$}2$, whereas the mean stellar mass decreases from $0.95$ to $0.57\,\mathrm{M}_{\odot}$ in the same regions, signifying clear evidence of mass segregation in the cluster. From an examination of the high-quality colour-magnitude data of the cluster and performing a series of Monte Carlo simulations we obtained a binary fraction of $f_{\rm bin}=34\pm12$ percent for stars with $0.70<m/\mathrm{M}_{\odot}<4.68$. This is significantly larger than the observed binary fraction, indicating that this open cluster contains a large population of unresolved binaries. Finally, we corrected the mass-function slopes for the effect of unresolved binaries and found low- and high-mass slopes of $\alpha_\mathrm{low}=0.89\pm0.11$ and $\alpha_\mathrm{high}=2.37\pm0.09$ and a total cluster mass of $352\,\mathrm{M}_{\odot}$.

Outer-halo globular clusters show large half-light radii and flat stellar mass functions, depleted in low-mass stars. Using N-body simulations of globular clusters on eccentric orbits within a Milky Way-like potential, we show how a cluster's half-mass radius and its mass function develop over time. The slope of the central mass function flattens proportionally to the amount of mass a cluster has lost, and the half-mass radius grows to a size proportional to the average strength of the tidal field. The main driver of these processes is mass segregation of dark remnants. We conclude that the extended, depleted clusters observed in the Milky Way must have had small half-mass radii in the past, and that they expanded due to the weak tidal field they spend most of their lifetime in. Moreover, their mass functions must have been steeper in the past but flattened significantly as a cause of mass segregation and tidal mass loss.

We present $N$-body simulations of globular clusters, exploring the effect of different galactic potentials on cluster sizes, $r_h$. For various galactocentric distances, $R_G$, we assess how cluster sizes change when we vary the virial mass and concentration of the host galaxy's dark-matter halo. We show that sizes of GCs are determined by the local galactic mass density rather than the virial mass of the host galaxy. We find that clusters evolving in the inner haloes of less concentrated galaxies are significantly more extended than those evolving in more concentrated ones, while the sizes of those orbiting in the outer halo are almost independent of concentration. Adding a baryonic component to our galaxy models does not change these results much, since its effect is only significant in the very inner halo. Our simulations suggest that there is a relation between $r_h$ and $R_G$, which systematically depends on the physical parameters of the halo. Hence, observing such relations in individual galaxies can put a new observational constraint on dark-matter halo characteristics. However, by varying the halo mass in a wide range of $10^{9}\leq M_{vir}/\msun \leq10^{13}$, we find that the $r_h-R_G$ relationship will be nearly independent of halo mass, if one assumes $M_{vir}$ and $c_{vir}$ as two correlated parameters, as is suggested by cosmological simulations.

We perform a series of direct $N$-body calculations to investigate the effect of residual gas expulsion from the gas-embedded progenitors of present-day globular clusters (GCs) on the stellar mass function (MF). Our models start either tidally filling or underfilling, and either with or without primordial mass segregation. We cover 100 Myr of the evolution of modeled clusters and show that the expulsion of residual gas from initially mass-segregated clusters leads to a significantly shallower slope of the stellar MF in the low- ($m\leq 0.50 M_\odot$) and intermediate-mass ($\simeq 0.50-0.85 M_\odot$) regime. Therefore, the imprint of residual gas expulsion and primordial mass segregation might be visible in the present-day MF. We find that the strength of the external tidal field, as an essential parameter, influences the degree of flattening, such that a primordially mass-segregated tidally-filling cluster with $r_h/r_t\geq 0.1$ shows a strongly depleted MF in the intermediate stellar mass range. Therefore, the shape of the present-day stellar MF in this mass range probes the birth place of clusters in the Galactic environment. We furthermore find that this flattening agrees with the observed correlation between the concentration of a cluster and its MF slope, as found by de Marchi et al.. We show that if the expansion through the residual gas expulsion in primordial mass segregated clusters is the reason for this correlation then GCs most probably formed in strongly fluctuating local tidal fields in the early proto-Milky Way potential, supporting the recent conclusion by Marks \& Kroupa.

In order to understand the orbital history of Galactic halo objects, such as globular clusters, authors usually assume a static potential for our Galaxy with parameters that appear at the present-day. According to the standard paradigm of galaxy formation, galaxies grow through a continuous accretion of fresh gas and a hierarchical merging with smaller galaxies from high redshift to the present day. This implies that the mass and size of disc, bulge, and halo change with time. We investigate the effect of assuming a live Galactic potential on the orbital history of halo objects and its consequences on their internal evolution. We numerically integrate backwards the equations of motion of different test objects located in different Galactocentric distances in both static and time-dependent Galactic potentials in order to see if it is possible to discriminate between them. We show that in a live potential, the birth of the objects, 13 Gyr ago, would have occurred at significantly larger Galactocentric distances, compared to the objects orbiting in a static potential. Based on the direct N-body calculations of star clusters carried out with collisional N-body code, NBODY6, we also discuss the consequences of the time-dependence of a Galactic potential on the early- and long-term evolution of star clusters in a simple way, by comparing the evolution of two star clusters embedded in galactic models, which represent the galaxy at present and 12 Gyr ago, respectively. We show that assuming a static potential over a Hubble time for our Galaxy as it is often done, leads to an enhancement of mass-loss, an overestimation of the dissolution rates of globular clusters, an underestimation of the final size of star clusters, and a shallower stellar mass function.

We use direct $N$-body calculations to investigate the impact of primordial mass segregation on the size scale and mass-loss rate of star clusters in a galactic tidal field. We run a set of simulations of clusters with varying degrees of primordial mass segregation at various galactocentric radii and show that, in primordially segregated clusters, the early, impulsive mass-loss from stellar evolution of the most massive stars in the innermost regions of the cluster leads to a stronger expansion than for initially non-segregated clusters. Therefore, models in stronger tidal fields dissolve faster due to an enhanced flux of stars over the tidal boundary. Throughout their lifetimes, the segregated clusters are more extended by a factor of about 2, suggesting that (at least) some of the very extended globular clusters in the outer halo of the Milky Way may have been born with primordial mass segregation. We finally derive a relation between star-cluster dissolution time, $T_{diss}$, and galactocentric radius, $R_G$, and show how it depends on the degree of primordial mass segregation.

We use direct $N$-body calculations to study the evolution of the unusually extended outer halo globular cluster Palomar 4 (Pal~4) over its entire lifetime in order to reproduce its observed mass, half-light radius, velocity dispersion and mass function slope at different radii. We find that models evolving on circular orbits, and starting from a non-mass segregated, canonical initial mass function (IMF) can reproduce neither Pal 4's overall mass function slope nor the observed amount of mass segregation. Including either primordial mass segregation or initially flattened IMFs does not reproduce the observed amount of mass segregation and mass function flattening simultaneously. Unresolved binaries cannot reconcile this discrepancy either. We find that only models with both a flattened IMF and primordial segregation are able to fit the observations. The initial (i.e. after gas expulsion) mass and half-mass radius of Pal~4 in this case are about 57000 M${\odot}$ and 10 pc, respectively. This configuration is more extended than most globular clusters we observe, showing that the conditions under which Pal~4 formed must have been significantly different from that of the majority of globular clusters. We discuss possible scenarios for such an unusual configuration of Pal~4 in its early years.

We present the first ever direct $N$-body computations of an old Milky Way globular cluster over its entire life time on a star-by-star basis. Using recent GPU hardware at Bonn University, we have performed a comprehensive set of $N$-body calculations to model the distant outer halo globular cluster Palomar 14 (Pal 14). By varying the initial conditions we aim at finding an initial $N$-body model which reproduces the observational data best in terms of its basic parameters, i.e. half-light radius, mass and velocity dispersion. We furthermore focus on reproducing the stellar mass function slope of Pal 14 which was found to be significantly shallower than in most globular clusters. While some of our models can reproduce Pal 14's basic parameters reasonably well, we find that dynamical mass segregation alone cannot explain the mass function slope of Pal 14 when starting from the canonical Kroupa initial mass function (IMF). In order to seek for an explanation for this discrepancy, we compute additional initial models with varying degrees of primordial mass segregation as well as with a flattened IMF. The necessary degree of primordial mass segregation turns out to be very high. This modelling has shown that the initial conditions of Pal 14 after gas expulsion must have been a half-mass radius of about 20 pc, a mass of about 50000 M$_{\odot}$, and possibly some mass segregation or an already established non-canonical IMF depleted in low-mass stars. Such conditions might be obtained by a violent early gas-expulsion phase from an embedded cluster born with mass segregation. Only at large Galactocentric radii are clusters likely to survive as bound entities the destructive gas-expulsion process we seem to have uncovered for Pal 14. In addition we compute a model with a 5% primordial binary fraction to test if such a population has an effect on the cluster's evolution.

Alternative gravitations of Milgrom (MOND), Moffat (MOG), and CDM scenarios all simulate rotation curves of spirals with reasonable details. They display significant disparities however in predicting the stellar mass-to-light ($M_*/L$) ratios of the galaxies. We maintain this feature could serve as a distinguishing factor between different alternative theories. We analyze the rotation curves of 46 low- and high-surface brightness galaxies and compare the resulting $M_*/L$s with the predictions of the Stellar Population Synthesis (SPS) scheme. The color-$M_*/L$ correlation obtained for MOND is consistent with predictions of SPS models. MOG does not show this consistency, and the $M_*/L$s of CDM model shows large dispersions. Furthermore, $M_*/L$ ratios of MOND with Bekenstein interpolating function favor Kroupa's initial mass function (IMF) of the SPS scheme, while those of MOND with standard and simple interpolating functions are consistent with Salpeter's IMF. Here is another indication to differentiate between different IMFs that are used in SPS context.

Rotation curves of spiral galaxies \emph{i}) fall off much less steeply than the Keplerian curves do, and \emph{ii}) have asymptotic speeds almost proportional to the fourth root of the mass of the galaxy, the Tully-Fisher relation. These features alone are sufficient for assigning a dark companion to the galaxy in an unambiguous way. In regions outside a spherical system, we design a spherically symmetric spacetime to accommodate the peculiarities just mentioned. Gravitation emerges in excess of what the observable matter can produce. We attribute the excess gravitation to a hypothetical, dark, perfect fluid companion to the galaxy and resort to the Tully-Fisher relation to deduce its density and pressure. The dark density turns out to be proportional to the square root of the mass of the galaxy and to fall off as $r^{-(2+\alpha)}, \alpha\ll 1$. The dark equation of state is barrotropic. For the interior of the configuration, we require the continuity of the total force field at the boundary of the system. This enables us to determine the size and the distribution of the interior dark density and pressure in terms of the structure of the observable matter. The formalism is nonlocal and nonlinear, and the density and pressure of the dark matter at any spacetime point turn out to depend on certain integrals of the baryonic matter over all or parts of the system in a nonlinear manner.

We model the dynamics of Magellanic Stream with the ram-pressure scenario in the logarithmic and power-law galactic halo models and construct numerically the past orbital history of Magellanic clouds and Magellanic Stream. The parameters of models include the asymptotic rotation velocity of spiral arms, halo flattening, core radius and rising or falling parameter of rotation curve. We obtain the best fit parameters of galactic models through the maximum likelihood analysis, comparing the high resolution radial velocity data of HI in Magellanic Stream with that of theoretical models. The initial condition of the Magellanic Clouds is taken from the six different values reported in the literature. We find that oblate and nearly spherical shape halos provide a better fit to the observation than the prolate halos. This conclusion is almost independent of choosing the initial conditions and is valid for both logarithmic and power-law models.