We report on the detection of pulsations of three pulsating subdwarf B stars observed by the TESS satellite and our results of mode identification in these stars based on an asymptotic period relation. SB 459 (TIC 067584818), SB 815 (TIC 169285097) and PG 0342+026 (TIC 457168745) have been monitored during single sectors resulting in 27 days coverage. These datasets allowed for detecting, in each star, a few tens of frequencies, which we interpreted as stellar oscillations. We found no multiplets, though we partially constrained mode geometry by means of period spacing, which recently became a key tool in analyses of pulsating subdwarf B stars. Standard routine that we have used allowed us to select candidates for trapped modes that surely bear signatures of non-uniform chemical profile inside the stars. We have also done statistical analysis using collected spectroscopic and asteroseismic data of previously known subdwarf B stars along with our three stars. Making use of high precision trigonometric parallaxes from the Gaia mission and spectral energy distributions we converted atmospheric parameters to stellar ones. Radii, masses and luminosities are close to their canonical values for extreme horizontal branch stars. In particular, the stellar masses are close to the canonical one of 0.47 M$_\odot$ for all three stars but uncertainties on the mass are large. The results of the analyses presented here will provide important constrains for asteroseismic modelling.

Cepheid masses derived from pulsations or binary dynamics are generally lower than those derived from stellar evolution models. Recent efforts have been dedicated to investigating the effects of abundances, mass loss, rotation, convection and overshooting prescriptions for modifying the evolution tracks to reduce or remove this Cepheid mass discrepancy. While these approaches are promising, either alone or in combination, more work is required to distinguish between possible solutions. Here we investigate nuclear reaction rate and opacity modifications on Cepheid evolution using the MESA code. We discuss the effects of opacity increases at envelope temperatures of 200,000-400,000 K proposed to explain the pulsation properties of hybrid main-sequence beta Cep/Slowly Pulsating B (SPB) variables which will evolve into Cepheids. We make use of the RSP nonlinear radial pulsation modeling capability in MESA to calculate periods and radial velocity amplitudes of Galactic Cepheids V1334 Cyg, Polaris, and delta Cep.

Cosmological observational analysis frequently assumes that the Universe is spatially flat. We aim to non-perturbatively check the conditions under which a flat or nearly flat expanding dust universe, including the LCDM model if interpreted as strictly flat, forbids the gravitational collapse of structure. We quantify spatial curvature at turnaround. We use the Hamiltonian constraint to determine the pointwise conditions required for an overdensity to reach its turnaround epoch in an exactly flat spatial domain. We illustrate this with a plane-symmetric, exact, cosmological solution of the Einstein equation, extending earlier work. More generally, for a standard initial power spectrum, we use the relativistic Zel'dovich approximation implemented in 'inhomog' to numerically estimate how much positive spatial curvature is required for turnaround to be allowed at typical epochs/length scales in almost-EdS and almost-LCDM models with inhomogeneous curvature. We find that gravitational collapse in a spatially exactly flat, irrotational, expanding, dust universe is relativistically forbidden pointwise. In the spatially flat plane-symmetric model considered here, pancake collapse is excluded both pointwise and in averaged domains. In an almost-EdS or LCDM model, the per-domain average curvature in collapsing domains almost always becomes strongly positive prior to turnaround, with the expansion-normalised curvature functional reaching $\Omega_{\cal R}^{\cal D} \sim -5$. We show analytically that a special case gives $\Omega_{\cal R}^{\cal D} = -5$ exactly (if normalised using the EdS expansion rate) at turnaround. An interpretation of LCDM as literally 3-Ricci flat would forbid structure formation. The difference between relativistic cosmology and a strictly flat LCDM model is fundamental in principle, but we find that the geometrical effect is weak.

Aim: Our aim is to obtain high-accuracy measurements of the physical and orbital parameters of two evolved eclipsing binary systems, and to use these measurements to study their evolutionary status. We also aim to derive distances to the systems by using a surface brightness - colour relation and compare these distances with the measurements provided by GAIA. Methods: We measured the physical and orbital parameters on both systems based on V-band and I-band photometry from OGLE, near-infrared photometry obtained with the NTT telescope and the instrument SOFI, as well as high-resolution spectra obtained at ESO 3.6m/HARPS and Clay 6.5/MIKE spectrographs. The light curves and radial-velocity curves were analysed with the Wilson-Devinney code. Results: We analysed two double-lined eclipsing binary systems OGLE-BLG-ECL-23903 and OGLE-BLG-ECL-296596 from the Optical Gravitational Lensing Experiment (OGLE) catalogue. Both systems have a configuration of two well-detached giant stars. The masses of the components ofOGLE-BLG-ECL-123903 are M_1= 2.045 $\pm$ 0.027 and M_2=2.074 $\pm$ 0.023 $M_\odot$ and the radii are R_1=9.540 $\pm$ 0.049 and R_2=9.052 $\pm$ 0.060 $R_\odot$. For OGLE-BLG-ECL-296596, the masses are M_1=1.093 $\pm$ 0.015 and M_2=1.125 $\pm$ 0.014 $M_\odot$, while the radii are R_1=18.06 $\pm$ 0.28 and R_2=29.80 $\pm$ 0.33 $R_\odot$. Evolutionary status was discussed based on the isochrones and evolutionary tracks from PARSEC and MESA codes. The ages of the systems were establishes to be around 1.3 Gyr for the OGLE-BLG-ECL-123903 and 7.7 Gyr for the OGLE-BLG-ECL-296596. We also determined the distance to both systems. For OGLE-BLG-ECL-123903 this is equal to d=2.95 $\pm$ 0.06 (stat.) $\pm$ 0.07 (syst.) kpc, while for the OGLE-BLG-ECL-296596 it is d=5.68 $\pm$ 0.07 (stat.) $\pm$ 0.14 (syst.) kpc. This is the first analysis of its kind for these unique evolved eclipsing binary systems.

This paper investigates the growth of the most massive cosmological objects. We utilize the Simsilun simulation, which is based on the approximation of the silent universe. In the limit of spatial homogeneity and isotropy the silent universes reduce to the standard FLRW models. We show that within the approximation of the silent universe the formation of the most massive cosmological objects differs from the standard background-dependent approaches. For objects with masses above $10^{15} M_\odot$ the effect of spatial curvature (overdense regions are characterized by positive spatial curvature) leads to measurable effects. The effect is analogous to the effect that the background cosmological model has on the formation of these objects (i.e. the higher matter density and spatial curvature the faster the growth of cosmic structures). We measure this by the means of the mass function and show that the mass function obtained from the Simsilun simulation has a higher amplitude at the high-mass end compared to standard mass function such as the Press-Schechter or the Tinker mass function. For comparison, we find that the expected mass of most massive objects using the Tinker mass function is $4.4^{+0.8}_{-0.6} \, \times 10^{15} M_\odot$, whereas for the Simsilun simulation is $6.3^{+1.0}_{-0.8} \, \times 10^{15} M_\odot$. We conclude that the nonlinear relativistic effects could affect the formation of the most massive cosmological objects, leading to a relativistic environment-dependence of the growth rate of the most massive clusters.

We present a determination of precise fundamental physical parameters of twenty detached, double- lined, eclipsing binary stars in the Large Magellanic Cloud (LMC) containing G- or early K-type giant stars. Eleven are new systems, the remaining nine are systems already analyzed by our team for which we present updated parameters. The catalogue results from our long-term survey of eclipsing binaries in the Magellanic Clouds suitable for high-precision determination of distances (the Araucaria project). The V-band brightnesses of the systems range from 15.4 mag to 17.7 mag and their orbital periods range from 49 days to 773 days. Six systems have favorable geometry showing total eclipses. The absolute dimensions of all eclipsing binary components are calculated with a precision of better than 3% and all systems are suitable for a precise distance determination. The measured stellar masses are in the range 1.4 to 4.6 M_sun and comparison with the MESA isochrones gives ages between 0.1 and 2.1 Gyr. The systems show some weak age-metallicity relation. Two systems have components with very different masses: OGLE LMC-ECL-05430 and OGLE LMC-ECL-18365. Neither system can be fitted by single stellar evolution isochrone, explained by a past mass transfer scenario in the case of ECL-18365 and a gravitational capture or a hierarchical binary merger scenario in the case of ECL-05430. The longest period system OGLE LMC SC9 230659 shows a surprising apsidal motion which shifts the apparent position of the eclipses. In one spectrum of OGLE LMC-ECL-12669 we noted a peculiar dimming of one of the components by 65% well outside of the eclipses. We interpret this observation as arising from an extremely rare occultation event as a foreground Galactic object covers only one component of an extragalactic eclipsing binary.

Time-series photometry from space in the ultraviolet can be presently done with only a few platforms, none of which is able to provide wide-field long-term high-cadence photometry. We present a concept of UVSat, a twin space telescope which will be capable to perform this kind of photometry, filling an observational niche. The satellite will host two telescopes, one for observations in the ultraviolet, the other for observations in the optical band. We also briefly show what science can be done with UVSat.

Self-consistent treatment of cosmological structure formation and expansion within the context of classical general relativity may lead to "extra" expansion above that expected in a structureless universe. We argue that in comparison to an early-epoch, extrapolated Einstein-de Sitter model, about 10-15% "extra" expansion is sufficient at the present to render superfluous the "dark energy" 68% contribution to the energy density budget, and that this is observationally realistic.

The current cosmological model ($\Lambda$CDM) with the underlying FLRW metric relies on the assumption of local isotropy, hence homogeneity of the Universe. Difficulties arise when one attempts to justify this model as an average description of the Universe from first principles of general relativity, since in general, the Einstein tensor built from the averaged metric is not equal to the averaged stress--energy tensor. In this context, the discrepancy between these quantities is called "cosmological backreaction" and has been the subject of scientific debate among cosmologists and relativists for more than $20$ years. Here we present one of the methods to tackle this problem, i.e. averaging the scalar parts of the Einstein equations, together with its application, the cosmological mass function of galaxy clusters.

We revise the evolutionary status of the B-type supergiant HD\,163899 based on the new determinations of the mass-luminosity ratio, effective temperature and rotational velocity as well as on the interpretation of the oscillation spectrum of the star. The observed value of the nitrogen-to-carbon abundance fixes the value of the rotation rate of the star. Now, more massive models are strongly preferred than those previously considered and it is very likely that the star is still in the main sequence stage. The rotationally induced mixing manifests as the nitrogen overabundance in the atmosphere, which agrees with our analysis of the HARPS spectra. Thus, HD\,163899 belongs probably to a group of evolved nitrogen-rich main-sequence stars.

In relativistic inhomogeneous cosmology, structure formation couples to average cosmological expansion. A conservative approach to modelling this assumes an Einstein--de Sitter model (EdS) at early times and extrapolates this forward in cosmological time as a "background model" against which average properties of today's Universe can be measured. This requires adopting an early-epoch--normalised background Hubble constant $H_1^{bg}$. Here, we show that the $\Lambda$CDM model can be used as an observational proxy to estimate $H_1^{bg}$ rather than choose it arbitrarily. We assume (i) an EdS model at early times; (ii) a zero dark energy parameter; (iii) bi-domain scalar averaging---division of the spatial sections into over- and underdense regions; and (iv) virialisation (stable clustering) of collapsed regions. We find $H_1^{bg}= 37.7 \pm 0.4$ km/s/Mpc (random error only) based on a Planck $\Lambda$CDM observational proxy. Moreover, since the scalar-averaged expansion rate is expected to exceed the (extrapolated) background expansion rate, the expected age of the Universe should be much less than $2/(3 H_1^{bg}) = 17.3$ Gyr. The maximum stellar age of Galactic Bulge microlensed low-mass stars (most likely: 14.7 Gyr; 68\% confidence: 14.0--15.0 Gyr) suggests an age about a Gyr older than the (no-backreaction) $\Lambda$CDM estimate.

We present a spectroscopic survey of known and candidate $\gamma$\,Doradus stars. The high-resolution, high signal-to-noise spectra of 52 objects were collected by five different spectrographs. The spectral classification, atmospheric parameters (\teff, $\log g$, $\xi$), $v\sin i$ and chemical composition of the stars were derived. The stellar spectral and luminosity classes were found between G0-A7 and IV-V, respectively. The initial values for \teff\ and \logg\ were determined from the photometric indices and spectral energy distribution. Those parameters were improved by the analysis of hydrogen lines. The final values of \teff, \logg\ and $\xi$ were derived from the iron lines analysis. The \teff\ values were found between 6000\,K and 7900\,K, while \logg\,values range from 3.8 to 4.5\,dex. Chemical abundances and $v\sin i$ values were derived by the spectrum synthesis method. The $v\sin i$ values were found between 5 and 240\,km\,s$^{-1}$. The chemical abundance pattern of $\gamma$\,Doradus stars were compared with the pattern of non-pulsating stars. It turned out that there is no significant difference in abundance patterns between these two groups. Additionally, the relations between the atmospheric parameters and the pulsation quantities were checked. A strong correlation between the $v\sin i$ and the pulsation periods of $\gamma$\,Doradus variables was obtained. The accurate positions of the analysed stars in the H-R diagram have been shown. Most of our objects are located inside or close to the blue edge of the theoretical instability strip of $\gamma$\,Doradus.

The general relativistic description of cosmological structure formation is an important challenge from both the theoretical and the numerical point of views. In this paper we present a brief prescription for a general relativistic treatment of structure formation and a resulting mass function on galaxy cluster scales in a highly generic scenario. To obtain this we use an exact scalar averaging scheme together with the relativistic generalization of Zel'dovich's approximation (RZA) that serves as a closure condition for the averaged equations.

Backreaction in the cosmological context is a longstanding problem that is especially important in the present era of precise cosmology. The standard model of a homogeneous background plus density perturbations is most probably oversimplified and is expected to fail to fully account for the near-future observations of sub-percent precision. From a theoretical point of view, the problem of backreaction is very complicated and deserves careful examination. Recently, Green and Wald claimed in a series of papers to have developed a formalism to properly describe the influence of density inhomogeneities on average properties of the Universe, i.e., the backreaction effect. A brief discussion of this framework is presented, focussing on its drawbacks and on misconceptions that have arisen during the "backreaction debate".

No. In a number of papers Green and Wald argue that the standard FLRW model approximates our Universe extremely well on all scales, except close to strong field astrophysical objects. In particular, they argue that the effect of inhomogeneities on average properties of the Universe (backreaction) is irrelevant. We show that this latter claim is not valid. Specifically, we demonstrate, referring to their recent review paper, that (i) their two-dimensional example used to illustrate the fitting problem differs from the actual problem in important respects, and it assumes what is to be proven; (ii) the proof of the trace-free property of backreaction is unphysical and the theorem about it fails to be a mathematically general statement; (iii) the scheme that underlies the trace-free theorem does not involve averaging and therefore does not capture crucial non-local effects; (iv) their arguments are to a large extent coordinate-dependent, and (v) many of their criticisms of backreaction frameworks do not apply to the published definitions of these frameworks. It is therefore incorrect to infer that Green and Wald have proven a general result that addresses the essential physical questions of backreaction in cosmology.

We present the new interpretation of the oscillation spectrum of HD 163899 based on the new determinations of the effective temperature, mass-luminosity ratio and rotational velocity. These new parameters strongly prefer the more massive models than previously considered. Now it is also possible that the star could be in the main sequence stage. Using the oscillation spectrum as a gauge, we intend to establish which stage of evolution corresponds better to HD 163899.

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.

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.

Evolutionary tracks and pulsational analysis of models with masses of 13-18 $M_\odot$ are presented. We address two important questions. The first one deals with one of the most unresolved problems in astrophysics, i.e., the existence of a blue loop after core helium ignition; the so called "to loop or not to loop" problem. We show that inward overshooting from the outer convective zone in the red giant phase is prerequisite for the development of the blue loop. Our second question concerns pulsational instability of models in the core helium burning phase. We present for the first time that models on the blue loop can have unstable modes driven by the $\kappa$ mechanism operating in the $Z-$bump. Contrary to post-main sequence models in the shell hydrogen burning phases, pulsational instability of the blue loop models depends mainly on effective temperature and metallicity is of secondary importance. Finally, we try to interpret the oscillation spectrum of the blue supergiant HD 163899, the only member of the SPBsg class, and to get some clue on the evolutionary status of the star.

We present results of pulsation analyses of B-type supergiant models with masses of $14 - 18 M_\odot$, considering evolutionary stages before and after helium core ignition. Using a non-adiabatic pulsation code, we compute instability domains for low degree modes. For selected models in these two evolutionary phases, we compare properties of pulsation modes. Significant differences are found in oscillation spectra and the kinetic energy density of pulsation modes.

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.

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.

As a response to the white paper call, we propose to turn Kepler to the South Ecliptic Pole (SEP) and observe thousands of large amplitude variables for years with high cadence in the frame of the Kepler-SEP Mission. The degraded pointing stability will still allow observing these stars with reasonable (probably better than mmag) accuracy. Long-term continuous monitoring already proved to be extremely helpful to investigate several areas of stellar astrophysics. Space-based missions opened a new window to the dynamics of pulsation in several class of pulsating variable stars and facilitated detailed studies of eclipsing binaries. The main aim of this mission is to better understand the fascinating dynamics behind various stellar pulsational phenomena (resonances, mode coupling, chaos, mode selection) and interior physics (turbulent convection, opacities). This will also improve the applicability of these astrophysical tools for distance measurements, population and stellar evolution studies. We investigated the pragmatic details of such a mission and found a number of advantages: minimal reprogramming of the flight software, a favorable field of view, access to both galactic and LMC objects. However, the main advantage of the SEP field comes from the large sample of well classified targets, mainly through OGLE. Synergies and significant overlap (spatial, temporal and in brightness) with both ground- (OGLE, LSST) and space-based missions (GAIA, TESS) will greatly enhance the scientific value of the Kepler-SEP mission. GAIA will allow full characterization of the distance indicators. TESS will continuously monitor this field for at least one year, and together with the proposed mission provide long time series that cannot be obtained by other means. If Kepler-SEP program is successful, there is a possibility to place one of the so-called LSST "deep-drilling" fields in this region.

We present results of the instability analysis of the post-main sequence massive star models against radial and nonradial pulsations. We confirm that both p- and g-modes can be excited by the $\kappa$-mechanism acting in the metal opacity bump. However, as opposed to the previous claims, we find that an intermediate convective zone (ICZ) related to the hydrogen burning shell is not necessary for excitation of g-modes. These modes can be reflected at a minimum of the Brunt-V\"ais\"al\"a frequency, located at the top of the chemical composition gradient region surrounding the radiative helium core. This minimum is associated with the change of actual temperature gradient from the adiabatic value in the semiconvective zone to the radiative value above it. Thus, the existence of pulsations at this evolutionary stage does not prove the existence of the convective zone but only some reflective layer. Finally, we show that no regular patterns can be expected in oscillation spectra of blue supergiant pulsators but there is a prospect for identification of the mode degree, $\ell$, from multicolour photometry.

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.

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\'e 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\zeta(3) G(M/L^3)\,x$, where $\zeta(3)$ is Ap\'ery'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$.