Context. The interpretation of polarised radiation emerging from a planetary atmosphere must rely on solutions to the vector Radiative Transport Equation (vRTE). Monte Carlo integration of the vRTE is a valuable approach for its flexible treatment of complex viewing and/or illumination geometries and because it can intuitively incorporate elaborate physics. Aims. We present a novel Pre-Conditioned Backward Monte Carlo (PBMC) algorithm for solving the vRTE and apply it to planetary atmospheres irradiated from above. As classical BMC methods, our PBMC algorithm builds the solution by simulating the photon trajectories from the detector towards the radiation source, i.e. in the reverse order of the actual photon displacements. Methods. We show that the neglect of polarisation in the sampling of photon propagation directions in classical BMC algorithms leads to unstable and biased solutions for conservative, optically-thick, strongly-polarising media such as Rayleigh atmospheres. The numerical difficulty is avoided by pre-conditioning the scattering matrix with information from the scattering matrices of prior (in the BMC integration order) photon collisions. Pre-conditioning introduces a sense of history in the photon polarisation states through the simulated trajectories. Results. The PBMC algorithm is robust and its accuracy is extensively demonstrated via comparisons with examples drawn from the literature for scattering in diverse media. Since the convergence rate for MC integration is independent of the integral's dimension, the scheme is a valuable option for estimating the disk-integrated signal of stellar radiation reflected from planets. Such a tool is relevant in the prospective investigation of exoplanetary phase curves. We lay out two frameworks for disk integration and, as an application, explore the impact of atmospheric stratification on planetary phase curves...

We cross-correlate the third-year WMAP data with galaxy samples extracted from the SDSS DR5 (SDSS5) covering 16% of the sky. These measurements confirm a positive cross-correlation, which is well fitted by the integrated Sachs-Wolfe (ISW) effect for flat LCDM models with a cosmological constant.