We perform (3+1)-dimensional simulations of primordial black hole (PBH) formation starting from the spheroidal super-horizon perturbations. We investigate how the ellipticity (prolateness or oblateness) affects the threshold of PBH formation in terms of the peak amplitude of curvature perturbation. We find that, in the case of the radiation-dominated universe, the effect of ellipticity on the threshold is negligibly small for large amplitude of perturbations expected for PBH formation.

We discuss the effect of local type non-Gaussianity on the abundance of primordial black holes (PBH) based on the peak theory. We provide the PBH formation criterion based on the so-called compaction function and use the peak theory statistics associated with the curvature perturbation with the local type non-Gaussianity. Providing a method to estimate the PBH abundance, we demonstrate the effects of non-Gaussianity. It is explicitly shown that the value of non-linear parameter $|f_{\rm NL}| \sim 1$ induces a similar effect to a few factors of difference in the amplitude of the power spectrum.

The production rate of primordial black holes is often calculated by considering a nearly Gaussian distribution of cosmological perturbations, and assuming that black holes will form in regions where the amplitude of such perturbations exceeds a certain threshold. A threshold $\zeta_{\rm th}$ for the curvature perturbation is somewhat inappropriate for this purpose, because it depends significantly on environmental effects, not essential to the local dynamics. By contrast, a threshold $\delta_{\rm th}$ for the density perturbation at horizon crossing seems to provide a more robust criterion. On the other hand, the density perturbation is known to be bounded above by a maximum limit $\delta_{\rm max}$, and given that $\delta_{\rm th}$ is comparable to $\delta_{\rm max}$, the density perturbation will be far from Gaussian near or above the threshold. In this paper, we provide a new plausible estimate for the primordial black hole abundance based on peak theory. In our approach, we assume that the curvature perturbation is given as a random Gaussian field with the power spectrum characterized by a single scale, while an optimized criterion for PBH formation is imposed, based on the locally averaged density perturbation. Both variables are related by the full nonlinear expression derived in the long-wavelength approximation of general relativity. We do not introduce a window function, and the scale of the inhomogeneity is introduced as a random variable in the peak theory. We find that the mass spectrum is shifted to larger mass scales by one order of magnitude or so, compared to a conventional calculation. The abundance of PBHs becomes significantly larger than the conventional one, by many orders of magnitude, mainly due to the optimized criterion for PBH formation and the removal of the suppresion associated with a window function.

We show that there is a class of spacetime curvature singularities which can be resolved with metric and matter field transformations. As an example, we consider an anisotropic power-law inflation model with gauge and scalar fields in which a space-like curvature singularity exists at the beginning of time. First, we provide a transformation of the metric to the flat geometry. The transformation is regular in the whole region of spacetime except for the singularity. In general, matter fields are still singular after such a metric transformation. However, we explicitly show that there is a case in which the singular behavior of the matter fields can be completely removed by a field re-definition. Since the action is invariant under any metric and matter field transformations, the regularity of the action at the original singularity is a necessary condition for the complete removal of a singularity.

Angular momentum plays very important roles in the formation of PBHs in the matter-dominated phase if it lasts sufficiently long. In fact, most collapsing masses are bounced back due to centrifugal force, since angular momentum significantly grows before collapse. As a consequence, most of the formed PBHs are rapidly rotating near the extreme value $a_{*}=1$, where $a_{*}$ is the nondimensional Kerr parameter at their formation. The smaller the density fluctuation $\sigma_{H}$ at horizon entry is, the stronger the tendency towards the extreme rotation. Combining the effect of angular momentum with that of anisotropy, we estimate the black hole production rate. We find that the production rate suffers from suppression dominantly due to angular momentum for a smaller value of $\sigma_{H}$, while due to anisotrpopy for a larger value of $\sigma_{H}$. We argue that matter domination significantly enhances the production of PBHs despite the suppression. If the matter-dominated phase does not last so long, the effect of the finite duration significantly suppresses PBH formation and weakens the tendency towards large spins. (abridged)

We investigate the validity of screening mechanism of the fifth force for chameleon field in highly inhomogeneous density profile. For simplicity, we consider a static and spherically symmetric system which is composed of concentric infinitely thin shells. We calculate the fifth force profile by using a numerical method for a relatively large Compton wavelength of the chameleon field. An approximate solution is also derived for the small Compton wavelength limit. Our results show that, if the thin-shell condition for the corresponding smoothed density profile is satisfied, the fifth force is safely screened outside the system irrespective of the configuration of the shells inside the system. In contrast to the screening outside the system, we find that the fifth force can be comparable to the Newtonian gravitational force inside the system. If the system is highly inhomogeneous, the chameleon field cannot trace the potential minimum varying with the density and repeats being kicked, climbing up and rolling down the potential even when the effective mass of the chameleon field is sufficiently large in the system on average. One should not feel complacent about the wellbehavedness of the fifth force field with an averaged density distribution when we consider inhomogeneous objects.

We derive an observational constraint on a spherical inhomogeneity of the void centered at our position from the angular power spectrum of the cosmic microwave background(CMB) and local measurements of the Hubble parameter. The late time behaviour of the void is assumed to be well described by the so-called $\Lambda$-Lema\^itre-Tolman-Bondi~($\Lambda$LTB) solution. Then, we restrict the models to the asymptotically homogeneous models each of which is approximated by a flat Friedmann-Lema\^itre-Robertson-Walker model. The late time $\Lambda$LTB models are parametrized by four parameters including the value of the cosmological constant and the local Hubble parameter. The other two parameters are used to parametrize the observed distance-redshift relation. Then, the $\Lambda$LTB models are constructed so that they are compatible with the given distance-redshift relation. Including conventional parameters for the CMB analysis, we characterize our models by seven parameters in total. The local Hubble measurements are reflected in the prior distribution of the local Hubble parameter. As a result of a Markov-Chains-Monte-Carlo analysis for the CMB temperature and polarization anisotropies, we found that the inhomogeneous universe models with vanishing cosmological constant are ruled out as is expected. However, a significant under-density around us is still compatible with the angular power spectrum of CMB and the local Hubble parameter.

We investigate primordial black hole formation in the matter-dominated phase of the Universe, where nonspherical effects in gravitational collapse play a crucial role. This is in contrast to the black hole formation in a radiation-dominated era. We apply the Zel'dovich approximation, Thorne's hoop conjecture, and Doroshkevich's probability distribution and subsequently derive the production probability $\beta_{0}$ of primordial black holes. The numerical result obtained is applicable even if the density fluctuation $\sigma$ at horizon entry is of the order of unity. For $\sigma\ll 1$, we find a semi-analytic formula $\beta_{0}\simeq 0.05556 \sigma^{5}$, which is comparable with the Khlopov-Polnarev formula. We find that the production probability in the matter-dominated era is much larger than that in the radiation-dominated era for $\sigma\lesssim 0.05$, while they are comparable with each other for $\sigma\gtrsim 0.05$. We also discuss how $\sigma$ can be written in terms of primordial curvature perturbations.

Spherically symmetric dust universe models with a positive cosmological constant $\Lambda$, known as $\Lambda$-Lema\^itre-Tolman-Bondi($\Lambda$LTB) models, are considered. We report a method to construct the $\Lambda$LTB model from a given distance-redshift relation observed at the symmetry center. The spherical inhomogeneity is assumed to be composed of growing modes. We derive a set of ordinary differential equations for three functions of the redshift, which specify the spherical inhomogeneity. Once a distance-redshift relation is given, with careful treatment of possible singular points, we can uniquely determine the model by solving the differential equations for each value of $\Lambda$. As a demonstration, we fix the distance-redshift relation as that of the flat $\Lambda$CDM model with $(\Omega^{\rm dis}_{\rm m0}, \Omega^{\rm dis}_{\rm \Lambda 0})=(0.3,0.7)$, where $\Omega^{\rm dis}_{\rm m0}$ and $\Omega^{\rm dis}_{\rm \Lambda 0}$ are the normalized matter density and the cosmological constant, respectively. Then, we construct the $\Lambda$LTB model for several values of $\Omega_{\rm \Lambda 0}:=\Lambda/(3H_0^2)$, where $H_0$ is the present Hubble parameter observed at the symmetry center. We obtain void structure around the symmetry center for $\Omega_{\Lambda 0}<\Omega^{\rm dis}_{\Lambda 0}$. We show the relation between the ratio $\Omega_{\Lambda0}/\Omega^{\rm dis}_{\Lambda 0}$ and the amplitude of the inhomogeneity.

The Square Kilometre Array will revolutionize pulsar studies with its wide field-of-view, wide-band observation and high sensitivity, increasing the number of observable pulsars by more than an order of magnitude. Pulsars are of interest not only for the study of neutron stars themselves but for their usage as tools for probing fundamental physics such as general relativity, gravitational waves and nuclear interaction. In this article, we summarize the activity and interests of SKA-Japan Pulsar Science Working Group, focusing on an investigation of modified gravity theory with the supermassive black hole in the Galactic Centre, gravitational-wave detection from cosmic strings and binary supermassive black holes, a study of the physical state of plasma close to pulsars using giant radio pulses and determination of magnetic field structure of Galaxy with pulsar pairs.

The so-called black hole shadow is a dark region which is expected to appear in a fine image of optical observation of black holes. It is essentially an absorption cross section of black hole, and the boundary of shadow is determined by unstable circular orbits of photons (UCOP). If there exists a compact object possessing UCOP but no black hole horizon, it can provide us with the same shadow image with black holes, and a detection of shadow image cannot be a direct evidence of black hole existence. Then, this paper examine whether or not such compact objects can exist under some suitable conditions. We investigate thoroughly the static spherical polytropic ball of perfect fluid with single polytrope index, and then investigate a representative example of the piecewise polytropic ball. Our result is that the spherical polytropic ball which we have investigated cannot possess UCOP, if the sound speed at center is subluminal (slower-than-light). This means that, if the polytrope treated in this paper is a good model of stellar matter in compact objects, the detection of shadow image is regarded as a good evidence of black hole existence. As a by-product, we have found the upper bound of the mass-to-radius radio (M/R) of polytropic ball with single index, M/R < 0.281, under the subluminal-sound-speed condition.

The discrepancy between the amplitudes of matter fluctuations inferred from Sunyaev-Zel'dovich (SZ) cluster number counts, the primary temperature, and the polarization anisotropies of the cosmic microwave background (CMB) measured by the Planck satellite can be reconciled if the local universe is embedded in an under-dense region as shown by Lee, 2014. Here using a simple void model assuming the open Friedmann-Robertson-Walker geometry and a Markov Chain Monte Carlo technique, we investigate how deep the local under-dense region needs to be to resolve this discrepancy. Such local void, if exists, predicts the local Hubble parameter value that is different from the global Hubble constant. We derive the posterior distribution of the local Hubble parameter from a joint fitting of the Planck CMB data and SZ cluster number counts assuming the simple void model. We show that the predicted local Hubble parameter value of $H_{\rm loc}=70.1\pm0.34~{\rm km\,s^{-1}Mpc^{-1}}$ is in better agreement with direct local Hubble parameter measurements, indicating that the local void model may provide a consistent solution to the cluster number counts and Hubble parameter discrepancies.

The ratio of total mass $M$ to surface radius $R$ of spherical perfect fluid ball has an upper bound, $M/R < B$. Buchdahl obtained $B = 4/9$ under the assumptions; non-increasing mass density in outward direction, and barotropic equation of states. Barraco and Hamity decreased the Buchdahl's bound to a lower value $B = 3/8$ $(< 4/9)$ by adding the dominant energy condition to Buchdahl's assumptions. In this paper, we further decrease the Barraco-Hamity's bound to $B \simeq 0.3636403$ $(< 3/8)$ by adding the subluminal (slower-than-light) condition of sound speed. In our analysis, we solve numerically Tolman-Oppenheimer-Volkoff equations, and the mass-to-radius ratio is maximized by variation of mass, radius and pressure inside the fluid ball as functions of mass density.

We construct cosmological long-wavelength solutions without symmetry in general gauge conditions compatible with the long-wavelength scheme. We then specify the relationship among the solutions in different time slicings. Applying this general framework to spherical symmetry, we derive the correspondence relation between long-wavelength solutions in the constant mean curvature slicing with conformally flat spatial coordinates and asymptotic quasihomogeneous solutions in the comoving gauge and compare the numerical results of PBH formation in these two different approaches. To discuss the PBH formation, it is convenient and conventional to use $\tilde{\delta}_{c}$, the value which the averaged density perturbation at threshold in the comoving slicing would take at horizon entry in the lowest-order long-wavelength expansion. We numerically find that within compensated models, the sharper the transition from the overdense region to the FRW universe is, the larger the $\tilde{\delta}_{c}$ becomes. We suggest that, for the equation of state $p=(\Gamma-1)\rho$, we can apply the analytic formulas for the minimum $\tilde{\delta}_{c, {\rm min}}\simeq [3\Gamma/(3\Gamma+2)]\sin^{2}\left[\pi\sqrt{\Gamma-1}/(3\Gamma-2)\right]$ and the maximum $\tilde{\delta}_{c, {\rm max}}\simeq 3\Gamma/(3\Gamma+2)$. As for the threshold peak value of the curvature variable $\psi_{0,c}$, we find that the sharper the transition is, the smaller the $\psi_{0,c}$ becomes. We analytically explain this feature. Using simplified models, we also analytically deduce an environmental effect that $\psi_{0,c}$ can be significantly larger (smaller) if the underlying density perturbation of much longer wavelength is positive (negative).

Usually the effects of isotropic inhomogeneities are not seriously taken into account in the determination of the cosmological parameters because of Copernican principle whose statement is that we do not live in the privileged domain in the universe. But Copernican principle has not been observationally confirmed yet in sufficient accuracy, and there is the possibility that there are non-negligible large-scale isotropic inhomogeneities in our universe. In this paper, we study the effects of the isotropic inhomogeneities on the determination of the cosmological parameters and show the probability that non-Copernican isotropic inhomogeneities mislead us into believing, for example, the phantom energy of the equation of state, $p=w\rho$ with $w<-1$, even in case that $w=-1$ is the true value.

Recently, the present authors studied perturbations in the Lemaitre-Tolman-Bondi cosmological model by applying the second-order perturbation theory in the dust Friedmann-Lemaitre-Robertson-Walker universe model. Before this work, the same subject was studied in some papers by analyzing linear perturbations in the Lemaitre-Tolman-Bondi cosmological model under the assumption proposed by Clarkson, Clifton and February, in which two of perturbation variables are negligible. However, it is a non-trivial issue in what situation the Clarkson-Clifton-February assumption is valid. In this paper, we investigate differences between these two approaches. It is shown that, in general, these two approaches are not compatible with each other. That is, in our perturbative procedure, the Clarkson-Clifton-February assumption is not valid at the order of our interest.

We consider a test of the Copernican Principle through observations of the large-scale structures, and for this purpose we study the self-gravitating system in a relativistic huge void universe model which does not invoke the Copernican Principle. If we focus on the the weakly self-gravitating and slowly evolving system whose spatial extent is much smaller than the scale of the cosmological horizon in the homogeneous and isotropic background universe model, the cosmological Newtonian approximation is available. Also in the huge void universe model, the same kind of approximation as the cosmological Newtonian approximation is available for the analysis of the perturbations contained in a region whose spatial size is much smaller than the scale of the huge void: the effects of the huge void are taken into account in a perturbative manner by using the Fermi-normal coordinates. By using this approximation, we derive the equations of motion for the weakly self-gravitating perturbations whose elements have relative velocities much smaller than the speed of light, and show the derived equations can be significantly different from those in the homogeneous and isotropic universe model, due to the anisotropic volume expansion in the huge void. We linearize the derived equations of motion and solve them. The solutions show that the behaviors of linear density perturbations are very different from those in the homogeneous and isotropic universe model.

We investigate the effect of small scale inhomogeneities on standard candle observations, such as type Ia supernovae (SNe) observations. Existence of the small scale inhomogeneities may cause a tension between SNe observations and other observations with larger diameter sources, such as the cosmic microwave background (CMB) observation. To clarify the impact of the small scale inhomogeneities, we use the Dyer-Roeder approach. We determined the smoothness parameter $\alpha(z)$ as a function of the redshift $z$ so as to compensate the deviation of cosmological parameters for SNe from those for CMB. The range of the deviation which can be compensated by the smoothness parameter $\alpha(z)$ satisfying $0\leq\alpha(z)\leq1$ is reported. Our result suggests that the tension may give us the information of the small scale inhomogeneities through the smoothness parameter.

Time evolution of a black hole lattice universe with a positive cosmological constant $\Lambda$ is simulated. The vacuum Einstein equations are numerically solved in a cubic box with a black hole in the center. Periodic boundary conditions on all pairs of opposite faces are imposed. Configurations of marginally trapped surfaces are analyzed. We describe the time evolution of not only black hole horizons, but also cosmological horizons. Defining the effective scale factor by using the area of a surface of the cubic box, we compare it with that in the spatially flat dust dominated FLRW universe with the same value of $\Lambda$. It is found that the behaviour of the effective scale factor is well approximated by that in the FLRW universe. Our result suggests that local inhomogeneities do not significantly affect the global expansion law of the universe irrespective of the value of $\Lambda$.

Based on a physical argument, we derive a new analytic formula for the amplitude of density perturbation at the threshold of primordial black hole formation in the universe dominated by a perfect fluid with the equation of state $p=w\rho c^{2}$ for $w\ge 0$. The formula gives $\delta^{\rm UH}_{H c}=\sin^{2}[\pi \sqrt{w}/(1+3w)]$ and $\tilde{\delta}_{c}=[3(1+w)/(5+3w)]\sin^{2}[\pi\sqrt{w}/(1+3w)]$, where $\delta^{\rm UH}_{H c}$ and $\tilde{\delta}_{c}$ are the amplitude of the density perturbation at the horizon crossing time in the uniform Hubble slice and the amplitude measure used in numerical simulations, respectively, while the conventional one gives $\delta^{\rm UH}_{H c}=w$ and $\tilde{\delta}_{c}=3w(1+w)/(5+3w)$. Our formula shows a much better agreement with the result of recent numerical simulations both qualitatively and quantitatively than the conventional formula. For a radiation fluid, our formula gives $\delta^{\rm UH}_{H c}=\sin^{2}(\sqrt{3}\pi/6)\simeq 0.6203$ and $\tilde{\delta}_{c}=(2/3)\sin^{2}(\sqrt{3}\pi/6)\simeq 0.4135$. We also discuss the maximum amplitude and the cosmological implications of the present result.

We study the network of Type-I cosmic strings using the field-theoretic numerical simulations in the Abelian-Higgs model. For Type-I strings, the gauge field plays an important role, and thus we find that the correlation length of the strings is strongly dependent upon the parameter \beta, the ratio between the self-coupling constant of the scalar field and the gauge coupling constant, namely, \beta=\lambda/2e^2. In particular, if we take the cosmic expansion into account, the network becomes densest in the comoving box for a specific value of \beta for \beta<1.

Time evolution of a black hole lattice universe is simulated. The vacuum Einstein equations in a cubic box with a black hole at the origin are numerically solved with periodic boundary conditions on all pairs of opposite faces. Defining effective scale factors by using the area of a surface and the length of an edge of the cubic box, we compare them with that in the Einstein-deSitter universe. It is found that the behaviour of the effective scale factors is well approximated by that in the Einstein-deSitter universe. Our result suggests that local inhomogeneities do not significantly affect the global expansion law of the universe even if the inhomogeneity is extremely nonlinear.

We study the two-point correlation function of density perturbations in a spherically symmetric void universe model which does not employ the Copernican principle. First we solve perturbation equations in the inhomogeneous universe model and obtain density fluctuations by using a method of non-linear perturbation theory which was adopted in our previous paper. From the obtained solutions, we calculate the two-point correlation function and show that it has a local anisotropy at the off-center position differently from those in homogeneous and isotropic universes. This anisotropy is caused by the tidal force in the off-center region of the spherical void. Since no tidal force exists in homogeneous and isotropic universes, we may test the inhomogeneous universe by observing statistical distortion of the two-point galaxy correlation function.

We consider the femto-lensing due to a cosmic string. If a cosmic string with the deficit angle $\Delta\sim 100$ [femto-arcsec] $\sim10^{-18}$ [rad] exists around the line of sight to a gamma-ray burst, we may observe characteristic interference patterns caused by gravitational lensing in the energy spectrum of the gamma-ray burst. This "femto-lensing" event was first proposed as a tool to probe small mass primordial black holes. In this paper, we propose use of the femto-lensing to probe cosmic strings with extremely small tension. Observability conditions and the event rate are discussed. Differences between the cases of a point mass and a cosmic string are presented.

We numerically construct an one-parameter family of initial data of an expanding inhomogeneous universe model which is composed of regularly aligned black holes with an identical mass. They are initial data for vacuum solutions of the Einstein equations. We call this universe model the "black hole universe" and analyze the structure of these initial data. We study the relation between the mean expansion rate of the 3-space, which corresponds to the Hubble parameter, and the mass density of black holes. The result implies that the same relation as that of the Einstein-de Sitter universe is realized in the limit of the large separation between neighboring black holes. The applicability of the cosmological Newtonian $N$-body simulation to the dark matter composed of black holes is also discussed. The deviation of the spatial metric of the cosmological Newtonian $N$-body system from that of the black hole universe is found to be smaller than about 1% in a region distant from the particles, if the separation length between neighboring particles is 20 times larger than their gravitational radius. By contrast, the deviation of the square of the Hubble parameter of the cosmological Newtonian $N$-body system from that of the black hole universe is about 20% for the same separation length.

We study the evolution of linear density perturbations in a large spherical void universe which accounts for the acceleration of the cosmic volume expansion without introducing dark energy. The density contrast of this void is not large within the light cone of an observer at the center of the void. Therefore, we describe the void structure as a perturbation with a dimensionless small parameter $\kappa$ in a homogeneous and isotropic universe within the region observable for the observer. We introduce additional anisotropic perturbations with a dimensionless small parameter $\epsilon$, whose evolution is of interest. Then, we solve perturbation equations up to order $\kappa \epsilon$ by applying second-order perturbation theory in the homogeneous and isotropic universe model. By this method, we can know the evolution of anisotropic perturbations affected by the void structure. We show that the growth rate of the anisotropic density perturbations in the large void universe is significantly different from that in the homogeneous and isotropic universe. This result suggests that the observation of the distribution of galaxies may give a strong constraint on the large void universe model.

We study the effect of weak lensing by cosmic (super-)strings on the anisotropies of cosmic microwave background (CMB). In developing a method to calculate the lensing convergence field due to strings, and thereby temperature and polarization angular power spectra of CMB, we clarify how the nature of strings, characterized by the intercommuting probability, can influence these CMB anisotropies. Assuming that the power spectrum is dominated by Poisson-distributed string segments, we find that the convergence spectrum peaks at low multipoles and is mostly contributed from strings located at relatively low redshifts. As the intercommuting probability decreases, the spectra of the convergence and hence the lensed temperature and polarizations are gained because the number density of string segments becomes larger. An observationally important feature of the string-induced CMB polarizations is that the string-lensed spectra decay more slowly on small scales compared with primordial scalar perturbations from standard inflation.

If we assume that we live in the center of a spherical inhomogeneous universe, we can explain the apparent accelerating expansion of the universe without introducing the unknown dark energy or modifying gravitational theory. Direct measurement of the cosmic acceleration can be a powerful tool in distinguishing $\Lambda$CDM and the inhomogeneous models. If $\Lambda$CDM is the correct model, we have shown that DECIGO/BBO has sufficient ability to detect the positive redshift drift of the source by observing gravitational waves from neutron star binaries for 5-10 years. This enables us to rule out any Lema\^itre-Tolman-Bondi (LTB) void model with monotonically increasing density profile. Furthermore, by detecting the positive redshift drift at $z\sim 0$, we can even rule out generic LTB models unless we allow unrealistically steep density gradient at $z\sim 0$. We also show that the measurement accuracy is slightly improved when we consider the joint search of DECIGO/BBO and the third generation Einstein Telescope. This test can be performed with GW observations alone without any reference to electromagnetic observations.

One possibility for explaining the apparent accelerating expansion of the universe is that we live in the center of a spherically inhomogeneous universe. Although current observations cannot fully distinguish $\Lambda$CDM and these inhomogeneous models, direct measurement of the acceleration of the universe can be a powerful tool in probing them. We have shown that, if $\Lambda$CDM is the correct model, DECIGO/BBO would be able to detect the positive redshift drift (which is the time evolution of the source redshift $z$) in 3--5 year gravitational wave (GW) observations from neutron-star binaries, which enables us to rule out any Lema\^itre-Tolman-Bondi (LTB) void model with monotonically increasing density profile. We may even be able to rule out any LTB model unless we allow unrealistically steep density profile at $z\sim 0$. This test can be performed with GW observations alone, without any reference to electromagnetic observations, and is more powerful than the redshift drift measurement using Lyman $\alpha$ forest.

The inverse problem with Lema\^itre-Tolman-Bondi (LTB) universe models is discussed. The LTB solution for the Einstein equations describes the spherically symmetric dust-filled spacetime. The LTB solution has two physical functional degrees of freedom of the radial coordinate. The inverse problem is constructing an LTB model requiring that the LTB model be consistent with selected important observational data. In this paper, we assume that the observer is at the center and consider the distance-redshift relation $\da$ and the redshift-space mass density $\mu$ as the selected important observational data. We give $\da$ and $\mu$ as functions of the redshift $z$. Then, we explicitly show that, for general functional forms of $\da(z)$ and $\mu(z)$, the regular solution does not necessarily exist in the whole redshift domain. We also show that the condition for the existence of the regular solution %in terms of $\da(z)$ and $\mu(z)$ is satisfied by the distance-redshift relation and the redshift-space mass density in $\Lambda$CDM models. Deriving regular differential equations for the inverse problem with the distance-redshift relation and the redshift-space mass density in $\Lambda$CDM models, we numerically solve them for the case $(\Omega_{\rm M0},\Omega_{\Lambda0})=(0.3,0.7)$. A set of analytic fitting functions for the resultant LTB universe model is given. How to solve the inverse problem with the simultaneous big-bang and a given function $\da(z)$ for the distance-redshift relation is provided in the Appendix.

It is known that the Meissner-like effect is seen in a magnetosphere without an electric current in black hole spacetime: no non-monopole component of magnetic flux penetrates the event horizon if the black hole is extreme. In this paper, in order to see how an electric current affects the Meissner-like effect, we study a force-free electromagnetic system in a static and spherically symmetric extreme black hole spacetime. By assuming that the rotational angular velocity of the magnetic field is very small, we construct a perturbative solution for the Grad-Shafranov equation, which is the basic equation to determine a stationary, axisymmetric electromagnetic field with a force-free electric current. Our perturbation analysis reveals that, if an electric current exists, higher multipole components may be superposed upon the monopole component on the event horizon, even if the black hole is extreme.

We study the redshift drift, i.e., the time derivative of the cosmological redshift in the Lema\^itre-Tolman-Bondi (LTB) solution in which the observer is assumed to be located at the symmetry center. This solution has often been studied as an anti-Copernican universe model to explain the acceleration of cosmic volume expansion without introducing the concept of dark energy. One of decisive differences between LTB universe models and Copernican universe models with dark energy is believed to be the redshift drift. The redshift drift is negative in all known LTB universe models, whereas it is positive in the redshift domain $z \lesssim 2$ in Copernican models with dark energy. However, there have been no detailed studies on this subject. In the present paper, we prove that the redshift drift of an off-center source is always negative in the case of LTB void models. We also show that the redshift drift can be positive with an extremely large hump-type inhomogeneity. Our results suggest that we can determine whether we live near the center of a large void without dark energy by observing the redshift drift.

We study the kinematic Sunyaev-Zel'dovich (kSZ) effect in a Lem\^itre-Tolman-Bondi (LTB) universe model whose distance-redshift relation agrees with that of the concordance $\Lambda$CDM model at redshifts $z\lesssim2$. This LTB universe model has a void with size comparable to the Hubble horizon scale. We first determine the decoupling epoch in this LTB universe model by an approximate analytical condition under a few simplified assumptions on the physical quantities at that epoch. Then we calculate the cosmic microwave background (CMB) anisotropy observed in the rest frame of clusters of galaxies which are assumed to be at rest in the spatial comoving coordinates of the LTB universe model. We find that the obtained temperature anisotropies are dominated by dipole, although there may exist higher multi-poles in general. We may interpret this dipole anisotropy as the drift velocity of a cluster of galaxies relative to the CMB rest frame. Hence it gives rise to the kSZ effect. We calculate this effect and compare it with observational data. We find that if we assume the conventional adiabatic perturbation scenario at the time of decoupling, the drift velocity of clusters of galaxies becomes unacceptably large. Conversely, this observational constraint may be relaxed by introducing a non-adiabatic (i.e., primordially isocurvature) component of inhomogeneities at the time of decoupling. However, our result indicates that the necessary isocurvature perturbation amplitude is very large.

We present a new analytical method to calculate the small angle CMB temperature angular power spectrum due to cosmic (super-)string segments. In particular, using our method, we clarify the dependence on the intercommuting probability $P$. We find that the power spectrum is dominated by Poisson-distributed string segments. The power spectrum for a general value of $P$ has a plateau on large angular scales and shows a power-law decrease on small angular scales. The resulting spectrum in the case of conventional cosmic strings is in very good agreement with the numerical result obtained by Fraisse et al.. Then we estimate the upper bound on the dimensionless tension of the string for various values of $P$ by assuming that the fraction of the CMB power spectrum due to cosmic (super-)strings is less than ten percents at various angular scales up to $\ell=2000$. We find that the amplitude of the spectrum increases as the intercommuting probability. As a consequence, strings with smaller intercommuting probabilities are found to be more tightly constrained.

Acoustic peaks in the spectrum of the cosmic microwave background in spherically symmetric inhomogeneous cosmological models are studied. At the photon-baryon decoupling epoch, the universe may be assumed to be dominated by non-relativistic matter, and thus we may treat radiation as a test field in the universe filled with dust which is described by the Lema\^itre-Tolman-Bondi (LTB) solution. First, we give an LTB model whose distance-redshift relation agrees with that of the concordance $\Lambda$CDM model in the whole redshift domain and which is well approximated by the Einstein-de Sitter universe at and before decoupling. We determine the decoupling epoch in this LTB universe by Gamow's criterion and then calculate the positions of acoustic peaks. Thus obtained results are not consistent with the WMAP data. However, we find that one can fit the peak positions by appropriately modifying the LTB model, namely, by allowing the deviation of the distance-redshift relation from that of the concordance $\Lambda$CDM model at $z>2$ where no observational data are available at present. Thus there is still a possibility of explaining the apparent accelerated expansion of the universe by inhomogeneity without resorting to dark energy if we abandon the Copernican principle. Even if we do not take this extreme attitude, it also suggests that local, isotropic inhomogeneities around us may seriously affect the determination of the density contents of the universe unless the possible existence of such inhomogeneities is properly taken into account.

We compute the one-point probability distribution function of small-angle cosmic microwave background temperature fluctuations due to curved cosmic (super-)strings with a simple model of string network by performing Monte Carlo simulations. Taking into account of the correlation between the curvature and the velocity of string segments, there appear non-Gaussian features, specifically non-Gaussian tails and a skewness, in the one-point pdf. The obtained sample skewness for the conventional field-theoretic cosmic strings is $g_1\approx -0.14$, which is consistent with the result reported by Fraisse et al.. We also discuss the dependence of the pdf on the intercommuting probability. We find that the standard deviation of the Gaussian part increases and non-Gaussian features are suppressed as the intercommuting probability decreases. For sufficiently small intercommuting probability, the skewness is given by $\lesssim$ (a\ few) $\times 10^{-2}$.

We compute analytically the small-scale temperature fluctuations of the cosmic microwave background from cosmic (super-)strings and study the dependence on the string intercommuting probability $P$. We develop an analytical model which describes the evolution of a string network and calculate the numbers of string segments and kinks in a horizon volume. Then we derive the probability distribution function (pdf) which takes account of finite angular resolution of observation. The resultant pdf consists of a Gaussian part due to frequent scatterings by long string segments and a non-Gaussian tail due to close encounters with kinks. The dispersion of the Gaussian part is reasonably consistent with that obtained by numerical simulations by Fraisse et al.. On the other hand, the non-Gaussian tail contains two phenomenological parameters which are determined by comparison with the numerical results for P=1. Extrapolating the pdf to the cases with $P<1$, we predict that the non-Gaussian feature is suppressed for small $P$.

We construct the Lema\^itre-Tolman-Bondi (LTB) dust universe whose distance-redshift relation is equivalent to that in the concordance $\Lambda$ cold dark matter ($\Lambda$CDM) cosmological model. In our model, the density distribution and velocity field are not homogeneous, whereas the big-bang time is uniform, which implies that the universe is homogeneous at its beginning. We also study the effects of local clumpiness in the density distribution as well as the effects of large-scale inhomogeneities on the distance-redshift relation, and show that these effects may reduce the amplitude of large-scale inhomogeneities necessary for having a distance-redshift relation that is the same as that of the concordance $\Lambda$CDM universe. We also study the temporal variation of the cosmological redshift and show that, by the observation of this quantity, we can distinguish our LTB universe model from the concordance $\Lambda$CDM model, even if their redshift-distance relations are equivalent to each other.

We study cosmology of the Einstein-Yang-Mills theory in ten dimensions with a quartic term in the Yang-Mills field strength. We obtain analytically a class of cosmological solutions in which the extra dimensions are static and the scale factor of the four-dimensional Friedmann-Lemaitre-Robertson-Walker metric is an exponential function of time. This means that the model can explain inflation. Then we look for solutions that describe dynamical compactification of the extra dimensions. The effective cosmological constant $\lambda_1$ in the four-dimensional universe is determined from the gravitational coupling, ten-dimensional cosmological constant, gauge coupling and higher derivative coupling. By numerical integration, the solution with $\lambda_1=0$ is found to behave as a matter-dominated universe which asymptotically approaches flat space-time, while the solution with a non-vanishing $\lambda_1$ approaches de Sitter space-time in the asymptotic future.

Lensing effects on light rays from point light sources, such like Type Ia supernovae, are simulated in a clumpy universe model. In our universe model, it is assumed that all matter in the universe takes the form of randomly distributed objects each of which has finite size and is transparent for light rays. Monte-Carlo simulations are performed for several lens models, and we compute probability distribution functions of magnification. In the case of the lens models that have a smooth density profile or the same degree of density concentration as the spherical NFW (Navarro-Frenk-White) lens model at the center, the so-called gamma distributions fit well the magnification probability distribution functions if the size of lenses is sufficiently larger than the Einstein radius. In contrast, the gamma distributions do not fit the magnification probability distribution functions in the case of the SIS (Singular Isothermal Sphere) lens model. We find, by using the power law cusp model, that the magnification probability distribution function is fitted well using the gamma distribution only when the slope of the central density profile is not very steep. These results suggest that we may obtain information about the slope of the central density profiles of dark matter halo from the lensing effect of Type Ia supernovae.

The distance-redshift relation determined by means of gravitational waves in the clumpy universe is simulated numerically by taking into account the effects of gravitational lensing. It is assumed that all of the matter in the universe takes the form of randomly distributed point masses, each of which has the identical mass $M_L$. Calculations are carried out in two extreme cases: $\lambda\gg GM_L/c^2$ and $\lambda\ll GM_L/c^2$, where $\lambda$ denotes the wavelength of gravitational waves. In the first case, the distance-redshift relation for the fully homogeneous and isotropic universe is reproduced with a small distance dispersion, whereas in the second case, the distance dispersion is larger. This result suggests that we might obtain information about the typical mass of lens objects through the distance-redshift relation gleaned through observation of gravitational waves of various wavelengths. In this paper, we show how to set limitations on the mass $M_L$ through the observation of gravitational waves in the clumpy universe model described above.