Cálculo didáctico de la edad del Universo y la

importancia de la constante cosmológica en un

modelo FRW

Rafael Hernández Jiménez1, Claudia Moreno2, Daniel Sánchez Guzmán3 y

Ricardo García Salcedo3
1
División de Ciencias e Ingenierías del Campus León de la Universidad de Guanajuato,
Loma del Bosque No. 103 Col. Lomas del Campestre C.P. 37150, León, Guanajuato,
México.
2
Centro Universitario de Ciencias Exactas e Ingeniería de la Universidad de
Guadalajara, Av. Revolución No. 1500, Sector Reforma, C.P. 44430, Guadalajara, Jal.,
México.
3
Centro de Investigación en Ciencia Aplicada y Tecnología Avanzada-Unidad Legaría
del Instituto Politécnico Nacional, Legaría #694. Col. Irrigación, C.P. 11500, México
D.F.

E-mail: fabulosorafa@fisica.ugto.mx; claudia.moreno@cucei.udg.mx;

dsanchezg@ipn.mx; rigarcias@ipn.mx

(Recibido el 2 de Febrero de 2011, aceptado el 11 de Septiembre de 2011)

Resumen
Partiendo del modelo cosmológico de Friedmann-Robertson-Walker (FRW), se puede encontrar una solución explícita
para el factor de escala a(t) utilizando la ecuación de estado de un fluido barotrópico y la ecuación de Friedmann. La
ecuación de estado describe varios tipos de materia (energía): materia no relativista, radiación y energía de vacío
(constante cosmológica), donde para cada componente existe una expresión explícita para a=a(t) que a su vez se
relaciona con una densidad de energía ρ=ρ(t). Se considera que cada etapa del Universo estuvo dominada por un tipo de
materia distinto y las expresiones que describen la evolución temporal de cada una de ellas son independientes unas de
otras. En este artículo se calcularán todos los casos para el tiempo transcurrido entre la famosa Gran Explosión y
nuestra época actual, para ello se utilizarán los datos experimentales tomados del artículo [1].

Palabras clave: Modelo cosmológico FRW, Edad del Universo, Observaciones cosmológicas.

Abstract
Starting with the Friedmann-Robertson-Walker (FRW) cosmological model, we can find an explicit solution for the
scale factor a(t) using the baryotropic equation of estat and the Friedmann equation. This barotropic equation describes
several types of matter (energy): Non-relativistic matter, radiation and vacuum energy (cosmological constant),
therefore exists an explicit expression for a=a(t) that at the same time relates to a density of energy ρ=ρ(t) for each
component. It is considered that each stage of the Universe was dominated by one type of material, and the expressions
that describe the time evolution are independent each other. By using experimental data from [1] the calculations of the
age of the Universe for each type of matter are carried out.

Keywords: Cosmological model FRW, Age of the Universe, Cosmological observations.

PACS: 01.55.+b, 04.20.Jb, 98.80.Jk, 98.80.Es ISSN 1870-9095

I. INTRODUCCIÓN (planetas) que sirvieron posteriormente a Johannes Kepler

para la formulación de las famosas leyes que gobiernan el
La cosmología estudia el nacimiento y evolución del movimiento planetario.
Universo. Antiguamente la cosmología estaba basada en las Actualmente, la cosmología está basada desde el punto
observaciones astronómicas realizadas a través de los de vista teórico en la Teoría de la Relatividad General
instrumentos más modernos de la época. Recordemos por (TRG) formulada por Albert Einstein (1915). La primera
ejemplo, la utilización de observatorios sin telescopios solución cosmológica a las ecuaciones de campo
utilizados por el astrónomo danés Tycho Brahe alrededor gravitacional de la TRG está basada en la idea de que el
de 1573, quien realizó innumerables observaciones de las Universo es muy parecido en todas las direcciones en que
posiciones de estrellas y de los llamados “errantes” se observe (isotropía), y que su composición es la misma en
Lat. Am. J. Phys. Educ. Vol. 5, No. 3, Sept. 2011 556 http://www.lajpe.org
Rafael Hernández Jiménez, Claudia Moreno, Daniel Sánchez Guzmán y Ricardo García Salcedo
todos lados (homogéneo) a escalas mayores a 100 que se encuentran de nosotros, están correlacionadas de tal
megaparsec1. Los argumentos anteriores dieron origen a lo forma que mientras más alejadas estuvieran de la Tierra,
que se conoce como el principio cosmológico [2-3], y el más rápidamente se alejaban de nosotros, este resultado se
modelo que resulta es conocido como la “Gran Explosión” conoce como la ley de Hubble [8].
o por su nombre en inglés “Big Bang” (BB) [4-5]; también Una vez que se conocieron los datos que aseguraban
es conocido como modelo de Friedmann-Robertson- que las galaxias se alejaban unas de otras, se entendió que
Walker. el Universo se está expandiendo [9]. Hubble encontró la
El modelo de FRW nos da una idea de la forma en la siguiente relación, conocida actualmente como ley de
que se originó el Universo, aunque no precisamente al Hubble:
tiempo cero, ya que para ello se requiere un periodo
conocido como inflación [6], el cual queda fuera del v  Hd , (1)
alcance de este artículo y que será tratado en un artículo
posteriormente. En este sentido, el modelo de FRW nos donde v es la velocidad de recesión entre las galaxias, d es
dice que el Universo comenzó con una densidad infinita de la distancia entre ellas y, H es la constante de Hubble. El
materia, en este caso constituida por materia ordinaria que valor para esta constante H fué determinada por Hubble, y
puede ser polvo no relativista (p=0) radiación km / s
electromagnética  p    o incluso lo que se conoce como tenia un valor de 500
Mpc
, es decir, por cada megaparsec
 3
constante cosmológica (p=-ρ) [4]. Al mismo tiempo, este de distancia la velocidad de recesión aumenta en 500km/s.
modelo nos predice tanto la cantidad como el tipo de La utilización de cada vez mejores métodos de medición de
átomos que surgieron de dicha densidad de energía distancias y velocidades ha llevado a que estos valores sean
(nucleosíntesis) [7]. cada vez más precisos, el valor más actual ha sido
Para los propósitos de este trabajo utilizaremos para la proporcionado por el satélite Wikilson Microwave
deducción de las ecuaciones dinámicas las ecuaciones de Anisotropy Probe (WMAP), el cual tiene un valor de
Einstein (EE), en donde consideraremos una métrica km / s
H 0  100 h , donde el subíndice 0 corresponde a
homogénea e isotrópica de curvatura constante y un tensor Mpc
de energía momento para un fluido perfecto [4]. Una vez considerar el valor actual (en esta época de la historia del
que sean obtenidas las ecuaciones dinámicas, estas se Universo en la que vivimos). El valor de la razón de
emplearan para calcular la edad del Universo en los expansión normalizada, h, es 0, 72±0, 03 [1].
diferentes tipos de materia mencionados previamente. La ley de Hubble es de gran importancia para tratar de
El presente artículo está organizado de la siguiente entender la historia del Universo, ya que ha dado pruebas
manera: En la sección II, se obtienen las ecuaciones contundentes de que el Universo se encuentra en una etapa
dinámicas del modelo cosmológico FRW a partir de la de expansión.
TRG. En la sección III se obtienen las soluciones a las Tomando en cuenta el principio cosmológico, vamos a
ecuaciones para ρ y a(t). La edad del Universo para los considerar una geometría que represente un espacio-tiempo
distintos tipos de materia se calcula en la sección IV, y de curvatura constante, homogéneo e isótropo. Esta
finalmente en la sección V se da una breve conclusión. geometría, la cual está representada por lo que se conoce
como métrica g  , es la de Robertson-Walker:
II. ECUACIONES BÁSICAS DEL MODELO
2  
COSMOLOGÍCO DE FRW ds  g  dx dx ,

2  dr 
 
2
(2)
r 2 d 2 sin 2  d 2  ,
2 2
En la década de 1920, Hubble se encontraba estudiando  c dt  a (t ) 
1kr
2
objetos cósmicos desde el Monte Wilson, los cuales se 
conocían como nebulosas o “universos isla”. Wilson
determinó la distancia a la galaxia de Andrómeda, 700.000 en la primera expresión un índice arriba y el mismo índice
años luz, mediante estrellas de brillo variable conocidas abajo implica suma sobre él,  ,  0, ..., 3, c es la velocidad
como cefeidas [8]. Asimismo, usando un método conocido de la luz, t es el tiempo cosmológico, a(t) es el factor de
como escala de distancias cósmicas [5], encontró las escala (magnitud que permite dar una medida de la
distancias a diferentes nebulosas que resultaron estar más evolución de las distancias entre dos puntos fijos en una
allá del tamaño de nuestra galaxia. Concluyó entonces que sección espacial dada), k es la curvatura de la sección
los “universos islas” están situados fuera de nuestra galaxia espacial dada y toma los valores 0,1 y ‐1 (espacio de
y se determinó que eran otras galaxias como la nuestra. curvatura cero, positiva y negativa), r la coordenada radial,
Debido a las observaciones de Wilson, Hubble θ la coordenada ángulo polar y φ la coordenada ángulo
determinó que la velocidad de recesión (velocidad con la azimutal. Visto de esta manera, el espacio‐tiempo admite
que se alejan de nosotros) de las galaxias y la distancia a la ser rebanado en cortes espaciales perpendiculares a la
dirección del tiempo cosmológico t, donde k es la curvatura
1 de cada corte espacial, homogéneo e isótropo. La forma de
El pársec o parsec (símbolo pc) es una unidad de longitud
equivalente a 3.26 años luz (3.0857x1016m).

Lat. Am. J. Phys. Educ. Vol. 5, No. 3, Sept. 2011 557 http://www.lajpe.org
 Cálculo didáctico de la edad del Universo y la importancia de la constante cosmológica en un modelo FRW
a(t) depende de las propiedades de la materia del universo, 
vacío como  vac  2 . La primera Ec. (6) es conocida
como lo veremos en la siguiente sección. 
El corrimiento al rojo es la relación entre la luz emitida como ecuación de Friedmann (EF), mientras que la segunda
por un objeto estelar y la luz de ese mismo objeto que de ellas es conocida como ecuación de Raychaudhuri (ER).
nosotros medimos, la cual es distinta debido a la expansión La ecuación de evolución del contenido material de nuestro
del Universo. En muchas ocasiones se utiliza el corrimiento Universo, la cual se conoce como ecuación de conservación
al rojo como parámetro de evolución del Universo. La de la energía, está dada por:
longitud de onda de la luz emitida se incrementa
proporcionalmente al factor de escala, cuyo efecto puede 
 T  0,
ser cuantificado por este corrimiento al rojo:

0 a0 donde es la derivada covariante, es decir, la derivada

1 z   . (3)
 a que es invariante ante transformaciones generales de
coordenadas. Sustituyendo (5) en la ecuación anterior se
La TRG de Ensitein fue formulada con el fin de generalizar obtiene que, para un fluido perfecto, la ecuación de
la teoría de la relatividad especial para sistemas inerciales, conservación es:
y considera sistemas de referencia no inerciales [7]. La
TRG esencialmente nos dice cómo es que el espacio-tiempo   3H    p   0. (7)
se deforma debido a la distribución de masa-energía
contenida en él, por lo que se considera una generalización De estas tres ecuaciones dinámicas, (6) y (7), dos son
de la teoría de la gravitación universal de Newton. independientes y la tercera se puede obtener de las otras
Las ecuaciones dinámicas de la TRG vienen dadas por dos. Estas representan el modelo conocido como FRW o de
las EE: la Gran Explosión y sus soluciones nos dicen cómo
evoluciona nuestro Universo. El modelo FRW también
1 2 puede obtenerse con argumentos que no requieren TRG [4-
R  Rg   g   4 T , (4)
2 c 5], por ejemplo. En la siguiente sección damos la solución
de estas ecuaciones.
Una cantidad útil para nuestros fines, es el parámetro de
donde k2=8πG, g  es el tensor métrico, R el tensor de
densidad cosmológica dado por:

curvatura de Ricci, R  g  R el escalar de curvatura de
Ricci,  es la constante cosmológica, c es la velocidad de la  2   vac
  2   , (8)
luz (en lo que sigue del texto se considerará con valor 1), G 3H 2
3H c c
es la constante de gravitación de Newton y, finalmente, el
tensor de energía-momento que representa la distribución que nos describe la densidad de materia total en el
de la materia-energía en el espacio-tiempo es: Universo, donde es conocida como la densidad critica:

T     p  u u  pg  , (5) 3H 2
c  . (9)
2
el cual en este caso, representa un fluido perfecto de
densidad ρ y de presión isotrópica p=p(ρ), ambas Actualmente, la densidad crítica tiene un valor de
cantidades son funciones que dependen solo de t. Además 26 2 3
c 0  1, 88  10 h kgm , donde h0.720.03 [1]. Esta
la 4‐velocidad local del fluido temporal está dada por
densidad corresponde aproximadamente a 3 átomos de
u  1,0,0,0  lo cual implica que T  diag (  , p, p, p ). hidrógeno por metro cúbico.
Sustituyendo (2) y (5) en las EE (4), solo se obtienen Tomando en cuenta la expresión (8), la ecuación de
dos ecuaciones independientes: Friedmann puede escribirse también como:

2  
2
a k k
H   
2
  ,   1  b     1  , (10)
a 3 a2 3 H 2a2
(6)
2 
  3 p   3 ,
a 
 donde    , corresponde a la densidad de energía de
a 6 3H 2
 2
donde H  H (t ) es el parámetro de Hubble, el punto sobre la constante cosmológica, y b  , corresponde a la
3H 2
la letra indica derivada con respecto a t y la constante
densidad de la materia no relativista y la radiación. El signo
cosmológica se relaciona con la densidad de energía del
de k está determinado por el valor de  y, por tanto, nos
indica a cuál de las tres geometrías de RW corresponde:
Lat. Am. J. Phys. Educ. Vol. 5, No. 3, Sept., 2011 558 http://www.lajpe.org
Rafael Hernández Jiménez, Claudia Moreno, Daniel Sánchez Guzmán y Ricardo García Salcedo
1. Si la densidad de materia    c 0 entonces se la expansión del Universo es acelerada. Tanto  como q
son parámetros adimensionales.
sigue que   1, y por lo tanto k=-1 lo cual representa
un Universo abierto desde el punto de vista geométrico.
2. Si la densidad de materia    c 0 entonces se III. SOLUCIONES DE LAS ECUACIONES
sigue que   1, y por lo tanto k  1 lo cual representa PARA  y a(t )
un Universo cerrado desde el punto de vista geométrico.
3. Si la densidad de materia    c 0 entonces se En esta sección se dan algunas soluciones para las
ecuaciones de evolución de una cosmología tipo FRW, Ecs.
sigue que   1, y por lo tanto k  0 el cual representa
(6) y (7). Estas soluciones pueden encontrarse en cualquier
un Universo plano desde el punto de vista geométrico. libro de cosmología, por ejemplo [7]; sin embargo,
El parámetro  es de suma importancia ya que nos podemos recomendar la referencia [4] para verlo de una
dice el tipo de materia que contiene el Universo. Así, si forma más didáctica.
tenemos varios tipos de materia o energía dentro de nuestro Recordemos que solo dos de las ecuaciones en (6) y (7)
Universo, entonces la expresión general de la Ec. (10) toma son independientes, obteniéndose de ellas un conjunto de
la forma: tres incógnitas (a(t), ρ(t), p(t)), por lo que debemos
especificar una ecuación adicional para que el sistema sea
k consistente. La ecuación a utilizar para tener un sistema
  1   i     1  , (11)
i H 2a2 completo, es una ecuación de estado que determina el típo
de fluido que contendrá el Universo, la cual está dada por:
donde la suma es sobre las diferentes especies de materia
que contiene el Universo. p   1  , (14)
A finales de la década de 1970, Vera Rubin y W. Kent
Ford, Jr. [10] descubrieron lo que ahora se conoce como donde  es conocido como parámetro barotrópico y puede
materia obscura, un tipo de materia que tiene efecto sobre el
tomar los siguientes valores 0 2 , por ejemplo, la
movimiento de las galaxias pero que no interactúa con
materia ordinaria. La densidad total de materia actual está radiación corresponde al valor  4/3 , la materia no
dada por: relativista tiene un valor de  1 y el valor para la energía
de vacío es  0 . En algunas ocasiones también se utiliza la
m  b  eo , (12) ecuación p , con  constante relacionada con  de la
siguiente forma:   1 .
donde  eo representa la densidad de materia obscura y  b Cuando se habla de materia no relativista nos referimos
al tipo de materia cuyas velocidades de translación están
la densidad de materia bariónica, que representa toda la muy por debajo de la velocidad de la luz. Este tipo de
materia que podemos observar y de la que estamos materia es la que nosotros vemos cotidianamente: planetas,
constituidos nosotros mismos. La discusión sobre la materia polvo estelar (sin incluir la radiación), nebulosas, galaxias,
obscura está fuera del alcance de este trabajo; sin embargo etc., de la misma forma, también se puede considerar la
ya se está trabajando para un próximo artículo. llamada materia obscura, que toma un papel importante a
Otro parámetro importante es el de deceleración nivel intergaláctico, y en la actualidad se le considera de
definido como: suma importancia para la descripción del movimiento de las
galaxias [10]. De manera observacional, la presión en este
aa H tipo de materia parece ser despreciable comparada con la
q  1  2 , (13)
a2 H densidad de energía, como lo hemos estado haciendo notar
a lo largo de este trabajo, es decir, p  0 , lo que implica
el cual mide la tasa de cambio de la velocidad de expansión que en la ecuación de estado (14).
del Universo. En principio es posible determinar el valor de Cuando hablamos de radiación nos referimos a los
q 0 de forma observacional. Por ejemplo, para un conjunto fotones que viajan desde cualquier rincón del Universo
idéntico de supernovas en una galaxia remota, la relación hasta nuestro planeta, toda radiación electromagnética: luz
entre su brillo y su corrimiento al rojo depende del valor del visible, ondas de radio, infrarrojo, ultravioleta, etc.
parámetro de deceleración. Aunque las medidas de esta Finalmente, la energía de vacio se considera como un
clase son muy difíciles de realizar e interpretar, las tipo de energía puramente cuántica, que existe incluso aún
observaciones recientes tienden a favorecer modelos en ausencia total de materia. Esta energía tiene una presión
cosmológicos que se aceleran [11], cuya explicación se negativa, lo que significa que si un fluido de presión
basa en los modelos conocidos como de energía obscura. positiva puede empujar las paredes del recipiente que lo
En observaciones recientes se encontró que contiene, en este caso tendría el efecto contrario, comprimir
el recipiente hacia dentro.
q0  0, 6138  0,1635 [1], lo que significa que actualmente

Lat. Am. J. Phys. Educ. Vol. 5, No. 3, Sept. 2011 559 http://www.lajpe.org
 Cálculo didáctico de la edad del Universo y la importancia de la constante cosmológica en un modelo FRW
Las soluciones a la Ec. (7) para la ecuación barotrópica  3p
(14) nos da la evolución de la densidad de energía: 0 1 
    3 
q0 =   0  1 . (20)
3 
2  2 
  a(t ) 
  , (15)
0  a0  Como podemos ver, el Universo en etapas tempranas de su
evolución estuvo dominado por la radiación, posteriormente
donde  0 y a0 son constantes de integración se fue enfriando y la materia no relativista comenzó a
dominar. Actualmente, el Universo está siendo dominado
correspondientes a los valores actuales de cada una de las por la energía de vacío (de acuerdo a las observaciones
variables. De esta forma, la densidad como función del [11]). A pesar del éxito de este modelo, tiene algunos
factor de escala para diferentes tipos de materia barotrópica problemas que se han resuelto si se toma en cuenta un
se expresa como: periodo de expansión acelerada en una etapa muy temprana
del Universo, conocida como inflación, para mayor
p  13    a 4 , información de esta etapa ver [6].
p 0   a 3 ,
, (16)
p     constante,
IV. CÁLCULO DE LA EDAD DEL UNIVERSO
para radiación electromagnética, materia no relativista y Desde el punto de vista observacional, podemos mencionar
energía de vacío respectivamente. que se ha determinado la edad de cúmulos globulares en la
La ecuación de Friedmann, primera Ec. (6), la podemos Vía Láctea en alrededor de 13, 5  2 mil millones de años
escribir sin pérdida de generalidad como: [12], lo anterior se ha realizado utilizando un método
independiente de la distancia. Por otro lado, si se utiliza un
a (t )   a (t ) .
2 2
(17) método de secuencia de enfriamiento de enanas blancas
[13-14], se obtiene una edad de 12, 7  0, 7 mil millones de
Al sustituir la densidad de energía obtenida en la Ec. (15) años en el cúmulo globular M4. Estos datos nos indican la
tenemos que el factor de escala como función de t es: cota mínima que se tiene para la edad del Universo:
2
t0  11  12 mil de millones de años. Los datos de 7 años de
a   t t0  3 , (18) la sonda WMAP determinaron, con un error de hasta el uno
por ciento, que la edad definitiva del Universo es de 13,73
y para los distintos tipos de materia que hemos considerado mil millones de años.
previamente se tiene que: Una forma muy sencilla de estimar la edad del Universo
es considerar que si velocidad de expansión es constante
desde su nacimiento, de tal forma que sustituyendo la ley de
a t t0    t 2 ,
1/2
Hubble (1) en la relación d=vt, se tiene que:
a t t0    t 2 ,
2/3
, (19)
ae Ht
  constante, t  H0 ,
1
(21)

para radiación, materia no relativista y energía de vacío parámetro que se conoce como el tiempo de Hubble. De las
respectivamente. observaciones del telescopio espacial Hubble y otros
La Ec. (19) nos dice que en el caso de Universo proyectos [15], el parámetro de Hubble está entre los
dominado por radiación o materia no relativista o, quizá, valores:
una combinación de ambos, la densidad de energía tiene
una singularidad al tiempo cero, es decir, la densidad de 1 1
energía es infinita en el famoso "Big Bang" o gran H0  9.776h , 0.64  h  0.80 , (22)
explosión, y conforme se va expandiendo el Universo esa
energía se va distribuyendo en toda su extensión. Como un en unidades de mil millones de años.
ejercicio para el lector, se pide encontrar las expresiones de Para calcular la edad del Universo tomando en cuenta el
la temperatura del Universo como función del tiempo y contenido material del Universo, tomemos la EF (6) junto
comprobar que mientras el Universo se expande (tiempo con la solución que se obtuvo para la ecuación de
cada vez más grande) su temperatura va disminuyendo. conservación (15):
El parámetro de deceleración (13) escrito para un fluido
H  H 0 0  a0 / a   0 k 0  a0 / a   ,
3
barotrópico toma la forma: 2 2 2
(23)

Lat. Am. J. Phys. Educ. Vol. 5, No. 3, Sept., 2011 560 http://www.lajpe.org
Rafael Hernández Jiménez, Claudia Moreno, Daniel Sánchez Guzmán y Ricardo García Salcedo
  Dependiendo del valor de la dendidad de energía  m , la
0
donde  1 ,   0 2 ,    0 2 y, en caso de que en
0 0

3H 0 3H 0 integral (29) puede ser distinta.

nuestro modelo contenga tanto radiación como materia, Una forma distinta de escribir la edad del Universo, es
   r   m , la ecuación de Friedmann toma la forma: como se propone en [16]. Para ello, despejamos a de la
definición del parámetro de deceleración (13) y lo
4 3 2 sustituimos en la EF (6):
H  H 0 r 0  a / a0  m0  a / a0  0 k 0  a / a0   . (24)
2 2

2
a 2  a 
k    H 0 12 q0  2 q0 0  , (30)
k 
0
Además,
H 0 2 a0 2
es el parámetro de densidad  0
a  a
asociado a la curvatura. Los subíndices 0 indican valores
donde   0,   0, y 2q0   m , además de considerar la
0
actuales. La Ec. (24) se obtuvo sustituyendo las relaciones
segunda de las Ecs. en (25).
4 3
r  a(t )  m  a(t )  La solución a esta Ec. (30) queda determinada por:
  ,   , (25)
r 0  a0  m 0  a0 
a

1 dx / a0
t0  , (31)
en la Ec. (23). H0 0 12q0  2q0 a0 / x
Ahora, podemos expresar H en términos del
corrimiento al rojo, Ec. (3), de la siguiente forma: donde x es una variable de integración y además a  0 en
t  0 . Si ahora tomamos w  x / a0 entonces (31) toma la
H  H 0 r 0 1 z  m0 1 z  0 k 0 1 z   .
2 2 4 3 2
(26) forma:

El siguiente paso consiste en derivar la Ec. (3) y sustituir la a


1 dw
Ec. (24) para escribir la edad del Universo como: t0  . (32)
H0 0 12q0  2q0 / w
t0 
 
dz
t0  dt  ,
0 0 H 1 z  Como podemos ver, tenemos dos expresiones equivalentes

(27) para el cálculo de la edad del Universo, Ecs. (29) y (32).

1 dz
 . Nuestro siguiente paso será tomar los valores que pudiera
H0 1 z  r 1 z   m 1 z    0  k 0 1 z 
0 4 0 3 2
0
tener  m y realizar la integral explícitamente para conocer
0

la edad del Universo.

Como podemos observar en la ecuación anterior, la escala
de tiempo para la edad del Universo es proporcional al
A.1 Universo plano  m  1
0
tiempo de Hubble; sin embargo, ahora aparecen
correcciones debidas a la presencia de algún tipo de
contenido material dentro del Universo. Este es un Universo cuyo contenido material está
En la siguiente sección, consideraremos casos enteramente dado por materia barotrópica, es decir,
particulares para el cálculo de la edad del Universo. 1
 m = 1 y q0 =
0
, por lo que la integral (29) queda como:
2
A. Época donde domina la materia no relativista
2 1 3/2
Consideremos un Universo lleno puramente de materia no t= H (1  z ) , (33)
3 0
relativista desde su nacimiento. Las soluciones a las
ecuaciones dinámicas son las siguientes:
entonces, la edad del Universo con z = 0 y los valores en
(22), está en el rango t0  8  10 mil millones de años, lo
p  0, a   t t0     t t0 
2/3 2
, ,
cual no satisface la cota mínima que muestran las
 0 (28) observaciones.
 r  0,    0,  k   m  1, q0 = m .
0 0 0 0

2
0
A.2 Universo cerrado  m  1
Entonces, la edad del Universo para este caso se calcula a
través de la siguiente integral, y usando las Ecs. (27) y (10):
Este caso no es, de acuerdo con las observaciones, el mejor

modelo para describir nuestro Universo, sin embargo, como

1 dz ejercicio haremos algunos cálculos que nos darán una idea
t0  . (29)
1 z  m0 z 1
2
H0 0 de la edad del Universo.

Lat. Am. J. Phys. Educ. Vol. 5, No. 3, Sept. 2011 561 http://www.lajpe.org
 Cálculo didáctico de la edad del Universo y la importancia de la constante cosmológica en un modelo FRW
La integral (31) en este caso, queda como: de donde observamos que H 0 t0  1, cuando  m  0, y
0

2
q0  sin   , cuando  m  1, lo cual nos da los límites
0
H 0 t0 
H 0t  , (34) 3
 
3
2 q0 1 correctos del modelo. Vemos que conforme decrece 0
m , la
edad del Universo va hacia el tiempo de Hubble. Las
si se realiza el cambio de variable [16]: observaciones [15] restringen la curvatura a un valor muy
cercano a cero, de tal forma que el Universo parece ser
 2q0 1 a plano, es decir  m 0 1, esta edad no es más grande que las
1  cos   . (35)
q0 a0 estrellas más viejas, por lo que tampoco es muy buen
modelo para explicar las observaciones.
Es muy fácil darse cuenta que las Ecs. (34) y (35) son las El problema del tiempo, puede ser resuelto si se
ecuaciones paramétricas de una cicloide para la solución considera un Universo que contenga tanto materia barionica
a (t ) . Una cicloide es la curva que genera un punto fijo en como constante cosmológica, suponiendo un modelo de
una circunferencia que gira sin deslizarse horizontalmente. curvatura cero. En este caso, la integral (27) queda como:
En esta curva a (t )  0 cuando   0 , y va creciendo hasta


que en θ=π alcanza el máximo y después decrece 1 dz
t0  ,
H0 1 z   m 1 z    0
3
nuevamente a(t)=0 cuando θ=2π. El tiempo Tm, cuando se 0 0

alcanza el valor máximo de a(t), se puede obtener mediante (39)

el criterio de la primera derivada, así como el valor del 2  1  0 
 ln   ,
factor de escala en ese momento a(Tm): 3  0   0 
 m 

 
3
Tm =  q0 H 0 2 q0 1 ,
donde  m    1 . En este caso, los valores asintóticos
0 0
(36)
1
a (Tm ) = 2 q0 (2 q0  1) a0 .
 m  0,
0
corresponden a que H 0 t0   cuando y
La edad actual del Universo se obtiene al considerar θ= θ0, 2
cuando  m  1. Por ejemplo, en un modelo
0
H 0 t0 
a=a0 en la Ec. (35) y sustituir en (34): 3
m  0.3,   0.7,
0 0
donde se tiene que
t  H0
1 q0 cos 1  q 1 1q 1 2 q 1 , (37) t0  0, 964 H 0
1
 13,1 mil millones de años si h0.72 . Este
   
3 0 0 0
2q0 1 valor, nos permite evitar la restricción impuesta
anteriormente, por lo que podemos concluir que la inclusión
de una constante cosmológica resuelve el problema de la
en cuya ecuación podemos dar valores adecuados para q0 y
edad del Universo.
h, y de esta forma obtener el valor actual de la edad de un Así como se ha calculado la edad del Universo para
Universo cerrado. Por ejemplo, si se consideran q0  2 y contenido de materia no relativista, también podemos
3 hacerlo para cuando existe un dominio de radiación a lo
h=0,72 entonces la edad de este Universo sería de t0  8, 53 largo de toda su evolución, es decir:
mil millones de años, el cual es un valor que queda por
a   t t0  ,    t t0 
1/2 2
debajo de la cota dada por la edad de algunas estrellas. ,
(40)
 r  1,    0,  m  0, q0 =  r
0 0 0 0
A.3 Universo abierto  m  1
0 ,

Este parece ser el caso más real debido a que las sin embargo, este cálculo lo dejaremos para el lector. Por
observaciones han mostrado que existen otros tipos de ahora, es suficiente mencionar que la edad para este modelo
elementos adicionales, además de la materia ordinaria, que es:
1 1
contiene el Universo. Entonces, tomando  m  1 y
0
t0 H , (41)
2 0
1
q0  , de tal forma que la integral (29) queda definida
2 el cual es un valor que no podría explicar la edad de las
como: estrellas [12-14].

1 m 0 1 1 0  B. Época donde domina energía de vacío, constante

H 0 t0 =  ln  m 
cosmológica,  0
 
1m0 2 1 0 , (38)
1 1m0 
3

Lat. Am. J. Phys. Educ. Vol. 5, No. 3, Sept., 2011 562 http://www.lajpe.org
Rafael Hernández Jiménez, Claudia Moreno, Daniel Sánchez Guzmán y Ricardo García Salcedo
Con anterioridad se ha introducido la constante cosmológica. Una vez obtenidas las soluciones, se procedió
cosmológica en un modelo donde tambien se contenia a utilizar la ecuación de Friedmann (6) para encontrar una
materia barotrópica e incluso materia obscura, y se obtuvo expresión general para la edad del Universo.
un resultado adecuado a las observaciones. En esta sección, Posteriormente, se analizaron los casos en donde domina la
y solo como análisis didáctico, se calcula la edad del materia no relativista, toda aquella que podemos observar
Universo cuando éste está dominado por la constante en sus diferentes opciones, dependiendo del valor del
cosmológica a lo largo de toda su evolución, es decir, un 0
parámetro densidad de energía  m ya sea mayor, menor o
Universo que se expande aceleradamente desde su origen.
Es interesante hacer la observación de que una vez que igual que cero.
0
Einstein publicó su TRG se dedicó a hacer consideraciones Cuando se considera el valor  m  1 se obtiene un
cosmológicas. Descubrió que era imposible construir un valor para la edad del Universo menor que el que presentan
modelo matemático de un Universo estático, coherente con las observaciones de estrellas en diversos cúmulos
el estado en ese momento de las observaciones (recordemos globulares de nuestra galaxia.
que la expansión del Universo se descubrio hasta finales de 0
la década de 1920 por parte de Hubble). Sin embargo, Cuando se tiene un valor  m  1 se observa que la
Einstein estaba convencido de que el Universo era estático, edad de este Universo queda igualmente por debajo de la
por lo que introdujo una constante, conocida ahora como cota de la edad de algunas estrellas.
constante cosmológica  . El objetivo de dicha constante 0
El último caso que falta es aquel donde  m  1 el cual
consistía en ofrecer una fuerza de atracción contra la
solución expansiva que se obtenía. Actualmente, a esta es el más interesante ya que permite que el Universo no
constante cosmológica se le da una interpretación distinta, solo contenga materia barionica o materia obscura, sino que
ya que ahora puede interpretarse como una presión negativa brinda la oportunidad de coexistencia de esta materia con
que contrarresta la fuerza de atracción y, además, puede alguna otra de tal forma que la densidad de energía total del
explicar la expansión acelerada que se observa [15], ésta Universo sea 1. Por lo que en este caso, si se considera que
constante tiene características parecidas a la energía del hay coexistencia de materia bariónica y/o obscura junto con
0
constante cosmológica tal que  m  0, 3 y    0, 7 , la
0
vacío.
Así, considerando el modelo de un Universo plano edad del Universo sería más cercana a la observada, es
 k 0 0 que contiene sólo energía del vacío  0 1 , se decir, t0  13,1 mil millones de años. Por tanto, aquí es
tiene que la edad actual de este universo es: razonable pensar en la existencia de dicha constante
cosmológica o energía de vacío en nuestro Universo.
1 1/2 1 1/2  Algunas veces esta constante cosmológica es reemplazada
t0  H 0 
2
ln  1/2  . (42)
1 
3 por modelos un poco más exóticos de energía obscura.
Finalmente se mencionan los resultados para un
Universo dominado completamente por radiación y otro por
Es interesante notar que a diferencia de los otros dos
constante cosmológica, en donde se observan resultados
modelos anteriores, un Universo con    0, 74 es más
0
que no son del todo adecuados, de acuerdo a las
1 observaciones.
viejo que H 0 , esto se debe a que la razón de expansión se
Un complemento a este trabajo, sería explorar algún
acelera. Y, por supuesto, en el límite cuando , modelo que pudiera conectar las cuatro etapas que
km / s
t0   . En particular, si    0, 74 y H 0  72 Mpc
, la actualmente se reconocen en la evolución del Universo: una
etapa de inflación (de expansión acelerada a tiempo
edad de este Universo dominado por la energía del vacío
temprano), una época de dominio de la radiación, otra etapa
sería de 15, 21 mil de millones de años que, comparado con
de dominio de materia no relativista y, finalmente, una
el tiempo de Hubble, este Universo sería más viejo que el etapa de coexistencia de materia no relativista y energía
tiempo de Hubble. obscura.

V. CONCLUSIONES AGRADECIMIENTOS
El modelo FRW nos da una idea de cómo surgió el Agradecemos los comentarios invaluables de los arbitros
Universo, aunque para ello habría que anexarle el llamado anónimos de este trabajo, así como a César Mora. Este
modelo inflacionario [6]. Este modelo nos dice que el trabajo fue parcialmente apoyado CONACyT México, a
Universo comenzó con una densidad infinita de materia través de los proyectos No. 49924-J y 105079; también a
(19), en este caso constituido por materia ordinaria que SNI-CONACyT. R. G-S. Agradece el apoyo de las becas
puede ser polvo no relativista o radiación electromagnética COFAA y EDI, así como de los proyectos del IPN SIP-
e incluso, constante cosmológica. En este artículo 20100610, y SIP-20110664. R. H. y C. M. agradecen a
revisamos las EE y sus soluciones para tres tipos distintos ProSNI-UdG.
de materia barotrópica que incluye a la constante

Lat. Am. J. Phys. Educ. Vol. 5, No. 3, Sept. 2011 563 http://www.lajpe.org
 Cálculo didáctico de la edad del Universo y la importancia de la constante cosmológica en un modelo FRW
REFERENCIAS [10] Rubin, V. C. y Ford, W. K., Rotation of the
Andromeda Nebula from a Spectroscopic Survey of
[1] Nakamura, K., et al., Review of Particle Physics, Emission Regions, Astrophys. J. 159, 379 (1970).
Journal of Physics G, Nuclear and particle physics 37, [11] Riess, A. G. et al., Supernova Limits on the Cosmic
075021 (2010). Equation of State, Astrophys. J. 509, 74-79 (1998) [astro-
[2] Kelvin, K., Wu, S., Ofer, L. & Rees, J. M., The large- ph/9806396]; S. Perlmutter et al., Measurements of Omega
scale smoothness of the Universe, Nature 397, 225-230 and Lambda from 42 High-Redshift Supernovae,
(1999), [arXiv:astro-ph/9804062v2]. Astrophys. J. 517, 565-586 (1999) [astro-ph/9812133].
[3] Ofer, L., How Smooth is the Universe on Large Scales?, [12] Jiménez, R. et al., Ages of globular clusters: A new
(1998) [arXiv:astro-ph/9807093v1]. approach, Mon. Not. R. Astron. Soc. 282, 926-942 (1996).
[4] García, S. R. y Moreno, C., Descripción de la evolución [13] Ritcher, H. et al., The lower main sequnce and mass
del Universo: Una presentación para alumnos function of the globular cluster Messier 4, Astrophys. J.
preuniversitarios, Lat. Am. J. Phys. Educ. 1, 95-100 574, L151-L154 (2002).
(2007). [14] Hansen, B. et al., The white dwarf cooling sequence of
[5] Mora, C., Deducción de los primeros modelos the globular cluster messier 4, Astrophys. J. 574, L155-
cosmológicos, Lat. Am. J. Phys. Educ. 2(2), 180-180 L158 (2002).
(2008). [15] Freedman, W. L. et al., Final results from the Hubble
[6] Guth, H. A., El universo inflacionario: La búsqueda de Space Telescope key project to measure the Hubble
una nueva teoría sobre los orígenes del cosmos, (Editorial constant, Astrophys. J. 553, 47-72 (2001); WMAP Collab.
Debate, Madrid, 1999). (D. N. Spigel, et al.), Wilkinson Microwave Anisotropy
[7] Carroll, S. M., Spacetime and Geometry: An Probe (WMAP) Three Year Results: Implications for
Introduction to General Relativity, (Addison-Wesley, San Cosmology, Astrophys. J. Suppl. 170, 377 (2007) [astro-
Francisco CA, 2003). Coles, P. y Lucchin, F., Cosmology. ph/0603449`]; SDSS Collab. M. Tegmark et al.
The origin and evolution of cosmic structure, 2a Ed. (John Cosmological parameters from SDSS and WMAP, Phys.
Wiley & Sons, LTD, Great Britain, 2002). Rev. D69, 103501 (2004); U. Seljak et al.,Cosmological
[8] Hubble, E., A Relation Between Distance nd Radial parameter analysis including SDSS Ly-alpha forest and
Velocity Among Extra-Galactic Nebulae, Proceedings de la galaxy bias: Constraints on the primordial spectrum of
Academia Nacional de Ciencias de EEUU 15, 168-173 fluctuations, neutrino mass and dark energy, Phys. Rev.
(1929): http://www.pnas.org/cgi/reprint/15/3/168. D71, 103515 (2005).
[9] Lineweaver, C. H. and Davis, T. M., Misconceptions [16] Islam, J. N., An introduction to Mathematical
about the Big Bang, Scientific American March 2005, 36- cosmology, 1a Ed. (Cambridge U. P., Cambridge, 1992).
45 (2005).

Lat. Am. J. Phys. Educ. Vol. 5, No. 3, Sept., 2011 564 http://www.lajpe.org