Independent Recovery of Vanishing Sources on POSS-I Photographic Plates Using Automated Source Detection and Cross-Epoch Matching

1. Problema e Contexto

O projeto VASCO (Vanishing and Appearing Sources during a Century of Observations) catalogou milhares de fontes transitórias candidatas comparando placas fotográficas históricas (como a POSS-I, de 1949–1957) com arquivos digitais modernos. Resultados notáveis incluem nove fontes pontuais que apareceram e desapareceram simultaneamente em uma única placa de abril de 1950 e um catálogo público de 5.399 candidatos a fontes desaparecidas.

A natureza desses transientes é debatida: alguns argumentam que são falhas de emulsão, enquanto outros sugerem flashes ópticos subsegundos. Mais controversamente, um estudo anterior (Bruehl & Villarroel, 2025) relatou uma associação estatística entre a ocorrência desses transientes e testes nucleares atmosféricos realizados durante o período da POSS-I.

O objetivo deste trabalho é realizar uma replicação independente desses resultados utilizando um pipeline totalmente automatizado e diferente, sem referência aos catálogos existentes de transientes, para validar a robustez das detecções e testar a associação temporal com testes nucleares.

2. Metodologia

O autor desenvolveu um pipeline automatizado que processa recortes digitais da POSS-I. As etapas principais são:

• Dados: Utilização de recortes do Digitized Sky Survey (DSS) do STScI, cobrindo três épocas: POSS-I Vermelho (1949–1957), POSS-I Azul (mesma época) e POSS-II Vermelho (1985–2000).
• Detecção e Filtragem:
  • Detecção de fontes via limiarização com sigma-clipping acima do fundo local.
  • Filtragem morfológica baseada no PSF (Função de Espalhamento de Ponto): são mantidas apenas fontes dentro de uma faixa específica de FWHM (largura total a meia altura) e com baixa elipticidade, rejeitando objetos estendidos e ruído.
  • O pipeline foi calibrado em um conjunto de referência de 20 casos, resultando em parâmetros de 8σ, FWHM de 0,7–1,5× a mediana e elipticidade < 0,3.
• Refinamento Astrométrico: Correção local afim para cada recorte, alinhando as 20–50 fontes mais brilhantes entre as épocas, reduzindo resíduos para 0,2–0,5 arcsec.
• Identificação de Transientes: Uma fonte na POSS-I Vermelho é marcada como candidata se não houver correspondência na POSS-I Azul e na POSS-II Red dentro de um raio de 8 arcsec.
• Verificação Pan-STARRS: Candidatos são cruzados com o Pan-STARRS DR1. Fontes com correspondência moderna dentro de 3 arcsec são rejeitadas como objetos persistentes que falharam no cruzamento devido a diferenças astrométricas ou morfológicas.
• Análise Temporal:
  • Binária: Verifica se houve observação de pelo menos um transiente em dias específicos relativos a testes nucleares.
  • Modelo GLM Binomial Negativo: Modela a contagem de candidatos por noite, usando a cobertura de patches noturna como termo de exposição, para controlar a variação na área do céu observada.

3. Principais Contribuições e Resultados

Validação e Desempenho

• Recuperação Independente: O pipeline recuperou 8/9 fontes conhecidas no campo de abril de 1950 e 3/3 no campo de julho de 1952.
• Taxa de Falsos Positivos: Em campos de controle aleatórios não-crowded, a taxa foi de ~0,2 por campo de 10 arcmin. Em campos da Via Láctea densos, a taxa aumentou devido à confusão de fontes (até 22 por campo).
• Comparação com Solano et al. (2022): Ao cruzar o catálogo gerado (2,85 milhões de candidatos) com o catálogo público de 5.399 fontes de Solano et al., o pipeline recuperou independentemente 3.450 fontes (63,9%), com uma separação angular mediana de 0,94 arcsec.
  • Das fontes não recuperadas, ~54% foram detectadas mas rejeitadas pelo filtro de PSF (indicando discordância morfológica), e ~46% não foram detectadas acima do limiar de 5σ.
  • Nenhuma das fontes não recuperadas (dentro da cobertura) tinha uma contraparte no Pan-STARRS DR1, validando a classificação de "ausência de arquivo".

Levantamento de Campo Cheio (Full-Sky)

• O pipeline varreu 122.991 posições de grade (30 arcmin) sobre a cobertura POSS-I, gerando um catálogo filtrado de 2,85 milhões de candidatos a fontes desaparecidas.
• O custo computacional foi de aproximadamente $50 USD.

Análise Temporal com Testes Nucleares

• Análise Binária (Dia de Calendário): Observou-se uma elevação na taxa de detecção no dia imediatamente seguinte aos testes nucleares ([+1]), com um Risco Relativo (RR) de 1,35 (IC 95% [0,91, 2,00]). No entanto, o resultado não foi estatisticamente significativo (p = 0,17, teste exato de Fisher bilateral). O autor nota que este resultado é sensível ao fato de que todas as 368 noites de observação no intervalo de estudo continham pelo menos um candidato, tornando o indicador binário dependente do cronograma de observação.
• Modelo Binomial Negativo: Ao modelar a contagem de candidatos por noite (controlando a cobertura do céu), o efeito desapareceu. O Incidence Rate Ratio (IRR) para a janela ±1 dia foi de 1,03 (IC 95% [0,89, 1,18]; p = 0,71). Nenhuma das janelas temporais testadas mostrou significância estatística.

4. Significado e Conclusões

• Validação Robusta: O estudo confirma que uma grande fração (63,9%) das fontes desaparecidas catalogadas anteriormente pode ser recuperada independentemente com alta precisão posicional, sugerindo que não são meros artefatos de processamento específicos de um único pipeline.
• Inconclusividade Temporal: Ao contrário da replicação do catálogo (que foi forte), a associação temporal com testes nucleares permanece inconclusiva com os dados atuais. O modelo estatístico mais rigoroso (contagem de candidatos por noite) não encontrou evidências de um efeito temporal.
• Limitações e Futuro: O catálogo atual de 2,85 milhões de fontes foi gerado com parâmetros preliminares (5σ), resultando em uma taxa de falsos positivos mais alta do que o ideal. O autor conclui que qualquer afirmação futura sobre a associação com testes nucleares exigirá:
  1. Um reprocessamento completo com os parâmetros calibrados (8σ) para reduzir a linha de base de candidatos.
  2. Controle explícito do cronograma de observação e agrupamento de placas.
  3. Janelas temporais pré-especificadas para evitar interpretações exploratórias excessivas.

Em resumo, o trabalho valida a existência de um catálogo robusto de fontes desaparecidas na POSS-I, mas não fornece evidência independente que suporte a hipótese de que essas fontes são causadas por testes nucleares atmosféricos, destacando a necessidade de filtragem mais rigorosa para tais análises temporais.