Measuring the local dark matter density

We examine systematic problems in determining the local matter density from the vertical motion of stars, i.e. the 'Oort limit'. Using collisionless simulations and a Monte Carlo Markov Chain technique, we determine the data quality required to detect local dark matter at its expected density. We find that systematic errors are more important than observational errors and apply our technique to Hipparcos data to reassign realistic error bars to the local dark matter density.Comment: 3 pages, 1 figure, to be published in "Hunting for the Dark: The Hidden Side of Galaxy Formation", Malta, 19-23 Oct. 2009, eds. V.P. Debattista & C.C. Popescu, AIP Conf. Se

A Semi-Analytic dynamical friction model that reproduces core stalling

We present a new semi-analytic model for dynamical friction based on Chandrasekhar's formalism. The key novelty is the introduction of physically motivated, radially varying, maximum and minimum impact parameters. With these, our model gives an excellent match to full N-body simulations for isotropic background density distributions, both cuspy and shallow, without any fine-tuning of the model parameters. In particular, we are able to reproduce the dramatic core-stalling effect that occurs in shallow/constant density cores, for the first time. This gives us new physical insight into the core-stalling phenomenon. We show that core stalling occurs in the limit in which the product of the Coulomb logarithm and the local fraction of stars with velocity lower than the infalling body tends to zero. For cuspy backgrounds, this occurs when the infalling mass approaches the enclosed background mass. For cored backgrounds, it occurs at larger distances from the centre, due to a combination of a rapidly increasing minimum impact parameter and a lack of slow moving stars in the core. This demonstrates that the physics of core-stalling is likely the same for both massive infalling objects and low-mass objects moving in shallow density backgrounds. We implement our prescription for dynamical friction in the direct summation code NBODY6 as an analytic correction for stars that remain within the Roche volume of the infalling object. This approach is computationally efficient, since only stars in the inspiralling system need to be evolved with direct summation. Our method can be applied to study a variety of astrophysical systems, including young star clusters orbiting near the Galactic Centre; globular clusters moving within the Galaxy; and dwarf galaxies orbiting within dark matter halos.Comment: 16 pages, 21 figures, Accepted for publication in MNRA

Methane cracking over cobalt molybdenum carbides

The catalytic behaviour of Co3Mo3C, Co6Mo6C, Co3Mo3N and Co6Mo6N for methane cracking has been studied to determine the relationship between the methane cracking activity and the chemical composition. The characterisation of post-reaction samples showed a complex phase composition with the presence of Co3Mo3C, α-Co and β-Mo2C as catalytic phases and the deposition of different forms of carbon during reaction

The mass distribution of the Fornax dSph: constraints from its globular cluster distribution

Uniquely among the dwarf spheroidal (dSph) satellite galaxies of the Milky Way, Fornax hosts globular clusters. It remains a puzzle as to why dynamical friction has not yet dragged any of Fornax's five globular clusters to the centre, and also why there is no evidence that any similar star cluster has been in the past (for Fornax or any other dSph). We set up a suite of 2800 N-body simulations that sample the full range of globular-cluster orbits and mass models consistent with all existing observational constraints for Fornax. In agreement with previous work, we find that if Fornax has a large dark-matter core then its globular clusters remain close to their currently observed locations for long times. Furthermore, we find previously unreported behaviour for clusters that start inside the core region. These are pushed out of the core and gain orbital energy, a process we call 'dynamical buoyancy'. Thus a cored mass distribution in Fornax will naturally lead to a shell-like globular cluster distribution near the core radius, independent of the initial conditions. By contrast, CDM-type cusped mass distributions lead to the rapid infall of at least one cluster within \Delta t = 1-2Gyr, except when picking unlikely initial conditions for the cluster orbits (\sim 2% probability), and almost all clusters within \Delta t = 10Gyr. Alternatively, if Fornax has only a weakly cusped mass distribution, dynamical friction is much reduced. While over \Delta t = 10Gyr this still leads to the infall of 1-4 clusters from their present orbits, the infall of any cluster within \Delta t = 1-2Gyr is much less likely (with probability 0-70%, depending on \Delta t and the strength of the cusp). Such a solution to the timing problem requires that in the past the globular clusters were somewhat further from Fornax than today; they most likely did not form within Fornax, but were accreted.Comment: 12 pages, 8 figures, 3 tables, submitted to MNRA

A semi-analytic dynamical friction model for cored galaxies

We present a dynamical friction model based on Chandrasekhar's formula that reproduces the fast inspiral and stalling experienced by satellites orbiting galaxies with a large constant density core. We show that the fast inspiral phase does not owe to resonance. Rather, it owes to the background velocity distribution function for the constant density cores being dissimilar from the usually-assumed Maxwellian distribution. Using the correct background velocity distribution function and the semi-analytic model from Petts et al. (2015), we are able to correctly reproduce the infall rate in both cored and cusped potentials. However, in the case of large cores, our model is no longer able to correctly capture core-stalling. We show that this stalling owes to the tidal radius of the satellite approaching the size of the core. By switching off dynamical friction when rt(r) = r (where rt is the tidal radius at the satellite's position) we arrive at a model which reproduces the N-body results remarkably well. Since the tidal radius can be very large for constant density background distributions, our model recovers the result that stalling can occur for Ms/Menc << 1, where Ms and Menc are the mass of the satellite and the enclosed galaxy mass, respectively. Finally, we include the contribution to dynamical friction that comes from stars moving faster than the satellite. This next-to-leading order effect becomes the dominant driver of inspiral near the core region, prior to stalling.Comment: 13 pages, 12 figures, resubmitted to MNRAS after responding to feedback from the refere

Enhancement of field renormalization in scalar theories via functional renormalization group

The flow equations of the Functional Renormalization Group are applied to the O(N)-symmetric scalar theory, for N=1 and N=4, in four Euclidean dimensions, d=4, to determine the effective potential and the renormalization function of the field in the broken phase. In our numerical analysis, the infrared limit, corresponding to the vanishing of the running momentum scale in the equations, is approached to obtain the physical values of the parameters by extrapolation. In the N=4 theory a non-perturbatively large value of the physical renormalization of the longitudinal component of the field is observed. The dependence of the field renormalization on the UV cut-off and on the bare coupling is also investigated.Comment: 20 pages, 7 figures. To appear in Physical Review

The influence of Massive Black Hole Binaries on the Morphology of Merger Remnants

Massive black hole (MBH) binaries, formed as a result of galaxy mergers, are expected to harden by dynamical friction and three-body stellar scatterings, until emission of gravitational waves (GWs) leads to their final coalescence. According to recent simulations, MBH binaries can efficiently harden via stellar encounters only when the host geometry is triaxial, even if only modestly, as angular momentum diffusion allows an efficient repopulation of the binary loss cone. In this paper, we carry out a suite of N-body simulations of equal-mass galaxy collisions, varying the initial orbits and density profiles for the merging galaxies and running simulations both with and without central MBHs. We find that the presence of an MBH binary in the remnant makes the system nearly oblate, aligned with the galaxy merger plane, within a radius enclosing 100 MBH masses. We never find binary hosts to be prolate on any scale. The decaying MBHs slightly enhance the tangential anisotropy in the centre of the remnant due to angular momentum injection and the slingshot ejection of stars on nearly radial orbits. This latter effect results in about 1% of the remnant stars being expelled from the galactic nucleus. Finally, we do not find any strong connection between the remnant morphology and the binary hardening rate, which depends only on the inner density slope of the remnant galaxy. Our results suggest that MBH binaries are able to coalesce within a few Gyr, even if the binary is found to partially erase the merger-induced triaxiality from the remnant.Comment: 16 pages, 13 figures, 4 tables; accepted for publication in MNRA

The Local Dark Matter Density from SDSS-SEGUE G-dwarfs

We derive the local dark matter density by applying the integrated Jeans equation method from Silverwood et al. (2016) to SDSS-SEGUE G-dwarf data processed and presented by B\"udenbender et al. (2015). We use the MultiNest Bayesian nested sampling software to fit a model for the baryon distribution, dark matter and tracer stars, including a model for the 'tilt term' that couples the vertical and radial motions, to the data. The $\alpha$-young population from B\"udenbender et al. (2015) yields the most reliable result of $\rho_{\rm DM} = 0.46^{+0.07}_{-0.09}\, {{\rm GeV\, cm}^{-3}} = 0.012^{+0.001}_{-0.002}\, {{\rm M}_\odot \, {\rm pc}^{-3}}$. Our analyses yield inconsistent results for the $\alpha$-young and $\alpha$-old data, pointing to problems in the tilt term and its modelling, the data itself, the assumption of a flat rotation curve, or the effects of disequilibria.Comment: 17 pages, 10 figures, submitted to MNRA