Thin shell wormholes are constructed by joining two asymptotically flat spacetimes along their inner boundaries. The junction conditions imposed on the spacetimes specify the equation of state of the matter called thin shell distributed along the joined boundaries. Barcelo and Visser (2000) reported that spherically symmetric thin shell wormholes have their shells, namely the wormhole throats, on the photon spheres if the wormholes are $Z_2$-symmetric across the throats and the shells are of pure tension. In this paper, first, we consider general joined spacetimes (JSTs) and show that any $Z_2$-symmetric pure-tensional JST (Z2PTJST) of $\Lambda$-vacuum has its shell on a photon surface, a generalized object of photon spheres, without assuming any other symmetries. The class of Z2PTJSTs also includes, for example, brane world models with the shells being the branes we live in. Second, we investigate the shell stability of Z2PTJSTs by analyzing the stability of the corresponding photon surfaces. Finally, applying the uniqueness theorem of photon spheres by Cederbaum (2015), we establish the uniqueness theorem of static wormholes of Z2PTJST.

In an accretion of fluid, its velocity may transit from subsonic to supersonic. The point at which such transition occurs is called sonic point and often mathematically special. We consider a steady-state and spherically symmetric accretion problem of ideal photon gas in general static spherically symmetric spacetime neglecting back reaction. Our main result is that the EOS of ideal photon gas leads to correspondence between its sonic point and the photon sphere of the spacetime in general situations. Despite of the dependence of the EOS on the dimension of spacetime, this correspondence holds for spacetimes of arbitrary dimensions.

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).