Prof Subir Sarkar

Abstract:

AbstractIn the late 1990’s, observations of two directionally-skewed samples of, in total, 93 Type Ia supernovae were analysed in the framework of the Friedmann–Lemaître–Robertson–Walker (FLRW) cosmology. Assuming these to be ‘standard(isable) candles’ it was inferred that the Hubble expansion rate is accelerating as if driven by a positive Cosmological Constant $$\varLambda $$ Λ in Einstein’s theory of gravity. This is still the only direct evidence for the ‘dark energy’ that is the dominant component of today’s standard $$\varLambda $$ Λ CDM cosmological model. Other data such as baryon acoustic oscillations (BAO) in the large-scale distribution of galaxies, temperature fluctuations in the cosmic microwave background (CMB), measurement of stellar ages, the rate of growth of structure, etc are all ‘concordant’ with this model but do not provide independent evidence for accelerated expansion. The recent discussions about whether the inferred acceleration is real rests on analysis of a larger sample of 740 SNe Ia which shows that these are not quite standard candles, and more importantly highlights the ‘corrections’ that are applied to analyse the data in the FLRW framework. The latter holds in the reference frame in which the CMB is isotropic, whereas observations are carried out in our heliocentric frame in which the CMB has a large dipole anisotropy. This is assumed to be of kinematic origin i.e. due to our non-Hubble motion driven by local inhomogeneity in the matter distribution which has grown under gravity from primordial density perturbations traced by the CMB fluctuations. The $$\varLambda $$ Λ CDM model predicts how this peculiar velocity should fall off as the averaging scale is raised and the universe becomes sensibly homogeneous. However observations of the local ‘bulk flow’ are inconsistent with this expectation and convergence to the CMB frame is not seen. Moreover, the kinematic interpretation implies a corresponding dipole in the sky distribution of high redshift quasars, which is rejected by observations at $$4.9\sigma $$ 4.9 σ . Hence the peculiar velocity corrections employed in supernova cosmology are inconsistent and discontinuous within the data. The acceleration of the Hubble expansion rate is in fact anisotropic at $$3.9\sigma $$ 3.9 σ and aligned with the bulk flow. Thus dark energy could be an artefact of analysing data assuming that we are idealised observers in an FLRW universe, when in fact the real universe is inhomogeneous and anisotropic out to distances large enough to impact on cosmological analyses.

Abstract:

We present a high-energy neutrino event generator, called LeptonInjector, alongside an event weighter, called LeptonWeighter. Both are designed for large-volume Cherenkov neutrino telescopes such as IceCube. The neutrino event generator allows for quick and flexible simulation of neutrino events within and around the detector volume, and implements the leading Standard Model neutrino interaction processes relevant for neutrino observatories: neutrino-nucleon deep-inelastic scattering and neutrino-electron annihilation. In this paper, we discuss the event generation algorithm, the weighting algorithm, and the main functions of the publicly available code, with examples. Program summary: Program Titles: LeptonInjector and LeptonWeighter CPC Library link to program files: https://doi.org/10.17632/662gkpjfd9.1 Developer's repository links: https://github.com/icecube/LeptonInjector and https://github.com/icecube/LeptonWeighter Licensing provisions: GNU Lesser General Public License, version 3. Programming Language: C++11 External Routines: • Boost • HDF5 • nuflux (https://github.com/icecube/nuflux) • nuSQuIDS (https://github.com/arguelles/nuSQuIDS) • Photospline (https://github.com/icecube/photospline) • SuiteSparse (https://github.com/DrTimothyAldenDavis/SuiteSparse) Nature of problem: LeptonInjector: Generate neutrino interaction events of all possible topologies and energies throughout and around a detector volume. LeptonWeighter: Reweight Monte Carlo events, generated by a set of LeptonInjector Generators, to any desired physical neutrino flux or cross section. Solution method: LeptonInjector: Projected ranges of generated leptons and the extent of the detector, in terms of column depth, are used to inject events in and around the detector volume. Event kinematics follow distributions provided in cross section files. LeptonWeighter: Event generation probabilities are calculated for each Generator, which are then combined into a generation weight and used to calculate an overall event weight.