# La Frontera a Alto Corrimiento al Rojo: Historia de la Formación de las Galaxias

## 1 Introducción

Ciertas cuestiones que han resultado en grandes avances en la historia de la astronomía y la cosmología tienen, además, un interés especial para la humanidad porque tratan sobre fenómenos totalmente imprescindibles para nuestra existencia en el universo: la fuente de energía del Sol, la formación de sistemas planetarios, y la síntesis de los elementos en el interior de las estrellas son ejemplos que acuden a la mente. La cosmología moderna, después de habernos encaminado por vez primera al estudio científico del origen del universo observable en su totalidad, en el marco de la teoría de la Gran Explosión, trae consigo también otra cuestión dentro de esa categoría: el origen de las fluctuaciones primordiales y la formación de las galaxias. El universo debe empezar con condiciones iniciales prácticamente homogéneas para resultar en el estado presente de homogeneidad a gran escala; pero la existencia de fluctuaciones de densidad primordiales, de una amplitud $`\delta \rho /\rho 10^5`$, es absolutamente necesaria para la formación de las galaxias. La gravedad puede amplificar y llevar las fluctuaciones iniciales al colapso no lineal, pero no puede crearlas. El proceso que generó esas fluctuaciones en el universo primitivo, sumido todavía en el misterio, dió lugar a la gran diversidad del universo no lineal, incluida la existencia de la vida, y evitó la continuación indefinida del universo lineal y homogéneo, conteniendo únicamente un mar de radiación y átomos de hidrógeno y helio.

Estamos actualmente en la era de exploración y descubrimiento hacia la frontera de alto corrimiento al rojo o, en otras palabras, de las mayores distancias des de las que es posible recibir mensajes en el universo. Dada la velocidad de la luz, las grandes distancias proporcionan a la cosmología la evidencia de la historia pasada de la formación de las galaxias. Durante las últimas tres décadas, nuestra visión del pasado del universo ha progresado continuamente con el descubrimiento de galaxias y núcleos activos a distancias cada vez mayores, el estudio de fondos de radiación cósmica a distintas frecuencias, y el análisis de espectros de absorción del hidrógeno intergaláctico interpuesto en la dirección de fuentes luminosas. Presentamos en este artículo de revisión un breve resumen del estado actual de la teoría y observaciones sobre la formación de las galaxias y la evolución del medio intergaláctico. No podemos hacer justicia en ese breve artículo al inmenso número de trabajos publicados en diversos temas de gran impacto sobre el estudio de la formación de las galaxias. Entre otros artículos de revisión, Ellis (1997) describe el estudio de galaxias débiles, y Rauch (1998) presenta una excelente exposición de observaciones y teorías del bosque Ly$`\alpha `$ de hidrógeno intergaláctico. Cabe recomendar asimismo varios artículos en libros de conferencias recientes de gran utilidad para ponerse al día en ese campo: véase, por ejemplo, Madau 1999, Steidel 1998a,b.

## 2 Teoría de la Materia Invisible Fría

De entre los distintos modelos propuestos de formación de galaxias, la teoría de la Materia Invisible Fría (que abreviaremos MIF; en inglés, Cold Dark Matter) ha resultado ser la de mayor éxito, y claramente favorecida por las observaciones. La teoría postula que la materia invisible, cuya presencia se deduce de las determinaciones dinámicas de la masa de galaxias y cúmulos (p.e., Trimble 1987), consiste en objetos o partículas “frías”, sin dispersión de velocidades inicial, y que las fluctuaciones primordiales de densidad son adiabáticas (es decir, manteniéndose constante la razón de la densidad de fotones, bariones y materia invisible), Gausianas, y con un espectro de potencias invariante de escala (véase por ejemplo, Blumenthal et al. 1984, Ostriker 1993). Generalmente, se supone que las fluctuaciones fueron generadas en un período de inflación por algún proceso que se mantuvo constante mientras el rango de escalas observable en el presente cruzaba el horizonte de acontecimientos, lo cual implica la invariancia de escala.

Inicialmente, la teoría MIF se consideró de forma casi exclusiva dentro del modelo cosmológico con densidad de materia crítica, $`\mathrm{\Omega }\rho /\rho _{crit}=1`$ (donde $`\rho `$ es la densidad media de materia y $`\rho _{crit}`$ es la densidad crítica). Desde los inicios de la cosmología observacional, las medidas de la densidad media del universo indicaron un valor $`\mathrm{\Omega }<1`$ (véase, p.e. Gott et al. 1974). Este resultado observacional se ha mantenido hasta la actualidad: la densidad media de luz en la banda B es en el presente $`\rho _B2\times 10^8hL_{}\mathrm{Mpc}^3`$ (Zucca et al. 1997 y referencias incluídas; utilizamos aquí$`H_0=100h\mathrm{km}\mathrm{s}^1\mathrm{Mpc}^1`$), y la razón masa-luminosidad en diversos sistemas colapsados, desde galaxias individuales a los mayores cúmulos, medida a radios suficientemente grandes, tiende generalmente a valores $`150hM/L_B500h`$ (p.e., Bahcall, Lubin, & Dorman 1995), resultando en una densidad de materia $`\overline{\rho }5\times 10^{10}h^2M_{}\mathrm{Mpc}^3`$, y $`\overline{\rho }/\rho _{crit}=\mathrm{\Omega }0.2`$.

El valor de $`\mathrm{\Omega }`$ así deducido a partir de la densidad luminosa del universo puede verse afectado por el efecto del sesgo de las galaxias. En principio, la mayor parte de la masa del universo podría ubicarse en los grandes vacíos, sin pertenecer a ninguna estructura virializada conteniendo galaxias donde la masa total pueda determinarse a través de los métodos dinámicos habituales. Esa fue la posibilidad bajo la cual el modelo Einstein-de Sitter ($`\mathrm{\Omega }=1`$) fue forzado a resguardarse de la evidencia observacional. Sin embargo, varios avances durante la última década han confirmado un valor de $`\mathrm{\Omega }0.3`$. En los cúmulos de galaxias de mayor masa, la fracción de la masa constituída por bariones deducida por la intensidad y el espectro de rayos X emitidos por el gas caliente es $`f_{bar}0.06h^{3/2}`$ (S. D. M. White et al. 1993; D. A. White & Fabian 1994). Dada la densidad de bariones deducida de la teoría de nucleosíntesis primordial, $`\mathrm{\Omega }_b0.019h^2`$ (p.e., Burles & Tytler 1998), y el hecho de que la fracción de masa bariónica en los cúmulos más masivos debe ser representativa de la fracción media en el universo (White et al. 1993), deducimos un valor $`\mathrm{\Omega }0.3h^{1/2}`$. Al mismo tiempo, las observaciones recientes de curvas de luz de supernovas tipo Ia deducen el mismo valor de $`\mathrm{\Omega }`$, concluyendo además que la expansión del universo está acelerando de la forma esperada en el modelo del universo con geometría espacial plana, donde la constante cosmológica proporciona la densidad de energía adicional necesaria para llegar a la densidad crítica (Perlmutter et al. 1998, Riess et al. 1998). Otros métodos de medir el valor de $`\mathrm{\Omega }`$ son, por lo general, consistentes con este resultado (por ejemplo, el valor $`\mathrm{\Omega }0.3`$ con una constante cosmológica $`\mathrm{\Lambda }=1\mathrm{\Omega }`$ también es favorecido por el valor de la constante de Hubble y la edad de las estrellas más antiguas en cúmulos globulares). Varias observaciones en el futuro próximo (tales como supernovas Tipo Ia y fluctuaciones en la radiación cósmica de fondo) deberán permitir medir $`\mathrm{\Omega }`$ con mayor precisión, y clarificar si la reciente aceleración del universo se debe a una constante cosmológica o al modelo más general de un campo escalar con presión negativa (p.e., Peebles & Vilenkin 1998 y referencias incluídas).

Una vez adoptamos el modelo cosmológico deducido a partir de esas observaciones, las predicciones de la tería MIF se ajustan bien a los datos observacionales sobre estructura a gran escala, tales como la función de correlación espacial de las galaxias, la abundancia de cúmulos de galaxias, la evolución en CR de esas cantidades, y las fluctuaciones en la radiación cósmica de fondo recientemente detectadas. Para obtener esas predicciones, es preciso calcular el espectro de potencias de las fluctuaciones primordiales.

El espectro de potencias de las fluctuaciones de densidad se obtiene en la teoría MIF calculando la evolución en el régimen lineal de esas fluctuaciones una vez cruzan el horizonte (para artículos de revisión, véase Efstathiou 1990, Bond 1996, Bertschinger 1996). En el límite de pequeñas escalas, las fluctuaciones entraron en el horizonte durante la época en que la densidad media del universo era dominada por radiación. Los bariones están entonces ligados a la radiación, y la presión radiativa resulta en una oscilación de las fluctuaciones, impidiendo el crecimiento gravitatorio. La materia invisible no interactúa con la radiación, y por lo tanto aumenta sus fluctuaciones de densidad. No obstante, ese crecimiento gravitatorio es muy lento incluso en la ausencia de acoplamiento con la radiación, cuando la radiación domina la densidad de energía.

Dado que la densidad de radiación en el presente es $`\mathrm{\Omega }_{rad}=aT_{rad}^4/(\rho _{crit}c^2)=2.47\times 10^5h^2`$ (donde $`T_{rad}=2.76`$ K es la temperatura de la radiación de fondo), la época de igualación de las densidades de materia y radiación corresponde al CR $`1+z=4.04\times 10^4\mathrm{\Omega }h^2`$. La longitud del horizonte en esta época era del orden $`cH^1(z)=cH_0^1\mathrm{\Omega }^{1/2}(1+z)^{3/2}`$, y por lo tanto la longitud comóvil era $`L_{ig}=cH_0^1\mathrm{\Omega }^{1/2}(1+z_{ig})^{1/2}15(\mathrm{\Omega }h)^1h^1\mathrm{Mpc}`$ (donde utilizamos la notación habitual para la constante de Hubble en el presente, $`H_0=100h\mathrm{km}\mathrm{s}^1\mathrm{Mpc}^1`$, y $`H(z)`$ es la constante de Hubble a CR $`z`$). A escalas comóviles mucho menores que el horizonte en la época de igualación, la amplitud de fluctuaciones de densidad debe conservar la forma de las fluctuaciones primordiales cuando emergen del horizonte, ya que las fluctuaciones empezaron a crecer solamente a partir de la época de igualación, y el crecimiento en la época anterior fue muy lento. Eso implica que la amplitud de fluctuaciones tiende a una constante para el caso más habitual de invariancia de escala, con $`n=1`$ (donde el espectro primordial es $`P(k)k^n`$). Pero en el límite de escalas mucho mayores, $`LL_{ig}`$, las fluctuaciones crecen en el régimen lineal proporcionalmente al factor de escala, $`a(1+z)^1`$, desde el momento en que entran en el horizonte; dado que la longitud comóvil del horizonte es proporcional a $`(1+z)^{1/2}`$, el factor total de crecimiento a una época fija debe ser proporcional a $`L^2`$.

Mostramos en la Figura 1 la dispersión en la fluctuación relativa de la masa total contenida dentro de una esfera de radio comóvil $`R`$, cuando la evolución lineal de fluctuaciones primordiales se extrapola hasta el presente. Conviene recordar aquí que la distribución de probabilidad de esta masa es Gausiana (dada la suposición de que el campo de densidades es Gausiano), con la dispersión mostrada en la Figura 1. La línea sólida es para el modelo $`\mathrm{\Omega }=0.3`$, $`\mathrm{\Lambda }=0.7`$, $`n=1`$, $`H_0=65\mathrm{km}\mathrm{s}^1\mathrm{Mpc}^1`$, y hemos normalizado las fluctuaciones a $`\sigma _8<\sigma ^2>^{1/2}(R=8h^1\mathrm{Mpc})=0.9`$. La línea de rayas es para el modelo $`\mathrm{\Omega }=1`$, $`\mathrm{\Lambda }=0`$, $`H_0=50\mathrm{km}\mathrm{s}^1\mathrm{Mpc}^1`$, y $`\sigma _8=0.55`$. Las curvas han sido calculadas con las fórmulas presentadas por Hu & Sugiyama (1996). Además del comportamiento esperado en el límite de pequeñas y grandes escalas, vemos también que el modelo de la constante cosmológica dispone de mayor potencia a grandes escalas, debido al mayor valor de la longitud característica $`L_{ig}(\mathrm{\Omega }h)^1`$.

La normalización del espectro de fluctuaciones, dada por el parámetro $`\sigma _8`$, debe considerarse como un parámetro ajustable de la teoría. Su valor ha sido escogido aquí para reproducir aproximadamente la abundancia de cúmulos de galaxias en el presente (Eke, Cole, & Frenk 1996). Desde que las fluctuaciones de intensidad de la radiación de fondo de microondas fueron detectadas por el satélite COBE, la normalización viene fijada independientemente por esas fluctuaciones (Bennett et al. 1996; Bunn & White 1997). Esas dos medidas de la amplitud concuerdan para el modelo con constante cosmológica con $`\mathrm{\Omega }0.3`$, mientras que son claramente conflictivas si se supone $`\mathrm{\Omega }=1`$. Es notable el hecho de que cuando se toma el valor de $`\mathrm{\Omega }`$ deducido de las observaciones de la geometría global del universo y la densidad promedio de materia, la teoría MIF predice correctamente la relación entre las fluctuaciones de la radiación de fondo y la abundancia de cúmulos de galaxias, y concuerda además con la correlación espacial de galaxias (Jing, Mo, & Börner 1998 y referencias incluídas), así como la evolución de la abundancia de cúmulos en CR (Bahcall & Fan 1998).

## 3 Colapso gravitatorio no-lineal: formación de halos

En la teoría MIF, la materia invisible empieza a formar objetos virializados (o “halos”) por colapso gravitatorio a partir de las escalas más pequeñas, cuando la fluctuación de densidad alcanza un valor $`1`$, entrando en el régimen no lineal. Luego, cuando las fluctuaciones a escalas mayores llegan al colapso, halos de mayor masa se forman por mersión de halos de masa menor colapsados anteriormente. Sólo podemos presentar aquí un resumen muy breve de la teoría de formación de halos y galaxias; para una discusión mucho más extensa, véase White (1996).

Un modelo analítico de gran utilidad para entender el proceso del colapso no lineal, llevando a la formación y mersión de halos, es el modelo de Press-Schechter (Press & Schechter 1974). El modelo se basa en la solución analítica del colapso de una perturbación de densidad constante con simetría esférica, en cuyo caso el colapso al punto central sucede cuando se alcanza el valor de la sobredensidad extrapolada linealmente $`\delta =\delta _c=3/5(3\pi /2)^{2/3}=1.686`$ (Peebles 1980). La distribución de la densidad promediada sobre una esfera de radio $`R`$ es una Gausiana con dispersión $`\sigma (R)`$ (Figura 1), y el modelo Press-Schechter consiste en suponer que la fracción de la masa del universo que ha colapsado en halos de masa mayor que $`M=(4\pi /3)\overline{\rho }R^3`$ es igual a $`\mathrm{erf}\{\delta _c/[2^{1/2}\sigma (R)]\}`$. Nótese que esa fracción es en realidad dos veces la fracción del volumen con sobredensidad superior a $`\delta _c`$ en las condiciones iniciales. Aunque la introducción del factor 2 puede justificarse mejor con un tratamiento en el espacio de Fourier (Bond et al. 1991), el modelo de Press-Schechter debe considerarse únicamente como una aproximación que se utiliza frecuentemente debido a su gran simplicidad, y especialmente al hecho de que sus predicciones se ajustan a resultados numéricos de la abundancia de halos con bastante precisión (p.e., Governato et al. 1998; a la práctica, el parámetro $`\delta _c`$ puede considerarse ajustable a los resultados numéricos, en vez de tomar el valor requerido para el colapso esférico).

Una vez obtenida la abundancia de halos de masa mayor que $`M`$, se obtiene por diferenciación la abundancia diferencial, y asignamos a cada escala $`R`$ una dispersión de velocidades para el halo obtenida también por el modelo esférico: $`\sigma ^2=(3\pi /2)^{2/3}[H(z_c)R]^2/2`$, donde $`H(z_c)`$ es la constante de Hubble en el momento del colapso. Básicamente, esa relación nos dice que la dispersión de velocidad de un cúmulo es del orden de la velocidad de expansión de Hubble sobre la región comóvil desde la cual el cúmulo ha colapsado, en el momento del colapso. La dispersión de velocidad de un halo (o, análogamente, el cociente de la masa sobre el radio) es la propiedad de un halo que determina las cantidades observables en galaxias y cúmulos, que son la dispersión de velocidad de las galaxias, la temperatura del gas difuso en equilibrio dinámico en el halo (determinada con espectros de rayos X), o la deflección gravitatoria de luz de fondo (medida en lentes gravitatorias).

En las Figuras 2(a,b), las tres líneas sólidas gruesas nos dan la dispersión de velocidades de un halo formado a partir de una fluctuación de densidad inicial igual a $`(1,2,3)\sigma `$, en función del CR, para los mismos dos modelos utilizados en la Figura 1 (para el modelo con constante cosmológica en la Fig. 2b, los valores de $`\delta _c`$ y de la dispersión de velocidades correspondiente a cada escala $`R`$ deben ser modificados debido al efecto de repulsión de la constante cosmológica, y el ritmo de crecimiento lineal de fluctuaciones es también distinto; véase p.e., Viana & Liddle 1996). De acuerdo con la distribución Gausiana, la fracción de masa que ha colapsado en halos de mayor dispersión de velocidades es, respectivamente, (64%,10%,0.6%). Tomemos, por ejemplo, el momento presente ($`z=0`$) para el modelo $`\mathrm{\Omega }=0.3`$ (Figura 2b). Ese modelo nos predice que un $`0.6\%`$ de la masa en el presente está en objetos colapsados con dispersión de velocidad mayor que $`1000\mathrm{km}\mathrm{s}^1`$. Esos halos se corresponden evidentemente con los cúmulos de galaxias más masivos, y se formaron a partir del colapso de raras fluctuaciones a escalas grandes, $`R15h^1\mathrm{Mpc}`$, donde $`3\sigma (R)=\delta _c`$. Las líneas de rayas nos dan el valor de la masa de los halos (con un factor 10 en masa separando líneas sucesivas), e indican que la masa de uno de esos cúmulos es $`2\times 10^{15}M_{}`$. La mayor parte de la masa del universo debe encontrarse en objetos formados a partir de fluctuaciones más habituales: una fluctuación $`1\sigma `$ produce un halo con dispersión de velocidades $`200\mathrm{km}\mathrm{s}^1`$, un valor típico de los pequeños grupos de galaxias, tales como el grupo Local, donde residen la mayoría de las galaxias en el universo. La masa indicada en la Figura 2 se refiere siempre a la masa total del halo en el momento de su formación. Por ejemplo, la galaxia de la Vía Láctea pudo formarse a $`z2`$, en un halo con $`\sigma 150\mathrm{km}\mathrm{s}^1`$ de masa $`M10^{12}M_{}`$, y en el presente el Grupo Local está colapsando en un halo mayor. Los dos modelos en la Figura 2 han sido normalizados para ajustar sus predicciones a la abundancia observada de cúmulos.

La dispersión de velocidades puede relacionarse con la temperatura del gas difuso en el halo cuando está en equilibrio hidorstático: $`kT/\mu =\sigma ^2`$. Mostramos también la temperatura en el eje vertical. Hemos utilizado aquí la masa media para el caso de gas primordial ionizado: $`\mu 0.6m_H10^{24}\mathrm{g}`$.

Una propiedad general de la teoría MIF es que, a medida que las fluctuaciones colapsan a escalas cada vez mayores, los halos deben mergerse y formar nuevos halos con mayor dispersión de velocidad; además, el ritmo de crecimiento de la masa y la dispersión de velocidad es mucho más rápido a alto CR que en el presente. Eso es una consecuencia directa de la forma del espectro de potencias en la Figura 1: a pequeñas escalas, la amplitud de fluctuaciones es prácticamente constante, y por lo tanto la masa de los halos aumenta rápidamente cuando esas fluctuaciones empiezan a colapsar. Pero a mayores escalas, la amplitud disminuye con mucha más rapidez con la escala, y por lo tanto las fluctuaciones a escalas grandes colapsan mucho más tarde.

## 4 Formación de galaxias

Una pregunta central en la teoría de formación de galaxias ha sido la siguiente: cuál es la causa de la diferencia entre una galaxia y un cúmulo de galaxias? Evidentemente, los cúmulos de galaxias han colapsado gravitatoriamente, pero deducimos a partir de observaciones de rayos X que la mayor parte de la materia bariónica está en forma de gas difuso y caliente distribuído en el halo. La masa total del gas difuso en un cúmulo masivo puede superar las $`10^{14}M_{}`$, pero no se forman nunca galaxias de masa tan grande.

Esa cuestión fue investigada por Binney (1977), Rees & Ostriker (1977), Silk (1977), y especialmente en el contexto de la teoría MIF por White & Rees (1978). En el colapso gravitatorio de un halo, la materia bariónica se calienta inicialmente en ondas de choque, hasta la temperatura necesaria para permanecer en equilibrio hidrostático. Luego, para poder condensar a una galaxia central, es preciso que el gas pierda su energía térmica por procesos de radiación. Si el tiempo necesario para que toda la energía térmica de los bariones sea radiada es corto comparado con el tiempo entre mersiones sucesivas de halos, entonces el gas va a poder concentrarse hacia el centro del halo, y el proceso de colapso gravitatorio puede proseguir hasta llegar a la formación de estrellas (despues de fragmentación de nubes de gas, y posiblemente de haber formado un disco debido a la conservación del momento angular; véase Fall & Efstathiou 1980, Fall & Rees 1985) o de un objeto masivo central, dando lugar a un núcleo activo. En cambio, si el tiempo de enfriamiento es demasiado largo, nuevas mersiones para formar halos de mayor masa van a calentar de nuevo el gas, manteniéndolo en estado difuso.

La condición que debemos analizar es, por consiguiente, que el tiempo de enfriamiento del gas sea igual a la edad del universo (que es aproximadamente el tiempo entre mersiones sucesivas). Los procesos de enfriamiento más importantes son la radiación bremsstrahlung, y la excitación y ionización colisional de iones. A temperaturas $`T10^5K`$, el enfriamiento por hidrógeno y helio domina, pero en el rango $`10^5K<T<10^7K`$ los iones de elementos pesados son importantes (p.e., Gaetz & Salpeter 1983). En las figuras 2(a,b), las líneas de puntos indican la dispersión de velocidad de halos en los que el tiempo de enfriamiento es igual a la edad del universo, para tres valores de la metalicidad del gas, $`(0,0.1,1)Z_{}`$. Por encima de las líneas de puntos, el tiempo de enfriamiento es demasiado largo y se predice que el gas caliente se mantiene en el halo, mientras que por debajo de las curvas el gas puede enfriarse y formar una galaxia. El ritmo de enfriamiento por radiación se ha calculado suponiendo que una fracción $`\mathrm{\Omega }_b/\mathrm{\Omega }`$ de la masa del halo está en forma de gas difuso, donde $`\mathrm{\Omega }_b=0.019h^2`$, de acuerdo con las mediciones de la abundancia de deuterio (Burles & Tytler 1998), y suponiendo el valor de la sobredensidad media de los halos en el modelo del colapso esférico (igual a $`18\pi ^2`$ para $`\mathrm{\Omega }=1`$). Si las galaxias pueden formarse en halos por debajo de las líneas de puntos, obtenemos la predicción de que hasta $`z4`$, el colapso de un halo debió resultar siempre en la formación de una galaxia central, pero a menores CR los halos más masivos pueden contener solamente galaxias formadas anteriormente, y la mayor parte del gas difuso no puede enfriarse. Mencionamos aquí que, incluso en los cúmulos de mayor masa, una parte del gas en la región central (con mayor densidad) puede radiar más rápidamente; la razón por la cual eso no da resultado a grandes tasas de formación estelar en las galaxias centrales de cúmulos masivos es una cuestión que permanece sin resolverse (Fabian, Nulsen, & Canizares 1994).

Vemos, pues, que la física de enfriamiento del gas proporciona una explicación satisfactoria del límite superior de la masa de las galaxias. Otra cuestión distinta es el límite inferior, y eso nos lleva a una transición fundamental en el universo para la historia de formación de las galaxias: la reionización del medio intergaláctico.

## 5 La reionización del universo

En la teoría de la Gran Explosión, la materia del universo debe formar átomos por vez primera cuando la temperatura de la radiación de fondo disminuye hasta $`T3000`$ K, a $`z10^3`$. La materia intergaláctica permanece posteriormente en estado atómico, hasta el momento en que el colapso no-lineal de las perturbaciones lleva a la formación de objetos que emiten radiación ionizante, y pueden ionizar de nuevo el gas difuso.

Las observaciones del espectro de fuentes luminosas a alto CR (generalmente, cuásares) demuestran que el medio intergaláctico fue ionizado anteriormente a $`z5`$. La luz de una fuente a longitudes de onda menores que la línea Ly$`\alpha `$ de hidrógeno atómico puede ser absorbida por un átomo a lo largo de su trayectoria a través del universo, en el punto en que la longitud de onda se ha corrido hasta coincidir con la línea Ly$`\alpha `$ . Si una parte importante de la densidad bariónica media del universo estuviera distribuída por el espacio en forma atómica, la profundidad óptica de absorción sería enorme ($`10^5`$ a $`z=3`$), con lo cual el flujo debería disminuir a cero abruptamente a la longitud de onda de Ly$`\alpha `$ (Gunn & Peterson 1965). Lo que se observa en realidad es la presencia de múltiples líneas de absorción que producen un decremento neto del flujo de solamente un $`30\%`$ a $`z=3`$. Eso implica que el medio intergaláctico está altamente ionizado, con una fracción neutra de $`10^5`$. Esa diminuta fracción del hidrógeno intergaláctico puede explicar satisfactoriamente las propiedades de las líneas de absorción observadas (denominadas usualmente como el “bosque Ly$`\alpha `$ ”), originadas en las variaciones de la densidad del gas que no ha colapsado todavía en halos virializados de gran densidad. Varios trabajos recientes, utilizando modelos analíticos y simulaciones numéricas, han mostrado que la teoría MIF predice de forma genérica que el medio intergaláctico ionizado debe dar lugar a ese bosque de líneas (véase McGill 1990; Bi, Börner, & Chu 1992; Bi 1993; Miralda-Escudé & Rees 1993; Cen et al. 1994; Zhang et al. 1995, 1998; Hernquist et al. 1996; Miralda-Escudé et al. 1996). El hecho de que la densidad del gas neutro es proporcional al cuadrado de la densidad del gas cuando se establece el equilibrio entre recombinación y fotoionización por el fondo cósmico ionizante es la causa de que un medio continuo da lugar a un espectro de absorción con la apariencia de “líneas” individuales, cada vez que se encuentra un máximo de la densidad del gas a lo largo de la línea de visual. El espectro de absorción tiene una escala de suavizamiento natural de $`20\mathrm{km}\mathrm{s}^1`$, debido a la dispersión de velocidad térmica del gas fotoionizado (con temperatura $`T2\times 10^4`$ K). El tamaño transversal de las estructuras de absorción Ly$`\alpha `$ debe ser del orden de su dispersión de velocidad multiplicada por la edad del universo, ya que tales estructuras no han tenido todavía tiempo de colapsar y de llegar a un equilibrio hidrostático; observaciones recientes del tamaño transversal en pares de cuásares (Bechtold et al. 1994; Dinshaw et al. 1994) confirman esa predicción.

Qué causó la reionización del universo? Dos mecanismos pueden ionizar el gas intergaláctico: fotoionización, o ionización colisional una vez el gas ha sido calentado por ondas de choque provenientes de alguna explosión de gran energía. La fotoionización es el método más eficiente para ionizar el medio de baja densidad, ya que requiere menos energía, y la radiación se transporta con gran eficacia a todas las regiones del espacio. Los dos tipos de fuentes de radiación ionizante conocidas que pueden formarse tan pronto como los primeros halos colapsan a pequeñas escalas son estrellas y núcleos activos (producidos por agujeros negros masivos en el centro de una galaxia). La cantidad de estrellas necesaria para ionizar el universo entero es solamente una parte muy pequeña de todas las estrellas que se han formado hasta el presente, y puede relacionarse fácilmente con la metalicidad media producida. Ya que tanto los elementos pesados provenientes de explosiones supernova como la radiación ionizante son producidas por estrellas masivas, la razón de esas dos cantidades está relativamente fijada: característicamente, una estrella de $`30M_{}`$ fusiona unas $`4M_{}`$ de hidrógeno en la secuencia principal, obteniendo una energía $`0.03M_{}c^2`$, de la cual $`0.01M_{}c^2`$ se emite en fotones ionizantes. Esa misma estrella produce unas $`5M_{}`$ de elementos pesados en la explosión supernova. Si los fotones son absorbidos por hidrógeno intergaláctico (necesitando una energía de $`20`$ eV para cada ionización, o una fracción $`2\times 10^8`$ de la masa-energía en reposo), la masa ionizada es de $`5\times 10^5M_{}`$, y por lo tanto la metalicidad media aumenta sólo hasta $`10^5`$ una vez se ha emitido un fotón ionizante para cada barión. Dado que la metalicidad media en el universo presente es mucho mayor, es evidente que la primera generación de estrellas puede fácilmente ionizar todo el universo (p.e., Couchman & Rees 1986).

Evidentemente, los cuásares pueden ser también las fuentes dominantes para la reionización, dado que la eficiencia en convertir la masa acretada por un agujero negro en radiación ionizante es generalmente mucho mayor que para estrellas. Las observaciones de la abundancia de cuásares y de la intensidad del fondo cósmico de radiación ionizante a $`z3`$ indican que los cuásares son probablemente los mayores contribuyentes de esa radiación cósmica, pero la naturaleza de las fuentes causantes de la reionización a mayor CR es todavía incierta (Miralda-Escudé & Ostriker 1990; Madau 1991, 1992, 1999; Haardt & Madau 1996; Rauch et al. 1997; Haiman & Loeb 1997, 1998; Miralda-Escudé, Haehnelt, & Rees 1999).

## 6 Las primeras galaxias

Las primeras galaxias donde pudieron formarse las primeras estrellas y cuásares surgieron, en la teoría MIF, en los primeros halos que colapsaron donde el gas pudo enfriarse. El hidrógeno atómico solamente empieza a radiar en líneas de excitación cuando la temperatura supera los $`10^4`$ K; a temperaturas inferiores, las colisiones con electrones térmicos no tienen nunca la energía necesaria para excitación a los niveles atómicos $`n=2`$. El hidrógeno molecular proporciona la única fuente de enfriamiento a menores temperaturas. En la materia primordial, en ausencia total de elementos pesados, el hidrógeno molecular se forma a partir del ión $`H^{}`$ (Saslaw & Zipoy 1967), el cual resulta a su vez de colisiones de hidrógeno con los electrones del residuo de ionización que permanece después de la época de recombinación (Peebles 1968). Ese proceso de formación de moléculas es muy ineficaz (en el presente, el hidrógeno molecular se forma en la superficie de granos de polvo interestelar), y sólo una pequeña fracción del hidrógeno forma moléculas. El enfriamiento resultante sólo es suficiente para resultar en la formación de galaxias cuando la temperatura de un halo supera los $`2000K`$ (Tegmark et al. 1997; Abel et al. 1998).

En las Figuras 2(a,b), la línea sólida a mayor CR indica la mínima temperatura para que un halo formado a CR $`z`$ pueda enfriar su gas y formar una galaxia, que hemos reproducido de Tegmark et al. (1997). Vemos que la teoría MIF, con los parámetros obtenidos de observaciones mencionadas antes, predice que las primeras estrellas pudieron formarse a $`z20`$, en halos de dispersión de velocidades de $`5\mathrm{km}\mathrm{s}^1`$ y masa total $`10^7M_{}`$. Esa primera generación de galaxias, formadas a través del enfriamiento por hidrógeno molecular, formó probablemente una cantidad muy pequeña de estrellas debido a la ineficacia de esta forma de enfriamiento (véase también Haiman, Rees, & Loeb 1997). En esta época, el ritmo de mersiones es muy elevado y la masa y dispersión de velocidad de los halos aumenta rápidamente con el tiempo; así pues, inmediatamente después de la formación de esas primeras estrellas, el enfriamiento atómico en halos de $`10^8M_{}`$ empieza a ser activo, dando lugar a una nueva generación de galaxias que debió reionizar el universo.

Una vez el medio intergaláctico es reionizado, el enfriamiento es suprimido por la fotoionización, que calienta el gas y disminuye la abundancia de átomos que pueden ser excitados colisionalmente. La línea sólida delgada a bajo CR nos da la temperatura en la cual el tiempo de enfriamiento es menor que el tiempo de Hubble para el gas ionizado (hemos adoptado aquí el modelo de Haardt & Madau 1996 para el fondo de radiación ionizante). Las galaxias pueden formarse solamente por encima de esta línea una vez la reionización se ha completado. En realidad, la fotoionización tiene también otro efecto: el calentamiento del gas durante la reionización puede producirse antes del colapso, y luego la temperatura del gas sube adiabáticamente cuando la densidad aumenta si el enfriamiento no es importante (recordemos que las líneas indicando un tiempo de enfriamiento igual a la edad del universo son para una sobredensidad de $`18\pi ^2178`$, el valor obtenido en el modelo esférico en el momento de la virialización, y que el tiempo de enfriamiento es mayor a menor densidad). Debido a este efecto, el gas en halos con $`\sigma 30\mathrm{km}\mathrm{s}^1`$ no puede en general enfriarse para formar galaxias (p.e., Thoul & Weinberg 1996). Otro efecto que puede disminuir la eficiencia de formación de galaxias en halos de baja dispersión de velocidad es el hecho de que la energía liberada por el propio proceso de formación estelar puede crear un viento galáctico con energía suficiente para calentar y expulsar el gas en el proceso de acreción (p.e., Dekel & Silk 1986). Ese proceso parece ser necesario para evitar un exceso de galaxias de baja luminosidad comparado con las observaciones (White & Frenk 1991; Navarro & Steinmetz 1997).

## 7 Observaciones de galaxias a alto CR

Grandes avances en cosmología observacional han tenido lugar recientemente, con el descubrimiento de gran número de galaxias a alto CR. Una de las técnicas de mayor éxito es la selección fotométrica de objetos débiles para detección de la brecha de Lyman (“Lyman break”; Guhathakurta, Tyson, & Majewski 1990; Steidel et al. 1996). Todas las galaxias, cuya luz resulta de la superposición de los espectros de muchas estrellas, muestran en su espectro una caída abrupta del flujo a la longitud de onda del límite de Lyman, como vemos en los modelos de espectros de distintas edades con ritmo constante de formación estelar en la Figura 3, que reproducimos de Bruzual & Charlot (1993). La brecha se produce en la atmósfera de las estrellas masivas, que producen la mayor parte de la radiación ultravioleta. Además, la presencia de hidrógeno atómico en el medio interestelar de la galaxia resulta generalmente en la absorción de la mayor parte de los fotones ionizantes, aumentando la amplitud de la brecha. A alto CR, la brecha de Lyman se corre hasta longitudes de onda en el visible, lo que permite seleccionar las galaxias a partir de sus colores. Por ejemplo, un objeto azul en V-R y muy rojo en B-V debe tener la brecha entre las bandas B y V, situándolo a $`z4`$.

El número de galaxias detectadas por este método, de magnitudes en el rango $`23B29`$, ha permitido la primera medición de la función de luminosidad y de la tasa global de formación estelar a alto CR, que resulta ser mucho mayor que la actual (Madau et al. 1996; Steidel et al. 1999). La cantidad total de estrellas formadas en estas galaxias, en el intervalo $`1z5`$, puede dar cuenta de la mayor parte de la población estelar vieja en el universo presente. Sin embargo, esas estimaciones están sujetas todavía a varias incertidumbres: en general, solamente la tasa de formación de estrellas masivas se puede deducir de esas observaciones, y existe la posibilidad de que la función de masa inicial fuera distinta a alto CR. Por otra parte, la tasa total de formación estelar podría ser mucho mayor si galaxias de baja luminosidad que permanecen por debajo de los límites de detección dominaran la emisión total, o si la mayor parte de la radiación ultravioleta es absorbida por polvo interestelar y reemitida en el infrarojo lejano (p.e., Calzetti 1999 y referencias incluídas).

Efectivamente, la radiación reemitida por polvo constituye otra forma de detección de galaxias a alto CR. Recientemente, la radiación de fondo en el infrarojo lejano, debida a la combinación de todas las galaxias, ha sido detectada por COBE (Fixsen et al. 1998; Hauser et al. 1998), y fuentes individuales se han descubierto con el nuevo detector SCUBA (Eales et al. 1998; Hughes et al. 1998; Smail et al. 1997), sugeriendo que gran parte de la radiación estelar a alto CR pudo ser procesada por polvo.

El origen de esas galaxias a alto CR puede comprenderse simplemente en la teoría MIF a partir de la Figura 4. Las líneas punteadas nos indican aquí la magnitud aparente de una galaxia (en el sistema $`AB`$, donde la magnitud es $`AB=48.62.5log_{10}(f_\nu )`$, y el flujo $`f_\nu `$ se expresa en unidades cgs), formada en un halo de dispersión de velocidad $`\sigma `$ en función de $`z`$, para un modelo de máxima luminosidad en que todos los bariones contenidos en el halo forman estrellas en un intervalo de tiempo igual a la mitad de la edad del universo en el momento del colapso. En otras palabras, el modelo representa la máxima eficiencia de formación estelar posible, donde toda la materia forma estrellas en un tiempo del orden del tiempo dinámico del halo de materia invisible (véase Miralda-Escudé & Rees 1998 para más detalles del modelo). Las galaxias de mayor masa pueden empezar a formarse a $`z4`$, cuando el gas puede enfriarse según la figura 2b. El flujo máximo de esas galaxias corresponde a $`AB22`$. La mayor parte de las galaxias de la brecha de Lyman son algo más débiles, como es de esperar cuando el ritmo de conversión de gas a estrellas es menos eficiente, y cuando una parte de la radiación ultravioleta es absorbida por polvo. Vemos también que las galaxias más masivas pueden formarse más fácilmente cuando el gas ha sido ya enriquecido, debido al mayor ritmo de enfriamiento a alta metalicidad.

Puesto que las galaxias detectadas a alto CR de mayor luminosidad están asociadas con las mayores fluctuaciones de densidad a las mayores escalas de colapso gravitatorio, su correlación espacial debería ser mucho mayor que la correlación de la masa, como se espera de los altos picos de densidad en un campo Gausiano (Kaiser 1984). Esa correlación ha sido detectada, y es en general consistente con las expectativas en la teoría MIF (Adelberger et al. 1998; Giavalisco et al. 1998; Kauffmann, Nusser, & Steinmetz 1997). Los estudios de la correlación de galaxias a alto CR abren un nuevo abanico enorme de posibilidades para observar la evolución de estructura a gran escala, que hemos solamente empezado a investigar.

A mayor $`z`$, la absorción del flujo ultravioleta a $`\lambda <1216\mathrm{\AA }`$ debido al bosque Ly$`\alpha `$ aumenta rápidamente, de tal forma que a $`z5`$, la caída de flujo en Ly$`\alpha `$ se convierte en la característica más importante para seleccionar objetos a este CR. Así pues, la técnica para encontrar galaxias a $`z>5`$ va a ser muy parecida a la de la brecha de Lyman, sustituyendo ésta por el hoyo de Gunn-Peterson (Gunn & Peterson 1965). Aunque inicialmente el término del hoyo de Gunn-Peterson se utilizó únicamente para referirse a la absorción en Ly$`\alpha `$ producida por el medio intergaláctico atómico antes de la reionización, el medio ionizado puede dar lugar también a un hoyo cuando la fracción neutra es suficientemente grande para que todas las partes del medio intergaláctico, incluso las de menor densidad en los vacíos, absorban esencialmente todo el flujo. Debido al aumento de la densidad media del gas en el universo, y al mayor ritmo de recombinación a alto CR, es de esperar que prácticamente todo el flujo a longitudes de onda menores que la línea Ly$`\alpha `$ sea absorbido a $`z6`$, incluso si la reionización ocurrió a mayor CR (Miralda-Escudé et al. 1999).

Finalmente, otra técnica importante de detección de galaxias a alto CR es mediante la línea de emisión de Ly$`\alpha `$ (Thompson et al. 1995; Thommes et al. 1998; Meisenheimer et al. 1998; Hu, Cowie, & McMahon 1998). En regiones de formación estelar, la mayor parte de la radiación ionizante emitida por estrellas jóvenes es generalmente absorbida por hidrógeno interestelar, y la energía se reemite en fotones de recombinación, siendo Ly$`\alpha `$ la línea más brillante. La búsqueda de líneas de emisión de objetos a alto CR tiene la ventaja de que el fondo de cielo puede reducirse observando solamente en una banda estrecha de longitudes de onda, especialmente cuando se seleccionan longitudes de onda de buena transparencia atmosférica (especialmente importante a $`z5`$, cuando la línea Ly$`\alpha `$ se corre hasta el infrarojo).

Es evidente a partir de las Figuras 2 y 4 que la dificultad en la detección de galaxias a CR progresivamente mayor aumenta rápidamente para $`z5`$, puesto que el flujo de las galaxias debe disminuir no sólo debido al mayor CR, sino a la baja luminosidad de las primeras galaxias, debida a su menor masa comparado con las galaxias actuales. Esa predicción de la teoría MIF sugiere que en el futuro, las fuentes detectadas a mayor CR podrían ser supernovas, las cuales debieron ocurrir tan pronto como las primeras estrellas se formaron (Miralda-Escudé & Rees 1997). La proyectada misión NGST podría detectar supernovas hasta $`z10`$ (Stockman & Mather 1997). Una posibilidad interesante que podría permitir acelerar el estudio observacional de las primeras estrellas es que los estallidos de rayos gamma, con sus brillantes contrapartidas ópticas (p.e., Metzger et al. 1997), sean un fenómeno asociado con estrellas masivas, ocurriendo por consiguiente a los mayores CR donde existían estrellas. Las contrapartidas ópticas serían detectables en observatorios terrestres, corridas al infrarojo, mientras que en el visual no habría contrapartida debido al hoyo de Gunn-Peterson y absorción por fotoionización. En el caso de los estallidos de rayos gamma, el reto observacional podría consistir en identificar una pequeña fracción de esos eventos a mayor CR que cualquier otra fuente conocida, y distinguirlos de estallidos a menor CR pero con un gran enrojecimiento debido a polvo interestelar cerca de la fuente, que pueden tener características fotométricas similares.

## 8 Conclusiones

La cosmología observacional está entrando en una nueva etapa de descubrimiento, con nuevas técnicas para la detección de galaxias débiles que empujan la frontera de alto CR hacia la época de la formación de las primeras galaxias. Al mismo tiempo, la medición precisa de las fluctuaciones en la radiación de fondo, la realización de nuevos escrutamientos de CR de galaxias, y la continuación de la búsqueda de supernovas a alto CR para la medición de la geometría cósmica, prometen fijar los parámetros del universo y del modelo cosmológico. Hemos visto que la teoría MIF reproduce con gran éxito las observaciones acumuladas hasta el presente sobre estructura a gran escala. Sin embargo, el estado actual de la teoría deja muchas preguntas abiertas: cuál es la naturaleza de la materia invisible? Existe realmente una “energía de vacío” que da cuenta de la densidad necesaria para alcanzar la densidad crítica? Cuál es la naturaleza de esta energía de vacío, cuál es su ecuación de estado, y por qué existe? Qué proceso generó las fluctuaciones de densidad? Qué determinó su amplitud? Son las fluctuaciones primordiales exactamente adiabáticas y Gausianas, y es su espectro perfectamente invariante de escala, o existen pequeñas desviaciones de esta simple hipótesis? Esas preguntas, que nos llevan al misterio de la época de inflación y al origen del universo, ocuparán probablemente el centro de interés en el futuro de la cosmología observacional, cuyo avance permitirá también descifrar la historia de la formación de galaxias, desde el colapso de las primeras estrellas hasta el presente.

Quisiera agradecer David Weinberg por proporcionar un código para calcular el espectro de potencias, y Gustavo Bruzual, Stephane Charlot y Max Tegmark por permitir la reproducción de resultados de sus trabajos en este artículo. Agradezco también diversas conversaciones con Martin Haehnelt, Martin Rees y David Weinberg.