Extensive air showers, originating from ultra-high energy cosmic rays, have been successfully measured through the use of arrays of water-Cherenkov detectors (WCDs). Sophisticated analyses exploiting WCD data have made it possible to demonstrate that shower simulations, based on different hadronic-interaction models, cannot reproduce the observed number of muons at the ground. The accurate knowledge of the WCD response to muons is paramount in establishing the exact level of this discrepancy. In this work, we report on a study of the response of a WCD of the Pierre Auger Observatory to atmospheric muons performed with a hodoscope made of resistive plate chambers (RPCs), enabling us to select and reconstruct nearly 600 thousand single muon trajectories with zenith angles ranging from 0$^\circ$ to 55$^\circ$. Comparison of distributions of key observables between the hodoscope data and the predictions of dedicated simulations allows us to demonstrate the accuracy of the latter at a level of 2%. As the WCD calibration is based on its response to atmospheric muons, the hodoscope data are also exploited to show the long-term stability of the procedure.

Starburst galaxies generate large-scale winds powered by the activity in the star-forming regions located in the galactic disks. Fragmentation of the disk produced by the outbreak of the wind results in the formation of clouds. Bowshocks caused by the supersonic outflow appear around such clouds. In this paper we discuss the acceleration of relativistic particles and the production of non-thermal radiation in such scenario. Cosmic rays accelerated at the bowshocks do not reach the highest energies, although the high-energy luminosity generated is significant. We show that up to ~10% of the gamma-ray emission in starbursts might come from these sources outside the galactic disks. Discrete X-ray sources with a power-law component are also expected.

Context. Detecting and characterizing substellar companions for which the luminosity, mass, and age can be determined independently is of utter importance to test and calibrate the evolutionary models due to uncertainties in their formation mechanisms. HD 19467 is a bright and nearby star hosting a cool brown dwarf companion detected with RV and imaging, making it a valuable object for such studies. Aims. We aim to further characterize the orbital, spectral, and physical properties of the HD 19467 system. Methods. We present new high-contrast imaging data with the SPHERE and NaCo instruments. We also analyze archival data from HARPS, NaCo, HIRES, UVES, and ASAS. We also use proper motion data of the star from Hipparcos and Gaia. Results. We refine the properties of the host star and derive an age of 8.0$^{+2.0}_{-1.0}$ Gyr based on isochrones, gyrochronology, and chemical and kinematic arguments. This estimate is slightly younger than previous estimates of ~9-11 Gyr. No orbital curvature is seen in the current imaging, RV, and astrometric data. From a joint fit of the data, we refine the orbital parameters for HD 19467B: period 398$^{+95}_{-93}$ yr, inclination 129.8$^{+8.1}_{-5.1}$ deg, eccentricity 0.56$\pm$0.09, longitude of the ascending node 134.8$\pm$4.5 deg, and argument of the periastron 64.2$^{+5.5}_{-6.3}$ deg. We assess a dynamical mass of 74$^{+12}_{-9}$ MJ. The fit with atmospheric models of the spectrophotometric data of HD 19467B indicates an atmosphere without clouds or with very thin clouds, an effective temperature of 1042$^{+77}_{-71}$ K, and a large surface gravity of 5.34$^{+0.08}_{-0.09}$ dex. The comparison to model predictions of the bolometric luminosity and dynamical mass of HD 19467B, assuming our system age estimate, indicates a better agreement with the Burrows et al. models; whereas the other evolutionary models used tend to underestimate its cooling rate.

We search for signals of magnetically-induced effects in the arrival directions of ultra-high-energy cosmic rays detected at the Pierre Auger Observatory. We apply two different methods. One is a search for sets of events that show a correlation between their arrival direction and the inverse of their energy, which would be expected if they come from the same point-like source, they have the same electric charge and their deflection is relatively small and coherent. We refer to these sets of events as "multiplets". The second method, called "thrust", is a principal axis analysis aimed to detect the elongated patterns in a region of interest. We study the sensitivity of both methods using a benchmark simulation and we apply them to data in two different searches. The first search is done assuming as source candidates a list of nearby active galactic nuclei and starburst galaxies. The second is an all-sky blind search. We report the results and we find no statistically significant features. We discuss the compatibility of these results with the indications on the mass composition inferred from data of the Pierre Auger Observatory.

Active galactic nuclei are supermassive black holes surrounded by an accretion disk, two populations of clouds, bipolar jets, and a dusty torus. The clouds move in Keplerian orbits at high velocities. In particular, the broad-line region (BLR) clouds have velocities ranging from $1000$ to $10000$ km s$^{-1}$. Given the extreme proximity of these clouds to the supermassive black hole, frequent collisions with the accretion disk should occur. The impact of BLR clouds onto the accretion disk can produce strong shock waves where particles might be accelerated. The goal of this work is to investigate the production of relativistic particles, and the associated non-thermal radiation in these events. In particular, we apply the model we develop to the Seyfert galaxy NGC 1068. We analyze the efficiency of diffusive shock acceleration in the shock of colliding clouds of the BLR with the accretion disk. We calculate the spectral energy distribution of photons generated by the relativistic particles and estimate the number of simultaneous impacts needed to explain the gamma radiation observed by the Fermi satellite in Seyfert galaxies. We find that is possible to understand the measured gamma emission in terms of the interaction of clouds with the disk if the hard X-ray emission of the source is at least obscured between $20\%$ and $40\%$. The total number of clouds contained in the BLR region might be between $3\times10^{8}$ and $6\times10^{8}$, which are values in good agreement with the observational evidence. The maximum energy achieved by the protons ($\sim$ PeV) in this context allows the production of neutrinos in the observing range of IceCube.

We present measurements of the large-scale cosmic-ray anisotropies in right ascension, using data collected by the surface detector array of the Pierre Auger Observatory over more than 14 years. We determine the equatorial dipole component, $\vec{d}_\perp$, through a Fourier analysis in right ascension that includes weights for each event so as to account for the main detector-induced systematic effects. For the energies at which the trigger efficiency of the array is small, the ``East-West'' method is employed. Besides using the data from the array with detectors separated by 1500 m, we also include data from the smaller but denser sub-array of detectors with 750 m separation, which allows us to extend the analysis down to $\sim 0.03$ EeV. The most significant equatorial dipole amplitude obtained is that in the cumulative bin above 8~EeV, $d_\perp=6.0^{+1.0}_{-0.9}$%, which is inconsistent with isotropy at the 6$\sigma$ level. In the bins below 8 EeV, we obtain 99% CL upper-bounds on $d_\perp$ at the level of 1 to 3 percent. At energies below 1 EeV, even though the amplitudes are not significant, the phases determined in most of the bins are not far from the right ascension of the Galactic center, at $\alpha_{\rm GC}=-94^\circ$, suggesting a predominantly Galactic origin for anisotropies at these energies. The reconstructed dipole phases in the energy bins above 4 EeV point instead to right ascensions that are almost opposite to the Galactic center one, indicative of an extragalactic cosmic ray origin.

The occurrence rate of long-period giant planets around young stars is highly uncertain since it is not only governed by the protoplanetary disc structure and planet formation process, but also reflects dynamical re-structuring processes after planet formation as well as possible capture of planets not formed in-situ. Direct imaging is currently the only feasible method to detect such wide-orbit planets and constrain their occurrence rate. We carry out a large L'-band high-contrast direct imaging survey for giant planets around young stars with protoplanetary or debris discs using the NACO instrument at the ESO Very Large Telescope on Cerro Paranal in Chile. We use very deep angular differential imaging observations with typically >60 deg field rotation, and employ a vector vortex coronagraph where feasible to achieve the best possible point source sensitivity down to an inner working angle of about 100mas. This paper introduces our NACO Imaging Survey for Planets around Young stars ("NACO-ISPY"), its goals and strategy, the target list, and data reduction scheme, and presents preliminary results from the first 2.5 survey years. We achieve a mean 5 sigma L' contrast of 6.4mag at 150mas and a background limit of 16.5mag at >1.5". Our detection probability is >50\% for companions with 8\,M$_{\rm Jup}$\ at semi-major axes 80-200au. It thus compares well to the detection space of other state-of-the-art high-contrast imaging surveys. We have contributed to the characterisation of two new planets originally discovered by VLT/SPHERE, but we have not yet independently discovered new planets around any of our target stars. We report the discovery of close-in low-mass stellar companions around four young stars and show L'-band scattered light images of the discs around eleven stars, six of which have never been imaged at L'-band before.

The combined effects of supernova explosions and stellar winds produce a hot bubble in the central regions of starburst galaxies. As the bubble expands, it can outbreak into the galactic halo driving a superwind that transports hot gas and fields to the intergalactic space. We present estimates of cosmic ray generation and gamma-ray emission in both this large-scale wind and the bow shocks created around the embedded clouds.

The Argentine Institute of Radio astronomy (IAR) is equipped with two single-dish 30mts radio antennas capable of performing daily observations of pulsars and radio transients in the southern hemisphere at 1.4 GHz. We aim to introduce to the international community the upgrades performed and to show that IAR observatory has become suitable for investigations in numerous areas of pulsar radio astronomy, such as pulsar timing arrays, targeted searches of continuous gravitational waves sources, monitoring of magnetars and glitching pulsars, and studies of short time scale interstellar scintillation. We refurbished the two antennas at IAR to achieve high-quality timing observations. We gathered more than $1\,000$ hours of observations with both antennas to study the timing precision and sensitivity they can achieve. We introduce the new developments for both radio telescopes at IAR. We present observations of the millisecond pulsar J0437$-$4715 with timing precision better than 1~$\mu$s. We also present a follow-up of the reactivation of the magnetar XTE J1810--197 and the measurement and monitoring of the latest (Feb. 1st. 2019) glitch of the Vela pulsar (J0835--4510). We show that IAR is capable of performing pulsar monitoring in the 1.4 GHz radio band for long periods of time with a daily cadence. This opens the possibility of pursuing several goals in pulsar science, including coordinated multi-wavelength observations with other observatories. In particular, observations of the millisecond pulsar J0437$-$4715 will increase the gravitational wave sensitivity of the NANOGrav array in their current blind spot. We also show IAR's great potential for studying targets of opportunity and transient phenomena such as magnetars, glitches, and fast-radio-burst sources.

With the Surface Detector array (SD) of the Pierre Auger Observatory we can detect neutrinos with energy between $10^{17}\,$eV and $10^{20}\,$eV from point-like sources across the sky, from close to the Southern Celestial Pole up to $60^\circ$ in declination, with peak sensitivities at declinations around $\sim -53^\circ$ and $\sim+55^\circ$, and an unmatched sensitivity for arrival directions in the Northern hemisphere. A search has been performed for highly-inclined air showers induced by neutrinos of all flavours with no candidate events found in data taken between 1 Jan 2004 and 31 Aug 2018. Upper limits on the neutrino flux from point-like steady sources have been derived as a function of source declination. An unrivaled sensitivity is achieved in searches for transient sources with emission lasting over an hour or less, if they occur within the field of view corresponding to the zenith angle range between $60^\circ$ and $~95^\circ$ where the SD of the Pierre Auger Observatory is most sensitive to neutrinos.

The determination of the primary energy of extensive air showers using the fluorescence detection technique requires an estimation of the energy carried away by particles that do not deposit all their energy in the atmosphere. This estimation is typically made using Monte Carlo simulations and thus depends on the assumed primary particle mass and on model predictions for neutrino and muon production. In this work we present a new method to obtain the invisible energy from events detected by the Pierre Auger Observatory. The method uses measurements of the muon number at ground level, and it allows us to reduce significantly the systematic uncertainties related to the mass composition and the high energy hadronic interaction models, and consequently to improve the estimation of the energy scale of the Observatory.

Neutrinos with energies above $10^{17}$ eV are detectable with the Surface Detector Array of the Pierre Auger Observatory. The identification is efficiently performed for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for Earth-skimming $\tau$ neutrinos with nearly tangential trajectories relative to the earth. No neutrino candidates were found in $\sim\,14.7$ years of data taken up to 31 August 2018. This leads to restrictive upper bounds on their flux. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high-energy neutrinos with an $E_\nu^{-2}$ spectrum in the energy range $1.0 \times 10^{17}~{\rm eV} - 2.5 \times 10^{19}~{\rm eV}$ is $E^2 {\rm d}N_\nu/{\rm d}E_\nu < 4.4 \times 10^{-9}~{\rm GeV~cm^{-2}~s^{-1}~sr^{-1}}$, placing strong constraints on several models of neutrino production at EeV energies and on the properties of the sources of ultra-high-energy cosmic rays.

H$\alpha$ is a powerful tracer of accretion and chromospheric activity, which has been detected in the case of young brown dwarfs and even recently in planetary mass companions (e.g. PDS70 b and c). H$\alpha$ detections and characterisation of brown dwarf and planet companions can further our knowledge of their formation and evolution, and expanding such a sample is therefore our primary goal. We used the Zurich IMaging POLarimeter (ZIMPOL) of the SPHERE instrument at the Very Large Telescope (VLT) to observe the known $38-72$M$_{\mathrm{J}}$ companion orbiting PZ Tel, obtaining simultaneous angular differential imaging observations in both continuum and narrow H$\alpha$ band. We detect H$\alpha$ emission from the companion, making this only the second H$\alpha$ detection of a companion using the SPHERE instrument. We used our newly added astrometric measurements to update the orbital analysis of PZ Tel B, and we used our photometric measurements to evaluate the H$\alpha$ line flux. Given the estimated bolometric luminosity, we obtained an H$\alpha$ activity (log$\mathrm{(L_{H\alpha}/L_{bol}})$) between $-4.16$ and $-4.31$. The H$\alpha$ activity of PZ Tel B is consistent with known average activity levels for M dwarf of the same spectral type. Given the absence of a known gaseous disk and the relatively old age of the system (24 Myr), we conclude that the H$\alpha$ emission around PZ Tel B is likely due to chromospheric activity.

Contributions of the Pierre Auger Collaboration to the 36th International Cosmic Ray Conference (ICRC 2019), 24 July - 1 August 2019, Madison, Wisconsin, USA.

The Southern Photometric Local Universe Survey (S-PLUS) is imaging ~9300 deg^2 of the celestial sphere in twelve optical bands using a dedicated 0.8 m robotic telescope, the T80-South, at the Cerro Tololo Inter-American Observatory, Chile. The telescope is equipped with a 9.2k by 9.2k e2v detector with 10 um pixels, resulting in a field-of-view of 2 deg^2 with a plate scale of 0.55"/pixel. The survey consists of four main subfields, which include two non-contiguous fields at high Galactic latitudes (8000 deg^2 at |b| > 30 deg) and two areas of the Galactic plane and bulge (for an additional 1300 deg^2). S-PLUS uses the Javalambre 12-band magnitude system, which includes the 5 u, g, r, i, z broad-band filters and 7 narrow-band filters centered on prominent stellar spectral features: the Balmer jump/[OII], Ca H+K, H-delta, G-band, Mg b triplet, H-alpha, and the Ca triplet. S-PLUS delivers accurate photometric redshifts (delta_z/(1+z) = 0.02 or better) for galaxies with r < 20 AB mag and redshift < 0.5, thus producing a 3D map of the local Universe over a volume of more than 1 (Gpc/h)^3. The final S-PLUS catalogue will also enable the study of star formation and stellar populations in and around the Milky Way and nearby galaxies, as well as searches for quasars, variable sources, and low-metallicity stars. In this paper we introduce the main characteristics of the survey, illustrated with science verification data highlighting the unique capabilities of S-PLUS. We also present the first public data release of ~336 deg^2 of the Stripe-82 area, which is available at http://datalab.noao.edu/splus.

Relative cross sections for $m$-fold photoionization ($m=1,\ldots,5$) of Fe$^{3+}$ by single photon absorption were measured employing the photon-ion merged-beams setup PIPE at the PETRA III synchrotron light source operated at DESY in Hamburg, Germany. The photon energies used spanned the range of $680-950\,\mathrm{eV}$, covering both the photoexcitation resonances from the $2p$ and $2s$ shells as well as the direct ionization from both shells. Multiconfiguration Dirac-Hartree-Fock (MCDHF) calculations were performed to simulate the total photoexcitation spectra. Good agreement was found with the experimental results. These computations helped to assign several strong resonance features to specific transitions. We also carried out Hartree-Fock calculations with relativistic extensions taking into account both photoexcitation and photoionization. Furthermore, we performed extensive MCDHF calculations of the Auger cascades that result when an electron is removed from the $2p$ and $2s$ shells of Fe$^{3+}$. Our theoretically predicted charge-state fractions are in good agreement with the experimental results, representing a substantial improvement over previous theoretical calculations. The main reason for the disagreement with the previous calculations is their lack of inclusion of slow Auger decays of several configurations that can only proceed when accompanied by de-excitation of two electrons. In such cases, this additional shake-down transition of a (sub-)valence electron is required to gain the necessary energy for the release of the Auger electron.

Context. The interaction between low-mass companions and the debris discs they reside in is still not fully understood. A debris disc can evolve due to self-stirring, a process in which planetesimals can excite their neighbours to the point of destructive collisions. In addition, the presence of a companion could further stir the disc (companion-stirring). Additional information is necessary to understand this fundamental step in the formation and evolution of a planetary system, and at the moment of writing only a handful of systems are known where a companion and a debris disc have both been detected and studied at the same time. Aims. Our primary goal is to augment the sample of these systems and to understand the relative importance between self-stirring and companion-stirring. Methods. In the course of the VLT/NaCo Imaging Survey for Planets around Young stars (ISPY), we observed HD 193571, an A0 debris disc hosting star at a distance of 68 pc with an age between 60 and 170 Myr. We obtained two sets of observations in L' band and a third epoch in H band using the GPI instrument at Gemini-South. Results. A companion was detected in all three epochs at a projected separation of 11 au (0.17 arcsec), and co-motion was confirmed through proper motion analysis. Given the inferred disc size of 120 au, the companion appears to reside within the gap between the host star and the disc. Comparison between the L' and H band magnitude and evolutionary tracks suggests a mass of 0.31 - 0.39 solar masses . Conclusions. We discovered a previously unknown M-dwarf companion around HD 193571, making it the third low-mass stellar object discovered within a debris disc. A comparison to self- and companion-stirring models suggests that the companion is likely responsible for the stirring of the disc.

The profile of the longitudinal development of showers produced by ultra-high energy cosmic rays carries information related to the interaction properties of the primary particles with atmospheric nuclei. In this work, we present the first measurement of the average shower profile in traversed atmospheric depth at the Pierre Auger Observatory. The shapes of profiles are well reproduced by the Gaisser-Hillas parametrization within the range studied, for E > 10^{17.8} eV. A detailed analysis of the systematic uncertainties is performed using 10 years of data and a full detector simulation. The average shape is quantified using two variables related to the width and asymmetry of the profile, and the results are compared with predictions of hadronic interaction models for different primary particles.

An overview of the multi-messenger capabilities of the Pierre Auger Observatory is presented. The techniques and performance of searching for Ultra-High Energy neutrinos, photons and neutrons are described. Some of the most relevant results are reviewed, such as stringent upper bounds that were placed to a flux of diffuse cosmogenic neutrinos and photons, bounds placed on neutrinos emitted from compact binary mergers that were detected by LIGO and Virgo during their first and second observing runs, as well as searches for high energy photons and neutrons from the Galactic center that constrain the properties of the putative Galactic PeVatron, observed by the H.E.S.S.\ collaboration. The observation of directional correlations between ultra-high energy cosmic rays and either high energy astrophysical neutrinos or specific source populations, weighted by their electromagnetic radiation, are also discussed. They constitute additional multi-messenger approaches aimed at identifying the sources of high energy cosmic rays.

Within the NaCo-ISPY exoplanet imaging program, we aim at detecting and characterizing the population of low-mass companions at wide separations ($\gtrsim$10AU), focusing in particular on young stars either hosting a known protoplanetary disk or a debris disk. R CrA is one of the youngest (1-3 Myr) and most promising objects in our sample because of two previous studies that suggested the presence of a close companion. Our aim is to directly image and characterize the companion for the first time. We observed R CrA twice with the NaCo instrument at VLT in the $L'$ filter with a one year time baseline in between. The high-contrast imaging data were reduced and analyzed, and in both datasets the companion candidate was detected. The companion is detected at a separation of $196.8\pm4.5$/$196.6\pm5.9$ mas ($18.7\pm1.3$/$18.7\pm1.4$ AU) and position angle of $134.7\pm0.5^\circ/133.7\pm0.7^\circ$ in the first/second epoch observation. We measure a contrast of $7.29\pm0.18$/$6.70\pm0.15$ mag with respect to the primary. Stellar proper motion study rejects the hypothesis of the signal being a background object. The companion candidate orbits in the clockwise direction and, if on a face-on circular orbit, its period is $\sim43-47$ yr. This value disagrees with the estimated orbital motion and therefore a face-on circular orbit may be excluded. Depending on the assumed age, extinction and brightness of the primary, the stellar companion has a mass between $0.10\pm0.02\, M_\odot$ and $1.03^{+0.20}_{-0.18}\,M_\odot$ range, if no contribution from circumsecondary material is taken into account. The presence of the companion needs to be taken into account when analyzing the complex circumstellar environment of R CrA.

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.

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.

Galactic all-sky maps at very disparate frequencies, like in the radio and $\gamma$-ray regime, show similar morphological structures. This mutual information reflects the imprint of the various physical components of the interstellar medium. We want to use multifrequency all-sky observations to test resolution improvement and restoration of unobserved areas for maps in certain frequency ranges. For this we aim to reconstruct or predict from sets of other maps all-sky maps that, in their original form, lack a high resolution compared to other available all-sky surveys or are incomplete in their spatial coverage. Additionally, we want to investigate the commonalities and differences that the ISM components exhibit over the electromagnetic spectrum. We build an $n$-dimensional representation of the joint pixel-brightness distribution of $n$ maps using a Gaussian mixture model and see how predictive it is: How well can one map be reproduced based on subsets of other maps? Tests with mock data show that reconstructing the map of a certain frequency from other frequency regimes works astonishingly well, predicting reliably small-scale details well below the spatial resolution of the initially learned map. Applied to the observed multifrequency data sets of the Milky Way this technique is able to improve the resolution of, e.g., the low-resolution Fermi LAT maps as well as to recover the sky from artifact-contaminated data like the ROSAT 0.855 keV map. The predicted maps generally show less imaging artifacts compared to the original ones. A comparison of predicted and original maps highlights surprising structures, imaging artifacts (fortunately not reproduced in the prediction), and features genuine to the respective frequency range that are not present at other frequency bands. We discuss limitations of this machine learning approach and ideas how to overcome them.

Eta Car is one of the most intriguing luminous blue variables in the Galaxy. Observations and models at different wavelengths suggest a central binary with a 5.54 yr period residing in its core. 2D and 3D radiative transfer and hydrodynamic simulations predict a primary with a dense and slow stellar wind that interacts with the faster and lower density wind of the secondary. The wind-wind collision scenario suggests that the secondary's wind penetrates the primary's wind creating a low-density cavity in it, with dense walls where the two winds interact. We aim to trace the inner ~5-50 au structure of Eta Car's wind-wind interaction, as seen through BrG and, for the first time, through the He I 2s-2p line. We have used spectro-interferometric observations with GRAVITY at the VLTI. Our modeling of the continuum allows us to estimate its FWHM angular size close to 2 mas and an elongation ratio of 1.06 +/- 0.05 over a PA = 130 +/- 20 deg. Our CMFGEN modeling helped us to confirm that the role of the secondary should be taken into account to properly reproduce the observed BrG and He I lines. Chromatic images across BrG reveal a southeast arc-like feature, possibly associated to the hot post-shocked winds flowing along the cavity wall. The images of He I 2s-2p served to constrain the 20 mas structure of the line-emitting region. The observed morphology of He I suggests that the secondary is responsible for the ionized material that produces the line profile. Both the BrG and the He I 2s-2p maps are consistent with previous hydrodynamical models of the colliding wind scenario. Future dedicated simulations together with an extensive interferometric campaign are necessary to refine our constraints on the wind and stellar parameters of the binary, which finally will help us predict the evolutionary path of Eta Car.

PDS 70b is a recently discovered and directly imaged exoplanet within the wide ($\gtrsim$40 au) cavity around PDS 70 (Keppler et al. 2018, M\"uller et al. 2018). Ongoing accretion onto the central star suggests that accretion onto PDS 70b may also be ongoing. We present the first high contrast images at H$\alpha$ (656 nm) and nearby continuum (643 nm) of PDS 70 utilizing the MagAO system. The combination of these filters allows for the accretion rate of the young planet to be inferred, as hot infalling hydrogen gas will emit strongly at H$\alpha$ over the optical continuum. We detected a source in H$\alpha$ at the position of PDS 70b on two sequential nights in May 2018, for which we establish a false positive probability of $<$0.1%. We conclude that PDS 70b is a young, actively accreting planet. We utilize the H$\alpha$ line luminosity to derive a mass accretion rate of $\dot M= 10^{-8\pm1}$ M$_{Jup}/yr$, where the large uncertainty is primarily due to the unknown amount of optical extinction from the circumstellar and circumplanetary disks. PDS 70b represents the second case of an accreting planet interior to a disk gap, and is among the early examples of a planet observed during its formation.

Young circumstellar disks are of prime interest to understand the physical and chemical conditions under which planet formation takes place. Only very few detections of planet candidates within these disks exist, and most of them are currently suspected to be disk features. In this context, the transition disk around the young star PDS 70 is of particular interest, due to its large gap identified in previous observations, indicative of ongoing planet formation. We aim to search for the presence of planets and search for disk structures indicative for disk-planet interactions and other evolutionary processes. We analyse new and archival near-infrared (NIR) images of the transition disk PDS 70 obtained with the VLT/SPHERE, VLT/NaCo and Gemini/NICI instruments in polarimetric differential imaging (PDI) and angular differential imaging (ADI) modes. We detect a point source within the gap of the disk at about 195 mas (about 22 au) projected separation. The detection is confirmed at five different epochs, in three filter bands and using different instruments. The astrometry results in an object of bound nature, with high significance. The comparison of the measured magnitudes and colours to evolutionary tracks suggests that the detection is a companion of planetary mass. We confirm the detection of a large gap of about 54 au in size within the disk in our scattered light images, and detect a signal from an inner disk component. We find that its spatial extent is very likely smaller than about 17 au in radius. The images of the outer disk show evidence of a complex azimuthal brightness distribution which may in part be explained by Rayleigh scattering from very small grains. Future observations of this system at different wavelengths and continuing astrometry will allow us to test theoretical predictions regarding planet-disk interactions, planetary atmospheres and evolutionary models.

Aims: We aim to characterize the orbital and atmospheric properties of PDS 70 b, which was first identified on May 2015 in the course of the SHINE survey with SPHERE, the extreme adaptive-optics instrument at the VLT. Methods: We obtained new deep SPHERE/IRDIS imaging and SPHERE/IFS spectroscopic observations of PDS 70 b. The astrometric baseline now covers 6 years which allows us to perform an orbital analysis. For the first time, we present spectrophotometry of the young planet which covers almost the entire near-infrared range (0.96 to 3.8 micrometer). We use different atmospheric models covering a large parameter space in temperature, log(g), chemical composition, and cloud properties to characterize the properties of the atmosphere of PDS 70 b. Results: PDS 70 b is most likely orbiting the star on a circular and disk coplanar orbit at ~22 au inside the gap of the disk. We find a range of models that can describe the spectrophotometric data reasonably well in the temperature range between 1000-1600 K and log(g) no larger than 3.5 dex. The planet radius covers a relatively large range between 1.4 and 3.7 R_jupiter with the larger radii being higher than expected from planet evolution models for the age of the planet of 5.4 Myr. Conclusions: This study provides a comprehensive dataset on the orbital motion of PDS 70 b, indicating a circular orbit and a motion coplanar with the disk. The first detailed spectral energy distribution of PDS 70 b indicates a temperature typical for young giant planets. The detailed atmospheric analysis indicates that a circumplanetary disk may contribute to the total planet flux.

Context: Runaway stars produce bowshocks that are usually observed at infrared (IR) wavelengths. Non-thermal radio emission has been detected so far only from the bowshock of BD+43{\deg}3654, whereas the detection of non-thermal radiation from these bowshocks at high energies remains elusive. Aims: We aim at characterising in detail the radio, X-ray, and gamma-ray emission from stellar bowshocks accounting for the structure of the region of interaction between the stellar wind and its environment. Methods: We develop a broadband-radiative, multi-zone model for stellar bowshocks that takes into account the spatial structure of the emitting region and the observational constraints. The model predicts the evolution and the emission of the relativistic particles accelerated and streaming together with the shocked flow. Results: We present broadband non-thermal spectral energy distributions for different scenarios, synthetic radio-cm synchrotron maps that reproduce the morphology of BD+43{\deg}3654, and updated predictions in X-ray and gamma-ray energy ranges. We also compare the results of the multi-zone model applied in this work with those of a refined one-zone model. Conclusions: A multi-zone model provides better constraints than a one-zone model on the relevant parameters, namely the magnetic field intensity and the amount of energy deposited in non-thermal particles. However, one-zone models can be improved by carefully characterising the intensity of the IR dust photon field and the escape rate of the plasma from the shocked region. Finally, comparing observed radio maps with those obtained from a multi-zone model enables constraints to be obtained on the direction of stellar motion with respect to the observer.

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}.

Starbursts are galaxies undergoing massive episodes of star formation. The combined effect of stellar winds from hot stars and supernova explosions creates a high-temperature cavity in the nuclear region of these objects. The very hot gas expands adiabatically and escapes from the galaxy creating a superwind which sweeps matter from the galactic disk. The superwind region in the halo is filled with a multi-phase gas with hot, warm, cool, and relativistic components. The shocks associated with the superwind of starbursts and the turbulent gas region of the bubble inflated by them might accelerate cosmic rays up to high energies. In this work we calculate the cosmic ray production associated with the superwind using parameters that correspond to the nearby southern starburst galaxy NGC 253, which has been suggested as a potential accelerator of ultra-high energy cosmic rays. We evaluate the efficiency of both diffusive shock acceleration (DSA) and stochastic diffusive acceleration (SDA) in the superwind of NGC 253. We estimate the distribution of both hadrons and leptons and calculate the corresponding spectral energy distributions of photons. We find that the strong mass load of the superwind, recently determined through ALMA observations, strongly attenuates the efficiency of DSA in NGC 253, whereas SDA is constrained by the age of the starburst. We conclude that NGC 253 and similar starbursts can only accelerate iron nuclei beyond $\sim10^{18}$ eV under very special conditions. If the central region of the galaxy harbors a starved supermassive black hole of $\sim10^6$ $M_{\odot}$, as suggested by some recent observations, a contribution in the range $10^{18}-10^{19}$ eV can be present for accretion rates $\dot{m}\sim10^{-3}$ in Eddington units. Only if very strong magnetic field amplification occurs close to the superwind shock energies of the order of 100 EeV might be achieved.

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.

High-velocity clouds (HVCs) are interstellar clouds of atomic hydrogen that do not partake of the Galactic rotation and have velocities of a several hundred kilometers per second. A considerable number of these clouds are falling down towards the Galactic disk. HVCs form large and massive complexes, so their collisions with the disk must release a great amount of energy into the interstellar medium. The cloud-disk interaction produces two shocks, one propagates through the cloud and the other through the disk; the properties of these shocks depend mainly on the cloud velocity and the disk-cloud density ratio. In this work we study the conditions necessary for these shocks to accelerate particles by diffusive shock acceleration and the produced non-thermal radiation. We analyze particle acceleration in both the cloud and disk shocks. Solving a time-dependent 2-D transport equation for both relativistic electrons and protons we obtain particle distributions and non-thermal spectral energy distributions. In a shocked cloud significant synchrotron radio emission is produced along with soft gamma rays. In the case of acceleration in the shocked disk, the non-thermal radiation is stronger; the gamma rays, of leptonic origin, might be detectable with current instruments. A large number of protons are injected into the Galactic interstellar medium, and locally exceed the cosmic-ray background. We conclude that under adequate conditions the contribution from HVC-disk collisions to the galactic population of relativistic particles and the associated extended non-thermal radiation might be important.

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.

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.

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

Absolute cross sections for m-fold photoionization (m=1,...,6) of Fe+ by a single photon were measured employing the photon-ion merged-beams setup PIPE at the PETRA III synchrotron light source, operated by DESY in Hamburg, Germany. Photon energies were in the range 680-920 eV which covers the photoionization resonances associated with 2p and 2s excitation to higher atomic shells as well as the thresholds for 2p and 2s ionization. The corresponding resonance positions were measured with an uncertainty of +- 0.2 eV. The cross section for Fe+ photoabsorption is derived as the sum of the individually measured cross-sections for m-fold ionization. Calculations of the Fe+ absorption cross sections have been carried out using two different theoretical approaches, Hartree-Fock including relativistic extensions and fully relativistic Multi-Configuration Dirac Fock. Apart from overall energy shifts of up to about 3 eV, the theoretical cross sections are in good agreement with each other and with the experimental results. In addition, the complex deexcitation cascades after the creation of inner-shell holes in the Fe+ ion have been tracked on the atomic fine-structure level. The corresponding theoretical results for the product charge-state distributions are in much better agreement with the experimental data than previously published configuration-average results. The present experimental and theoretical results are valuable for opacity calculations and are expected to pave the way to a more accurate determination of the iron abundance in the interstellar medium.