We herein report a theoretical study of the geometrical structure of porous dust aggregates formed via ballistic cluster-cluster aggregation (BCCA). We calculated the gyration radius $R_{\rm gyr}$ and the graph-based geodesic radius $R_{\rm geo}$ as a function of the number of constituent particles $N$. We found that $R_{\rm gyr} / r_{0} \sim N^{0.531 \pm 0.011}$ and $R_{\rm geo} / r_{0} \sim N^{0.710 \pm 0.013}$, where $r_{0}$ is the radius of constituent particles. Furthermore, we defined two constants that characterize the geometrical structure of fractal aggregates: $D_{\rm f}$ and $\alpha$. The definition of $D_{\rm f}$ and $\alpha$ are $N \sim {( R_{\rm gyr} / r_{0} )}^{D_{\rm f}}$ and ${R_{\rm geo}} / {r_{0}} \sim {\left( {R_{\rm gyr}} / {r_{0}} \right)}^{\alpha}$, respectively. Our study revealed that $D_{\rm f} \simeq 1.88$ and $\alpha \simeq 1.34$ for the clusters of the BCCA. In addition, we also studied the filling factor dependence of thermal conductivity of statically compressed fractal aggregates. From this study, we reveal that the thermal conductivity of statically compressed aggregates $k$ is given by $k \sim 2 k_{\rm mat} {( r_{\rm c} / r_{0} )} \phi^{(1 + \alpha) / (3 - D_{\rm f})}$, where $k_{\rm mat}$ is the material thermal conductivity, $r_{\rm c}$ is the contact radius of constituent particles, and $\phi$ is the filling factor of dust aggregates.

The poor stickiness of silicate dust grains is a major obstacle to the formation of rocky planetesimals. In this study, we examine the possibility that silicate grains with an organic mantle, which we call Organic-Mantled Grains (OMGs), form planetesimals through direct coagulation. Organic mantles are commonly found in interplanetary dust particles, and laboratory experiments show that they are softer than silicates, in particular in warm environments. This, combined with the theory of particle adhesion, implies that OMGs are stickier than bare silicate grains. Because organic mantles can survive up to 400 K, silicate grains inside the water snow line in protoplanetary disks can in principle hold such mantles. We construct a simple grain adhesion model to estimate the threshold collision velocity below which aggregates of OMGs can grow. The model shows that aggregates of 0.1 \micron-sized OMGs can overcome the fragmentation barrier in protoplanetary disks if the mantles are as thick as those in interplanetary dust particles and if the temperature is above $\sim$ 200 K. We use this adhesion model to simulate the global evolution of OMG aggregates in the inner part of a protoplanetary disk, demonstrating that OMG aggregates can indeed grow into planetesimals under favorable conditions. Because organic matter is unstable at excessively high temperatures, rocky planetesimal formation by the direct sticking of OMGs is expected to occur in a disk annulus corresponding to the temperature range $ \sim$ 200--400 K. The organic-rich planetesimals may grow into carbon-poor rocky planetesimals by accreting a large amount of carbon-poor chondrules.

Shock-wave heating within the solar nebula is one of the leading candidates for the source of chondrule-forming events. Here, we examine the possibility of compound chondrule formation via optically thin shock waves. Several features of compound chondrules indicate that compound chondrules are formed via the collisions of supercooled precursors. We evaluate whether compound chondrules can be formed via the collision of supercooled chondrule precursors in the framework of the shock-wave heating model by using semi-analytical methods and discuss whether most of the crystallized chondrules can avoid destruction upon collision in the post-shock region. We find that chondrule precursors immediately turn into supercooled droplets when the shock waves are optically thin and they can maintain supercooling until the condensation of evaporated fine dust grains. Owing to the large viscosity of supercooled melts, supercooled chondrule precursors can survive high-speed collisions on the order of $1\ {\rm km}\ {\rm s}^{-1}$ when the temperature is below $\sim 1400\ {\rm K}$. From the perspective of the survivability of crystallized chondrules, shock waves with a spatial scale of $\sim 10^{4}\ {\rm km}$ may be potent candidates for the chondrule formation mechanism. Based on our results from one-dimensional calculations, a fraction of compound chondrules can be reproduced when the chondrule-to-gas mass ratio in the pre-shock region is $\sim 2 \times 10^{-3}$, which is approximately half of the solar metallicity.

Understanding the heat transfer mechanism within dust aggregates is of great importance for many subjects in planetary science. We calculated the coordination number and the thermal conductivity through the solid network of compressed dust aggregates. We found a simple relationship between the coordination number and the filling factor and revealed that the thermal conductivity through the solid network of aggregates is represented by a power-law function of the filling factor and the coordination number.

We examined the excitations of planetesimals caused by the resonances of a giant planet in a protoplanetary gas disk. The highly excited planetesimals generate bow shocks, the mechanism of which results in chondrule formation, crystallization of silicate dust, and evaporation of icy planetesimals. The planetesimals beyond 2:1 resonance migrate owing to the gas drag and obtain the maximum eccentricity around 3:1 resonance, which is located at approximately half the planetary distance. The eccentricity depends on the parameters of the planetesimals and the Jovian planet, such as size and location, and gas density of the disk. The maximum relative velocity of a 100-km-sized planetesimal with respect to the gas disk reaches up to ~12 km/s in the case of Jupiter owing to secular resonance, which occurs because of the disk's gravity. We find that if a Jovian mass planet is located within 10 au, the planetesimals larger than 100 km gain sufficient velocity to cause the melting of chondrule precursors and crystallization of the silicate. The maximum velocity is higher for large planetesimals and eccentric planets. Planetesimals are trapped temporarily in the resonances and continue to have high speed over >1 Myr after the formation of a Jovian planet. This duration fits into the timescale of chondrule formation suggested by the isotopic data. The evaporation of icy planetesimals occurs when a Jovian planet is located within 15 au. This mechanism can be a new indicator of planet formation in exosystems if some molecules ejected from icy planetesimals are detected.

We report the results of detailed investigation of the Earth's transmission spectra during the lunar eclipse on UT 2011 December 10. The spectra were taken by using the High Dispersion Spectrograph (HDS) mounted on the Subaru 8.2 m telescope with unprecedented resolutions both in time and wavelength (300 s exposure time in umbra and 160,000 spectral resolution, respectively). In our penumbra and umbra data, we detected the individual absorption lines of O2 and H2O in transmission spectra and found that it became deeper as the eclipse became deeper. It indicates that the sunlight reaching the Moon passed through lower layers in the Earth's atmosphere with time because we monitored a given point on the Moon during the full eclipse duration. From the comparison between the observed and theoretically constructed transmission spectra, the lowest altitude at which the sunlight actually passed through the atmosphere is estimated to be about 10 km from the ground, which suggests the existence of sunlight blocking clouds below that altitude. Our result would be a test case for future investigations of atmospheric structure of Earth-like exoplanets via transmission spectroscopy including the refraction effect of the planetary atmosphere.

Highly porous dust aggregates can break through the radial drift barrier, but previous studies assumed disks in their later stage, where the disks have a very small mass and low temperature. In contrast, dust coagulation should begin in the very early stage such as the disk formation stage because the growth timescale of dust is shorter than the disk formation timescale if there is no process to suppress the collisional growth of dust. We investigate the possibility of planetesimal formation via direct collisional growth in the very early stage of a protoplanetary disk. We show that, in the very early stage of protoplanetary disks, icy dust aggregates suffer radial drift and deplete without forming planetesimal-sized objects. This is because as the disk temperature easily increases by the viscous heating in the disk formation stage, the area where the dust can break through the radial drift barrier is restricted only to the inside the snowline. This suggests that in the disk formation stage, icy planetesimal formation via direct collisional growth of dust is difficult.

Infrared Imaging Spectrograph (IRIS) is the first light instrument for the Thirty Meter Telescope (TMT) that consists of a near-infrared (0.84 to 2.4 micron) imager and integral field spectrograph (IFS) which operates at the diffraction-limit utilizing the Narrow-Field Infrared Adaptive Optics System (NFIRAOS). The imager will have a 34 arcsec x 34 arcsec field of view with 4 milliarcsecond (mas) pixels. The IFS consists of a lenslet array and slicer, enabling four plate scales from 4 mas to 50 mas, multiple gratings and filters, which in turn will operate hundreds of individual modes. IRIS, operating in concert with NFIRAOS will pose many challenges for the data reduction system (DRS). Here we present the updated design of the real-time and post-processing DRS. The DRS will support two modes of operation of IRIS: (1) writing the raw readouts sent from the detectors and performing the sampling on all of the readouts for a given exposure to create a raw science frame; and (2) reduction of data from the imager, lenslet array and slicer IFS. IRIS is planning to save the raw readouts for a given exposure to enable sophisticated processing capabilities to the end users, such as the ability to remove individual poor seeing readouts to improve signal-to-noise, or from advanced knowledge of the point spread function (PSF). The readout processor (ROP) is a key part of the IRIS DRS design for writing and sampling of the raw readouts into a raw science frame, which will be passed to the TMT data archive. We discuss the use of sub-arrays on the imager detectors for saturation/persistence mitigation, on-detector guide windows, and fast readout science cases (< 1 second).

The InfraRed Imaging Spectrograph (IRIS) is the first-light client instrument for the Narrow Field Infrared Adaptive Optics System (NFIRAOS) on the Thirty Meter Telescope (TMT). IRIS includes three natural guide star (NGS) On-Instrument Wavefront Sensors (OIWFS) to measure tip/tilt and focus errors in the instrument focal plane. NFIRAOS also has an internal natural guide star wavefront sensor, and IRIS and NFIRAOS must precisely coordinate the motions of their wavefront sensor positioners to track the locations of NGSs while the telescope is dithering (offsetting the telescope to cover more area), to avoid a costly re-acquisition time penalty. First, we present an overview of the sequencing strategy for all of the involved subsystems. We then predict the motion of the telescope during dithers based on finite-element models provided by TMT, and finally analyze latency and jitter issues affecting the propagation of position demands from the telescope control system to individual motor controllers.

$\mathit{Context.}$ The thermal conductivity of highly porous dust aggregates is a key parameter for many subjects in planetary science; however, it is not yet fully understood. $\mathit{Aims.}$ In this study, we investigate the thermal conductivity of fluffy dust aggregates with filling factors of less than $10^{-1}$. $\mathit{Methods.}$ We determine the temperature structure and heat flux of the porous dust aggregates calculated by $N$-body simulations of static compression in the periodic boundary condition. $\mathit{Results.}$ We derive an empirical formula for the thermal conductivity through the solid network $k_{\rm sol}$ as a function of the filling factor of dust aggregates $\phi$. The results reveal that $k_{\rm sol}$ is approximately proportional to ${\phi}^{2}$, and the thermal conductivity through the solid network is significantly lower than previously assumed. In light of these findings, we must reconsider the thermal histories of small planetary bodies.

Several pieces of evidence suggest that silicate grains in primitive meteorites are not interstellar grains but condensates formed in the early solar system. Moreover, the size distribution of matrix grains in chondrites implies that these condensates might be formed as nanometer-sized grains. Therefore, we propose a novel scenario for rocky planetesimal formation in which nanometer-sized silicate grains are produced by evaporation and recondensation events in early solar nebula, and rocky planetesimals are formed via aggregation of these nanograins. We reveal that silicate nanograins can grow into rocky planetesimals via direct aggregation without catastrophic fragmentation and serious radial drift, and our results provide a suitable condition for protoplanet formation in our solar system.

We present a novel model showing that compound chondrules are formed by collisions of supercooled droplets. This model reproduces two prominent observed features of compound chondrules: the nonporphyritic texture and the size ratio between two components.

We made near infrared multicolor imaging observations of a disk around Herbig Be star HD100546 using Gemini/NICI. K (2.2\,$\mu$m), H$_2$O ice (3.06\,$\mu$m), and L'(3.8\,$\mu$m) disk images were obtained and we found the 3.1\,$\mu$m absorption feature in the scattered light spectrum, likely due to water ice grains at the disk surface. We compared the observed depth of the ice absorption feature with the disk model based on \cite{Oka2012} including water ice photodesorption effect by stellar UV photons. The observed absorption depth can be explained by the both disk models with/without photodesorption effect within the measurement accuracy, but slightly favors the model with photodesorption effects, implying that the UV photons play an important role on the survival/destruction of ice grains at the Herbig Ae/Be disk surface. Further improvement on the accuracy of the observations of the water ice absorption depth is needed to constrain the disk models.

Motivated by recent detection of transiting high-density super-Earths, we explore the detectability of hot rocky super-Earths orbiting very close to their host stars. In the environment hot enough for their rocky surfaces to be molten, they would have the atmosphere composed of gas species from the magma oceans. In this study, we investigate the radiative properties of the atmosphere that is in the gas/melt equilibrium with the underlying magma ocean. Our equilibrium calculations yield Na, K, Fe, Si, SiO, O, and O$_2$ as the major atmospheric species. We compile the radiative-absorption line data of those species available in literature, and calculate their absorption opacities in the wavelength region of 0.1--100~$\mathrm{\mu m}$. Using them, we integrate the thermal structure of the atmosphere. Then, we find that thermal inversion occurs in the atmosphere because of the UV absorption by SiO. In addition, we calculate the ratio of the planetary to stellar emission fluxes during secondary eclipse, and find prominent emission features induced by SiO at 4~$\mathrm{\mu m}$ detectable by Spitzer, and those at 10 and 100~$\mathrm{\mu m}$ detectable by near-future space telescopes.

It is proposed that planetesimals perturbed by Jovian mean-motion resonances are the source of shock waves that form chondrules. It is considered that this shock-induced chondrule formation requires the velocity of the planetesimal relative to the gas disk to be on the order of > 7 km/s at 1 AU. In previous studies on planetesimal excitation, the effects of Jovian mean-motion resonance together with the gas drag were investigated, but the velocities obtained were at most 8 km/s in the asteroid belt, which is insufficient to account for the ubiquitous existence of chondrules. In this paper, we reexamine the effect of Jovian resonances and take into account the secular resonance in the asteroid belt caused by the gravity of the gas disk. We find that the velocities relative to the gas disk of planetesimals a few hundred kilometers in size exceed 12 km/s, and that this is achieved around the 3:1 mean-motion resonance. The heating region is restricted to a relatively narrow band between 1.5 AU and 3.5 AU. Our results suggest that chondrules were produced effectively in the asteroid region after Jovian formation. We also find that many planetesimals are scattered far beyond Neptune. Our findings can explain the presence of crystalline silicate in comets if the scattered planetesimals include silicate dust processed by shock heating.

We have discovered that Europa, Ganymede and Callisto are bright around 1.5 {\mu}m even when not directly lit by sunlight, based on observations from the Hubble Space Telescope and the Subaru Telescope. The observations were conducted with non-sidereal tracking on Jupiter outside of the field of view to reduce the stray light subtraction uncertainty due to the close proximity of Jupiter. Their eclipsed luminosity was $10^{-6}$-$10^{-7}$ of their uneclipsed brightness, which is low enough that this phenomenon has been undiscovered until now. In addition, Europa in eclipse was <1/10 of the others at 1.5 {\mu}m, a potential clue to the origin of the source of luminosity. Likewise, Ganymede observations were attempted at 3.6 {\mu}m by the Spitzer Space Telescope but it was not detected, suggesting a significant wavelength dependence. The reason why they are luminous even when in the Jovian shadow is still unknown, but forward-scattered sunlight by haze in the Jovian upper atmosphere is proposed as the most plausible candidate. If this is the case, observations of these Galilean satellites while eclipsed by the Jovian shadow provide us a new technique to investigate Jovian atmospheric composition, and investigating the transmission spectrum of Jupiter by this method is important for investigating the atmosphere of extrasolar giant planets by transit spectroscopy.

During their evolution, short-period exoplanets may lose envelope mass through atmospheric escape owing to intense XUV (X-ray and extreme ultraviolet) radiation from their host stars. Roche-lobe overflow induced by orbital evolution or intense atmospheric escape can also contribute to mass loss. To study the effects of mass loss on inner planet populations, we calculate the evolution of hot Jupiters considering mass loss of their envelopes and thermal contraction. Mass loss is assumed to occur through XUV-driven atmospheric escape and the following Roche-lobe overflow. The runaway effect of mass loss results in a dichotomy of populations: hot Jupiters that retain their envelopes and super Earths whose envelopes are completely lost. Evolution primarily depends on the core masses of planets and only slightly on migration history. In hot Jupiters with small cores (~10 Earth masses), runaway atmospheric escape followed by Roche-lobe overflow may create sub-Jupiter deserts, as observed in both mass and radius distributions of planetary populations. Comparing our results with formation scenarios and observed exoplanets populations, we propose that populations of closely orbiting exoplanets are formed by capturing of planets at/inside the inner edges of protoplanetary disks and subsequent evaporation of sub-Jupiters.

Photoevaporation by stellar ionizing radiation is believed to play an important role in the dispersal of disks around young stars. The mass loss model for dust-free disks developed by Hollenbach et al. is currently regarded as a conventional one and has been used in a wide variety of studies. However, the rate in this model was derived by the crude so-called 1+1D approximation of ionizing radiation transfer, which assumes that diffuse radiation propagates in a direction vertical to the disk. In this study, we revisit the photoevaporation of dust-free disks by solving the 2D axisymmetric radiative transfer for steady-state disks. Unlike that solved by the conventional model, we determine that direct stellar radiation is more important than the diffuse field at the disk surface. The radial density distribution at the ionization boundary is represented by the single power-law with an index -3/2 in contrast to the conventional double power-law. For this distribution, the photoevaporation rate from the entire disk can be written as a function of the ionizing photon emissivity, Phi_EUV, from the central star and the disk outer radius, r_d, as follows: Mdot_PE = 5.4 x 10^-5 x (Phi_EUV/10^49 sec^-1)^1/2 x (r_d/1000 AU)^1/2 Msun/yr. This new rate depends on the outer disk radius rather than on the gravitational radius as in the conventional model, caused by the enhanced contribution to the mass loss from the outer disk annuli. In addition, we discuss its applications to present-day as well as primordial star formation.

We investigate the effect of photodesorption on the snow line position at the surface of a protoplanetary disk around a Herbig Ae/Be star, motivated by the detection of water ice particles at the surface of the disk around HD142527 by Honda et al. For this aim, we obtain the density and temperature structure in the disk with a 1+1D radiative transfer and determine the distribution of water ice particles in the disk by the balance between condensation, sublimation, and photodesorption. We find that photodesorption induced by the far-ultraviolet radiation from the central star depresses the ice-condensation front toward the mid-plane and pushes the surface snow line outward significantly when the stellar effective temperature exceeds a certain critical value. This critical effective temperature depends on the stellar luminosity and mass, the water abundance in the disk, and the yield of photodesorption. We present an approximate analytic formula for the critical temperature. We separate Herbig Ae/Be stars into two groups on the HR diagram according to the critical temperature; one is the disks where photodesorption is effective and from which we may not find ice particles at the surface, and the other is the disks where photodesorption is not effective. We estimate the snow line position at the surface of the disk around HD142527 to be 100--300 AU which is consistent with the water ice detection at >140 AU in the disk. All results depend on the dust grain size by a complex way and this point requires more work in future.

Evolution of a snow line in an optically-thick protoplanetary disk is investigated with numerical simulations. The ice-condensing region in the disk is obtained by calculating the temperature and the density with the 1+1D approach. The snow line migrates as the mass accretion rate (\dot{M}) in the disk decreases with time. Calculations are carried out from an early phase with high disk accretion rates (\dot{M} \sim 10^{-7} M_sun/yr) to a later phase with low disk accretion rates (\dot{M} \sim 10^{-12} M_sun/yr) using the same numerical method. It is found that the snow line moves inward for \dot{M} > 10^{-10} M_sun/yr, while it gradually moves outward in the later evolution phase with \dot{M} < 10^{-10} M_sun/yr. In addition to the silicate opacity, the ice opacity is taken into consideration. In the inward migration phase, the additional ice opacity increases the distance of the snow line from the central star by a factor of 1.3 for dust grains < 10 micro-meter in size and 1.6 for > 100 micro-meter. It is inevitable that the snow line comes inside the Earth's orbit in the course of the disk evolution, if the alpha viscosity parameter is in a range 0.001-0.1, the dust-to-gas mass ratio is higher than a tenth of the solar abundance value, and the dust grains are smaller than 1 mm. The formation of water-devoid planetesimals in the terrestrial planet region seems to be difficult throughout the disk evolution, which imposes a new challenge to planet formation theory.

Several lines of evidence suggest that fine silicate crystals observed in primitive meteorite and interplanetary dust particles (IDPs) nucleated in a supersaturated silicate vapor followed by crystalline growth. We investigated evaporation of $\mu$m-sized silicate particles heated by a bow shock produced by a planetesimal orbiting in the gas in the early solar nebula and condensation of crystalline silicate from the vapor thus produced. Our numerical simulation of shock-wave heating showed that these {\mu}m-sized particles evaporated almost completely when the bow shock is strong enough to cause melting of chondrule precursor dust particles. We found that the silicate vapor cools very rapidly with expansion into the ambient unshocked nebular region; the cooling rate is estimated, for instance, to be as high as 2000 K s$^{-1}$ for a vapor heated by a bow shock associated with a planetesimal of radius 1 km. The rapid cooling of the vapor leads to nonequilibrium gas-phase condensation of dust at temperatures much lower than those expected from the equilibrium condensation. It was found that the condensation temperatures are lower by a few hundred K or more than the equilibrium temperatures. This explains the results of the recent experimental studies of condensation from a silicate vapor that condensation in such large supercooling reproduces morphologies similar to those of silicate crystals found in meteorites. Our results suggest strongly that the planetesimal bow shock is one of the plausible sites for formation of not only chondrules but also other cosmic crystals in the early solar system.

In order to assess the contribution of Lyman break galaxies (LBGs) and Lyman alpha emitters (LAEs) at redshifts 3<z<7 to the ionization of intergalactic medium (IGM), we investigate the escape fractions of ionizing photons from supernova-dominated primordial galaxies by solving the three-dimensional radiative transfer. The model galaxy is employed from an ultra-high-resolution chemodynamic simulation of a primordial galaxy by Mori & Umemura (2006),which well reproduces the observed properties of LAEs and LBGs. The total mass of model galaxy is 10^{11}M_sun. We solve not only photo-ionization but also collisional ionization by shocks. In addition, according to the chemical enrichment, we incorporate the effect of dust extinction, taking the size distributions of dust into account. As a result, we find that dust extinction reduces the escape fractions by a factor 1.5-8.5 in the LAE phase and by a factor 2.5-11 in the LBG phase, while the collisional ionization by shocks increases the escape fractions by a factor ~ 2. The resultant escape fractions are 0.07-0.47 in the LAE phase and 0.06-0.17 in the LBG phase. These results are well concordant with the recent estimations derived from the flux ratio at 1500 angstrom to 900 angstrom of LAEs and LBGs. Combining the resultant escape fractions with the luminosity functions of LAEs and LBGs, we find that high-z LAEs and LBGs can ionize the IGM at z=3-5. However, ionizing radiation from LAEs as well as LBGs falls short to ionize the IGM at z>6. That implies that additional ionization sources may required at z>6.

The temperature in the optically thick interior of protoplanetary discs is essential for the interpretation of millimeter observations of the discs, for the vertical structure of the discs, for models of the disc evolution and the planet formation, and for the chemistry in the discs. Since large icy grains have a large albedo even in the infrared, the effect of scattering of the diffuse radiation in the discs on the interior temperature should be examined. We have performed a series of numerical radiation transfer simulations including isotropic scattering by grains with various typical sizes for the diffuse radiation as well as for the incident stellar radiation. We also have developed an analytic model including isotropic scattering to understand the physics concealed in the numerical results. With the analytic model, we have shown that the standard two-layer approach is valid only for grey opacity (i.e. grain size $\ga10$ \micron) even without scattering. A three-layer interpretation is required for grain size $\la10$ \micron. When the grain size is 0.1--10 \micron, the numerical simulations show that isotropic scattering reduces the temperature of the disc interior. This reduction is nicely explained by the analytic three-layer model as a result of the energy loss by scatterings of the incident stellar radiation and of the warm diffuse radiation in the disc atmosphere. For grain size $\ga10$ \micron (i.e. grey scattering), the numerical simulations show that isotropic scattering does not affect the interior temperature. This is nicely explained by the analytic two-layer model; the energy loss by scattering in the disc atmosphere is exactly offset by the "green-house effect" due to scattering of the cold diffuse radiation in the interior.

We discuss how we obtain the spatial distribution of ice on the surface of the circumstellar disk around young stars. Ice in the disks plays a very important role in various issues, for instance, on the disk structure, on the planet formation, on the isotopic anomaly in meteorites, and on the origin of the sea on the Earth. Therefore, the spatially resolved observation of the condensation/sublimation front of ice, so-called ``snow line'' is strongly required. Here, we propose a new method for obtaining the spatially resolved ``snow line'' on the circumstellar disks by observing 3 \micron H$_2$O ice feature in the scattered light. Based on radiative transfer considerations, we show that the feature is clearly imprinted in the spectrum of the scattered light from both optically thick and thin circumstellar disks. We also show that the scattered light and the H$_2$O ice feature from protoplanetary disks are detectable and spatially resolvable with the current instruments through a $H_2O$ narrowband filter around 3 \micron. Finally, we present a diagnostics of disk dust properties on the $K-H_2O$ and $K-L'$ two color diagram.

We propose a new scenario for compound chondrule formation named as "fragment-collision model," in the framework of the shock-wave heating model. A molten cm-sized dust particle (parent) is disrupted in the high-velocity gas flow. The extracted fragments (ejectors) are scattered behind the parent and the mutual collisions between them will occur. We modeled the disruption event by analytic considerations in order to estimate the probability of the mutual collisions assuming that all ejectors have the same radius. We found that the estimated collision probability, which is the probability of collisions experienced by an ejector in one disruption event, can account for the observed fraction of compound chondrules. In addition, the model predictions are qualitatively consistent with other observational data (oxygen isotopic composition, textural types, and size ratios of constituents). Based on these results, we concluded that this new model can be one of the strongest candidates for the compound chondrule formation. It should be noted that all collisions do not necessarily lead to the compound chondrule formation. The formation efficiency and the future works which should be investigated in the forthcoming paper are also discussed.

Millimeter-sized, spherical silicate grains abundant in chondritic meteorites, which are called as chondrules, are considered to be a strong evidence of the melting event of the dust particles in the protoplanetary disk. One of the most plausible scenarios is that the chondrule precursor dust particles are heated and melt in the high-velocity gas flow (shock-wave heating model). We developed the non-linear, time-dependent, and three-dimensional hydrodynamic simulation code for analyzing the dynamics of molten droplets exposed to the gas flow. We confirmed that our simulation results showed a good agreement in a linear regime with the linear solution analytically derived by Sekiya et al. (2003). We found that the non-linear terms in the hydrodynamical equations neglected by Sekiya et al. (2003) can cause the cavitation by producing negative pressure in the droplets. We discussed that the fragmentation through the cavitation is a new mechanism to determine the upper limit of chondrule sizes. We also succeeded to reproduce the fragmentation of droplets when the gas ram pressure is stronger than the effect of the surface tension. Finally, we compared the deformation of droplets in the shock-wave heating with the measured data of chondrules and suggested the importance of other effects to deform droplets, for example, the rotation of droplets. We believe that our new code is a very powerful tool to investigate the hydrodynamics of molten droplets in the framework of the shock-wave heating model and has many potentials to be applied to various problems.

Chondrules are considered to have much information on dust particles and processes in the solar nebula. It is naturally expected that protoplanetary disks observed in present star forming regions have similar dust particles and processes, so study of chondrule formation may provide us great information on the formation of the planetary systems. Evaporation during chondrule melting may have resulted in depletion of volatile elements in chondrules. However, no evidence for a large degree of heavy-isotope enrichment has been reported in chondrules. In order to meet this observed constraint, the rapid heating rate at temperatures below the silicate solidus is required to suppress the isotopic fractionation. We have developed a new shock-wave heating model taking into account the radiative transfer of the dust thermal continuum emission and the line emission of gas molecules and calculated the thermal history of chondrules. We have found that optically-thin shock waves for the thermal continuum emission from dust particles can meet the rapid heating constraint, because the dust thermal emission does not keep the dust particles high temperature for a long time in the pre-shock region and dust particles are abruptly heated by the gas drag heating in the post-shock region. We have also derived the upper limit of optical depth of the pre-shock region using the radiative diffusion approximation, above which the rapid heating constraint is not satisfied. It is about 1 - 10.

Radiative transfer simulations are now at the forefront of numerical astrophysics. They are becoming crucial for an increasing number of astrophysical and cosmological problems; at the same time their computational cost has come to the reach of currently available computational power. Further progress is retarded by the considerable number of different algorithms (including various flavours of ray-tracing and moment schemes) developed, which makes the selection of the most suitable technique for a given problem a non-trivial task. Assessing the validity ranges, accuracy and performances of these schemes is the main aim of this paper, for which we have compared 11 independent RT codes on 5 test problems: (0) basic physics, (1) isothermal H II region expansion and (2) H II region expansion with evolving temperature, (3) I-front trapping and shadowing by a dense clump, (4) multiple sources in a cosmological density field. The outputs of these tests have been compared and differences analyzed. The agreement between the various codes is satisfactory although not perfect. The main source of discrepancy appears to reside in the multi-frequency treatment approach, resulting in different thicknesses of the ionized-neutral transition regions and different temperature structure. The present results and tests represent the most complete benchmark available for the development of new codes and improvement of existing ones. To this aim all test inputs and outputs are made publicly available in digital form.

The quasi-steady structure of super-critical accretion flows around a black hole is studied based on the two-dimensional radiation-hydrodynamical (2D-RHD) simulations. The super-critical flow is composed of two parts: the disk region and the outflow regions above and below the disk. Within the disk region the circular motion as well as the patchy density structure are observed, which is caused by Kelvin-Helmholtz instability and probably by convection. The mass-accretion rate decreases inward, roughly in proportion to the radius, and the remaining part of the disk material leaves the disk to form outflow because of strong radiation pressure force. We confirm that photon trapping plays an important role within the disk. Thus, matter can fall onto the black hole at a rate exceeding the Eddington rate. The emission is highly anisotropic and moderately collimated so that the apparent luminosity can exceed the Eddington luminosity by a factor of a few in the face-on view. The mass-accretion rate onto the black hole increases with increase of the absorption opacity (metalicity) of the accreting matter. This implies that the black hole tends to grow up faster in the metal rich regions as in starburst galaxies or star-forming regions.

A Compton camera is the most promising approach for gamma-ray detection in the energy region from several hundred keV to MeV, especially for application in high energy astrophysics. In order to obtain good angular resolution, semiconductor detectors such as silicon, germanium and cadmium telluride(CdTe) have several advantages over scintillation detectors, which have been used so far. Based on the recent advances of high resolution CdTe and silicon imaging detectors, we are working on a Si/CdTe Compton camera. We have developed 64-pixel CdTe detectors with a pixel size of 2mmx2mm and double-sided Si strip detectors(DSSDs) with a position resolution of 800 micron. As a prototype Si/CdTe Compton camera, we use a DSSD as a scatterer and two CdTe pixel detectors as an absorber. In order to verify its performance, we irradiate the camera with 100% linearly polarised 170keV gamma-rays and demonstrate the system works properly as a Compton camera. The resolution of the reconstructed scattering angle is 22 degrees(FWHM). Measurement of polarization is also reported. The polarimetric modulation factor is obtained to be 43%, which is consistent with the prediction of Monte Carlo simulations.

We are developing a Compton telescope based on high resolution Si and CdTe imaging devices in order to obtain a high sensitivity astrophysical observation in sub-MeV gamma-ray region. In this paper, recent results from the prototype Si/CdTe semiconductor Compton telescope are reported. The Compton telescope consists of a double-sided Si strip detector (DSSD) and CdTe pixel detectors, combined with low noise analog LSI, VA32TA. With this detector, we obtained Compton reconstructed images and spectra from line gamma-rays ranging from 81 keV up to 356 keV. The energy resolution is 3.8 keV and 7.9 keV at 122 keV and 356 keV, respectively, and the angular resolution is 9.9 degrees and 5.7 degrees at 122 keV and 356 keV, respectively.

Compton telescopes based on semiconductor technologies are being developed to explore the gamma-ray universe in an energy band 0.1--20 MeV, which is not well covered by the present or near-future gamma-ray telescopes. The key feature of such Compton telescopes is the high energy resolution that is crucial for high angular resolution and high background rejection capability. The energy resolution around 1 keV is required to approach physical limit of the angular resolution due to Doppler broadening. We have developed a low noise front-end ASIC (Application-Specific Integrated Circuit), VA32TA, to realize this goal for the readout of Double-sided Silicon Strip Detector (DSSD) and Cadmium Telluride (CdTe) pixel detector which are essential elements of the semiconductor Compton telescope. We report on the design and test results of the VA32TA. We have reached an energy resolution of 1.3 keV (FWHM) for 60 keV and 122 keV at 0 degree C with a DSSD and 1.7 keV (FWHM) with a CdTe detector.

A Semiconductor Multiple-Compton Telescope (SMCT) is being developed to explore the gamma-ray universe in an energy band 0.1--20 MeV, which is not well covered by the present or near-future gamma-ray telescopes. The key feature of the SMCT is the high energy resolution that is crucial for high angular resolution and high background rejection capability. We have developed prototype modules for a low noise Double-sided Silicon Strip Detector (DSSD) system which is an essential element of the SMCT. The geometry of the DSSD is optimized to achieve the lowest noise possible. A new front-end VLSI device optimized for low noise operation is also developed. We report on the design and test results of the prototype system. We have reached an energy resolution of 1.3 keV (FWHM) for 60 keV and 122 keV at 0 degree C.

A new method which enables us to estimate the inclination of Class I young stellar objects is proposed. Since Class I objects are not spherically symmetric, it is likely that the observed feature is sensitive to the inclination of the system. Thus, we construct a protostar model by carefully treating two-dimensional (2D) radiative transfer and radiative equilibrium. We show from the present 2D numerical simulations that the emergent luminosity L_SED,which is the frequency integration of spectral energy distribution (SED), depends strongly on the inclination of the system i, whereas the peak flux is insensitive to i. Based on this result, we introduce a novel indicator f_L, which is the ratio of L_SED to the peak flux, as a good measure for the inclination. By using f_L, we can determine the inclination regardless of the other physical parameters. The inclination would be determined by f_L within the accuracy of +- 5 degree, if the opening angle of bipolar outflows is specified by any other procedure. Since this spectrophotometric method is easier than a geometrical method or a full SED fitting method, this method could be a powerful tool to investigate the feature of protostars statistically with observational data which will be provided by future missions, such as SIRTF, ASTRO-F, and ALMA.

We explore the origin of the flat spectrum seen in some T Tauri stars by considering a three-component structure: a central star, a circumstellar disk, and a dusty halo. The radiative energy transport is faithfully treated by solving the angle- and frequency-dependent radiative transfer equation in two space dimensions assuming axisymmetry, and hence the radiative equilibrium temperature in the disk and halo is determined simultaneously. The disk is effectively heated by the scattering and reprocessing of stellar radiation through the halo. The large mid- to far-infrared excess originates from the photosphere of the warmed disk, resulting in a flat spectrum, as observed. The halo which we consider is observed as a compact reflection nebula, and is discriminated from extended, disk-like envelopes around flat-spectrum T Tauri stars. We show that the overall spectral shape of flat-spectrum T Tauri stars can be reproduced by the present disk-halo model.

Using 3D radiative transfer calculations on the reionization of an inhomogeneous universe, QSO absorption line systems are simulated and they are compared with observations of Ly$\alpha$ continuum depression at high redshifts. By this comparison, it is found that the metagalactic UV intensity decreases rapidly with $z$ at $z>4$ as $I_{21}=0.5\exp [3(4-z)]$, and the reionization must have taken place between $z=6$ and 10. Based on this time-dependence of UV background intensity, we explore the collapse of pregalactic clouds in the UV background, and find that the self-shielding is prominent above a mass scale as $ M_{BIF}=3.0\times 10^{11} M_\odot [(1+z_c)/5]^{-4.2}(I_{21}/0.5)^{0.6}. $ This mass scale predicts the bifurcation of galactic morphology, and by confrontation with observations it turns out that the bifurcation mass successfully discriminates between elliptical and spiral galaxies.

Full radiative transfer in 3D space has been solved to pursue the reionization history in an inhomogeneous universe. It has been shown that the reionization of an inhomogeneous universe is not a prompt event, but a fairly slow process. Also, QSO absorption line systems are simulated with using the results of radiative transfer calculations. Ly$\alpha$ continuum depression implies that the metagalactic UV intensity decreases rapidly with $z$ at $z>5$, and the reionization must have taken place between $z=6$ and 10. Finally, it is stressed that H$\alpha$ forest is a more powerful tool to probe the reionization history and the density fluctuations in the universe at $z>5$.