The typical accreting neutron star, Aquila X-1, was observed with Suzaku seven times in the decay phase of an outburst in 2007 September-October. Among them, the second to the fourth observations were performed 10 to 22 days after the outburst peak, when the source was in the hard state with a luminosity of 2x10^36 erg/sec. A unified spectral model for this type of objects approximately reproduced the 0.8--100~keV spectra obtained in these 3 observations. However, the spectra all exhibited an enigmatic hump-like excess around 30 keV, above the hard X-ray continuum which is interpreted as arising via Comptonization. The excess feature was confirmed to be significant against statistical and systematic uncertainties. It was successfully represented by a Gaussian centered at ~32 keV, with a width (sigma) of ~6 keV and an equivalent width of ~8.6 keV. Alternatively, the feature can also be explained by a recombination edge model, that produces a quasi-continuum above an edge energy of ~27 keV with an electron temperature of ~11 keV and an equivalent width of ~6.3 keV. These results are discussed in the context of atomic features of heavy elements synthesized via rapid-proton capture process during thermonuclear flashes.

We have analyzed X-ray data of SMC X-1 obtained with \textit{Suzaku}, and discovered a combination of an absorption line and an underlying broadened emission line centered at 6.4 keV in an observation performed on May 19, 2012. This absorption line is centered at $6.61_{-0.03}^{+0.02}$~keV with an absorption strength of $5.9_{-1.4}^{+1.7}$, naturally interpreted as an He$_{\alpha}$ resonance line of Fe at 6.7~keV that has a redshift of $4000_{-1300}^{+1400}\;\mathrm{km\;s^{-1}}$. Although \textit{Suzaku} observed this system ten times during 11 months in 2011--2012, the absorption feature has been seen only in a single observation when the NS was in a rising phase of the super-orbital modulation, which can be regarded as an egress from occultation by an extended accretion disk. We therefore attribute the line to a low density, highly ionized absorber in an accretion disk corona arising from the disk illuminated by the NS's intense X-rays. This interpretation also agrees with a discussion on the photoionization degree and the line depth.

The Hard X-ray Imager (HXI) onboard Hitomi (ASTRO-H) is an imaging spectrometer covering hard X-ray energies of 5-80 keV. Combined with the hard X-ray telescope, it enables imaging spectroscopy with an angular resolution of $1^\prime.7$ half-power diameter, in a field of view of $9^\prime\times9^\prime$. The main imager is composed of 4 layers of Si detectors and 1 layer of CdTe detector, stacked to cover wide energy band up to 80 keV, surrounded by an active shield made of BGO scintillator to reduce the background. The HXI started observations 12 days before the Hitomi loss, and successfully obtained data from G21.5$-$0.9, Crab and blank sky. Utilizing these data, we calibrate the detector response and study properties of in-orbit background. The observed Crab spectra agree well with a powerlaw model convolved with the detector response, within 5% accuracy. We find that albedo electrons in specified orbit strongly affect the background of Si top layer, and establish a screening method to reduce it. The background level over the full field of view after all the processing and screening is as low as the pre-flight requirement of $1$-$3\times10^{-4}$ counts s$^{-1}$ cm$^{-2}$ keV$^{-1}$.

The imaging and spectral performance of CdTe double-sided strip detectors (CdTe-DSDs) was evaluated for the ASTRO-H mission. The charcterized CdTe-DSDs have a strip pitch of 0.25 mm, an imaging area of 3.2 cm$\times$3.2 cm and a thickness of 0.75 mm. The detector was successfully operated at a temperature of $-20^\circ$C and with an applied bias voltage of 250 V. By using two-strip events as well as one-strip events for the event reconstruction, a good energy resolution of 2.0 keV at 59.5 keV and a sub-strip spatial resolution was achieved. The hard X-ray and gamma-ray response of CdTe-DSDs is complex due to the properties of CdTe and the small pixel effect. Therefore, one of the issues to investigate is the response of the CdTe-DSD. In order to investigate the spatial dependence of the detector response, we performed fine beam scan experiments at SPring-8, a synchrotron radiation facility. From these experiments, the depth structure of the electric field was determined as well as properties of carriers in the detector and successfully reproduced the experimental data with simulated spectra.

We promote the precise gamma-ray observation project Gamma-Ray Astro-Imager with Nuclear Emulsion (GRAINE), which uses balloon-borne emulsion gamma-ray telescopes. The emulsion telescope realizes observations with high angular resolution, polarization sensitivity, and large aperture area in the 0.01--100 GeV energy region. Herein, we report the data analysis of emulsion tracks and the first demonstration of gamma-ray imaging via an emulsion telescope by using the flight data from the balloon experiment performed in 2015 (GRAINE 2015). The emulsion films were scanned by the latest read-out system for a total area of 41 m$^2$ in three months, and then the gamma-ray event selection was automatically processed. Millions of electron-pair events are accumulated in the balloon-borne emulsion telescope. The emulsion telescope detected signals from a calibration source (gamma rays from the interaction of cosmic rays with an aluminum plate) with a high significance during the balloon observation and created a gamma-ray image consistent with the source size and the expected angular resolution in the energy range of 100--300 MeV. The flight performance obtained in the GRAINE 2015 experiment proves that balloon-borne emulsion telescope experiments with larger area are feasible while maintaining expected imaging performance.

Hard X-ray astronomical observatories in orbit suffer from a significant amount of background due to radioactivation induced by cosmic-ray protons and/or geomagnetically trapped protons. Within the framework of a full Monte Carlo simulation, we present modeling of in-orbit instrumental background dominated by radioactivation. To reduce the computation time required by straightforward simulations of delayed emissions from activated isotopes, we insert a semi-analytical calculation that converts production probabilities of radioactive isotopes by interaction of the primary protons into decay rates at measurement time of all secondary isotopes. Therefore, our simulation method is separated into three steps: (1) simulation of isotope production, (2) semi-analytical conversion to decay rates, and (3) simulation of decays of the isotopes at measurement time. This method is verified by a simple setup that has a CdTe semiconductor detector, and shows a 100-fold improvement in efficiency over the straightforward simulation. The simulation framework was tested against data measured with a CdTe sensor in the Hard X-ray Imager onboard the Hitomi X-ray Astronomy Satellite, which was put into a low Earth orbit with an altitude of 570 km and an inclination of 31 degrees, and thus experienced a large amount of irradiation from geomagnetically trapped protons during its passages through the South Atlantic Anomaly. The simulation is able to treat full histories of the proton irradiation and multiple measurement windows. The simulation results agree very well with the measured data, showing that the measured background is well described by the combination of proton-induced radioactivation of the CdTe detector itself and thick Bi4Ge3O12 scintillator shields, leakage of cosmic X-ray background and albedo gamma-ray radiation, and emissions from naturally contaminated isotopes in the detector system.

Lightning and thundercloud are the most dramatic natural particle accelerators on the Earth. Relativistic electrons accelerated by electric fields therein emit bremsstrahlung gamma rays, which have been detected at ground observations, by airborne detectors, and as terrestrial gamma-ray flashes (TGFs) from space. The energy of the gamma rays is sufficiently high to potentially invoke atmospheric photonuclear reactions 14N(gamma, n)13N, which would produce neutrons and eventually positrons via beta-plus decay of generated unstable radioactive isotopes, especially 13N. However, no clear observational evidence for the reaction has been reported to date. Here we report the first detection of neutron and positron signals from lightning with a ground observation. During a thunderstorm on 6 February 2017 in Japan, a TGF-like intense flash (within 1 ms) was detected at our monitoring sites 0.5-1.7 km away from the lightning. The subsequent initial burst quickly subsided with an exponential decay constant of 40-60 ms, followed by a prolonged line emission at about 0.511 megaelectronvolt (MeV), lasting for a minute. The observed decay timescale and spectral cutoff at about 10 MeV of the initial emission are well explained with de-excitation gamma rays from the nuclei excited by neutron capture. The centre energy of the prolonged line emission corresponds to the electron-positron annihilation, and hence is the conclusive indication of positrons produced after the lightning. Our detection of neutrons and positrons is unequivocal evidence that natural lightning triggers photonuclear reactions. No other natural event on the Earth is known to trigger photonuclear reactions. This discovery places lightning as only the second known natural channel on the Earth after the atmospheric cosmic-ray interaction, in which isotopes, such as 13C, 14C, and 15N, are produced.

Studies were made of the 1-70 keV persistent spectra of fifteen magnetars as a complete sample observed with Suzaku from 2006 to 2013. Combined with early NuSTAR observations of four hard X-ray emitters, nine objects showed a hard power-law emission dominating at $\gtrsim$10 keV with the 15--60 keV flux of $\sim$1-$11\times 10^{-11}$ ergs s$^{-1}$ cm$^{-2}$. The hard X-ray luminosity $L_{\rm h}$, relative to that of a soft-thermal surface radiation $L_{\rm s}$, tends to become higher toward younger and strongly magnetized objects. Updated from the previous study, their hardness ratio, defined as $\xi=L_{\rm h}/L_{\rm s}$, is correlated with the measured spin-down rate $\dot{P}$ as $\xi=0.62 \times (\dot{P}/10^{-11}\,{\rm s}\,{\rm s}^{-1})^{0.72}$, corresponding with positive and negative correlations of the dipole field strength $B_{\rm d}$ ($\xi \propto B_{\rm d}^{1.41}$) and the characteristic age $\tau_{\rm c}$ ($\xi \propto \tau_{\rm c}^{-0.68}$), respectively. Among our sample, five transients were observed during X-ray outbursts, and the results are compared with their long-term 1-10 keV flux decays monitored with Swift/XRT and RXTE/PCA. Fading curves of three bright outbursts are approximated by an empirical formula used in the seismology, showing a $\sim$10-40 d plateau phase. Transients show the maximum luminosities of $L_{\rm s}$$\sim$$10^{35}$ erg s$^{-1}$, which is comparable to those of the persistently bright ones, and fade back to $\lesssim$$10^{32}$ erg s$^{-1}$. Spectral properties are discussed in a framework of the magnetar hypothesis.

High-resolution X-ray spectroscopy with Hitomi was expected to resolve the origin of the faint unidentified E=3.5 keV emission line reported in several low-resolution studies of various massive systems, such as galaxies and clusters, including the Perseus cluster. We have analyzed the Hitomi first-light observation of the Perseus cluster. The emission line expected for Perseus based on the XMM-Newton signal from the large cluster sample under the dark matter decay scenario is too faint to be detectable in the Hitomi data. However, the previously reported 3.5 keV flux from Perseus was anomalously high compared to the sample-based prediction. We find no unidentified line at the reported high flux level. Taking into account the XMM measurement uncertainties for this region, the inconsistency with Hitomi is at a 99% significance for a broad dark-matter line and at 99.7% for a narrow line from the gas. We do not find anomalously high fluxes of the nearby faint K line or the Ar satellite line that were proposed as explanations for the earlier 3.5 keV detections. We do find a hint of a broad excess near the energies of high-n transitions of Sxvi (E=3.44 keV rest-frame) -- a possible signature of charge exchange in the molecular nebula and another proposed explanation for the unidentified line. While its energy is consistent with XMM pn detections, it is unlikely to explain the MOS signal. A confirmation of this interesting feature has to wait for a more sensitive observation with a future calorimeter experiment.

Clusters of galaxies are the most massive gravitationally-bound objects in the Universe and are still forming. They are thus important probes of cosmological parameters and a host of astrophysical processes. Knowledge of the dynamics of the pervasive hot gas, which dominates in mass over stars in a cluster, is a crucial missing ingredient. It can enable new insights into mechanical energy injection by the central supermassive black hole and the use of hydrostatic equilibrium for the determination of cluster masses. X-rays from the core of the Perseus cluster are emitted by the 50 million K diffuse hot plasma filling its gravitational potential well. The Active Galactic Nucleus of the central galaxy NGC1275 is pumping jetted energy into the surrounding intracluster medium, creating buoyant bubbles filled with relativistic plasma. These likely induce motions in the intracluster medium and heat the inner gas preventing runaway radiative cooling; a process known as Active Galactic Nucleus Feedback. Here we report on Hitomi X-ray observations of the Perseus cluster core, which reveal a remarkably quiescent atmosphere where the gas has a line-of-sight velocity dispersion of 164+/-10 km/s in a region 30-60 kpc from the central nucleus. A gradient in the line-of-sight velocity of 150+/-70 km/s is found across the 60 kpc image of the cluster core. Turbulent pressure support in the gas is 4% or less of the thermodynamic pressure, with large scale shear at most doubling that estimate. We infer that total cluster masses determined from hydrostatic equilibrium in the central regions need little correction for turbulent pressure.

Nowadays there is compelling evidence for the existence of dark matter in the Universe. A general consensus has been expressed on the need for a directional sensitive detector to confirm, with a complementary approach, the candidates found in conventional searches and to finally extend their sensitivity beyond the limit of neutrino-induced background. We propose here the use of a detector based on nuclear emulsions to measure the direction of WIMP-induced nuclear recoils. The production of nuclear emulsion films with nanometric grains is established. Several measurement campaigns have demonstrated the capability of detecting sub-micrometric tracks left by low energy ions in such emulsion films. Innovative analysis technologies with fully automated optical microscopes have made it possible to achieve the track reconstruction for path lengths down to one hundred nanometers and there are good prospects to further exceed this limit. The detector concept we propose foresees the use of a bulk of nuclear emulsion films surrounded by a shield from environmental radioactivity, to be placed on an equatorial telescope in order to cancel out the effect of the Earth rotation, thus keeping the detector at a fixed orientation toward the expected direction of galactic WIMPs. We report the schedule and cost estimate for a one-kilogram mass pilot experiment, aiming at delivering the first results on the time scale of six years.

The Suzaku data of the highly variable magnetar 1E 1547.0$-$5408, obtained during the 2009 January activity, were reanalyzed. The 2.07 s pulsation of the 15--40 keV emission detected with the HXD was found to be phase modulated, with a period of $36.0^{+4.5}_{-2.5}$ ks and an amplitude of $0.52 \pm 0.14$ s. The modulation waveform is suggested to be more square-wave like rather than sinusoidal. While the effect was confirmed with the 10--14 keV XIS data, the modulation amplitude decreased towards lower energies, becoming consistent with 0 below 4 keV. After the case of 4U 0142+61, this makes the 2nd example of this kind of behavior detected from magnetars. The effect can be interpreted as a manifestation of free precession of this magnetar, which is suggested to be oblately deformed under the presence of strong toroidal field of $\sim 10^{16}$ G.

Ages of the magnetar 1E 2259+586 and the associated supernova remnant CTB~109 were studied. Analyzing the Suzaku data of CTB~109, its age was estimated to be $\sim$14~kyr, which is much shorter than the measured characteristic age of 1E 2259+586, 230 kyr. This reconfirms the previously reported age discrepancy of this magnetar/remnant association, and suggests that the characteristic ages of magnetars are generally over-estimated as compared to their true ages. This discrepancy is thought to arise because the former are calculated without considering decay of the magnetic fields. This novel view is supported independently by much stronger Galactic-plane concentration of magnetars than other pulsars. The process of magnetic field decay in magnetars is mathematically modeled. It is implied that magnetars are much younger objects than previously considered, and can dominate new-born neutron stars.

As the endpoints of massive star evolution, neutron stars are enigmatic celestial objects characterized by extremely dense and exotic nuclear matter, magnetospheres with positrons (antimatter), rapid rotation and ultra-strong magnetic fields. Such an extreme environment has provided an accessible astrophysical laboratory to study physics under conditions unattainable on Earth and to tackle a range of fundamental questions related to: the aftermath of stellar evolution and the powerful explosions of massive stars, the equation of state and physics of some of the most exotic and magnetic stars in the Universe, the workings of the most powerful particle accelerators in our Galaxy and beyond, and the sources of gravitational waves that are yet to be detected. Recent observations revealed a great diversity of neutron stars, including ultra-strongly magnetized pulsars, referred to as "magnetars", and unusual types of accreting X-ray pulsars. In this white paper, we highlight the prospects of the upcoming X-ray mission, ASTRO-H, in studying these highly magnetized neutron stars.

Magnetars are a special type of neutron stars, considered to have extreme dipole magnetic fields reaching ~1e+11 T. The magnetar 4U 0142+61, one of prototypes of this class, was studied in broadband X-rays (0.5-70 keV) with the Suzaku observatory. In hard X-rays (15-40 keV), its 8.69 sec pulsations suffered slow phase modulations by +/-0.7 sec, with a period of ~15 hours. When this effect is interpreted as free precession of the neutron star, the object is inferred to deviate from spherical symmetry by ~1.6e-4 in its moments of inertia. This deformation, when ascribed to magnetic pressure, suggests a strong toroidal magnetic field, ~1e+12 T, residing inside the object. This provides one of the first observational approaches towards toroidal magnetic fields of magnetars.

Supernova remnants (SNRs) retain crucial information about both their parent explosion and circumstellar material left behind by their progenitor. However, the complexity of the interaction between supernova ejecta and ambient medium often blurs this information, and it is not uncommon for the basic progenitor type (Ia or core-collapse) of well-studied remnants to remain uncertain. Here we present a powerful new observational diagnostic to discriminate between progenitor types and constrain the ambient medium density of SNRs solely using Fe K-shell X-ray emission. We analyze all extant Suzaku observations of SNRs and detect Fe K alpha emission from 23 young or middle-aged remnants, including five first detections (IC 443, G292.0+1.8, G337.2-0.7, N49, and N63A). The Fe K alpha centroids clearly separate progenitor types, with the Fe-rich ejecta in Type Ia remnants being significantly less ionized than in core-collapse SNRs. Within each progenitor group, the Fe K alpha luminosity and centroid are well correlated, with more luminous objects having more highly ionized Fe. Our results indicate that there is a strong connection between explosion type and ambient medium density, and suggest that Type Ia supernova progenitors do not substantially modify their surroundings at radii of up to several parsecs. We also detect a K-shell radiative recombination continuum of Fe in W49B and IC 443, implying a strong circumstellar interaction in the early evolutionary phases of these core-collapse remnants.

Detecting the first electron pairs with nuclear emulsion allows a precise measurement of the direction of incident gamma-rays as well as their polarization. With recent innovations in emulsion scanning, emulsion analyzing capability is becoming increasingly powerful. Presently, we are developing a balloon-borne gamma-ray telescope using nuclear emulsion. An overview and a status of our telescope is given.

In this study, we are doing R&D for directional dark matter search with nuclear emulsion. First of all, higher resolution nuclear emulsion with fine silver halide crystals was developed in the production facility of emulsion at Nagoya university, and we confirmed that it can detect the expected nuclear recoil tracks. The readout of submicron tracks was required the new technology. We developed the expansion technique, and could readout the signal by shape analysis with optical microscopy. The two dimensional angular resolution is 36 degrees at the original track length of range from 150nm to 200nm with optical microscopy. Finally we demonstrated by using recoiled nuclei induced by 14.8MeV neutron, and confirmed the technique.Moreover, we developed the X-ray microscope system with SPring-8 as final check with higher resolution of selected candidate tracks with optical microscopy. The angular resolution was improved from 31 degrees with optical microscopy to 17degrees with X-ray microscopy at the track length of range from 150nm to 250nm. We are developing the practical system and planning for start of the test running with prototype detector.

We discuss the possibility to probe the QCD critical point during the dynamical black hole formation from a gravitational collapse of a massive star, where the temperature and the baryon chemical potential become as high as T ~ 90 MeV and $\mu_B$ ~ 1300 MeV. Comparison with the phase boundary in chiral effective models suggests that quark matter is likely to be formed before the horizon is formed. Furthermore, the QCD critical point may be probed during the black hole formation. The critical point is found to move in the lower temperature direction in asymmetric nuclear matter, and in some of the chiral models it is found to be in the reachable region during the black hole formation processes.

Improvements of in-orbit calibration of GSO scintillators in the Hard X-ray Detector on board Suzaku are reported. To resolve an apparent change of the energy scale of GSO which appeared across the launch for unknown reasons, consistent and thorough re-analyses of both pre-launch and in-orbit data have been performed. With laboratory experiments using spare hardware, the pulse height offset, corresponding to zero energy input, was found to change by ~0.5 of the full analog voltage scale, depending on the power supply. Furthermore, by carefully calculating all the light outputs of secondaries from activation lines used in the in-orbit gain determination, their energy deposits in GSO were found to be effectively lower, by several percent, than their nominal energies. Taking both these effects into account, the in-orbit data agrees with the on-ground measurements within ~5%, without employing the artificial correction introduced in the previous work (Kokubun et al. 2007). With this knowledge, we updated the data processing, the response, and the auxiliary files of GSO, and reproduced the HXD-PIN and HXD-GSO spectra of the Crab Nebula over 12-300 keV by a broken powerlaw with a break energy of ~110 keV.

We present the case for a dark matter detector with directional sensitivity. This document was developed at the 2009 CYGNUS workshop on directional dark matter detection, and contains contributions from theorists and experimental groups in the field. We describe the need for a dark matter detector with directional sensitivity; each directional dark matter experiment presents their project's status; and we close with a feasibility study for scaling up to a one ton directional detector, which would cost around $150M.

We investigate the detailed processes working in the drift of magnetic fields in molecular clouds. To the frictional force, whereby the magnetic force is transmitted to neutral molecules, ions contribute more than half only at cloud densities $n_{\rm H} < 10^4 {\rm cm}^{-3}$, and charged grains contribute more than 90% at $n_{\rm H} > 10^6 {\rm cm}^{-3}$. Thus grains play a decisive role in the process of magnetic flux loss. Approximating the flux loss time $t_B$ by a power law $t_B \propto B^{-\gamma}$, where $B$ is the mean field strength in the cloud, we find $\gamma \approx 2$, characteristic to ambipolar diffusion, only at $n_{\rm H} < 10^7 {\rm cm}^{-3}$. At higher densities, $\gamma$ decreases steeply with $n_{\rm H}$, and finally at $n_{\rm H} \approx n_{\rm dec} \approx {\rm a few} \times 10^{11} {\rm cm}^{-3}$, where magnetic fields effectively decouple from the gas, $\gamma << 1$ is attained, reminiscent of Ohmic dissipation, though flux loss occurs about 10 times faster than by Ohmic dissipation. Ohmic dissipation is dominant only at $n_{\rm H} > 1 \times 10^{12} {\rm cm}^{-3}$. While ions and electrons drift in the direction of magnetic force at all densities, grains of opposite charges drift in opposite directions at high densities, where grains are major contributors to the frictional force. Although magnetic flux loss occurs significantly faster than by Ohmic dissipation even at very high densities as $n_{\rm H} \approx n_{\rm dec}$, the process going on at high densities is quite different from ambipolar diffusion in which particles of opposite charges are supposed to drift as one unit.

We investigate where in protoplanetary disks magnetorotational instability operates, which can cause angular momentum transport in the disks. We investigate the spatial distribution of various charged particles and the unstable regions for a variety of models for protoplanetary disks taking into account the recombination of ions and electrons at grain surfaces, which is an important process in most parts of the disks. We find that for all the models there is an inner region which is magnetorotationally stable due to ohmic dissipation. This must make the accretion onto the central star non-steady. For the model of the minimum-mass solar nebula, the critical radius, inside of which the disk is stable, is about 20 AU, and the mass accretion rate just outside the critical radius is 10^{-7} - 10^{-6} M_{\odot} yr^{-1}. The stable region is smaller in a disk of lower column density. Dust grains in protoplanetary disks may grow by mutual sticking and may sediment toward the midplane of the disks. We find that the stable region shrinks as the grain size increases or the sedimentation proceeds. Therefore in the late evolutionary stages, protoplanetary disks can be magnetorotationally unstable even in the inner regions.

In this letter, we show that a massive black hole (MBH) which falls into the center of a galaxy in dynamical timescale leaves a weak cusp ($\rho \propto r^{-1/2}$) around it, which is in good agreement with the recent observations of luminous ellipticals by Hubble Space Telescope. Such event is a natural outcome of merging of two galaxies which have central MBHs. This is the only known mechanism to form weak cusps in luminous ellipticals. Therefore, the existence of the weak cusps indicates the central BHs of luminous ellipticals have fallen to the center from outside, most likely during a major merger event.

We investigated the dynamical reaction of the central region of galaxies to a falling massive black hole by N-body simulations. As the initial galaxy model, we used an isothermal King model and placed a massive black hole at around the half-mass radius of the galaxy. We found that the central core of the galaxy is destroyed by the heating due to the black hole and that a very weak density cusp ($\rho \propto r^{-\alpha}$, with $\alpha \sim 0.5$) is formed around the black hole. This result is consistent with recent observations of large elliptical galaxies with Hubble Space Telescope. The velocity of the stars becomes tangentially anisotropic in the inner region, while in the outer region the stars have radially anisotropic velocity dispersion. The radius of the weak cusp region is larger for larger black hole mass. Our result naturally explains the formation of the weak cusp found in the previous simulations of galaxy merging, and implies that the weak cusp observed in large elliptical galaxies may be formed by the heating process by sinking black holes during merging events.

We investigate the evolution of molecular abundance in quiescent protoplanetary disks which are presumed to be around weak-line T Tauri stars. In the region of surface density less than $10^2$ g cm$^{-2}$ (distance from the star $\gtrsim 10$ AU in the minimum- mass solar nebula), cosmic rays are barely attenuated even in the midplane of the disk and produce chemically active ions such as He$^+$ and H$_{3}^+$. Through reactions with these ions CO and N$_2$ are finally transformed into CO$_2$, NH$_3$, and HCN. In the region where the temperature is low enough for these products to freeze onto grains, considerable amount of carbon and nitrogen is locked up in the ice mantle and is depleted from the gas phase in a time scale $\lesssim 3\times 10^6$ yr. Oxidized (CO$_2$) ice and reduced (NH$_3$ and hydrocarbon) ice naturally coexist in this part of the disk. The molecular abundance both in the gas phase and in ice mantle varies significantly with the distance from the central star.