We study the formation of massive Population III binary stars using a newly developed radiation hydrodynamics code with the adaptive mesh refinement and adaptive ray-tracing methods. We follow the evolution of a typical primordial star-forming cloud obtained from a cosmological hydrodynamics simulation. Several protostars form as a result of disk fragmentation and grow in mass by the gas accretion, which is finally quenched by the radiation feedback from the protostars. Our code enables us, for the first time, to consider the feedback by both the ionizing and dissociating radiation from the multiple protostars, which is essential for self-consistently determining their final masses. At the final step of the simulation, we observe a very wide ($\gtrsim 10^4\,\mathrm{au}$) binary stellar system consisting of $60$ and $70\,M_\odot$ stars. One of the member stars also has two smaller mass ($10\,M_\odot$) companion stars orbiting at $200$ and $800\,\mathrm{au}$, making up a mini-triplet system. Our results suggest that massive binary or multiple systems are common among Population III stars.

The first hydrostatic core, the first quasi-hydrostatic object formed during the star formation process, is still the observational missing link between the prestellar and protostellar phases, mainly due to its short lifetime. Although we have not established a clear method to identify this rare object, recent theoretical studies predict that the first core has millimeter continuum emission and low-velocity outflow with a wide opening angle. An extensive continuum/outflow survey toward a large number of $"$starless$"$ cores in nearby star-forming regions works as a pathfinder. We observed 32 prestellar cores in Taurus with an average density of $\gtrsim$10$^5$ cm$^{-3}$ in 1.3 mm continuum and molecular lines using the ALMA-ACA stand-alone mode. Among the targets, MC35-mm centered at one of the densest $"$starless$"$ cores in Taurus has blueshifted/redshifted wings in the $^{12}$CO (2-1) line, indicating that there is deeply embedded object driving molecular outflow. The observed velocities and sizes of the possible outflow lobes are 2-4 km s$^{-1}$, and $\sim$2 $\times$10$^3$ au, respectively, and the dynamical time is calculated to be $\sim$10$^3$ yr. In addition to this, the core is one of the strongest N$_2$D$^{+}$ (3-2) emitter in our sample. All the observed signatures do not conflict with any of the theoretical predictions about the first hydrostatic core so far, and thus MC35-mm is unique as the only first-core candidate in the Taurus molecular cloud.

We have performed survey-type observations in 1 mm continuum and molecular lines toward dense cores (32 prestellar + 7 protostellar) with an average density of $\gtrsim$10$^5$ cm$^{-3}$ in the Taurus molecular clouds using the ALMA-ACA (Atacama Compact Array) stand-alone mode with an angular resolution of 6.$''$5 ($\sim$900 au). The primary purpose of this study is to investigate the innermost part of dense cores toward understanding the initial condition of star formation. In the protostellar cores, contributions from protostellar disks dominate the observed continuum flux with a range of 35-90\% except for the very low-luminosity object. For the prestellar cores, we have successfully confirmed continuum emission from dense gas with a density of $\gtrsim$3 $\times$10$^5$ cm$^{-3}$ toward approximately one-third of the targets. Thanks to the lower spatial frequency coverage with the ACA-7 m array, the detection rate is significantly higher than that of the previous surveys, which have 0 or 1 continuum detected sources among large number of starless samples using the ALMA Main array. The statistical counting method tells us that the lifetime of the prestellar cores until protostar formation therein approaches the free-fall time as the density increases. Among the prestellar cores, at least two targets have possible internal substructures, which are detected in continuum emission with the size scale of $\sim$1000 au if we consider the molecular line (C$^{18}$O and N$_2$D$^{+}$) distributions. These results suggest that small scale fragmentation/coalescence processes within a $\lesssim$0.1 pc scale region occur to determine the final core mass leading to protostars before the dynamical collapse whose central density range is $\sim$(0.3-1) $\times$10$^6$ cm$^{-3}$.

We report ALMA Cycle 4 observations of the Class I binary protostellar system L1551 IRS 5 in the 0.9-mm continuum emission, C18O (J=3-2), OCS (J=28-27), and four other Band 7 lines. At ~0.07" (= 10 au) resolution in the 0.9 mm emission, two circumstellar disks (CSDs) associated with the binary protostars are separated from the circumbinary disk (CBD). The CBD is resolved into two spiral arms, one connecting to the CSD around the northern binary source, Source N, and the other to Source S. As compared to the CBD in the neighboring protobinary system L1551 NE, the CBD in L1551 IRS 5 is more compact (r ~150 au) and the m=1 mode of the spirals found in L1551 NE is less obvious in L1551 IRS 5. Furthermore, the dust and molecular-line brightness temperatures of CSDs and CBD reach >260 K and >100 K, respectively, in L1551 IRS 5, much hotter than those in L1551 NE. The gas motions in the spiral arms are characterized by rotation and expansion. Furthermore, the transitions from the CBD to the CSD rotations at around the L2 and L3 Lagrangian points and gas motions around the L1 point are identified. Our numerical simulations reproduce the observed two spiral arms and expanding gas motion as a result of gravitational torques from the binary, transitions from the CBD to the CSD rotations, and the gas motion around the L1 point. The higher temperature in L1551 IRS 5 likely reflects the inferred FU-Ori event.

We study star cluster formation in various environments with different metallicities and column densities by performing a suite of three-dimensional radiation hydrodynamics simulations. We find that the photoionization feedback from massive stars controls the star formation efficiency (SFE) in a star-forming cloud, and its impact sensitively depends on the gas metallicity $Z$ and initial cloud surface density $\Sigma$. At $Z=1~Z_{\odot}$, for instance, the SFE increases from 0.02 to 0.3 with $\Sigma$ from $10$ to $300~{\rm M_{\odot} {\rm pc^{-2}}}$. In low-metallicity cases $10^{-2}- 10^{-1} Z_{\odot}$, star clusters form from atomic warm gases because the molecule formation time is not short enough with respect to the cooling or dynamical time. In addition, the whole cloud is disrupted more easily by expanding HII bubbles which have higher temperature owing to less efficient cooling. With smaller dust attenuation, the ionizing radiation feedback from nearby massive stars is stronger and terminate star formation in dense clumps. These effects result in inefficient star formation in low-metallicity environments: the SFE is reduced by a factor of $\sim 0.3$ for the cases with $Z=10^{-2}~Z_{\odot}$ regardless of $\Sigma$. Newborn star clusters are also gravitationally less bound. We further develop a new semi-analytical model that can reproduce the simulation results well, particularly the observed dependencies of the SFEs on the cloud surface densities and metallicities.

The Gravity Collaboration detected a near-infrared hotspot moving around Sgr A* during the 2018 July 22 flare. They fitted the partial loop the hotspot made on the sky with a circular Keplerian orbit of radius $\simeq7.5\,r_{\rm g}$ around the supermassive black hole (BH), where $r_{\rm g}$ is the gravitational radius. However, because the hotspot traversed the loop in a short time, models in which the hotspot tracks the motion of some fluid element tend to produce a best-fit trajectory smaller than the observed loop. This is true for a circular Keplerian orbit, even when BH spin is accounted for, and for motion along a RIAF streamline. A marginally bound geodesic suffers from the same problem; in addition, it is not clear what the origin of an object following the geodesic would be. The observed hotspot motion is more likely a pattern motion. Circular motion with $r\simeq12.5\,r_{\rm g}$ and a super-Keplerian speed $\simeq0.8\,c$ is a good fit. Such motion must be pattern motion because it cannot be explained by physical forces. The pattern speed is compatible with magnetohydrodynamic perturbations, provided that the magnetic field is sufficiently strong. Circular pattern motion of radius $\sim20\, r_{\rm g}$ on a plane above the BH is an equally good alternative; in this case, the hotspot may be caused by a precessing outflow interacting with a surrounding disk. As all our fits have relatively large radii, we cannot constrain the BH spin using these observations.

We have observed the Class I protostar L1489 IRS with the Atacama Millimeter/submillimeter Array (ALMA) in Band 6. The C$^{18}$O $J=$2-1 line emission shows flattened and non-axisymmetric structures in the same direction as its velocity gradient due to rotation. We discovered that the C$^{18}$O emission shows dips at a radius of ~200-300 au while the 1.3 mm continuum emission extends smoothly up to r~400 au. At the radius of the C$^{18}$O dips, the rotational axis of the outer portion appears to be tilted by ~15 degrees from that of the inner component. Both the inner and outer components with respect to the C$^{18}$O dips exhibit the $r^{-0.5}$ Keplerian rotation profiles until r~600 au. These results not only indicate that a Keplerian disk extends up to ~600 au but also that the disk is warped. We constructed a three dimensional warped disk model rotating at the Keplerian velocity, and demonstrated that the warped disk model reproduces main observed features in the velocity channel maps and the PV diagrams. Such a warped disk system can form by mass accretion from a misaligned envelope. We also discuss a possible disk evolution scenario based on comparisons of disk radii and masses between Class I and Class II sources.

The magnetic field structure of a star-forming Bok globule BHR 71 was determined based on near-infrared polarimetric observations of background stars. The magnetic field in BHR 71 was mapped from 25 stars. By using a simple 2D parabolic function, the plane-of-sky magnetic axis of the core was found to be $\theta_{\rm mag} = 125^{\circ} \pm 11^{\circ}$. The plane-of-sky mean magnetic field strength of BHR 71 was found to be $B_{\rm pos} = 8.8 - 15.0$ $\mu$G, indicating that the BHR 71 core is magnetically supercritical with $\lambda = 1.44 - 2.43$. Taking into account the effect of thermal/turbulent pressure and the plane-of-sky magnetic field component, the critical mass of BHR 71 was $M_{\rm cr} = 14.5-18.7$ M$_{\odot}$, which is consistent with the observed core mass of $M_{\rm core} \approx 14.7$ M$_{\odot}$ (Yang et al. 2017). We conclude that BHR 71 is in a condition close to a kinematically critical state, and the magnetic field direction lies close to the plane of sky. Since BHR 71 is a star-forming core, a significantly subcritical condition (i.e., the magnetic field direction deviating from the plane of sky) is unlikely, and collapsed from a condition close to a kinematically critical state. There are two possible scenarios to explain the curved magnetic fields of BHR 71, one is an hourglass-like field structure due to mass accumulation and the other is the Inoue \& Fukui (2013) mechanism, which proposes the interaction of the core with a shock wave to create curved magnetic fields wrapping around the core.

Some short gamma-ray bursts (SGRBs) show a longer lasting emission phase, called extended emission (EE) lasting $\sim 10^{2-3}\,\rm s$, as well as a plateau emission (PE) lasting $\sim10^{4-5}\,\rm s$. While a long-lasting activity of the central engines is a promising explanation for powering both emissions, their physical origin and their emission mechanisms are still uncertain. In this work, we study the properties of the EEs and their connection with the PEs. First, we constrain the minimal Lorentz factor $\Gamma$ of the outflows powering EEs, using compactness arguments and find that the outflows should be relativistic, $\Gamma\gtrsim10$. We propose a consistent scenario for the PEs, where the outflow eventually catches up with the jet responsible for the prompt emission, injecting energy into the forward shock formed by the prior jet, which naturally results in a PE. We also derive the radiation efficiency of EEs and the Lorentz factor of the outflow within our scenario for 10 well-observed SGRBs accompanied by both EE and PE. The efficiency has an average value of $\sim3\,\%$ but shows a broad distribution ranging from $\sim0.01$ to $\sim100\%$. The Lorentz factor is $\sim20-30$, consistent with the compactness arguments. These results suggest that EEs are produced by a slower outflow via more inefficient emission than the faster outflow which causes the prompt emission with a high radiation efficiency.

In this study, the detailed magnetic field structure of the dense protostellar core Barnard 335 (B335) was revealed based on near-infrared polarimetric observations of background stars to measure dichroically polarized light produced by magnetically aligned dust grains in the core. Magnetic fields pervading B335 were mapped using 24 stars after subtracting unrelated ambient polarization components, for the first time revealing that they have an axisymmetrically distorted hourglass-shaped structure toward the protostellar core. On the basis of simple two- and three-dimensional magnetic field modeling, magnetic inclination angles in the plane-of-sky and line-of-sight directions were determined to be $90^{\circ} \pm 7^{\circ}$ and $50^{\circ} \pm 10^{\circ}$, respectively. The total magnetic field strength of B335 was determined to be $30.2 \pm 17.7$ $\mu {\rm G}$. The critical mass of B335, evaluated using both magnetic and thermal/turbulent support against collapse, was determined to be $M_{\rm cr} = 3.37 \pm 0.94$ ${\rm M}_{\odot}$, which is identical to the observed core mass of $M_{\rm core}=3.67$ M$_{\odot}$. We thus concluded that B335 started its contraction from a condition near equilibrium. We found a linear relationship in the polarization versus extinction diagram, up to $A_V \sim 15$ mag toward the stars with the greatest obscuration, which verified that our observations and analysis provide an accurate depiction of the core.

A $\gamma$-ray source must have a limited optical depth to pair production. This simple condition, called compactness, implies that gamma-ray bursts (GRBs) must involve a highly relativistic motion ($\Gamma \gtrsim 100$) giving the first and most important clue on their nature. So far, this condition has been discussed under the assumption that the $\gamma$-ray sources are viewed on-axis, that is, by an observer within the beaming cone of the relativistic source. Recently, following the detection of the weak short GRB 170817A, an extensive interest arose in the possibility that some $\gamma$-ray sources are viewed off-axis. We generalize here the compactness formalism for an arbitrary viewing angle taking several possible opacity processes and $\gamma$-ray spectra into account. We find that for a given observables (peak luminosity, temporal variability, and spectra) the minimal Lorentz factor, $\Gamma_{\rm min}$, is obtained, as expected, for an on-axis observer. More remarkably we find that compactness dictates also a maximal viewing angle, $\theta_{\rm max} \simeq 1/2\Gamma_{\rm min}$. Our limit implies for regular GRBs a very small allowed viewing angle ($\lesssim10^{-2}\,\rm rad$), making it extremely unlikely that they are viewed off-axis. For GRB 170817A we confirm earlier results that rule out the possibility that the observed $\gamma$-rays were seen by an on-axis observer as a regular short GRB. The short GRB 150101B was also suggested to be an off-axis event. We show that its maximal viewing angle $\lesssim0.05\,\rm rad$, which is inconsistent with the off-axis model. Finally we show that for low luminosity GRBs, compactness does not exclude by itself an off-axis model, but when combined with other consideration this option is strongly disfavored.

von Kienlin et al. (2019) selected 11 short gamma-ray bursts (sGRBs) whose characteristics are similar to GRB 170817A. These bursts, like GRB 170817A, have a hard spike followed by a soft thermal tail. However, as their redshifts are unknown it is not clear if their luminosities are as low as that of GRB 170817A. Comparing the positions in the $\epsilon_{\rm p}$-$E_{\rm\gamma,iso}$ (spectral peak energy - isotropic-equivalent energy) plane and using compactness arguments to estimate the minimal Lorentz factor, $\Gamma$, we find that all the bursts in this sample are consistent with being regular sGRBs if they are located at $z\simeq0.3-3$. They are also consistent with being similar to GRB 170817A if they are located at $z\lesssim0.1$. Even in the latter case, the events must involve at least mildly relativistic ($\Gamma \gtrsim 2$) motion within the sources. We, further, find that at most one or two bursts in the sample are consistent with the cocoon shock-breakout model. Finally, we calculate the event rate of either off-axis emission from a jet core or from a jet-wing (surrounding the core). We find that the off-axis emission model as an origin of the sample is rejected as it predicts too small event rate. The wing model can be consistent with the observed rate but the model parameters cannot be constrained by the current observations.

The detailed magnetic field structure of the starless dense core CB81 (L1774, Pipe 42) in the Pipe Nebula was determined based on near-infrared polarimetric observations of background stars to measure dichroically polarized light produced by magnetically aligned dust grains in the core. The magnetic fields pervading CB81 were mapped using 147 stars and axisymmetrically distorted hourglass-like fields were identified. On the basis of simple 2D and 3D magnetic field modeling, the magnetic inclination angles in the plane-of-sky and line-of-sight directions were determined to be $4^{\circ} \pm 8^{\circ}$ and $20^{\circ} \pm 20^{\circ}$, respectively. The total magnetic field strength of CB81 was found to be $7.2 \pm 2.3$ $\mu{\rm G}$. Taking into account the effects of thermal/turbulent pressure and magnetic fields, the critical mass of CB81 was calculated to be $M_{\rm cr}=4.03 \pm 0.40$ M$_{\odot}$, which is close to the observed core mass of $M_{\rm core}=3.37 \pm 0.51$ M$_{\odot}$. We thus conclude that CB81 is in a condition close to the critical state. In addition, a spatial offset of $92''$ was found between the center of magnetic field geometry and the dust extinction distribution; this offset structure could not have been produced by self-gravity. The data also indicate a linear relationship between polarization and extinction up to $A_V \sim 30$ mag going toward the core center. This result confirms that near-infrared polarization can accurately trace the overall magnetic field structure of the core.

Observational data for the hourglass-like magnetic field toward the starless dense core FeSt 1-457 were compared with a flux freezing magnetic field model (Myers et al. 2018). Fitting of the observed plane-of-sky magnetic field using the flux freezing model gave a residual angle dispersion comparable with the results based on a simple three-dimensional parabolic model. The best-fit parameters for the flux freezing model were a line-of-sight magnetic inclination angle of $\gamma_{\rm mag} = 35^{\circ} \pm 15^{\circ}$ and a core center to ambient (background) density contrast of $\rho_{\rm c} / \rho_{\rm bkg} = 75$. The initial density for core formation ($\rho_0$) was estimated to be $\rho_{\rm c} / 75 = 4670$ cm$^{-3}$, which is about one order of magnitude higher than the expected density ($\sim 300$ cm$^{-3}$) for the inter-clump medium of the Pipe Nebula. FeSt 1-457 is likely to have been formed from the accumulation of relatively dense gas, and the relatively dense background column density of $A_V \simeq 5$ mag supports this scenario. The initial radius (core formation radius) $R_0$ and the initial magnetic field strength $B_0$ were obtained to be 0.15 pc ($1.64 R$) and $10.8-14.6$ $\mu$G, respectively. We found that the initial density $\rho_0$ is consistent with the mean density of the nearly critical magnetized filament with magnetic field strength $B_0$ and radius $R_0$. The relatively dense initial condition for core formation can be naturally understood if the origin of the core is the fragmentation of magnetized filaments.

The magnetic field structure, kinematical stability, and evolutionary status of the starless dense core Barnard 68 (B68) are revealed based on the near-infrared polarimetric observations of background stars, measuring the dichroically polarized light produced by aligned dust grains in the core. After subtracting unrelated ambient polarization components, the magnetic fields pervading B68 are mapped using 38 stars and axisymmetrically distorted hourglass-like magnetic fields are obtained, although the evidence for the hourglass field is not very strong. On the basis of simple 2D and 3D magnetic field modeling, the magnetic inclination angles on the plane-of-sky and in the line-of-sight direction are determined to be $47^{\circ} \pm 5^{\circ}$ and $20^{\circ} \pm 10^{\circ}$, respectively. The total magnetic field strength of B68 is obtained to be $26.1 \pm 8.7$ $\mu {\rm G}$. The critical mass of B68, evaluated using both magnetic and thermal/turbulent support, is $M_{\rm cr} = 2.30 \pm 0.20$ ${\rm M}_{\odot}$, which is consistent with the observed core mass of $M_{\rm core}=2.1$ M$_{\odot}$, suggesting nearly critical state. We found a relatively linear relationship between polarization and extinction up to $A_V \sim 30$ mag toward the stars with deepest obscuration. Further theoretical and observational studies are required to explain the dust alignment in cold and dense regions in the core.

During a winter thunderstorm on 2017 November 24, a strong burst of gamma rays with energies up to $\sim$10~MeV was detected coincident with a lightning discharge, by scintillation detectors installed at Kashiwazaki-Kariwa Nuclear Power Station at sea level in Japan. The burst had a sub-second duration, which is suggestive of photoneutron productions. The leading part of the burst was resolved into four intense gamma-ray bunches, each coincident with a low-frequency radio pulse. These bunches were separated by 0.7--1.5~ms, with a duration of $\ll$1~ms each. Thus, the present burst may be considered as a ``downward" terrestrial gamma-ray flash (TGF), which is analogous to up-going TGFs observed from space. Although the scintillation detectors were heavily saturated by these bunches, the total dose associated with them was successfully measured by ionization chambers, employed by nine monitoring posts surrounding the power plant. From this information and Monte Carlo simulations, the present downward TGF is suggested to have taken place at an altitude of 2500 $\pm$ 500~m, involving $8^{+8}_{-4} \times 10^{18}$ avalanche electrons with energies above 1~MeV. This number is comparable to those in up-going TGFs.

We performed a fluctuation analysis of the Hubble Extremely Deep Field (XDF) at four optical wavelength bands and found large fluctuations that are significantly brighter than those expected for ordinary galaxies. Good cross-correlations with flat spectra are found down to 0.2 arcsec, indicating the existence of a spatial structure even at the 0.2 arcsec scale. The detected auto and cross-correlations provide a lower limit of 24 nW m$^{-2}$ sr$^{-1}$ for the absolute sky brightness at 700-900 nm, which is consistent with previous observations. We searched for candidate objects to explain the detected large fluctuation using the catalog of the Hubble Ultra Deep Field (UDF), and found that the surface number density of faint compact objects (FCOs) rapidly increases toward the faint end. Radial profiles of FCOs are indistinguishable from the PSF, and the effective radius based on de Vaucouleur's law is estimated to be smaller than 0.02 arcsec. The SEDs of FCOs follow a power law at optical wavelengths, but show greater emission and structure at $\lambda$ > 1 $\mu$m. Assuming that the FCOs are the cause of the excess brightness and fluctuations, the faint magnitude limit is 34.9 mag for the F775W band, and the surface number density reaches $2.6 \times 10^3$ (arcsec)$^{-2}$. Recent gamma-ray observations require that the redshift of FCOs must be less than 0.1, if FCOs are the origin of the excess optical and infrared background. Assuming that FCOs consist of missing baryons, the mass and luminosity of a single FCO range from $10^{2}$ to 1$0^{3}$ solar units, and mass-to-luminosity ratio is significantly lower than 1.0 solar unit. The maximum effective radius of an FCO is 4.7 pc. These results and the good correlation between near-infrared and X-ray background indicate that FCOs could be powered by the gravitational energy associated with black holes.

The formation scenario of brown dwarfs is still unclear because observational studies to investigate its initial condition are quite limited. Our systematic survey of nearby low-mass star-forming regions using the Atacama Compact Array (aka Morita array) and the IRAM 30 m telescope in 1.2 mm continuum has identified a centrally concentrated starless condensation with a central H$_2$ volume density of $\sim$10$^6$ cm$^{-3}$, MC5-N, connected to a narrow (width $\sim$0.03 pc) filamentary cloud in the Taurus L1495 region. The mass of the core is $\sim$0.2-0.4 $M_{\odot}$, which is an order of magnitude smaller than typical low-mass prestellar cores. Taking into account a typical core to star formation efficiency for prestellar cores ($\sim$20%-40%) in nearby molecular clouds, brown dwarf(s) or very low-mass star(s) may be going to be formed in this core. We have found possible substructures at the high-density portion of the core, although much higher angular resolution observation is needed to clearly confirm them. The subsequent N$_2$H$^+$ and N$_2$D$^+$ observations using the Nobeyama 45 m telescope have confirmed the high-deuterium fractionation ($\sim$30%). These dynamically and chemically evolved features indicate that this core is on the verge of proto-brown dwarf or very low-mass star formation and is an ideal source to investigate the initial conditions of such low-mass objects via gravitational collapse and/or fragmentation of the filamentary cloud complex.

A new Riemann solver is presented for the ideal magnetohydrodynamics (MHD) equations with the so-called Boris correction. The Boris correction is applied to reduce wave speeds, avoiding an extremely small timestep in MHD simulations. The proposed Riemann solver, Boris-HLLD, is based on the HLLD solver. As done by the original HLLD solver, (1) the Boris-HLLD solver has four intermediate states in the Riemann fan when left and right states are given, (2) it resolves the contact discontinuity, Alfv\'en waves, and fast waves, and (3) it satisfies all the jump conditions across shock waves and discontinuities except for slow shock waves. The results of a shock tube problem indicate that the scheme with the Boris-HLLD solver captures contact discontinuities sharply and it exhibits shock waves without any overshoot when using the minmod limiter. The stability tests show that the scheme is stable when $|u| \lesssim 0.5c$ for a low Alfv\'en speed ($V_A \lesssim c$), where $u$, $c$, and $V_A$ denote the gas velocity, speed of light, and Alfv\'en speed, respectively. For a high Alfv\'en speed ($V_A \gtrsim c$), where the plasma beta is relatively low in many cases, the stable region is large, $|u| \lesssim (0.6-1) c$. We discuss the effect of the Boris correction on physical quantities using several test problems. The Boris-HLLD scheme can be useful for problems with supersonic flows in which regions with a very low plasma beta appear in the computational domain.

We investigate the gas structures around young binary stars by using three-dimensional numerical simulations. Each model exhibits circumstellar disks, spiral arms, and a circumbinary disk with an inner gap or cavity. The circumbinary disk has an asymmetric pattern rotating at an angular velocity of approximately one-fourth of the binary orbit of the moderate-temperature models. Because of this asymmetry, the circumbinary disk has a density bump and a vortex, both of which continue to exist until the end of our calculation. The density bump and vortex are attributed to enhanced angular momentum, which is promoted by the gravitational torque of the stars. In a hot model ($c \ge 2.0$), the asymmetry rotates considerably more slowly than in the moderate-temperature models. The cold models ($c \le 0.02$) exhibit eccentric circumbinary disks, the precession of which is approximated by a secular motion of the ballistic particles. The asymmetry in the circumbinary disk does not depend on the mass ratio, but it becomes less clear as the specific angular momentum of the infalling envelope increases. The relative accretion rate onto the stars is sensitive to the angular momentum of the infalling envelope. For envelopes with constant angular momentum, the secondary tends to have a higher accretion rate than the primary, except in very low angular momentum cases. For envelopes with a constant angular velocity, the primary has a higher accretion rate than the secondary because gas with low specific angular momentum falls along the polar directions.

The relationship between submillimeter (submm) dust emission polarization and near-infrared (NIR) $H$-band polarization produced by dust dichroic extinction was studied for the cold starless dense core FeSt 1-457. The distribution of polarization angles ($90^{\circ}$-rotated for submm) and degrees were found to be very different between at submm and NIR wavelengths. The mean polarization angles for FeSt 1-457 at submm and NIR wavelengths are $132.1^{\circ} \pm 22.0^{\circ}$ and $2.7^{\circ} \pm 16.2^{\circ}$, respectively. The correlation between $P_H$ and $A_V$ was found to be linear from outermost regions to relatively dense line of sight of $A_V \approx 25$ mag, indicating that NIR polarization reflects overall polarization (magnetic field) structure of the core at least in this density range. The flat $P_H/A_V$ versus $A_V$ correlations were confirmed, and the polarization efficiency was found to be comparable to the observational upper limit (Jones 1989). On the other hand, as reported by Alves et al., submm polarization degrees show clear linearly decreasing trend against $A_V$ from $A_V \approx 20$ mag to the densest center ($A_V \approx 41$ mag), appearing as "polarization hole" structure. The power law index for the $P_{\rm submm}$ versus $A_V$ relationship was obtained to be $\approx -1$, indicating that the alignment for the submm sensitive dust is lost. These very different polarization distributions at submm and NIR wavelengths suggest that (1) there is different radiation environment at these wavelengths or (2) submm-sensitive dust is localized or the combination of them.

Short gamma-ray bursts (sGRBs) are thought to be produced by binary NS mergers. While a sGRB requires a relativistic jet to break out of ejecta, the jet may be choked and fails to produce a successful sGRB. We propose a "delayed breakout" scenario where a late-time jet launched by a long-term engine activity can penetrate ejecta even if a prompt jet is choked. Observationally, such a late-time jet is supported by the long-lasting high-energy emissions in sGRBs. Solving the jet propagation in ejecta, we show that a typical late-time activity easily achieves the delayed breakout. This event shows not prompt $\gamma$-rays but long-time X-ray emissions for $\sim10^{2-3}\,\rm{s}$ or even $\sim10^{4-5}\,\rm{s}$. Some delayed events may be already detected as soft-long GRBs without supernova signatures. In an optimistic case, a few events coincident with gravitational-waves (GWs) are detected by the second-generation GW detectors every year. X-ray followups of merger events without $\gamma$-rays will be a probe of long-lasting engine activities in binary mergers.

The Multipurpose InfraRed Imaging System (MIRIS) performed the MIRIS Pa{\alpha} Galactic Plane Survey (MIPAPS), which covers the entire Galactic plane within the latitude range of -3{\deg} < b < +3{\deg} at Pa{\alpha} (1.87 um). We present the first result of the MIPAPS data extracted from the longitude range of l = 96.5{\deg}-116.3{\deg}, and demonstrate the data quality and scientific potential of the data by comparing them with H{\alpha} maps obtained from the INT Photometric H{\alpha} Survey (IPHAS) data. We newly identify 90 H II region candidates in the WISE H II region catalog as definite H II regions by detecting the Pa{\alpha} and/or H{\alpha} recombination lines, out of which 53 H II regions are detected at Pa{\alpha}. We also report the detection of additional 29 extended and 18 point-like sources at Pa{\alpha}. We estimate the E(B-V) color excesses and the total Lyman continuum luminosities for H II regions by combining the MIPAPS Pa{\alpha} and IPHAS H{\alpha} fluxes. The E(B-V) values are found to be systematically lower than those estimated from point stars associated with H II regions. Utilizing the MIPAPS Pa{\alpha} and IPHAS H{\alpha} images, we obtain an E(B-V) map for the entire region of the H II region Sh2-131 with an angular size of ~2.5{\deg}. The E(B-V) map shows not only numerous high-extinction filamentary features but also negative E(B-V) regions, indicating H{\alpha} excess. The H{\alpha} excess and the systematic underestimation of E(B-V) are attributed to light scattered by dust.

Macronovae (or kilonovae) are the optical and NIR counterparts of binary neutron star mergers. While the macronova in GW170817 was detected about 10 hours after the GW detection, future observations can possibly detect them within the first hour after the merger. Early- time macronovae are potentially powered by some mechanisms such as the beta-decay heating of the surviving free neutrons. In this paper, we propose that the polarimetric observation can be a useful tool to study the early macronova emissions. If free nucleons remain in the outermost layer of the ejecta, the electron scattering produces a larger polarization than that by the r-process element-rich ejecta. The degree of polarization can show a large value of $\sim3\%$ for the first $0.3-1$ hour for the free nucleon mass of $10^{-5}-10^{-4}\,M_{\odot}$. Quick polarimetric observations enable us to study not only the aspherical morphology of the ejecta but also the amount of the free nucleons in the ejecta, which is helpful to discriminate the emission mechanisms of the early macronovae.

We reanalyze the Imaging Photopolarimeter data from Pioneer 10 to study the zodiacal light in the B and R bands beyond Earth orbit, applying an improved method to subtract integrated star light (ISL) and diffuse Galactic light (DGL). We found that there exists a significant instrumental offset, making it difficult to examine the absolute sky brightness. Instead, we analyzed the differential brightness, i.e., the difference in sky brightness from the average at high ecliptic latitude, and compared with that expected from the model zodiacal light. At a heliocentric distance of r<2 au, we found a fairly good correlation between the J-band model zodiacal light and the residual sky brightness after subtracting the ISL and DGL. The reflectances of the interplanetary dust derived from the correlation study are marginally consistent with previous works. The zodiacal light is not significantly detectable at r>3 au, as previously reported. However, a clear discrepancy from the model is found at r=2.94 au which indicates the existence of a local dust cloud produced by the collision of asteroids or dust trail from active asteroids (or main-belt comets). Our result confirms that the main component of the zodiacal light (smooth cloud) is consistent with the model even beyond the earth orbit, which justifies the detection of the extragalactic background light after subtracting the zodiacal light based on the model.

We report ALMA Cycle 3 observations in CO isotopes toward a dense core, MC27/L1521F in Taurus, which is considered to be at an early stage of multiple star formation in a turbulent environment. Although most of the high-density parts of this core are considered to be as cold as $\sim$10 K, high-angular resolution ($\sim$20 au) observations in $^{12}$CO ($J$ = 3--2) revealed complex warm ($>$15--60 K) filamentary/clumpy structures with the sizes from a few tens of au to $\sim$1,000 au. The interferometric observations of $^{13}$CO and C$^{18}$O show that the densest part with arc-like morphologies associated with the previously identified protostar and condensations are slightly redshifted from the systemic velocity of the core. We suggest that the warm CO clouds may be consequences of shock heating induced by interactions among the different density/velocity components that originated from the turbulent motions in the core. However, such a small-scale and fast turbulent motion does not correspond to a simple extension of the line-width-size relation (i.e., Larson'{}s law), and thus the actual origin remains to be studied. The high-angular resolution CO observations are expected to be essential in detecting small-scale turbulent motions in dense cores and to investigate protostar formation therein.

The gravitational waves from the neutron star merger event GW170817 were accompanied by an unusually weak short GRB 170817A, by an optical/IR macronova/kilonova and by a long lasting radio to X-rays counterpart. While association of short GRBs with mergers was predicted a long time ago, the luminosity of this prompt {\gamma}-ray emission was weaker by a few orders of magnitude than all known previous sGRBs and it was softer than typical sGRBs. This raise the question whether the {\gamma}-rays that we have seen were a regular sGRB viewed off-axis. We revisit this question following recent refined analyses of the {\gamma}-ray signal and the VLBI observations that revealed the angular structure of the relativistic outflow: observing angle of $\sim\,20^\circ$, a narrow jet with core $\lesssim5^\circ$ and $E_{\rm iso}>10^{52}\,{\rm ergs}$. We show here that: (i) The region emitting the observed $\gamma$-rays must have been moving with a Lorentz factor $\Gamma\gtrsim5$; (ii) The observed {\gamma}-rays were not "off-axis" emission (viewing angle $>1/\Gamma$) emerging from the core of the jet, where a regular sGRB was most likely produced; (iii) The $\gamma$-ray emission region was either "on-axis" (at an angle $<1/\Gamma$) or if it was "off-axis" then the observing angle must have been small ($<5^\circ$) and the on-axis emission from this region was too faint and too hard to resemble a regular sGRB.

We report near-infrared (IR) observations of high Galactic latitude clouds to investigate diffuse Galactic light (DGL), which is starlight scattered by interstellar dust grains. The observations were performed at $1.1$ and $1.6\,\rm{\mu m}$ with a wide-field camera instrument, the Multi-purpose Infra-Red Imaging System (MIRIS) onboard the Korean satellite STSAT-3. The DGL brightness is measured by correlating the near-IR images with a far-IR $100\,\rm{\mu m}$ map of interstellar dust thermal emission. The wide-field observation of DGL provides the most accurate DGL measurement achieved to date. We also find a linear correlation between optical and near-IR DGL in the MBM32 field. To study interstellar dust properties in MBM32, we adopt recent dust models with or without $\rm{\mu m}$-sized very large grains and predict the DGL spectra, taking into account reddening effect of interstellar radiation field. The result shows that observed color of the near-IR DGL is closer to the model spectra without very large grains. This may imply that dust growth in the observed MBM32 field is not active owing to its low density of interstellar medium.

The gravitational wave event GW170817 from a binary neutron star (NS) merger is accompanied by electromagnetic counterparts, and the optical and near-infared emission is called a macronova (or kilonova). Although the radioactivity of synthesized r-process elements is widely discussed as an energy source, its decisive evidence is not clearly shown yet. We discuss a macronova powered by the central engine activities such as jet activities and X-rays from the matter fallback, and show that the engine model allows much broader parameter spaces, in particular smaller ejecta mass ($\sim10^{-4}-0.01\,\Msun$) than the r-process model. The blue and red macronovae are naturally explained by various combinations of the ejecta such as a cocoon and merger ejecta with the energy sources of jets and X-rays. The required energy injection is very similar to the X-ray excess observed in GRB 130603B with the power-law slope of $\sim-1.3$. The required lanthanoid fraction for the opacity can be also consistent with the Galactic one. Early or late multi-wavelength observations are crucial for revealing the central engine of short gamma-ray bursts and the r-process nucleosynthesis.

Almost all black holes (BHs) and BH candidates in our Galaxy have been discovered as soft X-ray transients, so-called X-ray novae. X-ray novae are usually considered to arise from binary systems. Here we propose that X-ray novae are also caused by isolated single BHs. We calculate the distribution of the accretion rate from interstellar matter to isolated BHs, and find that BHs in molecular clouds satisfy the condition of the hydrogen-ionizaiton disk instability, which results in X-ray novae. The estimated event rate is consistent with the observed one. Possible candidates include IGR J17454-2919, XTE J1908-094, and SAX J1711.6-3808. Near infrared photometric and spectroscopic follow-ups can exclude companion stars for a BH census in our Galaxy.

A full 3-dimensional compressible magnetohydrodynamic (MHD) simulation is conducted to investigate the thermal responses of a coronal loop to the dynamic dissipation processes of MHD waves. When the foot points of the loop are randomly and continuously forced, the MHD waves become excited and propagate upward. Then, a 1-MK temperature corona is produced naturally as the wave energy dissipates. The excited wave packets become non-linear just above the magnetic canopy, and the wave energy cascades into smaller spatial scales. Moreover, collisions between counter-propagating Alfv\'{e}n wave packets increase the heating rate, resulting in impulsive temperature increases. Our model demonstrates that the heating events in the wave-heated loops can be nanoflare-like in the sense that they are spatially localized and temporally intermittent.

We report results of the observations of 15 regions in several molecular lines for a statistical study of massive cluster-forming clumps. We identified 24 clumps based on the C18O (J=1-0) data obtained by the NRO 45 m telescope, and found that 16 of them are associated with young clusters. The clumps associated with clusters have a typical mass, radius, and molecular density of ~1 X 10^3 Mo, ~0.5 pc, ~1 X 10^5 cm^-3, respectively. We categorized the clumps and clusters into four types according to the spatial coincidence of gas and star density, and discussed their evolutions: Clumps without clusters (Type 1), clumps showing good correlations with clusters (Type 2), clumps showing poor correlations with clusters (Type 3), and clusters with no associated clumps (Type 4). Analyses of the velocity structures and the chemical compositions imply that the clump + cluster systems should evolve from Type 1 to Type 4. We found that some of the Type 2 clumps are infalling on the clump-scale to form clusters at the clump center, which should commonly occur in the beginning of cluster formation. Interestingly, all of the identified Type 1 clumps are likely to be older than the Type 2 clumps in terms of chemical compositions, suggesting that they have been gravitationally stable for a long time possibly being supported by the strong magnetic field of > 1 mG.Type 1 clumps younger than the observed Type 2 clumps should be very rare to find because of their short lifetime.

We present a photometric catalog for Spitzer Space Telescope warm mission observations of the North Ecliptic Pole (NEP; centered at $\rm R.A.=18^h00^m00^s$, $\rm Decl.=66^d33^m38^s.552$). The observations are conducted with IRAC in 3.6 $\mu$m and 4.5 $\mu$m bands over an area of 7.04 deg$^2$ reaching 1$\sigma$ depths of 1.29 $\mu$Jy and 0.79 $\mu$Jy in the 3.6 $\mu$m and 4.5 $\mu$m bands respectively. The photometric catalog contains 380,858 sources with 3.6 $\mu$m and 4.5 $\mu$m band photometry over the full-depth NEP mosaic. Point source completeness simulations show that the catalog is 80% complete down to 19.7 AB. The accompanying catalog can be utilized in constraining the physical properties of extra-galactic objects, studying the AGN population, measuring the infrared colors of stellar objects, and studying the extra-galactic infrared background light.

We report ALMA observations in 0.87 mm continuum and $^{12}$CO ($J$ = 3--2) toward a very low-luminosity ($<$0.1 $L_{\odot}$) protostar, which is deeply embedded in one of the densest core MC27/L1521F, in Taurus with an indication of multiple star formation in a highly dynamical environment. The beam size corresponds to $\sim$20 AU, and we have clearly detected blueshifted/redshifted gas in $^{12}$CO associated with the protostar. The spatial/velocity distributions of the gas show there is a rotating disk with a size scale of $\sim$10 AU, a disk mass of $\sim$10$^{-4}$ $M_{\odot}$ and a central stellar mass of $\sim$0.2 $M_{\odot}$. The observed disk seems to be detached from the surrounding dense gas, although it is still embedded at the center of the core whose density is $\sim$10$^{6}$ cm$^{-3}$. The current low-outflow activity and the very low luminosity indicate that the mass accretion rate onto the protostar is extremely low in spite of a very early stage of star formation. We may be witnessing the final stage of the formation of $\sim$0.2 $M_{\odot}$ protostar. However, we cannot explain the observed low luminosity with the standard pre-main-sequence evolutionary track unless we assume cold accretion with an extremely small initial radius of the protostar ($\sim$0.65 $R_\odot$). These facts may challenge our current understanding of the low mass star formation, in particular the mass accretion process onto the protostar and the circumstellar disk.

The protoplanetary disk around HL Tau is so far the youngest candidate of planet formation, and it is still embedded in a protostellar envelope with a size of thousands of au. In this work, we study the gas kinematics in the envelope and its possible influence on the embedded disk. We present our new ALMA cycle 3 observational results of HL Tau in the 13CO (2-1) and C18O (2-1) emission at resolutions of 0.8" (110 au), and we compare the observed velocity pattern with models of different kinds of gas motions. Both the 13CO and C18O emission lines show a central compact component with a size of 2" (280 au), which traces the protoplanetary disk. The disk is clearly resolved and shows a Keplerian motion, from which the protostellar mass of HL Tau is estimated to be 1.8+/-0.3 M$_\odot$, assuming the inclination angle of the disk to be 47 deg from the plane of the sky. The 13CO emission shows two arc structures with sizes of 1000-2000 au and masses of 3E-3 M$_\odot$ connected to the central disk. One is blueshifted and stretches from the northeast to the northwest, and the other is redshifted and stretches from the southwest to the southeast. We find that simple kinematical models of infalling and (counter-)rotating flattened envelopes cannot fully explain the observed velocity patterns in the arc structures. The gas kinematics of the arc structures can be better explained with three-dimensional infalling or outflowing motions. Nevertheless, the observed velocity in the northwestern part of the blueshifted arc structure is ~60-70% higher than the expected free-fall velocity. We discuss two possible origins of the arc structures: (1) infalling flows externally compressed by an expanding shell driven by XZ Tau and (2) outflowing gas clumps caused by gravitational instabilities in the protoplanetary disk around HL Tau.

Recent observations suggest that intensive molecular cloud collision can trigger massive star/cluster formation. The most important physical process caused by the collision is a shock compression. In this paper, the influence of a shock wave on the evolution of a molecular cloud is studied numerically by using isothermal magnetohydrodynamics (MHD) simulations with the effect of self-gravity. Adaptive-mesh-refinement and sink particle techniques are used to follow long-time evolution of the shocked cloud. We find that the shock compression of turbulent inhomogeneous molecular cloud creates massive filaments, which lie perpendicularly to the background magnetic field as we have pointed out in a previous paper. The massive filament shows global collapse along the filament, which feeds a sink particle located at the collapse center. We observe high accretion rate dot{M}_acc > 10^{-4} M_sun/yr that is high enough to allow the formation of even O-type stars. The most massive sink particle achieves M>50 M_sun in a few times 10^5 yr after the onset of the filament collapse.

Fragmentation of filaments into dense cores is thought to be an important step in forming stars. The bar-mode instability of spherically collapsing cores found in previous linear analysis invokes a possibility of re-fragmentation of the cores due to their ellipsoidal (prolate or oblate) deformation. To investigate this possibility, here we perform three-dimensional self-gravitational hydrodynamics simulations that follow all the way from filament fragmentation to subsequent core collapse. We assume the gas is polytropic with index \gamma, which determines the stability of the bar-mode. For the case that the fragmentation of isolated hydrostatic filaments is triggered by the most unstable fragmentation mode, we find the bar mode grows as collapse proceeds if \gamma < 1.1, in agreement with the linear analysis. However, it takes more than ten orders-of-magnitude increase in the central density for the distortion to become non-linear. In addition to this fiducial case, we also study non-fiducial ones such as the fragmentation is triggered by a fragmentation mode with a longer wavelength and it occurs during radial collapse of filaments and find the distortion rapidly grows. In most of astrophysical applications, the effective polytropic index of collapsing gas exceeds 1.1 before ten orders-of-magnitude increase in the central density. Thus, supposing the fiducial case of filament fragmentation, re-fragmentation of dense cores would not be likely and their final mass would be determined when the filaments fragment.

The Extragalactic Background Light (EBL) captures the total integrated emission from stars and galaxies throughout the cosmic history. The amplitude of the near-infrared EBL from space absolute photometry observations has been controversial and depends strongly on the modeling and subtraction of the Zodiacal light foreground. We report the first measurement of the diffuse background spectrum at 0.8-1.7 um from the CIBER experiment. The observations were obtained with an absolute spectrometer over two flights in multiple sky fields to enable the subtraction of Zodiacal light, stars, terrestrial emission, and diffuse Galactic light. After subtracting foregrounds and accounting for systematic errors, we find the nominal EBL brightness, assuming the Kelsall Zodiacal light model, is 42.7+11.9/-10.6 nW/m2/sr at 1.4 um. We also analyzed the data using the Wright Zodiacal light model, which results in a worse statistical fit to the data and an unphysical EBL, falling below the known background light from galaxies at <1.3 um. Using a model-independent analysis based on the minimum EBL brightness, we find an EBL brightness of 28.7+5.1/-3.3 nW/m2/sr at 1.4 um. While the derived EBL amplitude strongly depends on the Zodiacal light model, we find that we cannot fit the spectral data to Zodiacal light, Galactic emission, and EBL from solely integrated galactic light from galaxy counts. The results require a new diffuse component, such as an additional foreground or an excess EBL with a redder spectrum than that of Zodiacal light.

We investigate the formation of circumstellar disks and outflows subsequent to the collapse of molecular cloud cores with the magnetic field and turbulence. Numerical simulations are performed by using an adaptive mesh refinement to follow the evolution up to $\sim 1000$~yr after the formation of a protostar. In the simulations, circumstellar disks are formed around the protostars; those in magnetized models are considerably smaller than those in nonmagnetized models, but their size increases with time. The models with stronger magnetic field tends to produce smaller disks. During evolution in the magnetized models, the mass ratios of a disk to a protostar is approximately constant at $\sim 1-10$\%. The circumstellar disks are aligned according to their angular momentum, and the outflows accelerate along the magnetic field on the $10-100$~au scale; this produces a disk that is misaligned with the outflow. The outflows are classified into two types: a magneto-centrifugal wind and a spiral flow. In the latter, because of the geometry, the axis of rotation is misaligned with the magnetic field. The magnetic field has an internal structure in the cloud cores, which also causes misalignment between the outflows and the magnetic field on the scale of the cloud core. The distribution of the angular momentum vectors in a core also has a non-monotonic internal structure. This should create a time-dependent accretion of angular momenta onto the circumstellar disk. Therefore, the circumstellar disks are expected to change their orientation as well as their sizes in the long-term evolutions.

We report the ALMA Cycle 2 observations of the Class I binary protostellar system L1551 NE in the 0.9-mm continuum, C18O (3-2), 13CO (3-2), SO (7_8-6_7), and the CS (7-6) emission. At 0.18" (= 25 AU) resolution, ~4-times higher than that of our Cycle 0 observations, the circumbinary disk as seen in the 0.9-mm emission is shown to be comprised of a northern and a southern spiral arm, with the southern arm connecting to the circumstellar disk around Source B. The western parts of the spiral arms are brighter than the eastern parts, suggesting the presence of an m=1 spiral mode. In the C18O emission, the infall gas motions in the inter-arm regions and the outward gas motions in the arms are identified. These observed features are well reproduced with our numerical simulations, where gravitational torques from the binary system impart angular momenta to the spiral-arm regions and extract angular momenta from the inter-arm regions. Chemical differentiation of the circumbinary disk is seen in the four molecular species. Our Cycle 2 observations have also resolved the circumstellar disks around the individual protostars, and the beam-deconvolved sizes are 0.29" X 0.19" (= 40 X 26 AU) (P.A. = 144 deg) and 0.26" X 0.20" (= 36 X 27 AU) (P.A. = 147 deg) for Sources A and B, respectively. The position and inclination angles of these circumstellar disks are misaligned with that of the circumbinary disk. The C18O emission traces the Keplerian rotation of the misaligned disk around Source A.

We present near-infrared (0.8-1.8 $\mu$m) spectra of 105 bright (${m_{J}}$ $<$ 10) stars observed with the low resolution spectrometer on the rocket-borne Cosmic Infrared Background Experiment (CIBER). As our observations are performed above the earth's atmosphere, our spectra are free from telluric contamination, which makes them a unique resource for near-infrared spectral calibration. Two-Micron All Sky Survey (2MASS) photometry information is used to identify cross-matched stars after reduction and extraction of the spectra. We identify the spectral types of the observed stars by comparing them with spectral templates from the Infrared Telescope Facility (IRTF) library. All the observed spectra are consistent with late F to M stellar spectral types, and we identify various infrared absorption lines.

The first direct detections of gravitational waves (GWs) from black hole (BH) mergers, GW150914, GW151226 and LVT151012, give a robust lower limit $\sim 70000$ on the number of merged, highly-spinning BHs in our Galaxy. The total spin energy is comparable to all the kinetic energy of supernovae that ever happened in our Galaxy. The BHs release the spin energy to relativistic jets by accreting matter and magnetic fields from the interstellar medium (ISM). By considering the distributions of the ISM density, BH mass and velocity, we calculate the luminosity function of the BH jets, and find that they can potentially accelerate TeV-PeV cosmic-ray particles in our Galaxy with total power $\sim 10^{37\pm 3}$ erg s$^{-1}$ as PeVatrons, positron factories and/or unidentified TeV gamma-ray sources. Additional $\sim 300$ BH jet nebulae could be detectable by CTA (Cherenkov Telescope Array). We also argue that the accretion from the ISM can evaporate and blow away cold material around the BH, which has profound implications for some scenarios to predict electromagnetic counterparts to BH mergers.

Quiescent prominences host a diverse range of flows, including Rayleigh-Taylor instability driven upflows and impulsive downflows, and so it is no surprise that turbulent motions also exist. As prominences are believed to have a mean horizontal guide field, investigating any turbulence they host could shed light on the nature of MHD turbulence in a wide range of astrophysical systems. In this paper we have investigated the nature of the turbulent prominence motions using structure function analysis on the velocity increments estimated from H$\alpha$ Dopplergrams constructed with observational data from Hinode SOT. The pdf of the velocity increments shows that as we look at increasingly small spatial separations the distribution displays greater departure from a reference Gaussian distribution, hinting at intermittency in the velocity field. Analysis of the even order structure functions for both the horizontal and vertical separations showed the existence of two distinct regions displaying different exponents of the power law with the break in the power law at approximately 2000km. We hypothesise this to be a result of internal turbulence excited in the prominence by the dynamic flows of the system found at this spatial scale. We found that the scaling exponents of the p-th order structure functions for these two regions generally followed the p/2 (smaller scales) and p/4 (larger scales) laws that are the same as those predicted for weak MHD turbulence and Kraichnan-Iroshnikov turbulence respectively. However, the existence of the p/4 scaling at larger scales than the p/2 scaling is inconsistent with the increasing nonlinearity expected in MHD turbulence. Estimating the heating from the turbulent energy dissipation showed that this turbulent heating would be very inefficient, but that the mass diffusion through turbulence driven reconnection was of the order of $10^{10}$cm$^2$/s.

Supermassive black holes observed at high redshift $z\gtrsim6$ could grow from direct collapse black holes (DCBHs) with mass $\sim10^5\,M_{\odot}$, which result from the collapse of supermassive stars (SMSs). If a relativistic jet is launched from a DCBH, it can break out of the collapsing SMS and produce a gamma-ray burst (GRB). Although most of the GRB jets are off-axis from our line of sight, we show that the energy injected from the jet into a cocoon is huge $\sim10^{55-56}\,{\rm{erg}}$, so that the cocoon fireball is observed as ultra-luminous supernovae of $\sim10^{45-46}\rm{\,erg\,s^{-1}}$ for $\sim5000 [(1+z)/16] \rm{\,days}$. They are detectable by the future telescopes with near infrared bands, such as, $Euclid$, $WFIRST$, $WISH$, and $JWST$ up to $z\sim20$ and $\sim 100$ events per year, providing a direct evidence of the DCBH scenario.

We report the discovery of infalling motion with rotation of S235AB the massive cluster-forming clump (~10^3 Mo) in the S235 region. Our C18O observations with the 45m telescope at the Nobeyama Radio Observatory have revealed the elliptical shape of the clump. Position-velocity (PV) diagram taken along its major axis exhibits two well-defined peaks symmetrically located with respect to the clump center, which is similar to that found for a dynamically infalling envelope with rotation around a single protostar modeled by N. Ohashi and his collaborators, indicating that the cluster-forming clump is also collapsing by the self-gravity toward the clump center. With analogue to Ohashi's model, we made a simple model of an infalling, rotating clump to fit the observed data. Based on the inferred model parameters as well as results of earlier observations and simulations in the literature, we discuss structures of the clump such as the relation among the global mass infall rate (~10^-3 Mo/yr), formation of a compact core (with a mass and size of ~4 Mo and < 0.1 pc) at the center, and a massive star (~11 Mo) forming in the core.

We present observations at 7 mm that fully resolve the two circumstellar disks, and a reanalyses of archival observations at 3.5 cm that resolve along their major axes the two ionized jets, of the class I binary protostellar system L1551 NE. We show that the two circumstellar disks are better fit by a shallow inner and steep outer power-law than a truncated power-law. The two disks have very different transition radii between their inner and outer regions of $\sim$18.6 AU and $\sim$8.9 AU respectively. Assuming that they are intrinsically circular and geometrically thin, we find that the two circumstellar disks are parallel with each other and orthogonal in projection to their respective ionized jets. Furthermore, the two disks are closely aligned if not parallel with their circumbinary disk. Over an interval of $\sim$10 yr, source B (possessing the circumsecondary disk) has moved northwards with respect to and likely away from source A, indicating an orbital motion in the same direction as the rotational motion of their circumbinary disk. All the aforementioned elements therefore share the same axis for their angular momentum, indicating that L1551 NE is a product of rotationally-driven fragmentation of its parental core. Assuming a circular orbit, the relative disk sizes are compatible with theoretical predictions for tidal truncation by a binary system having a mass ratio of $\sim$0.2, in agreement with the reported relative separations of the two protostars from the center of their circumbinary disk. The transition radii of both disks, however, are a factor of $\gtrsim$1.5 smaller than their predicted tidally-truncated radii.

2.5-dimensional magnetohydrodynamic (MHD) simulations are performed with high spatial resolution in order to distinguish between competing models of the coronal heating problem. A single coronal loop powered by Alfv\'{e}n waves excited in the photosphere is the target of the present study. The coronal structure is reproduced in our simulations as a natural consequence of the transportation and dissipation of Alfv\'{e}n waves. Further, the coronal structure is maintained as the spatial resolution is changed from 25 to 3 km, although the temperature at the loop top increases with the spatial resolution. The heating mechanisms change gradually across the magnetic canopy at a height of 4 Mm. Below the magnetic canopy, both the shock and the MHD turbulence are dominant heating processes. Above the magnetic canopy, the shock heating rate reduces to less than 10 % of the total heating rate while the MHD turbulence provides significant energy to balance the radiative cooling and thermal conduction loss or gain. The importance of compressibility shown in the present study would significantly impact on the prospects of successful MHD turbulence theory in the solar chromosphere.

We present the results of ALMA observations of dust continuum emission and molecular rotational lines toward a dense core MC27 (aka L1521F) in Taurus, which is considered to be at a very early stage of star formation. The detailed column density distribution on size scales from a few tens AU to ~10,000 AU scale are revealed by combining the ALMA (12 m array + 7 m array) data with the published/unpublished single-dish data. The high angular resolution observations at 0.87 mm with a synthesized beam size of ~0."74 x 0."32 reveal that a protostellar source, MMS-1, is not spatially resolved and lacks associated gas emission, while a starless high-density core, MMS-2, has substructures both in dust and molecular emission. The averaged radial column density distribution of the inner part of MC27/L1521F (r $\lesssim$ 3000 AU) is N(H2) ~r$^{-0.4}$, clearly flatter than that of the outer part, ~r$^{-1.0}$. The complex velocity/spatial structure obtained with previous ALMA observations is located inside the inner flatter region, which may reflect the dynamical status of the dense core.

Either bulk rotation or local turbulence is widely invoked to drive fragmentation in collapsing cores so as to produce multiple star systems. Even when the two mechanisms predict different manners in which the stellar spins and orbits are aligned, subsequent internal or external interactions can drive multiple systems towards or away from alignment thus masking their formation process. Here, we demonstrate that the geometrical and dynamical relationship between the binary system and its surrounding bulk envelope provide the crucial distinction between fragmentation models. We find that the circumstellar disks of the binary protostellar system L1551 IRS 5 are closely parallel not just with each other but also with their surrounding flattened envelope. Measurements of the relative proper motion of the binary components spanning nearly 30 yr indicate an orbital motion in the same sense as the envelope rotation. Eliminating orbital solutions whereby the circumstellar disks would be tidally truncated to sizes smaller than are observed, the remaining solutions favor a circular or low-eccentricity orbit tilted by up to $\sim$25$^\circ$ from the circumstellar disks. Turbulence-driven fragmentation can generate local angular momentum to produce a coplanar binary system, but which bears no particular relationship with its surrounding envelope. Instead, the observed properties conform with predictions for rotationally-driven fragmentation. If the fragments were produced at different heights or on opposite sides of the midplane in the flattened central region of a rotating core, the resulting protostars would then exhibit circumstellar disks parallel with the surrounding envelope but tilted from the orbital plane as is observed.

The existence of black holes (BHs) of mass ~ 10^{9} M_sun at z > 6 is a big puzzle in astrophysics because even optimistic estimates of the accretion time are insufficient for stellar mass BHs of ~ 10 M_sun to grow into such supermassive BHs. A resolution of this puzzle might be the direct collapse of supermassive stars with mass M ~ 10^{5} M_sun into massive seed BHs. We find that if a jet is launched from the accretion disk around the central BH, the jet can break out the star because of the structure of the radiation pressure-dominated envelope. Such ultra-long gamma-ray bursts with duration of ~ 10^{4} - 10^{6} s and flux of 10^{-11} - 10^{-8} erg s^{-1} cm^{-2} could be detectable by Swift. We estimate an event rate of < 1yr^{-1}. The total explosion energy is > 10^{55} - 10^{56} erg. The resulting negative feedback delays the growth of the remnant BH by about 70 Myr or evacuates the host galaxy completely.

We reanalyze data of near-infrared background taken by Infrared Telescope in Space (IRTS) based on up-to-date observational results of zodiacal light, integrated star light and diffuse Galactic light. We confirm the existence of residual isotropic emission, which is slightly lower but almost the same as previously reported. At wavelengths longer than 2 {\mu}m, the result is fairly consistent with the recent observation with AKARI. We also perform the same analysis using a different zodiacal light model by Wright and detected residual isotropic emission that is slightly lower than that based on the original Kelsall model. Both models show the residual isotropic emission that is significantly brighter than the integrated light of galaxies.

Solar flares are an explosive phenomenon, where super-sonic flows and shocks are expected in and above the post-flare loops. To understand the dynamics of post-flare loops, a two-dimensional magnetohydrodynamic (2D MHD) simulation of a solar flare has been carried out. We found new shock structures in and above the post-flare loops, which were not resolved in the previous work by Yokoyama and Shibata 2001. To study the dynamics of flows along the reconnected magnetic field, kinematics and energetics of the plasma are investigated along selected field lines. It is found that shocks are crucial to determine the thermal and flow structures in the post-flare loops. On the basis of the 2D MHD simulation, we have developed a new post-flare loop model which we call the pseudo-2D MHD model. The model is based on the 1D MHD equations, where all the variables depend on one space dimension and all the three components of the magnetic and velocity fields are considered. Our pseudo-2D model includes many features of the multi-dimensional MHD processes related to magnetic reconnection (particularly MHD shocks), which the previous 1D hydrodynamic models are not able to include. We compare the shock formation and energetics of a specific field line in the 2D calculation with those in our pseudo-2D MHD model, and we found that they give similar results. This model will allow us to study the evolution of the post-flare loops in a wide parameter space without expensive computational cost and without neglecting important physics associated with magnetic reconnection.

We report spatial fluctuation analysis of the sky brightness in near-infrared from observations toward the north ecliptic pole (NEP) by the AKARI at 2.4 and 3.2 micron. As a follow up study of our previous work on the Monitor field of AKARI, we used NEP deep survey data, which covered a circular area of about 0.4 square degrees, in order to extend fluctuation analysis at angular scales up to 1000". We found residual fluctuation over the estimated shot noise at larger angles than the angular scale of the Monitor field. The excess fluctuation of the NEP deep field smoothly connects with that of the Monitor field at angular scales with a few hundreds arcseconds and extends without any significant variation to larger angular scales up to 1000". By comparing excess fluctuations at two wavelengths, we confirm a blue spectrum feature similar to the result of the Monitor field. We find that the result of this study is consistent with Spitzer Space Telescope observations at 3.6 micron. The origin of the excess fluctuation in the near-infrared background remains to be answered, but we could exclude zodiacal light, diffuse Galactic light, and unresolved faint galaxies at low-redshift based on the comparison with mid- and far-infrared brightness, ground based near-infrared images.

We report measurements of the Diffuse Galactic Light (DGL) spectrum in the near-infrared, spanning the wavelength range 0.95-1.65 {\mu}m by the Cosmic Infrared Background ExpeRiment (CIBER). Using the low-resolution spectrometer (LRS) calibrated for absolute spectro-photometry, we acquired long-slit spectral images of the total diffuse sky brightness towards four high-latitude fields spread over four sounding rocket flights. To separate the DGL spectrum from the total sky brightness, we correlated the spectral images with a 100 {\mu}m intensity map, which traces the dust column density in optically thin regions. The measured DGL spectrum shows no resolved features and is consistent with other DGL measurements in the optical and at near-infrared wavelengths longer than 1.8 {\mu}m. Our result implies that the continuum is consistently reproduced by models of scattered starlight in the Rayleigh scattering regime with a few large grains.

We investigated the formation of arc-like structures in the infalling envelope around protostars, motivated by the recent Atacama Large Millimeter/Submillimeter Array (ALMA) observations of the high-density molecular cloud core, MC27/L1521F. We performed self-gravitational hydrodynamical numerical simulations with an adaptive mesh refinement code. A filamentary cloud with a 0.1~pc width fragments into cloud cores because of perturbations due to weak turbulence. The cloud core undergoes gravitational collapse to form multiple protostars, and gravitational torque from the orbiting protostars produces arc structures extending up to a 1000~AU scale. As well as on a spatial extent, the velocity ranges of the arc structures, $\sim0.5\,\mathrm{km\,s}^{-1}$, are in agreement with the ALMA observations. We also found that circumstellar disks are often misaligned in triple system. The misalignment is caused by the tidal interaction between the protostars when they undergo close encounters because of a highly eccentric orbit of the tight binary pair.

Using self-gravitational hydrodynamical numerical simulations, we investigated the evolution of high-density turbulent molecular clouds swept by a colliding flow. The interaction of shock waves due to turbulence produces networks of thin filamentary clouds with a sub-parsec width. The colliding flow accumulates the filamentary clouds into a sheet cloud and promotes active star formation for initially high-density clouds. Clouds with a colliding flow exhibit a finer filamentary network than clouds without a colliding flow. The probability distribution functions (PDFs) for the density and column density can be fitted by lognormal functions for clouds without colliding flow. When the initial turbulence is weak, the column density PDF has a power-law wing at high column densities. The colliding flow considerably deforms the PDF, such that the PDF exhibits a double peak. The stellar mass distributions reproduced here are consistent with the classical initial mass function with a power-law index of $-1.35$ when the initial clouds have a high density. The distribution of stellar velocities agrees with the gas velocity distribution, which can be fitted by Gaussian functions for clouds without colliding flow. For clouds with colliding flow, the velocity dispersion of gas tends to be larger than the stellar velocity dispersion. The signatures of colliding flows and turbulence appear in channel maps reconstructed from the simulation data. Clouds without colliding flow exhibit a cloud-scale velocity shear due to the turbulence. In contrast, clouds with colliding flow show a prominent anti-correlated distribution of thin filaments between the different velocity channels, suggesting collisions between the filamentary clouds.

Extragalactic background light (EBL) anisotropy traces variations in the total production of photons over cosmic history, and may contain faint, extended components missed in galaxy point source surveys. Infrared EBL fluctuations have been attributed to primordial galaxies and black holes at the epoch of reionization (EOR), or alternately, intra-halo light (IHL) from stars tidally stripped from their parent galaxies at low redshift. We report new EBL anisotropy measurements from a specialized sounding rocket experiment at 1.1 and 1.6 micrometers. The observed fluctuations exceed the amplitude from known galaxy populations, are inconsistent with EOR galaxies and black holes, and are largely explained by IHL emission. The measured fluctuations are associated with an EBL intensity that is comparable to the background from known galaxies measured through number counts, and therefore a substantial contribution to the energy contained in photons in the cosmos.

We report results of molecular line observations carried out toward a massive dense core in the Cyg OB 7 molecular cloud. The core has an extraordinarily large mass ($\sim1.1 \times 10^4$ $M_\odot$) and size ($\sim2 \times 5$ pc$^2$), but there is no massive young star forming therein. We observed this core in various molecular lines such as C$^{18}$O($J=1-0$) using the 45m telescope at Nobeyama Radio Observatory. We find that the core has an elongated morphology consisting of several filaments and core-like structures. The filaments are massive ($10^2-10^3$ $M_\odot$), and they are apparently colliding against each other. Some candidates of YSOs are distributed around their intersection, suggesting that the collisions of the filaments may have influenced on their formation. To understand the formation and evolution of such colliding filaments, we performed numerical simulations using the adaptive mesh refinement (AMR) technique adopting the observed core parameters (e.g., the mass and size) as the initial conditions. Results indicate that the filaments are formed as seen in other earlier simulations for small cores in literature, but we could not reproduce the collisions of the filaments simply by assuming the large initial mass and size. We find that the collisions of the filaments occur only when there is a large velocity gradient in the initial core in a sense to compress it. We suggest that the observed core was actually compressed by an external effect, e.g., shocks of nearby supernova remnants including HB21 which has been suggested to be interacting with the Cyg OB 7 molecular cloud.

We report the ALMA observation of the Class I binary protostellar system L1551 NE in the 0.9-mm continuum, C18O (3-2), and 13CO (3-2) lines at a ~1.6 times higher resolution and a ~6 times higher sensitivity than those of our previous SMA observations, which revealed a r ~300 AU-scale circumbinary disk in Keplerian rotation. The 0.9-mm continuum shows two opposing U-shaped brightenings in the circumbinary disk, and exhibits a depression between the circumbinary disk and the circumstellar disk of the primary protostar. The molecular lines trace non-axisymmetric deviations from Keplerian rotation in the circumbinary disk at higher velocities relative to the systemic velocity, where our previous SMA observations could not detect the lines. In addition, we detect inward motion along the minor axis of the circumbinary disk. To explain the newly-observed features, we performed a numerical simulation of gas orbits in a Roche potential tailored to the inferred properties of L1551 NE. The observed U-shaped dust features coincide with locations where gravitational torques from the central binary system are predicted to impart angular momentum to the circumbinary disk, producing shocks and hence density enhancements seen as a pair of spiral arms. The observed inward gas motion coincides with locations where angular momentum is predicted to be lowered by the gravitational torques. The good agreement between our observation and model indicates that gravitational torques from the binary stars constitute the primary driver for exchanging angular momentum so as to permit infall through the circumbinary disk of L1551 NE.

We are now investigating and studying a small satellite mission HiZ-GUNDAM for future observation of gamma-ray bursts (GRBs). The mission concept is to probe "the end of dark ages and the dawn of formation of astronomical objects", i.e. the physical condition of early universe beyond the redshift z > 7. We will consider two kinds of mission payloads, (1) wide field X-ray imaging detectors for GRB discovery, and (2) a near infrared telescope with 30 cm in diameter to select the high-z GRB candidates effectively. In this paper, we explain some requirements to promote the GRB cosmology based on the past observations, and also introduce the mission concept of HiZ-GUNDAM and basic development of X-ray imaging detectors.

Starless dense cores eventually collapse dynamically, which forms protostars inside them, and the physical properties of the cores determine the nature of the forming protostars. We report ALMA observations of dust continuum emission and molecular rotational lines toward MC27 or L1521F, which is considered to be very close to the first protostellar core phase. We found a few starless high-density cores, one of which has a very high density of ~10^7 cm^-3, within a region of several hundred AU around a very low-luminosity protostar detected by Spitzer. A very compact bipolar outflow with a dynamical timescale of a few hundred years was found toward the protostar. The molecular line observation shows several cores with an arc-like structure, possibly due to the dynamical gas interaction. These complex structures revealed in the present observations suggest that the initial condition of star formation is highly dynamical in nature, which is considered to be a key factor in understanding fundamental issues of star formation such as the formation of multiple stars and the origin of the initial mass function of stars.

We have performed a 2.5 dimensional magnetohydrodynamic simulation that resolves the propagation and dissipation of Alfven waves in the solar atmosphere. Alfvenic fluctuations are introduced on the bottom boundary of the extremely large simulation box that ranges from the photosphere to far above the solar wind acceleration region. Our model is ab initio in the sense that no corona and no wind are assumed initially.The numerical experiment reveals the quasi-steady solution that has the transition from the cool to the hot atmosphere and the emergence of the high speed wind. The global structure of the resulting hot wind solution fairly well agree with the coronal and the solar wind structure inferred from observations. The purpose of this study is to complement the previous paper by Matsumoto & Suzuki (2012) and describe the more detailed results and the analysis method. These results include the dynamics of the transition region and the more precisely measured heating rate in the atmosphere. Particularly, the spatial distribution of the heating rate helps us to interpret the actual heating mechanisms in the numerical simulation.Our estimation method of heating rate turned out to be a good measure for dissipation of Alfven waves and low beta fast waves.

We first obtained the spectrum of the diffuse Galactic light (DGL) at general interstellar space in 1.8-5.3 um wavelength region with the low-resolution prism spectroscopy mode of the AKARI Infra-Red Camera (IRC) NIR channel. The 3.3 um PAH band is detected in the DGL spectrum at Galactic latitude |b| < 15 deg, and its correlations with the Galactic dust and gas are confirmed. The correlation between the 3.3 um PAH band and the thermal emission from the Galactic dust is expressed not by a simple linear correlation but by a relation with extinction. Using this correlation, the spectral shape of DGL at optically thin region (5 deg < |b| < 15 deg) was derived as a template spectrum. Assuming that the spectral shape of this template spectrum is uniform at any position, DGL spectrum can be estimated by scaling this template spectrum using the correlation between the 3.3 um PAH band and the thermal emission from the Galactic dust.

The Extragalactic Background Light (EBL) as an integrated light from outside of our Galaxy includes information of the early universe and the Dark Ages. We analyzed the spectral data of the astrophysical diffuse emission obtained with the low-resolution spectroscopy mode on the AKARI Infra-Red Camera (IRC) in 1.8-5.3 um wavelength region. Although the previous EBL observation in this wavelength region is restricted to the observations by DIRBE and IRTS, this study adds a new independent result with negligible contamination of Galactic stars owing to higher sensitivity for point sources. Other two major foreground components, the zodiacal light (ZL) and the diffuse Galactic light (DGL), were subtracted by taking correlations with ZL brightness estimated by the DIRBE ZL model and with the 100 um dust thermal emission, respectively. The isotropic emission was obtained as EBL, which shows significant excess over integrated light of galaxies at <4 um. The obtained EBL is consistent with the previous measurements by IRTS and DIRBE.

The formation of a circumstellar disk in collapsing cloud cores is investigated with three-dimensional magnetohydrodynamic simulations. We prepare four types of initial cloud having different density profiles and calculate their evolution with or without a sink. To investigate the effect of magnetic dissipation on disk formation, the Ohmic dissipation is considered in some models. Calculations show that disk formation is very sensitive to both the initial cloud configuration and the sink treatment. The disk size considerably differs in clouds with different density profiles even when the initial clouds have almost the same mass-to-flux ratio. Only a very small disk (\sim 10 AU in size) appears in clouds with a uniform density profile, whereas a large disk (\sim 100 AU in size) forms in clouds with a Bonnor-Ebert density profile. In addition, a large sink accretion radius numerically impedes disk formation during the main accretion phase and tends to foster the misleading notion that disk formation is completely suppressed by magnetic braking. The protostellar outflow is also greatly affected by the sink properties. A sink accretion radius of \lesssim 1 AU and sink threshold density of \gtrsim 10^ 13 cm^-3 are necessary for investigating disk formation during the main accretion phase.

We have developed a near-infrared spectrometer designed to measure the absolute intensity of the Solar 854.2 nm CaII Fraunhofer line, scattered by interplanetary dust, in the Zodiacal light spectrum. Based on the known equivalent line width in the Solar spectrum, this measurement can derive the Zodiacal brightness, testing models of the Zodiacal light based on morphology that are used to determine the extragalactic background light in absolute photometry measurements. The spectrometer is based on a simple high-resolution tipped filter placed in front of a compact camera with wide-field refractive optics to provide the large optical throughput and high sensitivity required for rocket-borne observations. We discuss the instrument requirements for an accurate measurement of the absolute Zodiacal light brightness, the measured laboratory characterization, and the instrument performance in flight.

We present the near- and mid-infrared zodiacal light spectrum obtained with the AKARI Infra-Red Camera (IRC). A catalog of 278 spectra of the diffuse sky covering a wide range of Galactic and ecliptic latitudes was constructed. The wavelength range of this catalog is 1.8-5.3 {\mu}m with wavelength resolution of \lambda /\Delta \lambda ~20. Advanced reduction methods specialized for the slit spectroscopy of diffuse sky spectra are developed for constructing the spectral catalog. Based on the comparison analysis of the spectra collected in different seasons and ecliptic latitudes, we confirmed that the spectral shape of the scattered component and the thermal emission component of the zodiacal light in our wavelength range does not show any dependence on location and time, but relative brightness between them varies with location. We also confirmed that the color temperature of the zodiacal emission at 3-5 {\mu}m is 300+/-10 K at any ecliptic latitude. This emission is expected to be originated from sub-micron dust particles in the interplanetary space.

We investigate mass losses via stellar winds from sun-like main sequence stars with a wide range of activity levels. We perform forward-type magnetohydrodynamical numerical experiments for Alfven wave-driven stellar winds with a wide range of the input Poynting flux from the photosphere. Increasing the magnetic field strength and the turbulent velocity at the stellar photosphere from the current solar level, the mass loss rate rapidly increases at first owing to the suppression of the reflection of the Alfven waves. The surface materials are lifted up by the magnetic pressure associated with the Alfven waves, and the cool dense chromosphere is intermittently extended to 10 -- 20 % of the stellar radius. The dense atmospheres enhance the radiative losses and eventually most of the input Poynting energy from the stellar surface escapes by the radiation. As a result, there is no more sufficient energy remained for the kinetic energy of the wind; the stellar wind saturates in very active stars, as observed in Wood et al. The saturation level is positively correlated with B_{r,0}f_0, where B_{r,0} and f_0 are the magnetic field strength and the filling factor of open flux tubes at the photosphere. If B_{r,0}f_0 is relatively large >~ 5 G, the mass loss rate could be as high as 1000 times. If such a strong mass loss lasts for ~ 1 billion years, the stellar mass itself is affected, which could be a solution to the faint young sun paradox. We derive a Reimers-type scaling relation that estimates the mass loss rate from the energetics consideration of our simulations. Finally, we derive the evolution of the mass loss rates, \dot{M} t^{-1.23}, of our simulations, combining with an observed time evolution of X-ray flux from sun-like stars, which is shallower than \dot{M} t^{-2.33+/-0.55} in Wood et al.(2005).

We report the first three-dimensional radiation magnetohydrodynamic (RMHD) simulations of protostellar collapse with and without Ohmic dissipation.We take into account many physical processes required to study star formation processes, including a realistic equation of state. We follow the evolution from molecular cloud cores until protostellar cores are formed with sufficiently high resolutions without introducing a sink particle. The physical processes involved in the simulations and adopted numerical methods are described in detail.We can calculate only about one year after the formation of the protostellar cores with our direct three-dimensional RMHD simulations because of the extremely short timescale in the deep interior of the formed protostellar cores, but successfully describe the early phase of star formation processes. The thermal evolution and the structure of the first and second (protostellar) cores are consistent with previous one-dimensional simulations using full radiation transfer, but differ considerably from preceding multi-dimensional studies with the barotropic approximation. The protostellar cores evolve virtually spherically symmetric in the ideal MHD models because of efficient angular momentum transport by magnetic fields, but Ohmic dissipation enables the formation of the circumstellar disks in the vicinity of the protostellar cores as in previous MHD studies with the barotropic approximation. The formed disks are still small (less than 0.35 AU) because we simulate only the earliest evolution. We also confirm that two different types of outflows are naturally launched by magnetic fields from the first cores and protostellar cores in the resistive MHD models.

We investigate molecular evolution from a molecular cloud core to a first hydrostatic core in three spatial dimensions. We perform a radiation hydrodynamic simulation in order to trace fluid parcels, in which molecular evolution is investigated, using a gas-phase and grain-surface chemical reaction network. We derive spatial distributions of molecular abundances and column densities in the core harboring the first core. We find that the total of gas and ice abundances of many species in a cold era (10 K) remain unaltered until the temperature reaches ~500 K. The gas abundances in the warm envelope and the outer layer of the first core (T < 500 K) are mainly determined via the sublimation of ice-mantle species. Above 500 K, the abundant molecules, such as H2CO, start to be destroyed, and simple molecules, such as CO, H2O and N2 are reformed. On the other hand, some molecules are effectively formed at high temperature; carbon-chains, such as C2H2 and cyanopolyynes, are formed at the temperature of >700 K. We also find that large organic molecules, such as CH3OH and HCOOCH3, are associated with the first core (r < 10 AU). Although the abundances of these molecules in the first core stage are comparable or less than in the protostellar stage (hot corino), reflecting the lower luminosity of the central object, their column densities in our model are comparable to the observed values toward the prototypical hot corino, IRAS 16293-2422. We propose that these large organic molecules can be good tracers of the first cores.

We present the smoothness of the mid-infrared sky from observations by the Japanese infrared astronomical satellite AKARI. AKARI monitored the north ecliptic pole (NEP) during its cold phase with nine wavebands covering from 2.4 to 24 micron, out of which six mid-infrared bands were used in this study. We applied power spectrum analysis to the images in order to search for the fluctuation of the sky brightness. Observed fluctuation is explained by fluctuation of photon noise, shot noise of faint sources, and Galactic cirrus. The fluctuations at a few arcminutes scales at short mid-infrared wavelengths (7, 9, and 11 micron) are largely caused by the diffuse Galactic light of the interstellar dust cirrus. At long mid-infrared wavelengths (15, 18, and 24 micron), photon noise is the dominant source of fluctuation over the scale from arcseconds to a few arcminutes. The residual fluctuation amplitude at 200 arcseconds after removing these contributions is at most 1.04 +/- 0.23 nW m^-2 sr^-1 or 0.05% of the brightness at 24 micron and at least 0.47 +/- 0.14 nW m^-2 sr^-1 or 0.02% at 18 micron. We conclude that the upper limit of the fluctuation in the zodiacal light towards the NEP is 0.03% of the sky brightness, taking 2{\sigma} error into account.

We have developed and characterized an imaging instrument to measure the spatial properties of the diffuse near-infrared extragalactic background light in a search for fluctuations from z > 6 galaxies during the epoch of reionization. The instrument is part of the Cosmic Infrared Background Experiment (CIBER), designed to observe the extragalactic background light above the Earth's atmosphere during a suborbital sounding rocket flight. The imaging instrument incorporates a 2x2 degree field of view, to measure fluctuations over the predicted peak of the spatial power spectrum at 10 arcminutes, and 7"x7" pixels, to remove lower redshift galaxies to a depth sufficient to reduce the low-redshift galaxy clustering foreground below instrumental sensitivity. The imaging instrument employs two cameras with \Delta \lambda / \lambda ~0.5 bandpasses centered at 1.1 and 1.6 microns to spectrally discriminate reionization extragalactic background fluctuations from local foreground fluctuations. CIBER operates at wavelengths where the electromagnetic spectrum of the reionization extragalactic background is thought to peak, and complements fluctuations measurements by AKARI and Spitzer at longer wavelengths. We have characterized the instrument in the laboratory, including measurements of the sensitivity, flat-field response, stray light performance, and noise properties. Several modifications were made to the instrument following a first flight in 2009 February. The instrument performed to specifications in subsequent flights in 2010 July and 2012 March, and the scientific data are now being analyzed.

Solar spicules are the fundamental magnetic structures in the chromosphere and considered to play a key role in channelling the chromosphere and corona. Recently, it was suggested by De Pontieu et al. that there were two types of spicules with very different dynamic properties, which were detected by space- time plot technique in the Ca ii H line (3968 A) wavelength from Hinode/SOT observations. 'Type I' spicule, with a 3-7 minute lifetime, undergoes a cycle of upward and downward motion; in contrast, 'Type II' spicule fades away within dozens of seconds, without descending phase. We are motivated by the fact that for a spicule with complicated 3D motion, the space-time plot, which is made through a slit on a fixed position, could not match the spicule behavior all the time and might lose its real life story. By revisiting the same data sets, we identify and trace 105 and 102 spicules in quiet sun (QS) and coronal hole (CH), respectively, and obtain their statistical dynamic properties. First, we have not found a single convincing example of 'Type II' spicules. Secondly, more than 60% of the identified spicules in each region show a complete cycle, i.e., majority spicules are 'Type I'. Thirdly, the lifetime of spicules in QS and CH are 148 s and 112 s, respectively, but there is no fundamental lifetime difference between the spicules in QS and CH reported earlier. Therefore, the suggestion of coronal heating by 'Type II' spicules should be taken with cautions. Subject headings: Sun: chromosphere Sun:transition region Sun:corona

The solar wind emanates from the hot and tenuous solar corona. Earlier studies using 1.5 dimensional simulations show that Alfv\'{e}n waves generated in the photosphere play an important role in coronal heating through the process of non-linear mode conversion. In order to understand the physics of coronal heating and solar wind acceleration together, it is important to consider the regions from photosphere to interplanetary space as a single system. We performed 2.5 dimensional, self-consistent magnetohydrodynamic simulations, covering from the photosphere to the interplanetary space for the first time. We carefully set up the grid points with spherical coordinate to treat the Alfv\'{e}n waves in the atmosphere with huge density contrast, and successfully simulate the solar wind streaming out from the hot solar corona as a result of the surface convective motion. The footpoint motion excites Alfv\'{e}n waves along an open magnetic flux tube, and these waves traveling upwards in the non-uniform medium undergo wave reflection, nonlinear mode conversion from Alfv\'{e}n mode to slow mode, and turbulent cascade. These processes leads to the dissipation of Alfv\'{e}n waves and acceleration of the solar wind. It is found that the shock heating by the dissipation of the slow mode wave plays a fundamental role in the coronal heating process whereas the turbulent cascade and shock heating drive the solar wind.

Absolute spectrophotometric measurements of diffuse radiation at 1 \mu m to 2 \mu m are crucial to our understanding of the radiative content of the Universe from nucleosynthesis since the epoch of reionization, the composition and structure of the Zodiacal dust cloud in our solar system, and the diffuse galactic light arising from starlight scattered by interstellar dust. The Low Resolution Spectrometer (LRS) on the rocket-borne Cosmic Infrared Background Experiment (CIBER) is a \lambda / \Delta \lambda \sim 15-30 absolute spectrophotometer designed to make precision measurements of the absolute near-infrared sky brightness between 0.75 \mu m < \lambda < 2.1 \mu m. This paper presents the optical, mechanical and electronic design of the LRS, as well as the ground testing, characterization and calibration measurements undertaken before flight to verify its performance. The LRS is shown to work to specifications, achieving the necessary optical and sensitivity performance. We describe our understanding and control of sources of systematic error for absolute photometry of the near-infrared extragalactic background light.

The Cosmic Infrared Background Experiment (CIBER) is a suite of four instruments designed to study the near infrared (IR) background light from above the Earth's atmosphere. The instrument package comprises two imaging telescopes designed to characterize spatial anisotropy in the extragalactic IR background caused by cosmological structure during the epoch of reionization, a low resolution spectrometer to measure the absolute spectrum of the extragalactic IR background, and a narrow band spectrometer optimized to measure the absolute brightness of the Zodiacal light foreground. In this paper we describe the design and characterization of the CIBER payload. The detailed mechanical, cryogenic, and electrical design of the system are presented, including all system components common to the four instruments. We present the methods and equipment used to characterize the instruments before and after flight, and give a detailed description of CIBER's flight profile and configurations. CIBER is designed to be recoverable and has flown twice, with modifications to the payload having been informed by analysis of the first flight data. All four instruments performed to specifications during the second flight, and the scientific data from this flight are currently being analyzed.

We report a search for fluctuations of the sky brightness toward the north ecliptic pole (NEP) with the Japanese infrared astronomical satellite AKARI,at 2.4, 3.2, and 4.1 \mum. We obtained circular maps with 100 diameter field of view which clearly show a spatial structure on scale of a few hundred arcsec. A power spectrum analysis shows that there is a significant excess fluctuation at angular scales larger than 10000 that can't be explained by zodiacal light, diffuse Galactic light, shot noise of faint galaxies or clustering of low redshift galaxies. These results are consistent with observations at 3.6 and 4.5 \mum by NASA's Spitzer Space Telescope. The fluctuating component observed at large angular scales has a blue stellar spectrum which is similar to that of the spectrum of the excess isotropic emission observed with IRTS. A significant spatial correlation between wavelength bands was found, and the slopes of the linear correlations is consistent with the spectrum of the excess fluctuation. These findings indicate that the detected fluctuation could be attributed to the first stars of the universe, i.e. pop. III stars. The observed fluctuation provides an important constraints on the era of the first stars.

Star formation efficiency controlled by the protostellar outflow in a single cloud core is investigated by three-dimensional resistive MHD simulations. Starting from the prestellar cloud core, the star formation process is calculated until the end of the main accretion phase. In the calculations, the mass of the prestellar cloud is parameterized. During the star formation, the protostellar outflow is driven by the circumstellar disk. The outflow extends also in the transverse direction until its width becomes comparable to the initial cloud scale, and thus, the outflow has a wide opening angle of >40 degrees. As a result, the protostellar outflow sweeps up a large fraction of the infalling material and ejects it into the interstellar space. The outflow can eject at most over half of the host cloud mass, significantly decreasing star formation efficiency. The outflow power is stronger in clouds with a greater initial mass. Thus, the protostellar outflow effectively suppresses star formation efficiency in a massive cloud. The outflow weakens significantly and disappears in several free-fall timescales of the initial cloud after the cloud begins to collapse. The natal prestellar core influences the lifetime and size of the outflow. At the end of the main accretion phase, a massive circumstellar disk comparable in mass to the protostar remains. Calculations show that typically, ~30% of the initial cloud mass is converted into the protostar and ~20% remains in the circumstellar disk, while ~40% is ejected into the interstellar space by the protostellar outflow. Therefore, a single cloud core typically has a star formation efficiency of 30-50%.

Using resistive magnetohydrodynamics simulation, we investigate circumstellar disk formation in a strongly magnetized cloud. As the initial state, an isolated cloud core embedded in a low-density interstellar medium with a uniform magnetic field is adopted. The cloud evolution is calculated until almost all gas inside the initial cloud falls onto either the circumstellar disk or a protostar, and a part of the gas is ejected into the interstellar medium by the protostellar outflow driven by the circumstellar disk. In the early main accretion phase, the disk size is limited to \sim 10 AU because the angular momentum of the circumstellar disk is effectively transferred by both magnetic braking and the protostellar outflow. In the later main accretion phase, however, the circumstellar disk grows rapidly and exceeds 100 AU by the end of the main accretion phase. This rapid growth of the circumstellar disk is caused by the depletion of the infalling envelope, while magnetic braking is effective when the infalling envelope is more massive than the circumstellar disk. The infalling envelope cannot brake the circumstellar disk when the latter is more massive than the former. In addition, the protostellar outflow weakens and disappears in the later main accretion phase, because the outflow is powered by gas accretion onto the circumstellar disk. Although the circumstellar disk formed in a magnetized cloud is considerably smaller than that in an unmagnetized cloud, a circumstellar disk exceeding 100 AU can form even in a strongly magnetized cloud.

An implicit method for the ohmic dissipation is proposed. The proposed method is based on the Crank-Nicolson method and exhibits second-order accuracy in time and space. The proposed method has been implemented in the SFUMATO adaptive mesh refinement (AMR) code. The multigrid method on the grids of the AMR hierarchy converges the solution. The convergence is fast but depends on the time step, resolution, and resistivity. Test problems demonstrated that decent solutions are obtained even at the interface between fine and coarse grids. Moreover, the solution obtained by the proposed method shows good agreement with that obtained by the explicit method, which required many time steps. The present method reduces the number of time steps, and hence the computational costs, as compared with the explicit method.

A first adiabatic core is a transient object formed in the early phase of star formation. The observation of a first core is believed to be difficult because of its short lifetime and low luminosity. On the basis of radiation hydrodynamic simulations, we propose a novel theoretical model of first cores, Exposed Long-lifetime First core (ELF). In the very low-mass molecular core, the first core evolves slowly and lives longer than 10,000 years because the accretion rate is considerably low. The evolution of ELFs is different from that of ordinary first cores because radiation cooling has a significant effect there. We also carry out radiation transfer calculation of dust-continuum emission from ELFs to predict their observational properties. ELFs have slightly fainter but similar SEDs to ordinary first cores in radio wavelengths, therefore they can be observed. Although the probabilities that such low mass cores become gravitationally unstable and start to collapse are low, we still can expect that a considerable number of ELFs can be formed because there are many low-mass molecular cloud cores in star-forming regions that can be progenitors of ELFs.

The formation and evolution of a circumstellar disk in magnetized cloud cores is investigated from prestellar core stage until sim 10^4 yr after protostar formation. In the circumstellar disk, fragmentation first occurs due to gravitational instability in a magnetically inactive region, and substellar-mass objects appear. The substellar-mass objects lose their orbital angular momenta by gravitational interaction with the massive circumstellar disk and finally fall onto the protostar. After this fall, the circumstellar disk increases its mass by mass accretion and again induces fragmentation. The formation and falling of substellar-mass objects are repeated in the circumstellar disk until the end of the main accretion phase. In this process, the mass of fragments remain small, because the circumstellar disk loses its mass by fragmentation and subsequent falling of fragments before it becomes very massive. In addition, when fragments orbit near the protostar, they disturb the inner disk region and promote mass accretion onto the protostar. The orbital motion of substellar-mass objects clearly synchronizes with the time variation of the accretion luminosity of the protostar. Moreover, as the objects fall, the protostar shows a strong brightening for a short duration. The intermittent protostellar outflows are also driven by the circumstellar disk whose magnetic field lines are highly tangled owing to the orbital motion of fragments. The time-variable protostellar luminosity and intermittent outflows may be a clue for detecting planetary-mass objects in the circumstellar disk.

Ultraviolet emission from the first generation of stars in the Universe ionized the intergalactic medium in a process which was completed by z~6; the wavelength of these photons has been redshifted by (1+z) into the near infrared today and can be measured using instruments situated above the Earth's atmosphere. First flying in February 2009, the Cosmic Infrared Background Experiment (CIBER) comprises four instruments housed in a single reusable sounding rocket borne payload. CIBER will measure spatial anisotropies in the extragalactic IR background caused by cosmological structure from the epoch of reionization using two broadband imaging instruments, make a detailed characterization of the spectral shape of the IR background using a low resolution spectrometer, and measure the absolute brightness of the Zodical light foreground with a high resolution spectrometer in each of our six science fields. This paper presents the scientific motivation for CIBER and details of its first two flights, including a review of the published scientific results from the first flight and an outlook for future reionization science with CIBER data.

We investigate protostellar collapse of molecular cloud cores by numerical simulations, taking into account turbulence and magnetic fields. By using the adaptive mesh refinement technique, the collapse is followed over a wide dynamic range from the scale of a turbulent cloud core to that of the first core. The cloud core is lumpy in the low density region owing to the turbulence, while it has a smooth density distribution in the dense region produced by the collapse. The shape of the dense region depends mainly on the mass of the cloud core; a massive cloud core tends to be prolate while a less massive cloud core tends to be oblate. In both cases, anisotropy of the dense region increases during the isothermal collapse. The minor axis of the dense region is always oriented parallel to the local magnetic field. All the models eventually yield spherical first cores supported mainly by the thermal pressure. Most of turbulent cloud cores exhibit protostellar outflows around the first cores. These outflows are classified into two types, bipolar and spiral flows, according to the morphology of the associated magnetic field. Bipolar flow often appears in the less massive cloud core. The rotation axis of the first core is oriented parallel to the local magnetic field for bipolar flow, while the orientation of the rotation axis from the global magnetic field depends on the magnetic field strength. In spiral flow, the rotation axis is not aligned with the local magnetic field.

The formation and evolution of the circumstellar disk in the collapsing molecular cloud is investigated from the prestellar stage resolving both the molecular cloud core and the protostar itself. In the collapsing cloud, the first adiabatic core appears prior to the protostar formation. Reflecting the thermodynamics of the collapsing gas, the first core is much more massive than the protostar. When the molecular cloud has no angular momentum, the first core falls onto the protostar and disappears a few years after the protostar formation. On the other hand, when the molecular cloud has an angular momentum, the first core does not disappear even after the protostar formation, and directly evolves into the circumstellar disk with a Keplerian rotation. There are two paths for the formation of the circumstellar disk. When the initial cloud has a considerably small rotational energy, two nested disks appear just after the protostar formation. During the early main accretion phase, the inner disk increases its size and merges with the outer disk (i.e. first core) to form a single circumstellar disk with a Keplerian rotation. On the other hand, when the molecular cloud has a rotational energy comparable to observations, a single Keplerian disk that corresponds to the first core already exists prior to the protostar formation. In such a cloud, the first core density gradually increases, maintaining the Keplerian rotation and forms the protostar inside it. Thus, the protostar is born in the Keplerian disk. In other words, a massive disk already exists before the protostar formation. In each case, the protostar at its formation is already surrounded by a massive circumstellar disk. The circumstellar disk is about 10-100 times more massive than the protostar in the main accretion disk. Such disks are favourable sites for the formation of binary companions and gas-giant planets.

Interplanetary dust (IPD) scatters solar radiation which results in the zodiacal light that dominates the celestial diffuse brightness at optical and near-infrared wavelengths. Both asteroid collisions and cometary ejections produce the IPD, but the relative contribution from these two sources is still unknown. The Low Resolution Spectrometer (LRS) onboard the Cosmic Infrared Background Experiment (CIBER) observed the astrophysical sky spectrum between 750 and 2100 nm over a wide range of ecliptic latitude. The resulting zodiacal light spectrum is redder than the solar spectrum, and shows a broad absorption feature, previously unreported, at approximately 900 nm, suggesting the existence of silicates in the IPD material. The spectral shape of the zodiacal light is isotropic at all ecliptic latitudes within the measurement error. The zodiacal light spectrum, including the extended wavelength range to 2500 nm using IRTS data, is qualitatively similar to the reflectance of S-type asteroids. This result can be explained by the proximity of S-type asteroidal dust to Earth's orbit, and the relativily high albedo of asteridal dust compared with cometary dust.

Studies of the structure and evolution of protoplanetary disks are important for understanding star and planet formation. Here, we present the direct image of an interacting binary protoplanetary system. Both circumprimary and circumsecondary disks are resolved in the near-infrared. There is a bridge of infrared emission connecting the two disks and a long spiral arm extending from the circumprimary disk. Numerical simulations show that the bridge corresponds to gas flow and a shock wave caused by the collision of gas rotating around the primary and secondary stars. Fresh material streams along the spiral arm, consistent with the theoretical scenarios where gas is replenished from a circummultiple reservoir.

We use resistive magnetohydrodynamical simulations with the nested grid technique to study the formation of protoplanetary disks around protostars from molecular cloud cores that provide the realistic environments for planet formation. We find that gaseous planetary-mass objects are formed in the early evolutionary phase by gravitational instability in regions that are decoupled from the magnetic field and surrounded by the injection points of the magnetohydrodynamical outflows during the formation phase of protoplanetary disks. Magnetic decoupling enables massive disks to form and these are subject to gravitational instability, even at ~ 10 AU. The frequent formation of planetary-mass objects in the disk suggests the possibility of constructing a hybrid planet formation scenario, where the rocky planets form later under the influence of the giant planets in the protoplanetary disk.

A cooperative observation with Hida observatory and Hinode satellite was performed on an emerging flux region. The successive Ca II K spectro-heliograms of the emerging flux region were taken by the Domeless Solar Telescope of Hida observatory. Hinode observed the emerging flux region with Ca II H and Fe I Stokes IQUV filtergrams. In this study, detailed dynamics and temporal evolution of the magnetic flux emergence was studied observationally. The event was first detected in the photospheric magnetic field signals. 3 minutes later, the horizontal expansion of the dark area was detected. And then, 7 minutes later than the horizontal expansion, the emerging loops were detected with the maximal rise speed of 2.1 km/s at chromospheric heights. The observed dynamics of emerging magnetic flux from the photosphere to the upper chromosphere is well consistent with the results of previous simulation works. The gradual rising phase of flux tubes with a weak magnetic strength was confirmed by our observation.

We have derived the temporal power spectra of the horizontal velocity of the solar photosphere. The data sets for 14 quiet regions observed with the Gband filter of Hinode/SOT are analyzed to measure the temporal fluctuation of the horizontal velocity by using the local correlation tracking (LCT) method. Among the high resolution (~0.2") and seeing-free data sets of Hinode/SOT, we selected the observations whose duration is longer than 70 minutes and cadence is about 30 s. The so-called k-{\omega} diagrams of the photospheric horizontal velocity are derived for the first time to investigate the temporal evolution of convection. The power spectra derived from k-omega diagrams typically have a double power law shape bent over at a frequency of 4.7 mHz. The power law index in the high frequency range is -2.4 while the power law index in the low frequency range is -0.6. The root mean square of the horizontal speed is about 1.1 km/s when we use a tracer size of 0.4" in LCT method. Autocorrelation functions of intensity fluctuation, horizontal velocity, and its spatial derivatives are also derived in order to measure the correlation time of the stochastic photospheric motion. Since one of possible energy sources of the coronal heating is the photospheric convection, the power spectra derived in the present study will be of high value to quantitatively justify various coronal heating models.

We investigate gas accretion onto a protoplanet, by considering the thermal effect of gas in three-dimensional hydrodynamical simulations, in which the wide region from a protoplanetary gas disk to a Jovian radius planet is resolved using the nested-grid method. We estimate the mass accretion rate and growth timescale of gas giant planets. The mass accretion rate increases with protoplanet mass for M_p < M_cri, while it becomes saturated or decreases for M_p > M_cri, where M_cri = 0.036 M_Jup (a_p/1AU)^0.75, and M_Jup and a_p are the Jovian mass and the orbital radius, respectively. The growth timescale of a gas giant planet or the timescale of the gas accretion onto the protoplanet is about 10^5 yr, that is two orders of magnitude shorter than the growth timescale of the solid core. The thermal effects barely affect the mass accretion rate because the gravitational energy dominates the thermal energy around the protoplanet. The mass accretion rate obtained in our local simulations agrees quantitatively well with those obtained in global simulations with coarser spatial resolution. The mass accretion rate is mainly determined by the protoplanet mass and the property of the protoplanetary disk. We find that the mass accretion rate is correctly calculated when the Hill or Bondi radius is sufficiently resolved. Using the oligarchic growth of protoplanets, we discuss the formation timescale of gas giant planets.

We studied spicular jets over a plage area and derived their dynamic characteristics using Hinode Solar Optical Telescope (SOT) high-resolution images. The target plage region was near the west limb of the solar disk. This location permitted us to study the dynamics of spicular jets without the overlapping effect of spicular structures along the line of sight. In this work, to increase the ease with which we can identify spicules on the disk, we applied the image processing method `MadMax' developed by Koutchmy et al. (1989). It enhances fine, slender structures (like jets), over a diffuse background. We identified 169 spicules over the target plage. This sample permits us to derive statistically reliable results regarding spicular dynamics. The properties of plage spicules can be summarized as follows: (1) In a plage area, we clearly identified spicular jet features. (2) They were shorter in length than the quiet region limb spicules, and followed ballistic motion under constant deceleration. (3) The majority (80%) of the plage spicules showed the cycle of rise and retreat, while 10% of them faded out without a complete retreat phase. (4) The deceleration of the spicule was proportional to the velocity of ejection (i.e. the initial velocity).

We perform a three-dimensional nested-grid radiation magneto-hydrodynamics (RMHD) simulation with self-gravity to study the early phase of the low-mass star formation process from a rotating molecular cloud core to a first adiabatic core just before the second collapse begins. Radiation transfer is handled with the flux-limited diffusion approximation, operator-splitting and implicit time-integrator. In the RMHD simulation, the outer region of the first core attains a higher entropy and the size of first core is larger than that in the magnetohydrodynamics simulations with the barotropic approximation. Bipolar molecular outflow consisting of two components is driven by magnetic Lorentz force via different mechanisms, and shock heating by the outflow is observed. Using the RMHD simulation we can predict and interpret the observed properties of star-forming clouds, first cores and outflows with millimeter/submillimeter radio interferometers, especially the Atacama Large Millimeter/submillimeter Array (ALMA).

We have performed MHD simulations of Alfven wave propagation along an open flux tube in the solar atmosphere. In our numerical model, Alfven waves are generated by the photospheric granular motion. As the wave generator, we used a derived temporal spectrum of the photospheric granular motion from G-band movies of Hinode/SOT. It is shown that the total energy flux at the corona becomes larger and the transition region height becomes higher in the case when we use the observed spectrum rather than white/pink noise spectrum as the wave generator. This difference can be explained by the Alfven wave resonance between the photosphere and the transition region. After performing Fourier analysis on our numerical results, we have found that the region between the photosphere and the transition region becomes an Alfven wave resonant cavity. We have confirmed that there are at least three resonant frequencies, 1, 3 and 5 mHz, in our numerical model. Alfven wave resonance is one of the most effective mechanisms to explain the dynamics of the spicules and the sufficient energy flux to heat the corona.

The formation and evolution of the circumstellar disk in unmagnetized molecular clouds is investigated using three-dimensional hydrodynamic simulations from the prestellar core until the end of the main accretion phase. In collapsing clouds, the first (adiabatic) core with a size of ~10AU forms prior to the formation of the protostar. At its formation, the first core has a thick disk-like structure, and is mainly supported by the thermal pressure. After the protostar formation, it decreases the thickness gradually, and becomes supported by the centrifugal force. We found that the first core is a precursor of the circumstellar disk. This indicates that the circumstellar disk is formed before the protostar formation with a size of ~10AU, which means that no protoplanetary disk smaller than <10AU exists. Reflecting the thermodynamics of the collapsing gas, at the protostar formation epoch, the circumstellar disk has a mass of ~0.01-0.1 solar mass, while the protostar has a mass of ~10^-3 solar mass. Thus, just after the protostar formation, the circumstellar disk is about 10-100 times more massive than the protostar. Even in the main accretion phase that lasts for ~10^5yr, the circumstellar disk mass dominates the protostellar mass. Such a massive disk is unstable to gravitational instability, and tends to show fragmentation. Our calculations indicate that the planet or brown-dwarf mass object may form in the circumstellar disk in the main accretion phase. In addition, the mass accretion rate onto the protostar shows strong time variability that is caused by the perturbation of proto-planets and/or the spiral arms in the circumstellar disk. Such variability provides a useful signature for detecting the planet-sized companion in the circumstellar disk around very young protostars.

Star formation in relic HII regions of the first stars is investigated using magneto-hydrodynamical simulations with a nested grid method that covers 10 orders of magnitude in spatial scale and 20 orders of magnitude in density contrast. Due to larger fraction of H_2 and HD molecules, its prestellar thermal evolution is considerably different from that in the first star formation. Reflecting the difference, two hydrostatic cores appear in a nested manner: a protostar is enclosed by a transient hydrostatic core, which appears during the prestellar collapse. If the initial natal core rotates fast at a rate with rotational to gravitational energy ratio \beta_0 > 0.01-0.1, the transient hydrostatic core fragments to \sim 10 M_\odot sub-cores at density \sim 10^9 cm^-3. With smaller rotation energy, fragmentation occurs at higher density while a single protostar forms without fragmentation if rotation is extremely slow with \beta_0 < 10^-6 -10^-5. If magnetic field is present, these threshold values of \beta_0 is boosted owing to angular momentum transport by the magnetic breaking. Magnetic field also drives the protostellar outflows. With strong magnetic field, two distinct outflows are observed: The slower one emanates from the transient hydrostatic core, while the faster one from the protostar. These flows may affect the final stellar mass by ejecting some of masses in the initial core, and also may play some role in driving and maintenance of interstellar turbulence in young galaxies.

Brown dwarf formation and star formation efficiency are studied using a nested grid simulation that covers five orders of magnitude in spatial scale (10^4 - 0.1AU). Starting with a rotating magnetized compact cloud with a mass of 0.22 M_sun, we follow the cloud evolution until the end of main accretion phase. Outflow of about 5 km/s emerges about 100 yr before the protostar formation and does not disappear until the end of the calculation. The mass accretion rate declines from 10^-6 M_sun/yr to 10^-8 - 10^-12 M_sun/yr in a short time (about 10^4 yr) after the protostar formation. This is because (1) a large fraction of mass is ejected from the host cloud by the protostellar outflow and (2) the gas escapes from the host cloud by the thermal pressure. At the end of the calculation, 74% (167 M_Jup) of the total mass (225 M_Jup) is outflowing from the protostar, in which 34% (77 M_Jup) of the total mass is ejected by the protostellar outflow with supersonic velocity and 40% (90 M_Jup) escapes with subsonic velocity. On the other hand, 20% (45 M_Jup) is converted into the protostar and 6% (13 M_Jup) remains as the circumstellar disk. Thus, the star formation efficiency is epsilon = 0.2. The resultant protostellar mass is in the mass range of brown dwarfs. Our results indicate that brown dwarfs can be formed in compact cores in the same manner as hydrogen-burning stars, and the magnetic field and protostellar outflow are essential in determining the star formation efficiency and stellar mass.

The fragmentation process in collapsing clouds with various metallicities is studied using three-dimensional nested-grid hydrodynamics. Initial clouds are specified by three parameters: cloud metallicity, initial rotation energy and initial cloud shape. For different combinations of these parameters, we calculate 480 models in total and study cloud evolution, fragmentation conditions, orbital separation and binary frequency. For the cloud to fragment during collapse, the initial angular momentum must be higher than a threshold value, which decreases with decreasing metallicity. Although the exact fragmentation conditions depend also on the initial cloud shape, this dependence is only modest. Our results indicate a higher binary frequency in lower-metallicity gas. In particular, with the same median rotation parameter as in the solar neighbourhood, a majority of stars are born as members of binary/multiple systems for < 10^-4 Z_sun. With initial mass < 0.1 M_sun, if fragments are ejected in embryo from the host clouds by multi-body interaction, they evolve to substellar-mass objects. This provides a formation channel for low-mass stars in zero- or low-metallicity environments.

The protostellar outflows have indispensable role in the formation of single stars, because they carry off the excess angular momentum from the centre of the shrinking gas cloud, and permits further collapse to form a star. On the other hand, a significant fraction of stars is supposed to be born as binaries with circumbinary disk that are frequently observed. Here, we investigate the evolution of a magnetized rotating cloud using three-dimensional resistive MHD nested-grid code, and show that the outflow is driven by the circumbinary disk and has an important role even in the binary formation. After the adiabatic core formation in the collapsing cloud core, the magnetic flux is significantly removed from the centre of the cloud by the Ohmic dissipation. Since this removal makes the magnetic braking ineffective, the adiabatic core continuously acquires the angular momentum to induce fragmentation and subsequent binary formation. The magnetic field accumulates in the circumbinary disk where the removal and accretion of magnetic field are balanced, and finally drives circumbinary outflow. This result explains the spectacular morphology of some specific young stellar objects such as L1551 IRS5. We can infer that most of the bipolar molecular outflows observed by low density tracers (i.e., CO) would correspond to circumbinary or circum-multiple outflows found in this report, since most of the young stellar objects are supposed to be binaries or multiples.

We investigate the driving mechanism of outflows and jets in star formation process using resistive MHD nested grid simulations. We found two distinct flows in the collapsing cloud core: Low-velocity outflows (sim 5 km/s) with a wide opening angle, driven from the first adiabatic core, and high-velocity jets (sim 50 km/s) with good collimation, driven from the protostar. High-velocity jets are enclosed by low-velocity outflow. The difference in the degree of collimation between the two flows is caused by the strength of the magnetic field and configuration of the magnetic field lines. The magnetic field around an adiabatic core is strong and has an hourglass configuration. Therefore, the low-velocity outflow from the adiabatic core are driven mainly by the magnetocentrifugal mechanism and guided by the hourglass-like field lines. In contrast, the magnetic field around the protostar is weak and has a straight configuration owing to Ohmic dissipation in the high-density gas region. Therefore, high-velocity jet from the protostar are driven mainly by the magnetic pressure gradient force and guided by straight field lines. Differing depth of the gravitational potential between the adiabatic core and the protostar cause the difference of the flow speed. Low-velocity outflows correspond to the observed molecular outflows, while high-velocity jets correspond to the observed optical jets. We suggest that the protostellar outflow and the jet are driven by different cores (the first adiabatic core and protostar), rather than that the outflow being entrained by the jet.

The Cosmic Infrared Background ExpeRiment (CIBER) is a rocket-borne absolute photometry imaging and spectroscopy experiment optimized to detect signatures of first-light galaxies present during reionization in the unresolved IR background. CIBER-I consists of a wide-field two-color camera for fluctuation measurements, a low-resolution absolute spectrometer for absolute EBL measurements, and a narrow-band imaging spectrometer to measure and correct scattered emission from the foreground zodiacal cloud. CIBER-I was successfully flown on February 25th, 2009 and has one more planned flight in early 2010. We propose, after several additional flights of CIBER-I, an improved CIBER-II camera consisting of a wide-field 30 cm imager operating in 4 bands between 0.5 and 2.1 microns. It is designed for a high significance detection of unresolved IR background fluctuations at the minimum level necessary for reionization. With a FOV 50 to 2000 times largerthan existing IR instruments on satellites, CIBER-II will carry out the definitive study to establish the surface density of sources responsible for reionization.

(Brief Summary) What is the total radiative content of the Universe since the epoch of recombination? The extragalactic background light (EBL) spectrum captures the redshifted energy released from the first stellar objects, protogalaxies, and galaxies throughout cosmic history. Yet, we have not determined the brightness of the extragalactic sky from UV/optical to far-infrared wavelengths with sufficient accuracy to establish the radiative content of the Universe to better than an order of magnitude. Among many science topics, an accurate measurement of the EBL spectrum from optical to far-IR wavelengths, will address: What is the total energy released by stellar nucleosynthesis over cosmic history? Was significant energy released by non-stellar processes? Is there a diffuse component to the EBL anywhere from optical to sub-millimeter? When did first stars appear and how luminous was the reionization epoch? Absolute optical to mid-IR EBL spectrum to an astrophysically interesting accuracy can be established by wide field imagingat a distance of 5 AU or above the ecliptic plane where the zodiacal foreground is reduced by more than two orders of magnitude.

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.

We show three-dimensional magnetohydrodynamical simulations of core collapse supernova in which the progenitor has magnetic fields inclined to the rotation axis. The simulations employed a simple empirical equation of state in which the pressure of degenerate gas is approximated by piecewise polytropes for simplicity. Neither energy loss due to neutrino is taken into account for simplicity. The simulations start from the stage of dynamical collapse of an iron core. The dynamical collapse halts at $ t $ = 189 ms by the pressure of high density gas and a proto-neutron star (PNS) forms. The evolution of PNS was followed about 40 milli-seconds in typical models. When the initial rotation is mildly fast and the initial magnetic fields are mildly strong, bipolar jets are launched from an upper atmosphere ($ r \sim 60 {\rm km} $) of the PNS. The jets are accelerated to $ \sim 3 \times 10 ^4 $ km s$^{-1}$, which is comparable to the escape velocity at the foot point. The jets are parallel to the initial rotation axis. Before the launch of the jets, magnetic fields are twisted by rotation of the PNS. The twisted magnetic fields form torus-shape multi-layers in which the azimuthal component changes alternately. The formation of magnetic multi-layers is due to the initial condition in which the magnetic fields are inclined with respect to the rotation axis. The energy of the jet depends only weakly on the initial magnetic field assumed. When the initial magnetic fields are weaker, the time lag is longer between the PNS formation and jet ejection. It is also shown that the time lag is related to the Alfv\'en transit time. Although the nearly spherical prompt shock propagates outward in our simulations, it is ...

Jet driving and fragmentation process in collapsing primordial cloud are studied using three-dimensional MHD nested grid simulations. Starting from a rotating magnetized spherical cloud with the number density of n=10^3 cm^-3, we follow the evolution of the cloud up to the stellar density n=10^22 cm^-3. We calculate 36 models parameterizing the initial magnetic and rotational energies (\gamma_0, \beta_0). In the collapsing primordial clouds, the cloud evolutions are characterized by the ratio of the initial rotational to magnetic energy, \gamma_0/\beta_0. The Lorentz force significantly affects the cloud evolution when \gamma_0 > \beta_0, while the centrifugal force is more dominant than the Lorentz force when \beta_0 > \gamma_0. When the cloud rotates rapidly with angular velocity of \Omega_0 > 10^-17 (n/10^3 cm^-3)^2/3 s^-1 and \beta_0 > \gamma_0, fragmentation occurs before the protostar is formed, but no jet appears after the protostar formation. On the other hand, a strong jet appears after the protostar formation without fragmentation when the initial cloud has the magnetic field of B_0 > 10^-9 (n/10^3 cm^-3)^2/3 G and \gamma_0 > \beta_0. Our results indicate that proto-Population III stars frequently show fragmentation and protostellar jet. Population III stars are therefore born as binary or multiple stellar systems, and they can drive strong jets, which disturb the interstellar medium significantly, as well as in the present-day star formation, and thus they may induce the formation of next generation stars.

We investigate the accretion of angular momentum onto a protoplanet system using three-dimensional hydrodynamical simulations. We consider a local region around a protoplanet in a protoplanetary disk with sufficient spatial resolution. We describe the structure of the gas flow onto and around the protoplanet in detail. We find that the gas flows onto the protoplanet system in the vertical direction crossing the shock front near the Hill radius of the protoplanet, which is qualitatively different from the picture established by two-dimensional simulations. The specific angular momentum of the gas accreted by the protoplanet system increases with the protoplanet mass. At Jovian orbit, when the protoplanet mass M_p is M_p < 1 M_J, where M_J is Jovian mass, the specific angular momentum increases as j \propto M_p. On the other hand, it increases as j \propto M_p^2/3 when the protoplanet mass is M_p > 1 M_J. The stronger dependence of the specific angular momentum on the protoplanet mass for M_p < 1 M_J is due to thermal pressure of the gas. The estimated total angular momentum of a system of a gas giant planet and a circumplanetary disk is two-orders of magnitude larger than those of the present gas giant planets in the solar system. A large fraction of the total angular momentum contributes to the formation of the circumplanetary disk. We also discuss the satellite formation from the circumplanetary disk.

High resolution imaging observation of a sunspot umbra was done with Hinode Solar Optical Telescope (SOT). Filtergrams in wavelengths of blue and green continuum were taken during three consecutive days. The umbra consisted of a dark core region, several diffuse components and numerous umbral dots. We derived basic properties of umbral dots (UDs), especially their temperatures, lifetimes, proper motions, spatial distribution and morphological evolution. Brightness of UDs is confirmed to depend on the brightness of their surrounding background. Several UDs show fission and fusion. Thanks to the stable condition of space observation, we could first follow the temporal behavior of these events. The derived properties of internal structure of the umbra are discussed in viewpoint of magnetoconvection in a strong magnetic field.

The fragmentation process of primordial-gas cores during prestellar collapse is studied using three-dimensional nested-grid hydrodynamics. Starting from the initial central number density of n \sim10^3 cm^-3, we follow the evolution of rotating spherical cores up to the stellar density n \simeq 10^{22} cm^-3. An initial condition of the cores is specified by three parameters: the ratios of the rotation and thermal energies to the gravitational energy (\beta_0, and \alpha_0, respectively), and the amplitude of the bar-mode density perturbation (A_\phi). Cores with rotation \beta_0 > 10^{-6} are found to fragment during the collapse. The fragmentation condition hardly depends on either the initial thermal energy \alpha_0 or amplitude of bar-mode perturbation A_\phi. Since the critical rotation parameter for fragmentation is lower than that expected in the first star formation, binaries or multiples are also common for the first stars.

A new numerical code, called SFUMATO, for solving self-gravitational magnetohydrodynamics (MHD) problems using adaptive mesh refinement (AMR) is presented. A block-structured grid is adopted as the grid of the AMR hierarchy. The total variation diminishing (TVD) cell-centered scheme is adopted as the MHD solver, with hyperbolic cleaning of divergence error of the magnetic field also implemented. The self-gravity is solved by a multigrid method composed of (1) full multigrid (FMG)-cycle on the AMR hierarchical grids, (2) V-cycle on these grids, and (3) FMG-cycle on the base grid. The multigrid method exhibits spatial second-order accuracy, fast convergence, and scalability. The numerical fluxes are conserved by using a refluxing procedure in both the MHD solver and the multigrid method. The several tests are performed indicating that the solutions are consistent with previously published results.

We observed small scale magnetic flux emergence in a sunspot moat region by the Solar Optical Telescope (SOT) aboard the Hinode satellite. We analyzed filtergram images observed in the wavelengths of Fe 6302 angstrom, G-band and Ca II H. In Stokes I images of Fe 6302 angstrom, emerging magnetic flux were recognized as dark lanes. In G-band, they showed their shapes almost the same as in Stokes I images. These magnetic flux appeared as dark filaments in Ca II H images. Stokes V images of Fe 6302 angstrom showed pairs of opposite polarities at footpoints of each filament. These magnetic concentrations are identified to correspond to bright points in G-band/Ca II H images. From the analysis of time-sliced diagrams, we derived following properties of emerging flux, which are consistent with the previous works. (1) Two footpoints separate each other at a speed of 4.2 km/s during the initial phase of evolution and decreases to about 1 km/s in 10 minutes later. (2) Ca II H filaments appear almost simultaneously with the formation of dark lanes in Stokes I in the observational cadence of 2 minutes. (3) The lifetime of the dark lanes in Stokes I and G-band is 8 minutes, while that of Ca filament is 12 minutes. An interesting phenomena was observed that an emerging flux tube expands laterally in the photosphere with a speed of 3.8 km/s. Discussion on the horizontal expansion of flux tube will be given with refernce to previous simulation studies.

Fragmentation and binary formation processes are studied using three-dimensional resistive MHD nested grid simulations. Starting with a Bonnor-Ebert isothermal cloud rotating in a uniform magnetic field, we calculate the cloud evolution from the molecular cloud core (n=10^4 cm^-3) to the stellar core (n \simeq 10^22 cm^-3). We calculated 147 models with different initial magnetic, rotational, and thermal energies, and the amplitudes of the non-axisymmetric perturbation. In a collapsing cloud, fragmentation is mainly controlled by the initial ratio of the rotational to the magnetic energy, regardless of the initial thermal energy and amplitude of the non-axisymmetric perturbation. When the clouds have large rotational energies in relation to magnetic energies, fragmentation occurs in the low-density evolution phase (10^12 cm^-3 < n < 10^15 cm^-3) with separations of 3-300 AU. Fragments that appeared in this phase are expected to evolve into wide binary systems. On the other hand, fragmentation does not occur in the low-density evolution phase, when initial clouds have large magnetic energies in relation to the rotational energies. In these clouds, fragmentation only occurs in the high-density evolution phase (n > 10^17 cm^-3) after the clouds experience significant reduction of the magnetic field owing to Ohmic dissipation in the period of 10^12 cm^-3 < n < 10^15 cm^-3. Fragments appearing in this phase have separations of < 0.3 AU, and are expected to evolve into close binary systems. As a result, we found two typical fragmentation epochs, which cause different stellar separations. Although these typical separations are disturbed in the subsequent gas accretion phase, we might be able to observe two peaks of binary separations in extremely young stellar groups.

The Infrared Camera (IRC) is one of the two instruments on board the AKARI satellite. In addition to deep imaging from 1.8-26.5um for the pointed observation mode of the AKARI, it has a spectroscopic capability in its spectral range. By replacing the imaging filters by transmission-type dispersers on the filter wheels, it provides low-resolution (lambda/d_lambda ~ 20-120) spectroscopy with slits or in a wide imaging field-of-view (approximately 10'X10'). The IRC spectroscopic mode is unique in space infrared missions in that it has the capability to perform sensitive wide-field spectroscopic surveys in the near- and mid-infrared wavelength ranges. This paper describes specifications of the IRC spectrograph and its in-orbit performance.

We present the results of optical identifications for 257 mid-infrared sources detected with a deep 15um survey over approximately 80 arcmin^2 area in the AKARI performance verification field near the North Ecliptic Pole. The 15um fluxes of the sources range from 1 mJy down to 40 uJy, approximately a half of which are below 100 uJy. Optical counterparts were searched for within a 2-3 arcsec radius in both the BVRi'z' catalog generated by using the deep Subaru/Suprime-cam field which covers one-third of the performance verification field, and the g'r'i'z' catalog based on observations made with MegaCam at CFHT. We found B-R and R-z' colours of sources with successful optical identifications are systematically redder than that of the entire optical sample in the same field. Moreover, approximately 40% of the 15um sources show colours R-L15>5, which cannot be explained by the spectral energy distribution (SED) of normal quiescent spiral galaxies, but are consistent with SEDs of redshifted (z>1) starburst or ultraluminous infrared galaxies. This result indicates that the fraction of the ultraluminous infrared galaxies in our faint 15um sample is much larger than that in our brighter 15um sources, which is consistent with the evolving mid-infrared luminosity function derived by recent studies based on the Spitzer 24um deep surveys. Based on an SED fitting technique, the nature of the faint 15um sources is further discussed for a selected number of sources with available K_s-band data.

We report the first results of a near- and mid- infrared deep survey with the Infrared Camera (IRC) onboard AKARI in the performance verification phase. Simultaneous observations by the NIR, MIR-S and MIR-L channels of the IRC with effective integration times of 4529, 4908, and 4417 seconds at 3, 7, and 15 micron, covering 86.0, 70.3, and 77.3 arcmin^2 area, detected 955, 298 and 277 sources, respectively. The 5 sigma detection limits of the survey are 6.0, 31.5 and 71.2 micro Jy and the 50% completeness limit are 24.0, 47.5, and 88.1 micro Jy at 3, 7, and 15 micron, respectively. The observation is limited by source confusion at 3 micron. We have confirmed the turnover in the 15 micron differential source counts around 400 micro Jy, previously detected by surveys with the Infrared Space Observatory. The faint end of 15 micron raw source counts agree with the results from the deep surveys in the GOODS fields carried out with the Spitzer IRS peak up imager and the predictions of current galaxy evolution models. These results indicate that deep surveys with comprehensive wavelength coverage at mid-infrared wavelength are very important to investigate the evolution of infrared galaxies at high redshifts.

We present the observations of the reflection nebulae IC4954 and IC4955 region with the Infrared Camera (IRC) and the Far-Infrared Surveyor (FIS) on board the infrared astronomical satellite AKARI during its performance verification phase. We obtained 7 band images from 7 to 160um with higher spatial resolution and higher sensitivities than previous observations. The mid-infrared color of the S9W (9um) and L18W (18um) bands shows a systematic variation around the exciting sources. The spatial variation in the mid-infrared color suggests that the star-formation in IC4954/4955 is progressing from south-west to north-east. The FIS data also clearly resolve two nebulae for the first time in the far-infrared. The FIS 4-band data from 65um to 160um allow us to correctly estimate the total infrared luminosity from the region, which is about one sixth of the energy emitted from the existing stellar sources. Five candidates for young stellar objects have been detected as point sources for the first time in the 11um image. They are located in the red S9W to L18W color regions, suggesting that current star-formation has been triggered by previous star-formation activities. A wide area map of the size of about 1 x 1 (deg^2) around the IC4954/4955 region was created from the AKARI mid-infrared all-sky survey data. Together with the HI 21cm data, it suggests a large hollow structure of a degree scale, on whose edge the IC4954/4955 region has been created, indicating star formation over three generations in largely different spatial scales.

We describe the flight performance of the cryogenic system of the infrared astronomical satellite AKARI, which was successfully launched on 2006 February 21 (UT). AKARI carries a 68.5 cm telescope together with two focal plane instruments, Infrared Cameras (IRC) and Far Infrared Surveyor (FIS), all of which are cooled down to cryogenic temperature to achieve superior sensitivity. The AKARI cryogenic system is a unique hybrid system, which consists of cryogen (liquid helium) and mechanical coolers (2-stage Stirling coolers). With the help of the mechanical coolers, 179 L (26.0 kg) of super-fluid liquid helium can keep the instruments cryogenically cooled for more than 500 days. The on-orbit performance of the AKARI cryogenics is consistent with the design and pre-flight test, and the boil-off gas flow rate is as small as 0.32 mg/s. We observed the increase of the major axis of the AKARI orbit, which can be explained by the thrust due to thermal pressure of vented helium gas.

AKARI, the first Japanese satellite dedicated to infrared astronomy, was launched on 2006 February 21, and started observations in May of the same year. AKARI has a 68.5 cm cooled telescope, together with two focal-plane instruments, which survey the sky in six wavelength bands from the mid- to far-infrared. The instruments also have the capability for imaging and spectroscopy in the wavelength range 2 - 180 micron in the pointed observation mode, occasionally inserted into the continuous survey operation. The in-orbit cryogen lifetime is expected to be one and a half years. The All-Sky Survey will cover more than 90 percent of the whole sky with higher spatial resolution and wider wavelength coverage than that of the previous IRAS all-sky survey. Point source catalogues of the All-Sky Survey will be released to the astronomical community. The pointed observations will be used for deep surveys of selected sky areas and systematic observations of important astronomical targets. These will become an additional future heritage of this mission.

The evolution of the magnetic field and angular momentum in the collapsing cloud core is studied using three-dimensional resistive MHD nested grid simulations. Starting with a Bonnor-Ebert isothermal cloud rotating in a uniform magnetic field, we calculate the cloud evolution from the molecular cloud core (n=10^4 cm^-3) to the stellar core (n \simeq 10^22 cm^-3). The magnetic field strengths at the center of the clouds converge to a certain value as the clouds collapse, when the clouds have the same angular momenta but different strengths of the magnetic fields at the initial state. For 10^12 cm^-3 < n < 10^16 cm^-3, Ohmic dissipation considerably removes the magnetic field from the collapsing cloud core, and the magnetic field lines, which are strongly twisted for n <10^12 cm^-3, is de-collimated. The magnetic field lines are twisted and amplified again for nc > 10^16 cm^-3, because the magnetic field is recoupled with the warm gas. Finally, protostars at their formation epoch have 0.1-1kG of the magnetic fields, which are comparable to observations. The magnetic field strength of protostar slightly depends on the angular momentum of the host cloud. The protostar formed from the slowly rotating cloud core has a stronger magnetic field. The evolution of the angular momentum is closely related to the evolution of the magnetic field. The angular momentum in the collapsing cloud is removed by the magnetic effect. The formed protostars have 0.1-2 days of the rotation period at their formation epoch, which are slightly shorter than the observation. This indicates that the further removal mechanism of the angular momentum such as interaction between the protostar and disk, wind gas or jet is important in further evolution of the protostar.

We report the detection of an H-alpha emission line in the low resolution spectrum of a quasar, RX J1759.4+6638, at a redshift of 4.3 with the InfraRed Camera (IRC) onboard the AKARI. This is the first spectroscopic detection of an H-alpha emission line in a quasar beyond z=4. The overall spectral energy distribution (SED) of RX J1759.4+6638 in the near- and mid-infrared wavelengths agrees with a median SED of the nearby quasars and the flux ratio of F(Ly-alpha)/F(H-alpha) is consistent with those of previous reports for lower-redshift quasars.

We report preliminary results of AKARI observations of two globular clusters, NGC104 and NGC362. Imaging data covering areas of about 10x10 arcmin^2 centered on the two clusters have been obtained with InfraRed Camera (IRC) at 2.4, 3.2, 4.1, 7.0, 9.0, 11.0, 15.0, 18.0 and 24.0 mu. We used F11/F2 and F24/F7 flux ratios as diagnostics of circumstellar dust emission. Dust emissions are mainly detected from variable stars obviously on the asymptotic giant branch, but some variable stars that reside below the tip of the first-ascending giant branch also show dust emissions. We found eight red sources with F24/F7 ratio greater than unity in NGC362. Six out of the eight have no 2MASS counterparts. However, we found no such source in NGC104.

The Infrared Camera (IRC) is one of two focal-plane instruments on the AKARI satellite. It is designed for wide-field deep imaging and low-resolution spectroscopy in the near- to mid-infrared (1.8--26.5um) in the pointed observation mode of AKARI. IRC is also operated in the survey mode to make an all-sky survey at 9 and 18um. It comprises three channels. The NIR channel (1.8--5.5um) employs a 512 x 412 InSb array, whereas both the MIR-S (4.6--13.4um) and MIR-L (12.6--26.5um) channels use 256 x 256 Si:As impurity band conduction arrays. Each of the three channels has a field-of-view of about 10' x 10' and are operated simultaneously. The NIR and MIR-S share the same field-of-view by virtue of a beam splitter. The MIR-L observes the sky about $25' away from the NIR/MIR-S field-of-view. IRC gives us deep insights into the formation and evolution of galaxies, the evolution of planetary disks, the process of star-formation, the properties of interstellar matter under various physical conditions, and the nature and evolution of solar system objects. The in-flight performance of IRC has been confirmed to be in agreement with the pre-flight expectation. This paper summarizes the design and the in-flight operation and imaging performance of IRC.

The driving mechanism of jets and outflows in star formation process is studied using resistive MHD nested grid simulations. We calculated cloud evolution from the molecular cloud core to the stellar core. In the collapsing cloud core, we found two distinct flows: Low-velocity flows (sim 5 km/s) with a wide opening angle, driven from the adiabatic core, and high-velocity flows (sim 30 km/s) with good collimation, driven from the protostar. High-velocity flows are enclosed by low-velocity flows after protostar formation. The difference in the degree of collimation between the two flows is caused by the strength of the magnetic field and configuration of the magnetic field lines. The magnetic field around an adiabatic core is strong and has an hourglass configuration; therefore, flows from the adiabatic core are driven mainly by the magnetocentrifugal mechanism and guided by the hourglass-like field lines. In contrast, the magnetic field around the protostar is weak and has a straight configuration owing to Ohmic dissipation in the high-density gas region. Therefore, flows from the protostar are driven mainly by the magnetic pressure gradient force and guided by straight field lines. Differing depth of the gravitational potential between the adiabatic core and the protostar cause the difference of the flow speed. Low-velocity flows correspond to the observed molecular outflows, while high-velocity flows correspond to the observed optical jets. We suggest that the outflow and the jet are driven by different cores, rather than that the outflow being entrained by the jet.

We present an initial analysis of AKARI 18um-selected galaxies using all 9 photometric bands at 2-24um available in the InfraRed Camera (IRC), in order to demonstrate new capabilities of AKARI cosmological surveys. We detected 72 sources at 18um in an area of 50.2 arcmin^2 in the AKARI/IRC monitor field towards the North Ecliptic Pole (NEP). From this sample, 25 galaxies with probable redshifts z>~ 0.5 are selected with a single colour cut (N2-N3>0.1) for a detailed SED analysis with ground-based BVRi'z'JK data. Using an SED radiative transfer model of starbursts covering the wavelength range UV -- submm, we derive photometric redshifts from the optical-MIR SEDs of 18um-selected galaxies. From the best-fit SED models, we show that the IRC all-band photometry is capable of tracing the steep rise in flux at the blue side of the PAH 6.2um emission feature. This indicates that the IRC all-band photometry is useful to constrain the redshift of infrared galaxies, specifically for dusty galaxies with a less prominent 4000A break. Also, we find that the flux dip between the PAH 7.7 and 11.2um emission feature is recognizable in the observed SEDs of galaxies at z~1. By using such a colour anomaly due to the PAH and silicate absorption features, unique samples of ULIRGs at z~1, `silicate-break' galaxies, can be constructed from large cosmological surveys of AKARI towards the NEP, i.e. the NEP-Deep and NEP-Wide survey. This pilot study suggests the possibility of detecting many interesting galaxy properties in the NEP-Deep and Wide surveys, such as a systematic difference in SEDs between high- and low-z ULIRGs, and a large variation of the PAH inter-band strength ratio in galaxies at high redshifts. [abridged]

We have conducted deep JHKs imaging polarimetry of a ~8' x 8' area of the NGC 2071 star forming region. Our polarization data have revealed various infrared reflection nebulae (IRNe) associated with the central IR young star cluster NGC2071IR and identified their illuminating sources. There are at least 4 IRNe in NGC2071IR and several additional IRNe are identified around nearby young stars in the same field-of-view. Each illuminating source coincides with a known near-IR source except for IRS3, which is only a part of IRN2 and is illuminated by the radio source 1c. Aperture polarimetry of each cluster source is used to detect unresolved circumstellar disk/outflow systems. Aperture polarimetry of the other point-like sources within the field is made in this region for the first time. The magnetic field structures (from ~1 pc down to \~0.1 pc) are derived using both aperture polarimetry of the point-like sources and imaging polarimetry of the shocked H2 emission that is seen as the dominant knotty nebulae in the Ks band image; they are both of dichroic origin and the derived field directions are consistent with each other. The magnetic field direction projected on the sky is also consistent with that inferred from the 850 micron thermal continuum emission polarimetry of the central 0.2 pc region, but running roughly perpendicular (~75 degrees) to the direction of the large scale outflow. We argue that the field strength is too weak to align the outflow in the large scale field direction via magnetic braking.

We investigate outflows driven by a giant protoplanet using three-dimensional MHD nested grid simulations. We consider a local region around the protoplanet in the protoplanetary disk, and calculate three models: (a) unmagnetized disk model, (b) magnetized disk model having magnetic field azimuthally parallel to the disk, and (c) magnetic field perpendicular to the disk. Outflows with velocities, at least, 10 km/s are driven by the protoplanets in both magnetized disk models, while outflow does not appear in unmagnetized disk model. Tube-like outflows along the azimuthal direction of the protoplanetary disk appear in model with magnetic field being parallel to the disk. In this model, the magnetically dominated regions (i.e., density gap) are clearly contrasted from other regions and spiral waves appear near the protoplanet. On the other hand, in model with magnetic field being perpendicular to the disk, outflows are driven by a protoplanet with cone-like structure just as seen in the outflow driven by a protostar. Magnetic field lines are strongly twisted near the protoplanet and the outflows have well-collimated structures in this model.These outflows can be landmarks for searching exo-protoplanets in their formation stages. Our results indicate that the accretion rate onto the protoplanet tend to have a larger value than that expected from previous hydrodynamical calculations, since a fraction of the angular momentum of circum-planetary disk is removed by outflows, enhanced non-axisymmetric patterns caused by magnetic field, and magnetic braking. Possible implications for observation are also briefly discussed.

The protostellar jets driven by the formation of the first stars are studied by using MHD nested grid simulations. Starting from a slowly rotating spherical cloud of 5.1 times 10^4 Msun permeated by a uniform magnetic field, we follow the evolution from the central number density n = 10^3 cm^-3 to n simeq 10^23 cm^-3. We calculate four models that differ in initial magnetic field strengths and angular velocities. In all models, protostars of simeq 10^-3 Msun are formed at n simeq 10^22 cm^-3 in accordance with one-dimensional calculations. By this epoch, the magnetic flux density is amplified by 10 orders of magnitude from the initial value. Consequently, the formed protostar possesses the magnetic field of \sim 10^6 G that is much larger than the flux density of the present counterparts, reflecting the fact that the dissipation of a magnetic field is ineffective in primordial gas clouds. If the initial magnetic field B > 10^-9 (n/10^3 cm^-3)^2/3 G, the protostellar jet is launched and its velocities reaches \sim 70 km s^-1 by the time the protostellar mass becomes (4-6) times 10^-3 Msun, and a fraction (3-10%) of the accreting matter is blown off from the central region. Owing to the interaction of these ejecta with surrounding matter, expanding bow shocks are created at both heads of the jet. If this jet continues to sweep out the surrounding gas that otherwise accretes onto the central star or circumstellar disk, the final mass of the first star can be substantially reduced. In addition, dense post-shock regions behind the bow shocks are expected to promote the chemical reactions (formation of H_2 and HD), and this provides possible environments for subsequent low-mass star formation in the early universe.

The stellar core formation and high speed jets driven by the formed core are studied by using three-dimensional resistive MHD nested grid simulations. Starting with a Bonnor-Ebert isothermal cloud rotating in a uniform magnetic field, we calculate the cloud evolution from the molecular cloud core (n = 10^6 cm^-3, r_c = 4.6 times 10^4 AU) to the stellar core (n \simeq 10^23 cm^-3, r_c \simeq 1 solar radius). We resolve cloud structure over 7 orders of magnitude in spatial extent and over 17 orders of magnitude in density contrast. For comparison, we calculate two models: resistive and ideal MHD models. Both models have the same initial condition, but the former includes dissipation process of magnetic field while the latter does not. The magnetic fluxes in resistive MHD model are extracted from the first core during 10^12 cm^-3 < n < 10^16 cm^-3 by Ohmic dissipation. Magnetic flux density of the formed stellar core (n \simeq 10^20 cm^-3) in resistive MHD model is two orders of magnitude smaller than that in ideal MHD model. Since magnetic braking is less effective in resistive MHD model, rapidly rotating stellar core (the second core) is formed. After stellar core formation, the magnetic field of the core is largely amplified both by magneto-rotational instability and the shearing motion between the stellar core and ambient medium. As a consequence, high speed (simeq 45 km,s^-1) jets are driven by the second core, which results in strong mass ejection. A cocoon-like structure around the second core also forms with clear bow shocks.

We studied the collapse of rotating molecular cloud cores with inclined magnetic fields, based on three-dimensional numerical simulations.The numerical simulations start from a rotating Bonnor-Ebert isothermal cloud in a uniform magnetic field. The magnetic field is initially taken to be inclined from the rotation axis. As the cloud collapses, the magnetic field and rotation axis change their directions. When the rotation is slow and the magnetic field is relatively strong, the direction of the rotation axis changes to align with the magnetic field, as shown earlier by Matsumoto & Tomisaka. When the magnetic field is weak and the rotation is relatively fast, the magnetic field inclines to become perpendicular to the rotation axis. In other words, the evolution of the magnetic field and rotation axis depends on the relative strength of the rotation and magnetic field. Magnetic braking acts to align the rotation axis and magnetic field, while the rotation causes the magnetic field to incline through dynamo action. The latter effect dominates the former when the ratio of the angular velocity to the magnetic field is larger than a critical value \Omega_0/ B_0 > 0.39 G^1/2 c_s^-1, where B_0, \Omega_0, G, and c_s^-1 denote the initial magnetic field, initial angular velocity, gravitational constant, and sound speed, respectively. When the rotation is relatively strong, the collapsing cloud forms a disk perpendicular to the rotation axis and the magnetic field becomes nearly parallel to the disk surface in the high density region. A spiral structure appears due to the rotation and the wound-up magnetic field in the disk.

We estimate the polarized thermal dust emission from MHD simulations of protostellar collapse and outflow formation in order to investigate alignment of outflows with magnetic fields. The polarization maps indicate that alignment of an outflow with the magnetic field depends on the field strength inside the cloud core; the direction of the outflow, projected on the plane of the sky, is aligned preferentially with the mean polarization vector for a cloud core with a magnetic field strength of 80 microgauss, while it does not tend to be aligned for 50 microgauss as long as the 1000 AU scale is considered. The direction of the magnetic field at the cloud center is probed by the direction of the outflow. In addition, the magnetic field at the cloud center can be revealed by ALMA even when the source is embedded deeply in the envelope. The Chandrasekhar-Fermi formula is examined using the polarization maps, indicating that the field strength predicted by the formula should be corrected by a factor of 0.24 - 0.44. The correction factor has a tendency to be lower for a cloud core with a weaker magnetic field.

We are developing a rocket-borne instrument (the Cosmic Infrared Background ExpeRiment, or CIBER) to search for signatures of primordial galaxy formation in the cosmic near-infrared extra-galactic background. CIBER consists of a wide-field two-color camera, a low-resolution absolute spectrometer, and a high-resolution narrow-band imaging spectrometer. The cameras will search for spatial fluctuations in the background on angular scales from 7 arcseconds to 2 degrees over a range of angular scales poorly covered by previous experiments. CIBER will determine if the fluctuations reported by the IRTS arise from first-light galaxies or have a local origin. In a short rocket flight CIBER has sensitivity to probe fluctuations 100 times fainter than IRTS/DIRBE. By jointly observing regions of the sky studied by Spitzer and ASTRO-F, CIBER will build a multi-color view of the near-infrared background, accurately assessing the contribution of local (z = 1-3) galaxies to the observed background fluctuations, allowing a deep and comprehensive survey for first-light galaxy background fluctuations. The low-resolution spectrometer will search for a redshifted Lyman cutoff feature between 0.8 - 2.0 microns. The high-resolution spectrometer will trace zodiacal light using the intensity of scattered Fraunhofer lines, providing an independent measurement of the zodiacal emission and a new check of DIRBE zodiacal dust models. The combination will systematically search for the infrared excess background light reported in near-infrared DIRBE/IRTS data, compared with the small excess reported at optical wavelengths.

Subsequent to Paper I, the evolution and fragmentation of a rotating magnetized cloud are studied with use of three-dimensional MHD nested-grid simulations. After the isothermal runaway collapse, an adiabatic gas forms a protostellar first core at the center of the cloud. When the isothermal gas is stable for fragmentation in a contracting disk, the adiabatic core often breaks into several fragments. Conditions for fragmentation and binary formation are studied. All the cores which show fragmentations are geometrically thin, as the diameter-to-thickness ratio is larger than 3. Two patterns of fragmentation are found. (1) When a thin disk is supported by centrifugal force, the disk fragments through a ring configuration (ring fragmentation). This is realized in a fast rotating adiabatic core as Omega >0.2 tau_ff^-1, where Omega and tau_ff represent the angular rotation speed and the free-fall time of the core, respectively. (2) On the other hand, the disk is deformed to an elongated bar in the isothermal stage for a strongly magnetized or rapidly rotating cloud. The bar breaks into 2 - 4 fragments (bar fragmentation). Even if a disk is thin, the disk dominated by the magnetic force or thermal pressure is stable and forms a single compact body. In either ring or bar fragmentation mode, the fragments contract and a pair of outflows are ejected from the vicinities of the compact cores. The orbital angular momentum is larger than the spin angular momentum in the ring fragmentation. On the other hand, fragments often quickly merge in the bar fragmentation, since the orbital angular momentum is smaller than the spin angular momentum in this case. Comparison with observations is also shown.

We discuss evolution of the magnetic flux density and angular velocity in a molecular cloud core, on the basis of three-dimensional numerical simulations, in which a rotating magnetized cloud fragments and collapses to form a very dense optically thick core of > 5 times 10 ^10 cm^-3 . As the density increases towards the formation of the optically thick core, the magnetic flux density and angular velocity converge towards a single relationship between the two quantities. If the core is magnetically dominated its magnetic flux density approaches 1.5 (n/5 times 10^10 cm^-3)^1/2 mG, while if the core is rotationally dominated the angular velocity approaches 2.57 times 10^-3, (n/5 times 10^10 cm^-3)^1/2 yr^-1, where n is the density of the gas. We also find that the ratio of the angular velocity to the magnetic flux density remains nearly constant until the density exceeds 5 times 10^10 cm^-3. Fragmentation of the very dense core and emergence of outflows from fragments are shown in the subsequent paper.

We searched for near-infrared extragalactic background light (EBL) in the data from the Near-InfraRed Spectrometer (NIRS) on the Infrared Telescope in Space (IRTS). After subtracting the contribution of faint stars and the zodiacal component based on modeling, a significant isotropic emission is obtained in the wavelength bands from 1.4 micron to 4.0 micron. The in-band flux amounts to ~ 35 nW m-2 sr-1 which is too bright to be explained by the integrated light from faint galaxies. Significant fluctuations in sky brightness were also detected which can not be explained by fluctuations due to faint stars, zodiacal components and normal galaxies. The excess fluctuation amounts to ~ 1/4 of the excess emission over the integrated light of galaxies. A two-point correlation analysis shows that IRTS/NIRS data has an angular scale of fluctuations of a few degrees. Recent WMAP results of CMB polarization indicate that the reionization of the Universe occurred at z ~ 17 or earlier. The observed near infrared background light could be redshifted UV light from the first generation of massive stars (Pop.III stars) which caused the reionization of the Universe. The spectral jump around 1 micron over the optical EBL suggests that massive star formation terminated at z~9. The observed fluctuations, however, are considerably larger than the theoretical predictions and require a new scenario for the Pop.III era.

This is the first attempt to construct detailed X-ray spectra of clusters of galaxies from the results of high-resolution hydrodynamic simulations and simulate X-ray observations in order to study velocity fields of the intracluster medium (ICM). The hydrodynamic simulations are based on the recently proposed tsunami model, in which cluster cores are affected by bulk motions of the ICM and turbulence is produced. We note that most other solutions of the cooling flow problem also involve the generation of turbulence in cluster cores. From the mock X-ray observations with Astro-E2 XRS, we find that turbulent motion of the ICM in cluster cores could be detected with the satellite. The Doppler shifts of the metal lines could be used to discriminate among turbulence models. The gas velocities measured through the mock observations are consistent with the line-emission weighted values inferred directly from hydrodynamic simulations.

We explore the possibility for the formation of Population III binaries. The collapse of a rotating cylinder is simulated with a three-dimensional, high-resolution nested grid, assuming the thermal history of primordial gas. The simulations are done with dimensionless units, and the results are applicable to low-mass as well as massive systems by scaling with the initial density. We find that if the initial angular momentum is as small as $\beta \approx 0.1$, where $\beta$ is the ratio of centrifugal force to pressure force, then the runaway collapse of the cloud stops to form a rotationally-supported disk. After the accretion of the envelope, the disk undergoes a ring instability, eventually fragmenting into a binary. If the initial angular momentum is relatively large, a bar-type instability arises, resulting in the collapse into a single star through rapid angular momentum transfer. The present results show that a significant fraction of Pop III stars are expected to form in binary systems, even if they are quite massive or less massive. The cosmological implications of Population III binaries are briefly discussed.

The collapse of slowly rotating molecular cloud cores threaded by magnetic fields is investigated by high-resolution numerical simulation. Outflow formation in the collapsing cloud cores is also followed. In the models examined, the cloud core and parent cloud rotate rigidly and are initially threaded by a uniform magnetic field. The simulations show that the cloud core collapses along the magnetic field lines. The magnetic field in the dense region of the cloud core rotates faster than that of the parent cloud as a consequence of spin-up of the central region during the collapse. The cloud core exhibits significant precession of the rotation axis, magnetic field, and disk orientation, with precession highest in the models with low initial field strength ($\lesssim 20 \mu {\rm G}$). Precession in models with initial fields of $\sim 40 \mu {\rm G}$ is suppressed by strong magnetic braking. Magnetic braking transfers angular momentum form the central region and acts more strongly on the component of angular momentum oriented perpendicular to the magnetic field. After the formation of an adiabatic core, outflow is ejected along the local magnetic field lines. Strong magnetic braking associated with the outflow causes the direction of angular momentum to converge with that of the local magnetic field, resulting in the convergence of the local magnetic field, angular momentum, outflow, and disk orientation by the outflow formation phase. The magnetic field of a young star is inclined at an angle of no more than $30^\circ$ from that of the parent cloud at initial field strengths of $\sim 20 \mu {\rm G}$, while at an initial field strength of $\sim 40 \mu {\rm G}$, the magnetic field of the young star is well aligned with that of the parent cloud.

Based on high-resolution two-dimensional hydrodynamic simulations, we show that the bulk gas motions in a cluster of galaxies, which are naturally expected during the process of hierarchical structure formation of the universe, have a serous impact on the core. We found that the bulk gas motions represented by acoustic-gravity waves create local but strong turbulence, which reproduces the complicated X-ray structures recently observed in cluster cores. Moreover, if the wave amplitude is large enough, they can suppress the radiative cooling of the cores. Contrary to the previous studies, the heating is operated by the turbulence, not weak shocks. The turbulence could be detected in near-future space X-ray missions such as ASTRO-E2.

(Abridged) Recent cosmic microwave background anisotropy results from the Wilkinson Microwave Anisotropy Probe suggest that the universe was reionized at a redshift around 20. Such an early reionization could arise through the ionizing radiation emitted by metal-free population III stars at redshifts of 10 and higher. We discuss infrared background (IRB) surface brightness spatial fluctuations from such a generation of early star formation. We show that the spatial clustering of these stars at tens of arcminute scales generates a contribution to the angular power spectrum of the IRB anisotropies at the same angular scales. This excess can be potentially detected when resolved foreground galaxies out to a redshift of a few is removed from the clustering analysis. Additionally, assuming a population III stellar spectrum, we suggest that the clustering excess related to the first generation of stars can be separated from brightness fluctuations resulting from other foreground sources and galaxies using multifrequency observations in the wavelength range of ~ 1 to 5 microns. In addition to removing the low-redshift galaxy population, the multifrequency data are essential to account for certain foreground contaminants such zodiacal light, which, if varying spatially over degree scales, can be a significant source of confusion for the proposed study. Using various instruments, we study the extent to which spatial fluctuations of the IRB can be studied in the near future.

This is the first paper about the fragmentation and mass outflow in the molecular cloud by using three-dimensional MHD nested-grid simulations. The binary star formation process is studied paying particular attention to the fragmentation of a rotating magnetized molecular cloud. We assume an isothermal rotating and magnetized cylindrical cloud in hydrostatic balance. Non-axisymmetric as well as axisymmetric perturbations are added to the initial state and the subsequent evolutions are studied. The evolution is characterized by three parameters: the amplitude of the non-axisymmetric perturbations, the rotation speed, and the magnetic field strength. As a result, it is found that non-axisymmetry hardly evolves in the early phase, but begins to grow after the gas contracts and forms a thin disk. Disk formation is strongly promoted by the rotation speed and the magnetic field strength. There are two types of fragmentation: fragmentation from a ring and that from a bar. Thin adiabatic cores fragment if a thickness is smaller than 1/4 of the radius. For the fragments to survive, they should be formed in a heavily elongated barred core or a flat round disk. In the models showing fragmentation, outflows from respective fragments are found as well as those driven by the rotating bar or the disk.

We present results of a deep mid-infrared survey in the SSA13 field with the Infrared Space Observatory (ISO). In order to probe the near-infrared light at high redshifts, we surveyed the field with the broad band LW2 (5-8.5um) filter of the mid-infrared camera ISOCAM. Adopting a highly redundant imaging strategy for the 23 hour observation and carefully treating gradual changes in the detector responsivity caused by a very high rate of cosmic ray impacts, we succeeded in reaching an 80% completeness limit of 16uJy in the central 7 arcmin2 region. Utilizing the signal-to-noise ratio map, we detected 65 sources down to 6uJy in the 16 arcmin2 field. Integral galaxy number counts at 6.7um are then derived, reaching 1.3e4/deg2 at the faint limit with a slope of -1.6 between 13uJy and 130uJy. Integrating individual sources in this flux range, the resolved fraction of the extragalactic background light at 6.7um is estimated to be 0.56nW/m2/sr. These results, which reach a flux limit three times fainter than those in the Hubble Deep Fields, are in fairly good agreement with a model prediction by Franceschini et al. (1997). Finally, we discuss the relation of distant massive E/S0 galaxies to the faint 6.7um galaxy population.

The fragmentation of molecular cloud cores a factor of 1.1 denser than the critical Bonnor-Ebert sphere is examined though three-dimensional numerical simulations. A nested grid is employed to resolve fine structure down to 1 AU while following the entire structure of the molecular cloud core of radius 0.14 pc. A total of 225 models are shown to survey the effects of initial rotation speed, rotation law, and amplitude of bar mode perturbation. The simulations show that the cloud fragments whenever the cloud rotates sufficiently slowly to allow collapse but fast enough to form a disk before first-core formation. The latter condition is equivalent to $\Omega_0 t_{\rm ff} \gtrsim 0.05$, where $\Omega_0$ and $t_{\rm ff}$ denote the initial central angular velocity and the freefall time measured from the central density, respectively. Fragmentation is classified into six types: disk-bar, ring-bar, satellite, bar, ring, and dumbbell types according to the morphology of collapse and fragmentation. When the outward decrease in initial angular velocity is more steep, the cloud deforms from spherical at an early stage. The cloud deforms into a ring only when the bar mode m = 2 perturbation is very minor. The ring fragments into two or three fragments via ring-bar type fragmentation and into at least three fragments via ring type fragmentation. When the bar mode is significant, the cloud fragments into two fragments via either bar or dumbbell type fragmentation. These fragments eventually merge due to their low angular momenta, after which several new fragments form around the merged fragment via satellite type fragmentation.

We present a Monte-Carlo technique to probe the cosmic Infrared background radiation (CIB) with faint source confusion in deep, wide-field 4 micron images. We report a detection of surface brightness fluctuations using an ensemble of 5000 pixels. The brightness fluctuations do not show correlated structure and appear to be dominated by galaxies. We estimate a lower limit (95% C.L.) on CIB by modeling the number-count behavior of contributing sources to the detected fluctuations.

We review the next generation Japanese infrared space mission, ASTRO-F. ASTRO-F will be the first survey of the entire sky at infrared wavelengths since the IRAS mission almost 20 years ago. ASTRO-F will survey the entire sky in 4 far-infrared bands from 50-200microns and 2 mid-infrared bands at 9 and 20microns to sensitivities of 10-1000 times deeper than the IRAS satellite at angular resolutions of 25-45arcsec (c.f. IRAS 2-5arcmins). ASTRO-F can be considered a SUPER-IRAS. Using the galaxy evolution model of Pearson (2001) we produce expected numbers of sources under 3 different cosmological world models. We predict that ASTRO-F will detect of the order of 10's millions of sources in the far-infrared wavelength bands, most of which will be dusty LIG/ULIGs of which as many as half will lie at redshifts greater than unity. We produce number-redshift distributions, flux-redshift and colour-colour diagrams for the survey and discuss various segregation and photometric redshift techniques. Furthermore, we investigate the large scale structure scales that will be accessed by ASTRO-F, discovering that ASTRO-F and SIRTF-SWIRE probe both different scales and redshift domains and concluding that the 2 missions will supplement rather than supplant one another.

We present a numerical method for solving the Poisson equation on a nested grid. The nested grid consists of uniform grids having different grid spacing and is designed to cover the space closer to the center with a finer grid. Thus our numerical method is suitable for computing the gravity of a centrally condensed object. It consists of two parts: the difference scheme for the Poisson equation on the nested grid and the multi-grid iteration algorithm. It has three advantages: accuracy, fast convergence, and scalability. First it computes the gravitational potential of a close binary accurately up to the quadraple moment, even when the binary is resolved only in the fine grids. Second residual decreases by a factor of 300 or more by each iteration. We confirmed experimentally that the iteration converges always to the exact solution of the difference equation. Third the computation load of the iteration is proportional to the total number of the cells in the nested grid. Thus our method gives a good solution at the minimum expense when the nested grid is large. The difference scheme is applicable also to the adaptive mesh refinement in which cells of different sizes are used to cover a domain of computation.

We report on optical and near-infrared (NIR) follow-up spectroscopy of faint far-infrared (FIR) sources found in our deep FIR survey by Kawara et al.

We present basic observational strategies for ASTRO-F (also known as the Imaging Infra Red Surveyor (IRIS)) to be launched in 2004 by the Japanese Institute of Space and Astronautical Science (ISAS). We examine 2 survey scenarios, a deep ~1sq.deg. survey reaching sensitivities an order of magnitude below than the deepest surveys performed by ISO in the mid-IR, and a shallow \~18sq.deg mid-infrared (7-25um in 6 bands) covering an area greater than the entire area covered by all ISO mid-IR surveys. Using 2 cosmological models the number of galaxies predicted for each survey is calculated. The first model uses an enhancement of the pure luminosity evolution model of Pearson & Rowan-Robinson while the new models incorporate a strongly evolving ULIG component. For the deep survey, between 20,000-30,000 galaxies should be detected in the shortest wavebands and ~5000 in the longest (25um) band. The shallow survey would be expected to detect of the order of 100,000 - 150,000 sources.

We investigate stability of Toomre-Hayashi model for a self-gravitating, rotating gas disk. The rotation velocity, $ v_\phi $, and sound speed, $ c _{\rm s} $, are spatially constant in the model. We show that the model is unstable against an axisymmetric perturbation irrespectively of the values of $ v_\phi $ and $ c_{\rm s} $. When the ratio, $ v_\phi / c_{\rm s} $, is large, the model disk is geometrically thin and unstable against an axisymmetric perturbation having a short radial wavelength, i.e., unstable against fragmentation. When $ v_\phi / c_{\rm s} $ is smaller than 1.20, it is unstable against total collapse. Toomre-Hayashi model of $ 1.7 \la v_\phi / c_{\rm s} \la 3 $ was thought to be stable against axisymmetric perturbations in earlier studies in which only radial motion was taken into account. Thus, the instability of Toomre-Hayashi model having a medium $ v_\phi / c_{\rm s} $ is due to not change in the surface density but that in the disk height. We also find that the singular isothermal perturbation is stable against non-spherical peturbations.

We present the results of ground-based imaging spectroscopy of the [Ne II] 12.8 micron line emitted from the ultracompact (UC) H II regions; W51d, G45.12+0.13, G35.20-1.74 and Monoceros R2, with 2arcsec spatial resolution. We found that the overall distribution of the [Ne II] emission is generally in good agreement with the radio (5 or 15 GHz) VLA distribution for each source. The Ne+ abundance ([Ne+/H+]) distributions are also derived from the [Ne II] and the radio maps. As for G45.12+0.13 and W51d, the Ne+ abundance decreases steeply from the outer part of the map toward the radio peak. On the other hand, the Ne+ abundance distributions of G35.20-1.74 and Mon R2 appear rather uniform. These results can be interpreted by the variation of ionizing structures of neon, which is primarily determined by the spectral type of the ionizing stars. We have evaluated the effective temperature of the ionizing star by comparing the Ne+ abundance averaged over the whole observed region with that calculated by H II region models based on recent non-LTE stellar atmosphere models: 39,100 (+1100, -500) K (O7.5V-O8V) for W51d, 37,200 (+1000, -700) K (O8V-O8.5V) for G45.12+0.13, 35,000-37,600 (+1500, -600) K (O8V-O9V) for G35.20-1.74, and < 34,000 K (< B0V) for Mon R2. These effective temperatures are consistent with those inferred from the observed Ne+ abundance distributions.

We investigate the characteristics of FIR brightness fluctuations at 90 micron and 170 micron in the Lockman Hole, which were surveyed with ISOPHOT aboard the Infrared Space Observatory(ISO). We first calculated the angular correlation function of each field and then its Fourier transform(the angular Power Spectral Density: PSD) over the spatial frequency range of f=0.05 - 1 arcmin^-1. The PSDs are found to be rather flat at low spatial frequencies (f <= 0.1 arcmin^-1), slowly decreasing toward higher frequencies. These spectra are unlike the power-law ones seen in the IR cirrus fluctuations, and are well explained by randomly distributed point sources. Furthermore, point-to-point comparison between 90 micron and 170 micron brightness shows a linear correlation between them, and the slope of the linear fit is much shallower than that expected from the IR cirrus color, and is consistent with the color of galaxies at low or moderate redshift(z < 1). We conclude that the brightness fluctuations in the Lockman Hole are not caused by the IR cirrus, but are most likely due to faint star-forming galaxies. We also give the constraints on the galaxy number counts down to 35 mJy at 90 micron and 60 mJy at 170 micron, which indicate the existence of a strong evolution down to these fluxes in the counts. The galaxies responsible for the fluctuations also significantly contribute to the cosmic infrared background radiation.

We present a search for near-infrared (3.5-5 micron) emission from baryonic dark matter in the form of low-mass stars and/or brown dwarfs in the halo of the nearby edge-on spiral galaxy NGC 5907. The observations were made using a 256 by 256 InSb array with a pixel scale of 17" at the focus of a liquid-helium-cooled telescope carried above the Earth's atmosphere by a sounding rocket. In contrast to previous experiments which have detected a halo around NGC 5907 in the V, R, I, J and K bands at galactic radii 6kpc < r < 10kpc, our search finds no evidence for emission from a halo at 10kpc < r < 30kpc. Assuming a halo mass density scaling as r^(-2), which is consistent with the flat rotation curves that are observed out to radii of 32kpc, the lower limit of the mass-to-light ratio at 3.5-5 microns for the halo of NGC 5907 is 250 (2 sigma) in solar units. This is comparable to the lower limit we have found previously for NGC 4565 (Uemizu et al. 1998). Based on recent models, our non-detection implies that hydrogen- burning stars contribute < 15% of the mass of the dark halo of NGC 5907. Our results are consistent with the previous detection of extended emission at r < 10kpc if the latter is caused by a stellar population that has been ejected from the disk because of tidal interactions. We conclude that the dark halo of NGC 5907, which is evident from rotation curves that extend far beyond 10kpc, is not comprised of hydrogen burning stars.

Results of a narrow-band imaging for the redshifted [OII]3727 emission around a quasar at z=4.7 obtained with the Subaru telescope and CISCO (a Cassegrain near infrared camera) are presented. A significant emission line is detected in the narrow-band H_2 (v=1-0 S(1)) filter at a location 2.4" northwest from the quasar, where the presence of a companion has been reported in Lyman alpha emission and the rest-frame UV continuum. We identify this line as [OII]3727 emission and confirm that the source really is a companion at z=4.7. The [OII]3727 flux from the companion is estimated to be 2.5 x 10^{-17} erg s^{-1} cm^{-2}. If the companion is a star forming object, the inferred star formation rate is as high as 45-230 M_{solar} yr^{-1} even without assuming the extinction correction. This value is higher than those derived from the Lyman alpha emission or from the UV continuum. Thus, provided that the difference is caused by dust extinction, the extinction corrected star formation rate is calculated to be 45 to 2300 M_{solar} yr^{-1} depending on the assuming extinction curves.

We investigate stability of a gravitationally collapsing iron core against non-spherical perturbation. The gravitationally collapsing iron core is approximated by a similarity solution for dynamically collapsing polytropic gas sphere. We find that the similarity solution is unstable against non-spherical perturbations. The perturbation grows in proportion to $ (t - t_0) ^{-\sigma} $ while the the central density increases in proportion to $ (t - t_0) ^{-2} $. The growth rate is $ \sigma = 1/3 + \ell (\gamma - 4/3) $, where $ \gamma $ and $ \ell $ denote the polytropic index and the parameter $ \ell $ of the spherical harmonics, $ Y_{\ell} ^m (\theta, \phi) $, respectively. The growing perturbation is dominated by vortex motion. Thus it excites global convection during the collapse and may contribute to material mixing in a type II supernova.

We discuss stability of dynamically collapsing gas spheres. We use a similarity solution for a dynamically collapsing sphere as the unperturbed state. In the similarity solution the gas pressure is approximated by a polytrope of $ P = K \rho ^\gamma $. We examine three types of perturbations: bar ($ \ell = 2$) mode, spin-up mode, and Ori-Piran mode. When $ \gamma < 1.097 $, it is unstable against bar-mode. It is unstable against spin-up mode for any $ \gamma $. When $ \gamma < 0.961 $, the similarity solution is unstable against Ori-Piran mode. The unstable mode grows in proportion to $ | t - t_0 | ^{-\sigma} $ while the central density increases in proportion to $ \rho_c \propto (t - t_0) ^{-2} $ in the similarity solution. The growth rate, $ \sigma $ is obtained numerically as a function of $ \gamma $ for bar mode and Ori-Piran mode. The growth rate of the bar mode is larger for a smaller $ \gamma $. The spin-up mode has the growth rate of $ \sigma = 1/3 $ for any $ \gamma $.

We report the results of a search for neutrino-induced particle cascades using a deep ocean water Cherenkov detector. The effective mass of the detector, a string of seven 40 cm diameter photomultipliers at 5.2 m spacing, is found through simulation analysis to be surprisingly large: greater than 1 megaton of water at incident neutrino energies of 1 PeV. We find no evidence for neutrino-induced cascades in 18.6 hours of observation. Although the limit implied by this observation is the strongest yet for predictions of active galatic nuclei (AGN) neutrinos at energies above 100 TeV, perhaps the more intriguing result is that the power of these techniques can be exploited to test these AGN models in a relatively short time.

The dynamical instability of a self-gravitating magnetized filamentary cloud was investigated while taking account of rotation around its axis. The filamentary cloud of our model is supported against self-gravity in part by both a magnetic field and rotation. The density distribution in equilibrium was assumed to be a function of the radial distance from the axis, $ \rho _0 (r) \, = \, \rho _{\rm c} \, ( 1 \, + \, r ^2 / 8 H ^2 ) ^{-2} $, where $ \rho _{\rm c} $ and $ H $ are model parameters specifying the density on the axis and the length scale, respectively; the magnetic filed was assumed to have both longitudinal ($z$-) and azimuthal ($\varphi$-) components with a strength of $ B _0 (r) \, \propto \, \sqrt{ \rho _0 (r) } $. The rotation velocity was assumed to be $ v _{0\varphi} \, = \, \Omega _{\rm c} \, r \, (1 \, + \, r ^2 / 8 H ^2 ) ^{-1/2} $. We obtained the growth rate and eigenfunction numerically for (1) axisymmetric $ ( m \, = \, 0 ) $ perturbations imposed on a rotating cloud with a longitudinal magnetic field, (2) non-axisymmetric $ ( m \, = \, 1 ) $ perturbations imposed on a rotating cloud with a longitudinal magnetic filed, and (3) axisymmetric perturbations imposed on a rotating cloud with a helical magnetic field. The fastest growing perturbation is an axisymmetric one for all of the model clouds studied. Its wavelength is $ \lambda _{\rm max} \, \le \, 11.14 \, H $ for a non-rotating cloud without a magnetic field, and is shorter