Learning the fundamental mid-infrared spectral components of galaxies with non-negative matrix factorization

The mid-infrared (MIR) spectra observed with the Spitzer Infrared Spectrograph (IRS) provide a valuable data set for untangling the physical processes and conditions within galaxies. This paper presents the first attempt to blindly learn fundamental spectral components of MIR galaxy spectra, using non-negative matrix factorization (NMF). NMF is a recently developed multivariate technique shown to be successful in blind source separation problems. Unlike the more popular multivariate analysis technique, principal component analysis, NMF imposes the condition that weights and spectral components are non-negative. This more closely resembles the physical process of emission in the MIR, resulting in physically intuitive components. By applying NMF to galaxy spectra in the Cornell Atlas of Spitzer/IRS sources, we find similar components amongst different NMF sets. These similar components include two for active galactic nucleus (AGN) emission and one for star formation. The first AGN component is dominated by fine structure emission lines and hot dust, the second by broad silicate emission at 10 and 18 μm. The star formation component contains all the polycyclic aromatic hydrocarbon features and molecular hydrogen lines. Other components include rising continuums at longer wavelengths, indicative of colder grey-body dust emission. We show an NMF set with seven components can reconstruct the general spectral shape of a wide variety of objects, though struggle to fit the varying strength of emission lines. We also show that the seven components can be used to separate out different types of objects. We model this separation with Gaussian mixtures modelling and use the result to provide a classification tool. We also show that the NMF components can be used to separate out the emission from AGN and star formation regions and define a new star formation/AGN diagnostic which is consistent with all MIR diagnostics already in use but has the advantage that it can be applied to MIR spectra with low signal-to-noise ratio or with limited spectral range. The seven NMF components and code for classification are available at https://github.com/pdh21/NMF_software/

Warm Dust and Spatially Variable PAH Emission in the Dwarf Starburst Galaxy NGC 1705

We present Spitzer observations of the dwarf starburst galaxy NGC 1705 obtained as part of SINGS. The galaxy morphology is very different shortward and longward of ~5 microns: short-wavelength imaging shows an underlying red stellar population, with the central super star cluster (SSC) dominating the luminosity; longer-wavelength data reveals warm dust emission arising from two off-nuclear regions offset by ~250 pc from the SSC. These regions show little extinction at optical wavelengths. The galaxy has a relatively low global dust mass (~2E5 solar masses, implying a global dust-to-gas mass ratio ~2--4 times lower than the Milky Way average). The off-nuclear dust emission appears to be powered by photons from the same stellar population responsible for the excitation of the observed H Alpha emission; these photons are unassociated with the SSC (though a contribution from embedded sources to the IR luminosity of the off-nuclear regions cannot be ruled out). Low-resolution IRS spectroscopy shows moderate-strength PAH emission in the 11.3 micron band in the eastern peak; no PAH emission is detected in the SSC or the western dust emission complex. There is significant diffuse 8 micron emission after scaling and subtracting shorter wavelength data; the spatially variable PAH emission strengths revealed by the IRS data suggest caution in the interpretation of diffuse 8 micron emission as arising from PAH carriers alone. The metallicity of NGC 1705 falls at the transition level of 35% solar found by Engelbracht and collaborators; the fact that a system at this metallicity shows spatially variable PAH emission demonstrates the complexity of interpreting diffuse 8 micron emission. A radio continuum non-detection, NGC 1705 deviates significantly from the canonical far-IR vs. radio correlation. (Abridged)Comment: ApJ, in press; please retrieve full-resolution version from http://www.astro.wesleyan.edu/~cannon/pubs.htm

The Nature of Infrared Emission in the Local Group Dwarf Galaxy NGC 6822 As Revealed by Spitzer

We present Spitzer imaging of the metal-deficient (Z ~30% Z_sun) Local Group dwarf galaxy NGC 6822. On spatial scales of ~130 pc, we study the nature of IR, H alpha, HI, and radio continuum emission. Nebular emission strength correlates with IR surface brightness; however, roughly half of the IR emission is associated with diffuse regions not luminous at H alpha (as found in previous studies). The global ratio of dust to HI gas in the ISM, while uncertain at the factor of ~2 level, is ~25 times lower than the global values derived for spiral galaxies using similar modeling techniques; localized ratios of dust to HI gas are about a factor of five higher than the global value in NGC 6822. There are strong variations (factors of ~10) in the relative ratios of H alpha and IR flux throughout the central disk; the low dust content of NGC 6822 is likely responsible for the different H alpha/IR ratios compared to those found in more metal-rich environments. The H alpha and IR emission is associated with high-column density (> ~1E21 cm^-2) neutral gas. Increases in IR surface brightness appear to be affected by both increased radiation field strength and increased local gas density. Individual regions and the galaxy as a whole fall within the observed scatter of recent high-resolution studies of the radio-far IR correlation in nearby spiral galaxies; this is likely the result of depleted radio and far-IR emission strengths in the ISM of this dwarf galaxy.Comment: ApJ, in press; please retrieve full-resolution version from http://www.astro.wesleyan.edu/~cannon/pubs.htm

Clinicopathologic predictors of finding additional inguinal lymph node metastases in penile cancer patients following positive dynamic sentinel node biopsy:a European multicentre evaluation

OBJECTIVES: Following tumour positive sentinel lymph node biopsy (+DSNB), completion inguinal lymph node dissection (ILND) is negative in 84-89% of basins. Thus, ILND after +DSNB may be considered overtreatment resulting in substantial morbidity. This study aimed to develop a predictive model for additional inguinal lymph node metastases (LNM) at ILND following +DSNB using DSNB characteristics to identify a patient group in which ILND might be omitted. PATIENTS AND METHODS: A retrospective study of 407 inguinal basins with a +DSNB of penile cancer patients who underwent subsequent ILND from seven European centres. From the histopathology reports, the number of positive and negative lymph nodes, extranodal extension and size of the metastasis were recorded. Using bootstrapped logistic regression, variables were selected for the clinical prediction model based on the optimisation of Akaike's information criterion. The area under the curve (AUC) of the receiver operating characteristic (ROC) curve was calculated for the resulting model. Decision curve analysis (DCA) was used to evaluate the clinical utility of the model. RESULTS: 64 (16%) of +DSNB harboured additional LNM at ILND. The number of positive nodes at +DSNB (odds ratio (OR) 2.19; 95% confidence interval (CI) 1.17-4.00; p=0.01) and the largest metastasis size in mm (OR 1.06; 95%CI 1.03-1.10; p=0.001) were selected for the clinical prediction model. The AUC was 0.67 (95%CI 0.60-0.74). The DCA showed no clinical benefit of using the clinical prediction model. CONCLUSION: A small but clinically important group of basins harbour additional LNM at completion ILND following +DSNB. While DSNB characteristics were associated with additional LNM, they did not improve the selection of basins in which ILND could be omitted. Thus, completion ILND remains necessary in all basins with a +DSNB

Wavy Strings: Black or Bright?

Recent developments in string theory have brought forth a considerable interest in time-dependent hair on extended objects. This novel new hair is typically characterized by a wave profile along the horizon and angular momentum quantum numbers $l,m$ in the transverse space. In this work, we present an extensive treatment of such oscillating black objects, focusing on their geometric properties. We first give a theorem of purely geometric nature, stating that such wavy hair cannot be detected by any scalar invariant built out of the curvature and/or matter fields. However, we show that the tidal forces detected by an infalling observer diverge at the `horizon' of a black string superposed with a vibration in any mode with $l \ge 1$. The same argument applied to longitudinal ($l=0$) waves detects only finite tidal forces. We also provide an example with a manifestly smooth metric, proving that at least a certain class of these longitudinal waves have regular horizons.Comment: 45 pages, latex, no figure

Vertical and horizontal distribution of regional new particle formation events in Madrid

The vertical profile of new particle formation (NPF) events was studied by comparing the aerosol size number distributions measured aloft and at surface level in a suburban environment in Madrid, Spain, using airborne instruments. The horizontal distribution and regional impact of the NPF events was investigated with data from three urban, urban background, and suburban stations in the Madrid metropolitan area. Intensive regional NPF episodes followed by particle growth were simultaneously recorded at three stations in and around Madrid during a field campaign in July 2016. The urban stations presented larger formation rates compared to the suburban station. Condensation and coagulation sinks followed a similar evolution at all stations, with higher values at urban stations. However, the total number concentration of particles larger than 2.5 nm was lower at the urban station and peaked around noon, when black carbon (BC) levels are at a minimum. The vertical soundings demonstrated that ultrafine particles (UFPs) are formed exclusively inside the mixed layer. As convection becomes more effective and the mixed layer grows, UFPs are detected at higher levels. The morning soundings revealed the presence of a residual layer in the upper levels in which aged particles (nucleated and grown on previous days) prevail. The particles in this layer also grow in size, with growth rates significantly smaller than those inside the mixed layer. Under conditions with strong enough convection, the soundings revealed homogeneous number size distributions and growth rates at all altitudes, which follow the same evolution at the other stations considered in this study. This indicates that UFPs are detected quasi-homogenously in an area spanning at least 17 km horizontally. The NPF events extend over the full vertical extension of the mixed layer, which can reach as high as 3000 m in the area, according to previous studies. On some days a marked decline in particle size (shrinkage) was observed in the afternoon, associated with a change in air masses. Additionally, a few nocturnal nucleation-mode bursts were observed at the urban stations, for which further research is needed to elucidate their origin.Peer reviewe

Cosmic Ray Protons Accelerated at Cosmological Shocks and Their Impact on Groups and Clusters of Galaxies

We investigate the production of cosmic ray (CR) protons at cosmological shocks by performing, for the first time, numerical simulations of large scale structure formation that include directly the acceleration, transport and energy losses of the high energy particles. CRs are injected at shocks according to the thermal leakage model and, thereafter, accelerated to a power-law distribution as indicated by the test particle limit of the diffusive shock acceleration theory. The evolution of the CR protons accounts for losses due to adiabatic expansion/compression, Coulomb collisions and inelastic p-p scattering. Our results suggest that CR protons produced at shocks formed in association with the process of large scale structure formation could amount to a substantial fraction of the total pressure in the intra-cluster medium. Their presence should be easily revealed by GLAST through detection of gamma-ray flux from the decay of neutral pions produced in inelastic p-p collisions of such CR protons with nuclei of the intra-cluster gas. This measurement will allow a direct determination of the CR pressure contribution in the intra-cluster medium. We also find that the spatial distribution of CR is typically more irregular than that of the thermal gas because it is more influenced by the underlying distribution of shocks. This feature is reflected in the appearance of our gamma-ray synthetic images. Finally, the average CR pressure distribution appears statistically slightly more extended than the thermal pressure.Comment: 16 pages, 7 figures, ApJ in pres