We examine the uncertainty of the calculation of the atmospheric neutrino flux and present a way to reduce it using accurately measured atmospheric muon flux. Considering the difference of the hadronic interaction model and the real one as a variation of hadronic interaction, we find a quantitative estimation method for the error of the atmospheric neutrino flux calculation from the residual of the reconstruction of the atmospheric muon flux observed in a precision experiment, by the study of atmospheric neutrino and muon fluxes response to the variation of hadronic interaction. However, the efficiencty of this method is largely dependent on the observation site of the atmospheric muon flux, as the relation of the error of the atmospheric neutrino flux calculation and the residual of the reconstruction of the atmospheric muon flux is also largely dependent on the muon observation site, especially for the low energy neutrinos. We calculate several observation sites, near Kamioka at sea level, same but 2770m a.s.l.., Hanle India (4500m a.s.l.), and at Balloon altitude ($\sim$ 32km). Then we estimate how stringently can the atmospheric muon reduce the error in the calculation of the atmospheric neutrino flux. We also discuss on the source of error which is difficult to reduce by only the observation of atmospheric muon.

Comet 21P/Giacobini-Zinner (hereafter, comet 21P/G-Z) is a Jupiter-family comet and a parent comet of the October Draconids meteor shower. If meteoroids originating from a Jupiter-family comet contain complex organic molecules, such as amino acids, they are essential pieces of the puzzle regarding the origin of life on Earth. We observed comet 21P/G-Z in the mid-infrared wavelength region using the Cooled Mid-infrared Camera and Spectrometer (COMICS) on the 8.2 m Subaru Telescope on UT 2005 July 5. Here, we report the unidentified infrared (UIR) emission features of comet 21P/G-Z, which are likely due to complex organic molecules (both aliphatic and aromatic hydrocarbons), and the thermal emission from amorphous/crystalline silicates and amorphous carbon grains in its mid-infrared low-resolution spectrum. The UIR features at ~8.2 micron, ~8.5 micron, and ~11.2 micron found in the spectrum of comet 21P/G-Z could be attributed to polycyclic aromatic hydrocarbons (or hydrogenated amorphous carbons) contaminated by N- or O-atoms, although part of the feature at ~11.2 micron comes from crystalline olivine. The other feature at ~9.2 micron might originate from aliphatic hydrocarbons. Comet 21P/G-Z is enriched in complex organic molecules. Considering that the derived mass fraction of crystalline silicates in comet 21P/G-Z is typical of comets, we propose that the comet originated from a circumplanetary disk of giant planets (similar to Jupiter and Saturn) where was warmer than the typical comet-forming region (5-30 au from the Sun) and was suitable for the formation of complex organic molecules. Comets from circumplanetary disks might be enriched in complex organic molecules, such as comet 21P/G-Z, and may have provided pre-biotic molecules to ancient Earth by direct impact or meteor showers.

In this paper, we analyse the upper limit fluxes of sub-millimeter ortho-H$_{2}$$^{16}$O 321 GHz, para-H$_{2}$$^{18}$O 322 GHz, and HDO 335 GHz lines from the protoplanetary disk around the Herbig Ae star HD 163296, using the Atacama Large Millimeter/Submillimeter Array (ALMA). These water lines are considered to be the best candidate sub-millimeter lines to locate the position of the H$_{2}$O snowline, on the basis of our previous model calculations. We compare the upper limit fluxes with the values calculated by our models with dust emission included, and we constrain the line emitting region and the dust opacity from the observations. We conclude that, if the outer edge of the region with high water vapor abundance and if the position of the water snowline are beyond 8 au also, the mm dust opacity will have a value larger than 2.0 cm$^{2}$ g$^{-1}$. In addition, the position of the water snowline will be inside 20 au, if the mm dust opacity is 2.0 cm$^{2}$ g$^{-1}$. Future observations of the dust continuum emission at higher angular resolution and sub-millimeter water lines with longer observation time are required to clarify the detailed structures and the position of the H$_{2}$O snowline in the disk midplane.

Dust grains of crystalline silicate, which are rarely presented in interstellar space, were found in cometary nuclei. These crystalline silicates are thought to have formed by annealing of amorphous silicate grains or direct condensation of gaseous materials near the Sun in the solar nebula, and incorporated into cometary nuclei in the cold comet-forming region after radial transportation of grains in the solar nebula. Abundances of the crystalline silicate dust grains were therefore expected to be smaller farther from the Sun. We aim to better understand the formation mechanism of minerals incorporated into comet 17P/Holmes based on its mineral abundances. To derive the mineral composition of comet 17P/Holmes, we applied a thermal emission model for cometary dust grains to mid-infrared spectra of comet 17P/Holmes taken with the Cooled Mid-Infrared Camera and Spectrograph (COMICS) mounted on the Subaru Telescope a few days later the great outburst in October 2007. The resulting mass fraction of crystalline silicate, f_ cry, and an olivine-to-pyroxene abundance ratio, f_OP, are f_cry = 0.31 +/- 0.03 and f_OP = 1.20 +0.16/-0.12, respectively. Based on a simple consideration of the mixing of dust grains originating in both the interstellar medium and solar nebula, the minerals of 17P/Holmes formed by non-equilibrium condensation. This result is consistent with theoretical and experimental predictions for vaporization and condensation of olivine in the solar nebula.

In this paper, we extend the results presented in our former papers (Notsu et al. 2016, 2017) on using ortho-$\mathrm{H_2}$$^{16}\mathrm{O}$ line profiles to constrain the location of the $\mathrm{H_2O}$ snowline in T Tauri and Herbig Ae disks, to include sub-millimeter para-$\mathrm{H_2}$$^{16}\mathrm{O}$ and ortho- and para-$\mathrm{H_2}$$^{18}\mathrm{O}$ lines. Since the number densities of the ortho- and para-H$_{2}$$^{18}$O molecules are about 560 times smaller than their $^{16}$O analogues, they trace deeper into the disk than the ortho-H$_{2}$$^{16}$O lines (down to $z=0$, i.e., the midplane). Thus these H$_{2}$$^{18}$O lines are potentially better probes of the position of the H$_{2}$O snowline at the disk midplane, depending on the dust optical depth. The values of the Einstein $A$ coefficients of sub-millimeter candidate water lines tend to be lower (typically $<$$10^{-4}$ s$^{-1}$) than infrared candidate water lines (Notsu et al. 2017). Thus in the sub-millimeter candidate water line cases, the local intensity from the outer optically thin region in the disk is around $10^{4}$ times smaller than that in the infrared candidate water line cases. Therefore, in the sub-millimeter lines, especially H$_{2}$$^{18}$O and para-H$_{2}$$^{16}$O lines with relatively lower upper state energies ($\sim$ a few 100K) can also locate the position of the $\mathrm{H_2O}$ snowline. We also investigate the possibility of future observations with ALMA to identify the position of the water snowline. There are several candidate water lines that trace the hot water vapor inside the $\mathrm{H_2O}$ snowline in ALMA Bands $5-10$.

Transitional disks around the Herbig Ae/Be stars are fascinating targets in the contexts of disk evolution and also planet formation. Oph IRS 48 is one of such Herbig Ae stars, which shows an inner dust cavity and azimuthally lopsided large dust distribution. We present new images of Oph IRS 48 at eight mid-infrared (MIR) wavelengths from 8.59 to 24.6\,$\mu {\rm m}$ taken with the COMICS mounted on the 8.2\,m Subaru Telescope. The N-band (7 to 13\,$\mu {\rm m}$) images show that the flux distribution is centrally peaked with a slight spatial extent, while the Q-band (17 to 25\,$\mu {\rm m}$) images show asymmetric double peaks (east and west). Using 18.8 and 24.6\,$\mu$m images, we derived the dust temperature at both east and west peaks to be 135$\pm$22 K. Thus, the asymmetry may not be attributed to a difference in the temperature. % thus other reason is necessary to explain the asymmetry of east and west peaks in Q-band. Comparing our results with previous modeling works, we conclude that the inner disk is aligned to the outer disk. A shadow casted by the optically thick inner disk has a great influence on the morphology of MIR thermal emission from the outer disk.

Observationally measuring the location of the H$_{2}$O snowline is crucial for understanding the planetesimal and planet formation processes, and the origin of water on Earth. In disks around Herbig Ae stars ($T_{\mathrm{*}}\sim$ 10,000K, $M_{\mathrm{*}}\gtrsim$ 2.5$M_{\bigodot}$), the position of the H$_{2}$O snowline is further from the central star compared with that around cooler, and less massive T Tauri stars. Thus, the H$_{2}$O emission line fluxes from the region within the H$_{2}$O snowline are expected to be stronger. In this paper, we calculate the chemical composition of a Herbig Ae disk using chemical kinetics. Next, we calculate the H$_{2}$O emission line profiles, and investigate the properties of candidate water lines across a wide range of wavelengths (from mid-infrared to sub-millimeter) that can locate the position of the H$_{2}$O snowline. Those line identified have small Einstein $A$ coefficients ($\sim 10^{-6} -10^{-3}$ s$^{-1}$) and relatively high upper state energies ($\sim$ 1000K). The total fluxes tend to increase with decreasing wavelengths. We investigate the possibility of future observations (e.g., ALMA, SPICA/SMI-HRS) to locate the position of the H$_{2}$O snowline. Since the fluxes of those identified lines from Herbig Ae disks are stronger than those from T Tauri disks, the possibility of a successful detection is expected to increase for a Herbig Ae disk.

Traces of 2-3 Myr old 60Fe were recently discovered in a manganese crust and in lunar samples. We have found that this signal is extended in time and is present in globally distributed deep-sea archives. A second 6.5-8.7 Myr old signature was revealed in a manganese crust. The existence of the Local Bubble hints to a recent nearby supernova-activity starting 13 Myr ago. With analytical and numerical models generating the Local Bubble, we explain the younger 60Fe-signature and thus link the evolution of the solar neighborhood to terrestrial anomalies.

Inside the H$_{2}$O snowline of protoplanetary disks, water evaporates from the dust-grain surface into the gas phase, whereas it is frozen out on to the dust in the cold region beyond the snowline. H$_{2}$O ice enhances the solid material in the cold outer part of a disk, which promotes the formation of gas-giant planet cores. We can regard the H$_{2}$O snowline as the surface that divides the regions between rocky and gaseous giant planet formation. Thus observationally measuring the location of the H$_{2}$O snowline is crucial for understanding the planetesimal and planet formation processes, and the origin of water on Earth. In this paper, we find candidate water lines to locate the H$_{2}$O snowline through future high-dispersion spectroscopic observations. First, we calculate the chemical composition of the disk and investigate the abundance distributions of H$_{2}$O gas and ice, and the position of the H$_{2}$O snowline. We confirm that the abundance of H$_{2}$O gas is high not only in the hot midplane region inside the H$_{2}$O snowline but also in the hot surface layer of the outer disk. Second, we calculate the H$_{2}$O line profiles and identify those H$_{2}$O lines which are promising for locating the H$_{2}$O snowline: the identified lines are those which have small Einstein $A$ coefficients and high upper state energies. The wavelengths of the candidate H$_{2}$O lines range from mid-infrared to sub-millimeter, and they overlap with the regions accessible to ALMA and future mid-infrared high dispersion spectrographs (e.g., TMT/MICHI, SPICA).

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.

We investigate the dust and gas distribution in the disk around HD 142527 based on ALMA observations of dust continuum, 13CO(3-2), and C18O(3-2) emission. The disk shows strong azimuthal asymmetry in the dust continuum emission, while gas emission is more symmetric. In this paper, we investigate how gas and dust are distributed in the dust-bright northern part of the disk and in the dust-faint southern part. We construct two axisymmetric disk models. One reproduces the radial profiles of the continuum and the velocity moments 0 and 1 of CO lines in the north and the other reproduces those in the south. We have found that the dust is concentrated in a narrow ring having ~50AU width (in FWHM; w_d=30AU in our parameter definition) located at ~170-200AU from the central star. The dust particles are strongly concentrated in the north. We have found that the dust surface density contrast between the north and south amounts to ~70. Compared to the dust, the gas distribution is more extended in the radial direction. We find that the gas component extends at least from ~100AU to ~250AU from the central star, and there should also be tenuous gas remaining inside and outside of these radii. The azimuthal asymmetry of gas distribution is much smaller than dust. The gas surface density differs only by a factor of ~3-10 between the north and south. Hence, gas-to-dust ratio strongly depends on the location of the disk: ~30 at the location of the peak of dust distribution in the south and ~3 at the location of the peak of dust distribution in the north. Despite large uncertainties, the overall gas-to-dust ratio is inferred to be ~10-30, indicating that the gas depletion may have already been under way.

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

In this paper, we extend the calculation of the atmospheric neutrino flux~\cite{hkkm2004,hkkms2006,hkkm2011} to the sites in polar and tropical regions. In our earliest full 3D-calculation~\cite{hkkm2004}, we used DPMJET-III~\cite{dpm} for the hadronic interaction model above 5~GeV, and NUCRIN~\cite{nucrin} below 5~GeV. We modified DPMJET-III as in Ref.~\cite{hkkms2006} to reproduce the experimental muon spectra better, mainly using the data observed by BESS group~\cite{BESSTeVpHemu}. In a recent work~\cite{hkkm2011}, we introduced JAM interaction model for the low energy hadronic interactions. JAM is a nuclear interaction model developed with PHITS (Particle and Heavy-Ion Transport code System)~\cite{phits}. In Ref.~\cite{hkkm2011}, we could reproduce the observed muon flux at the low energies at balloon altitude with DPMJET-III above 32 GeV and JAM below that better than the combination of DPMJET-III above 5~GeV and NUCRIN below that. Besides the interaction model, we have also improved the calculation scheme according to the increase of available computational power, such as the "virtual detector correction" introduced in Ref.~\cite{hkkms2006} and the optimization of it in Ref.~\cite{hkkm2011}. The statistics of the Monte Carlo simulation is also improved at every step of the work.

Coronagraphic imagery of the circumstellar disk around HD 169142 in H-band polarized intensity (PI) with Subaru/HiCIAO is presented. The emission scattered by dust particles at the disk surface in 0.2" <= r <= 1.2", or 29 <= r <= 174 AU, is successfully detected. The azimuthally-averaged radial profile of the PI shows a double power-law distribution, in which the PIs in r=29-52 AU and r=81.2-145 AU respectively show r^{-3}-dependence. These two power-law regions are connected smoothly with a transition zone (TZ), exhibiting an apparent gap in r=40-70 AU. The PI in the inner power-law region shows a deep minimum whose location seems to coincide with the point source at \lambda = 7 mm. This can be regarded as another sign of a protoplanet in TZ. The observed radial profile of the PI is reproduced by a minimally flaring disk with an irregular surface density distribution or with an irregular temperature distribution or with the combination of both. The depletion factor of surface density in the inner power-law region (r< 50 AU) is derived to be <= 0.16 from a simple model calculation. The obtained PI image also shows small scale asymmetries in the outer power-law region. Possible origins for these asymmetries include corrugation of the scattering surface in the outer region, and shadowing effect by a puffed up structure in the inner power-law region.

We imaged circumstellar disks around 22 Herbig Ae/Be stars at 25 \mu m using Subaru/COMICS and Gemini/T-ReCS. Our sample consists of equal numbers of objects belonging to the two categories defined by Meeus et al. (2001); 11 group I (flaring disk) and II (at disk) sources. We find that group I sources tend to show more extended emission than group II sources. Previous studies have shown that the continuous disk is hard to be resolved with 8 meter class telescopes in Q-band due to the strong emission from the unresolved innermost region of the disk. It indicates that the resolved Q-band sources require a hole or gap in the disk material distribution to suppress the contribution from the innermost region of the disk. As many group I sources are resolved at 25 \mu m, we suggest that many, not all, group I Herbig Ae/Be disks have a hole or gap and are (pre-)transitional disks. On the other hand, the unresolved nature of many group II sources at 25 \mu m supports that group II disks have continuous at disk geometry. It has been inferred that group I disks may evolve into group II through settling of dust grains to the mid-plane of the proto-planetary disk. However, considering growing evidence for the presence of a hole or gaps in the disk of group I sources, such an evolutionary scenario is unlikely. The difference between groups I and II may reflect different evolutionary pathways of protoplanetary disks.

We report the first independent, second-epoch (re-)detection of a directly-imaged protoplanet candidate. Using $L^\prime$ high-contrast imaging of HD 100546 taken with the Near-Infrared Coronagraph and Imager (NICI) on Gemini South, we recover `HD 100546 b' with a position and brightness consistent with the original VLT/NaCo detection from Quanz et al, although data obtained after 2013 will be required to decisively demonstrate common proper motion. HD 100546 b may be spatially resolved, up to $\approx$ 12-13 AU in diameter, and is embedded in a finger of thermal IR bright, polarized emission extending inwards to at least 0.3". Standard hot-start models imply a mass of $\approx$ 15 $M_{J}$. But if HD 100546 b is newly formed or made visible by a circumplanetary disk, both of which are plausible, its mass is significantly lower (e.g. 1--7 $M_{J}$). Additionally, we discover a thermal IR-bright disk feature, possibly a spiral density wave, at roughly the same angular separation as HD 100546 b but 90 degrees away. Our interpretation of this feature as a spiral arm is not decisive, but modeling analyses using spiral density wave theory implies a wave launching point exterior to $\approx$ 0.45" embedded within the visible disk structure: plausibly evidence for a second, hitherto unseen wide-separation planet. With one confirmed protoplanet candidate and evidence for 1--2 others, HD 100546 is an important evolutionary precursor to intermediate-mass stars with multiple super-jovian planets at moderate/wide separations like HR 8799.

We present a new method of analysis for determining the surface geometry of five protoplanetary disks observed with near-infrared imaging polarimetry using Subaru-HiCIAO. Using as inputs the observed distribution of polarized intensity (PI), disk inclination, assumed properties for dust scattering, and other reasonable approximations, we calculate a differential equation to derive the surface geometry. This equation is numerically integrated along the distance from the star at a given position angle. We show that, using these approximations, the local maxima in the PI distribution of spiral arms (SAO 206462, MWC 758) and rings (2MASS J16042165-2130284, PDS 70) are associated with local concave-up structures on the disk surface. We also show that the observed presence of an inner gap in scattered light still allows the possibility of a disk surface that is parallel to the light path from the star, or a disk that is shadowed by structures in the inner radii. Our analysis for rings does not show the presence of a vertical inner wall as often assumed in studies of disks with an inner gap. Finally, we summarize the implications of spiral and ring structures as potential signatures of ongoing planet formation.

We report ALMA observations of dust continuum, 13CO J=3--2, and C18O J=3--2 line emission toward a gapped protoplanetary disk around HD 142527. The outer horseshoe-shaped disk shows the strong azimuthal asymmetry in dust continuum with the contrast of about 30 at 336 GHz between the northern peak and the southwestern minimum. In addition, the maximum brightness temperature of 24 K at its northern area is exceptionally high at 160 AU from a star. To evaluate the surface density in this region, the grain temperature needs to be constrained and was estimated from the optically thick 13CO J=3--2 emission. The lower limit of the peak surface density was then calculated to be 28 g cm-2 by assuming a canonical gas-to-dust mass ratio of 100. This finding implies that the region is locally too massive to withstand self-gravity since Toomre's Q <~1--2, and thus, it may collapse into a gaseous protoplanet. Another possibility is that the gas mass is low enough to be gravitationally stable and only dust grains are accumulated. In this case, lower gas-to-dust ratio by at least 1 order of magnitude is required, implying possible formation of a rocky planetary core.

*Context The evolution of young massive protoplanetary disks toward planetary systems is expected to include the formation of gaps and the depletion of dust and gas. *Aims A special group of flaring disks around Herbig Ae/Be stars do not show prominent silicate emission features. We focus our attention on four key Herbig Ae/Be stars to understand the structural properties responsible for the absence of silicate feature emission. *Methods We investigate Q- and N-band images taken with Subaru/COMICS, Gemini South/T-ReCS and VLT/VISIR. Our radiative transfer modeling solutions require a separation of inner- and outer- disks by a large gap. From this we characterize the radial density structure of dust and PAHs in the disk. *Results The inner edge of the outer disk has a high surface brightness and a typical temperature between ~100-150 K and therefore dominates the emission in the Q-band. We derive radii of the inner edge of the outer disk of 34, 23, 30 and 63 AU for HD97048, HD169142, HD135344B and Oph IRS 48 respectively. For HD97048 this is the first detection of a disk gap. The continuum emission in the N-band is not due to emission in the wings of PAHs. This continuum emission can be due to VSGs or to thermal emission from the inner disk. We find that PAH emission is not always dominated by PAHs on the surface of the outer disk. *Conclusions. The absence of silicate emission features is due to the presence of large gaps in the critical temperature regime. Many, if not all Herbig disks with Spectral Energy Distribution (SED) classification `group I' are disks with large gaps and can be characterized as (pre-) transitional. An evolutionary path from the observed group I to the observed group II sources seems no longer likely. Instead, both might derive from a common ancestor.

We present the calculation of the atmospheric neutrino fluxes for the neutrino experiments proposed at INO, South Pole and Pyh\"asalmi. Neutrino fluxes have been obtained using ATMNC, a simulation code for cosmic ray in the atmosphere. Even using the same primary flux model and the interaction model, the calculated atmospheric neutrino fluxes are different for the different sites due to the geomagnetic field. The prediction of these fluxes in the present paper would be quite useful in the experimental analysis.

We present high resolution H-band polarized intensity (PI; FWHM = 0."1: 14 AU) and L'-band imaging data (FWHM = 0."11: 15 AU) of the circumstellar disk around the weak-lined T Tauri star PDS 70 in Centaurus at a radial distance of 28 AU (0."2) up to 210 AU (1."5). In both images, a giant inner gap is clearly resolved for the first time, and the radius of the gap is ~70 AU. Our data show that the geometric center of the disk shifts by ~6 AU toward the minor axis. We confirm that the brown dwarf companion candidate to the north of PDS 70 is a background star based on its proper motion. As a result of SED fitting by Monte Carlo radiative transfer modeling, we infer the existence of an optically thick inner disk at a few AU. Combining our observations and modeling, we classify the disk of PDS 70 as a pre-transitional disk. Furthermore, based on the analysis of L'-band imaging data, we put an upper limit mass of companions at ~30 to ~50MJ within the gap. Taking account of the presence of the large and sharp gap, we suggest that the gap could be formed by dynamical interactions of sub-stellar companions or multiple unseen giant planets in the gap.

We present H-band polarimetric imagery of UX Tau A taken with HiCIAO/AO188 on the Subaru Telescope. UX Tau A has been classified as a pre-transitional disk object, with a gap structure separating its inner and outer disks. Our imagery taken with the 0.15 (21 AU) radius coronagraphic mask has revealed a strongly polarized circumstellar disk surrounding UX Tau A which extends to 120 AU, at a spatial resolution of 0.1 (14 AU). It is inclined by 46 \pm 2 degree as the west side is nearest. Although SED modeling and sub-millimeter imagery suggested the presence of a gap in the disk, with the inner edge of the outer disk estimated to be located at 25 - 30 AU, we detect no evidence of a gap at the limit of our inner working angle (23 AU) at the near-infrared wavelength. We attribute the observed strong polarization (up to 66 %) to light scattering by dust grains in the disk. However, neither polarization models of the circumstellar disk based on Rayleigh scattering nor Mie scattering approximations were consistent with the observed azimuthal profile of the polarization degrees of the disk. Instead, a geometric optics model of the disk with nonspherical grains with the radii of 30 micron meter is consistent with the observed profile. We suggest that the dust grains have experienced frequent collisional coagulations and have grown in the circumstellar disk of UX Tau A.

The disk around the Herbig Ae star HD\,169142 was imaged and resolved at 18.8 and 24.5\,$\mu$m using Subaru/COMICS. We interpret the observations using a 2D radiative transfer model and find evidence for the presence of a large gap. The MIR images trace dust that emits at the onset of the strong rise in the spectral energy distribution (SED) at 20\,$\mu$m, therefore are very sensitive to the location and characteristics of the inner wall of the outer disk and its dust. We determine the location of the wall to be 23$^{+3}_{-5}$\,AU from the star. An extra component of hot dust must exist close to the star. We find that a hydrostatic optically thick inner disk does not produce enough flux in the NIR and an optically thin geometrically thick component is our solution to fit the SED. Considering the recent findings of gaps and holes in a number of Herbig Ae/Be group I disks, we suggest that such disk structures may be common in group I sources. Classification as group I should be considered a support for classification as a transitional disk, though improved imaging surveys are needed to support this speculation.

We present high-resolution, H-band, imaging observations, collected with Subaru/HiCIAO, of the scattered light from the transitional disk around SAO 206462 (HD 135344B). Although previous sub-mm imagery suggested the existence of the dust-depleted cavity at r~46AU, our observations reveal the presence of scattered light components as close as 0.2" (~28AU) from the star. Moreover, we have discovered two small-scale spiral structures lying within 0.5" (~70AU). We present models for the spiral structures using the spiral density wave theory, and derive a disk aspect ratio of h~0.1, which is consistent with previous sub-mm observations. This model can potentially give estimates of the temperature and rotation profiles of the disk based on dynamical processes, independently from sub-mm observations. It also predicts the evolution of the spiral structures, which can be observable on timescales of 10-20 years, providing conclusive tests of the model. While we cannot uniquely identify the origin of these spirals, planets embedded in the disk may be capable of exciting the observed morphology. Assuming that this is the case, we can make predictions on the locations and, possibly, the masses of the unseen planets. Such planets may be detected by future multi-wavelengths observations.

Transitional circumstellar disks around young stellar objects have a distinctive infrared deficit around 10 microns in their Spectral Energy Distributions (SED), recently measured by the Spitzer Infrared Spectrograph (IRS), suggesting dust depletion in the inner regions. These disks have been confirmed to have giant central cavities by imaging of the submillimeter (sub-mm) continuum emission using the Submillimeter Array (SMA). However, the polarized near-infrared scattered light images for most objects in a systematic IRS/SMA cross sample, obtained by HiCIAO on the Subaru telescope, show no evidence for the cavity, in clear contrast with SMA and Spitzer observations. Radiative transfer modeling indicates that many of these scattered light images are consistent with a smooth spatial distribution for micron-sized grains, with little discontinuity in the surface density of the micron-sized grains at the cavity edge. Here we present a generic disk model that can simultaneously account for the general features in IRS, SMA, and Subaru observations. Particularly, the scattered light images for this model are computed, which agree with the general trend seen in Subaru data. Decoupling between the spatial distributions of the micron-sized dust and mm-sized dust inside the cavity is suggested by the model, which, if confirmed, necessitates a mechanism, such as dust filtration, for differentiating the small and big dust in the cavity clearing process. Our model also suggests an inwardly increasing gas-to-dust-ratio in the inner disk, and different spatial distributions for the small dust inside and outside the cavity, echoing the predictions in grain coagulation and growth models.

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.

We present a study of the radial distribution of dust species in young brown dwarf disks. Our work is based on a compositional analysis of the 10 and 20 micron silicate emission features for brown dwarfs in the Taurus-Auriga star-forming region. A fundamental finding of our work is that brown dwarfs exhibit stronger signs of dust processing in the cold component of the disk, compared to the higher mass T Tauri stars in Taurus. For nearly all of our targets, we find a flat disk structure, which is consistent with the stronger signs of dust processing observed in these disks. For the case of one brown dwarf, 2M04230607, we find the forsterite mass fraction to be a factor of ~3 higher in the outer disk compared to the inner disk region. Simple large-scale radial mixing cannot account for this gradient in the dust chemical composition, and some local crystalline formation mechanism may be effective in this disk. The relatively high abundance of crystalline silicates in the outer cold regions of brown dwarf disks provides an interesting analogy to comets. In this context, we have discussed the applicability of the various mechanisms that have been proposed for comets on the formation and the outward transport of high-temperature material. We also present Chandra X-ray observations for two Taurus brown dwarfs, 2M04414825 and CFHT-BD-Tau 9. We find 2M04414825, which has a ~12% crystalline mass fraction, to be more than an order of magnitude brighter in X-ray than CFHT-BD-Tau 9, which has a ~35% crystalline mass fraction. Combining with previous X-ray data, we find the inner disk crystalline mass fractions to be anti-correlated with the X-ray strength.

We present the SPICA Coronagraphic Instrument (SCI), which has been designed for a concentrated study of extra-solar planets (exoplanets). SPICA mission provides us with a unique opportunity to make high contrast observations because of its large telescope aperture, the simple pupil shape, and the capability for making infrared observations from space. The primary objectives for the SCI are the direct coronagraphic detection and spectroscopy of Jovian exoplanets in infrared, while the monitoring of transiting planets is another important target. The specification and an overview of the design of the instrument are shown. In the SCI, coronagraphic and non-coronagraphic modes are applicable for both an imaging and a spectroscopy. The core wavelength range and the goal contrast of the coronagraphic mode are 3.5--27$\mu$m, and 10$^{-6}$, respectively. Two complemental designs of binary shaped pupil mask coronagraph are presented. The SCI has capability of simultaneous observations of one target using two channels, a short channel with an InSb detector and a long wavelength channel with a Si:As detector. We also give a report on the current progress in the development of key technologies for the SCI.

We present mid-infrared narrow-band images of the Orion BN/KL region, and N-band low-resolution spectra of IRc2 and the nearby radio source "I." The distributions of the silicate absorption strength and the color temperature have been revealed with a sub-arcsecond resolution. The detailed structure of the 7.8 micron/12.4 micron color temperature distribution was resolved in the vicinity of IRc2. A mid-infrared counterpart to source I has been detected as a large color temperature peak. The color temperature distribution shows an increasing gradient from IRc2 toward source I, and no dominant temperature peak is seen at IRc2. The spectral energy distribution of IRc2 could be fitted by a two-temperature component model, and the "warmer component" of the infrared emission from IRc2 could be reproduced by scattering of radiation from source I. IRc2 itself is not self-luminous, but is illuminated and heated by an embedded luminous young stellar object located at source I.

We investigate gravitational lensing effects of an extrasolar planet transiting its host star. We focus on the `rising spikes' of the light curve just before and after the transit, which is a peculiar feature of the gravitational lensing, and find that it could be a novel observable for determining physical parameters. Detectability of such an effect is also discussed.

We report high-resolution 1.6 $\micron$ polarized intensity ($PI$) images of the circumstellar disk around the Herbig Ae star AB Aur at a radial distance of 22 AU ($0."15$) up to 554 AU (3.$"$85), which have been obtained by the high-contrast instrument HiCIAO with the dual-beam polarimetry. We revealed complicated and asymmetrical structures in the inner part ($\lesssim$140 AU) of the disk, while confirming the previously reported outer ($r$ $\gtrsim$200 AU) spiral structure. We have imaged a double ring structure at $\sim$40 and $\sim$100 AU and a ring-like gap between the two. We found a significant discrepancy of inclination angles between two rings, which may indicate that the disk of AB Aur is warped. Furthermore, we found seven dips (the typical size is $\sim$45 AU or less) within two rings as well as three prominent $PI$ peaks at $\sim$40 AU. The observed structures, including a bumpy double ring, a ring-like gap, and a warped disk in the innermost regions, provide essential information for understanding the formation mechanism of recently detected wide-orbit ($r$ $>$20 AU) planets.

We present the calculation of the atmospheric neutrino fluxes with an interaction model named JAM, which is used in PHITS (Particle and Heavy-Ion Transport code System). The JAM interaction model agrees with the HARP experiment a little better than DPMJET-III. After some modifications, it reproduces the muon flux below 1~GeV/c at balloon altitudes better than the modified-DPMJET-III which we used for the calculation of atmospheric neutrino flux in previous works. Some improvements in the calculation of atmospheric neutrino flux are also reported.

The recent findings of gas giant planets around young A-type stars suggest that disks surrounding Herbig Ae/Be stars will develop planetary systems. An interesting case is HD142527, for which previous observations revealed a complex circumstellar environment and an unusually high ratio of infrared to stellar luminosity. Its properties differ considerably from other Herbig Ae/Be stars. This suggests that the disk surrounding HD142527 is in an uncommon evolutionary stage. We aim for a better understanding of the geometry and evolutionary status of the circumstellar material around the Herbig Ae/Be star HD142527. We map the composition and spatial distribution of the dust around HD142527. We analyze SEST and ATCA millimeter data, VISIR N and Q-band imaging and spectroscopy. We gather additional relevant data from the literature. We use the radiative transfer code MCMax to construct a model of the geometry and density structure of the circumstellar matter, which fits all of the observables satisfactorily. We find that the disk of HD142527 has three geometrically distinct components separated by a disk gap running from 30 to 130 AU. There is a geometrically flat inner disk running from 0.3 AU up to 30 AU; an optically thin halo-like component of dust in the inner disk regions; and a massive self-shadowed outer disk running from 130 AU up to 200 AU. We derived a total dust mass in small grains of 10^-3 M0 and a vertical height of the inner wall of the outer disk of h = 60 AU. ....

We propose a new analytic scaling of the cut-off frequency of synchrotron radiation from active galactic nucleus (AGN) jets that are nonuniformly filled with many filaments. The theoretical upper limit is provided independent of magnetic intensity, spectral index, coherence and correlation length of filamentary turbulence, etc., such that \nu_c\simeq 6\times 10^{20}\delta[(r-1)/r]^{4/3}(b/10^{-4}) Hz, where \delta, r and b are the Doppler beaming factor, shock-compression ratio and energy-density ratio of the perturbed/local mean magnetic field of the filaments, respectively. Combining our results with observational data for 18 extragalactic sources, a constraint on the filament correlation length is found, in order to give the number scaling of filaments. The results suggest that, in particular, the jets of compact BL Lacs possess a large number of filaments with transverse size scale smaller than the emission-region size. The novel concept of the quantization of flowing plasma is suggested.

Blazars have been regarded as one of the most powerful sources of the highest energy cosmic rays and also their byproducts, neutrinos. Provided that a magnetized filamentary system is established in a blazar jet as well, we could apply the mechanism of multi-stage diffusive shock acceleration to a feasible TeV emitter, Mrk 501 to evaluate the achievable maximum energy of protons. Taking conceivable energy restriction into account systematically, it seems adequate to say that EeV-protons are produced at this site by our present model. We also estimate neutrino fluxes generated by these accelerated protons and discuss the detectability based on an updated kilometre-scale telescope such as IceCube.

The disk around AB Aur was imaged and resolved at 24.6\,$\mu$m using the Cooled Mid-Infrared Camera and Spectrometer on the 8.2m Subaru Telescope. The gaussian full-width at half-maximum of the source size is estimated to be 90 $\pm$ 6 AU, indicating that the disk extends further out at 24.6\,$\mu$m than at shorter wavelengths. In order to interpret the extended 24.6\,$\mu$m image, we consider a disk with a reduced surface density within a boundary radius $R_c$, which is motivated by radio observations that suggest a reduced inner region within about 100 AU from the star. Introducing the surface density reduction factor $f_c$ for the inner disk, we determine that the best match with the observed radial intensity profile at 24.6\,$\mu$m is achieved with $R_c$=88 AU and $f_c$=0.01. We suggest that the extended emission at 24.6\,$\mu$m is due to the enhanced emission from a wall-like structure at the boundary radius (the inner edge of the outer disk), which is caused by a jump in the surface density at $R_c$. Such reduced inner disk and geometrically thick outer disk structure can also explain the more point-like nature at shorter wavelengths. We also note that this disk geometry is qualitatively similar to a pre-transitional disk, suggesting that the AB Aur disk is in a pre-transitional disk phase.

We present H- and Ks-band imaging data resolving the gap in the transitional disk around LkCa 15, revealing the surrounding nebulosity. We detect sharp elliptical contours delimiting the nebulosity on the inside as well as the outside, consistent with the shape, size, ellipticity, and orientation of starlight reflected from the far-side disk wall, whereas the near-side wall is shielded from view by the disk's optically thick bulk. We note that forward-scattering of starlight on the near-side disk surface could provide an alternate interpretation of the nebulosity. In either case, this discovery provides confirmation of the disk geometry that has been proposed to explain the spectral energy distributions (SED) of such systems, comprising an optically thick outer disk with an inner truncation radius of ~46 AU enclosing a largely evacuated gap. Our data show an offset of the nebulosity contours along the major axis, likely corresponding to a physical pericenter offset of the disk gap. This reinforces the leading theory that dynamical clearing by at least one orbiting body is the cause of the gap. Based on evolutionary models, our high-contrast imagery imposes an upper limit of 21 Jupiter masses on companions at separations outside of 0.1" and of 13 Jupiter masses outside of 0.2". Thus, we find that a planetary system around LkCa 15 is the most likely explanation for the disk architecture.

The SPace Infrared telescope for Cosmology and Astrophysics (SPICA) is a proposed mid-to-far infrared (4-200 um) astronomy mission, scheduled for launch in 2017. A single, 3.5m aperture telescope would provide superior image quality at 5-200 um, and its very cold (~5 K) instrumentation would provide superior sensitivity in the 25-200 um wavelength regimes. This would provide a breakthrough opportunity for studies of exoplanets, protoplanetary and debris disk, and small solar system bodies. This paper summarizes the potential scientific impacts for the proposed instrumentation.

We evaluate the achievable maximum energy of nuclei diffusively accelerated by shock wave in the jet of Cen A, based on an updated model involving the stochastic magnetic fields that are responsible for recent synchrotron X-ray measurements. For the maximum energy analysis, conceivable energy constraints from spatiotemporal scales are systematically considered for the jet-wide including discrete X-ray knots. We find that in the inner region within ~1 arcmin from galactic core, which includes knots AX and BX, proton and iron nucleus can be accelerated to 10^{19}-10^{20} and 10^{21} eV (10-100 EeV and ZeV) ranges, respectively. The upper cutoff energy of the very energetic neutrinos produced via photopion interaction is also provided. These are essential for identifying the acceleration site of the ultra-high energy cosmic ray detected in the Pierre Auger Observatory, which signifies the arrival from nearby galaxies including Cen A.

We made mid-infrared observations of the 10Msun Herbig Be star HD200775 with the Cooled Mid-Infrared Camera and Spectrometer (COMICS) on the 8.2m Subaru Telescope. We discovered diffuse emission of an elliptical shape extended in the north-south direction inabout 1000AU radius around unresolved excess emission. The diffuse emission is perpendicular to the cavity wall formed by the past outflow activity and is parallel to the projected major axis of the central close binary orbit. The centers of the ellipse contours of the diffuse emission are shifted from the stellar position and the amount of the shift increases as the contour brightness level decreases. The diffuse emission is well explained in all of geometry, size, and configuration by an inclined flared disk where only its surface emits the mid-infrared photons. Our results give the first well-resolved infrared disk images around a massive star and strongly support that HD200775 is formed through the disk accretion. The disk survives the main accretion phase and shows a structure similar to that around lower-mass stars with 'disk atmosphere'. At the same time, the disk also shows properties characteristic to massive stars such as photoevaporation traced by the 3.4mm free-free emission and unusual silicate emission with a peak at 9.2micron, which is shorter than that of many astronomical objects. It provides a good place to compare the disk properties between massive and lower-mass stars.

Massive young stellar objects (MYSO) are surrounded by massive dusty envelopes. Our aim is to establish their density structure on scales of ~1000 AU, i.e. a factor 10 increase in angular resolution compared to similar studies performed in the (sub)mm. We have obtained diffraction-limited (0.6") 24.5 micron images of 14 well-known massive star formation regions with Subaru/COMICS. The images reveal the presence of discrete MYSO sources which are resolved on arcsecond scales. For many sources, radiative transfer models are capable of satisfactorily reproducing the observations. They are described by density powerlaw distributions (n(r) ~ r^(-p)) with p = 1.0 +/-0.25. Such distributions are shallower than those found on larger scales probed with single-dish (sub)mm studies. Other sources have density laws that are shallower/steeper than p = 1.0 and there is evidence that these MYSOs are viewed near edge-on or near face-on, respectively. The images also reveal a diffuse component tracing somewhat larger scale structures, particularly visible in the regions S140, AFGL 2136, IRAS 20126+4104, Mon R2, and Cep A. We thus find a flattening of the MYSO envelope density law going from ~10 000 AU down to scales of ~1000 AU. We propose that this may be evidence of rotational support of the envelope (abridged).

We present diffraction limited (0.6") 24.5micron Subaru/COMICS images of the red supergiant mu Cep. We report the detection of a circumstellar nebula, that was not detected at shorter wavelengths. It extends to a radius of at least 6" in the thermal infrared. On these angular scales, the nebula is roughly spherical, in contrast, it displays a pronounced asymmetric morphology closer in. We simultaneously model the azimuthally averaged intensity profile of the nebula and the observed spectral energy distribution using spherical dust radiative transfer models. The models indicate a constant mass-loss process over the past 1000 years, for mass-loss rates a few times 10^(-7) Msun/yr. This work supports the idea that at least part of the asymmetries in shells of evolved massive stars and supernovae may be due to the mass-loss process in the red supergiant phase.

Hierarchical turbulent structure constituting a jet is considered to reproduce energy-dependent variability in blazars, particularly, the correlation between X- and gamma-ray light curves measured in the TeV blazar Markarian 421. The scale-invariant filaments are featured by the ordered magnetic fields that involve hydromagnetic fluctuations serving as electron scatterers for diffusive shock acceleration, and the spatial size scales are identified with the local maximum electron energies, which are reflected in the synchrotron spectral energy distribution (SED) above the near-infrared/optical break. The structural transition of filaments is found to be responsible for the observed change of spectral hysteresis.

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.

Strong, anisotropic turbulence reflecting magnetized filaments is considered, to model the diffusive acceleration of particles by shock waves in active galactic nucleus jets. We address that at knot A of the nearby M87 jet, the shock involving the filamentary turbulence can accelerate an iron nucleus to zetta-eV (ZeV; 10^{21} eV) ranges. A smaller value of the particle diffusion coefficient is found to be essential to achieve a ZeV cosmic-ray accelerator, a "Zevatron."

This paper presents the effects of electron-positron pair production on the linear growth of the resistive hose instability of a filamentary beam that could lead to snake-like distortion. For both the rectangular radial density profile and the diffuse profile reflecting the Bennett-type equilibrium for a self-collimating flow, the modified eigenvalue equations are derived from a Vlasov-Maxwell equation. While for the simple rectangular profile, current perturbation is localized at the sharp radial edge, for the realistic Bennett profile with an obscure edge, it is non-locally distributed over the entire beam, removing catastrophic wave-particle resonance. The pair production effects likely decrease the betatron frequency, and expand the beam radius to increase the resistive decay time of the perturbed current; these also lead to a reduction of the growth rate. It is shown that, for the Bennett profile case, the characteristic growth distance for a preferential mode can exceed the observational length-scale of astrophysical jets. This might provide the key to the problem of the stabilized transport of the astrophysical jets including extragalactic jets up to Mpc ($\sim 3\times 10^{24}$ cm) scales.

Mid-infrared (MIR) images of the Herbig Ae star HD 142527 were obtained at 18.8 and 24.5 micron with the Subaru/COMICS. Bright extended arc-like emission (outer disk) is recognized at r=0.85" together with a strong central source (inner disk) and a gap around r=0.6" in the both images. Thermal emission of the eastern side is much brighter than that of the western side in the MIR. We estimate the dust size as a few micron from the observed color of the extended emission and the distance from the star. The dust temperature T and the optical depth tau of the MIR emitting dust are also derived from the two images as T=82+/-1K, tau=0.052+/-0.001 for the eastern side and T=85+/-3K, tau=0.018+/-0.001 for the western side. The observed asymmetry in the brightness can be attributed to the difference in the optical depth of the MIR emitting dust. To account for the present observations, we propose an inclined disk model, in which the outer disk is inclined along the east-west direction with the eastern side being in the far side and the inner rim of the outer disk in the eatern side is exposed directly to us. The proposed model can successfully account for the MIR observations as well as near-infrared (NIR) images of the scattering light, in which the asymmetry is seen in the opposite sense and the forward scattering light (near side -- western side) is brighter.

Using the ``modified DPMJET-III'' model explained in the previous paper, we calculate the atmospheric neutrino flux. The calculation scheme is almost the same as HKKM04 \cite{HKKM2004}, but the usage of the ``virtual detector'' is improved to reduce the error due to it. Then we study the uncertainty of the calculated atmospheric neutrino flux summarizing the uncertainties of individual components of the simulation. The uncertainty of $K$-production in the interaction model is estimated by modifying FLUKA'97 and Fritiof 7.02 so that they also reproduce the atmospheric muon flux data correctly, and the calculation of the atmospheric neutrino flux with those modified interaction models. The uncertainties of the flux ratio and zenith angle dependence of the atmospheric neutrino flux are also studied.

We have studied the hadronic interaction for the calculation of the atmospheric neutrino flux by summarizing the accurately measured atmospheric muon flux data and comparing with simulations. We find the atmospheric muon and neutrino fluxes respond to errors in the $\pi$-production of the hadronic interaction similarly, and compare the atmospheric muon flux calculated using the HKKM04 code with experimental measurements. The $\mu^++\mu^-$ data show good agreement in the 1$\sim$30 GeV/c range, but a large disagreement above 30 GeV/c. The the $\mu^+/\mu^-$ ratio shows sizable differences at lower and higher momenta for opposite directions. As the disagreements are considered to be due to assumptions in the hadronic interaction model, we try to improve it phenomenologically based on the quark parton model. The improved interaction model reproduces the observed muon flux data well. The calculation of the atmospheric neutrino flux will be reported in the following paper

We present the results of mid-infrared N-band spectroscopy of the Herbig Ae/Be system MWC1080 using the Cooled Mid-Infrared Camera and Spectrometer (COMICS) on board the 8m Subaru Telescope. The MWC1080 has a geometry such that the diffuse nebulous structures surround the central Herbig B0 type star. We focus on the properties of polycyclic aromatic hydrocarbons (PAHs) and PAH-like species, which are thought to be the carriers of the unidentified infrared (UIR) bands in such environments. A series of UIR bands at 8.6, 11.0, 11.2, and 12.7um is detected throughout the system and we find a clear increase in the UIR 11.0um/11.2um ratio in the vicinity of the central star. Since the UIR 11.0um feature is attributed to a solo-CH out-of-plane wagging mode of cationic PAHs while the UIR 11.2um feature to a solo-CH out-of-plane bending mode of neutral PAHs, the large 11.0um/11.2um ratio directly indicates a promotion of the ionization of PAHs near the central star.

We propose a modified version of the X-ray spectral index and an intrinsic cutoff frequency of inverse Compton radiation from the brightest knot of the M87 jet, in conjunction with an application of the new conceptions of injection and diffusive shock acceleration (DSA) of electrons in magnetized filamentary plasma to the specified source. The drop of the X-ray flux density in a transitive frequency region is associated with the interplay of ordinary synchrotron cooling and weaker magnetic fields concomitant with the smaller scale filaments that allow the electron injection, while the radio-optical synchrotron continuum is dominantly established by the major electrons that are quasi-secularly bound to larger filaments. With reference to, particularly, the updated external Compton model, we demonstrate that in the Klein-Nishina regime fading inverse Comptonization, the injected electrons can be stochastically energized up to a Lorentz factor as high as $5\times 10^{10}$ in the temporal competition with diffuse synchrotron cooling; this value is larger than that attainable for a simple DSA scenario based on the resonant scattering diffusion of the gyrating electrons bound to a supposed magnetic field homogeneously pervading the entire knot. The upper limits of the photon frequency boosted via conceivable inverse Compton processes are predicted to be of the common order of $\sim 10^{30}$ Hz. The variability of the broadband spectrum is also discussed in comparison to the features of a blazar light curve. The present scenario of a peta-eV (PeV; $10^{15}$ eV) electron accelerator, the "Pevatron," might provide some guidance for exploring untrod hard X-ray and gamma-ray bands in forthcoming observations.

We expand the off-resonant scattering theory for particle diffusion in magnetized current filaments that can be typically compared to astrophysical jets, including active galactic nucleus jets. In a high plasma beta region where the directional bulk flow is a free-energy source for establishing turbulent magnetic fields via current filamentation instabilities, a novel version of quasi-linear theory to describe the diffusion of test particles is proposed. The theory relies on the proviso that the injected energetic particles are not trapped in the small-scale structure of magnetic fields wrapping around and permeating a filament but deflected by the filaments, to open a new regime of the energy hierarchy mediated by a transition compared to the particle injection. The diffusion coefficient derived from a quasi-linear type equation is applied to estimating the timescale for the stochastic acceleration of particles by the shock wave propagating through the jet. The generic scalings of the achievable highest energy of an accelerated ion and electron, as well as of the characteristic time for conceivable energy restrictions, are systematically presented. We also discuss a feasible method of verifying the theoretical predictions. The strong, anisotropic turbulence reflecting cosmic filaments might be the key to the problem of the acceleration mechanism of the highest energy cosmic rays exceeding 100 EeV (10^{20} eV), detected in recent air shower experiments.

Cosmic ray particles are more rapidly accelerated in oblique shocks, with the magnetic field inclined with respect to the shock normal direction, than in parallel shocks, as a result of mirror reflection at the shock surface and slower diffusion in the shock normal direction. We investigate quantitatively how these effects contribute to reducing the acceleration time over the whole range of magnetic field inclinations. It is shown that, for quasi-perpendicular inclination, the mirror effect plays a remarkable role in reducing the acceleration time; whereas, at relatively small inclination, the anisotropic diffusion effect is dominant in reducing that time. These results are important for a detailed understanding of the mechanism of particle acceleration by an oblique shock in space and heliosphereic plasmas.

We estimate the highest energy of proton diffusively accelerated by shock in knot A1 of the jet in luminous nearby quasar 3C 273. Referring to the recent polarization measurements using very long baseline interferometry (VLBI), we consider the shock propagation across magnetic field lines, namely, configuration of the oblique shock. For larger inclination of the field lines, the effects of particle reflection at the shock front are more pronounced, to significantly increase acceleration efficiency. The quasiperpendicular shock turns out to be needed for safely achieving the proton acceleration to the energy above 100 EeV (10^20 eV) in a parameter domain reflecting conceivable energy restrictions.

We have observed the 8-13 $\mu$m spectrum (R$\sim$250) of the Vega-like star candidate HD145263 using Subaru/COMICS. The spectrum of HD145263 shows the broad trapezoidal silicate feature with the shoulders at 9.3 $\mu$m and 11.44 $\mu$m, indicating the presence of crystalline silicate grains. This detection implies that crystalline silicate may also be commonly present around Vega-like stars. The 11.44 $\mu$m feature is slightly shifted to a longer wavelength compared to the usual 11.2-3 $\mu$m crystalline forsterite feature detected toward Herbig Ae/Be stars and T Tauri stars. Although the peak shift due to the effects of the grain size can not be ruled out, we suggest that Fe-bearing crystalline olivine explains the observed peak wavelength fairly well. Fe-bearing silicates are commonly found in meteorites and most interplanetary dust particles, which originate from planetesimal-like asteroids. According to studies of meteorites, Fe-bearing silicate must have been formed in asteroidal planetesimals, supporting the scenario that dust grains around Vega-like stars are of planetesimal origin, if the observed 11.44 $\mu$m peak is due to Fe-bearing silicates.

We have revised the calculation of the flux of atmospheric neutrinos based on a 3-dimensional scheme with the realistic IGRF geomagnetic model. The primary flux model has been revised, based on the AMS and BESS observations, and the interaction model updated to DPMJET-III. With a fast simulation code and computer system, the statistical errors in the Monte Carlo study are negligible. We estimate the total uncertainty of the atmospheric neutrino flux prediction is reduced to $\lesssim$ 10 % below 10 GeV. The `3-dimensional effects' are found to be almost the same as the study with the dipole magnetic field, but the muon curvature effect remains up to a few 10 GeV for horizontal directions. The uncertainty of the absolute normalization of the atmospheric neutrino is still large above 10 GeV due to the uncertainty of the primary cosmic ray flux above 100 GeV. However, the zenith angle variation is not affected by these uncertainties.

We report centimeter continuum observations with the VLA and the VLBA as well as mid-infrared observations with COMICS/SUBARU toward the components of the YLW 15 very young binary system, VLA 1 and VLA 2. The centimeter emission of the two components traces partially thick free-free emission, likely due to collimated, ionized winds. VLA 1 is an embedded protostar, undetected in the near-IR, and possibly in the Class 0 to Class I transition and powering a Herbig-Haro outflow. Its mid-IR emission appears slightly resolved with a diameter of ~16 AU, possibly tracing circumstellar material from both the envelope and the disk. VLA 2 is a typical Class I object, unresolved in the mid-IR, and is the responsible of the strong X-ray emission associated with YLW 15. The expected centimeter ''peri-stellar'' emission associated with the X-ray emission is not detected with the VLBA at 6 cm likely due to the high optical depth of the free-free emission. Strikingly, the near to mid-IR properties of YLW 15 suggest that VLA 1 is a more embedded YSO, or alternatively, less luminous than VLA 2, whereas orbital proper motions of this binary system indicate that VLA 1 is more massive than VLA 2. This result is apparently against the expected evolutionary scenario, where one expects that the more massive YSO in a binary system is the more evolved and luminous YSO. Finally, the nearby source YLW 16A is detected with the VLA, its position coincides with reported near-IR and X-ray sources.

Precise knowledge of the hadronic interaction between primary cosmic rays and atmospheric nuclei is very important and fundamental to study atmospheric neutrinos and their oscillations. We studied hadronic interaction models using the data of primary and secondary cosmic rays observed by BESS experiments. By comparing the observed spectra with the ones calculated by several interaction models, we found DPMJET-III is most favored among them.

I theoretically found eigenmodes and growth rates of relativistic current filamentation instability in collisional regimes, deriving a generalized dispersion relation from self-consistent beam-Maxwell equations. For symmetrically counterstreaming, fully relativistic electron currents, the collisional coupling between electrons and ions creates the unstable modes of growing oscillation and wave, which stand out for long-wavelength perturbations. In the stronger collisional regime, the growing oscillatory mode tends to be dominant for all wavelengths. In the collisionless limit, those modes vanish, while maintaining another purely growing mode that exactly coincides with a standard relativistic Weibel mode. It is also shown that the effects of electron-electron collisions and thermal spread lower the growth rate of the relativistic Weibel instability. The present mechanisms of filamentation dynamics are essential for transport of homogeneous electron beam produced by the interaction of high power laser pulses with plasma.

A cosmic ray detection system, consisting of standardized detector stations connected through the Internet, is described. The system can be used for detecting air showers that arrive over a wide area with correlated time. The data at each site are exchanged in (quasi) real time, which makes the system suitable for displaying cosmic ray arrivals at museumns and schools.

We propose a novel mathematical method to construct an exact polytropic sphere in self-gravitating hydrostatic equilibrium, improving the non-linear Poisson equation. The central boundary condition for the present equation requires a ratio of gas pressure to total one at the centre, which is uniquely identified by the whole mass and molecular weight of the system. The special solution derived from the Lane-Emden equation can be reproduced. This scheme is now available for modelling the molecular cloud cores in interstellar media. The mass-radius relation of the first core is found to be consistent with the recent results of radiation hydrodynamic simulations.

We analyze the data of the gravitational microlensing survey carried out by by the MOA group during 2000 towards the Galactic Bulge (GB). Our observations are designed to detect efficiently high magnification events with faint source stars and short timescale events, by increasing the the sampling rate up to 6 times per night and using Difference Image Analysis (DIA). We detect 28 microlensing candidates in 12 GB fields corresponding to 16 deg^2. We use Monte Carlo simulations to estimate our microlensing event detection efficiency, where we construct the I-band extinction map of our GB fields in order to find dereddened magnitudes. We find a systematic bias and large uncertainty in the measured value of the timescale $t_{\rm Eout}$ in our simulations. They are associated with blending and unresolved sources, and are allowed for in our measurements. We compute an optical depth tau = 2.59_{-0.64}^{+0.84} \times 10^{-6} towards the GB for events with timescales 0.3<t_E<200 days. We consider disk-disk lensing, and obtain an optical depth tau_{bulge} = 3.36_{-0.81}^{+1.11} \times 10^{-6}[0.77/(1-f_{disk})] for the bulge component assuming a 23% stellar contribution from disk stars. These observed optical depths are consistent with previous measurements by the MACHO and OGLE groups, and still higher than those predicted by existing Galactic models. We present the timescale distribution of the observed events, and find there are no significant short events of a few days, in spite of our high detection efficiency for short timescale events down to t_E = 0.3 days. We find that half of all our detected events have high magnification (>10). These events are useful for studies of extra-solar planets.

We have carried out mid-infrared N-band spectroscopic observations of the T Tauri star Hen 3-600A in the TW Hydra association with the COMICS on the 8.2m Subaru Telescope and found structured features in its spectrum. These structured features are well explained by a combination of crystalline forsterite, crystalline enstatite, silica and glassy olivine grains. Among intermediate-mass young stellar objects (YSOs), crystalline silicates have already been detected, but no firm detection has been reported so far for low-mass YSOs such as T Tauri stars. This is the first clear detection of crystalline silicates in low-mass YSOs and shows that the crystallization event occurs even in the protoplanetary disk of low-mass YSOs in the T Tauri phase. The physical processes leading to the inferred dust composition in the Hen3-600A system may be analogous to those occured in the early epoch of the Solar system.

A novel model for collimation and transport of electron-positron-ion jets is presented. Analytical results show that the filamentary structures can be sustained by self-induced toroidal magnetic fields permeating through the filaments, whose widths significantly expand in the pair-dominant regimes. The magnetic field strength reflects a characteristic of equipartition of excess kinetic energy of the jets. It is also shown that growth of the hoselike instability is strongly suppressed. Essential features derived from this model are consistent with recent results observed by using very long baseline telescopes.

A search for extra-solar planets was carried out in three gravitational microlensing events of high magnification, MACHO 98-BLG-35, MACHO 99-LMC-2, and OGLE 00-BUL-12. Photometry was derived from observational images by the MOA and OGLE groups using an image subtraction technique. For MACHO 98-BLG-35, additional photometry derived from the MPS and PLANET groups was included. Planetary modeling of the three events was carried out in a super-cluster computing environment. The estimated probability for explaining the data on MACHO 98-BLG-35 without a planet is <1%. The best planetary model has a planet of mass ~(0.4-1.5) X 10^-5 M_Earth at a projected radius of either ~1.5 or ~2.3 AU. We show how multi-planet models can be applied to the data. We calculated exclusion regions for the three events and found that Jupiter-mass planets can be excluded with projected radii from as wide as about 30 AU to as close as around 0.5 AU for MACHO 98-BLG-35 and OGLE 00-BUL-12. For MACHO 99-LMC-2, the exclusion region extends out to around 10 AU and constitutes the first limit placed on a planetary companion to an extragalactic star. We derive a particularly high peak magnification of ~160 for OGLE 00-BUL-12. We discuss the detectability of planets with masses as low as Mercury in this and similar events.

One hundred and forty six long-period red variable stars in the Large Magellanic Cloud (LMC) from the three year MOA project database were analysed. A careful periodic analysis was performed on these stars and a catalogue of their magnitudes, colours, periods and amplitudes is presented. We convert our blue and red magnitudes to $K$ band values using 19 oxygen-rich stars. A group of red short-period stars separated from the Mira sequence has been found on a (log P, K) diagram. They are located at the short period side of the Mira sequence consistent with the work of Wood and Sebo (1996). There are two interpretations for such stars; a difference in pulsation mode or a difference in chemical composition. We investigated the properties of these stars together with their colour, amplitude and periodicity. We conclude that they have small amplitudes and less regular variability. They are likely to be higher mode pulsators. A large scatter has been also found on the long period side of the (log P, K) diagram. This is possibly a systematic spread given that the blue band of our photometric system covers both standard B and V bands and affects carbon-rich stars.

We describe observations carried out by the MOA group of the Galactic Bulge during 2000 that were designed to detect efficiently gravitational microlensing of faint stars in which the magnification is high and/or of short duration. These events are particularly useful for studies of extra-solar planets and faint stars. Approximately 17 degrees square were monitored at a sampling rate of up to 6 times per night. The images were analysed in real-time using a difference imaging technique. Twenty microlensing candidates were detected, of which 8 were alerted to the microlensing community whilst in progress. Approximately half of the candidates had high magnifications (>~10), at least one had very high magnification (>~50), and one exhibited a clear parallax effect. The details of these events are reported here, together with details of the on-line difference imaging technique. Some nova-like events were also observed and these are described, together with one asteroid.

Recent work by the MOA gravitational microlensing group is briefly described, including (i) the current observing strategy, (ii) use of a high-speed parallel computer for analysis of results by inverse ray shooting, (iii) analysis of the light curve of event OGLE-2000-BUL12 in terms of extra-solar planets, and (iv) the MOA alert system using difference imaging.

We present observations of microlensing event MACHO-98-BLG-35 which reached a peak magnification factor of almost 80. These observations by the Microlensing Planet Search (MPS) and the MOA Collaborations place strong constraints on the possible planetary system of the lens star and show intriguing evidence for a low mass planet with a mass fraction $4\times 10^{-5} \leq \epsilon \leq 2\times 10^{-4}$. A giant planet with $\epsilon = 10^{-3}$ is excluded from 95% of the region between 0.4 and 2.5 $R_E$ from the lens star, where $R_E$ is the Einstein ring radius of the lens. This exclusion region is more extensive than the generic "lensing zone" which is $0.6 - 1.6 R_E$. For smaller mass planets, we can exclude 57% of the "lensing zone" for $\epsilon = 10^{-4}$ and 14% of the lensing zone for $\epsilon = 10^{-5}$. The mass fraction $\epsilon = 10^{-5}$ corresponds to an Earth mass planet for a lensing star of mass $\sim 0.3 \msun$. A number of similar events will provide statistically significant constraints on the prevalence of Earth mass planets. In order to put our limits in more familiar terms, we have compared our results to those expected for a Solar System clone averaging over possible lens system distances and orientations. We find that such a system is ruled out at the 90% confidence level. A copy of the Solar System with Jupiter replaced by a second Saturn mass planet can be ruled out at 70% confidence. Our low mass planetary signal (few Earth masses to Neptune mass) is significant at the $4.5\sigma$ confidence level. If this planetary interpretation is correct, the MACHO-98-BLG-35 lens system constitutes the first detection of a low mass planet orbiting an ordinary star without gas giant planets.

A review of the MOA (Microlensing Observations in Astrophysics) project is presented. MOA is a collaboration of approximately 30 astronomers from New Zealand and Japan established with the aim of finding and detecting microlensing events towards the Magellanic Clouds and the Galactic bulge, which may be indicative of either dark matter or of planetary companions. The observing program commenced in 1995, using very wide band blue and red filters and a nine-chip mosaic CCD camera. As a by-product of these observations a large database of CCD photometry for 1.4 million stars towards both LMC and SMC has been established. In one preliminary analysis 576 bright variable stars were confirmed, nearly half of them being Cepheids. Another analysis has identified large numbers of blue variables, and 205 eclipsing binaries are included in this sample. In addition 351 red variables (AGB stars) have been found. Light curves have been obtained for all these stars. The observations are carried out on a 61-cm f/6.25 telescope at Mt John University Observatory where a new larger CCD camera was installed in 1998 July. From this latitude (44 S) the Magellanic Clouds can be monitored throughout the year.

We present analysis of MACHO Alert 95-30, a dramatic gravitational microlensing event towards the Galactic bulge whose peak magnification departs significantly from the standard point-source microlensing model. Alert 95-30 was observed in real-time by the Global Microlensing Alert Network (GMAN), which obtained densely sampled photometric and spectroscopic data throughout the event. We interpret the light-curve ``fine structure'' as indicating transit of the lens across the extended face of the source star. This signifies resolution of a star several kpc distant. We find a lens angular impact parameter theta_{min}/theta_{source} = 0.715 +/- 0.003. This information, along with the radius and distance of the source, provides an additional constraint on the lensing system. Spectroscopic and photometric data indicate the source is an M4 III star of radius 61 +/- 12 Rsun, located on the far side of the bulge at 9 kpc. We derive a lens angular velocity, relative to the source, of 21.5 +/- 4.9 km/s/kpc, where the error is dominated by uncertainty in the source radius. Likelihood analysis yields a median lens mass of 0.67{+2.53}{-0.46} Msun, located with 80% probability in the Galactic bulge at a distance of 6.93{+1.56}{-2.25} kpc. If the lens is a main-sequence star, we can include constraints on the lens luminosity. This modifies our estimates to M_lens = 0.53{+0.52}{-0.35} Msun and D_lens = 6.57{+0.99}{-2.25} kpc. Spectra taken during the event show that the absorption line equivalent widths of H alpha and the TiO bands near 6700 A vary, as predicted for microlensing of an extended source. This is most likely due to center-to-limb variation in the stellar spectral lines. These data demonstrate the feasibility of using microlensing limb crossings as a tool to probe stellar atmospheres directly.