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.

The evolutionary paths among galaxy types on the red sequence at 0.3 < Z < 1.5

We have studied the main evolutionary paths among the galaxy types residing on the massive end of the Red Sequence and nearby locations on the Green Valley during the last ∼9 Gyr. The morphological and star formation properties of a sample of these galaxies at 0 . 3 < z < 1 .5 with stellar masses M_∗ > 5 × 10^10 M_⊙ have been analysed. We present direct observational evidence for the first time of the existence of two main evolutionary paths among the different red galaxy types since z ∼ 1 .5, which provide some clues on the nature of the processes that have governed the assembly of present-day massive quiescent galaxies. The results are in excellent agreement with the hierarchical evolutionary framework proposed in the Eliche-Moral et al. (2010) model. Data from SHARDS (one of the ESO/GTC Large Programmes approved in 2009A) will complement and improve the present findings, shedding some light into many of the still unsettled questions concerning the migration of galaxies from the Blue Cloud to the Red Sequence at z < 1 .5.

The circumstellar environment and evolutionary status of the supergiant B[e] star Wd1-9

Context. Historically, supergiant (sg)B[e] stars have been difficult to include in theoretical schemes for the evolution of massive OB stars. Aims. The location of Wd1-9 within the coeval starburst cluster Westerlund 1 means that it may be placed into a proper evolutionary context and we therefore aim to utilise a comprehensive multiwavelength dataset to determine its physical properties and consequently its relation to other sgB[e] stars and the global population of massive evolved stars within Wd1. Methods. Multi-epoch R- and I-band VLT/UVES and VLT/FORS2 spectra are used to constrain the properties of the circumstellar gas, while an ISO-SWS spectrum covering 2.45−45μm is used to investigate the distribution, geometry and composition of the dust via a semi-analytic irradiated disk model. Radio emission enables a long term mass-loss history to be determined, while X-ray observations reveal the physical nature of high energy processes within the system. Results. Wd1-9 exhibits the rich optical emission line spectrum that is characteristic of sgB[e] stars. Likewise its mid-IR spectrum resembles those of the LMC sgB[e] stars R66 and 126, revealing the presence of equatorially concentrated silicate dust, with a mass of ~10−4M⊙. Extreme historical and ongoing mass loss (≳ 10−4M⊙yr−1) is inferred from the radio observations. The X-ray properties of Wd1-9 imply the presence of high temperature plasma within the system and are directly comparable to a number of confirmed short-period colliding wind binaries within Wd1. Conclusions. The most complete explanation for the observational properties of Wd1-9 is that it is a massive interacting binary currently undergoing, or recently exited from, rapid Roche-lobe overflow, supporting the hypothesis that binarity mediates the formation of (a subset of) sgB[e] stars. The mass loss rate of Wd1-9 is consistent with such an assertion, while viable progenitor and descendent systems are present within Wd1 and comparable sgB[e] binaries have been identified in the Galaxy. Moreover, the rarity of sgB[e] stars - only two examples are identified from a census of ~ 68 young massive Galactic clusters and associations containing ~ 600 post-Main Sequence stars - is explicable given the rapidity (~ 104yr) expected for this phase of massive binary evolution.

The Gaia-ESO Survey: processing FLAMES-UVES spectra

The Gaia-ESO Survey is a large public spectroscopic survey that aims to derive radial velocities and fundamental parameters of about 105 Milky Way stars in the field and in clusters. Observations are carried out with the multi-object optical spectrograph FLAMES, using simultaneously the medium-resolution (R ~ 20 000) GIRAFFE spectrograph and the high-resolution (R ~ 47 000) UVES spectrograph. In this paper we describe the methods and the software used for the data reduction, the derivation of the radial velocities, and the quality control of the FLAMES-UVES spectra. Data reduction has been performed using a workflow specifically developed for this project. This workflow runs the ESO public pipeline optimizing the data reduction for the Gaia-ESO Survey, automatically performs sky subtraction, barycentric correction and normalisation, and calculates radial velocities and a first guess of the rotational velocities. The quality control is performed using the output parameters from the ESO pipeline, by a visual inspection of the spectra and by the analysis of the signal-to-noise ratio of the spectra. Using the observations of the first 18 months, specifically targets observed multiple times at different epochs, stars observed with both GIRAFFE and UVES, and observations of radial velocity standards, we estimated the precision and the accuracy of the radial velocities. The statistical error on the radial velocities is σ ~ 0.4 km s-1 and is mainly due to uncertainties in the zero point of the wavelength calibration. However, we found a systematic bias with respect to the GIRAFFE spectra (~0.9 km s-1) and to the radial velocities of the standard stars (~0.5 km s-1) retrieved from the literature. This bias will be corrected in the future data releases, when a common zero point for all the set-ups and instruments used for the survey is be established.

A VLT/FLAMES survey for massive binaries in Westerlund 1. IV. Wd1-5 – binary product and a pre-supernova companion for the magnetar CXOU J1647-45?

Context. The first soft gamma-ray repeater was discovered over three decades ago, and was subsequently identified as a magnetar, a class of highly magnetised neutron star. It has been hypothesised that these stars power some of the brightest supernovae known, and that they may form the central engines of some long duration gamma-ray bursts. However there is currently no consenus on the formation channel(s) of these objects. Aims. The presence of a magnetar in the starburst cluster Westerlund 1 implies a progenitor with a mass ≥40 M⊙, which favours its formation in a binary that was disrupted at supernova. To test this hypothesis we conducted a search for the putative pre-SN companion. Methods. This was accomplished via a radial velocity survey to identify high-velocity runaways, with subsequent non-LTE model atmosphere analysis of the resultant candidate, Wd1-5. Results. Wd1-5 closely resembles the primaries in the short-period binaries, Wd1-13 and 44, suggesting a similar evolutionary history, although it currently appears single. It is overluminous for its spectroscopic mass and we find evidence of He- and N-enrichement, O-depletion, and critically C-enrichment, a combination of properties that is difficult to explain under single star evolutionary paradigms. We infer a pre-SN history for Wd1-5 which supposes an initial close binary comprising two stars of comparable (~ 41 M⊙ + 35 M⊙) masses. Efficient mass transfer from the initially more massive component leads to the mass-gainer evolving more rapidly, initiating luminous blue variable/common envelope evolution. Reverse, wind-driven mass transfer during its subsequent WC Wolf-Rayet phase leads to the carbon pollution of Wd1-5, before a type Ibc supernova disrupts the binary system. Under the assumption of a physical association between Wd1-5 and J1647-45, the secondary is identified as the magnetar progenitor; its common envelope evolutionary phase prevents spin-down of its core prior to SN and the seed magnetic field for the magnetar forms either in this phase or during the earlier episode of mass transfer in which it was spun-up. Conclusions. Our results suggest that binarity is a key ingredient in the formation of at least a subset of magnetars by preventing spin-down via core-coupling and potentially generating a seed magnetic field. The apparent formation of a magnetar in a Type Ibc supernova is consistent with recent suggestions that superluminous Type Ibc supernovae are powered by the rapid spin-down of these objects.

ATCA HI observations of the peculiar galaxy IC 2554

ATCA H I and radio continuum observations of the peculiar southern galaxy IC 2554 and its surroundings reveal typical signatures of an interacting galaxy group. We detected a large H I cloud between IC 2554 and the elliptical galaxy NGC 3136B. The gas dynamics in IC 2554 itself, which is sometimes described as a colliding pair, are surprisingly regular, whereas NGC 3136B was not detected. The H I cloud, which emerges from IC 2554 as a large arc-shaped plume, has a size of similar to30 kpc, larger than that of IC 2554. The total H I mass of the IC 2554 system is similar to2 x 10(9) M., one-third of which resides in the H I cloud. It is possible that tidal interaction between IC 2554 and NGC 3136B caused this spectacular H I cloud, but the possibility of IC 2554 being a merger remnant is also discussed. We also detected H I gas in the nearby galaxies ESO 092-G009 and RKK 1959 and an associated H I cloud, ATCA J1006-6710. Together they have an H I mass of similar to4.6 x 10(8) M-.. Another new H I source, ATCA J1007-6659, with an H I mass of only similar to2.2 x 10(7) M. was detected roughly between IC 2554 and ESO 092-G009 and corresponds to a face-on low surface brightness dwarf galaxy. Star formation is evident only in the galaxy IC 2554 with a rate of similar to4 M. yr(-1).

Galaxy threshing and ultra-compact dwarfs in the Fornax cluster

We have discovered a new type of galaxy in the Fornax Cluster: 'ultra-compact' dwarfs (UCDs). The UCDs are unresolved in ground-based imaging and have spectra typical of old stellar systems. Although the UCDs resemble overgrown globular clusters, based on VLT UVES echelle spectroscopy, they appear to be dynamically distinct systems with higher internal velocity dispersions and M/L ratios for a given luminosity than Milky Way or M31 globulars. Our preferred explanation for their origin is that they are the remnant nuclei of dwarf elliptical galaxies which have been tidally stripped, or 'threshed' by repeated encounters with the central cluster galaxy, NGC1399. If correct, then tidal stripping of nucleated dwarfs to form UCDs may, over a Hubble time, be an important source of the plentiful globular cluster population in the halo of NGC1399, and, by implication, other cD galaxies. In this picture, the dwarf elliptical halo contents, up to 99% of the original dwarf luminosity, contribute a significant fraction of the populations of intergalactic stars, globulars, and gas in galaxy clusters.

At the vigintennial of the Butcher-Oemler effect

In their study of the evolution of galaxies within clusters, Butcher and Oemler discovered evidence for a strong evolution in star-formation rate with redshift. Later studies confirmed this conclusion and uncovered several aspects of the effect: photometric, spectroscopic, and morphological. This article reviews a broad sample of these works and discusses selection effects, biases, and driving mechanisms that might be responsible for the changes in star-formation rate.

Australia Telescope Compact Array H I observations of the NGC 6845 galaxy group

We present the results of Australia Telescope Compact Array (ATCA) H i line and 20-cm radio continuum observations of the galaxy quartet NGC 6845. The H i emission extends over all four galaxies but can only be associated clearly with the two spiral galaxies, NGC 6845A and B, which show signs of strong tidal interaction. We derive a total H i mass of at least 1.8 x 10(10) M-., most of which is associated with NGC 6845A, the largest galaxy of the group. We investigate the tidal interaction between NGC 6845A and B by studying the kinematics of distinct H i components and their relation to the known H ii regions. No H i emission is detected from the two lenticular galaxies, NGC 6845C and D. A previously uncatalogued dwarf galaxy, ATCA J2001-4659, was detected 4.4 arcmin NE from NGC 6845B and has an H i mass of similar to5 x 10(8) M-.. No H i bridge is visible between the group and its newly detected companion. Extended 20-cm radio continuum emission is detected in NGC 6845A and B as well as in the tidal bridge between the two galaxies. We derive star formation rates of 15-40 M-. yr(-1).