Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Friedmann Equation

  1. Dec 27, 2008 #1
    Hey! I'd like to see an easy way to arrive to Friedmann equation by using Newton's Mechanics. I've seen many ways but they skip many steps so I don't understand the whole.

  2. jcsd
  3. Dec 27, 2008 #2
    Friedmann equations...

    As stated by Wikipedia, the Friedmann equations are derived from General Relativity, not from Newtonian mechanics. If you are interested in how these equations were derived, I recommend reading Friedmann's actual papers themselves, they are listed in reference 1 below, at the end of page in Reference as well as most of the important equations.

    http://en.wikipedia.org/wiki/Friedmann_equations" [Broken]
    Last edited by a moderator: May 3, 2017
  4. Dec 27, 2008 #3
    Re: Friedmann equations...

    Yes, but for people who doesn't have the level:tongue2: there is another way to get them, you can arrive by using Newtonian Mechanics, and when finally you have the constant then you have to compare with the real one and then you substitute your constant by the one that appears in the original one. That's the way I was asking for.

    Last edited by a moderator: May 3, 2017
  5. Dec 27, 2008 #4


    User Avatar
    Staff Emeritus
    Science Advisor
    Gold Member

  6. Dec 28, 2008 #5
    Friedmann Equation derived from Newton mechanics...

    The Friedmann_equation derived from Newton mechanics.

    This model describes a open finite Universe relativistic harmonic oscillator with a positive non-zero cosmological constant:
    [tex]\boxed{\Lambda = \frac{1}{dt^2} = 10^{-35} \; \text{s}^{-2}}[/tex]

    Harmonic oscillator force:
    [tex]F = -kx[/tex]

    Newtonian cosmological constant force:
    [tex]F_{\Lambda} = m \frac{d^2 r}{dt^2} = \frac{\Lambda m r}{3}[/tex]

    Harmonic oscillator potential energy:
    [tex]U = \pm \frac{dF \cdot dL}{2} = - \frac{kx^2}{2} = \frac{\Lambda m r^2}{6}[/tex]

    Co-moving coordinates:
    [tex]\boxed{r = a(t) x}[/tex]

    Kinetic energy:
    [tex]E_k = \frac{m \dot{r}^2}{2} = \frac{m \dot{a}^2 x^2}{2}[/tex]

    [tex]\boxed{E_k = \frac{m \dot{a}^2 x^2}{2}}[/tex]

    cosmological constant potential energy:
    [tex]U_{\Lambda} = \frac{\Lambda m r^2}{6}} = \frac{\Lambda m a^2 x^2}{6}}[/tex]

    [tex]\boxed{U_{\Lambda} = \frac{\Lambda m a^2 x^2}{6}}[/tex]

    Particle of mass m and radius r from a uniform expanding medium of density ρ.
    The total mass within radius r:
    [tex]m = \frac{4 \pi \rho r^3}{3}[/tex]

    Newton's universal law of gravitation:
    [tex]F = - \frac{G m^2}{r^2} = - \frac{4 \pi G \rho m r}{3}[/tex]

    Particle gravitational potential energy:
    [tex]V = - dF \cdot dL = - \frac{4 \pi G \rho m r^2}{3} = - \frac{4 \pi G \rho m a^2 x^2}{3}[/tex]

    [tex]\boxed{V = - \frac{4 \pi G \rho m a^2 x^2}{3}}[/tex]

    Relativistic harmonic oscillator potential energy:
    [tex]\boxed{U = - mc^2 \left( \frac{kx^2}{2} - \frac{\Lambda a^2 x^2}{6} \right)}[/tex]

    Particle energy conservation:
    [tex]U = E_k + V[/tex]

    Integration by substitution:
    [tex]- mc^2 \left( \frac{kx^2}{2} - \frac{\Lambda a^2 x^2}{6} \right) = \frac{m \dot{a}^2 x^2}{2} - \frac{4 \pi G \rho m a^2 x^2}{3}[/tex]

    Factor out m:
    [tex]- c^2 \left( \frac{kx^2}{2} - \frac{\Lambda a^2 x^2}{6} \right) = \frac{\dot{a}^2 x^2}{2} - \frac{4 \pi G \rho a^2 x^2}{3}[/tex]

    Multiply completely through by: [tex]\frac{2}{a^2 x^2}[/tex]
    [tex]- \frac{kc^2}{a^2} + \frac{\Lambda c^2}{3} = \left( \frac{\dot{a}}{a} \right)^2 - \frac{8 \pi G \rho}{3}[/tex]

    Solve for Friedmann equation:
    [tex]\boxed{H^2 = \left( \frac{\dot{a}}{a} \right)^2 = \frac{8 \pi G \rho}{3} - \frac{kc^2}{a^2} + \frac{\Lambda c^2}{3}}[/tex]

    http://zebu.uoregon.edu/lambda2.html" [Broken]
    http://en.wikipedia.org/wiki/Cosmological_constant" [Broken]
    http://en.wikipedia.org/wiki/Scale_factor_(universe)" [Broken]
    http://en.wikipedia.org/wiki/Friedmann_equations" [Broken]
    Last edited by a moderator: May 3, 2017
  7. Jan 1, 2009 #6
    Compton wavelength and general relativity...

    Does the Compton wavelength obey general relativity and special relativity?

    Universe quantum radius:
    [tex]\boxed{r_u = \overline{\lambda_u}}[/tex]

    [tex]m_u[/tex] - Universe mass

    Universe Compton wavelength:
    [tex]r_u = \overline{\lambda_u} = \frac{\hbar}{m_u c}[/tex]
    [tex]\boxed{r_u = \frac{\hbar}{m_u c}}[/tex]

    Universe quantum age:
    [tex]t_u = \frac{r_u}{c} = \frac{\hbar}{m_u c^2}[/tex]

    [tex]\boxed{t_u = \frac{\hbar}{m_u c^2}}[/tex]

    Universe quantum sphere average density:
    [tex]\rho_u(t_u) = \frac{m_u}{V_u} = \frac{3 m_u}{4 \pi r_u^3} = \frac{3 m_u}{4 \pi} \left( \frac{m_u c}{\hbar} \right)^3 = \frac{3 c^3 m_u^4}{4 \pi \hbar^3}[/tex]

    [tex]\boxed{\rho_u(t_u) = \frac{3 c^3 m_u^4}{4 \pi \hbar^3}}[/tex]

    Friedmann equation:
    [tex]H^2 = \left( \frac{\dot{a}}{a} \right)^2 = \frac{8 \pi G \rho_u(t_u)}{3} - \frac{kc^2}{a^2} + \frac{\Lambda c^2}{3}[/tex]

    Integration by substitution:
    [tex]H^2 = \left( \frac{\dot{a}}{a} \right)^2 = \frac{8 \pi G}{3} \left( \frac{3 c^3 m_u^4}{4 \pi \hbar^3}\right) - \frac{kc^2}{a(t_u)^2} + \frac{\Lambda(t_u) c^2}{3} = \frac{2 G c^3 m_u^4}{\hbar^3} - \frac{kc^2}{a(t_u)^2} + \frac{\Lambda(t_u) c^2}{3}[/tex]

    Quantum Compton-Friedmann equation:
    [tex]\boxed{H^2 = \left( \frac{\dot{a}}{a} \right)^2 = \frac{2 G c^3 m_u^4}{\hbar^3} - \frac{kc^2}{a(t_u)^2} + \frac{\Lambda(t_u) c^2}{3}}[/tex]

    http://en.wikipedia.org/wiki/Compton_wavelength" [Broken]
    http://en.wikipedia.org/wiki/Friedmann_equations" [Broken]
    http://en.wikipedia.org/wiki/General_relativity" [Broken]
    http://en.wikipedia.org/wiki/Special_relativity" [Broken]
    Last edited by a moderator: May 3, 2017
  8. Jan 4, 2009 #7
    Friedmann scale factor...

    Current Universe radius:
    [tex]\boxed{r_u(t_0) = 4.4 \cdot 10^{26} \; \text{m}}[/tex]

    Friedmann scale factor:
    [tex]\boxed{a(t) = \frac{r_u(t)}{r_u(t_0)}}[/tex]

    [tex]a(t_u) = \frac{r_u(t_u)}{r_u(t_0)} = \frac{\hbar}{m_u c r_u(t_0)}[/tex]

    [tex]\boxed{a(t_u) = \frac{\hbar}{m_u c r_u(t_0)}}[/tex]

    Quantum Compton-Friedmann equation:
    [tex]H^2 = \left( \frac{\dot{a}}{a} \right)^2 = \frac{2 G c^3 m_u^4}{\hbar^3} - \frac{kc^2}{a(t_u)^2} + \frac{\Lambda(t_u) c^2}{3}[/tex]

    Integration by substitution:
    [tex]H^2 = \left( \frac{\dot{a}}{a} \right)^2 = \frac{2 G c^3 m_u^4}{\hbar^3} - kc^2 \left( \frac{m_u c r_u(t_0)}{\hbar} \right)^2 + \frac{\Lambda(t_u) c^2}{3} = \frac{2 G c^3 m_u^4}{\hbar^3} - k \left( \frac{m_u c^2 r_u(t_0)}{\hbar} \right)^2 + \frac{\Lambda c^2}{3}[/tex]

    Quantum Compton-Friedmann equation:
    [tex]\boxed{H^2 = \left( \frac{\dot{a}}{a} \right)^2 = \frac{2 G c^3 m_u^4}{\hbar^3} - k \left( \frac{m_u c^2 r_u(t_0)}{\hbar} \right)^2 + \frac{\Lambda c^2}{3}}[/tex]

    [tex]\Lambda = \frac{1}{dt^2} = 10^{-35} \; \text{s}^{-2}[/tex]

    http://en.wikipedia.org/wiki/Scale_factor_(universe)" [Broken]
    http://en.wikipedia.org/wiki/Universe" [Broken]
    Last edited by a moderator: May 3, 2017
  9. Jan 5, 2009 #8


    User Avatar
    Gold Member

    Orion1: Thought your presentation of the equations in posts #5, #6 & #7 was excellent, but wanted to try to better understand the rationale behind some of the steps by way of clarification.

    In post #5:
    You start with a definition of Lambda = 1/s^2, but in the final Friedmann equation, the units have switched to 1/m^2 by virtue of including c^2. I know a lot of people like using natural units, but it can sometimes be misleading when trying to correlate units in the derivation.

    The 3rd equation was interesting and I followed the link that explained its premise. If you were to model a classical circular orbit, which accounted for the Lambda force, would the following form be correct?

    [1] [tex] \frac {mv^2}{r} + \frac{\Lambda c^2 m r}{3} = \frac{GMm}{r^2}[/tex]

    I am assuming that in an expanding universe, Lambda would act in opposition to the force of gravity. I raise this point because presumably it is important when we switch from force to energy in the later equations.

    As another point of clarification, would equating Lambda to dark energy imply that the effective mass (M), related to the energy density contained in the volume defined by [r], would increase with the expansion of the universe, if the density energy remains constant under expansion?

    Finally, I am also trying to understand how the concept of a SHM oscillator is being introduced. Without damping, an oscillator would just maintain the conservation of energy between kinetic and potential energy, but I shall also speculatively introduce the Lambda energy terms such that we might write:

    [2] [tex]dE_T = (dE_K + dE_\Lambda) + (-dE_P) = 0 [/tex]

    I realise this is a dubious conservation of energy assumption in an expanding universe and I am using the differential form simply because the equation ignores the energy put into the system to start it oscillating in the first place. However, if I do a very quick and dirty derivation based on the premise of equation [2], but bearing in mind the energy associated with Lambda:

    [3] [tex]E_K = E_P - E_\Lambda[/tex]

    In this form, the kinetic energy only relates to radial velocity and I am going to ignore that complexity of comoving coordinates by assuming [k=0] from the outset. I now substitute for Newtonian kinetic and gravitational energy plus your Lambda energy, but highlight the difference in sign:

    [4] [tex]1/2 mv^2 = \frac{GMm}{r} - \frac{\Lambda c^2 m r^2}{6} [/tex]

    I can divide through by (m) and substitute for (M=density*volume):

    [5] [tex] \frac{v^2}{2} = \frac{4 \pi G \rho r^2}{3} - \frac{\Lambda c^2 r^2}{6} [/tex]

    Leading to:

    [6] [tex] H^2 = \frac{v^2}{r^2} = \frac{8 \pi G \rho}{3} - \frac{\Lambda c^2 }{3} [/tex]

    I know the sign is wrong, but wanted to discuss this point as you have clearly been looking at these equations in some detail and may be able to put me straight on a couple of issues that have been puzzling me. If I initially ignore Lamba, equation [6] would reduce to the form

    [7] [tex] H^2 = \frac{v^2}{r^2} = \frac{8 \pi G \rho}{3} [/tex]

    Given that this equation seems to only relate the velocity [v] of our unit mass (m) to a gravitational term, why would [v] be described as an expansive velocity as normally associated with the Hubble constant (H)?

    If we re-introduce Lamba as an expansive energy (negative pressure); wouldn’t it act in opposition to the gravitational potential in an expanding universe?

    I will defer comment on #6 and #7 as this post is already a bit long, but would appreciate any clarification you might be able to offer on the points raised. Thanks.
  10. Jan 5, 2009 #9
    Classical Newtonian harmonic oscillator...

    Affirmative, the units have changed in the relativistic harmonic oscillator potential energy well.

    Non-relativistic classical Newtonian harmonic oscillator potential energy:
    [tex]U = U_k + U_{\Lambda} = - \frac{kx^2}{2} + \frac{\Lambda m r^2}{6}} = - \frac{kx^2}{2} + \frac{\Lambda m a^2 x^2}{6} = - \left( \frac{kx^2}{2} - \frac{\Lambda m a^2 x^2}{6} \right)[/tex]

    [tex]\boxed{U = - \left( \frac{kx^2}{2} - \frac{\Lambda m a^2 x^2}{6} \right)} \; \; \; \; \; \; \boxed{\Lambda = \frac{1}{dt^2}}[/tex]

    Relativistic harmonic oscillator potential energy:
    [tex]\boxed{U = - mc^2 \left( \frac{kx^2}{2} - \frac{\Lambda a^2 x^2}{6} \right)} \; \; \; \; \; \; \boxed{\Lambda = \frac{1}{dL^2}}[/tex]

    When the harmonic oscillator switches from a classical newtonian harmonic oscillator to a relativistic harmonic oscillator, the Systeme International units change for the cosmological constant.

    Newton-Einstein relative unit conversion:
    [tex]\boxed{\frac{\Lambda_N}{\Lambda_E} = \frac{dL^2}{dt^2} = c^2} \; \; \; \; \; \; \boxed{\frac{k_N}{k_E} = dF \cdot dL = \frac{mc^2}{a}} \; \; \; \; \; \; \boxed{\frac{\rho_N}{\rho_E} = 2}[/tex]


    Newton's universal law of gravitation:
    [tex]F_g = m \frac{d^2 r}{dt^2} = - \frac{G m^2}{r^2}[/tex]

    Newtonian centripetal force:
    [tex]F_c = m \frac{d^2 r}{dt^2} = \frac{mv^2}{r}[/tex]

    Harmonic oscillator force:
    [tex]F_k = m \frac{d^2 r}{dt^2} = - kx[/tex]

    Newtonian cosmological constant force:
    [tex]F_{\Lambda} = m \frac{d^2 r}{dt^2} = \frac{\Lambda m r}{3} \; \; \; \; \; \; \boxed{\Lambda = \frac{1}{dt^2}}[/tex]

    Newtonian force theorem:
    [tex]\boxed{F_g + F_c = F_k + F_{\Lambda}}[/tex]

    [tex]\boxed{F_g = -F_c + F_k + F_{\Lambda}}[/tex]

    [tex]F = m \frac{d^2 r}{dt^2} = - \frac{Gm^2}{r^2} = - \frac{mv^2}{r} - kx + \frac{\Lambda m r}{3}[/tex]

    [tex]\boxed{\frac{Gm^2}{r^2} = \frac{mv^2}{r} + kx - \frac{\Lambda m r}{3}}[/tex]

    Multiply completely through by [tex]\frac{r}{2}[/tex]:
    [tex]\frac{Gm^2}{2r} = \frac{mv^2}{2} + \frac{k r x}{2} - \frac{\Lambda m r^2}{6}[/tex]

    Kinetic energy:
    [tex]\boxed{E_k = \frac{mv^2}{2} = \frac{Gm^2}{2r} - \frac{k r x}{2} + \frac{\Lambda m r^2}{6}}[/tex]

    Multiply completely through by [tex]\frac{2}{m}[/tex]:

    [tex]v^2 = \frac{G}{r} \left( \frac{4 \pi \rho r^3}{3} \right) - \frac{k r x}{m} + \frac{\Lambda r^2}{3} = \frac{4 \pi G \rho r^2}{3} - \frac{k r x}{m} + \frac{\Lambda r^2}{3}[/tex]

    [tex]\boxed{v^2 = \frac{4 \pi G \rho r^2}{3} - \frac{k r x}{m} + \frac{\Lambda r^2}{3}} \; \; \; \; \; \; \boxed{a = \frac{r}{x}}[/tex] - co-moving coordinates

    Non-relativistic Newtonian Hubble parameter equation:
    [tex]\boxed{H^2 = \frac{v^2}{r^2} = \frac{4 \pi G \rho}{3} - \frac{k}{ma} + \frac{\Lambda}{3}} \; \; \; \; \; \; \boxed{\Lambda = \frac{1}{dt^2}}[/tex]
    Last edited: Jan 6, 2009
  11. Jan 9, 2009 #10
    Friedmann Equation...

    This 'Hubble-like' parameter is only a factor in the Friedmann equations relative only to the cosmological constant:

    [tex]\boxed{\Lambda = 1.112 \cdot 10^{-52} \; \text{m}^{-2}}[/tex]

    Friedmann equation:
    [tex]H^2 = \left( \frac{\dot{a}}{a} \right)^2 = \frac{8 \pi G \rho}{3} - \frac{kc^2}{a^2} + \frac{\Lambda c^2}{3}[/tex]

    The Hubble parameter:
    [tex]\boxed{H_0 = \sqrt{\frac{8 \pi G \rho}{3} - \frac{kc^2}{a^2} + \frac{\Lambda c^2}{3}}}[/tex]

    [tex]\boxed{H_0 = 2.26 \cdot 10^{-}^{18} \; \text{s}^{-}^{1}}[/tex]

    [tex]H_{\Lambda}^2 = \frac{\Lambda c^2}{3}[/tex]

    [tex]H_{\Lambda} = c \sqrt{\frac{\Lambda}{3}}[/tex]

    [tex]\boxed{H_{\Lambda} = 1.825 \cdot 10^{-18} \; \text{s}^{-1}}[/tex]

    http://en.wikipedia.org/wiki/Friedmann_equations" [Broken]
    Last edited by a moderator: May 3, 2017
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook