Aims: This paper has two main objectives: (1) To determine the intrinsic properties of 16 faint and mostly unstudied open clusters in the poorly known sector of the Galaxy at 270$^\circ-$300$^\circ$, to probe the Milky Way structure in future investigations. (2) To address previously reported systematics in Gaia DR2 parallaxes by comparing the cluster distances derived from photometry with those derived from parallaxes. Methods: Deep UBVI photometry of 16 open clusters was carried out. Observations were reduced and analyzed in an automaticway using the ASteCA package to get individual distances, reddening, masses, ages and metallicities. Photometric distances were compared to those obtained from a Bayesian analysis of Gaia DR2 parallaxes. Results: Ten out of the 16 clusters are true or highly probable open clusters. Two of them are quite young and follow the trace of the Carina Arm and the already detected warp. The rest of the clusters are placed in the interarm zone between the Perseus and Carina Arms as expected for older objects. We found that the cluster van den Berg-Hagen 85 is 7.5$\times$10$^9$ yrs old becoming then one of the oldest open cluster detected in our Galaxy so far. The relationship of these ten clusters with the Galaxy structure in the solar neighborhood is discussed. The comparison of distances from photometry and parallaxes data, in turn, reveals a variable level of disagreement. Conclusions: Various zero point corrections for Gaia DR2 parallax data recently reported were considered for a comparison between photometric and parallax based distances. The results tend to improve with some of these corrections. Photometric distance analysis suggest an average correction of $\sim$+0.026 mas (to be added to the parallaxes). The correction may have a more intricate distance dependency, but addressing that level of detail will require a larger cluster sample.

We present an analysis of the binary central star of the planetary nebula NGC 2346 based on archival data from the International Ultraviolet Explorer (IUE), and new low- and high-resolution optical spectra (3700 - 7300{\AA}). By including in the spectral analysis the contribution of both stellar and nebular continuum, we reconciled long-time discrepant UV and optical diagnostics and derive $E(B-V)=0.18\pm0.01$. We re-classified the companion star as A5IV by analyzing the wings of the Balmer absorption lines in the high-resolution ($R=67\,000$) optical spectra. Using the distance to the nebula of 1400 pc from Gaia DR2, we constructed a photoionization model based on abundances and line intensities derived from the low-resolution optical spectra, and obtained a temperature of $T_{\rm eff}=130\,000$K and a luminosity $L=170$L$_\odot$ for the ionizing star, consistent with the UV continuum. This analysis allows us to better characterize the binary system's evolution. We conclude that the progenitor star of NGC 2346 has experienced a common envelope phase, in which the companion star has accreted mass and evolved off the main-sequence.

We present a polarimetric investigation of the protoplanetary nebula Frosty Leo performed with the Submillimeter Array. We were able to detect, in the low continuum level (peak at 14.4 mJy beam$^{-1}$), a marginal polarization at $\sim2.6\sigma$. The molecular line investigation based on the CO $J=3\rightarrow2$ emission shows a peak emission of 68.1 Jy beam$^{-1}$ km s$^{-1}$ and the polarization detection in this CO line is also marginal, with a peak at $\sim3.8\sigma$. In both cases, it was therefore not possible to use the electric vector maps (B-field) to accurately trace the magnetic field (B-field) within the PPN. The spatio-kinematic modelling realised with the different velocity channel maps indicates three main structures: a distorted torus accompanied by a bipolar outflow or jet aligned with its axis and a flattened spherical "cap". The comparison of the CO polarization segments with our model suggests that the polarized emission probably arises in the first two components.

We present a detailed study of the large-scale anisotropies of cosmic rays with energies above 4 EeV measured using the Pierre Auger Observatory. For the energy bins [4,8] EeV and $E\geq 8$ EeV, the most significant signal is a dipolar modulation in right ascension at energies above 8 EeV, as previously reported. In this paper we further scrutinize the highest-energy bin by splitting it into three energy ranges. We find that the amplitude of the dipole increases with energy above 4 EeV. The growth can be fitted with a power law with index $\beta=0.79\pm 0.19$. The directions of the dipoles are consistent with an extragalactic origin of these anisotropies at all the energies considered. Additionally we have estimated the quadrupolar components of the anisotropy: they are not statistically significant. We discuss the results in the context of the predictions from different models for the distribution of ultrahigh-energy sources and cosmic magnetic fields.

The parallax expansion and kinematics of a nova shell can be used to assess its age and distance, and to investigate the interaction of the ejecta with the circumstellar medium. These are key to understand the expansion and dispersal of the nova ejecta in the Galaxy. Multi-epoch images and high-dispersion spectroscopic observations of the recently discovered classical nova shell IPHASXJ210204.7+471015 around a nova-like system have been used to derive a present day expansion rate of 0.100 arcsec/yr and an expansion velocity of 285 km/s. These data are combined to obtain a distance of 600 pc to the nova. The secular expansion of the nova shell place the event sometime between 1850 and 1890, yet it seems to have been missed at that time. Despite its young age, 130-170 yrs, we found indications that the ejecta has already experienced a noticeable deceleration, indicating the interaction of this young nova shell with the surrounding medium.

With the Auger Engineering Radio Array (AERA) of the Pierre Auger Observatory, we have observed the radio emission from 561 extensive air showers with zenith angles between 60$^\circ$ and 84$^\circ$. In contrast to air showers with more vertical incidence, these inclined air showers illuminate large ground areas of several km$^2$ with radio signals detectable in the 30 to 80\,MHz band. A comparison of the measured radio-signal amplitudes with Monte Carlo simulations of a subset of 50 events for which we reconstruct the energy using the Auger surface detector shows agreement within the uncertainties of the current analysis. As expected for forward-beamed radio emission undergoing no significant absorption or scattering in the atmosphere, the area illuminated by radio signals grows with the zenith angle of the air shower. Inclined air showers with EeV energies are thus measurable with sparse radio-antenna arrays with grid sizes of a km or more. This is particularly attractive as radio detection provides direct access to the energy in the electromagnetic cascade of an air shower, which in case of inclined air showers is not accessible by arrays of particle detectors on the ground.

An in-situ calibration of a logarithmic periodic dipole antenna with a frequency coverage of 30 MHz to 80 MHz is performed. Such antennas are part of a radio station system used for detection of cosmic ray induced air showers at the Engineering Radio Array of the Pierre Auger Observatory, the so-called Auger Engineering Radio Array (AERA). The directional and frequency characteristics of the broadband antenna are investigated using a remotely piloted aircraft (RPA) carrying a small transmitting antenna. The antenna sensitivity is described by the vector effective length relating the measured voltage with the electric-field components perpendicular to the incoming signal direction. The horizontal and meridional components are determined with an overall uncertainty of 7.4^{+0.9}_{-0.3} % and 10.3^{+2.8}_{-1.7} % respectively. The measurement is used to correct a simulated response of the frequency and directional response of the antenna. In addition, the influence of the ground conductivity and permittivity on the antenna response is simulated. Both have a negligible influence given the ground conditions measured at the detector site. The overall uncertainties of the vector effective length components result in an uncertainty of 8.8^{+2.1}_{-1.3} % in the square root of the energy fluence for incoming signal directions with zenith angles smaller than 60{\deg}.

We present a combined fit of a simple astrophysical model of UHECR sources to both the energy spectrum and mass composition data measured by the Pierre Auger Observatory. The fit has been performed for energies above $5 \cdot 10^{18}$ eV, i.e.~the region of the all-particle spectrum above the so-called "ankle" feature. The astrophysical model we adopted consists of identical sources uniformly distributed in a comoving volume, where nuclei are accelerated through a rigidity-dependent mechanism. The fit results suggest sources characterized by relatively low maximum injection energies, hard spectra and heavy chemical composition. We also show that uncertainties about physical quantities relevant to UHECR propagation and shower development have a non-negligible impact on the fit results.

A new analysis of the dataset from the Pierre Auger Observatory provides evidence for anisotropy in the arrival directions of ultra-high-energy cosmic rays on an intermediate angular scale, which is indicative of excess arrivals from strong, nearby sources. The data consist of 5514 events above 20 EeV with zenith angles up to 80 deg recorded before 2017 April 30. Sky models have been created for two distinct populations of extragalactic gamma-ray emitters: active galactic nuclei from the second catalog of hard Fermi-LAT sources (2FHL) and starburst galaxies from a sample that was examined with Fermi-LAT. Flux-limited samples, which include all types of galaxies from the Swift-BAT and 2MASS surveys, have been investigated for comparison. The sky model of cosmic-ray density constructed using each catalog has two free parameters, the fraction of events correlating with astrophysical objects and an angular scale characterizing the clustering of cosmic rays around extragalactic sources. A maximum-likelihood ratio test is used to evaluate the best values of these parameters and to quantify the strength of each model by contrast with isotropy. It is found that the starburst model fits the data better than the hypothesis of isotropy with a statistical significance of 4.0 sigma, the highest value of the test statistic being for energies above 39 EeV. The three alternative models are favored against isotropy with 2.7-3.2 sigma significance. The origin of the indicated deviation from isotropy is examined and prospects for more sensitive future studies are discussed.

With its bright and wide equatorial waist seen almost edge-on (the butterfly body) and the faint and broad bipolar extensions (the butterfly wings), NGC650-1 is the archetypical example of bipolar planetary nebula (PN) with butterfly morphology. We present here deep high-resolution broad- and narrow-band optical images that expose the rich and intricate fine structure of this bipolar PN, with small-scale bubble-like features and collimated outflows. A SHAPE spatio-kinematical model indicates that NGC650-1 has a broad central torus with an inclination angle of 75 degrees with respect to the line of sight, whereas that of the bipolar lobes, which are clearly seen in the position-velocity maps, is 85 degrees. Large field of view deep images show, for first time, an arc-like diffuse envelope in low- and high-excitation emission lines located up to 180 arcsec towards the East-Southeast of the central star, well outside the main nebula. This morphological component is confirmed by Spitzer MIPS and WISE infrared imaging, as well as by long-slit low- and high-dispersion optical spectroscopic observations. HST images of NGC650-1 obtained at two different epochs (14 yrs apart), reveal the proper motion of the central star along this direction. We propose that this motion of the star through the interstellar medium compresses the remnant material of a slow Asymptotic Giant Branch wind, producing this bow-shock-like feature.

We present a new method for probing the hadronic interaction models at ultra-high energy and extracting details about mass composition. This is done using the time profiles of the signals recorded with the water-Cherenkov detectors of the Pierre Auger Observatory. The profiles arise from a mix of the muon and electromagnetic components of air-showers. Using the risetimes of the recorded signals we define a new parameter, which we use to compare our observations with predictions from simulations. We find, firstly, inconsistencies between our data and predictions over a greater energy range and with substantially more events than in previous studies. Secondly, by calibrating the new parameter with fluorescence measurements from observations made at the Auger Observatory, we can infer the depth of shower maximum for a sample of over 81,000 events extending from 0.3 EeV to over 100 EeV. Above 30 EeV, the sample is nearly fourteen times larger than currently available from fluorescence measurements and extending the covered energy range by half a decade. The energy dependence of the average depth of shower maximum is compared to simulations and interpreted in terms of the mean of the logarithmic mass. We find good agreement with previous work and extend the measurement of the mean depth of shower maximum to greater energies than before, reducing significantly the statistical uncertainty associated with the inferences about mass composition.

AMIGA (Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory designed to extend its energy range of detection and to directly measure the muon content of the cosmic ray primary particle showers. The array will be formed by an infill of surface water-Cherenkov detectors associated with buried scintillation counters employed for muon counting. Each counter is composed of three scintillation modules, with a 10 m$^2$ detection area per module. In this paper, a new generation of detectors, replacing the current multi-pixel photomultiplier tube (PMT) with silicon photo sensors (aka. SiPMs), is proposed. The selection of the new device and its front-end electronics is explained. A method to calibrate the counting system that ensures the performance of the detector is detailed. This method has the advantage of being able to be carried out in a remote place such as the one where the detectors are deployed. High efficiency results, i.e. 98 % efficiency for the highest tested overvoltage, combined with a low probability of accidental counting ($\sim$2 %), show a promising performance for this new system.

Contributions of the Pierre Auger Collaboration to the 35th International Cosmic Ray Conference (ICRC 2017), 12-20 July 2017, Bexco, Busan, Korea.

We present a novel method to measure precisely the relative spectral response of the fluorescence telescopes of the Pierre Auger Observatory. We used a portable light source based on a xenon flasher and a monochromator to measure the relative spectral efficiencies of eight telescopes in steps of 5 nm from 280 nm to 440 nm. Each point in a scan had approximately 2 nm FWHM out of the monochromator. Different sets of telescopes in the observatory have different optical components, and the eight telescopes measured represent two each of the four combinations of components represented in the observatory. We made an end-to-end measurement of the response from different combinations of optical components, and the monochromator setup allowed for more precise and complete measurements than our previous multi-wavelength calibrations. We find an overall uncertainty in the calibration of the spectral response of most of the telescopes of 1.5% for all wavelengths; the six oldest telescopes have larger overall uncertainties of about 2.2%. We also report changes in physics measureables due to the change in calibration, which are generally small.

Cosmic rays are atomic nuclei arriving from outer space that reach the highest energies observed in nature. Clues to their origin come from studying the distribution of their arrival directions. Using $3 \times 10^4$ cosmic rays above $8 \times 10^{18}$ electron volts, recorded with the Pierre Auger Observatory from a total exposure of 76,800 square kilometers steradian year, we report an anisotropy in the arrival directions. The anisotropy, detected at more than the 5.2$\sigma$ level of significance, can be described by a dipole with an amplitude of $6.5_{-0.9}^{+1.3}$% towards right ascension $\alpha_{d} = 100 \pm 10$ degrees and declination $\delta_{d} = -24_{-13}^{+12}$ degrees. That direction indicates an extragalactic origin for these ultra-high energy particles.

We report a multi-resolution search for anisotropies in the arrival directions of cosmic rays detected at the Pierre Auger Observatory with local zenith angles up to $80^\circ$ and energies in excess of 4 EeV ($4 \times 10^{18}$ eV). This search is conducted by measuring the angular power spectrum and performing a needlet wavelet analysis in two independent energy ranges. Both analyses are complementary since the angular power spectrum achieves a better performance in identifying large-scale patterns while the needlet wavelet analysis, considering the parameters used in this work, presents a higher efficiency in detecting smaller-scale anisotropies, potentially providing directional information on any observed anisotropies. No deviation from isotropy is observed on any angular scale in the energy range between 4 and 8 EeV. Above 8 EeV, an indication for a dipole moment is captured; while no other deviation from isotropy is observed for moments beyond the dipole one. The corresponding $p$-values obtained after accounting for searches blindly performed at several angular scales, are $1.3 \times 10^{-5}$ in the case of the angular power spectrum, and $2.5 \times 10^{-3}$ in the case of the needlet analysis. While these results are consistent with previous reports making use of the same data set, they provide extensions of the previous works through the thorough scans of the angular scales.

A search for ultra-high energy photons with energies above 1 EeV is performed using nine years of data collected by the Pierre Auger Observatory in hybrid operation mode. An unprecedented separation power between photon and hadron primaries is achieved by combining measurements of the longitudinal air-shower development with the particle content at ground measured by the fluorescence and surface detectors, respectively. Only three photon candidates at energies 1 - 2 EeV are found, which is compatible with the expected hadron-induced background. Upper limits on the integral flux of ultra-high energy photons of 0.027, 0.009, 0.008, 0.008 and 0.007 km$^{-2}$ sr$^{-1}$ yr$^{-1}$ are derived at 95% C.L. for energy thresholds of 1, 2, 3, 5 and 10 EeV. These limits bound the fractions of photons in the all-particle integral flux below 0.1%, 0.15%, 0.33%, 0.85% and 2.7%. For the first time the photon fraction at EeV energies is constrained at the sub-percent level. The improved limits are below the flux of diffuse photons predicted by some astrophysical scenarios for cosmogenic photon production. The new results rule-out the early top-down models $-$ in which ultra-high energy cosmic rays are produced by, e.g., the decay of super-massive particles $-$ and challenge the most recent super-heavy dark matter models.

Simultaneous measurements of air showers with the fluorescence and surface detectors of the Pierre Auger Observatory allow a sensitive search for EeV photon point sources. Several Galactic and extragalactic candidate objects are grouped in classes to reduce the statistical penalty of many trials from that of a blind search and are analyzed for a significant excess above the background expectation. The presented search does not find any evidence for photon emission at candidate sources, and combined $p$-values for every class are reported. Particle and energy flux upper limits are given for selected candidate sources. These limits significantly constrain predictions of EeV proton emission models from non-transient Galactic and nearby extragalactic sources, as illustrated for the particular case of the Galactic center region.

We present the first part of an investigation on the planetary nebula (PN) PC 22 which focuses on the use of deep imaging and high resolution echelle spectroscopy to perform a detailed morpho-kinematical analysis. PC 22 is revealed to be a multipolar PN emitting predominantly in [O III] and displaying multiple non-symmetric outflows. Its central region is found to be also particularly inhomogeneous with a series of low ionization structures (knots) located on the path of the outflows. The morpho-kinematical model obtained with SHAPE indicates that i) the de-projected velocities of the outflows are rather large, > 100 km/s, while the central region has expansion velocities in the range ~25 to ~45 km/s following the "Wilson effect", ii) the majority of the measured structures share similar inclination, ~100 degrees, i.e. they are coplanar, and iii) all outflows and lobes are coeval (within the uncertainties). All these results make us to suggest that PC 22 is an evolved starfish PN. We propose that the mechanism responsible for the morphology of PC 22 consists of a wind-shell interaction, where the fast post-AGB wind flows through a filamentary AGB shell with some large voids.

Aims. To produce an homogeneous catalog of astrophysical parameters of 239 resolved star clusters located in the Small and Large Magellanic Clouds, observed in the Washington photometric system. Methods. The cluster sample was processed with the recently introduced Automated Stellar Cluster Analysis ($\texttt{ASteCA}$) package, which ensures both an automatized and a fully reproducible treatment, together with a statistically based analysis of their fundamental parameters and associated uncertainties. The fundamental parameters determined with this tool for each cluster, via a color-magnitude diagram (CMD) analysis, are: metallicity, age, reddening, distance modulus, and total mass. Results. We generated an homogeneous catalog of structural and fundamental parameters for the studied cluster sample, and performed a detailed internal error analysis along with a thorough comparison with values taken from twenty-six published articles. We studied the distribution of cluster fundamental parameters in both Clouds, and obtained their age-metallicity relationships. Conclusions. The $\texttt{ASteCA}$ package can be applied to an unsupervised determination of fundamental cluster parameters; a task of increasing relevance as more data becomes available through upcoming surveys.

On September 14, 2015 the Advanced LIGO detectors observed their first gravitational-wave (GW) transient GW150914. This was followed by a second GW event observed on December 26, 2015. Both events were inferred to have arisen from the merger of black holes in binary systems. Such a system may emit neutrinos if there are magnetic fields and disk debris remaining from the formation of the two black holes. With the surface detector array of the Pierre Auger Observatory we can search for neutrinos with energy above 100 PeV from point-like sources across the sky with equatorial declination from about -65 deg. to +60 deg., and in particular from a fraction of the 90% confidence-level (CL) inferred positions in the sky of GW150914 and GW151226. A targeted search for highly-inclined extensive air showers, produced either by interactions of downward-going neutrinos of all flavors in the atmosphere or by the decays of tau leptons originating from tau-neutrino interactions in the Earth's crust (Earth-skimming neutrinos), yielded no candidates in the Auger data collected within $\pm 500$ s around or 1 day after the coordinated universal time (UTC) of GW150914 and GW151226, as well as in the same search periods relative to the UTC time of the GW candidate event LVT151012. From the non-observation we constrain the amount of energy radiated in ultrahigh-energy neutrinos from such remarkable events.

We present narrow-, broad-band, and WISE archive images, and high- and intermediate-resolution long-slit spectra of K 4-37, a planetary nebula that has never been analyzed in detail. Although K 4-37 appears bipolar, the morphokinematical analysis discloses the existence of three distinct axes and additional particular directions in the object, indicating that K 4-37 is a multi-axis planetary nebula that has probably been shaped by several bipolar outflows at different directions. A 4-6 M$_{\odot}$ main-sequence progenitor is estimated from the derived high nebular He and N abundances, and very high N/O abundance ratio ($\sim$2.32). The general properties are compatible with K 4-37 being a highly evolved planetary nebula located at $\sim$14 kpc. The WISE image at 22 $\mu$m reveals K 4-37 to be surrounded by a large ($\sim$13$\times$8 pc$^2$) elliptical detached shell probably related to material ejected from the AGB progenitor. The observed elliptical morphology suggests deformation of an originally spherical AGB shell by the ISM magnetic field or by the influence of a companion. We compare K 4-37 and NGC 6309 and found remarkable similarities in their physical structure but noticeably different chemical abundances that indicate very different progenitor mass. This strongly suggests that, irrespective of the initial mass, their (presumably binary) central stars have shared a very similar mass ejection history.

We report a first measurement for ultra-high energy cosmic rays of the correlation between the depth of shower maximum and the signal in the water Cherenkov stations of air-showers registered simultaneously by the fluorescence and the surface detectors of the Pierre Auger Observatory. Such a correlation measurement is a unique feature of a hybrid air-shower observatory with sensitivity to both the electromagnetic and muonic components. It allows an accurate determination of the spread of primary masses in the cosmic-ray flux. Up till now, constraints on the spread of primary masses have been dominated by systematic uncertainties. The present correlation measurement is not affected by systematics in the measurement of the depth of shower maximum or the signal in the water Cherenkov stations. The analysis relies on general characteristics of air showers and is thus robust also with respect to uncertainties in hadronic event generators. The observed correlation in the energy range around the `ankle' at $\lg(E/{\rm eV})=18.5-19.0$ differs significantly from expectations for pure primary cosmic-ray compositions. A light composition made up of proton and helium only is equally inconsistent with observations. The data are explained well by a mixed composition including nuclei with mass $A > 4$. Scenarios such as the proton dip model, with almost pure compositions, are thus disfavoured as the sole explanation of the ultrahigh-energy cosmic-ray flux at Earth.

Ultrahigh energy cosmic ray air showers probe particle physics at energies beyond the reach of accelerators. Here we introduce a new method to test hadronic interaction models without relying on the absolute energy calibration, and apply it to events with primary energy 6-16 EeV (E_CM = 110-170 TeV), whose longitudinal development and lateral distribution were simultaneously measured by the Pierre Auger Observatory. The average hadronic shower is 1.33 +- 0.16 (1.61 +- 0.21) times larger than predicted using the leading LHC-tuned models EPOS-LHC (QGSJetII-04), with a corresponding excess of muons.

We present a search for ultra-relativistic magnetic monopoles with the Pierre Auger Observatory. Such particles, possibly a relic of phase transitions in the early universe, would deposit a large amount of energy along their path through the atmosphere, comparable to that of ultrahigh-energy cosmic rays (UHECRs). The air shower profile of a magnetic monopole can be effectively distinguished by the fluorescence detector from that of standard UHECRs. No candidate was found in the data collected between 2004 and 2012, with an expected background of less than 0.1 event from UHECRs. The corresponding 90% confidence level (C.L.) upper limits on the flux of ultra-relativistic magnetic monopoles range from $10^{-19}$ (cm$^{2}$ sr s)$^{-1}$ for a Lorentz factor $\gamma=10^9$ to $2.5 \times10^{-21}$ (cm$^{2}$ sr s)$^{-1}$ for $\gamma=10^{12}$. These results - the first obtained with a UHECR detector - improve previously published limits by up to an order of magnitude.

We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principle calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.