^{1}

^{*}

^{2}

We propose a mathematical model for determining the probability of meteorite origin, impacting the earth. Our method is based on axioms similar to both the complex networks and emergent gravity. As a consequence, we are able to derive a link between complex networks and Newton’s gravity law, and as a possible application of our model we discuss several aspects of the Bacubirito meteorite. In particular, we analyze the possibility that the origin of this meteorite may be alpha Centauri system. Moreover, we find that in order for the Bacubirito meteorite to come from alpha Cen and be injected into our Solar System, its velocity must be reduced one order of magnitude of its ejected scape velocity from alpha Cen. There are several ways how this could happened, for example through collision with the Oort cloud objects (located outside the boundary of our Solar System), and/or through collisions within the Solar meteorites belt (located between Mars and Jupiter). We also argue that it may be interesting to study the Bacubirito meteorite from the perspective of the recently discovered Oumuamua object.

In this work we are interested in answering the question whether meteorites found on the earth actually came from the closest star Alpha Centauri (alpha Cen) system [

Alpha Centaury system consists of three stars: the pair alpha Cen A and alpha Cen B and Proxima Cen [

Recently, astronomers have reported an Earth-size planet orbiting one of the nearest stars in our galaxy; alpha Cen B [

It is worth mentioning that both alpha Cen A and B are high metallicity stars, which would promote the existence of circumstellar discs with a high fraction of solid materials (meteorites) at early times (see Refs. [

If circumstellar discs of alpha Cen system is capable of forming planets it is naturally to assume that there are also analogue asteroid belt, Kuiper belt or Oort cloud. From this point of view one must assume the possibility that some asteroids found on earth come from the alpha Cen system. Of course, the process can be vice versa in the sense that there has been interchange of asteroids between the alpha Cen system and the solar system. At present, roughly speaking the Oort cloud starts approximately at distance of 1 light-year and ends at 2 light-year from the sun. So one may assume that the analogue Oort cloud of the alpha Cen system starts also a distance of 1 light-year and ends at 2 light-year (from the alpha Centaury system). Considering that the average distance from the sun and alpha Centaury system is of 4 light-year it seems reasonable to assume that both Oort clouds are in interaction. If this is the case at present times presumably in the past the evolution of both systems must be of great importance in the sense that part of the matter from both systems must be interchanged. So, it is expected that some asteroids trapped in the solar system actually come from alpha Cen system. Consequently, although the meteorite origin probability found on earth is higher in the case of structure near to Earth such as the asteroid belt the probability is not necessarily zero if one assumes an alpha Cen origin of some meteorites. In this work we propose a mathematical model that could answer somehow this question. Such a model corresponds to complex network adapted to gravitational phenomena, and as a possible application of our gravitational complex network model, we consider the Bacubirito meteorite [

Finally, it is interesting to mention that a future interstellar spacecraft, including a fleet of StarChip spacecraft is currently being developed for a flyby mission by the Breakthrough Starshot project [

This work is organized as follows. In Section 2 we comment about complex networks. In Section 3, we obtain an expression for gravitational complex networks. In Section 4, we explore if in the case of the Bacubirito meteorite we can apply our formalism. In Section 6 we mention that the possibility that some meteorite find on earth may have interstellar origin as the Oumuamua object [

It is known that random networks with complex topology describe a wide range of systems in Nature. Surprisingly, recent advances in this scenario show that most large networks can be described by mean-field method applied to a system with scale-free features (see Refs. [

P ( k ) ~ 1 k γ , (1)

where P ( k ) is the probability that a vertex in the network is connected to k other vertices and γ is a numerical scale-free parameter so called “connectivity distribution exponent”. Random networks with complex topology are based in two principles:

1) Growth: starting with small number of vertices v 0 , at every time step t one adds a new vertex with e < v 0 edges that will be connected to the vertices already present in the system.

2) Preferential attachment: When choosing the vertices to which the new vertex connects, one assumes that the probability Π ( k i ) that a new vertex will be connected to vertex i depends on the connectivity (node degree) k i of that vertex and is given by

Π ( k i ) = k i ∑ j = 1 v 0 + t − 1 k j . (2)

(The reacher becomes reacher.) Observe that the sum in (2) goes over all vertices in the system except the new one. Assuming that k i is continuous parameter, one can assume that the variation of k i with respect the time is proportional to this probability

∂ k i ∂ t = e Π ( k i ) , (3)

where e is proportional constant. Thus, considering (2) we have

∂ k i ∂ t = e k i ∑ j = 1 v 0 + t − 1 k j . (4)

It is possible to show that

∑ j = 1 v 0 + t − 1 k j = 2 e t . (5)

Therefore, one gets equation

∂ k i ∂ t = k i 2 t , (6)

which has the following solution (given the condition k i ( t i ) = e )

k i ( t ) = e ( t t i ) 1 / 2 . (7)

Then, using this expression, the probability that a vertex has connectivity k i smaller than k can be written as

P ( k i ( t ) < k ) = P ( t i > e 2 t k 2 ) = 1 − P ( t i ≤ e 2 t k 2 ) = 1 − e 2 t k 2 ( v 0 + t ) , (8)

where we have assumed that the probability density for t i is P ( t i ) = 1 / ( v 0 + t ) . Making the differentiation of this expression with respect to k, one obtain the probability that a vertex in the network is connected to k other vertices

P ( k ) = ∂ P ( k i ( t ) < k ) ∂ k = α 1 k 3 , with α = ( 2 e 2 t v 0 + t ) . (9)

Comparing this expression with (1), one sees that in this model the free-scaling parameter becomes γ = 3 .

The idea of connecting gravity with networks has been of great interest through the years (see Ref. [

P ~ 1 k 3 , (10)

and

F ~ 1 r 2 , (11)

was considered. Consequently the possible relation between the radio r and the connectivity k was established,

r ~ k 3 / 2 . (12)

In fact, the expression (10) can be generalized in the form

P ~ 1 k γ , (13)

where, as it was mentioned in section 1, γ is just a free-scale parameter called the connectivity distribution exponent. It turns out that the scale-free parameter γ is a model dependent. For instance, in the observed networks the values in the range 2 ≤ γ ≤ 3 [

P ~ 1 k 2 (14)

and therefore one can make the identification r ~ k to obtain

P ~ G M m r 2 , (15)

where the constant of proportionality must have units inverse of force units. This expression can be interpreted as “the probability that object of mass m is connected to other object of mass M is inversely proportional to the square distance between the two masses”. Thus, from of point of view of complex networks the Newton gravitational law is the emergent probabilistic expression (A7) (see Appendix), which can be used to estimate the probability for a meteorite to impact the Earth from a given location, that is, to determine meteorite origin impacting the Earth.

The Bacubirito meteorite [

At present, it is not known the origin of the Bacubirito meteorite. Of course, one should expect that its most probable origin is the asteroid belt, the Kuiper belt or the Oort cloud. However, since it is considered anomalous iron meteorite we would like to leave open the possibility that its origin is the alpha Cen system. According to the Ref. [

P Alfa P Kuiper = ( r K r A ) 2 ≈ 10 − 10 , (16)

where M = 5.972 × 10 24 kg , m = 22 × 10 3 kg , r A = 4.13 × 10 16 m , and r K = 1.5 × 10 11 m . Moreover, alpha Cen A is one of the closest stellar systems and has spectral type very similar to the sun. There are observation results [

It is known that in order to understand the evolution of the Milky Way one uses the chemical composition of stars. Using this data, one concludes that neighbor stars are not isolated system but rather the vicinity of one another are affected by the same astrophysical events. So it is likely that alpha Cen system and the sun interact in different forms during their evolution formation. In particular one should expect that the matter interchange between both systems was a very possible scenario. In particular, in different time periods some asteroids of alpha Cen system could reach the solar system and vice versa. In this context, the study of the comparison of chemical abundance in both neighbors stars and different meteorites found of earth can be of great importance.

We would like also to analyze the possibility that the scape velocity play an important physical role in our search for a meteorite from alpha Cen enter to the solar system. As we mentioned Alpha Cen system is mainly made up of alpha Cen A, alpha Centauri B and Proxima Centauri which is a red dwarf. They masses are respectively m A = 1.1 M ⊙ , m B = 0.907 M ⊙ , and m P = 0.123 M ⊙ , where M ⊙ = 1.98855 × 10 30 kg . So, the total mass of the systems is M = 2.13 M ⊙ . Now, the scape velocity of a body associated with the system of mass M is given by

v = 2 G M r , (17)

where G is the gravitational constant ( G = 6.67 × 10 − 11 m 3 ⋅ kg − 1 ⋅ s − 2 ), M is the mass of the system, and r the center of mass of the body-system. If r ⊙ is the center of mass of a body-Solar system, the ratio between the scape velocities for both system is (total mass of our Solar System is dominated by the mass of the Sun)

v v ⊙ = M M ⊙ r ⊙ r . (18)

Therefore, one has that this comparison ratio between alfa Cen and our Solar System is

v v ⊙ ≈ 1.46, (19)

where it has been assumed that r ≈ r ⊙ one can guess that this object in alpha Centauri must be located between alpha Cen A and alpha Cen B where the non linear resonances can produce instabilities in the objects and can be ejected from the system, where the distance from alpha Cen A to alpha Cen B is about the distance from our Sun to Pluto, between 4.4 to 7.4 billion Km. In this way, any object escaping from Alpha Centauri in direction to our Solar System will have a speed such that the Solar System will not be able freely captured it into its system, unless directly hit the sun of a planet; in particular the earth. Thus, a direct hit of meteorite of the size of the Bacubirito meteorite from alpha Cen will produce a crater bigger than size left by the meteorite which made the Arizona crater (about one Km of diameter). However, there is not such a crater on the place where Bacubirito was found. Therefore, its enter energy in Earth atmosphere had to be much smaller. Let us see with some details: taking r ≈ 3 AU = 4.5 × 10 8 km the scape speed of this object coming out of Alpha Centauri is about v ≈ 35.4 km / s (making a trip for 37,000 years to get our Solar System, covering a distance of 41.3 × 10 12 km ). Therefore, its kinetic energy is approximately 10.5 × 10 12 J ≈ 2.6 × 10 3 Ton ( TNT ) which is equivalent to one Kton atomic bomb. Therefore, in order for the Bacubirito meteorite to come from alpha Cen, its velocity must be reduced about one order of magnitude of its scape speed. There are several ways this energy can be reduced, for example through collision with the Oort sloud objects (outside the plates orbits of our Solar System) and/or meteorites belt (between Mars and Jupiter) of our Solar System.

While we were preparing and refining the present article for publication we became aware of the surprisingly discovered, on 2017 October 19, of the first interstellar object called 1I/2017 U1 (Oumuamua) (see Refs. [

This work was partially supported by PROFAPI-UAS/2012. J. A. Nieto would like to thank J. H. Abundis, C. Añorve and E. Terán for helpful comments.

The authors declare no conflicts of interest regarding the publication of this paper.

López, G.V. and Nieto, J.A. (2018) Alpha Centauri System and Meteorites Origin. Journal of Applied Mathematics and Physics, 6, 2370-2381. https://doi.org/10.4236/jamp.2018.611199

Gravitational information theory

In order to explain this development let us recall how the Newton’s law of gravity can be obtained from information theory (see Ref. [^{2}. Moreover, the information is measure by bits. Thus, one assumes that the number of bits N storage in a sphere is proportional to the area A, that is

N = A l p 2 , (A1)

where l p is the Planck’s length

l p = G ℏ c 3 . (A2)

In addition, one assumes the following two basic conditions for small Δ x displacement and two conditions:

1) Force F, entropy S and temperature T relation: F Δ x = T Δ S .

2) Entropy S and Δ x relation: Δ S = 2π λ c Δ x , with λ c = ℏ m c being the Compton wave length.

3) The equipartition rule for the energy: E = 1 2 N T , with k B = 1 .

4) The rest mass equation expression: E = M c 2 .

From 1) and 2) one has

F = T Δ S Δ x = 4 π T λ c . (A3)

But from 3) one sees that

T = 2 E N , (A4)

and using 4), one finds that F becomes

F = 4 π λ c M c 2 N = 4 π ℏ M m c 3 N . (A5)

This expression can also be written as

F = 4 π G M m N ℏ G / c 3 = 4 π G M m N l p 2 . (A6)

Thus, using (A1) and (A2) one gets the formula

F = G M m r 2 , (A7)

which is the familiar Newton’s gravitatio law. Here, M denotes the mass enclosed by a spherical screen S^{2} (see also Ref. [