next up previous
Next: Scheme of integration Up: PM Code Previous: Introduction

Equations and dimensionless variables

The equations of motion of particles in expanding coordinates, which are used in our PM code, were presented by Kates et al.(1991). Different numerical effects (including resolution) were discussed in Klypin et al.(1996). We use comoving coordinates ofr particles x = x (t) , which are related to proper coordinates by r = a(t)x , where a(t) = (1 + z)- 1 is the expansion parameter. Instead of using peculiar velocity v pec = a$\dot{\mathbf}$x we write the equations for particle momenta p :


  \begin{figure}
% latex2html id marker 49
\begin{equation}
\mathbf{p}=a^2\dot\mathbf{x}, \quad \mathbf{v}_{\rm pec}=\mathbf{p}/a\end{equation}\end{figure}

This choice of ``velocity'' simplifies the equations of motion by removing a few terms with $\dot{a}$/a . It is also convenient to change the time variable from time t to the expansion parameter a . The equations governing the motion of particles are:

 
$\displaystyle{d\mathbf{p}\over da}$ = - $\displaystyle{\nabla\phi\over\dot a}$$\displaystyle{d\mathbf{x}\over da}$ = $\displaystyle{\mathbf{p}\over \dot a a^2}$ (2)

 
$\displaystyle\nabla^{2}_{}$$\displaystyle\phi$ = 4$\displaystyle\pi$G$\displaystyle\Omega_{\rm m}^{}$(t)a 2$\displaystyle\rho_{\rm cr}^{}$(t)$\displaystyle\delta$ = 4$\displaystyle\pi$G$\displaystyle\Omega_{0}^{}$$\displaystyle\rho_{\rm cr,0}^{}$$\displaystyle{\delta\over a}$$\displaystyle\delta$ $\displaystyle\equiv$ $\displaystyle{\rho(\mathbf{x})-\rho_b\over\rho_b}$, (3)

 
$\displaystyle\dot{a}$$\displaystyle\sqrt{a}$ = H0$\displaystyle\sqrt{\Omega_0 
 +\Omega_{\rm curv,0}a
 +\Omega_{\Lambda,0}a^3}$$\displaystyle\Omega_{0}^{}$ + $\displaystyle\Omega_{\rm curv,0}^{}$ + $\displaystyle\Omega_{\Lambda,0}^{}$ = 1, (4)

where $\Omega_{0}^{}$ = $\Omega_{m}^{}$(z = 0) , $\Omega_{\rm curv,0}^{}$ , and $\Omega_{\Lambda,0}^{}$ are the density of the matter, effective density of the curvature, and the cosmological constant in units of the critical density at z = 0 . The curvature contribution is positive for negative curvature.

Dimensionless variables (shown with tildas below) are defined by introducing the length of a cell of the grid x0 and by measuring the time in units of 1/H0 :


  \begin{figure}
% latex2html id marker 89
\begin{equation}
\mathbf{x} = x_0\tilde\mathbf{x}, \quad t=\tilde t/H_0,\end{equation}\end{figure}


  \begin{figure}
% latex2html id marker 101
\begin{equation}
\mathbf{v}_{\rm pec}=...
 ...0)\tilde\mathbf{p}/a, 
 \quad\phi=\tilde\phi(x_0H_0)^2\end{equation}\end{figure}

 
$\displaystyle\rho$ = $\displaystyle{\tilde\rho\over a^3}$$\displaystyle{3H_0^2\over 8\pi G}$$\displaystyle\Omega_{0}^{}$ (7)

Equations (2-5) can be rewritten in terms of dimensionless variables:

 
$\displaystyle{d\tilde\mathbf{p}\over da}$ = - F(a)$\displaystyle\tilde{\nabla}$$\displaystyle\tilde{\phi}$$\displaystyle{d\tilde\mathbf{x}\over da}$ = F(a)$\displaystyle{\tilde\mathbf{p}\over a^2}$ (8)

 
$\displaystyle\tilde{\nabla}^{2}_{}$$\displaystyle\tilde{\phi}$ = $\displaystyle{3\over 2}$$\displaystyle{\Omega_0\over a}$($\displaystyle\tilde{\rho}$ - 1), (9)

where

 
F(a) $\displaystyle\equiv$ H0/$\displaystyle\dot{a}$ = $\displaystyle\left({ 
 \Omega_0 +\Omega_{\rm curv,0}a+\Omega_{\Lambda,0}a^3}
 \over a \right)^$-1/2 . (10)

Equations (8 - 9) are solved numerically by the PM code.

If L is the length of the computational box at z = 0 , Ngrid is the number of grid cells in one direction, and Nrow is the number of particles in one direction, which contribute a fraction $\Omega_{0}^{}$ of the critical density, then the transformations from dimensionless variables given by the code to dimensional variables are given by

 
x = x0$\displaystyle\tilde{x}$x0 = $\displaystyle{L\over N_{\rm grid}}$ = 7.8kpc $\displaystyle\left({L_{\rm Mpc}\over N_{\rm grid}/128}\right),$ (11)

 
vpec = (x0H0)$\displaystyle{\tilde p\over a}$ = 0.781$\displaystyle{\rm km\over s}$ $\displaystyle\cdot$ $\displaystyle{\tilde p\over a}$ $\displaystyle\cdot$ $\displaystyle{L_{\rm Mpc}h\over N_{\rm grid}/128}$, (12)

 
Mass = Nparticles $\displaystyle\cdot$ m1m1 = $\displaystyle\Omega_{0}^{}$$\displaystyle\rho_{\rm cr,0}^{}$$\displaystyle\left({L\over N_{\rm row}}\right)^$3 = 1.32 $\displaystyle\cdot$ 105($\displaystyle\Omega_{0}^{}$h 2)$\displaystyle\left({L_{\rm Mpc}\over N_{\rm row}/128}\right)^$3 (13)


next up previous
Next: Scheme of integration Up: PM Code Previous: Introduction
Jon Holtzman
12/9/1997