1 / 25

The Baryon Induced Transformation of CDM Halos

The Baryon Induced Transformation of CDM Halos. Mario G. Abadi Universidad Nacional de Córdoba, CONICET Argentina In collaboration with Julio Navarro and the Canadian Computational Cosmology Collaboration (C4), University of Victoria, Canada LENAC Latin-american Workshop

gagan
Download Presentation

The Baryon Induced Transformation of CDM Halos

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. The Baryon Induced Transformation of CDM Halos Mario G. Abadi Universidad Nacional de Córdoba, CONICET Argentina In collaboration with Julio Navarro and the Canadian Computational Cosmology Collaboration (C4), University of Victoria, Canada LENAC Latin-american Workshop October 29 to November 1 Guaruja Brazil

  2. How dark matter halos are transformed by baryons? Profile (Theory vs Simulations) Shape (Sphericity vs Triaxiality) Orientation (Halo vs Disk)

  3. Introduction

  4. Profile: Theory Blumenthal et al. (1986) Contraction of a dark matter halo in response to condensation of baryons in its center. The cooling of gas in the centers of dark matter halos is expected to lead to a more concentrated dark matter distribution. The response of dark matter to the condensation of baryons is usually calculated using the standard model of adiabatic contraction

  5. Adiabatic contraction Spherically simetric dark matter halo with circular orbits Conservation of the adiabatic invariant M(r)r I=Integral of q dp = Integral of p dq If p is the circular velocity, then dq=r dtheta I=Integral of v r dtheta using v^2=GM/r I=Int(GM/r)^(1/2) r dtheta = (GM/r)^(1/2) Int(dtheta)/2pi=(GM/r)^(1/2) I=GMr=r^2 v^2

  6. Profile: Simulations Gnedin et al (2004) High-resolution cosmological simulations which include gas dynamics, radiative cooling, and star formation Particle orbits in the halos are highly eccentric Dissipation of gas indeed increases the density of dark matter and steepens its radial profile in the inner regions of halos compared to the case without cooling Simple modification of the assumed invariant from M(r)r to M(r_av)r, where r and r_av are the current and orbit-averaged particle positions

  7. Profile: Simulations Abadi et al. (2006) The adiabatic contraction model C4 Numerical simulations Simulations of 13 galactic dark matter halos With and without gas High and low resolution One halo also at super high and super low resolution

  8. Profile: Simulations Abadi et al. (2006) The abadiatic contraction model C4 Numerical simulations Simulations of 13 galactic dark matter halos With and without gas High and low resolution One halo also at super high and super low resolution

  9. Dark matter halo: with and without gas

  10. Circular Velocity Profile Circular velocity and density profiles are increased by the presence of the baryons Different models give different dark matter profiles in the inner parts

  11. Contracted dark matter halos Polinomial fit to infer contracted dark matter profiles from non-contracted dark matter + gas profiles

  12. Inverse model • Total mass of the disk Mdisk • Contribution of the exponential disk to the total velocity • Vdisk^2(r)=2 G Mdisk/Rdisk x^2 (I0(x) K0(x)-I1(x) K1(x)) • where x=r/2/Rdisk • Contribution of the (contracted) dark matter halo to the total velocity Vdark^2=Vrot^2-Vdisk^2 • Invert the model in order to obtain the circular velocity of the non-contracted (i.e. without gas) dark matter (only) halo at r=2.2Rdisk • Assuming the shape of the dark matter halo (i.e. a concentration parameter “c” for the NFW fit), you have the non-contracted dark matter halo density profile and its Vvirial and Vmaximum

  13. Application to the Milky Way Mdisk=6.0x10^10 solar masses Rdisk=3.5 kpc Vrot(2.2Rdisk)=220 km/sec Mvir=1.9, 1.0 and 0.3 10^12 solar masses Mvir > Mdisk/f_b=0.4 10^12 solar masses Ours: Vmax=188, Vvir=156 km/s Gnedin: Vmax=155, Vvir=129 km/s Standard: Vmax=107, Vvir= 89 km/s

  14. Application Semianalytic models

  15. Application to other galaxies • Main observables for galaxies: x-band surface brightness profile (photometry) and rotation curve (kinematic) • Obtain x-band scalelength Rdisk of and exponential disk, rotational velocity Vrot(r=2.2Rdisk) and also x-band total luminosity L • Assume a stellar mass-to-light ratio M/L (REM: depends on the color index) • It is possible to obtain the disk total mass Mdisk • Go back 3 slides: “Inverse model” • Compare Vrot vs Vmax

  16. Luminosities and Scalelenghts • Courteau (1996, 1997) 306 galaxies with luminosities and scalelengths in Kent-Thuan-Gunn system r-band. Also luminosities in Johnson B-band from RC3 taken from NED • r=R+0.94 (Jorgensen 1994) • Courteau et al (2000) only 36 galaxies with absolute magnitudes and scalelenghts in "Landolt" (is interchangeable with Johnson) I-band and SDSS colors (g-r) and (g-i) • i=I+0.62 and g=V+0.49 (Courteau et al 2006) • M-Msun=-2.5Log(L/Lsun)

  17. Mass-to-light ratio guesstimation from Bell & de Jong (2001) • Compute R=r-0.354 (Jorgensen 1994), then compute B-R, thenM/LR=aR+bR(B-R) • Compute i=I+0.62, then g=i+(g-i), then V=g-0.49, then (V-I), then M/LI=aI+bI (V-I)

  18. Mass-to-light ratio guesstimation from Bell et al (2003) • M/Li=ai+bi(g-r) and/or M/Li=ai+bi(g-i)

  19. The baryonic Tully-Fisher relation • B-R colors (white dots) • g-i colors (green dots)

  20. The baryonic Tully-Fisher relation • B-R colors (white dots) • g-r colors (green dots)

  21. The baryonic Tully-Fisher relation • B-R colors (white dots) • V-I colors (green dots)

  22. Application to UGC 5794 Mdisk=7.4x10^9 solar masses Rdisk=2.5 kpc Vrot(2.2Rdisk)=128 km/sec Mvir=1.9, 1.8 and 1.6 10^12 solar masses Mvir > Mdisk/f_b=5.4 10^10 solar masses

  23. Vrot vs Vmax: Galaxies “with” disk • Vrot = Vmax corresponds to semianalytic models that simultaneously reproduce both the Tully-Fisher and the luminosity function • There are differences between the 3 models

  24. Vrot vs Vmax: Galaxies “without” disk • Vrot = Vmax corresponds to semianalytic models that simultaneously reproduce both the Tully-Fisher and the luminosity function • There are no differences between the 3 models

  25. Conclusions • A new model for the contraction of dark matter halos • Previous models (probably?) overestimate the amount of contraction (Bower & Benson: reason for small disk in semianalytic models?) • Nice rotation curve for the Milky Way • Agreement with semianalytic models • Pending: application to other galaxies with more or less reliable M/L ratios

More Related