O objetivo desta dissertação é a busca de anisotropias em grandes escalas angulares no fluxo de raios cósmicos acima de 4 × 10^{18} eV e o confronto dos resultados obtidos com as previsões dos principais modelos astrofísicos na região de ultra altas energias. Com esse propósito foram analisados os dados de chuveiros atmosféricos extensos registrados pelo Observatório Pierre Auger no período de 01 de janeiro de 2004 até 31 de agosto de 2016, compreendendo um total de 108 480 eventos detetados com ângulos zenitais entre 0 e 80 graus acima de 4 × 10^{18} eV. Efeitos sistemáticos locais como a influência da variação das condições atmosféricas e a influência do campo geomagnético nos estimadores de energia das cascatas atmosféricas foram corrigidos cuidadosamente para evitar a introdução de sinais espúrios de anisotropias na análise. Especial ênfase foi dedicada à construção de um mapa bidimensional na esfera celeste com os valores da exposição para cada direção no céu, calculado com uma precisão de 1% a partir da informação da área de coleção efetiva no Observatório Pierre Auger a cada minuto. Da mesma forma, a distribuição das direções de chegada dos eventos foi registrada em mapas bidimensionais na esfera celeste que, após corrigidos pela exposição, são transformados nos mapas de fluxo usados na análise quantitativa para a busca de anisotropias. Neste trabalho procurou-se padrões dipolares de anisotropias, os quais foram quantificados usando o método de expansão multipolar em harmônicos esféricos. Diferentemente do método bidimensional de Rayleigh usado na maioria dos estudos deste tipo, o método aqui utilizado fornece de maneira direta as três componentes cartesianas do dipolo, por meio das quais a reconstrução da amplitude e orientação espacial do dipolo é feita de maneira mais natural. A busca de anisotropias foi realizada em dois intervalos de energia independentes, assim foram analisados por separado os eventos com energias entre 4 × 10^{18} e 8 × 10^{18} eV e de energias maiores que 8 × 10^{18} eV. No primeiro intervalo não foi encontrado nenhum sinal de anisotropia estatisticamente significativo. No segundo intervalo foi encontrada uma anisotropia dipolar de amplitude d = (7.9 ± 1.1)% com coordenadas equatoriais (_d , _d ) = (103 ± 10 , 38 ± 7 ) graus e uma probabilidade de ocorrência ao acaso de 6.90 × 10^{8} , correspondendo a uma significância estatística acima do limiar de descoberta de 5 . Tal amplitude e direção são melhor explicadas assumindo-se uma origem extragaláctica para os raios cósmicos com energias acima de 8 × 10^{18} eV.

The aim of this dissertation is the search for large scale anisotropies in the cosmic ray flux above 4 × 10^{18} eV and the comparison of the results with the predictions of the main astrophysical models in the ultra high energy region. With this purpose we analyzed the data of extensive air showers detected by the Pierre Auger Observatory in the period of 1st January 2004 to 31st August 2016 comprising a total of 108 480 events with zenith angles between 0 and 80 degrees above the threshold of 4 × 10^{18} eV. Local systematic effects such as the influence of the variation of the atmospheric conditions and the effect of the geomagnetic field in the energy estimators were properly corrected to avoid the introduction of spurious anisotropy signals in the analysis. Special emphasis was placed on the construction of a two-dimensional map on the celestial sphere with the exposure values for each direction on the sky, calculated with an accuracy of 1% from the information of the effective collection area of the observatory at every minute. Likewise, the arrival directions distribution of the events was recorded in two-dimensional maps on the celestial sphere that after correcting by the exposure were transformed into the flux maps to be used in the quantitative analysis for the anisotropies search. In this work we search for dipolar-like anisotropy patterns, which were quantified by performing a spherical harmonics decomposition of the cosmic ray flux map. Unlike the traditional two-dimensional Rayleigh method used in the majority of the studies of this kind, the method used here directly provides the three cartesian components of the dipole, allowing the reconstruction of its amplitude and orientation in a more natural way. The anisotropy search was made in two independent energy bins, so the distributions of the events with energies between 4 × 10^{18} and 8 × 10^{18} eV and those with energies greater than 8 × 10^{18} eV were analyzed separately. In the first energy bin no statistically significant departure from isotropy was found. In the second energy bin a dipole anisotropy of amplitude d = (7.9 ± 1.1)% with equatorial coordinates (_d , _d ) = (103 ± 10 , 38 ± 7 ) degrees and a chance probability of 6.90 × 10^{8} was found, corresponding to a statistical significance above the discovery threshold of 5 . Such amplitude and direction for the reconstructed dipole is better explained if an extragalactic origin is assumed for the bulk of cosmic rays above 8 × 10^{18} eV.