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

The Mice at play in the CALIFA survey. A case study of a gas-rich major merger between first passage and coalescence

We present optical integral field spectroscopy (IFS) observations of the Mice, a major merger between two massive (≳10^11 M_⊙) gas-rich spirals NGC 4676A and B, observed between first passage and final coalescence. The spectra provide stellar and gas kinematics, ionised gas properties, and stellar population diagnostics, over the full optical extent of both galaxies with ~1.6 kpc spatial resolution. The Mice galaxies provide a perfect case study that highlights the importance of IFS data for improving our understanding of local galaxies. The impact of first passage on the kinematics of the stars and gas has been significant, with strong bars most likely induced in both galaxies. The barred spiral NGC 4676B exhibits a strong twist in both its stellar and ionised gas disk. The edge-on disk galaxy NGC 4676A appears to be bulge free, with a strong bar causing its “boxy” light profile. On the other hand, the impact of the merger on the stellar populations has been minimal thus far. By combining the IFS data with archival multiwavelength observations we show that star formation induced by the recent close passage has not contributed significantly to the total star formation rate or stellar mass of the galaxies. Both galaxies show bicones of high ionisation gas extending along their minor axes. In NGC 4676A the high gas velocity dispersion and Seyfert-like line ratios at large scaleheight indicate a powerful outflow. Fast shocks (vs ~ 350 km s^-1) extend to ~6.6 kpc above the disk plane. The measured ram pressure (P/k = 4.8 × 10^6 K cm^-3) and mass outflow rate (~8−20 M_⊙ yr^-1) are similar to superwinds from local ultra-luminous infrared galaxies, although NGC 4676A only has a moderate infrared luminosity of 3 × 10^10 L_⊙. Energy beyond what is provided by the mechanical energy of the starburst appears to be required to drive the outflow. Finally, we compare the observations to mock kinematic and stellar population maps extracted from a hydrodynamical merger simulation. The models show little enhancement in star formation during and following first passage, in agreement with the observations. We highlight areas where IFS data could help further constrain the models.

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.

No surviving evolved companions of the progenitor of SN 1006

Type Ia supernovae are thought to occur when a white dwarf made of carbon and oxygen accretes sufficient mass to trigger a thermonuclear explosion(1). The accretion could be slow, from an unevolved (main-sequence) or evolved (subgiant or giant) star(2,3) (the single-degenerate channel), or rapid, as the primary star breaks up a smaller orbiting white dwarf(3,4) (the double-degenerate channel). A companion star will survive the explosion only in the single-degenerate channel(5). Both channels might contribute to the production of type Ia supernovae(6,7), but the relative proportions of their contributions remain a fundamental puzzle in astronomy. Previous searches for remnant companions have revealed one possible case for SN 1572 (refs 8, 9), although that has been questioned(10). More recently, observations have restricted surviving companions to be small, main-sequence stars(11-13), ruling out giant companions but still allowing the single-degenerate channel. Here we report the results of a search for surviving companions of the progenitor of SN 1006 (ref. 14). None of the stars within 4 arc minutes of the apparent site of the explosion is associated with the supernova remnant, and we can firmly exclude all giant and subgiant stars from being companions of the progenitor. In combination with previous results, our findings indicate that fewer than 20 per cent of type Ia supernovae occur through the single-degenerate channel.

Effect of magnetic activity saturation in chromospheric flux-flux relationships

We present a homogeneous study of chromospheric and coronal flux–flux relationships using a sample of 298 late-type dwarf active stars with spectral types F to M. The chromospheric lines were observed simultaneously in each star to avoid spread as a result of long-term variability. Unlike other works, we subtract the basal chromospheric contribution in all the spectral lines studied. For the first time, we quantify the departure of dMe stars from the general relations. We show that dK and dKe stars also deviate from the general trend. Studying the flux–colour diagrams, we demonstrate that the stars deviating from the general relations are those with saturated X-ray emission and we show that these stars also present saturation in the Hα line. Using several age spectral indicators, we show that these are younger stars than those following the general relationships. The non-universality of flux–flux relationships found in this work should be taken into account when converting between fluxes in different chromospheric activity indicators.

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.

Stars and brown dwarfs in the σ Orionis cluster (Research Note) III. OSIRIS/GTC low-resolution spectroscopy of variable sources

Context. Although many studies have been performed so far, there are still dozens of low-mass stars and brown dwarfs in the young σ Orionis open cluster without detailed spectroscopic characterisation. Aims. We look for unknown strong accretors and disc hosts that were undetected in previous surveys. Methods. We collected low-resolution spectroscopy (R ~ 700) of ten low-mass stars and brown dwarfs in σ Orionis with OSIRIS at the Gran Telescopio Canarias under very poor weather conditions. These objects display variability in the optical, infrared, Hα, and/or X-rays on time scales of hours to years. We complemented our spectra with optical and near-/mid-infrared photometry. Results. For seven targets, we detected lithium in absorption, identified Hα, the calcium doublet, and forbidden lines in emission, and/or determined spectral types for the first time. We characterise in detail a faint, T Tauri-like brown dwarf with an 18 h-period variability in the optical and a large Hα equivalent width of –125  ±  15 Å, as well as two M1-type, X-ray-flaring, low-mass stars, one with a warm disc and forbidden emission lines, the other with a previously unknown cold disc with a large inner hole. Conclusions. New unrevealed strong accretors and disc hosts, even below the substellar limit, await discovery among the list of known σ Orionis stars and brown dwarfs that are variable in the optical and have no detailed spectroscopic characterisation yet.

Emission-line galaxy surveys with GTC: The Hα Universe

Emission line galaxies are the most easily detected and studied objects in the high redshift Universe. They are being used to trace the evolution of critical observables of the Universe such as Star Formation Rate densities, starburst properties and abundances. Most of the research is being done using [OII]3727 and UV lines, but Hα is still one the best tracers for Star Formation Rate and physical properties of current star-forming galaxies. As a complementary contribution to studies of galaxy evolution, our team has focused into a long-term project to study the population of Hα-selected star-forming galaxies of the Universe at different redshifts. In 1995 we first determined the local Hα luminosity function, and from it the Star Formation Rate density (SFRd) of the local Universe. We then, using narrow-band imaging in the optical, extended this measurement to z ≃ 0.24 and z ≃ 0.4. Working in the near-infrared, GTC will be a very powerful tool to study the evolution of the Hα emission-line galaxy populations at different redshifts. We will both quantify the SFRd evolution and characterize the star-forming galaxy populations by directly comparing the same observables at all redshifts up to z ≃ 2.5.

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.