Spectral Leakage and Masking Effects in the Measurement of Hyperuniformity

Resumen Técnico: Fuga Espectral y Efectos de Enmascaramiento en la Medición de la Hiperuniformidad

Planteamiento del Problema
Los sistemas hiperuniformes desordenados (DHU) se caracterizan por la supresión completa de las fluctuaciones de densidad en longitudes de onda infinitas ($\lim_{|k|\to 0} S(k) = 0$) mientras carecen de un orden de largo alcance convencional. La detección y clasificación de estos sistemas dependen de la medición precisa del comportamiento de escalamiento del factor de estructura en números de onda pequeños, $S(k) \sim k^\alpha$. Sin embargo, en la práctica, las mediciones suelen estar restringidas a subsistemas finitos o implican observaciones incompletas (por ejemplo, regiones enmascaradas) debido a limitaciones experimentales o restricciones de simulación. Estudios previos han señalado que el muestreo con ventanas finitas afecta los exponentes medidos, pero faltaba un marco teórico unificado que cuantificara cómo tanto las ventanas de observación compactas como las máscaras binarias espacialmente correlacionadas distorsionan el factor de estructura observado ($S_{obs}(k)$). Esta distorsión conlleva el riesgo de identificar erróneamente el orden hiperuniform genuino como artefactos o de oscurecer el verdadero exponente hiperuniform $\alpha$.

Metodología
Los autores desarrollan un marco teórico unificado basado en el teorema de la convolución en el espacio de Fourier. La premisa central es que cualquier procedimiento de observación que multiplique el campo de densidad $\rho(r)$ por una función de observación $W(r)$ (que representa una ventana finita) o una máscara binaria $M(r)$ en el espacio real, resulta en una convolución del factor de estructura intrínseco $S(k)$ con la densidad espectral de la función de observación en el espacio recíproco.

1. Efectos de Ventana Finita: El artículo deriva una fórmula de convolución exacta donde el factor de estructura observado $S_{obs}(k)$ es la convolución del $S(k)$ intrínseco con el espectro de potencia de la función de ventana $|\tilde{W}(k)|^2$. Los autores analizan dos regímenes asintóticos:
   • Régimen de Fuga Espectral ($k \lesssim 1/L$): Donde $L$ es el tamaño lineal de la ventana.
   • Régimen Intrínseco ($1/L \ll k \ll q_c$): Donde $q_c$ es la escala de cruce donde se mantiene la ley de potencia $S(k) \sim k^\alpha$.
2. Efectos de Enmascaramiento: El marco se extiende a máscaras binarias espacialmente correlacionadas. El espectro observado se expresa como una convolución de $S(k)$ con la densidad espectral de la máscara $\chi_M(k)$. Los autores analizan casos específicos: máscaras deterministas, adelgazamiento de Bernoulli independiente (no correlacionado) y máscaras con correlaciones a gran escala (máscaras de Debye).
3. Verificación Numérica: Las predicciones teóricas se validan mediante simulaciones numéricas de configuraciones de puntos hiperuniformes desordenados en $\mathbb{R}^2$ con diversos exponentes de hiperuniformidad ($\alpha = 0.5$ a $3.0$) y sistemas hiperuniformes sigilosos (stealthy). Estos sistemas son sometidos a ventanas cuadradas finitas y tanto a máscaras de Bernoulli como de Debye.

Contribuciones Clave y Resultados

• Fuga Cuadrática Universal: Para sistemas hiperuniformes genéricos con $S(k) \sim k^\alpha$, los autores demuestran que las ventanas de observación finitas inducen un término de fuga cuadrática universal en números de onda suficientemente pequeños ($k \lesssim 1/L$). Esto se manifiesta como un escalamiento aparente de $k^2$ ($S_{obs}(k) \approx A_0 + A_2 k^2$) que es independiente del exponente $\alpha$ real. En consecuencia, el exponente intrínseco $\alpha$ solo puede extraerse de manera fiable dentro del régimen intermedio $1/L \ll k \ll q_c$.
• Dependencia del Exponente $\alpha$: Se demuestra que la precisión en la estimación de $\alpha$ depende fuertemente de su valor real. Para sistemas con exponentes grandes ($\alpha \gtrsim 2$), la señal intrínseca decae rápidamente a números de onda pequeños, lo que la hace altamente susceptible a la contaminación por el piso de fuga y al ruido aditivo. Esto conduce a un sesgo sistemático a la baja en los exponentes estimados, a menos que se utilicen ventanas de observación muy grandes o un promedio de ensamble extensivo.
• Hiperuniformidad Sigilosa (Stealthy): En los sistemas hiperuniformes sigilosos, donde el $S(k)$ intrínseco posee un hueco espectral ($S(k)=0$ para $|k|<K$), toda la potencia observada a números de onda pequeños proviene enteramente del mecanismo de convolución. El escalamiento observado es, por tanto, un reflejo directo de la densidad espectral de la ventana y no de las fluctuaciones de densidad intrínsecas.
• Efectos de Enmascaramiento:
  • Máscaras No Correlacionadas (Bernoulli): Resultan en una atenuación del espectro padre más un término de ruido blanco aditivo. No alteran significativamente el exponente de escalamiento $\alpha$ medido dentro de la tolerancia numérica.
  • Máscaras Correlacionadas (Debye): Las máscaras con potencia sustancial a números de onda pequeños (largas longitudes de correlación) alteran drásticamente el comportamiento de $S_{obs}(k)$ a números de onda bajos. El exponente medido puede ser completamente diferente del sistema original, pudiendo destruir la firma de la hiperuniformidad o inducir un escalamiento artificial.
• Limitaciones del Problema Inverso: El artículo establece que recuperar el exponente intrínseco a partir de mediciones de ventana finita es un problema mal condicionado si el rango de escalamiento intrínseco es estrecho, si $\alpha$ es grande, o si los pisos de fuga/ruido son significativos. Para sistemas sigilosos, la señal a números de onda pequeños es enteramente inducida por la observación y no puede interpretarse como una ley de potencia intrínseca.

Significado
El artículo proporciona criterios cuantitativos para distinguir el orden hiperuniform genuino de los artefactos inducidos por la medición. Al establecer la relación de convolución entre los factores de estructura observados e intrínsecos, el marco ofrece directrices prácticas para:

1. Diseño Experimental y Numérico: Determinar el tamaño necesario de la ventana de observación $L$ para acceder al régimen intrínseco ($1/L \ll k$) y evitar el régimen dominado por la fuga.
2. Interpretación de Datos: Reconocer que el escalamiento aparente de $k^2$ en los números de onda más bajos accesibles es a menudo un artefacto del muestreo finito y no una propiedad física del sistema.
3. Estrategias de Corrección: Resaltar que la extracción fiable de $\alpha$ requiere un control cuidadoso de la función de observación o el uso de deconvolución/ajuste basado en modelos, particularmente cuando se trata de máscaras correlacionadas o sistemas sigilosos.

Los autores concluyen que, si bien las mediciones de ventana finita pueden caracterizar el comportamiento efectivo dependiente de la escala, no pueden por sí mismas establecer el límite estricto de $k \to 0$ sin considerar las propiedades espectrales específicas de la función de observación. Este trabajo contribuye a una caracterización más robusta y reproducible de los estados hiperuniformes en sistemas físicos, biológicos e ingenieriles.