Projection effects in star-forming regions: I. Nearest-neighbour statistics and observational biases

Resumo Técnico: Efeitos de Projeção em Regiões de Formação Estelar I. Estatística de Vizinhos Mais Próximos e Vieses Observacionais

Enunciado do Problema
As estrelas se formam através da fragmentação hierárquica de nuvens moleculares em núcleos densos, filamentos e subaglomerados. Uma medida observável primária usada para restringir a física dessa fragmentação é o espaçamento característico entre núcleos densos, frequentemente comparado a escalas teóricas como o comprimento de Jeans ou o comprimento sônico. No entanto, as medições observacionais são, inerentemente, projeções bidimensionais de estruturas tridimensionais. Embora um fator de desprojeção geométrica simples de $4/\pi \approx 1,27$ seja frequentemente aplicado sob a suposição de isotropia, este artigo argumenta que tal fator é inadequado. A abordagem padrão falha ao não considerar dois vieses concorrentes: (1) a projeção não apenas encurta os separações, mas também "reconfigura" a rede de vizinhos mais próximos (NN - nearest-neighbour), criando conexões 2D artificiais entre fontes que estão amplamente separadas em 3D; e (2) a resolução angular finita introduz a fusão de feixes (beam blending), que funde vizinhos próximos e infla os separações aparentes. Sem um arcabouço de correção robusto, as inferências sobre escalas de fragmentação e os processos físicos dominantes (térmico vs. turbulento) podem ser sistematicamente enviesadas.

Metodologia
Os autores utilizam uma hierarquia de modelos numéricos de brinquedo controlados para quantificar os vieses de projeção através de um amplo espaço de parâmetros.

• Modelos: Eles geram distribuições 3D de "núcleos" pontuais usando duas geometrias principais: (i) perfis esféricos uniformes e centralmente concentrados, e (ii) distribuições fractais hierárquicas (prescrição de Cartwright & Whitworth) com dimensões fractais variando ($D \approx 1,6–2,5$) para simular subestruturas clumpadas e filamentares.
• Processo de Simulação: Para cada configuração, a distribuição 3D é rotacionada e projetada em um plano 2D. Os autores então computam o grafo NN tanto em 3D quanto em 2D.
• Efeitos de Resolução: Para mimetizar limitações observacionais, eles aplicam uma etapa de "fusão de feixe" onde pontos dentro de um raio de feixe são fundidos em um único centroide. Isso é parametrizado pelo Alcance Dinâmico Espacial (SDR), definido como a razão entre o Campo de Visão (FoV) e a FWHM (largura total à meia altura) do feixe.
• Análise: O estudo varia sistematicamente o tamanho da amostra ($N$), a dimensão fractal e o SDR. Eles comparam comprimentos de aresta NN 3D intrínsecos ($\ell_{3D}$) contra comprimentos NN 2D projetados ($\ell_{2D}$), analisando tanto a razão média $C = \langle \ell_{3D} \rangle / \langle \ell_{2D} \rangle$ quanto a sobreposição topológica (similaridade de Jaccard) entre as redes 3D e 2D.

Resultos Principais

1. Falha do Fator $4/\pi$: O fator geométrico padrão de $1,27$ subestima significativamente o espaçamento 3D real. Para casos bem amostrados e bem resolvidos, a razão $C$ é substancialmente maior.
2. Reconfiguração de Rede: A projeção não apenas encurta as conexões existentes; ela reorganiza fundamentalmente a rede NN. Em distribuições uniformes com $N=200$, apenas $\sim20\%$ das conexões NN 3D sobrevivem à projeção. As conexões restantes são substituídas por novos vizinhos 2D muito mais curtos, elevando a razão média $C$ para $\sim2,5$ em casos idealizados de alta resolução.
3. Dependência do Tamanho da Amostra: O viés escala com o número de objetos. Para amostras pequenas ($N \lesssim 10$), a razão é menor ($\sim1,5$), mas aumenta sistematicamente com $N$, seguindo uma lei de potência consistente com a diferença de escala dimensional entre os separações médias 3D ($N^{-1/3}$) e 2D ($N^{-1/2}$).
4. Dependência de Resolução (Fusão de Feixe): A resolução finita atua como um viés concorrente. Quando o SDR é baixo ($\lesssim 10$), a fusão de feixe funde núcleos distintos, removendo os separações intrínsecas mais curtas. Isso suprime o encurtamento induzido pela projeção, levando $C$ de volta para a unidade ($C \approx 1$).
5. Fórmula de Correção Empírica: Os autores derivam um fator de correção empírico simples $C(N, \text{SDR})$ que captura a dependência conjunta de tamanho de amostra e resolução:
   $$C(N, \text{SDR}) = C_\infty \left[1 - \exp(-\text{SDR}/S_0)\right] \left(\frac{N}{100}\right)^\beta$$
   Parâmetros ajustados são $C_\infty = 1,94$, $S_0 = 21,8$ e $\beta = 0,173$.
   • Regime I (Baixo $N$ ou Baixo SDR): $C \approx 1$. A fusão de feixe domina, apagando diferenças topológicas.
   • Regime II (Intermediário): $C$ sobe rapidamente conforme a reconfiguração de projeção domina.
   • Regime III (Alto $N$, Alto SDR): $C$ satura próximo a $\sim2,0$. Aqui, os separações 3D reais são tipicamente o dobro dos espaçamentos 2D observados.

Significância e Aplicação
O artigo fornece a primeira calibração sistemática do viés de projeção em separações de núcleos, indo além da suposição simplista de $4/\pi$. Os autores demonstram que aplicar um fator multiplicativo dependente de resolução e tamanho de amostra permite que observadores estimem os separações 3D intrínsecos a partir de dados 2D.

• Impacto na Interpretação: Aplicar este arcabouço aos dados do levantamento ASHES revela que, embora os separações 2D não corrigidos pareçam consistentes com o comprimento de Jeans térmico, aplicar uma correção de tamanho de amostra ($C(N)$) sugere separações 3D maiores. No entanto, incorporar a correção de resolução ($C(N, \text{SDR})$) reduz novamente o espaçamento 3D inferido, destacando que o acordo aparente com escalas teóricas em mapas marginalmente resolvidos pode ser um artefato da fusão de feixe, em vez de uma realidade física.
• Limitações e Ressalvas: Os autores declaram explicitamente que este trabalho é um "primeiro passo" usando modelos de brinquedo idealizados. Eles deliberadamente omitiram complexidades como limites de sensibilidade, confusão de fundo, campos de visão incompletos e morfologias não uniformes. Consequentemente, as correções derivadas devem ser tratadas como aproximadas, com incertezas sistemáticas de $\sim30–40\%$ devido à morfologia e orientação. O fator de correção baseado em $N_{obs}$ observado é um limite inferior conservador para a correção real.
• Ferramentas: Os autores disponibilizaram o pacote corespacing3d para facilitar a aplicação destas correções.

Em conclusão, este trabalho estabelece que os efeitos de projeção em regiões de formação estelar são uma convolução de encurtamento geométrico, reconfiguração de rede e fusão induzida pela resolução. Um fator global único é insuficiente; a inferência precisa da física de fragmentação requer correções que levem em conta tanto o número de núcleos detectados quanto o alcance dinâmico espacial da observação.