Recently, the ALICE Collaboration reported an enhancement of the yield ratio of strange and multi-strange hadrons to charged pions as a function of multiplicity at mid-rapidity in proton-proton, proton-lead, lead-lead, and xenon-xenon scattering. ALICE observations provide a strong indication that a quark-gluon plasma is partly formed in high multiplicity events of both small and large colliding systems. Motivated by ALICE's results, we propose a new test for hadronic interaction models used for analyzing ultra-high-energy-cosmic-ray (UHECR) collisions with air nuclei. The test is grounded in the almost equal column-energy density in UHECR-air collisions and lead-lead collisions at the LHC. We applied the test to post-LHC event generators describing hadronic phenomena of UHECR scattering and show that these QCD Monte Carlo-based codes must be retuned to accommodate the strangeness enhancement relative to pions observed in LHC data.

IceCube has observed 80 astrophysical neutrino candidates in the energy range 0.02 < E_\nu/PeV < 2. Deep inelastic scattering of these neutrinos with nucleons on Antarctic ice sheet probe center-of-mass energies $\sqrt{s} \sim$ 1 TeV. By comparing the rates for two classes of observable events, any departure from the benchmark (perturbative QCD) neutrino-nucleon cross section can be constrained. Using the projected sensitivity of South Pole next generation neutrino telescope we show that this facility will provide a unique probe of strong interaction dynamics. In particular, we demonstrate that the high-energy high-statistics data sample to be recorded by IceCube-Gen2 in the very near future will deliver a direct measurement of the neutrino-nucleon cross section at $\sqrt{s} \sim 1$ TeV, with a precision comparable to perturbative QCD informed by HERA data. We also use IceCube data to extract the neutrino-nucleon cross section at $\sqrt{s} \sim 1$ TeV through a likelihood analysis, considering (for the first time) both the charged-current and neutral-current contributions as free parameters of the likelihood function.

The Antarctic Impulsive Transient Antenna (ANITA) experiment has observed two air shower events with energy $\sim 500~{\rm PeV}$ emerging from the Earth with exit angles $\sim 30^\circ$ above the horizon. As was immediately noted by the ANITA Collaboration, these events (in principle) could originate in the atmospheric decay of an upgoing $\tau$-lepton produced through a charged current interaction of a $\nu_\tau$ inside the Earth. However, the relatively steep arrival angles of these perplexing events create tension with the standard model (SM) neutrino-nucleon interaction cross section. Deepening the conundrum, the IceCube neutrino telescope and the Pierre Auger Observatory with substantially larger exposures to cosmic $\nu_\tau$'s in this energy range have not observed any events. This lack of observation implies that the messenger particle (MP) giving rise to ANITA events must produce an air shower event rate at least a factor of 40 larger than that produced by a flux of $\tau$-neutrinos to avoid conflicts with the upper limits reported by the IceCube and the Pierre Auger collaborations. In addition, the sensitivity of ANITA to MP-induced events must be comparable to or larger than those of IceCube and Auger to avoid conflict with the non-observation of any signal at these facilities. Beyond SM interpretations of ANITA events can be classified according to whether the MPs: (i) live inside the Earth, (ii) originate in neutrino-nucleon collisions inside the Earth, (iii) come from cosmological distances. In this communication we investigate the positive and negative facets of these three classes of models.

We find that the ratio r_{mu e} of the muon to the electromagnetic component of an extended air shower at the ground level provides an indirect measure of the depth X_{max} of the shower maximum. This result, obtained with the air-shower code AIRES, is independent of the hadronic model used in the simulation. We show that the value of r_{mu e} in a particular shower discriminates its proton or iron nature with a 98% efficiency. We also show that the eventual production of forward heavy quarks inside the shower may introduce anomalous values of r_{mu e} in isolated events.

The total photoproduction cross section at ultra-high energies is obtained using a model based on QCD minijets and soft-gluon resummation and the ansatz that infrared gluons limit the rise of total cross sections. This cross section is introduced into the Monte Carlo system AIRES to simulate extended air-showers initiated by cosmic ray photons. The impact of the new photoproduction cross section on common shower observables, especially those related to muon production, is compared with previous results.

A model for the total photoproduction cross section based on the ansatz that resummation of infrared gluons limits the rise induced by QCD minijets in all the total cross-sections, is used to simulate extended air showers initiated by cosmic rays with the AIRES simulation program. The impact on common shower observables, especially those related with muon production, is analysed and compared with the corresponding results obtained with previous photoproduction models.

Very energetic charm and bottom hadrons may be produced in the upper atmosphere when a primary cosmic ray or the leading hadron in an extensive air shower collide with a nucleon. At $E\approx 10^8$ GeV their decay length becomes of the order of 10 km, implying that they tend to interact in the air instead of decaying. Since the inelasticity in these collisions is much smaller than the one in proton and pion collisions, there could be rare events where a heavy-hadron component transports a significant amount of energy deep into the atmosphere. We have developed a module for the detailed simulation of these processes and have included it in a new version of the air shower simulator AIRES. We study the frequency, the energy distribution and the depth of charm and bottom production, as well as the depth and the energy distribution of these quarks when they decay. As an illustration, we consider the production and decay of tau leptons (from $D_s$ decays) and the lepton flux at PeV energies from a 30 EeV proton primary. The proper inclusion of charm and bottom hadrons in AIRES opens the possibility to search for air-shower observables that are sensitive to heavy quark effects.

Recent data from the Auger Observatory on low energy secondary cosmic ray particles are analyzed to study temporal correlations together with data on the daily sunspot numbers and neutron monitor data. Standard spectral analysis demonstrates that the available data shows 1/f^{\beta} fluctuations with {\beta} approximately 1 in the low frequency range. All data behave like Brownian fluctuations in the high frequency range. The existence of long-range correlations in the data was confirmed by detrended fluctuation analysis. The real data confirmed the correlation between the scaling exponent of the detrended analysis and the exponent of the spectral analysis.

A comparative analysis of the secondary particles output of the main hadronic interaction packages used in simulations of extensive air showers is presented. Special attention is given to the study of events with very energetic leading secondary particles, including diffractive interactions.

We emphasize the inelasticity distribution of events detected at the IceCube neutrino telescope as an important tool for revealing new physics. This is possible because the unique energy resolution at this facility allows to separately assign the energy fractions for emergent muons and taus in neutrino interactions. As a particular example, we explore the possibility of probing second and third generation leptoquark parameter space (coupling and mass). We show that production of leptoquarks with masses \agt 250 GeV and diagonal generation couplings of O(1) can be directly tested if the cosmic neutrino flux is at the Waxman-Bahcall level.

The variational method with constraints recently developed by Verkley and Gerkema to describe maximum-entropy atmospheric profiles is generalized to ideal gases but with temperature-dependent specific heats. In so doing, an extended and non standard potential temperature is introduced that is well suited for tackling the problem under consideration. This new formalism is successfully applied to the atmosphere of Venus. Three well defined regions emerge in this atmosphere up to a height of $100 km$ from the surface: the lowest one up to about $35 km$ is adiabatic, a transition layer located at the height of the cloud deck and finally a third region which is practically isothermal.

Monthly averages of the sunspot number visible on the sun, observed from 1749, Zurich Observatory, and from 1848 other observatories, have been analyzed. This time signal presents a frequency power spectra with a clear $1/f^\alpha$ behavior with $\alpha\simeq 0.8 \pm 0.2$. The well known cycle of approximately 11 years, clearly present in the spectrum, does not produce a sensible distortion of that behavior. The eventual characterization of the sunspot time series as a fractal is analyzed by means of the detrended fluctuation analysis. The jump size distribution of the signal is also studied.

An analysis of the influence of the Landau-Migdal-Pomeranchuk (LPM) effect on the development of air showers initiated by astroparticles is presented. The theory of Migdal is studied and compared with other theoretical methods, particularly the Blankenbecler and Drell approach. By means of realistic computer simulations and using algorithms that emulate Migdal's theory, including also the so-called dielectric suppression, we study the behavior of the relevant observables in the case of ultra-high energy primaries. We find that the LPM effect can significantly modify the development of high energy electromagnetic showers in certain cases.

An analysis of the influence of the Landau-Migdal-Pomeranchuck (LPM) effect on the development of air showers initiated by astroparticles, is presented. By means of computer simulations using algorithms that emulate Migdal's theory, including also the so-called dielectric suppression, we study the behaviour of the relevant observables in the case of ultra high energy primaries. We find that the LPM effect can significantly modify the development of high energy electromagnetic showers.