Source author record

Boudewijn F. Roukema

Boudewijn F. Roukema appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

42works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

42 published item(s)

preprint2022arXiv

Toward Long-Term and Archivable Reproducibility

Analysis pipelines commonly use high-level technologies that are popular when created, but are unlikely to be readable, executable, or sustainable in the long term. A set of criteria is introduced to address this problem: Completeness (no execution requirement beyond a minimal Unix-like operating system, no administrator privileges, no network connection, and storage primarily in plain text); modular design; minimal complexity; scalability; verifiable inputs and outputs; version control; linking analysis with narrative; and free and open source software. As a proof of concept, we introduce "Maneage" (Managing data lineage), enabling cheap archiving, provenance extraction, and peer verification that has been tested in several research publications. We show that longevity is a realistic requirement that does not sacrifice immediate or short-term reproducibility. The caveats (with proposed solutions) are then discussed and we conclude with the benefits for the various stakeholders. This article is itself a Maneage'd project (project commit 54e4eb2).

preprint2020arXiv

[not Rp] Reproducibility of 'Poincare dodecahedral space parameter estimates'

Is a scientific research paper based on (i) public, online observational data files and (ii) providing free-licensed software for reproducing its results easy to reproduce by the same author a decade later? This paper attempts to reproduce a cosmic topology observational paper published in 2008 and satisfying both criteria (i) and (ii). The reproduction steps are defined formally in a free-licensed git repository package "0807.4260" and qualitatively in the current paper. It was found that the effort in upgrading the Fortran 77 code at the heart of the software, interfaced with a C front end, and originally compiled with g77, in the content of the contemporary gfortran compiler, risked being too great to be justified on any short time scale. In this sense, the results of RBG08 are not as reproducible as they appeared to be, despite both (i) data availability and (ii) free-licensing and public availability of the software. The software and a script to reproduce the steps of this incomplete reproduction are combined in a new git repository named 0807.4260, following the ArXiv identity code (arXiv:0807.4260) of RBG08.

preprint2020arXiv

Cosmological signatures of torsion and how to distinguish torsion from the dark sector

Torsion is a non-Riemannian geometrical extension of general relativity that allows including the spin of matter and the twisting of spacetime. Cosmological models with torsion have been considered in the literature to solve problems of either the very early (high redshift $z$) or the present-day Universe. This paper focuses on distinguishable observational signatures of torsion that could not be otherwise explained with a scalar field in pseudo-Riemannian geometry. We show that when torsion is present, the cosmic duality relation between the angular diameter distance, $D_{\mathrm A}$, and the luminosity distance, $D_{\mathrm L}$, is broken. We show how the deviation described by the parameter $η= D_{\mathrm L}/[ D_{\mathrm A}(1+z)^2] -1 $ is linked to torsion and how different forms of torsion lead to special-case parametrisations of $η$, including $η_0 z$, $η_0 z/(1+z)$, and $η_0 \ln (1+z)$. We also show that the effects of torsion could be visible in low-redshift data, inducing biases in supernovae-based $H_0$ measurements. We also show that torsion can impact the Clarkson-Bassett-Lu (CBL) function ${\cal C}(z) = 1 + H^2 (D D'' - D'^2) + H H' D D'$, where $D$ is the transverse comoving distance. If $D$ is inferred from the luminosity distance, then, in general non-zero torsion models, ${\cal C}(z) \ne 0$. For pseudo-Riemannian geometry, the Friedmann-Lemaitre-Robertson-Walker (FLRW) metric has ${\cal C}(z) \equiv 0$; thus, measurement of the CBL function could provide another diagnostic of torsion.

preprint2016arXiv

The comptonization parameter from simulations of single-frequency, single-dish, dual-beam, cm-wave observations of galaxy clusters and mitigating CMB confusion using the Planck sky survey

Systematic effects in dual-beam, differential, radio observations of extended objects are discussed in the context of the One Centimeter Receiver Array (OCRA). We use simulated samples of Sunyaev-Zel'dovich (SZ) galaxy clusters at low ($z<0.4$) and intermediate ($0.4<z<1.0$) redshifts to study the implications of operating at a single frequency (30 GHz) on the accuracy of extracting SZ flux densities and of reconstructing comptonization parameters with OCRA. We analyze dependences on cluster mass, redshift, observation strategy, and telescope pointing accuracy. Using $Planck$ data to make primary cosmic microwave background (CMB) templates, we test the feasibility of mitigating CMB confusion effects in observations of SZ profiles at angular scales larger than the separation of the receiver beams.

preprint2015arXiv

Evidence for an environment-dependent shift in the baryon acoustic oscillation peak

The Friedmann-Lemaitre-Robertson-Walker (FLRW) metric assumes comoving spatial rigidity of metrical properties. The curvature term in comoving coordinates is environment-independent and cannot evolve. In the standard model, structure formation is interpreted accordingly: structures average out on the chosen metrical background, which remains rigid in comoving coordinates despite nonlinear structure growth. The latter claim needs to be tested, since it is a hypothesis that is not derived using general relativity. We introduce a test of the comoving rigidity assumption by measuring the two-point auto-correlation function on comoving scales---assuming FLRW comoving spatial rigidity---in order to detect shifts in the baryon acoustic oscillation (BAO) peak location for Large Red Galaxy (LRG) pairs of the Sloan Digital Sky Survey Data Release 7. In tangential directions, subsets of pairs overlapping with superclusters or voids show the BAO peak. The tangential BAO peak location for overlap with Nadathur & Hotchkiss superclusters is $4.3\pm1.6$ Mpc/h less than that for LRG pairs unselected for supercluster overlap, and $6.6\pm2.8$ Mpc/h less than that of the complementary pairs. Liivamagi et al. superclusters give corresponding differences of $3.7\pm2.9$ Mpc/h and $6.3\pm2.6$ Mpc/h, respectively. We have found moderately significant evidence (Kolmogorov--Smirnov tests suggest very significant evidence) that the BAO peak location for supercluster-overlapping pairs is compressed by about 6% compared to that of the complementary sample, providing a potential challenge to FLRW models and a benchmark for predictions from models based on an averaging approach that leaves the spatial metric a priori unspecified.

preprint2015arXiv

Is the baryon acoustic oscillation peak a cosmological standard ruler?

In the standard model of cosmology, the Universe is static in comoving coordinates; expansion occurs homogeneously and is represented by a global scale factor. The baryon acoustic oscillation (BAO) peak location is a statistical tracer that represents, in the standard model, a fixed comoving-length standard ruler. Recent gravitational collapse should modify the metric, rendering the effective scale factor, and thus the BAO standard ruler, spatially inhomogeneous. Using the Sloan Digital Sky Survey, we show to high significance (P < 0.001) that the spatial compression of the BAO peak location increases as the spatial paths' overlap with superclusters increases. Detailed observational and theoretical calibration of this BAO peak location environment dependence will be needed when interpreting the next decade's cosmological surveys.

preprint2014arXiv

Simplicity in cosmology: add virialisation, remove $Λ$, keep classical GR

Present-day extragalactic observations are mostly rather well-modelled by a general-relativistic model, the $Λ$ CDM model. The model appears to surpass the limits of known physics by requiring that the Universe be dominated by "dark energy". However, the model sacrifices physical simplicity in favour of applied mathematical simplicity. A physically simpler, general-relativistic alternative to the $Λ$ CDM model is described here, along with preliminary observational checks. Thus, it will be argued that extragalactic observations such as the distance-modulus--redshift relation of type Ia supernovae are well-modelled within classical general relativity, without the addition of "new physics".

preprint2013arXiv

A first-digit anomaly in the 2009 Iranian presidential election

A local bootstrap method is proposed for the analysis of electoral vote-count first-digit frequencies, complementing the Benford's Law limit. The method is calibrated on five presidential-election first rounds (2002--2006) and applied to the 2009 Iranian presidential-election first round. Candidate K has a highly significant (p< 0.15%) excess of vote counts starting with the digit 7. This leads to other anomalies, two of which are individually significant at p\sim 0.1%, and one at p\sim 1%. Independently, Iranian pre-election opinion polls significantly reject the official results unless the five polls favouring candidate A are considered alone. If the latter represent normalised data and a linear, least-squares, equal-weighted fit is used, then either candidates R and K suffered a sudden, dramatic (70%\pm 15%) loss of electoral support just prior to the election, or the official results are rejected (p\sim 0.01%).

preprint2013arXiv

Dark energy with rigid voids versus relativistic voids alone

The standard model of cosmology is dominated - at the present epoch - by dark energy. Its voids are rigid and Newtonian within a relativistic background. The model prevents them from becoming hyperbolic. Observations of rapid velocity flows out of voids are normally interpreted within the standard model that is rigid in comoving coordinates, instead of allowing the voids' density parameter to drop below critical and their curvature to become negative. Isn't it time to advance beyond nineteenth century physics and relegate dark energy back to the "no significant evidence" box?

preprint2013arXiv

Deep redshift topological lensing: strategies for the T^3 candidate

The 3-torus (T^3) FLRW model better fits the nearly zero large-scale auto-correlation of the WMAP CMB sky maps than the infinite flat model. The T^3 model's WMAP parameters imply approximately equal-redshift topological lensing at z \sim 6. We investigate observational strategies for rejecting the T^3 solution or providing candidate topologically lensed galaxy pairs. T^3 holonomies are applied to (i) existing z \sim 6 observations and (ii) simulated observations, creating multiply connected catalogues. Corresponding simply connected catalogues are generated. Each catalogue is analysed using a successive filter method and collecting matched quadruples. Quadruple statistics between the multiply and simply connected catalogues are compared. The expected rejection of the hypothesis, or detection of candidate topologically lensed galaxies, is possible at a significance of 5% for a pair of T^3 axis-centred northern and southern surveys if photometric redshift accuracy is σ(\zphot) < 0.01 for a pair of nearly complete 100 deg^2 surveys with a total of > 500 galaxies over 4.3 < z < 6.6, or for a pair of 196 deg^2 surveys with > 400 galaxies and σ(\zphot) < 0.02 over 4<z<7. Dropping the maximum time interval in a pair from Δt =1 Gyr/h to Δt =0.1 Gyr/h yields a requirement of σ(\zphot) < 0.005 or σ(\zphot) < 0.01, respectively. Millions of z \sim 6 galaxies will be observed over fields of these sizes during the coming decades, implying much stronger constraints. The question is not if the hypothesis will be rejected or confirmed, it is when.

preprint2013arXiv

Topological implications of inhomogeneity

The approximate homogeneity of spatial sections of the Universe is well supported observationally, but the inhomogeneity of the spatial sections is even better supported. Here, we consider the implications of inhomogeneity in dust models for the connectedness of spatial sections at early times. We consider a non-global Lemaitre-Tolman-Bondi (LTB) model designed to match observations, a more general, heuristic model motivated by the former, and two specific, global LTB models. We propose that the generic class of solutions of the Einstein equations projected back in time from the spatial section at the present epoch includes subclasses in which the spatial section evolves (with increasing time) smoothly (i) from being disconnected to being connected, or (ii) from being simply connected to being multiply connected, where the coordinate system is comoving and synchronous. We show that (i) and (ii) each contain at least one exact solution. These subclasses exist because the Einstein equations allow non-simultaneous big bang times. The two types of topology evolution occur over time slices that include significantly post-quantum epochs if the bang time varies by much more than a Planck time. In this sense, it is possible for cosmic topology evolution to be "mostly" classical.

preprint2013arXiv

Virialisation-induced curvature as a physical explanation for dark energy

The geometry of the dark energy and cold dark matter dominated cosmological model (LambdaCDM) is commonly assumed to be given by a Friedmann-Lemaitre-Robertson-Walker (FLRW) metric, i.e. it assumes homogeneity in the comoving spatial section. The homogeneity assumption fails most strongly at (i) small distance scales and (ii) recent epochs, implying that the FLRW approximation is most likely to fail at these scales. We use the virialisation fraction to quantify (i) and (ii), which approximately coincide with each other on the observational past light cone. For increasing time, the virialisation fraction increases above 10% at about the same redshift (sim 1-3) at which Omega_Lambda grows above 10% (approx 1.8). Thus, instead of non-zero Omega_Lambda, we propose an approximate, general-relativistic correction to the matter-dominated (Omega_m =1, Omega_Lambda=0), homogeneous metric (Einstein de Sitter, EdS). A low-redshift effective matter-density parameter of Omega_m^eff(0) = 0.26 pm 0.05 is inferred. Over redshifts 0 < z < 3, the distance modulus of the virialisation-corrected EdS model approximately matches the LambdaCDM distance modulus. This rough approximation assumes "old physics" (general relativity), not "new physics". Thus, pending more detailed calculations, we strengthen the claim that "dark energy" should be considered as an artefact of emerging average curvature in the void-dominated Universe, via a novel approach that quantifies the relation between virialisation and average curvature evolution.

preprint2013arXiv

Virialization-induced curvature versus dark energy

The concordance model is successful in explaining numerous observable phenomena at the price of introducing an exotic source of unknown origin: dark energy. Dark energy dominance occurs at recent epochs, when we expect most cosmological structures to have already formed, and thus, when the error induced by forcing the homogeneous FLRW metric onto the data is expected to be the most significant. We propose a way to quantify the impact of deviations from homogeneity on the evolution of cosmological parameters. Using a multi-scale partitioning approach and the virialization fraction estimated from numerical simulations in an Einstein-de Sitter model, we obtain an observationally realistic distance modulus over redshifts 0 < z < 3 by a relativistic correction of the FLRW metric.

preprint2012arXiv

A new look at the long-period eclipsing binary V383 Sco

V383Sco was discovered to be an eclipsing binary at the beginning of the XX century. This system has one of the longest orbital periods known (13.5yr) and was initially classified as a zet_Aur-type variable. It was then forgotten for decades. This study provides a detailed look at the V383Sco, using new data obtained around the last eclipse in 2007/8. There was a suspicion that this system could be similar to eclipsing systems with extensive dusty disks like EECep and eps_Aur. This and other, alternative hypotheses are considered. The ASAS-3 VI light curves have been used to examine photometric changes. Low-(LRS) and high-res.(HRS) spectra have been used for spectral classification, to analyse line profiles, as well as to determine the reddening, radial velocities (RVs) and distance. The SED was analysed. Using original numerical code, we performed a simplified model of the eclipse, taking into account the pulsations of one of the components. The LRS shows traces of molecular bands, characteristic of an M-type supergiant. The presence of this star in the system is confirmed by SED, by a strong dependence of the eclipse depth on the photometric bands, and by pulsational changes. The presence of a low excitation nebula around the system has been inferred from [OI] 6300A emission. Analysis of the RVs, reddening, and P-L relation for Mira-type stars imply a distance to the V383Sco of 8.4+-0.6 kpc. The distance to the nearby V381Sco is 6.4+-0.8 kpc. The very different and oppositely directed RVs of these systems (89.8 vs -178.8 km/s) seem to be in agreement with a bulge/bar kinematic model of the Galactic centre and inconsistent with purely circular motion. We have found evidence for the presence of a pulsating M-type supergiant in the V383Sco which periodically obscures the much more luminous F0I-type star, causing the deep (possibly total) eclipses which vary in duration and shape.

preprint2012arXiv

A relativistic model of the topological acceleration effect

It has previously been shown heuristically that the topology of the Universe affects gravity, in the sense that a test particle near a massive object in a multiply connected universe is subject to a topologically induced acceleration that opposes the local attraction to the massive object. This effect distinguishes different comoving 3-manifolds, potentially providing a theoretical justification for the Poincaré dodecahedral space observational hypothesis and a dynamical test for cosmic topology. It is necessary to check if this effect occurs in a fully relativistic solution of the Einstein equations that has a multiply connected spatial section. A Schwarzschild-like exact solution that is multiply connected in one spatial direction is checked for analytical and numerical consistency with the heuristic result. The T$^1$ (slab space) heuristic result is found to be relativistically correct. For a fundamental domain size of $L$, a slow-moving, negligible-mass test particle lying at distance $x$ along the axis from the object of mass $M$ to its nearest multiple image, where $GM/c^2 \ll x \ll L/2$, has a residual acceleration away from the massive object of $4ζ(3) G(M/L^3)\,x$, where $ζ(3)$ is Apéry's constant. For $M \sim 10^14 M_\odot$ and $L \sim 10$ to $20\hGpc$, this linear expression is accurate to $\pm10%$ over $3\hMpc \ltapprox x \ltapprox 2\hGpc$. Thus, at least in a simple example of a multiply connected universe, the topological acceleration effect is not an artefact of Newtonian-like reasoning, and its linear derivation is accurate over about three orders of magnitude in $x$.

preprint2012arXiv

Cosmic topology affects dynamics

The role of global topology in the dynamics of the Universe is poorly understood. Along with observational programmes for determining the topology of the Universe, some small theoretical steps have recently been made. Heuristic Newtonian-like arguments suggest a topological acceleration effect that differs for differing spatial sections. A relativistic spacetime solution solution shows that the effect is not just a Newtonian artefact.

preprint2012arXiv

Topological acceleration in relativistic cosmology

Heuristic approaches in cosmology bypass more difficult calculations that would more strictly agree with the standard Einstein equation. These give us the well-known Friedmann-Lemaitre-Robertson-Walker (FLRW) models, and, more recently, the feedback effect of the global topology of spatial sections on the acceleration of test particles. Forcing the FLRW heuristic model on observations leads to dark energy, which, pending fully relativistic calculations, is best interpreted as an artefact. Could topological acceleration also be an artefact of using a heuristic approach? A multiply connected exact solution of the Einstein equation shows that topological acceleration is present in at least one fully relativistic case---it is not an artefact of Newtonian-like thinking.

preprint2011arXiv

The size of the Universe according to the Poincare dodecahedral space hypothesis

One of the FLRW models that best fits the WMAP sky maps of the CMB is the Poincare dodecahedral space. The optimal fit of this model to WMAP data was recently found using an optimal cross-correlation method, but the systematic error in the estimate of the matched-circle angular radius α, or equivalently, the (comoving) size of the Universe 2\rinj (twice the injectivity radius), might be much higher than the random error. In order to increase the falsifiability of the model, it would be useful to reduce the uncertainty in this estimate and to estimate the fraction of the sky where multiply imaged gravitationally bound objects should potentially be detectable. "Matched discs" are defined in order to describe a useful subset of multiply imaged objects. The cross-correlation method at \ltapprox 1 \hGpc is applied to WMAP 7-year data in order to improve the estimate of α. The improved matched-circle radius estimate is α= 23 \pm 1.4 deg, where the uncertainty represents systematic error dependent on the choices of galactic mask and all-sky map. This is equivalent to 2\rinj = 18.2\pm 0.5 \hGpc for matter density parameter Omega_m=0.28\pm 0.02. The lowest redshift of multiply imaged objects is z=106\pm18. Multiply imaged high overdensity (rare) peaks visible during 200>z>106 should be present in matched discs of radius 14.8\pm2.3 deg. The accuracy in the matched circle radius estimate is considerably improved by using the higher resolution signal. The predicted matched discs (over 200>z>106) project to about 20% of the full sky. Since any object located exactly in the discs would be multiply imaged at equal redshifts, evolutionary effects would be small for objects that are nearly located in the discs.

preprint2010arXiv

A measure on the set of compact Friedmann-Lemaitre-Robertson-Walker models

Compact, flat Friedmann-Lemaitre-Robertson-Walker (FLRW) models have recently regained interest as a good fit to the observed cosmic microwave background temperature fluctuations. However, it is generally thought that a globally, exactly-flat FLRW model is theoretically improbable. Here, in order to obtain a probability space on the set F of compact, comoving, 3-spatial sections of FLRW models, a physically motivated hypothesis is proposed, using the density parameter Omega as a derived rather than fundamental parameter. We assume that the processes that select the 3-manifold also select a global mass-energy and a Hubble parameter. The inferred range in Omega consists of a single real value for any 3-manifold. Thus, the obvious measure over F is the discrete measure. Hence, if the global mass-energy and Hubble parameter are a function of 3-manifold choice among compact FLRW models, then probability spaces parametrised by Omega do not, in general, give a zero probability of a flat model. Alternatively, parametrisation by the injectivity radius r_inj ("size") suggests the Lebesgue measure. In this case, the probability space over the injectivity radius implies that flat models occur almost surely (a.s.), in the sense of probability theory, and non-flat models a.s. do not occur.

preprint2010arXiv

Galactic Plane image sharpness as a check on cosmic microwave background mapmaking

The largest uncollapsed inhomogeneity in the observable Universe is statistically represented in the quadrupole signal of the cosmic microwave background (CMB) sky maps as observed by the Wilkinson Microwave Anisotropy Probe (WMAP). The constant temporal offset of -25.6 ms between the timestamps of the spacecraft attitude and observational data records in the time-ordered data (TOD) of the WMAP observations was suspected to imply that previously derived all-sky CMB maps are erroneous, and that the quadrupole is in large part an artefact. The optimal focussing of bright objects in the Galactic Plane plays a key role in showing that no error occurred at the step of mapmaking from the calibrated TOD. Instead, the error had an effect when the uncalibrated TOD were calibrated. Estimates of the high-latitude quadrupole based on the wrongly calibrated WMAP maps are overestimated by about 15-60%.

preprint2010arXiv

On the suspected timing error in WMAP map-making

A large fraction of the previously estimated WMAP CMB quadrupole signal would be an artefact of incorrect Doppler dipole subtraction if the hypothesis of a small timing interpolation error were correct. Observations of bright foreground objects constitute part of the time-ordered-data (TOD). Scans of an object in different directions should be shifted by the would-be timing error, causing a blurring effect. Three half-years of the calibrated, filtered WMAP TOD are compiled individually for the four W band differencing assemblies (DA's), with no masking of bright objects, giving 12 maps for each timing offset. Percentiles of the temperature-fluctuation distribution in each map at HEALPix resolution N_side=2048 are used to determine the dependence of all-sky image sharpness on the timing offset. In the W band, the 99.999% percentile, i.e. the temperature fluctuation in the approx 503-rd brightest pixel, is the least noisy percentile. Using this statistic, the hypothesis that a -25.6 ms offset relative to the timing adopted by the WMAP collaboration gives a focus at least as sharp as the uncorrected timing is rejected at 4.6σsignificance. The Q and V band maps also reject the -25.6 ms offset hypothesis at high statistical significance. The requirement that the correct choice of timing offset must maximise image sharpness implies that the hypothesis of a timing error in the WMAP collaboration's compilation of the WMAP calibrated, filtered TOD is rejected at high statistical significance in each of the Q, V and W wavebands. However, the hypothesis that a timing error was applied during calibration of the raw TOD, leading to a dipole-induced difference signal, is not excluded by this method.

preprint2010arXiv

On the suspected timing-offset-induced calibration error in the Wilkinson microwave anisotropy probe time-ordered data

In the time-ordered data (TOD) files of the WMAP CMB observations, there is an undocumented timing offset of -25.6 ms between the spacecraft attitude and radio flux density timestamps. If the offset induced an error during calibration of the raw TOD, then estimates of the WMAP CMB quadrupole might be substantially in error. A timing error during calibration would not only induce an artificial quadrupole-like signal in the mean sky map, it would also add variance per pixel. This variance would be present in the calibrated TOD. Low-resolution map-making as a function of timing offset should show a minimum variance for the correct timing offset. Three years of the calibrated, filtered WMAP 3-year TOD are compiled into sky maps at HEALPix resolution N_side=8, individually for each of the K, Ka, Q, V and W band differencing assemblies (DA's), as a function of timing offset. The median per map of the temperature fluctuation variance per pixel is calculated and minimised against timing offset. Minima are clearly present. The timing offsets that minimise the median variance are -38 \pm 8 ms (K, Ka), -27 \pm 3 ms (Q), -43 \pm 8 ms (V), and -47 \pm 194 ms (W), i.e. an average of -30 \pm 3 ms, where the WMAP collaboration's preferred offset is 0 \pm 1.7 ms. A non-parametric bootstrap analysis rejects the latter at a significance of 99.999%. The hypothesis of a -25.6 ms offset, suggested by Liu, Xiong & Li from the TOD file timing offset, is consistent with these minima. It is difficult to avoid the conclusion that the WMAP calibrated TOD and inferred maps are wrongly calibrated. CMB quadrupole estimates (3/pi)C_2 based on the incorrectly calibrated TOD are overestimated by roughly 64 \pm 6% (KQ85 mask) to 94 \pm 10% (KQ75 mask). Ideally, the WMAP map-making pipelines should be redone starting from the uncalibrated TOD and using the -25.6 ms timing offset correction.

preprint2010arXiv

Which FLRW comoving 3-manifold is preferred observationally and theoretically?

The lack of structure greater than 10 Gpc/h in Wilkinson Microwave Anisotropy Probe (WMAP) observations of the cosmic microwave background (CMB) favours compact Friedmann-Lemaitre-Robertson-Walker (FLRW) models of the Universe. The present best candidates based on observations are the Poincare dodecahedral space S^3/I^* and the 3-torus T^3. The residual gravity effect favours the Poincare space, while a measure space argument where the density parameter is a derived parameter favours flat spaces almost surely.

preprint2009arXiv

A residual acceleration effect due to an inhomogeneity

The perturbed Friedmann-Lemaitre-Robertson-Walker models allow many different possibilities for the 3-manifold of the comoving spatial section of the Universe. It used to be thought that global properties of the spatial section, including the topology of the space, have no feedback effect on dynamics. However, an elementary, weak-limit calculation shows that in the presence of a density perturbation, a gravitational feedback effect that is algebraically similar to dark energy does exist. Moreover, the effect differs between different 3-spaces. Among the well-proportioned spaces, the effect disappears down to third order in several cases, and down to fifth order for the Poincare dodecahedral space S^3/I^*. The Poincare space, that which also is preferred in many observational analyses, is better-balanced than the other spaces.

preprint2009arXiv

Does gravity prefer the Poincare dodecahedral space?

The missing fluctuations problem in cosmic microwave background observations is naturally explained by well-proportioned small universe models. Among the well-proportioned models, the Poincare dodecahedral space is empirically favoured. Does gravity favour this space? The residual gravity effect is the residual acceleration induced by weak limit gravity from multiple topological images of a massive object on a nearby negligible mass test object. At the present epoch, the residual gravity effect is about a million times weaker in three of the well-proportioned spaces than in ill-proportioned spaces. However, in the Poincare space, the effect is 10,000 times weaker still, i.e. the Poincare space is about 10^{10} times "better balanced" than ill-proportioned spaces. Both observations and weak limit dynamics select the Poincare space to be special.

preprint2009arXiv

Some spaces are more equal than others

It has generally been thought that in perturbed Friedmann-Lemaitre-Robertson-Walker models of the Universe, global topology should not have any feedback effects on dynamics. However, a weak-field limit heuristical argument, assuming a finite particle horizon for the transmission of gravitational signals, shows that a residual acceleration effect can occur. The nature of this effect differs algebraically between different constant curvature 3-manifolds. This potentially provides a selection mechanism for the 3-manifold of comoving space.

preprint2009arXiv

The residual gravity acceleration effect in the Poincare dodecahedral space

In a flat space, the global topology of comoving space can induce a weak acceleration effect similar to dark energy. Does a similar effect occur in the case of the Poincare dodecahedral space S^3/I^*? Does the effect distinguish the Poincare space from other well-proportioned spaces? The residual acceleration effect in the Poincare space is studied here using a massive particle and a nearby test particle of negligible mass, in S^3 embedded in R^4. The weak limit gravitational attraction on a test particle at distance r is set \propto [r_C \sin(r/r_C)]^{-2}, where r_C = curvature radius, in order to satisfy Stokes' theorem. A finite particle horizon large enough to include the adjacent topological images of the massive particle is assumed. The regular, flat, 3-torus T^3 is re-examined, and two other well-proportioned spaces, S^3/T^* and S^3/O^*, are also studied. The residual gravity effect occurs in all four cases. In a perfectly regular 3-torus of side length L_a, and in S^3/T^* and S^3/O^*, the highest order term in the residual acceleration is the third order term in the Taylor expansion of r/L_a (3-torus), or r/r_C, respectively. However, the Poincare dodecahedral space is unique among the four spaces. The third order cancels, leaving the fifth order term \sim \pm 300 (r/r_C)^5 as the most significant. Not only are three of the four perfectly regular well-proportioned spaces better balanced than most other multiply connected spaces in terms of the residual acceleration effect by a factor of about a million (setting r/L_a = r/r_C \sim 10^{-3}), but the fourth of these spaces is about 10^4 times better balanced than the other three. This is the Poincare dodecahedral space. Is this unique dynamical property of the Poincare space a clue towards a theory of cosmic topology?

preprint2009arXiv

There was movement that was stationary, for the four-velocity had passed around

Is the Doppler interpretation of galaxy redshifts in a Friedmann-Lemaitre-Robertson-Walker (FLRW) model valid in the context of the approach to comoving spatial sections pioneered by de Sitter, Friedmann, Lemaitre and Robertson, i.e. according to which the 3-manifold of comoving space is characterised by both its curvature and topology? Holonomy transformations for flat, spherical and hyperbolic FLRW spatial sections are proposed. By quotienting a simply-connected FLRW spatial section by an appropriate group of holonomy transformations, the Doppler interpretation in a non-expanding Minkowski space-time, obtained via four-velocity parallel transport along a photon path, is found to imply that an inertial observer is receding from herself at a speed greater than zero, implying contradictory world-lines. The contradiction in the multiply-connected case occurs for arbitrary redshifts in the flat and spherical cases, and for certain large redshifts in the hyperbolic case. The link between the Doppler interpretation of redshifts and cosmic topology can be understood physically as the link between parallel transport along a photon path and the fact that the comoving spatial geodesic corresponding to a photon's path can be a closed loop in an FLRW model of any curvature. Closed comoving spatial loops are fundamental to cosmic topology.

preprint2006arXiv

WMAP First Year Sky Map: Hints of Poincare Dodecahedral Topology

Several analyses of the cosmic microwave background map made by the satellite WMAP suggest that the global shape of a spatial section of the Universe is that of a Poincare dodecahedral space. A summary of some of these analyses and a description of independent tests which should be able to either increase confidence in the hypothesis or else refute it to extremely high significance will be presented.

preprint2005arXiv

A constraint on any topological lensing hypothesis in the spherical case: it must be a root of the identity

Three-dimensional catalogues of objects at cosmological distances can potentially yield candidate topologically lensed pairs of sets of objects, which would be a sign of the global topology of the Universe. In the spherical case, a necessary condition, which does not exist for either null or negative curvature, can be used to falsify such hypotheses, without needing to loop through a list of individual spherical 3-manifolds. This condition is that the isometry between the two sets of objects must be a root of the identity isometry in the covering space S^3. This enables numerical falsification of topological lensing hypotheses without needing to assume any particular spherical 3-manifold. By embedding S^3 in R^4, this condition can be expressed as the requirement that M^n = I for an integer n, where M is the matrix representation of the hypothesised lensing isometry and I is the identity. Moreover, this test becomes even simpler with the requirement that the two rotation angles, theta, phi, corresponding to the given isometry, satisfy 2π/ θ, 2π/ ϕ\in Z. The calculation of this test involves finding the two eigenplanes of the matrix M. A GNU General Public Licence numerical package, called eigenplane, is made available at http://cosmo.torun.pl/GPLdownload/eigen/ for finding the rotation angles and eigenplanes of an arbitrary isometry M of S^3.

preprint2005arXiv

Cosmic crystallography using short-lived objects - active galactic nuclei

Cosmic crystallography is based on the principle that peaks in the pair separation histogram (PSH) of objects in a catalogue should be induced by the high number of topologically lensed pairs that are separated by Clifford translations, in excess to ``random'' pairs of objects. Here we present modifications of this method that successively improve the signal-to-noise ratio by removing a large part of the noise and then false signals induced by selection effects. Given the transient nature of the most readily available tracer objects, active galactic nuclei (AGNs), the former is possible because a natural filter for removing many of the noise pairs is available: when counting pairs of objects in order to create PSHs, only those with nearly identical redshifts need to be counted. This redshift filter (a maximum value of $Δz/z = 0.005$) was applied to a compilation of AGN catalogues. Further noise was removed by applying a second filter, a maximum angle $Δθ=0.075$ rad, and a minimum number of pairs $\protect\npairs=3$ to find each ``bunch of pairs'' (BoP) where the {\em vectors} (in Euclidean comoving space) defined by pairs are required to be nearly equal, whereas in the PSH only the {\em lengths} must be nearly equal. These filters reveal significant signals, which, however, are due to selection effects. A third filter, a minimum length $\Lselec=150$ {\hMpc} between the (parallel) vectors in a BoP, is found to effectively remove these selection effect pairs. After application of these successive filters, no significant topological signal was found.

preprint2004arXiv

A Hint of Poincaré Dodecahedral Topology in the WMAP First Year Sky Map

Luminet et al. (2003) suggested that WMAP data are better matched by a Poincaré dodecahedral FLRW model of global geometry, rather than by an infinite flat model. The analysis by Cornish et al. (2003) for angular radii 25-90 degrees failed to support this. Here, a matched circles analysis specifically designed to detect dodecahedral patterns of matched circles is performed over angular radii in the range 1-40 degrees on the one-year WMAP ILC map, using a correlation statistic and an rms difference statistic. Extreme value distributions of these statistics are calculated for left-handed and right-handed 36 degree `screw motions' (Clifford translations) when matching circles and for a zero (unphysical) rotation. The most correlated circles appear for circle radii of 11\pm1 degrees, for the left-handed screw motion, but not for the right-handed one, nor for the zero rotation. The favoured six dodecahedral face centres in galactic coordinates are (l,b)= (252, +65), (51, +51), (144,+38), (207,+10), (271,+3), (332,+25) and their opposites. The six pairs of circles_independently_ each favour a circle angular radius of 11\pm1 degrees. Whether or not these six circle pairs centred on dodecahedral faces match via a 36 degree rotation only due to unexpected statistical properties of the WMAP ILC map, or whether they match due to global geometry, it is clear that the WMAP ILC map has some unusual statistical properties which mimic a potentially interesting cosmological signal. The software for reading the WMAP data and for carrying out this analysis are released under the GNU General Public License.

preprint2000arXiv

Tangential Large Scale Structure as a Standard Ruler: Curvature Parameters from Quasars

Several observational analyses suggest that matter is spatially structured at $\approx 130h^{-1}Mpc$ at low redshifts. This peak in the power spectrum provides a standard ruler in comoving space which can be used to compare the local geometry at high and low redshifts, thereby constraining the curvature parameters. It is shown here that this power spectrum peak is present in the observed quasar distribution at $z\sim 2$: qualitatively, via wedge diagrams which clearly show a void-like structure, and quantitatively, via one-dimensional Fourier analysis of the quasars' tangential distribution. The sample studied here contains 812 quasars. The method produces strong constraints (68% confidence limits) on the density parameter $Ω_0$ and weaker constraints on the cosmological constant $λ_0$, which can be expressed by the relation $ Ω_0 = (0.24 \pm0.15) + (0.10\pm0.08) λ_0 $. Independently of $λ_0$ (in the range $λ_0 \in [0,1]$), the constraint is $0.1 < Ω_0 < 0.45$. Combination of the present results with SN Type Ia results yields $Ω_0 = (0.30\pm0.11) + (0.57\pm0.11) (λ_0-0.7),$ $0.55 < λ_0 < 0.95,$ (68% confidence limits). This strongly supports the possibility that the observable universe satisfies a nearly flat, perturbed Friedmann-Lema\^ıtre-Robertson-Walker model, independently of any cosmic microwave background observations.

preprint2000arXiv

The Topology of the Universe

The Hilbert-Einstein equations are insufficient to describe the geometry of the Universe, as they only constrain a local geometrical property: curvature. A global knowledge of the geometry of space, if possible, would require measurement of the topology of the Universe. Since the subject was discussed in 1900 by Schwarzschild, observational attempts to measure global topology have been rare for most of this century, but have accelerated in the 1990's due to the rapidly increasing amount of observations of non-negligible fractions of the observational sphere. A brief review of basic concepts of cosmic topology and of the rapidly growing gamut of diverse and complementary observational strategies for measuring the topology of the Universe is provided here.

preprint1999arXiv

Constraining Curvature Parameters via Topology

If the assumption that physical space has a trivial topology is dropped, then the Universe may be described by a multiply connected Friedmann-Lema\^ıtre model on a sub-horizon scale. Specific candidates for the multiply connected space manifold have already been suggested. How precisely would a significant detection of multiple topological images of a single object, or a region on the cosmic microwave background, (due to photons arriving at the observer by multiple paths which have crossed the Universe in different directions), constrain the values of the curvature parameters $Ω_0$ and $λ_0$? The way that the constraints on $Ω_0$ and $λ_0$ depend on the redshifts of multiple topological images and on their radial and tangential separations is presented and calculated. The tangential separations give the tighter constraints: multiple topological images of known types of astrophysical objects at redshifts $z \ltapprox 3$ would imply values of $Ω_0$ and $λ_0$ preciser than $\sim 1%$ and $\sim 10%$ respectively. Cosmic microwave background `spots' identified with lower redshift objects by the Planck or MAP satellites would provide similar precision. This method is purely geometrical: no dynamical assumptions (such as the virial theorem) are required and the constraints are independent of the Hubble constant, $H_0.$

preprint1999arXiv

Evolution of the correlation function as traced by the HDF galaxies

An initially highly biased value of the spatial correlation function, which decreases with time up to a transition redshift z_t ~ 1-2, was predicted theoretically at least as early as 1993, and shown to be consistent with the HDF-N angular correlation function estimates in 1997. Observational analyses are presented here which show (i) an HDF-N estimate of the correlation function amplitude of galaxies selected at z ~ 2 (by the UV drop-{\em in} technique), which supports the estimate z_t ~ 2 and (ii) an HDF-N estimate of the correlation function amplitude of galaxies selected at z ~ 3.7 by using photometric redshifts, which suggests that the correlation function amplitude evolves as (1+z)^{2\pm3.5} during epochs earlier than z_t.

preprint1999arXiv

Galaxy Clustering at z ~ 2 and Halo Radii

The amplitude of the angular two-point galaxy correlation function w(θ) for galaxies at z~2 is estimated for galaxies in the Hubble Deep Field by using a U < 27 complete sub-sample. (i) It is confirmed that the amplitude of the correlation can be corrected for the integral constraint without having to make assumptions about the shape of the correlation function and by avoiding the introduction of linear error terms. The estimate using this technique is w(5'') = 0.10 \pm 0.09. (ii) If the biases introduced in faint galaxy selection due to obscuration by large objects are not corrected for by masking areas around them, then the estimate would be w(5'') =0.16\pm 0.07. (iii) The effective (3-D) galaxy pair separation at 5'' and this redshift range is ~ 25-250 /h kpc, so the correction to the spatial correlation function ξ(r) due to exclusion of overlapping galaxy dark matter haloes should be considered. For clustering stable in proper units in an Ω=1,λ=0 universe, our w(5\arcs) estimate (a) implies a present-day correlation length of r_0 ~ 2.6^{+1.1}_{-1.7}/h Mpc if halo overlapping is ignored, but (b) for a present-day correlation length of r_0=5.5/h Mpc implies that a typical halo exclusion radius is r_halo=70^{+420}_{-30}/h kpc. (iv) The decreasing correlation period (DCP) of a high initial bias in the spatial correlation function is not detected at this redshift. For an Ω=1,λ=0 universe and (proper) stable clustering, possible detections of the DCP in other work would imply that ξat redshifts greater than z_t = 1.7\pm0.9 would be [(1+z)/(1+z_t)]^{2.1\pm3.6} times higher than at z_t, which is consistent with our lack of a detection at z ~ 2.

preprint1998arXiv

Three-dimensional Topology-Independent Methods to Look for Global Topology

The space-like hypersurface of the Universe at the present cosmological time is a three-dimensional manifold. A non-trivial global topology of this space-like hypersurface would imply that the apparently observable universe (the sphere of particle horizon radius) could contain several images of the single, physical Universe. Recent three-dimensional techniques for constraining and/or detecting this topology are reviewed. Initial applications of these techniques using X-ray bright clusters of galaxies and quasars imply (weak) candidates for a non-trivial topology.

preprint1997arXiv

Constraining Cosmological Topology via Highly Luminous X-ray Clusters

The topology of the observable Universe is not yet known. The most significant observational sign of a non-trivial topology would be multiple images (``ghosts'') of a single object at (in general) different sky positions and redshifts. It is pointed out that the previous search by Gott (1980) for ghost images of the Coma cluster can be extended by using highly X-ray luminous clusters of galaxies. This is likely to be more efficient than with other astrophysical objects viewable on these scales since (1) X-ray clusters would be at least as easy to identify if viewed from other angles as any other objects and (2) the X-ray emitting thermally heated gas is likely to be simpler than for other objects. Possibilities that the highly luminous cluster RXJ~1347.5-1145 ($z=0\.45$) has a ``ghost image'' at lower redshift are analysed. It is noted that RXJ~1347.5-1145, the Coma cluster and the cluster CL~09104+4109 form nearly a right angle ($\approx 88°$) with arms of nearly identical length ($970h^{-1}$ and $960h^{-1}$ Mpc respectively) for $Ω_0=1, λ_0=0$ curvature ($h\equiv H_0/100 kms^{-1} Mpc^{-1}$). This is a clue that the three clusters could be ghost images of one and the same cluster, for a hypertoroidal topology. However, several arguments are presented that this relation is not physical.

preprint1997arXiv

Merging History Trees of Dark Matter Haloes: a Tool for Exploring Galaxy Formation Models

A method of deriving and using merging history trees of dark matter galaxy haloes directly from pure gravity N-body simulations is presented. This combines the full non-linearity of N-body simulations with the flexibility of the semi-analytical approach. Merging history trees derived from power-law initial perturbation spectrum simulations (for indices $n=-2$ and $n=0$) by Warren et al. (1992) are shown. As an example of a galaxy formation model, these are combined with evolutionary stellar population synthesis, via simple scaling laws for star formation rates, showing that if most star formation occurs during merger-induced bursts, then a nearly flat faint-end slope of the galaxy luminosity function may be obtained in certain cases. Interesting properties of hierarchical halo formation are noted: (1) In a given model, merger rates may vary widely between individual haloes, and typically 20%-30% of a halo's mass may be due to infall of uncollapsed material. (2) Small mass haloes continue to form at recent times: as expected, the existence of young, low redshift, low metallicity galaxies (e.g., Izotov et al. 1997) is consistent with hierarchical galaxy formation models. (3) For $n=-2,$ the halo spatial correlation function can have a very high initial bias due to the high power on large scales.

preprint1997arXiv

On the Effective Spatial Separations in the Clustering of Faint Galaxies

Several recent measurements have been made of the angular correlation function $w(θ,m)$ of faint galaxies in deep surveys (e.g., Hubble Deep Field, HDF). Are the measured correlations indicative of gravitational growth of primordial perturbations or of the relationship between galaxies and (dark matter dominated) galaxy haloes? A first step in answering this question is to determine the typical spatial separations of galaxies whose spatial correlations, $ξ(r,z),$ contribute most to w. The median spatial separation of galaxy pairs contributing to a fraction $p$ of $w$ for a galaxy survey is denoted by r_{eff} and compared with the perpendicular distance, $r_\perp,$ at the median redshift, z_{median}, of the galaxies. Over a wide range in $ξ$ growth rates $ε$ and in $z_{median}, \r_eff(50%) \ltsim r_\perp$ and $\r_eff(90%) \sim 4 r_\perp.$ Values of r_eff for typical redshift distributions indicate that many $w$ measurements correspond to spatial correlations at comoving length scales well below 1 h^-1 Mpc. For $Ω_0=1$ and $λ_0=0$, the correlation signal at 4\arcsec predominant in the Villumsen et al. (1996) estimates of $w$ for faint HDF galaxies corresponds to r_eff(50%)(HDF) \approx 40 h^-1 kpc; other cosmologies and angles up to 10\arcsec can increase this to r_eff(50%)(HDF) $\ltsim 200 h^-1 kpc. The proper separations are (1+z) times smaller. These scales are small: the faint galaxy $w$ measurements are at scales where halo and/or galaxy existence, let alone interactions, may modify the spatial correlation function. These measurements could maybe be used to probe the radial extent of haloes at high redshift.