Is there anybody out there?
Abstract
The Fermi paradox is the discrepancy between the strong likelihood of alien intelligent life emerging (under a wide variety of assumptions) and the absence of any visible evidence for such emergence. We use this intriguing unlikeness to derive an upper limit on the fraction of living intelligent species that develop communication technology . indicates average over all the multiple manners civilizations can arise, grow, and develop such technology, starting at any time since the formation of our Galaxy in any location inside it. Following Drake, we factorize as the product of the fractions in which: (i) life arises, (ii) intelligence develops, and (iii) communication technology is developed. This averaging procedure must be regarded as a crude approximation because the characteristics of the initial conditions in a planet and its surroundings may affect the three phenomena with high complexity. In this approximation, the number of communicating intelligent civilizations that exist in the Galaxy at any given time is found to be , where is the average production rate of potentially habitable rocky planets with a longlasting () ecoshell and is the length of time that a typical civilization communicates. We estimate the production rate of exoplanets in the habitable zone and using recent determinations of the rate of gammaray bursts (GRBs) and their luminosity function, we calculate the probability that a lifethreatening (lethal) GRB could make a planet inhospitable to life, yielding . Our current measurement of then implies at the 95%C.L., where we have taken such that propagation distances of Galactic scales (), ensuring that any advanced civilization living in the Milky Way would be able to communicate with us.
Is there anybody out there?
\FullConference35th International Cosmic Ray Conference  ICRC217 
1020 July, 2017
Bexco, Busan, Korea
On May 25, 1961 President Kennedy’s announcement to put a man on the moon and bring him back safely before the end of the decade set the advent of human exploration of space for NASA, culminating to the landing on the Moon on July 16, 1969. It is difficult to believe that this is the only time such an event has ever happened in the history of the universe. On the other hand, if there is alien life capable of pulling off such a feat we must ask, as Fermi did, where is everybody? [1].
By adopting the starting point of a first approximation of the answer, we can write the number of intelligent civilizations in our galaxy at any given time capable of releasing detectable signals of their existence into space using a quite simple functional form,
(1) 
where is the average rate of star formation, is the fraction of stars with planetary systems, is the number of planets (per solar system) with a longlasting () ecoshell, is the fraction of suitable planets on which life actually appears, is the fraction of living species that develop intelligence, is the fraction of intelligent species with communications technology, and is the length of time such civilizations release detectable signals into space (i.e. the lifetime of the communicative phase) [2].
Following [3] we separate (1) into its astrophysical and biotechnological factors
(2) 
where represents the production rate of habitable planets with longlasting ecoshell (determined through astrophysics) and represents the product of all chemical, biological and technological factors leading to the development of a technological civilization. indicates average over all the multiple manners civilizations can arise, grow, and develop such technology, starting at any time since the formation of our Galaxy in any location inside it. This averaging procedure must be regarded as a crude approximation because the characteristics of the initial conditions in a planet and its surroundings may affect , , and with high complexity. In this work we estimate the production rate of exoplanets in the habitable zone and the rate of planetary catastrophes which could threaten the evolution of life on the surface of these worlds. Armed with these estimates we use our current measurement of to set an upper limit on .
The star formation rate in the Galaxy is estimated to be [4, 5]. This estimate has been derived assuming the Kroupa initial mass function (IMF) [6, 7]. The shape of this IMF is lognormallike and exhibits a peak around [8], suggesting there are roughly 2 stars per . Altogether, this yields . Now, only 10% of these stars are appropriate for harboring habitable planets. This is because the mass of the star to be sufficiently longlived (with main sequence lifetimes larger than 4.5 Gyr) and to possess circumstellar habitable zones outside the tidally locked region [9].
Next, in line with our stated plan, we estimate a rough probability that a habitable planet will survive and remain in a habitable zone to present day. Gammaray bursts (GRBs) are shortlived, luminous explosions, thought to originate from relativistic plasma launched at the deaths of massive stars. The widely accepted interpretation of GRB phenomenology is that the observable effects are due to the dissipation of the kinetic energy of a relativistically expanding fireball [15]. The physical conditions in the dissipation region produce a heavy flux of photons with energies above about 100 keV. It has been suggested that a nearby (Galactic) GRB may destroy the ozone layer, possibly making it damaging to life on Earth.
GRBs are generally divided in two groups according to their duration: long ) and short (). Short GRBs are weaker and hence their life threatening effect is negligible [22]. Throughout we consider only long GRBs. We want to estimate the expected number of GRBs that can terminate life in a planet that is situated at a Galactic radius . To this end we should take into account the following considerations:
 The luminosityrate function

, which measures the number of GRBs with a luminosity in a small range around a given value occurring per unit time and volume, is given by
(3) where , , , , and [23]. Here, is the volumetric rate of GRBs at . We consider a fiducial value of [23]. To accommodate the metallicity bias determined in [24] we follow [22] and take correction of a factor . It has been noted in [25] that such a low metallicity correction factor yields an upper limit on the volumetric rate of long GRBs . To derive our upper bound on we will adopt , since the larger the number of GRBs, the smaller the planets with longlasting ecoshell, and therefore the larger the value of . In general, any GRB regardless of its luminosity could terminate life if it is close enough to a planet. This is taken into account in our next point.
 The fluence

measures the amount of energy per unit area arriving to a planet from a distant GRB. If the distance between the GRB and the planet is , an isotropic emission and conservation of energy implies
(4) where is the total energy released by the GRB.
The effects that a copious flux of gamma rays may cause on the Earth atmosphere have been studied in [19, 20]. A fluence of could cause a depletion of roughly 68% of the ozone layer on a time scale of a month, whereas fluences of and could cause depletions of about 91% and 98%, respectively. This implies that a fluence of could cause some damage to life, while will wipe out nearly the whole atmosphere causing a catastrophic life extinction event. Note, however, that the complete removal of all life on an Earthlike planet is a very unlikely event [26]. The critical fluence gives the limit on acceptable fluence on a planet. Following [22], in our calculations we take as our fiducial life threatening critical fluence.
 The hazard distance

of a GRB, characterized by its total energy , measures the length for which the fluence is higher than the critical fluence . Any planet within this distance will have life terminated. For a fixed GRB energy and critical fluence, we use (4) to obtain the hazard distance
(5) where in the second rendition we have assumed that the total GRB energy is roughly the average ( half) of the peak flux [22]. A good but rough estimate of follows from the assumption that the typical duration of long GRBs is .
 The fraction of hazardous galaxy

measures the fraction of the total galactic mass that is within the hazard distance , for any point at radius from the Galactic center. The fraction of galactic mass that is contained within a surface element at radius is given by,
(6) such that , with [27]. As illustrated in Fig. 1, to calculate the fraction of hazardous galaxy around a point with radius , one has to integrate (6) in a circular region of radius centred at ,
(7) Defining the unit vectors and perpendicular to , we can write and . Since ,
(8) Defining , we have , where runs from to , so that
(9)  The rate of lifethreatening GRBs
Taking into account all of these considerations we now turn to estimate the expected number of GRBs that can terminate life in a planet that is situated at a Galactic radius . As it turns out is . This means that for time scales we expect to have (on average) a small number of GRB events. This kind of estimate is then well suited to Poisson statistics. The Poisson distribution is a discrete probability distribution for the counts of events that occur randomly in a given interval of time (or space). Of particular interest here, the probability of having GRBs when the expected average is is given by . The probability of having or more GRBs is . The average number of GRBs during a time is , and therefore
(11) 
In Fig. 1 we show probability contours of at least one GRB having occurred in the past time with enough flux to produce significant life extinction on Earth. Since the probabilities get too high in most of the parameter space of interest, the legend does not show the probability (which would be too close to ), but instead the parameter defined by
(12) 
(12) gives the usual correspondence between probabilities and the standard deviation for a normal distribution, which leads to . Some value correspondences between and the probabilities in are , , and . The circles indicate selected values of and . Our estimates are in good agreement with those given in Table II of [22]. These findings seem to indicate that a nearby GRB may have caused one of the five greatest mass extinctions on Earth. However, it is important to remind the reader that the estimates shown in Fig. 1 rely on an upper limit of the volumetric rate of long GRBs; namely, .
Coming back to our calculation, we can use (11) to compute the probability of having at least one lethal GRB, for a critical fluency and a fixed lookback time , as a function of the distance . Life has been evolving on Earth for close to 4 Gyr [30, 31], but complex life is well under 1 Gyr old, and intelligent life is only a Myr old at most. In what follows we adopt and 4 Gyr as critical time intervals for life evolution [32]. In Fig. 2 we show the probability of having one or more lethal GRBs for a critical fluency and and 4, as a function of distance to the Galactic center.
All we need to do now is add the components together to arrive at the production rate of habitable planets with a longlasting ecoshell,
(13) 
It is of interest to study how depends on the different parameters involved. According to (3), (10) and (11), the dependence of in and can be grouped in the same functional dependence. Introducing , one can rewrite , where and . We can then rewrite (13) as , where
(14) 
The value of is shown in Fig. 3 as a function of for different values of . Note that for , with and , (13) leads to and , respectively.
Finally, to determine the upper bound on we must decide on the possible minimum . Herein we consider such that propagation distances of Galactic scales (). This would provide enough time to receive electromagnetic (and/or highenergy neutrino [33]) signals from any advanced civilization living in the Milky Way which is trying to communicate with us.
As of today, the nonobservation of evidence of advanced civilizations implies that models of predicting are excluded at the 95% C.L. [34].
In summary, in this paper we have derived an upper bound on the average fraction of living intelligent species that develop communication technology: at the 95% C.L. Future observations
could help to tighten this bound. In particular, a new arsenal of data will certainly provide an ideal testing ground to improve our understanding about: (i) the occurrence of exoplanets in the habitable zone, (ii) the early star formation rate models, and (iii) the GRB phenomenology. The past few years have witnessed the discovery of more and more rocky planets that are larger and heftier than Earth. Finding the Earthtwins is a higher order challenge, because these smaller planets produce fainter signals and hence only a few have been discovered. Technology to detect and image Earthlike planets has been developed for use of the next generation space telescopes. The Transiting Exoplanet Survey Satellite (TESS) is NASA’s next step in the search for planets outside of our solar system, including those that could support life. The NASA roadmap will subsequently continue with the launch of the James Webb Space Telescope (JWST) [37] and perhaps the proposed Wide Field Infrared Survey Telescope  Astrophysics Focused Telescope Assets (WFIRSTAFTA) early in the next decade [38]. The ability to detect alien life may still be years or more away, but the quest is underway.
We thank Michael Unger for discussion and Pink Floyd for inspiration. This research was supported by the U.S. NSF (Grant PHY1620661) and NASA (Grant NNX13AH52G).
Footnotes
 The habitable zone is the orbital range around a star within which surface liquid water could be sustained. Since water is essential for life as we know it, the search for biosignature gases naturally focuses on planets located in the habitable zone of their host stars. The habitable zone of the solar system looks like a ring around the Sun. Rocky planets with an orbit within this ring may have liquid water to support life. The habitable zone around a single star looks similar to the habitable zone in our Solar System. The only difference is the size of the ring. If the star is bigger than the Sun it has a wider zone, if the star is smaller it has a narrower zone. It might seem that the bigger the star the better. However, the biggest stars have relatively short lifespans, so the life around them probably would not have enough time to evolve. The habitable zones of small stars face a different problem. Besides being narrow they are relatively close to the star. A hypothetical planet in such a region would be tidally locked. That means that one half of it would always face the star and be extremely hot, while the opposite side would always be facing away and freezing. Such conditions are not very favorable for life.
 Ozone (O) forms a kind of layer in the stratosphere, which normally prevents about 90% of the solar UVB () radiation from reaching the Earth’s surface. UVB is extremely damaging to most organisms, particularly since it easily damages DNA and proteins.
 If a corresponding hypothesis test is performed, the confidence level (C.L.) is the complement of the level of statistical significance, e.g, a 95% confidence interval reflects a significance level of 0.05.
References
 E. M. Jones, Where is everybody?, An account of Fermi’s Question, LA10311MS, UC34b, CIC14 Report Collection, March 1985; Physics Today 38, p.11 August 1985, (without full correspondence). http://www.fas.org/sgp/othergov/doe/lanl/la10311ms.pdf
 F. D. Drake, The radio search for intelligent extraterrestrial life, in Current aspects of exobiology, (Eds. G. Mamikunian and M. H. Briggs, Elsevier Ltd., 1965), p.323. ISBN: 9781483200477
 N. Prantzos, A joint analysis of the Drake equation and the Fermi paradox, Int. J. Astrobiol. 12, 246 (2013) doi:10.1017/S1473550413000037 [arXiv:1301.6411].
 T. C. Licquia and J. A. Newman, Improved estimates of the Milky Way’s stellar mass and star formation rate from hierarchical bayesian metaanalysis, Astrophys. J. 806, 96 (2015) doi:10.1088/0004637X/806/1/96 [arXiv:1407.1078].
 L. Chomiuk and M. S. Povich, Toward a unification of star formation rate determinations in the Milky Way and other galaxies, Astron. J. 142, 197 (2011) doi:10.1088/00046256/142/6/197 [arXiv:1110.4105].
 P. Kroupa, The initial mass function of stars: Evidence for uniformity in variable systems, Science 295, 82 (2002) doi:10.1126/science.1067524 [astroph/0201098].
 P. Kroupa and C. Weidner, Galactic  field IMFs of massive stars, Astrophys. J. 598, 1076 (2003) doi:10.1086/379105 [astroph/0308356].
 C. J. Lada, Star formation in the local Milky Way, arXiv:1508.02711.
 S. Seager, Exoplanet habitability, Science, 340, 577 (2013) doi: 10.1126/science.1232226
 W. J. Borucki et al. [Kepler Collaboration], Kepler planetdetection mission: Introduction and first results, Science 327, 977 (2010). doi:10.1126/science.1185402
 W. J. Borucki et al. [Kepler Collaboration], Characteristics of planetary candidates observed by Kepler, II: Analysis of the first four months of data, Astrophys. J. 736, 19 (2011) doi:10.1088/0004637X/736/1/19 [arXiv:1102.0541].
 N. M. Batalha et al. [Kepler Collaboration], Planetary candidates observed by Kepler, III: Analysis of the first 16 months of data, Astrophys. J. Supp. 204, 24 (2013) doi:10.1088/00670049/204/2/24 [arXiv:1202.5852].
 C. D. Dressing and D. Charbonneau, The occurrence rate of small planets around small stars, Astrophys. J. 767, 95 (2013) doi:10.1088/0004637X/767/1/95 [arXiv:1302.1647].
 R. k. Kopparapu, A revised estimate of the occurrence rate of terrestrial planets in the habitable zones around Kepler Mdwarfs, Astrophys. J. 767, L8 (2013) doi:10.1088/20418205/767/1/L8 [arXiv:1303.2649].
 T. Piran, Gammaray bursts and the fireball model, Phys. Rept. 314, 575 (1999) doi:10.1016/S03701573(98)001276 [astroph/9810256].
 M. A. Ruderman, Possible consequences of nearby supernova explosions for atmospheric ozone and terrestrial life, Science 185, 1079 (1974).
 S. E. Thorsett, Terrestrial implications of cosmological gammaray burst models, Astrophys. J. 444, L53 (1995) doi:10.1086/187858 [astroph/9501019].
 A. Dar, A. Laor and N. J. Shaviv, Life extinctions by cosmic ray jets, Phys. Rev. Lett. 80, 5813 (1998) doi:10.1103/PhysRevLett.80.5813 [astroph/9705008].
 B. C. Thomas, C. H. Jackman, A. L. Melott, C. M. Laird, R. S. Stolarski, N. Gehrels, J. K. Cannizzo and D. P. Hogan, Terrestrial ozone depletion due to a Milky Way gammaray burst, Astrophys. J. 622, L153 (2005) doi:10.1086/429799 [astroph/0411284].
 B. C. Thomas et al., Gammaray bursts and the Earth: Exploration of atmospheric, biological, climatic and biogeochemical effects, Astrophys. J. 634, 509 (2005) doi:10.1086/496914 [astroph/0505472].
 A. Melott et al., Did a gammaray burst initiate the late Ordovician mass extinction?, Int. J. Astrobiol. 3, 55 (2004) doi:10.1017/S1473550404001910 [astroph/0309415].
 T. Piran and R. Jimenez, Possible role of gammaray bursts on life extinction in the Universe, Phys. Rev. Lett. 113, 231102 (2014) doi:10.1103/PhysRevLett.113.231102 [arXiv:1409.2506].
 D. Wanderman and T. Piran, The luminosity function and the rate of Swift’s gammaray bursts, Mon. Not. Roy. Astron. Soc. 406, 1944 (2010) doi:10.1111/j.13652966.2010.16787.x [arXiv:0912.0709].
 R. Jimenez and T. Piran, Reconciling the gammaray burst rate and star formation histories, Astrophys. J. 773, 126 (2013) doi:10.1088/0004637X/773/2/126 [arXiv:1303.4809].
 M. G. Gowanlock, Astrobiological effects of gammaray bursts in the Milky Way galaxy, Astrophys. J. 832, 38 (2016) doi:10.3847/0004637X/832/1/38 [arXiv:1609.09355].
 D. Sloan, R. Alves Batista, and A. Loeb, The resilience of life to astrophysical events, Scientific Reports (to be published).
 J. Bovy and H. W. Rix, A direct dynamical measurement of the Milky Way’s disk surface density profile, disk scale length, and dark matter profile at 4 kpc R 9 kpc, Astrophys. J. 779, 115 (2013) doi:10.1088/0004637X/779/2/115 [arXiv:1309.0809].
 Y. Li and B. Zhang, Can life survive gammaray bursts in the highredshift universe?, Astrophys. J. 810, 41 (2015) doi:10.1088/0004637X/810/1/41 [arXiv:1507.05966 [astroph.HE]].
 P. J. McMillan, Mass models of the Milky Way, Mon. Not. Roy. Astron. Soc. 414, 2446 (2011) doi:10.1111/j.13652966.2011.18564.x [arXiv:1102.4340 [astroph.GA]].
 S. J. Mojzsis, G. Arrhenius, K. D. McKeegan, T. M. Harrison, A. P. Nutman, and C. R. L. Friend, Evidence for life on Earth before 3,800 million years ago, Nature 384, 55 (1996) doi:10.1038/384055a0.
 M. S. Dodd, D. Papineau, T. Grenne, J. F. Slack, M. Rittner, F. Pirajno, J. O’Neil, and C. T. S. Little, Evidence for early life in Earth’s oldest hydrothermal vent precipitates, Nature 543, 60 (2017) doi:10.1038/nature21377.
 A. Loeb, R. A. Batista and D. Sloan, Relative likelihood for life as a function of cosmic time, JCAP 1608, 040 (2016) doi:10.1088/14757516/2016/08/040 [arXiv:1606.08448 [astroph.CO]].
 J. G. Learned, S. Pakvasa and A. Zee, Galactic neutrino communication, Phys. Lett. B 671, 15 (2009) doi:10.1016/j.physletb.2008.11.057 [arXiv:0805.2429].
 G. J. Feldman and R. D. Cousins, A Unified approach to the classical statistical analysis of small signals, Phys. Rev. D 57, 3873 (1998) doi:10.1103/PhysRevD.57.3873 [physics/9711021].
 R. C. Kennicutt, Jr. and N. J. Evans II, Star formation in the Milky Way and nearby galaxies, Ann. Rev. Astron. Astrophys. 50, 531 (2012) doi:10.1146/annurevastro081811125610 [arXiv:1204.3552].
 https://tess.gsfc.nasa.gov
 https://www.jwst.nasa.gov
 https://wfirst.gsfc.nasa.gov