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.

High-speed impact ejecta at velocities comparable to the impact velocity are expected to contribute to material transport between planetary bodies and deposition of ejecta far from the impact crater. We investigated the behavior of high-speed ejecta produced at angles of 45 and 90 degrees, using both experimental and numerical methods. The experimental system developed at the Planetary Exploration Research Center of Chiba Institute of Technology (Japan) allowed us to observe the initial growth of the ejecta. We succeeded in imaging high-speed ejecta at 0.2 $\mathrm{{\mu}}$s intervals for impacts of polycarbonate projectiles of 4.8 mm diameter onto a polycarbonate plate at an impact velocity of ~4 km s$^{-1}$. Smoothed particle hydrodynamics (SPH) simulations of various numerical resolutions were conducted for the same impact conditions as pertaining to the experiments. We compared the morphology and velocities of the ejecta for the experiments and simulations, and we confirmed a close match for high-resolution simulations (with $\geq10^6$ SPH particles representing the projectile). According to the ejecta velocity distributions obtained from our high-resolution simulations, the ejection velocities of the high-speed ejecta for oblique impacts are much greater than those for vertical impacts. The translational motion of penetrating projectiles parallel to the target surface in oblique impacts could cause long-term, sustained acceleration at the root of the ejecta.

We study effect of magnetic field on massive dense core formation in colliding unequal molecular clouds by performing magnetohydrodynamic simulations with sub-parsec resolution (0.015 pc) that can resolve the molecular cores. Initial clouds with the typical gas density of the molecular clouds are immersed in various uniform magnetic fields. The turbulent magnetic fields in the clouds consistent with the observation by Crutcher et al. (2010) are generated by the internal turbulent gas motion before the collision, if the uniform magnetic field strength is 4.0 $\mu$G. The collision speed of 10 km s$^{-1}$ is adopted, which is much larger than the sound speeds and the Alfv\'{e}n speeds of the clouds. We identify gas clumps with gas densities greater than 5$\times$10$^{-20}$ g cm$^{-3}$ as the dense cores and trace them throughout the simulations to investigate their mass evolution and gravitational boundness. We show that a greater number of massive, gravitationally bound cores are formed in the strong magnetic field (4.0 $\mu$G) models than the weak magnetic field (0.1 $\mu$G) models. This is partly because the strong magnetic field suppresses the spatial shifts of the shocked layer that should be caused by the nonlinear thin shell instability. The spatial shifts promote formation of low-mass dense cores in the weak magnetic field models. The strong magnetic fields also support low-mass dense cores against gravitational collapse. We show that the numbers of massive, gravitationally bound cores formed in the strong magnetic field models are much larger than the isolated, non-colliding cloud models, which are simulated for comparison. We discuss the implications of our numerical results on massive star formation.

We report on the detection of the [CII] 157.7 $\mu$m emission from the Lyman break galaxy (LBG) MACS0416_Y1 at z = 8.3113, by using the Atacama Large Millimeter/submillimeter Array (ALMA). The luminosity ratio of [OIII] 88 $\mu$m (from previous campaigns) to [CII] is 9.31 $\pm$ 2.6, indicative of hard interstellar radiation fields and/or a low covering fraction of photo-dissociation regions. The emission of [CII] is cospatial to the 850 $\mu$m dust emission (90 $\mu$m rest-frame, from previous campaigns), however the peak [CII] emission does not agree with the peak [OIII] emission, suggesting that the lines originate from different conditions in the interstellar medium. We fail to detect continuum emission at 1.5 mm (160 $\mu$m rest-frame) down to 18 $\mu$Jy (3$\sigma$). This nondetection places a strong limit on the dust spectrum, considering the 137 $\pm$ 26 $\mu$Jy continuum emission at 850 $\mu$m. This suggests an unusually warm dust component (T $>$ 80 K, 90% confidence limit), and/or a steep dust-emissivity index ($\beta_{\rm dust}$ $>$ 2), compared to galaxy-wide dust emission found at lower redshifts (typically T $\sim$ 30 - 50 K, $\beta_{\rm dust}$ $\sim$ 1 - 2). If such temperatures are common, this would reduce the required dust mass and relax the dust production problem at the highest redshifts. We therefore warn against the use of only single-wavelength information to derive physical properties, recommend a more thorough examination of dust temperatures in the early Universe, and stress the need for instrumentation that probes the peak of warm dust in the Epoch of Reionization.

We succeeded in observing two large spicules simultaneously with the Atacama Large Millimeter/submillimeter Array (ALMA), the Interface Region Imaging Spectrograph (IRIS), and the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory. One is a spicule seen in the IRIS Mg II slit-jaw images and AIA 304\AA\ images (MgII/304A spicule). The other one is a spicule seen in the 100GHz images obtained with ALMA (100GHz spicule). Although the 100GHz spicule overlapped with the MgII/304A spicule in the early phase, it did not show any corresponding structures in the IRIS Mg II and AIA 304A images after the early phase. It suggests that the spicules are individual events and do not have a physical relationship. To obtain the physical parameters of the 100GHz spicule, we estimate the optical depths as a function of temperature and density using two different methods. One is using the observed brightness temperature by assuming a filling factor, and the other is using an emission model for the optical depth. As a result of comparing them, the kinetic temperature of the plasma and the number density of ionized hydrogens in the 100GHz spicule are ~6800 K and 2.2 x 10^10 cm^-3. The estimated values can explain the absorbing structure in the 193A image, which appear as a counterpart of the 100GHz spicule. These results suggest that the 100GHz spicule presented in this paper is classified to a macrospicule without a hot sheath in former terminology.

We investigated the off-limb spicules observed in the Mg II h and k lines by IRIS in a solar polar coronal hole. We analyzed the large dataset of obtained spectra to extract quantitative information about the line intensities, shifts, and widths. The observed Mg II line profiles are broad and double-peaked at lower altitudes, broad but flat-topped at middle altitudes, and narrow and single-peaked with the largest Doppler shifts at higher altitudes. We use 1D non-LTE vertical slab models (i.e. models which consider departures from Local Thermodynamic Equilibrium) in single-slab and multi-slab configurations to interpret the observations and to investigate how a superposition of spicules along the line of sight (LOS) affects the synthetic Mg II line profiles. The used multi-slab models are either static, i.e. without any LOS velocities, or assume randomly assigned LOS velocities of individual slabs, representing the spicule dynamics. We conducted such single-slab and multi-slab modeling for a broad set of model input parameters and showed the dependence of the Mg II line profiles on these parameters. We demonstrated that the observed line widths of the h and k line profiles are strongly affected by the presence of multiple spicules along the LOS. We later showed that the profiles obtained at higher altitudes can be reproduced by single-slab models representing individual spicules. We found that the multi-slab model with a random distribution of the LOS velocities ranging from -25 to 25 km s$^{-1}$ can well reproduce the width and the shape of Mg II profiles observed at middle altitudes.

Dry lakebeds might constitute large volatile reservoirs on Mars. Hypervelocity impacts onto ancient dry lakebeds would have affected the volatile distribution on Mars. We developed a new experimental method to investigate the response of evaporitic minerals (halite and gypsum) to impact shocks in an open system. This technique does not result in chemical contamination from the operation of the gas gun. The technique is termed the two-valve method and the gun system is located in the Planetary Exploration Research Center, Chiba Institute of Technology, Japan. We detected the vaporization of halite at 31 GPa and devolatilization from gypsum at 11 GPa, suggesting that impact-induced volatile release from dry lakebeds has periodically occurred throughout Martian history. The vaporization of halite deposits might have enhanced the production of perchlorates, which are found globally on Mars. The water loss from gypsum possibly explains the coexisting types of Ca-sulfates found in Gale Crater.

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 new ALMA observations and physical properties of a Lyman Break Galaxy at z=7.15. Our target, B14-65666, has a bright ultra-violet (UV) absolute magnitude, $M_{\rm UV}\approx-22.4$, and has been spectroscopically identified in Ly$\alpha$ with a small rest-frame equivalent width of $\approx4$ \AA. Previous HST image has shown that the target is comprised of two spatially separated clumps in the rest-frame UV. With ALMA, we have newly detected spatially resolved [OIII] 88 $\mu$m, [CII] 158 $\mu$m, and their underlying dust continuum emission. In the whole system of B14-65666, the [OIII] and [CII] lines have consistent redshifts of $7.1520\pm0.0003$, and the [OIII] luminosity, $(34.4\pm4.1)\times10^{8}L_{\rm \odot}$, is about three times higher than the [CII] luminosity, $(11.0\pm1.4)\times10^{8}L_{\rm \odot}$. With our two continuum flux densities, the dust temperature is constrained to be $T_{\rm d}\approx50-60$ K under the assumption of the dust emissivity index of $\beta_{\rm d}=2.0-1.5$, leading to a large total infrared luminosity of $L_{\rm TIR}\approx1\times10^{12}L_{\rm \odot}$. Owing to our high spatial resolution data, we show that the [OIII] and [CII] emission can be spatially decomposed into two clumps associated with the two rest-frame UV clumps whose spectra are kinematically separated by $\approx200$ km s$^{-1}$. We also find these two clumps have comparable UV, infrared, [OIII], and [CII] luminosities. Based on these results, we argue that B14-65666 is a starburst galaxy induced by a major-merger. The merger interpretation is also supported by the large specific star-formation rate (defined as the star-formation rate per unit stellar mass), sSFR$=260^{+119}_{-57}$ Gyr$^{-1}$, inferred from our SED fitting. Probably, a strong UV radiation field caused by intense star formation contributes to its high dust temperature and the [OIII]-to-[CII] luminosity ratio.

We present the Atacama Large Millimeter/submillimeter Array (ALMA) detection of the [O III] 88 $\mu$m line and rest-frame 90 $\mu$m dust continuum emission in a Y-dropout Lyman break galaxy (LBG), MACS0416_Y1, lying behind the Frontier Field cluster MACS J0416.1-2403. This [O III] detection confirms the LBG with a spectroscopic redshift of $z = 8.3118 \pm 0.0003$, making this object one of the furthest galaxies ever identified spectroscopically. The observed 850 $\mu$m flux density of $137 \pm 26$ $\mu$Jy corresponds to a de-lensed total infrared (IR) luminosity of $L_{\rm IR} = (1.7 \pm 0.3) \times 10^{11} L_{\odot}$ if assuming a dust temperature of $T_{\rm dust} = 50$ K and an emissivity index of $\beta = 1.5$, yielding a large dust mass of $4 \times 10^6 M_{\odot}$. The ultraviolet-to-far IR spectral energy distribution modeling where the [O III] emissivity model is incorporated suggests the presence of a young ($\tau_{\rm age} \approx 4$ Myr), star-forming (SFR $\approx 60 M_{\odot}$ yr$^{-1}$), moderately metal-polluted ($Z \approx 0.2 Z_{\odot}$) stellar component with a mass of $M_{\rm star} = 3 \times 10^8 M_{\odot}$. An analytic dust mass evolution model with a single episode of star-formation does not reproduce the metallicity and dust mass in $\tau_{\rm age} \approx 4$ Myr, suggesting a pre-existing evolved stellar component with $M_{\rm star} \sim 3 \times 10^9 M_{\odot}$ and $\tau_{\rm age} \sim 0.3$ Gyr as the origin of the dust mass.

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 results from a pilot observation of nearby (~20 Mpc) galaxies with mass similar to that of the Milky Way (MW) to address the missing satellite problem. This is the first paper from an on-going project to address the problem with a statistical sample of galaxies outside of the Local Group (LG) without employing an assumption that the LG is a typical halo in the Universe. Thanks to the close distances of our targets, dwarf galaxies around them can be identified as extended, diffuse galaxies. By applying a surface brightness cut together with a careful visual screening to remove artifacts and background contamination, we construct a sample of dwarf galaxies. The luminosity function (LF) of one of the targets is broadly consistent with that of the MW, but the other has a more abundant dwarf population. Numerical simulations by Okamoto (2013) seem to overpredict the number of dwarfs on average, while more recent predictions from Copernicus Complexio are in a better agreement. In both observations and simulations, there is a large diversity in the LFs, demonstrating the importance of addressing the missing satellite problem with a statistically representative sample. We also characterize the projected spatial distributions of the satellites and do not observe strong evidence for alignments around the central galaxies. Based on this successful pilot observation, we are carrying out further observations to increase the sample of nearby galaxies, which we plan to report in our future paper.

We report the results of the Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of the solar chromosphere on the southern polar limb. Coordinated observations with the Interface Region Imaging Spectrograph (IRIS) are also conducted. ALMA provided unprecedented high spatial resolution in the millimeter band ($\approx$ 2.0 arcsec) at 100 GHz frequency with a moderate cadence (20 s). The results are as follows: (1) The ALMA 100 GHz images show saw-tooth patterns on the limb, and a comparison with SDO/AIA 171\AA\ images shows a good correspondence of the limbs with each other. (2) The ALMA 100 GHz movie shows a dynamic thorn-like structure elongating from the saw-tooth patterns on the limb, with lengths reaching at least 8 arcsec, thus suggesting jet-like activity in the ALMA microwave range. These ALMA jets are in good correspondence with IRIS jet clusters. (3) A blob ejection event is observed. By comparing with the IRIS Mg II slit-jaw images, the trajectory of the blob is located along the spicular patterns.

A fundamental quest of modern astronomy is to locate the earliest galaxies and study how they influenced the intergalactic medium a few hundred million years after the Big Bang. The abundance of star-forming galaxies is known to decline from redshifts of about 6 to 10, but a key question is the extent of star formation at even earlier times, corresponding to the period when the first galaxies might have emerged. Here we present spectroscopic observations of MACS1149-JD1, a gravitationally lensed galaxy observed when the Universe was less than four per cent of its present age. We detect an emission line of doubly ionized oxygen at a redshift of $9.1096\pm0.0006$, with an uncertainty of one standard deviation. This precisely determined redshift indicates that the red rest-frame optical colour arises from a dominant stellar component that formed about 250 million years after the Big Bang, corresponding to a redshift of about 15. Our results indicate the it may be possible to detect such early episodes of star formation in similar galaxies with future telescopes.

We performed coordinated observations of AR 12205, which produced a C-class flare on 2014 November 11, with the Interface Region Imaging Spectrograph (IRIS) and the Domeless Solar Telescope (DST) at Hida Observatory. Using spectral data in the Si IV 1403 \AA, C II 1335 \AA, and Mg II h and k lines from IRIS and the Ca II K, Ca II 8542 \AA, and H$\alpha$ lines from DST, we investigated a moving flare kernel during the flare. In the Mg II h line, the leading edge of the flare kernel showed the intensity enhancement in the blue wing, and the smaller intensity of the blue-side peak (h2v) than that of the red-side one (h2r). The blueshift lasted for 9-48 s with a typical speed of 10.1 $\pm$ 2.6 km s$^{-1}$ and it was followed by the high intensity and the large redshift with a speed of up to 51 km s$^{-1}$ detected in the Mg II h line. The large redshift was a common property for all six lines but the blueshift prior to it was found only in the Mg II lines. A cloud modeling of the Mg II h line suggests that the blue wing enhancement with such peak difference can be caused by a chromospheric-temperature (cool) upflow. We discuss a scenario in which an upflow of cool plasma is lifted up by expanding hot plasma owing to the deep penetration of non-thermal electrons into the chromosphere. Furthermore, we found that the blueshift persisted without any subsequent redshift in the leading edge of the flare kernel during its decaying phase. The cause of such long-lasting blueshift is also discussed.

Sunspots are the most notable structure on the solar surface with strong magnetic fields. The field is generally strongest in a dark area (umbra), but sometimes stronger fields are found in non-dark regions such as a penumbra and a light bridge. The formation mechanism of such strong fields outside umbrae is still puzzling. Here we report clear evidence of the magnetic field of 6,250 G, which is the strongest field among Stokes I profiles with clear Zeeman splitting ever observed on the Sun. The field was almost parallel to the solar surface and located in a bright region sandwiched by two opposite-polarity umbrae. Using a time series of spectral datasets, we discussed the formation process of the super-strong field and suggested that this strong field region was generated as a result of compression of one umbra pushed by the horizontal flow from the other umbra, like the subduction of the Earth's crust in plate tectonics.

We explore the effect of different galactic disc environments on the properties of star-forming clouds through variations in the background potential in a set of isolated galaxy simulations. Rising, falling and flat rotation curves expected in halo dominated, disc dominated and Milky Way-like galaxies were considered, with and without an additional two-arm spiral potential. The evolution of each disc displayed notable variations that are attributed to different regimes of stability, determined by shear and gravitational collapse. The properties of a typical cloud were largely unaffected by the changes in rotation curve, but the production of small and large cloud associations was strongly dependent on this environment. This suggests that while differing rotation curves can influence where clouds are initially formed, the average bulk properties are effectively independent of the global environment. The addition of a spiral perturbation made the greatest difference to cloud properties, successfully sweeping the gas into larger, seemingly unbound, extended structures and creating large arm-interarm contrasts.

We investigate the impact of radiation pressure on spatial dust distribution inside H$_\mathrm{II}$ regions using one-dimensional radiation hydrodynamic simulations, which include absorption and re-emission of photons by dust. In order to investigate grain size effects as well, we introduce two additional fluid components describing large and small dust grains in the simulations. Relative velocity between dust and gas strongly depends on the drag force. We include collisional drag force and coulomb drag force. We find that, in a compact H$_\mathrm{II}$ region, a dust cavity region is formed by radiation pressure. Resulting dust cavity sizes (~0.2 pc) agree with observational estimates reasonably well. Since dust inside an H$_\mathrm{II}$ region is strongly charged, relative velocity between dust and gas is mainly determined by the coulomb drag force. Strength of the coulomb drag force is about 2-order of magnitude larger than that of the collisional drag force. In addition, in a cloud of mass $10^5$ $M_{\odot}$, we find that the radiation pressure changes the grain size distribution inside H$_\mathrm{II}$ regions. Since large (0.1 $\mathrm{\mu m}$) dust grains are accelerated more efficiently than small (0.01 $\mathrm{\mu m}$) grains, the large to small grain mass ratio becomes smaller by an order of magnitude compared with the initial one. Resulting dust size distributions depend on the luminosity of the radiation source. The large and small grain segregation becomes weaker when we assume stronger radiation source, since dust grain charges become larger under stronger radiation and hence coulomb drag force becomes stronger.

Hypervelocity ejection of material by impact spallation is considered a plausible mechanism for material exchange between two planetary bodies. We have modeled the spallation process during vertical impacts over a range of impact velocities from 6 to 21 km/s using both grid- and particle-based hydrocode models. The Tillotson equations of state, which are able to treat the nonlinear dependence of density on pressure and thermal pressure in the strongly shocked matter, were used to study the hydrodynamic and thermodynamic response after impacts. The effects of material strength and gravitational acceleration were not considered. A two-dimensional time-dependent pressure field within a 1.5-fold projectile radius from the impact point was investigated in cylindrical coordinates to address the generation of spalled material. A resolution test was also performed to reject ejected materials with peak pressures that were too low due to artificial viscosity. The relationship between ejection velocity veject and peak pressure Ppeak was also derived. Our approach shows that late stage acceleration in an ejecta curtain occurs due to the compressible nature of the ejecta, resulting in an ejection velocity that can be higher than the ideal maximum of the resultant particle velocity after passage of a shock wave. We also calculate the ejecta mass that can escape from a planet like Mars (i.e., veject higher than 5 km/s) that matches the petrographic constraints from Martian meteorites, and which occurs when Ppeak from 30-50 GPa. Although the mass of such ejecta is limited to from 0.1-1 percent of the projectile mass in vertical impacts, this is sufficient for spallation to have been a plausible mechanism for the ejection of Martian meteorites. Finally, we propose that impact spallation is a plausible mechanism for the generation of tektites.

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 have investigated the impact of photoionization and radiation pressure on a dusty star-forming cloud using one-dimensional radiation hydrodynamic simulations, which include absorption and re-emission of photons by dust. We find that, in a cloud of mass 10^5 Msun and radius 17 pc, the effect of radiation pressure is negligible when star formation efficiency is 2%. The importance of radiation pressure increases with increasing star formation efficiency or an increasing dust-to-gas mass ratio. The net effect of radiation feedback, however, becomes smaller with the increasing dust-to-gas mass ratio, since the absorption of ultra-violet photons by dust grains suppresses photoionization and hence photoheating.

Massive early-type galaxies have higher metallicities and higher ratios of $\alpha$ elements to iron than their less massive counterparts. Reproducing these correlations has long been a problem for hierarchical galaxy formation theory, both in semi-analytic models and cosmological hydrodynamic simulations. We show that a simulation in which gas cooling in massive dark haloes is quenched by radio-mode active galactic nuclei (AGNs) feedback naturally reproduces the observed trend between $\alpha$/Fe and the velocity dispersion of galaxies, $\sigma$. The quenching occurs earlier for more massive galaxies. Consequently, these galaxies complete their star formation before $\alpha$/Fe is diluted by the contribution from type Ia supernovae. For galaxies more massive than $\sim 10^{11}~M_\odot$ whose $\alpha$/Fe correlates positively with stellar mass, we find an inversely correlated mass-metallicity relation. This is a common problem in simulations in which star formation in massive galaxies is quenched either by quasar- or radio-mode AGN feedback. The early suppression of gas cooling in progenitors of massive galaxies prevents them from recapturing enriched gas ejected as winds. Simultaneously reproducing the [$\alpha$/Fe]-$\sigma$ relation and the mass-metallicity relation is, thus, difficult in the current framework of galaxy formation.

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.

Fine-structure dynamics in solar prominences holds critical clues to understanding their physical nature of significant space-weather implications. We report evidence of rotational motions of horizontal helical threads in two active-region prominences observed by the \emph{Hinode} and/or \emph{IRIS} satellites at high resolution. In the first event, we found transverse motions of brightening threads at speeds up to 55~km~s$^{-1}$ seen in the plane of the sky. Such motions appeared as sinusoidal space--time trajectories with a typical period of $\sim$390~s, which is consistent with plane-of-sky projections of rotational motions. Phase delays at different locations suggest propagation of twists along the threads at phase speeds of 90--270~km~s$^{-1}$. At least 15 episodes of such motions occurred in two days, none associated with any eruption. For these episodes, the plane-of-sky speed is linearly correlated with the vertical travel distance, suggestive of a constant angular speed. In the second event, we found Doppler velocities of 30--40~km~s$^{-1}$ in opposite directions in the top and bottom portions of the prominence, comparable to the plane-of-sky speed. The moving threads have about twice broader line widths than stationary threads. These observations, when taken together, provide strong evidence for rotations of helical prominence threads, which were likely driven by unwinding twists triggered by magnetic reconnection between twisted prominence magnetic fields and ambient coronal fields.

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.

The physical properties and elemental abundances of the interstellar medium in galaxies during cosmic reionization are important for understanding the role of galaxies in this process. We report the Atacama Large Millimeter/submillimeter Array detection of an oxygen emission line at a wavelength of 88 micrometers from a galaxy at an epoch about 700 million years after the Big Bang. The oxygen abundance of this galaxy is estimated at about one-tenth that of the Sun. The non-detection of far-infrared continuum emission indicates a deficiency of interstellar dust in the galaxy. A carbon emission line at a wavelength of 158 micrometers is also not detected, implying an unusually small amount of neutral gas. These properties might allow ionizing photons to escape into the intergalactic medium.

We have performed very large and high resolution cosmological hydrodynamic simulations in order to investigate detectability of nebular lines in the rest-frame ultraviolet (UV) to optical wavelength range from galaxies at $z > 7$. We find that the expected line fluxes are very well correlated with apparent UV magnitudes. The C IV $1549 \rm \AA$ and the C III] $1909 \rm \AA$ lines of galaxies brighter than $26~ \rm AB$ are detectable with current facilities such as the VLT/XShooter and the Keck/MOSFIRE. Metal lines such as C IV $1549 \rm \AA$, C III] $1909 \rm \AA$, [O II] $3727 \rm \AA$ and [O III] $4959/5007 \rm \AA$ can be good targets for the spectroscopic observation with Thirty Meter Telescope (TMT), European Extremely Large Telescope (E-ELT), Giant Magellan Telescope (GMT) and James Webb Space Telescope (JWST). H$\alpha$ and H$\beta$ lines are also expected to be detectable with these telescopes. We also predict detectability of the nebular lines for $z > 10$ galaxies that will be found with JWST, Wide-Field Infrared Survey Telescope (WFIRST) and First Light And Reionization Explorer (FLARE) $(11 \leq z \leq 15$). We conclude that the C IV $1549 \rm \AA$, C III] $1909 \rm \AA$, [O III] $4959/5007 \rm \AA$ and H$\beta$ lines even from $z \sim 15$ galaxies can be strong targets for TMT, ELT and JWST. We also find that the magnification by gravitational lensing is of great help to detect such high-$z$ galaxies. According to our model, the C III] $1909 \rm \AA$ line in $z > 9$ galaxy candidates is detectable even using the current facilities.

Transverse magnetohydrodynamic (MHD) waves have been shown to be ubiquitous in the solar atmosphere and can in principle carry sufficient energy to generate and maintain the Sun's million-degree outer atmosphere or corona. However, direct evidence of the dissipation process of these waves and subsequent heating has not yet been directly observed. Here we report on high spatial, temporal, and spectral resolution observations of a solar prominence that show a compelling signature of so-called resonant absorption, a long hypothesized mechanism to efficiently convert and dissipate transverse wave energy into heat. Aside from coherence in the transverse direction, our observations show telltale phase differences around 180 degrees between transverse motions in the plane-of-sky and line-of-sight velocities of the oscillating fine structures or threads, and also suggest significant heating from chromospheric to higher temperatures. Comparison with advanced numerical simulations support a scenario in which transverse oscillations trigger a Kelvin-Helmholtz instability (KHI) at the boundaries of oscillating threads via resonant absorption. This instability leads to numerous thin current sheets in which wave energy is dissipated and plasma is heated. Our results provide direct evidence for wave-related heating in action, one of the candidate coronal heating mechanisms.

Transverse magnetohydrodynamic (MHD) waves are ubiquitous in the solar atmosphere and may be responsible for generating the Sun's million-degree outer atmosphere. However, direct evidence of the dissipation process and heating from these waves remains elusive. Through advanced numerical simulations combined with appropriate forward modeling of a prominence flux tube, we provide the observational signatures of transverse MHD waves in prominence plasmas. We show that these signatures are characterized by thread-like substructure, strong transverse dynamical coherence, an out-of-phase difference between plane-of-the-sky motions and LOS velocities, and enhanced line broadening and heating around most of the flux tube. A complex combination between resonant absorption and Kelvin-Helmholtz instabilities (KHI) takes place in which the KHI extracts the energy from the resonant layer and dissipates it through vortices and current sheets, which rapidly degenerate into turbulence. An inward enlargement of the boundary is produced in which the turbulent flows conserve the characteristic dynamics from the resonance, therefore guaranteeing detectability of the resonance imprints. We show that the features described in the accompanying paper (Okamoto et al. 2015) through coordinated Hinode and IRIS observations match well the numerical results.

We present a new numerical scheme to solve the transfer of diffuse radiation on three-dimensional mesh grids which is efficient on processors with highly parallel architecture such as recently popular GPUs and CPUs with multi- and many-core architectures. The scheme is based on the ray-tracing method and the computational cost is proportional to $N_{\rm m}^{5/3}$ where $N_{\rm m}$ is the number of mesh grids, and is devised to compute the radiation transfer along each light-ray completely in parallel with appropriate grouping of the light-rays. We find that the performance of our scheme scales well with the number of adopted CPU cores and GPUs, and also that our scheme is nicely parallelized on a multi-node system by adopting the multiple wave front scheme, and the performance scales well with the amount of the computational resources. As numerical tests to validate our scheme and to give a physical criterion for the angular resolution of our ray-tracing scheme, we perform several numerical simulations of the photo-ionization of neutral hydrogen gas by ionizing radiation sources without the "on-the-spot" approximation, in which the transfer of diffuse radiation by radiative recombination is incorporated in a self-consistent manner.

We investigate the evolution of two bars formed in fully self-consistent hydrodynamic simulations of the formation of Milky Way-mass galaxies. One galaxy shows higher central mass concentration and has a longer and stronger bar than the other at $z = 0$. The stronger bar evolves by transferring its angular momentum mainly to the dark halo. Consequently the rotation speed of the bar decreases with time, while the amplitude of the bar increases with time. These features qualitatively agree with the results obtained by idealized simulations. The pattern speed of the stronger bar largely goes up and down within a half revolution in its early evolutionary stage. These oscillations occur when the bar is misaligned with the $m = 4$ mode Fourier component. These oscillations correlate with the oscillations in the triaxilality of the dark matter halo, but differently from the way identified by idealized simulations. The amplitude of the weaker bar does not increase despite the fact that its rotation slows down with time.This result contradicts what is expected from idealized simulations and is caused by the decline of the central density associated with the mass loss and feedback from the stellar populations. The amplitude of the weaker bar is further weakens by the angular momentum injection by the interactions with stellar clumps in the disk. In the both galaxies, the bars are terminated around the 4:1 resonance.

We have investigated effects of dust attenuation on quasar luminosity functions using a semi-analytic galaxy formation model combined with a large cosmological N-body simulation. We estimate the dust attenuation of quasars self-consistently with that of galaxies by considering the dust in their host bulges. We find that the luminosity of the bright quasars is strongly dimmed by the dust attenuation, about 2 mag in the B-band. Assuming the empirical bolometric corrections for active galactic nuclei (AGNs) by Marconi et al., we find that this dust attenuation is too strong to explain the B-band and X-ray quasar luminosity functions simultaneously. We consider two possible mechanisms that weaken the dust attenuation. As such a mechanism, we introduce a time delay for AGN activity, that is, gas fueling to a central black hole starts some time after the beginning of the starburst induced by a major merger. The other is the anisotropy in the dust distribution. We find that in order to make the dust attenuation of the quasars negligible, either the gas accretion into the black holes has to be delayed at least three times the dynamical timescale of their host bulges or the dust covering factor is as small as 0.1.

Spectroscopic observations of prominence eruptions associated with coronal mass ejections (CMEs), although relatively rare, can provide valuable plasma and 3D geometry diagnostics. We report the first observations by the Interface Region Imaging Spectrograph (IRIS) mission of a spectacular fast CME/prominence eruption associated with an equivalent X1.6 flare on 2014 May 9. The maximum plane-of-sky and Doppler velocities of the eruption are 1200 and 460 km/s, respectively. There are two eruption components separated by ~200 km/s in Doppler velocity: a primary, bright component and a secondary, faint component, suggesting a hollow, rather than solid, cone-shaped distribution of material. The eruption involves a left-handed helical structure undergoing counter-clockwise (viewed top-down) unwinding motion. There is a temporal evolution from upward eruption to downward fallback with less-than-free-fall speeds and decreasing nonthermal line widths. We find a wide range of Mg II k/h line intensity ratios (less than ~2 expected for optically-thin thermal emission): the lowest ever-reported median value of 1.17 found in the fallback material and a comparably high value of 1.63 in nearby coronal rain and intermediate values of 1.53 and 1.41 in the two eruption components. The fallback material exhibits a strong ($> 5 \sigma$) linear correlation between the k/h ratio and the Doppler velocity as well as the line intensity. We demonstrate that Doppler dimming of scattered chromospheric emission by the erupted material can potentially explain such characteristics.

We present cosmological hydrodynamic simulations performed to study evolution of galaxy population. The simulations follow timed release of mass, energy, and metals by stellar evolution and employ phenomenological treatments of supernova feedback, pre-supernova feedback modeled as feedback by radiation pressure from massive stars, and quenching of gas cooling in large halos. We construct the fiducial model so that it reproduces the observationally estimated galaxy stellar mass functions and the relationships between the galaxy stellar mass and the host halo mass from $z = 4$ to 0. We find that the fiducial model constructed this way naturally explains the cosmic star formation history, the galaxy downsizing, and the star formation rate and metallicity of the star-forming galaxies. The simulations without the quenching of the gas cooling in large halos overproduce massive galaxies at $z < 2$ and fail to reproduce galaxy downsizing. The simulations that do not employ the radiation pressure feedback from young stars predict too strong redshift evolution of the mass-metallicity relation. Furthermore, the slope of the relation becomes too steep at low redshift without the radiation pressure feedback. The metallicity dependence in the radiation pressure feedback is a key to explain the observed mass-metallicity relation. These facts indicate that these two processes in addition to supernova feedback are essential for galaxy evolution. Our simple phenomenological model is suitable to construct a mock galaxy sample to study physical properties of observed galaxy populations.

We have performed a large cosmological hydrodynamics simulation tailored to the deep survey with the Hubble Space Telescope made in 2012, the so-called UDF12 campaign. After making a light-cone output, we have applied the same color selection criteria as the UDF12 campaign to select galaxies from our simulation, and then, have examined the physical properties of them as a proxy of the real observed UDF12 galaxies at $z > 7$. As a result, we find that the halo mass is almost linearly proportional to the observed ultraviolet (UV) luminosity ($4 \times 10^{11}~{\rm M_{\odot}}$ at $M_{\rm UV} = -21$). The dust attenuation and UV slope $\beta$ well correlates with the observed UV luminosity, which is consistent with observations quantitatively. The star formation rate (SFR) is also linearly proportional to the stellar mass and the specific SFR shows only a weak dependency on the mass. We also find an increasing star formation history with a time-scale of $\sim100$ Myr in the high-$z$ galaxies. An average metallicity weighted by the Lyman continuum luminosity reaches up to $>0.1$ Solar even at $z \sim 10$, suggesting a rapid metal enrichment. We also expect $\geq 0.1$ mJy at 350 GHz of the dust thermal emission from the galaxies with $H_{160} \leq 27$, which can be detectable with the Atacama Large Milimetre-submilimetre Array. The galaxies selected by the UDF12 survey contribute to only $52--12\%$ of the cosmic SFR density from $z \sim 7$ to $z \sim 10$, respectively. The James Webb Space Telescope will push the detection fraction up to $77--72\%$.

The Interface Region Imaging Spectrograph (IRIS) small explorer spacecraft provides simultaneous spectra and images of the photosphere, chromosphere, transition region, and corona with 0.33-0.4 arcsec spatial resolution, 2 s temporal resolution and 1 km/s velocity resolution over a field-of-view of up to 175 arcsec x 175 arcsec. IRIS was launched into a Sun-synchronous orbit on 27 June 2013 using a Pegasus-XL rocket and consists of a 19-cm UV telescope that feeds a slit-based dual-bandpass imaging spectrograph. IRIS obtains spectra in passbands from 1332-1358, 1389-1407 and 2783-2834 Angstrom including bright spectral lines formed in the chromosphere (Mg II h 2803 Angstrom and Mg II k 2796 Angstrom) and transition region (C II 1334/1335 Angstrom and Si IV 1394/1403 Angstrom). Slit-jaw images in four different passbands (C II 1330, Si IV 1400, Mg II k 2796 and Mg II wing 2830 Angstrom) can be taken simultaneously with spectral rasters that sample regions up to 130 arcsec x 175 arcsec at a variety of spatial samplings (from 0.33 arcsec and up). IRIS is sensitive to emission from plasma at temperatures between 5000 K and 10 MK and will advance our understanding of the flow of mass and energy through an interface region, formed by the chromosphere and transition region, between the photosphere and corona. This highly structured and dynamic region not only acts as the conduit of all mass and energy feeding into the corona and solar wind, it also requires an order of magnitude more energy to heat than the corona and solar wind combined. The IRIS investigation includes a strong numerical modeling component based on advanced radiative-MHD codes to facilitate interpretation of observations of this complex region. Approximately eight Gbytes of data (after compression) are acquired by IRIS each day and made available for unrestricted use within a few days of the observation.

We investigate the potential use of nebular emission lines in the rest-frame far-infrared (FIR) for determining spectroscopic redshift of z>8 galaxies with the Atacama Large Millimeter/sub-millimeter Array (ALMA). After making a line emissivity model as a function of metallicity, especially for the [O III] 88 micron line which is likely to be the strongest FIR line from H II regions, we predict the line fluxes from high-z galaxies based on a cosmological hydrodynamics simulation of galaxy formation. Since the metallicity of galaxies reaches at ~0.2 Zsun even at z>8 in our simulation, we expect the [O III] 88 micron line as strong as 1.3 mJy for 27 AB objects, which is detectable at a high significance by <1 hour integration with ALMA. Therefore, the [O III] 88 micron line would be the best tool to confirm the spectroscopic redshifts beyond z=8.

We investigated four major solar flare events that occurred in active regions NOAA 10930 (December 13 and 14, 2006) and NOAA 11158 (February 13 and 15, 2011) by using data observed by the Solar Optical Telescope (SOT) onboard the Hinode satellite. To reveal the trigger mechanism of solar flares, we analyzed the spatio-temporal correlation between the detailed magnetic field structure and the emission image of the Ca H line at the central part of flaring regions for several hours prior to the onset of flares. We observed in all the flare events that the magnetic shear angle in the flaring regions exceeded 70 degrees, as well as that characteristic magnetic disturbances developed at the centers of flaring regions in the pre-flare phase. These magnetic disturbances can be classified into two groups depending on the structure of their magnetic polarity inversion lines; the so-called "Opposite-Polarity" and "Reversed-Shear" magnetic field recently proposed by our group, although the magnetic disturbance in one event of the four samples is too subtle to clearly recognize the detailed structure. The result suggests that some major solar flares are triggered by rather small magnetic disturbances. We also show that the critical size of the flare-trigger field varies among flare events and briefly discuss how the flare-trigger process depends on the evolution of active regions.

(Abridged) We have developed a numerical software library for collisionless N-body simulations named "Phantom-GRAPE" which highly accelerates force calculations among particles by use of a new SIMD instruction set extension to the x86 architecture, AVX, an enhanced version of SSE. In our library, not only the Newton's forces, but also central forces with an arbitrary shape f(r), which has a finite cutoff radius r_cut (i.e. f(r)=0 at r>r_cut), can be quickly computed. Using an Intel Core i7--2600 processor, we measure the performance of our library for both the forces. In the case of Newton's forces, we achieve 2 x 10^9 interactions per second with 1 processor core, which is 20 times higher than the performance of an implementation without any explicit use of SIMD instructions, and 2 times than that with the SSE instructions. With 4 processor cores, we obtain the performance of 8 x 10^9 interactions per second. In the case of the arbitrarily shaped forces, we can calculate 1 x 10^9 and 4 x 10^9 interactions per second with 1 and 4 processor cores, respectively. The performance with 1 processor core is 6 times and 2 times higher than those of the implementations without any use of SIMD instructions and with the SSE instructions. These performances depend weakly on the number of particles. It is good contrast with the fact that the performance of force calculations accelerated by GPUs depends strongly on the number of particles. Substantially weak dependence of the performance on the number of particles is suitable to collisionless N-body simulations, since these simulations are usually performed with sophisticated N-body solvers such as Tree- and TreePM-methods combined with an individual timestep scheme. Collisionless N-body simulations accelerated with our library have significant advantage over those accelerated by GPUs, especially on massively parallel environments.

We adapt a modern scheme of smoothed particle hydrodynamics (SPH) to our tree N-body/SPH galactic chemodynamics code GCD+. The applied scheme includes imple- mentations of the artificial viscosity switch and artificial thermal conductivity pro- posed by Morris & Monaghan (1997), Rosswog & Price (2007) and Price (2008), to model discontinuities and Kelvin-Helmholtz instabilities more accurately. We first present hydrodynamics test simulations and contrast the results to runs undertaken without artificial viscosity switch or thermal conduction. In addition, we also explore the different levels of smoothing by adopting larger or smaller smoothing lengths, i.e. a larger or smaller number of neighbour particles, Nnb. We demonstrate that the new version of GCD+ is capable of modelling Kelvin-Helmholtz instabilities to a simi- lar level as the mesh code, Athena. From the Gresho vortex, point-like explosion and self-similar collapse tests, we conclude that setting the smoothing length to keep the number of neighbour particles as high as Nnb~58 is preferable to adopting smaller smoothing lengths. We present our optimised parameter sets from the hydrodynamics tests.

More than half of nearby disc galaxies have pseudobulges instead of classical bulges that are though to be end-products of galaxy mergers. Pseudobulges are presumed to develop overtime as a result of secular evolution of galaxy discs. We report simulations of galaxy formation, in which two disc galaxies with disky pseudobulges have formed. Based on the profile decomposition, the bulge-to-total mass ratio of the simulated galaxies is 0.6 for one galaxy and 0.3 for the other. We find that the main formation mechanism of the pseudobulges in our simulations is not the secular evolution of discs but high-redshift starbursts. The progenitors of the pseudobulges form as high-redshift discs with small scale lengths by rapid supply of low angular momentum gas. By redshift 2, before the main disc formation, pseudobulge formation has largely completed in terms of mass. The secular evolution such as bar instability accounts for about 30 % of the bulge mass for one galaxy and only 13 % for the other but does affect the final shape and kinematic properties of the pseudobulges.

We use large cosmological Smoothed-Particle-Hydrodynamics simulations to study the formation and evolution of sub-millimetre galaxies (SMGs). In our previous work, we studied the statistical properties of ultra-violet selected star-forming galaxies at high redshifts. We populate the same cosmological simulations with SMGs by calculating the reprocess of stellar light by dust grains into far-infrared to millimetre wavebands in a self-consistent manner. We generate light-cone outputs to compare directly the statistical properties of the simulated SMGs with available observations. Our model reproduces the submm source number counts and the clustering amplitude. We show that bright SMGs with flux $S > 1$ mJy reside in halos with mass of $\sim 10^{13} M_{\odot}$ and have stellar masses greater than $10^{11}\sim \rm M_{\odot}$. The angular cross-correlation between the SMGs and Lyman-$\alpha$ emitters is significantly weaker than that between the SMGs and Lyman-break galaxies. The cross-correlation is also weaker than the auto-correlation of the SMGs. The redshift distribution of the SMGs shows a broad peak at $z \sim 2$, where Bright SMGs contribute significantly to the global cosmic star formation rate density. Our model predicts that there are hundreds of SMGs with $S > 0.1$ mJy at $z > 5$ per 1 square degree field. Such SMGs can be detected by ALMA.

Bulges can be classified into classical and pseudobulges; the former are considered to be end products of galactic mergers and the latter to form via secular evolution of galactic disks. Observationally, bulges of disk galaxies are mostly pseudobulges, including the Milky Way's. We here show, by using self-consistent cosmological simulations of galaxy formation, that the formation of pseudobulges of Milky Way-sized disk galaxies has mostly completed before disk formation; thus the main channel of pseudobulge formation is not secular evolution of disks. Our pseudobulges form by rapid gas supply at high-redshift and their progenitors would be observed as high-redshift disks.

We compare the results of thirteen cosmological gasdynamical codes used to simulate the formation of a galaxy in the LCDM structure formation paradigm. The various runs differ in their hydrodynamical treatment (SPH, moving-mesh and AMR) but share the same initial conditions and adopt their latest published model of cooling, star formation and feedback. Despite the common halo assembly history, we find large code-to-code variations in the stellar mass, size, morphology and gas content of the galaxy at z=0, due mainly to the different implementations of feedback. Compared with observation, most codes tend to produce an overly massive galaxy, smaller and less gas-rich than typical spirals, with a massive bulge and a declining rotation curve. A stellar disk is discernible in most simulations, though its prominence varies widely from code to code. There is a well-defined trend between the effects of feedback and the severity of the disagreement with observation. Models that are more effective at limiting the baryonic mass of the galaxy come closer to matching observed galaxy scaling laws, but often to the detriment of the disk component. Our conclusions hold at two different numerical resolutions. Some differences can also be traced to the numerical techniques: more gas seems able to cool and become available for star formation in grid-based codes than in SPH. However, this effect is small compared to the variations induced by different feedback prescriptions. We conclude that state-of-the-art simulations cannot yet uniquely predict the properties of the baryonic component of a galaxy, even when the assembly history of its host halo is fully specified. Developing feedback algorithms that can effectively regulate the mass of a galaxy without hindering the formation of high-angular momentum stellar disks remains a challenge.

Ultraluminous infrared galaxies (ULIRGs) with multiple ($\ge 3$) nuclei are frequently observed. It has been suggested that these nuclei are produced by multiple major mergers of galaxies. The expected rate of such mergers is, however, too low to reproduce the observed number of ULIRGs with multiple nuclei. We have performed high-resolution simulations of the merging of two gas-rich disk galaxies. We found that extremely massive and compact star clusters form from the strongly disturbed gas disks after the first or second encounter between the galaxies. The mass of such clusters reaches $\sim 10^8 M_{\odot}$, and their half-mass radii are $20-30 \rm{pc}$. Since these clusters consist of young stars, they appear to be several bright cores in the galactic central region ($\sim \rm{kpc}$). The peak luminosity of these clusters reaches $\sim 10%$ of the total luminosity of the merging galaxy. These massive and compact clusters are consistent with the characteristics of the observed multiple nuclei in ULIRGs. Multiple mergers are not necessary to explain multiple nuclei in ULIRGs.

We present two types of numerical prescriptions that accelerate the radiative transfer calculation around point sources within a three-dimensional Cartesian grid by using the oct-tree structure for the distribution of radiation sources. In one prescription, distant radiation sources are grouped as a bright extended source when the group's angular size, $\theta_{\rm s}$, is smaller than a critical value, $\theta_{\rm crit}$, and radiative transfer is solved on supermeshes whose angular sizes are similar to that of the group of sources. The supermesh structure is constructed by coarse-graining the mesh structure. With this method, the computational time scales with $N_{\rm m} \log(N_{\rm m}) \log(N_{\rm s})$ where $N_{\rm m}$ and $N_{\rm s}$ are the number of meshes and that of radiation sources, respectively. While this method is very efficient, it inevitably overestimates the optical depth when a group of sources acts as an extended powerful radiation source and affects distant meshes. In the other prescription, a distant group of sources is treated as a bright point source ignoring the spatial extent of the group and the radiative transfer is solved on the meshes rather than the supermeshes. This prescription is simply a grid-based version of {\scriptsize START} by Hasegawa & Umemura and yields better results in general with slightly more computational cost ($\propto N_{\rm m}^{4/3} \log(N_{\rm s})$) than the supermesh prescription. Our methods can easily be implemented to any grid-based hydrodynamic codes and are well-suited to the adaptive mesh refinement methods.

We present a high-performance N-body code for self-gravitating collisional systems accelerated with the aid of a new SIMD instruction set extension of the x86 architecture: Advanced Vector eXtensions (AVX), an enhanced version of the Streaming SIMD Extensions (SSE). With one processor core of Intel Core i7-2600 processor (8 MB cache and 3.40 GHz) based on Sandy Bridge micro-architecture, we implemented a fourth-order Hermite scheme with individual timestep scheme (Makino and Aarseth, 1992), and achieved the performance of 20 giga floating point number operations per second (GFLOPS) for double-precision accuracy, which is two times and five times higher than that of the previously developed code implemented with the SSE instructions (Nitadori et al., 2006b), and that of a code implemented without any explicit use of SIMD instructions with the same processor core, respectively. We have parallelized the code by using so-called NINJA scheme (Nitadori et al., 2006a), and achieved 90 GFLOPS for a system containing more than N = 8192 particles with 8 MPI processes on four cores. We expect to achieve about 10 tera FLOPS (TFLOPS) for a self-gravitating collisional system with N 105 on massively parallel systems with at most 800 cores with Sandy Bridge micro-architecture. This performance will be comparable to that of Graphic Processing Unit (GPU) cluster systems, such as the one with about 200 Tesla C1070 GPUs (Spurzem et al., 2010). This paper offers an alternative to collisional N-body simulations with GRAPEs and GPUs.

We investigate the formation and evolution of satellite galaxies using smoothed particle hydrodynamics (SPH) simulations of a Milky Way(MW)-like system, focussing on the best resolved examples, analogous to the classical MW satellites. Comparing with a pure dark matter simulation, we find that the condensation of baryons has had a relatively minor effect on the structure of the satellites' dark matter halos. The stellar mass that forms in each satellite agrees relatively well over three levels of resolution (a factor of ~64 in particle mass) and scales with (sub)halo mass in a similar way in an independent semi-analytical model. Our model provides a relatively good match to the average luminosity function of the MW and M31. To establish whether the potential wells of our satellites are realistic, we measure their masses within observationally determined half-light radii, finding that the most massive examples have somewhat higher mass-to-light ratios than those derived for the MW dSphs from stellar kinematic data. A statistical test yields a ~9 percent probability that the simulated and observationally derived distributions of masses are consistent. Our results may suggest that either the MW halo is less massive than assumed in our simulations (~1.4e12 M_sun) or that there is substantial scatter in the satellite luminosity function or distribution of mass-to-light ratios at fixed host halo mass. Alternatively, feedback processes not properly captured by our simulations may have reduced the central densities of (sub)halos, or the subhalos may have initially formed with lower concentrations as would be the case, for example, if the dark matter were made of warm, rather than cold particles.

Alfv\'enic waves are thought to play an important role in coronal heating and acceleration of solar wind. Here we investigated the statistical properties of Alfv\'enic waves along spicules (jets that protrude into the corona) in a polar coronal hole using high cadence observations of the Solar Optical Telescope (SOT) onboard \emph{Hinode}. We developed a technique for the automated detection of spicules and high-frequency waves. We detected 89 spicules, and found: (1) a mix of upward propagating, downward propagating, as well as standing waves (occurrence rates of 59%, 21%, and 20%, respectively). (2) The phase speed gradually increases with height. (3) Upward waves dominant at lower altitudes, standing waves at higher altitudes. (4) Standing waves dominant in the early and late phases of each spicule, while upward waves were dominant in the middle phase. (5) In some spicules, we find waves propagating upward (from the bottom) and downward (from the top) to form a standing wave in the middle of the spicule. (6) The medians of the amplitude, period, and velocity amplitude were 55 km, 45 s, and 7.4 km s$^{-1}$, respectively. We speculate that upward propagating waves are produced near the solar surface (below the spicule) and downward propagating waves are caused by reflection of (initially) upward propagating waves off the transition region at the spicule top. The mix of upward and downward propagating waves implies that exploiting these waves to perform seismology of the spicular environment requires careful analysis and may be problematic.

We use very large cosmological Smoothed-Particle-Hydrodynamics simulations to study the properties of high redshift Lyman-alpha emitters (LAEs). We identify star-forming galaxies at z=3.1 in a cosmological volume of 100 h^-1Mpc on a side. We develop a phenomenological model of absorption, scattering and escape of Lyman-alpha photons on the assumption that the clumpiness of the inter-stellar medium in a galaxy is correlated with the larger scale substructure richness. The radiative transfer effect proposed by Neufeld (1991) allows a large fraction of Lyman-alpha photons to escape from a clumpy galaxy even if it contains a substantial amount of dust. Our model reproduces, for the first time, all of the following observed properties of LAEs at z=3.1: the angular correlation function, ultra-violet and Lyman-alpha luminosity functions, and the equivalent width distribution. A simple model that takes only dust absorption into account fails in matching the observational data, suggesting that the kind of effect we consider is needed. Our model also predicts a bimodal age distribution for LAEs. There are old, massive and dusty LAEs, similar to recently found high redshift LAEs. The large LAEs have escape fractions of Lyman-alpha photons of f_esc ~ 0.05-0.1.

We studied the formation process of star clusters using high-resolution N-body/smoothed particle hydrodynamcs simulations of colliding galaxies. The total number of particles is 1.2x10^8 for our high resolution run. The gravitational softening is 5 pc and we allow gas to cool down to \sim 10 K. During the first encounter of the collision, a giant filament consists of cold and dense gas found between the progenitors by shock compression. A vigorous starburst took place in the filament, resulting in the formation of star clusters. The mass of these star clusters ranges from 10^{5-8} Msun. These star clusters formed hierarchically: at first small star clusters formed, and then they merged via gravity, resulting in larger star clusters.

Continuous observations were performed of a quiescent prominence with the Solar Optical Telescope (SOT) on board the /emph{Hinode} satellite on 2006 December 23--24. A peculiar slowly-rising column of $/sim10^{4}$ K plasma develops from the lower atmosphere during the observations. The apparent ascent speed of the column is 2 km s$^{-1}$, while the fine structures of the column exhibit much faster motion of up to 20 km s$^{-1}$. The column eventually becomes a faint low-lying prominence. Associated with the appearance of the column, an overlying coronal cavity seen in the X-ray and EUV moves upward at $/sim$5 km s$^{-1}$. We discuss the relationship between these episodes, and suggest that they are due to the emergence of a helical flux rope that undergoes reconnection with lower coronal fields, possibly carrying material into the coronal cavity. Under the assumption of the emerging flux scenario, the lower velocity of 2 km s$^{-1}$ and the higher one of 20 km s$^{-1}$ in the column are attributed to the rising motion of the emerging flux and to the outflow driven by magnetic reconnection between the emerging flux and the pre-existing coronal field, respectively. The present paper gives a coherent explanation of the enigmatic phenomenon of the rising column with the emergence of the helical rope, and its effect on the corona. We discuss the implications that the emergence of such a helical rope has on the dynamo process in the convection zone.

We investigate the properties of satellite galaxies in cosmological N-body/SPH simulations of galaxy formation in Milky Way-sized haloes. Because of their shallow potential wells, satellite galaxies are very sensitive to heating processes which affect their gas content. Their properties can therefore be used to constrain the nature of feedback processes that regulate galaxy formation. In our simulations, we assume that all the energy produced by supernovae is used as kinetic energy to drive galactic winds. Several of our simulations produce bright, disc-dominated galaxies. We find that wind models in which the wind speed, v_w, is proportional to local velocity dispersion of dark matter, sigma, (and thus the wind mass-loading, eta_w \propto sigma^{-2}) have episodic star formation histories, reproduce the observed satellite luminosity function reasonably well (down to M_v=-7) and match the luminosity-metallicity relation observed in the Local Group satellites. By contrast, models that assume a constant wind speed overproduce faint satellites and predict an incorrect luminosity-metallicity relation. Our simulations therefore suggest that the feedback processes that operate on the scale of satellite galaxies should generate galactic outflows whose intensity varies inversely with the depth of the potential.

We investigate the magnitude and internal alignment of the angular momentum of cold dark matter haloes in simulations with and without baryons. We analyse the cumulative angular momentum profiles of hundreds of thousands of haloes in the Millennium simulation and in a smaller, but higher resolution, simulation, in total spanning 5 orders of magnitude in mass. For haloes of a given mass, the median specific angular momentum increases as j(<r) proportional to r. The direction of the vector varies considerably with radius: the median angle between the inner (< 0.25 Rvir) and total (< Rvir) angular momentum vectors is about 25degr. To investigate how baryons affect halo spin, we use another high resolution simulation, which includes gas cooling, star formation and feedback. This simulation produces a sample of galaxies with a realistic distribution of disc-to-total ratios. The formation of the galaxy spins up the dark matter within 0.1 Rvir such that the specific halo angular momentum increases by approx 50% in the median. The dark matter angular momentum becomes better aligned, but there remains a broad distribution of (mis-)alignments between the halo and the central galaxy, with a median angle between their angular momenta of ~ 30degr. Galaxies have a range of orientations relative to the shape of the halo: half of them have their minor axes misaligned by more than 45degr. Finally, we align a sample of haloes according to the orientation of their galaxies and stack the projected mass distributions. Although the individual haloes are aspherical, galaxy--halo misalignments produce a stacked mass distribution that is indistinguishable from circular. If the misalignment found in our simulations is realistic, it will be extremely difficult for weak lensing studies to measure the shape of dark haloes using this technique.

We study the formation histories and present-day structure of satellite galaxies formed in a high resolution hydrodynamic simulation of a Milky Way-like galaxy. The simulated satellites span nearly 4 orders of magnitude in luminosity but have a very similar mass within their inner 600 pc, ~ 3 10^7 M_solar, with very little scatter. This result is in agreement with the recent measurements for dwarf spheroidal galaxies (dSphs) in the Milky Way by Strigari et al. In our simulations a preferred mass scale arises naturally from the effects of the early reionisation of gas. These impose a sharp threshold of ~ 12 km/s on the circular velocity of haloes which can cool gas and make stars. At the present day, subhaloes that host satellites as luminous as the classical Milky Way dwarfs (L_V > 2.6 10^5 L_solar), have typically grown to have circular velocities of $> 20 km/s. There are, however, subhaloes with similar circular velocities today which were, nevertheless, below threshold at reionisation and thus remain dark. Star formation in above-threshold haloes is truncated when the halo is accreted into the main galaxy progenitor. Thus, many properties of today's dwarf satellites such as their luminosity and star formation rate are set by their accretion time.

The formation and evolution process and magnetic configuration of solar prominences remain unclear. In order to study the formation process of prominences, we examine continuous observations of a prominence in NOAA AR 10953 with the Solar Optical Telescope on the \emph{Hinode} satellite. As reported in our previous Letter, we find a signature suggesting that a helical flux rope emerges from below the photosphere under a pre-existing prominence. Here we investigate more detailed properties and photospheric indications of the emerging helical flux rope, and discuss their relationship to the formation of the prominence. Our main conclusions are: (1) A dark region with absence of strong vertical magnetic fields broadens and then narrows in Ca \textsc{ii} H-line filtergrams. This phenomenon is consistent with the emergence of the helical flux rope as photospheric counterparts. The size of the flux rope is roughly 30,000 km long and 10,000 km wide. The width is larger than that of the prominence. (2) No shear motion or converging flows are detected, but we find diverging flows such as mesogranules along the polarity inversion line. The presence of mesogranules may be related to the emergence of the helical flux rope. (3) The emerging helical flux rope reconnects with magnetic fields of the pre-existing prominence to stabilize the prominence for the next several days. We thus conjecture that prominence coronal magnetic fields emerge in the form of helical flux ropes that contribute to the formation and maintenance of the prominence.

We investigated the evolution of interacting disk galaxies using high-resolution $N$-body/SPH simulations, taking into account the multiphase nature of the interstellar medium (ISM). In our high-resolution simulations, a large-scale starburst occurred naturally at the collision interface between two gas disks at the first encounter, resulting in the formation of star clusters. This is consistent with observations of interacting galaxies. The probability distribution function (PDF) of gas density showed clear change during the galaxy-galaxy encounter. The compression of gas at the collision interface between the gas disks first appears as an excess at $n_{\rm H} \sim 10{\rm cm^{-3}}$ in the PDF, and then the excess moves to higher densities ($n_{\rm H} \gtrsim 100{\rm cm^{-3}}$) in a few times $10^7$ years where starburst takes place. After the starburst, the PDF goes back to the quasi-steady state. These results give a simple picture of starburst phenomena in galaxy-galaxy encounters.

The heating of the solar chromosphere and corona is a long-standing puzzle in solar physics. Hinode observations show the ubiquitous presence of chromospheric anemone jets outside sunspots in active regions. They are typically 3 to 7 arc seconds = 2000 to 5000 kilometers long and 0.2 to 0.4 arc second = 150 to 300 kilometers wide, and their velocity is 10 to 20 kilometers per second. These small jets have an inverted Y-shape, similar to the shape of x-ray anemone jets in the corona. These features imply that magnetic reconnection similar to that in the corona is occurring at a much smaller spatial scale throughout the chromosphere and suggest that the heating of the solar chromosphere and corona may be related to small-scale ubiquitous reconnection.

Hinode discovered a beautiful giant jet with both cool and hot components at the solar limb on 2007 February 9. Simultaneous observations by the Hinode SOT, XRT, and TRACE 195 satellites revealed that hot (5x10^6 K) and cool (10^4 K) jets were located side by side and that the hot jet preceded the associated cool jet (1-2 minutes). A current-sheet-like structure was seen in optical (Ca IIH), EUV (195A), and soft X-ray emissions, suggesting that magnetic reconnection is occurring in the transition region or upper chromosphere. Alfven waves were also observed with Hinode SOT. These propagated along the jet at velocities of 200 km/s with amplitudes (transverse velocity) of 5-15 km/s and a period of 200 s. We performed two-dimensional MHD simulation of the jets on the basis of the emerging flux-reconnection model, by extending Yokoyama and Shibata's model. We extended the model with a more realistic initial condition (10^6 K corona) and compared our model with multiwavelength observations. The improvement of the coronal temperature and density in the simulation model allowed for the first time the reproduction of the structure and evolution of both the cool and hot jets quantitatively, supporting the magnetic reconnection model. The generation and the propagation of Alfven waves are also reproduced self-consistently in the simulation model.

We perform cosmological hydrodynamic simulations to determine to what extent galaxies lose their gas due to photoheating from an ionizing background. We find that the characteristic mass at which haloes on average have lost half of their baryons is Mc ~ 6.5 x 10^9 Msun/h at z = 0, which corresponds to a circular velocity of 25 km/s. This is significantly lower than the filtering mass obtained by the linear theory, which is often used in semianalytical models of galaxy formation. We demonstrate it is the gas temperature at the virial radius which determines whether a halo can accrete gas. A simple model that follows the merger history of the dark matter progenitors, and where gas accretion is not allowed when this temperature is higher than the virial temperature of the halo, reproduces the results from the simulation remarkably well. This model can be applied to any reionization history, and is easy to incorporate in semianalytical models.

We investigate the evolution of angular momentum in simulations of galaxy formation in a cold dark matter universe. We analyse two model galaxies generated in the N-body/hydrodynamic simulations of Okamoto et al. Starting from identical initial conditions, but using different assumptions for the baryonic physics, one of the simulations produced a bulge-dominated galaxy and the other one a disc-dominated galaxy. The main difference is the treatment of star formation and feedback, both of which were designed to be more efficient in the disc-dominated object. We find that the specific angular momentum of the disc-dominated galaxy tracks the evolution of the angular momentum of the dark matter halo very closely: the angular momentum grows as predicted by linear theory until the epoch of maximum expansion and remains constant thereafter. By contrast, the evolution of the angular momentum of the bulge-dominated galaxy resembles that of the central, most bound halo material: it also grows at first according to linear theory, but 90% of it is rapidly lost as pre-galactic fragments, into which gas had cooled efficiently, merge, transferring their orbital angular momentum to the outer halo by tidal effects. The disc-dominated galaxy avoids this fate because the strong feedback reheats the gas which accumulates in an extended hot reservoir and only begins to cool once the merging activity has subsided. Our analysis lends strong support to the classical theory of disc formation whereby tidally torqued gas is accreted into the centre of the halo conserving its angular momentum.

We performed 3-dimensional N-body/SPH simulations to study how mass resolution and other model parameters such as the star formation efficiency parameter, C* and the threshold density, nth affect structures of the galactic gaseous/stellar disk in a static galactic potential. We employ 10^6 - 10^7 particles to resolve a cold and dense (T < 100 K & n_H > 100 cm^{-3}) phase. We found that structures of the ISM and the distribution of young stars are sensitive to the assumed nth. High-nth models with nth = 100 cm^{-3} yield clumpy multi-phase features in the ISM. Young stars are distributed in a thin disk of which half-mass scale height is 10 - 30 pc. In low-nth models with nth = 0.1 cm^{-3}, the stellar disk is found to be several times thicker, and the gas disk appears smoother than the high-nth models. A high-resolution simulation with high-nth is necessary to reproduce the complex structure of the gas disk. The global properties of the model galaxies in low-nth models, such as star formation histories, are similar to those in the high-nth models when we tune the value of C* so that they reproduce the observed relation between surface gas density and surface star formation rate density. We however emphasize that high-nth models automatically reproduce the relation, regardless of the values of C*. The ISM structure, phase distribution, and distributions of young star forming region are quite similar between two runs with values of C* which differ by a factor of 15. We also found that the timescale of the flow from n_H ~1 cm^{-3} to n_H > 100 cm^{-3} is about 5 times as long as the local dynamical time and is independent of the value of C*. The use of a high-nth criterion for star formation in high-resolution simulations makes numerical models fairy insensitive to the modelling of star formation. (Abridged)

Solar prominences are cool 10$^4$ Kelvin plasma clouds supported in the surrounding 10$^6$ Kelvin coronal plasma by as-yet undetermined mechanisms. Observations from \emph{Hinode} show fine-scale threadlike structures oscillating in the plane of the sky with periods of several minutes. We suggest these transverse magnetohydrodynamic waves may represent Alfv\'en waves propagating on coronal magnetic field lines and these may play a role in heating the corona.

Continuous observations were obtained of active region 10953 with the Solar Optical Telescope (SOT) on board the \emph{Hinode} satellite during 2007 April 28 to May 9. A prominence was located over the polarity inversion line (PIL) in the south-east of the main sunspot. These observations provided us with a time series of vector magnetic fields on the photosphere under the prominence. We found four features: (1) The abutting opposite-polarity regions on the two sides along the PIL first grew laterally in size and then narrowed. (2) These abutting regions contained vertically-weak, but horizontally-strong magnetic fields. (3) The orientations of the horizontal magnetic fields along the PIL on the photosphere gradually changed with time from a normal-polarity configuration to a inverse-polarity one. (4) The horizontal-magnetic field region was blueshifted. These indicate that helical flux rope was emerging from below the photosphere into the corona along the PIL under the pre-existing prominence. We suggest that this supply of a helical magnetic flux into the corona is associated with evolution and maintenance of active-region prominences.

In this paper, we present a new implementation of feedback due to active galactic nuclei (AGN) in cosmological simulations of galaxy formation. We assume that a fraction of jet energy, which is generated by an AGN, is transferred to the surrounding gas as thermal energy. Combining a theoretical model of mass accretion onto black holes with a multiphase description of star-forming gas, we self-consistently follow evolution of both galaxies and their central black holes. The novelty in our model is that we consider two distinct accretion modes: standard radiatively efficient thin accretion disks and radiatively inefficient accretion flows which we will generically refer to as RIAFs; motivated by theoretical models for jet production in accretion disks, we assume that only the RIAF is responsible for the AGN feedback. We find that, after an initial episode of bursting star formation, the accretion rate onto the central black hole drops so that the accretion disk switches to a RIAF structure. At this point, the feedback from the AGN becomes efficient and slightly suppresses star formation in the galactic disk and almost completely halts star formation in the bulge. As a result, the nucleus becomes a stochastically fuelled low-luminosity AGN (Seyfert galaxy) with recurrent short-lived episodes of activity after the star bursts. Our model predicts several properties of the low-luminosity AGN including the bolometric luminosity, jet powers, the effect on kpc-scale of the radio jet and the AGN lifetime, which are in broad agreement with observations of Seyfert galaxies and their radio activity. We also find that the mass ratios between the central black hole and the the host spheroid at z = 0 are ~10^{-3} regardless of the strength of either supernova feedback or AGN feedback. (abridged)

I review the current status of theoretical studies of galaxy formation. I outline the importance of the physics of baryonic component in galaxy formation by showing results obtained by using two major tools, semi-analytical approaches and cosmological simulations. In particular, I emphasis on roles of feedback in galaxy formation and discuss whether apparent conflictions between the standard theory of structure formation, the cold dark matter model, and observations can be solved by the feedback. I also discuss future prospects in numerical simulations of galaxy formation.

On 2005 September 13 a filament eruption accompanied by a halo CME occurred in the most flare-productive active region NOAA 10808 in Solar Cycle 23. Using multi-wavelength observations before the filament eruption on Sep. 13th, we investigate the processes leading to the catastrophic eruption. We find that the filament slowly ascended at a speed of 0.1km/s over two days before the eruption. During slow ascending, many small flares were observed close to the footpoints of the filament, where new magnetic elements were emerging. On the basis of the observational facts we discuss the triggering mechanism leading to the filament eruption. We suggest the process toward the eruption as follows: First, a series of small flares played a role in changing the topology of the loops overlying the filament. Second, the small flares gradually changed the equilibrium state of the filament and caused the filament to ascend slowly over two days. Finally, a C2.9 flare that occurred when the filament was close to the critical point for loss of equilibrium directly led to the catastrophic filament eruption right after itself.

We describe a numerical implementation of star formation in disk galaxies, in which the conversion of cooling gas to stars in the multiphase interstellar medium is governed by the rate at which molecular clouds are formed and destroyed. In the model, clouds form from thermally unstable ambient gas and get destroyed by feedback from massive stars and thermal conduction. Feedback in the ambient phase cycles gas into a hot galactic fountain or wind. We model the ambient gas hydrodynamically using smoothed particle hydrodynamics (SPH). However, we cannot resolve the Jeans mass in the cold and dense molecular gas and, therefore, represent the cloud phase with ballistic particles that coagulate when colliding. We show that this naturally produces a multiphase medium with cold clouds, a warm disk, hot supernova bubbles and a hot, tenuous halo. Our implementation of this model is based on the Gadget N-Body code. We illustrate the model by evolving an isolated Milky Way-like galaxy and study the properties of a disk formed in a rotating spherical collapse. Many observed properties of disk galaxies are reproduced well, including the molecular cloud mass spectrum, the molecular fraction as a function of radius, the Schmidt law, the stellar density profile and the appearance of a galactic fountain.

We investigate the properties of satellite galaxies formed in N-body/SPH simulations of galaxy formation in the $\Lambda$CDM cosmology. The simulations include the main physical effects thought to be important in galaxy formation and, in several cases, produce realistic spiral discs. In total, a sample of 9 galaxies of luminosity comparable to the Milky Way was obtained. At magnitudes brighter than the resolution limit, $M_V=-12$, the luminosity function of the satellite galaxies in the simulations is in excellent agreement with data for the Local Group. The radial number density profile of the model satellites, as well as their gas fractions also match observations very well. In agreement with previous N-body studies, we find that the satellites tend to be distributed in highly flattened configurations whose major axis is aligned with the major axis of the (generally triaxial) dark halo. In 2 out of 3 systems with sufficiently large satellite populations, the satellite system is nearly perpendicular to the plane of the galactic disc, a configuration analogous to that observed in the Milk Way. The discs themselves are perpendicular to the minor axis of their host halos in the inner parts, and the correlation between the orientation of the galaxy and the shape of the halo persists even out to the virial radius. However, in one case the disc's minor axis ends up, at the virial radius, perpendicular to the minor axis of the halo. The angular momenta of the galaxies and their host halo tend to be well aligned.

We investigate chemical evolution in Milky Way-like galaxies based on the cold dark matter model in which cosmic structures form via hierarchical merging. We introduce chemical enrichment due to type Ia supernovae (SNe Ia) into the Mitaka semi-analytic galaxy formation model developed by Nagashima & Yoshii. For the first time we derive distributions of stellar metallicities and their ratios in Milky Way-like galaxies treating chemical enrichment due to SNe Ia in a hierarchical galaxy formation model self-consistently. As a first attempt, we assume all SNe Ia to have the same lifetime, and assume instantaneous recycling for type II supernovae (SNe II). We find that our model reproduces well the metal abundance ratio [O/Fe] against [Fe/H] and the {iron metallicity distribution function} in the solar neighborhood. This means that the so-called G-dwarf problem is resolved by the hierarchical formation of galaxies, and a gas infall term introduced in traditional monolithic collapse models to solve this problem is well explained by the mixture of some physical processes such as hierarchical merging of dark halos, gas cooling, energy feedback and injection of gas and metals into hot gas due to supernovae. Our model predicts more oxygen-enhanced stars in bulges at [Fe/H] $\simeq 0$ than in disks. This trend seems to be supported by recent observations while they have still uncertainties. More data in number and accuracy will provide independent and important constraints on galaxy formation. (abridged)

Tidal disruption of dark matter halos around proto-globular clusters in a halo of a small galaxy is studied in the context of the hierarchical clustering scenario by using semi-cosmological N-body/SPH simulations assuming the standard cold dark matter model ($\Omega_0 = 1$). Our analysis on formation and evolution of the galaxy and its substructures archives until $z = 2.0$. In such a high-redshift universe, the Einstein-de Sitter universe is still a good approximation for a recently favored $\Lambda$-dominated universe, and then our results does not depend on the choice of cosmology. In order to resolve small gravitationally-bound clumps around galaxies and consider radiative cooling below $T = 10^4 K$, we adopt a fine mass resolution ($m_{\rm SPH} = 1.12 \times 10^3 \Msun$). Because of the cooling, each clump immediately forms a `core-halo' structure which consists of a baryonic core and a dark matter halo. The tidal force from the host galaxy mainly strips the dark matter halo from clumps and, as a result, theses clumps get dominated by baryons. Once a clump is captured by the host halo, its mass drastically decreases each pericenter passage. At $z = 2$, more than half of the clumps become baryon dominated systems (baryon mass/total mass $> 0.5$). Our results support the tidal evolution scenario of the formation of globular clusters and baryon dominated dwarf galaxies in the context of the cold dark matter universe.

We have performed hydrodynamic simulations of galaxy formation in a LCDM universe. We have followed galaxy formation in a dark matter halo, chosen to have a relatively quiet recent merger history, using different models for star formation and feedback. In all cases, we have adopted a multi-phase description of the interstellar medium and modelled star formation in quiescent and burst modes. We have explored two triggers for starbursts: strong shocks and high gas density, allowing for the possibility that stars in the burst may form with a top-heavy initial mass function. We find that the final morphology of the galaxy is extremely sensitive to the modelling of star formation and feedback. Starting from identical initial conditions, galaxies spanning the wide range of Hubble types, with B-band disc-to-total luminosity ratios ranging from 0.25 to 0.9, can form in the same dark matter halo. Models in which starbursts are induced by high gas density (qualitatively similar to models in which feedback is produced by AGN) generate energetic winds and result in galaxies with an early-type morphology. Models in which the starbursts are induced by strong shocks lead to extended discs. In this case, the feedback associated with the bursts suppresses the collapse of baryons in small halos, helping to create a reservoir of hot gas that is available for cooling after $z = 1$, following the bulk of the dynamical activity that builds up the halo. This gas then cools to form a large, young stellar disc.

We investigate the metal enrichment of elliptical galaxies in the framework of hierarchical models of galaxy formation. The semi-analytical model we use in this Letter, which has been used to study the metal enrichment of the intracluster medium (ICM) by Nagashima et al., includes the effects of flows of gas and metals both into and out of galaxies and the processes of metal enrichment due to both type Ia and type II supernovae. We adopt a solar neighbourhood IMF for star formation in discs, but consider models in which starbursts have either a solar neighbourhood IMF or a top-heavy IMF. We find that the alpha-element abundance in ellipticals is consistent with observed values only if the top-heavy IMF is used. This result is consistent with our previous study on the metal enrichment of the ICM. We also discuss the abundance ratio of alpha elements to iron as a function of velocity dispersion and metallicity. We find that models with a top-heavy IMF match the alpha/Fe ratios observed in typical L_* ellipticals, but none of the models reproduce the observed increase of alpha/Fe with velocity dispersion.

Seven cosmological hydrodynamic simulations of disk galaxy formation are analyzed to determine the alignment of the disk within the dark matter halo and the internal structure of the halo. We find that the orientation of the outer halo, beyond ~0.1 r_vir, is unaffected by the presence of the disk. In contrast, the inner halo is aligned such that the halo minor axis aligns with the disk axis. The relative orientation of these two regions of the halo are uncorrelated. The alignment of the disk and inner halo appears to take place simultaneously through their joint evolution. The disconnect between these two regions of the halo should be taken into account when modelling tidal streams in the halos of disk galaxies and when calculating intrinsic alignments of disk galaxies based on the properties of dark matter halos.

We investigate environmental effects on evolution of bright cluster galaxies in a Lambda-dominated cold dark matter universe using a combination of dissipationless N-body simulations and a semi-analytic galaxy formation model. We incorporate effects of ram-pressure stripping (RPS) and minor merger-induced small starburst (minor burst) into our model. By considering minor burst, observed morphology-radius relation is successfully reproduced. When we do not consider minor burst, the RPS hardly increases the intermediate B/T population. In addition, the RPS and minor burst are not important for colours or star formation rates of galaxies in the cluster core if star formation time-scale is properly chosen, because the star formation is sufficiently suppressed by consumption of the cold gas. We also find that SF in bulge-dominated galaxies is mainly terminated by starburst induced by major mergers in all environments.

We investigate the ability to constrain oscillatory features in the primordial power spectrum using current and future cosmic microwave background observations. In particular, we study the observability of an oscillation arising from imprints of physics at the cut-off energy scale. We perform a likelihood analysis on the WMAP data set, and find that the current data set constrains the amplitude of the oscillations to be less than 0.77 at 2-sigma, consistent with a power spectrum without oscillations. In addition, we investigate the fundamental limitations in the measurement of oscillation parameters by studying the constraints from a cosmic variance limited experiment. We find that such an experiment is capable of constraining the amplitude of such oscillations to be below 0.005, implying that reasonable models with cut-off energy scales Lambda>200 H_infl are unobservable through the microwave background.

We study the physical limitations placed on CMB temperature and polarization measurements of the initial power spectrum by geometric projection, acoustic physics, gravitational lensing and the joint fitting of cosmological parameters. Detailed information on the spectrum is greatly assisted by polarization information and localized to the acoustic regime k = 0.02-0.2 Mpc^{-1} with a fundamental resolution of Delta k/k>0.05. From this study we construct principal component based statistics, which are orthogonal to cosmological parameters including the initial amplitude and tilt of the spectrum, that best probe deviations from scale-free initial conditions. These statistics resemble Fourier modes confined to the acoustic regime and ultimately can yield ~50 independent measurements of the power spectrum features to percent level precision. They are straightforwardly related to more traditional parameterizations such as the the running of the tilt and in the future can provide many statistically independent measurements of it for consistency checks. Though we mainly consider physical limitations on the measurements, these techniques can be readily adapted to include instrumental limitations or other sources of power spectrum information.

We investigate the evolution of angular momentum in SPH simulations of galaxy formation, paying particular attention to artificial numerical effects. We find that a cold gas disc forming in an ambient hot gas halo receives a strong hydrodynamic torque from the hot gas. By splitting the hydrodynamic force into artificial viscosity and pressure gradients, we find that the angular momentum transport is caused not by the artificial viscosity but by the pressure gradients. Using simple test simulations of shear flows, we conclude that the pressure gradient-based viscosity can be divided into two components: one due to the noisiness of SPH and the other to ram pressure. The former is problematic even with very high resolution because increasing resolution does not reduce the noisiness. On the other hand, the ram pressure effect appears only when a cold gas disc or sheet does not contain enough particles. In such a case, holes form in the disc or sheet, and then ram pressure from intra-hole hot gas, causes significant deceleration. In simulations of galactic disc formation, star formation usually decreases the number of cold gas particles, and hole formation leads to the fragmentation of the disc. To circumvent these problem, we modify the SPH algorithm, decoupling the cold from the hot gas phases, i.e. inhibiting the hydrodynamic interaction between cold and hot particles. This, a crude modelling of a multi-phase fluid in SPH cosmological simulations, leads to the formation of smooth extended cold gas discs and to better numerical convergence. The decoupling is applicable in so far as the self-gravitating gas disc with negligible external pressure is a good approximation for a cold gas disc. (abridged)

Gravitational lensing of the microwave background by the intervening dark matter mainly arises from large-angle fluctuations in the projected gravitational potential and hence offers a unique opportunity to study the physics of the dark sector at large scales. Studies with surveys that cover greater than a percent of the sky will require techniques that incorporate the curvature of the sky. We lay the groundwork for these studies by deriving the full sky minimum variance quadratic estimators of the lensing potential from the CMB temperature and polarization fields. We also present a general technique for constructing these estimators, with harmonic space convolutions replaced by real space products, that is appropriate for both the full sky limit and the flat sky approximation. This also extends previous treatments to include estimators involving the temperature-polarization cross-correlation and should be useful for next generation experiments in which most of the additional information from polarization comes from this channel due to sensitivity limitations.

We investigate how cluster morphology is affected by the cosmological constant in low-density universes. Using high-resolution cosmological N-body/SPH simulations of flat (\Omega_0 = 0.3, \lambda_0 = 0.7, \Lambda CDM) and open (\Omega_0 = 0.3, \lambda_0 = 0, OCDM) cold dark matter universes, we calculate statistical indicators to quantify the irregularity of the cluster morphologies. We study axial ratios, center shifts, cluster clumpiness, and multipole moment power ratios as indicators for the simulated clusters at z=0 and 0.5. Some of these indicators are calculated for both the X-ray surface brightness and projected mass distributions. In \Lambda CDM all these indicators tend to be larger than those in OCDM at z=0. This result is consistent with the analytical prediction of Richstone, Loeb, & Turner, that is, clusters in \Lambda CDM are formed later than in OCDM, and have more substructure at z=0. We make a Kolmogorov-Smirnov test on each indicator for these two models. We then find that the results for the multipole moment power ratios and the center shifts for the X-ray surface brightness are under the significance level (5%). We results also show that these two cosmological models can be distinguished more clearly at z=0 than z = 0.5 by these indicators.

We investigate environmental effects on evolution of bright cluster galaxies in a $\Lambda$ dominated cold dark matter universe using a combination of dissipationless N-body simulations and a semi-analytic galaxy formation model. The N-body simulations enable us to calculate orbits of galaxies in simulated clusters. Therefore we can incorporate stripping of cold gas from galactic disks by ram pressure from intracluster medium into our model. In this paper we study how ram pressure stripping (RPS) and small starburst induced by a minor merger affect colors, star formation rates, and morphologies of cluster galaxies. We find that the RPS is not important for colors and SFRs of galaxies in the cluster core if star formation time-scale is properly chosen, because the star formation is sufficiently suppressed by consumption of the cold gas in the disks. Then observed color and SFR gradients can be reproduced without the RPS. The small starburst triggered by a minor merger hardly affects the SFRs and colors of the galaxies as well. We also examine whether these two processes can resolve the known problem that the major merger-driven bulge formation scenario predict too few galaxies of intermediate bulge-to-total luminosity ratio (B/T) in clusters. When the minor burst is taken into account, the intermediate B/T population is increased and the observed morphology gradients in clusters are successfully reproduced. When the minor burst is considered, the RPS also plays an important role in formation of the intermediate B/T galaxies. We present redshift evolution of morphological fractions predicted by our models. The predicted number ratios of the intermediate B/T galaxies to the bulge-dominated galaxies show nearly flat or slightly increasing trends with increasing redshift.

We develop the formalism necessary to study four-point functions of the cosmic microwave background (CMB) temperature and polarization fields. We determine the general form of CMB trispectra, with the constraints imposed by the assumption of statistical isotropy of the CMB fields, and derive expressions for their estimators, as well as their Gaussian noise properties. We apply these techniques to initial non-Gaussianity of a form motivated by inflationary models. Due to the large number of four-point configurations, the sensitivity of the trispectra to initial non-Gaussianity approaches that of the temperature bispectrum at high multipole moment. These trispectra techniques will also be useful in the study of secondary anisotropies induced for example by the gravitational lensing of the CMB by the large scale structure of the universe.

Weak gravitational lensing by the intervening large-scale structure of the Universe induces high-order correlations in the cosmic microwave background (CMB) temperature and polarization fields. We construct minimum variance estimators of the intervening mass distribution out of the six quadratic combinations of the temperature and polarization fields. Polarization begins to assist in the reconstruction when E-mode mapping becomes possible on degree-scale fields, i.e. for an experiment with a noise level of ~40 uK-arcmin and beam of ~7', similar to the Planck experiment; surpasses the temperature reconstruction at ~26 uK-arcmin and 4'; yet continues to improve the reconstruction until the lensing B-modes are mapped to l ~ 2000 at ~0.3 uK-arcmin and 3'. Ultimately, the correlation between the E and B modes can provide a high signal-to-noise mass map out to multipoles of L ~ 1000, extending the range of temperature-based estimators by nearly an order of magnitude. We outline four applications of mass reconstruction: measurement of the linear power spectrum in projection to the cosmic variance limit out to L ~ 1000 (or wavenumbers 0.002 < k < 0.2 in h/Mpc), cross-correlation with cosmic shear surveys to probe the evolution of structure tomographically, cross-correlation of the mass and temperature maps to probe the dark energy, and the separation of lensing and gravitational wave B-modes.

We present a model to investigate the formation and evolution of cluster galaxies using cosmological N-body simulations. The N-body simulations are used to make merging history trees of dark halos. Gas cooling, star formation, supernova feedback and mergers of galaxies within dark halos are included by using simple prescriptions taken from semi-analytic models of galaxy formation. In this {\it Letter} we exam the merger-driven bulge formation model and represent the morphology-density relation of cluster galaxies at z = 0. We find that this morphological evolution model can explain the distribution of elliptical galaxies in the clusters well and cannot reproduce the distribution of S0 galaxies. This result suggests that the elliptical galaxies are mainly formed by the major mergers, while, in the S0 formation, the processes other than major mergers play important role.

Using cosmological N-body simulations of critical (SCDM) and open ($\Omega = 0.3$, OCDM) cold dark matter models, we investigate evolution of cluster galaxies. From our numerical simulation, we construct merging history trees of the galaxies. By following their merging history, we show major merger fractions of the galaxies in cluster forming regions is roughly proportional to (1+z)^{4.5} at low redshifts (z < 2) and has a steep peak at $z \simeq 2.5$ and $z \simeq 3$ in SCDM and OCDM, respectively. We also show that the cluster galaxies are affected by tidal interaction after the clusters formed. Because the formation redshift of the cluster in SCDM, z_{form} = 0.15, is much more recent than that of the cluster in OCDM, z_{form} = 1.6, the cluster galaxies in SCDM show more rapid evolution by tidal interactions from z = 0.5 than those in OCDM.

We investigate effects of time evolution of a rich cluster of galaxies on its member galactic halos in the standard cold dark matter (SCDM) universe using high resolution N-body simulations. We identify several hundred galactic halos within virial radius of our simulated cluster. We also find that a large number of halos have been tidally disrupted at z = 0. Therefore we improve a method of deriving merging history trees of galaxies taking account of tidally stripped galaxies.