The Event Horizon Telescope Consortium is on the verge to provide the first telescopic image of massive black holes SgrA* and M87* surrounded by accretion disks, at a resolution scale comparable to the size of their event horizons. Well before this remarkable achievement made possible by VLBI radio astronomy, many researchers used the computer to reconstruct how a black hole surrounded by luminous material would look from close-up views. The images must experience extraordinary optical deformations due to the deflection of light rays produced by the strong curvature of the space-time in the vicinity. General relativity allows the calculation of such effects, both on a surrounding accretion disk and on the background star field. This article is an exhaustive and illustrated review of the numerical work on black hole imaging done during the first thirty years of its history.

Astronomical observations are about to deliver the very first telescopic image of the massive black hole lurking at the Galactic Center. The mass of data collected in one night by the Event Horizon Telescope network, exceeding everything that has ever been done in any scientific field, should provide a recomposed image during 2018. All this, forty years after the first numerical simulations done by the present author.

In the last decade, the study of the overall shape of the universe, called Cosmic Topology, has become testable by astronomical observations, especially the data from the Cosmic Microwave Background (hereafter CMB) obtained by WMAP and Planck telescopes. Cosmic Topology involves both global topological features and more local geometrical properties such as curvature. It deals with questions such as whether space is finite or infinite, simply-connected or multi-connected, and smaller or greater than its observable counterpart. A striking feature of some relativistic, multi-connected small universe models is to create multiples images of faraway cosmic sources. While the last CMB (Planck) data fit well the simplest model of a zero-curvature, infinite space model, they remain consistent with more complex shapes such as the spherical Poincare Dodecahedral Space, the flat hypertorus or the hyperbolic Picard horn. We review the theoretical and observational status of the field.

I give a critical review of the holographic hypothesis, which posits that a universe with gravity can be described by a quantum field theory in fewer dimensions. I first recall how the idea originated from considerations on black hole thermodynamics and the so-called information paradox that arises when Hawking radiation is taken into account. String Quantum Gravity tried to solve the puzzle using the AdS/CFT correspondence, according to which a black hole in a 5-D anti-de Sitter space is like a flat 4-D field of particles and radiation. Although such an interesting holographic property, also called gauge/gravity duality, has never been proved rigorously, it has impulsed a number of research programs in fields as diverse as nuclear physics, condensed matter physics, general relativity and cosmology. I finally discuss the pros and cons of the holographic conjecture, and emphasizes the key role played by black holes for understanding quantum gravity and possible dualities between distant fields of theoretical physics.

What is the shape of the Universe? Is it finite or infinite ? Is space multi-connected to create ghost images of faraway cosmic sources? After a "dark age" period, the field of cosmic topology has now become one of the major concerns in astronomy and cosmology, not only from theorists but also from observational astronomers. Here I give a personal account of the spectacular progress in the field since the beginning of the 1990's, when I started to work in it.

I discuss some aspects of the use of computers in Relativity, Astrophysics and Cosmology. For each section I provide two examples representative of the field, including gravitational collapse, black hole imagery, supernovae explosions, star-black hole tidal interactions, N-body cosmological simulations and detection of cosmic topology.

This is an editorial note to accompany printing as a Golden Oldie in the Journal of General Relativity and Gravitation of the fundamental article by Georges Lema\^itre first published in French in 1927, in which the author provided the first explanation of the observations on the recession velocities of galaxies as a natural consequence of dynamical cosmological solutions of Einstein's field equations, and discovered the so-called Hubble law. We analyze in detail the scientific contents of this outstanding work, we describe how it remained unread or poorly appreciated until 1930, and we list and explain the differences between the 1927 and 1931 versions. Indeed the English translation published in 1931 in MNRAS was not perfectly faithful to the original text - it was updated. As it turned out very recently, the updates were done by Lema\^itre himself, but the discrepancies between the two texts caused a temporary stir among historians. Our new translation - given in the Appendix - follows the 1927 version exactly.

This is an editorial note to accompany reprinting as a Golden Oldie in the Journal of General Relativity and Gravitation of the famous note by Georges Lemaitre on the quantum birth of the universe, published in Nature in 1931. We explain why this short (457 words) article can be considered to be the true "Charter" of the modern Big Bang theory.

Aims: We investigate the stellar pancake mechanism during which a solar-type star is tidally flattened within its orbital plane passing close to a 10^6 solar masses black hole. We simulate the relativistic orthogonal compression process and follow the associated shock waves formation. Methods: We consider a one-dimensional hydrodynamical stellar model moving in the relativistic gravitational field of a non-rotating black hole. The model is numerically solved using a Godunov-type shock-capturing source-splitting method in order to correctly reproduce the shock waves profiles. Results: Simulations confirm that the space-time curvature can induce several successive orthogonal compressions of the star which give rise to several strong shock waves. The shock waves finally escape from the star and repeatedly heat up the stellar surface to high energy values. Such a shock-heating could interestingly provide a direct observational signature of strongly disruptive star - black hole encounters through the emission of hard X or soft gamma-ray bursts. Timescales and energies of such a process are consistent with some observed events such as GRB 970815.

In 2007 an international conference engaged a reflection on the present conditions for sciences and scientific practice, to initiate a dialogue between science, philosophy, epistemology and sociology. Starting from many examples taken from the history of sciences and the present status of astrophysics, cosmology and fundamental physics, I analyze various issues such as the possible limits of theoretical and technological advances, the paradox of non-verifiability in string theory, incompleteness, the role of imagination in scientific research along with educational aspects, the epistemological value of the Ockham's razor principle for constraining theories, and our undersanding of man's place in Creation.

What is the shape of the Universe? Is it curved or flat, finite or infinite ? Is space "wrapped around" to create ghost images of faraway cosmic sources? We review how tessellations allow to build multiply-connected 3D Riemannian spaces useful for cosmology. We discuss more particularly the proposal of a finite, positively curved, dodecahedral space for explaining some puzzling features of the cosmic microwave background radiation, as revealed by the 2003-2006 WMAP data releases.

We interest in the case of a main-sequence star deeply penetrating within the tidal radius of a massive black hole. We focus on the compression phase leading to a so-called pancake configuration of the star at the instant of maximal compression. The aim is to study the tidal compression process paying particular attention to the development of shock waves;to deduce reliable estimates of the thermodynamical quantities involved in the pancake star; and to solve a controversy about whether or not thermonuclear reactions can be triggered in the core of a tidally compressed star. We have set up a one-dimensional hydrodynamical model well-adapted to the geometry of the problem. Based on the high-resolution shock-capturing Godunov-type approach, it allows to study the compression phase undergone by the star in the direction orthogonal to its orbital plane. We show the existence of two regimes depending on whether shock waves develop before or after the instant of maximal core compression. In both cases we confirm high compression and heating factors in the stellar core able to trigger a thermonuclear explosion. Moreover, we show that the shock waves carry outwards a brief but very high peak of temperature from the centre to the surface of the star. We tentatively conclude that the phenomenon could give rise to hard electromagnetic radiation, to be compared to some X-ray flares already observed in some galactic nuclei harbouring massive black holes. Finally, we estimate that the rate of pancake stars should be about $10^{-5}$ per galaxy per year. If generated in hard X- or $\gamma$-ray band, several events of this kind per year should be detectable within the full observable universe.

The full three-year WMAP results (WMAP3) reinforce the absence of large-angle correlations at scales greater than 60 degrees. The Poincare dodecahedral space (PDS) model model, which may naturally explain such features, thus remains a plausible cosmological model, despite recent controversy about whether matched circle searches would or would not push the topology beyond the horizon. Here, we have used new eigenmode calculations of the dodecahedral space to predict the cosmic microwave background (CMB) temperature fluctuations in such models, with an improved angular resolution. We have simulated CMB maps and confirmed the expected presence of matching circles. For a set of plausible cosmological parameters, we have derived the angular power spectrum of the CMB up to large wavenumbers. Comparison with the WMAP3 observations confirms a remarkable fit with a PDS model, for a value $\Omega_0 = 1.018$ of the average total energy density.

We provide an epistemological analysis of the developments of relativistic cosmology from 1917 to 2006, based on the seminal articles by Einstein, de Sitter, Friedmann, Lemaitre, Hubble, Gamow and other main historical figures of the field. It appears that most of the ingredients of the present-day standard cosmological model, such as the accelation of the expansion due to a repulsive dark energy, the interpretation of the cosmological constant as vacuum energy or the possible non-trivial topology of space, had been anticipated by Lemaitre, although his papers remain desperately unquoted.

General relativity does not allow one to specify the topology of space, leaving the possibility that space is multiply rather than simply connected. We review the main mathematical properties of multiply connected spaces, and the different tools to classify them and to analyse their properties. Following their mathematical classification, we describe the different possible muticonnected spaces which may be used to construct Friedmann-Lemaitre universe models. Observational tests concern the distribution of images of discrete cosmic objects or more global effects, mainly those concerning the Cosmic Microwave Background. According to the 2003-2006 WMAP data releases, various deviations from the flat infinite universe model predictions hint at a possible non-trivial topology for the shape of space. In particular, a finite universe with the topology of the Poincar\'e dodecahedral spherical space fits remarkably well the data and is a good candidate for explaining both the local curvature of space and the large angle anomalies in the temperature power spectrum. Such a model of a small universe, whose volume would represent only about 80% the volume of the observable universe, offers an observational signature in the form of a predictable topological lens effect on one hand, and rises new issues on the physics of the early universe on the other hand.

Multiply connected space sections of the universe on a scale smaller than the horizon size can leave an imprint on cosmic microwave background polarization maps, in such a way that the so-called ``circles-in-the-sky'' method can be used to detect or constrain the topology. We investigate some specific cases, namely toroidal and sixth-turn spaces, in order to show the influence of topology on CMB polarization. The correlation between matched points happens to be always positive and higher than 75% regardless of the angular scale and of the cosmological parameters, except for reionization. This figure is better than what occurs in temperature maps, but is achieved only in the absence of noise. It is only slightly reduced by the filtering scheme.

Conventional thinking says the universe is infinite. But it could be finite and relatively small, merely giving the illusion of a greater one, like a hall of mirrors. Recent astronomical measurements add support to a finite space with a dodecahedral topology.

What is the shape of space is a long-standing question in cosmology. In this talk I review recent advances in cosmic topology since it has entered a new era of experimental tests. High redshift surveys of astronomical sources and accurate maps of the Cosmic Microwave Background radiation (CMB) are beginning to hint at the shape of the universe, or at least to limit the wide range of possibilities. Among those possibilites are surprising "wrap around" universe models in which space, whatever its curvature, may be smaller than the observable universe and generate topological lensing effects on a detectable cosmic scale. In particular, the recent analysis of CMB data provided by the WMAP satellite suggests a finite universe with the topology of the Poincare dodecahedral spherical space. Such a model of a "small universe", the volume of which would represent only about 80% the volume of the observable universe, offers an observational signature in the form of a predictable topological lens effect on one hand, and rises new issues on the early universe physics on the other hand.

A widespread myth asserts that all small universe models suppress the CMB quadrupole. In actual fact, some models suppress the quadrupole while others elevate it, according to whether their low-order modes are weak or strong relative to their high-order modes. Elementary geometrical reasoning shows that a model's largest dimension determines the rough value ell_min at which the CMB power spectrum ell(ell + 1)C_ell/(2pi) effectively begins; for cosmologically relevant models, ell_min < 4. More surprisingly, elementary geometrical reasoning shows that further reduction of a model's smaller dimensions -- with its largest dimension held fixed -- serves to elevate modes in the neighborhood of ell_min relative to the high-ell portion of the spectrum, rather than suppressing them as one might naively expect. Thus among the models whose largest dimension is comparable to or less than the horizon diameter, the low-order C_ell tend to be relatively weak in well-proportioned spaces (spaces whose dimensions are approximately equal in all directions) but relatively strong in oddly-proportioned spaces (spaces that are significantly longer in some directions and shorter in others). We illustrate this principle in detail for the special cases of rectangular 3-tori and spherical spaces. We conclude that well-proportioned spaces make the best candidates for a topological explanation of the low CMB quadrupole observed by COBE and WMAP.

Starting from an anisotropic flat cosmological model(Bianchi type $I$), we show that conditions leading to isotropisation fall into 3 classes, respectively 1, 2, 3. We look for necessary conditions such that a Bianchi type $I$ model reaches a stable isotropic state due to the presence of several massive scalar fields minimally coupled to the metric with a perfect fluid for class 1 isotropisation. The conditions are written in terms of some functions $\ell$ of the scalar fields. Two types of theories are studied. The first one deals with scalar tensor theories resulting from extra-dimensions compactification, where the Brans-Dicke coupling functions only depend on their associated scalar fields. The second one is related to the presence of complex scalar fields. We give the metric and potential asymptotical behaviours originating from class 1 isotropisation. The results depend on the domination of the scalar field potential compared to the perfect fluid energy density. We give explicit examples showing that some hybrid inflation theories do not lead to isotropy contrary to some high-order theories, whereas the most common forms of complex scalar fields undergo a class 3 isotropisation, characterised by strong oscillations of the $\ell$ functions.

This article investigates the signature of the seventeen multi-connected flat spaces in cosmic microwave background (CMB) maps. For each such space it recalls a fundamental domain and a set of generating matrices, and then goes on to find an orthonormal basis for the set of eigenmodes of the Laplace operator on that space. The basis eigenmodes are expressed as linear combinations of eigenmodes of the simply connected Euclidean space. A preceding work, which provides a general method for implementing multi-connected topologies in standard CMB codes, is then applied to simulate CMB maps and angular power spectra for each space. Unlike in the 3-torus, the results in most multi-connected flat spaces depend on the location of the observer. This effect is discussed in detail. In particular, it is shown that the correlated circles on a CMB map are generically not back-to-back, so that negative search of back-to-back circles in the WMAP data does not exclude a vast majority of flat or nearly flat topologies.

Cosmology's standard model posits an infinite flat universe forever expanding under the pressure of dark energy. First-year data from the Wilkinson Microwave Anisotropy Probe (WMAP) confirm this model to spectacular precision on all but the largest scales (Bennett {\it et al.}, 2003 ; Spergel {\it et al.}, 2003). Temperature correlations across the microwave sky match expectations on scales narrower than $60^{\circ}$, yet vanish on scales wider than $60^{\circ}$. Researchers are now seeking an explanation of the missing wide-angle correlations (Contaldi {\it et al.}, 2003 ; Cline {\it et al.}, 2003). One natural approach questions the underlying geometry of space, namely its curvature (Efstathiou, 2003) and its topology (Tegmark {\it et al.}, 2003). In an infinite flat space, waves from the big bang would fill the universe on all length scales. The observed lack of temperature correlations on scales beyond $60^{\circ}$ means the broadest waves are missing, perhaps because space itself is not big enough to support them. Here we present a simple geometrical model of a finite, positively curved space -- the Poincar\'e dodecahedral space -- which accounts for WMAP's observations with no fine-tuning required. Circle searching (Cornish, Spergel and Starkman, 1998) may confirm the model's topological predictions, while upcoming Planck Surveyor data may confirm its predicted density of $\Omega_0 \simeq 1.013 > 1$. If confirmed, the model will answer the ancient question of whether space is finite or infinite, while retaining the standard Friedmann-Lema\^\i{}tre foundation for local physics.

General relativity does not allow one to specify the topology of space, leaving the possibility that space is multi-- rather than simply--connected. We review the main mathematical properties of multi--connected spaces, and the different tools to classify them and to analyse their properties. Following the mathematical classification, we describe the different possible muticonnected spaces which may be used to construct universe models. We briefly discuss some implications of multi--connectedness for quantum cosmology, and its consequences concerning quantum field theory in the early universe. We consider in details the properties of the cosmological models where space is multi--connected, with emphasis towards observable effects. We then review the analyses of observational results obtained in this context, to search for a possible signature of multi--connectedness, or to constrain the models. They may concern the distribution of images of cosmic objects like galaxies, clusters, quasars,..., or more global effects, mainly those concerning the Cosmic Microwave Background, and the present limits resulting from them.

This article investigates the computation of the eigenmodes of the Laplacian operator in multi-connected three-dimensional spherical spaces. General mathematical results and analytical solutions for lens and prism spaces are presented. Three complementary numerical methods are developed and compared with our analytic results and previous investigations. The cosmological applications of these results are discussed, focusing on the cosmic microwave background (CMB) anisotropies. In particular, whereas in the Euclidean case too small universes are excluded by present CMB data, in the spherical case there will always exist candidate topologies even if the total energy density parameter of the universe is very close to unity.

This article gives the construction and complete classification of all three-dimensional spherical manifolds, and orders them by decreasing volume, in the context of multiconnected universe models with positive spatial curvature. It discusses which spherical topologies are likely to be detectable by crystallographic methods using three-dimensional catalogs of cosmic objects. The expected form of the pair separation histogram is predicted (including the location and height of the spikes) and is compared to computer simulations, showing that this method is stable with respect to observational uncertainties and is well suited for detecting spherical topologies.

We investigate to what extent the cosmic crystallographic methods aimed to detect the topology of the universe using catalogues of cosmic objects would be damaged by various observational uncertainties. We find that the topological signature is robust in the case of Euclidean spaces, but is very fragile in the case of compact hyperbolic spaces. Comparing our results to the presently accepted range of values for the curvature parameters, the best hopes for detecting space topology rest on elliptical space models.

Multi-connected universe models with space identification scales smaller than the size of the observable universe produce topological images of cosmic sources. We generalise to locally hyperbolic spaces the crystallographic method, aimed to detect the topology from three-dimensional catalogs of cosmological objects. Our new method is based on the construction of a {\it collecting-correlated-pair} technique which enhances the topological signature and can make it detectable. The main idea is that in multi-connected universes, equal distances occur more often than by chance. We present an idealised version of this method as well as numerical simulations, we discuss the statistical relevance of the expected signature and we show how the extraction of a topological signal may also lead to a precise determination of the cosmological parameters. A more realistic version of the method which takes account of uncertainties in the position and velocity data is then discussed. We apply our technique to presently available data, namely a quasar catalog, and discuss the significance of the result. We show how the improved crystallographic method both competes with the two--dimensional methods based on cosmic microwave background analyses, and suffers from the same drawback: the shortage of present observational data.

``One could imagine that as a result of enormously extended astronomical experience, the entire universe consists of countless identical copies of our Milky Way, that the infinite space can be partitioned into cubes each containing an exactly identical copy of our Milky Way. Would we really cling on to the assumption of infinitely many identical repetitions of the same world? >... We would be much happier with the view that these repetitions are illusory, that in reality space has peculiar connection properties so that if we leave any one cube through a side, then we immediately reenter it through the opposite side.'' (Schwarzschild 1900, translation 1998) Developments in the theoretical and observational sides of cosmic topology were slow for most of the century, but are now progressing rapidly, at the scale of most interest which is 1-10 h^-1 Gpc rather than 10kpc. The historical, mathematical and observational sides of this subject are briefly reviewed in this course.

If the assumption that physical space has a trivial topology is dropped, then the Universe may be described by a multiply connected Friedmann-Lema\^{\i}tre model on a sub-horizon scale. Specific candidates for the multiply connected space manifold have already been suggested. How precisely would a significant detection of multiple topological images of a single object, or a region on the cosmic microwave background, (due to photons arriving at the observer by multiple paths which have crossed the Universe in different directions), constrain the values of the curvature parameters $\Omega_0$ and $\lambda_0$? The way that the constraints on $\Omega_0$ and $\lambda_0$ depend on the redshifts of multiple topological images and on their radial and tangential separations is presented and calculated. The tangential separations give the tighter constraints: multiple topological images of known types of astrophysical objects at redshifts $z \ltapprox 3$ would imply values of $\Omega_0$ and $\lambda_0$ preciser than $\sim 1%$ and $\sim 10%$ respectively. Cosmic microwave background `spots' identified with lower redshift objects by the Planck or MAP satellites would provide similar precision. This method is purely geometrical: no dynamical assumptions (such as the virial theorem) are required and the constraints are independent of the Hubble constant, $H_0.$

Universe models with compact spatial sections smaller than the observable universe produce a topological lens effect. Given a catalog of cosmic sources, we estimate the number of topological images in locally hyperbolic and locally elliptic spaces, as a function of the cosmological parameters, of the volume of the spatial sections and of the catalog depth. Next we apply the crystallographic method, aimed to detect a topological signal in the 3D distance histogram between images, to compact hyperbolic models. Numerical calculations in the Weeks manifold allows us to check the absence of crystallographic signature of topology, due to the fact that the number of copies of the fundamental domain in the observable covering space is low and that the points are not moved the same distance by the holonomies of space.

In the first part I set out some unexplored historical material about the early development of cosmic topology. In the second part I briefly comment new developments in the field since the Lachieze-Rey & Luminet report (1995), both from a theoretical and an observational point of view.

Our understanding of space and time is probed to its depths by black holes. These objects, which appear as a natural consequence of general relativity, provide a powerful analytical tool able to examine macroscopic and microscopic properties of the universe. This introductory article presents in a pictorial way the basic concepts of black hole's theory, as well as a description of the astronomical sites where black holes are suspected to lie, namely binary X-ray sources and galactic nuclei.

We assume that the Universe has a non trivial topology whose compact spatial sections have a volume significantly smaller than the horizon volume. By a topological lens effect, such a "folded" space configuration generates multiple images of cosmic sources, e.g. clusters of galaxies. We present a simple and powerful method to unveil non ambiguous observational effects, independently of the sign of the curvature and of the topological type.