Diffusive shock acceleration operating at expanding supernova remnant shells is by far the most popular model for the origin of galactic cosmic rays. Despite the general consensus received by the model, an unambiguous and conclusive proof of the supernova remnant hypothesis is still missing. In this context, the recent developments in gamma ray astronomy provide us with precious insights into the problem of the origin of galactic cosmic rays, since production of gamma rays is expected both during the acceleration of cosmic rays at supernova remnant shocks and during their subsequent propagation in the interstellar medium. In particular, the recent detection of a number of supernova remnants at TeV energies nicely fits with the model, but it still does not constitute a conclusive proof of it, mainly due to the difficulty of disentangling the hadronic and leptonic contributions to the observed gamma ray emission. The main goal of my research is to search for an unambiguous and conclusive observational test for proving (or disproving) the idea that supernova remnants are the sources of galactic cosmic rays with energies up to (at least) the cosmic ray knee. Our present comprehension of the mechanisms of particle acceleration at shocks and of the propagation of cosmic rays in turbulent magnetic fields encourages beliefs that such a conclusive test might come from future observations of supernova remnants and of the Galaxy in the almost unexplored domain of multi-TeV gamma rays.

Thesis (Ph.D. (Space Physics)--North-West University, Potchefstroom Campus, 2008.

The Voyager 1 spacecraft is now about 25 AU beyond the heliospheric termination shock and soon it should encounter the outer boundary of the heliosphere, the heliopause. This is set to be at 120 AU in the modulation model used for this study. This implies that Voyager 1, and soon afterwards also Voyager 2, should be able to measure the heliopause spectrum, to be interpreted as the lowest possible local interstellar spectrum, for low energy galactic electrons (1 MeV to 120 MeV). This could give an answer to a long outstanding question about the spectral shape (energy dependence) of the galactic electron spectrum at these low energies. These in situ electron observations from Voyager 1, until the year 2010 when it was already beyond 112 AU, are used for a comparative study with a comprehensive three dimensional numerical model for the solar modulation of galactic electrons from the inner to the outer heliosphere. A locally developed steady state modulation model which numerically solves the relevant heliospheric transport equation is used to compute and study modulated electron spectra from Earth up to the heliopause. The issue of the spectral shape of the local interstellar spectrum at these low energies is specifically addressed, taking into account modulation in the inner heliosheath, up to the heliopause, including the effects of the transition of the solar wind speed from supersonic to subsonic in the heliosheath. Modulated electron spectra from the inner to the outer heliosphere are computed, together with radial and latitudinal profiles, focusing on 12 MeV electrons. This is compared to Voyager 1 observations for the energy range 6–14 MeV. A heliopause electron spectrum is computed and presented as a new plausible local interstellar spectrum from 30 GeV down to 10 MeV. The comparisons between model predictions and observations from Voyager 1 and at Earth (e.g. from the PAMELA mission and from balloon flights) and in the inner heliosphere (e.g. from the Ulysses mission) are made. This enables one to make conclusions about diffusion theory applicable to electrons in the heliosphere, in particular the rigidity dependence of diffusion perpendicular and parallel to the local background solar magnetic field. A general result is that the rigidity dependence of both parallel and perpendicular diffusion coefficients needs to be constant below P < 0.4 GV and only be allowed to increase above this rigidity to assure compatibility between the modeling and observations at Earth and especially in the outer heliosphere. A modification in the radial dependence of the diffusion coefficients in the inner heliosheath is required to compute realistic modulation in this region. With this study, estimates of the intensity of low energy galactic electrons at Earth can be made. A new local interstellar spectrum is computed for these low energies to improve understanding of the modulation galactic electrons as compared to previous results described in the literature. / Thesis (M.Sc. (Physics))--North-West University, Potchefstroom Campus, 2012.

Thesis (Ph.D. (Space Physics)--North-West University, Potchefstroom Campus, 2008.

The Voyager 1 spacecraft is now about 25 AU beyond the heliospheric termination shock and soon it should encounter the outer boundary of the heliosphere, the heliopause. This is set to be at 120 AU in the modulation model used for this study. This implies that Voyager 1, and soon afterwards also Voyager 2, should be able to measure the heliopause spectrum, to be interpreted as the lowest possible local interstellar spectrum, for low energy galactic electrons (1 MeV to 120 MeV). This could give an answer to a long outstanding question about the spectral shape (energy dependence) of the galactic electron spectrum at these low energies. These in situ electron observations from Voyager 1, until the year 2010 when it was already beyond 112 AU, are used for a comparative study with a comprehensive three dimensional numerical model for the solar modulation of galactic electrons from the inner to the outer heliosphere. A locally developed steady state modulation model which numerically solves the relevant heliospheric transport equation is used to compute and study modulated electron spectra from Earth up to the heliopause. The issue of the spectral shape of the local interstellar spectrum at these low energies is specifically addressed, taking into account modulation in the inner heliosheath, up to the heliopause, including the effects of the transition of the solar wind speed from supersonic to subsonic in the heliosheath. Modulated electron spectra from the inner to the outer heliosphere are computed, together with radial and latitudinal profiles, focusing on 12 MeV electrons. This is compared to Voyager 1 observations for the energy range 6–14 MeV. A heliopause electron spectrum is computed and presented as a new plausible local interstellar spectrum from 30 GeV down to 10 MeV. The comparisons between model predictions and observations from Voyager 1 and at Earth (e.g. from the PAMELA mission and from balloon flights) and in the inner heliosphere (e.g. from the Ulysses mission) are made. This enables one to make conclusions about diffusion theory applicable to electrons in the heliosphere, in particular the rigidity dependence of diffusion perpendicular and parallel to the local background solar magnetic field. A general result is that the rigidity dependence of both parallel and perpendicular diffusion coefficients needs to be constant below P < 0.4 GV and only be allowed to increase above this rigidity to assure compatibility between the modeling and observations at Earth and especially in the outer heliosphere. A modification in the radial dependence of the diffusion coefficients in the inner heliosheath is required to compute realistic modulation in this region. With this study, estimates of the intensity of low energy galactic electrons at Earth can be made. A new local interstellar spectrum is computed for these low energies to improve understanding of the modulation galactic electrons as compared to previous results described in the literature. / Thesis (M.Sc. (Physics))--North-West University, Potchefstroom Campus, 2012.

Après plus d'un siècle d'étude des rayons cosmiques, des questions persistent à propos de leur origine. À l'échelle de notre galaxie, les vestiges de supernovæ sont d'excellents candidats en tant que sources de rayons cosmiques et il est aujourd'hui avéré que ces ondes de choc générées par l'explosion d'étoiles en fin de vie accélèrent des électrons jusqu'au TeV. Cependant, des preuves tangibles manquent toujours pour confirmer que les protons, qui représentent ~90% du rayonnement cosmique, y sont également accélérés. Ces preuves sont à chercher entre 100 MeV et 10 GeV où la forme spectrale de l'émission gamma nous renseigne sur la nature des particules accélérées.Ce domaine en énergie est observé par le satellite Fermi qui scrute le ciel à haute énergie grâce au Large Area Telescope (LAT), un instrument capable de détecter les rayons gamma d'énergies comprises entre 20 MeV et plus de 300 GeV. Avec la dernière version des données en date (Pass 8) rendue publique en 2015, la surface efficace a été sensiblement accrue, notamment à haute énergie. L'objet de cette thèse a été de mettre à profit les nouvelles performances du LAT pour étudier le rayonnement provenant de jeunes vestiges de supernovæ en coquille afin d'évaluer leur capacité à accélérer des particules à haute énergie. Au total, quatre vestiges de ce type ont été étudiés en détail : RCW 86, HESS J1731-347, SN 1006 et RX J1713.7-3946.L'analyse morphologique et spatiale de ces sources nous a menés à d'importants résultats : HESS J1731-347 et SN 1006 ont été détectés pour la première fois au GeV, nous avons mis en évidence une morphologie en coquille pour RCW 86 qui était précédemment détectée comme une source ponctuelle et RX J1713.7-3946 a révélé une forme spectrale plus complexe que ce qui était alors connu. Grâce à la modélisation de l'émission non-thermique de ces sources, en s'aidant des données au TeV obtenues par H.E.S.S., nous avons montré que leur émission gamma est très largement dominée par le rayonnement Compton inverse des électrons de haute énergie diffusant sur les champs de photons ambiants. Cependant, des modifications spectrales ont également été trouvées dans certaines parties de SN 1006 et RX J1713.7-3946 où l'onde de choc est en interaction avec des régions denses en matière, traçant la possible accélération de protons en leur sein. / After more than a century of cosmic-ray studies, questions about their origin remain unanswered. In our galaxy, supernova remnants are excellent candidates to be the sources of cosmic rays and it is now certain that these shock waves created by the explosion of dying stars accelerate electrons up to TeV energies. However, undeniable proof is still missing to confirm that protons, which represent ~90% of the cosmic radiation, are also accelerated. Such proof can be searched for between 100 MeV and 10 GeV, where the spectral shape of the gamma-ray emission changes according to the nature of the particles that are accelerated.This energy range is covered by the Fermi spacecraft which observes the high energy sky with the Large Area telescope (LAT), an instrument detecting gamma-rays from 20 MeV to more than 300 GeV. With the latest version of the data, released in 2015 (Pass 8), the effective area has been greatly improved, especially at high energy. The goal of this thesis was to take advantage of the improved performance of the LAT to study the radiation coming from young shell-type supernova remnants and to evaluate their ability to accelerate particles to high energy. Four such remnants were studied in detail: RCW 86, HESS J1731-347, SN 1006 and RX J1713.7-3946.Spatial and spectral analyses of these sources yielded important results: HESS J1731-347 and SN 1006 were identified for the first time in the GeV range, we detected a shell-like morphology for RCW 86 which was previously seen as a point source, and RX J1713.7-3946 revealed an unexpected spectral shape. Broadband modeling of the non-thermal emission of these remnants, using in particular TeV data obtained with H.E.S.S., showed that their gamma-ray emission is dominated by the inverse Compton scattering of electrons on ambient photon fields. However, we also found spectral modifications in some parts of SN 1006 and RX J1713.7-3946 where the shock is interacting with dense regions that could trace the acceleration of protons.

No estudo dos raios cósmicos de ultra alta energia utilizando detectores de superfície a energia da partícula primária é estimada pela distribuição lateral (LDF - Lateral Distribution Function), que descreve a amplitude do sinal das estações em função da distância ao centro do chuveiro. Entretanto, com exceção da estimativa do centro do chuveiro, não se utiliza a LDF para obter mais nenhuma informação sobre o chuveiro, talvez porque ela não possua uma parametrização que a descreva completamente, especialmente para os chuveiros com energias mais altas. As primeiras interações dos raios cósmicos com a atmosfera são determinantes para o desenvolvimento dos chuveiros atmosféricos extensos. Tais interações dependem, entre outras coisas, da composição química dos raios cósmicos. Diferenças nessas interações podem causar alterações nas flutuações da distribuição lateral. Através de simulações dos chuveiros com diferentes partículas primárias pode ser possível estimar a composição dos raios cósmicos de ultra alta energia, comparando as flutuações das distribuições laterais de eventos reais com as de eventos simulados. Uma das grandezas relevantes para a flutuação da LDF é a incerteza do sinal das estações. O framework de análise do Observatório Pierre Auger aplica uma correção parametrizada empiricamente à incerteza do sinal das estações. Neste trabalho foi apresentada uma justificativa estatística para esta correção, que está relacionada à distribuição do sinal não ser uma Poisson, mas uma composição de processos com diferentes distribuições. Para a realização deste trabalho foi gerada uma biblioteca de chuveiros produzidos por dois simuladores de chuveiros atmosféricos, AIRES e CORSIKA, sendo que o AIRES utilizou o SIBYLL como modelo de interação hadrônica enquanto o CORSIKA utilizou o EPOS. Foram produzidos chuveiros iniciados por prótons e núcleos de ferro com os dois simuladores, e suas distribuições angulares foram consideradas isotrópicas. A distribuição de energia desses eventos segue uma lei de potência e varia entre 1 e 200 EeV. Utilizando a função de Nishimura, Kamata e Greisen (NKG) como parametrização para a LDF se observa resíduos sistematicamente positivos nas estações mais distantes do centro do chuveiro, que também têm o sinal mais próximo do trigger. Uma das hipóteses levantadas em outros trabalhos para esse comportamento é que ele estaria relacionado à influência das estações silenciosas, mas este trabalho mostra que o uso dessas estações tem pouca influência na flutuação da LDF. Na verdade esse efeito é causado porque as parametrizações da LDF no Offline não consideram que os sinais das estações possuem um corte devido o trigger, isto é, que a função de densidade de probabilidade que descreve o sinal real não é a mesma que descreve o sinal que é observado. Neste trabalho é proposta uma correção para as parametrizações da LDF que é implementada no Offline. Como resultado desta correção é observada uma redução significativa dos resíduos que eram sistematicamente positivos e que, após a correção, passam a ser compatíveis com zero. Neste trabalho foram realizadas três análises independentes para comparar os eventos reais com os simulados, das quais duas não dependem diretamente do ajuste da LDF e também não são sensíveis à energia. Elas permitem uma comparação entre os sinais supondo uma relação simples entre eles. No primeiro caso supõe-se que a diferença no sinal é devido à componente muônica do chuveiro e no segundo supõe-se que os dois conjuntos de eventos comparados são bem descritos por funções NKG mas com parâmetros S1000 diferentes. A terceira análise utiliza os resíduos dos ajustes da LDF e permite observar a composição em função da energia da partícula primária. Essa última análise foi realizada utilizando a função NKG com e sem correção do efeito do trigger. As diferentes análises utilizadas para estimar a composição dos raios cósmicos apresentaram resultados consistentes entre si, apesar das limitações encontradas em algumas delas. Todos esses indicadores de composição da partícula primária obtidos pelo detector de superfície são consistentes com os resultados obtidos pelas análises de Xmax do detector de fluorescência, reforçando a tese de que a composição dos raios cósmicos é predominantemente de próton entre 1 e 10 EeV e entre próton e ferro para energias acima de aproximadamente 10 EeV. / The energy of ultra high energy cosmic rays can be estimated from the lateral distribution function (LDF) of the shower as measured by surface detectors. The LDF describes the particle density as a function of the distance from the shower center. However, with the exception of the position of the shower center, no other information is extracted from it, may because it does not have a parametrization or an analytic function that describes it completely. The first interactions of cosmic rays with the atmosphere are decisive for the development of the extensive air showers. Such interactions, among other things, depend on the chemical composition of comic rays. Differences in these interactions can cause changes in the fluctuation shape of lateral distribution. Through simulations of showers with different primary particles it may be possible to estimate the composition of ultra high energy cosmic rays comparing the fluctuation shape of the lateral distributions of real events with those from simulated ones. One of the quantities relevant to the fluctuation of the LDF signal is the uncertainty of the stations. The analysis framework of the Pierre Auger Observatory applies a correction to the signal uncertainty of the signal. The parameterization of this correction is obtained empirically. In this work a statistical justification for this correction is proposed and is related to distribution of the signal which is not Poisson, but a composition of processes with different distributions. For this work a library of showers using two simulators of air showers, AIRES and CORSIKA, was produced. The showers simulated with the AIRES used SIBYLL as a hadronic interaction model while COSIKA used EPOS. Showers initiated by protons and iron nuclei with the two simulators were produced, and their angular distribution was considered isotropic. The energy distribution of these events follows a power law and ranges from 1 to 200 EeV. Using the Nishimura, Kamata and Greisen (NKG) function as a parameterization for the LDF, one obtains residues that are systematically positive at stations further from the center of the shower. These stations have a signalclose to the trigger level. One of the hypothesis raised in other works for this behavior is that it is related to the influence of the silent stations, but this work shows that their use has little impact on the fluctuation shape of the LDF. In fact, this effect is caused because the parametrizations of LDF in the Offline ignore that the signals of the stations have a cut due to the trigger, ie, the probability density function that describes the real signal is not the same that describes the observed signal. This work proposes a correction to the parameterizations of the LDF and implements it in the Offline. As a result of this correction, the residues, which were always positive, are significantly reduced and compatible with zero. In this study three independent analysis were performed to compare real and simulated events, two of them not dependent directly of the LDF fit and also not sensitive to the primary particle energy. They allow a comparison between the signals assuming a simple relationship between them. The first case assumes that the difference in signal is due to the muonic component of the shower and the second assumes that the two compared sets of events are well described by NKG functions but with different S1000. The third analysis uses the residues of the LDF fits and is able to observe the composition of as a function of primary particle energy. This last analysis was performed using the NKG function with and without correction of the trigger effect. The different analysis used to estimate the composition of cosmic rays showed results consistent, despite the limitations found in some of them. The primary particle composition obtained from the surface detectors in this work is consistent with the results derived from the elongation rate measured by the fluorescence detectors, supporting the hypotesis that the composition of cosmic rays is predominantly proton becoming heavier for energies above 10 EeV.

Relativity is one of the most important and well tested theories and Lorentz invariance is one of its pillars. Lorentz invariance violation (LIV), however, has been discussed in several quantum gravity and high energy models. For this reason, it is crucial to test it. Several tests, both terrestrial and astrophysical, have been performed in the last years and provide limits on the violation. This work takes part in these efforts and discuss the possibility of testing LIV with ultra-high energy cosmic rays (UHECRs). The effects of LIV in their propagation and the resulting changes in the spectrum of UHECRs are obtained and compared to the experimental data from the Pierre Auger Observatory. An analytical calculation for the inelasticity in the laboratory frame with LIV of any a + b &rarr; c + d interaction is presented and used to obtain the phase space and the energy losses of the pion production for protons, the photodisintegration for nuclei and the pair production for photons with LIV. A parametrization for the threshold energy of the photodisintegration with LIV is also proposed. The main effect seen is a decrease in the phase space and a resulting decrease in the energy loss. These changes have been implemented in Monte Carlo propagation codes and the resulting spectra of protons, nuclei and photons on Earth have been obtained and fitted to the data from the Pierre Auger Observatory. It is shown that upper limits on the photon LIV coefficient can be derived from the upper limits on the photon flux from the Pierre Auger Observatory. / Relatividade é uma das mais importantes e bem testadas teorias e a invariância de Lorentz é um de seus pilares. A violação da invariância de Lorentz (VIL), todavia, tem sido discutida em diversos modelos de gravidade quântica e altas energias. Por tal motivo, é crucial testá-la. Diversos testes, tanto terrestres quanto astrofísicos, foram realizados nos últimos anos e fornecem limites na violação. Este trabalho se insere nesses esforços e discute a possibilidade de testar VIL com raios cósmicos de altíssima energia. Os efeitos da VIL em sua propagação e as consequentes mudanças no espectro de raios cósmicos de altíssima energia são obtidos e comparados com os dados experimentais do Observatório Pierre Auger. Um cálculo analítico para a inelasticidade no referencial do laboratório com VIL para qualquer interação da forma a + b &rarr; c + d é apresentado e usado para obter o espaço de fase e as perdas de energia para a produção de píons para prótons, a fotodesintegração para núcleos e a produção de pares para fótons com VIL. Uma parametrização para o limiar de energia da fotodesintegração com VIL também é proposta. O principal efeito observado é uma diminuição no espaço de fase e uma consequente diminuição nas perdas de energia. Tais mudanças foram implementadas em códigos de Monte Carlo para a propagação e os espectros resultantes para prótons, núcleos e fótons na Terra foram obtidos e ajustados aos dados do Observatório Pierre Auger. É mostrado que limites superiores nos coeficientes de VIL para o fóton podem ser deduzidos dos limites superiores para o fluxo de fótons do Observatório Pierre Auger.

Made available in DSpace on 2016-03-15T19:35:57Z (GMT). No. of bitstreams: 1 Edith Tueros Cuadros.pdf: 4357925 bytes, checksum: 7a9a7c9700e9ab1b8e24a65ab089b198 (MD5) Previous issue date: 2016-02-02 / Fundação de Amparo a Pesquisa do Estado de São Paulo / Cosmic rays are strongly influenced by solar, geomagnetic and atmospheric phenomena. CARPET detector, conceived for cosmic rays observation with energies in the range between 105 - 1012 eV, is an important tool for the study of these phenomena. The Earth s atmosphere conditions are also affected by changes in the cosmic rays flux, therefore, cosmic rays characterization is important to define physical and chemical conditions of our atmosphere. To characterize the cosmic rays flux variations, detected on the ground, prior elimination of atmospheric pressure and temperature effects on ground level is needed, thus, data recorded by meteorological instruments on CASLEO were used for that corrections. To eliminate the effect of temperature through the whole vertical atmosphere it was applied the integral and the mass-average temperature method by using vertical temperature profiles. Both methods were tested using CARPET-TEL data for the year 2009, this data were previously corrected by pressure influences. The mass-average temperature method shows a better response when comparing the corrected CARPET date with neutron monitor observations .The whole cosmic rays flux data analysis, for CARPET-TEL data corrected by integral method (for the period 2006/04/01 - 2014/06/30), shows an anti-correlation with sunspot number and a clear seasonal variation after 2008. Two Forbush decreases (FD), were detected by CARPET which were produced by geo-effective CMEs. Both FD onsets coincided with the interplanetary shock instant, which preceded intense geomagnetic storms. With these results, we can reinforce that CARPET is an important tool to study long and short term cosmic rays behavior, because it has a similar response to experiments that operates at atmosphere particle energy ranges. / Os raios cósmicos são fortemente influenciados pelos fenômenos solares, geomagnéticos e atmosféricos. O detector CARPET, concebido para a observação dos raios cósmicos com energias na faixa compreendida entre 105 e 1012 eV, é uma ferramenta importante para o estudo desses fenômenos. As condições da nossa atmosfera são também afetadas pelas variações no fluxo de raios cósmicos que chegam a Terra, portanto sua caracterização é importante para se definir as condições físicas e químicas da mesma. Para caracterizar as variações do fluxo de raios cósmicos a partir de observações no solo é necessário a eliminação prévia da influência dos efeitos da pressão e temperatura atmosférica na superfície, para tal foram utilizados dados meteorológicos obtidos no CASLEO. Para a eliminação do efeito da temperatura ao longo de toda atmosfera foram utilizados o método integral e o método de temperatura ponderada pela massa que utilizam perfis verticais de temperatura. Os dois métodos foram testados para o ano de 2009 com os dados do canal TEL do CARPET previamente corrigidos pela pressão, sendo que o método de temperatura ponderada pela massa apresentou uma resposta melhor ao se comparar os dados corrigidos com os dados observados com monitores de nêutrons. A análise de todos os dados do fluxo raios cósmicos do canal TEL do CARPET (1/04/2006 a 30/06/2014) corrigidos pelo método integral, mostra a anti-correlação com o número de manchas solares, e uma variação sazonal pronunciada após 2008. Foram detectados dois decréscimos Forbush (FD) produzidos por CMEs geo-efetivos. O início de ambos FD coincidiu com o instante do choque interplanetário, que antecedeu tempestades geomagnéticas intensas. Com estes resultados reforçamos que o CARPET é uma ferramenta importante para estudo das modulações de raios cósmicos de longo e curto prazo porque apresenta uma resposta semelhante com os experimentos que fazem observações em outras faixas de energia.

Ultra-High Energy Cosmic Rays (UHECR) can only be measured by the detection of Extensive Air Showers (EAS) created by the interaction of the cosmic ray particle with an atmospheric nuclei. The inference of some of the properties of UHECR, like their mass composition, is only possible by the comparison of measurements of EAS observables to predictions from Monte Carlo simulations. The most important source of uncertainties on the description of EAS by the simulations is the modeling of hadronic interactions. For many years it has been known that the hadronic interaction models fail on predicting the EAS observables related to their muonic component. The most evident manifestation of that is called muon deficit problem due to the fact that the number of muons in EAS with energies above 1018 eV predicted by simulations is smaller than the observed ones. The aim of this thesis is to approach this problem in three distinct fronts. First, a method is developed to interpret measurements of number of muons in terms of cosmic rays composition in despite of the muon deficit problem. Second, an EAS observable which is sensitive to the muon energy spectrum at ground and, consequently, can be used to constrain hadronic interaction models is proposed and tested. Third and final, the muon production in air showers is studied through measurements of hadron production spectra in pion-carbon interactions. / Raios Cósmicos Ultra Energéticos (Ultra-High Energy Cosmic Rays, UHECR) somente podem ser medidos através da detecção dos Chuveiros Atmosféricos Extensos (Extensive Air Showers, EAS) criados pela interação do raio cósmico primário com núcleos atmoféricos. A inferência de algumas propriedados dos UHECRs, como a composição de massa, é possível somente através da comparação entre medidas de observáveis dos EASs com predições geradas por simulações de Monte Carlo. A fonte de incerteza mais importante na descrição de EAS por simulações é a modelagem das interações hadrônicas. Por muitos anos é sabido que os modelos de interação hadrônica falham na predição de observáveis dos EASs relacionados a sua componente muônica. A manifestação mais evidente disso é chamada problema do déficit de múons devido ao fato que o número de múons em chuveiros com energias acima de 1018 eV predito por simulações é menor que os observados. O objetivo desta tese é abordar este problema através de três frentes. Primeiramente, um método é desenvolvido para interpretar as medidas do número de múons em termos de composição de raios cósmicos considerando o problema do déficit de múons. Segundo, a proposta e o teste de um observável que é sensível ao espectro de energia dos múons na superfície e, consequentemente, pode ser usado para discriminar entre os modelos de interação hadrônica. Por último, a produção de múons em chuveiros é estudada através de medidas do espectro de produção de hádrons em interações do tipo píon-carbono.