Using two-dimensional simulations, we compute the torque and rate of work (power) on a low-mass gravitational body, with softening length $R_{\rm soft}$, embedded in a gaseous disk when its orbit is eccentric and retrograde with respect to the disk. We explore orbital eccentricities $e$ between $0$ and $0.6$. We find that the power has its maximum at $e\simeq 0.25(h/0.05)^{2/3}$, where $h$ is the aspect ratio of the disk. We show that the power and the torque converge to the values predicted in the local (non-resonant) approximation of the dynamical friction (DF) when $R_{\rm soft}$ tends to zero. For retrograde inspirals with mass ratios $\lesssim 5\times 10^{-4}$ embedded in disks with $h\geq 0.025$, our simulations suggest that (i) the rate of inspiral barely depends on the orbital eccentricity and (ii) the local approximation provides the value of this inspiral rate within a factor of $1.5$. The implications of the results for the orbital evolution of extreme mass-ratio inspirals are discussed.

We study the dynamical evolution of Jupiter and Saturn embedded in a gaseous, solar-nebula-type disc by means of hydrodynamics simulations with the FARGO2D1D code. We study the evolution for different initial separations of the planets' orbits, Delta a_SJ , to investigate whether they become captured in mean motion resonance (MMR) and the direction of the subsequent migration of the planet (inward or outward). We also provide an assessment of the planet's orbital dynamics at different epochs of Saturn's growth. We find that the evolution of initially compact orbital configurations is dependent on the value of Delta a_SJ . This implies that an evolution as that proposed in the Grand Tack model depends on the precise initial orbits of Jupiter and Saturn and on the timescales for their formation. Capture in the 1:2 MMR and inward or (nearly) stalled migration are highly favoured. Within its limits, our work suggests that the reversed migration, associated with the resonance capture of Jupiter and Saturn, may be a low probability evolutionary scenario, so that other planetary systems with giant planets are not expected to have experienced a Grand Tack-like evolutionary path.

We study the tidal interaction between a low-mass companion (e.g., a protoplanet or a black hole) in orbit about a central mass, and the accretion disk within which it is submerged. We present results for a companion on a coplanar orbit with eccentricity e between 0.1 and 0.6. For these eccentricities, dynamical friction arguments in its local approximation, that is, ignoring differential rotation and the curvature of the orbit, provide simple analytical expressions for the rates of energy and angular momentum exchange between the disk and the companion. We examine the range of validity of the dynamical friction approach by conducting a series of hydrodynamical simulations of a perturber with softening radius R_soft embedded in a two-dimensional disk. We find close agreement between predictions and the values in simulations provided that R_soft is chosen sufficiently small, below a threshold value Rtilde_soft, which depends on the disk parameters and on eccentricity. We give Rtilde_soft for both razor-thin disks and disks with a finite scaleheight. For point-like perturbers, the local approximation is valid if the accretion radius is smaller than Rtilde_soft. This condition imposes an upper value on the mass of the perturber.

An extended gravitational object embedded in a parent system comprised of gas and collisionless particles may undergo both dynamical friction (DF) and mass loss by tidal forces. If the object is compact enough, it can increase its mass through accretion of material from the surrounding medium. We extend the classical linear analysis of DF on a constant-mass body in a gaseous medium to the case where its mass changes with time. We show that the structure of the wake may differ significantly from the constant-mass case. For instance, the front-back symmetry of density about subsonic constant-mass perturbers is broken down for variable-mass perturbers. The density wake keeps a memory of the past mass history of the perturber. For dissolving perturbers, the density field is more dense than expected using the instantaneous mass of the perturber in the classical formula. As a consequence, the instantaneous-mass approximation underestimates the drag force for mass-losing perturbers and overestimates it for mass-gaining perturbers. We present cases in which the percentage error in the drag force using the instantaneous-mass approximation is greater than 50%.

We evaluate the torque acting on a gravitational perturber on a retrograde circular orbit in the midplane of a gaseous disk. We assume that the mass of this satellite is so low it weakly disturbs the disk (type I migration). The perturber may represent the companion of a binary system with a small mass ratio. We compare the results of hydrodynamical simulations with analytic predictions. Our two-dimensional (2D) simulations indicate that the torque acting on a perturber with softening radius $R_{\rm soft}$ can be accounted for by a scattering approach if $R_{\rm soft}<0.3H$, where $H$ is defined as the ratio between the sound speed and the angular velocity at the orbital radius of the perturber. For $R_{\rm soft}>0.3H$, the torque may present large and persistent oscillations, but the resultant time-averaged torque decreases rapidly with increasing $R_{\rm soft}/H$, in agreement with previous analytical studies. We then focus on the torque acting on small-size perturbers embedded in full three-dimensional (3D) disks and argue that the density waves propagating at distances $\lesssim H$ from the perturber contribute significantly to the torque because they transport angular momentum. We find a good agreement between the torque found in 3D simulations and analytical estimates based on ballistic orbits. We compare the radial migration timescales of prograde versus retrograde perturbers. For a certain range of the perturber's mass and aspect ratio of the disk, the radial migration timescale in the retrograde case may be appreciably shorter than in the prograde case. We also provide the smoothing length required in 2D simulations in order to account for 3D effects.

We study the relativistic Bondi-Hoyle-Lyttleton accretion onto a Schwarzschild black hole (BH), which is surrounded by rigid and small, randomly distributed, bodies. These bodies are idealized representations of substructure like stars passing close to the BH, bubbles created by stellar winds or cold clumps.We explore cases where the filling factor of these bodies is small. The flow is assumed to be adiabatic and move supersonically towards the black hole. The interaction with these rigid obstacles transforms ram pressure of the flow into thermal pressure through bow shocks, slowing down the flow and making the accreting gas turbulent. As a consequence, although the flow reaches a statistically-steady state, the accretion rate presents some variability. For a flow Mach number at infinity of 4, a few of these objects (5 - 10) are enough to increase the accretion rate about 50% over the accretion rate without bodies, even though the gas is adiabatic and the filling factor of the obstacles is small.

We study gap formation in gaseous protoplanetary discs by a Jupiter mass planet. The planet's orbit is circular and inclined relative to the midplane of the disc. We use the impulse approximation to estimate the gravitational tidal torque between the planet and the disc, and infer the gap profile. For low-mass discs, we provide a criterion for gap opening when the orbital inclination is $\leq 30^{\circ}$. Using the FARGO3D code, we simulate the disc response to an inclined massive planet. The dependence of the depth and width of the gap obtained in the simulations on the inclination of the planet is broadly consistent with the scaling laws derived in the impulse approximation. Although we mainly focus on planets kept on fixed orbits, the formalism permits to infer the temporal evolution of the gap profile in cases where the inclination of the planet changes with time. This study may be useful to understand the migration of massive planets on inclined orbit, because the strength of the interaction with the disc depends on whether a gap is opened or not.

We analyse the rotation curves and gravitational stability of a sample of six bulgeless galaxies for which detailed images reveal no evidence for strong bars. We explore two scenarios: Newtonian dark matter models and MOdified Newtonian Dynamics (MOND). By adjusting the stellar mass-to-light ratio, dark matter models can match simultaneously both the rotation curve and bar-stability requirements in these galaxies. To be consistent with stability constraints, in two of these galaxies, the stellar mass-to-light ratio is a factor of ~1.5-2 lower than the values suggested from galaxy colours. In contrast, MOND fits to the rotation curves are poor in three galaxies, perhaps because the gas tracer contains noncircular motions. The bar stability analysis provides a new observational test to MOND. We find that most of the galaxies under study require abnormally-high levels of random stellar motions to be bar stable in MOND. In particular, for the only galaxy in the sample for which the line-of-sight stellar velocity dispersion has been measured (NGC 6503), the observed velocity dispersion is not consistent with MOND predictions because it is far below the required value to guarantee bar stability. Precise measurements of mass-weighted velocity dispersions in (unbarred and bulgeless) spiral galaxies are crucial to test the consistency of MOND.

Dynamical estimates of the mass surface density at the solar radius can be made up to a height of 4 kpc using thick disk stars as tracers of the potential. We investigate why different Jeans estimators of the local surface density lead to puzzling and conflicting results. Using the Jeans equations, we compute the vertical (F_z) and radial (F_R) components of the gravitational force, as well as Gamma(z), defined as the radial derivative of V_c^2, with V_c^{2}= -RF_R. If we assume that the thick disk does not flare and that all the components of the velocity dispersion tensor of the thick disk have a uniform radial scalelength of 3.5 kpc, Gamma takes implausibly large negative values, when using the currently available kinematical data of the thick disk. This implies that the input parameters or the model assumptions must be revised. We have explored, using a simulated thick disk, the impact of the assumption that the scale lengths of the density and velocity dispersions do not depend on the vertical height z above the midplane. In the lack of any information about how these scale radii depend on z, we define a different strategy. By using a parameterized Galactic potential, we find that acceptable fits to F_z, F_R and Gamma are obtained for a flaring thick disk and a spherical dark matter halo with a local density larger than 0.0064 M_sun pc^{-3}. Disk-like dark matter distributions might be also compatible with the current data of the thick disk. A precise measurement of Gamma at the midplane could be very useful to discriminate between models.

The cold dark matter (CDM) model has two unsolved issues: simulations overpredict the satellite abundance around the Milky Way (MW) and it disagrees with observations of the central densities of dwarf galaxies which prefer constant density (core) profiles.One alternative explanation known as the scalar field dark matter (SFDM) model, assumes that the dark matter is a scalar field of mass($\sim 10^{-22}$ eV/$c^2$); this model can reduce the overabundance issue due to the lack of halo formation below a mass scale of $\sim 10^8$M$_{\odot}$ and successfully fits the density distribution in dwarfs. One of the attractive features of the model is predicting core profiles in halos, although the determination of the core sizes is set by fitting the observational data. We perform \textit{N}-body simulations to explore the influence of tidal forces over a stellar distribution embedded in a SFDM halo orbiting a MW-like SFDM host halo with a disk. Our simulations intend to test the viability of SFDM as an alternative model by comparing the tidal effects that result in this paradigm with those obtained in CDM for similar mass halos. We found that galaxies in subhalos with core profiles and high central densities survive for 10 Gyr. The same occurs for galaxies in low density subhalos located far from the host disk influence, whereas satellites in low density DM halos and in tight orbits can eventually be stripped of stars. We conclude that SFDM shows consistency with results from CDM for dwarf galaxies, but naturally offer a possibility to solve the missing satellite problem.

Using time-dependent linear theory, we investigate the morphology of the gravitational wake induced by a binary, whose center of mass moves at velocity Vcm against a uniform background of gas. For simplicity, we assume that the binary's components are on circular orbits about their common center of mass. The consequences of dynamical friction is twofold. First, gas dynamical friction may drag the binary's center of mass and cause the binary to migrate. Second, drag forces also induce a braking torque, which causes the orbits of the binary components to shrink. We compute the drag forces acting on one component of the binary due to the gravitational interaction with its own wake. We show that the dynamical friction force responsible to decelerate the binary's center of mass is smaller than it is in the point-mass case because of the loss of gravitational focusing. We show that the braking internal torque depends on the Mach numbers of each binary component about their center of mass, and also on the Mach number of the center of mass of the binary. In general, the internal torque decreases with increasing the velocity of the binary relative to the ambient gas cloud. However, this is not always the case. We also mention the relevance of our results on the period distribution of binaries.

We provide constraints on the inclination angle of the H\,{\sc i} disk of the dwarf irregular galaxy Holmberg II (Ho II) from stability analysis of the outer gaseous disk. We point out that a mean inclination angle of 27 degrees and thus a flat circular velocity of ~60 km/s, is required to have a level of gravitational stability similar to that found in other galaxies. Adopting this inclination angle, we find that Ho II lies on the right location in the baryonic Tully-Fisher relation. Moreover, for this inclination, its rotation curve is consistent with MOND. However, the corresponding analysis of the stability under MOND indicates that this galaxy could be problematic for MOND because its outer parts are marginally unstable in this gravity theory. We urge MOND simulators to study numerically the non-linear stability of gas-rich dwarf galaxies since it may provide a new key test for MOND.

The cold dark matter model predicts cuspy dark matter halos. However, it has been found that, in some low-mass galaxies, cored dark halos provide a better description of their internal dynamics. Here we give constraints on the dark halo profile in the Sextans dwarf spheroidal galaxy by studying the longevity of two cold kinematic substructures detected in this galaxy. We perform N-body simulations of a stellar clump in the Sextans dwarf galaxy, including a live dark matter halo and the main stellar component. We find that, if the dark halo is cuspy, stellar clumps orbiting with semi-major axis ~400 pc are disrupted in ~5 Gyr, even if the clump is initially as compact stellar cluster with a radius of r_c=5 pc. Stellar clusters in an initial orbit with semi-major axis \leq 250 pc may survive to dissolution but their orbits decay towards the center by dynamical friction. In contrast, the stellar clumps can persist for a Hubble time within a cored dark matter halo, even if the initial clump's radius is as extended as r_c=80 pc. We also study the evolution of the clump in the MONDian context. In this scenario, we find that even an extended stellar clump with radius r_c=80 pc survives for a Hubble time, but an unrealistic value for the stellar mass-to-light ratio of 9.2 is needed.

Using axisymmetrical numerical simulations, we revisit the gravitational drag felt by a gravitational Plummer sphere with mass M and core radius Rs, moving at constant velocity V0 through a background homogeneous medium of adiabatic gas. Since the potential is non-diverging, there is no gas removal due to accretion. When Rs is larger than the Bondi radius RB, the perturbation is linear at every point and the drag force is well fitted by the time-dependent Ostriker's formula with r_{min}= 2.25Rs, where r_{min} is the minimum impact parameter in the Coulomb logarithm. In the deep nonlinear supersonic regime (Rs<< RB), the minimum radius is no longer related with Rs but with RB. We find r_min=3.3mach^{-2.5}RB, for Mach numbers of the perturber between $1.5$ and $4$, although r_{min} = 2\mach^{-2}RB=2GM/V0^{2} also provides a good fit at mach>2. As a consequence, the drag force does not depend sensitively on the nonlinearity parameter RB/Rs, for RB/Rs-values larger than a certain critical value. We show that our generalized Ostriker's formula for the drag force is more accurate than the formula suggested by Kim & Kim (2009).

We quantify the effects of magnetic fields, cosmic rays and gas pressure on the rotational velocity of HI gas in the Milky Way, at galactic distances between Rsun and 2Rsun. The magnetic field is modelled by two components; a mainly azimuthal magnetic component and a small-scale tangled field. We construct a range of plausible axisymmetric models consistent with the strength of the total magnetic field as inferred from radio synchrotron data. In a realistic Galactic disk, the pressure by turbulent motions, cosmic rays and the tangled turbulent field provide radial support to the disk. Large-scale (ordered) magnetic fields may or may not provide support to the disk, depending on the local radial gradient of the azimuthal field. We show that for observationally constrained models, magnetic forces cannot appreciably alter the tangential velocity of HI gas within a galactic distance of 2Rsun.

We have carried out a search for gas-rich dwarf galaxies that have lower rotation velocities in their outskirts than MOdified Newtonian Dynamics (MOND) predicts, so that the amplitude of their rotation curves cannot be fitted by arbitrarily increasing the mass-to-light ratio of the stellar component or by assuming additional undetected matter. With presently available data, the gas-rich galaxies UGC 4173, Holmberg II, ESO 245-G05, NGC 4861 and ESO 364-G029 deviate most from MOND predictions and, thereby, provide a sample of promising targets in testing the MOND framework. In the case of Holmberg II and NGC 4861, we find that their rotation curves are probably inconsistent with MOND, unless their inclinations and distances differ significantly from the nominal ones. The galaxy ESO 364-G029 is a promising target because its baryonic mass and rotation curve are similar to Holmberg II but presents a higher inclination. Deeper photometric and HI observations of ESO 364-G029, together with further decreasing systematic uncertainties, may provide a strong test to MOND.

Through numerical simulations, we study the dissolution timescale of the Ursa Minor cold stellar clump, due to the combination of phase-mixing and gravitational encounters with compact dark substructures in the halo of Ursa Minor. We compare two scenarios; one where the dark halo is made up by a smooth mass distribution of light particles and one where the halo contains 10% of its mass in the form of substructures (subhalos). In a smooth halo, the stellar clump survives for a Hubble time provided that the dark matter halo has a big core. In contrast, when the point-mass dark substructures are added, the clump survives barely for \sim 1.5 Gyr. These results suggest a strong test to the \Lambda-cold dark matter scenario at dwarf galaxy scale.

We consider the hypothesis that galactic dark matter is composed of ultra-light scalar particles and use internal properties of dwarf spheroidal galaxies to establish a preferred range for the mass m of these bosonic particles. We re-investigate the problem of the longevity of the cold clump in Ursa Minor and the problem of the rapid orbital decay of the globular clusters in Fornax and dwarf ellipticals. Treating the scalar field halo as a rigid background gravitational potential and using N-body simulations, we have explored how the dissolution timescale of the cold clump in Ursa Minor depends on m. It is demonstrated that for masses in the range 0.3x10^-22 eV < m <1x10^-22 eV, scalar field dark halos without self-interaction would have cores large enough to explain the longevity of the cold clump in Ursa Minor and the wide distribution of globular clusters in Fornax, but small enough to make the mass of the dark halos compatible with dynamical limits. It is encouraging to see that this interval of m is consistent with that needed to suppress the overproduction of substructure in galactic halos and is compatible with the acoustic peaks of cosmic microwave radiation. On the other hand, for self-interacting scalar fields with coupling constant l, values of m^4/l <= 0.55x10^3 eV^4 are required to account for the properties of the dark halos of these dwarf spheroidal galaxies.

The dynamical friction force experienced by a massive gravitating body moving through a gaseous medium is modified by sufficiently strong large-scale magnetic fields. Using linear perturbation theory, we calculate the structure of the wake generated by, and the gravitational drag force on, a body traveling in a straight-line trajectory in a uniformly magnetized medium. The functional form of the drag force as a function of the Mach number (V_0/c_s, where V_0 is the velocity of the body and c_s the sound speed) depends on the strength of the magnetic field and on the angle between the velocity of the perturber and the direction of the magnetic field. In particular, the peak value of the drag force is not near Mach number 1 for a perturber moving in a sufficiently magnetized medium. As a rule of thumb, we may state that for supersonic motion, magnetic fields act to suppress dynamical friction; for subsonic motion, magnetic fields tend to enhance dynamical friction. For perturbers moving along the magnetic field lines, the drag force at some subsonic Mach numbers may be stronger than it is at supersonic velocities. We also mention the relevance of our findings to black hole coalescence in galactic nuclei.

We explore a ballistic orbit model to infer the gravitational drag force on an accreting point mass M, such as a black hole, moving at a hypersonic velocity v_{0} through a gaseous environment of density \rho_{0}. The streamlines blend in the flow past the body and transfer momentum to it. The total drag force acting on the body, including the nonlinear contribution of those streamlines with small impact parameter that bend significantly and pass through a shock, can be calculated by imposing conservation of momentum. In this fully analytic approach, the ambiguity in the definition of the lower cut-off distance $r_{\rm min}$ in calculations of the effect of dynamical friction is removed. It turns out that $r_{\rm min}=\sqrt{e}GM/2v_{0}^{2}$. Using spherical surfaces of control of different sizes, we carry out a successful comparison between the predicted drag force and the one obtained from a high resolution, axisymmetric, isothermal flow simulation. We demonstrate that ballistic models are reasonably successful in accounting for both the accretion rate and the gravitational drag.

We compare the basic properties and kinematics of two gas-rich dwarf galaxies: KK246 and Holmberg II (HoII). HoII is 20 times more luminous in the blue-band than KK246 and its HI mass is a factor of 6 higher than in KK246. However, the amplitudes of the rotation curves (at the last measured point) of both galaxies are very similar, of about 40 km/s at a galactocentric radius of 7 kpc. This fact is challenging for modified theories of gravity that predict a one-to-one relation between gravity at any radius and the enclosed baryonic mass in galaxies. In particular, MOdified Newtonian Dynamics (MOND) predicts an asymptotic flat velocity of 60 km/s in HoII. Since the baryonic mass of HoII is dominated by the gas component, MOND overboosts its rotation speed even if the mass of the stellar disk is taken negligibly small. We conclude the rotation curve of HoII is probably inconsistent with MOND, unless the inclination and distance are both fine-tuned to very unlikely values.

Estimates of the dynamical surface mass density at the solar Galactocentric distance are commonly derived assuming that the disk is in vertical equilibrium with the Galactic potential. This assumption has recently been called into question, based on the claim that the ratio between the kinetic and the gravitational energy in such solutions is a factor of 3 larger than required if Virial equilibrium is to hold. Here we show that this ratio between energies was overerestimated and that the disk solutions are likely to be in Virial equilibrium after all. We additionally demonstrate, using one-dimensional numerical simulations, that the disks are indeed in equilibrium. Hence, given the uncertainties, we find no reason to cast doubt on the steady-state solutions which are traditionally used to measure the matter density of the disk.

We study the final equilibrium states of the Parker instability arising from an initially unstable cylindrical equilibrium configuration of gas in the presence of a radial gravitational field and a longitudinal magnetic field. The aim of this work is to compare the properties of the nonlinear final equilibria with those found in a system with Cartesian geometry. Maps of the density and magnetic field lines, when the strength of the gravitational field is constant, are given in both geometries. In the axisymmetric model, the magnetic field tends to expand in radius, forming magnetic arcades, while knots of gas are formed because the plasma drains radially and strangulates the magnetic field lines, leading to the formation of magnetic bottlenecks. We find that the magnetic buoyancy and the drainage of gas along field lines are less efficient under axial symmetry than in a Cartesian atmosphere. As a consequence, the column density enhancement arising in gas condensations in the axially-symmetric model is smaller than in Cartesian geometry. The magnetic-to-gas pressure ratio in the final state takes more extreme values in the Cartesian model. Models with non-uniform radial gravity are also discussed.

The survival of unbound density substructure against orbital mixing imposes strong constraints on the slope of the underlying gravitational potential and provides a new test on modified gravities. Here we investigate whether the interpretation that the stellar clump in Ursa Minor (UMi) dwarf spheroidal galaxy is a `dynamical fossil' is consistent with Modified Newtonian dynamics (MOND). For UMi mass models inferred by fitting the velocity dispersion profile, the stellar clump around the second peak of UMi is erased very rapidly, within 1.25 Gyr (6.5 orbits), even with the inclusion of self-gravity. We find that the clump can hardly survive for more than 2 Gyr even under more generous conditions. Alternative scenarios which could lead to a kinematically cold clump are discussed but, so far, none of them were found to be fully satisfactory. Our conclusion is that the cold clump in UMi poses a challenge for both LambdaCDM and MOND.

From various cosmological, astrophysical and terrestrial requirements, we derive conservative upper bounds on the present-day fraction of the mass of the Galactic dark matter (DM) halo in charged massive particles (CHAMPs). If dark matter particles are neutral but decay lately into CHAMPs, the lack of detection of heavy hydrogen in sea water and the vertical pressure equilibrium in the Galactic disc turn out to put the most stringent bounds. Adopting very conservative assumptions about the recoiling velocity of CHAMPs in the decay and on the decay energy deposited in baryonic gas, we find that the lifetime for decaying neutral DM must be > (0.9-3.4)x 10^3 Gyr. Even assuming the gyroradii of CHAMPs in the Galactic magnetic field are too small for halo CHAMPs to reach Earth, the present-day fraction of the mass of the Galactic halo in CHAMPs should be < (0.4-1.4)x 10^{-2}. We show that redistributing the DM through the coupling between CHAMPs and the ubiquitous magnetic fields cannot be a solution to the cuspy halo problem in dwarf galaxies.

The well-established correlations between the mass of massive black holes (BHs) in the nuclei of most studied galaxies and various global properties of their hosting galaxy lend support to the idea that dwarf galaxies and globular clusters could also host a BH in their centers. Direct kinematic detection of BHs in dwarf spheroidal (dSph) galaxies are seriously hindered by the small number of stars inside the gravitational influence region of the BH. The aim of this Letter is to establish an upper dynamical limit on the mass of the putative BH in the Ursa Minor (UMi) dSph galaxy. We present direct N-body simulations of the tidal disruption of the dynamical fossil observed in UMi, with and without a massive BH. We find that the observed substructure is incompatible with the presence of a massive BH of (2-3)x10^4 Msun within the core of UMi. These limits are consistent with the extrapolation of the M_{BH}-sigma relation to the M_{BH}<10^6 Msun regime. We also show that the BH may be off-center with respect to the center of symmetry of the whole galaxy.

There exist many physical processes that may contribute to the driving of turbulence in galactic disks. Some of them could drive turbulence even in the absence of star formation. For example, hydrodynamic (HD) or magnetohydrodynamic (MHD) instabilities, frequent mergers of small satellite clumps, ram pressure, or infalling gas clouds. In this work we present numerical simulations to study the interaction of compact high velocity clouds (CHVC) with the outskirts of magnetized gaseous disks. With our numerical simulations we show that the rain of small HVCs onto the disk is a potential source of random motions in the outer parts of HI disks.

We know that the galactic magnetic field possesses a random component in addition to the mean uniform component, with comparable strength of the two components. This random component is considered to play important roles in the evolution of the interstellar medium (ISM). In this work we present numerical simulations associated with the interaction of the supersonic flows located at high latitude in our Galaxy (High Velocity Clouds, HVC) with the magnetized galactic ISM in order to study the effect that produces a random magnetic field in the evolution of this objects.

We improve earlier Galactic bounds that can be placed on the fraction of dark matter in charged elemental particles (CHAMPs). These constraints are of interest for CHAMPs whose mass is too large for them to have seen through their electromagnetic interaction with ordinary matter, and whose gyroradius in the galactic magnetic field is too small for halo CHAMPs to reach Earth. If unneutralized CHAMPs in that mass range are well mixed in the halo, they can at most make up a fraction < (3-7) x 10^{-3} of the mass of the Galactic halo. CHAMPs might still be a solution to the cuspy halo problem if they decay to neutral dark matter but a fine-tuning is required. We also discuss the case where CHAMPs do not populate a spherical halo.

We calculate the structure of a wake generated by, and the dynamical friction force on, a gravitational perturber travelling through a gaseous medium of uniform density and constant background acceleration g_ext, in the context of Modified Newtonian Dynamics (MOND). The wake is described as a linear superposition of two terms. The dominant part displays the same structure as the wake generated in Newtonian gravity scaled up by a factor mu^{-1}(g_ext/a_0), where a_{0} is the constant MOND acceleration and mu the interpolating function. The structure of the second term depends greatly on the angle between g_{ext} and and the velocity of the perturber. We evaluate the dynamical drag force numerically and compare our MOND results with the Newtonian case. We mention the relevance of our calculations to orbit evolution of globular clusters and satellites in a gaseous proto-galaxy. Potential differences in the X-ray emission of gravitational galactic wakes in MOND and in Newtonian gravity with a dark halo are highlighted.

In many astrophysical situations, as in the coalescence of supermassive black hole pairs at gas rich galactic nuclei, the dynamical friction experienced by an object is a combination of its own wake as well as the wakes of its companions. Using a semi-analytic approach, we investigate the composite wake due to, and the resulting drag forces on, double perturbers that are placed at the opposite sides of the orbital center and move on a circular orbit in a uniform gaseous medium. The circular orbit makes the wake of each perturber asymmetric, creating an overdense tail at the trailing side. The tail not only drags the perturber backward but it also exerts a positive torque on the companion. For equal-mass perturbers, the positive torque created by the companion wake is, on average, a fraction ~40-50% of the negative torque created by its own wake, but this fraction may be even larger for perturbers moving subsonically. This suggests that the orbital decay of a perturber in a double system, especially in the subsonic regime, can take considerably longer than in isolation. We provide the fitting formulae for the forces due to the companion wake and discuss our results in light of recent numerical simulations for mergers of binary black holes.

The outskirts of galaxies are a very good laboratory for testing the nature of the gravitational field at low accelerations. By assuming that the neutral hydrogen gas is in hydrostatic equilibrium in the gravitational potential of the host galaxy, the observed flaring of the gas layer can be used to test modified gravities. For the first time we construct a simple framework to derive the scaleheight of the neutral hydrogen gas disc in the MOND scenario and apply this to the Milky Way. It is shown that using a constant gas velocity dispersion of ~9 km/s, MOND is able to give a very good fit to the observed HI flaring beyond a galactocentric distance of 17 kpc up to the last measured point (~40 kpc). Between 10 and 16 kpc, however, the observed scaleheight is about 40% more than what MOND predicts for the standard interpolating function and 70% for the form suggested by Famaey & Binney. Given the uncertainties in the non-thermal pressure support by cosmic rays and magnetic fields, MOND seems to be a plausible alternative to dark matter in explaining the Milky Way flaring. Studying the flaring of extended HI discs in external edge-on galaxies may be a promising approach to assess the viability of MOND.

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

High--resolution 2D--MHD numerical simulations have been carried out to investigate the effects of continuing infall of clumpy gas in extended HI galactic disks. Given a certain accretion rate, the response of the disk depends on its surface gas density and temperature. For Galactic conditions at a galactocentric distance of ~20 kpc, and for mass accretion rates consistent with current empirical and theoretical determinations in the Milky Way, the rain of compact high velocity clouds onto the disk can maintain transonic turbulent motions in the warm phase (~2500 K) of HI. Hence, the HI line width is expected to be ~6.5 km/s for a gas layer at 2500 K, if infall were the only mechanism of driving turbulence. Some statistical properties of the resulting forcing flow are shown in this Letter. The radial dependence of the gas velocity dispersion is also discussed.

We establish new dynamical constraints on the mass and abundance of compact objects in the halo of dwarf spheroidal galaxies. In order to preserve kinematically cold the second peak of the Ursa Minor dwarf spheroidal (UMi dSph) against gravitational scattering, we place upper limits on the density of compact objects as a function of their assumed mass. The mass of the dark matter constituents cannot be larger than 1000 solar masses at a halo density in UMi's core of 0.35 solar masses/pc^3. This constraint rules out a scenario in which dark halo cores are formed by two-body relaxation processes. Our bounds on the fraction of dark matter in compact objects with masses >3000 solar masses improve those based on dynamical arguments in the Galactic halo. In particular, objects with masses $\sim 10^{5}$ solar masses can comprise no more than a halo mass fraction $\sim 0.01$. Better determinations of the velocity dispersion of old overdense regions in dSphs may result in more stringent constraints on the mass of halo objects. For illustration, if the preliminary value of 0.5 km/s for the secondary peak of UMi is confirmed, compact objects with masses above $\sim 100$ solar masses could be excluded from comprising all its dark matter halo.

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

In this paper we consider the viability of these non-gravitational torques to take the gas off the plane. We show that magnetically generated warps are clearly flawed because they would wrap up into a spiral in less than two or three galactic rotations. The inclusion of any magnetic diffusivity to dilute the wrapping effect, causes the amplitude of the warp to damp. We also consider observational consequences of the accretion of an intergalactic plane-parallel flow at infinity. We have computed the amplitude and warp asymmetry in the accretion model, for a disc embedded in a flattened dark matter halo, including self-consistently the contribution of the modes with azimuthal wavenumbers m=0 and m=1. This model predicts quite asymmetric warps, maximum gas displacements on the two sides in the ratio 3:2 for the preferred Galactic parameters,and the presence of a fraction ~3.5% of U-shaped warps, at least. The azimuthal dependence of the moment transfer by the ram pressure would produce a strong asymmetry in the thickness of the HI layer and asymmetric density distributions in z, in conflict with observational data for the warp in our Galaxy and in external galaxies. The amount of accretion that is required to explain the Galactic warp would give gas scaleheights in the far outer disc that are too small. We conclude that accretion of a flow with no net angular momentum, cannot be the main and only cause of warps.

Self-interacting dark matter has been proposed as a hypothesis to explain the shallow central slopes of the density profiles of dark matter halos in galaxies. In order to be consistent with observational studies at scales of galaxy clusters, the cross section should scale inversely with the velocity of collision. In this paper we consider the mass density profile of the halo of the low surface brightness (LSB) galaxy NGC 5963 to place an upper limit on the dark matter cross section for collisions with velocities $\sim 150$ km s$^{-1}$, i.e. at the low velocity regime. After calibrating against cosmological simulations, we found that the large inferred dark matter concentration and central dark matter density in NGC 5963 are inconsistent with an effective collisional cross section per unit of mass >0.2 cm^2/g. Corrections were applied in order to account for reduction of the core by the adiabatic contraction caused by cooling baryons. Our limits that involve a number of simplifying, but always conservative, assumptions, exclude the last permitted interval for velocity-dependent cross sections to explain the flat density core in LSB galaxies. Implications for the nature of dark matter are also discussed.

The density profiles of dark halos are usually inferred from the rotation curves of disk galaxies based on the assumption that the gas is a good tracer of the gravitational potential of the galaxies. Some authors have suggested that magnetic pinching forces could alter significantly the rotation curves of spiral galaxies. In contrast to other studies which have concentrated in the vertical structure of the disk, here we focus on the problem of magnetic confinement in the radial direction to bound the magnetic effects on the HI rotation curves. It is shown that azimuthal magnetic fields can hardly speed up the HI disk of galaxies as a whole. In fact, based on virial constraints we show that the contribution of galactic magnetic fields to the rotation curves cannot be larger than ~10 km/s at the outermost point of HI detection, if the galaxies did not contain dark matter at all, and up to 20 km/s in the conventional dark halo scenario. The procedure to estimate the maximum effect of magnetic fields is general and applicable to any particular galaxy disk. The inclusion of the surface terms, namely the intergalactic (thermal, magnetic or ram) pressure, does not change our conclusions. Other problems related with the magnetic alternative to dark halos are highlighted. The relevance of magnetic fields in the cuspy problem of dark halos is also discussed.

We speculate that the dark halos of dwarf galaxies and low surface brightness galaxies soften their central cusps by the decay of a fraction of cold dark matter (CDM) particles to a stable particle with a recoiling velocity of a few tens km s$^{-1}$, after they have driven the formation of galactic halos. This process, however, does not necessarily produce a significant reduction of the central dark matter density of satellite dwarf spheroidals like Draco or Fornax. It is shown that the recovered distribution of concentration parameters $c$ for the initial (before decay) Navarro-Frenk-White halos, is in good agreement with CDM predictions. Other interesting potentials of unstable CDM are highlighted.

(Abridged) We discuss the nonlinear development of the isobaric mode of thermal instability (TI) in the context of the atomic interstellar medium (ISM), in order to assess the ability of TI to establish a well-segregated multi-phase structure in the turbulent ISM. The key parameter is the ratio of the cooling time to the dynamical crossing time $\eta$. The isobaric density perturbations of sizes $\gtrsim 15$ pc in media with mean density $\sim 1$ cm$^{-3}$ develop inflow motions with Mach number larger than 0.5 and a shock that propagates outwards from the condensation and brings the surrounding medium out of thermal equilibrium. The time for the dynamical transient state to subside ranges from 4 to 30 Myr for initial density perturbations of 20% and sizes 3 to 75 pc. By the time the condensations have formed, a substantial fraction of the mass is still traversing the unstable range and, therefore, clouds formed by TI should be bounded by accreting gas traversing the unstable range, rather than by sharp transitions to the stable warm phase. In the presence of random velocity forcing the condensation process can be suppressed for arbitrarily long times if the forcing causes a moderate rms Mach number ($\gtrsim 0.3$) and extends to small enough scales that $\eta>1$. We suggest that these mechanisms may be at the origin of the relatively large amounts of gas mass in the unstable regime. We remark that in the (stable) warm diffuse medium, $\eta>1$ for velocity perturbations of scales up to several pc. The flow's response to them is thus nearly adiabatic, and relatively weakly compressible, consistent with observations that suggest a nearly Kolmogorov power spectrum in this medium.

We discuss the temperature distribution in a two-dimensional, thermally unstable numerical simulation of the warm and cold gas in the Galactic disk, including the magnetic field, self-gravity, the Coriolis force, stellar energy injection and a realistic cooling function. We find that ~50% of the turbulent gas mass has temperatures in what would be the thermally unstable range if thermal instability were to be considered alone. This appears to be a consequence of there being many other forces at play than just thermal pressure. We also point out that a bimodal temperature pdf is a simple consequence of the form of the interstellar cooling function and is not necessarily a signature of discontinuous phase transitions.

The dynamical friction experienced by a body moving in a gaseous medium is different from the friction in the case of a collisionless stellar system. Here we consider the orbital evolution of a gravitational perturber inside a gaseous sphere using three-dimensional simulations, ignoring however self-gravity. The results are analysed in terms of a `local' formula with the associated Coulomb logarithm taken as a free parameter. For forced circular orbits, the asymptotic value of the component of the drag force in the direction of the velocity is a slowly varying function of the Mach number in the range 1.0-1.6. The dynamical friction timescale for free decay orbits is typically only half as long as in the case of a collisionless background, which is in agreement with E.C. Ostriker's recent analytic result. The orbital decay rate is rather insensitive to the past history of the perturber. It is shown that, similar to the case of stellar systems, orbits are not subject to any significant circularization. However, the dynamical friction timescales are found to increase with increasing orbital eccentricity for the Plummer model, whilst no strong dependence on the initial eccentricity is found for the isothermal sphere.