We study shapes and alignments of 45 dark matter (DM) haloes and their brightest cluster galaxies (BCGs) using a sample of 39 massive clusters from Hubble Frontier Field (HFF), Cluster Lensing And Supernova survey with Hubble (CLASH), and Reionization Lensing Cluster Survey (RELICS). We measure shapes of the DM haloes by strong gravitational lensing, whereas BCG shapes are derived from their light profiles in Hubble Space Telescope images. Our measurements from a large sample of massive clusters presented here provide new constraints on dark matter and cluster astrophysics. We find that DM haloes are on average highly elongated with the mean ellipticity of $0.482\pm 0.028$, and position angles of major axes of DM haloes and their BCGs tend to be aligned well with the mean value of alignment angles of $22.2\pm 3.9$ deg. We find that DM haloes in our sample are on average more elongated than their BCGs with the mean difference of their ellipticities of $0.11\pm 0.03$. In contrast, the Horizon-AGN cosmological hydrodynamical simulation predicts on average similar ellipticities between DM haloes and their central galaxies. While such a difference between the observations and the simulation may well be explained by the difference of their halo mass scales, other possibilities include the bias inherent to strong lensing measurements, limited knowledge of baryon physics, or a limitation of cold dark matter.

It is known observationally that the major axes of galaxy clusters and their brightest cluster galaxies are roughly aligned with each other. To understand the origin of the alignment, we identify 40 cluster-sized dark matter (DM) haloes with masses higher than $5\times10^{13}~M_{\odot}$ and their central galaxies (CGs) at $z\approx 0$ in the Horizon-AGN cosmological hydrodynamical simulation. We trace the progenitors at 50 different epochs between $0<z<5$. We then fit their shapes and orientations with a triaxial ellipsoid model. While the orientations of both DM haloes and CGs change significantly due to repeated mergers and mass accretions, their relative orientations are well aligned at each epoch even at high redshifts, $z>1$. The alignment becomes tighter with cosmic time; the major axes of the CGs and their host DM haloes at present are aligned on average within $\sim 30^{\circ}$ in the three dimensional space and $\sim 20^{\circ}$ in the projected plane. The orientations of the major axes of DM haloes on average follow one of the eigen-vectors of the surrounding tidal field that corresponds to the {\it slowest collapsing} (or even stretching) mode, and the alignment with the tidal field also becomes tighter. This implies that the orientations of CGs and DM haloes at the present epoch are largely imprinted in the primordial density field of the Universe, whereas strong dynamical interactions such as mergers are important to explain their mutual alignment at each epoch.

We present multi-wavelength studies of optically-defined merging clusters, based on the Hyper Suprime-Cam Subaru Strategic Program. Luminous red galaxies, tracing cluster mass distributions, enable to identify cluster subhalos at various merging stages, and thus make a homogeneous sample of cluster mergers, which is unbiased with respect to the merger boost of the intracluster medium (ICM). We define, using a peak-finding method, merging clusters with multiple-peaks and single clusters with single-peaks from the CAMIRA cluster catalog. Stacked weak-lensing analysis indicates that our sample of the merging clusters is categorized into major mergers. The average halo concentration for the merging clusters is $\sim70\%$ smaller than that of the single-peak clusters, which agrees well with predictions of numerical simulations. The spatial distribution of subhalos is less centrally concentrated than the mass distribution of the main halo. The fractions of red galaxies in the merging clusters are not higher than those of the single-peak clusters. We find a signature of the merger boost of the ICM from stacked Planck Sunyaev-Zeldovich effect and ROSAT X-ray luminosity, but not in optical richness. The stacked X-ray surface brightness distribution, aligned with the main-subhalo pairs of low redshift and massive clusters, shows that the central gas core is elongated along the merger axis and overall gas distribution is misaligned by $\sim60$ deg. The homogeneous, unbiased sample of cluster mergers and multi-wavelength follow-up studies provide a unique opportunity to make a complete picture of merger physics over the whole process.

While various observations measured ellipticities of galaxy clusters and alignments between orientations of the brightest cluster galaxies and their host clusters, there are only a handful of numerical simulations that implement realistic baryon physics to allow direct comparisons with those observations. Here we investigate ellipticities of galaxy clusters and alignments between various components of them and the central galaxies in the state-of-the-art cosmological hydrodynamical simulation Horizon-AGN, which contains dark matter, stellar, and gas components in a large simulation box of $(100 h^{-1}$ Mpc$)^3$ with high spatial resolution ($\sim1$ kpc). We estimate ellipticities of total matter, dark matter, stellar, gas surface mass density distributions, X-ray surface brightness, and the Compton $y$-parameter of the Sunyaev-Zel'dovich effect, as well as alignments between these components and the central galaxies for 120 projected images of galaxy clusters with masses $M_{200}>5\times10^{13}M_{\odot}$. Our results indicate that the distributions of these components are well aligned with the major-axes of the central galaxies, with the root mean square value of differences of their position angles of $\sim 20^\circ$, which vary little from inner to the outer regions. We also estimate alignments of these various components with total matter distributions, and find tighter alignments than those for central galaxies with the root mean square value of $\sim 15^\circ$. We compare our results with previous observations of ellipticities and position angle alignments and find reasonable agreements. The comprehensive analysis presented in this paper provides useful prior information for analyzing stacked lensing signals as well as designing future observations to study ellipticities and alignments of galaxy clusters.

We perform image stacking analysis of Sloan Digital Sky Survey (SDSS) photometric galaxies over the AKARI Far-Infrared Surveyor (FIS) maps at 65{\mu}m, 90{\mu}m, and 140{\mu}m. The resulting image profiles are decomposed into the central galaxy component (single term) and the nearby galaxy component (clustering term), as a function of the r-band magnitude, m_r of the central galaxy. We find that the mean far-infrared (FIR) flux of a galaxy with magnitude m_r is well fitted with f^s_{90{\mu}m}=13*10^{0.306(18-m_r)}[mJy]. The FIR amplitude of the clustering term is consistent with that expected from the angular-correlation function of the SDSS galaxies, but galaxy morphology dependence needs to be taken into account for a more quantitative conclusion. We also fit the spectral energy distribution of stacked galaxies at 65{\mu}m, 90{\mu}m, and 140{\mu}m, and derive a mean dust temperature of ~30K. This is consistent with the typical dust temperature of galaxies that are FIR luminous and individually detected.

Recently, a novel class of models for inflation has been found in which the inflationary dynamics is driven solely by (non-canonical) kinetic terms rather than by potential terms. As an obvious extension, we show that a scalar field with non-canonical kinetic terms alone behaves like an energy component which is time-varying and has negative pressure presently, i.e. quintessence. We present a model which has a constant equation of state, that is, a ``kinetic'' counterpart of the Ratra-Peebles model of a quintessence field with a potential term. We make clear the structure of the phase plane and show that the quintessential solution is a late-time attractor. We also give a model for the ``phantom'' component which has an equation of state with $w=p/\rho <-1$.

The formation probability of multiple winding topological defects is calculated by phase and flux distribution analysis based on the Kibble mechanism. The core size of defects is taken into account so that when it is much larger than the correlation length of the Higgs field, high winding configuration can be realized. For example, if the self coupling constant of the Higgs field is smaller than $4\times10^{-4}$, the topological inflation may occur at the grand unification energy scale by the multiple winding string produced during the GUT phase transition.