Optical evidence for the unification of active galactic nuclei and quasi-stellar objects.

There is a variety of optical evidence for some unification of different types of active galactic nuclei and quasi-stellar objects (QSOs). The case is very strong for the unification of at least some Seyfert galaxies, where polarization data show that the type assigned to the Seyfert galaxy must depend on viewing direction. It has been proposed that Fanaroff-Riley type 2 (FR2) radio galaxies are quasars seen in a direction from which the quasar is obscured, and there is some limited direct evidence for this picture. The broad absorption line QSOs may be normal QSOs seen from a special direction. Some of the sources observed to have high luminosities in the far infrared could be obscured QSOs and active nuclei. Mergers and interactions are likely to play an important role in nuclear activity, and active galaxies and QSOs could change their apparent types through these encounters followed by subsequent evolution.

Probing active galactic nuclei with H2O megamasers.

We describe the characteristics of the rapidly rotating molecular disk in the nucleus of the mildly active galaxy NGC4258. The morphology and kinematics of the disk are delineated by the point-like watervapor emission sources at 1.35-cm wavelength. High angular resolution [200 microas where as is arcsec, corresponding to 0.006 parsec (pc) at 6.4 million pc] and high spectral resolution (0.2 km.s-1 or nu/Deltanu = 1.4 x 10(6)) with the Very-Long-Baseline Array allow precise definition of the disk. The disk is very thin, but slightly warped, and is viewed nearly edge-on. The masers show that the disk is in nearly perfect Keplerian rotation within the observable range of radii of 0.13-0.26 pc. The approximately random deviations from the Keplerian rotation curve among the high-velocity masers are approximately 3.5 km.s-1 (rms). These deviations may be due to the masers lying off the midline by about +/-4 degrees or variations in the inclination of the disk by +/-4 degrees. Lack of systematic deviations indicates that the disk has a mass of <4 x 10(6) solar mass (M[symbol: see text]). The gravitational binding mass is 3.5 x 10(7) M[symbol: see text], which must lie within the inner radius of the disk and requires that the mass density be >4 x 10(9) M[symbol: see text].pc-3. If the central mass were in the form of a star cluster with a density distribution such as a Plummer model, then the central mass density would be 4 x 10(12) M[symbol: see text].pc-3. The lifetime of such a cluster would be short with respect to the age of the galaxy [Maoz, E. (1995) Astrophys. J. Lett. 447, L91-L94]. Therefore, the central mass may be a black hole. The disk as traced by the systemic velocity features is unresolved in the vertical direction, indicating that its scale height is <0.0003 pc (hence the ratio of thickness to radius, H/R, is <0.0025). For a disk in hydrostatic equilibrium the quadrature sum of the sound speed and Alfven velocity is <2.5 km.s-1, so that the temperature of the disk must be <1000 K and the toroidal magnetic field component must be <250 mG. If the molecular mass density in the disk is 10(10) cm-3, then the disk mass is approximately 10(4) M[symbol: see text], and the disk is marginally stable as defined by the Toomre stability parameter Q (Q = 6 at the inner edge and 1 at the outer edge). The inward drift velocity is predicted to be <0.007 km.s-1, for a viscosity parameter of 0.1, and the accretion rate is <7 x 10(-5) M[symbol: see text].yr-1. At this value the accretion would be sufficient to power the nuclear x-ray source of 4 x 10(40) ergs-1 (1 erg = 0.1 microJ). The volume of individual maser components may be as large as 10(46) cm3, based on the velocity gradients, which is sufficient to supply the observed luminosity. The pump power undoubtedly comes from the nucleus, perhaps in the form of x-rays. The warp may allow the pump radiation to penetrate the disk obliquely [Neufeld, D. A. & Maloney, P. R. (1995) Astrophys. J. Lett. 447, L17-L19]. A total of 15 H2O megamasers have been identified out of >250 galaxies searched. Galaxy NGC4258 may be the only case where conditions are optimal to reveal a well-defined nuclear disk. Future measurement of proper motions and accelerations for NGC4258 will yield an accurate distance and a more precise definition of the dynamics of the disk

Mass-metallicity relation explored with CALIFA I. Is there a dependence on the star-formation rate?

We studied the global and local ℳ-Z relation based on the first data available from the CALIFA survey (150 galaxies). This survey provides integral field spectroscopy of the complete optical extent of each galaxy (up to 2−3 effective radii), with a resolution high enough to separate individual H II regions and/or aggregations. About 3000 individual H II regions have been detected. The spectra cover the wavelength range between [OII]3727 and [SII]6731, with a sufficient signal-to-noise ratio to derive the oxygen abundance and star-formation rate associated with each region. In addition, we computed the integrated and spatially resolved stellar masses (and surface densities) based on SDSS photometric data. We explore the relations between the stellar mass, oxygen abundance and star-formation rate using this dataset. We derive a tight relation between the integrated stellar mass and the gas-phase abundance, with a dispersion lower than the one already reported in the literature (σ_Δlog (O/H) = 0.07 dex). Indeed, this dispersion is only slightly higher than the typical error derived for our oxygen abundances. However, we found no secondary relation with the star-formation rate other than the one induced by the primary relation of this quantity with the stellar mass. The analysis for our sample of ~3000 individual H II   regions confirms (i) a local mass-metallicity relation and (ii) the lack of a secondary relation with the star-formation rate. The same analysis was performed with similar results for the specific star-formation rate. Our results agree with the scenario in which gas recycling in galaxies, both locally and globally, is much faster than other typical timescales, such like that of gas accretion by inflow and/or metal loss due to outflows. In essence, late-type/disk-dominated galaxies seem to be in a quasi-steady situation, with a behavior similar to the one expected from an instantaneous recycling/closed-box model.

DUst around NEarby Stars. The survey observational results

Context. Debris discs are a consequence of the planet formation process and constitute the fingerprints of planetesimal systems. Their solar system counterparts are the asteroid and Edgeworth-Kuiper belts. Aims. The DUNES survey aims at detecting extra-solar analogues to the Edgeworth-Kuiper belt around solar-type stars, putting in this way the solar system into context. The survey allows us to address some questions related to the prevalence and properties of planetesimal systems. Methods. We used Herschel/PACS to observe a sample of nearby FGK stars. Data at 100 and 160 μm were obtained, complemented in some cases with observations at 70 μm, and at 250, 350 and 500 μm using SPIRE. The observing strategy was to integrate as deep as possible at 100 μm to detect the stellar photosphere. Results. Debris discs have been detected at a fractional luminosity level down to several times that of the Edgeworth-Kuiper belt. The incidence rate of discs around the DUNES stars is increased from a rate of ~12.1% ± 5% before Herschel to ~20.2% ± 2%. A significant fraction (~52%) of the discs are resolved, which represents an enormous step ahead from the previously known resolved discs. Some stars are associated with faint far-IR excesses attributed to a new class of cold discs. Although it cannot be excluded that these excesses are produced by coincidental alignment of background galaxies, statistical arguments suggest that at least some of them are true debris discs. Some discs display peculiar SEDs with spectral indexes in the 70–160 μm range steeper than the Rayleigh-Jeans one. An analysis of the debris disc parameters suggests that a decrease might exist of the mean black body radius from the F-type to the K-type stars. In addition, a weak trend is suggested for a correlation of disc sizes and an anticorrelation of disc temperatures with the stellar age.

Pipe3D, a pipeline to analyze integral field spectroscopy DATA: II. Analysis sequence and CALIFA dataproducts

We present PIPE3D, an analysis pipeline based on the FIT3D fitting tool, developed to explore the properties of the stellar populations and ionized gas of integral field spectroscopy (IFS) data. PIPE3D was created to provide coherent, simple to distribute, and comparable dataproducts, independently of the origin of the data, focused on the data of the most recent IFU surveys (e.g., CALIFA, MaNGA, and SAMI), and the last generation IFS instruments (e.g., MUSE). In this article we describe the different steps involved in the analysis of the data, illustrating them by showing the dataproducts derived for NGC 2916, observed by CALIFA and P-MaNGA. As a practical example of the pipeline we present the complete set of dataproducts derived for the 200 datacubes that comprises the V500 setup of the CALIFA Data Release 2 (DR2), making them freely available through the network. Finally, we explore the hypothesis that the properties of the stellar populations and ionized gas of galaxies at the effective radius are representative of the overall average ones, finding that this is indeed the case.

Pipe3D, a pipeline to analyse integral Field Spectroscopy data: I. New fitting philosophy of FIT3D

We present an improved version of FIT3D, a fitting tool for the analysis of the spectroscopic properties of the stellar populations and the ionized gas derived from moderate resolution spectra of galaxies. This tool was developed to analyze integral field spectroscopy data and it is the basis of PIPE3D, a pipeline used in the analysis of CALIFA, MaNGA, and SAMI data. We describe the philosophy and each step of the fitting procedure. We present an extensive set of simulations in order to estimate the precision and accuracy of the derived parameters for the stellar populations and the ionized gas. We report on the results of those simulations. Finally, we compare the results of the analysis using FIT3D with those provided by other widely used packages, and we find that the parameters derived by FIT3D are fully compatible with those derived using these other tools.

Low-resolution spectroscopy and spectral energy distributions of selected sources towards σ Orionis

Aims. We study in detail nine sources in the direction of the young σ Orionis cluster, which is considered to be a unique site for studying stellar and substellar formation. The nine sources were selected because of their peculiar properties, such as extremely-red infrared colours or excessively strong Hα emission for their blue optical colours. Methods. We acquired high-quality, low-resolution spectroscopy (R ∼ 500) of the nine targets with ALFOSC at the Nordic Optical Telescope. We also re-analysed [24]-band photometry from MIPS/Spitzer and compiled the highest quality photometric dataset available at the ViJHK_s passbands and the four IRAC/Spitzer channels, for constructing accurate spectral energy distributions between 0.55 and 24 μm. Results. The nine targets were classified into: one Herbig Ae/Be star with a scattering edge-on disc; two G-type stars; one X-ray flaring, early-M, young star with chromospheric Hα emission; one very low-mass, accreting, young spectroscopic binary; two young objects at the brown-dwarf boundary with the characteristics of classical T Tauri stars; and two emission-line galaxies, one undergoing star formation, and another whose spectral energy distribution is dominated by an active galactic nucleus. We also discovered three infrared sources associated with overdensities in a cold cloud of the cluster centre. Conclusions. Low-resolution spectroscopy and spectral energy distributions are a vital tool for measuring the physical properties and evolution of young stars and candidates in the σ Orionis cluster.

Caught in the act: gas and stellar velocity dispersions in a fast Quenching compact star-forming galaxy at z ~ 1.7

We present Keck I MOSFIRE spectroscopy in the Y and H bands of GDN-8231, a massive, compact, star-forming galaxy at a redshift of z ~ 1.7. Its spectrum reveals both Hα and [Nii] emission lines and strong Balmer absorption lines. The Hα and Spitzer MIPS 24 μm fluxes are both weak, thus indicating a low star-formation rate of SFR≲5-10 M_⨀ yr−1. This, added to a relatively young age of ~700 Myr measured from the absorption lines, provides the first direct evidence for a distant galaxy being caught in the act of rapidly shutting down its star formation. Such quenching allows GDN-8231 to become a compact, quiescent galaxy, similar to three other galaxies in our sample, by z ~ 1.5. Moreover, the color profile of GDN-8231 shows a bluer center, consistent with the predictions of recent simulations for an early phase of inside-out quenching. Its line-of-sight velocity dispersion for the gas, σ_LOG^gas = 127 ± 32 km s^−1, is nearly 40% smaller than that of its stars, σ_LOG^* = 215 ± 35 km s^−1. High-resolution hydro-simulations of galaxies explain such apparently colder gas kinematics of up to a factor of ~1.5 with rotating disks being viewed at different inclinations and/or centrally concentrated star-forming regions. A clear prediction is that their compact, quiescent descendants preserve some remnant rotation from their star-forming progenitors.

A search for anisotropy in the arrival directions of ultra high energy cosmic rays recorded at the Pierre Auger Observatory

Observations of cosmic rays arrival directions made with the Pierre Auger Observatory have previously provided evidence of anisotropy at the 99% CL using the correlation of ultra high energy cosmic rays (UHECRs) with objects drawn from the Veron-Cetty Veron catalog. In this paper we report on the use of three catalog independent methods to search for anisotropy. The 2pt-L, 2pt+ and 3pt methods, each giving a different measure of self-clustering in arrival directions, were tested on mock cosmic ray data sets to study the impacts of sample size and magnetic smearing on their results, accounting for both angular and energy resolutions. If the sources of UHECRs follow the same large scale structure as ordinary galaxies in the local Universe and if UHECRs are deflected no more than a few degrees, a study of mock maps suggests that these three method can efficiently respond to the resulting anisotropy with a P-value = 1.0% or smaller with data sets as few as 100 events. using data taken from January 1, 2004 to July 31, 2010 we examined the 20, 30, ... , 110 highest energy events with a corresponding minimum energy threshold of about 49.3 EeV. The minimum P-values found were 13.5% using the 2pt-L method, 1.0% using the 2pt+ method and 1.1% using the 3pt method for the highest 100 energy events. In view of the multiple (correlated) scans performed on the data set, these catalog-independent methods do not yield strong evidence of anisotropy in the highest energy cosmic rays.

Episodic star formation in a group of LAEs at z=5.07

We are undertaking a search for high-redshift low-luminosity Lyman Alpha sources in the SHARDS (Survey for High-z Absorption Red and Dead Sources) survey. Among the pre-selected Lyman Alpha sources two candidates were spotted, located 3.19 arcsec apart, and tentatively at the same redshift. Here, we report on the spectroscopic confirmation with Gran Telescopio Canarias of the Lyman Alpha emission from this pair of galaxies at a confirmed spectroscopic redshifts of z=5.07. Furthermore, one of the sources is interacting/merging with another close companion that looks distorted. Based on the analysis of the spectroscopy and additional photometric data, we infer that most of the stellar mass of these objects was assembled in a burst of star formation 100 Myr ago. A more recent burst (2 Myr old) is necessary to account for the measured Lyman Alpha flux. We claim that these two galaxies are good examples of Lyman Alpha sources undergoing episodic star formation. Besides, these sources very likely constitute a group of interacting Lyman Alpha emitters (LAEs).

Investigating evidence for different black hole accretion modes since redshift z ∼ 1

Chandra data in the COSMOS, AEGIS-XD and 4 Ms Chandra Deep Field South are combined with multiwavelength photometry available in those fields to determine the rest-frame U − V versus V − J colours of X-ray AGN hosts in the redshift intervals 0.1 < z < 0.6 (mean z¯=0.40) and 0.6 < z < 1.2 (mean z¯=0.85). This combination of colours provides an effective and least model-dependent means of separating quiescent from star-forming, including dust reddened, galaxies. Morphological information emphasizes differences between AGN populations split by their U − V versus V − J colours. AGN in quiescent galaxies consist almost exclusively of bulges, while star-forming hosts are equally split between early- and late-type hosts. The position of AGN hosts on the U − V versusV − J diagram is then used to set limits on the accretion density of the Universe associated with evolved and star-forming systems independent of dust induced biases. It is found that most of the black hole growth at z≈ 0.40 and 0.85 is associated with star-forming hosts. Nevertheless, a non-negligible fraction of the X-ray luminosity density, about 15–20 per cent, at both z¯=0.40 and 0.85, is taking place in galaxies in the quiescent region of the U − V versus V − J diagram. For the low-redshift sub-sample, 0.1 < z < 0.6, we also find tentative evidence, significant at the 2σ level, that AGN split by their U − V and V − J colours have different Eddington ratio distributions. AGN in blue star-forming hosts dominate at relatively high Eddington ratios. In contrast, AGN in red quiescent hosts become increasingly important as a fraction of the total population towards low Eddington ratios. At higher redshift, z > 0.6, such differences are significant at the 2σ level only for sources with Eddington ratios ≳ 10^− 3. These findings are consistent with scenarios in which diverse accretion modes are responsible for the build-up of supermassive black holes at the centres of galaxies. We compare these results with the predictions of theGALFORM semi-analytic model for the cosmological evolution of AGN and galaxies. This model postulates two black hole fuelling modes, the first is linked to star formation events and the second takes place in passive galaxies. GALFORM predicts that a substantial fraction of the black hole growth at z < 1 is associated with quiescent galaxies, in apparent conflict with the observations. Relaxing the strong assumption of the model that passive AGN hosts have zero star formation rate could bring those predictions in better agreement with the data.

[C II] and 12^CO(1-0) emission maps in HLSJ091828.6+514223: a strongly lensed interacting system AT z = 5.24

We present Submillimeter Array [C II] 158 μm and Karl G. Jansky Very Large Array 12^CO(1-0) line emission maps for the bright, lensed, submillimeter source at z = 5.2430 behind A 773: HLSJ091828.6+514223 (HLS0918). We combine these measurements with previously reported line profiles, including multiple 12^CO rotational transitions, [C I], water, and [N II], providing some of the best constraints on the properties of the interstellar medium in a galaxy at z > 5. HLS0918 has a total far-infrared (FIR) luminosity L_FIR(8–1000 μm) = (1.6 ± 0.1) × 10^14 L_☉ μ^–1, where the total magnification μ_total = 8.9 ± 1.9, via a new lens model from the [C II] and continuum maps. Despite a HyLIRG luminosity, the FIR continuum shape resembles that of a local LIRG. We simultaneously fit all of the observed spectral line profiles, finding four components that correspond cleanly to discrete spatial structures identified in the maps. The two most redshifted spectral components occupy the nucleus of a massive galaxy, with a source-plane separation <1 kpc. The reddest dominates the continuum map (demagnified L_FIR, component = (1.1 ± 0.2) × 10^13 L_☉) and excites strong water emission in both nuclear components via a powerful FIR radiation field from the intense star formation. A third star-forming component is most likely a region of a merging companion (ΔV ~ 500 km s^–1) exhibiting generally similar gas properties. The bluest component originates from a spatially distinct region and photodissociation region analysis suggests that it is lower density, cooler, and forming stars less vigorously than the other components. Strikingly, it has very strong [N II] emission, which may suggest an ionized, molecular outflow. This comprehensive view of gas properties and morphology in HLS0918 previews the science possible for a large sample of high-redshift galaxies once ALMA attains full sensitivity.

An extended Herschel drop-out source in the center of AS1063: a normal dusty galaxy at z=6.1 or SZ substructures?

In the course of our 870 μm APEX/LABOCA follow-up of the Herschel Lensing Survey we have detected a source in AS1063 (RXC J2248.7-4431) that has no counterparts in any of the Herschel PACS/SPIRE bands, it is a Herschel “drop-out” with S_870/S_500 ≥ 0.5. The 870 μm emission is extended and centered on the brightest cluster galaxy, suggesting either a multiply imaged background source or substructure in the Sunyaev-Zel’dovich increment due to inhomogeneities in the hot cluster gas of this merging cluster. We discuss both interpretations with emphasis on the putative lensed source. Based on the observed properties and on our lens model we find that this source may be the first submillimeter galaxy (SMG) with a moderate far-infrared (FIR) luminosity (L_FIR < 10^12 L_⊙) detected so far at z > 4. In deep HST observations we identified a multiply imaged z ~ 6 source and measured its spectroscopic redshift to be z = 6.107 with VLT/FORS. This source may be associated with the putative SMG, but it is most likely offset spatially by 10−30 kpc and they may be interacting galaxies. With a FIR luminosity in the range [5−15] × 10^11 L_⊙ corresponding to a star formation rate in the range [80−260] M_⊙ yr^-1, this SMG would be more representative of the z > 4 dusty galaxies than the extreme starbursts detected so far. With a total magnification of ~25 it would open a unique window to the normal dusty galaxies at the end of the epoch of reionization.

The merger history of massive spheroids since z ∼ 1 is size independent

Using a compilation of 379 massive (stellar mass M ≳ 10^11 M_⊙) spheroid-like galaxies from the near-infrared Palomar/DEEP-2 survey, we investigated, up to z ∼ 1, whether the presence of companions depends on the size of the host galaxy. We explored the presence of companions for mass ratios with respect to the central massive galaxy down to 1 : 10 and 1 : 100, and within projected distances of 30, 50 and 100 kpc of these objects. We found evidence that these companions are equally distributed around both compact and extended massive spheroid-like galaxies. This suggests that, at least since z ∼ 1, the merger activity in these objects is nearly homogeneous across the whole population and that the merger history is not affected by the size of the host galaxy. Our results could indicate that compact and extended massive spheroid-like galaxies are increasing in size at the same rate.

Deep observation of the NGC1275 region with MAGIC: search of diffuse gamma-ray emission from cosmic rays in the Perseus cluster

Clusters of galaxies are expected to be reservoirs of cosmic rays (CRs) that should produce diffuse γ-ray emission due to their hadronic interactions with the intra-cluster medium. The nearby Perseus cool-core cluster, identified as the most promising target to search for such an emission, has been observed with the MAGIC telescopes at very-high energies (VHE, E ≥ 100 GeV) for a total of 253 hr from 2009 to 2014. The active nuclei of NGC 1275, the central dominant galaxy of the cluster, and IC 310, lying at about 0.6º from the centre, have been detected as point-like VHE γ-ray emitters during the first phase of this campaign. We report an updated measurement of the NGC 1275 spectrum, which is described well by a power law with a photon index Γ = 3.6 ± 0.2_(stat) ± 0.2_(syst) between 90 GeV and 1200 GeV. We do not detect any diffuse γ-ray emission from the cluster and so set stringent constraints on its CR population. To bracket the uncertainties over the CR spatial and spectral distributions, we adopt different spatial templates and power-law spectral indexes α. For α = 2.2, the CR-to-thermal pressure within the cluster virial radius is constrained to be ≤ 1 − 2%, except if CRs can propagate out of the cluster core, generating a flatter radial distribution and releasing the CR-to-thermal pressure constraint to ≤ 20%. Assuming that the observed radio mini-halo of Perseus is generated by secondary electrons from CR hadronic interactions, we can derive lower limits on the central magnetic field, B_(0), that depend on the CR distribution. For α = 2.2, B_(0) ≥ 5 − 8 µG, which is below the ∼25 µG inferred from Faraday rotation measurements, whereas for α ≤ 2.1, the hadronic interpretation of the diffuse radio emission contrasts with our γ-ray flux upper limits independently of the magnetic field strength.