Dear all:
I want to get the \psi and \phi which are the potentials of Newtonian Gauge from CAMB. How can I do that ? Thanks in advance.
Cheers,
Gongbo Zhao
CAMB: Output the \psi and \phi

 Posts: 69
 Joined: January 04 2005
 Affiliation: ICG, Portsmouth
 Contact:

 Posts: 1435
 Joined: September 23 2004
 Affiliation: University of Sussex
 Contact:
Re: CAMB: Output the psi and phi
Equations relating \Psi and \Phi to the synchronous gauge quantities used by CAMB are given in e.g. Eqs. A1,A2,etc of astroph/0212248. Usually it is easiest to use the Poisson equation
[tex]k^2\Phi = \kappa S^2(\delta\rho + 3Hq/k)/2[/tex]
The total matter perturbation ([tex]\kappa S^2\delta\rho[/tex]) and heat fluxe (\kappa\S^2\q) are called dgrho and dgq in CAMB's equations.f90. You could write them out in the fderivs or output routine.
Note however that the result for \Phi is numerically unstable on large scales, so if you require an accurate result on large scales (and an analytical one won't do) you may have to evolve the \Phi equation separately. (The early time series solution for \Psi is given in the notes).
Note also that I have seen every possible different sign convention for both [tex]\Psi[/tex] and [tex]\Phi[/tex]!
[tex]k^2\Phi = \kappa S^2(\delta\rho + 3Hq/k)/2[/tex]
The total matter perturbation ([tex]\kappa S^2\delta\rho[/tex]) and heat fluxe (\kappa\S^2\q) are called dgrho and dgq in CAMB's equations.f90. You could write them out in the fderivs or output routine.
Note however that the result for \Phi is numerically unstable on large scales, so if you require an accurate result on large scales (and an analytical one won't do) you may have to evolve the \Phi equation separately. (The early time series solution for \Psi is given in the notes).
Note also that I have seen every possible different sign convention for both [tex]\Psi[/tex] and [tex]\Phi[/tex]!