- #1
fab13
- 318
- 6
Hello,
i wrote a Matlab program which calculate the apparent angle of an anisotropy with a size equal to "380.000 light-years" , i.e "0.1166 Mpc".
For this, i have used the following equation defining the apparent angle ( in arcminute unit) as a function of the redshift "z" with "Danisotropy=size of the anisotropy" :
[tex] \theta= \dfrac{180*60}{\pi}\,D_{anisotropy}\,\big[\dfrac{c}{H_{0}\,(1+z)} \ ,\int_{z_{begin}}^{z} \,(\Omega_{m}\,(1+z)^{3}+\Omega_{\Lambda}+\Omega_{k}\,(1+z)^{2})^{-1/2}\,dz\big]^{-1} [/tex]
For [tex]\Omega_{m}=0.3[/tex] and [tex]\Omega_{\Lambda}=0.7[/tex], i get the figure in attachment. For "z=1100", I have an angle apparent equal to about 30 arcminutes for Danisotropy=380.000 ly , i.e 0.5 degrees.
This value is relatively quite far from the first peak in CMB power spectrum which is localized at "l=200", i.e "theta= 1 degree".
Does the relation between multipole "l" and "theta" angle take the form : "theta(°)=200/l" ?
If it does, then my calculation is the value of one of the others peaks which are localized for smaller angles ?
If anyone could explain me to which peak corresponds this result,
Thanks in advance.
ps: here is the Matlab code :
i wrote a Matlab program which calculate the apparent angle of an anisotropy with a size equal to "380.000 light-years" , i.e "0.1166 Mpc".
For this, i have used the following equation defining the apparent angle ( in arcminute unit) as a function of the redshift "z" with "Danisotropy=size of the anisotropy" :
[tex] \theta= \dfrac{180*60}{\pi}\,D_{anisotropy}\,\big[\dfrac{c}{H_{0}\,(1+z)} \ ,\int_{z_{begin}}^{z} \,(\Omega_{m}\,(1+z)^{3}+\Omega_{\Lambda}+\Omega_{k}\,(1+z)^{2})^{-1/2}\,dz\big]^{-1} [/tex]
For [tex]\Omega_{m}=0.3[/tex] and [tex]\Omega_{\Lambda}=0.7[/tex], i get the figure in attachment. For "z=1100", I have an angle apparent equal to about 30 arcminutes for Danisotropy=380.000 ly , i.e 0.5 degrees.
This value is relatively quite far from the first peak in CMB power spectrum which is localized at "l=200", i.e "theta= 1 degree".
Does the relation between multipole "l" and "theta" angle take the form : "theta(°)=200/l" ?
If it does, then my calculation is the value of one of the others peaks which are localized for smaller angles ?
If anyone could explain me to which peak corresponds this result,
Thanks in advance.
ps: here is the Matlab code :
Code:
function dist_diam_ang
c=3*10^(5);
H0=71;
Omega2_m=0.3;
Omega2_red=0.7;
Omega2_k=1-Omega2_m-Omega2_red;
z_begin=0.1;
z_final=1100;
size1=(380000/3.26)*10^(-6); % Megaparsec
inter=1;
j=1;
for z=z_begin:inter:z_final
theta1(j)=(180*60/(pi))*size1*(c/(H0*(1+z))*(quad(@(x)myfunc(x,Omega2_m,Omega2_red,Omega2_k),z_begin,z)))^(-1);
j=j+1;
end
plot(z_begin:inter:z_final,theta1);
ylim([0 60]);
xlabel('Redshift z');
ylabel('Apparent Angle (arcminute)');
%-------------------%
function y = myfunc(x,Omega_m,Omega_red,Omega_k)
y=(Omega_m*(1+x).^(3)+Omega_red+Omega_k*(1+x).^(2)).^(-1/2);
end
%-------------------%
end