Search

In this era of spatially resolved observations of planet forming disks with ALMA and large ground-based telescopes such as the VLT, Keck and Subaru, we still lack statistically relevant information on the quantity and composition of the material that is building the planets, such as the total disk gas mass, the ice content of dust, and the state of water in planetesimals. SPICA is an infrared space mission concept developed jointly by JAXA and ESA to address these questions. The key unique capabilities of SPICA that enable this research are (1) the wide spectral coverage 10-220 micron, (2) the high line detection sensitivity of (1-2) 10-19 W m-2 with R~2000-5000 in the far-IR (SAFARI) and 10-20 W m-2 with R~29000 in the mid-IR (SMI, spectrally resolving line profiles), (3) the high far-IR continuum sensitivity of 0.45 mJy (SAFARI), and (4) the observing efficiency for point source surveys. This paper details how mid- to far-IR infrared spectra will be unique in measuring the gas masses and water/ice content of disks and how these quantities evolve during the planet forming period. These observations will clarify the crucial transition when disks exhaust their primordial gas and further planet formation requires secondary gas produced from planetesimals. The high spectral resolution mid-IR is also unique for determining the location of the snowline dividing the rocky and icy mass reservoirs within the disk and how the divide evolves during the build-up of planetary systems. Infrared spectroscopy (mid- to far-IR) of key solid state bands is crucial for assessing whether extensive radial mixing, which is part of our Solar System history, is a general process occurring in most planetary systems and whether extrasolar planetesimals are similar to our Solar System comets/asteroids. ... (abbreviated), Comment: accepted for publication in PASA

In JT-60U, simultaneous excitation of n = 1 helical cores (HCs) and m/n = 2/1 Tearing Modes (TMs) was observed [T. Bando et al., Plasma Phys. Control. Fusion 61 115014 (2019)]. In this paper, we have investigated the excitation mechanism of n = 1 HCs with m/n = 2/1 TMs based on the experimental observations and a simple quasi-linear MHD model. In the previous study, it was reported that a "coupling" on the phase of the MHD mode is observed between n = 1 HCs and m/n = 2/1 TMs. In this study, it is found that the coupling is observed with the mode frequency from several Hz to 6 kHz. This indicates that the resistive wall and the plasma control system do not induce the coupling because the both time scales are different from the mode frequency. In addition, n = 1 HCs appear to be the non-resonant mode from the two observations: n = 1 HCs do not rotate with the plasma around the q = 1 surface in the core and the coupling is also observed even when qmin > 1. It is also observed that the electron fluctuation due to an n = 1 HC in the core region disappears with the stabilization of an m/n = 2/1 neoclassical tearing mode by electron cyclotron current drive, implying that n = 1 HCs are driven by m/n = 2/1 TMs. This perspective, n = 1 HCs are driven by m/n = 2/1 TMs, is supported by the observation that the saturated amplitude of the m/n = 1/1 component of the radial displacement in the core is smaller than that of the m/n = 2/1 component. Finally, we revisit a quasi-linear MHD model where the m/n = 1/1 HC is induced directly by the sideband of the current for the m/n = 2/1 TM, which allows to excite the non-resonant m/n = 1/1 mode. The model also describes the characteristic of the coupling, fm/n=1/1(HC) = 2fm/n=2/1(TM).

We analyze the effects of the passing energetic particles on the resistive ballooning modes (RBM) and the energetic particle driven modes in JT-60SA plasma, which leads to the prediction of the stability in N-NBI heated plasma. The analysis is performed using the code FAR3d that solves the reduced MHD equations describing the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moments for the energetic particle (EP) species assuming an averaged Maxwellian EP distribution fitted to the slowing down distribution, including the effect of the acoustic modes. The simulations show the possible destabilization of a 3/2-4/2 TAE with a frequency (f) of 115 kHz, a 6/4-7/4 TAE with f=98 kHz and a 6/4 or 7/4 BAE with f=57 kHz in the ITER-like inductive scenario. If the energetic particle beta increases, Beta induced AEs (BAE), Toroidal AEs (TAE) and Elliptical AEs (EAE) are destabilized between the inner-middle plasma region, leading to the overlapping of AE of different toroidal families. If these instabilities coexist in the non-linear saturation phase the EP transport could be enhanced leading to a lower heating efficiency. For a hypothetical configuration based on the ITER-like inductive scenario but an center peaked EP profile, the EP beta threshold increases and several BAEs are destabilized in the inner plasma region, indicating an improved AE stability with respect to the off-axis peaked EP profile. In addition, the analysis of a hypothetical JT-60SA scenario with a resonant q=1 in the inner plasma region shows the destabilization of fishbones-like instabilities by the off-axis peaked EP profile. Also, the EPs have an stabilizing effect on the RBM, stronger as the population of EP with low energies (below 250 keV) increases at the plasma pedestal.

A unific helical field as geometric structure in an internal space is proposed which has the ability to serve attributes inherent to point particles. Transition between our observational space and the warped infinitesimal space is considered in an operational fashion. It is seen, that rotational eigenvalue equation satisfied by the vector field equivalent to Gromeka-Beltrami flow familiar to plasma astrophysics provides a spatio-unifier that sustains complex orthogonal coupling between rotational, internal coordinate space and angular momentum space. Self-consistent normalization of the rotational coordinate owing to the unifier is compared to renormalization of quantum electrodynamics. A numerical value comparable to the electromagnetic coupling constant in the nonrelativistic limit is derived from theoretical analysis involving the rotational eigenvalue that charged leptons should internally refer to. It is found that chiral asymmetry of the helical eigenflows can be reflected in electroweak symmetry breaking. This study might provide a key to settle geometrical problem of preceding theories accommodated by extra-dimensions of space., Comment: 19 pages; major revisions, results unchanged

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., Comment: 14 pages, 30 fugures

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., Comment: accepted in ApJ

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., Comment: 7 figures, Accepted by ApJ

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., Comment: 25 pages, 16 figures

*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., Comment: 20 pages, 42 Figures

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., Comment: 12Pages,9Figs

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., Comment: accepted by APJL

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., Comment: 20 pages, 8 figures, and 1 table. accepted to PASJ

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., Comment: 18 pages, 5 figures, accepted to ApJ

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., Comment: 41 pages (single column), 1 table, 10 figures, ApJ accepted

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., Comment: 8 pages, 5 figures, ApJL in press, typo corrected

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., Comment: Accepted in MNRAS

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., Comment: 22 pages, 10 figures

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., Comment: 12 pages, 3 figures

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. ...., Comment: 13 pages, 17 figures. Accepted by A&A 14 January 2011

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., Comment: 17 pages, 3 Postscript figures

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., Comment: 10 pages, 4 figures, accepted for publication in ApJ

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., Comment: 5 pages, 4 figures, accepted for publication in ApJ Letters. Minor change to Figure 4

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., Comment: 18 pages, 8 figures, accepted by The Astrophysical Journal

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)., Comment: 21 pages, accepted for A&A

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., Comment: Accepted for publication in ApJ Letters

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

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., Comment: 6 pages, 2 figures

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., Comment: 31 pages, 19 figures, submitted to Phys. Rev. D

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., Comment: 30 pages, 6 figures. Accepted in Astronomical Journal

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., Comment: 4 pages, 5 figures, a revise of the previous work, the proceeding of ICRC2003

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., Comment: 10 pages, 9 figures

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.310). These events are useful for studies of extra-solar planets., Comment: 65 pages and 30 figures, accepted for publication in ApJ. A systematic bias and uncertainty in the optical depth measurement has been quantified by simulations

We observed atmospheric gamma-rays around 10 GeV at balloon altitudes (15~25 km) and at a mountain (2770 m a.s.l). The observed results were compared with Monte Carlo calculations to find that an interaction model (Lund Fritiof1.6) used in an old neutrino flux calculation was not good enough for describing the observed values. In stead, we found that two other nuclear interaction models, Lund Fritiof7.02 and dpmjet3.03, gave much better agreement with the observations. Our data will serve for examining nuclear interaction models and for deriving a reliable absolute atmospheric neutrino flux in the GeV region., Comment: 9 pages 12 figures

Atmospheric neutrinos produced by cosmic-ray interactions in the atmosphere are of interest for several reasons. As a beam for studies of neutrino oscillations they cover a range of parameter space hitherto unexplored by accelerator neutrino beams. The atmospheric neutrinos also constitute an important background and calibration beam for neutrino astronomy and for the search for proton decay and other rare processes. Here we review the literature on calculations of atmospheric neutrinos over the full range of energy, but with particular attention to the aspects important for neutrino oscillations. Our goal is to assess how well the properties of atmospheric neutrinos are known at present., Comment: 68 pages, 26 figures. With permission from the Annual Review of Nuclear & Particle Science. Final version of this material is scheduled to appear in the Annual Review of Nuclear & Particle Science Vol. 52, to be published in December 2002 by Annual Reviews (http://annualreviews.org)

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., Comment: 19 pages, 19 figures, accepted for publication in MNRAS

A 3-dimensional calculation of atmospheric neutrinos flux is presented, and the results are compared with those of a 1-dimensional one. In this study, interaction and propagation of particles is treated in a 3-dimensional way including the curvature of charged particles due to the geomagnetic field, which is assumed to be a dipole field. The purpose of this paper is limited to the comparison of calculation schemes. The updated flux value with new interaction model and primary flux model will be reported in a separate paper. Except for nearly horizontal directions, the flux is very similar to the result of 1 dimensional calculations. However, for near-horizontal directions an enhancement of the neutrino flux is seen even at energies as high as 1 GeV. The production height of neutrinos is lower than the prediction by 1-dimensional calculation for near-horizontal directions, and is a little higher for near-vertical directions. However, the difference is not evident except for near-horizontal directions., Comment: 22 pages, 15figures

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., Comment: 14 pages, 16 embedded postscript figures, 3 PNG figures, revised version accepted by MNRAS

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., Comment: 14 pages, 9 embedded postscript figures, 2 PNG figures, replaced with revised version accepted for publication in MNRAS, full resolution version available at http://www.phys.canterbury.ac.nz/~physib/papers