I. INTRODUCCIÓN
Luego de su propagación en el espacio intergaláctico, gran cantidad de astropartículas atraviesan los campos magnéticos en la heliósfera y la magnetósfera terrestre. Estas finalmente arriban a la atmósfera terrestre e interactúan con los núcleos de la materia que la componen. En estas interacciones, pueden generarse hasta miles de millones de partículas secundarias, que continúan interactuando y decayendo, resultando en una “cascada” o “lluvia” de secundarios. A estas se las denominan genéricamente como EAS (por sus siglas en inglés de Extensive Air Showers). Muchos de estos secundarios arriban a la superficie terrestre y son una componente importante de la radiación natural de nuestro ambiente. Experimentos orientados a la búsqueda de materia oscura o al estudio de la física de neutrinos, requieren de un blindaje para reducir el ruido intrínseco producido por esta radiación de origen extraterrestre. Hoy en día se encuentran en funcionamiento más de una docena de laboratorios subterráneos, todos ellos en el hemisferio norte, construidos para albergar experimentos que deben ser realizados en ambientes con bajos niveles de radiación.
ANDES1 (por sus siglas en inglés de Agua Negra Deep Experiment Site) es un proyecto propuesto por un grupo de científicos latinoamericanos en 2012 y se trata de la construcción de un laboratorio subterráneo para albergar experimentos de materia oscura y detección de neutrinos además de estudios en otras áreas como biología y geología. Este laboratorio será instalado dentro del túnel de Agua Negra en los Andes a una profundidad de 1700 m bajo roca y albergará experimentos de colaboraciones de todo el mundo de diversas disciplinas como física, geología, sismología y estudios del medio ambiente. En el caso de ANDES, además de su profundidad, cuenta con el beneficio de su ubicación en el hemisferio sur, la cual es valiosa por diversos motivos, por ejemplo, en búsquedas indirectas de materia oscura ayudaría a eliminar el ruido de fondo estacional observado anualmente en modulaciones de experimentos en el hemisferio norte .2 Por lo tanto, es necesario estimar con precisión el fondo esperado de muones presente en el interior del laboratorio subterráneo ANDES.
Una estimación previa de este flujo fue realizada por miembros de la colaboración a cargo del proyecto ANDES. En ese estudio se concluyó que el flujo esperado sería de 3. Sin embargo, al momento de realizar esta estimación el equipo de ANDES aún no contaba con los estudios estructurales ni geológicos de la composición y distribución de roca de la montaña que albergará a ANDES y que resultan necesarios para una estimación más precisa.
El objetivo de este trabajo es actualizar, refinar los modelos y perfeccionar los métodos de cálculo utilizados para realizar una mejor estimación del fondo de muones en ANDES. El laboratorio se encontrará ubicado dentro de la formación geológica denominada Grupo Choiyoi, compuesta en su mayoría por riolitas porfídicas y dacita, ambas rocas volcánicas formadas esencialmente por cuarzo, plagioclasas, feldespata y biotita. Esta información fue aportada por un estudio geológico preliminar de la estructura de la montaña4 a partir del cual se planteó el modelo de caminos en la roca que se muestra en la Fig. 1.
II. MÉTODOS
Haciendo uso de los software de simulación CORSIKA5 y ARTI6 (ampliamente validados para este tipo de estudios7) se realizaron simulaciones de flujo de secundarios a dos alturas diferentes, el pico y la base de la montaña.
Luego, a partir de la definición de opacidad de la roca:
se utilizó esta magnitud para estimar la atenuación del flujo de secundarios de muy alta energía, con especial énfasis en los muones. Para esto se estudió su propagación en las direcciones indicadas en la Fig. 1 dentro de la montaña. Asimismo, se supuso simetría acimutal, y se diseñó una función de transferencia. Para esto se utilizó la aproximación de desaceleración continua8 (CSDA por sus siglas en inglés) ampliamente utilizada en experimentos como IceCube9. En ella se supone que la tasa de pérdida de energía de las partículas que pasan a través de un material es igual a su poder de frenado . Este está definido en cada porción diferencial del material. En particular, la pérdida de energía para los muones de altas energías puede ser descripta como:
siendo la pérdida de energía electrónica y la contribución en escala de energías proveniente de procesos radiativos -Bremsstrahlung, producción de pares e interacciones fotonucleares- y E la energía total.
Simulaciones de lluvias atmosféricas extendidas
Las simulaciones del flujo de secundarios sobre el sitio de ANDES fueron realizadas con el programa basado en el método de Montecarlo CORSIKA10 (versión 7.5600). A su vez se hizo uso de ARTI11 que consta de un paquete completo que convierte a CORSIKA, que esá hecho para simular lluvias monocromáticas, en un programa capaz de simular el flujo de fondo generado en un intervalo de tiempo. Este paquete fue creado en el marco del proyecto LAGO12.
De esta forma, se simuló el flujo de partículas secundarias a dos alturas distintas que denominamos “pico”: a 5400 m s.n.m (metros sobre el nivel del mar) y “base”: a 3700 m s.n.m. sobre el sitio de ANDES, utilizando la correspondiente selección de parámetros, a saber: modelo hadrónico a altas energías: QGSJETII-04 (Quark Gluon String Model with JETs); dirección de arribo del primario: entre y para el ángulo cenital y el ángulo acimutal entre - y ; límites de energía del primario: 1 TeV - 10 TeV; altura del sitio: base (4700 m s.n.m) y pico (5300 m s.n.m) de la montaña; tiempo de flujo: dos simulaciones equivalentes a tres meses de flujo, para cada altura con el fin de tener una buena cobertura estadística en la fracción de altas energías teniendo en cuenta la fuerte dependencia del flujo de primarios con la energía (); límite inferior de energía de secundarios: 800 GeV a fin de optimizar el tiempo de simulación; tipo de detector: detector volumétrico (el flujo sigue una distribución del tipo ).
En primer lugar, se integró el flujo de todos los primarios en el rango de (protones) a (hierros) en un tiempo total de 90 días y se aproximó el espectro de inyección de cada especie por una función de la forma según trabajos previos6. En la Fig. 2 (izquierda) puede observarse la comparativa de dos histogramas de energía que corresponden a la simulación sobre la base y sobre el pico respectivamente. A su vez, en la Fig. 2 (derecha) se encuentra la misma gráfica comparativa pero sólo para el flujo de muones sobre cada sitio. El flujo total resulta un 14.3 mayor en el pico que en la base, pero en el flujo de muones de alta energía esta diferencia decrece a un 2, siendo igualmente mayor el flujo en el pico. Esto indica que la tasa de pérdida de energía de los muones en su recorrido a lo largo de la atmósfera resulta significativamente menor que la de la mayoría de las otras partículas que componen la lluvia8. En el rango de energías considerado para los secundarios, GeV, el flujo de muones sobre el pico resultaría de aproximadamente díamsr y sobre la base de díamsr. Además de los muones, las otras partículas que forman parte de la lluvia atmosférica extendida, por ejemplo las que conforman la componente hadrónica, también podrían resultar una fuente de señal de ruido. Sin embargo, este trabajo explora la componente muónica por ser estos los que tienen mayor probabilidad de llegar al laboratorio debido a que son las que tienen una menor tasa relativa de pérdida de energía8. En cuanto a otras partículas, por ejemplo, piones, podrían decaer en muones en el interior de la roca. De hecho, en la Fig. 2 de la izquierda que representa el flujo de todas las partículas, se puede ver la presencia de un flujo reducido de piones en energías cercanas a GeV pero este flujo resulta 4 órdenes menor al de muones, con lo cual en este estudio no se los considera.
Dirección de arribo de los muones
Conociendo la dirección de los muones con respecto al eje z, es decir, su ángulo cenital, se consideró como “dispersión” la diferencia entre el ángulo cenital del primario que lo generó y el del muón. Este estudio fue de importancia para el objetivo principal de este trabajo por dos motivos: en primer lugar, el ángulo de incidencia de cada muón determinaría el camino a recorrer por la partícula en el modelo de la montaña y así también su pérdida de energía. En segundo lugar, la dispersión de los muones con respecto a su partícula primaria da una idea de la capacidad de deflexión de los muones en su camino a través de la roca.
A partir de la distribución de lo que se denominó “dispersión” se calcularon el valor medio y la desviación estándar de la misma obteniendo y . Estos valores permiten estimar que los muones altamente energéticos no se desvían en grandes ángulos con respecto a su primario, es decir que en su transcurso a través de la atmósfera mantienen una dirección aproximadamente constante.
Energía mínima de penetración y porcentaje de pérdida de energía
La distancia mínima de cobertura que presentará el laboratorio son 1700 m. A partir de este hecho y en base a los valores tabulados de poder de frenado de muones en roca,8 se estimó la energía mínima necesaria que requiere una partícula para atravesar 1700 m de roca estándar g cm. Se encontró esta energía en un valor alrededor de los 3 TeV.
Considerando que el poder de frenado depende de la energía inicial de la partícula resultó necesario conocer esta energía mínima de penetración al laboratorio para configurar las simulaciones de EAS en CORSIKA, en particular para establecer un límite inferior para la energía de las partículas secundarias y así optimizar las simulaciones.
Además, se estudió la variación de pérdida de energía en función de la densidad del material para tres partículas de distinta energía inicial: 800 GeV, 1000 GeV y 3000 GeV. En la Tabla 1 pueden verse los resultados de este estudio. Lo que se observa es una mayor pérdida de energía cuanto mayor es la energía inicial de la partícula, lo cual concuerda con los esperados según los estudios realizados en trabajos previos 8 para energías de estos órdenes, donde la curva de la aproximación de desaceleración continua es creciente.
Transferencia en la aproximación de desaceleración continua
Teniendo ya un modelo aproximado de la estructura montañosa que el flujo sobre la ladera debería atravesar para llegar al laboratorio se procedió a programar en lenguaje Python una función de transferencia de flujo para los muones. Los pasos del algoritmo fueron los siguientes: se tomaba un muon de las simulaciones y según su ángulo cenital se seleccionaba el camino que recorre dentro de la montaña. Luego, el muon seleccionado avanza 1 cm o 1 m según la derivada de la función del poder de frenado en la roca y se calcula su pérdida de energía. Entonces, si la energía del muon llega a 0 se lo descarta y se toma el siguiente muon, si no, avanza nuevamente. Para cada muon se repite este proceso hasta que su energía llega a 0 o hasta que alcanza la profundidad del laboratorio. Un diagrama de flujo del algoritmo explicado puede encontrarse en la Fig. 3.
Opacidad
A partir de la definición (1) es posible estudiar la opacidad direccional de la cobertura rocosa para así determinar la función de transferencia de los muones a lo largo de la montaña. De esta forma el flujo resultante en el laboratorio será una combinación de la opacidad de la montaña y del flujo de partículas que llegan a la superficie. Se calculó el mapa de opacidad en función del ángulo cenital de la montaña para un observador ubicado en el laboratorio, tal como se muestra en la Fig. 4.
En este mapa puede observarse que a partir de los ángulos cercanos a la vertical la opacidad decrece hasta llegar a un mínimo en aproximadamente los 50. A partir de este ángulo comienza a aumentar nuevamente, encontrándose en un valor máximo entre los 70 y los 90. Además, se calculó el mapa de opacidad para una montaña con la misma geometría pero considerando que estuviera formada en su totalidad por roca estándar (g cm). Este estudio se realizó con el fin de caracterizar la influencia de las distintas densidades. El mapa resultante puede verse en la Fig. 5.
III. RESULTADOS Y DISCUSIÓN
Distribución Angular
A partir de la identificación de la energía y dirección de propagación de cada secundario transferido, se determinó la distribución angular del flujo de muones esperado en el laboratorio. En la Fig. 6 puede verse un histograma con los mismos intervalos angulares que se usaron para la función de transferencia. En esta figura puede observarse una mayor incidencia de partículas desde el pico cuyo flujo se atenúa a medida que aumenta el ángulo cenital.
Este resultado concuerda con lo esperado según el mapa de opacidad presentado en la sección anterior. Al tener en cuenta la geometría de la montaña, cuyo radio de base es mayor que el doble de su altura y que los estratos más densos se encuentran en la base de la montaña, es esperable que desde estas direcciones se observe una atenuación pronunciada del flujo de muones.
Flujo total de muones esperado en el laboratorio
Una vez realizada la transferencia en roca, se obtuvo la cantidad de partículas que ingresan al laboratorio etiquetadas cada una de ellas con su energía resultante al final de la transferencia, es decir, la energía con la que llegarían al laboratorio. Luego, a partir de estos datos se realizaron histogramas de la distribución de energías de estas partículas. En la Fig. 7 se encuentra el total del flujo, donde puede observarse que las partículas abarcan un rango de energías desde el orden de 1 GeV hasta el orden de los 1000 GeV encontrándose el valor medio para estas energías en 243 GeV.
Una vez que se normalizó el flujo total por los 90 días de simulación y por unidad de ángulo sólido, se obtuvo un flujo de muones por día .
Este resultado resulta relevante a la hora de evaluar la atenuación ocasionada por la montaña dentro del laboratorio, que se esperaba que fuera del orden de 10-71. En primer lugar, el corte de energía de los secundarios en 800 GeV y el hecho de que se trabajó únicamente con muones reduce la cantidad de partículas a un orden de 105, y en segundo lugar, los muones con posibilidades de alcanzar el laboratorio son del orden de 1 en 100 en total, este cálculo confirma las estimaciones iniciales.
IV. CONCLUSIONES
En este trabajo se realizaron estudios acerca del blindaje de astropartículas en el laboratorio subterráneo ANDES, los cuales implicaron: dos simulaciones de tres meses de flujo cada una utilizando los Softwares CORSIKA y ARTI para los sitios designados como el pico y la base de la montaña bajo la cual se encontrará el laboratorio. Y además un análisis de la variación de la pérdida de energía de tres muones de energías iniciales de 800 GeV, 1000 GeV y 3000 GeV. Se observó que para una misma energía el efecto de la densidad sobre la pérdida de energía está entre un 3% y un 4% cada 0.5 g cm-3.
Partiendo del modelo geológico de la montaña, se estimó la función de transferencia direccional del flujo de muones en su propagación por el interior de la montaña en la aproximación CSDA. De esta forma, se determinó que el flujo esperado de muones en el interior del laboratorio subterráneo ANDES será de (1.47±0.02) por día.