## Recent Postings from Galaxies

### What produces the far-infrared/submm emission in the most luminous QSOs?

The AGN. I examine the average spectral energy distributions (SEDs) of two samples of the most powerful, unobscured QSOs at 2<z<3.5, with rest frame optical luminosities in the 46.2<log vLv(5100Ang)<47.4 range, corresponding to the tail of the 2<z<4 QSO optical luminosity function. I find that the AGN could potentially account for the entire broad-band emission from the UV to the submm, on the basis that the SEDs of these sources are similar to the intrinsic AGN SEDs derived for lower power, lower redshift QSOs. Although this does not preclude substantial star-formation in their host galaxies, I find that the AGN dominates the total infrared luminosity, removing the necessity for a star-forming component in the far-IR/submm. I argue that the origin of the far-IR/submm emission in such powerful QSOs includes a small contribution from the AGN torus, but is predominantly linked to dust at kpc-scales heated by the AGN. The latter component accounts for at least 5-10 per cent of the bolometric AGN luminosity and has an implied dust mass of the order of 10^8 solar masses.

### Excess C/O and C/H in outer protoplanetary disk gas

The compositions of nascent planets depend on the compositions of their birth disks. In particular, the elemental compositions of Gas Giant gaseous envelopes depend on the elemental composition of the disk gas from which the envelope is accreted. Previous models demonstrated that sequential freeze-out of O and C-bearing volatiles in disks will result in an supersolar C/O ratios and subsolar C/H ratios in the gas between water and CO snowlines. This result does not take into account, however, the expected grain growth and radial drift of pebbles in disks, and the accompanying re-distribution of volatiles from the outer to the inner disk. Using a toy model we demonstrate that when drift is considered, CO is enhanced between the water and CO snowline, resulting in both supersolar C/O and C/H ratios in the disk gas in the Gas Giant formation zone. This result appears robust to the details of the disk model as long as there is substantial pebble drift across the CO snowline, and the efficiency of CO vapor diffusion is limited. Gas Giants that accrete their gaseous envelopes exterior to the water snowline and do not experience substantial core-envelope mixing, may thus present both superstellar C/O and C/H ratios in their atmospheres. Pebble drift will also affect the nitrogen and noble gas abundances in the planet forming zones, which may explain some of Jupiter's peculiar abundance patterns.

### Another look at the size of the low-surface brightness galaxy VCC 1661 in the Virgo Cluster

We present new wide-field images of the low-surface brightness Virgo Cluster dwarf galaxy VCC 1661. The extant literature lists a broad range of radii for this object, covering a factor of more than four, depending on the filters used and the details of the analyses. While some studies find a radius typical of other Virgo dwarfs and note the normality of this object, any larger spatial extent, taken at face value, would render this galaxy the largest dwarf in the Virgo Cluster samples. Confirmation of a large extent of dwarf galaxies has often led to the discovery of tidal tails and would then, also in VCC 1661, indicate a severe state of tidal disruption. Given the importance of galactic sizes for assessing tidal interactions of the satellites with their hosts, we thus combine our surface brightness profile with data from the literature to investigate further the nature of this galaxy. However, our new characteristic radius for VCC 1661 of $r_e=24.1$"$\pm7.7$" and the previously noted smooth appearance of its isophotes are fully consistent with the remainder of the ACSVCS dwarf galaxy population without any need to invoke tidal perturbations.

### A centrally heated dark halo for our Galaxy

We construct a new family of models of our Galaxy in which dark matter and disc stars are both represented by distribution functions that are analytic functions of the action integrals of motion. The potential that is self-consistently generated by the dark matter, stars and gas is determined, and parameters in the distribution functions are adjusted until the model is compatible with observational constraints on the circular-speed curve, the vertical density profile of the stellar disc near the Sun, the kinematics of nearly 200 000 giant stars within 2 kpc of the Sun, and estimates of the optical depth to microlensing of bulge stars. We find that the data require a dark halo in which the phase-space density is approximately constant for actions |J| \lesssim 140 kpc km ^-1. In real space these haloes have core radii ~ 2 kpc.

### GMRT observations of IC 711 - The longest head-tail radio galaxy known

We present low-frequency, Giant Metrewave Radio Telescope (GMRT) observations at 240, 610 and 1300 MHz of IC 711, a narrow angle tail (NAT) radio galaxy. The galaxy lies in Abell cluster 1314 (redshift ~ 0.034) and has a long radio tail of total angular extent ~17 arcmin, corresponding to a projected linear size of ~700 kpc. This makes it the longest head-tail radio galaxy known. The objectives of the GMRT observations were to investigate the diffuse-emission of the long tail structure at low frequencies. The radio structure, especially initial ~10 arcmin of tail being a long straight feature, does not seem to be consistent with a simple circular motion around the cluster center, as previously suggested in the literature. Two sharp bends after the straight section of the tail cast doubt on the prevailing idea in the literature that the long narrow tails represent trails left behind by the fast moving parent optical galaxy with respect to the cluster medium, as the optical galaxy could not have undergone such sharp bends in its path, under any conceivable gravitational influence of some individual galaxy or of the overall cluster gravitational potential. In fact the tail does not seem to have been influenced by the gravitational field of any of the cluster-member galaxy. The radio tail shows a break in the spectrum. We derive an expression for the minimum energy in the case of a spectral break, in order to do the minimum energy calculation in diffuse tail regions of IC 711.

### The galaxy counterpart of the high-metallicity and 16 kpc impact parameter DLA towards Q0918+1636 - a challenge to galaxy formation models?

The quasar Q0918+1636 (z=3.07) has an intervening high-metallicity Damped Lyman-alpha Absorber (DLA) along the line of sight, at a redshift of z=2.58. The DLA is located at a large impact parameter of 16.2 kpc, and has an almost solar metallicity. It is shown, that a novel type of cosmological galaxy formation models, invoking a new SNII feedback prescription, the Haardt & Madau (2012) UVB field and explicit treatment of UVB self-shielding, can reproduce the observed characteristics of the DLA. UV radiation from young stellar populations in the galaxy, in particular in the photon energy range 10.36-13.61 eV (relating to Sulfur II abundance), are also considered in the analysis. It is found that a) for L~L* galaxies (at z=2.58), about 10% of the sight-lines through the galaxies at impact parameter 16.2 kpc will display a Sulfur II column density N(SII)$>$ 10$^{15.82}$ cm$^{-2}$ (the observed value for the DLA), and b) considering only cases where a near-solar metallicity will be detected at 16.2 kpc impact parameter, the probability distribution of galaxy SFR peaks near the value observed for the DLA galaxy counterpart of ~27 Msun/yr. It is argued, that the bulk of the alpha-elements, like Sulfur, traced by the high metal column density, b=16.2 kpc absorption lines, were produced by evolving young stars in the inner galaxy, and later transported outward by galactic winds.

### ATCA detections of massive molecular gas reservoirs in dusty, high-z radio galaxies

Observations using the 7 mm receiver system on the Australia Telescope Compact Array have revealed large reservoirs of molecular gas in two high-redshift radio galaxies: HATLAS J090426.9+015448 (z = 2.37) and HATLAS J140930.4+003803 (z = 2.04). Optically the targets are very faint, and spectroscopy classifies them as narrow-line radio galaxies. In addition to harbouring an active galactic nucleus the targets share many characteristics of sub-mm galaxies. Far-infrared data from Herschel-ATLAS suggest high levels of dust (>10^9 M_solar) and a correspondingly large amount of obscured star formation (~1000 M_solar / yr). The molecular gas is traced via the J = 1-0 transition of 12CO, its luminosity implying total H_2 masses of (1.7 +/- 0.3) x 10^11 and (9.5 +/- 2.4) x 10^10 ({\alpha}_CO/0.8) M_solar in HATLAS J090426.9+015448 and HATLAS J140930.4+003803 respectively. Both galaxies exhibit molecular line emission over a broad (~1000 km/s) velocity range, and feature double-peaked profiles. We interpret this as evidence of either a large rotating disk or an on-going merger. Gas depletion timescales are ~100 Myr. The 1.4 GHz radio luminosities of our targets place them close to the break in the luminosity function. As such they represent `typical' z > 2 radio sources, responsible for the bulk of the energy emitted at radio wavelengths from accretion-powered sources at high redshift, and yet they rank amongst the most massive systems in terms of molecular gas and dust content. We also detect 115 GHz rest-frame continuum emission, indicating a very steep high-radio-frequency spectrum, possibly classifying the targets as compact steep spectrum objects.

### The X-ray luminosity-temperature relation of a complete sample of low-mass galaxy clusters

We present \Chandra\ observations of 23 galaxy groups and low-mass galaxy clusters at $0.03<z<0.15$ with a median temperature of ~2keV. The sample is a statistically complete flux-limited subset of the 400 deg$^2$ survey. We investigated the scaling relation between X-ray luminosity ($L$) and temperature ($T$), taking selection biases fully into account. The logarithmic slope of the bolometric \LT\ relation was found to be $3.29\pm0.33$, consistent with values typically found for samples of more massive clusters. In combination with other recent studies of the \LT\ relation we show that there is no evidence for the slope, normalisation, or scatter of the \LT\ relation of galaxy groups being different than that of massive clusters. The exception to this is that in the special case of the most relaxed systems, the slope of the core-excised \LT\ relation appears to steepen from the self-similar value found for massive clusters to a steeper slope for the lower mass sample studied here. Thanks to our rigorous treatment of selection biases, these measurements provide a robust reference against which to compare predictions of models of the impact of feedback on the X-ray properties of galaxy groups.

### The Mass-Discrepancy Acceleration Relation: a Natural Outcome of Galaxy Formation in CDM halos

We analyze the total and baryonic acceleration profiles of a set of well-resolved galaxies identified in the EAGLE suite of hydrodynamic simulations. Our runs start from the same initial conditions but adopt different subgrid models for stellar and AGN feedback, resulting in diverse populations of galaxies by the present day. Some of them reproduce observed galaxy scaling relations, while others do not. However, regardless of the feedback implementation, all of our galaxies follow closely a simple relationship between the total and baryonic acceleration profiles, consistent with recent observations of rotationally supported galaxies. The relation has small scatter: different feedback processes -- which produce different galaxy populations -- mainly shift galaxies along the relation, rather than perpendicular to it. Furthermore, galaxies exhibit a single characteristic acceleration, $g_{\dagger}$, above which baryons dominate the mass budget, as observed. These observations have been hailed as evidence for modified Newtonian dynamics but can be accommodated within the standard cold dark matter paradigm.

### The SELGIFS data challenge: generating synthetic observations of CALIFA galaxies from hydrodynamical simulations

In this work we present a set of synthetic observations that mimic the properties of the Integral Field Spectroscopy (IFS) survey CALIFA, generated using radiative transfer techniques applied to hydrodynamical simulations of galaxies in a cosmological context. The simulated spatially-resolved spectra include stellar and nebular emission, kinematic broadening of the lines, and dust extinction and scattering. The results of the radiative transfer simulations have been post-processed to reproduce the main properties of the CALIFA V500 and V1200 observational setups. The data has been further formatted to mimic the CALIFA survey in terms of field of view size, spectral range and sampling. We have included the effect of the spatial and spectral Point Spread Functions affecting CALIFA observations, and added detector noise after characterizing it on a sample of 20 galaxies. The simulated datacubes are suited to be analyzed by the same algorithms used on real IFS data. In order to provide a benchmark to compare the results obtained applying IFS observational techniques to our synthetic datacubes, and test the calibration and accuracy of the analysis tools, we have computed the spatially-resolved properties of the simulations. Hence, we provide maps derived directly from the hydrodynamical snapshots or the noiseless spectra, in a way that is consistent with the values recovered by the observational analysis algorithms. Both the synthetic observations and the product datacubes are public and can be found in the collaboration website http://astro.ft.uam.es/selgifs/data_challenge/.

### Does turbulence determine the initial mass function? [Replacement]

We test the hypothesis that the initial mass function (IMF) is determined by the density probability distribution function (PDF) produced by supersonic turbulence. We compare 14 simulations of star cluster formation in 50 solar mass molecular cloud cores where the initial turbulence contains either purely solenoidal or purely compressive modes, in each case resolving fragmentation to the opacity limit to determine the resultant IMF. We find statistically indistinguishable IMFs between the two sets of calculations, despite a factor of two difference in the star formation rate and in the standard deviation of $\log(\rho)$. This suggests that the density PDF, while determining the star formation rate, is not the primary driver of the IMF.

### Does turbulence determine the initial mass function?

We test the hypothesis that the initial mass function (IMF) is determined by the density probability distribution function (PDF) produced by supersonic turbulence. We compare 14 simulations of star cluster formation in 50 solar mass molecular cloud cores where the initial turbulence contains either purely solenoidal or purely compressive modes, in each case resolving fragmentation to the opacity limit to determine the resultant IMF. We find statistically indistinguishable IMFs between the two sets of calculations, despite a factor of two difference in the star formation rate and in the standard deviation of $\log(\rho)$. This suggests that the density PDF, while determining the star formation rate, is not the primary driver of the IMF.

### The Tarantula Massive Binary Monitoring project: II. A first SB2 orbital and spectroscopic analysis for the Wolf-Rayet binary R145

We present the first SB2 orbital solution and disentanglement of the massive Wolf-Rayet binary R145 (P = 159d) located in the Large Magellanic Cloud. The primary was claimed to have a stellar mass greater than 300Msun, making it a candidate for the most massive star known. While the primary is a known late type, H-rich Wolf-Rayet star (WN6h), the secondary could not be so far unambiguously detected. Using moderate resolution spectra, we are able to derive accurate radial velocities for both components. By performing simultaneous orbital and polarimetric analyses, we derive the complete set of orbital parameters, including the inclination. The spectra are disentangled and spectroscopically analyzed, and an analysis of the wind-wind collision zone is conducted. The disentangled spectra and our models are consistent with a WN6h type for the primary, and suggest that the secondary is an O3.5 If*/WN7 type star. We derive a high eccentricity of e = 0.78 and minimum masses of M1 sin^3 i ~ M2 sin^3 i ~ 13 +- 2 Msun, with q = M2 / M1 = 1.01 +- 0.07. An analysis of emission excess stemming from a wind-wind collision yields a similar inclination to that obtained from polarimetry (i = 39 +- 6deg). Our analysis thus implies M1 = 53^{+40}_{-20} and M2 = 54^{+40}_{-20} Msun, excluding M1 > 300Msun. A detailed comparison with evolution tracks calculated for single and binary stars, as well as the high eccentricity, suggest that the components of the system underwent quasi-homogeneous evolution and avoided mass-transfer. This scenario would suggest current masses of ~ 80 Msun and initial masses of Mi,1 ~ 105 and Mi,2 ~ 90Msun, consistent with the upper limits of our derived orbital masses, and would imply an age of ~2.2 Myr.

### Galactic rotation in Gaia DR1

The spatial variations of the velocity field of local stars provide direct evidence of Galactic differential rotation. The local divergence, shear, and vorticity of the velocity field---the traditional Oort constants---can be measured based purely on astrometric measurements and in particular depend linearly on proper motion and parallax. I use data for 304,267 main-sequence stars from the Gaia DR1 Tycho-Gaia Astrometric Solution to perform a local, precise measurement of the Oort constants at a typical heliocentric distance of 230 pc. The pattern of proper motions for these stars clearly displays the expected effects from differential rotation. I measure the Oort constants to be: A = 15.3+/-0.4+/- 0.5 (syst.) km/s/kpc, B = -11.9+/-0.4 km/s/kpc, C = -3.2+/-0.4 km/s/kpc and K = -3.3+/-0.6 km/s/kpc, where the systematic uncertainty in A is due to its uncertain correction for the asymmetric drift. These measurements imply local values for the rotational frequency of \Omega_0 = 27.1+/-0.5+/- 0.5 (syst.) km/s/kpc, for the slope of the rotation curve of -3.4+/-0.6+/-0.5 (syst.) km/s/kpc, and for the circular velocity of 219+/-4+/-4 (syst.) km/s.

### The deep OB star population in Carina from the VST Photometric H$\alpha$ Survey (VPHAS+)

Massive OB stars are critical to the ecology of galaxies, and yet our knowledge of OB stars in the Milky Way, fainter than $V \sim 12$, remains patchy. Data from the VST Photometric H$\alpha$ Survey (VPHAS+) permit the construction of the first deep catalogues of blue excess-selected OB stars, without neglecting the stellar field. A total of 14900 candidates with 2MASS cross-matches are blue-selected from a 42 square-degree region in the Galactic Plane, capturing the Carina Arm over the Galactic longitude range $282^{\circ} \lesssim \ell \lesssim 293^{\circ}$. Spectral energy distribution fitting is performed on these candidates' combined VPHAS+ $u,g,r,i$ and 2MASS $J,H,K$ magnitudes. This delivers: effective temperature constraints, statistically separating O from early-B stars; high-quality extinction parameters, $A_0$ and $R_V$ (random errors typically $< 0.1$). The high-confidence O-B2 candidates number 5915 and a further 5170 fit to later B spectral type. Spectroscopy of 276 of the former confirms 97% of them. The fraction of emission line stars among all candidate B stars is 7--8% . Greyer ($R_V > 3.5$) extinction laws are ubiquitous in the region, over the distance range 2.5--3 kpc to $\sim$10~kpc. Near prominent massive clusters, $R_V$ tends to rise, with particularly large and chaotic excursions to $R_V \sim 5$ seen in the Carina Nebula. The data reveal a hitherto unnoticed association of 108 O-B2 stars around the O5If$+$ star LSS 2063 ($\ell = 289.77^{\circ}$, $b = -1.22^{\circ}$). Treating the OB star scale-height as a constant within the thin disk, we find an orderly mean relation between extinction ($A_0$) and distance in the Galactic longitude range, $287.6^{\circ} < \ell < 293.5^{\circ}$, and infer the subtle onset of thin-disk warping. A halo around NGC 3603, roughly a degree in diameter, of $\sim$500 O-B2 stars with $4 < A_0 (\rm{mag}) < 7$ is noted.

### The inner structure of early-type galaxies in the Illustris simulation

Early-type galaxies provide unique tests for the predictions of the cold dark matter cosmology and the baryonic physics assumptions entering models for galaxy formation. In this work, we use the Illustris simulation to study correlations of three main properties of early-type galaxies, namely, the stellar orbital anisotropies, the central dark matter fractions and the central radial density slopes, as well as their redshift evolution since $z=1.0$. We find that lower-mass galaxies or galaxies at higher redshift tend to be bluer in rest-frame colour, have higher central gas fractions, and feature more tangentially anisotropic orbits and steeper central density slopes than their higher-mass or lower-redshift counterparts, respectively. The projected central dark matter fraction within the effective radius shows no significant mass dependence but positively correlates with galaxy effective radii due to the aperture effect. The central density slopes obtained in the simulation by combining strong lensing measurements with single aperture kinematics are found to be shallower than the true density slopes. We identify systematic biases in this measurement due to two common modelling assumptions, isotropic stellar orbital distributions and power-law density profiles. We also compare the properties of early-type galaxies in Illustris to those from the SLACS, SL2S and BOSS surveys, finding in general broad agreement but also some tension, which appears to be mostly caused by too large galaxy sizes in Illustris.

### The segregation of baryons and dark matter during halo assembly

The standard galaxy formation theory assumes that baryons and dark matter are initially well-mixed before becoming segregated due to radiative cooling. We use non-radiative hydrodynamical simulations to explicitly examine this assumption and find that baryons and dark matter can also be segregated during the build-up of the halo. As a result, baryons in many haloes do not originate from the same Lagrangian region as the dark matter. When using the fraction of corresponding dark matter and gas particles in the initial conditions (the "paired fraction") as a proxy of the dark matter and gas segregation strength of a halo, on average about $25$ percent of the baryonic and dark matter of the final halo are segregated in the initial conditions. A consequence of this effect is that the baryons and dark matter of the same halo initially experience different tidal torques and thus their angular momentum vectors are often misaligned. This is at odds with the assumption of the standard galaxy formation model, and challenges the precision of some semi-analytical approaches which utilize dark matter halo merger trees to infer properties of gas associated to dark matter haloes.

### Gamma Rays From Dark Matter Subhalos Revisited: Refining the Predictions and Constraints

Utilizing data from the ELVIS and Via Lactea-II simulations, we characterize the local dark matter subhalo population, and use this information to refine the predictions for the gamma-ray fluxes arising from annihilating dark matter in this class of objects. We find that the shapes of nearby subhalos are significantly altered by tidal effects, and are generally not well described by NFW density profiles, instead prefering power-law profiles with an exponential cutoff. From the subhalo candidates detected by the Fermi Gamma-Ray Space Telescope, we place limits on the dark matter annihilation cross section that are only modestly weaker than those based on observations of dwarf galaxies. We also calculate the fraction of observable subhalos that are predicted to be spatially extended at a level potentially discernible to Fermi.

### Gamma Rays From Dark Matter Subhalos Revisited: Refining the Predictions and Constraints [Cross-Listing]

Utilizing data from the ELVIS and Via Lactea-II simulations, we characterize the local dark matter subhalo population, and use this information to refine the predictions for the gamma-ray fluxes arising from annihilating dark matter in this class of objects. We find that the shapes of nearby subhalos are significantly altered by tidal effects, and are generally not well described by NFW density profiles, instead prefering power-law profiles with an exponential cutoff. From the subhalo candidates detected by the Fermi Gamma-Ray Space Telescope, we place limits on the dark matter annihilation cross section that are only modestly weaker than those based on observations of dwarf galaxies. We also calculate the fraction of observable subhalos that are predicted to be spatially extended at a level potentially discernible to Fermi.

### The Grism Lens-Amplified Survey from Space (GLASS) X. Sub-kpc resolution gas-phase metallicity maps at cosmic noon behind the Hubble Frontier Fields cluster MACS1149.6+2223 [Replacement]

(Abridged) We combine deep HST grism spectroscopy with a new Bayesian method to derive maps of gas-phase metallicity, nebular dust extinction, and star-formation rate for 10 star-forming galaxies at high redshift ($1.2<z<2.3$). Exploiting lensing magnification by the foreground cluster MACS1149.6+2223, we reach sub-kpc spatial resolution and push the stellar mass limit associated with such high-z spatially resolved measurements below $10^8M_\odot$ for the first time. Our maps exhibit diverse morphologies, indicative of various effects such as efficient radial mixing from tidal torques, rapid accretion of low-metallicity gas, etc., which can affect the gas and metallicity distributions in individual galaxies. Based upon an exhaustive sample of all existing sub-kpc metallicity gradients at high-z, we find that predictions given by analytical chemical evolution models assuming a relatively extended star-formation profile in the early disk formation phase can explain the majority of observed gradients, without involving galactic feedback or radial outflows. We observe a tentative correlation between stellar mass and metallicity gradient, consistent with the downsizing galaxy formation picture that more massive galaxies are more evolved into a later phase of disk growth, where they experience more coherent mass assembly at all radii and thus show shallower metallicity gradients. In addition, we compile a sample of homogeneously cross-calibrated integrated metallicity measurements spanning three orders of magnitude in stellar mass at $z\sim1.8$. We use this sample to study the mass-metallicity relation (MZR) and test the fundamental metallicity relation (FMR). The slope of the observed MZR can rule out the momentum-driven wind model at 3-$\sigma$ confidence level. We find no significant offset with respect to the FMR, taking into account the intrinsic scatter and measurement uncertainties.

### The Grism Lens-Amplified Survey from Space (GLASS) X. Sub-kpc resolution gas-phase metallicity maps at cosmic noon behind the Hubble Frontier Fields cluster MACS1149.6+2223

(Abridged) We combine deep HST grism spectroscopy with a new Bayesian method to derive maps of gas-phase metallicity, nebular dust extinction, and star-formation rate for 10 star-forming galaxies at high redshift ($1.2<z<2.3$). Exploiting lensing magnification by the foreground cluster MACS1149.6+2223, we reach sub-kpc spatial resolution and push the stellar mass limit associated with such high-z spatially resolved measurements below $10^8M_\odot$ for the first time. Our maps exhibit diverse morphologies, indicative of various effects such as efficient radial mixing from tidal torques, rapid accretion of low-metallicity gas, etc., which can affect the gas and metallicity distributions in individual galaxies. Based upon an exhaustive sample of all existing sub-kpc metallicity gradients at high-z, we find that predictions given by analytical chemical evolution models assuming a relatively extended star-formation profile in the early disk formation phase can explain the majority of observed gradients, without involving galactic feedback or radial outflows. We observe a tentative correlation between stellar mass and metallicity gradient, consistent with the downsizing galaxy formation picture that more massive galaxies are more evolved into a later phase of disk growth, where they experience more coherent mass assembly at all radii and thus show shallower metallicity gradients. In addition, we compile a sample of homogeneously cross-calibrated integrated metallicity measurements spanning three orders of magnitude in stellar mass at $z\sim1.8$. We use this sample to study the mass-metallicity relation (MZR) and test the fundamental metallicity relation (FMR). The slope of the observed MZR can rule out the momentum-driven wind model at 3-$\sigma$ confidence level. We find no significant offset with respect to the FMR, taking into account the intrinsic scatter and measurement uncertainties.

### Magnetic flux density from the relative circular motion of stars and partially ionized gas in the Galaxy mid-plane vicinity

Observations suggest a slower stellar rotation relative to gas rotation in the outer part of the Milky Way Galaxy. This difference could be attributed to an interaction with the interstellar magnetic field. In a simple model, fields of order 10 micro Gauss are then required, consistently with the observed values. This coincidence suggests a tool for estimating magnetic fields in spiral galaxies. A North-South asymmetry in the rotation of gas in the Galaxy could be of magnetic origin too.

### The $\Lambda$CDM simulations of Keller and Wadsley do not account for the MOND mass-discrepancy-acceleration relation

Keller and Wadsley (2016) have smugly suggested, recently, that the end of MOND may be in view. This is based on their claim that their highly-restricted sample of $\Lambda$CDM-simulated galaxies are "consistent" with the observed MOND mass-discrepancy-acceleration relation (MDAR), in particular, with its recent update by McGaugh et al. (2016), based on the SPARC sample. From this they extrapolate to "$\Lambda$CDM is fully consistent" with the MDAR. I explain why these simulated galaxies do not show that $\Lambda$CDM accounts for the MDAR. a. Their sample of simulated galaxies contains only 18 high-mass galaxies, within a narrow range of one order of magnitude in baryonic mass, at the very high end of the observed, SPARC sample, which spans 4.5 orders of magnitude in mass. More importantly, the simulated sample has none of the low-mass, low-acceleration galaxies -- abundant in SPARC -- which encapsulate the crux and the nontrivial aspects of the predicted and observed MDAR. The low-acceleration part of the simulated MDAR is achieved, rather trivially, from the flattish-velocity-curve regions of the simulated high-mass galaxies. b. Half of the simulated galaxies have "wrong" rotation curves that differ greatly from any observed ones. This, does not prevent these wrong galaxies from lying on the observed MDAR (for trivial reasons, again). They, in fact, define the high-acceleration branch of the simulated MDAR. c. To boot, even if $\Lambda$CDM were made "consistent" with the MDAR through the elaborate adjustments that the simulations allow, this would not obviate MOND, which predicts much more than the MDAR.

### The $\Lambda$CDM simulations of Keller and Wadsley do not account for the MOND mass-discrepancy-acceleration relation [Cross-Listing]

Keller and Wadsley (2016) have smugly suggested, recently, that the end of MOND may be in view. This is based on their claim that their highly-restricted sample of $\Lambda$CDM-simulated galaxies are "consistent" with the observed MOND mass-discrepancy-acceleration relation (MDAR), in particular, with its recent update by McGaugh et al. (2016), based on the SPARC sample. From this they extrapolate to "$\Lambda$CDM is fully consistent" with the MDAR. I explain why these simulated galaxies do not show that $\Lambda$CDM accounts for the MDAR. a. Their sample of simulated galaxies contains only 18 high-mass galaxies, within a narrow range of one order of magnitude in baryonic mass, at the very high end of the observed, SPARC sample, which spans 4.5 orders of magnitude in mass. More importantly, the simulated sample has none of the low-mass, low-acceleration galaxies -- abundant in SPARC -- which encapsulate the crux and the nontrivial aspects of the predicted and observed MDAR. The low-acceleration part of the simulated MDAR is achieved, rather trivially, from the flattish-velocity-curve regions of the simulated high-mass galaxies. b. Half of the simulated galaxies have "wrong" rotation curves that differ greatly from any observed ones. This, does not prevent these wrong galaxies from lying on the observed MDAR (for trivial reasons, again). They, in fact, define the high-acceleration branch of the simulated MDAR. c. To boot, even if $\Lambda$CDM were made "consistent" with the MDAR through the elaborate adjustments that the simulations allow, this would not obviate MOND, which predicts much more than the MDAR.

### Radial acceleration relation from screening of fifth forces

We show that the radial acceleration relation for rotationally-supported galaxies may be explained, in the absence of cold dark matter, by a non-minimally coupled scalar field, whose fifth forces are partially screened on galactic scales by the symmetron mechanism. If realised in nature, this effect could have a significant impact on the inferred density of dark matter halos.

### Radial acceleration relation from screening of fifth forces [Cross-Listing]

We show that the radial acceleration relation for rotationally-supported galaxies may be explained, in the absence of cold dark matter, by a non-minimally coupled scalar field, whose fifth forces are partially screened on galactic scales by the symmetron mechanism. If realised in nature, this effect could have a significant impact on the inferred density of dark matter halos.

### Radial acceleration relation from screening of fifth forces [Cross-Listing]

We show that the radial acceleration relation for rotationally-supported galaxies may be explained, in the absence of cold dark matter, by a non-minimally coupled scalar field, whose fifth forces are partially screened on galactic scales by the symmetron mechanism. If realised in nature, this effect could have a significant impact on the inferred density of dark matter halos.

### Seyfert 1 Composite Spectrum using SDSS Legacy Survey Data

We present a rest-frame composite spectrum for Seyfert 1 galaxies using spectra obtained from the DR12 release of the Sloan Digital Sky Survey (SDSS). The spectrum is constructed by combining data from a total of 10,112 galaxies, spanning a redshift range of 0 to 0.793. We produce an electronic table of the median and geometric mean composite Seyfert 1 spectrum. We measure the spectral index of the composite spec- trum, and compare it with that of the composite quasar spectrum. We also measure the flux and width of the strong emission lines present in the composite spectrum. We compare the entire spectrum with the quasar spectrum in the context of the AGN unification model. The two composite spectra match extremely well in the blue part of the spectrum, while there is an offset in flux in the red portion of the spectrum.

### Galactic winds and the origin of large-scale magnetic fields

Observations of dwarf galaxies suggest the presence of large-scale magnetic fields. However the size and slow rotation of these galaxies appear insufficient to support a mean-field dynamo action to excite such fields. Here we suggest a new mechanism to explain large-scale magnetic fields in galaxies that are too small to support mean-field dynamo action. The key idea is that we do not identify large-scale and mean magnetic fields. In our scenario the the magnetic structures originate from a small-scale dynamo which produces small-scale magnetic field in the galactic disc and a galactic wind that transports this field into the galactic halo where the large turbulent diffusion increases the scale and order of the field. As a result, the magnetic field becomes large-scale; however its mean value remains vanishing in a strict sense. We verify the idea by numerical modelling of two distinct simplified configurations, a thin disc model using the no-$z$ approximation, and an axisymmetric model using cylindrical $r,z$ coordinates. Each of these allows reduction of the problem to two spatial dimensions. Taken together, the models support the proposition that the general trends will persist in a fully 3D model. We demonstrate that a pronounced large-scale pattern can develop in the galactic halo for a wide choice of the dynamo governing parameters. We believe that our mechanism can be relevant to explaining the presence of the fields observed in the halos of dwarf galaxies. We emphasize that detailed modelling of the proposed scenario needs 3D simulations, and adjustment to the specific dynamo governing parameters of dwarf galaxies.

### Water delivery from cores to disks: deuteration as a probe of the prestellar inheritance of H2O

We investigate the delivery of regular and deuterated forms of water from prestellar cores to circumstellar disks. We adopt a semi-analytical axisymmetric two-dimensional collapsing core model with post-processing gas-ice astrochemical simulations, in which a layered ice structure is considered. The physical and chemical evolutions are followed until the end of the main accretion phase. When mass averaged over the whole disk, a forming disk has a similar H2O abundance and HDO/H2O abundance ratio as their precollapse values (within a factor of 2), regardless of time in our models. Consistent with previous studies, our models suggest that interstellar water ice is delivered to forming disks without significant alteration. On the other hand, the local vertically averaged H2O ice abundance and HDO/H2O ice ratio can differ more, by up to a factor of several, depending on time and distance from a central star. Key parameters for the local variations are the fluence of stellar UV photons en route into the disk and the ice layered structure, the latter of which is mostly established in the prestellar stages. We also find that even if interstellar water ice is destroyed by stellar UV and (partly) reformed prior to disk entry, the HDO/H2O ratio in reformed water ice is similar to the original value. This finding indicates that some caution is needed in discussions on the prestellar inheritance of H2O based on comparisons between the observationally derived HDO/H2O ratio in clouds/cores and that in disks/comets. Alternatively, we propose that the ratio of D2O/HDO to HDO/H2O better probes the prestellar inheritance of H2O. It is also found that icy organics are more enriched in deuterium than water ice in forming disks. The differential deuterium fractionation in water and organics is inherited from the prestellar stages.

### GOTHIC: Gravitational oct-tree code accelerated by hierarchical time step controlling

The tree method is a widely implemented algorithm for collisionless $N$-body simulations in astrophysics well suited for GPU(s). Adopting hierarchical time stepping can accelerate $N$-body simulations; however, it is infrequently implemented and its potential remains untested in GPU implementations. We have developed a Gravitational Oct-Tree code accelerated by HIerarchical time step Controlling named \texttt{GOTHIC}, which adopts both the tree method and the hierarchical time step. The code adopts some adaptive optimizations by monitoring the execution time of each function on-the-fly and minimizes the time-to-solution by balancing the measured time of multiple functions. Results of performance measurements with realistic particle distribution performed on NVIDIA Tesla M2090, K20X, and GeForce GTX TITAN X, which are representative GPUs of the Fermi, Kepler, and Maxwell generation of GPUs, show that the hierarchical time step achieves a speedup by a factor of around 3--5 times compared to the shared time step. The measured elapsed time per step of \texttt{GOTHIC} is 0.30~s or 0.44~s on GTX TITAN X when the particle distribution represents the Andromeda galaxy or the NFW sphere, respectively, with $2^{24} =$~16,777,216 particles. The averaged performance of the code corresponds to 10--30\% of the theoretical single precision peak performance of the GPU.

### A Submillimeter Continuum Survey of Local Dust-Obscured Galaxies

We conduct a 350 micron dust continuum emission survey of 17 dust-obscured galaxies (DOGs) at z = 0.05-0.08 with the Caltech Submillimeter Observatory (CSO). We detect 14 DOGs with S_350 = 114-650 mJy and S/N > 3. By including two additional DOGs with submillimeter data in the literature, we are able to study dust contents for a sample of 16 local DOGs that consists of 12 bump and 4 power-law types. We determine their physical parameters with a two-component modified blackbody function model. The derived dust temperatures are in the range 57-122 K and 22-35 K for the warm and cold dust components, respectively. The total dust mass and the mass fraction of warm dust component are 3-34$\times10^{7} M_\odot$ and 0.03-2.52%, respectively. We compare these results with those of other submillimeter-detected infrared luminous galaxies. The bump DOGs, the majority of the DOG sample, show similar distributions of dust temperatures and total dust mass to the comparison sample. The power-law DOGs show a hint of smaller dust masses than other samples, but need to be tested with a larger sample. These findings support that the reason why DOGs show heavy dust obscuration is not an overall amount of dust content, but probably the spatial distribution of dust therein.

### Radiative transfer modelling of parsec-scale dusty warped discs

Warped discs have been found on (sub-)parsec scale in some nearby Seyfert nuclei, identified by their maser emission. Using dust radiative transfer simulations we explore their observational signatures in the infrared in order to find out whether they can partly replace the molecular torus. Strong variations of the brightness distributions are found, depending on the orientation of the warp with respect to the line of sight. Whereas images at short wavelengths typically show a disc-like and a point source component, the warp itself only becomes visible at far-infrared wavelengths. A similar variety is visible in the shapes of the spectral energy distributions. Especially for close to edge-on views, the models show silicate feature strengths ranging from deep absorption to strong emission for variations of the lines of sight towards the warp. To test the applicability of our model, we use the case of the Circinus galaxy, where infrared interferometry has revealed a highly elongated emission component matching a warped maser disc in orientation and size. Our model is for the first time able to present a physical explanation for the observed dust morphology as coming from the AGN heated dust. As opposed to available torus models, a warped disc morphology produces a variety of silicate feature shapes for grazing lines of sight, close to an edge-on view. This could be an attractive alternative to a claimed change of the dust composition for the case of the nearby Seyfert 2 galaxy NGC 1068, which harbours a warped maser disc as well.

### On the nature and correction of the spurious S-wise spiral galaxy winding bias in Galaxy Zoo 1

The Galaxy Zoo 1 catalog displays a bias towards the S-wise winding direction in spiral galaxies which has yet to be explained. The lack of an explanation confounds our attempts to verify the Cosmological Principle, and has spurred some debate as to whether a bias exists in the real universe. The bias manifests not only in the obvious case of trying to decide if the universe as a whole has a winding bias, but also in the more insidious case of selecting which galaxies to include in a winding direction survey. While the former bias has been accounted for in a previous image-mirroring study, the latter has not. Furthermore, the bias has never been {\em corrected} in the GZ1 catalog, as only a small sample of the GZ1 catalog was re-examined during the mirror study. We show that the existing bias is a human {\em selection} effect rather than a human chirality bias. In effect, the excess S-wise votes are spuriously "stolen" from the elliptical and edge-on-disk categories, not the Z-wise category. Thus, when selecting a set of spiral galaxies by imposing a threshold $T$ so that $\max(P_S,P_Z) > T$ or $P_S+P_Z>T$, we spuriously select more S-wise than Z-wise galaxies. We show that when a provably unbiased machine selects which galaxies are spirals independent of their chirality, the S-wise surplus vanishes, even if humans are still used to determine the chirality. Thus, when viewed across the entire GZ1 sample (and by implication, the Sloan catalog), the winding direction of arms in spiral galaxies as viewed from Earth is consistent with the flip of a fair coin.

### The Local Bubble in the interstellar medium and the origin of the low energy cosmic rays

An analysis of the energy spectra of cosmic rays and particularly the precise data from the AMS-02 experiment support the view about the important role of the Local Bubble in the nearby interstellar medium. It is suggested that the bulk of cosmic rays below about 200 GV of rigidity (momentum/charge ratio) comes from the modest number of supernova remnants in the Local Bubble which appear to have occurred some $10^6$ years ago and contributed to its formation. At higher rigidities the contribution from a 'Local Source', a single supernova remnant generated some $10^5$ years ago seems to dominate up to, at least 1000 GV.

### Lick Indices and Spectral Energy Distribution Analysis based on an M31 Star Cluster Sample: Comparisons of Methods and Models

Application of fitting techniques to obtain physical parameters---such as ages, metallicities, and $\alpha$-element to iron ratios---of stellar populations is an important approach to understand the nature of both galaxies and globular clusters (GCs). In fact, fitting methods based on different underlying models may yield different results, and with varying precision. In this paper, we have selected 22 confirmed M31 GCs for which we do not have access to previously known spectroscopic metallicities. Most are located at approximately one degree (in projection) from the galactic center. We performed spectroscopic observations with the 6.5 m MMT telescope, equipped with its Red Channel Spectrograph. Lick/IDS absorption-line indices, radial velocities, ages, and metallicities were derived based on the $\rm EZ\_Ages$ stellar population parameter calculator. We also applied full spectral fitting with the ULySS code to constrain the parameters of our sample star clusters. In addition, we performed $\chi^2_{\rm min}$ fitting of the clusters' Lick/IDS indices with different models, including the Bruzual & Charlot models (adopting Chabrier or Salpeter stellar initial mass functions and 1994 or 2000 Padova stellar evolutionary tracks), the GALEV, and the Thomas et al. models. For comparison, we collected their $UVBRIJK$ photometry from the Revised Bologna Catalogue (v.5) to obtain and fit the GCs' spectral energy distributions (SEDs). Finally, we performed fits using a combination of Lick/IDS indices and SEDs. The latter results are more reliable and the associated error bars become significantly smaller than those resulting from either our Lick/IDS indices-only or our SED-only fits.

### New fully empirical calibrations of strong-line metallicity indicators in star forming galaxies

We derive new empirical calibrations for strong-line diagnostics of gas phase metallicity in local star forming galaxies by uniformly applying the Te method over the full metallicity range probed by the Sloan Digital Sky Survey (SDSS). To measure electron temperatures at high metallicity, where the auroral lines needed are not detected in single galaxies, we stacked spectra of more than 110,000 galaxies from the SDSS in bins of log[O II]/H$\beta$ and log[O III]/H$\beta$. This stacking scheme does not assume any dependence of metallicity on mass or star formation rate, but only that galaxies with the same line ratios have the same oxygen abundance. We provide calibrations which span more than 1 dex in metallicity and are entirely defined on a consistent absolute Te metallicity scale for galaxies. We apply our calibrations to the SDSS sample and find that they provide consistent metallicity estimates to within 0.05 dex.

### The RESOLVE Survey Atomic Gas Census and Environmental Influences on Galaxy Gas Reservoirs

We present the HI mass inventory for the RESOLVE survey, a volume-limited, multi-wavelength census of >1500 z=0 galaxies spanning diverse environments and complete in baryonic mass down to dwarfs of 10^9 Msun. This first 21cm data release provides robust detections or strong upper limits (1.4M_HI < 5 to 10% of stellar mass M_stars) for 94% of RESOLVE. We examine global atomic gas-to-stellar mass ratios (G/S) in relation to galaxy environment using several metrics: group dark matter halo mass M_h , central/satellite designation, relative mass density of the cosmic web, and distance to nearest massive group. We find that at fixed M_stars, satellites have decreasing G/S with increasing M_h starting clearly at M_h = 10^12 Msun, suggesting the presence of starvation and/or stripping mechanisms associated with halo gas heating in intermediate-mass groups. The analogous relationship for centrals is uncertain because halo abundance matching builds in relationships between central G/S, stellar mass, and halo mass, which depend on the integrated group property used as a proxy for halo mass (stellar or baryonic mass). On larger scales G/S trends are less sensitive to the abundance matching method. At fixed M_h < 10^12 Msun, the fraction of gas-poor centrals increases with large-scale structure density. In overdense regions, we identify a rare population of gas-poor centrals in low-mass (M_h < 10^11.4 Msun) halos primarily located within 1.5 times the virial radius of more massive (M_h > 10^12 Msun) halos, suggesting that gas stripping and/or starvation may be induced by interactions with larger halos or the surrounding cosmic web. We find that the detailed relationship between G/S and environment varies when we examine different subvolumes of RESOLVE independently, which we suggest may be a signature of assembly bias.

### The second closest gamma-ray burst: sub-luminous GRB 111005A with no supernova in a super-solar metallicity environment

We report the detection of the radio afterglow of a long gamma-ray burst (GRB) 111005A at 5-345 GHz, including the very long baseline interferometry observations with the positional error of 0.2 mas. The afterglow position is coincident with the disk of a galaxy ESO 580-49 at z= 0.01326 (~1" from its center), which makes GRB 111005A the second closest GRB known to date, after GRB 980425. The radio afterglow of GRB 111005A was an order of magnitude less luminous than those of local low-luminosity GRBs, and obviously than those of cosmological GRBs. The radio flux was approximately constant and then experienced an unusually rapid decay a month after the GRB explosion. Similarly to only two other GRBs, we did not find the associated supernovae (SN), despite deep near- and mid-infrared observations 1-9 days after the GRB explosion, reaching ~20 times fainter than other SNe associated with GRBs. Moreover, we measured twice solar metallicity for the GRB location. The low gamma-ray and radio luminosities, rapid decay, lack of a SN, and super-solar metallicity suggest that GRB 111005A represents a different rare class of GRBs than typical core-collapse events. We modelled the spectral energy distribution of the GRB 111005A host finding that it is a dwarf, moderately star-forming galaxy, similar to the host of GRB 980425. The existence of two local GRBs in such galaxies is still consistent with the hypothesis that the GRB rate is proportional to the cosmic star formation rate (SFR) density, but suggests that the GRB rate is biased towards low SFRs. Using the far-infrared detection of ESO 580-49, we conclude that the hosts of both GRBs 111005A and 980425 exhibit lower dust content than what would be expected from their stellar masses and optical colours.

### The overlooked role of stellar variability in the extended main sequence of LMC intermediate-age clusters

Intermediate-age star clusters in the Large Magellanic Cloud show extended main sequence turn offs (MSTOs), which are not consistent with a canonical single stellar population. These broad turn offs have been interpreted as evidence for extended star formation and/or stellar rotation. Since most of these studies use single frames per filter to do the photometry, the presence of variable stars near the MSTO in these clusters has remained unnoticed and their impact totally ignored. We model the influence of Delta Scuti using synthetic CMDs, adding variable stars following different levels of incidence and amplitude distributions. We show that Delta Scuti observed at a single phase will produce a broadening of the MSTO without affecting other areas of a CMD like the upper MS or the red clump; furthermore, the amount of spread introduced correlates with cluster age as observed. This broadening is constrained to ages ~ 1-3 Gyr when the MSTO area crosses the instability strip, which is also consistent with observations. Variable stars cannot explain bifurcated MSTOs or the extended MSTOs seen in some young clusters, but they can make an important contribution to the extended MSTOs in intermediate-age clusters.

### Acceleration in Modified Gravity (MOG) and the Mass-Discrepancy Baryonic Relation [Cross-Listing]

The equation of motion in the generally covariant modified gravity (MOG) theory leads for weak gravitational fields and the non-relativistic limit to a modification of the Newtonian gravitational acceleration law, expressed in terms of two parameters $\alpha$ and $\mu$. The parameter $\alpha$ determines the strength of the gravitational field and $\mu$ is the effective mass of the vector field $\phi_\mu$, coupled with gravitational strength to baryonic matter. The MOG acceleration law for weak field gravitation and non-relativistic particles has been demonstrated to fit a wide range of galaxies, galaxy clusters and the Bullet Cluster and Train Wreck Cluster mergers. We demonstrate that the MOG acceleration law for a point mass source is in agreement with the McGaugh et al., correlation between the radial acceleration traced by galaxy rotation curves and the distribution of baryonic matter for the SPARC sample of 153 rotationally supported spiral and irregular galaxies.

### Acceleration in Modified Gravity (MOG) and the Mass-Discrepancy Baryonic Relation

The equation of motion in the generally covariant modified gravity (MOG) theory leads for weak gravitational fields and the non-relativistic limit to a modification of the Newtonian gravitational acceleration law, expressed in terms of two parameters $\alpha$ and $\mu$. The parameter $\alpha$ determines the strength of the gravitational field and $\mu$ is the effective mass of the vector field $\phi_\mu$, coupled with gravitational strength to baryonic matter. The MOG acceleration law for weak field gravitation and non-relativistic particles has been demonstrated to fit a wide range of galaxies, galaxy clusters and the Bullet Cluster and Train Wreck Cluster mergers. We demonstrate that the MOG acceleration law for a point mass source is in agreement with the McGaugh et al., correlation between the radial acceleration traced by galaxy rotation curves and the distribution of baryonic matter for the SPARC sample of 153 rotationally supported spiral and irregular galaxies.

### Diffuse Interstellar Bands and the Ultraviolet Extinction Curves: The Missing Link Revisited

A large number of interstellar absorption features at ~ 4000\AA\ -- 1.8 {\mu}m, known as the "diffuse interstellar bands" (DIBs), remains unidentified. Most recent works relate them to large polycyclic aromatic hydrocarbon (PAH) molecules or ultrasmall carbonaceous grains which are also thought to be responsible for the 2175 \AA\ extinction bump and/or the far ultraviolet (UV) extinction rise at $\lambda^{-1} > 5.9\ {\mu}m^{-1}$. Therefore, one might expect some relation between the UV extinction and DIBs. Such a relationship, if established, could put important constraints on the carrier of DIBs. Over the past four decades, whether DIBs are related to the shape of the UV extinction curves has been extensively investigated. However, the results are often inconsistent, partly due to the inconsistencies in characterizing the UV extinction. Here we re-examine the connection between the UV extinction curve and DIBs. We compile the extinction curves and the equivalent widths of 40 DIBs along 97 slightlines. We decompose the extinction curve into three Drude-like functions composed of the visible/near-infrared component, the 2175 \AA\ bump, and the far-UV extinction at $\lambda^{-1} > 5.9\ {\mu}m^{-1}$. We argue that the wavelength-integrated far-UV extinction derived from this decomposition technique best measures the strength of the far-UV extinction. No correlation is found between the far-UV extinction and most (~90\%) of the DIBs. We have also shown that the color excess E(1300-1700), the extinction difference at 1300 \AA\ and 1700 \AA\ often used to measure the strength of the far-UV extinction, does not correlate with DIBs. Finally, we confirm the earlier findings of no correlation between the 2175 \AA\ bump and DIBs or between the 2175 \AA\ bump and the far-UV extinction.

### The Megamaser Cosmology Project.IX. Black hole masses for three maser galaxies

As part of the Megamaser Cosmology Project (MCP), we present VLBI maps of nuclear water masers toward five galaxies. The masers originate in sub-parsec circumnuclear disks. For three of the galaxies, we fit Keplerian rotation curves to estimate their supermassive black hole (SMBH) masses, and determine (2.9 $\pm$ 0.3) $\times~10^{6}M_\odot$ for J0437+2456, (1.7 $\pm$ 0.1) $\times~10^{7}M_\odot$ for ESO 558$-$G009, and (1.1 $\pm$ 0.2) $\times~10^{7}M_\odot$ for NGC 5495. In the other two galaxies, Mrk 1029 and NGC 1320, the geometry and dynamics are more complicated and preclude robust black hole mass estimates. Including our new results, we compiled a list of 15 VLBI-confirmed disk maser galaxies with robust SMBH mass measurements. With this sample, we confirm the empirical relation of $R_{out} \propto 0.3 M_{SMBH}$ reported in Wardle & Yusef-Zadeh (2012). We also find a tentative correlation between maser disk outer radii and WISE luminosity. We find no correlations of maser disk size with X-ray 2-10 keV luminosity or [O III] luminosity.

### Two Groups of Red Giants with Distinct Chemical Abundances in the Bulge Globular Cluster NGC 6553 Through the Eyes of APOGEE

Multiple populations revealed in globular clusters (GCs) are important windows to the formation and evolution of these stellar systems. The metal-rich GCs in the Galactic bulge are an indispensable part of this picture, but the high optical extinction in this region has prevented extensive research. In this work, we use the high resolution near-infrared (NIR) spectroscopic data from APOGEE to study the chemical abundances of NGC 6553, which is one of the most metal-rich bulge GCs. We identify ten red giants as cluster members using their positions, radial velocities, iron abundances, and NIR photometry. Our sample stars show a mean radial velocity of $-0.14\pm5.47$ km s$^{-1}$, and a mean [Fe/H] of $-0.15\pm 0.05$. We clearly separate two populations of stars in C and N in this GC for the first time. NGC 6553 is the most metal-rich GC where the multiple stellar population phenomenon is found until now. Substantial chemical variations are also found in Na, O, and Al. However, the two populations show similar Si, Ca, and iron-peak element abundances. Therefore, we infer that the CNO, NeNa, and MgAl cycles have been activated, but the MgAl cycle is too weak to show its effect on Mg. Type Ia and Type II supernovae do not seem to have significantly polluted the second generation stars. Comparing with other GC studies, NGC 6553 shows similar chemical variations as other relatively metal-rich GCs. We also confront current GC formation theories with our results, and suggest possible avenues for improvement in the models.

### Gas and stellar spiral arms and their offsets in the grand-design spiral galaxy M51

Theoretical studies on the response of interstellar gas to a gravitational potential disc with a quasi-stationary spiral arm pattern suggest that the gas experiences a sudden compression due to standing shock waves at spiral arms. This mechanism, called a galactic shock wave, predicts that gas spiral arms move from downstream to upstream of stellar arms with increasing radius inside a corotation radius. In order to investigate if this mechanism is at work in the grand-design spiral galaxy M51, we have measured azimuthal offsets between the peaks of stellar mass and gas mass distributions in its two spiral arms. The stellar mass distribution is created by the spatially resolved spectral energy distribution fitting to optical and near infrared images, while the gas mass distribution is obtained by high-resolution CO and HI data. For the inner region (r < 150"), we find that one arm is consistent with the galactic shock while the other is not. For the outer region, results are less certain due to the narrower range of offset values, the weakness of stellar arms, and the smaller number of successful offset measurements. The results suggest that the nature of two inner spiral arms are different, which is likely due to an interaction with the companion galaxy.

### The Broadband Spectral Variability of Holmberg IX X-1

We present the results from four new broadband X-ray observations of the extreme ultraluminous X-ray source Holmberg IX X-1 ($L_{\rm{X}} > 10^{40}$ erg s$^{-1}$), performed by the $Suzaku$ and $NuSTAR$ observatories in coordination. Combined with the two prior observations coordinating $XMM$-$Newton$, $Suzaku$ and $NuSTAR$, we now have broadband observations of this remarkable source from six separate epochs. Two of these new observations probe lower fluxes than seen previously, allowing us to extend our knowledge of the broadband spectral variability exhibited by Holmberg IX X-1. The broadband spectra are well fit by two thermal blackbody components, which dominate the emission below 10 keV, as well as a steep ($\Gamma \sim 3.5$) powerlaw tail which dominates above $\sim$15 keV. Remarkably, while the 0.3-10.0 keV flux varies by a factor of $\sim$3 between all these epochs, the 15-40 keV flux varies by only $\sim$20%. Although the spectral variability is strongest in the $\sim$1-10 keV band, the broadband variability is not restricted to only one of the thermal components; both are required to vary when all epochs are considered. We also re-visit the search for iron absorption features, leveraging the high-energy $NuSTAR$ data to improve our sensitivity to extreme velocity outflows in light of the ultra-fast outflow recently detected in NGC1313 X-1. We find that iron absorption from a similar outflow along our line of sight can be ruled out in this case. We discuss these results in the context of models for super-Eddington accretion that invoke a funnel-like geometry for the inner flow, and propose a scenario in which we have an almost face-on view of a funnel that expands to larger radii with increasing flux, resulting in an increasing degree of geometrical collimation for the emission from the intermediate temperature regions.

### Structure and Physical Conditions in the Huygens Region of the Orion Nebula

HST images, MUSE maps of emission-lines, and an atlas of high velocity resolution emission-line spectra have been used to establish for the firrst time correlations of the electron temperature, electron density, radial velocity, turbulence, and orientation within the main ionization front of the nebula. From the study of the combined properties of multiple features, it is established that variations in the radial velocity are primarily caused by the photo-evaporating ionization front being viewed at different angles. There is a progressive increase of the electron temperature and density with decreasing distance from the dominant ionizing star Theta1 Ori C. The product of these characteristics (NexTe) is the most relevant parameter in modeling a blister-type nebula like the Huygens Region, where this quantity should vary with the surface brightness in Halpha. Several lines of evidence indicate that small-scale structure and turbulence exists down to the level of our resolution of a few arcseconds. Although photo-evaporative ow must contribute at some level to the well-known non-thermal broadening of the emission lines, comparison of quantitative predictions with the observed optical line widths indicate that it is not the major additive broad- ening component. Derivation of Te values for H+ from radio+optical and optical-only ionized hydro- gen emission showed that this temperature is close to that derived from [Nii] and that the transition from the well-known at extinction curve that applies in the Huygens Region to a more normal steep extinction curve occurs immediately outside of the Bright Bar feature of the nebula.

### Scaling Laws of Passive-Scalar Diffusion in the Interstellar Medium

Passive scalar mixing (metals, molecules, etc.) in the turbulent interstellar medium (ISM) is critical for abundance patterns of stars and clusters, galaxy and star formation, and cooling from the circumgalactic medium. However, the fundamental scaling laws remain poorly understood (and usually unresolved in numerical simulations) in the highly supersonic, magnetized, shearing regime relevant for the ISM.We therefore study the full scaling laws governing passive-scalar transport in idealized simulations of supersonic MHD turbulence, including shear. Using simple phenomenological arguments for the variation of diffusivity with scale based on Richardson diffusion, we propose a simple fractional diffusion equation to describe the turbulent advection of an initial passive scalar distribution. These predictions agree well with the measurements from simulations, and vary with turbulent Mach number in the expected manner, remaining valid even in the presence of a large-scale shear flow (e.g. rotation in a galactic disk). The evolution of the scalar distribution is not the same as obtained using simple, constant "effective diffusivity" as in Smagorinsky models, because the scale-dependence of turbulent transport means an initially Gaussian distribution quickly develops highly non-Gaussian tails. We also emphasize that these are mean scalings that only apply to ensemble behaviors (assuming many different, random scalar injection sites): individual Lagrangian "patches" remain coherent (poorly-mixed) and simply advect for a large number of turbulent flow-crossing times.

### You need to log in to vote

The blog owner requires users to be logged in to be able to vote for this post.

Alternatively, if you do not have an account yet you can create one here.