I'm having difficulties in calculating the angular size of the causally connected regions on the cosmic microwave background (CMB), as seen from Earth today. I read in several documents that this angle is of about $1^{\circ}$, but most authors are giving only crude hand waving arguments about that number. See for example that page (look at the last two paragraphs):
https://ned.ipac.caltech.edu/level5/Sept02/Kinney/Kinney4_2.html
I'm trying to reproduce that value by explicit calculations from the standard FLRW metric, in the case of a spatially flat geometry ($k = 0$): \begin{equation}\tag{1} ds^2 =dt^2 - a^2(t) (dx^2 + dy^2 + dz^2). \end{equation} At observation time $t_{obs}$ (today: $t_{obs} \approx 13,8~\mathrm{Gyears}$), the proper distance from a given source (emitting light at time $t_{em} \approx 300~000~\mathrm{years}$) is given by $ds^2 = 0$ (light-like spacetime intervall): \begin{equation}\tag{2} \mathcal{D}(t_{obs}, t_{em}) = a(t_{obs}) \int_{t_{em}}^{t_{obs}} \frac{1}{a(t)} \, dt. \end{equation} For example, this distance to the CMB surface, in the case of a dust universe, is found with the scale factor $a(t) \propto t^{2/3}$: \begin{equation}\tag{3} \mathcal{D} = 3 \, (\, t_{obs} - t_{obs}^{2/3} \, t_{em}^{1/3}). \end{equation} This gives $\mathcal{D} \approx 40,2~\mathrm{Gly}$ (more accurate models with radiation gives about $42$ or $45~\mathrm{Gly}$).
Now, the causally correlated regions on the CMB sphere should have a proper radius of (considering the dust only universe): \begin{equation}\tag{4} R_{causal} = a(t_{em}) \int_0^{t_{em}} \frac{1}{a(t)} \, dt = 3 \, t_{em}, \end{equation} i.e. $R_{causal} \approx 9 \times 10^5 ~ \mathrm{ly}$. As seen from Earth, the angular size of a causal patch should have an angular size $\alpha_{causal}$ of: \begin{equation}\tag{5} \alpha_{causal} = 2 \arctan{\Big( \frac{R_{causal}}{\mathcal{D}} \Big)} \approx 0.003^{\circ}. \end{equation} Of course, this is much too short, and I'm probably doing a naive calculation. I don't know where I'm making a mistake.
How should I fix the angular size (5)?
EDIT: Apparently, the right formula fixing (5) is the following (the factor 2 is to get the full angular diameter, and not just the angular radius of the causal patch): \begin{equation}\tag{6} \alpha_{causal} = 2 \arctan{\Big( \frac{\displaystyle{\int_{0}^{t_{em}} \frac{1}{a(t)} \, dt}}{\displaystyle{\int_{t_{em}}^{t_{ob}} \frac{1}{a(t)} \, dt}} \Big)}, \end{equation} but I don't understand why the angle is found by the ratio of the comoving lenghts instead of the proper lenghts.