La cosmología moderna se apoya en el análisis de grandes cartografiados de galaxias para entender la estructura a gran escala del universo. Un aspecto crucial es el modelado del espectro de potencia de los trazadores sesgados (galaxias), que no reflejan perfectamente la distribución de la materia oscura. Durante décadas, se han desarrollado plantillas perturbativas en marcos eulerianos y lagrangianos para el modelo cosmológico estándar $\Lambda$CDM. Sin embargo, para ir más allá de $\Lambda$CDM y explorar teorías de gravedad modificada, se requieren herramientas más sofisticadas que puedan manejar los regímenes no lineales de manera precisa.

Este trabajo aborda la implementación de la expansión perturbativa sesgada, dentro del esquema de sesgo lagrangiano local, en el marco de la Teoría de Campo Efectiva Híbrida (HEFT). HEFT combina la teoría de perturbaciones con simulaciones de materia oscura para modelar el régimen no lineal. Los investigadores se centraron en la gravedad $f(R)$, una teoría de gravedad modificada que presenta un crecimiento dependiente de la escala y un efecto de "camuflaje" (chameleon screening), lo que la convierte en un escenario particularmente desafiante para los cálculos de la teoría de perturbaciones lagrangiana y para la generación de simulaciones numéricas precisas.

Los autores presentan una descripción detallada de los elementos necesarios para calcular analíticamente los espectros de potencia sesgados con correcciones de bucle. Estas predicciones analíticas se comparan con los resultados de simulaciones completamente no perturbativas, validando el enfoque. Finalmente, proponen una estrategia para extender emuladores existentes basados en HEFT para $\Lambda$CDM, como \texttt{bacco} y \texttt{Aemulus}, a cosmologías que van más allá del modelo estándar, abriendo la puerta a un análisis más robusto de datos de futuros cartografiados de galaxias en el contexto de teorías de gravedad modificada.