Radial Distributions of Surface Mass Density and Mass-to-Luminosity Ratio in Spiral Galaxies

Similar documents
Rotation curve decomposition for size mass relations of bulge, disk, and dark halo components in spiral galaxies

arxiv:astro-ph/ v1 20 Nov 1996

Distance Measurement of Spiral Arms in the Galactic-Center Direction Using Radial-Velocity Gradients

Estimates of the Enclosed Mass and its Distribution. for several Spiral Galaxies. Abstract

Estimates of the Enclosed Mass and its Distribution. for several Spiral Galaxies. Abstract

Rotation curves of spiral galaxies

Star Formation Law in the Milky Way

Does Low Surface Brightness Mean Low Density? W.J.G. de Blok. Kapteyn Astronomical Institute. P.O. Box AV Groningen. The Netherlands.

Astronomy 330 Lecture 7 24 Sep 2010

arxiv:astro-ph/ v1 14 Nov 2003

Three-Dimensional Distribution of the ISM in the Milky Way Galaxy: IV. 3D Molecular Fraction and Galactic-Scale HI-to-H 2 Transition

Surface Brightness of Spiral Galaxies

Unified Rotation Curve of the Galaxy Decomposition into de Vaucouleurs Bulge, Disk, Dark Halo, and the 9-kpc Rotation Dip

CENTRAL ROTATION CURVES OF SPIRAL GALAXIES Y. SOFUE,1 Y. TUTUI,1 M. HONMA,1 A. TOMITA,2 T. TAKAMIYA,1 J. KODA,1 AND Y. TAKEDA1

Modelling mass distribution of the Milky Way galaxy using Gaia s billionstar

A Newtonian Model for Spiral Galaxy Rotation Curves. Abstract

Spatial distribution of stars in the Milky Way

Lecture Three: Observed Properties of Galaxies, contd.! Hubble Sequence. Environment! Globular Clusters in Milky Way. kpc

The motions of stars in the Galaxy

Contents 1. INTRODUCTION ROTATION CURVES Measurements of Galactic Rotation Rotation curve of the Galaxy

Three-Dimensional Distribution of the ISM in the Milky Way Galaxy: I. The H I Disk

Extra Gal 9: Photometry & Dynamics

The cosmic distance scale

Direct comparisons between VERA, VLBA and EVN astrometry results and an analytic gas dynamics model

The Correlation Between Supermassive Black Hole Mass and the Structure of Ellipticals and Bulges

The Virgo High-Resolution CO Survey II. Rotation Curves and Dynamical Mass Distributions

Normal Galaxies ASTR 2120 Sarazin

arxiv:astro-ph/ v1 1 Jan 2003

THE GALACTIC BULGE AS SEEN BY GAIA

The dynamics of face-on galaxies in MOND

On the Keplerian Rotation Curves of Galaxies

Scaling Relations of late-type galaxies

View of the Galaxy from within. Lecture 12: Galaxies. Comparison to an external disk galaxy. Where do we lie in our Galaxy?

Cours d astrophysique III :

Stellar Dynamics and Structure of Galaxies

ASTRO 310: Galactic & Extragalactic Astronomy Prof. Jeff Kenney

Physics of Galaxies 2016 Exercises with solutions batch I

AS1001:Extra-Galactic Astronomy

ASTR 200 : Lecture 22 Structure of our Galaxy

ASTRO 310: Galac/c & Extragalac/c Astronomy Prof. Jeff Kenney. Class 17 Mar 30, 2016 Starlight Distribu/ons in Disk Galaxies

Laws of Galactic Rotation. Stacy McGaugh Case Western Reserve University

Components of Galaxies Gas The Importance of Gas

Dark Matter. ASTR 333/433 Spring 2018 T R 4:00-5:15pm Sears 552

Scale Length of the Galactic Thin Disk

Gravitation and Dark Matter

If we see a blueshift on one side and a redshift on the other, this is a sign of rotation.

Unified Rotation Curve of the Galaxy Decomposition into de Vaucouleurs Bulge, Disk, Dark Halo, and the 9-kpc Rotation Dip

Thom et al. (2008), ApJ

ASTRON 449: Stellar (Galactic) Dynamics. Fall 2014

The Milky Way Galaxy. Some thoughts. How big is it? What does it look like? How did it end up this way? What is it made up of?

Surface Photometry Quantitative description of galaxy morphology. Hubble Sequence Qualitative description of galaxy morphology

A new mechanism for the formation of PRGs

Analyzing Spiral Galaxies Observed in Near-Infrared

distribution of mass! The rotation curve of the Galaxy ! Stellar relaxation time! Virial theorem! Differential rotation of the stars in the disk

Blue Compact Dwarfs:

arxiv: v1 [astro-ph.ga] 19 Mar 2017

Milky Way S&G Ch 2. Milky Way in near 1 IR H-W Rixhttp://online.kitp.ucsb.edu/online/galarcheo-c15/rix/

Whittle : EXTRAGALACTIC ASTRONOMY 5. SPIRAL GALAXIES

Galaxies. CESAR s Booklet

The Milky Way s rotation curve out to 100 kpc and its constraint on the Galactic mass distribution

Homework 1. Astronomy 202a. Fall 2009

Halo Gas Velocities Using Multi-slit Spectroscopy

Einführung in die Astronomie II

Number of Stars: 100 billion (10 11 ) Mass : 5 x Solar masses. Size of Disk: 100,000 Light Years (30 kpc)

arxiv: v2 [astro-ph.ga] 23 Nov 2017

The dark and visible matter content of low surface brightness disk galaxies Blok, W. J. G. de; McGaugh, S. S.

Correlation between kinematics and Hubble Sequence in disk galaxies?

The Milky Way - Chapter 23

Galaxies: The Nature of Galaxies

The origin of the steep vertical stellar distribution in the Galactic disc

Lecture Three: Observed Properties of Galaxies, contd. Longair, chapter 3 + literature. Monday 18th Feb

A galaxy is a self-gravitating system composed of an interstellar medium, stars, and dark matter.

Using globular clusters to test gravity in the weak acceleration regime

5.1 Circular Velocities and Rotation Curves

Astr 5465 Feb. 5, 2018 Kinematics of Nearby Stars

Is the exponential distribution a good approximation of dusty galactic disks?

Astronomy 114. Lecture 29: Internal Properties of Galaxies. Martin D. Weinberg. UMass/Astronomy Department

Our View of the Milky Way. 23. The Milky Way Galaxy

The origin of lopsidedness in galaxies

Low Surface Brightness Galaxies Erwin de Blok. Encyclopedia of Astronomy & Astrophysics P. Murdin

Active Galactic Nuclei-I. The paradigm

Galaxy Rotation Curve of Triangulum Galaxy (M-33) By Using New Quantum Gravity Theory Ravindra S. Mundase 23/11/2017

Summary: Nuclear burning in stars

EVOLUTION OF DUST EXTINCTION AND SUPERNOVA COSMOLOGY

An Introduction to Galaxies and Cosmology

Stars, Galaxies & the Universe Lecture Outline

Astronomy 330 Lecture Oct 2010

Kinematics of the Ionized Halo of NGC 5775

The Milky Way Galaxy (ch. 23)

H I in Galactic Disks

Galaxy Morphology. - a description of the structure of galaxies

More on Galaxy Classifcation

Density and Mass Distribution of the Spiral Galaxy NGC 3672 تمت دراسة توزيع الكثافة و الكتلة للمجرة الحلزونية NGC 3672

The Milky Way. Mass of the Galaxy, Part 2. Mass of the Galaxy, Part 1. Phys1403 Stars and Galaxies Instructor: Dr. Goderya

Lecture Five: The Milky Way: Structure

This week at Astro 3303

ASTRONOMY 460: PROJECT INTRO - GALACTIC ROTATION CURVE

arxiv: v1 [astro-ph.ga] 1 Feb 2015

Visible Matter. References: Ryden, Introduction to Cosmology - Par. 8.1 Liddle, Introduction to Modern Cosmology - Par. 9.1

Today in Astronomy 142: the Milky Way

Transcription:

Publ. Astron. Soc. Japan (204) 00(0), doi:.93/pasj/xxx000 Radial Distributions of Surface Mass Density and Mass-to-Luminosity Ratio in Spiral Galaxies Yoshiaki SOFUE Institute of Astronomy, The University of Tokyo, Mitaka, Tokyo 8-005, Japan E-mail: sofue@ioa.s.u-tokyo.ac.jp Received ; Accepted Abstract We present radial profiles of the surface mass density (SMD) in spiral galaxies directly calculated using rotation curves (RC) by two approximations of flat-disk (SMD-F) and spherical mass distribution (SMD-S). The SMDs are combined with surface brightness (SB) using photometric data to derive radial variations of the mass-to-luminosity ratio (ML). It is found that ML has generally a central peak or a plateau, and decreases to a local minimum at R 0. 0.2h, where R is the radius and h is the scale radius of optical disk. The ML ratio, then, increases rapidly till 0.5h, and is followed by gradual rise till 2h, remaining at around ML 2 in w band (infrared λ 3.4 µm) and [M L ] in r-band (λ6-7500 A). Beyond this radius, ML steeply increases toward the observed edges at R 5h, attaining values as high as ML 20 in w and 2 [M L ] in r-band, indicative of dominant dark matter. The general properties of the ML distributions will be useful to constrain cosmological formation models of spiral galaxies. The radial profiles of the RC, SMD, and ML are available in pdf/eps figures and machine-readable tables as an archival atlas at URL, http://www.ioa.s.u-tokyo.ac.jp/ sofue/smd208/ and as the supplementary data on PASJ home page. Key words: galaxies: rotation curve galaxies: kinematics and dynamics galaxies: mass distribution galaxies: spiral galaxies: structure Introduction The ultimate purpose to study rotation curves (RC) of spiral galaxies is to derive their mass distributions (e.g., Sofue and Rubin ; Sofue 207). By definitiaon, RC represents the circular velocity of the disk plane, and hence treats only axisymmetric quantities as a function of the radius. Accordingly, the derived mass distribution is such integrated in the rotation axis direction, which is the surface mass density (SMD) as projected on the galactic plane. Volume density may be estimated by dividing SMD by the vertical scale height, which needs modeling of the vertical density distribution. On the other hand, SMD has a great advantage that it can be directly compared with photometric observations, which also yields only projected surface brightness (SB) onto the galactic plane. The purpose of this paper is to compile rotation curves of spiral galaxies as many as possible, and calculate SMDs using the direct method as described below. The obtained SMDs are compared with optical and infrared SB in order to obtain radial distributions of the mass-to-luminosity ratio (ML). Based on the SMD and ML for many galaxies, we discuss statistical properties about the mass distribution in galaxies. There have been two major methods to derive the SMD (Sofue 207 for review). One is the deconvolution method, in which a rotation curve is decomposed into several components such as the bulge, disk and dark halo. Each component is expressed by analytic functions of radius including two free parameters (mass and scale radius). Hence, we determine the six parameters of three different types of non-linear functions by the least-squares fitting in a six dimensional parameter space. Finally the sum of thus fitted functions is vertically integrated to reach SMD. Threfore, this method requires sophysticated treatments of data and considerable computing time, so that the c 204. Astronomical Society of Japan.

2 Publications of the Astronomical Society of Japan, (204), Vol. 00, No. 0 method is seldom to be used for SMD, despite that the method is common in studies of galactic structure. On the other hand, this method has an advantage that it is free from the edge effects as described below. An alternative method is the direct method, in which the SMD is directly calculated using the observed rotation curve without assuming any functional forms (Takamiya and Sofue 0). Instead, two extreme approximations are made either that the galaxy is a flat (thin) disk or a sphere. The flat disk approximation results in under-estimation for spheroidal mass, whereas sphere over-estimates disky mass. Namely, the sphere approximation gives better estimation in the central bulge, while the outer disk is better represented by flat-disk. However, the difference of the two cases is at most a factor of.6, and we consider that the true value exists in between. In sphere approximation, we must be cautious about an edge effect caused by the limited extent of observed RC, by which the SMD near the end radius is significantly underestimated. This occurs due to the cylindrical (vertical) projection of spherical mass of finite radius onto the galactic plane. Similarly, although far more slightly, the flat-disk overestimates the mass near the edge due to artificial gravitational attraction caused by cutting off of the disk mass at a finite radius. Considering that the edge effect is slighter in flat disk and that spiral galaxies are more flat-disk like than sphere, we rely more on the direct method in this paper. although we present the results for both approximations. In the analyses, we use the rotation curves compiled by Sofue (206), which we refer to as S-sample, and optical rotation curves obtained by Courteau (996; 997) referred to as C-sample. The S-sample RCs contain those taken from Sofue et al. (999; Nearby galaxy RC atlas); Sofue et al. (3; Virgo galaxy CO line survey); Márquez et al. (4; Isolated galaxy survey); de Blok et al. (8; THINGS survey); Garrido et al. (5; GHASP survey); Noordermeer et al. (7; Early type spiral survey); Swaters et al. (9: Low-SB galaxy survey); and Martinsson et al.(203; DiskMass survey). Several more RCs for individual galaxies were taken from the literature cited in Sofue (206). All RCs used in this paper are shown in figure. For analyses of the mass-to-luminosity ratio (ML), we employ the following archival photometric data for surface brightness of the galaxies: (a) Optical photometry in r-band (λ6800 A (6-7500 A)) observed by Courteau (996, 997), which overlaps with C-sample. (b) Mid-infraraed photometry in w- band (λ3.4 µm (2.8-3.9µm)) by Pilyugin et al. (208), who used photometric maps from the Wide-field Infrared Survey Explorer (WISE; Wright et al. 209), and will be referred to as P-sample. We abbreviate the quantities as SMD-F for SMD calculated by flat-disk approximation; SMD-S for SMD by spherical; ML- (a) (b) Vrot (km/s) 500 300 0 20 30 40 50 Fig.. (a) Rotation curves from compilation by Sofue (206), and (b) optical RCs by Courteau (996, 997). F for ML that is equal to SMD-F divided by surface brightness (SB); ML-S for ML equal to SMD-S divided by SB. The ML in w-band (λ 3.4µm) is corrected for inclination angle of the galactic disk, and ML in r-band (λ 6800 A) is corrected for inclination as well as the extinction. These will be abbreviated as ML-w and ML-r, respectively. The SMD, SB, and ML are measured in units of [M pc 2 ], [L pc 2 ] and [M /L ], respectively. Atlas of full figures and machine-readable tables of RCs, calculated SMDs, and MLs for all analyzed galaxies are available at URL, http://www.ioa.s.u-tokyo.ac.jp/ sofue/smd208/ and on the PASJ home page as supplementary data. 2 Surface Mass Density In this section we present SMDs calculated for rotation curves of about two hundred galaxies from the catalog of Sofue (206), and about three hundred galaxies from the catalog of Courteau et al (996, 997). The surface mass densities, SMD-F for flatdisk approximation and SMD-S for sphere, are calculated by the following methods (Takamiya and Sofue 0), 2. Spherical approximation By spherical approximation, the SMD is represented by SMD S = Σ S(R) = 2π r dm(r) dr, () r 2 R 2 dr where M(r) = R rv (r)2. The spheroidal component in the central G

Publications of the Astronomical Society of Japan, (204), Vol. 00, No. 0 3 region is well represented by this method. However, in the outer disk regions, this equation under-estimates the mass due to the edge effect near the end of the observed rotation curve. Note that SMD is a cylindrical projection of mass on the disk, but a sphere of finite radius does not cover the region away from the disk plane. 2.2 Flat-disk approximation The SMD of a flat thin disk, Σ D(R), is calculated by solving the Poisson s equation assuming a disk of negligible thickness (Freeman 970; Binney&Tremaine 987), and is given by SMD F = Σ F(R) = π 2 G R R 0 ( ) dv 2 dr x ( ) x K dx + R R ( ) dv 2 dr x ( ) R dx K x x.(2) Here, K is the complete elliptic integral, which becomes very large when x R. Generally, this method gives a more reasonable estimates to the mass of a spiral galaxy. (a) 00 0 0 5 5 20 25 30 0 00 0 2.3 SMD plotted against radius Figure 2 shows calculated SMD-F and SMD-S for the Milky Way and M3 together with RC. The radial variations of SMDs are similar to each other, revealing highly concentrated central bulge, exponentially decreasing disk, and slowly decreasing outskirts representing the dark halo. In the Appendix by figures A. and A.2 we present samples of individual plots. Results for all galaxies are presented as an archival atlas at URL. Figure 3 shows SMD-F and SMD-S calculated for rotation curves of S sample galaxies (Sofue 206). Figures 4 and 5 show SMD-F and SMD-S for P and C sample galaxies, where SMDs are plotted against R as well as against the normalized radius, R/h, where h is the scale radius given in Pilyugin et al. (204) and Courteau (996, 997). Note that the plots for SMD-S show a cut-off behaviors at their end radii because of the edge effect during the integration in equation, whereas SMD-F exhibits more mild behavior by equation 2. We will use SMD-F in the following ML analyses. In the figures, we also present radial distributions of the surface brightness (SB=µ) in w (Pilyugin et al. 204) and r-bands (Courteau 996, 997). Here, SBs are corrected for inclination angle of the galactic disk and dust extinction by applying the method described in the next section. We also plot the same against normalized radius by the scale radius, R/h. The distributions, particularly SB, in the outer disks show more ordered and constant gradient behavior in the plots against R/h. In figure 6 we plot SMD-S against SMD-F for the same galaxies, where the data have been under-sampled by a factor (b) 0 5 5 20 25 30 0 Fig. 2. (a) Radial profiles of SMD-F (thick line), SMD-S (dashed line) and RC (thin line) for the Milky Way, and (b) M3. SMD-F (Msun sqpc) SMD-S (Msun/sqpc) 5 4 3 2 Flat disk 0 0 5 5 20 25 30 35 R kpc 5 4 3 2 Sphere 0 0 5 5 20 25 30 35 R kpc Fig. 3. Surface mass densities for flat-disk (SMD-F) (top) and sphere (SMD- S) (bottom) approximations plotted against radius R for S sample galaxies (Sofue 206).

4 Publications of the Astronomical Society of Japan, (204), Vol. 00, No. 0 4 SMD-F SMD-S SB w-band 4 5 w-band, corrected for inc. 3 2 SMD-S (Msun/sqpc) 3 2 Mu-w (Lsun/sqpc) 4 3 2 0 20 30 R kpc 0 20 30 R kpc 0 0 5 5 20 25 30 35 R kpc 4 4 5 w-band, corrected for inc. 3 2 SMD-S (Msun/sqpc) 3 2 Mu-w (Lsun/sqpc) 4 3 2 0 2 3 4 5 6 0 2 3 4 5 6 0 0 2 3 4 5 6 Fig. 4. [Upper panels] SMD-F, SMD-S, and SB (= µ) in w-band for S/P sample galaxies (3.4µm; Pilyugin et al. 204) corrected for i plotted against R, and [lower panels] same but plotted against R/h (bottom). of 5. The SMD-S is systematically greater than SMD-F by a factor of π/2. This is because of the spherical approximation of the mass in the SMD-S, where the mass located at large distances from the galactic plane less contributes to the radial force to balance the centrifugal force by the galactic rotation. The apparent under-estimation of SMD-S at low densities is due to the edge effect caused by the limited maximum radius of the rotation curve. 3 Mass-to-Luminosity Ratio Observed w-band (3.4 µm) SB (= µ) from Pilyugin et al. (204) in units of the solar luminosity in w band as [L arcsec 2 ] is adopted for the w-band ML analysis. Observed data of the r-band (6-7A) surface brightness from Courteau (996, 997) in unit of [mag arcsec 2 ] were converted to r-magnitude to luminosity by assuming that the visual absolute magnitude of the Sun is V = 4.83 mag, and the color index between the r and v band is V R = 0.452 (Tanriver et al. 206). This yields the absolute magnitude of the r-band absolute magnitude of the Sun to be m (r) = 4.38 mag, which was used to calculate the surface brightness SB (= µ) in r-band. SMD-S (Msun/sqpc) 5 4 3 2 0 0 2 3 4 5 Fig. 6. Direct SMD-S plotted against SMD-F for all analyzed galaxies. The former generally exceeds the latter by a factor of.6. Under estimation by SMD-S due to the edge-effect (see section ) appear in the low-density regions corresponding to outer ends of rotation curves.

Publications of the Astronomical Society of Japan, (204), Vol. 00, No. 0 5 4 SMD-F SMD-S SB r-band 4 4 r-band, corrected for inc., extinc. Flat disk Sphere 3 2 SMD-S (Msun/sqpc) 3 2 Mu-r (Lsun/sqpc) 3 2 0 20 30 0 20 30 0 0 20 30 4 Flat disk 4 Sphere 4 r-band, corrected for inc., extinc. 3 2 SMD-S (Msun/sqpc) 3 2 Mu-r (Lsun/sqpc) 3 2 0 2 3 4 5 6 0 2 3 4 5 6 0 0 2 3 4 5 6 Fig. 5. [Upper panels] SMD-F, SMD-S, and SB (= µ) in r-band (6-7A) for C-sample galaxies (Courteau 996; 997) corrected for i and dust extinction plotted against R, and [lower panels] same but against R/h (bottom). Note the ordered behaviors of SB in the plots against R/h. 3. Correction for inclination and interstellar extinction Correction for the inclination of the disk and the interstellar extinction was applied by (cos i) + ( cos i)e R/h SB = µ = µ obs, (3) ( + e τ )/2 where, µ obs is the observed surface brightness at radius R, h is the scale radius of the disk when it is fitted by an exponential function, and i is the inclination angle. Figure 7 illustrates the absorption and inclination condition. The inclination correction represented by the numerator of equation 3 was so applied that the central region inside h has a spherical brightness profile and the outer disk has an elliptical profile corresponding to the disk s inclination i. Thus, the outer disk brightness was corrected to a face-on value multiplying by cos i, whereas the central region is assumed to be spherical, dominated by the bulge brightness, yielding no or less inclination correction. The extinction correction is approximated by the denominator of equation 3. The optical depth τ along the line of sight through the inclined disk at R on the major axis is approximated by τ = τ 0 e R/h sec i. (4) We assume that the interstellar absorption occurs by dust layer distributed near the galactic plane, which is assumed to be sufficiently thinner than the stellar disk. Hence, light from the further side of the disk plane is absorbed, while the light from stars in front of the disk is not absorbed. The factor τ 0 in equation 4 is the vertical optical depth of the galactic disk in the galactic center. To estimate τ 0 in equation 4, we referred to the extinction law in the solar vicinity of the Milky Way, where the r-band extinction is approximately given by (Cardelli et al. 989; Gordon et al. 3; Predehl and Schmitt 995) A r = 0.8A v 0.8 mag [.79 2 H cm 2 ] (5). Taking the local H density to be 0.8 H cm 3 (Sofue 207b), we obtain A r =. mag kpc. Assuming that the thickness of the gas disk to be pc, the perpendicular optical depth of the galactic disk near the Sun can be approximated by τ = 0.3(n/n 0) = 0.3 e (r R 0)/h MW, (6)

6 Publications of the Astronomical Society of Japan, (204), Vol. 00, No. 0 5 Mu-w (Lsun/sqpc) 4 3 2 (a) 0 2 3 4 4 Mu-r (Lsun/sqpc) 3 2 (b) 0 2 3 4 Fig. 8. (a) SB (= µ) vs SMD-F for S/P sample in w-band, and (b) same for C sample in r-band. Fig. 7. Correction for the inclination and dust absorption. where n and n 0 are the gas density in the disk at galacto-centric distance r and r = R 0 = 8 kpc, respectively, which are assumed to be expressed by an exponential function of the radius, and h MW is the scale radius of the Galaxy. For h = 3.5 kpc, we obtain the optical depth at the Galactic Center to be τ 0 (r) = 0.996. We, then, generalize this value to the galactic centers in the analyzed galaxies, and express their optical depth by equation 4. The optical depth in the mid infrared w-band is assumed to be negligible, τ 0(w) 0. These approximations may yield uncertainty in the resulting SB (= µ) values by a factor of 2, mainly due to the too much simplified extinction correction, that may be variable with radii and depend on galaxy types. Thus, we must remember that the ML ratio in this paper includes such uncertainty in addition to the uncertainty arising from the SMD-F and SMD-S selection, which has systematic uncertainty of a factor for.5. However, these uncertainties are systematic varying from a galaxy to another, but not sensitive to the radius within a single galaxy. Therefore, the ML ratio diagrams presented in this paper may be used to discuss the general property of the radial variation, but may not be used for analysis of their absolute values, which may be uncertain within a factor of 3. 3.2 Plot of SB vs SMD Applying the above inclination and extinction corrections, we obtained radial distributions of SB perpendicular to the galactic planes in unit of [L pc 2 ]. The thus calculated radial profiles of SB for P and C sample galaxies are shown in figures 4 and 5. Figure 8 shows the obtained SB values in w and r-bands plotted against SMD-F for P and C sample galaxies. It is trivially shown that SB increases with SMD in all the galaxies, although the correlation is not necessarily linear.

Publications of the Astronomical Society of Japan, (204), Vol. 00, No. 0 7 3.3 Radial variation of ML The mass-to-luminosity ratio at a radius R is calculated by ML(R) = SMD(R) SB(R), (7) where the value has the dimension of [M /L ] in w or r band. Figure 9 shows radial variations of the ML ratio in NGC 97 in w-band from P sample and UGC 4 in r-band from C-sample galaxies. In figures A.3 and A.4 in Appendix we show more examples of individual ML-F and ML-S for P- and C-sample galaxies. All results are presented as an archive atlas at URL. In figure we plot all obtained ML profiles in w and r- bands plotted against radius R and R/h. In the right panels we enlarge the inner region at R/h 2. In the plots against R/h, the radial distributions of ML are similar to each other among the plotted galaxies. The middle panels show that the global ML increases monotonically from the disk toward the halo. However, as the right panels indicate, ML also increases toward the centers, attaining a central peak. Although the resolutions are not sufficient to detect and discuss the central active regions and massive objects in detail, it seems that the central regions at R/h 0. are dominated by high concentration of mass over the luminosity peaks. Leaving the center, ML decreases to local minimum at R/h 0.. It then increases rapidly till R/h 0.5. Beyond this radius ML remains nearly constant at ML 2 in w and in r-band, and gradually increases within the bright optical disk until R/h 2. Beyond this radius, ML starts to increase rapidly, and reaches values as large as several tens in w to a hundred in r-band in the outermost regions, indicating that the halo is dominated by dark matter. Figure shows averged ML-F of the C and S/P sample galaxies after Gaussian averaging in each bin of R/h with the interval of d(r/h) = 0.. Bars are standard errors of the mean value in each bin, which are smaller for C samples because of the larger number of data points in each bin. Note that the standard deviations are much larger as trivially seen in the original plots in figure. The general characteristics of the ML profiles as described above are equally recognized in these averaged plots, while local structures are much smoothed. 3.4 ML vs SMD Figure 2 shows plots of ML against SMD, and ML against SB. The ML generally decreases both with the SMD and SB. These inverse correlations manifest not only the dependence of the ML with the radius in individual galaxies, but also implies that the higher is the SB and the higher is the SMD, the smaller is the ML, regardless positions in galaxies. 00 ML, Mu(Lsun/sqpc), 0 (a) 0. 0 5 5 20 25 30 00 (b) 0 0 5 5 20 25 30 Fig. 9. (a) ML-F (thick increasing line), ML-S (thick dash), SB (= µ) (thick decreasing line), SMD-F (thin decreasing line) and SMD-S (thin dash) for NGC 97 in w-band from S/P sample, and (b) same for UGC 4 in r-band from C sample. 4 Discussion We calculated the SMD-F and SMD-S using compiled rotation curves for about five hundred galaxies. SMDs were compared with surface brightness SB (= µ) in w- and r-bands to obtain radial variations of ML ratios. The SB in r-band was corrected for inclination angle, and that in r-band was corrected both for inclination and dust extinction. The obtained results are presented as an archival atlas at URL. We summarize the results as follows. Galaxies show similar SMDs to each other, having a sharply concentrated central bulge, exponentially decreasing disk, and largely extended outskirt representing the massive halo. SB profiles are also similar to each other. However, the outskirt of SB decreases faster than SMD, obeying the exponentially decreasing function. The similarities of the radial profiles of SB and ML among galaxies are more pronounced, when they are plotted against

8 Publications of the Astronomical Society of Japan, (204), Vol. 00, No. 0 ML-F in w-band, S/P sample 2 2 ML-w (Msun/Lsun) 0 ML-w (Msun/Lsun) 0 ML-w (Msun/Lsun) 0-0 5 5 20 25 30 35 R kpc - 0 2 3 4 5 6-0.0 0.5.0.5 2.0 2 2 ML-F in r-band, C sample 2 ML-r (Msun/Lsun) ML-r (Msun/Lsun) ML-r (Msun/Lsun) 0 0 0 20 30 0 2 3 4 5 6 0 0.0 0.5.0.5 2.0 Fig.. ML-F against R (left panels) and R/h (middle) for P sample (w-band) (upper panels) and C sample (r-band) galaxies (lower), and the same but enlarged for the inner region at R/h 2 (right panels). 2 ML (Msun/lsun) C sample, r-band ML (Msun/Lsun) C sample, r-band 0 P/S sample, w-band 0 P/S sample, w-band (a) 0 2 3 4 5 (b) 0.0 0.5.0.5 2.0 Fig.. (a) Gaussian-running averaged ML-F in each bin of d(r/h) = 0. plotted against R/h for P/S and C sample galaxies. Bars are standard errors of the mean in each bin. Smaller errors for C sample are due to larger number of data points. (b) Same, but enlarge the central region with samller bins of d(r/h) = 0.02 for r-band, and 0.05 for w-band.

Publications of the Astronomical Society of Japan, (204), Vol. 00, No. 0 9 2 w-band 2 ML-w (Msun/Lsun) 0 ML-w (Msun/Lsun) 0-2 3 4 5-0 2 3 4 5 Mu-w (Lsun/sqpc) 2 r-band 2 ML-r (Msun/Lsun) ML-r (Msun/Lsun) 0 2 3 4 0 0 2 3 4 Mu-r (Lsun/sqpc) Fig. 2. ML-F in w-band and r-band against SMD-F and SB (= µ). normalized radius, R/h. ML profiles are also similar to each other among the galaxies. It shows a central peak or a plateau, indicative of a highly concentrated mass near the nucleus. ML, then, decreases to a local minimum at R/h 0. to 0.2. It, then, increases steeply till R/h 0.5 and is followed by gradual rise till R 2h remaining around ML 2 in w and in r-band. Beyond there, ML steeply increases toward observed edge, attaining large values, ML 20 in w and about a hundred in r-band at R 5h, representing dark matter dominant halos. ML ratios in w-band have generally smaller values than those in r-band, manifesting larger SBs in w-band than in r-band. We stress that the here obtained general characteristics of the radial distributions of the SMD and ML ratio will give a strong constraint on the formation and evolutionary scenarios of galaxies in the expanding universe (e.g., Behroozi et al.203; Miller et al. 204). Particularly, the radial variation of ML will be useful to constrain the evolution models of individual galaxies, giving stronger condition than the global relation among the decomposed masses of bulge, disk and dark halos. Acknowledgments: The author thanks Prof. S. Courteau for the data base of the optical rotation curves and r-band photometry of spiral galaxies, and Prof. L. S. Pilyugin for the w- band photometric data. The data analysis was carried out on the open use data analysis computer system at the Astronomy Data Center, ADC, of the National Astronomical Observatory of Japan (NAOJ). References Amado, P. J., & Byrne, P. B. 996, Irish Astronomical Journal, 23,

Publications of the Astronomical Society of Japan, (204), Vol. 00, No. 0 Behroozi, P. S., Wechsler, R. H., & Conroy, C. 203, ApJ, 770, 57 Binney & Tremaine 987 in Galactic Dynamics. Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 989, ApJ, 345, 245 Courteau S., 997, AJ, 4, 2402 Courteau S., 996, ApJS, 3, 363 de Blok, W. J. G., Walter, F., Brinks, E., et al. 8, AJ, 36, 2648 Freeman K. C., 970, ApJ, 60, 8 Gordon, K. D., Clayton, G. C., Misselt, K. A., Landolt, A. U., & Wolff, M. J. 3, ApJ, 594, 279 Garrido, O., Marcelin, M., Amram, P., et al. 5, MNRAS, 362, 27 Márquez I., Masegosa J., Moles M., Varela J., Bettoni D., Galletta G., 2, A&A, 393, 389 Martinsson, T. P. K., Verheijen, M. A. W., Westfall, K. B., et al. 203, AA, 557, AA3 Miller S. H., Ellis R. S., Newman A. B., Benson A., 204, ApJ, 782, 5 Noordermeer, E., van der Hulst, J. M., Sancisi, R., Swaters, R. S., & van Albada, T. S. 7a, MNRAS, 376, 53 Pilyugin L. S., Grebel E. K., Zinchenko I. A., Kniazev A. Y., 204, AJ, 48, 34 Pilyugin L. S., Grebel E. K., Kniazev A. Y., 204, AJ, 47, 3 Predehl, P., & Schmitt, J. H. M. M. 995, AA, 293, 889 Sofue Y., Koda J., Nakanishi H., Onodera S., 3, PASJ, 55, 59 Sofue, Y., Tutui, Y., Honma, M., et al. 999, ApJ, 523, 36 Sofue, Y. 206, PASJ, 68, 2 Sofue, Y. 207, PASJ, 69, R Sofue, Y. 207b, MNRAS, 468, 4030 Sofue Y., Rubin V.,, ARA&A, 39, 37 Swaters, R. A., Sancisi, R., van Albada, T. S., & van der Hulst, J. M. 9, AA, 493, 87 Takamiya, T., & Sofue, Y. 0, ApJ, 534, 670 Tanriver, M., & Özeren, F. F. 206, 9th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun (CS9), 54 Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 20, AJ, 40, 868-88 Appendix Atlas of RC, SMD and ML In this appendix, we present samples of the calculated SMD and ML based on the compiled RC as described in the text in figures A., A.2, A.3, and A.4. Plots for all analyzed galaxies and machine-readable tables of the results are available as an archival atlas at URL.

Publications of the Astronomical Society of Japan, (204), Vol. 00, No. 0 00 MW NGC 224 NGC 253 IC 342 00 00 00 0 0 0 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 00 NGC 97 NGC 365 NGC 808 NGC 2403 00 00 00 0 0 0 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 00 NGC 598 NGC 660 NGC 89 NGC 68 00 00 00 0 0 0 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 00 NGC 284 NGC 2855 NGC 2903 NGC 303 00 00 00 0 0 0 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 00 NGC 3034 NGC 3079 NGC 398 NGC 3628 00 00 00 0 0 0 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 Fig. A.. Examples of SMD-F/S (thick/dashed lines) and RC (thin lines) from 2 S-sample galaxies (Sofue 206). See the URL for all galaxies

2 Publications of the Astronomical Society of Japan, (204), Vol. 00, No. 0 00 U4 U U 66 U 79 00 00 00 0 0 0 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 00 U 95 U 2 U 229 U 272 00 00 00 0 0 0 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 00 U 299 U 3 U 346 U 533 00 00 00 0 0 0 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 00 U 540 U 562 U 565 U 624 00 00 00 0 0 0 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 00 U 633 U 673 U 732 U 783 00 00 00 0 0 0 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 0 5 5 20 25 30 0 Fig. A.2. Examples of SMD-F/S (thick/dashed lines) and RC (thin lines) from 289 C-sample galaxies (Courteau et al. 996, 997).

Publications of the Astronomical Society of Japan, (204), Vol. 00, No. 0 3 NGC 224 NGC 253 NGC 300 NGC 925 00 ML, Mu(Lsun/sqpc), 0 00 ML, Mu(Lsun/sqpc), 0 00 ML, Mu(Lsun/sqpc), 0 00 ML, Mu(Lsun/sqpc), 0 0. 0 5 5 20 25 30 0. 0 5 5 20 25 30 0. 0 5 5 20 25 30 0. 0 5 5 20 25 30 NGC 68 NGC 97 NGC 365 NGC 2403 00 ML, Mu(Lsun/sqpc), 0 00 ML, Mu(Lsun/sqpc), 0 00 ML, Mu(Lsun/sqpc), 0 00 ML, Mu(Lsun/sqpc), 0 0. 0 5 5 20 25 30 0. 0 5 5 20 25 30 0. 0 5 5 20 25 30 0. 0 5 5 20 25 30 NGC 2903 NGC 284 NGC 303 NGC 398 00 ML, Mu(Lsun/sqpc), 0 00 ML, Mu(Lsun/sqpc), 0 00 ML, Mu(Lsun/sqpc), 0 00 ML, Mu(Lsun/sqpc), 0 0. 0 5 5 20 25 30 0. 0 5 5 20 25 30 0. 0 5 5 20 25 30 0. 0 5 5 20 25 30 NGC 4254 NGC 4303 NGC 432 NGC 4395 00 ML, Mu(Lsun/sqpc), 0 00 ML, Mu(Lsun/sqpc), 0 00 ML, Mu(Lsun/sqpc), 0 00 ML, Mu(Lsun/sqpc), 0 0. 0 5 5 20 25 30 0. 0 5 5 20 25 30 0. 0 5 5 20 25 30 0. 0 5 5 20 25 30 NGC 450 NGC 4535 NGC 463 NGC 4654 00 ML, Mu(Lsun/sqpc), 0 00 ML, Mu(Lsun/sqpc), 0 00 ML, Mu(Lsun/sqpc), 0 00 ML, Mu(Lsun/sqpc), 0 0. 0 5 5 20 25 30 0. 0 5 5 20 25 30 0. 0 5 5 20 25 30 0. 0 5 5 20 25 30 Fig. A.3. Examples of w-band ML-F/S (increasing thick/dashed lines), w-band µ (SB) corrected for inclination and extinction (decreasing thick lines), and SMD-F/S (decreasing thin/dash) from 3 P-sample galaxies (Pilyugin et al. 204. See the URL for all galaxies.

4 Publications of the Astronomical Society of Japan, (204), Vol. 00, No. 0 00 U4 U U 66 U 79 00 00 00 0 0 0 0 0 5 5 20 25 30 0 5 5 20 25 30 0 5 5 20 25 30 0 5 5 20 25 30 00 U 95 U 2 U 229 U 272 00 00 00 0 0 0 0 0 5 5 20 25 30 0 5 5 20 25 30 0 5 5 20 25 30 0 5 5 20 25 30 00 U 299 U 3 U 346 U 533 00 00 00 0 0 0 0 0 5 5 20 25 30 0 5 5 20 25 30 0 5 5 20 25 30 0 5 5 20 25 30 00 U 540 U 562 U 565 U 624 00 00 00 0 0 0 0 0 5 5 20 25 30 0 5 5 20 25 30 0 5 5 20 25 30 0 5 5 20 25 30 00 U 633 U 673 U 732 U 783 00 00 00 0 0 0 0 0 5 5 20 25 30 0 5 5 20 25 30 0 5 5 20 25 30 0 5 5 20 25 30 Fig. A.4. Examples of r-band ML-F/S (increasing thick/dashed lines), r-band µ (SB) corrected for inclination and extinction (decreasing thick lines), and SMD- F/S (decreasing thin/dash) from 289 C-sample galaxies (Courteau s 996, 997). See URL for all galaxies.