Analysis of a cryolava flow-like feature on Titan

This paper reports on the analysis of the highest spatial resolution hyperspectral images acquired by the Visual and Infrared Mapping Spectrometer (VIMS) onboard the Cassini spacecraft during its prime mission. A bright area matches a flow-like feature coming out of a caldera-like feature observed in Synthetic Aperture Radar (SAR) data recorded by the Cassini radar experiment [Lopes et al., 2007. Cryovolcanic features on Titan's surface as revealed by the Cassini Titan Radar Mapper. Icarus 186, 395-412, doi:10.1016/j.icarus.2006.09.006]. In this SAR image, the flow extends about 160 km east of the caldera. The contrast in brightness between the flow and the surroundings progressively vanishes, suggesting alteration or evolution of the composition of the cryolava during the lifetime of the eruptions. Dunes seem to cover part of this flow on its eastern end. We analyze the different terrains using the Spectral Mixing Analysis (SMA) approach of the Multiple-Endmember Linear Unmixing Model (MELSUM, Combe et al., 2008). The study area can be fully modeled by using only two types of terrains. Then, the VIMS spectra are compared with laboratory spectra of known materials in the relevant atmospheric windows (from 1 to 2.78 μm). We considered simple molecules that could be produced during cryovolcanic events, including H₂O, CO₂ (using two different grain sizes), CH₄ and NH₃. We find that the mean spectrum of the cryoflow-like feature is not consistent with pure water ice. It can be best fitted by linear combinations of spectra of the candidate materials, showing that its composition is compatible with a mixture of H₂O, CH₄ and CO_2.     

Primordial non-Gaussianities in the intergalactic medium

We present results from the first high-resolution hydrodynamical simulations of non-Gaussian cosmological models. We focus on the statistical properties of the transmitted Lyman-α flux in the high-redshift intergalactic medium. Imprints of non-Gaussianity are present and are larger at high redshifts. Differences larger than 20 per cent at   z > 3  in the flux probability distribution function for high-transmissivity regions (voids) are expected for values of the non-linearity parameter   f _NL=±100  when compared to a standard Λ cold dark matter cosmology with   f _NL= 0  . We also investigate the one-dimensional flux bispectrum: at the largest scales (corresponding to tens of Mpc), we expect deviations in the flux bispectrum up to 20 per cent at   z ∼ 4  (for   f _NL=±100  ), significantly larger than deviations of ∼3 per cent in the flux power spectrum. We briefly discuss possible systematic errors that can contaminate the signal. Although challenging, a detection of non-Gaussianities in the interesting regime of scales and redshifts probed by the Lyman-α forest could be possible with future data sets.     

A fast Bayesian approach to discrete object detection in astronomical data sets – PowellSnakes I

A new fast Bayesian approach is introduced for the detection of discrete objects immersed in a diffuse background. This new method, called PowellSnakes, speeds up traditional Bayesian techniques by (i) replacing the standard form of the likelihood for the parameters characterizing the discrete objects by an alternative exact form that is much quicker to evaluate; (ii) using a simultaneous multiple minimization code based on Powell's direction set algorithm to locate rapidly the local maxima in the posterior and (iii) deciding whether each located posterior peak corresponds to a real object by performing a Bayesian model selection using an approximate evidence value based on a local Gaussian approximation to the peak. The construction of this Gaussian approximation also provides the covariance matrix of the uncertainties in the derived parameter values for the object in question. This new approach provides a speed up in performance by a factor of '100' as compared to existing Bayesian source extraction methods that use Monte Carlo Markov chain to explore the parameter space, such as that presented by Hobson & McLachlan. The method can be implemented in either real or Fourier space. In the case of objects embedded in a homogeneous random field, working in Fourier space provides a further speed up that takes advantage of the fact that the correlation matrix of the background is circulant. We illustrate the capabilities of the method by applying to some simplified toy models. Furthermore, PowellSnakes has the advantage of consistently defining the threshold for acceptance/rejection based on priors which cannot be said of the frequentist methods. We present here the first implementation of this technique (version I). Further improvements to this implementation are currently under investigation and will be published shortly. The application of the method to realistic simulated Planck observations will be presented in a forthcoming publication.     

Strong lensing of submillimetre galaxies: a tracer of foreground structure?

The steep source counts and negative K -corrections of bright submillimetre galaxies (SMGs) suggest that a significant fraction of those observed at high flux densities may be gravitationally lensed, and that the lensing objects may often lie at redshifts above 1, where clusters of galaxies are difficult to detect through other means. In this case, follow-up of bright SMGs may be used to identify dense structures along the line-of-sight. Here, we investigate the probability for SMGs to experience strong lensing, using the latest N -body simulations and observed source flux and redshift distributions. We find that almost all high-redshift sources with a flux density above 100 mJy will be lensed if they are not relatively local galaxies. We also give estimates of the fraction of sources experiencing strong lensing as a function of observed flux density. This has implications for planning follow-up observations for bright SMGs discovered in future surveys with the Submillimetre Common-User Bolometer Array 2 and other instruments. The largest uncertainty in these calculations is the maximum allowed lensing amplification, which is dominated by the presently unknown spatial extent of SMGs.     

Clustering of luminous red galaxies – II. Small-scale redshift-space distortions

This is the second paper of a series where we study the clustering of luminous red galaxies (LRG) in the recent spectroscopic Sloan Digital Sky Survey (SDSS) data release, DR6, which has 75 000 LRG covering over  1 Gpc³  h ⁻³  for  0.15 < z < 0.47  . Here, we focus on modelling redshift-space distortions in  ξ(σ, π)  , the two-point correlation in separate line-of-sight and perpendicular directions, at small scales and in the line-of-sight. We show that a simple Kaiser model for the anisotropic two-point correlation function in redshift space, convolved with a distribution of random peculiar velocities with an exponential form, can describe well the correlation of LRG on all scales. We show that to describe with accuracy the so-called 'fingers-of-God' (FOG) elongations in the radial direction, it is necessary to model the scale dependence of both bias b and the pairwise rms peculiar velocity σ₁₂ with the distance. We show how both quantities can be inferred from the  ξ(σ, π)  data. From   r ≃ 10 Mpc  h ⁻¹  to   r ≃ 1 Mpc  h ⁻¹  , both the bias and σ₁₂ are shown to increase by a factor of 2: from   b = 2  to 4 and from  σ₁₂= 400  to  800 km s⁻¹  . The latter is in good agreement, within a 5 per cent accuracy in the recovered velocities, with direct velocity measurements in dark matter simulations with  Ω_m= 0.25  and  σ₈= 0.85  .     

The study of topology of the Universe using multipole vectors

We study a multipole vector-based decomposition of cosmic microwave background data in order to search for signatures of a multiconnected topology of the universe. Using 10⁶ simulated maps, we analyse the multipole vector distribution on the sky for the lowest order multipoles together with the probability distribution function of statistics based on the sum of the dot products of the multipole vectors for both the simply connected flat universe and universes with the topology of a 3 torus. The estimated probabilities of obtaining lower values for these statistics as compared to the 5-yr Wilkinson Microwave Anisotropy Probe data indicate that the observed alignment of the quadrupole and octopole is statistically favoured in a 3-torus topology where at least one dimension of the fundamental domain is significantly shorter than the diameter of the observable Universe, as compared to the usual standard simply connected universe. However, none of the obtained results is able to clearly rule out the latter (at more than 97 per cent confidence level). Multipole vector statistics do not appear to be very sensitive to the signatures of a 3-torus topology if the shorter dimension of the domain becomes comparable to the diameter of the observable Universe. Unfortunately, the signatures are also significantly diluted by the integrated Sachs–Wolfe effect.     

The bending of light and lensing in modified gravity

Our modified gravity theory (MOG) was used successfully in the past to explain a range of astronomical and cosmological observations, including galaxy rotation curves, the cosmic microwave background acoustic peaks, and the galaxy mass power spectrum. MOG was also used successfully to explain the unusual features of the Bullet Cluster  1E0657−558  without exotic dark matter. In the present work, we derive the relativistic equations of motion in the spherically symmetric field of a point source in MOG and, in particular, we derive equations for light bending and lensing. Our results also have broader applications in the case of extended distributions of matter, and they can be used to validate the Bullet Cluster results and provide a possible explanation for the merging clusters in Abell 520.     

Electron Densities in Quiescent Prominences Derived from Eclipse Observations

We develop a diagnostic tool for determination of the electron densities in solar prominences using eclipse data. The method is based on analysis of the hydrogen Balmer-line intensities (namely Hα and Hβ) and the white-light emission due to Thomson scattering on the prominence electrons. Our approach represents a generalization of the ratio method already used by Koutchmy, Lebecq, and Stellmacher (Astron. Astrophys. 119, 261, 1983). In this paper we use an extended grid of non-LTE prominence models of Gouttebroze, Heinzel, and Vial (Astron. Astrophys. Suppl. Ser. 99, 513, 1993) and derive various useful relations between prominence radiation properties and electron densities. Simultaneously, an effective geometrical thickness of the prominence can also be obtained. As an example we apply our general technique to original eclipse data of Koutchmy, Lebecq, and Stellmacher (Astron. Astrophys. 119, 261, 1983). Finally, we use our results to determine the color of prominences as it should be seen during total eclipses.