Parallel Version of CORSIKA Code with Cherenkov Option for SPHERE-3 Project

Este artigo descreve o desenvolvimento de uma versão multithread do código CORSIKA com opção de luz Cherenkov, executada no supercomputador Lomonosov-2, para otimizar a geração de eventos de chuveiros atmosféricos extensos e evitar que processos sejam interrompidos por limites de tempo na fila de processamento.

O essencial

Imagine que você está tentando prever o tempo, mas em vez de nuvens e chuva, você está tentando entender o que acontece quando uma pedra gigante do espaço (um raio cósmico) bate na atmosfera da Terra.

Esse é o trabalho dos cientistas do projeto SPHERE-3. Eles querem saber do que são feitas essas pedras cósmicas (se são de ferro, hidrogênio, etc.) para entender melhor o universo.

O problema é que, quando essas pedras batem no ar, elas criam uma "explosão" de outras partículas menores, como uma cascata de dominós caindo. Para estudar isso, eles usam um programa de computador chamado CORSIKA, que simula essa explosão.

Aqui está o resumo do que os autores fizeram, explicado de forma simples:

1. O Problema: A Corrida Contra o Relógio

Imagine que você precisa desenhar uma paisagem complexa, mas o computador que você usa tem uma regra estranha: se você demorar mais de 10 horas para desenhar, ele apaga tudo e você começa do zero.

Para simular as explosões de raios cósmicos mais energéticos (aquelas muito poderosas), o programa original demorava 20 horas para simular apenas uma única explosão. Como o computador "desligava" a simulação antes de terminar, os cientistas ficavam sem dados. Era como tentar encher uma banheira com um balde, mas a torneira fecha antes de encher.

2. A Solução: A Equipe de Resgate (Multithreading)

Os cientistas decidiram mudar a estratégia. Em vez de uma única pessoa (o processador do computador) fazer todo o trabalho, eles criaram uma equipe.

Eles modificaram o código do programa para que ele trabalhe em duas etapas:

• Etapa 1: O Chefe (Thread Principal)
  O "Chefe" entra na simulação e guia a partícula principal (a pedra cósmica) até que ela perca um pouco de sua energia. Ele faz isso sozinho, rápido e sem dividir a atenção. Ele para quando a partícula principal ainda tem cerca de 2% da sua energia original.
  Analogia: É como um maestro que organiza a orquestra e dá o primeiro batimento, mas não toca todos os instrumentos.

• Etapa 2: A Equipe (Threads Secundários)
  Quando o Chefe para, ele olha para a "bagunça" de partículas que foram criadas até aquele momento. Ele pega essa pilha de partículas e divide o trabalho entre 10 amigos (processadores) ao mesmo tempo.
  Analogia: Imagine que o Chefe tem uma pilha enorme de caixas para entregar. Em vez de ele carregar todas sozinho, ele divide a pilha em 10 montes menores e entrega um para cada um dos 10 ajudantes. Todos entregam suas caixas ao mesmo tempo.

3. O Desafio da Divisão Justa

Dividir o trabalho não é tão fácil quanto parece. Algumas partículas são "pesadas" (como um raio gama gigante) e levam muito tempo para serem processadas. Se o Chefe der uma dessas partículas "pesadas" para um ajudante e as outras partículas para os outros, o ajudante com a partícula pesada vai demorar muito, enquanto os outros ficam parados esperando.

Os cientistas criaram um algoritmo inteligente para tentar equilibrar a carga, garantindo que ninguém fique ocioso por muito tempo, embora ainda haja momentos em que alguns computadores ficam esperando.

4. Os Resultados: Mais Rápido e Preciso

O resultado dessa nova abordagem foi impressionante:

• Velocidade: O que antes levava 20 horas para ser feito, agora leva apenas 7,5 horas. O programa ficou quase 3 vezes mais rápido.
• Precisão: Eles verificaram se a "divisão de trabalho" estragou os dados. A resposta foi não. As simulações feitas pela equipe são idênticas às feitas pelo "Chefe" sozinho. É como se você tivesse 3 pessoas desenhando a mesma paisagem em 3 partes diferentes; quando juntam, a imagem final é perfeita e não tem erros.

Conclusão

Essa nova versão do programa é como transformar um carro de corrida de um único assento em um trem de alta velocidade com vários vagões. Agora, os cientistas podem gerar milhões de simulações de raios cósmicos para o projeto SPHERE-3, o que ajudará a responder a perguntas antigas sobre a composição da matéria no universo, sem que o computador "desista" no meio do caminho.

Em resumo: Eles pegaram um programa lento, dividiram o trabalho em equipe e conseguiram terminar a tarefa antes que o relógio acabasse.

Resumo técnico

Aqui está um resumo técnico detalhado do artigo em português:

Título: Versão Paralela do Código CORSIKA com Opção de Luz Cherenkov para o Projeto SPHERE-3

1. O Problema

O projeto SPHERE-3 visa estudar a composição de massa dos raios cósmicos primários (PCR) na faixa de energia de $10^{15}$a$10^{17}$ eV (1 PeV a 100 PeV). Para isso, é necessário simular vastas quantidades de chuveiros atmosféricos extensos (EAS) que incluam a geração de luz Cherenkov.

• Limitação Computacional: O código original CORSIKA-7 não suporta processamento paralelo de sub-cascadas com geração de luz Cherenkov, exigindo o uso de um único núcleo de CPU por evento.
• Tempo de Execução: Em energias acima de 70 PeV, a simulação de um único evento excede o limite de tempo da fila de processamento do supercomputador Lomonosov-2 (onde o projeto é executado).
• Consequência: Muitos eventos são "matados" (cancelados) antes de serem concluídos, inviabilizando a geração da base de dados necessária para a calibração e otimização do detector.

2. Metodologia

Os autores desenvolveram uma versão multithread (multithreading) do código CORSIKA modificada, utilizando uma estratégia de dois estágios para distribuir a carga computacional:

• Modificações de Saída: O código foi adaptado para gerar arquivos binários compactos contendo distribuições espaciais e temporais da luz Cherenkov em diferentes altitudes (nível do solo e 500m, 1000m, 1500m), eliminando a necessidade de arquivos de saída padrão de partículas secundárias.
• Algoritmo de Dois Estágios:
  1. Fase Sequencial (Mestre): O thread principal (master) realiza as inicializações padrão e rastreia a partícula primária e sua descendente mais energética (chamada de "líder") através da atmosfera. O rastreamento do líder continua até que sua energia caia abaixo de um limiar de 2% da energia primária ou até que apareça um fóton de alta energia na pilha (stack).
  2. Fase Paralela (Escravos): Uma vez atingido o limiar, o thread mestre distribui a pilha de partículas secundárias restantes entre threads escravos (workers).
• Algoritmo de Distribuição de Carga (Stack Distribution):
  • A pilha de partículas é ordenada por energia crescente.
  • O objetivo é equilibrar a energia total atribuída a cada thread escravo ($e_p = E_{total} / n_{slaves}$).
  • Partículas são agrupadas em sub-pilhas até que a energia acumulada exceda o alvo.
  • Desafio Identificado: A indivisibilidade das partículas e a presença de fótons de alta energia (que podem representar até 50% da energia total) podem causar desbalanceamento, deixando alguns threads ociosos. Uma melhoria proposta é ajustar o alvo de energia para $E_{total} / (n_{slaves} + 1)$.

3. Contribuições Chave

• Desenvolvimento de Código Paralelo: Criação da primeira versão multithread do CORSIKA capaz de gerar luz Cherenkov, permitindo o uso eficiente de múltiplos núcleos de CPU.
• Otimização de Memória e I/O: Redução significativa do tamanho dos arquivos de evento (de ~6 GB para <1 GB comprimido) e otimização do uso de memória (7 GB para o processo mestre + 450 MB por escravo).
• Validação Física: Implementação de testes rigorosos para garantir que a paralelização não introduzisse viéses sistemáticos nas distribuições físicas dos chuveiros.

4. Resultados

Os testes foram realizados em um servidor local (AMD Ryzen 9 5950X, 16 núcleos) e validados contra a versão sequencial original.

• Desempenho (Aceleração):
  • Para prótons de $10^{17}$ eV, o tempo de simulação caiu de 20 horas (sequencial) para 7,5 horas (paralelo).
  • O fator de aceleração ($S = T_{seq} / T_{par}$) variou entre 2,2 e 3,6, dependendo da energia e da composição nuclear (próton vs. ferro).
  • A versão paralela permitiu processar eventos que antes seriam cancelados por limite de tempo na fila do supercomputador.
• Validação Física:
  • Distribuição Lateral (LDF): As distribuições transversais de fótons Cherenkov entre as versões sequencial e paralela foram idênticas, com diferenças apenas no nível das flutuações estatísticas.
  • Rendimento de Fótons: A discrepância no número médio total de fótons foi de 1–4% para núcleos de ferro e 1–8% para prótons. Essa pequena variação é atribuída à variabilidade intrínseca dos chuveiros de prótons e a flutuações de amostragem, não a erros sistemáticos do algoritmo.

5. Significância e Conclusão

• Viabilidade do Projeto SPHERE-3: A nova versão do código torna viável a geração da vasta base de dados de eventos necessária para o projeto SPHERE-3, que utiliza um telescópio aéreo para observar a luz Cherenkov refletida na neve do Lago Baikal.
• Eficiência Computacional: A solução supera as limitações de tempo de fila do supercomputador Lomonosov-2, permitindo simulações de alta energia que antes eram impraticáveis.
• Limitações e Futuro: O principal gargalo restante é a distribuição não uniforme de partículas entre os threads devido à natureza estocástica das interações de partículas de alta energia. Trabalhos futuros focarão em otimizar o algoritmo de particionamento da pilha e explorar a aceleração via GPU para fragmentos específicos do código.

Em resumo, o trabalho apresenta uma solução técnica robusta que permite a simulação eficiente de chuveiros atmosféricos de ultra-alta energia com luz Cherenkov, garantindo a continuidade e o sucesso do projeto de pesquisa de raios cósmicos SPHERE-3.