Next Article in Journal
The Effect of Higher Dimensional QCD Operators on the Spectroscopy of Bottom-Up Holographic Models
Next Article in Special Issue
Our Peculiar Motion Inferred from Number Counts of Mid Infra Red AGNs and the Discordance Seen with the Cosmological Principle
Previous Article in Journal
Perturbations of the Gravitational Energy in the TEGR: Quasinormal Modes of the Schwarzschild Black Hole
Previous Article in Special Issue
Quantum Physics Literacy Aimed at K12 and the General Public
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gravitational Interaction in the Chimney Lattice Universe †

1
Department of Mathematics and Physics, North Carolina Central University, 1801 Fayetteville St., Durham, NC 27707, USA
2
Department of Physics, Istanbul Technical University, Maslak, 34469 Istanbul, Turkey
3
Institute of Theoretical Astrophysics, University of Oslo, P.O. Box 1029 Blindern, N-0315 Oslo, Norway
4
Astronomical Observatory, Odessa I.I. Mechnikov National University, Dvoryanskaya St. 2, 65082 Odessa, Ukraine
*
Author to whom correspondence should be addressed.
This paper is an extended version from the proceeding paper: Eingorn, M.; McLaughlin, A., II; Canay, E.; Brilenkov, M.; Zhuk, A. Gravitation in the space with chimney topology. In Proceedings of the 1st Electronic Conference on Universe, Online, 22–28 February 2021.
Universe 2021, 7(4), 101; https://doi.org/10.3390/universe7040101
Submission received: 30 March 2021 / Revised: 10 April 2021 / Accepted: 12 April 2021 / Published: 15 April 2021

Abstract

:
We investigate the influence of the chimney topology T × T × R of the Universe on the gravitational potential and force that are generated by point-like massive bodies. We obtain three distinct expressions for the solutions. One follows from Fourier expansion of delta functions into series using periodicity in two toroidal dimensions. The second one is the summation of solutions of the Helmholtz equation, for a source mass and its infinitely many images, which are in the form of Yukawa potentials. The third alternative solution for the potential is formulated via the Ewald sums method applied to Yukawa-type potentials. We show that, for the present Universe, the formulas involving plain summation of Yukawa potentials are preferable for computational purposes, as they require a smaller number of terms in the series to reach adequate precision.
PACS:
04.25.Nx—Post-Newtonian approximation; perturbation theory; related approximations; 98.80.Jk—Mathematical and relativistic aspects of cosmology

1. Introduction

The shape of the space, whether it is positively curved, negatively curved, or flat, and whether there is a limit to the size of the Universe are all among essential topics of contemporary debate in theoretical physics and cosmology. Spatial topology of the Universe, its function at the very early stages of evolution (in the quantum gravity regime), and in the later process of large scale structure formation are quite interesting questions yet to be answered. General Relativity does not favor any particular topology; hence, on theoretical grounds, the space might be simply connected, in agreement with concordance cosmology, or equally as well, multiply connected. It is worth noting that some “exotic” non-simply connected spacetimes are timelike geodesically incomplete, since they have singularities [1]. Hence, such topologies are not viable.
If the Universe is multiply connected, it may have a finite volume and yet be negatively curved or flat [2]. The current available data cannot reveal the finiteness of its volume if the Universe covers a much wider region than the observable sector. However, a rather smaller volume points at the possibility of finding observational indications of its topological features [3]. For instance, a photon can travel plenty of times across the volume of multiply connected space and, thus, generate multiple images of the emitting source as a signature [4,5]. Spaces with toroidal topology in one to three dimensions may be presented as common examples of multiply connected spaces. To this class belongs the three-torus T × T × T , chimney T × T × R , and slab T × R × R topologies.
There are various comprehensive studies in the literature on potential indicators of the shape of the space [6,7,8,9,10], and the majority of research is focused on their relation to the Cosmic Microwave Background (CMB) data. Indeed, there exists a very appealing conjectural relation which suggests that CMB anomalies in large angular scale observations, e.g., the suppression of the quadrupole moment and the quadrupole and octopole alignment, are imprints of spatial topology [11,12]. The weak wide-angle temperature correlations in the CMB can be also explained, e.g., with the help of dodecahedral topology of the Universe [13]. In the present work we study the chimney topology, which admits a single infinite axis subject to interpretations such as the preferred direction of the quadrupole and octopole alignment and the commonly named “axis of evil” [14] (see [15] for additional observable signatures of a preferred axis).
In connection with the investigation of possible topological imprints in CMB observations, Planck 2013 data [2] place the constraint R i > 0.71 χ rec on the radius of the largest sphere that may be inscribed in the topological domain for a flat Universe with the equal-sided chimney topology. The parameter χ rec specifies the distance from the recombination surface and it is of the same order with the particle horizon, that is, χ rec 14 Gpc . For the toroidal topologies, the former restrictions on the size of the Universe from the seven- and nine-year WMAP temperature map analyses are presented in [14,16]. The smallest possible size of the fundamental topological domain for flat space, according to the seven-year WMAP results, is d = 2 R LSS cos α min 27.9 Gpc [12], where R LSS stands for the distance from the last scattering surface.
In this paper, we study the chimney topology T × T × R in terms of the gravitational characteristics of the Universe, manifested in the shape of the gravitational potential and force. In the cosmological setting, the inhomogeneous gravitational field is sourced by fluctuations in the matter density [17] and, as expected, in the Newtonian limit, the potential satisfies the Poisson equation. The form of the gravitational potential in the case of toroidal topologies was previously studied in [18]. Particularly, the authors have shown that there exists no physically justified nontrivial solution of the Poisson equation for the T × T × R model. On the other hand, by employing the perturbed Einstein equations from the very beginning, one automatically includes the essential relativistic effects in the formulation and, for the gravitational potential, obtains a Helmholtz-type equation instead of the Poisson one [19,20,21]. Quite remarkably, as we show in the present work for the chimney topology, it then becomes possible to obtain exact solutions of this equation that are nontrivial and physically meaningful. Herein, we derive distinct expressions for the gravitational potential and force through alternative methods and point out the particular solutions appearing in the form of summed Yukawa potentials as ready-to-use notable sources for numerical computations. It is worth mentioning that, in the above outlined approach, we make no presumptions regarding the spatial distribution of gravitating bodies.
The outline of the paper is as follows. In Section 2, following [22], we introduce the main equations and derive alternative expressions for the gravitational potential using distinct methods, now including the Ewald technique. Subsequently, in Section 3.1, extending the results of [22], we compare these expressions in view of their usefulness for numerical computations. In Section 3.2, we obtain the gravitational force expressions for each form of the potential. We briefly review the results of our work in the concluding Section 4.

2. Methods

The Model and Basic Equations

It is well known that the gravitational potential Φ , created by fluctuations in the matter density, is defined by scalar perturbations of the metric coefficients [23] and that in the framework of General Relativity, it satisfies the linearized Einstein equation (see, e.g., [24,25]). Ignoring peculiar velocities, in the case of the Λ CDM cosmological model, this equation reads [19,20,21]
Δ Φ 0 3 κ ρ ¯ c 2 2 a Φ 0 = κ c 2 2 a ρ ρ ¯ ,
where κ 8 π G N / c 4 (with G N and c being the Newtonian gravitational constant and the speed of light, respectively), a is the scale factor, while Δ represents the Laplace operator in comoving coordinates. Here, ρ and ρ ¯ = const are the comoving mass density and its averaged value, respectively. As we operate within the Λ CDM model, matter is pressureless, and we consider it in the form of discrete point-like gravitating bodies with masses m n to represent, e.g., galaxies. Therefore, the comoving mass density
ρ = n m n δ ( r r n ) .
The 0 subscript of Φ in Equation (1) refers to the fact that peculiar velocities have been disregarded (see also [26]).
The shifted gravitational potential
Φ ^ 0 Φ 0 1 3
fulfils the equation
Δ Φ ^ 0 a 2 λ 2 Φ ^ 0 = κ c 2 2 a ρ ,
where the screening length [19]
λ 3 κ ρ ¯ c 2 2 a 3 1 / 2 .
The presence of the term Φ 0 in Equation (1) (consequently, Φ ^ 0 in Equation (4)) results in the Yukawa-type cutoff of the potential with the characteristic length λ . The term Φ enters as a summand into energy-momentum fluctuations generating metric perturbations [19].
In what follows, the overhat indicates that the gravitational potential is shifted. A significant bonus of working with the shifted potential is that it is now possible to employ the superposition principle in solving Equation (4): once we find a solution for a single particle that is located at the center of Cartesian coordinates, we may immediately generalize it for a collection of particles at random positions.
We consider the space with chimney topology T 1 × T 2 × R , where the tori T 1 and T 2 have periods l 1 and l 2 along, e.g., the x- and y-axes, respectively. Hence, each gravitating body has its images positioned away from the original point in multiples of periods l 1 and l 2 along the corresponding axes. Now, let us place a particle with mass m at the center of Cartesian coordinates. For the above indicated topology, the delta functions δ ( x ) and δ ( y ) may be presented as
δ ( x ) = 1 l 1 k 1 = + cos 2 π k 1 l 1 x , δ ( y ) = 1 l 2 k 2 = + cos 2 π k 2 l 2 y ,
which implicitly include the contribution from the images of the particle. Consequently, Equation (4) for this particle reads
Δ Φ ^ 0 a 2 λ 2 Φ ^ 0 = κ c 2 2 a m l 1 l 2 k 1 = + k 2 = + cos 2 π k 1 l 1 x cos 2 π k 2 l 2 y δ z ,
so, we are motivated to consider the solution
Φ ^ 0 = k 1 = + k 2 = + C k 1 k 2 ( z ) cos 2 π k 1 l 1 x cos 2 π k 2 l 2 y ,
where the coefficients C k 1 k 2 ( z ) satisfy the equation
k 1 = + k 2 = + C k 1 k 2 ( z ) 4 π 2 k 1 2 l 1 2 + k 2 2 l 2 2 C k 1 k 2 ( z ) a 2 λ 2 C k 1 k 2 ( z ) κ c 2 2 a m l 1 l 2 δ z × cos 2 π k 1 l 1 x cos 2 π k 2 l 2 y = 0 .
Using the condition d 2 | z | / d z 2 = 2 δ ( z ) , we can easily obtain the explicit expressions for the coefficients C k 1 k 2 ( z ) , so that the shifted gravitational potential for the selected particle and all its images eventually reads
Φ ^ 0 = κ c 2 4 a m l 1 l 2 k 1 = + k 2 = + 4 π 2 k 1 2 l 1 2 + k 2 2 l 2 2 + a 2 λ 2 1 / 2 × exp 4 π 2 k 1 2 l 1 2 + k 2 2 l 2 2 + a 2 λ 2 | z | cos 2 π k 1 l 1 x cos 2 π k 2 l 2 y .
The above expression has the correct behavior in the Newtonian limit in the neighborhood of the considered particle, where it is no longer possible to distinguish between different (infinite and periodic) axes. For such regions, the summations in (10) may be replaced by the integrals:
Φ ^ 0 κ c 2 m 4 a + d k x + d k y 4 π 2 k x 2 + k y 2 + a 2 λ 2 1 / 2 × exp 4 π 2 k x 2 + k y 2 + a 2 λ 2 | z | cos 2 π k x x + k y y
for k x k 1 / l 1 and k y k 2 / l 2 . Introducing the vectors k = ( k x , k y ) , r = ( x , y ) with the absolute values k = k x 2 + k y 2 and r = x 2 + y 2 , and assuming an angle φ between them, we obtain
Φ ^ 0 κ c 2 m 4 a 0 + k d k 4 π 2 k 2 + a 2 λ 2 1 / 2 exp 4 π 2 k 2 + a 2 λ 2 | z | × 0 2 π d φ cos 2 π k r cos φ = π κ c 2 m 2 a 0 + k d k 4 π 2 k 2 + a 2 λ 2 1 / 2 exp 4 π 2 k 2 + a 2 λ 2 | z | J 0 ( 2 π k r ) = G N m c 2 1 Z 2 + R 2 exp 1 λ Z 2 + R 2 G N m c 2 1 Z 2 + R 2 ,
where Z = a z and R = a r represent the physical distances, and the last integration is performed by using the formula 2.12.10(10) of [27].
Evidently, for a system of randomly positioned gravitating bodies, we have
Φ 0 = 1 3 + Φ ^ 0 m n m n ; x , y , z x x n , y y n , z z n .
For linear fluctuations, the averaged value of this expression is equal to zero, as it should be (also see [28]). Indeed,
0 l 1 d x 0 l 2 d y + d z Φ ^ 0 = κ c 2 4 a m λ a + exp a λ | z | d z = κ c 2 2 a m λ 2 a 2 exp a λ z 0 + = κ c 2 2 a m λ 2 a 2 = 1 3 m ρ ¯ ,
and, hence, the spatial average of the total gravitational potential is
Φ 0 ¯ = 1 3 1 3 m ρ ¯ · N l 1 l 2 L z = 1 3 1 3 = 0 , m N l 1 l 2 L z = ρ ¯ .
For the sake of simplicity, here we have considered the particular configuration in which all N bodies in the volume V = l 1 l 2 L z are assigned identical masses m.
Equation (4) is of Helmholtz type and we can likewise solve it by considering the contribution of periodic images. In this case, the resulting expression consists of summed Yukawa potentials attributed to each one of them:
Φ ^ 0 = κ c 2 m 8 π a k 1 = + k 2 = + 1 ( x k 1 l 1 ) 2 + ( y k 2 l 2 ) 2 + z 2 × exp a ( x k 1 l 1 ) 2 + ( y k 2 l 2 ) 2 + z 2 λ .
As we have noted previously, the peculiar motion of gravitating bodies is disregarded in Equation (1) and, consequently, in (4). Nevertheless, the significance of such a contribution has recently been pointed out in [29], where the authors have also shown that peculiar velocities may be effectively restored by employing the effective screening length λ eff (given by the formula (41) of [29]) instead of the screening length λ in Equations (1) and (4). Specifically, in the matter-dominated epoch, the two quantities λ eff and λ are related to one another as λ eff = 3 / 5 λ . Returning to our formulation, the effect of peculiar motion is included by replacing λ with λ eff in the formulas (10) and (16), which yields
Φ ˜ cos κ c 2 8 π a m l 1 Φ ^ cos = k 1 = + k 2 = + k 1 2 + k 2 2 + 1 4 π 2 λ ˜ eff 2 1 / 2 × exp 4 π 2 k 1 2 + k 2 2 + 1 λ ˜ eff 2 | z ˜ | cos 2 π k 1 x ˜ cos 2 π k 2 y ˜
and
Φ ˜ exp κ c 2 8 π a m l 1 Φ ^ exp = k 1 = + k 2 = + 1 x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 × exp x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 λ ˜ eff .
For simpler demonstration, we have assumed l 1 = l 2 = l and introduced the rescaled quantities
x = x ˜ l , y = y ˜ l , z = z ˜ l , λ eff = λ ˜ eff a l .
Two alternative solutions are labeled with the subscripts “cos” and “exp” in Equations (17) and (18). Now that the peculiar velocities are included in the calculations, the 0 subscript is omitted in the new formulas.
There is also a third way to express the gravitational potential for the given topology. Indeed, Yukawa-type interactions that are subject to periodic boundary conditions can be formulated via Ewald sums, so that the expression for the potential consists of two rapidly converging series, one in each of the real and Fourier spaces. The technique is commonly employed while modeling particle interactions in plasma and colloids, and, in such a context, the corresponding potential for quasi two-dimensional systems, i.e., three-dimensional systems with two-dimensional periodicity, has previously been derived in [30,31]. Being implemented in the cosmological setting considered in our paper, the discussed expression for the “Yukawa–Ewald” potential reads
Φ ˜ mix κ c 2 8 π a m l 1 Φ ^ mix = k 1 = + k 2 = + D x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 ; α ; λ ˜ eff 2 x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 + π cos 2 π k 1 x ˜ + k 2 y ˜ F 4 π 2 ( k 1 2 + k 2 2 ) + λ ˜ eff 2 ; z ˜ ; α 4 π 2 ( k 1 2 + k 2 2 ) + λ ˜ eff 2 ,
where
D x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 ; α ; λ ˜ eff exp x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 λ ˜ eff × erfc α x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 + 1 2 α λ ˜ eff + exp x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 λ ˜ eff × erfc α x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 1 2 α λ ˜ eff
and
F 4 π 2 ( k 1 2 + k 2 2 ) + λ ˜ eff 2 ; z ˜ ; α exp 4 π 2 ( k 1 2 + k 2 2 ) + λ ˜ eff 2 z ˜ erfc 4 π 2 ( k 1 2 + k 2 2 ) + λ ˜ eff 2 2 α + α z ˜ + exp 4 π 2 ( k 1 2 + k 2 2 ) + λ ˜ eff 2 z ˜ erfc 4 π 2 ( k 1 2 + k 2 2 ) + λ ˜ eff 2 2 α α z ˜ .
In these formulas, erfc is the complementary error function and the free parameter α , as indicated in [31], is to be chosen in such a way that a balanced interplay of computational cost and satisfactory precision is achieved. For definiteness, we set α equal to λ ˜ eff .
In the forthcoming section, we will compare three expressions and present the optimum formula in view of its efficiency in use for numerical analysis.

3. Results

3.1. Gravitational Potentials

Formulas (17), (18) and (20) describe the gravitational potential due to a point-like body, with mass m, placed at ( x , y , z ) = ( 0 , 0 , 0 ) , and by the accompanying images placed at ( x , y , z ) = ( k 1 l , k 2 l , 0 ) , where k 1 , 2 = 0 , ± 1 , ± 2 , All three forms of the rescaled potential are composed of infinite series. Hence, for any desired precision, one needs to determine the minimum number of terms that are needed to numerically calculate the potential. The criterion that we use to specify this number n is the following: the ratio | exact Φ ˜ approximate Φ ˜ | / | exact Φ ˜ | should be less than 0.001. This defines the order of accuracy in our analysis. Evidently, for each form of the potential the number n can be different, so we denote these as n exp , n cos and n mix , correspondingly. The formula requiring the smallest number of terms to define Φ ˜ (i.e., to calculate “approximate Φ ˜ ”) up to the adopted accuracy is clearly the best alternative for numerical computation purposes. In this connection, we are interested in comparing (17), (18) and (20) here. Because the formulas include double series, the accompanying numbers n are to be ascribed the smallest possible number of combinations ( k 1 , k 2 ) that can provide the necessary precision. We find these by listing the summands in increasing order for k 1 2 + k 2 2 and assigning to n the total number of terms that are included in the list eventually. This procedure of generating a sequence of combinations ( k 1 , k 2 ) and finding n is performed using Mathematica [32].
Table 1 and Table 2 show the outputs for eight selected points. As to the adopted accuracy, for all n n exp , the approximate value of Φ ˜ exp (calculated by (18)) differs from the exact value by less than one tenth of a percent. In both tables, the exact value Φ ˜ is calculated from (18) for n n exp . The quantities n cos and n mix indicate the numbers of terms in formulas (17) and (20), respectively, which one needs to keep in order to obtain the same values of the potential at the selected points with the same precision as attained by using (18). In the n cos column, the dash reflects either incorrect outputs that are produced because of complications in the computational process, or the fact that unreasonably large number of summands is necessary. Because our results depend on the ratio of λ eff to the physical size a l of the periods of tori, i.e., λ ˜ eff = λ eff / ( a l ) , we present the results that were obtained for small and large values of λ ˜ eff separately in Table 1 and Table 2, which include λ ˜ eff = 0.01 , 0.1 and λ ˜ eff = 1 , 3 , respectively.
According to these tables, both expressions (18) and (20) seem preferable for numerical calculations in the case λ ˜ eff < 1 since n exp , n mix n cos , although Equation (20) is, of course, much more complicated than Equation (18), and the computation of its every single summand takes longer. However, for λ ˜ eff 1 , the Yukawa–Ewald formula (20) alone becomes superior to the remaining two and the distinction grows as λ ˜ eff becomes larger. According to Planck 2013 data [2], the lower limit on the periods of tori (in the case of chimney topology) is ∼20 Gpc . Meanwhile, the current value of the effective cosmological screening length, as indicated in [29], is 2.6 Gpc. Thus, the region that is defined by λ ˜ eff < 1 depicts the observable Universe and, here, as we have just discussed, Equation (18) is more convenient for numerical analysis.
Concluding this section, we also present Figure 1, Figure 2, Figure 3 and Figure 4, demonstrating the shape of the rescaled potential Φ ˜ for the same values of λ ˜ eff as those picked for Table 1 and Table 2. To plot these figures (using Mathematica [32]), we use (18) for n n exp .

3.2. Gravitational Forces

It is also interesting to study the forces (per unit mass) associated with the alternative forms of the gravitational potential derived in the previous section. We intend to consider the projections of these forces on the x- and z-axes. Owing to the symmetry of the model, the x and y projections are similar. We calculate the gravitational forces for the same points as for the potentials and, among these, naturally, we only investigate the points at which projections on the axis of interest are nonzero. In this connection, the points A 1 , A 2 , A 3 , C 1 , C 2 and the points A 3 , B 3 are omitted for the x- and z-components, respectively. The accuracy of force calculations is of the same level as that of the potentials. Here, once again, we compare the number of terms needed to achieve this accuracy, but now for three different forms of the gravitational force presentation.

3.2.1. x-Component of the Gravitational Force

From (17), (18) and (20), we derive three alternative expressions for the x-component of the rescaled force:
x ˜ Φ ˜ cos = 2 π k 1 = + k 2 = + k 1 2 + k 2 2 + 1 4 π 2 λ ˜ eff 2 1 / 2 × exp 4 π 2 k 1 2 + k 2 2 + 1 λ ˜ eff 2 | z ˜ | k 1 sin 2 π k 1 x ˜ cos 2 π k 2 y ˜ ,
x ˜ Φ ˜ exp = k 1 = + k 2 = + exp x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 λ ˜ eff × x ˜ k 1 x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 3 / 2 + x ˜ k 1 λ ˜ eff x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 ,
x ˜ Φ ˜ mix = 1 2 k 1 = + k 2 = + x ˜ k 1 D x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 ; α ; λ ˜ eff x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 3 / 2 + x ˜ k 1 x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 · C exp x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 λ ˜ eff + x ˜ k 1 x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 · C + exp x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 λ ˜ eff + 4 π 2 k 1 sin 2 π k 1 x ˜ + k 2 y ˜ F 4 π 2 k 1 2 + k 2 2 + λ ˜ eff 2 ; z ˜ ; α 4 π 2 k 1 2 + k 2 2 + λ ˜ eff 2 ,
where
C = C x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 ; α ; λ ˜ eff 2 α π exp α x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 1 2 α λ ˜ eff 2 ± 1 λ ˜ eff erfc α x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 1 2 α λ ˜ eff .
We present the results of the calculations that were performed in Mathematica [32] in Table 3 and Table 4 for λ ˜ eff = 0.01 , 0.1 and λ ˜ eff = 1 , 3 , respectively. As in the case for the gravitational potential, a straightforward analysis shows that formulas (24) and (25) that are related to the Yukawa and Yukawa–Ewald potentials, respectively, are preferable over (23) for the physically significant case λ ˜ eff < 1 (although the structure of the expression (25) is again much more complicated compared to (24)). Meanwhile, when λ ˜ eff 1 , the Yukawa–Ewald force becomes superior. In both tables, Equation (24) was employed for n n exp while computing the values of the rescaled x-component Φ ˜ x .
Additionally, we present Figure 5, Figure 6, Figure 7 and Figure 8, demonstrating the x-component of the rescaled force Φ ˜ x for the same values of λ ˜ eff as those picked for Table 3 and Table 4. To plot these figures (using Mathematica [32]), we employ the formula (24) for n n exp .

3.2.2. z-Component of the Gravitational Force

For the z-component of the rescaled force, three alternative formulas are:
z ˜ Φ ˜ cos = 2 π k 1 = + k 2 = + exp 4 π 2 k 1 2 + k 2 2 + 1 λ ˜ eff 2 z ˜ × cos 2 π k 1 x ˜ cos 2 π k 2 y ˜ ,
where, for simplicity, z ˜ > 0 ,
z ˜ Φ ˜ exp = k 1 = + k 2 = + exp x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 λ ˜ eff × z ˜ x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 3 / 2 + z ˜ λ ˜ eff x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2
and
z ˜ Φ ˜ mix = 1 2 k 1 = + k 2 = + z ˜ D x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 ; α ; λ ˜ eff x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 3 / 2 + z ˜ x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 · C exp x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 λ ˜ eff + z ˜ x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 · C + exp x ˜ k 1 2 + y ˜ k 2 2 + z ˜ 2 λ ˜ eff 2 π cos 2 π k 1 x ˜ + k 2 y ˜ 1 4 π 2 k 1 2 + k 2 2 + λ ˜ eff 2 × F exp z 4 π 2 k 1 2 + k 2 2 + 1 λ ˜ eff 2 + F + exp z 4 π 2 k 1 2 + k 2 2 + 1 λ ˜ eff 2 ,
where
F = F 4 π 2 k 1 2 + k 2 2 + λ ˜ eff 2 ; z ˜ ; α ± 2 α π exp 4 π 2 k 1 2 + k 2 2 + λ ˜ eff 2 2 α α z ˜ 2 4 π 2 k 1 2 + k 2 2 + λ ˜ eff 2 erfc 4 π 2 k 1 2 + k 2 2 + λ ˜ eff 2 2 α α z ˜ ,
and C are given by (26).
Now, we employ these formulas to calculate the nonzero z-components of the gravitational force at the previously selected set of points and, again, for the desired precision. The results that were obtained in Mathematica [32] are presented in Table 5 and Table 6, which show that while λ ˜ eff < 1 , depicting well the observational restrictions, two formulas (28) and (29) are favorable (as before, the latter is much more cumbersome). On the other hand, for λ ˜ eff 1 , the Yukawa–Ewald formula (29) gives the best results. Herein, the quantity Φ ˜ z stands for the z-component of the rescaled force, calculated from Equation (28) for n n exp . We depict the behavior of this component in Figure 9 and Figure 10 for the section y = 0 . Obviously, the projection of the gravitational force on the z-axis is absent for the section z = 0 due to the symmetry of the model.

4. Conclusions

In this work, we have studied how the chimney topology T × T × R of the Universe affects the form of the gravitational potential and, consequently, that of the gravitational force. In this connection, we have proposed three alternative forms for each of the solutions. One of them (see Equation (17)) relies on the Fourier expansion of the delta functions into series while using periodicity in two toroidal dimensions in the model. The second one (see Equation (18)) follows from the summation of solutions of the Helmholtz equation, each in the form of the Yukawa potential, for a source mass and all of its periodic images. Finally, the third form of the potential (see Equation (20)) is formulated via the Ewald sums for Yukawa potentials. Subsequently, we have presented three alternative forms of the gravitational force (see Equations (23)–(25) and (27)–(29) for the x- and z-components, respectively) derived from the potential expressions.
In all three alternative forms, the screening length λ ˜ eff serves as a crucial parameter, as it specifies the distance (from the source or the periodic images) where the gravitational potential undergoes exponential cutoff. This fact is most clearly demonstrated in the case where the solution takes on the form of summed Yukawa potentials (18). The observational data show that this parameter should be less than 1 ( λ ˜ eff < 1 ) in today’s Universe.
One of the main goals of this work was to reveal which of the obtained alternative formulas would serve better as a tool to be employed in numerical calculations. Namely, to show which formula would require less terms in the series to reach the desired precision. Our calculations have demonstrated that, for both the gravitational potentials and forces, two expressions involving plain summations of Yukawa potentials are preferable for the physically significant case λ ˜ eff < 1 . However, in the case λ ˜ eff 1 , the Yukawa–Ewald presentation stands out as the best alternative.
Additionally, we have produced Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10 for λ ˜ eff = 0.01 , 0.1 , 1 , 3 to provide graphical demonstration of the gravitational potentials and force projections.

Author Contributions

Conceptualization, M.E.; methodology, M.E., E.C. and A.Z.; formal analysis, M.E., A.M.II, E.C., M.B. and A.Z.; investigation, M.E., A.M.II, E.C., M.B. and A.Z.; writing—original draft preparation, A.Z.; writing—review and editing, M.E. and E.C.; visualization, M.E. and A.M.II; supervision, M.E. and A.Z.; project administration, M.E.; funding acquisition, M.E. All authors have read and agreed to the published version of the manuscript.

Funding

The work of M.E. and A.M.II was supported by the National Science Foundation HRD Award number 1954454.

Conflicts of Interest

The authors have no conflicts of interest to declare that are relevant to the content of this article.

References

  1. Gannon, D. Singularities in nonsimply connected space-times. J. Math. Phys. 1975, 16, 2364. [Google Scholar] [CrossRef]
  2. Ade, P.A.R.; et al. [Planck Collaboration] Planck 2013 results. XXVI. Background geometry and topology of the Universe. Astron. Astrophys. 2014, 571, A26. [Google Scholar]
  3. Luminet, J.-P. The shape and topology of the Universe. arXiv 2008, arXiv:0802.2236. [Google Scholar]
  4. Roukema, B.F. Topology of the Universe: Background and recent observational approaches. Pramana 1999, 53, 945–949. [Google Scholar] [CrossRef] [Green Version]
  5. Piechocki, W. Topology of the Universe. arXiv 1999, arXiv:gr-qc/9910055. [Google Scholar]
  6. Bielewicz, P.; Riazuelo, A. The study of topology of the universe using multipole vectors. Mon. Not. R. Astron. Soc. 2009, 396, 609–623. [Google Scholar] [CrossRef] [Green Version]
  7. Bielewicz, P.; Banday, A.J.; Gorski, K.M. Constraining the topology of the Universe using the polarized CMB maps. Mon. Not. R. Astron. Soc. 2012, 421, 1064–1072. [Google Scholar] [CrossRef] [Green Version]
  8. Vaudrevange, P.M.; Starkman, G.D.; Cornish, N.J.; Spergel, D.N. Constraints on the topology of the Universe: Extension to general geometries. Phys. Rev. D 2012, 86, 083526. [Google Scholar] [CrossRef] [Green Version]
  9. Fabre, O.; Prunet, S.; Uzan, J.-P. Topology beyond the horizon: How far can it be probed? Phys. Rev. D 2015, 92, 043003. [Google Scholar] [CrossRef] [Green Version]
  10. Bucher, M.; et al. [Planck Collaboration] Planck 2015 results. XVIII. Background geometry and topology. Astron. Astrophys. 2016, 594, A18. [Google Scholar]
  11. Bielewicz, P.; Banday, A.J. Constraints on the topology of the Universe derived from the 7-year WMAP data. Mon. Not. R. Astron. Soc. 2011, 412, 2104–2110. [Google Scholar] [CrossRef] [Green Version]
  12. Bielewicz, P.; Banday, A.J.; Gorski, K.M. Constraints on the topology of the Universe derived from the 7-year WMAP CMB data and prospects of constraining the topology using CMB polarization maps. In Proceedings of the XLVIIth Rencontres de Moriond, La Tuile, Italy, 3–10 March 2012; Auge, E., Dumarchez, J., Tran Thanh Van, J., Eds.; ARISF: Paris, France, 2012; p. 91. [Google Scholar]
  13. Luminet, J.P.; Weeks, J.R.; Riazuelo, A.; Lehoucq, R.; Uzan, J.P. Dodecahedral space topology as an explanation for weak wide-angle temperature correlations in the cosmic microwave background. Nature 2003, 425, 593–595. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Aslanyan, G.; Manohar, A.V. The topology and size of the Universe from the Cosmic Microwave Background. J. Cosmol. Astropart. Phys. 2012, 06, 003. [Google Scholar] [CrossRef] [Green Version]
  15. Floratos, E.G.; Leontaris, G.K. On topological modifications of Newton’s law. J. Cosmol. Astropart. Phys. 2012, 2012, 24. [Google Scholar] [CrossRef]
  16. Aslanyan, G.; Manohar, A.V.; Yadav, A.P.S. The topology and size of the Universe from CMB temperature and polarization data. J. Cosmol. Astropart. Phys. 2013, 2013, 9. [Google Scholar] [CrossRef]
  17. Peebles, P.J.E. The Large-Scale Structure of the Universe; Princeton University Press: Princeton, NJ, USA, 1980. [Google Scholar]
  18. Brilenkov, M.; Eingorn, M.; Zhuk, A. Lattice Universe: Examples and problems. Eur. Phys. J. C 2015, 75, 217. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Eingorn, M. First-order cosmological perturbations engendered by point-like masses. Astrophys. J. 2016, 825, 84. [Google Scholar] [CrossRef] [Green Version]
  20. Eingorn, M.; Kiefer, C.; Zhuk, A. Scalar and vector perturbations in a universe with discrete and continuous matter sources. J. Cosmol. Astropart. Phys. 2016, 2016, 32. [Google Scholar] [CrossRef] [Green Version]
  21. Eingorn, M.; Kiefer, C.; Zhuk, A. Cosmic screening of the gravitational interaction. Int. J. Mod. Phys. D 2017, 26, 1743012. [Google Scholar] [CrossRef] [Green Version]
  22. Eingorn, M.; McLaughlin II, A.; Canay, E.; Brilenkov, M.; Zhuk, A. Gravitation in the space with chimney topology. In Proceedings of the 1st Electronic Conference on Universe, online, 22–28 February 2021. [Google Scholar]
  23. Landau, L.D.; Lifshitz, E.M. The Classical Theory of Fields; Course of Theoretical Physics Series, V.2; Pergamon Press: Oxford, UK, 2000. [Google Scholar]
  24. Mukhanov, V.F. Physical Foundations of Cosmology; Cambridge University Press: Cambridge, UK, 2005. [Google Scholar]
  25. Gorbunov, D.S.; Rubakov, V.A. Introduction to the Theory of the Early Universe: Cosmological Perturbations and Inflationary Theory; World Scientific: Singapore, 2011. [Google Scholar]
  26. Eingorn, M. Cosmological law of universal gravitation. Int. J. Mod. Phys. D 2017, 26, 1750121. [Google Scholar] [CrossRef] [Green Version]
  27. Prudnikov, A.P.; Brychkov, Y.A.; Marichev, O.I. Integrals and Series, Volume 2. Special Functions; Gordon and Breach Science Publishers: New York, NY, USA, 1992. [Google Scholar]
  28. Eingorn, M.; Brilenkov, M.; Vlahovic, B. Zero average values of cosmological perturbations as an indispensable condition for the theory and simulations. Eur. Phys. J. C 2015, 75, 381. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Canay, E.; Eingorn, M. Duel of cosmological screening lengths. Phys. Dark Universe 2020, 29, 100565. [Google Scholar] [CrossRef] [Green Version]
  30. Mazars, M. Ewald sums for Yukawa potentials in quasi-two-dimensional systems. J. Chem. Phys. 2007, 126, 056101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Mazars, M. Yukawa potentials in systems with partial periodic boundary conditions I: Ewald sums for quasi-two dimensional systems. Mol. Phys. 2007, 105, 1909–1925. [Google Scholar] [CrossRef] [Green Version]
  32. Wolfram Research, Inc. Mathematica; Version 11.3; Wolfram Research, Inc.: Champaign, IL, USA, 2018. [Google Scholar]
Figure 1. Φ ˜ = G N m / ( c 2 a l ) 1 Φ ^ for λ ˜ eff = 0.01 for the sections z = 0 (left panel) and y = 0 (right panel).
Figure 1. Φ ˜ = G N m / ( c 2 a l ) 1 Φ ^ for λ ˜ eff = 0.01 for the sections z = 0 (left panel) and y = 0 (right panel).
Universe 07 00101 g001
Figure 2. Φ ˜ = G N m / ( c 2 a l ) 1 Φ ^ for λ ˜ eff = 0.1 for the sections z = 0 (left panel) and y = 0 (right panel).
Figure 2. Φ ˜ = G N m / ( c 2 a l ) 1 Φ ^ for λ ˜ eff = 0.1 for the sections z = 0 (left panel) and y = 0 (right panel).
Universe 07 00101 g002
Figure 3. Φ ˜ = G N m / ( c 2 a l ) 1 Φ ^ for λ ˜ eff = 1 for the sections z = 0 (left panel) and y = 0 (right panel).
Figure 3. Φ ˜ = G N m / ( c 2 a l ) 1 Φ ^ for λ ˜ eff = 1 for the sections z = 0 (left panel) and y = 0 (right panel).
Universe 07 00101 g003
Figure 4. Φ ˜ = G N m / ( c 2 a l ) 1 Φ ^ for λ ˜ eff = 3 for the sections z = 0 (left panel) and y = 0 (right panel).
Figure 4. Φ ˜ = G N m / ( c 2 a l ) 1 Φ ^ for λ ˜ eff = 3 for the sections z = 0 (left panel) and y = 0 (right panel).
Universe 07 00101 g004
Figure 5. x-component of the rescaled force Φ ˜ x for λ ˜ eff = 0.01 for the sections z = 0 (left panel) and y = 0 (right panel).
Figure 5. x-component of the rescaled force Φ ˜ x for λ ˜ eff = 0.01 for the sections z = 0 (left panel) and y = 0 (right panel).
Universe 07 00101 g005
Figure 6. x-component of the rescaled force Φ ˜ x for λ ˜ eff = 0.1 for the sections z = 0 (left panel) and y = 0 (right panel).
Figure 6. x-component of the rescaled force Φ ˜ x for λ ˜ eff = 0.1 for the sections z = 0 (left panel) and y = 0 (right panel).
Universe 07 00101 g006
Figure 7. x-component of the rescaled force Φ ˜ x for λ ˜ eff = 1 for the sections z = 0 (left panel) and y = 0 (right panel).
Figure 7. x-component of the rescaled force Φ ˜ x for λ ˜ eff = 1 for the sections z = 0 (left panel) and y = 0 (right panel).
Universe 07 00101 g007
Figure 8. x-component of the rescaled force Φ ˜ x for λ ˜ eff = 3 for the sections z = 0 (left panel) and y = 0 (right panel).
Figure 8. x-component of the rescaled force Φ ˜ x for λ ˜ eff = 3 for the sections z = 0 (left panel) and y = 0 (right panel).
Universe 07 00101 g008
Figure 9. z-component of the rescaled force Φ ˜ z for λ ˜ eff = 0.01 and λ ˜ eff = 0.1 (left and right panels, respectively).
Figure 9. z-component of the rescaled force Φ ˜ z for λ ˜ eff = 0.01 and λ ˜ eff = 0.1 (left and right panels, respectively).
Universe 07 00101 g009
Figure 10. z-component of the rescaled force Φ ˜ z for λ ˜ eff = 1 and λ ˜ eff = 3 (left and right panels, respectively).
Figure 10. z-component of the rescaled force Φ ˜ z for λ ˜ eff = 1 and λ ˜ eff = 3 (left and right panels, respectively).
Universe 07 00101 g010
Table 1. Potentials Φ ˜ as well as the numbers n exp , n cos and n mix at a selection of points for λ ˜ eff = 0.01 and λ ˜ eff = 0.1 in the left and right tables, respectively.
Table 1. Potentials Φ ˜ as well as the numbers n exp , n cos and n mix at a selection of points for λ ˜ eff = 0.01 and λ ˜ eff = 0.1 in the left and right tables, respectively.
x ˜ y ˜ z ˜ Φ ˜ n exp n cos n mix x ˜ y ˜ z ˜ Φ ˜ n exp n cos n mix
A 1 0.500.5 5.524 × 10 31 210072 A 1 0.500.5 2.418 × 10 3 7407
A 2 0.500.1 2.810 × 10 22 22 A 2 0.500.1 2.398 × 10 2 68086
A 3 0.500 7.715 × 10 22 22 A 3 0.500 2.700 × 10 2 44
B 1 0.100.5 1.405 × 10 22 11871 B 1 0.100.5 1.203 × 10 2 4284
B 2 0.100.1 5.101 × 10 6 121191 B 2 0.100.1 1.719 13801
B 3 0.100 4.540 × 10 4 11 B 3 0.100 3.679 11
C 1 000.5 3.857 × 10 22 12361 C 1 000.5 1.353 × 10 2 4374
C 2 000.1 4.540 × 10 4 114791 C 2 000.1 3.679 14901
Table 2. Potentials Φ ˜ as well as the numbers n exp , n cos and n mix in a selection of points for λ ˜ eff = 1 and λ ˜ eff = 3 in the left and right tables, respectively.
Table 2. Potentials Φ ˜ as well as the numbers n exp , n cos and n mix in a selection of points for λ ˜ eff = 1 and λ ˜ eff = 3 in the left and right tables, respectively.
x ˜ y ˜ z ˜ Φ ˜ n exp n cos n mix x ˜ y ˜ z ˜ Φ ˜ n exp n cos n mix
A 1 0.500.5 3.783 174915 A 1 0.500.5 15.93 141876
A 2 0.500.1 5.067 16322915 A 2 0.500.1 17.60 13791209
A 3 0.500 5.153 16315 A 3 0.500 17.71 13779
B 1 0.100.5 3.990 1711013 B 1 0.100.5 16.14 141186
B 2 0.100.1 9.405 13316411 B 2 0.100.1 22.00 12901389
B 3 0.100 12.34 12310 B 3 0.100 24.96 12429
C 1 000.5 4.014 1701313 C 1 000.5 16.17 141087
C 2 000.1 12.30 1233579 C 2 000.1 24.91 12432869
Table 3. Numerical values of the x-component of the rescaled force Φ ˜ x Φ ˜ / x ˜ as well as the numbers n exp , n cos and n mix for points B 1 , B 2 and B 3 for λ ˜ eff = 0.01 and λ ˜ eff = 0.1 in the left and right tables, respectively.
Table 3. Numerical values of the x-component of the rescaled force Φ ˜ x Φ ˜ / x ˜ as well as the numbers n exp , n cos and n mix for points B 1 , B 2 and B 3 for λ ˜ eff = 0.01 and λ ˜ eff = 0.1 in the left and right tables, respectively.
x ˜ y ˜ z ˜ Φ ˜ x n exp n cos n mix x ˜ y ˜ z ˜ Φ ˜ x n exp n cos n mix
B 1 0.100.5 2.810 × 10 21 12631 B 1 0.100.5 2.783 × 10 2 5545
B 2 0.100.1 3.862 × 10 4 124481 B 2 0.100.1 20.75 15921
B 3 0.100 4.994 × 10 2 11 B 3 0.100 73.57 11
Table 4. Numerical values of the x-component of the rescaled force Φ ˜ x Φ ˜ / x ˜ as well as the numbers n exp , n cos and n mix for points B 1 , B 2 and B 3 for λ ˜ eff = 1 and λ ˜ eff = 3 in the left and right tables, respectively.
Table 4. Numerical values of the x-component of the rescaled force Φ ˜ x Φ ˜ / x ˜ as well as the numbers n exp , n cos and n mix for points B 1 , B 2 and B 3 for λ ˜ eff = 1 and λ ˜ eff = 3 in the left and right tables, respectively.
x ˜ y ˜ z ˜ Φ ˜ x n exp n cos n mix x ˜ y ˜ z ˜ Φ ˜ x n exp n cos n mix
B 1 0.100.5 4.730 × 10 1 1303821 B 1 0.100.5 4.920 × 10 1 8623821
B 2 0.100.1 34.65 205539 B 2 0.100.1 34.88 7755213
B 3 0.100 99.14 198 B 3 0.100 99.49 349
Table 5. Numerical values of the z-component of the rescaled force Φ ˜ z Φ ˜ / z ˜ as well as the numbers n exp , n cos and n mix for points A 1 , A 2 , B 1 , B 2 , C 1 and C 2 for λ ˜ eff = 0.01 and λ ˜ eff = 0.1 in the left and right tables, respectively.
Table 5. Numerical values of the z-component of the rescaled force Φ ˜ z Φ ˜ / z ˜ as well as the numbers n exp , n cos and n mix for points A 1 , A 2 , B 1 , B 2 , C 1 and C 2 for λ ˜ eff = 0.01 and λ ˜ eff = 0.1 in the left and right tables, respectively.
x ˜ y ˜ z ˜ Φ ˜ z n exp n cos n mix x ˜ y ˜ z ˜ Φ ˜ z n exp n cos n mix
A 1 0.500.5 3.962 × 10 29 210702 A 1 0.500.5 1.946 × 10 2 6476
A 2 0.500.1 5.620 × 10 21 22 A 2 0.500.1 5.620 × 10 2 216472
B 1 0.100.5 1.405 × 10 20 11871 B 1 0.100.5 1.407 × 10 1 2332
B 2 0.100.1 3.862 × 10 4 122281 B 2 0.100.1 20.75 16491
C 1 000.5 3.935 × 10 20 12401 C 1 000.5 1.620 × 10 1 3443
C 2 000.1 4.994 × 10 2 116201 C 2 000.1 73.58 17221
Table 6. Numerical values of the z-component of the rescaled force Φ ˜ z Φ ˜ / z ˜ as well as the numbers n exp , n cos and n mix for points A 1 , A 2 , B 1 , B 2 , C 1 and C 2 for λ ˜ eff = 1 and λ ˜ eff = 3 in the left and right tables, respectively.
Table 6. Numerical values of the z-component of the rescaled force Φ ˜ z Φ ˜ / z ˜ as well as the numbers n exp , n cos and n mix for points A 1 , A 2 , B 1 , B 2 , C 1 and C 2 for λ ˜ eff = 1 and λ ˜ eff = 3 in the left and right tables, respectively.
x ˜ y ˜ z ˜ Φ ˜ z n exp n cos n mix x ˜ y ˜ z ˜ Φ ˜ z n exp n cos n mix
A 1 0.500.5 3.571 852115 A 1 0.500.5 5.072 4442111
A 2 0.500.1 1.673 6490015 A 2 0.500.1 2.037 33186313
B 1 0.100.5 5.045 742013 B 1 0.100.5 6.593 397199
B 2 0.100.1 35.48 154445 B 2 0.100.1 36.04 574448
C 1 000.5 5.241 732613 C 1 000.5 6.795 3922412
C 2 000.1 99.97 86784 C 2 000.1 100.7 216777
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Eingorn, M.; McLaughlin, A., II; Canay, E.; Brilenkov, M.; Zhuk, A. Gravitational Interaction in the Chimney Lattice Universe. Universe 2021, 7, 101. https://doi.org/10.3390/universe7040101

AMA Style

Eingorn M, McLaughlin A II, Canay E, Brilenkov M, Zhuk A. Gravitational Interaction in the Chimney Lattice Universe. Universe. 2021; 7(4):101. https://doi.org/10.3390/universe7040101

Chicago/Turabian Style

Eingorn, Maxim, Andrew McLaughlin, II, Ezgi Canay, Maksym Brilenkov, and Alexander Zhuk. 2021. "Gravitational Interaction in the Chimney Lattice Universe" Universe 7, no. 4: 101. https://doi.org/10.3390/universe7040101

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop