The spatial clustering of active galactic nuclei (AGNs) is considered to be one of the important diagnostics for the understanding of the underlying processes behind their activities complementary to measurements of the luminosity function (LF). We analyse the AGN clustering from a recent semi-analytic model performed on a large cosmological $N$-body simulation covering a cubic gigaparsec comoving volume. We have introduced a new time-scale of gas accretion on to the supermassive black holes to account for the loss of the angular momentum on small scales, which is required to match the faint end of the observed X-ray LF. The large simulation box allows us accurate determination of the auto-correlation function of the AGNs. The model prediction indicates that this time-scale plays a significant role in allowing massive haloes to host relatively faint population of AGNs, leading to a higher bias factor for those AGNs. The model predictions are in agreement with observations of X-ray selected AGNs in the luminosity range $10^{41.5}~\mathrm{erg} \ \mathrm{s}^{-1} \leq L_{2-10\mathrm{keV}} \leq 10^{44.5}~\mathrm{erg} \ \mathrm{s}^{-1}$, with the typical host halo mass of $10^{12.5-13.5} h^{-1}\,{\rm M}_{\odot}$ at $z \lesssim 1$. This result shows that the observational clustering measurements impose an independent constraint on the accretion time-scale complementary to the LF measurements. Moreover, we find that not only the effective halo mass corresponding to the overall bias factor, but the extended shape of the predicted AGN correlation function shows remarkable agreement with those from observations. Further observational efforts towards the low luminosity end at $z \sim 1$ would give us stronger constraints on the triggering mechanisms of AGN activities through their clustering.

We show the significance of the super-Eddington accretion for the cosmic growth of supermassive black holes (SMBHs) with a semi-analytical model for galaxy and black hole evolution. The model explains various observed properties of galaxies and active galactic nuclei at a wide redshift range. By tracing the growth history of individual SMBHs, we find that the fraction of the SMBH mass acquired during the super-Eddington accretion phases to the total SMBH mass becomes larger for less massive black holes and at higher redshift. Even at z = 0, SMBHs with > 1e+9 Msun have acquired more than 50% of their mass by super-Eddington accretions, which is apparently inconsistent with classical Soltan's argument. However, the mass-weighted radiation efficiency of SMBHs with > 1e+8 Msun obtained with our model, is about 0.08 at z = 0, which is consistent with Soltan's argument within the observational uncertainties. We, therefore, conclude that Soltan's argument cannot reject the possibility that SMBHs are grown mainly by super-Eddington accretions.

We construct a model of H$\alpha$ emitters (HAEs) based on a semi-analytic galaxy formation model, the New Numerical Galaxy Catalog ($\nu^2$GC). In this paper, we report our estimate for the field variance of the HAE distribution. By calculating the H$\alpha$ luminosity from the star-formation rate of galaxies, our model well reproduces the observed H$\alpha$ luminosity function (LF) at $z=0.4$. The large volume of the $\nu^2$GC makes it possible to examine the spatial distribution of HAEs over a region of (411.8 Mpc)$^3$ in the comoving scale. The surface number density of $z=0.4$ HAEs with $L_{\rm H\alpha} \geq 10^{40}$ erg s$^{-1}$ is 308.9 deg$^{-2}$. We have confirmed that the HAE is a useful tracer for the large-scale structure of the Universe because of their significant overdensity ($>$ 5$\sigma$) at clusters and the filamentary structures. The H$\alpha$ LFs within a survey area of $\sim$2 deg$^2$ (typical for previous observational studies) show a significant field variance up to $\sim$1 dex. Based on our model, one can estimate the variance on the H$\alpha$ LFs within given survey areas.

We present the clustering properties of low-$z$ $(z\leq1.4)$ galaxies selected by the Hyper Suprime-Cam Subaru Strategic Program Wide layer over $145$ deg$^{2}$. The wide-field and multi-wavelength observation yields $5,064,770$ galaxies at $0.3\leq z\leq1.4$ with photometric redshifts and physical properties. This enables the accurate measurement of angular correlation functions and subsequent halo occupation distribution (HOD) analysis allows the connection between baryonic properties and dark halo properties. The fraction of less-massive satellite galaxies at $z\lesssim1$ is found to be almost constant at $\sim20\%$, but it gradually decreases beyond $M_{\star} \sim 10^{10.4}h^{-2}M_{\odot}$. However, the abundance of satellite galaxies at $z>1$ is quite small even for less-massive galaxies due to the rarity of massive centrals at high-$z$. This decreasing trend is connected to the small satellite fraction of Lyman break galaxies at $z>3$. The stellar-to-halo mass ratios at $0.3\leq z\leq1.4$ are almost consistent with the predictions obtained using the latest empirical model; however, we identify small excesses from the theoretical model at the massive end. The pivot halo mass is found to be unchanged at $10^{11.9-12.1}h^{-1}M_{\odot}$ at $0.3\leq z\leq1.4$, and we systematically show that $10^{12}h^{-1}M_{\odot}$ is a universal pivot halo mass up to $z\sim5$ that is derived using only the clustering/HOD analyses. Nevertheless, halo masses with peaked instantaneous baryon conversion efficiencies are much smaller than the pivot halo mass regardless of a redshift, and the most efficient stellar-mass assembly is thought to be in progress in $10^{11.0-11.5}h^{-1}M_{\odot}$ dark haloes.

Galaxy imaging surveys enable us to map the cosmic matter density field through weak gravitational lensing analysis. The density reconstruction is compromised by a variety of noise originating from observational conditions, galaxy number density fluctuations, and intrinsic galaxy properties. We propose a deep-learning approach based on generative adversarial networks (GANs) to reduce the noise in the weak lensing map under realistic conditions. We perform image-to-image translation using conditional GANs in order to produce noiseless lensing maps using the first-year data of the Subaru Hyper Suprime-Cam (HSC) survey. We train the conditional GANs by using 30000 sets of mock HSC catalogs that directly incorporate observational effects. We show that an ensemble learning method with GANs can reproduce the one-point probability distribution function (PDF) of the lensing convergence map within a $0.5-1\sigma$ level. We use the reconstructed PDFs to estimate a cosmological parameter $S_{8} = \sigma_{8}\sqrt{\Omega_{\rm m0}/0.3}$, where $\Omega_{\rm m0}$ and $\sigma_{8}$ represent the mean and the scatter in the cosmic matter density. The reconstructed PDFs place tighter constraint, with the statistical uncertainty in $S_8$ reduced by a factor of $2$ compared to the noisy PDF. This is equivalent to increasing the survey area by $4$ without denoising by GANs. Finally, we apply our denoising method to the first-year HSC data, to place $2\sigma$-level cosmological constraints of $S_{8} < 0.777 \, ({\rm stat}) + 0.105 \, ({\rm sys})$ and $S_{8} < 0.633 \, ({\rm stat}) + 0.114 \, ({\rm sys})$ for the noisy and denoised data, respectively.

We perform an ensemble of $N$-body simulations with $2048^3$ particles for 101 flat $w$CDM cosmological models sampled based on a maximin-distance Sliced Latin Hypercube Design. By using the halo catalogs extracted at multiple redshifts in the range of $z=[0,1.48]$, we develop Dark Emulator, which enables fast and accurate computations of the halo mass function, halo-matter cross-correlation, and halo auto-correlation as a function of halo masses, redshift, separations and cosmological models, based on the Principal Component Analysis and the Gaussian Process Regression for the large-dimensional input and output data vector. We assess the performance of the emulator using a validation set of $N$-body simulations that are not used in training the emulator. We show that, for typical halos hosting CMASS galaxies in the Sloan Digital Sky Survey, the emulator predicts the halo-matter cross correlation, relevant for galaxy-galaxy weak lensing, with an accuracy better than $2\%$ and the halo auto-correlation, relevant for galaxy clustering correlation, with an accuracy better than $4\%$. We give several demonstrations of the emulator. It can be used to study properties of halo mass density profiles such as the mass-concentration relation and splashback radius for different cosmologies. The emulator outputs can be combined with an analytical prescription of halo-galaxy connection such as the halo occupation distribution at the equation level, instead of using the mock catalogs, to make accurate predictions of galaxy clustering statistics such as the galaxy-galaxy weak lensing and the projected correlation function for any model within the $w$CDM cosmologies, in a few CPU seconds.

The distributions of the pairwise line-of-sight velocity between galaxies and their host clusters are segregated according to the galaxy's colour and morphology. We investigate the velocity distribution of red-spiral galaxies, which represents a rare population within galaxy clusters. We find that the probability distribution function of the pairwise line-of-sight velocity $v_{\rm{los}}$ between red-spiral galaxies and galaxy clusters has a dip at $v_{\rm{los}} = 0$, which is a very odd feature, at 93\% confidence level. To understand its origin, we construct a model of the phase space distribution of galaxies surrounding galaxy clusters in three-dimensional space by using cosmological $N$-body simulations. We adopt a two component model that consists of the infall component, which corresponds to galaxies that are now falling into galaxy clusters, and the splashback component, which corresponds to galaxies that are on their first (or more) orbit after falling into galaxy clusters. We find that we can reproduce the distribution of the line-of-sight velocity of red-spiral galaxies with the dip with a very simple assumption that red-spiral galaxies reside predominantly in the infall component, regardless of the choice of the functional form of their spatial distribution. Our results constrain the quenching timescale of red-spiral galaxies to a few Gyrs, and the radius where the morphological transformation is effective as $r \sim 0.2 h^{-1} \rm{Mpc}$.

We show the Eddington ratio distributions of supermassive black holes at a wide redshift range (0 < z < 8) obtained with a semi-analytic model of galaxy formation. The distribution is broadly consistent with observational estimates at low redshift. We find that the growth rate of black holes at higher redshift is more likely to exceed the Eddington limit because the typical gas fraction of the host galaxies is higher at higher redshift. We also find that the super- Eddington growth is more common for less massive supermassive black holes, supporting an idea that supermassive black holes have been formed via super-Eddington accretion. These results indicate the "slowing down" of cosmic growth of supermassive black holes: the growth of supermassive black holes with a higher Eddington ratio peaks at higher redshift. We also show the effect of the sample selection on the shape of the Eddington ratio distribution functions and find that shallower observations will miss active galactic nuclei with not only the smaller but also higher Eddington ratios.

We present the latest results of a semi-analytic model of galaxy formation, "New Numerical Galaxy Catalogue", which is combined with large cosmological N-body simulations. This model can reproduce statistical properties of galaxies at z < 6.0. We focus on the properties of active galactic nuclei (AGNs) and supermassive black holes, especially on the accretion timescale onto black holes. We find that the number density of AGNs at z < 1.5 and at hard X-ray luminosity 10^{ 44 }< erg/s is underestimated compared with recent observational estimates when we assume the exponentially decreasing accretion rate and the accretion timescale which is proportional to the dynamical time of the host halo or the bulge, as is often assumed in semi-analytic models. We show that to solve this discrepancy, the accretion timescale of such less luminous AGNs instead should be a function of the black hole mass and the accreted gas mass. This timescale can be obtained from a phenomenological modelling of the gas angular momentum loss in the circumnuclear torus and/or the accretion disc. Such models predict a longer accretion timescale for less luminous AGNs at z < 1.0 than bright QSOs whose accretion timescale would be 10^{ 7-8 } yr. With this newly introduced accretion timescale, our model can explain the observed luminosity functions of AGNs at z < 6.0.

We present morphological properties of dusty star-forming galaxies at z=1-3 determined with high-resolution (FWHM~0"19) Atacama Large Milllimeter/submilimeter Array (ALMA) 1-mm band maps of our ASAGAO survey covering a 26-arcmin^2 area in GOODS-S. In conjunction with the ALMA archival data, the present sample consists of 42 ALMA sources with a wide rest-frame far-infrared (FIR) luminosity L_FIR range of ~10^11-10^13 Lo. To obtain an average rest-frame FIR profile, we perform individual measurements and careful stacking of the ALMA sources using the uv-visibility method that includes positional-uncertainty and smoothing-effect evaluations through Monte-Carlo simulations. We find that the dusty star-forming galaxies have the average FIR-wavelength Sersic index and effective radius of n_FIR=1.2+/-0.2 and R_e,FIR=1.0-1.3 kpc, respectively, additionally with a point source at the center, indicative of the existence of AGN. The average FIR profile agrees with a morphology of an exponential-disk clearly distinguished from a spheroidal profile (Sersic index of 4). We also examine the rest-frame optical Sersic index n_opt and effective radius R_e,opt with the deep Hubble Space Telescope (HST) images. Interestingly, we obtain n_opt=0.9+/-0.3 (~n_FIR) and R_e,opt=3.2+/-0.6 kpc (>R_e,FIR), suggesting that the FIR-emitting disk is embedded within a larger stellar disk. The rest-frame UV and FIR data of HST and ALMA provide us a radial surface density profile of the total star-formation rate (SFR), where the FIR SFR dominates over the UV SFR at the center. Under the simple assumption of a constant SFR, a compact stellar distribution found in z~1-2 compact quiescent galaxies (cQGs) is well reproduced, while a spheroidal stellar morphology of cQGs (n_opt=4) cannot, suggestive of other important mechanisms such as dynamical dissipation.

We present a statistical study of the redshift evolution of the cluster galaxy population over a wide redshift range from 0.1 to 1.1, using $\sim 1900$ optically-selected CAMIRA clusters from $\sim 232$~deg$^2$ of the Hyper Suprime-Cam (HSC) Wide S16A data. Our stacking technique with a statistical background subtraction reveals color-magnitude diagrams of red-sequence and blue cluster galaxies down to faint magnitudes of $m_z\sim 24$. We find that the linear relation of red-sequence galaxies in the color-magnitude diagram extends down to the faintest magnitudes we explore with a small intrinsic scatter $\sigma_{\rm int}(g-r)<0.1$. The scatter does not evolve significantly with redshift. The stacked color-magnitude diagrams are used to define red and blue galaxies in clusters for studying their radial number density profiles without resorting to photometric redshifts of individual galaxies. We find that red galaxies are significantly more concentrated toward cluster centers and blue galaxies dominate the outskirt of clusters. We explore the fraction of red galaxies in clusters as a function of redshift, and find that the red fraction decreases with increasing distances from cluster centers. The red fraction exhibits a moderate decrease with increasing redshift. The radial number density profiles of cluster member galaxies are also used to infer the location of the steepest slope in the three dimensional galaxy density profiles. For a fixed threshold in richness, we find little redshift evolution in this location.

Super-Eddington mass accretion has been suggested as an efficient mechanism to grow supermassive black holes (SMBHs). We investigate the imprint left by the radiative efficiency of the super-Eddington accretion process on the clustering of quasars using a new semi-analytic model of galaxy and quasar formation based on large-volume cosmological $N$-body simulations. Our model includes a simple model for the radiative efficiency of a quasar, which imitates the effect of photon trapping for a high mass accretion rate. We find that the model of radiative efficiency affects the relation between the quasar luminosity and the quasar host halo mass. The quasar host halo mass has only weak dependence on quasar luminosity when there is no upper limit for quasar luminosity. On the other hand, it has significant dependence on quasar luminosity when the quasar luminosity is limited by its Eddington luminosity. In the latter case, the quasar bias also depends on the quasar luminosity, and the quasar bias of bright quasars is in agreement with observations. Our results suggest that the quasar clustering studies can provide a constraint on the accretion disc model.

We examine the clustering of quasars over a wide luminosity range, by utilizing 901 quasars at $\overline{z}_{\rm phot}\sim3.8$ with $-24.73<M_{\rm 1450}<-22.23$ photometrically selected from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) S16A Wide2 date release and 342 more luminous quasars at $3.4<z_{\rm spec}<4.6$ having $-28.0<M_{\rm 1450}<-23.95$ from the Sloan Digital Sky Survey (SDSS) that fall in the HSC survey fields. We measure the bias factors of two quasar samples by evaluating the cross-correlation functions (CCFs) between the quasar samples and 25790 bright $z\sim4$ Lyman Break Galaxies (LBGs) in $M_{\rm 1450}<-21.25$ photometrically selected from the HSC dataset. Over an angular scale of \timeform{10.0"} to \timeform{1000.0"}, the bias factors are $5.93^{+1.34}_{-1.43}$ and $2.73^{+2.44}_{-2.55}$ for the low and high luminosity quasars, respectively, indicating no luminosity dependence of quasar clustering at $z\sim4$. It is noted that the bias factor of the luminous quasars estimated by the CCF is smaller than that estimated by the auto-correlation function (ACF) over a similar redshift range, especially on scales below \timeform{40.0"}. Moreover, the bias factor of the less-luminous quasars implies the minimal mass of their host dark matter halos (DMHs) is $0.3$-$2\times10^{12}h^{-1}M_{\odot}$, corresponding to a quasar duty cycle of $0.001$-$0.06$.

We present measurements of the clustering properties of a sample of infrared (IR) bright dust-obscured galaxies (DOGs). Combining 125 deg$^2$ of wide and deep optical images obtained with the Hyper Suprime-Cam on the Subaru Telescope and all-sky mid-IR (MIR) images taken with Wide-Field Infrared Survey Explorer, we have discovered 4,367 IR-bright DOGs with $(i - [22])_{\rm AB}$ $>$ 7.0 and flux density at 22 $\mu$m $>$ 1.0 mJy. We calculate the angular autocorrelation function (ACF) for a uniform subsample of 1411 DOGs with 3.0 mJy $<$ flux (22 $mu$m) $<$ 5.0 mJy and $i_{\rm AB}$ $<$ 24.0. The ACF of our DOG subsample is well-fit with a single power-law, $\omega (\theta)$ = (0.010 $\pm$ 0.003) $\theta^{-0.9}$, where $\theta$ in degrees. The correlation amplitude of IR-bright DOGs is larger than that of IR-faint DOGs, which reflects a flux-dependence of the DOG clustering, as suggested by Brodwin et al. (2008). We assume that the redshift distribution for our DOG sample is Gaussian, and consider 2 cases: (1) the redshift distribution is the same as IR-faint DOGs with flux at 22 $\mu$m $<$ 1.0 mJy, mean and sigma $z$ = 1.99 $\pm$ 0.45, and (2) $z$ = 1.19 $\pm$ 0.30, as inferred from their photometric redshifts. The inferred correlation length of IR-bright DOGs is $r_0$ = 12.0 $\pm$ 2.0 and 10.3 $\pm$ 1.7 $h^{-1}$ Mpc, respectively. IR-bright DOGs reside in massive dark matter halos with a mass of $\log [\langle M_{\mathrm{h}} \rangle / (h^{-1} M_{\odot})]$ = 13.57$_{-0.55}^{+0.50}$ and 13.65$_{-0.52}^{+0.45}$ in the two cases, respectively.

We present a new cosmological galaxy formation model, $\nu^2$GC, as an updated version of our previous model $\nu$GC. We adopt the so-called "semi-analytic" approach, in which the formation history of dark matter halos is computed by ${\it N}$-body simulations, while the baryon physics such as gas cooling, star formation and supernova feedback are simply modeled by phenomenological equations. Major updates of the model are as follows: (1) the merger trees of dark matter halos are constructed in state-of-the-art ${\it N}$-body simulations, (2) we introduce the formation and evolution process of supermassive black holes and the suppression of gas cooling due to active galactic nucleus (AGN) activity, (3) we include heating of the intergalactic gas by the cosmic UV background, and (4) we tune some free parameters related to the astrophysical processes using a Markov chain Monte Carlo method. Our ${\it N}$-body simulations of dark matter halos have unprecedented box size and mass resolution (the largest simulation contains 550 billion particles in a 1.12 Gpc/h box), enabling the study of much smaller and rarer objects. The model was tuned to fit the luminosity functions of local galaxies and mass function of neutral hydrogen. Local observations, such as the Tully-Fisher relation, size-magnitude relation of spiral galaxies and scaling relation between the bulge mass and black hole mass were well reproduced by the model. Moreover, the model also well reproduced the cosmic star formation history and the redshift evolution of rest-frame ${\it K}$-band luminosity functions. The numerical catalog of the simulated galaxies and AGNs is publicly available on the web.

We explore the effect of varying the mass of the seed black hole on the resulting black hole mass - bulge mass relation at z ~ 0, using a semi-analytic model of galaxy formation combined with large cosmological N-body simulations. We constrain our model by requiring the observed properties of galaxies at z ~ 0 are reproduced. In keeping with previous semi-analytic models, we place a seed black hole immediately after a galaxy forms. When the mass of the seed is set at 10^5 M_sun, we find that the model results become inconsistent with recent observational results of the black hole mass - bulge mass relation for dwarf galaxies. In particular, the model predicts that bulges with ~ 10^9 M_sun harbour larger black holes than observed. On the other hand, when we employ seed black holes with 10^3 M_sun, or randomly select their mass within a 10^(3-5) M_sun range, the resulting relation is consistent with observation estimates, including the observed dispersion. We find that to obtain stronger constraints on the mass of seed black holes, observations of less massive bulges at z ~ 0 are a more powerful comparison than the relations at higher redshifts.

We present the evolution of dark matter halos in six large cosmological N-body simulations, called the $\nu^2$GC (New Numerical Galaxy Catalog) simulations on the basis of the LCDM cosmology consistent with observational results obtained by the Planck satellite. The largest simulation consists of $8192^3$ (550 billion) dark matter particles in a box of $1.12 \, h^{-1} \rm Gpc$ (a mass resolution of $2.20 \times 10^{8} \, h^{-1} M_{\odot}$). Among simulations utilizing boxes larger than $1 \, h^{-1} \rm Gpc$, our simulation yields the highest resolution simulation that has ever been achieved. A $\nu^2$GC simulation with the smallest box consists of eight billions particles in a box of $70 \, h^{-1} \rm Mpc$ (a mass resolution of $3.44 \times 10^{6} \, h^{-1} M_{\odot}$). These simulations can follow the evolution of halos over masses of eight orders of magnitude, from small dwarf galaxies to massive clusters. Using the unprecedentedly high resolution and powerful statistics of the $\nu^2$GC simulations, we provide statistical results of the halo mass function, mass accretion rate, formation redshift, and merger statistics, and present accurate fitting functions for the Planck cosmology. By combining the $\nu^2$GC simulations with our new semi-analytic galaxy formation model, we are able to prepare mock catalogs of galaxies and active galactic nuclei, which will be made publicly available in the near future.

We perform dry merger simulations to investigate the role of dry mergers in the size growth of early-type galaxies in high density environments. We replace the virialized dark matter haloes obtained by a large cosmological $N$-body simulation with $N$-body galaxy models consisting of two components, a stellar bulge and a dark matter halo, which have higher mass resolution than the cosmological simulation. We then re-simulate nine cluster forming regions, whose masses range from 1e+14 Msun to 5e+14 Msun. Masses and sizes of stellar bulges are also assumed to satisfy the stellar mass--size relation of high-z compact massive early-type galaxies. We find that dry major mergers considerably contribute to the mass and size growth of central massive galaxies. One or two dry major mergers double the average stellar mass and quadruple the average size between $z=2$ and $z=0$. These growths favorably agree with observations. Moreover, the density distributions of our simulated central massive galaxies grow from the inside-out, which is consistent with recent observations. The mass--size evolution is approximated as R propto M_{*}^{alpha}, with alpha ~ 2.24. Most of our simulated galaxies are efficiently grown by dry mergers, and their stellar mass--size relations match the ones observed in the local Universe. Our results show that the central galaxies in the cluster haloes are potential descendants of high-z (z ~ 2-3) compact massive early-type galaxies. This conclusion is consistent with previous numerical studies which investigate the formation and evolution of compact massive early-type galaxies.

We investigate clustering properties of quasars using a new version of our semi-analytic model of galaxy and quasar formation with state-of-the-art cosmological N-body simulations. In this study, we assume that a major merger of galaxies triggers cold gas accretion on to a supermassive black hole and quasar activity. Our model can reproduce the downsizing trend of the evolution of quasars. We find that the median mass of quasar host dark matter haloes increases with cosmic time by an order of magnitude from z=4 (a few 1e+11 Msun) to z=1 (a few 1e+12 Msun), and depends only weakly on the quasar luminosity. Deriving the quasar bias through the quasar--galaxy cross-correlation function in the model, we find that the quasar bias does not depend on the quasar luminosity, similar to observed trends. This result reflects the fact that quasars with a fixed luminosity have various Eddington ratios and thus have various host halo masses that primarily determine the quasar bias. We also show that the quasar bias increases with redshift, which is in qualitative agreement with observations. Our bias value is lower than the observed values at high redshifts, implying that we need some mechanisms that make quasars inactive in low-mass haloes and/or that make them more active in high-mass haloes.

We report the discovery of a new ultra-bright submillimeter galaxy (SMG) behind the Small Magellanic Cloud (SMC). This SMG is detected as a 43.3+-8.4 mJy point source (MM J01071-7302, hereafter MMJ0107) in the 1.1 mm continuum survey of the SMC by AzTEC on the ASTE telescope. MMJ0107 is also detected in the radio (843 MHz), Herschel/SPIRE, Spitzer MIPS 24 {\mu}m, all IRAC bands, Wide-field Infrared Survey Explorer, and near-infrared (J, H, KS). We find an optical (U, B, V) source, which might be the lensing object, at a distance of 1.4 arcsec from near-infrared and IRAC sources. Photometric redshift estimates for the SMG using representative spectral energy distribution templates show the redshifts of 1.4-3.9. We estimate total far-infrared luminosity of (0.3-2.2)x10^14 {\mu}^-1 L_sun and a star formation rate of 5600-39, 000 {\mu}^-1 M_sun yr^-1, where {\mu} is the gravitational magnification factor. This apparent extreme star formation activity is likely explained by a highly magnified gravitational lens system.

Recent observations show evidence that high-z (z\sim 2 - 3) early-type galaxies (ETGs) are more compact than those with comparable mass at z\sim 0. Such a size evolution is most likely explained by the `Dry Merger Sceanario'. However, previous studies based on this scenario are not able to consistantly explain both the properties of the high-z compact massive ETGs and the local ETGs. We investigate the effect of multiple sequential dry minor mergers on the size evolution of the compact massive ETGs. From an analysis of the Millennium Simulation Database, we show that such minor (stellar mass ratio $M_{2}/M_{1} < 1/4$) mergers are extremely common during hierarchical structure formation. We perform N-body simulations of sequential minor mergers with parabolic and head-on orbits, including a dark matter component and a stellar component. Typical mass ratios of the minor mergers are $1/20 < M_{2}/M_{1} < 1/10$. We show that sequential minor mergers of compact satellite galaxies are the most efficient at promoting size growth and decreasing the velocity dispersion of the compact massive ETGs in our simulations. The change of stellar size and density of the merger remnants is consistent with recent observations. Furthermore, we construct the merger histories of candidates for the high-z compact massive ETGs using the Millennium Simulation Database, and estimate the size growth of the galaxies by the dry minor merger scenario. We can reproduce the mean size growth factor between $z=2$ and $z=0$, assuming the most efficient size growth obtained during sequential minor mergers in our simulations. However, we note that our numerical result is only valid for merger histories with typical mass ratios between 1/20 and 1/10 with parabolic and head-on orbits, and that our most efficient size growth efficiency is likely to an upper limit.