In previous papers it was shown that in a quasi-spherical Szekeres (QSS) metric, impulses of gamma radiation can arise that have several properties in common with the observed gamma-ray bursts. This happens when the bang-time function $t_B(r)$ has a gate-shaped hump around the origin of the QSS region. The gamma rays arise along two preferred directions of the QSS geometry (coincident with dipole extrema when axially symmetric, otherwise unrelated). In these directions, the rays of the relic radiation are blueshifted rather than redshifted. The blueshift is generated in a thin region between the Big Bang (BB) and the extremum-redshift hypersurface (ERH). But the Szekeres models can describe the real Universe only forward in time from the last-scattering hypersurface (LSH) because the matter in them has zero pressure. The ERH is tangent to the BB at the origin, so in a neighbourhood thereof the ERH lies earlier than the LSH and no blueshift is generated in the physical region. The question thus arose whether the BB and ERH can be ``unglued'' if the QSS region has no origin, but the areal radius function $\Phi$ has a local maximum or minimum somewhere. In the present paper it is demonstrated that this is indeed the case. If the hump in $t_B(r)$ is centered around the minimum of $\Phi$, then the BB and ERH in general do not coincide there and a stronger blueshift is generated on rays passing nearby. It follows that a lower and narrower hump on the BB set can generate sufficient blueshift to move the initial frequencies of the relic radiation to the gamma range. These facts are demonstrated by numerical calculations in an explicit example of a QSS region.

This paper is a follow-up on two previous ones, in which properties of blueshifted rays were investigated in Lema\^{\i}tre - Tolman (L-T) and quasispherical Szekeres (QSS) spacetimes. In the present paper, an axially symmetric QSS deformation is superposed on such an L-T background that was proved, in the first paper, to mimic several properties of gamma-ray bursts. The present model makes $z$ closer to $-1$ than in the background L-T spacetime, and, as implied by the second paper, strong blueshifts exist in it only along two opposite directions. The QSS region is matched into a Friedmann background. The Big Bang (BB) function $t_B(r)$, which is constant in the Friedmann region, has a gate-shaped hump in the QSS region. Since a QSS island generates stronger blueshifts than an L-T island, the BB hump can be made lower - then it is further removed from the observer and implies a smaller observed angular radius of the source. Consequently, more sources can be fitted into the sky - all these facts are confirmed by numerical computations. Null geodesics reaching present observers from different directions relative to the BB hump are numerically calculated. Patterns of redshift across the image of the source and along the rays are displayed.

In Lema\^{\i}tre -- Tolman (L--T) models, light rays emitted \textit{radially} at the Big Bang (BB) at such radial coordinates $r$ where the bang-time function $t_B(r)$ has $\dril {t_B} r \neq 0$ reach every observer with \textit{infinite blueshift}, $z = -1$. Consequently, there exist rays, emitted soon after the BB, that will reach later observers with finite blueshift ($-1 < z < 0$). But in spacetimes without symmetry there are no radial directions. The question thus arises whether blueshifts can exist at all in the Szekeres models that contain L--T as a limit, but in general have no symmetry. The aim of the present paper is to show that strong blueshifts can be generated in quasi-spherical Szekeres (QSS) models. It is shown that in an axially symmetric QSS model, infinite blueshift can appear only on axial rays, which intersect every space orthogonal to the dust flow on the symmetry axis. In an exemplary QSS model it is numerically shown that if such a ray is emitted from the Big Bang where $\dril {t_B} r \neq 0$, then indeed observers see it with $z \approx -1$. Rays emitted shortly after the BB and running close to the symmetry axis will reach the observers with a strong blueshift, too. Then, in a toy QSS model that has no symmetry, it was shown by numerical calculations that two null lines exist such that rays in their vicinity have redshift profiles similar to those in a vicinity of the axial rays in the axially symmetric case. This indicates that rays generating infinite blueshifts exist in general QSS spacetimes and are concentrated around two directions.

It is shown that the basic observed properties of the gamma-ray bursts (GRBs) are accounted for if one assumes that the GRBs arise by blueshifting the emission radiation of hydrogen and helium generated during the last scattering epoch. The blueshift generator for a single GRB is a region with a nonconstant bang-time function $t_B(r)$ (described by a Lema\^{\i}tre -- Tolman (L--T) exact solution of Einstein's equations) matched into a homogeneous and isotropic (Friedmann) background. Blueshift visible to the present observer arises \textit{only on those rays that are emitted radially in an L--T region}. The paper presents three L--T models with different Big Bang profiles, adapted for the highest and the lowest end of the GRB frequency range. The models account for: (1) The observed frequency range of the GRBs; (2) Their limited duration; (3) The afterglows; (4) Their hypothetical collimation into narrow jets; (5) The large distances to their sources; (6) The multitude of the observed GRBs. Properties (2), (3) and (6) are accounted for only qualitatively. With a small correction of the parameters of the model, the implied perturbations of the CMB radiation will be consistent with those actually caused by the GRBs. A complete model of the Universe would consist of many L--T regions with different $t_B(r)$ profiles, matched into the same Friedmann background. This paper is meant to be an initial exploration of the possibilities offered by models of this kind; the actual fitting of all parameters to observational results requires fine-tuning of several interconnected variables and is left for a separate study.

A hypothesis is proposed that the gamma-ray bursts (GRBs) may arise by blueshifting the emission radiation of hydrogen and helium generated during the last scattering epoch. The blueshift mechanism is provided by such a Lema\^itre -- Tolman (L--T) model, in which the bang-time function $t_B(r)$ is not everywhere constant. Blueshift arises on \textit{radial} rays that are emitted over regions where $\dril{t_B} r \neq 0$. The paper presents an L--T model adapted for this purpose and shows how it accounts for the observed properties of the GRBs, some properties are accounted for only qualitatively.

In the Lema\^{\i}tre -- Tolman (L--T) models that have nonconstant bang-time function $t_B(r)$, light emitted close to those points of the Big Bang where $\dril {t_B} r \neq 0$ is blueshifted at the receiver position. The blueshifted rays are expected to perturb the temperature of the cosmic microwave background (CMB) radiation along the lines of sight of the present central observer. It is shown here that, in a general L--T model, the $t_B(r)$ can be chosen so that the blueshift-generating region is hidden before the recombination time, where the L--T model does not apply. The rest of the paper is devoted to investigating blueshifts in one specific L--T model, called L--T$(t_B)$ -- the one that duplicates the luminosity distance vs. redshift relation of the $\Lambda$CDM model using nonconstant $t_B(r)$ alone. The location of the blueshift-generating region in the L--T$(t_B)$ spacetime is determined. Profiles of redshift/blueshift along several rays intersecting the past light cone of the present central observer are calculated. The L--T$(t_B)$ model matched to Friedmann is considered, and profiles of redshift/blueshift in such a composite model are calculated. The requirement of invisibility of blueshifts makes the L--T$(t_B)$ model astrophysically unacceptable if it should apply back to the recombination time, but does not "rule out" a general L--T model -- it only puts limits on $\dril {t_B} r$.

It is shown how certain observations interpreted in the background of the Friedmann model with $\Lambda < 0 = k$ (the $\Lambda$CDM model) can be re-interpreted using the $\Lambda = 0$ Lema\^{\i}tre - Tolman (L-T) model so as to do away with the "dark energy". The purpose of the paper is to clarify the underlying geometrical relations by doing the calculations as much as possible analytically or by very simple numerical programs. In the first part of the paper (fictitious) observations of the distribution of expansion velocity along the past light cone of the observer are considered. It is shown that the whole past light cone of the $\Lambda$CDM observer can be reproduced in the L-T model with $\Lambda = 0 = E$. This is a geometric exercise that has the advantage of being free of numerical complications. In the second part, the luminosity distance - redshift relation of the $\Lambda$CDM model is duplicated using the L-T model with $-k = 2E/r^2 =$ constant $> 0$. The value of $k$ and the function $t_B(r)$ are determined by the $\Lambda$CDM parameters. General properties of this L-T model are described. Difficulties of carrying the numerical calculations through the apparent horizon are presented in detail and mostly solved. The second model is a counterexample to the claim that an L-T model mimicking $\Lambda$CDM must contain a void around the center - it has a peak of density at $R = 0$.

The notion of an apparent horizon (AH) in a collapsing object can be carried over from the Lema\^{\i}tre -- Tolman (L--T) to the quasispherical Szekeres models in three ways: 1. Literally by the definition -- the AH is the boundary of the region, in which every bundle of null geodesics has negative expansion scalar. 2. As the locus, at which null lines that are as nearly radial as possible are turned toward decreasing areal radius $R$. These lines are in general nongeodesic. The name "absolute apparent horizon" (AAH) is proposed for this locus. 3. As the boundary of a region, where null \textit{geodesics} are turned toward decreasing $R$. The name "light collapse region" (LCR) is proposed for this region (which is 3-dimensional in every space of constant $t$); its boundary coincides with the AAH. The AH and AAH coincide in the L--T models. In the quasispherical Szekeres models, the AH is different from (but not disjoint with) the AAH. Properties of the AAH and LCR are investigated, and the relations between the AAH and the AH are illustrated with diagrams using an explicit example of a Szekeres metric. It turns out that an observer who is already within the AH is, for some time, not yet within the AAH. Nevertheless, no light signal can be sent through the AH from the inside. The analogue of the AAH for massive particles is also considered.

In this article we present a model of formation of a galaxy with a black hole in the center. It is based on the Lema\^{\i}tre--Tolman solution and is a refinement of an earlier model. The most important improvement is the choice of the interior geometry of the black hole allowing for the formation of Gyrs old black holes. Other refinements are the use of an arbitrary Friedmann model as the background (unperturbed) initial state and the adaptation of the model to an arbitrary density profile of the galaxy. Our main interest was the M87 galaxy (NGC 4486), which hosts a supermassive black hole of mass $3.2\cdot 10^{9}M_{\odot}$. It is shown that for this particular galaxy, within the framework of our model and for the initial state being a perturbation of the $\Lambda$CDM model, the age of the black hole can be up to 12.7 Gyrs. The dependence of the model on the chosen parameters at the time of last scattering was also studied. The maximal age of the black hole as a function of the $\Omega_m$ and $\Omega_\Lambda$ parameters for the M87 galaxy can be 3.717 or 12.708 Gyr.

The maximally extended Reissner--Nordstr\"om (RN) manifold with $e^2 < m^2$ begs for attaching a material source to it that would preserve the infinite chain of asymptotically flat regions and evolve through the wormhole between the RN singularities. So far, the attempts were discouraging. Here we try one more possible source -- a solution found by Ruban in 1972 that is a charged generalisation of an inhomogeneous Kantowski--Sachs-type dust solution. It can be matched to the RN solution, and the matching surface must stay all the time between the two RN event horizons. However, shell crossings do not allow even half a cycle of oscillation between the maximal and the minimal size.

Advantages of inhomogeneous cosmological models that are exact solutions of Einstein's equations over linearised perturbations of homogeneous models are presented. Examples of effects that can be described in the inhomogeneous ones are given: the non-repeatable light paths, the observed anisotropies in the cosmic microwave background, the redshift drift and the maximum diameter distance. Criticisms of inhomogeneous models that are based on misunderstandings or fallacious reasonings are pointed out and corrected; these include the "weak singularity", the positivity of deceleration "theorem", the "pathology" of redshift behaviour at the "critical point" and the alleged necessity of the bang time to be constant.

Recently, inhomogeneous generalisations of the Friedmann-Lemaitre-Robertson-Walker cosmological models have gained interest in the astrophysical community and are more often employed to study cosmological phenomena. However, in many papers the inhomogeneous cosmological models are treated as an alternative to the FLRW models. In fact, they are not an alternative, but an exact perturbation of the latter, and are gradually becoming a necessity in modern cosmology. The assumption of homogeneity is just a first approximation introduced to simplify equations. So far this assumption is commonly believed to have worked well, but future and more precise observations will not be properly analysed unless inhomogeneities are taken into account. This paper reviews recent developments in the field and shows the importance of an inhomogeneous framework in the analysis of cosmological observations.

Conditions for the existence of repeatable light paths (RLPs) in the shearfree normal cosmological models are investigated. It is found that in the conformally nonflat models the only RLPs are radial null geodesics (in the spherical case) and their analogues in the plane- and hyperbolically symmetric cases. In the conformally flat Stephani models, there exist special spherically-, plane- and hyperbolically symmetric subcases, in which all null geodesics are RLPs. They are slightly more general than the Friedmann -- Lema\^{\i}tre -- Robertson -- Walker (FLRW) models of the corresponding symmetries: their curvature index function $k(t)$ and the scale factor $R(t)$ are expressed through a single function of time. In addition to that, there exist special cases of the Stephani solution in which some of the null geodesics are RLPs. All these special cases are identified.

Lema\^itre - Tolman (L--T) toy models with a central observer have been used to study the effect of large scale inhomogeneities on the SN Ia dimming. Claims that a giant void is mandatory to explain away dark energy in this framework are currently dominating. Our aim is to show that L-T models exist that reproduce a few features of the $\Lambda$CDM model, but do not contain the giant cosmic void. We propose to use two sets of data - the angular diameter distance together with the redshift-space mass-density and the angular diameter distance together with the expansion rate - both defined on the past null cone as functions of the redshift. We assume that these functions are of the same form as in the $\Lambda$CDM model. Using the Mustapha-Hellaby-Ellis algorithm, we numerically transform these initial data into the usual two L-T arbitrary functions and solve the evolution equation to calculate the mass distribution in spacetime. For both models, we find that the current density profile does not exhibit a giant void, but rather a giant hump. However, this hump is not directly observable, since it is in a spacelike relation to a present observer. The alleged existence of the giant void was a consequence of the L-T models used earlier because their generality was limited a priori by needless simplifying assumptions, like, for example, the bang-time function being constant. Instead, one can feed any mass distribution or expansion rate history on the past light cone as initial data to the L-T evolution equation. When a fully general L-T metric is used, the giant void is not implied.