Bar properties and photometry of barred low surface brightness disc galaxies
Wesley Peters, Rachel Kuzio de Naray

TL;DR
This study analyzes the properties of bars in low surface brightness galaxies, comparing them to high surface brightness galaxies, and introduces a new method for measuring bar length that accounts for seeing effects.
Contribution
It provides new measurements of bar properties in LSBs, compares them with HSBs, and introduces an improved bar length measurement technique based on azimuthal light profiles.
Findings
Bars in LSBs are shorter and weaker than in HSBs.
Bar pattern speeds in LSBs are comparable to HSBs, mostly 'fast'.
Barred LSBs are slightly brighter and have shallower disc scale lengths than unbarred LSBs.
Abstract
We measure the bar properties (length, strength, and corotation radius) of a sample of barred low surface brightness (LSBs) galaxies and compare to previous results for both LSBs and high surface brightness galaxy (HSBs) samples. In addition, we present new optical - and -band surface brightness profiles, magnitudes, and colours. We find that bars in LSBs are shorter and weaker when compared with those in HSBs galaxies. Based on analysis of four different bar length measures on simulated galaxy images, we find our bar length measure based on azimuthal light profiles is the most accurate out of the four tested here when seeing effects are taken into account. We also find that bars in LSBs have comparable relative bar pattern speeds to those in HSBs, the majority of which we find to be `fast' (i.e. 1.4). In general, we find that our barred LSBs have…
| Galaxy | R.A. | Dec. | Date | Bands | Seeing | P.A. | D | Morph. | |
|---|---|---|---|---|---|---|---|---|---|
| (J2000) | (J2000) | () | (∘) | (∘) | (Mpc) | ||||
| UGC 628 | 01:00:51.9 | +19:28:33 | 7 August 2016 | 1.74, 1.50 | -42.80 | 58.20 | 86 | SBc | |
| LEDA 135682 | 03:39:34.2 | -00:30:43 | 29 October 2017 | 1.06, 0.92 | 42.66 | 32.75 | 104 | SBbc | |
| LEDA 135684 | 03:41:44.5 | -02:00:09 | 29 October 2017 | 1.12, 0.91 | 30.07 | 44.35 | 46 | SBd | |
| LEDA 135693 | 03:53:26.3 | +00:50:32 | 29 October 2017 | 0.88, 0.87 | 37.84 | 22.12 | 154 | SB(r)b | |
| UGC 2925 | 04:01:02.5 | -00:43:03 | 29 October 2017 | 1.00, 0.79 | -36.64 | 31.79 | 57 | SB(r)c | |
| F563-V2 | 08:53:03.8 | +18:26:09 | 24, 28 February 2017 | 1.50, 1.54 | -32.00 | 29 | 60 | Irr | |
| F568-1 | 10:26:06.3 | +22:26:01 | 24, 28 February 2017 | 1.43, 1.08 | -86.00 | 24.90 | 92 | SBbc | |
| F568-3 | 10:27:20.2 | +22:14:24 | 24, 28 February 2017 | 1.59, 1.15 | -11.40 | 39.60 | 83 | SBcd | |
| LEDA 135782 | 11:28:29.6 | +00:08:40 | 13 May 2018 | 1.27, 1.21 | 5.27 | 35.91 | 208 | SBc | |
| UGC 8066 | 12:57:00.3 | +01:01:43 | 13 May 2018 | 1.11, 0.91 | -8.46 | 50.90 | 40 | SB(r)d | |
| LEDA 135867 | 14:48:56.4 | -00:43:38 | 13 May 2018 | 1.32, 1.25 | 73.40 | 38.30 | 119 | SBd | |
| F602-1 | 22:34:46.0 | +22:33:48 | 25 August 2017 | 0.84, 0.97 | -6.00 | 28.58 | 105 | SBc | |
| PGC 70352 | 23:03:21.5 | +01:53:12 | 25 August 2017 | 1.11, 1.10 | 60.90 | 31.09 | 73 | SBd | |
| ASK 25131 | 23:12:21.0 | -01:05:42 | 25 August 2017 | 1.04, 0.83 | -2.96 | 38.00 | 107 | SBcd | |
| ISI96 2329-0204 | 23:31:40.8 | -01:48:01 | 25 August 2017 | 0.87, 0.74 | 20.29 | 31.43 | … | SBcd |
| Galaxy | ||||
|---|---|---|---|---|
| UGC 628 | 11.63 0.68 | 14.36 0.68 | 16.96 0.46 | 11.21 0.92 |
| LEDA 135682 | 2.05 0.2 | 2.73 0.26 | 2.90 0.23 | 2.05 0.46 |
| LEDA 135684 | 2.87 0.27 | 9.01 0.86 | 6.30 0.46 | 7.07 0.92 |
| LEDA 135693 | 7.91 0.76 | 8.70 0.83 | 6.10 0.46 | 5.70 0.46 |
| UGC 2925 | 11.60 1.11 | 12.76 1.22 | 9.10 0.46 | 8.44 0.92 |
| F563-V2 | 5.02 0.91 | 7.96 0.46 | 6.65 0.46 | |
| F568-1 | 8.44 0.91 | 10.26 0.91 | 7.76 0.46 | 4.37 0.46 |
| F568-3 | 7.75 0.68 | 16.64 0.68 | 13.96 0.46 | 8.93 0.92 |
| LEDA 135782 | 3.11 0.3 | 3.77 0.36 | 3.50 0.46 | 2.66 0.68 |
| UGC 8066 | 6.07 0.58 | 15.75 1.50 | 16.50 0.34 | 7.07 0.68 |
| LEDA 135867 | 6.77 0.65 | 8.19 0.78 | 4.39 0.34 | |
| F602-1 | 1.19 0.27 | 2.80 0.27 | 4.20 0.34 | 4.11 0.34 |
| PGC 70352 | 5.28 0.50 | 7.03 0.67 | 5.02 0.68 | |
| ASK 25131 | 3.28 0.31 | 9.36 0.89 | 8.70 0.34 | 7.76 0.92 |
| ISI96 2329-0204 | 1.64 0.16 | 2.64 0.25 | 2.05 0.46 |
| Galaxy | Sb | |||||
|---|---|---|---|---|---|---|
| (arcsec) | (kpc) | (arcsec) | (kpc) | |||
| UGC 628 | 11.21 0.92 | 4.67 0.38 | 13.96 0.46 | 5.81 0.19 | 1.25 0.11 | 0.32 |
| LEDA 135682 | 2.05 0.46 | 1.04 0.23 | 2.70 0.34 | 1.36 0.17 | 1.32 0.34 | 0.07 |
| LEDA 135684 | 7.07 0.92 | 1.57 0.20 | 7.90 0.46 | 1.75 0.10 | 1.12 0.16 | 0.30 |
| LEDA 135693 | 5.70 0.46 | 4.25 0.69 | 7.30 0.46 | 5.45 0.34 | 1.28 0.22 | 0.36 |
| UGC 2925 | 8.44 0.92 | 2.34 0.26 | 12.50 0.46 | 3.47 0.13 | 1.48 0.17 | 0.18 |
| F563-V2 | 6.65 0.46 | 1.95 0.13 | 15.86 0.46 | 4.64 0.13 | 2.38 0.18 | 0.29 |
| F568-1 | 4.37 0.46 | 1.94 0.20 | 5.86 0.46 | 2.61 0.20 | 1.34 0.18 | 0.13 |
| F568-3 | 8.93 0.46 | 3.61 0.37 | 10.06 0.46 | 4.07 0.19 | 1.13 0.13 | 0.25 |
| LEDA 135782 | 2.66 0.68 | 2.69 0.69 | 2.80 0.34 | 2.83 0.34 | 1.05 0.30 | 0.11 |
| UGC 8066 | 7.07 0.68 | 1.39 0.13 | 7.50 0.34 | 1.47 0.07 | 1.06 0.11 | 0.38 |
| LEDA 135867 | 4.39 0.34 | 2.53 0.20 | 5.40 0.34 | 3.11 0.20 | 1.23 0.12 | 0.37 |
| F602-1 | 4.11 0.34 | 2.09 0.17 | 5.30 0.34 | 2.70 0.17 | 1.29 0.14 | 0.21 |
| PGC 70352 | 5.02 0.68 | 1.79 0.24 | 6.60 0.46 | 2.35 0.16 | 1.32 0.20 | 0.13 |
| ASK 25131 | 7.76 0.92 | 4.02 0.48 | 10.60 0.46 | 5.49 0.24 | 1.37 0.17 | 0.49 |
| ISI96 2329-0204 | 2.05 0.46 | 2.80 0.28 | 1.37 0.34 | 0.18 |
| Galaxy | Band | ||||||
| (mag arcsec-2) | (arcsec) | (kpc) | (mag) | (mag) | (mag) | ||
| UGC 628 | B | 22.36 | 13.87 | 5.78 | 15.35 | -19.32 | 0.158 |
| I | 20.15 | 8.33 | 3.47 | 14.25 | -20.42 | 0.065 | |
| LEDA 135682 | B | 22.31 | 10.09 | 5.10 | 15.49 | -19.60 | 0.380 |
| I | 20.97 | 7.27 | 3.67 | 14.86 | -20.23 | 0.158 | |
| LEDA 135684 | B | 22.37 | 10.61 | 2.13 | 15.61 | -17.69 | 0.344 |
| I | 21.58 | 9.58 | 1.40 | 15.05 | -18.25 | 0.143 | |
| LEDA 135693 | B | 21.42 | 6.30 | 4.70 | 15.51 | -20.43 | 1.086 |
| I | 20.17 | 5.72 | 4.27 | 14.47 | -21.47 | 0.451 | |
| UGC 2925 | B | 20.71 | 10.18 | 2.83 | 13.86 | -19.93 | 1.781 |
| I | 18.75 | 8.99 | 2.50 | 12.19 | -21.60 | 0.739 | |
| F563-V2 | B | 21.95 | 7 | 2.05 | 15.87 | -18.18 | |
| I | |||||||
| F568-1 | B | 23.65 | 12.86 | 5.72 | 16.23 | -18.58 | |
| I | 22.32 | 12.17 | 5.41 | 15.02 | -19.79 | ||
| F568-3 | B | 22.79 | 10.76 | 4.35 | 15.91 | -18.70 | |
| I | 21.04 | 8.37 | 3.38 | 14.70 | -19.91 | ||
| LEDA 135782 | B | 22.22 | 4.58 | 4.63 | 17.16 | -19.43 | 0.094 |
| I | 20.41 | 3.35 | 3.38 | 16.01 | -20.58 | 0.039 | |
| UGC 8066 | B | 23.34 | 14.04 | 2.76 | 16.10 | -16.94 | 0.063 |
| I | 21.83 | 11.34 | 2.23 | 15.07 | -17.97 | 0.026 | |
| LEDA 135867 | B | 22.38 | 4.92 | 2.84 | 17.19 | -18.19 | 0.179 |
| I | 21.07 | 4.15 | 2.39 | 16.24 | -19.14 | 0.074 | |
| F602-1 | B | 22.47 | 11.87 | 4.92 | 15.24 | -19.87 | 0.169 |
| I | 20.92 | 9.66 | 4.65 | 14.14 | -20.97 | 0.070 | |
| PGC 70352 | B | 21.58 | 7.64 | 2.72 | 15.33 | -19.00 | 0.261 |
| I | 20.18 | 6.64 | 2.37 | 14.24 | -20.09 | 0.108 | |
| ASK 25131 | B | 23.09 | 9.14 | 4.74 | 16.55 | -18.59 | 0.135 |
| I | 21.16 | 6.49 | 3.36 | 15.36 | -19.78 | 0.056 | |
| ISI96 2329-0204 | B | 23.14 | 9.38 | 16.48 | 0.162 | ||
| I | 21.62 | 9.04 | 15.02 | 0.067 |
| Galaxy | bar | disc | area |
|---|---|---|---|
| UGC 628 | 1.49 | 0.98 | 1.08 |
| LEDA 135682 | 1.61 | 0.54 | 0.55 |
| LEDA 135684 | 0.93 | 0.45 | 0.49 |
| LEDA 135693 | 1.52 | 0.86 | 0.90 |
| UGC 2925 | 1.94 | 1.53 | 1.55 |
| F563-V2 | 1.64 | 1.57 | |
| F568-1 | 1.57 | 1.32 | |
| F568-3 | 1.58 | 1.29 | |
| LEDA 135782 | 1.61 | 1.05 | 1.09 |
| UGC 8066 | 1.52 | 1.04 | 1.06 |
| LEDA 135867 | 1.46 | 0.87 | 0.96 |
| F602-1 | 1.57 | 1.13 | 1.14 |
| PGC 70352 | 1.57 | 1.03 | 1.05 |
| ASK 25131 | 1.53 | 1.20 | 1.24 |
| ISI96 2329-0204 | 1.64 | 1.41 | 1.42 |
| averages | 1.540.20 | 1.010.31 | 1.110.31 |
| Galaxy | ||||||
|---|---|---|---|---|---|---|
| (Jy km/s) | () | () | () | |||
| UGC 628 | 4.37 | 9.88 | 9.81 | 0.50 | 9.51 | 0.77 |
| LEDA 135684 | 7.52 | 9.57 | 8.94 | 0.21 | 8.26 | 0.97 |
| UGC 2925 | 4.66 | 9.56 | 10.28 | 1.01 | 10.28 | 0.21 |
| F563-V2 | 3.57 | 9.49 | 8.75 | 1.02 | 8.76 | 0.88 |
| F568-1 | 1.91 | 9.58 | 9.56 | 0.71 | 9.41 | 0.67 |
| F568-3 | 2.31 | 9.58 | 9.60 | 0.68 | 9.43 | 0.66 |
| LEDA 135782 | 1.45 | 10.17 | 9.87 | 0.50 | 9.57 | 0.85 |
| UGC 8066 | 4.06 | 9.20 | 8.83 | 0.54 | 8.56 | 0.86 |
| F602-1 | 1.94 | 9.70 | 10.03 | 0.54 | 9.76 | 0.55 |
| PGC 70352 | 3.38 | 9.63 | 9.68 | 0.48 | 9.36 | 0.72 |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Bar properties and photometry of barred low surface brightness disc galaxies
Wesley Peters1 and Rachel Kuzio de Naray1
1Department of Physics & Astronomy, Georgia State University, PO Box 5060, Atlanta, GA 30302-5060, USA E-mail: [email protected]: [email protected]
(Accepted XXX. Received YYY; in original form ZZZ)
Abstract
We measure the bar properties (length, strength, and corotation radius) of a sample of barred low surface brightness (LSBs) galaxies and compare to previous results for both LSBs and high surface brightness galaxy (HSBs) samples. In addition, we present new optical B- and I-band surface brightness profiles, magnitudes, and colours. We find that bars in LSBs are shorter and weaker when compared with those in HSBs galaxies. Based on analysis of four different bar length measures on simulated galaxy images, we find our bar length measure based on azimuthal light profiles is the most accurate out of the four tested here when seeing effects are taken into account. We also find that bars in LSBs have comparable relative bar pattern speeds to those in HSBs, the majority of which we find to be ‘fast’ (i.e. 1.4). In general, we find that our barred LSBs have comparable central surface brightnesses and total magnitudes to unbarred LSBs, albeit slightly brighter. Our barred LSBs fall on a shallower disc scale length (h) vs relation than unbarred LSBs. Finally, we find our barred LSBs to be just as gas rich as unbarred LSBs (i.e. ).
keywords:
galaxies: structure — galaxies:photometry
††pubyear: 2018††pagerange: Bar properties and photometry of barred low surface brightness disc galaxies–Bar properties and photometry of barred low surface brightness disc galaxies
1 Introduction
Low surface brightness galaxies (LSBs) are disc galaxies typically defined as having central surface brightness mag arcsec*-2* (e.g. McGaugh et al., 1995b) or mag arcsec*-2* (e.g. Impey & Bothun, 1997). Although LSBs occupy a large portion of the galaxy population, they are often biased against in large surveys due to their faint nature (McGaugh et al., 1995b; O’Neil & Bothun, 2000; Trachternach et al., 2006). LSBs are not exclusively low-mass galaxies, but instead span all mass ranges (McGaugh et al., 1995a; O’Neil & Bothun, 2000). LSBs are also HI rich (i.e. ) when compared with HSBs (de Blok et al., 1996; Pahwa & Saha, 2018), suggesting they may be relatively unevolved systems.
LSBs are vital to our galaxy formation and evolution theories. In particular, LSBs serve as ideal dark matter laboratories, as analysis of their rotation curves have indicated they are dominated by dark matter at nearly all radii (e.g. de Blok et al., 1996; de Blok & McGaugh, 1997; McGaugh et al., 2000; Swaters et al., 2003; Haghi et al., 2018). As such, the dynamics of unbarred LSBs are particularly well understood. However, given the inherent noncircular motions of barred LSBs, these galaxies and their bars are less well studied. Nonetheless, these systems could shed further light on the nature of dark matter.
Due to their gas rich nature (e.g. McGaugh & de Blok, 1997; Cervantes Sodi & Sánchez García, 2017), dark matter domination (e.g. de Blok & McGaugh, 1997; Bothun et al., 1997), and isolated environments (e.g Bothun et al., 1993; Mo et al., 1994; Rosenbaum & Bomans, 2004; Du et al., 2015; Honey et al., 2018), LSBs are expected to be quite stable against bar formation. Indeed, the bar fraction for LSBs has historically been quite low, 8, (Mihos et al., 1997; Honey et al., 2016) when compared with the HSB bar fraction, 60 (e.g. Eskridge et al., 2000; Marinova & Jogee, 2007). However, large-scale LSB studies have begun increasing this bar fraction to 20 (see Cervantes Sodi, 2017; Pahwa & Saha, 2018). Finally, while LSBs are relatively isolated, there is recent work to suggest that galaxy isolation has less of an effect on bar formation than previously thought (Casteels et al., 2013; Lin et al., 2014; Zana et al., 2018).
In a previous paper (Peters & Kuzio de Naray, 2018, hereafter Paper I), we measured the photometric bar properties of four well-studied LSBs: UGC 628, F563-V2, F568-1, and F568-3. We found that the bars in these galaxies had comparable lengths and relative bar pattern speeds to those in HSBs, albeit weaker. Most interestingly, we found that three of the four galaxies were hosts to ‘fast’ bars, contrary to expectations for objects in non-rotating dark matte haloes.
Here, we use the same methods as in Paper I for 11 new targets that are comparably less well studied, i.e. no existing surface brightness profiles, mass modeling, or rotation curves. We have used the same measurement techniques as those in Paper I, with an additional bar length measure used here. Finally, we have created mock galaxy images in order to fully test our various bar length measurements, specifically to quantify how accurately we can determine the other two parameters we are interested in: strength and corotation radius.
We have also obtained broadband - and -band photometry of our full sample, and derived surface brightness profiles, magnitudes, and colours. In addition, we have used published HI fluxes to obtain HI mass estimates, and population synthesis models to obtain stellar mass estimates in order to examine the gas fractions of a subset of our sample.
The organization of this paper is as follows. We detail our sample, data acquisition and reduction in Sec. 2. Our methods and results for the bar properties of our sample are presented in Sec. 3, and we describe and show results of our mock galaxy image procedure in Sec. 4. We outline our methods for measuring photometry and present the results in Sec. 5. We use published HI fluxes and obtain stellar mass estimates for a subset of our sample in Sec. 6 to examine gas fractions. Finally, our discussion and conclusions are in Sec. 7.
2 Sample and Data
In this section we briefly outline our sample selection and observations of barred LSBs. We also detail our data reduction.
2.1 Sample
As discussed in Paper I, we have assembled a sample of 15 barred LSBs drawn from the catalogs of Schombert et al. (1992) and Impey et al. (1996). Not much is known about the 11 new targets in this sample, with few of the galaxies being observed since appearing in Impey et al. (1996). We show the entire observed sample in Table 1. Here, we list the R.A. (col. 1) and Dec. (col. 2) coordinates, date of observation (col. 3), photometric bands observed (col. 4), average seeing for observations (col. 5), derived position angles (col. 6) and inclinations (col. 7).
In col. 8 we list distances (Mpc). Only one of our galaxies, UGC 628, has a redshift independent distance (Tully et al., 2016). We estimated the distances for the remainder of the galaxies by averaging the derived spectroscopic redshift distances on the NASA Extragalactic Database (NED)111The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration., assuming = 73 km s*-1* Mpc*-1*. There are no available spectra or distance measurements for ISI96 2329-0204. For those galaxies with distances less than 100 Mpc, these distance estimates likely have large errors. However, we use the distances here to obtain estimates of physical bar lengths and disc-scale lengths from our arcsec measurements, allowing the physical values to be scaled later on if a redshift independent distance is determined at a later time.
In col. 9 we list the morphologies of the sample, following the same method as for HSBs (i.e. de Vaucouleurs, 1959). As with LSBs in general, our targets are all late- to very late-type galaxies. With the exception of F563-V2, all our galaxies have a spiral structure, some much more pronounced than others.
2.2 Observations and Reduction
To observe our targets, we used the ARCTIC imager on the 3.5-m telescope at Apache Point Observatory222Based on observations obtained with the Apache Point Observatory 3.5-meter telescope, which is owned and operated by the Astrophysical Research Consortium.. ARCTIC has a field of view of 7.57.5. We used ARCTIC in single read out mode with 22 binning, giving a plate scale of 0.228 arcsec pix*-1*. We obtained - and -band images of 13 new galaxies. In addition, we obtained V- and R-band images for a subset of four galaxies.
We observed each galaxy in all bands with 3600 sec. exposures, with 15 dithering between each image. In addition, we used a circular dithering pattern for the I-band images in order to create a master fringe pattern for reduction. We reduced the data using standard packages and routines in IRAF333IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation: CCDPROC, IMCOMBINE, etc.. For photometric calibration, stars near our targets from Landolt (1992) were observed to obtain photometric zero points. These zero points were found to be accurate to 0.2 mags by checking with stars in the field from Pan-STARRS1 (Chambers et al., 2016) and TASS Mark III (Richmond et al., 2000). Sky values were determined by subtracting the average value of six 5050 pixel star-free boxes in each image (see Paper I and Schombert & McGaugh, 2014).
As with Paper I, we deproject our galaxy images for our analysis. To do this, we use the IRAF task ELLIPSE to determine the disc position angle and inclination of each galaxy in each photometric band. We then take the average value between all available bands to be the P.A. and inclination of the galaxy, with the largest deviation from the average assigned as the error. Errors in P.A. and inclination are typically 5*∘* and 2*∘*, respectively. The deprojection is then performed using GEOTRAN to rotate and ‘stretch’ the image. Our derived P.A. and inclinations are listed in Table 1.
Our final, reduced, on-sky B- and I-band images for our new galaxies are presented in Fig. 1. Here, images of the same photometric band have been scaled the same way for easy comparison (except for the I-band image of UGC 2925 which has a larger scaling range), and every image uses asinh scaling to best show the structure in each galaxy. Each image is 400400 pixels, or 1.521.52.
3 Bar Properties
Here we describe the methods used to characterise each of the bar properties for the full sample. In addition, we present and discuss the results in relation to previous LSB and HSB results.
3.1 Methods
We follow the methods outlined in Paper I for determining the bar length, strength, and corotation radius but we also employ a new bar length measure with the aim of improving the elliptical isophotal bar length measure for our sample. For a more detailed description of each of these methods, see Sec. 3 in Paper I. All of our measurements use the deprojected -band galaxy images, since bars are stellar features and are more prominent and well-defined in redder photometric bands. For the sake of example, we will use LEDA 135782 as a walkthrough for each meaasurement. The remaining galaxies are shown in the Appendix figures.
3.1.1 Bar Lengths
Since the bar strength and corotation radius are dependent on the bar radius, it is important that the bar radius be accurately determined. Therefore, as in Paper I, we use multiple methods for measuring the bar length: (1) elliptical isophotes (e.g., Wozniak et al., 1995; Aguerri et al., 1998, 2000a, 2000b, 2009); (2) Fourier analysis of azimuthal light profiles (e.g. Ohta et al., 1990; Aguerri et al., 1998, 2000a); (3) azimuthal behaviour of the bar (Ohta et al., 1990). We briefly describe each method below. We select the best measurement of the bar length by visually comparing the various measures on the deprojected galaxy images (see Paper I). For a few targets, the bars are no longer than 2 (8.77 pixels), so some methods are unable to find a bar length.
Elliptical Isophotes
When fitting the light distribution of a galaxy with elliptical isophotes, a bar can be found by analyzing the behaviour of both the eccentricity and position angles of the ispohotes. For example, a traditional approach is to use the radius of maximum ellipticity as the bar length (e.g. Wozniak et al., 1995).
However, the bar length can be underestimated when using the radius of maximum ellipticity (Michel-Dansac & Wozniak, 2006). In addition, because the bar strength and relative bar pattern speed depend on the bar length, choosing the most accurate measure of the bar length is crucial. Due to this, we have employed an additional bar length measure using the elliptical isophotes that was not used in Paper I. This method defines the bar length to be where the position angle of the isophotes diverges a certain amount away from the position angle of the bar (Wozniak et al., 1995; Aguerri et al., 2009). We therefore have two measures using the elliptical isophotes: (1) , the radius of maximum ellipticity, and (2) , the radius where the P.A. of the isophotes differs by more than 5 degrees from the value at .
The deprojected radial plots of eccentricity and P.A. for LEDA 135782 are shown in Fig. 2. In each panel, the vertical dashed line indicates the bar length from the respective method. Using the radius of maximum eccentricity, we find a bar length of , and the bar length from the discontinuous P.A. is . The errors for both these methods come from the radial spacing between the points. We find that discontinuous P.A. radius method fails for one of our galaxies, F563-V2.
The bar lengths for each galaxy are listed in Columns 1 and 2 in Table 2. Plots similar to Fig. 2 for the remaining galaxies are in shown in the online Appendix in Fig. 1.
Fourier Analysis
The Fourier analysis is done by decomposing our deprojected images into azimuthal light profiles , and then decomposing these via a Fourier transform:
[TABLE]
The Fourier coefficients are
[TABLE]
[TABLE]
and the amplitudes are
[TABLE]
[TABLE]
In order to determine the bar length, we use the bar/interbar () Fourier intensities method from Aguerri et al. (2000a), where the bar region is defined as
[TABLE]
The last radius at which the above is satisfied is taken to be the bar radius. The bar and interbar intensities are defined as
[TABLE]
[TABLE]
While Aguerri et al. (2000a) state Equation 6 describes the bar region better than , both our results from Paper I and the results from Aguerri et al. (2009) show that this method can tend to overestimate the true bar length. Regardless, we use this method to determine how effective (or not) it is when applied to LSBs.
The bar/interbar Fourier intensity plot for LEDA 135782 is shown in Fig. 3, where the dashed horizontal line denotes Equation 6. We only consider data near the bar region as including radii in the disc can bias the method due to spiral arms, hence why the profile ends near 4. Using this method, we find a Fourier bar length of . The error for this method comes from the radial spacing of the azimuthal light profiles.
LEDA 135782 has interesting Fourier amplitudes, shown in Fig. 4. Here we see strong (deep blue) modes in the bar and disc region, which is not surprising for a barred spiral. However, we also see strong (deep red) and (yellow) modes in the disc. When looking at Fig. 1, LEDA 135782 appears to have three spiral arms, consistent with the strong odd modes.
We find that the strength of the spiral arm pattern greatly affects the accuracy of this method. For instance, for those galaxies with a strong dual-arm pattern (see UGC 628 and F568-1 in Paper I and LEDA 135682 in Fig. 1), the bar length is consistently overpredicted. The method finds the bar length more easily when the spiral arms are weaker.
In addition, we find the Fourier method fails to find a bar length for three of our galaxies, most likely due to the spiral arms: LEDA 135867, PGC 70352, and ISI96 2329-0204 (middle center, middle right, and bottom right in Fig. 1). When examining the Fourier amplitudes for these three galaxies, we find even modes that are stronger in the disc than in the bar, preventing the method from accurately determining a bar region.
The bar lengths for each galaxy using the Fourier Analysis method are listed in Column 3 of Table 2. Figures similar to Fig. 3 and Fig. 4 for the remaining galaxies are shown in the online Appendix in Fig. 2 and Fig. 3.
Azimuthal Light Profile
Finally, we use azimuthal light profiles to characterize the azimuthal behaviour of the bar. Following the procedure outlined in Paper I, we determine the azimuthal centroid of the bar and track either where the humps in the azimuthal light profiles no longer trace a constant azimuthal centroid, or where the azimuthal profiles begin to trace the spiral arms at the end of the bar. To do this, we make the following assumptions: (1) the bar is constant in azimuthal angle, and (2) the (two) spiral arms, if present, are separated by 180*∘* azimuthally. While these two assumptions may not always be true (e.g. the bar may have internal structure that causes it to be warped), they allow us to analyze the behaviour of the bar in order to determine the length.
Due to the different morphologies of our sample and sizes of our galaxies on the sky, we varied the azimuthal spacing and starting radius when constructing the azimuthal light profiles, as opposed to leaving all galaxies with the same radial spacing of 0.46 and starting radii of 2 as was the case in Paper I. This allowed bars that were quite short, 2, to be measured.
The -band and -band azimuthal light profiles for LEDA 135782 are shown in Fig. 5. Here, the profiles are plotted every 0.7 for clarity, and the colour scheme indicates the radial position with redder colours at smaller radii. We can see the clear presence of a bar in both bands, as indicated by the black arrows at 90*∘* and 270*∘* in the inner radii, as well as a three arm pattern present in the outer radii (green profiles) in both bands.
We show the azimuthal information for the bar in LEDA 135782 in Fig. 6, deriving the azimuthal positions and errors of the humps by fitting gaussians to the profiles. Here we show the azimuthal positions of the bar ‘humps’ (top left), the azimuthal difference between the two humps (top right), the azimuthal difference from the bar centroid (bottom left), and the ADU intensity of the bar humps (bottom right). The dashed vertical line in the bottom left denotes the bar length, or where the azimuthal difference from the bar centroid has diverged. LEDA 135782 is a prime case of a very small bar (on the sky), which can pose a challenge when using this method. However, it is clear from looking at the bottom left panel that the first three points are relatively constant, best seen in the open yellow points, and show a dramatic change after the dashed line. We find a bar length of 2.66. The errors are determined based on the behaviour of the azimuthal light profiles, as discussed in Paper I. We assign the larger errors due to the small bar in this galaxy.
For LEDA 135682 and ISI96 2329-0204, the bar is quite small on the sky (2) and hard to measure with this method. In order to get as accurate a measurement as possible, we decreased the starting radius of our azimuthal light profiles to 1 for each galaxy, and lowered the radial spacing for LEDA 135682 to 1.5 pixels (0.34) and to 1 pixel (0.228) for ISI96 2329-0204. To ensure we maintained at least one pixel in each azimuthal division, we increased the size of the azimuthal bins to 9*∘* for both galaxies. This is a loss of resolution on these azimuthal light profiles, but the azimuthal locations of the bar remained unchanged.
The bar lengths for each galaxy measured using azimuthal light profiles are listed in Column 4 of Table 2. Figures similar to Fig. 5 and Fig. 6 for the remaining galaxies are shown in the online Appendix in Fig. 4 and Fig. 5 respectively.
3.1.2 Bar Strengths
Instead of using the same bar strength measure as in Paper I (i.e. from Aguerri et al., 2000a), we modify the method to include the higher order even components, as these can be significant (Ohta et al., 1990):
[TABLE]
where are the even modes, and is the bar radius from the azimuthal method. We note that our measure is only a lower limit, as we cannot probe down to , since this would result in azimuthal bins with no pixels when constructing the azimuthal light profiles.
3.1.3 Corotation Radii
We use the method put forth in Puerari & Dottori (1997) to measure the corotation radius () of the bar, the radius where disc orbits are equal to the pattern speed of the bar. This method determines via the intersection of phase profiles of the Fourier transforms of the B- and I-band images. The phase is given by
[TABLE]
The idea behind this method is based on using the two different photometric bands, and , as proxies for two different stellar populations, young and old. We take the first phase intersection after the bar length to be the corotation radius. Corotation should not occur within the bar region (Contopoulos, 1980; Bureau et al., 1999), and multiple intersections after the bar may be representative of the pattern speed of the disc (Puerari & Dottori, 1997).
Again, we use LEDA 135782 as an example for this process. The phase profiles for LEDA 135782 are shown in Fig. 7. Here we see there are multiple phase intersections near the bar region: 2.2, 2.8, and 3.1. Using the azimuthal bar length from Table 2, this means the bar corotation radius is at 2.8, the first intersection after the bar radius. In addition, we find a phase intersection farther out, at 6.8. This is likely a disc corotation radius.
Plots similar to Fig. 7 for the remaining galaxies are shown in the online Appendix in Fig. 6.
3.2 Results
Here, we discuss the three bar properties (length, strength, corotation radius) for the entire sample. We do not discuss each galaxy individually, but we do focus on interesting or problematic galaxies when applicable (for more detailed discussion on these methods and individual results on our first sample, see Paper I). The bar lengths from each of the four techniques for each galaxy are listed in Table 2, and our final derived bar properties are listed in Table 3.
3.2.1 Bar Lengths
The bar lengths measured using each technique for each galaxy are listed in Table 2. We also include the bar length measurements from Paper I, as well as apply the measure to the four galaxies from Paper I. We determine the final bar length for all galaxies by visually plotting the various bar length measures over the deprojected -band images (see Fig. 8, Fig. 7 in the online Appendix, and Paper I). We report these values as in Table 3. In general, we find an average bar length of 2.5 kpc for our sample.
We plot the four different bar length measures over the deprojected -band image of LEDA 135782 in Fig. 8. Here we see that the azimuthal bar length (blue) is the best measure of the bar in this galaxy, as the other three extend into the spiral arms. We find that this holds for the other galaxies in our sample (see Fig. 7 in the online Appendix), suggesting that the azimuthal method is the most accurate of the four used here, consistent with our findings in Paper I.
With this in mind, in Fig. 9 we show the various measures in relation to the azimuthal method. Here, denotes the difference between either (red triangles), (open black squares), or (blue circles) and . We discuss each measure in relation to below.
We find that the two bar length measures based on the behaviour of the elliptical isophotes do not produce consistent results. equally over and under predicts the bar length compared to . This is because the measure can be biased towards longer values due to to a very elliptical bar or the presence of spiral arms, as well as biased towards shorter bars due to any highly elliptical feature in the inner regions.
We find that almost consistently overpredicts the bar length, often quite extremely (i.e. 2). This is not too surprising, given that this measure is dependent on . Since can only be larger than , this method can only obtain an accurate measurement if underpredicts the bar length.
In general, the Fourier method overpredicts the bar length compared to . However, we find better results for this full sample when compared with our results from Paper I. In fact, appears to be a better predictor of bar length than or when looking at Fig. 9 (at least when ).
The azimuthal method does not rely on prior knowledge about the size of the bar, or the bar region. The three previous methods, however, can be biased towards shorter or longer values depending on how they are used. For example, if there is inner structure within the bar, the radius of maximum ellipticity () can be pushed inwards towards shorter values (see LEDA 135684, F602-1 and ASK 25131 in Table 2). This can in turn bias , as this is determined via a change in P.A. relative to the value at . To mitigate this effect, one can begin fitting isophotes at a larger initial radius, but this can lead to missing structure and requires knowledge about the bar length, possibly introducing bias. Our various bar length measures will be explored in more detail in Sec. 4.
In order to obtain final bar lengths, we have visually examined the bar lengths over the deprojected -band images for all galaxies (see Fig. 8, Fig. 25, and Paper I). We have selected the best bar length from this examination and list them as the bar length in Table 3. In all cases, the best bar length measure comes from our azimuthal method.
3.2.2 Bar Strengths
Using Equation 9 and as , we find bar strengths that range from 0.07 (LEDA 135682) to 0.49 (ASK 25131), with an average value of 0.25. As with Paper I, we examined the Fourier amplitudes of each galaxy and decreased the starting location of the azimuthal light profiles in order to gauge the confidence of our bar strength measure. We found that this did not significantly change the values, suggesting that our bar strengths, while lower limits, are accurate. We find that accounting for the higher order even modes gives a better picture of the bar strength, as some galaxies have significant contribution from these modes. For example, we find a bar strength of 0.25 for F568-3 here, and a bar strength of 0.19 in Paper I where we only accounted for the mode.
3.2.3 Corotation Radii
We find that the majority of the bar corotation radii in our sample occur very close to the end of the bars in our sample. This is consistent with results from HSBs (see e.g. Pérez et al., 2012; Aguerri et al., 2015; Sierra et al., 2015), and from the LSBs in Paper I. We address a problem galaxy, ISI96 2329-0204, below.
Special Note on ISI96 2329-0204
Measuring a corotation radius for ISI96 2329-0204 proves to be quite difficult. For this galaxy, we have BVRI images, and thus create phase profiles for all four bands (see Fig 24 in the Appendix). When examining the phase intersections, we find that the B-band behaves quite differently from the other three bands. So much so that it does not intersect the I-band phase profile until 11. This would leave us with a relative bar pattern speed that is 10, not a realistic number as this is more likely indicating the location of a disc corotation. This deviation from the other three bands is most pronounced at radii greater than 3, near the end of the bar.
Because this bar is quite small on the sky, we decreased the radial spacing of the azimuthal light profiles to one pixel, or 0.228. This did not produce better results, and no phase intersection between the B- and I-bands are observed at radii less than 11.
To add to this conundrum, the B-band phase profile intersects the R-band phase profile at 2.8 and 3, both at the end of the bar. Whether or not this is the radius of corotation is not entirely clear, as the method from Puerari & Dottori (1997) needs photometric bands separated by a large enough wavelength range. However, there has been previous work that uses ‘bluer’ bands than I- to determine the corotation radius (see Sierra et al., 2015). With this in mind, we will take the phase intersection between the B- and R-bands at 2.8 to be the corotation radius of the bar, as using the phase intersections between B- and I- is most likely not due to the bar.
3.3 Discussion of Bar Properties
Our final bar properties for the whole sample are listed in Table 3. Our sample has an average bar length of 2.5 kpc. This is shorter than what is typical for HSBs. When looking at surveys of HSBs, bars in these galaxies tend to be in the range kpc 5 kpc (see e.g. Erwin, 2005; Marinova & Jogee, 2007; Aguerri et al., 1998, 2009).
With the corotation radii of our sample, we can calculate the relative bar pattern speeds for our sample, , also listed in Table 3. Bars are considered ‘fast’ rotators if and ‘slow’ rotators if (Athanassoula, 1992; Elmegreen et al., 1996; Debattista & Sellwood, 2000). While historically the almost unanimous result for HSBs are fast bars (e.g. Pérez et al., 2012; Aguerri et al., 2015; Sierra et al., 2015), recent work has found an increasing number of ‘slow’ bars (Font et al., 2017; Guo et al., 2019). The results for dark matter dominated galaxies, LSBs specifically, are much more unclear due to the lack of pattern speed measurements.
Centrally dense, nonrotating, dark matter halos are expected to dynamically slow down bars over time (Weinberg, 1985; Debattista & Sellwood, 2000). Indeed, very slow bars in dark matter dominated galaxies have been reported in the literature (Bureau et al., 1999; Chemin & Hernandez, 2009; Banerjee et al., 2013). In order to examine possible trends between and galaxy properties, Guo et al. (2019) used a sample of 53 barred galaxies and found a significant number of very slow bars. Interestingly, they found no correlation between a very slow bar and dark matter content. Our previous results in Paper I indicated that three out of the four galaxies we analyzed were hosts to fast bars. Here we find that 13 out of the 15 galaxies in this work are hosts to fast bars.
Noteably, only one galaxy in our sample has a very slow bar (i.e. ), F563-V2, and one galaxy, UGC 2925, has a slow bar that is comparable to previous results on HSBs (Aguerri et al., 1998; Sierra et al., 2015; Font et al., 2017; Guo et al., 2019). The remainder of the bars in our sample are fast. While this may seem to indicate that the bars in our sample have not been slowed down, dark matter halos are not static and can have significant angular momentum, arising from tidal exchanges during the formation of the galaxy (see Peebles, 1969). This in turn means that the angular momentum exchange between the disc and the halo can be mitigated, sometimes significantly (Long et al., 2014). Thus, halos cannot serve as a pure angular momentum sink, and the angular momentum of the halo, often characterised by the halo spin parameter , can dictate important aspects of a bar.
For example, Fujii et al. (2019) found that for high halos, bars actually speed up, while also becoming shorter and weaker. This is especially of note, as LSB discs are expected to form in high halos (e.g. Dalcanton et al., 1997; Jimenez et al., 1998; Kim & Lee, 2013). Thus, without knowing , we cannot truly make a definitive statement on the implications of our relative bar pattern speeds.
We show versus bar strength in Fig. 10, including the HSBs from Aguerri et al. (1998), as this work uses the same measure as used here for corotation radius and bar strength. This is an updated plot of Fig. 14 from Paper I. Here we also show the fit to both LSBs and HSBs, excluding the outlier of F563-V2 (the point far above the rest) and find with a scatter of 0.13. This is close to identical to the relation from Paper I, with a near identical scatter, and is consistent with the lack of relation found by Guo et al. (2019). Our relation is plotted in Fig. 10 as the solid line, with the shaded region indicating the scatter.
When looking at Fig. 10, we find that LSBs and HSBs form a continuum, with the LSBs clustered on the weaker strength end (average ) transitioning into the stronger HSBs (0.37 for the HSBs shown here). This also then translates to surface brightness, with bars being fast rotators across all surface brightnesses.
As both the strength and relative bar pattern speed are dependent on the bar length, it is crucial that this property be measured as accurately as possible. For this reason, we have employed the four different measurement techniques here and have found that our method performs the best for all the galaxies in our sample. However, in order to be as confident in this method as possible, we have tested the four methods used here on simulated galaxy images consisting of a disc and bar. This is therefore an idealized case for each method, as no spiral arms or disc features can affect the measurement. We spend the entirety of the next section discussing this process. Because the azimuthal bar length performs so well on our data, we expect this method to perform just as well on the simulated images. However, the lack of spiral arms or disc features could affect the accuracy of the other three methods as well.
4 Analysis of Mock Galaxy Images
In order to determine how well we can determine the bar properties of our galaxies, we constructed fake galaxy images with randomised structural and photometric parameters. We used observational parameters similar to those of our observed sample and those of the ARCTIC imager in order to compare with our observations. We broadly follow the procedure in Aguerri et al. (2009).
4.1 Creating the Images
For each galaxy, we assume two components: an exponential disc and a Ferrers bar (see Laurikainen et al., 2005). We do not consider a bulge component since LSBs do not typically have a bulge (Pahwa & Saha, 2018), and because our sample lacks significant bulges. In addition, we do not include spiral arms in our images in order to provide an ideal scenario for each bar length fitting method. The intensity profile of an exponential disc is given by
[TABLE]
where is the central disc intensity and is the disc scale length. The intensity profile of a Ferrers bar is given by
[TABLE]
where is the central bar intensity, is the length of the bar, and is the bar shape parameter. The intensity is defined to be zero for . The radial coordinate is given by
[TABLE]
where is the ellipticity of the bar, and is a parameter that controls the shape of the bar isophotes (Athanassoula et al., 1990). Pure elliptical isophotes have , boxy isophotes have , and discy isophotes have .
In order to construct observationally accurate galaxies for our purposes, we selected galaxy properties that match the general LSB population, with bar values similar to our results from Paper I and from Sec. 3.2. We begin by randomly assigning a disc central surface brightness
[TABLE]
and bar central surface brightness
[TABLE]
with the condition that the central bar surface brightness may not exceed the disc central surface brightness by more than 1 magnitude. This condition was chosen in order to closely match the photometric properties of our sample (see Sec. 5). We then convert these to intensity values via
[TABLE]
where is the exposure time, is the plate scale of the detector, and is the photometric zeropoint. We used values that matched our observations to obtain images similar to our observed galaxies: sec., arcsec pix*-1*, and . In order to produce images similar to our observations, we select distances
[TABLE]
For the other parameters, such as bar ellipticity and disc inclination, we broadly follow the methods laid out in Aguerri et al. (2009). For the bar length, we select from a range representative of our results in Sec. 3.2. We also force the bar length to be shorter than the disc scale length in order to prevent constructing very unrealistic galaxies (i.e. a bar length of 15 and a disc scale length of 3 is not considered realistic).
After the galaxy images were created, we added noise to each pixel in order to closely match our observations. The gain and read-out noise of ARCTIC are 1.98 e*-* ADU*-1* and 3.8 e*-1* respectively. Sky noise was added in assuming a 600 sec. exposure in B. Finally, to account for seeing, we convolved the images with a 2D Moffat PSF
[TABLE]
where is the seeing parameter, and is typically taken to be 3.5 (see Trujillo et al., 2001). To test a wide range of observations, we randomly select our seeing to be between 0.8 and 2, consistent with the average seeing values of our observations (Table 1).
Examples of the images we created are shown in Fig. 11. These images were randomly selected, and are representative of the types of galaxy images we have created. In total we have created 200 images. Histograms of bar lengths, ellipticity, and seeing are shown in Fig. 12. Here we can see that we are sampling a large range of bar lengths, but that we have a large over sampling of shorter bar lengths. This is because the bar lengths in our sample are quite short (Sec. 3.2) and we are hoping to gauge how well the various methods work for our sample. We see that we are sampling a large range of bar ellipticities, . The two peaks in the seeing histogram arise from the distributions around ‘good’ and ‘bad’ seeing, which we have defined here to be 1 and 1.8 respectively.
4.2 Measurements of Mock Images
To maintain consistency, we follow the same reduction process as our real observations. We determine the sky value by using the box method of Schombert & McGaugh (2014) and use ELLIPSE and GEOTRAN to deproject our images. We then use our four methods described in Sec. 3.1 to determine the bar length.
As opposed to our real galaxies (see Sec. 3.2), we do not have concerns over spiral arms influencing our isophotal bar length measures, as we have explicitly not included them. Therefore, we have an ideal scenario for each measurement technique.
Due to the large number of images tested (200), we automated the azimuthal light profile measurement by standardizing each set of profiles. We did this by finding the azimuthal centroid of the bar, usually located at the maximum value of the profile. We then adjusted the entire set of profiles so that the bar centroid is aligned with either 90*∘* or 270*∘, depending on which is closest to the centroid. It is then assumed that the two humps in the azimuthal light profiles will remain at 90∘* and 270*∘*. As there are no spiral arms in our images, we instead use the relative intensity of the bar hump to the disc to determine where the bar ends.
4.3 Results of Measurements
We show the comparison between the true and measured bar length for our mock images in Fig. 13: (top left), (top right), (bottom left), and (bottom right). Blue points denote those images with less than 1 seeing, and red points denote images with greater than 1 seeing. In each of these subfigures the solid line denotes unity, and the dashed lines are fits to the blue points (i.e. those with good seeing). It is clear from Fig. 13 that the seeing in each image greatly increases the scatter of the measurement.
We find that the radius of maximum ellipticity (top left panel) does a poor job of measuring the bar length, often measuring a length 50 shorter than the true length. This is an expected result, especially for Ferrers bars (Michel-Dansac & Wozniak, 2006; Aguerri et al., 2009). The fit for this measure is . Even though (top right panel) depends on , we find that it does an excellent job of finding the bar length, only really failing at larger radii (i.e. ). The fit for this measure is . Our bar length measure using the azimuthal light profiles (bottom left), , performs the best out of the four methods tested, albeit only slightly better than . The fit for this measure is . Lastly, we find that under predicts the bar length, contrary to our results using real data (see UGC 628 and F568-1 in Paper I and LEDA 135682 in Fig. 1, for example) and previous works (Aguerri et al., 2009), suggesting that spiral arms can heavily bias this method to longer bar lengths. The fit for this measure is .
Based on this modeling, we find that the bar length measure based on the behaviour of the azimuthal light profiles to be the best measure out of the four used here. While almost performs as well, the dependence on for this measure is concerning, and can result in wrong measurements. That performs so much better here than with our real data is surprising, although this could simply be that underpredicts the bar length so significantly here. In addition, the change in P.A. used here is simply following what has been used before (i.e. from Aguerri et al., 2009). Contrary to this, there is strong motivation for the azimuthal method, as it traces the location of the bar and any spiral arms in the galaxy, only being limited by the resolution and seeing. We do, however, note the large amount of scatter at smaller bar lengths for . In additon, we see that this technique fails to find an accurate bar length for three bars, seen as the wildly innacurate points in the bottom left panel of Fig. 13. When examining these images, we found that these failures were caused by a combination of factors. First, all three images contain bars and discs that have central surface brightnesses within 0.4 mag arcsec*-2* of each other, making it somewhat hard to distinguish the two components in the azimuthal light profiles. Second, the position angles of the bars and discs are all within 5*∘* of each other, which causes the bar to become slightly lost in the inclined disc. Finally, two of the failures have poor seeing, which further amplifies the issues caused by the first two points.
This suggests that our bar properties reported in Sec. 3.2 are accurate. More importantly, this means that the bars in our LSBs are indeed fast (but see Sec. 7).
5 Photometry
In this section, we detail our methods for measuring the photometry, and present our results. The surface brightness profiles, magnitudes, and colours for the 11 targets not in Paper I are new measurements: as to the best of our knowledge no such measurements have been previously published.
5.1 Methods
Here we detail how we construct surface brightness profiles, as well as obtaining total magnitudes and colours.
5.1.1 Surface Brightness Profiles
We first obtain intensity profiles () of our galaxies using the IRAF routine ELLIPSE. This method of constructing surface brightness profiles for LSBs has been used succesfully in the literature (e.g. Wittmann et al., 2017). We then convert these to surface brightness values via:
[TABLE]
where is the exposure time in seconds, is the plate scale of the instrument, and is the photometric zeropoint of the observation (Sec. 2.2). We then fit for the central surface brightness () and disc scale length (), starting beyond the bar region, (see Sec. 3.2).
In order to correct for the likely low internal dust extinction in our galaxies (e.g. Bothun et al., 1997; Matthews & Wood, 2001; Wyder et al., 2009; Hinz et al., 2007; Honey et al., 2016), we assume the discs are optically thin slabs. The correction for this assumption is simply:
[TABLE]
where is the observed central surface brightness, is the face-on central surface brightness, and a and b are the semi-major and -minor axes respectively (Table 1). This correction increases with the inclination, making magnitudes brighter for highly inclined discs.
5.1.2 Magnitudes
Due to the faint nature of LSBs, using a typical magnitude is not an adequate measure of the total luminosity. Instead, we follow de Blok et al. (1995) and use the total magnitude, , of the disc, found via
[TABLE]
where is the inclination of the disc (Table 1). Integrating out to infinity involves extrapolating the disc scale length, which can introduce error. To counter this, de Blok et al. (1995) also use another magnitude measure, , which takes the entire data of the galaxy disc into account. As their data probe down to 28 mag arcsec*-2* in , almost all the light of the galaxy is being probed. Our observations are not quite as deep, down to 26 mag arcsec*-2* in for most targets, so we do not use here.
5.1.3 Colours
Colour profiles were obtained by subtracting the - and -band surface brightness profiles from each other.
We obtain three different total colours for each galaxy:
a bar colour consisting of the average colour within the bar region () 2. 2.
a disc colour consisting of the average colour outside the bar region () 3. 3.
an area-weighted colour within the 25.5 -mag arcsec*-2* isophote (e.g. McGaugh & Bothun, 1994; de Blok et al., 1995) ()
We measure these three colours in the same manner. First, we create a () colour map for each galaxy by rebinning each image by a factor of 2 to reduce the noise, and subtracting the new -band image from the -band image, obtaining a colour in each pixel (see de Blok et al., 1995; Schombert et al., 2011). We finally obtain colours by taking the mean value in each of the defined regions, rejecting divergent pixels and those with large errors.
5.2 Results
We list central surface brightnesses, magnitudes, disc scale lengths and other information for our whole sample in Table 4. All surface brightnesses and magnitudes reported are corrected for Galactic extinction using values from Schlafly & Finkbeiner (2011), also listed in Table 4. Data for F563-V2, F568-1, and F568-3 are taken from McGaugh & Bothun (1994) and de Blok et al. (1995). Although LSBs are thought to have low dust content, we correct for any internal dust by assuming the discs are optically thin slabs (see Sec. 5.1.1).
With the exception of UGC 628, F563-V2, F568-1, and F568-3, the surface brightness profiles, magnitudes, and colours are all new. In addition, this is one of the few samples of exclusively barred LSB photometry.
5.2.1 Surface Brightness Profiles
Our B- and I-band surface brightness profiles and () colour profiles are shown in Fig. 14. We do not have photometrically calibrated ARCTIC data for F563-V2, F568-1, or F568-3. The majority of our galaxies exhibit very nearly Freeman Type I profiles (Freeman, 1970), or pure exponential disks. For some galaxies, the bar is not noticeable in these surface brightness profiles, LEDA 135782 or PGC 70352 for example. For others, it very clearly dominates the inner light profile, as in the case of LEDA 135693 or ISI96 2329-0204.
The mean -band and -band central surface brightnesses are 22.39 and 20.28 mag arcsec*-2* respectively, consistent with the findings in Zhong et al. (2008), who looked at roughly 12000 LSBs from the SDSS catalog. We find a mean B-band disc scale length of 3.68 kpc and a mean I-band disc scale length of 2.03 kpc. We find no correlation between scale length and central surface brightness, consistent with previous works (McGaugh & Bothun, 1994; de Blok et al., 1995; Zhong et al., 2008; Pahwa & Saha, 2018), as shown in the top panel of Fig. 15. Here, black points are our data, blue squares are de Blok et al. (1995), and red triangles are McGaugh & Bothun (1994), both of which use comparable methods for both barred and unbarred LSBs. We can see that roughly one dex in scale length spans roughly 4 dex in central surface brightness.
5.2.2 Magnitudes
The mean absolute -band and -band magnitudes are and respectively. We find a relation between absolute magnitude and disc scale length, shown as the longer dashed line in the bottom panel of Fig. 15. The relation is given by , with a scatter of (shown as the shaded gray region in Fig. 15). This is slightly shallower than the relation for the general LSB population in Zhong et al. (2008). We also fit for the unbarred LSBs in Fig. 15 from McGaugh & Bothun (1994) and de Blok et al. (1995) and found the relation , shown as the shorter dashed line, indicating a significantly different slope between the two LSB populations.
In the Fig. 16 we show a comparison of the B-band central surface brightness () and absolute magnitude () for our barred sample (black circles) with the unbarred LSBs in McGaugh & Bothun (1994) (red triangles) and de Blok et al. (1995) (blue squares). We find that our barred LSBs are noticeably brighter than their unbarred counterparts, extending off of the trend of the unbarred LSBs.
5.2.3 Colours
Here we discuss our radial colour profiles, three different total colour measures, and our colour maps.
Radial Profiles
We show the radial () profiles in the bottom panels of Fig. 14. We can see that our galaxies are rather blue, with () values of roughly 1 in the disc region outside the bar. In addition, we find that most colour profiles are more red and constant within the bar region, consistent with a primarily stellar feature. All show a clear trend towards bluer values with increasing radius, as expected for LSBs (de Blok et al., 1995).
Total Colours
Our three total colour measures (see Sec 5.1.3) are shown in Table 5. Color data for F563-V2 are taken from McGaugh & Bothun (1994) and data for F568-1 and F568-3 are taken from de Blok et al. (1995) and McGaugh & de Blok (1997) respectively. Average colours are also listed at the bottom of the table.
We find our bar colours (average value of 1.540.20) are comparable to the nuclear colours found in McGaugh & Bothun (1994) (average value of 1.52) and de Blok et al. (1995) (average value of 1.47). In these works, the nuclear region was defined to be the color within a 5 aperature for all galaxies. Here, our bar color is the color within the bar region for each galaxy, which is dependent on each individual galaxy.
Our disc colours are noticeably bluer than the bar colours, with an average value of 1.010.31. We find that our area colour is often very close to the value of the disc colour, with an average value of 1.110.31. Since the bars in our sample are quite small, it is not surprising that the area colours are heavily weighted towards the larger area of the bluer disc.
Colour Maps
We show the () maps of the galaxies in our sample with calibrated photometry used to determine the total colours in Sec. 5.2.3 in Fig. 17. Here, we use a colour bar ranging (with the exception of UGC 2925). Spiral arms are present in these maps as the slightly bluer (whiter) band in some of our maps (LEDA 135693 and F602-1, for example). In addition, a few galaxies show HII regions as white blobs in spiral arms (PGC 70352 and ISI962329-0204, for example). The bars in our galaxies are clear as the redder (darker) regions in the centres of the maps.
All of our galaxies show a stark contrast between the bar and disc regions, with an almost immediate shift from redder colours to much bluer, consistent with our radial colour plots and total colours. For those galaxies with more tenuous disc structure (LEDA 135867 or ASK 25131, for example), the colour maps are not as clear or defined, quickly getting lost in the noise.
Interestingly, some galaxies show bluer regions within the bar: noteably F602-1, ASK 25131, and PGC 70352. This is not a result of the images being misaligned when combining, but instead a real feature. When examining the - and -band images individually (Fig. 1), it is noticeable that the bar in the -band is noticeably different than in the -band for these galaxies. Most noticeable is F602-1, which has a very narrow bar in compared with .
6 Gas Fraction
We have shown that bars in LSBs are characteristically shorter and weaker than those in HSBs (Sec. 3.3). It is thought that high gas content in galaxies prohibits bar formation, as well as forms shorter and weaker bars (e.g. Mayer & Wadsley, 2004; Cervantes Sodi & Sánchez García, 2017). Therefore, we wish to examine the gas fractions of the galaxies in our sample in order to see if this explains the short and weak nature of the bars found.
To obtain gas fractions of our sample, we require estimates of both the gas and stellar mass. We use the HyperLeda online data base444http://leda.univ-lyon1.fr/ (Makarov et al., 2014) to obtain 21cm HI fluxes (). Out of our whole sample, ten galaxies have published HI magnitudes, or , which we convert into fluxes. We then convert these fluxes to HI masses via the relation from Haynes & Giovanelli (1984):
[TABLE]
where D is the distance in Mpc (Table 1). Although LSBs are thought to generally contain very little molecular gas (Mihos et al., 1999), we obtain the gas mass by .
With our -band magnitudes and () colours it is possible to obtain a stellar mass via a mass-to-light ratio, which requires using stellar population models. We first use our -band absolute magnitudes (Table 4) to determine the -band luminosity for each galaxy. We then use the () area colour (Table 5) to obtain an -band stellar mass-to-light ratio using the appropriate equation from Into & Portinari (2013):
[TABLE]
Because we do not have an -band absolute magnitude for F563-V2, we take the apparent magnitude from HyperLeda, and use the distance reported in Table 1. Finally, we obtain gas fractions via
[TABLE]
In Table 6 we report the HI fluxes (), HI masses (), -band luminosity (), stellar mass-to-light ratios (), stellar masses () for our sample, and gas fractions ().
The mean HI mass for our galaxies is , consistent with the total LSB spiral population in Honey et al. (2018). The mean stellar mass for our galaxies is . We find that the majority of our galaxies are gas rich (), with only one galaxy being gas poor, UGC 2925. Specifically, the majority of our sample have . This is consistent with barred LSBs (Pahwa & Saha, 2018) and LSBs in general (de Blok et al., 1996).
In Fig. 18 we show as a function . The solid line is the fit , slightly steeper than the relation for general LSB spirals in Honey et al. (2018), (the dashed line in Fig. 18). The scatter about our relation, , is shown as the shaded region in Fig. 18.
7 Discussion and Conclusions
We have measured the photometric bar properties (length, strength, and corotation radius) of fifteen barred LSB disc galaxies assembled from Schombert et al. (1992) and Impey et al. (1996) using optical - and -band photometry. We have found that bars in LSBs are shorter and weaker than those in HSBs, consistent with numerical simulations. Our mean bar length and strength are 2.5 kpc and 0.19 respectively. In order to determine the best measure of the bar length, we have created and analyzed fake galaxy images consisting of an exponential disc and a bar. We have found that our azimuthal bar length technique performs best out of the four tested here, and that the measure, while dependent on the inaccurate measure, also does an excellent job at measuring the bar length.
Interestingly, we find that the corotation radius of the bar falls very close to the end of the bar for our sample, implying that bars in LSBs are fast rotators. In addition, the only very slow bar is located in the only galaxy without any clear spiral structure, F563-V2. Recent mergers are thought to decrease the bar pattern speed (Gerin et al., 1990), but LSBs are known to be more isolated than HSBs (Bothun et al., 1993; Mo et al., 1994; Rosenbaum & Bomans, 2004; Du et al., 2015; Honey et al., 2018). In addition, there are no known companions to F563-V2.
As discussed in Paper I, the dark matter halos should dynamically slow down any bar that may be present (Weinberg, 1985; Debattista & Sellwood, 2000), although this only applies to nonrotating halos (see Sec. 3.3 why this is likely not the case for LSBs). While there have been reports of slow bars found in dark matter dominated objects (see Bureau et al., 1999; Chemin & Hernandez, 2009; Banerjee et al., 2013), our results imply that bars in LSBs are generally fast. The reason for this remains unclear, and requires further detailed spectroscopic analysis in order to obtain bar pattern speeds. However, the faint nature and rough discs of LSBs makes application of the only direct bar pattern speed measurement (Tremaine & Weinberg, 1984) prohibatively difficult.
In addition, it is possible that can give misleading information regarding the actual rotational speed of the bar. Font et al. (2017) found that by examining the bar pattern speeds () of a large sample of galaxies, some bars that were considered fast (i.e. ) were actually some of the slowest rotators. In order to determine this, they examined the ‘normalised’ bar pattern speed which they defined as with being the disc pattern speed. By setting a deliminator of as a means of separating fast and slow rotators, they found that bars with could have , meaning the bar rotates at a comparable rate to the disc. Because of this, the authors conclude that bars have indeed been slowed down due to the dark matter halos, despite the prevalance of bars with , and suggest as a better means of separating fast and slow rotators.
The galaxies in our sample have only been shallowly observed before. We have therefore presented - and -band surface brightness profiles, magnitudes, and colours of our sample. We have found that barred LSBs are brighter than unbarred LSBs, and fall on a shallower vs relation. The disc scale lengths of our sample are identical to those of unbarred LSBs. Our sample is also quite blue, having a mean () area colour of 1.11, while the bar region is significantly redder, with a mean colour of 1.54. Due to the small bars in our galaxies, it is not surprising that our area colours are heavily weighted towards the bluer disc.
Finally, we have used available 21cm HI fluxes to determine the HI masses and population synthesis models to get stellar masses of our sample and found our galaxies to be quite gas rich. Our sample falls on the same vs relation as the general LSB spiral population. It is thought that gas rich, low luminosity galaxies with high halo spin form short and weak bars (Cervantes-Sodi et al., 2013; Cervantes Sodi, 2017), consistent with our findings here.
We list our major conclusions based on our observed sample here.
Bars in LSBs are shorter and weaker than those in HSBs, with an average length and strength of 2.5 kpc and 0.19 respectively. 2. 2.
Bars in LSBs are fast rotators (), with the corotation radius occuring close to the end of the bar. Barred LSBs also show multiple - and -band phase intersections, possibly indicating disc corotation radii. 3. 3.
Barred LSBs are slightly brighter than unbarred LSBs, with average central surface brightnesses of mag arcsec*-2* and mag arcsec*-2*, and average absolute magnitudes of and . 4. 4.
Barred LSBs fall on a slightly shallower vs. relation than unbarred LSBs. 5. 5.
Barred LSBs have HI masses that are nearly identical to the general spiral LSB population, and are just as gas rich as unbarred LSBs (). 6. 6.
In order to fully probe the nature of the dark matter halos of barred LSBs, extensive spectroscopy is required in order to obtain bar pattern speeds and halo spin parameters ().
Acknowledgements
The authors would like to thank the anonymous referee for their very helpful comments in assisting with the discussions in this paper. In addition, would like to thank Russet McMillan for her help taking images to create a master I-band fringe pattern for the ARCTIC imager. We would also like to thank Dr. Stacy McGaugh for discussions on LSB photometry, especially morphological classifications.
This research has made use of the NASA/ IPAC Infrared Science Archive, which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This research made use of Astropy, a community-developed core Python package for Astronomy (Astropy Collaboration et al., 2013; The Astropy Collaboration et al., 2018). We acknowledge the usage of the HyperLeda database (http://leda.univ-lyon1.fr).
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Aguerri et al. (1998) Aguerri, J. A. L., Beckman, J. E., & Prieto, M. 1998, AJ, 116, 2136
- 2Aguerri et al. (2000 a) Aguerri, J. A. L., Muñoz-Tuñón, C., Varela, A. M., & Prieto, M. 2000, A&A, 361, 841
- 3Aguerri et al. (2000 b) Aguerri, J. A. L., Varela, A. M., Prieto, M., & Muñoz-Tuñón, C. 2000, AJ, 119, 1638
- 4Aguerri et al. (2009) Aguerri, J. A. L., Méndez-Abreu, J., & Corsini, E. M. 2009, A&A, 495, 491
- 5Aguerri et al. (2015) Aguerri, J. A. L., Méndez-Abreu, J., Falcón-Barroso, J., et al. 2015, A&A, 576, A 102
- 6Astropy Collaboration et al. (2013) Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A 33
- 7The Astropy Collaboration et al. (2018) The Astropy Collaboration, Price-Whelan, A. M., Sipőcz, B. M., et al. 2018, ar Xiv:1801.02634
- 8Athanassoula et al. (1990) Athanassoula, E., Morin, S., Wozniak, H., et al. 1990, MNRAS, 245, 130
