• ### Particle acceleration in the superwinds of starburst galaxies(1801.06483)

April 6, 2018 astro-ph.HE
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.
• ### The Effect of Toroidal Magnetic Fields on Solar Oscillation Frequencies(1801.07932)

Jan. 24, 2018 astro-ph.SR
Solar oscillation frequencies change with the level of magnetic activity. Localizing subsurface magnetic field concentrations in the Sun with helioseismology will help us to understand the solar dynamo. Because the magnetic fields are not considered in standard solar models, adding them to the basic equations of stellar structure changes the eigenfunctions and eigenfrequencies. We use quasi-degenerate perturbation theory to calculate the effect of toroidal magnetic fields on solar oscillation mean multiplet frequencies for six field configurations. In our calculations, we consider both the direct effect of the magnetic field, which describes the coupling of modes, and the indirect effect, which accounts for changes in stellar structure due to the magnetic field. We limit our calculations to self-coupling of modes. We find that the magnetic field affects the multiplet frequencies in a way that depends on the location and the geometry of the field inside the Sun. Comparing our theoretical results with observed shifts, we find that strong tachocline fields cannot be responsible for the observed frequency shifts of p modes over the solar cycle. We also find that part of the surface effect in helioseismic oscillation frequencies might be attributed to magnetic fields in the outer layers of the Sun. The theory presented here is also applicable to models of solar-like stars and their oscillation frequencies.
• ### The Direct Effect of Toroidal Magnetic Fields on Stellar Oscillations: An Analytical Expression for the General Matrix Element(1709.02454)

Sept. 7, 2017 astro-ph.SR
Where is the solar dynamo located and what is its modus operandi? These are still open questions in solar physics. Helio- and asteroseismology can help answer them by enabling us to study solar and stellar internal structures through global oscillations. The properties of solar and stellar acoustic modes are changing with the level of magnetic activity. However, until now, the inference on subsurface magnetic fields with seismic measures has been very limited. The aim of this paper is to develop a formalism to calculate the effect of large-scale toroidal magnetic fields on solar and stellar global oscillation eigenfunctions and eigenfrequencies. If the Lorentz force is added to the equilibrium equation of motion, stellar eigenmodes can couple. In quasi-degenerate perturbation theory, this coupling, also known as the direct effect, can be quantified by the general matrix element. We present the analytical expression of the matrix element for a superposition of subsurface zonal toroidal magnetic field configurations. The matrix element is important for forward calculations of perturbed solar and stellar eigenfunctions and frequency perturbations. The results presented here will help to ascertain solar and stellar large-scale subsurface magnetic fields, and their geometric configuration, strength, and their change over the course of activity cycles.
• ### Inversions for Deep Solar Meridional Flow Using Spherical Born Kernels(1707.08803)

July 27, 2017 astro-ph.SR
The solar meridional flow is a crucial ingredient in modern dynamo theory. Seismic estimates of this flow have, however, been contradictory in deeper layers below about $0.9\,R_\odot$. Results from time-distance helioseismology have so far been obtained using the ray approximation. Here, we perform inversions using the Born approximation. The initial result is similar to the result previously obtained by Jackiewicz et al. (2015) using ray kernels while using the same set of GONG data and the SOLA inversion technique. However, we show that the assumption of uncorrelated measurements used in earlier studies may lead to inversion errors being underestimated by a factor of about two to four. In a second step, refined inversions are performed using the full covariance matrix and a regularization for cross-talk. As the results are found to depend on the threshold used in the singular value decomposition, they were obtained for a medium threshold ($10^{-7}-10^{-5}$, about 50% of the values used) and a threshold lower by a factor of 10 (about 70% of the values used). The result obtained with the medium threshold is again similar to the original, with less latitudinal variation. However, using the lower threshold, the inverted flow in the southern hemisphere shows two or three cells stacked radially depending on the associated radial flows. Both the single-cell and the multi-cell profiles are consistent with the measured travel times. All our results confirm a shallow return flow at about $0.9\,R_\odot$.
• ### A new look at sunspot formation using theory and observations(1704.04062)

April 13, 2017 astro-ph.SR
Sunspots are of basic interest in the study of the Sun. Their relevance ranges from them being an activity indicator of magnetic fields to being the place where coronal mass ejections and flares erupt. They are therefore also an important ingredient of space weather. Their formation, however, is still an unresolved problem in solar physics. Observations utilize just 2D surface information near the spot, but it is debatable how to infer deep structures and properties from local helioseismology. For a long time, it was believed that flux tubes rising from the bottom of the convection zone are the origin of the bipolar sunspot structure seen on the solar surface. However, this theory has been challenged, in particular recently by new surface observation, helioseismic inversions, and numerical models of convective dynamos. In this article we discuss another theoretical approach to the formation of sunspots: the negative effective magnetic pressure instability. This is a large-scale instability, in which the total (kinetic plus magnetic) turbulent pressure can be suppressed in the presence of a weak large-scale magnetic field, leading to a converging downflow, which eventually concentrates the magnetic field within it. Numerical simulations of forced stratified turbulence have been able to produce strong super-equipartition flux concentrations, similar to sunspots at the solar surface. In this framework, sunspots would only form close to the surface due to the instability constraints on stratification and rotation. Additionally, we present some ideas from local helioseismology, where we plan to use the Hankel analysis to study the pre-emergence phase of a sunspot and to constrain its deep structure and formation mechanism.
• ### Validation of Spherical Born Approximation Sensitivity Functions for Measuring Deep Solar Meridional Flow(1703.03700)

March 10, 2017 astro-ph.SR
Accurate measurements of deep solar meridional flow are of vital interest for understanding the solar dynamo. In this paper, we validate a recently developed method for obtaining sensitivity functions (kernels) for travel-time measurements to solar interior flows using the Born approximation in spherical geometry, which is expected to be more accurate than the classical ray approximation. Furthermore, we develop a numerical approach to efficiently compute a large number of kernels based on the separability of the eigenfunctions into their horizontal and radial dependence. The validation is performed using a hydrodynamic simulation of linear wave propagation in the Sun, which includes a standard single-cell meridional flow profile. We show that, using the Born approximation, it is possible to accurately model observational quantities relevant for time-distance helioseismology such as the mean power spectrum, disc-averaged cross-covariance functions, and travel times in the presence of a flow field. In order to closely match the model to observations, we show that it is beneficial to use mode frequencies and damping rates which were extracted from the measured power spectrum. Furthermore, the contribution of the radial flow to the total travel time is found to reach 20% of the contribution of the horizontal flow at travel distances over $40^\circ$. Using the Born kernels and a 2D SOLA inversion of travel times, we can recover most features of the input meridional flow profile. The Born approximation is thus a promising method for inferring large-scale solar interior flows.
• ### Stellar magnetic activity and variability of oscillation parameters - An investigation of 24 solar-like stars observed by Kepler(1611.02029)

Nov. 7, 2016 astro-ph.SR
Context. The Sun and solar-like stars undergo activity cycles for which the underlying mechanisms are not well understood. The oscillations of the Sun are known to vary with its activity cycle and these changes provide diagnostics on the conditions below the photosphere. Kepler has detected oscillations in hundreds of solar-like stars but as of yet, no widespread detection of signatures of magnetic activity cycles in the oscillation parameters of these stars have been reported. Aims. We analyse the photometric short cadence Kepler time series of a set of 24 solar-like stars, which were observed for at least 960 days each, with the aim to find signatures of stellar magnetic activity in the oscillation parameters. Methods. We analyse the temporal evolution of oscillation parameters by measuring mode frequency shifts, changes in the height of the p-mode envelope, as well as granulation time scales. Results. For 23 of the 24 investigated stars, we find significant frequency shifts in time. We present evidence for magnetic activity in six of them. We find that the amplitude of the frequency shifts decreases with stellar age and rotation period. For the most prominent example, KIC 8006161, we find that, similar to the solar case, frequency shifts are smallest for the lowest and largest for the highest p-mode frequencies. Conclusions. These findings show that magnetic activity can be routinely observed in the oscillation parameters for solar-like stars, which opens up the possibility to place the solar activity cycle in the context of other stars by asteroseismology.
• 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.
• ### Verification of the helioseismic Fourier-Legendre analysis for meridional flow measurements(1606.05202)

June 16, 2016 astro-ph.SR
Measuring the Sun's internal meridional flow is one of the key issues of helioseismology. Using the Fourier-Legendre analysis is a technique for addressing this problem. We validate this technique with the help of artificial helioseismic data. The analysed data set was obtained by numerically simulating the effect of the meridional flow on the seismic wave field in the full volume of the Sun. In this way, a 51.2-hour long time series was generated. The resulting surface velocity field is then analyzed in various settings: Two $360^\circ \times 90^\circ$ halfspheres, two $120^\circ \times 60^\circ$ patches on the front and farside of the Sun (North and South, respectively) and two $120^\circ \times 60^\circ$ patches on the northern and southern frontside only. We compare two possible measurement setups: observations from Earth and from an additional spacecraft on the solar farside, and observations from Earth only, in which case the full information of the global solar oscillation wave field was available. We find that, with decreasing observing area, the accessible depth range decreases: the $360^\circ \times 90^\circ$ view allows us to probe the meridional flow almost to the bottom of the convection zone, while the $120^\circ \times 60^\circ$ view means only the outer layers can be probed. These results confirm the validity of the Fourier-Legendre analysis technique for helioseismology of the meridional flow. Furthermore these flows are of special interest for missions like Solar Orbiter that promises to complement standard helioseismic measurements from the solar nearside with farside observations.
• ### Experimental observation of $\beta$-delayed neutrons from $^{9}$Li as a way to study short-pulse laser-driven deuteron production(1605.05702)

A short-pulse laser-driven deuteron beam is generated in the relativistic transparency regime and aimed at a beryllium converter to generate neutrons at the TRIDENT laser facility. These prompt neutrons have been used for active interrogation to detect nuclear materials, the first such demonstration of a laser-driven neutron source. During the experiments, delayed neutrons from $^9$Li decay was observed. It was identified by its characteristic half-life of 178.3 ms. Production is attributed to the nuclear reactions $^9$Be(d,2p)$^9$Li and $^9$Be(n,p)$^9$Li inside the beryllium converter itself. These reactions have energy thresholds of 18.42 and 14.26 MeV respectively, and we estimate the (d,2p) reaction to be the dominant source of $^9$Li production. Therefore, only the higher-energy portion of the deuteron spectrum contributes to the production of the delayed neutrons. It was observed that the delayed-neutron yield decreases with increasing distance between the converter and the deuteron source. This behavior is consistent with deuteron production with energy greater than $\sim$20 MeV within a cone with a half-angle greater than 40$^{\circ}$. Prompt-neutron time-of-flight measurements at varying separation between the converter and the laser target indicate that the fast deuteron population above threshold is severely depleted on axis out to $\sim$20$^{\circ}$. These measurements are consistent with emission of the fast deuterons (i.e., above 10 MeV/nucleon) in a ring-like fashion around the central axis. Such an inferred ring-like structure is qualitatively consistent with a documented signature of the breakout afterburner (BOA) laser-plasma ion acceleration mechanism. The measurement of $\beta$-delayed neutrons from $^9$Li decay could provide an important new diagnostic tool for the study of the features of the deuteron production mechanism in a non-intrusive way.
• Theano is a Python library that allows to define, optimize, and evaluate mathematical expressions involving multi-dimensional arrays efficiently. Since its introduction, it has been one of the most used CPU and GPU mathematical compilers - especially in the machine learning community - and has shown steady performance improvements. Theano is being actively and continuously developed since 2008, multiple frameworks have been built on top of it and it has been used to produce many state-of-the-art machine learning models. The present article is structured as follows. Section I provides an overview of the Theano software and its community. Section II presents the principal features of Theano and how to use them, and compares them with other similar projects. Section III focuses on recently-introduced functionalities and improvements. Section IV compares the performance of Theano against Torch7 and TensorFlow on several machine learning models. Section V discusses current limitations of Theano and potential ways of improving it.
• ### Sensitivity Kernels for Flows in Time-Distance Helioseismology: Extension to Spherical Geometry(1604.03803)

April 14, 2016 astro-ph.SR
We extend an existing Born approximation method for calculating the linear sensitivity of helioseismic travel times to flows from Cartesian to spherical geometry. This development is necessary for using the Born approximation for inferring large-scale flows in the deep solar interior. In a first sanity check, we compare two $f-$mode kernels from our spherical method and from an existing Cartesian method. The horizontal and total integrals agree to within 0.3 %. As a second consistency test, we consider a uniformly rotating Sun and a travel distance of 42 degrees. The analytical travel-time difference agrees with the forward-modelled travel-time difference to within 2 %. In addition, we evaluate the impact of different choices of filter functions on the kernels for a meridional travel distance of 42 degrees. For all filters, the sensitivity is found to be distributed over a large fraction of the convection zone. We show that the kernels depend on the filter function employed in the data analysis process. If modes of higher harmonic degree ($90\lesssim l \lesssim 170$) are permitted, a noisy pattern of a spatial scale corresponding to $l\approx 260$ appears near the surface. When mainly low-degree modes are used ($l\lesssim70$), the sensitivity is concentrated in the deepest regions and it visually resembles a ray-path-like structure. Among the different low-degree filters used, we find the kernel for phase-speed filtered measurements to be best localized in depth.
• ### Recent Developments in Helioseismic Analysis Methods and Solar Data Assimilation(1603.04742)

March 15, 2016 astro-ph.SR
We review recent advances and results in enhancing and developing helioseismic analysis methods and in solar data assimilation. In the first part of this paper we will focus on selected developments in time-distance and global helioseismology. In the second part, we review the application of data assimilation methods on solar data. Relating solar surface observations as well as helioseismic proxies with solar dynamo models by means of the techniques from data assimilation is a promising new approach to explore and to predict the magnetic activity cycle of the Sun.
• Neutrinos in the cosmic ray flux with energies near 1 EeV and above are detectable with the Surface Detector array of the Pierre Auger Observatory. We report here on searches through Auger data from 1 January 2004 until 20 June 2013. No neutrino candidates were found, yielding a limit to the diffuse flux of ultra-high energy neutrinos that challenges the Waxman-Bahcall bound predictions. Neutrino identification is attempted using the broad time-structure of the signals expected in the SD stations, and is efficiently done for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for "Earth-skimming" neutrino interactions in the case of tau neutrinos. In this paper the searches for downward-going neutrinos in the zenith angle bins $60^\circ-75^\circ$ and $75^\circ-90^\circ$ as well as for upward-going neutrinos, are combined to give a single limit. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high energy neutrinos with an $E^{-2}$ spectrum in the energy range $1.0 \times 10^{17}$ eV - $2.5 \times 10^{19}$ eV is $E_\nu^2 dN_\nu/dE_\nu < 6.4 \times 10^{-9}~ {\rm GeV~ cm^{-2}~ s^{-1}~ sr^{-1}}$.
• A measurement of the cosmic-ray spectrum for energies exceeding $4{\times}10^{18}$ eV is presented, which is based on the analysis of showers with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger Observatory between 1 January 2004 and 31 December 2013. The measured spectrum confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$ eV, the "ankle", the flux can be described by a power law $E^{-\gamma}$ with index $\gamma=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by a smooth suppression region. For the energy ($E_\text{s}$) at which the spectral flux has fallen to one-half of its extrapolated value in the absence of suppression, we find $E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$ eV.
• ### The FUNK search for Hidden Photon Dark Matter in the eV range(1510.05869)

Oct. 20, 2015 hep-ph, hep-ex, physics.ins-det
We give a brief update on the search for Hidden Photon Dark Matter with FUNK. The experiment uses a large spherical mirror, which, if Hidden Photon Dark Matter exists in the accessible mass and coupling parameter range, would yield an optical signal in the mirror's center in an otherwise dark environment. After a test run with a CCD, preparations for a run with a low-noise PMT are under way and described in this proceedings.
• ### Determination of fundamental asteroseismic parameters using the Hilbert transform(1504.08066)

April 30, 2015 astro-ph.SR, astro-ph.IM
Context. Solar-like oscillations exhibit a regular pattern of frequencies. This pattern is dominated by the small and large frequency separations between modes. The accurate determination of these parameters is of great interest, because they give information about e.g. the evolutionary state and the mass of a star. Aims. We want to develop a robust method to determine the large and small frequency separations for time series with low signal-tonoise ratio. For this purpose, we analyse a time series of the Sun from the GOLF instrument aboard SOHO and a time series of the star KIC 5184732 from the NASA Kepler satellite by employing a combination of Fourier and Hilbert transform. Methods. We use the analytic signal of filtered stellar oscillation time series to compute the signal envelope. Spectral analysis of the signal envelope then reveals frequency differences of dominant modes in the periodogram of the stellar time series. Results. With the described method the large frequency separation $\Delta\nu$ can be extracted from the envelope spectrum even for data of poor signal-to-noise ratio. A modification of the method allows for an overview of the regularities in the periodogram of the time series.
• ### Hidden Photon Dark Matter Search with a Large Metallic Mirror(1410.0200)

Oct. 1, 2014 hep-ph, hep-ex, physics.ins-det
If Dark Matter is composed of hidden-sector photons that kinetically mix with photons of the visible sector, then Dark Matter has a tiny oscillating electric field component. Its presence would lead to a small amount of visible radiation being emitted from a conducting surface, with the photon frequency given approximately by the mass of the hidden photon. Here, we report on experimental efforts that have started recently to search for such hidden photon Dark Matter in the (sub-)eV regime with a prototype mirror for the Auger fluorescence detector at the Karlsruhe Institute for Technology.
• ### An upper bound from helioseismology on the stochastic background of gravitational waves(1401.6888)

March 28, 2014 gr-qc, astro-ph.SR
The universe is expected to be permeated by a stochastic background of gravitational radiation of astrophysical and cosmological origin. This background is capable of exciting oscillations in solar-like stars. Here we show that solar-like oscillators can be employed as giant hydrodynamical detectors for such a background in the muHz to mHz frequency range, which has remained essentially unexplored until today. We demonstrate this approach by using high-precision radial velocity data for the Sun to constrain the normalized energy density of the stochastic gravitational-wave background around 0.11 mHz. These results open up the possibility for asteroseismic missions like CoRoT and Kepler to probe fundamental physics.
• ### On the feasibility of employing solar-like oscillators as detectors for the stochastic background of gravitational waves(1401.6883)

Jan. 27, 2014 gr-qc, astro-ph.SR
We present a hydrodynamic model that describes excitation of linear stellar oscillations by a stochastic background of gravitational waves (SBGW) of astrophysical and cosmological origin. We find that this excitation mechanism is capable of generating solar g-mode amplitudes close to or comparable with values expected from excitation by turbulent convection, which is considered to be the main driving force for solar-like oscillations. A method is presented that places direct upper bounds on the SBGW in a frequency range, in which the SBGW is expected to contain rich astrophysical information. Employing estimates for solar g-mode amplitudes, the proposed method is demonstrated to have the potential to compete with sensitivities reached by gravitational wave experiments in other frequency ranges.
• Contributions of the Pierre Auger Collaboration to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013
• ### The energy spectrum of cosmic rays at the highest energies(1306.6138)

June 26, 2013 astro-ph.IM, astro-ph.HE
One of several working groups established for this workshop was charged with examining results and methods associated with the UHECR energy spectrum. We summarize the results of our discussions, which include a better understanding of the analysis choices made by groups and their motivation. We find that the energy spectra determined by the larger experiments are consistent in normalization and shape after energy scaling factors are applied. Those scaling factors are within systematic uncertainties in the energy scale, and we discuss future work aimed at reducing these systematics.
• ### Disentangling the Air Shower Components Using Scintillation and Water Cherenkov Detectors(1208.2167)

Aug. 10, 2012 astro-ph.HE
We consider a ground array of scintillation and water Cherenkov detectors with the purpose of determining the muon content of air showers. The different response characteristics of these two types of detectors to the components of the air shower provide a way to infer their relative contributions. We use a detailed simulation to estimate the impact of parameters, such as scintillation detector size, in the determination of the size of the muon component.
• ### Mass Composition Sensitivity of an Array of Water Cherenkov and Scintillation Detectors(1109.1692)

Sept. 8, 2011 astro-ph.IM, astro-ph.HE
We consider a hybrid array composed of scintillation and water Cherenkov detectors designed to measure the cosmic ray primary mass composition at energies of about 1 EeV. We have developed a simulation and reconstruction chain to study the theoretical performance of such an array. In this work we investigate the sensitivity of mass composition observables in relation to the geometry of the array. The detectors are arranged in a triangular grid with fixed 750 m spacing and the configuration of the scintillator detectors is optimized for mass composition sensitivity. We show that the performance for composition determination can be compared favorably to that of Xmax measurements after the difference in duty cycles is considered.