Dust continuum observation is one of the best methods to constrain the properties of protoplanetary disks. Recent theoretical studies have suggested that the dust scattering at the millimeter wavelength potentially reduces the observed intensity, which results in an underestimate in the dust mass. We investigate whether the dust scattering indeed reduces the observed continuum intensity by comparing the ALMA archival data of the TW Hya disk at Band 3, 4, 6, 7 and 9 to models obtained by radiative transfer simulations. We find that the model with scattering by 300 ${\rm \mu m}$-sized grains well reproduces the observed SED of the central part of the TW Hya disk while the model without scattering is also consistent within the errors of the absolute fluxes. To explain the intensity at Band 3, the dust surface density needs to be $\sim$ 10 ${\rm g\,cm^{-2}}$ at 10 au in the model with scattering, which is 26 times more massive than previously predicted. The model without scattering needs 2.3 times higher dust mass than the model with scattering because it needs lower temperature. At Band 7, scattering reduces the intensity by $\sim$ 35% which makes the disk looks optically thin even though it is optically thick. Our study suggests the TW Hya disk is still capable of forming cores of giant planets at where the current solar system planets exist.

We are developing an ultra-wideband spectroscopic instrument, DESHIMA (DEep Spectroscopic HIgh-redshift MApper), based on the technologies of an on-chip filter-bank and Microwave Kinetic Inductance Detector (MKID) to investigate dusty star-burst galaxies in the distant universe at millimeter and submillimeter wavelength. An on-site experiment of DESHIMA was performed using the ASTE 10-m telescope. We established a responsivity model that converts frequency responses of the MKIDs to line-of-sight brightness temperature. We estimated two parameters of the responsivity model using a set of skydip data taken under various precipitable water vapor (PWV, 0.4-3.0 mm) conditions for each MKID. The line-of-sight brightness temperature of sky is estimated using an atmospheric transmission model and the PWVs. As a result, we obtain an average temperature calibration uncertainty of $1\sigma=4$%, which is smaller than other photometric biases. In addition, the average forward efficiency of 0.88 in our responsivity model is consistent with the value expected from the geometrical support structure of the telescope. We also estimate line-of-sight PWVs of each skydip observation using the frequency response of MKIDs, and confirm the consistency with PWVs reported by the Atacama Large Millimeter/submillimeter Array.

Ultra-wideband 3D imaging spectrometry in the millimeter-submillimeter (mm-submm) band is an essential tool for uncovering the dust-enshrouded portion of the cosmic history of star formation and galaxy evolution. However, it is challenging to scale up conventional coherent heterodyne receivers or free-space diffraction techniques to sufficient bandwidths ($\geq$1 octave) and numbers of spatial pixels (>$10^2$). Here we present the design and first astronomical spectra of an intrinsically scalable, integrated superconducting spectrometer, which covers 332-377 GHz with a spectral resolution of $F/\Delta F \sim 380$. It combines the multiplexing advantage of microwave kinetic inductance detectors (MKIDs) with planar superconducting filters for dispersing the signal in a single, small superconducting integrated circuit. We demonstrate the two key applications for an instrument of this type: as an efficient redshift machine, and as a fast multi-line spectral mapper of extended areas. The line detection sensitivity is in excellent agreement with the instrument design and laboratory performance, reaching the atmospheric foreground photon noise limit on sky. The design can be scaled to bandwidths in excess of an octave, spectral resolution up to a few thousand and frequencies up to $\sim$1.1 THz. The miniature chip footprint of a few $\mathrm{cm^2}$ allows for compact multi-pixel spectral imagers, which would enable spectroscopic direct imaging and large volume spectroscopic surveys that are several orders of magnitude faster than what is currently possible.

We perform simulations of the dust and gas disk evolution to investigate the observational features of a dust-pileup at the dead-zone inner edge. We show that the total mass of accumulated dust particles is sensitive to the turbulence strength in the dead zone, $\alpha_{\rm dead}$, because of the combined effect of turbulence-induced particle fragmentation (which suppresses particle radial drift) and turbulent diffusion. For a typical critical fragmentation velocity of silicate dust particles of $1~{\rm m~s^{-1}}$, the stress to pressure ratio $\alpha_{\rm dead}$ needs to be lower than $3 \times 10^{-4}$ for dust trapping to operate. The obtained dust distribution is postprocessed using the radiative transfer code RADMC-3D to simulate infrared scattered-light images of the inner part of protoplanetary disks with a dust pileup. We find that a dust pileup at the dead-zone inner edge, if present, casts a shadow extending out to $\sim 10~{\rm au}$. In the shadowed region the temperature significantly drops, which in some cases yields even multiple water snow lines. We also find that even without a dust pileup at the dead-zone inner edge, the disk surface can become thermally unstable, and the excited waves can naturally produce shadows and ring-like structures in observed images. This mechanism might account for the ring-like structures seen in the scattered-light images of some disks, such as the TW Hya disk.

The total amount of dust (or "metallicity") and the dust distribution in protoplanetary disks are crucial for planet formation. Dust grains radially drift due to gas--dust friction, and the gas is affected by the feedback from dust grains. We investigate the effects of the feedback from dust grains on the viscous evolution of the gas, taking into account the vertical dust settling. The feedback from the grains pushes the gas outward. When the grains are small and the dust-to-gas mass ratio is much smaller than unity, the radial drift velocity is reduced by the feedback effect but the gas still drifts inward. When the grains are sufficiently large or piled-up, the feedback is so effective that forces the gas flows outward. Although the dust feedback is affected by dust settling, we found that the 2D approximation reasonably reproduces the vertical averaged flux of gas and dust. We also performed the 2D two-fluid hydrodynamic simulations to examine the effect of the feedback from the grains on the evolution of the gas disk. We show that when the feedback is effective, the gas flows outward and the gas density at the region within $\sim 10\ \mbox{AU}$ is significantly depleted. As a result, the dust-to-gas mass ratio at the inner radii may significantly excess unity, providing the environment where planetesimals are easily formed via, e.g., streaming instability. We also show that a simplified 1D model well reproduces the results of the 2D two-fluid simulations, which would be useful for future studies.

We analytically derive the expressions for the structure of the inner region of protoplanetary disks based on the results from the recent hydrodynamical simulations. The inner part of a disk can be divided into four regions: dust-free region with gas temperature in the optically thin limit, optically thin dust halo, optically thick condensation front and the classical optically thick region in order from the inside. We derive the dust-to-gas mass ratio profile in the dust halo using the fact that partial dust condensation regulates the temperature to the dust evaporation temperature. Beyond the dust halo, there is an optically thick condensation front where all the available silicate gas condenses out. The curvature of the condensation surface is determined by the condition that the surface temperature must be nearly equal to the characteristic temperature $\sim 1200{\,\rm K}$. We derive the mid-plane temperature in the outer two regions using the two-layer approximation with the additional heating by the condensation front for the outermost region. As a result, the overall temperature profile is step-like with steep gradients at the borders between the outer three regions. The borders might act as planet traps where the inward migration of planets due to gravitational interaction with the gas disk stops. The temperature at the border between the two outermost regions coincides with the temperature needed to activate magnetorotational instability, suggesting that the inner edge of the dead zone must lie at this border. The radius of the dead-zone inner edge predicted from our solution is $\sim$ 2-3 times larger than that expected from the classical optically thick temperature.

In the Solar System, interplanetary dust particles (IDPs) originating mainly from asteroid collisions and cometary activities drift to the Earth orbit due to the Poynting-Robertson drag. We analyzed the thermal emission from IDPs that was observed by the first Japanese infrared astronomical satellite, AKARI. The observed surface brightness in the trailing direction of the Earth orbit is 3.7% greater than that in the leading direction in the $9{\rm \mu m}$ band and 3.0% in the $18{\rm \mu m}$ band. In order to reveal dust properties causing the leading-trailing surface brightness asymmetry, we numerically integrated orbits of the Sun, the Earth, and a dust particle as a restricted three-body problem including radiation from the Sun. The initial orbits of particles are determined according to the orbits of main-belt asteroids or Jupiter-family comets. The orbital trapping in mean motion resonances results in a significant leading-trailing asymmetry so that intermediate sized dust (~10-100${\rm \mu m}$) produces a greater asymmetry than the zodiacal light has. The leading-trailing surface brightness difference integrated over the size distribution of the asteroidal dust is obtained to be the values of 27.7% and 25.3% in the $9{\rm \mu m}$ and $18{\rm \mu m}$ bands, respectively. In contrast, the brightness difference for cometary dust is calculated as the values of 3.6% and 3.1% in the $9{\rm \mu m}$ and $18{\rm \mu m}$ bands, respectively, if the maximum dust radius is set to be $s_{\rm max} = 3000{\rm \mu m}$. Taking into account these values and their errors, we conclude that the contribution of asteroidal dust to the zodiacal infrared emission is less than ~10%, while cometary dust of the order of 1 mm mainly accounts for the zodiacal light in infrared.

We describe a novel GaAs/AlGaAs double-quantum-well device for the infrared photon detection, called Charge-Sensitive Infrared Phototransistor (CSIP). The principle of CSIP detector is the photo-excitation of an intersubband transition in a QW as an charge integrating gate and the signal amplification by another QW as a channel with very high gain, which provides us with extremely high responsivity (10^4 -- 10^6 A/W). It has been demonstrated that the CSIP designed for the mid-infrared wavelength (14.7 um) has an excellent sensitivity; the noise equivalent power (NEP) of 7x10^-19 W/rHz with the quantum efficiency of ~2%. Advantages of the CSIP against the other highly sensitive detectors are, huge dynamic range of >10^6, low output impedance of 10^3 -- 10^4 Ohms, and relatively high operation temperature (>2K). We discuss possible applications of the CSIP to FIR photon detection covering 35 -- 60 um waveband, which is a gap uncovered with presently available photoconductors.

The accretion-powered pulsar Her X-1 was observed with Suzaku twice in its main-on state, on 2005 October 5-6 and 2006 March 29-30, for a net exposure of 30.5 ks and 34.4 ks, respectively. In the 2005 and 2006 observations, the source was detected at an average 10-30 keV intensity of 290 mCrab and 230 mCrab, respectively. The intrinsic pulse period was measured on both occasions at 1.23776 s by HXD-PIN, after barycentric and binary corrections. The pulse phase-averaged spectra in the energy range above 10 keV are well fitted by ``Negative and Positive power-law times EXponential (NPEX)'' model, multiplied by a fundamental cyclotron resonance scattering feature at ~36 keV which appears very significantly in the HXD-PIN data. The resonance profiles were reproduced successfully by the Lorentzian type scattering cross section, rather than by a Gaussian type alternative. The pulse phase-averaged HXD-GSO data, covering 50-120 keV, are featureless. However, in a differential spectrum between the pulse-decay phase and off-pulse phase, the second harmonic cyclotron resonance was detected in the GSO data at ~73 keV, with a depth of 1.6_-0.7^+0.9. This makes Her X-1 a 6th pulsar with established second harmonic resonance. Implications of these results are briefly discussed.