Source author record

John Dubinski

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

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

17 published item(s)

preprint2015arXiv

Mergers in Galaxy Groups. II. The Fundamental Plane of Elliptical Galaxies

Observations consistently show that elliptical galaxies follow a tight "fundamental plane" scaling relation between size, mean surface brightness and velocity dispersion, with the form $R \propto σ^{a}μ^{b}$. This relation not only has very small (<0.05 dex) intrinsic scatter, but also has significantly different coefficients from the expected virial scaling (a "tilt"). We analyze hundreds of simulations of elliptical galaxies formed from mergers of spiral galaxies in groups to determine if the fundamental plane can emerge from multiple, mostly minor and hierarchical collisionless mergers. We find that these simulated ellipticals lie on a similar fundamental plane with $a \approx 1.7$ and $b \approx 0.3$. The scatter about this plane is not larger than observed, while the tilt is in the correct sense, although $a$ is larger than for typical observations. This supports the idea that collisionless mergers can contribute significantly to the tilt of the fundamental plane. The tilt is mainly driven by a mass-dependent dark matter fraction, such that more massive galaxies have larger dark matter fractions within $R_{e}$. We further discuss the origin of this mass-dependent dark matter fraction and its compatibility with strong lensing observations, as well as the links between the fundamental plane, dynamical masses and the virial theorem.

preprint2013arXiv

Mergers in Galaxy Groups - I. Structure and Properties of Elliptical Remnants

We present collisionless simulations of dry mergers in groups of three to twenty-five galaxies to test the hypothesis that elliptical galaxies form at the centers of such groups. Mock observations of the central remnants confirm their similarity to ellipticals, despite having no dissipational component. We vary the profile of the original spiral's bulge and find that ellipticals formed from spirals with exponential bulges have too low Sersic indices. Mergers of spirals with de Vaucouleurs (classical) bulges produce remnants with larger Sersic indices correlated with luminosity, as with SDSS ellipticals. Exponential bulge mergers are better fits to faint ellipticals, whereas classical bulge mergers better match luminous ellipticals. Similarly, luminous ellipticals are better reproduced by remnants undergoing many (>5) mergers, and fainter ellipticals by those with fewer mergers. The remnants follow tight size-luminosity and velocity dispersion-luminosity (Faber-Jackson) relations (<0.12 dex scatter), demonstrating that stochastic merging can produce tight scaling relations if the merging galaxies also follow tight scaling relations. The slopes of the size-luminosity and Faber-Jackson relations are close to observations but slightly shallower in the former case. Both relations' intercepts are offset - remnants are too large but have too low dispersions at fixed luminosity. Some remnants show substantial (v/σ> 0.1) rotational support, although most are slow rotators and few are very fast rotators (v/σ> 0.5). These findings contrast with previous studies concluding that dissipation necessary to produce ellipticals from binary mergers of spirals. Multiple, mostly minor and dry mergers can produce bright ellipticals, whereas significant dissipation could be required to produce faint, rapidly-rotating ellipticals.

preprint2012arXiv

PAndAS in the mist: The stellar and gaseous mass within the halos of M31 and M33

Large scale surveys of the prominent members of the Local Group have provided compelling evidence for the hierarchical formation of massive galaxies, revealing a wealth of substructure that is thought to be the debris from ancient and on-going accretion events. In this paper, we compare two extant surveys of the M31-M33 subgroup of galaxies; the Pan-Andromeda Archaeological Survey (PAndAS) of the stellar structure, and a combination of observations of the HI gaseous content, detected at 21cm. Our key finding is a marked lack of spatial correlation between these two components on all scales, with only a few potential overlaps between stars and gas.The paucity of spatial correlation significantly restricts the analysis of kinematic correlations, although there does appear to the HI kinematically associated with the Giant Stellar Stream where it passes the disk of M31. These results demonstrate that that different processes must significantly influence the dynamical evolution of the stellar and HI components of substructures, such as ram pressure driving gas away from a purely gravitational path. Detailed modelling of the offset between the stellar and gaseous substructure will provide a determination of the properties of the gaseous halo of M31 and M33.

preprint2012arXiv

The Fundamental Plane of Galaxy Group Mergers

We present a series of hundreds of collisionless simulations of galaxy group mergers. These simulations are designed to test whether the properties of elliptical galaxies - including the key fundamental plane scaling relation, morphology and kinematics - can be simultaneously reproduced by dry multiple mergers in galaxy groups. Preliminary results indicate that galaxy group mergers can produce elliptical remnants lying on a tilted fundamental plane, even without a central dissipational component from a starburst. This suggests that multiple mergers in groups are an alternate avenue for the formation of elliptical galaxies which could well dominate for luminous ellipticals.

preprint2010arXiv

The photometric properties of a vast stellar substructure in the outskirts of M33

We have surveyed $\sim40$sq.degrees surrounding M33 with CFHT MegaCam in the g and i filters, as part of the Pan-Andromeda Archaeological Survey. Our observations are deep enough to resolve the top 4mags of the red giant branch population in this galaxy. We have previously shown that the disk of M33 is surrounded by a large, irregular, low-surface brightness substructure. Here, we quantify the stellar populations and structure of this feature using the PAndAS data. We show that the stellar populations of this feature are consistent with an old population with $<[Fe/H]>\sim-1.6$dex and an interquartile range in metallicity of $\sim0.5$dex. We construct a surface brightness map of M33 that traces this feature to $μ_V\simeq33$mags\,arcsec$^{-2}$. At these low surface brightness levels, the structure extends to projected radii of $\sim40$kpc from the center of M33 in both the north-west and south-east quadrants of the galaxy. Overall, the structure has an "S-shaped" appearance that broadly aligns with the orientation of the HI disk warp. We calculate a lower limit to the integrated luminosity of the structure of $-12.7\pm0.5$mags, comparable to a bright dwarf galaxy such as Fornax or AndII and slightly less than $1\$ of the total luminosity of M33. Further, we show that there is tentative evidence for a distortion in the distribution of young stars near the edge of the HI disk that occurs at similar azimuth to the warp in HI. The data also hint at a low-level, extended stellar component at larger radius that may be a M33 halo component. We revisit studies of M33 and its stellar populations in light of these new results, and we discuss possible formation scenarios for the vast stellar structure. Our favored model is that of the tidal disruption of M33 in its orbit around M31.

preprint2009arXiv

PAndAS' cubs: discovery of two new dwarf galaxies in the surroundings of the Andromeda and Triangulum galaxies

We present the discovery of two new dwarf galaxies, Andromeda XXI and Andromeda XXII, located in the surroundings of the Andromeda and Triangulum galaxies (M31 and M33). These discoveries stem from the first year data of the Pan-Andromeda Archaeological Survey (PAndAS), a photometric survey of the M31/M33 group conducted with the Megaprime/MegaCam wide-field camera mounted on the Canada-France-Hawaii Telescope. Both satellites appear as spatial overdensities of stars which, when plotted in a color-magnitude diagram, follow metal-poor, [Fe/H]=-1.8, red giant branches at the distance of M31/M33. Andromeda XXI is a moderately bright dwarf galaxy (M_V=-9.9+/-0.6), albeit with low surface brightness, emphasizing again that many relatively luminous M31 satellites still remain to be discovered. It is also a large satellite, with a half-light radius close to 1 kpc, making it the fourth largest Local Group dwarf spheroidal galaxy after the recently discovered Andromeda XIX, Andromeda II and Sagittarius around the Milky Way, and supports the trend that M31 satellites are larger than their Milky Way counterparts. Andromeda XXII is much fainter (M_V=-6.5+/-0.8) and lies a lot closer in projection to M33 than it does to M31 (42 vs. 224 kpc), suggesting that it could be the first Triangulum satellite to be discovered. Although this is a very exciting possibility in the context of a past interaction of M33 with M31 and the fate of its satellite system, a confirmation will have to await a good distance estimate to confirm its physical proximity to M33. Along with the dwarf galaxies found in previous surveys of the M31 surroundings, these two new satellites bring the number of dwarf spheroidal galaxies in this region to 20.

preprint2009arXiv

The Horizon Run N-body Simulation: Baryon Acoustic Oscillations and Topology of Large Scale Structure of the Universe

In support of the new Sloan III survey, which will measure the baryon oscillation scale using the luminous red galaxies (LRGs), we have run the largest N-body simulation to date using $4120^3 = 69.9$ billion particles, and covering a volume of $(6.592 h^{-1} {\rm Gpc})^3$. This is over 2000 times the volume of the Millennium Run, and corner-to-corner stretches all the way to the horizon of the visible universe. LRG galaxies are selected by finding the most massive gravitationally bound, cold dark matter subhalos, not subject to tidal disruption, a technique that correctly reproduces the 3D topology of the LRG galaxies in the Sloan Survey.We have measured the covariance function, power spectrum, and the 3D topology of the LRG galaxy distribution in our simulation and made 32 mock surveys along the past lightcone to simulate the Sloan III survey. Our large N-body simulation is used to accurately measure the non-linear systematic effects such as gravitational evolution, redshift space distortion, past light cone space gradient, and galaxy biasing, and to calibrate the baryon oscillation scale and the genus topology. For example, we predict from our mock surveys that the baryon acoustic oscillation peak scale can be measured with the cosmic variance-dominated uncertainty of about 5% when the SDSS-III sample is divided into three equal volume shells, or about 2.6% when a thicker shell with $0.4<z<0.6$ is used. We find that one needs to correct the scale for the systematic effects amounting up to 5.2% to use it to constrain the linear theories. And the uncertainty in the amplitude of the genus curve is expected to be about 1% at 15 $h^{-1}$Mpc scale. We are making the simulation and mock surveys publicly available.

preprint2009arXiv

The remnants of galaxy formation from a panoramic survey of the region around M31

In hierarchical cosmological models, galaxies grow in mass through the continual accretion of smaller ones. The tidal disruption of these systems is expected to result in loosely bound stars surrounding the galaxy, at distances that reach $10 - 100$ times the radius of the central disk. The number, luminosity and morphology of the relics of this process provide significant clues to galaxy formation history, but obtaining a comprehensive survey of these components is difficult because of their intrinsic faintness and vast extent. Here we report a panoramic survey of the Andromeda galaxy (M31). We detect stars and coherent structures that are almost certainly remnants of dwarf galaxies destroyed by the tidal field of M31. An improved census of their surviving counterparts implies that three-quarters of M31's satellites brighter than $M_V < -6$ await discovery. The brightest companion, Triangulum (M33), is surrounded by a stellar structure that provides persuasive evidence for a recent encounter with M31. This panorama of galaxy structure directly confirms the basic tenets of the hierarchical galaxy formation model and reveals the shared history of M31 and M33 in the unceasing build-up of galaxies.

preprint2009arXiv

Warps and Bars from the External Tidal Torques of Tumbling Dark Halos

The dark matter halos in $Λ$CDM cosmological simulations are triaxial and highly flattened. In many cases, these triaxial equilibria are also tumbling slowly, typically about their short axes, with periods of order a Hubble time. Halos may therefore exert a slowly changing external torque on spiral galaxies that can affect their dynamical evolution in interesting ways. We examine the effect of the external torques exerted by a tumbling quadrupolar tidal field on the evolution of spiral galaxies using N-body simulations with realistic, disk galaxy models. We measure the amplitude of the external quadrupole moments of dark halos in cosmological simulations and use these to force disk galaxy models in a series of N-body experiments for a range of pattern speeds. We find that the torques are strong enough to induce long lived transient warps in disks similar to those observed in real spirals and also induce the bar instability at later times in some galaxy models that are otherwise stable for long periods of time in isolation. We also observe forced spiral structure near the edge of the disk where normally self-gravity is too weak to be responsible for such structure. This overlooked influence of dark halos may well be responsible for many of the peculiar aspects of disk galaxy dynamics.

preprint2008arXiv

Anatomy of the Bar Instability in Cuspy Dark Matter Halos

We examine the bar instability in galactic models with an exponential disk and a cuspy dark matter (DM) halo with a Navarro-Frenk-White (NFW) cosmological density profile. We construct equilibrium models from a 3-integral composite distribution function that are subject to the bar instability. We generate a sequence of models with a range of mass resolution from 1.8K to 18M particles in the disk and 10K to 100M particles in the halo along with a multi-mass model with an effective resolution of ~10^10 particles. We describe how mass resolution affects the bar instability, including its linear growth phase, the buckling instability, pattern speed decay through the resonant transfer of angular momentum to the DM halo, and the possible destruction of the halo cusp. Our higher resolution simulations show a converging spectrum of discrete resonance interactions between the bar and DM halo orbits. As the pattern speed decays, orbital resonances sweep through most of the DM halo phase space and widely distribute angular momentum among the halo particles. The halo does not develop a flat density core and preserves the cusp, except in the region dominated by gravitational softening. The formation of the bar increases the central stellar density and the DM is compressed adiabatically increasing the halo central density by 1.7X. Overall, the evolution of the bar displays a convergent behavior for halo particle numbers between 1M and 10M particles, when comparing bar growth, pattern speed evolution, the DM halo density profile and a nonlinear analysis of the orbital resonances. Higher resolution simulations clearly illustrate the importance of discrete resonances in transporting the angular momentum from the bar to the halo.

preprint2008arXiv

Bars in Cuspy Dark Halos

We examine the bar instability in models with an exponential disk and a cuspy NFW-like dark matter (DM) halo inspired by cosmological simulations. Bar evolution is studied as a function of numerical resolution in a sequence of models spanning 10K to 100M DM particles - including a multi-mass model with an effective resolution of 10G. The goal is to find convergence in dynamical behaviour. We characterize the bar growth, the buckling instability, pattern speed decay through resonant transfer of angular momentum, and possible destruction of the DM halo cusp. Overall, most characteristics converge in behaviour in detail for halos containing more than 10M particles. Notably, the formation of the bar does not destroy the density cusp in this case. These higher resolution simulations clearly illustrate the importance of discrete resonances in transporting angular momentum from the bar to the halo.

preprint2008arXiv

Dynamical Blueprints for Galaxies

We present an axisymmetric, equilibrium model for late-type galaxies which consists of an exponential disk, a Sersic bulge, and a cuspy dark halo. The model is specified by a phase space distribution function which, in turn, depends on the integrals of motion. Bayesian statistics and the Markov Chain Monte Carlo method are used to tailor the model to satisfy observational data and theoretical constraints. By way of example, we construct a chain of 10^5 models for the Milky Way designed to fit a wide range of photometric and kinematic observations. From this chain, we calculate the probability distribution function of important Galactic parameters such as the Sersic index of the bulge, the disk scale length, and the disk, bulge, and halo masses. We also calculate the probability distribution function of the local dark matter velocity dispersion and density, two quantities of paramount significance for terrestrial dark matter detection experiments. Though the Milky Way models in our chain all satisfy the prescribed observational constraints, they vary considerably in key structural parameters and therefore respond differently to non-axisymmetric perturbations. We simulate the evolution of twenty-five models which have different Toomre Q and Goldreich-Tremaine X parameters. Virtually all of these models form a bar, though some, more quickly than others. The bar pattern speeds are ~ 40 - 50 km/s/kpc at the time when they form and then decrease, presumably due to coupling of the bar with the halo. Since the Galactic bar has a pattern speed ~50 km/s/kpc we conclude that it must have formed recently.

preprint2003arXiv

GOTPM: A Parallel Hybrid Particle-Mesh Treecode

We describe a parallel, cosmological N-body code based on a hybrid scheme using the particle-mesh (PM) and Barnes-Hut (BH) oct-tree algorithm. We call the algorithm GOTPM for Grid-of-Oct-Trees-Particle-Mesh. The code is parallelized using the Message Passing Interface (MPI) library and is optimized to run on Beowulf clusters as well as symmetric multi-processors. The gravitational potential is determined on a mesh using a standard PM method with particle forces determined through interpolation. The softened PM force is corrected for short range interactions using a grid of localized BH trees throughout the entire simulation volume in a completely analogous way to P$^3$M methods. This method makes no assumptions about the local density for short range force corrections and so is consistent with the results of the P$^3$M method in the limit that the treecode opening angle parameter, $θ\to 0$. (abridged)

preprint1997arXiv

Parallel TreeSPH

We describe PTreeSPH, a gravity treecode combined with an SPH hydrodynamics code designed for massively parallel supercomputers having distributed memory. Our computational algorithm is based on the popular TreeSPH code of Hernquist & Katz (1989). PTreeSPH utilizes a domain decomposition procedure and a synchronous hypercube communication paradigm to build self-contained subvolumes of the simulation on each processor at every timestep. Computations then proceed in a manner analogous to a serial code. We use the Message Passing Interface (MPI) communications package, making our code easily portable to a variety of parallel systems. PTreeSPH uses individual smoothing lengths and timesteps, with a communication algorithm designed to minimize exchange of information while still providing all information required to accurately perform SPH computations. We have additionally incorporated cosmology, periodic boundary conditions with forces calculated using a quadrupole Ewald summation method, and radiative cooling and heating from a parameterized ionizing background following Katz, Weinberg & Hernquist (1996). The addition of other physical processes, such as star formation, is straightforward. A cosmological simulation from z=49 to z=2 with 64^3 gas particles and 64^3 dark matter particles requires ~6000 node-hours on a Cray T3D, with a communications overhead of ~10% and is load balanced to a ~90% level. When used on the new Cray T3E, this code will be capable of performing cosmological hydrodynamical simulations down to z=0 with ~2x10^6 particles, or to z=2 with ~10^7 particles, in a reasonable amount of time. Even larger simulations will be practical in situations where the matter is not highly clustered or when periodic boundaries are not required.

preprint1995arXiv

Nearly Self-Consistent Disk-Bulge-Halo Models for Galaxies

We describe methods for setting up self-consistent disk-bulge-halo galaxy models. The bulge and halo distribution functions (\df) are functions of $E$ and $L_z$ only. The halo's flattening and rotation can be specified. The disk \df\ is a function of $E$ and $L_z$ and a third ``integral'', $E_z$, the vertical energy, which is approximately conserved in a warm disk with vertical extent. The models also have finite extent making them suitable for N-body simulation. A simulation of a sample model shows that in practice the models are very close to equilibrium making them ideal for experiments on instabilities in galactic disks. We also present a sequence of models closely resembling the Milky Way mass distribution with 5 exponential scale radii and varying halo mass and radial extent.

preprint1993arXiv

Lowered Evans Models: Analytic Distribution Functions of Oblate Halo Potentials

We describe an analytic distribution function of a finite, oblate stellar system that is useful for the practical modelling of dark halos. The function is determined by lowering Evans's (1993) distribution function of a flattened, cored isothermal system in analogy to the lowering of the singular, isothermal sphere in the definition of the King (1966) models. We derive analytic expressions of the density, maximal streaming velocity and velocity dispersion profiles as a function of the potential. As for the King models, the potential must be calculated numerically. We also present a recipe for generating N-body realizations of this distribution function and examine the stability in three models with dimensionless spins $λ=$0.0, 0.05 and 0.18 using N-body simulations with 50,000 particles. The $λ=0.18$ model is unstable to the formation of a triaxial bar within $\sim 5$ King radii while the other models appear stable. We conclude that the slowly rotating systems are useful for modelling flattened dark halos.