BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

J. A. Sellwood

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

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COastro-phastro-ph.IM

Catalog footprint

What is connected

38works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

38 published item(s)

preprint2022arXiv

The ultra-diffuse galaxy AGC 114905 needs dark matter

Recent 21 cm line observations of the ultra-diffuse galaxy AGC~114905 indicate a rotating disc largely supported against gravity by orbital motion, as usual. Remarkably, this study has revealed that the form and amplitude of the HI rotation curve is completely accounted for by the observed distribution of baryonic matter, stars and neutral gas, implying that no dark halo is required. It is surprising to find a DM-free galaxy for a number of reasons, one being that a bare Newtonian disk having low velocity dispersion would be expected to be unstable to both axi- and non-axisymmetric perturbations that would change the structure of the disc on a dynamical timescale, as has been known for decades. We present $N$-body simulations of the DM-free model, and one having a low-density DM halo, that confirm this expectation: the disc is chronically unstable to just such instabilities. Since it is unlikely that a galaxy that is observed to have a near-regular velocity pattern would be unstable, our finding calls into question the suggestion that the galaxy may lack, or have little, dark matter. We also show that if the inclination of this near face-on system has been substantially overestimated, the consequent increased amplitude of the rotation curve would accommodate a halo massive enough for the galaxy to be stable.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Spiral Instabilities in N-body Simulations: Emergence from noise II

An earlier paper presented the potentially significant discovery that disturbances in simplified simulations of a stellar disc model that was predicted to be stable in linear theory grew to large amplitude over a long period of time. The ultimate appearance of true instabilities was attributed to non-linear scattering by a succession of collective waves excited by shot noise from the finite number of particles. The paper concluded that no finite number of particles, however large, could mimic a smooth disc. As this surprising finding has been challenged as an artifact of the numerical scheme employed, we here present a new calculation of the same model using a different grid geometry that confirms the original behaviour.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Three Mechanisms for Bar Thickening

We present simulations of bar-unstable stellar discs in which the bars thicken into box/peanut shapes. Detailed analysis of the evolution of each model revealed three different mechanisms for thickening the bars. The first mechanism is the well-known buckling instability, the second is the vertical excitation of bar orbits by passage through the 2:1 vertical resonance, and the third is a gradually increasing fraction of bar orbits trapped into this resonance. Since bars in many galaxies may have formed and thickened long ago, we have examined the models for fossil evidence in the velocity distribution of stars in the bar, finding a diagnostic to discriminate between a bar that had buckled from the other two mechanisms.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Bar instability in disk-halo systems

We show that the exponential growth rate of a bar in a stellar disk is substantially greater when the disk is embedded in a live halo than in a rigid one having the same mass distribution. We also find that the vigor of the instability in disk-halo systems varies with the shape of the halo velocity ellipsoid. Disks in rigid halos that are massive enough to be stable by the usual criteria, quickly form bars in isotropic halos and much greater halo mass is needed to avoid a strong bar; thus stability criteria derived for disks in rigid halos do not apply when the halo is responsive. The study presented here is of an idealized family of models with near uniform central rotation and that lack an extended halo; we present more realistic models with extended halos in a companion paper. The puzzle presented by the absence of strong bars in some galaxies having gently rising inner rotation curves is compounded by the results presented here.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

DiskFit: a code to fit simple non-axisymmetric galaxy models either to photometric images or to kinematic maps

This posting announces public availability of version 1.2 of the DiskFit software package developed by the authors, which may be used to fit simple non-axisymmetric models either to images or to velocity fields of disk galaxies. Here we give an outline of the capability of the code and provide the link to downloading executables, the source code, and a comprehensive on-line manual. We argue that in important respects the code is superior to rotcur for fitting kinematic maps and to galfit for fitting multi-component models to photometric images.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Made-To-Measure Models of Self-Similar Triaxial Halos with Steep Inner Density Gradients

We use the Made-to-Measure method to construct N-body realizations of self-similar, triaxial ellipsoidal halos having cosmologically realistic density profiles. Our implementation parallels previous work with a few numerical refinements, but we show that orbital averaging is an intrinsic feature of the force of change equation and argue that additional averaging or smoothing schemes are redundant. We present models having the Einasto radial mass profile that range from prolate to strongly triaxial. We use a least-squares polynomial fit to the expansion coefficients to obtain an analytical representation of the particle density from which we derive density contours and eccentricity profiles more efficiently than by the usual particle smoothing techniques. We show that our N-body realizations both retain their shape in unconstrained evolution and recover it after large amplitude perturbations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Relaxation in N-body simulations of spherical systems

I present empirical measurements of the rate of relaxation in N-body simulations of stable spherical systems and distinguish two separate types of relaxation: energy diffusion that is largely independent of particle mass, and energy exchange between particles of differing masses. While diffusion is generally regarded as a Fokker-Planck process, it can equivalently be viewed as the consequence of collective oscillations that are driven by shot noise. Empirical diffusion rates scale as N^{-1} in inhomogeneous models, in agreement with Fokker-Planck predictions, but collective effects cause relaxation to scale more nearly as N^{-1/2} in the special case of a uniform sphere. I use four different methods to compute the gravitational field, and a 100-fold range in the numbers of particles in each case. I find the rate at which energy is exchanged between particles of differing masses does not depend at all on the force determination method, but I do find the energy diffusion rate is marginally lower when a field method is used. The relaxation rate in 3D is virtually independent of the method used because it is dominated by distant encounters; any method to estimate the gravitational field that correctly captures the contributions from distant particles must also capture their statistical fluctuations and the collective modes they drive.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The RINGS Survey: High-Resolution H-alpha Velocity Fields of Nearby Spiral Galaxies with the SALT Fabry-Perot

We have obtained high-spatial-resolution spectrophotometric data on several nearby spiral galaxies with the Southern African Large Telescope (SALT) Fabry-Pérot interferometer on the Robert Stobie Spectrograph (RSS) as a part of the RSS Imaging spectroscopy Nearby Galaxy Survey (RINGS). We have successfully reduced two tracks of Fabry-Pérot data for the galaxy NGC 2280 to produce a velocity field of the H-alpha line of excited hydrogen. We have modeled these data with the DiskFit modeling software and found these models to be in excellent agreement both with previous measurements in the literature and with our lower-resolution HI velocity field of the same galaxy. Despite this good agreement, small regions exist where the difference between the H-alpha and HI velocities is larger than would be expected from typical dispersions. We investigate these regions of high velocity difference and offer possible explanations for their existence.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

GALAXY package for N-body simulation

This posting announces public availability of the GALAXY software package developed by the author over the past 40 years. It is a highly efficient code for the evolution of (almost) isolated, collisionless stellar systems, both disk-like and ellipsoidal. In addition to the N-body code galaxy, which offers eleven different methods to compute the gravitational accelerations, the package also includes sophisticated set-up and analysis software. This paper gives an outline of the contents of the package and provides links to the source code and a comprehensive on-line manual. While not as versatile as tree codes, the particle-mesh methods in this package are shown, for certain restricted applications, to be between 50 and 200 times faster than a widely-used tree code.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Re-interpretation of "Bar slowdown and the distribution of dark matter in barred galaxies" by Athanassoula

Athanassoula (2014) has claimed that measurements of the ratio of corotation radius to bar length in galaxies do not place a constraint on the disk to halo mass ratio. Her conclusion was based on a series of simulations published by Athanassoula et al. (2013). Here we show that these results are, in fact, consistent with previous work on the slow down of bars due to dynamical friction because gas inflow rearranges the disk mass before the bar settles. It therefore remains true that a bar rotating sufficiently fast that corotation is not far beyond the bar end requires a near maximum disk.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Smoothing Rotation Curves and Mass Profiles

We show that spiral activity can erase pronounced features in disk galaxy rotation curves. We present simulations of growing disks, in which the added material has a physically motivated distribution, as well as other examples of physically less realistic accretion. In all cases, attempts to create unrealistic rotation curves were unsuccessful because spiral activity rapidly smoothed away features in the disk mass profile. The added material was redistributed radially by the spiral activity, which was itself provoked by the density feature. In the case of a ridge-like feature in the surface density profile, we show that two unstable spiral modes develop, and the associated angular momentum changes in horseshoe orbits remove particles from the ridge and spread them both inwards and outwards. This process rapidly erases the density feature from the disk. We also find that the lack of a feature when transitioning from disk to halo dominance in the rotation curves of disk galaxies, the so called "disk-halo conspiracy", could also be accounted for by this mechanism. We do not create perfectly exponential mass profiles in the disk, but suggest that this mechanism contributes to their creation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Formation of Spiral Galaxies: Adiabatic Compression with Young's Algorithm and the Relation of Dark Matter Haloes to Their Primordial Antecedents

We utilize Young's algorithm to model the adiabatic compression of the dark matter haloes of galaxies in the THINGS survey to determine the relationship between the halo fit to the rotation curve and the corresponding primordial halo prior to compression. Young's algorithm conserves radial action and angular momentum, resulting in less halo compression than more widely utilized approximations. We find that estimates of the parameters of the NFW haloes fit to the current dark matter distribution systematically overestimate the concentration and underestimate the virial velocity of the corresponding primordial halo. It is the latter that is predicted by dark matter simulations; so accounting for compression is a necessary step for evaluating whether massive galaxies are consistent with dark matter-only simulations. The inferred primordial haloes broadly follow the c-V200 relation expected in a LCDM cosmogony, but often scatter to lower concentrations. We are unable to obtain fits at all for those galaxies whose current dark matter haloes are poorly described by the NFW form. We thus find a mixed bag: some galaxies are reasonably well described by adiabatic compression within a primordial NFW halo, while others require additional mechanisms that reduce the density of dark matter below the primordial initial condition.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The RINGS Survey I: Halpha and HI Velocity Maps of Galaxy NGC 2280

Precise measurements of gas kinematics in the disk of a spiral galaxy can be used to estimate its mass distribution. The Southern African Large Telescope (SALT) has a large collecting area and field of view, and is equipped with a Fabry-Perot interferometer that can measure gas kinematics in a galaxy from the Halpha line. To take advantage of this capability, we have constructed a sample of 19 nearby spiral galaxies, the RSS Imaging and Spectroscopy Nearby Galaxy Survey (RINGS), as targets for detailed study of their mass distributions and have collected much of the needed data. In this paper, we present velocity maps produced from Halpha Fabry-Perot interferometry and HI aperture synthesis for one of these galaxies, NGC 2280, and show that the two velocity measurements are generally in excellent agreement. Minor differences can mostly be attributed to the different spatial distributions of the excited and neutral gas in this galaxy, but we do detect some anomalous velocities in our Halpha velocity map of the kind that have previously been detected in other galaxies. Models produced from our two velocity maps agree well with each other and our estimates of the systemic velocity and projection angles confirm previous measurements of these quantities for NGC 2280.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Transient spirals as superposed instabilities

We present evidence that recurrent spiral activity, long manifested in simulations of disk galaxies, results from the super-position of a few transient spiral modes. Each mode lasts between 5 and 10 rotations at its corotation radius where its amplitude is greatest. The scattering of stars as each wave decays takes place over narrow ranges of angular momentum, causing abrupt changes to the impedance of the disk to subsequent traveling waves. Partial reflections of waves at these newly created features, allows new standing-wave instabilities to appear that saturate and decay in their turn, scattering particles at new locations, creating a recurring cycle. The spiral activity causes the general level of random motion to rise, gradually decreasing the ability of the disk to support further activity unless the disk contains a dissipative gas component from which stars form on near-circular orbits. We also show that this interpretation is consistent with the behavior reported in other recent simulations with low mass-disks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Relaxation in N-body simulations of disk galaxies

I use N-body simulations with two mass species of particles to demonstrate that disk galaxy simulations are subject to collisional relaxation at a higher rate than is widely assumed. Relaxation affects the vertical thickness of the disk most strongly, and drives the velocity ellipsoid to a moderately flattened shape similar to that observed for disk stars in the solar neighborhood. The velocity ellipsoid in simulations with small numbers of particles quickly approaches this shape, but shot noise also dominates the in-plane behavior. Simulations with higher, but reachable, numbers of particles relax slowly enough to be considered collisionless, allowing the in-plane dispersions to rise due to spiral activity without heating the vertical motions. Relaxation may have affected many previously published simulations of the formation and evolution of galaxy disks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Secular Evolution in Disk Galaxies

Disk galaxies evolve over time through processes that may rearrange both the radial mass profile and the metallicity distribution within the disk. This review of such slow changes is largely, though not entirely, restricted to internally-driven processes that can be distinguished from evolution driven by galaxy interactions. It both describes our current understanding of disk evolution, and identifies areas where more work is needed. Stellar disks are heated through spiral scattering, which increases random motion components in the plane, while molecular clouds redirect some fraction of the random energy into vertical motion. The recently discovered process of radial migration at the corotation resonance of a transient spiral mode does not alter the underlying structure of the disk, since it neither heats the disk nor causes it to spread, but it does have a profound effect on the expected distribution of metallicities among the disk stars. Bars in disks are believed to be major drivers of secular evolution through interactions with the outer disk and with the halo. Once the material that makes up galaxy disks is converted into stars, their overall angular momentum distribution cannot change by much, but that of the gas is generally far more liable to rearrangement, allowing rings and pseudo-bulges to form. While simulations are powerful tools from which we have learned a great deal, those of disks may suffer from collisional relaxation that requires some results to be interpreted with caution.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Test for Radial Mixing Using Local Star Samples

We use samples of local main-sequence stars to show that the radial gradient of [Fe/H] in the thin disk of the Milky Way decreases with mean effective stellar temperature. Many of these stars are visiting the solar neighborhood from the inner and outer Galaxy. We use the angular momentum of each star about the Galactic center to determine the guiding center radius and to eliminate the effects of epicyclic motion, which would otherwise blur the estimated gradients. We interpret the effective temperature as a proxy for mean age, and conclude that the decreasing gradient is consistent with the predictions of radial mixing due to transient spiral patterns. We find some evidence that the trend of decreasing gradient with increasing mean age breaks to a constant gradient for samples of stars whose main-sequence life-times exceed the likely age of the thin disk.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An Algorithm for Preferential Selection of Spectroscopic Targets in LEGUE

We describe a general target selection algorithm that is applicable to any survey in which the number of available candidates is much larger than the number of objects to be observed. This routine aims to achieve a balance between a smoothly-varying, well-understood selection function and the desire to preferentially select certain types of targets. Some target-selection examples are shown that illustrate different possibilities of emphasis functions. Although it is generally applicable, the algorithm was developed specifically for the LAMOST Experiment for Galactic Understanding and Exploration (LEGUE) survey that will be carried out using the Chinese Guo Shou Jing Telescope. In particular, this algorithm was designed for the portion of LEGUE targeting the Galactic halo, in which we attempt to balance a variety of science goals that require stars at fainter magnitudes than can be completely sampled by LAMOST. This algorithm has been implemented for the halo portion of the LAMOST pilot survey, which began in October 2011.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

LAMOST Experiment for Galactic Understanding and Exploration (LEGUE) The survey science plan

We describe the current plans for a spectroscopic survey of millions of stars in the Milky Way galaxy using the Guo Shou Jing Telescope (GSJT, formerly the Large Area Multi-Object Spectroscopic Telescope - LAMOST). The survey will obtain spectra for 2.5 million stars brighter than $r<19$ during dark/grey time, and 5 million stars brighter than $r<17$ or $J<16$ on nights that are moonlit or have low transparency. The survey will begin in fall of 2012, and will run for at least four years. The telescope design constrains the optimal declination range for observations to $10^\circ<δ<50^\circ$, and site conditions lead to an emphasis on stars in the direction of the Galactic anticenter. The survey is divided into three parts with different target selection strategies: disk, anticenter, and spheroid. The resulting dataset will be used to study the merger history of the Milky Way, the substructure and evolution of the disks, the nature of the first generation of stars through identification of the lowest metallicity stars, and star formation through study of open clusters and the OB associations. Detailed design of the LEGUE survey will be completed after a review of the results of the pilot survey in summer 2012.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Radial Migration in Galactic Thick Discs

We present a study of the extent to which the Sellwood & Binney radial migration of stars is affected by their vertical motion about the midplane. We use both controlled simulations in which only a single spiral mode is excited, as well as slightly more realistic cases with multiple spiral patterns and a bar. We find that rms angular momentum changes are reduced by vertical motion, but rather gradually, and the maximum changes are almost as large for thick disc stars as for those in a thin disc. We find that particles in simulations in which a bar forms suffer slightly larger angular momentum changes than in comparable cases with no bar, but the cumulative effect of multiple spiral events still dominates. We have determined that vertical action, and not vertical energy, is conserved on average during radial migration.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Searching for Non-axisymmetries in NGC 6503: A Weak End-on Bar

The isolation, simple apparent structure, and low luminosity of the nearby spiral galaxy NGC 6503 make it an ideal candidate for an in-depth kinematic and photometric study. We introduce a new publicly available code, DiskFit, that implements procedures for fitting non-axisymmetries in either kinematic or photometric data. We use DiskFit to analyze new Halpha and CO velocity field data as well as HI kinematics from Greisen et al. to search for non-circular motions in the disc of NGC 6503. We find NGC 6503 to have remarkably regular gas kinematics that are well-described by rotation only. We also use DiskFit and a new Ks-band image of NGC 6503 to constrain photometric models of the disc, bar and bulge. We find the galaxy to be photometrically dominated by the disc. We find NGC 6503 to contain a faint bar and an exponential bulge which together contribute only ~5% of the total galaxy light. The combination of our kinematic and photometric DiskFit models suggest that NGC 6503 contains a weak, end-on bar that may have produced its Type II surface brightness profile but is unlikely to be responsible for its strong sigma-drop.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Spiral Instabilities in N-body Simulations I: Emergence from Noise

The origin of spiral patterns in galaxies is still not fully understood. Similar features also develop readily in N-body simulations of isolated cool, collisionless disks, yet even here the mechanism has yet to be explained. In this series of papers, I present a detailed study of the origin of spiral activity in simulations in the hope that the mechanism that causes the patterns is also responsible for some of these features galaxies. In this first paper, I use a suite of highly idealized simulations of a linearly stable disk that employ increasing numbers of particles. While the amplitudes of initial non-axisymmetric features scale as the inverse square-root of the number of particles employed, the final amplitude of the patterns is independent of the particle number. I find that the amplitudes of non-axisymmetric disturbances grow in two distinct phases: slow growth occurs when the relative overdensity is below ~2%, but above this level the amplitude rises more rapidly. I show that all features, even of very low amplitude, scatter particles at the inner Lindblad resonance, changing the distribution of particles in the disk in such a way as to foster continued growth. Stronger scattering by larger amplitude waves provokes a vigorous instability that is a true linear mode of the modified disk.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamics of Disks and Warps

This chapter reviews theoretical work on the stellar dynamics of galaxy disks. All the known collective global instabilities are identified, and their mechanisms described in terms of local wave mechanics. A detailed discussion of warps and other bending waves is also given. The structure of bars in galaxies, and their effect on galaxy evolution, is now reasonably well understood, but there is still no convincing explanation for their origin and frequency. Spiral patterns have long presented a special challenge, and ideas and recent developments are reviewed. Other topics include scattering of disk stars and the survival of thin disks.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Velocity-space substructure from nearby RAVE and SDSS stars

We extract a sample of disc stars within 200 pc of the Sun from the RAVE and SDSS surveys. Distances are estimated photometrically and proper motions are from ground-based data. We show that the velocity-space substructure first revealed in the Geneva-Copenhagen sample is also present in this completely independent sample. We also evaluate action-angle variables for these stars and show that the Hyades stream stars in these data are again characteristic of having been scattered at a Lindblad resonance. Unfortunately, analysis of such local samples can determine neither whether it is an inner or an outer Lindblad resonance, nor the multiplicity of the pattern.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Recent Lindblad Resonance in the Solar Neighbourhood

Stars in the solar neighbourhood do not have a smooth distribution of velocities. Instead, the distribution of velocity components in the Galactic plane manifests a great deal of kinematic substructure. Here I present an analysis in action-angle variables of the Geneva-Copenhagen survey of ~14,000 nearby F & G dwarfs with distances and full space motions. I show that stars in the so-called "Hyades stream" have both angle and action variables characteristic of their having been scattered at an inner Lindblad resonance of a rotating disturbance potential. This discovery seems to favour spiral patterns as recurrent, short-lived instabilities.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

New Developments in Spiral Structure Theory

After a short review of the principal theories of spiral structure in galaxies, I describe two new developments. First, it now seems clear that linear theory cannot yield a full description for the development of spiral patterns because N-body simulations suggest that non-linear effects manifest themselves at a relative overdensity of ~2%, which is well below the believed spiral amplitudes in galaxies. Second, I summarize the evidence that some stars in the solar neighborhood have been scattered at an inner Lindblad resonance. This evidence strongly supports a picture of spirals as a recurring cycle of transient instabilities, each caused by resonant scattering by a previous wave.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Spiral instabilities in N-body simulations

N-body simulations of disc galaxies that display recurrent transient spiral patterns are comparatively easy to construct, but are harder to understand. In this paper, I summarise the evidence from such experiments that the spiral patterns result from a recurrent spiral instability cycle. Each wave starts as rapidly growing, small-amplitude instability caused by a deficiency of particles at a particular angular momentum. The resulting large-amplitude wave creates, through resonant scattering, the conditions needed to precipitate a new instability.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Lifetimes of Spiral Patterns in Disc Galaxies

The rate of internally-driven evolution of galaxy discs is strongly affected by the lifetimes of the spiral patterns they support. Evolution is much faster if the spiral patterns are recurrent short-lived transients rather than long-lived, quasi-steady features. As rival theories are still advocated based on these two distinct hypotheses, I review the evidence that bears on the question of the lifetimes of spiral patterns in galaxies. Observational evidence from external galaxies is frustratingly inconclusive, but the velocity distribution in the solar neighbourhood is more consistent with the transient picture. I present simulations of galaxy models that have been proposed to support quasi-steady, two-arm spiral modes that in fact evolve quickly due to multi-arm instabilities. I also show that all simulations to date manifest short-lived patterns, despite claims to the contrary. Thus the transient hypothesis is favoured by both numerical results and the velocity distribution in the solar neighbourhood.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Stochasticity in N-body Simulations of Disc Galaxies

We demonstrate that the chaotic nature of N-body systems can lead to macroscopic variations in the evolution of collisionless simulations containing rotationally supported discs. The unavoidable stochasticity that afflicts all simulations generally causes mild differences between the evolution of similar models but, in order to illustrate that this is not always true, we present a case that shows extreme bimodal divergence. The divergent behaviour occurs in two different types of code and is independent of all numerical parameters. We identify and give explicit illustrations of several sources of stochasticity, and also show that macroscopic variations in the evolution can originate from differences at the round-off error level. We obtain somewhat more consistent results from simulations in which the halo is set up with great care compared with those started from more approximate equilibria, but we have been unable to eliminate diverging behaviour entirely because the main sources of stochasticity are intrinsic to the disc. We show that the divergence is only temporary and that halo friction is merely delayed, for a substantial time in some cases. We argue that the delays are unlikely to arise in real galaxies, and that our results do not affect dynamical friction constraints on halo density. Stochastic variations in the evolution are inevitable in all simulations of disc-halo systems, irrespective of how they were created, although their effect is generally far less extreme than we find here. The possibility of divergent behaviour complicates comparison of results from different workers.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Dark Matter Densities during the Formation of the First Stars and in Dark Stars

The first stars in the universe form inside $\sim 10^6 M_\odot$ dark matter (DM) haloes whose initial density profiles are laid down by gravitational collapse in hierarchical structure formation scenarios. During the formation of the first stars in the universe, the baryonic infall compresses the dark matter further. The resultant dark matter density is presented here, using an algorithm originally developed by Young to calculate changes to the profile as the result of adiabatic infall in a spherical halo model; the Young prescription takes into account the non-circular motions of halo particles. The density profiles obtained in this way are found to be within a factor of two of those obtained using the simple adiabatic contraction prescription of Blumenthal et al. Our results hold regardless of the nature of the dark matter or its interactions and rely merely on gravity. If the dark matter consists of weakly interacting massive particles, which are their own antiparticles, their densities are high enough that their annihilation in the first protostars can indeed provide an important heat source and prevent the collapse all the way to fusion. In short, a ``Dark Star'' phase of stellar evolution, powered by DM annihilation, may indeed describe the first stars in the universe.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Dark Matter Density in Disk Galaxies

I show that the predicted densities of the inner dark matter halos in LCDM models of structure formation appear to be higher than estimates from real galaxies and constraints from dynamical friction on bars. This inconsistency would not be a problem for the LCDM model if physical processes that are omitted in the collisionless collapse simulations were able to reduce the dark matter density in the inner halos. I review the mechanisms proposed to achieve the needed density reduction.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Modeling Non-Circular Motions in Disk Galaxies: A Bar in NGC 2976

We give a brief description of a new model for non-circular motions in disk galaxy velocity fields, that does not invoke epicycles. We assume non-circular motions to stem from a bar-like or oval distortion to the potential, as could arise from a triaxial halo or a bar in the mass distribution of the baryons. We apply our model to the high-quality CO and Halpha kinematics of NGC 2976 presented by Simon et al. 2003; it fits the data as well as their model with unrealistic radial flows, but yields a steeper rotation curve. Our analysis and other evidence suggests that NGC 2976 hosts a bar, implying a large baryonic contribution to the potential and thus limiting the allowed dark matter halo density.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Modeling Non-Circular Motions in Disk Galaxies: Application to NGC 2976

We present a new procedure to fit non-axisymmetric flow patterns to 2-D velocity maps of spiral galaxies. We concentrate on flows caused by bar-like or oval distortions to the total potential that may arise either from a non-axially symmetric halo or a bar in the luminous disk. We apply our method to high-quality CO and Halpha data for the nearby, low-mass spiral NGC 2976 previously obtained by Simon et al., and find that a bar-like model fits the data at least as well as their model with large radial flows. We find supporting evidence for the existence of a bar in the baryonic disk. Our model suggests that the azimuthally averaged central attraction in the inner part of this galaxy is larger than estimated by these authors. It is likely that the disk is also more massive, which will limit the increase to the allowed dark halo density. Allowance for bar-like distortions in other galaxies may either increase or decrease the estimated central attraction.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Anomalously Weak Dynamical Friction in Halos

A bar rotating in a pressure-supported halo generally loses angular momentum and slows down due to dynamical friction. Valenzuela & Klypin report a counter-example of a bar that rotates in a dense halo with little friction for several Gyr, and argue that their result invalidates the claim by Debattista & Sellwood that fast bars in real galaxies require a low halo density. We show that it is possible for friction to cease for a while should the pattern speed of the bar fluctuate upward. The reduced friction is due to an anomalous gradient in the phase-space density of particles at the principal resonance created by the earlier evolution. The result obtained by Valenzuela & Klypin is probably an artifact of their adaptive mesh refinement method, but anyway could not persist in a real galaxy. The conclusion by Debattista & Sellwood still stands.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Bar-halo Friction in Galaxies I: Scaling Laws

It has been known for some time that rotating bars in galaxies slow due to dynamical friction against the halo. However, recent attempts to use this process to place constraints on the dark matter density in galaxies and possibly also to drive dark matter out of the center have been challenged. This paper uses simplified numerical experiments to clarify several aspects of the friction mechanism. I explicitly demonstrate the Chandrasekhar scaling of the friction force with bar mass, halo density, and halo velocity dispersion. I present direct evidence that exchanges between the bar and halo orbits at major resonances are responsible for friction and study both individual orbits and the net changes at these resonances. I also show that friction alters the phase space density of particles in the vicinity of a major resonance, which is the reason the magnitude of the friction force depends on the prior evolution. I demonstrate that bar slow down can be captured correctly in simulations having modest spatial resolution and practicable numbers of particles. Subsequent papers in this series delineate the dark matter density that can be tolerated in halos of different density profiles.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Dynamical Friction and the Distribution of Dark Matter in Barred Galaxies

We use fully self-consistent N-body simulations of barred galaxies to show that dynamical friction from a dense dark matter halo dramatically slows the rotation rate of bars. Our result supports previous theoretical predictions for a bar rotating within a massive halo. On the other hand, low density halos, such as those required for maximum disks, allow the bar to continue to rotate at a high rate. There is somewhat meager observational evidence indicating that bars in real galaxies do rotate rapidly and we use our result to argue that dark matter halos must have a low central density in all high surface brightness disk galaxies, including the Milky Way. Bars in galaxies that have larger fractions of dark matter should rotate slowly, and we suggest that a promising place to look for such candidate objects is among galaxies of intermediate surface brightness.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Instabilities of a family of oblate stellar spheroids

We have examined the stability of a sequence of oblate elliptical galaxy models having the Stackel form suggested by Kuz'min & Kutuzov. We have employed the 2-integral DFs given by Dejonghe & de Zeeuw for which flattened non-rotating models are characterized by counter-streaming motion and are radially cool; we introduce net rotation in some models by changing the sign of the z-component of angular momentum for a fraction of the particles. We have found that all non-rotating and slowly rotating members of this sequence rounder than E7 are stable, and that even maximally rotating models rounder than E4 are stable. In the absence of strong rotation, the most disruptive instability, and the last to be stabilized by increasing thickness, is a lopsided (m=1) mode. This instability appears to be driven by counter-rotation in radially cool models. Its vigour is lessened as rotation is increased, but it remains strong even in models with net angular momentum 90% of that of a maximally rotating model before finally disappearing in maximally rotating models. Strongly rotating models are more unstable to bar-forming modes which afflict maximally rotating models with c/a >~ 0.5, but this mode is quickly stabilized by moderate fractions of counter-rotating particles. Bending instabilities appear not to be very important; they are detectable in the inner parts of the flatter models, but are less vigorous and more easily stabilized than the lop-sided or bar modes in every case. We briefly discuss the possible relevance of the lop-sided instability to the existence of many lop-sided disk galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Secular Evolution in Barred Galaxies

A strong bar rotating within a massive halo should lose angular momentum to the halo through dynamical friction, as predicted by Weinberg. We have conducted fully self-consistent, numerical simulations of barred galaxy models with a live halo population and find that bars are indeed braked very rapidly. Specifically, we find that the bar slows sufficiently within a few rotation periods that the distance from the centre to co-rotation is more than twice the semi-major axis of the bar. Observational evidence (meagre) for bar pattern speeds seems to suggest that this ratio typically lies between 1.2 to 1.5 in real galaxies. We consider, a number of possible explanations for this discrepancy between theoretical prediction and observation, and conclude that no conventional alternative seems able to account for it.

0 citations0 reviewsNo code linkedOpen access