We have developed a method to make a spectral-line-based survey of hot cores, which represent an important stage of high-mass star formation, and applied the method to the data of the FUGIN (FOREST Unbiased Galactic plane Imaging survey with the Nobeyama 45-m telescope) survey. First, we select hot core candidates by searching the FUGIN data for the weak hot core tracer lines (HNCO and CH$_3$CN) by stacking, and then we conduct follow-up pointed observations on these candidates in C$^{34}$S, SO, OCS, HC$_3$N, HNCO, CH$_3$CN, and CH$_3$OH $J=2-1$ and $J=8-7$ lines to confirm and characterize them. We applied this method to the $l = 10^\circ-20^\circ$ portion of the FUGIN data and identified 22 "Hot Cores" (compact sources with more than two significant detection of the hot core tracer lines, i.e., SO, OCS, HC$_3$N, HNCO, CH$_3$CN, or CH$_3$OH $J=8-7$ lines) and 14 "Dense Clumps" (sources with more than two significant detection of C$^{34}$S, CH$_3$OH $J=2-1$, or the hot core tracer lines). The identified Hot Cores are found associated with signposts of high-mass star formation such as ATLASGAL clumps, WISE HII regions, and Class II methanol masers. For those associated with ATLASGAL clumps, their bolometric luminosity to clump mass ratios are consistent with the star formation stages centered at the hot core phase. The catalog of FUGIN Hot Cores provides a useful starting point for further statistical studies and detailed observations of high-mass star forming regions.

We have searched for exotic neutrino-electron interactions that could be produced by a neutrino millicharge, by a neutrino magnetic moment, or by dark photons using solar neutrinos in the XMASS-I liquid xenon detector. We observed no significant signals in 711 days of data. We obtain an upper limit for neutrino millicharge of 5.4$\times$10$^{-12} e$ at 90\% confidence level assuming all three species of neutrino have common millicharge. We also set flavor dependent limits assuming the respective neutrino flavor is the only one carrying a millicharge, $7.3 \times 10^{-12} e$ for $\nu_e$, $1.1 \times 10^{-11} e$ for $\nu_{\mu}$, and $1.1 \times 10^{-11} e$ for $\nu_{\tau}$. These limits are the most stringent yet obtained from direct measurements. We also obtain an upper limit for the neutrino magnetic moment of 1.8$\times$10$^{-10}$ Bohr magnetons. In addition, we obtain upper limits for the coupling constant of dark photons in the $U(1)_{B-L}$ model of 1.3$\times$10$^{-6}$ if the dark photon mass is 1$\times 10^{-3}$ MeV$/c^{2}$, and 8.8$\times$10$^{-5}$ if it is 10 MeV$/c^{2}$.

We report X-ray spectroscopic results for four giant solar flares occurred on 2005 September 7 (X17.0), 2005 September 8 (X5.4), 2005 September 9 (X6.2), and 2006 December 5 (X9.0), obtained from Earth albedo data with the X-ray imaging spectrometer (XIS) onboard Suzaku. The good energy resolution of the XIS (FWHM$\sim$100 eV) enables us to separate a number of line-like features and detect the underlying continuum emission. These features include Si He$\alpha$, Si Ly$\alpha$, S He$\alpha$, S Ly$\alpha$, Ar He$\alpha$, and Ca He$\alpha$ originating from solar flares as well as fluorescent Ar K$\alpha$ and Ar K$\beta$ from the Earth atmosphere. Absolute elemental abundances (X/H) averaged over the four flares are obtained to be $\sim$2.0 (Ca), $\sim$0.7 (Si), $\sim$0.3 (S), and $\sim$0.9 (Ar) at around flare peaks. This abundance pattern is similar to those of active stars' coronae showing inverse first ionization potential (i-FIP) effects, i.e., elemental abundances decrease with decreasing FIP with a turnover at the low end of FIP. The abundances are almost constant during the flares, with an exception of Si which increases by a factor of $\sim$2 in the decay phase. The evolution of the Si abundance is consistent with a picture that the i-FIP plasma originates from the chromosphere evaporation and then mixes with the surrounding low-FIP biased materials. Flare-to-flare abundance varied by a factor of 2, agreeing with past observations of solar flares. Finally, we emphasize that Earth albedo data acquired by X-ray astronomy satellites like Suzaku and XRISM can significantly contribute to studies of solar physics.

A search for Weakly Interacting Massive Particles (WIMPs) was conducted with the single-phase liquid-xenon detector XMASS through inelastic scattering in which $^{129}$Xe nuclei were excited, using an exposure ($\rm 327\; kg \times 800.0 \; days$) 48 times larger than that of our previous study. The inelastic excitation sensitivity was improved by detailed evaluation of background, event classification based on scintillation timing that distinguished $\gamma$-rays and $\beta$-rays, and simultaneous fitting of the energy spectra of $\gamma$-like and $\beta$-like samples. No evidence of a WIMP signal was found. Thus, we set the upper limits of the inelastic channel cross section at 90\% confidence level, for example, $4.1\times 10^{-39} \;{\rm cm^2}$ for a $200\; {\rm GeV}/c^2$ WIMP. This result provides the most stringent limits on the SD WIMP-neutron interaction and is better by a factor of 7.7 at $200\;{\rm GeV}/c^2$ than the existing experimental limit.

In recent years, the outskirts of galaxy clusters have emerged as one of the new frontiers and unique laboratories for studying the growth of large scale structure in the universe. Modern cosmological hydrodynamical simulations make firm and testable predictions of the thermodynamic and chemical evolution of the X-ray emitting intracluster medium. However, recent X-ray and Sunyaev-Zeldovich effect observations have revealed enigmatic disagreements with theoretical predictions, which have motivated deeper investigations of a plethora of astrophysical processes operating in the virialization region in the cluster outskirts. Much of the physics of cluster outskirts is fundamentally different from that of cluster cores, which has been the main focus of X-ray cluster science over the past several decades. A next-generation X-ray telescope, equipped with sub-arcsecond spatial resolution over a large field of view along with a low and stable instrumental background, is required in order to reveal the full story of the growth of galaxy clusters and the cosmic web and their applications for cosmology.

XMASS-I is a single-phase liquid xenon detector whose purpose is direct detection of dark matter. To achieve the low background requirements necessary in the detector, a new model of photomultiplier tubes (PMTs), R10789, with a hexagonal window was developed based on the R8778 PMT used in the XMASS prototype detector. We screened the numerous component materials for their radioactivity. During development, the largest contributions to the reduction of radioactivity came from the stem and the dynode support. The glass stem was exchanged to the Kovar alloy one and the ceramic support were changed to the quartz one. R10789 is the first model of Hamamatsu Photonics K. K. that adopted these materials for low background purposes and provided a groundbreaking step for further reductions of radioactivity in PMTs. Measurements with germanium detectors showed 1.2$\pm$0.3 mBq/PMT of $^{226}$Ra, less than 0.78 mBq/PMT of $^{228}$Ra, 9.1$\pm$2.2 mBq/PMT of $^{40}$K, and 2.8$\pm$0.2 mBq/PMT of $^{60}$Co. In this paper, the radioactive details of the developed R10789 are described together with our screening methods and the components of the PMT.

A search for dark matter using an underground single-phase liquid xenon detector was conducted at the Kamioka Observatory in Japan, particularly for Weakly Interacting Massive Particles (WIMPs). We have used 705.9 live days of data in a fiducial volume containing 97 kg of liquid xenon at the center of the detector. The event rate in the fiducial volume after the data reduction was ${\rm (4.2 \pm 0.2) \times 10^{-3} \, day^{-1}kg^{-1} keV_{ee}^{-1}}$ at ${\rm 5 \, keV_{ee}}$, with a signal efficiency of ${\rm 20\%}$. All the remaining events are consistent with our background evaluation, mostly of the "mis-reconstructed events" originated from $^{210}$Pb in the copper plates lining the detector's inner surface. The obtained upper limit on a spin-independent WIMP-nucleon cross section was ${\rm 2.2 \times 10^{-44} \, cm^{2}}$ for a WIMP mass of ${\rm 60 \, GeV/c^{2}}$ at the $90\%$ confidence level, which was the most stringent limit among results from single-phase liquid xenon detectors.

A search for dark matter (DM) with mass in the sub-GeV region (0.32-1 GeV) was conducted by looking for an annual modulation signal in XMASS, a single-phase liquid xenon detector. Inelastic nuclear scattering accompanied by bremsstrahlung emission was used to search down to an electron equivalent energy of 1 keV. The data used had a live time of 2.8 years (3.5 years in calendar time), resulting in a total exposure of 2.38 ton-years. No significant modulation signal was observed and 90% confidence level upper limits of $1.6 \times 10^{-33}$ cm$^2$ at 0.5 GeV was set for the DM-nucleon cross section. This is the first experimental result of a search for DM mediated by the bremsstrahlung effect. In addition, a search for DM with mass in the multi-GeV region (4-20 GeV) was conducted with a lower energy threshold than previous analysis of XMASS. Elastic nuclear scattering was used to search down to a nuclear recoil equivalent energy of 2.3 keV, and upper limits of 2.9 $\times$10$^{-42}$ cm$^2$ at 8 GeV was obtained.

High-resolution spectroscopy of the core of the Perseus Cluster of galaxies, using the $Hitomi$ satellite above 2 keV and the $XMM$-$Newton$ Reflection Grating Spectrometer at lower energies, provides reliable constraints on the abundances of O, Ne, Mg, Si, S, Ar, Ca, Cr, Mn, Fe, and Ni. Accounting for all known systematic uncertainties, the Ar/Fe, Ca/Fe, and Ni/Fe ratios are determined with a remarkable precision of less than 10%, while the constraints on Si/Fe, S/Fe, and Cr/Fe are at the 15% level, and Mn/Fe is measured with a 20% uncertainty. The average biases in determining the chemical composition using archival CCD spectra from $XMM$-$Newton$ and $Suzaku$ range typically from 15-40%. A simple model in which the enrichment pattern in the Perseus Cluster core and the proto-solar nebula are identical gives a surprisingly good description of the high-resolution X-ray spectroscopy results, with $\chi^2=10.7$ for 10 d.o.f. However, this pattern is challenging to reproduce with linear combinations of existing supernova nucleosynthesis calculations, particularly given the precise measurements of intermediate $\alpha$-elements enabled by $Hitomi$. We discuss in detail the degeneracies between various supernova progenitor models and explosion mechanisms, and the remaining uncertainties in these theoretical models. We suggest that including neutrino physics in the core-collapse supernova yield calculations may improve the agreement with the observed pattern of $\alpha$-elements in the Perseus Cluster core. Our results provide a complementary benchmark for testing future nucleosynthesis calculations required to understand the origin of chemical elements.

Hidden photons and axion-like particles are candidates for cold dark matter if they were produced non-thermally in the early universe. We conducted a search for both of these bosons using 800 live-days of data from the XMASS detector with 327 kg of liquid xenon in the fiducial volume. No significant signal was observed, and thus we set constraints on the $\alpha' / \alpha$ parameter related to kinetic mixing of hidden photons and the coupling constant $g_{Ae}$ of axion-like particles in the mass range from 40 to 120 keV/$c^2$, resulting in $\alpha' / \alpha < 6 \times 10^{-26}$ and $g_{Ae} < 4 \times 10^{-13}$. These limits are the most stringent derived from both direct and indirect searches to date.

We present results of a search for unidentified line emission and absorption signals in the 2-12keV energy band of spectra extracted from Perseus Cluster core region observations obtained with the 5eV energy resolution Hitomi Soft X-ray Spectrometer. No significant unidentified line emission or absorption is found. Line flux upper limits (1 sigma per resolution element) vary with photon energy and assumed intrinsic width, decreasing from 100 photons cm$^{-2}$ s$^{-1}$ sr $^{-1}$ at 2keV to $<10$ cm$^{-2}$ s$^{-1}$ sr $^{-1}$ over most of the 5-10 keV energy range for a Gaussian line with Doppler broadening of 640 km/s. Limits for narrower and broader lines have a similar energy dependence and are systematically smaller and larger, respectively. These line flux limits are used to constrain the decay rate of hypothetical dark matter candidates. For the sterile neutrino decay rate, new constraints over the the mass range of 4-24 keV with mass resolution better than any previous X-ray analysis are obtained. Additionally, the accuracy of relevant thermal spectral models and atomic data are evaluated. The Perseus cluster spectra may be described by a composite of multi-temperature thermal and AGN power-law continua. Superposed on these, a few line emission signals possibly originating from unmodeled atomic processes (including Si XIV and Fe XXV) are marginally detected and tabulated. Comparisons with previous X-ray upper limits and future prospects for dark matter searches using high-energy resolution spectroscopy are discussed.

We present the results from the Hitomi Soft Gamma-ray Detector (SGD) observation of the Crab nebula. The main part of SGD is a Compton camera, which in addition to being a spectrometer, is capable of measuring polarization of gamma-ray photons. The Crab nebula is one of the brightest X-ray / gamma-ray sources on the sky, and, the only source from which polarized X-ray photons have been detected. SGD observed the Crab nebula during the initial test observation phase of Hitomi. We performed the data analysis of the SGD observation, the SGD background estimation and the SGD Monte Carlo simulations, and, successfully detected polarized gamma-ray emission from the Crab nebula with only about 5 ks exposure time. The obtained polarization fraction of the phase-integrated Crab emission (sum of pulsar and nebula emissions) is (22.1 $\pm$ 10.6)% and, the polarization angle is 110.7$^o$ + 13.2 / $-$13.0$^o$ in the energy range of 60--160 keV (The errors correspond to the 1 sigma deviation). The confidence level of the polarization detection was 99.3%. The polarization angle measured by SGD is about one sigma deviation with the projected spin axis of the pulsar, 124.0$^o$ $\pm$0.1$^o$.

The X-ray Integral Field Unit (X-IFU) is the high resolution X-ray spectrometer of the ESA Athena X-ray observatory. Over a field of view of 5' equivalent diameter, it will deliver X-ray spectra from 0.2 to 12 keV with a spectral resolution of 2.5 eV up to 7 keV on ~5 arcsecond pixels. The X-IFU is based on a large format array of super-conducting molybdenum-gold Transition Edge Sensors cooled at about 90 mK, each coupled with an absorber made of gold and bismuth with a pitch of 249 microns. A cryogenic anti-coincidence detector located underneath the prime TES array enables the non X-ray background to be reduced. A bath temperature of about 50 mK is obtained by a series of mechanical coolers combining 15K Pulse Tubes, 4K and 2K Joule-Thomson coolers which pre-cool a sub Kelvin cooler made of a 3He sorption cooler coupled with an Adiabatic Demagnetization Refrigerator. Frequency domain multiplexing enables to read out 40 pixels in one single channel. A photon interacting with an absorber leads to a current pulse, amplified by the readout electronics and whose shape is reconstructed on board to recover its energy with high accuracy. The defocusing capability offered by the Athena movable mirror assembly enables the X-IFU to observe the brightest X-ray sources of the sky (up to Crab-like intensities) by spreading the telescope point spread function over hundreds of pixels. Thus the X-IFU delivers low pile-up, high throughput (>50%), and typically 10 eV spectral resolution at 1 Crab intensities, i.e. a factor of 10 or more better than Silicon based X-ray detectors. In this paper, the current X-IFU baseline is presented, together with an assessment of its anticipated performance in terms of spectral resolution, background, and count rate capability. The X-IFU baseline configuration will be subject to a preliminary requirement review that is scheduled at the end of 2018.

An annual modulation signal due to the Earth orbiting around the Sun would be one of the strongest indications of the direct detection of dark matter. In 2016, we reported a search for dark matter by looking for this annual modulation with our single-phase liquid xenon XMASS-I detector. That analysis resulted in a slightly negative modulation amplitude at low energy. In this work, we included more than one year of additional data, which more than doubles the exposure to 800 live days with the same 832 kg target mass. When we assume weakly interacting massive particle (WIMP) dark matter elastically scattering on the xenon target, the exclusion upper limit for the WIMP-nucleon cross section was improved by a factor of 2 to 1.9$\times$10$^{-41}$cm$^2$ at 8 GeV/c$^2$ at 90\% confidence level with our newly implemented data selection through a likelihood method. For the model-independent case, without assuming any specific dark matter model, we obtained more consistency with the null hypothesis than before with a $p$-value of 0.11 in the 1$-$20 keV energy region. This search probed this region with an exposure that was larger than that of DAMA/LIBRA. We also did not find any significant amplitude in the data for periodicity with periods between 50 and 600 days in the energy region between 1 to 6 keV.

We present results from the Hitomi X-ray observation of a young composite-type supernova remnant (SNR) G21.5$-$0.9, whose emission is dominated by the pulsar wind nebula (PWN) contribution. The X-ray spectra in the 0.8-80 keV range obtained with the Soft X-ray Spectrometer (SXS), Soft X-ray Imager (SXI) and Hard X-ray Imager (HXI) show a significant break in the continuum as previously found with the NuSTAR observation. After taking into account all known emissions from the SNR other than the PWN itself, we find that the Hitomi spectra can be fitted with a broken power law with photon indices of $\Gamma_1=1.74\pm0.02$ and $\Gamma_2=2.14\pm0.01$ below and above the break at $7.1\pm0.3$ keV, which is significantly lower than the NuSTAR result ($\sim9.0$ keV). The spectral break cannot be reproduced by time-dependent particle injection one-zone spectral energy distribution models, which strongly indicates that a more complex emission model is needed, as suggested by recent theoretical models. We also search for narrow emission or absorption lines with the SXS, and perform a timing analysis of PSR J1833$-$1034 with the HXI and SGD. No significant pulsation is found from the pulsar. However, unexpectedly, narrow absorption line features are detected in the SXS data at 4.2345 keV and 9.296 keV with a significance of 3.65 $\sigma$. While the origin of these features is not understood, their mere detection opens up a new field of research and was only possible with the high resolution, sensitivity and ability to measure extended sources provided by an X-ray microcalorimeter.

We present the result of the in-flight calibration of the effective area of the Soft X-ray Spectrometer (SXS) onboard the Hitomi X-ray satellite using an observation of the Crab nebula. We corrected for the artifacts when observing high count rate sources with the X-ray microcalorimeter. We then constructed a spectrum in the 0.5-20 keV band, which we modeled with a single power-law continuum attenuated by an interstellar extinction. We evaluated the systematic uncertainty upon the spectral parameters by various calibration items. In the 2-12 keV band, the SXS result is consistent with the literature values in flux (2.20 $\pm$ 0.08) $\times$10$^{-8}$ erg s$^{-1}$ cm$^{-2}$ with a 1$\sigma$ statistical uncertainty) but is softer in the power-law index (2.19 $\pm$ 0.11). The discrepancy is attributable to the systematic uncertainty of about $+$6/$-$7% and $+$2/$-$5% respectively for the flux and the power-law index. The softer spectrum is affected primarily by the systematic uncertainty of the Dewar gate valve transmission and the event screening.

The present paper investigates the temperature structure of the X-ray emitting plasma in the core of the Perseus cluster using the 1.8--20.0 keV data obtained with the Soft X-ray Spectrometer (SXS) onboard the Hitomi Observatory. A series of four observations were carried out, with a total effective exposure time of 338 ks and covering a central region $\sim7'$ in diameter. The SXS was operated with an energy resolution of $\sim$5 eV (full width at half maximum) at 5.9 keV. Not only fine structures of K-shell lines in He-like ions but also transitions from higher principal quantum numbers are clearly resolved from Si through Fe. This enables us to perform temperature diagnostics using the line ratios of Si, S, Ar, Ca, and Fe, and to provide the first direct measurement of the excitation temperature and ionization temperature in the Perseus cluster. The observed spectrum is roughly reproduced by a single temperature thermal plasma model in collisional ionization equilibrium, but detailed line ratio diagnostics reveal slight deviations from this approximation. In particular, the data exhibit an apparent trend of increasing ionization temperature with increasing atomic mass, as well as small differences between the ionization and excitation temperatures for Fe, the only element for which both temperatures can be measured. The best-fit two-temperature models suggest a combination of 3 and 5 keV gas, which is consistent with the idea that the observed small deviations from a single temperature approximation are due to the effects of projection of the known radial temperature gradient in the cluster core along the line of sight. Comparison with the Chandra/ACIS and the XMM-Newton/RGS results on the other hand suggests that additional lower-temperature components are present in the ICM but not detectable by Hitomi SXS given its 1.8--20 keV energy band.

The Hitomi SXS spectrum of the Perseus cluster, with $\sim$5 eV resolution in the 2-9 keV band, offers an unprecedented benchmark of the atomic modeling and database for hot collisional plasmas. It reveals both successes and challenges of the current atomic codes. The latest versions of AtomDB/APEC (3.0.8), SPEX (3.03.00), and CHIANTI (8.0) all provide reasonable fits to the broad-band spectrum, and are in close agreement on best-fit temperature, emission measure, and abundances of a few elements such as Ni. For the Fe abundance, the APEC and SPEX measurements differ by 16%, which is 17 times higher than the statistical uncertainty. This is mostly attributed to the differences in adopted collisional excitation and dielectronic recombination rates of the strongest emission lines. We further investigate and compare the sensitivity of the derived physical parameters to the astrophysical source modeling and instrumental effects. The Hitomi results show that an accurate atomic code is as important as the astrophysical modeling and instrumental calibration aspects. Substantial updates of atomic databases and targeted laboratory measurements are needed to get the current codes ready for the data from the next Hitomi-level mission.

We present Hitomi observations of N132D, a young, X-ray bright, O-rich core-collapse supernova remnant in the Large Magellanic Cloud (LMC). Despite a very short observation of only 3.7 ks, the Soft X-ray Spectrometer (SXS) easily detects the line complexes of highly ionized S K and Fe K with 16-17 counts in each. The Fe feature is measured for the first time at high spectral resolution. Based on the plausible assumption that the Fe K emission is dominated by He-like ions, we find that the material responsible for this Fe emission is highly redshifted at ~800 km/s compared to the local LMC interstellar medium (ISM), with a 90% credible interval of 50-1500 km/s if a weakly informative prior is placed on possible line broadening. This indicates (1) that the Fe emission arises from the supernova ejecta, and (2) that these ejecta are highly asymmetric, since no blue-shifted component is found. The S K velocity is consistent with the local LMC ISM, and is likely from swept-up ISM material. These results are consistent with spatial mapping that shows the He-like Fe concentrated in the interior of the remnant and the S tracing the outer shell. The results also show that even with a very small number of counts, direct velocity measurements from Doppler-shifted lines detected in extended objects like supernova remnants are now possible. Thanks to the very low SXS background of ~1 event per spectral resolution element per 100 ks, such results are obtainable during short pointed or slew observations with similar instruments. This highlights the power of high-spectral-resolution imaging observations, and demonstrates the new window that has been opened with Hitomi and will be greatly widened with future missions such as the X-ray Astronomy Recovery Mission (XARM) and Athena.

We report a Hitomi observation of IGR J16318-4848, a high-mass X-ray binary system with an extremely strong absorption of N_H~10^{24} cm^{-2}. Previous X-ray studies revealed that its spectrum is dominated by strong fluorescence lines of Fe as well as continuum emission. For physical and geometrical insight into the nature of the reprocessing material, we utilize the high spectroscopic resolving power of the X-ray microcalorimeter (the soft X-ray spectrometer; SXS) and the wide-band sensitivity by the soft and hard X-ray imager (SXI and HXI) aboard Hitomi. Even though photon counts are limited due to unintended off-axis pointing, the SXS spectrum resolves Fe K{\alpha_1} and K{\alpha_2} lines and puts strong constraints on the line centroid and width. The line width corresponds to the velocity of 160^{+300}_{-70} km s^{-1}. This represents the most accurate, and smallest, width measurement of this line made so far from any X-ray binary, much less than the Doppler broadening and shift expected from speeds which are characteristic of similar systems. Combined with the K-shell edge energy measured by the SXI and HXI spectra, the ionization state of Fe is estimated to be in the range of Fe I--IV. Considering the estimated ionization parameter and the distance between the X-ray source and the absorber, the density and thickness of the materials are estimated. The extraordinarily strong absorption and the absence of a Compton shoulder component is confirmed. These characteristics suggest reprocessing materials which are distributed in a narrow solid angle or scattering primarily with warm free electrons or neutral hydrogen.

The origin of the narrow Fe-K{\alpha} fluorescence line at 6.4 keV from active galactic nuclei has long been under debate; some of the possible sites are the outer accretion disk, the broad line region, a molecular torus, or interstellar/intracluster media. In February-March 2016, we performed the first X-ray microcalorimeter spectroscopy with the Soft X-ray Spectrometer (SXS) onboard the Hitomi satellite of the Fanaroff-Riley type I radio galaxy NGC 1275 at the center of the Perseus cluster of galaxies. With the high energy resolution of ~5 eV at 6 keV achieved by Hitomi/SXS, we detected the Fe-K{\alpha} line with ~5.4 {\sigma} significance. The velocity width is constrained to be 500-1600 km s$^{-1}$ (FWHM for Gaussian models) at 90% confidence. The SXS also constrains the continuum level from the NGC 1275 nucleus up to ~20 keV, giving an equivalent width ~20 eV of the 6.4 keV line. Because the velocity width is narrower than that of broad H{\alpha} line of ~2750 km s$^{-1}$, we can exclude a large contribution to the line flux from the accretion disk and the broad line region. Furthermore, we performed pixel map analyses on the Hitomi/SXS data and image analyses on the Chandra archival data, and revealed that the Fe-K{\alpha} line comes from a region within ~1.6 kpc from the NGC 1275 core, where an active galactic nucleus emission dominates, rather than that from intracluster media. Therefore, we suggest that the source of the Fe-K{\alpha} line from NGC 1275 is likely a low-covering fraction molecular torus or a rotating molecular disk which probably extends from a pc to hundreds pc scale in the active galactic nucleus system.

In theories with the large extra dimensions beyond the standard 4-dimensional spacetime, axions could propagate in such extra dimensions, and acquire Kaluza-Klein (KK) excitations. These KK axions are produced in the Sun and could solve unexplained heating of the solar corona. While most of the solar KK axions escape from the solar system, a small fraction is gravitationally trapped in orbits around the Sun. They would decay into two photons inside a terrestrial detector. The event rate is expected to modulate annually depending on the distance from the Sun. We have searched for the annual modulation signature using $832\times 359$ kg$\cdot$days of XMASS-I data. No significant event rate modulation is found, and hence we set the first experimental constraint on the KK axion-photon coupling of $4.8 \times 10^{-12}\, \mathrm{GeV}^{-1}$ at 90% confidence level for a KK axion number density of $\bar{n}_\mathrm{a} = 4.07 \times 10^{13}\, \mathrm{m}^{-3}$, the total number of extra dimensions $n = 2$, and the number of extra dimensions $\delta = 2$ that axions can propagate in.

Extending the earlier measurements reported in Hitomi collaboration (2016, Nature, 535, 117), we examine the atmospheric gas motions within the central 100~kpc of the Perseus cluster using observations obtained with the Hitomi satellite. After correcting for the point spread function of the telescope and using optically thin emission lines, we find that the line-of-sight velocity dispersion of the hot gas is remarkably low and mostly uniform. The velocity dispersion reaches maxima of approximately 200~km~s$^{-1}$ toward the central active galactic nucleus (AGN) and toward the AGN inflated north-western `ghost' bubble. Elsewhere within the observed region, the velocity dispersion appears constant around 100~km~s$^{-1}$. We also detect a velocity gradient with a 100~km~s$^{-1}$ amplitude across the cluster core, consistent with large-scale sloshing of the core gas. If the observed gas motions are isotropic, the kinetic pressure support is less than 10\% of the thermal pressure support in the cluster core. The well-resolved optically thin emission lines have Gaussian shapes, indicating that the turbulent driving scale is likely below 100~kpc, which is consistent with the size of the AGN jet inflated bubbles. We also report the first measurement of the ion temperature in the intracluster medium, which we find to be consistent with the electron temperature. In addition, we present a new measurement of the redshift to the brightest cluster galaxy NGC~1275.

We present a joint X-ray, optical and weak-lensing analysis for X-ray luminous galaxy clusters selected from the MCXC (Meta-Catalog of X-Ray Detected Clusters of Galaxies) cluster catalog in the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) survey field with S16A data, As a pilot study of our planned series papers, we measure hydrostatic equilibrium (H.E.) masses using XMM-Newton data for four clusters in the current coverage area out of a sample of 22 MCXC clusters. We additionally analyze a non-MCXC cluster associated with one MCXC cluster. We show that H.E. masses for the MCXC clusters are correlated with cluster richness from the CAMIRA catalog (Oguri et al. 2017), while that for the non-MCXC cluster deviates from the scaling relation. The mass normalization of the relationship between the cluster richness and H.E. mass is compatible with one inferred by matching CAMIRA cluster abundance with a theoretical halo mass function. The mean gas mass fraction based on H.E. masses for the MCXC clusters is $\langle f_{\rm gas} \rangle = 0.125\pm0.012$ at spherical overdensity $\Delta=500$, which is $\sim80-90$ percent of the cosmic mean baryon fraction, $\Omega_b/\Omega_m$, measured by cosmic microwave background experiments. We find that the mean baryon fraction estimated from X-ray and HSC-SSP optical data is comparable to $\Omega_b/\Omega_m$. A weak-lensing shear catalog of background galaxies, combined with photometric redshifts, is currently available only for three clusters in our sample. Hydrostatic equilibrium masses roughly agree with weak-lensing masses, albeit with large uncertainty. This study demonstrates that further multiwavelength study for a large sample of clusters using X-ray, HSC-SSP optical and weak lensing data will enable us to understand cluster physics and utilize cluster-based cosmology.

Thanks to its high spectral resolution (~5 eV at 6 keV), the Soft X-ray Spectrometer (SXS) on board Hitomi enables us to measure the detailed structure of spatially resolved emission lines from highly ionized ions in galaxy clusters for the first time. In this series of papers, using the SXS we have measured the velocities of gas motions, metallicities and the multi-temperature structure of the gas in the core of the Perseus cluster. Here, we show that when inferring physical properties from line emissivities in systems like Perseus, the resonant scattering (RS) effect should be taken into account. In the Hitomi waveband, RS mostly affects the FeXXV He$\alpha$ line ($w$) - the strongest line in the spectrum. The flux measured by Hitomi in this line is suppressed by a factor ~1.3 in the inner ~30 kpc, compared to predictions for an optically thin plasma; the suppression decreases with the distance from the center. The $w$ line also appears slightly broader than other lines from the same ion. The observed distortions of the $w$ line flux, shape and distance dependence are all consistent with the expected effect of the resonant scattering in the Perseus core. By measuring the ratio of fluxes in optically thick ($w$) and thin (FeXXV forbidden, He$\beta$, Ly$\alpha$) lines, and comparing these ratios with predictions from Monte Carlo radiative transfer simulations, the velocities of gas motions have been obtained. The results are consistent with the direct measurements of gas velocities from line broadening described elsewhere in this series, although the systematic and statistical uncertainties remain significant. Further improvements in the predictions of line emissivities in plasma models, and deeper observations with future X-ray missions will enable RS measurements to provide powerful constraints on the amplitude and anisotropy of clusters gas motions.

To search for giant X-ray pulses correlated with the giant radio pulses (GRPs) from the Crab pulsar, we performed a simultaneous observation of the Crab pulsar with the X-ray satellite Hitomi in the 2 -- 300 keV band and the Kashima NICT radio observatory in the 1.4 -- 1.7 GHz band with a net exposure of about 2 ks on 25 March 2016, just before the loss of the Hitomi mission.The timing performance of the Hitomi instruments was confirmed to meet the timing requirement and about 1,000 and 100 GRPs were simultaneously observed at the main and inter-pulse phases, respectively, and we found no apparent correlation between the giant radio pulses and the X-ray emission in either the main or inter-pulse phases.All variations are within the 2 sigma fluctuations of the X-ray fluxes at the pulse peaks, and the 3 sigma upper limits of variations of main- or inter- pulse GRPs are 22\% or 80\% of the peak flux in a 0.20 phase width, respectively, in the 2 -- 300 keV band.The values become 25\% or 110\% for main or inter-pulse GRPs, respectively, when the phase width is restricted into the 0.03 phase.Among the upper limits from the Hitomi satellite, those in the 4.5-10 keV and the 70-300 keV are obtained for the first time, and those in other bands are consistent with previous reports.Numerically, the upper limits of main- and inter-pulse GRPs in the 0.20 phase width are about (2.4 and 9.3) $\times 10^{-11}$ erg cm$^{-2}$, respectively. No significant variability in pulse profiles implies that the GRPs originated from a local place within the magnetosphere and the number of photon-emitting particles temporally increases.However, the results do not statistically rule out variations correlated with the GRPs, because the possible X-ray enhancement may appear due to a $>0.02$\% brightening of the pulse-peak flux under such conditions.

The Crab nebula originated from a core-collapse supernova (SN) explosion observed in 1054 A.D. When viewed as a supernova remnant (SNR), it has an anomalously low observed ejecta mass and kinetic energy for an Fe-core collapse SN. Intensive searches were made for a massive shell that solves this discrepancy, but none has been detected. An alternative idea is that the SN1054 is an electron-capture (EC) explosion with a lower explosion energy by an order of magnitude than Fe-core collapse SNe. In the X-rays, imaging searches were performed for the plasma emission from the shell in the Crab outskirts to set a stringent upper limit to the X-ray emitting mass. However, the extreme brightness of the source hampers access to its vicinity. We thus employed spectroscopic technique using the X-ray micro-calorimeter onboard the Hitomi satellite. By exploiting its superb energy resolution, we set an upper limit for emission or absorption features from yet undetected thermal plasma in the 2-12 keV range. We also re-evaluated the existing Chandra and XMM-Newton data. By assembling these results, a new upper limit was obtained for the X-ray plasma mass of <~ 1Mo for a wide range of assumed shell radius, size, and plasma temperature both in and out of the collisional equilibrium. To compare with the observation, we further performed hydrodynamic simulations of the Crab SNR for two SN models (Fe-core versus EC) under two SN environments (uniform ISM versus progenitor wind). We found that the observed mass limit can be compatible with both SN models if the SN environment has a low density of <~ 0.03 cm-3 (Fe core) or <~ 0.1 cm-3 (EC) for the uniform density, or a progenitor wind density somewhat less than that provided by a mass loss rate of 10-5 Mo yr-1 at 20 km s-1 for the wind environment.

We report on the results of a $Suzaku$ observation of the plasma in the filament located between the two massive clusters of galaxies Abell 399 and Abell 401. Abell 399 ($z$=0.0724) and Abell 401 ($z$=0.0737) are expected to be in the initial phase of a cluster merger. In the region between the two clusters, we find a clear enhancement in the temperature of the filament plasma from 4 keV (expected value from a typical cluster temperature profile) to $kT\sim$6.5 keV. Our analysis also shows that filament plasma is present out to a radial distance of 15' (1.3 Mpc) from a line connecting the two clusters. The temperature profile is characterized by an almost flat radial shape with $kT\sim$6-7 keV within 10' or $\sim$0.8 Mpc. Across $r$=8'~from the axis, the temperature of the filament plasma shows a drop from 6.3 keV to 5.1 keV, indicating the presence of a shock front. The Mach number based on the temperature drop is estimated to be ${\cal M}\sim$1.3. We also successfully determined the abundance profile up to 15' (1.3 Mpc), showing an almost constant value ($Z$=0.3 solar) at the cluster outskirt. We estimated the Compton $y$-parameter to be $\sim$14.5$\pm1.3\times10^{-6}$, which is in agreement with $Planck$'s results (14-17$\times10^{-6}$ on the filament). The line of sight depth of the filament is $l\sim$1.1 Mpc, indicating that the geometry of filament is likely a pancake shape rather than cylindrical. The total mass of the filamentary structure is $\sim$7.7$\times10^{13}~\rm M_{\odot}$. We discuss a possible interpretation of the drop of X-ray emission at the rim of the filament, which was pushed out by the merging activity and formed by the accretion flow induced by the gravitational force of the filament.

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.

We report $Suzaku$ X-ray observations of the dark subhalo associated with the merging group of NGC 4839 in the Coma cluster. The X-ray image exhibits an elongated tail toward the southwest. The X-ray peak shifts approximately $1'$ away from the weak-lensing mass center toward the opposite direction of the Coma cluster center. We investigated the temperature, normalization, pressure, and entropy distributions around the subhalo. Excluding the X-ray tail, the temperature beyond the truncation radius is 8-10$~\rm keV$, which is two times higher than that of the subhalo and the X-ray tail. The pressure is nearly uniform excluding southern part of the subhalo at two times of the truncation radius. We computed the gas mass within the truncation radius and the X-ray tail. While the gas fraction within the truncation radius is about 5 times smaller than that of regular groups, the gas mass in the subhalo and the X-ray tail to weak-lensing mass ratio is consistent with that of regular groups. Assuming the infall velocity, $2000~\rm km~s^{-1}$, the ram pressure is 1.4 times greater than gravitational force per unit area. Assuming the Kelvin-Helmholtz instabilities, the total lost mass is approximately $3\times10^{11}~M_{\odot}$. If this gas had originally been within the truncation radius, the gas mass fraction of the subhalo would have been comparable with those of regular groups before infalling to the Coma cluster.

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.

We have measured an annual parallax of the Mira variable R~Ursae~Majoris (R~UMa) with the VLBI exploration for Radio Astronomy (VERA). From the monitoring VLBI observations spanning about two years, we detected H$_2$O maser spots in the LSR velocities ranges from 37 to 42 km\,s$^{-1}$. We derived an annual parallax of 1.97$\pm$0.05\,mas, and it gives a corresponding distance of 508$\pm$13\,pc. The VLBI maps revealed 72 maser spots distributed in $\sim$110 au area around an expected stellar position. Circumstellar kinematics of the maser spots were also revealed by subtracting a systemic motion in the Hipparcos catalog from proper motions of each maser spots derived from our VLBI observations. Infrared photometry is also conducted to measure a $K$ band apparent magnitude, and we obtained a mean magnitude of $m_K$ = 1.19$\pm$0.02\,mag. Using the trigonometric distance, the $m_K$ is converted to a $K$ band absolute magnitude of $M_K = -$7.34$\pm$0.06\,mag. This result gives a much more accurate absolute magnitude of R~UMa than previously provided. We solved a zero-point of $M_K - \log P$ relation for the Galactic Mira variables and obtained a relation of $M_K = -$3.52 $\log P$ + (1.09 $\pm$ 0.14). Other long period variables including red supergiants, whose distances were determined from astrometric VLBI, were also compiled to explore the different sequences of $M_K - \log P$ relation.

We perform a weak-lensing study of the nearby cool-core galaxy clusters, Hydra A ($z=0.0538$) and A478 ($z=0.0881$), of which brightest cluster galaxies (BCGs) host powerful activities of active galactic nuclei (AGNs). For each cluster, the observed tangential shear profile is well described either by a single Navarro--Frenk--White model or a two-component model including the BCG as an unresolved point mass. For A478, we determine the BCG and its host-halo masses from a joint fit to weak-lensing and stellar photometry measurements. We find that the choice of initial mass functions (IMFs) can introduce a factor of two uncertainty in the BCG mass, whereas the BCG host halo mass is well constrained by data. We perform a joint analysis of weak-lensing and stellar kinematics data available for the Hydra A cluster, which allows us to constrain the central mass profile without assuming specific IMFs. We find that the central mass profile ($r<300$kpc) determined from the joint analysis is in excellent agreement with those from independent measurements, including dynamical masses estimated from the cold gas disk component, X-ray hydrostatic total mass estimates, and the central stellar mass estimated based on the Salpeter IMF. The observed dark-matter fraction around the BCG for Hydra A is found to be smaller than those predicted by adiabatic contraction models, suggesting the importance of other physical processes, such as the the AGN feedback and/or dissipationless mergers.

Looking at the region connecting two clusters is a promising way to identify and study the Warm-Hot Intergalactic Medium. Observations show that the spectrum of the bridge between A3556 and A3558 has a stronger soft X-ray emission than the nearby region. Suzaku observations could not discriminate the origin of the extra emission. In this work we analyze a dedicated Chandra observation of the same target to identify point sources and characterize the background emission in the bridge. We find that the count number of the point sources is much higher than average field population (using CDFS~4~Ms as a reference). Moreover, the shape of the cumulative distribution resembles that of galaxy distribution suggesting that the point sources are galaxies in a filament. The Suzaku extra emission is well explained by the high abundance of point sources identified by Chandra. Furthermore, we used optical/IR observations of point sources in the same field to estimate the density of the putative filament as rho~150 rho_b$, below Suzaku sensitivity.

We observed three massive subhalos in the Coma cluster with {\it Suzaku}. These subhalos, labeled "ID 1", "ID 2", and "ID 32", were detected with a weak-lensing survey using the Subaru/Suprime-Cam (Okabe et al. 2014a), and are located at the projected distances of 1.4 $r_{500}$, 1.2 $r_{500}$, and 1.6 $r_{500}$ from the center of the Coma cluster, respectively. The subhalo "ID 1" has a compact X-ray excess emission close to the center of the weak-lensing mass contour, and the gas mass to weak-lensing mass ratio is about 0.001. The temperature of the emission is about 3 keV, which is slightly lower than that of the surrounding intracluster medium (ICM) and that expected for the temperature vs. mass relation of clusters of galaxies. The subhalo "ID 32" shows an excess emission whose peak is shifted toward the opposite direction from the center of the Coma cluster. The gas mass to weak-lensing mass ratio is also about 0.001, which is significantly smaller than regular galaxy groups. The temperature of the excess is about 0.5 keV and significantly lower than that of the surrounding ICM and far from the temperature vs. mass relation of clusters. However, there is no significant excess X-ray emission in the "ID 2" subhalo. Assuming an infall velocity of about 2000 $\rm km~s^{-1}$, at the border of the excess X-ray emission, the ram pressures for "ID 1" and "ID 32" are comparable to the gravitational restoring force per area. We also studied the effect of the Kelvin-Helmholtz instability to strip the gas. Although we found X-ray clumps associated with the weak-lensing subhalos, their X-ray luminosities are much lower than the total ICM luminosity in the cluster outskirts.

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions developed by the Institute of Space and Astronautical Science (ISAS), with a planned launch in 2015. The ASTRO-H mission is equipped with a suite of sensitive instruments with the highest energy resolution ever achieved at E > 3 keV and a wide energy range spanning four decades in energy from soft X-rays to gamma-rays. The simultaneous broad band pass, coupled with the high spectral resolution of Delta E < 7 eV of the micro-calorimeter, will enable a wide variety of important science themes to be pursued. ASTRO-H is expected to provide breakthrough results in scientific areas as diverse as the large-scale structure of the Universe and its evolution, the behavior of matter in the gravitational strong field regime, the physical conditions in sites of cosmic-ray acceleration, and the distribution of dark matter in galaxy clusters at different redshifts.

The next generation X-ray observatory ASTRO-H will open up a new dimension in the study of galaxy clusters by achieving for the first time the spectral resolution required to measure velocities of the intracluster plasma, and extending at the same time the spectral coverage to energies well beyond 10 keV. This white paper provides an overview of the capabilities of ASTRO-H for exploring gas motions in galaxy clusters including their cosmological implications, the physics of AGN feedback, dynamics of cluster mergers as well as associated high-energy processes, chemical enrichment of the intracluster medium, and the nature of missing baryons and unidentified dark matter.

In asteroseismology, the observed time series often suffers from incomplete time coverage due to gaps. The presence of periodic gaps may generate spurious peaks in the power spectrum that limit the analysis of the data. Various methods have been developed to deal with gaps in time series data. However, it is still important to improve these methods to be able to extract all the possible information contained in the data. In this paper, we propose a new approach to handle the problem, the so-called inpainting method. This technique, based on a sparsity prior, enables to judiciously fill-in the gaps in the data, preserving the asteroseismic signal, as far as possible. The impact of the observational window function is reduced and the interpretation of the power spectrum is simplified. This method is applied both on ground and space-based data. It appears that the inpainting technique improves the oscillation modes detection and estimation. Additionally, it can be used to study very long time series of many stars because its computation is very fast. For a time series of 50 days of CoRoT-like data, it allows a speed-up factor of 1000, if compared to methods of the same accuracy.

We report properties of the intracluster medium (ICM) in Abell~1246 to the virial radius ($r_{200}$) and further outside as observed with Suzaku. The ICM emission is clearly detected to $r_{200}$, and we derive profiles of electron temperature, density, entropy, and cluster mass based on the spectral analysis. The temperature shows variation from $\sim 7$ keV at the central region to $\sim 2.5$ keV around $r_{200}$. The total mass in $r_{500}$ is $(4.3 \pm 0.4) \times 10^{14}~M_{\odot}$, assuming hydrostatic equilibrium. At $r>r_{500}$, the hydrostatic mass starts to decline and we, therefore, employ the total mass within $r_{200}$ based on weak-lens mass profile obtained from a sample of lower mass clusters. This yields the gas mass fraction at $r_{200}$ consistent with the cosmic baryon fraction, i.e. $\sim 17$%. The entropy profile indicates a flatter slope than that of the numerical simulation, particularly in $r>r_{500}$. These tendencies are similar to those of other clusters observed with Suzaku. We detect no significant ICM emission outside of $r_{200}$, and $2\sigma$ upper limits of redshifted OVII and OVIII line intensities are constrained to be less than 2.9 and $5.6\times 10^{-7}$ photons cm$^{-2}$ s$^{-1}$ arcmin$^{-2}$, respectively. The OVII line upper limit indicates $n_{\rm H}< 4.7\times 10^{-5}$ cm$^{-3}$ ($Z/0.2~Z_{\odot}$)$^{-1/2}$ ($L/20~{\rm Mpc}$)$^{-1/2}$, which corresponds to an overdensity, $\delta<160$ ($Z/0.2~Z_{\odot}$)$^{-1/2}$ ($L/20~{\rm Mpc}$)$^{-1/2}$.

We conduct a joint X-ray and weak-lensing study of four relaxed galaxy clusters (Hydra A, A478, A1689 and A1835) observed by both Suzaku and Subaru out to virial radii, with an aim to understand recently-discovered unexpected feature of the ICM in cluster outskirts. We show that the average hydrostatic-to-lensing total mass ratio for the four clusters decreases from \sim 70% to \sim 40% as the overdensity contrast decreases from 500 to the virial value.The average gas mass fraction from lensing total mass estimates increases with cluster radius and agrees with the cosmic mean baryon fraction within the virial radius, whereas the X-ray-based gas fraction considerably exceeds the cosmic values due to underestimation of the hydrostatic mass. We also develop a new advanced method for determining normalized cluster radial profiles for multiple X-ray observables by simultaneously taking into account both their radial dependence and multivariate scaling relations with weak-lensing masses. Although the four clusters span a range of halo mass, concentration, X-ray luminosity and redshift, we find that the gas entropy, pressure, temperature and density profiles are all remarkably self-similar when scaled with the lensing M_200 mass and r_200 radius.The entropy monotonically increases out to \sim 0.5r_200 following the accretion shock heating model K(r)\propto r^1.1, and flattens at \simgt 0.5r_200.The universality of the scaled entropy profiles indicates that the thermalization mechanism over the entire cluster region (>0.1r_200) is controlled by gravitation in a common to all clusters, although the heating efficiency in the outskirts needs to be modified from the standard law.The bivariate scaling functions of the gas density and temperature reveal that the flattening of the outskirts entropy profile is caused by the steepening of the temperature, rather than the flattening of the gas density.

We conducted phase referencing VLBI observations of the Mira variable T~Lepus (T~Lep) using VERA, from 2003 to 2006. The distance to the source was determined from its annual parallax which was measured to be 3.06$\pm$0.04 mas, corresponding to a distance of 327$\pm$4\,pc. Our observations revealed the distribution and internal kinematics of H$_2$O masers in T~Lep, and we derived a source systemic motions of 14.60$\pm$0.50 mas yr$^{-1}$ and $-$35.43$\pm$0.79 mas yr$^{-1}$ in right ascension and declination, respectively. We also determined a LSR velocity of $V_\mathrm{LSR}^{\ast} = -$27.63 km s$^{-1}$. Comparison of our result with an image recently obtained from the VLTI infrared interferometer reveals a linear scale picture of the circumstellar structure of T~Lep. Analysis of the source systemic motion in the Galacto-centric coordinate frame indicates a large peculiar motion, which is consistent with the general characteristics of AGB stars. This source makes a contribution to the calibration of the period-luminosity relation of Galactic Mira variables. From the compilation of data for nearby Mira variables found in the literature, whose distances were derived from astrometric VLBI observations, we have calibrated the Galactic Mira period-luminosity relation to a high degree of accuracy.

The central cD galaxy of the Hydra A cluster has one of the most powerful active galactic nuclei (AGNs) in the nearby Universe (z <0.2). We report on the discovery of a dust lane in the cD galaxy using Subaru telescope. The i'-band image shows the existence of a dark band of the size of 3.6"x 0.7" (4 kpc x 0.8 kpc), which appears to be quite similar to the dust lane observed in Centaurus A. The morphology indicates that the cold disk that seen as the dust lane is almost edge-on and rotates around the AGN. Since the minor axis of the dust lane is nearly parallel to the radio jets emerging from the AGN, the disk is probably feeding its gas into the central black hole. From the absorption, we estimate the hydrogen column density of the lane is N_H=2.0 x 10^21 cm^-2, and the mass of the disk is ~8 x 10^7 M_sun. The column density is consistent with constraints obtained from Chandra X-ray observations. The age of the disk is >~ 4 x 10^7 yr. The position angle of the disk and the galaxy's photometric axis are misaligned, which may imply that the cold gas in the disk is brought via galaxy mergers. Our observations may indicate that the supply of cold gas by galaxy mergers is required for the most intensive feedback from AGNs.

We present the results from Suzaku observations of the Perseus galaxy cluster, which is relatively close, the brightest in the X-ray sky and a relaxed object with a cool core. A number of exposures of central regions and offset pointing with the X-ray Imaging Spectrometer cover a region within radii of 20'-30'. The central data are used to evaluate the instrumental energy-scale calibration with accuracy confirmed to within around 300 km/s, by the spatial and temporal variation of the instruments. These deep and well-calibrated data are used to measure X-ray redshifts of the intracluster medium. A hint of gas bulk motion, with radial velocity of about -(150-300) km/s, relative to the main system was found at 2-4 arcmin (45-90kpc) west of the cluster center, where an X-ray excess and a cold front were found previously. No other velocity structure was discovered. Over spatial scales of 50-100kpc and within 200kpc radii of the center, the gas-radial-velocity variation is below 300 km/s, while over scales of 400 kpc within 600 kpc radii, the variation is below 600 km/s. These X-ray redshift distributions are compared spatially with those of optical member galaxies for the first time in galaxy clusters. Based on X-ray line widths gas turbulent velocities within these regions are also constrained within 1000-3000 km/s. These results of gas dynamics in the core and larger scales in association with cluster merger activities are discussed and future potential of high-energy resolution spectroscopy with ASTRO-H is considered.

We studied the distributions of metal abundances and metal-mass-to-light ratios in the intracluster medium (ICM) of four galaxy groups, MKW 4, HCG 62, the NGC 1550 group, and the NGC 5044 group, out to 0.5 r180 observed with Suzaku. The Fe abundance decreases with radius, and about 0.2-0.4 solar beyond 0.1 r180. At a given radius in units of r180, the Fe abundance in the ICM of the four galaxy groups were consistent or smaller than those of clusters of galaxies. The Mg/Fe and Si/Fe ratios in the ICM are nearly constant at the solar ratio out to 0.5 r180. We also studied systematic uncertainties in the derived metal abundances comparing the results from two versions of atomic data for astrophysicists (ATOMDB) and single- and two temperature model fits. Since the metals have been synthesized in galaxies, we collected K-band luminosities of galaxies from Two Micron All Sky Survey catalogue (2MASS) and calculated the integrated iron-mass-to-light-ratios (IMLR), or the ratios of the iron mass in the ICM to light from stars in galaxies. The groups with smaller gas mass to light ratios have smaller IMLR values and the IMLR inversely correlated with the entropy excess. Based on these abundance features, we discussed the past history of metal enrichment process in groups of galaxies.

The cluster of galaxies MS 1512.4+3647 (z=0.372) was observed with Suzaku for 270 ks. Besides the Fe abundance, the abundances of Mg, Si, S, and Ni are separately determined for the first time in a medium redshift cluster (z>0.3). The derived abundance pattern of MS 1512.4+3647 is consistent with those of nearby clusters, suggesting that the system has similar contributions from supernovae (SNe) Ia and SNe II to nearby clusters. The number ratio of SNe II to SNe Ia is $\sim3$. The estimated total numbers of both SNe II and SNe Ia against gas mass indicate similar correlation with those for the nearby clusters. The abundance results of MS 1512.4+3647 is consistent with the standard scenario that the SN II rate history roughly follows the star-formation history which has a peak at 1<z<2 and then declines by about one order of magnitude toward $z\sim0$. The similar number of SNe Ia to the nearby clusters suggests that the SN Ia rate declines steeply from z=0.37 to z=0 and/or SN Ia explosions occurred predominantly at larger redshifts.

Rotation is thought to drive cyclic magnetic activity in the Sun and Sun-like stars. Stellar dynamos, however, are poorly understood owing to the scarcity of observations of rotation and magnetic fields in stars. Here, inferences are drawn on the internal rotation of a distant Sun-like star by studying its global modes of oscillation. We report asteroseismic constraints imposed on the rotation rate and the inclination of the spin axis of the Sun-like star HD 52265, a principal target observed by the CoRoT satellite that is known to host a planetary companion. These seismic inferences are remarkably consistent with an independent spectroscopic observation (rotational line broadening) and with the observed rotation period of star spots. Furthermore, asteroseismology constrains the mass of exoplanet HD 52265b. Under the standard assumption that the stellar spin axis and the axis of the planetary orbit coincide, the minimum spectroscopic mass of the planet can be converted into a true mass of 1.85 (+0.52,-0.42) M_Jupiter, which implies that it is a planet, not a brown dwarf.

We investigate the observational constraints on the inhomogeneous big-bang nucleosynthesis that Matsuura et al. suggested the possibility of the heavy element production beyond ${}^7$Li in the early universe. From the observational constraints on light elements of ${}^4$He and D, possible regions are found on the plane of the volume fraction of the high density region against the ratio between high-and low-density regions. In these allowed regions, we have confirmed that the heavy elements beyond Ni can be produced appreciably, where $p$- and/or $r$-process elements are produced well simultaneously.

Based on a scenario of the inhomogeneous big-bang nucleosynthesis (IBBN), we investigate the detailed nucleosynthesis that includes the production of heavy elements beyond Li-7. From the observational constraints on light elements of He4 and D for the baryon-to-photon ratio given by WMAP, possible regions found on the plane of the volume fraction of the high density region against the ratio between high- and low-density regions. In these allowed regions, we have confirmed that the heavy elements beyond Fe can be produced appreciably, where p- and/or r-process elements are produced well simultaneously compared to the solar system abundances. We suggest that recent observational signals such as He4 overabundance in globular clusters and high metallicity abundances in quasars could be partly due to the results of IBBN. Possible implications are given for the formation of the first generation stars

We present the analysis of 21 bright X-ray knots in the Cassiopeia A supernova remnant from observations spanning 10 yr. We performed a comprehensive set of measurements to reveal the kinematic and thermal state of the plasma in each knot, using a combined analysis of two high energy resolution High Energy Transmission Grating (HETG) and four medium energy resolution Advanced CCD Imaging Spectrometer (ACIS) sets of spectra. The ACIS electron temperature estimates agree with the HETG-derived values for approximately half of the knots studied, yielding one of the first comparisons between high resolution temperature estimates and ACIS-derived temperatures. We did not observe the expected spectral evolutionpredicted from the ionization age and density estimates for each knotin all but three of the knots studied. The incompatibility of these measurements with our assumptions has led us to propose a dissociated ejecta model, with the metals unmixed inside the knots, which could place strong constraints on supernova mixing models.

We present results of four-pointing Suzaku X-ray observations (total ~200 ks) of the intracluster medium (ICM) in the Abell 1835 galaxy cluster (kT ~ 8 keV, z = 0.253) out to the virial radius (r_vir ~ 2.9 Mpc) and beyond. Faint X-ray emission from the ICM out to r_vir is detected. The temperature gradually decreases with radius from ~8 keV in the inner region to ~2 keV at r_vir. The entropy profile is shown to flatten beyond r_500, in disagreement with the r_1.1 dependence predicted from the accretion shock heating model. The thermal pressure profile in the range 0.3r_500 < r < r_vir agrees well with that obtained from the stacked Sunyaev-Zel'dovich effect observations with the Planck satellite. The hydrostatic mass profile in the cluster outskirts (r_500 < r < r_vir) falls well short of the weak lensing one derived from Subaru/Suprime-Cam observations, showing an unphysical decrease with radius. The gas mass fraction at r_vir defined with the lensing total mass agrees with the cosmic baryon fraction from the WMAP 7-year data. All these results indicate, rather than the gas-clumping effect, that the bulk of the ICM in the cluster outskirts is far from hydrostatic equilibrium and infalling matter retained some of its kinetic energy. Finally, combining with our recent Suzaku and lensing analysis of Abell 1689, a cluster of similar mass, temperature, and redshift, we show that the cluster temperature distribution in the outskirts is significantly correlated with the galaxy density field in the surrounding large-scale environment at (1-2)r_vir.

We studied the distributions of Si, Fe, and Ni in the intracluster medium (ICM) of the Coma cluster, one of the largest clusters in the nearby universe, using XMM-Newton data up to 0.5 r180 and Suzaku data of the central region up to 0.16 r180. Using the flux ratios of Ly alpha of H-like Si and 7.8 keV blend to K alpha of He-like Fe, the abundance ratios of Si to Fe and Ni to Fe of the ICM were derived using APEC model v2.0.1. The Si/Fe ratio in the ICM of the Coma cluster shows no radial gradient. The emission weighted averages of the Si/Fe ratio in the ICM within 0.0--0.2 r180, 0.2--0.5 r180, and 0.0--0.5 r180 are 0.97 +- 0.11, 1.05 +- 0.36 and 0.99 +- 0.13, respectively, in solar units using the solar abundance of Lodders (2003). These values are close to those of smaller clusters and groups of galaxies. Using the Suzaku data of the central region, the derived Ni/Fe ratio of the ICM is 0.6--1.5 in solar units, according to the same solar abundance table. The systematic difference in the derived abundance ratios by different plasma codes are about 10%. Therefore, for the ICM in the Coma cluster, the abundance pattern of Si, Fe, and Ni is consistent with the same mixture of the yields of SN II and SN Ia in our Galaxy. Within 0.5 r180}, the cumulative iron-mass-to-light ratio increases with radius, and its radial profile is similar to those of relaxed smaller clusters with cD galaxies at their center. Considering the observed Si/Fe ratio, the cumulative metal-mass-to-light ratios at 0.5 r180 are compared with theoretical expectations.

We analyzed XMM-Newton data of the Perseus cluster out to $\sim$1 Mpc, or approximately half the virial radius. Using the flux ratios of Lyalpha lines of H-like Si and S to Kalpha line of He-like Fe, the abundance ratios of Si/Fe and S/Fe of the intracluster medium (ICM) were derived using the APEC plasma code v2.0.1. The temperature dependence of the line ratio limits the systematic uncertainty in the derived abundance ratio. The Si/Fe and S/Fe in the ICM of the Perseus cluster show no radial gradient. The emission-weighted averages of the Si/Fe and S/Fe ratios outside the cool core are 0.91 +- 0.08 and 0.93 +- 0.10, respectively, in solar units according to the solar abundance table of Lodders (2003). These ratios indicate that most Fe was synthesized by supernovae Ia. We collected K-band luminosities of galaxies and calculated the ratio of Fe and Si mass in the ICM to K-band luminosity, iron-mass-to-light ratio (IMLR) and silicon-mass-to-light ratio (SMLR). Within $\sim$1 Mpc, the cumulative IMLR and SMLR increase with radius. Using Suzaku data for the northwest and east directions, we also calculated the IMLR out to $\sim$ 1.8 Mpc, or about the virial radius. We constrained the SMLR out to this radius and discussed the slope of the initial mass function of stars in the cluster. Using the cumulative IMLR profile, we discuss the past supernova Ia rate.

By implementing widely-used equations of state (EOS) from Lattimer & Swesty (LS) and H. Shen et al. (SHEN) in core-collapse supernova simulations, we explore possible impacts of these EOS on the post-bounce dynamics prior to the onset of neutrino-driven explosions. Our spherically symmetric (1D) and axially symmetric (2D) models are based on neutrino radiation hydrodynamics including spectral transport, which is solved by the isotropic diffusion source approximation. We confirm that in 1D simulations neutrino-driven explosions cannot be obtained for any of the employed EOS. Impacts of the EOS on the post-bounce hydrodynamics are more clearly visible in 2D simulations. In 2D models of a 15 M_sun progenitor using the LS EOS, the stalled bounce shock expands to increasingly larger radii, which is not the case using the SHEN EOS. Keeping in mind that the omission of the energy drain by heavy-lepton neutrinos in the present scheme could facilitate explosions, we find that 2D models of an 11.2 M_sun progenitor produce neutrino-driven explosions for all the EOS under investigation. Models using the LS EOS are slightly more energetic compared to those with the SHEN EOS. The more efficient neutrino heating in the LS models coincides with a higher electron antineutrino luminosity and a larger mass that is enclosed within the gain region. The models based on the LS EOS also show a more vigorous and aspherical downflow of accreting matter to the surface of the protoneutron star (PNS). The accretion pattern is essential for the production and strength of outgoing pressure waves, that can push in turn the shock to larger radii and provide more favorable conditions for the explosion. [abbreviated]

We present analyses to determine the fundamental parameters of the Galaxy based on VLBI astrometry of 52 Galactic maser sources obtained with VERA, VLBA and EVN. We model the Galaxy's structure with a set of parameters including the Galaxy center distance R_0, the angular rotation velocity at the LSR Omega_0, mean peculiar motion of the sources with respect to Galactic rotation (U_src, V_src, W_src), rotation-curve shape index, and the V component of the Solar peculiar motions V_sun. Based on a Markov chain Monte Carlo method, we find that the Galaxy center distance is constrained at a 5% level to be R_0 = 8.05 +/- 0.45 kpc, where the error bar includes both statistical and systematic errors. We also find that the two components of the source peculiar motion U_src and W_src are fairly small compared to the Galactic rotation velocity, being U_src = 1.0 +/- 1.5 km/s and W_src = -1.4 +/- 1.2 km/s. Also, the rotation curve shape is found to be basically flat between Galacto-centric radii of 4 and 13 kpc. On the other hand, we find a linear relation between V_src and V_sun as V_src = V_sun -19 (+/- 2) km/s, suggesting that the value of V_src is fully dependent on the adopted value of V_sun. Regarding the rotation speed in the vicinity of the Sun, we also find a strong correlation between Omega_0 and V_sun. We find that the angular velocity of the Sun, Omega_sun, which is defined as Omega_sun = Omega_0 + V_sun/R_0, can be well constrained with the best estimate of Omega_sun = 31.09 +/- 0.78 km/s/kpc. This corresponds to Theta_0 = 238 +/- 14 km/s if one adopts the above value of R_0 and recent determination of V_sun ~ 12 km/s.

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions initiated by the Institute of Space and Astronautical Science (ISAS). ASTRO-H will investigate the physics of the high-energy universe via a suite of four instruments, covering a very wide energy range, from 0.3 keV to 600 keV. These instruments include a high-resolution, high-throughput spectrometer sensitive over 0.3-2 keV with high spectral resolution of Delta E < 7 eV, enabled by a micro-calorimeter array located in the focal plane of thin-foil X-ray optics; hard X-ray imaging spectrometers covering 5-80 keV, located in the focal plane of multilayer-coated, focusing hard X-ray mirrors; a wide-field imaging spectrometer sensitive over 0.4-12 keV, with an X-ray CCD camera in the focal plane of a soft X-ray telescope; and a non-focusing Compton-camera type soft gamma-ray detector, sensitive in the 40-600 keV band. The simultaneous broad bandpass, coupled with high spectral resolution, will enable the pursuit of a wide variety of important science themes.

We report Suzaku observations of the northern half of the Hydra A cluster out to ~1.4 Mpc, reaching the virial radius. This is the first Suzaku observations of a medium-size (kT ~3 keV) cluster out to the virial radius. Two observations were conducted, north-west and north-east offsets, which continue in a filament direction and a void direction of the large-scale structure of the Universe, respectively. The X-ray emission and distribution of galaxies elongate in the filament direction. The temperature profiles in the two directions are mostly consistent with each other within the error bars and drop to 1.5 keV at 1.5 r_500. As observed by Suzaku in hot clusters, the entropy profile becomes flatter beyond r_500, in disagreement with the r^1.1 relationship that is expected from accretion shock heating models. When scaled with the average intracluster medium (ICM) temperature, the entropy profiles of clusters observed with Suzaku are universal and do not depend on system mass. The hydrostatic mass values in the void and filament directions are in good agreement, and the Navarro, Frenk, and White universal mass profile represents the hydrostatic mass distribution up to ~ 2 r_500. Beyond r_500, the ratio of gas mass to hydrostatic mass exceeds the result of the Wilkinson microwave anisotropy probe, and at r_100, these ratios in the filament and void directions reach 0.4 and 0.3, respectively. We discuss possible deviations from hydrostatic equilibrium at cluster outskirts. We derived radial profiles of the gasmass- to-light ratio and iron-mass-to-light ratio out to the virial radius. Within r_500, the iron-mass-to-light ratio of the Hydra A cluster was compared with those in other clusters observed with Suzaku.

We report the first observations of metal distributions in the intracluster medium of the Centaurus cluster up to ~0.17 r_180 with Suzaku. Radial profiles of the O, Mg, Si, S, Ar, Ca, and Fe were determined at the outer region of the cluster, and their variations appear to be similar to each other. Within the cool core region (r<0.045 r_ 180), all the metal distributions sharply increased toward the center. In the central region (r<0.015 r_180), the abundances of Si, S, Ar, Ca, and Fe were 1.5--1.8 solar, while those of O and Mg were approximately 1 solar. The derived abundance ratios of O and Mg to Fe were slightly lower than those of a set of other clusters. In contrast, the calculated mass-to-light ratios (MLRs) for O, Mg, and Fe were larger than those of the other clusters. For the outer region of the cool core (r>0.07 r_180), all the abundances were almost constant at 0.5 solar. The derived MLRs were comparable to those of the other clusters. This suggests that the cD galaxy of the Centaurus cluster efficiently supplies more Fe than the other clusters.

We present the preliminary results of a frequency and line-profile analysis of the CoRoT gamma Dor candidate HD171834. The data consist of 149 days of CoRoT light curves and a ground-based dataset of more than 1400 high-resolution spectra, obtained with six different instruments. Low-amplitude frequencies between 0 and 5 c/d, dominated by a frequency near 0.96 c/d and several of its harmonics, are detected. These findings suggest that HD171834 is not a mere gamma Dor pulsator and that stellar activity plays an important role in its variable behaviour. Based on CoRoT space data and on ground-based observations with ESO Telescopes at the La Silla Observatory under the ESO Large Programmes ESO LP 178.D-0361 and ESO LP 182.D-0356 (FEROS/2.2m and HARPS/3.6m), and data collected with FOCES/2.2m at the Centro Astronomico Hispano Aleman at Calar Alto, SOPHIE/1.93m at Observatoire de Haute Provence, FIES/NOT at Observatorio del Roque de los Muchachos, and HERCULES/1.0m at Mount John University Observatory.

We investigated the spatial correlation between the arrival directions of the highest energy cosmic rays (HECRs) detected by the Akeno Giant Air Shower Array (AGASA) with energies above $4 \times 10^{19}$ eV and the positions of galaxies observed by the Sloan Digital Sky Survey (SDSS) within $z = 0.024$. We systematically tested the dependence of the correlation on the redshift ranges and properties of the galaxies, i.e., absolute luminosity, color, and morphology, to understand where HECR sources are and what objects are HECR sources. In the systematic survey, we found potential signals of the positive correlation at small angular scale ($<10^{\circ}$) with the (non-penalized) chance probability less than 5% in intermediate redshift ranges. Then, we estimated penalized probabilities to compensate the trial effects of angular scan, and found that the strongest correlation is produced by early-type galaxies in $0.012 \leq z < 0.018$ at 90% C.L. The possible origin of HECRs which correlating galaxies imply is also discussed.

Suzaku observed a central region and five offset regions within 0.2 r180 in the Fornax cluster, a nearby poor cluster, and XMM-Newton mapped the cluster with 15 pointings out to 0.3 r180. The distributions of O, Mg, Si, S, and Fe in the intracluster medium (ICM) were studied with Suzaku, and those of Fe and temperature were studied with XMM. The temperature of the ICM gradually decreases with radius from 1.3 keV at 0.04 r180 to 1 keV at 0.2-0.3 r180. If the new solar abundances of Lodders et al. (2003) and a single-temperature plasma model are adopted, O, Mg, Si, S, and Fe show similar abundances: 0.4-0.6 solar within 0.02-0.2 r180. This Fe abundance is similar to those at 0.1-0.2 r180 in rich clusters and other groups of galaxies. At 0.2-0.3 r180, the Fe abundance becomes 0.2-0.3 solar. A two-temperature plasma model yields ICM abundances that are higher by a factor of 1.2-1.5, but gives similar abundance ratios among O, Mg, Si, S, and Fe. The northern region has a lower ICM temperature and higher brightness and Fe abundance, whereas the southern region has a higher ICM temperature and lower brightness and Fe abundance. These results indicate that the cD galaxy may have traveled from the north because of recent dynamical evolution. The cumulative oxygen- and iron-mass-to-light ratios within 0.3 r180 are more than an order of magnitude lower than those of rich clusters and some relaxed groups of galaxies. Past dynamical evolution might have hindered the strong concentration of hot gas in the Fornax cluster's central region. Scatter in the IMLR and similarity in the element abundances in the ICM of groups and clusters of galaxies indicate early metal synthesis.

Suzaku performed observations of 3 regions in and around the Shapley supercluster: a region located between A3558 and A3556, at ~0.9 times the virial radii of both clusters, and two other regions at 1{\deg}and 4{\deg}away from the first pointing. The 4{\deg}-offset observation was used to evaluate the Galactic foreground emission. We did not detect significant redshifted Oxygen emission lines (O VII and O VIII) in the spectra of all three pointings, after subtracting the contribution of foreground and background emission. An upper limit for the redshifted O VIII Ka line intensity of the warm-hot intergalactic medium (WHIM) is 1.5 \times 10^-7 photons s^-1 cm^-2 arcmin^-2, which corresponds to an overdensity of ~380 (Z/0.1 Z_solar)^{-1/2} (L/3 Mpc)^{-1/2}, assuming T=3\times10^6 K. We found excess continuum emission in the 1{\deg}-offset and on-filament regions, represented by thermal models with kT ~1 keV and ~2 keV, respectively. The redshifts of both 0 and that of the supercluster (0.048) are consistent with the observed spectra. The ~1 keV emission can be also fitted with Ne-rich Galactic (zero redshift) thin thermal emission. Radial intensity profile of 2 keV component suggests contribution from A3558 and A3556, but with significant steepening of the intensity slope in the outer region of A3558. Finally, we summarized the previous Suzaku search for the WHIM and discussed the feasibility of constraining the WHIM. An overdensity of < 400 can be detectable using O VII and O VIII emission lines in a range of 1.4\times10^6 K < T < 5\times10^6 K or a continuum emission in a relatively high temperature range T > 5\times10^6 K with the Suzaku XIS. The non detection with Suzaku suggests that typical line-of-sight average overdensity is < 400.

The Coma cluster was observed with an X-ray Imaging Spectrometer (XIS) onboard Suzaku in six pointings, including the central X-ray peak region, the 14'west offset region, 30'and 34'north-west offset regions, and 44'and 60'south-west offset regions. Owing to its lower background level, Suzaku has better sensitivity to Fe K-alpha lines than other satellites. Using precise Fe line measurements, we studied the temperature structure, possible bulk motions, and iron abundance distributions in the intracluster medium (ICM). The observed spectra were well-represented by a single-temperature model, and two- or three- temperature model did not improve chi-square substantially. The temperature, derived from K-alpha line ratios of H-like and He-like Fe, agree with those derived from the single-temperature model. Because the line ratio is a steep function of temperature, the consistency supports the accuracy of temperature measurements conducted with Suzaku. Within the 34'region, the redshift derived from the central energy of the He-like Fe line is consistent with that from optical observations, within a calibration error of 18 eV or 818 km/s in the line of sight. This value is smaller than the sound velocity of ICM, which is 1500 km/s. The central energy of Fe lines at the 44'offset region around the NGC 4839 subcluster is also consistent with those within the 34'region. These results on the temperature and velocity structure suggest that the core of the cluster is in a relaxed state, and non-thermal electrons relevant to the radio halo are accelerated by intracluster turbulence rather than large-scale shocks. Fe abundance is almost constant at 0.4 solar within the 34'region, and decreases with radius. This value is slightly lower than those of other clusters, which means the gas have been mixed well during a past merger associated with the growth of the cluster.

We observed outer regions of a bright cluster of galaxies A2142 with Suzaku. Temperature and brightness structures were measured out to the virial radius ($r_{200}$) with good sensitivity. We confirmed the temperature drop from 9 keV around the cluster center to about 3.5 keV at $r_{200}$, with the density profile well approximated by the $\beta$ model with $\beta = 0.85$. Within $0.4\r_{200}$, the entropy profile agrees with $r^{1.1}$, as predicted by the accretion shock model. The entropy slope becomes flatter in the outer region and negative around $r_{200}$. These features suggest that the intracluster medium in the outer region is out of thermal equilibrium. Since the relaxation timescale of electron-ion Coulomb collision is expected to be longer than the elapsed time after shock heating at $r_{200}$, one plausible reason of the low entropy is the low electron temperature compared to that of ions. Other possible explanations would be gas clumpiness, turbulence and bulk motions of ICM\@. We also searched for a warm-hot intergalactic medium around $r_{200}$ and set an upper limit on the oxygen line intensity. Assuming a line-of-sight depth of 2 Mpc and oxygen abundance of 0.1 solar, the upper limit of an overdensity is calculated to be 280 or 380, depending on the foreground assumption.

By performing a series of one- and two-dimensional (1-, 2D) hydrodynamic simulations with spectral neutrino transport, we study possible impacts of collective neutrino oscillations on the dynamics of core-collapse supernovae. To model the spectral swapping which is one of the possible outcome of the collective neutrino oscillations, we parametrize the onset time when the spectral swap begins, the radius where the spectral swap occurs, and the threshold energy above which the spectral interchange between heavy-lepton neutrinos and electron/anti-electron neutrinos takes place, respectively. By doing so, we systematically study how the neutrino heating enhanced by the spectral swapping could affect the shock evolution as well as the matter ejection. We also investigate the progenitor dependence using a suite of progenitor models (13, 15, 20, and 25 $M_\odot$). We find that there is a critical heating rate induced by the spectral swapping to trigger explosions, which significantly differs between the progenitors. The critical heating rate is generally smaller for 2D than 1D due to the multidimensionality that enhances the neutrino heating efficiency. For the progenitors employed in this paper, the final remnant masses are estimated to range in 1.1-1.5$M_\odot$. For our 2D model of the $15M_\odot$ progenitor, we find a set of the oscillation parameters that could account for strong supernova explosions ($\sim 10^{51}$ erg), simultaneously leaving behind the remnant mass close to $\sim 1.4 M_\odot$.

The star HD 52265 is a G0V metal-rich exoplanet-host star observed in the seismology field of the CoRoT space telescope from November 2008 to March 2009. The satellite collected 117 days of high-precision photometric data on this star, showing that it presents solar-like oscillations. HD 52265 was also observed in spectroscopy with the Narval spectrograph at the same epoch. We characterise HD 52265 using both spectroscopic and seismic data. The fundamental stellar parameters of HD 52265 were derived with the semi-automatic software VWA, and the projected rotational velocity was estimated by fitting synthetic profiles to isolated lines in the observed spectrum. The parameters of the observed p modes were determined with a maximum-likelihood estimation. We performed a global fit of the oscillation spectrum, over about ten radial orders, for degrees l=0 to 2. We also derived the properties of the granulation, and analysed a signature of the rotation induced by the photospheric magnetic activity. Precise determinations of fundamental parameters have been obtained: Teff = 6100 +- 60 K, log g = 4.35 +- 0.09, [M/H] = 0.19 +- 0.05, as well as vsini = 3.6 +0.3 -1.0 km/s. We have measured a mean rotation period P_rot = 12.3 +- 0.15 days, and find a signature of differential rotation. The frequencies of 31 modes are reported in the range 1500-2550 micro-Hz. The large separation exhibits a clear modulation around the mean value <Dnu> = 98.3 +- 0.1 micro-Hz. Mode widths vary with frequency along an S-shape with a clear local maximum around 1800 micro-Hz. We deduce lifetimes ranging between 0.5 and 3 days for these modes. Finally, we find a maximal bolometric amplitude of about 3.96 +- 0.24 ppm for radial modes.

ORIGIN is a proposal for the M3 mission call of ESA aimed at the study of metal creation from the epoch of cosmic dawn. Using high-spectral resolution in the soft X-ray band, ORIGIN will be able to identify the physical conditions of all abundant elements between C and Ni to red-shifts of z=10, and beyond. The mission will answer questions such as: When were the first metals created? How does the cosmic metal content evolve? Where do most of the metals reside in the Universe? What is the role of metals in structure formation and evolution? To reach out to the early Universe ORIGIN will use Gamma-Ray Bursts (GRBs) to study their local environments in their host galaxies. This requires the capability to slew the satellite in less than a minute to the GRB location. By studying the chemical composition and properties of clusters of galaxies we can extend the range of exploration to lower redshifts (z ~ 0.2). For this task we need a high-resolution spectral imaging instrument with a large field of view. Using the same instrument, we can also study the so far only partially detected baryons in the Warm-Hot Intergalactic Medium (WHIM). The less dense part of the WHIM will be studied using absorption lines at low redshift in the spectra for GRBs.

We have measured the annual parallax of the H2O maser source associated with an infrared dark cloud MSXDC G034.43+00.24 from the observations with VERA (VLBI Exploration of Radio Astrometry). The parallax is 0.643 +/- 0.049 mas, corresponding to the distance of 1.56 +0.12/-0.11 kpc. This value is less than the half of the previous kinematic distance of 3.7 kpc. We revise the core mass estimates of MSXDC G034.43+00.24, based on virial masses, LTE masses and dust masses and show that the core masses decrease from the previous estimations of ~1000 Mo to hundreds of Mo. The spectral type derived from the luminosity also changes from O9.5 to B1 in the case of MM1. This spectral type is still consistent with that of the massive star. The radial velocity derived from the flat rotation model is smaller than the observed velocity, which corresponds to the peculiar motion of ~40 km/s in the line-of-sight direction.

We present the results of astrometric observations with VERA toward the H2O maser sources in IRAS 05137+3919, which is thought to be located in the far outer Galaxy. We have derived the parallax of \pi = 0.086 +/- 0.027 mas, which corresponds to the source distance of D=11.6+5.3-2.8 kpc. Although the parallax measurement is only 3-sigma level and thus the distance uncertainty is considerably large, we can strongly constrain the minimum distance to this source, locating the source at the distance from the Sun greater than 8.3 kpc (or 16.7 kpc from the Galaxy's center) at 90% confidence level. Our results provide an astrometric confirmation that this source is located in the far outer Galaxy beyond 15 kpc from the Galaxy center, indicating that IRAS 05137+3919 is one of the most distant star-forming regions from the Galaxy center.

We report on the astrometric observations of the 22 GHz H2O masers in the high mass star-forming region G5.89-0.39 with VERA (VLBI Exploration of Radio Astrometry). Newly derived distance of 1.28^{+0.09}_{-0.08} kpc is the most precise and significantly nearer than previous values. We revised physical parameters and reconsidered nature of G5.89-0.39 based on the new distance as follows. (1) The ionizing star of the ultra compact (UC) HII region is a late O-type (O8 - 8.5) zero age main sequence (ZAMS) star, consistent with previously established limits based on its infrared spectral line emission. (2) Crescent-like maser alignment at the position of the O type ZAMS star may trace accretion disk (or its remnant), which suggests that the star is still young and before complete evaporation of circumstellar materials. (3) Although the revised mass for the east-west outflow has been reduced, it still quite large (100 Msun) which indicates that a significant fraction of the mass is entrained material and that the dynamical age significantly underestimates the actual outflow age. Our newly-derived distance emphasizes that G5.89-0.39 is one of the nearest targets to investigate ongoing high-mass star formation and evolution in a compact cluster containing a young O-type star.

Since the detection of the asymptotic properties of the dipole gravity modes in the Sun, the quest to find individual gravity modes has continued. An extensive and deeper analysis of 14 years of continuous GOLF/SoHO observational data, unveils the presence of a pattern of peaks that could be interpreted as individual dipole gravity modes in the frequency range between 60 and 140 microHz, with amplitudes compatible with the latest theoretical predictions. By collapsing the power spectrum we have obtained a quite constant splitting for these patterns in comparison to regions where no g modes were expected. Moreover, the same technique applied to simultaneous VIRGO/SoHO data unveils some common signals between the power spectra of both instruments. Thus, we are able to identify and characterize individual g modes with their central frequencies, amplitudes and splittings allowing to do seismic inversions of the rotation profile inside the solar core. These results open a new light on the physics and dynamics of the solar deep core.

The helioseismic Global Oscillation at Low Frequency (GOLF) and the Variability of solar Irradiance and Gravity Oscillations (VIRGO) instruments onboard SoHO, have been observing the Sun continuously for the last 14 years. In this preliminary work, we characterize the acoustic modes over the entire p-mode range in both, Doppler velocity and luminosity, with a special care for the low-frequency modes taking advantage of the stability and the high duty cycle of space observations.

We studied the high temperature plasma in the direction of the Sculptor supercluster at z=0.108 with Suzaku. Suzaku carried out four observations in the supercluster: namely, A2811, A2811 offset, A2804, A2801 regions in 2005 Nov.--Dec., including the regions beyond the virial radii of these clusters. The study needed precise background estimation because the measured intensity of the redshifted lines, especially those from oxygen, were strongly affected by the the Galactic emission. The spectra taken in the regions outside of the virial radii of the member clusters were used as the background which included both the Galactic and Cosmic X-ray Background (CXB) components. We also used the background data which were taken near the Sculptor supercluster. Temperature and metal abundance profiles were determined to the virial radii of the member clusters, and then we searched for the oxygen line emission in the region outside of the virial radii of the clusters. As a result, the temperature of the clusters decreased toward the virial radii, and the spectral fits for the filament region did not require extra component other than the Galactic and CXB components. We constrained the intensities of O VII and O VIII lines to be less than 8.1 and 5.1 photons cm^-2 s^-1 arcmin^-2, respectively, as 2-sigma upper limits. The intensity of O VII indicates n_H < 1.6e-5 cm^-3 (Z/0.1 Z_solar)^-1/2 (L/25 Mpc)^-1/2, which corresponds to an over density, delta < 60 (Z/0.1 Z_solar)^-1/2 (L/25 Mpc)^-1/2.

The primary science goal of the Kepler Mission is to provide a census of exoplanets in the solar neighborhood, including the identification and characterization of habitable Earth-like planets. The asteroseismic capabilities of the mission are being used to determine precise radii and ages for the target stars from their solar-like oscillations. Chaplin et al. (2010) published observations of three bright G-type stars, which were monitored during the first 33.5 days of science operations. One of these stars, the subgiant KIC 11026764, exhibits a characteristic pattern of oscillation frequencies suggesting that it has evolved significantly. We have derived asteroseismic estimates of the properties of KIC 11026764 from Kepler photometry combined with ground-based spectroscopic data. We present the results of detailed modeling for this star, employing a variety of independent codes and analyses that attempt to match the asteroseismic and spectroscopic constraints simultaneously. We determine both the radius and the age of KIC 11026764 with a precision near 1%, and an accuracy near 2% for the radius and 15% for the age. Continued observations of this star promise to reveal additional oscillation frequencies that will further improve the determination of its fundamental properties.

By performing axisymmetric hydrodynamic simulations of core-collapse supernovae with spectral neutrino transport based on the isotropic diffusion source approximation scheme, we support the assumption that the neutrino-heating mechanism aided by the standing accretion shock instability and convection can initiate an explosion of a 13 $M_{\odot}$ star. Our results show that bipolar explosions are more likely to be associated with models which include rotation. We point out that models, which form a north-south symmetric bipolar explosion, can lead to larger explosion energies than for the corresponding unipolar explosions.

We studied the temperature and metal abundance distributions of the intra-cluster medium (ICM) in a group of galaxies NGC 1550 observed with Suzaku. The NGC 1550 is classified as a fossil group, which have few bright member galaxies except for the central galaxy. Thus, such a type of galaxy is important to investigate how the metals are enriched to the ICM. With the Suzaku XIS instruments, we directly measured not only Si, S, and Fe lines but also O and Mg lines and obtained those abundances to an outer region of ~0.5 r_180 for the first time, and confirmed that the metals in the ICM of such a fossil group are indeed extending to a large radius. We found steeper gradients for Mg, Si, S, and Fe abundances, while O showed almost flat abundance distribution. Abundance ratios of alpha-elements to Fe were similar to those of the other groups and poor clusters. We calculated the number ratio of type II to type Ia supernovae for the ICM enrichment to be 2.9 +- 0.5 within 0.1 r_180, and the value was consistent with those for the other groups and poor clusters observed with Suzaku. We also calculated metal mass-to-light ratios (MLRs) for Fe, O and Mg with B-band and K-band luminosities of the member galaxies of NGC 1550. The derived MLRs were comparable to those of NGC 5044 group in the r<0.1 r_180 region, while those of NGC 1550 are slightly higher than those of NGC 5044 in the outer region.

The propagation trajectories of the highest energy cosmic rays (HECRs) are deflected by not only intergalactic magnetic field but also Galactic magnetic field (GMF). These magnetic fields can weaken the positive correlation between the arrival directions of HECRs and the positions of their sources. In order to explore the effect of GMF on the expected correlation, we simulate the arrival distribution of protons with the energy above $6 \times 10^{19}$ eV taking several GMF models into account, and then test the correlation between the protons and their sources assumed in the simulation. The dependence of the correlation signals on GMF models are also investigated. The correlation can be observed by accumulating $\sim 200$ protons in a half hemisphere. Typical angular scale at which the positive signal of the correlation is maximized depends on the spiral component of GMF models. That angular scale is $\sim 5^o$ for bisymmetric spiral (BS) GMF models and $\sim 7^o$ for axisymmetric spiral (AS) GMF models if the number density of HECR sources, $n_s$, is $\sim 10^{-4}$ Mpc$^{-3}$. An additional vertical (dipole) component of GMF affects these angular scale by $0.5^o$ - $1^o$. The difference between the correlation signal for the BS models and that for the AS models is prominent in the northern sky. Significance of the positive correlation depends on source distribution. The probability that the number of simulated HECR events correlating with sources is smaller than the number of random events correlating with the same sources by chance is much less than $10^{-3}$ ($\sim 3\sigma$) in almost all the source distributions with $n_s = 10^{-4}$ Mpc$^{-3}$ under 200 protons detection, but $\sim 10\%$ of source distributions predicts the chance probability more than $10^{-3}$ in the AS GMF model. In addition, we also briefly discuss the effect of GMF for heavy-nuclei dominated composition.

We report on the first asteroseismic analysis of solar-type stars observed by Kepler. Observations of three G-type stars, made at one-minute cadence during the first 33.5d of science operations, reveal high signal-to-noise solar-like oscillation spectra in all three stars: About 20 modes of oscillation can clearly be distinguished in each star. We discuss the appearance of the oscillation spectra, including the presence of a possible signature of faculae, and the presence of mixed modes in one of the three stars.

In Helio- and asteroseismology, it is important to have continuous, uninterrupted, data sets. However, seismic observations usually contain gaps and we need to take them into account. In particular, if the gaps are not randomly distributed, they will produce a peak and a series of harmonics in the periodogram that will destroy the stellar information. An interpolation of the data can be good solution for this problem. In this paper we have studied an interpolation method based on the so-called 'inpainting' algorithms. To check the algorithm, we used both VIRGO and CoRoT satellite data to which we applied a realistic artificial window of a real CoRoT observing run to introduce gaps. Next we compared the results with the original, non-windowed data. Therefore, we were able to optimize the algorithm by minimizing the difference between the power spectrum density of the data with gaps and the complete time series. In general, we find that the power spectrum of the inpainted time series is very similar to the original, unperturbed one. Seismic inferences obtained after interpolating the data are the same as in the original case.

The CoRoT mission is in its third year of observation and the data from the second long run in the galactic centre direction are being analysed. The solar-like oscillating stars that have been observed up to now have given some interesting results, specially concerning the amplitudes that are lower than predicted. We present here the results from the analysis of the star HD 170987.The goal of this research work is to characterise the global parameters of HD 170987. We look for global seismic parameters such as the mean large separation, maximum amplitude of the modes, and surface rotation because the signal-to-noise ratio in the observations do not allow us to measure individual modes. We also want to retrieve the stellar parameters of the star and its chemical composition.We have studied the chemical composition of the star using ground-based observations performed with the NARVAL spectrograph. We have used several methods to calculate the global parameters from the acoustic oscillations based on CoRoT data. The light curve of the star has been interpolated using inpainting algorithms to reduce the effect of data gaps. We find power excess related to p modes in the range [400 - 1200]muHz with a mean large separation of 55.2+-0.8muHz with a probability above 95% that increases to 55.9 +-0.2muHz in a higher frequency range [500 - 1250] muHz and a rejection level of 1%. A hint of the variation of this quantity with frequency is also found. The rotation period of the star is estimated to be around 4.3 days with an inclination axis of i=50 deg +20/-13. We measure a bolometric amplitude per radial mode in a range [2.4 - 2.9] ppm around 1000 muHz. Finally, using a grid of models, we estimate the stellar mass, M=1.43+-0.05 Msun, the radius, R=1.96+-0.046 Rsun, and the age ~2.4 Gyr.

Structure of hot dense matter at subnuclear densities is investigated by quantum molecular dynamics (QMD) simulations. We analyze nucleon distributions and nuclear shapes using two-point correlation functions and Minkowski functionals to determine the phase-separation line and to classify the phase of nuclear matter in terms of the nuclear structure. Obtained phase diagrams show that the density of the phase boundaries between the different nuclear structures decreases with increasing temperature due to the thermal expansion of nuclear matter region. The critical temperature for the phase separation is $\agt 6$ MeV for the proton fraction $x=0.5$ and $\agt 5$ MeV for $x=0.3$. Our result suggests the existence of "spongelike" phases with negative Euler characteristic in addition to the simple "pasta" phases in supernova cores until $T \alt 3$ MeV.

We examine the model dependence of the phase diagram of inhomogeneous nulcear matter in supernova cores using the quantum molecular dynamics (QMD). Inhomogeneous matter includes crystallized matter with nonspherical nuclei - "pasta" phases - and the liquid-gas phase separating nuclear matter. Major differences between the phase diagrams of the QMD models can be explained by the energy of pure neutron matter at low densities and the saturation density of asymmetric nuclear matter. We show the density dependence of the symmetry energy is also useful to understand uncertainties of the phase diagram. We point out that, for typical nuclear models, the mass fraction of the pasta phases in the later stage of the collapsing cores is higher than 10-20 %.

We present X-ray observations of the northern outskirts of the relaxed galaxy cluster A1413 with Suzaku, whose XIS instrument has the low intrinsic background needed to make measurements of these low surface brightness regions. We excise 15 point sources superimposed on the image above a flux of $1\times 10^{-14}$ \fluxunit (2--10keV) using XMM-Newton and Suzaku images of the cluster. We quantify all known systematic errors as part of our analysis, and show our statistical errors encompasses them for the most part. Our results extend previous measurements with Chandra and XMM-Newton, and show a significant temperature drop to about 3keV at the virial radius, $r_{200}$. Our entropy profile in the outer region ($> 0.5 r_{200}$) joins smoothly onto that of XMM-Newton, and shows a flatter slope compared with simple models, similar to a few other clusters observed at the virial radius. The integrated mass of the cluster at the virial radius is approximately $7.5\times10^{14}M_{\odot}$ and varies by about 30% depending on the particular method used to measure it.

We present preliminary asteroseismic results from Kepler on three G-type stars. The observations, made at one-minute cadence during the first 33.5d of science operations, reveal high signal-to-noise solar-like oscillation spectra in all three stars: About 20 modes of oscillation may be clearly distinguished in each star. We discuss the appearance of the oscillation spectra, use the frequencies and frequency separations to provide first results on the radii, masses and ages of the stars, and comment in the light of these results on prospects for inference on other solar-type stars that Kepler will observe.

Over the past 15 years, our knowledge of the interior of the Sun has tremendously progressed by the use of helioseismic measurements. However, to go further in our understanding of the solar core, we need to measure gravity (g) modes. Thanks to the high quality of the Doppler-velocity signal measured by GOLF/SoHO, it has been possible to unveil the signature of the asymptotic properties of the solar g modes, thus obtaining a hint of the rotation rate in the core. However, the quest for the detection of individual g modes is not yet over. In this work, we apply the latest theoretical developments to guide our research using GOLF velocity time series. In contrary to what was thought till now, we are maybe starting to identify individual low-frequency g modes...

In supernova cores, nuclear "pasta" phases such as triangular lattice of rod-like nuclei and layered structure of slab-like nuclei are considered to exist. However, it is still unclear whether or not they are actually formed in collapsing supernova cores. Using {\it ab-initio} numerical simulations called the Quantum Molecular Dynamics (QMD), we here solve this problem by demonstrating that a lattice of rod-like nuclei is formed from a bcc lattice by compression. We also find that, in the transition process, the system undergoes zigzag configuration of elongated nuclei, which are formed by a fusion of two original spherical nuclei.

The paper by Sato and Kobayashi in 1977 studied the cosmological effects of a massive neutrino and obtained constraints on its properties. This paper initiated many studies to use cosmology as a laboratory of particle physics or to use particle physics to explore the very early universe.

We investigate correlation between the arrival directions of ultra-high-energy cosmic rays (UHECRs) and the large-scale structure (LSS) of the Universe by using statistical quantities which can find the angular scale of the correlation. The Infrared Astronomical Satellite Point Source Redshift Survey (IRAS PSCz) catalog of galaxies is adopted for LSS. We find a positive correlation of the highest energy events detected by the Pierre Auger Observatory (PAO) with the IRAS galaxies inside $z=0.018$ within the angular scale of $\sim 15^{\circ}$. This positive correlation observed in the southern sky implies that a significant fraction of the highest energy events comes from nearby extragalactic objects. We also analyze the data of the Akeno Giant Air Shower Array (AGASA) which observed the northern hemisphere, but the obvious signals of positive correlation with the galaxy distribution are not found. Since the exposure of the AGASA is smaller than the PAO, the cross-correlation in the northern sky should be tested using a larger number of events detected in the future. We also discuss the correlation using the all-sky combined data sets of both the PAO and AGASA, and find a significant correlation within $\sim 8^{\circ}$. These angular scales can constrain several models of intergalactic magnetic field. These cross-correlation signals can be well reproduced by a source model in which the distribution of UHECR sources is related to the IRAS galaxies.

The Suzaku X-ray satellite observed the nearby spiral galaxy NGC 4258 for a total good exposure time of 100 ks. We present an analysis of the Suzaku XIS data, in which we confirm that the 0.5--2 keV spectra of the interstellar medium (ISM) are well-represented by a two-temperature model. The cool and hot ISM temperatures are 0.23+-0.02 and 0.59 +-0.01 keV, respectively. Suzaku's excellent spectral sensitivity enables us to measure the metal abundances of O, Ne, Mg, Si and Fe of the ISM for the first time. The resultant abundance pattern of O, Mg, Si, and Fe is consistent with that of the new solar abundance table of Lodders (2003), rather than Anders & Grevesse (1989). This suggests that the metal enrichment processes of NGC 4258 and of our Galaxy are similar.

Context: Measurements of intracluster gas temperatures out to large radii are important for the use of clusters for precision cosmology and for studies of cluster physics. Previous attempts to measure robust temperatures at cluster virial radii failed. Aims: The goal of this work is to measure the temperature profile of the very relaxed galaxy cluster Abell 2204 out to large radii, possibly reaching the virial radius. Methods: Taking advantage of its low particle background due to its low-Earth orbit, Suzaku data are used to measure the outer temperature profile of Abell 2204. These data are combined with Chandra and XMM-Newton data of the same cluster in order to make the connection to the inner regions, unresolved by Suzaku, and to determine the smearing due to Suzaku's PSF. Results: The temperature profile of Abell 2204 is determined from 10 kpc to 1800 kpc, close to an estimate of r200 (the approximation to the virial radius). The temperature rises steeply from below 4 keV in the very center up to more than 8 keV in the intermediate range and then decreases again to about 4 keV at the largest radii. Varying the measured particle background normalization artificially by +-10 percent does not change the results significantly. Predictions for outer temperature profiles based on hydrodynamic simulations show good agreement. In particular, we find the observed temperature profile to be slightly steeper but consistent with a drop of a factor of 0.6 from 0.3 r200 to r200, as predicted by simulations. Conclusions: Temperature measurements up to the virial radius seem feasible with Suzaku, when a careful analysis of the different background components and the effects of the PSF is performed. The result obtained here indicates that numerical simulations capture the intracluster gas physics well in cluster outskirts.

There are two promising scenarios that explain the ankle, which is a dip in the spectrum of cosmic rays at $\sim 10^{19}$ eV. A scenario interprets the ankle as the transition from Galactic to extragalactic cosmic rays ({\it ankle-transition scenario}), while the other is that the dip caused by pair production on the cosmic microwave background radiation ({\it proton-dip scenario}). In this paper, we consider whether cosmogenic neutrinos can be a clue to judge which scenario is favored. We calculated the fluxes of cosmogenic neutrinos following these scenarios with plausible physical parameter sets, and found several important features as follows. First of all, the neutrino flux at $\sim 10^{20}$ eV becomes much higher in the ankle-transition scenario as long as the maximum energy of the cosmic rays at sources is sufficiently high. On the other hand, the neutrino spectrum has a characteristic peak at $\sim 10^{16}$ eV in the proton-dip scenario on the condition that extragalactic protons significantly contribute to the observed cosmic rays down to $10^{17}$ eV. Thus, we conclude cosmogenic neutrinos should give us a clue to judge which scenario is favored, unless these features are masked by the neutrino background coming from possible, powerful neutrino sources such as AGNs and GRBs. We also found an interesting feature that the neutrino flux at $\sim 10^{18}$ eV depends only on the cosmological evolution of the cosmic ray sources. That means cosmogenic neutrinos with the energy bring us information on the cosmological evolution of the sources of ultra-high energy cosmic rays. Finally, we compare the fluxes of cosmogenic neutrinos with the expected sensitivity curves of several neutrino detectors, and conclude the detection of cosmogenic neutrinos in the near future is promising.

We perform a series of two-dimensional magnetohydrodynamic core-collapse simulations of rapidly rotating and strongly magnetized massive stars. To study the properties of magnetic explosions for a longer time stretch of postbounce evolution, we develop a new code under the framework of special relativity including a realistic equation of state with a multiflavor neutrino leakage scheme. Our results show the generation of the magnetically-dominated jets in the two ways. One is launched just after the core-bounce in a prompt way and another is launched at $ \sim 100 $ ms after the stall of the prompt shock. We find that the shock-revival occurs when the magnetic pressure becomes strong, due to the field wrapping, enough to overwhelm the ram pressure of the accreting matter. The critical toroidal magnetic fields for the magnetic shock-revival are found to be universal of $\sim 10^{15}\mathrm{G}$ behind the jets. We point out that the time difference before the shock-revival has a strong correlation with the explosions energies. Our results suggest that the magnetically dominated jets are accompanied by the formation of the magnetars. Since the jets are mildly relativistic, we speculate that they might be the origin of some observed X-ray flashes.

Under the fossil field hypothesis of the origin of magnetar magnetic fields, the magnetar inherits its magnetic field from its progenitor. We show that during the supernova of such a progenitor, protons may be accelerated to \sim 10^4 GeV as the supernova shock propagates in the magnetic stellar envelope. Inelastic nuclear collisions of these protons produce a flash of high-energy neutrinos arriving a few hours after thermal (10 MeV) neutrinos. The neutrino flash is characterized by energies up to O(100) GeV and durations seconds to hours, depending on the progenitor: those from smaller Type Ibc progenitors are typically shorter in duration and reach higher energies compared to those from larger Type II progenitors. A Galactic Type Ib supernova leaving behind a magnetar remnant will yield up to \sim 160 neutrino induced muon events in Super-Kamiokande, and up to \sim 7000 in a km^3 class detector such as IceCube, providing a means of probing supernova models and the presence of strong magnetic fields in the stellar envelope.

The metallicity distribution in the intracluster medium of the NGC 5044 group was studied up to 0.3 r_180 using the XIS instrument on board the Suzaku satellite. Abundances of O, Mg, Si, S, and Fe were measured with high accuracy. The region within a radius of 0.05 r_180 from the center shows approximately solar abundances of Mg, Si, S, and Fe, while the O/Fe ratio is about 0.5--0.6 in solar units. In the outer region, the Fe abundance gradually drops to 0.3 solar. Radial abundance profiles of Mg, Si and S are similar to that of Fe, while that of O seems to be flatter. At r>0.05 r_180, the mass density profile of O differs from that of Fe, showing a shoulder-like structure that traces the luminosity density profile of galaxies. The mass-to-light ratios for O and Fe in NGC 5044 are one of the largest among groups of galaxies, but they are still smaller than those in rich clusters. These abundance features probably reflect the metal enrichment history of this relaxed group hosting a giant elliptical galaxy in the center.

We studied the temperature and abundance distributions of intra-cluster medium (ICM) in the Abell 262 cluster of galaxies observed with Suzaku. Abell 262 is a bright, nearby poor cluster with ICM temperature of ~2 keV, thus providing useful information about the connection of ICM properties between groups and clusters of galaxies. The observed spectrum of the central region was well-represented by two temperature models, and the spectra for the outer regions were described by single temperature model. With the XIS instrument, we directly measured not only Si, S, and Fe lines but also O and Mg lines and obtained those abundances to an outer region of ~0.43 r_180 for the first time. We found steep gradients for Mg, Si, S, and Fe abundances, while O showed almost flat abundance distribution. Abundance ratios of alpha-elements to Fe were found to be similar to those of the other clusters and groups. We calculated the number ratio of type II to type Ia supernovae for the ICM enrichment to be 3.0 +- 0.6 within 0.1 r_180, and the value was consistent with those for other clusters and groups. We also calculated metal mass-to-light ratios (MLRs) for Fe, O and Mg (IMLR, OMLR, MMLR) with B-band and K-band luminosities of the member galaxies of Abell 262. The derived MLRs were comparable to those for other clusters with kT = 3-4 keV.

Suzaku observation of the edge-on spiral galaxy NGC 4631 confirmed its X-ray halo extending out to about 10 kpc from the galactic disk. The XIS spectra yielded the temperature and metal abundance for the disk and the halo regions. The observed abundance pattern for O, Ne, Mg, Si and Fe is consistent with the metal yield from type II supernovae, with an O mass of about 1E6 Msun contained in the halo. These features imply that metal-rich gas produced by type II supernova is brought into the halo region very effectively, most likely through a galactic wind. Temperature and metal abundance may be affected by charge exchange and dust. An upper limit for the hard X-ray flux was obtained, corresponding to a magnetic field higher than 0.5 uG.

We present a distance measurement for the semiregular variable S Crateris (S Crt) based on its annual parallax. With the unique dual beam system of the VLBI Exploration for Radio Astrometry (VERA) telescopes, we measured the absolute proper motion of a water maser spot associated with S Crt, referred to the quasar J1147-0724 located at an angular separation of 1.23$^{\circ}$. In observations spanning nearly two years, we have detected the maser spot at the LSR velocity of 34.7 km s$^{-1}$, for which we measured the annual parallax of 2.33$\pm$0.13 mas corresponding to a distance of 430$^{+25}_{-23}$ pc. This measurement has an accuracy one order of magnitude better than the parallax measurements of HIPPARCOS. The angular distribution and three-dimensional velocity field of maser spots indicate a bipolar outflow with the flow axis along northeast-southwest direction. Using the distance and photospheric temperature, we estimate the stellar radius of S Crt and compare it with those of Mira variables.

We present results of multi-epoch VLBI observations with VERA (VLBI Exploration of Radio Astrometry) of the 22 GHz H$_{2}$O masers associated with a young stellar object (YSO) IRAS 22198+6336 in a dark cloud L1204G. Based on the phase-referencing VLBI astrometry, we derive an annual parallax of IRAS 22198+6336 to be 1.309$\pm$0.047 mas, corresponding to the distance of 764$\pm$27 pc from the Sun. Although the most principal error source of our astrometry is attributed to the internal structure of the maser spots, we successfully reduce the errors in the derived annual parallax by employing the position measurements for all of the 26 detected maser spots. Based on this result, we reanalyze the spectral energy distribution (SED) of IRAS 22198+6336 and find that the bolometric luminosity and total mass of IRAS 22198+6336 are 450$L_{\odot}$ and 7$M_{\odot}$, respectively. These values are consistent with an intermediate-mass YSO deeply embedded in the dense dust core, which has been proposed to be an intermediate-mass counterpart of a low-mass Class 0 source. In addition, we obtain absolute proper motions of the H$_{2}$O masers for the most blue-shifted components. We propose that the collimated jets aligned along the east-west direction are the most plausible explanation for the origin of the detected maser features.

We report astrometric observations of H2O masers around the red supergiant VY Canis Majoris (VY CMa) carried out with VLBI Exploration of Radio Astrometry (VERA). Based on astrometric monitoring for 13 months, we successfully measured a trigonometric parallax of 0.88 +/- 0.08 mas, corresponding to a distance of 1.14 +0.11/-0.09 kpc. This is the most accurate distance to VY CMa and the first one based on an annual parallax measurement. The luminosity of VY CMa has been overestimated due to a previously accepted distance. With our result, we re-estimate the luminosity of VY CMa to be (3 +/- 0.5) x 10^5 L_sun using the bolometric flux integrated over optical and IR wavelengths. This improved luminosity value makes location of VY CMa on the Hertzsprung-Russel (HR) diagram much closer to the theoretically allowable zone (i.e. the left side of the Hayashi track) than previous ones, though uncertainty in the effective temperature of the stellar surface still does not permit us to make a final conclusion.

We estimate the local number density of sources of ultra-high-energy cosmic rays (UHECRs) based on the statistical features of their arrival direction distribution. We calculate the arrival distributions of protons above $10^{19}$ eV taking into account their propagation process in the Galactic magnetic field and a structured intergalactic magnetic field, and statistically compare those with the observational result of the Pierre Auger Observatory. The anisotropy in the arrival distribution at the highest energies enables us to estimate the number density of UHECR sources as $\sim 10^{-4} {\rm Mpc}^{-3}$ assuming the persistent activity of UHECR sources. We compare the estimated number density of UHECR sources with the number densities of known astrophysical objects. This estimated number density is consistent with the number density of Fanaroff-Reily I galaxies. We also discuss the reproducability of the observed {\it isotropy} in the arrival distribution above $10^{19}$ eV. We find that the estimated source model cannot reproduce the observed isotropy. However, the observed isotropy can be reproduced with the number density of $10^{-2}$-$10^{-3} {\rm Mpc}^{-3}$. This fact indicates the existence of UHECR sources with a maximum acceleration energy of $\sim 10^{19}$ eV whose number density is an order of magnitude more than that injecting the highest energy cosmic rays.

This is a reply report to astro-ph/0604264. We studied heavy element production in high baryon density region in early universe astro-ph/0507439. However it is claimed in astro-ph/0604264 that small scale but high baryon density region contradicts observations for the light element abundance or in order not to contradict to observations high density region must be so small that it cannot affect the present heavy element abundance. In this paper we study big bang nucleosynthesis in high baryon density region and show that in certain parameter spaces it is possible to produce enough amount of heavy element without contradiction to CMB and light element observations.

We have used the Japanese VLBI array VERA to perform high-precision astrometry of an H2O maser source in the Galactic star-forming region NGC 281 West, which has been considered to be part of a 300-pc superbubble. We successfully detected a trigonometric parallax of 0.355+/-0.030 mas, corresponding to a source distance of 2.82+/-0.24 kpc. Our direct distance determination of NGC 281 has resolved the large distance discrepancy between previous photometric and kinematic studies; likely NGC 281 is in the far side of the Perseus spiral arm. The source distance as well as the absolute proper motions were used to demonstrate the 3D structure and expansion of the NGC 281 superbubble, ~650 pc in size parallel to the Galactic disk and with a shape slightly elongated along the disk or spherical, but not vertically elongated, indicating the superbubble expansion may be confined to the disk. We estimate the expansion velocity of the superbubble as ~20 km/s both perpendicular to and parallel to the Galactic disk with a consistent timescale of ~20 Myr.

We study the effects of rotation on the neutrino emission from Population III (Pop III) stars by performing a series of two-dimensional rotational collapse simulations of Pop III stellar cores. Our results show that rotation enhances the neutrino luminosities and the average energies of emitted neutrinos. This is because the thermalized inner core, which is the dominant neutrino source from Pop III stars, can be enlarged, due to rotational flattening, enough to extend the inner core outside the neutrinospheres. This is in sharp contrast to the case of spherical collapse, in which the case of inner core shrinks deeper inside the neutrinospheres before black hole formation, which hinders the efficient neutrino emission. In the case of rotational core-collapse, the emitted neutrino energies are found to become larger in the vicinity near the pole than the ones near the equatorial plane. These factors make the emergent neutrino spectrum broader and harder than the spherical collapse case. By computing the overall neutrino signals produced by the ensemble of individual rotating Pop III stars, we find that the amplitudes of the relic neutrinos, depending on their star formation rates, can dominate over the contributions from ordinary core-collapse supernovae below a few MeV. A detection of this signal could be an important tool to probe star formation history in the early universe.

Temperature and abundance distributions of the intra-cluster medium (ICM) in the NGC 507 group of galaxies were studied with Suzaku. Observed concentric annular spectra were well-represented by a two temperature model for ICM, and we found steeper abundance gradients for Mg, Si, S, and Fe compared with O in the central region. Abundance ratios of alpha-elements to iron were found to be similar to those in other groups and poor clusters. We calculated metal mass-to-light ratios for Fe, O and Mg (IMLR, OMLR, MMLR) for NGC 507, and values for different systems were compared. Hotter and richer systems tend to show higher values of IMLR, OMLR, and MMLR. OMLR and MMLR were measured to an outer region for the first time with Suzaku, while IMLR was consistent with that with ASCA. We also looked into 2-dimensional map of the hardness ratio, but found no significant deviation from the circular symmetry.

We investigate the deflections of UHE protons by Galactic magnetic field(GMF) using four conventional GMF models in order to discuss the positional correlation between the arrival distribution of UHECRs and their sources. UHE protons coming from the direction around the Galactic center are highly deflected above $8^{\circ}$ by the dipole magnetic field during their propagation in Galactic space. However, in bisymmetric spiral field models, there are directions with the deflection angle below $1^{\circ}$. One of these directions is toward Centaurus A, the nearest radio-loud active galactic nuclei that is one of possible candidates of UHECR sources. On the other hand, UHE protons arriving from the direction of the anti-Galactic center are less deflected, especially in bisymmetric spiral field models. Thus, the northern hemisphere, not including the Galactic center, is suitable for the studies of correlation with sources. The dependence on model parameters is also investigated. The deflection angles of UHE protons are dependent on the pitch angle of the spiral field. We also investigate distortion of the supergalactic plane by GMF. Since the distortion in the direction around Galactic center strongly depends on the GMF model, we can obtain information on GMF around Galactic center if Pierre Auger Observatory finds the significant positional correlation around the supergalactic plane.

We present results of 120 ks observation of a compact group of galaxies HCG~62 ($z=0.0145$) with Suzaku XIS and HXD-PIN\@. The XIS spectra for four annular regions were fitted with two temperature {\it vapec} model with variable abundance, combined with the foreground Galactic component. The Galactic component was constrained to have a common surface brightness among the four annuli, and two temperature {\it apec} model was preferred to single temperature model. We confirmed the multi-temperature nature of the intra-group medium reported with Chandra and XMM-Newton, with a doughnut-like high temperature ring at radii 3.3--6.5$'$ in a hardness image. We found Mg, Si, S, and Fe abundances to be fairly robust. We examined the possible ``high-abundance arc'' at $\sim 2'$ southwest from the center, however Suzaku data did not confirm it. We suspect that it is a misidentification of an excess hot component in this region as the Fe line. Careful background study showed no positive detection of the extended hard X-rays previously reported with ASCA, in 5--12 keV with XIS and 12--40 keV with HXD-PIN, although our upper limit did not exclude the ASCA result. There is an indication that the X-ray intensity in $r<3.3'$ region is $70\pm 19$% higher than the nominal CXB level (5--12 keV), and Chandra and Suzaku data suggest that most of this excess could be due to concentration of hard X-ray sources with an average photon index of $\Gamma=1.38\pm 0.06$. Cumulative mass of O, Fe and Mg in the group gas and the metal mass-to-light ratio were derived and compared with those in other groups. Possible role of AGN or galaxy mergers in this group is also discussed.

We discuss the possibility of accurately estimating the source number density of ultra-high-energy cosmic rays (UHECRs) using small-scale anisotropy in their arrival distribution. The arrival distribution has information on their source and source distribution. We calculate the propagation of UHE protons in a structured extragalactic magnetic field (EGMF) and simulate their arrival distribution at the Earth using our previously developed method. The source number density that can best reproduce observational results by Akeno Giant Air Shower Array is estimated at about $10^{-5} {\rm Mpc}^{-3}$ in a simple source model. Despite having large uncertainties of about one order of magnitude, due to small number of observed events in current status, we find that more detection of UHECRs in the Auger era can sufficiently decrease this so that the source number density can be more robustly estimated. 200 event observation above $4 \times 10^{19} {\rm eV}$ in a hemisphere can discriminate between $10^{-5}$ and $10^{-6} {\rm Mpc}^{-3}$. Number of events to discriminate between $10^{-4}$ and $10^{-5} {\rm Mpc}^{-3}$ is dependent on EGMF strength. We also discuss the same in another source model in this paper.

We investigate the possibility that near future observations of ultra-high-energy cosmic rays (UHECRs) can unveil their local source distribution, which reflects the observed local structures if their origins are astrophysical objects. In order to discuss this possibility, we calculate the arrival distribution of UHE protons taking into account their propagation process in intergalactic space i.e. energy losses and deflections by extragalactic magnetic field (EGMF). For a realistic simulation, we construct and adopt a model of a structured EGMF and UHECR source distribution, which reproduce the local structures actually observed around the Milky Way. The arrival distribution is compared statistically to their source distribution using correlation coefficient. We specially find that UHECRs above $10^{19.8}$eV are best indicators to decipher their source distribution within 100 Mpc, and detection of about 500 events on all the sky allows us to unveil the local structure of UHE universe for plausible EGMF strength and the source number density. This number of events can be detected by five years observation by Pierre Auger Observatory.

We report on the results of multi-epoch VLBI observations with VERA (VLBI Exploration of Radio Astrometry) of the 22 GHz H2O masers associated with the young stellar object SVS 13 in the NGC 1333 region. We have carried out phase-referencing VLBI astrometry and measured an annual parallax of the maser features in SVS 13 of 4.25+/-0.32 mas, corresponding to the distance of 235+/-18 pc from the Sun. Our result is consistent with a photometric distance of 220 pc previously reported. Even though the maser features were detectable only for 6 months, the present results provide the distance to NGC 1333 with much higher accuracy than photometric methods. The absolute positions and proper motions have been derived, revealing that the H2O masers with the LSR (local standard of rest) velocities of 7-8 km s-1 are most likely associated with VLA4A, which is a radio counterpart of SVS 13. The origin of the observed proper motions of the maser features are currently difficult to attribute to either the jet or the rotating circumstellar disk associated with VLA4A, which should be investigated through future high-resolution astrometric observations of VLA4A and other radio sources in NGC 1333.

We have performed high-precision astrometry of H2O maser sources in Galactic star forming region Sharpless 269 (S269) with VERA. We have successfully detected a trigonometric parallax of 189+/-8 micro-arcsec, corresponding to the source distance of 5.28 +0.24/-0.22 kpc. This is the smallest parallax ever measured, and the first one detected beyond 5 kpc. The source distance as well as proper motions are used to constrain the outer rotation curve of the Galaxy, demonstrating that the difference of rotation velocities at the Sun and at S269 (which is 13.1 kpc away from the Galaxy's center) is less than 3%. This gives the strongest constraint on the flatness of the outer rotation curve and provides a direct confirmation on the existence of large amount of dark matter in the Galaxy's outer disk.

We studied the properties of the intra-cluster medium (ICM) in two clusters of galaxies (AWM 7 and Abell 1060) and two groups (HCG 62 and NGC 507) with the X-ray Observatory Suzaku. Based on spatially resolved energy spectra, we measured for the first time precise cumulative ICM metal masses within 0.1 and ~0.3 r_{180}. Comparing our results with supernova nucleosynthesis models, the number ratio of type II (SNe II) to type Ia (SNe Ia) is estimated to be ~3.5, assuming the metal mass in the ICM is represented by the sum of products synthesized in SNe Ia and SNe II. Normalized by the K-band luminosities of present galaxies, and including the metals in stars, the integrated number of past SNe II explosions is estimated to be close to or somewhat higher than the star formation rate determined from Hubble Deep Field observations.

We carried out 3 observations of the cluster of galaxies AWM 7, for the central region and 20'-east and 20'-west offset regions, with Suzaku. Temperature and abundance profiles are measured out to 27'~ 570 /h_70 kpc, which corresponded to ~0.35 r_180. The temperature of the intra-cluster medium (ICM) slightly decreases from 3.8 keV at the center to 3.4 keV in ~0.35 r_180 region, indicating a flatter profile than those in other nearby clusters. Abundance ratio of Si to Fe is almost constant in our observation, while Mg to Fe ratio increases with radius from the cluster center. O to Fe ratio in the west region shows increase with radius, while that in the east region is almost flat, though the errors are relatively large. These features suggest that the enrichment process is significantly different between products of type II supernovae (O and Mg) and those by type Ia supernovae (Si and Fe). We also examined positional shift of the central energy of He-like Fe-Ka line, in search of possible rotation of the ICM. The 90% upper limit for the line-of-sight velocity difference was derived to be v ~ 2000 km/s, suggesting that the ellipticity of AWM 7 is rather caused by a recent directional infall of the gas along the large-scale filament.

We simulate the arrival distribution of ultra-high-energy (UHE) protons by following their propagation processes in several strengths of a structured extragalactic magnetic field (EGMF). Comparing our result to observational one by Akeno Giant Air Shower Array, we constrain the number density of UHE cosmic ray sources with the small-scale anisotropy. As a result, the source number density is $\sim 10^{-5} {\rm Mpc}^{-3}$ with uncertainty of about an order of magnitude due to the small number of observed events. This hardly depends on our structured EGMF strength. We also investigate future prospects for this approarch. The near future observations, such as Pierre Auger Observatory, can distinguish $10^{-6} {\rm Mpc}^{-3}$ accurately from the more source densities. Number of events to discriminate between $10^{-4}$ and $10^{-5} {\rm Mpc}^{-3}$ is dependent on the EGMF strength.

We estimate the gravitational wave (GW) background from Population III (Pop III) stars using the results from our hydrodynamic simulations (Suwa et al. 2007). We calculate gravitational waveforms from matter motions and anisotropic neutrino emissions for single Pop III stars. We find that the GW amplitudes from matter motions are dominant until just after bounce, but those from neutrinos dominate later on at frequencies below $\sim 10$ Hz in the GW spectrum. Computing the overall signal produced by the ensemble of such Pop III stars, we find that the resultant density parameter of the GW background peaks at the amplitude of $\Omega_\mathrm{GW}\sim 10^{-10}$ in the frequency interval $\sim 1-10$ Hz. We show that such signals, depending on the formation rate of Pop III stars, can be within the detection limits of future planned interferometers such as DECIGO and BBO in the frequency interval of $\sim 0.1-1$ Hz. Our results suggest that the detection of the GW background from Pop III stars can be an important tool to supply the information about the star formation history in the early universe.

We present the initial results of multi-epoch VLBI observations of the 22 GHz H2O masers in the Orion KL region with VERA (VLBI Exploration of Radio Astrometry). With the VERA dual-beam receiving system, we have carried out phase-referencing VLBI astrometry and successfully detected an annual parallax of Orion KL to be 2.29+/-0.10 mas, corresponding to the distance of 437+/-19 pc from the Sun. The distance to Orion KL is determined for the first time with the annual parallax method in these observations. Although this value is consistent with that of the previously reported, 480+/-80 pc, which is estimated from the statistical parallax method using proper motions and radial velocities of the H2O maser features, our new results provide the much more accurate value with an uncertainty of only 4%. In addition to the annual parallax, we have detected an absolute proper motion of the maser feature, suggesting an outflow motion powered by the radio source I along with the systematic motion of source I itself.

We report on absolute proper-motion measurements of an H2O maser source in the NGC 281 West molecular cloud, which is located ~320 pc above the Galactic plane and is associated with an HI loop extending from the Galactic plane. We have conducted multi-epoch phase-referencing observations of the maser source with VERA (VLBI Exploration of Radio Astrometry) over a monitoring period of 6 months since May 2006. We find that the H2O maser features in NGC 281 West are systematically moving toward the southwest and further away from the Galactic plane with a vertical velocity of ~20-30 km/s at its estimated distance of 2.2-3.5 kpc. Our new results provide the most direct evidence that the gas in the NGC 281 region on the HI loop was blown out from the Galactic plane, most likely in a superbubble driven by multiple or sequential supernova explosions in the Galactic plane.

Nuclear ``pasta'', nonspherical nuclei in dense matter, is predicted to occur in collapsing supernova cores. We show how pasta phases affect the neutrino transport cross section via weak neutral current using several nuclear models. This is the first calculation of the neutrino opacity of the phases with rod-like and slab-like nuclei taking account of finite temperature effects, which are well described by the quantum molecular dynamics. We also show that pasta phases can occupy 10--20% of the mass of supernova cores in the later stage of the collapse.

We perform a series of two-dimensional magnetorotational core-collapse simulations of Pop III stars. Changing the initial distributions of rotation and magnetic fields prior to collapse in a parametric manner, we compute 19 models. By so doing, we systematically investigate how rotation and magnetic fields affect the collapse dynamics and explore how the properties of the black-hole formations and neutrino emissions could be affected. As for the microphysics, we employ a realistic equation of state and approximate the neutrino transfer by a multiflavour leakage scheme. With these computations, we find that the jet-like explosions are obtained by the magnetodriven shock waves if the initial magnetic field is as large as $10^{12}$ G. We point out that the black-hole masses at the formation decrease with the initial field strength, on the other hand, increase with the initial rotation rates. As for the neutrino properties, we point out that the degree of the differential rotation plays an important role to determine which species of the neutrino luminosity is more dominant than the others. Furthermore, we find that the stronger magnetic fields make the peak neutrino luminosities smaller, because the magnetic pressure acts to halt the collapse in the central regions, leading to the suppression of the releasable gravitational binding energies.

We carried out observations of the central and 20' east offset regions of the cluster of galaxies Abell 1060 with Suzaku. Spatially resolved X-ray spectral analysis has revealed temperature and abundance profiles of Abell 1060 out to 27' ~ 380; /h_70 kpc, which corresponded to ~ 0.25; r_180. Temperature decrease of the intra cluster medium from 3.4 keV at the center to 2.2 keV in the outskirt region are clearly observed. Abundances of Si, S and Fe also decrease by more than 50% from the center to the outer, while Mg shows fairly constant abundance distribution at ~ 0.7 solar within r < 17'. O shows lower abundance of ~ 0.3 solar in the central region (r~ 6'), and indicates a similar feature with Mg, however it is sensitive to the estimated contribution of the Galactic components of kT_1 ~ 0.15 keV and kT_2 ~ 0.7 keV in the outer annuli (r ~ 13'). Systematic effects due to the point spread function tails, contamination on the XIS filters, instrumental background, cosmic and/or Galactic X-ray background, and the assumed solar abundance tables are carefully examined. Results on temperature and abundances of Si, S, and Fe are consistent with those derived by XMM-Newton at r < 13'. Formation and metal enrichment process of the cluster are discussed based on the present results.

We have developed a framework for the Monte-Carlo simulation of the X-Ray Telescopes (XRT) and the X-ray Imaging Spectrometers (XIS) onboard Suzaku, mainly for the scientific analysis of spatially and spectroscopically complex celestial sources. A photon-by-photon instrumental simulator is built on the ANL platform, which has been successfully used in ASCA data analysis. The simulator has a modular structure, in which the XRT simulation is based on a ray-tracing library, while the XIS simulation utilizes a spectral "Redistribution Matrix File" (RMF), generated separately by other tools. Instrumental characteristics and calibration results, e.g., XRT geometry, reflectivity, mutual alignments, thermal shield transmission, build-up of the contamination on the XIS optical blocking filters (OBF), are incorporated as completely as possible. Most of this information is available in the form of the FITS (Flexible Image Transport System) files in the standard calibration database (CALDB). This simulator can also be utilized to generate an "Ancillary Response File" (ARF), which describes the XRT response and the amount of OBF contamination. The ARF is dependent on the spatial distribution of the celestial target and the photon accumulation region on the detector, as well as observing conditions such as the observation date and satellite attitude. We describe principles of the simulator and the ARF generator, and demonstrate their performance in comparison with in-flight data.

The mass spectrum of the primordial black holes formed by density perturbation in the radiation-dominated era of the Randall-Sundrum type-2 cosmology is given. The spectrum coincides with standard four-dimensional one on large scales but the deviation is apparent on smaller scales. The mass spectrum is initially softer than standard four-dimensional one, while after accretion during the earliest era it becomes harder than that. We also show expected extragalactic diffuse photon background spectra varying the initial perturbation power-law power spectrum and give constraints on the blue spectra and/or the reheating temperature. The most recent observations on the small scale density perturbation from WMAP, SDSS and Lyman-\alpha Forest are used. What we get are interpreted as constraints on the smaller scale inflation on the brane connected to the larger one at the scale of Lyman-\alpha Forest. If we set the bulk curvature radius to be 0.1 mm and assume the reheating temperature is higher than 10^6 GeV, the scalar spectral index from the smaller scale inflation is constrained to be n \lesssim 1.3. Typically, the constraints are tighter than standard four-dimensional one, which is also revised by us using the most recent observations.

Two-dimensional hydrodynamic simulations are performed to investigate explosive nucleosynthesis in a collapsar using the model of MacFadyen and Woosley (1999). It is shown that 56Ni is not produced in the jet of the collapsar sufficiently to explain the observed amount of a hypernova when the duration of the explosion is \sim 10 sec, which is considered to be the typical timescale of explosion in the collapsar model. Even though a considerable amount of 56Ni is synthesized if all explosion energy is deposited initially, the opening angles of the jets become too wide to realize highly relativistic outflows and gamma-ray bursts in such a case. From these results, it is concluded that the origin of 56Ni in hypernovae associated with GRBs is not the explosive nucleosynthesis in the jet. We consider that the idea that the origin is the explosive nucleosynthesis in the accretion disk is more promising. We also show that the explosion becomes bi-polar naturally due to the effect of the deformed progenitor. This fact suggests that the 56Ni synthesized in the accretion disk and conveyed as outflows are blown along to the rotation axis, which will explain the line features of SN 1998bw and double peaked line features of SN 2003jd. Some fraction of the gamma-ray lines from 56Ni decays in the jet will appear without losing their energies because the jet becomes optically thin before a considerable amount of 56Ni decays as long as the jet is a relativistic flow. We show that abundance of nuclei whose mass number \sim 40 in the ejecta depends sensitively on the energy deposition rate. So it may be determined by observations of chemical composition in metal poor stars which model is the proper one as a model of a gamma-ray burst accompanied by a hypernova.

We perform two-dimensional, axisymmetric, magnetohydrodynamic simulations of the collapse of a rotating star of 40 Msun and in the light of the collapsar model of gamma-ray burst. Considering two distributions of angular momentum, up to \sim 10^{17} cm^2/s, and the uniform vertical magnetic field, we investigate the formation of an accretion disk around a black hole and the jet production near the hole. After material reaches to the black hole with the high angular momentum, the disk is formed inside a surface of weak shock. The disk becomes in a quasi-steady state for stars whose magnetic field is less than 10^{10} G before the collapse. We find that the jet can be driven by the magnetic fields even if the central core does not rotate as rapidly as previously assumed and outer layers of the star has sufficiently high angular momentum. The magnetic fields are chiefly amplified inside the disk due to the compression and the wrapping of the field. The fields inside the disk propagate to the polar region along the inner boundary near the black hole through the Alfv{\'e}n wave, and eventually drive the jet. The quasi-steady disk is not an advection-dominated disk but a neutrino cooling-dominated one. Mass accretion rates in the disks are greater than 0.01 Msun/sec with large fluctuations. The disk is transparent for neutrinos. The dense part of the disk, which locates near the hole, emits neutrino efficiently at a constant rate of < 8 \times 10^{51} erg/s. The neutrino luminosity is much smaller than those from supernovae after the neutrino burst.

We present a new astrophysical site of the big bang nucleosynthesis (BBN) that are very peculiar compared with the standard BBN. Some models of the baryogenesis suggest that very high baryon density regions were formed in the early universe. On the other hand, recent observations suggest that heavy elements already exist in high red-shifts and the origin of these elements become a big puzzle. Motivated by these, we investigate BBN in very high baryon density regions. BBN proceeds in proton-rich environment, which is known to be like the p-process. However, by taking very heavy nuclei into account, we find that BBN proceeds through both the p-process and the r-process simultaneously. P-nuclei such as 92Mo, 94Mo, 96Ru, 98Ru whose origin is not well known are also synthesized.

We report a current status of our radiation-magnetohydrodynamic code for the study of core-collapse supernovae. In this contribution, we discuss the accuracy of our newly developed numerical code by presenting the test problem in a static background model. We also present the application to the spherically symmetric core-collapse simulations. Since close comparison with the previously published models is made, we are now applying it for the study of magnetorotational core-collapse supernovae.

We present numerical simulations on propagation of Ultra-High Energy Cosmic Rays (UHECRs) above $10^{19}$ eV in a structured extragalactic magnetic field (EGMF) and simulate their arrival distributions at the earth. We use the IRAS PSCz catalogue in order to construct a model of the EGMF and source models of UHECRs, both of which reproduce the local structures observed around the Milky Way. We also consider modifications of UHECR arrival directions by the galactic magnetic field. We follow an inverse process of their propagation from the earth and record the trajectories. This enables us to calculate only trajectories of UHECRs arriving at the earth, which saves the CPU time. From these trajectories and our source models, we construct arrival distributions of UHECRs and calculate the harmonic amplitudes and the two point correlation functions of them. We estimate number density of sources which reproduces the Akeno Ground Air Shower Array (AGASA) observation best. As a result, we find that the most appropriate number density of the sources is $\sim 5 \times 10^{-6}$ Mpc$^{-3}$. This constrains the source candidates of UHECRs. We also demonstrate skymaps of their arrival distribution with the event number expected by future experiments and examine how the EGMF affects their arrival distribution. A main result is diffusion of clustering events which are obtained from calculations in the absence of the EGMF. This tendency allows us to reproduce the observed two point correlation function better.

We present results from an XMM-Newton observation of the head-tail radio galaxy IC 310 located in the southwest region of the Perseus cluster. The spectrum is well-fitted by an absorbed power-law model with a photon index of $2.50 \pm 0.02$ with no significant absorption excess. The X-ray image shows a point-like emission at IC 310 without any signs of a structure correlated with the radio halo tail. The temperature of the intracluster medium surrounding IC 310 declines as a function of distance from the cluster center, from $ kT \sim 6$ keV in the northeast corner of the field of view to about 3 keV in the southwest region. Although we do not find any sharp edges in the surface brightness profile, a brightness excess over a smooth $\beta$ model by about 20% is seen. The temperature also rises by about 10% in the same region. This indicates that the IC 310 region is a subcluster probably infalling into the Perseus cluster, and the gas in front of IC 310 towards the Perseus cluster is likely to be compressed by the large-scale motion, which supports the view that the IC 310 system is undergoing a merger.

Core-collapse supernovae are among the most energetic explosions in the universe marking the catastrophic end of massive stars. In spite of rigorous studies for several decades, we still don't understand the explosion mechanism completely. Since they are related to many astrophysical phenomena such as nucleosynthesis, gamma-ray bursts and acceleration of cosmic rays, understanding of their physics has been of wide interest to the astrophysical community. In this article, we review recent progress in the study of core-collapse supernovae focusing on the explosion mechanism, supernova neutrinos, and the gravitational waves. As for the explosion mechanism, we present a review paying particular attention to the roles of multidimensional aspects, such as convection, rotation, and magnetic fields, on the neutrino heating mechanism. Next, we discuss supernova neutrinos, which is a powerful tool to probe not only deep inside of the supernovae but also intrinsic properties of neutrinos. For this purpose, it is necessary to understand neutrino oscillation which has been established recently by a lot of experiments. Gravitational astronomy is now also becoming reality. We present an extensive review on the physical foundations and the emission mechanism of gravitational waves in detail, and discuss the possibility of their detections.

We report on the multi-epoch observations of H2O maser emission in star forming region OH 43.8-0.1 carried out with VLBI Exploration of Radio Astrometry (VERA). The large-scale maser distributions obtained by single-beam VLBI mapping reveal new maser spots scattered in area of 0.7 x 1.0 arcsec, in addition to a `shell-like' structure with a scale of 0.3 x 0.5 arcsec which was previously mapped by Downes et al.(1979). Proper motions are also obtained for 43 spots based on 5-epoch monitoring with a time span of 281 days. The distributions of proper motions show a systematic outflow in the north-south direction with an expansion velocity of ~8 km/s, and overall distributions of maser spots as well as proper motions are better represented by a bipolar flow plus a central maser cluster with a complex structure, rather than a shell with a uniform expansion such as those found in Cep A R5 and W75N VLA2. The distance to OH 43.8-0.1 is also estimated based on the statistical parallax, yielding D = 2.8 +/- 0.5 kpc. This distance is consistent with a near kinematic distance and rules out a far kinematics distance (~9 kpc), and the LSR velocity of OH 43.8-0.1 combined with the distance provides a constraint on the flatness of the galactic rotation curve, that there is no systematic difference in rotation speeds at the Sun and at the position of OH 43.8-0.1, which is located at the galacto-centric radius of ~6.3 kpc.

We investigate $r$-process nucleosynthesis during the magnetohydrodynamical (MHD) explosion of supernova in a massive star of 13 $M_{\odot}$. Contrary to the case of the spherical explosion, jet-like explosion due to the combined effects of the rotation and magnetic field lowers the electron fraction significantly inside the layers above the iron core. We find that the ejected material of low electron fraction responsible for the $r$-process comes out from the silicon rich layer of the presupernova model. This leads to the production up to the third peak in the solar $r$-process elements. We examine whether the fission affects the $r$-process paths by using the full nuclear reaction network with both the spontaneous and $\beta$-delayed fission included. Moreover, we pay particular attention how the mass formula affects the $r$-process peaks with use of two mass formulae. It is found that both formulae can reproduce the global abundance pattern up to the third peak though detailed distributions are rather different. We point out that there are variations in the $r$-process nucleosynthesis if the MHD effects play an important role in the supernova explosion.

We investigate cosmic-ray antiprotons emitted from the galactic primordial black holes in the Randall-Sundrum type-2 braneworld. The recent results of the BESS antiproton observation implies the existence of exotic primary sub-GeV antiprotons, one of whose most probable origin is PBHs in Our Galaxy. We show that the magnitude of antiproton flux from PBHs in the RS braneworld is proportional to negative power of the AdS radius, and immediately find that a large extra-dimension can relax upper-limits on the abundance of the Galactic PBHs. If actually there are more PBHs than the known upper-limit obtained in the pure 4D case, they set a lower bound on the size of the extra dimension above at least 10^{20} times 4D Planck-length to avoid inconsistency. On completion of the numerical studies, we show that these constraints on the AdS radius is comparable to those obtained from the diffuse photon background by some of the authors in the previous paper. Moreover, in the low accretion-rate case, only antiprotons can constrain the braneworld. We show that we will detect signatures of the braneworld as a difference between the flux of the antiprotons predicted in 4D and 5D by future observations in sub-GeV region with a few percent precision.

We study the impact of possible high density baryonic bubbles on the early formed QSO, IGM, and metal poor stars. Such bubbles could be created, under certain conditions, in Affleck-Dine model of baryogenesis and may occupy a relatively small fraction of space, while the dominant part of the cosmological volume has the normal observed baryon-to-photon ratio $\eta = 6\cdot 10^{-10}$. The value of $\eta$ in the bubbles, could be much larger than the usually accepted one (it might be even close to unity) without contradicting the existing data on light element abundances and the observed angular spectrum of CMBR. We find upper bounds on $\eta$ by comparing heavy elements' abundances produced in BBN and those of metal poor stars. We conclude that $\eta$ should be smaller than $10^{-5}$ in some metal poor star regions.

Calculations of equilibrium properties of dense matter predict that at subnuclear densities nuclei can be rodlike or slablike. To investigate whether transitions between phases with non-spherical nuclei can occur during the collapse of a star, we perform quantum molecular dynamic simulations of the compression of dense matter. We have succeeded in simulating the transitions between rodlike and slablike nuclei and between slablike nuclei and cylindrical bubbles. Our results strongly suggest that non-spherical nuclei can be formed in the inner cores of collapsing stars.

Present and future observations of supernova relic neutrinos (SRNs), i.e., a cosmological neutrino background from past core-collapse supernova explosions, potentially give us useful information concerning various fields of astrophysics, cosmology and particle physics. We review recent progress of theoretical and observational studies of SRNs, particularly focusing on the detectability and also on implications for cosmic star formation history and neutrino physics.

We perform a series of two-dimensional magnetohydrodynamic simulations of supernova cores. Since the distributions of the angular momentum and the magnetic fields of strongly magnetized stars are quite uncertain, we systematically change the combinations of the strength of the angular momentum, the rotations law, the degree of differential rotation, and the profiles of the magnetic fields to construct the initial conditions. By so doing, we estimate how the rotation-induced anisotropic neutrino heating are affected by the strong magnetic fields through parity-violating effects and first investigate how the north-south asymmetry of the neutrino heating in a strongly magnetized supernova core could be. As for the microphysics, we employ a realistic equation of state based on the relativistic mean field theory and take into account electron captures and the neutrino transport via the neutrino leakage scheme. With these computations, we find that the parity-violating corrections reduce $ \lesssim 0.5 %$ of the neutrino heating rate than that without the magnetic fields in the vicinity of the north pole of a star, on the other hand, enhance about $ \lesssim 0.5 %$ in the vicinity of the south pole. If the global asymmetry of the neutrino heating in the both of the poles develops in the later phases, the newly born neutron star might be kicked toward the north pole in the subsequent time.

We perform a series of two-dimensional magnetohydrodynamic simulations of the rotational core-collapse of a magnetized massive star. We employ a realistic equation of state and take into account the neutrino cooling by the so-called leakage scheme. In this study we systematically investigate how the strong magnetic field and the rapid rotation affect the propagation of the shock waves. Our results show that in the case of the strong initial poloidal magnetic field, the toroidal magnetic field amplified by the differential rotation, becomes strong enough to generate a tightly collimated shock wave along the rotational axis. On the other hand, in the case of the weak initial magnetic field, although the differential rotation amplifies toroidal magnetic field over the long rotational period, the launched shock wave is weak and the shape of it becomes wider. The former case is expected to be accompanied by the formation of the so-called magnetar. Our models with rapid rotation and strong magnetic field can create a nozzle formed by the collimated shock wave. This might be the analogous situation of the collapsar that is plausible for the central engine of the Gamma-Ray Bursts.

We report on a strong outburst of the W49N H2O maser observed with VERA. Single-dish monitoring with VERA 20 m telescopes detected a strong outburst of the maser feature at V_LSR = -30.7 km/s in 2003 October. The outburst had a duration of ~100 days and a peak intensity of 7.9 x 10^4 Jy, being one of the strongest outbursts in W49N observed so far. VLBI observations with the VERA array were also carried out near to the maximum phase of the outburst, and the outburst spot was identified in the VLBI map. While the map was in good agreement with previous studies, showing three major concentrations of maser spots, we found a newly formed arc-like structure in the central maser concentration, which may be a shock front powered by a forming star or a star cluster. The outburst spot was found to be located on the arc-like structure, indicating a possible connection of the present outburst to a shock phenomenon.

Recently, the High Resolution Fly's Eye (HiRes) experiment claims that there is no small scale anisotropy in the arrival distribution of ultra-high energy cosmic rays (UHECRs) above $E>10^{19}$ eV contrary to the Akeno Giant Air Shower Array (AGASA) observation. In this paper, we discuss the statistical significance of this discrepancy between the two experiments. We calculate arrival distribution of UHECRs above $10^{19}$ eV predicted by the source models constructed using the Optical Redshift Survey galaxy sample. We apply the new method developed by us for calculating arrival distribution in the presence of the galactic magnetic field. The great advantage of this method is that it enables us to calculate UHECR arrival distribution with lower energy ($\sim 10^{19}$ eV) than previous studies within reasonable time by following only the trajectories of UHECRs actually reaching the earth. It has been realized that the small scale anisotropy observed by the AGASA can be explained with the source number density $\sim 10^{-5 \sim -6}$ Mpc$^{-3}$ assuming weak extragalactic magnetic field ($B \le 1$ nG). We find that the predicted small scale anisotropy for this source number density is also consistent with the current HiRes data. We thus conclude that the statement by the HiRes experiment that they do not find small scale anisotropy in UHECR arrival distribution is not statistically significant at present. We also show future prospect of determining the source number density with increasing amount of observed data.

We perform a series of two-dimensional, axisymmetric, magnetohydrodynamic simulations of the rotational collapse of a supernova core. In order to calculate the waveforms of the gravitational wave, we derive the quadrupole formula including the contributions from the electromagnetic fields. Recent stellar evolution calculations imply that the magnetic fields of the toroidal components are much stronger than those of the poloidal ones at the presupernova stage. Thus, we systematically investigate the effects of the toroidal magnetic fields on the amplitudes and waveforms. Furthermore, we employ the two kinds of the realistic equation of states, which are often used in the supernova simulations. Then, we investigate the effects of the equation of states on the gravitational wave signals. With these computations, we find that the peak amplitudes are lowered by an order of 10% for the models with the strongest toroidal magnetic fields. However, the peak amplitudes are mostly within sensitivity range of laser interferometers such as TAMA and the first LIGO for a source at a distance of 10 kpc. Furthermore, we point out that the amplitudes of second peaks are still within the detection limit of the first LIGO for the source, although the characteristics of second peaks are reduced by the magnetic fields. We stress the importance of the detection, since it will give us information about the angular momentum distribution of massive evolved stars. When we compare the gravitational waves from the two realistic equation of states, significant differences are not found, except that the typical frequencies of the gravitational wave become slightly higher for the softer equation of state.

We investigate the spectral shape of a high-energy diffuse photon emitted by evaporating primordial black holes (PBHs) in the Randall-Sundrum type II (RS2) braneworld. In their braneworld scenario, the nature of small PBHs is drastically modified from the ordinary four-dimensional case for the following two reasons. (i) dropping Hawking temperature, which equivalently lengthens the lifetime of the individual PBH due to the change of space-time topology and (ii) the effective increase of the total amount of PBHs caused by accretion during the earliest part of the radiation-dominated epoch, the brane high-energy phase. From studies of the expected spectral shape and its dependence on braneworld parameters, we obtain two qualitatively distinctive possibilities of constraints on the braneworld PBHs from the observations of diffuse high-energy photon background. If the efficiency of accretion in the high-energy phase exceeds a critical value, the existence of the extra dimension gives a more stringent upper bound on the abundance of PBHs than the 4D case and a small length scale for the extra dimension is favored. On the contrary, in the case below the critical accretion efficiency, we find that the constraint on the PBH abundance can be relaxed by a few orders of magnitude in exchange for the existence of the large extra dimension; its size may be even bounded in the region above 10^{19} times 4D Planck length scale provided the rest mass energy density of the PBHs relative to energy density of radiation is actually larger than 10^{-27} (4D upper bound) at their formation time. The above analytical studies are also confirmed numerically, and an allowed region for braneworld parameters and PBH abundance is clearly obtained.

Structure of cold and hot dense matter at subnuclear densities is investigated by quantum molecular dynamics (QMD) simulations. Obtained phase diagrams show that the density of the phase boundaries between the different nuclear structures decreases with increasing temperature due to the thermal expansion of nuclear matter region.

Resonant spin-flavour (RSF) conversions of supernova neutrinos, which is induced by the interaction between the nonzero neutrino magnetic moment and supernova magnetic fields, are studied for both normal and inverted mass hierarchy. As the case for the pure matter-induced neutrino oscillation (Mikheyev--Smirnov--Wolfenstein (MSW) effect), we find that the RSF transitions are strongly dependent on the neutrino mass hierarchy as well as the value of $\theta_{13}$. Flavour conversions are solved numerically for various neutrino parameter sets, with presupernova profile calculated by Woosley and Weaver. In particular, it is very interesting that the RSF-induced $\nu_\rme\to\bar\nu_\rme$ transition occurs, if the following conditions are all satisfied: the value of $\mu_\nu B$ ($\mu_\nu$ is the neutrino magnetic moment, and $B$ is the magnetic field strength) is sufficiently strong, the neutrino mass hierarchy is inverted, and the value of $\theta_{13}$ is large enough to induce adiabatic MSW resonance. In this case, the strong peak due to original $\nu_\rme$ emitted from neutronization burst would exist in time profile of the neutrino events detected at the Super-Kamiokande detector. If this peak were observed in reality, it would provide fruitful information on the neutrino properties. On the other hand, characters of the neutrino spectra are also different between the neutrino models, but we find that there remains degeneracy among several models. Dependence on presupernova models is also discussed.

Structure of cold nuclear matter at subnuclear densities for the proton fraction $x=0.5$, 0.3 and 0.1 is investigated by quantum molecular dynamics (QMD) simulations. We demonstrate that the phases with slablike and rodlike nuclei, etc. can be formed dynamically from hot uniform nuclear matter without any assumptions on nuclear shape, and also systematically analyze the structure of cold matter using two-point correlation functions and Minkowski functionals. In our simulations, we also observe intermediate phases, which has complicated nuclear shapes. It is found out that these phases can be characterized as those with negative Euler characteristic. Our result implies the existence of these kinds of phases in addition to the simple ``pasta'' phases in neutron star crusts and supernova inner cores. In addition, we investigate the properties of the effective QMD interaction used in the present work to examine the validity of our results. The resultant energy per nucleon $\epsilon_{n}$ of the pure neutron matter, the proton chemical $\mu_{p}^{(0)}$ in pure neutron matter and the nuclear surface tension $E_{\rm surf}$ are generally reasonable in comparison with other nuclear interactions.

We present numerical simulations on the propagation of UHE protons with energies of $(10^{19.5}-10^{22})$ eV in extragalactic magnetic fields over 1 Gpc. We use the ORS galaxy sample, which allow us to accurately quantify the contribution of nearby sources to the energy spectrum and the arrival distribution, as a source model. We calculate three observable quantities, cosmic ray spectrum, harmonic amplitude, and two point correlation function from our data of numerical simulations. With these quantities, we compare the results of our numerical calculations with the observation. We show that the three observable quantities including the GZK cutoff of the energy spectrum can be reproduced in the case that the number fraction $\sim 10^{-1.7}$ of the ORS galaxies more luminous than -20.5 mag is selected as UHECR sources. In terms of the source number density, this constraint corresponds to $10^{-6}$ Mpc$^{-3}$. However, since mean number of sources within the GZK sphere is only $\sim 0.5$ in this case, the AGASA 8 events above $10^{20.0}$ eV, which do not constitute the clustered events with each other, can not be reproduced. On the other hand, if the cosmic ray flux measured by the HiRes, which is consistent with the GZK cutoff, is correct and observational features about the arrival distribution of UHECRs are same as the AGASA, our source model can explain both the arrival distribution and the flux at the same time. Thus, we conclude that large fraction of the AGASA 8 events above $10^{20}$ eV might originate in the topdown scenarios, or that the cosmic ray flux measured by the HiRes experiment might be better. We also discuss the origin of UHECRs below $10^{20.0}$ eV through comparisons between the number density of astrophysical source candidates and our result ($\sim 10^{-6}$ Mpc$^{-3}$).

We present a new method for calculating arrival distribution of UHECRs including modifications by the galactic magnetic field. We perform numerical simulations of UHE anti-protons, which are injected isotropically at the earth, in the Galaxy and record the directions of velocities at the earth and outside the Galaxy for all of the trajectories. We then select some of them so that the resultant mapping of the velocity directions outside the Galaxy of the selected trajectories corresponds to a given source location scenario, applying Liouville's theorem. We also consider energy loss processes of UHE protons in the intergalactic space. Applying this method to our source location scenario which is adopted in our recent study and can explain the AGASA observation above 4 \times 10^{19} eV, we calculate the arrival distribution of UHECRs including lower energy (E>10^19 eV) ones. We find that our source model can reproduce the large-scale isotropy and the small-scale anisotropy on UHECR arrival distribution above 10^19 eV observed by the AGASA. We also demonstrate the UHECR arrival distribution above 10^19 eV with the event number expected by future experiments in the next few years. The interesting feature of the resultant arrival distribution is the arrangement of the clustered events in the order of their energies, reflecting the directions of the galactic magnetic field. This is also pointed out by Alvarez-Muniz et al.(2002). This feature will allow us to obtain some kind of information about the composition of UHECRs and the magnetic field with increasing amount of data.

We have done a series of two-dimensional hydrodynamic simulations of the rotational collapse of a supernova core in axisymmetry. We have employed a realistic equation of state (EOS) and taken into account electron captures and neutrino transport by the so-called leakage scheme. It is an important progress to apply the realistic EOS coupled with the microphysics to 2-D simulations for computing gravitational radiation in rotational core collapse. We have used the quadrupole formula to calculate the amplitudes and the waveforms of gravitational wave assuming the Newtonian gravity. From these computations, we have extended the conventional category of the gravitational waveforms. Our results have shown that the peak amplitudes of gravitational wave are mostly within the sensitivity range of the laser interferometers such as TAMA and first LIGO for a source at a distance of 10 kpc. Furthermore we have found that the amplitudes of the second peaks are within the detection limit of first LIGO for the source and first pointed out the importance of the detections, since they will give us the information as to the angular momentum distribution of evolved massive stars.

We perform 2-dimensional relativistic hydrodynamical simulations in the context of collapsar model. Calculations of explosive nucleosynthesis are also accomplished. We investigate the influence of the structure of the progenitor and energy deposition rate on the resulting explosive nucleosynthesis, assuming that 56Ni is mainly synthesized in the jet launched by the neutrino heating. We show the amount of 56Ni is very sensitive to the energy deposition rate. Thus we conclude that it is quite natural not to detect an underlying supernova in some X-ray afterglows as in GRB 010921. We also point out the possibility that the relative abundance of the elements with intermediate mass number such as Si and S in the X-ray afterglow of GRB 011211 may be naturally explained if the energy deposition rate at the central engine is relatively long because little amount of 56Ni should be synthesized under such an environment. If this discussion is true, there should be correlation between the line features in the X-ray afterglow and duration of the GRB. It should be noted that the duration of GRB 011211 is 270 seconds, making it the longest burst ever observed by Beppo-SAX although it suffers from the effect of red-shift (z_host=2.14), which supports our conclusion. Our results also suggest that the type I collapsar model in which the energy deposition rate is relatively low (\dot{E} \sim 10^{51} erg/s) might have difficulty in reproducing the observed amount of 56Ni in a hypernova such as SN 1998bw. This means that the mechanism of the central engine of a hypernova accompanying GRB may be constrained by the discussion of explosive nucleosynthesis.

We study the resonant spin-flavor (RSF) conversion of supernova neutrinos, which is induced by the interaction between the nonzero neutrino magnetic moment and the supernova magnetic fields, and its dependence on presupernova models. As the presupernova models, we adopt the latest ones by Woosley, Heger, and Weaver, and, further, models with both solar and zero metallicity are investigated. Since the $(1-2Y_e)$ profile of the new presupernova models, which is responsible for the RSF conversion, suddenly drops at the resonance region, the completely adiabatic RSF conversion is not realized, even if $\mu_\nu B_0= (10^{-12}\mu_B)(10^{10} G)$, where $B_0$ is the strength of the magnetic field at the surface of the iron core. In particular for the model with zero metallicity, the conversion is highly nonadiabatic in the high energy region, reflecting the $(1-2Y_e)$ profile of the model. In calculating the flavor conversion, we find that the shock wave propagation, which changes density profiles drastically, is a much more severe problem than it is for the pure Mikheyev-Smirnov-Wolfenstein (MSW) conversion case. This is because the RSF effect occurs at a far deeper region than the MSW effect. To avoid the uncertainty concerning the shock propagation, we restrict our discussion to 0.5 s after the core bounce (and for more conservative discussion, 0.25 s), during which the shock wave is not expected to affect the RSF region. We also evaluate the energy spectrum at the Super-Kamiokande detector for various models using the calculated conversion probabilities, and find that it is very difficult to obtain useful information on the supernova metallicities and magnetic fields or on the neutrino magnetic moment from the supernova neutrino observation. Future prospects are also discussed.

We present the results of first dual-beam observations with VERA (VLBI Exploration of Radio Astrometry). The first dual-beam observations of a pair of H2O maser sources W49N and OH43.8-0.1 have been carried out on 2002 May 29 and July 23, and fringes of H2O maser lines at 22 GHz have been successfully detected. While the residual fringe phases of both sources show rapid variations over 360 degree due to the atmospheric fluctuation, the differential phase between the two sources remains constant for 1 hour with r.m.s. of 8 degree, demonstrating that the atmospheric phase fluctuation is removed effectively by the dual-beam phase referencing. The analysis based on Allan standard deviation reveals that the differential phase is mostly dominated by white phase noise, and the coherence function calculated from the differential phase shows that after phase referencing the fringe visibility can be integrated for arbitrarily long time. These results demonstrate VERA's high capability of phase referencing, indicating that VERA is a promising tool for phase referencing VLBI astrometry at 10 microarcsec level accuracy.

Recent deep Chandra LETG+HRC-S observations suggest the possibility that RX J1856.5--3754 is a compact star whose radiation radius is 3.8--8.2 km. In this paper, we systematically calculate the mass-radius relations of quark stars within the bag model. Assuming that RX J1856.5--3754 is a pure quark star, we derive an upper limit on its mass for various sets of the bag-model parameters.

We predict the arrival distribution of UHECRs above $4 \times 10^{19}$ eV with the event number expected by future experiments in the next few years. We perform event simulations with the source model which is adopted in our recent study and can explain the current AGASA observation. At first, we calculate the harmonic amplitude and the two point correlation function for the simulated event sets. We find that significant anisotropy on large angle scale will be observed when $\sim 10^3$ cosmic rays above $4 \times 10^{19}$ eV are detected by future experiments. The statistics of the two point correlation function will also increase. The angle scale at which the events have strong correlation with each other corresponds to deflection angle of UHECR in propagating in the EGMF, which in turn can be determined by the future observations. We further investigate the relation between the number of events clustered at a direction and the distance of their sources. Despite the limited amount of data, we find that the C2 triplet events observed by the AGASA may originate from the source within 100 Mpc. Merger galaxy Arp 299 (NGC 3690 + IC 694) is the best candidate for their source. If data accumulate, the UHECR sources within $\sim 100$ Mpc can be identified from observed event clusterings significantly. This will provide some kinds of information about poorly known parameters which influence the propagation of UHECRs, such as extragalactic and galactic magnetic field, chemical composition of observed cosmic rays. Also, we will reveal their origin with our method to identify the sources of UHECR. Finally, we predict the arrival distribution of UHECRs above $10^{20}$ eV, which is expected to be observed if the current HiRes spectrum is correct, and discuss their statistical features and implications.

We investigate the flux of supernova relic neutrinos (SRN) for several neutrino oscillation models with parameters inferred from recent experimental results. In the calculation, we adopt the realistic {\it time-dependent} supernova density profile, which is very different from the static progenitor profile owing to shock propagation. The Earth matter effect is also included appropriately using realistic density profile of the Earth. As a result, these two effects are found to induce the flux difference by a few % in the detection energy range ($E_\nu > 19.3$ MeV). We also set 90% C.L. upper limit on SRN flux for various oscillation models using the recently released observational result by Super-Kamiokande (SK). In the near future, further reduced upper limit is actually expected, when the current SK data are reanalyzed using some technique to reduce background events against SRN signals. It is expected that the reduced upper limit is sufficient to provide useful implications for neutrino oscillation.

The effect of the shock propagation on neutrino oscillation in supernova is studied paying attention to evolution of average energy of $\nu_{e}$ and $\bar{\nu}_{e}$. We show that the effect appears as a decrease in average $\nu_{e}$ (in case of inverted mass hierarchy, $\bar{\nu}_{e}$) energy at stellar surface as the shock propagates. It is found that the effect is significant 2 seconds after bounce if $3 \times 10^{-5} < \sin^{2}{\theta_{13}} < 10^{-2}$.

If Ultra-High Energy Cosmic Rays (UHECRs) are originated from nearby galaxies, modeling of the distribution of nearby galaxies is important to an accurate estimate the source number density of UHECRs. We investigate uncertainty of the selection function of the Optical Redshift Survey (ORS), which we used to construct a source model of UHECRs. The investigation is based on a comparison of numbe counts of ORS galaxies with those of the spectroscopic sample of the Sloan Digital Sky Survey (SDSS) Early Data Release (EDR). We carefully count galaxies in the same absolute magnitude bin from the two samples. We find a slight systematic overestimate of the ORS counts outside 5000 km s$^{-1}$ by about a factor of 2. We revise the selection function of the ORS assuming that the SDSS counts are correct. Our revision is based on the absorption given in the ORS catalog as well as that computed from Schlegel et al. (1998), which is systematically larger than the former by $A_B \sim 0.1$ mag in the region of low absorption. It is found that introduction of Schlegel et al.'s absorption changes one of the parameters of the ORS selection function by more than 10%. The revision should be taken into account in the future analysis of the source number density of UHECRs based on the ORS. Using the revised selection function, we determine the global structure of the Local Supercluster (LSC) with a source model of UHECRs, that is, a number-density model consisting of a uniform spherical halo and an exponential disk. We find that the revision is insignificant in terms of the structure of the LSC. However, the revised selection function will be useful to other studies such as peculiar velocity and correlation function.

We investigate resonant spin-flavor (RSF) conversions of supernova neutrinos which are induced by the interaction of neutrino magnetic moment and supernova magnetic fields. From the formulation which includes all three-flavor neutrinos and anti-neutrinos, we give a new crossing diagram that includes not only ordinary MSW resonance but also magnetically-induced RSF effect. With the diagram, it is found that four conversions occur in supernova, two are induced by the RSF effect and two by the pure MSW. We also numerically calculate neutrino conversions in supernova matter, using neutrino mixing parameters inferred from recent experimental results and a realistic supernova progenitor model. The results indicate that until 0.5 seconds after core bounce, the RSF-induced $\bar{\nu}_e \leftrightarrow \nu_\tau$ transition occurs efficiently (adiabatic resonance), when $\mu_\nu \agt 10^{-12} \mu_B (B_0 / 5 \times 10^{9} \mathrm G)^{-1}$, where $B_0$ is the strength of the magnetic field at the surface of iron core. We also evaluate the energy spectrum as a function of $\mu_\nu B_0$ at the SuperKamiokande detector and the Sudbury Neutrino Observatory using the calculated conversion probabilities, and find that the spectral deformation might have possibility to provide useful information on neutrino magnetic moment as well as magnetic field strength in supernovae.

Structure of cold dense matter at subnuclear densities is investigated by quantum molecular dynamics (QMD) simulations. We succeeded in showing that the phases with slab-like and rod-like nuclei etc. can be formed dynamically from hot uniform nuclear matter without any assumptions on nuclear shape. We also observe intermediate phases, which has complicated nuclear shapes. Geometrical structures of matter are analyzed with Minkowski functionals, and it is found out that intermediate phases can be characterized as ones with negative Euler characteristic. Our result suggests the existence of these kinds of phases in addition to the simple ``pasta'' phases in neutron star crusts.

We investigate the flux and the event rate of the supernova relic neutrino background (SRN) at the SuperKamiokande detector for various neutrino oscillation models with parameters inferred from recent experimental results. A realistic model of neutrino emission from supernova explosions and several models of the cosmic star formation history are adopted in the calculation. The number flux over entire energy range is found to be $11-15 \mathrm{cm^{-2}s^{-1}}$. We discuss the detection possibility of SRN at SuperKamiokande, comparing this SRN flux with other background neutrinos in more detail than previous studies. Even though there is no energy window in which SRN is dominant, we might detect it as the distortion of the other background event. We found in the energy range $17-25 \mathrm{MeV}$ the expected event rate at SuperKamiokande $0.4-0.8 ~\mathrm{yr^{-1}}$. In this range, ten-year observation might enable us to detect SRN signal (at one sigma level) in the case of LMA solar neutrino solution. We also investigate event rate at SNO and KamLAND. Although we can find energy window, the expected event rate is rather small (0.03 yr$^{-1}$ for SNO, 0.1 yr$^{-1}$ for KamLAND).

We estimate the flux of the GRB neutrino background and compute the event rate at SK and TITAND in the collapsar model, assuming that GRB formation rate is proportional to the star formation rate. We find that the flux and the event rate depend sensitively on the mass-accretion rate. Although the detection of signals from GRBs seems to be difficult by SK, we find that we can detect them by TITAND (optimistically about 20 events per year in the energy range where the energy of anti-electron neutrino is larger than 50 MeV). Therefore, we show that we will obtain the informations on the mass-accretion rate of collapsar by observing the GRB neutrino background, which, in turn, should give us informations on the total explosion energy of GRBs.

We perform a detailed study of the Earth matter effects on supernova neutrinos with neutrino oscillation parameter LMA and small $\theta_{13}$. The Earth effects show significant dependences on the neutrino path length inside the Earth and the value of $\Delta m^{2}_{12}$. We investigate rather optimistically a possibility that we can probe the value of $\Delta m^{2}_{12}$ by the Earth effects. We assume that $\theta_{12}$ and the direction of the supernova are known with enough accuracy and that the resonance that occurs at higher density in supernova envelope is completely nonadiabatic. Further the neutrino spectra before neutrinos go through the Earth are assumed to be known. Then we show that making use of these dependences, we can obtain implication for the value of $\Delta m^{2}_{12}$ by comparing the observed energy spectrum to the predicted one. When SK detects neutrinos from supernova at 10kpc which traveled through the Earth (nadir angle $<$ 80 degree), $\Delta m^{2}_{12}$ can be determined with an accuracy of $\sim 10%$. In much of the neutrino-detection-time-$\Delta m^{2}_{12}$ plane, $\Delta m^{2}_{12}$ might be determined with an accuracy equal to or better than $\pm 0.5 \times 10^{-5} {\rm eV}^{2}$.

Structure of cold dense stellar matter at subnuclear densities is investigated with quantum molecular dynamics (QMD). We succeeded in showing that the phases with slab-like and rod-like nuclei etc. can be formed dynamically without any assumptions on the nuclear shape. Our result suggests the existence of these kinds of phases in neutron star crusts.

Formation and decay of rod-like and slab-like nuclei, as may be encountered in the outer parts of neutron stars, are examined within a zero-temperature stability analysis with respect to perturbations inducing fission and proton clustering. The nonspherical nuclei are found to be stable against these perturbations, implying the possibility that the size of the stellar region containing such nuclei exceeds the equilibrium prediction.

The effects of three-flavor neutrino oscillation on the supernova neutrino spectrum are studied. We calculate the expected event rate and energy spectra, and their time evolution at the Superkamiokande (SK) and the Sudbury Neutrino Observatory (SNO), by using a realistic neutrino burst model based on numerical simulations of supernova explosions. We also employ a realistic density profile based on a presupernova model for the calculation of neutrino conversion probability in supernova envelopes. These realistic models and numerical calculations allow us to quantitatively estimate the effects of neutrino oscillation in a more realistic way than previous studies. We then found that the degeneracy of the solutions of the solar neutrino problem can be broken by the combination of the SK and SNO detections of a future Galactic supernova.

We have performed the detailed numerical simulations on the propagation of the UHE protons in the energy range $E=(10^{19.5} - 10^{22.0}$) eV in the relatively strong extra-galactic magnetic field with strength $B= (10, 100)$ nG within about 40 Mpc. In this case, the deflection angles of UHECRs become so large that the no counterparts problem is simply solved. As for the source distribution, we assumed that it is proportional to the number distribution of galaxies within the GZK sphere. We have found many clusters, which mean the small-scale anisotropy, in our simulations. It has been also shown that the observed energy spectrum is well reproduced in our models without any fine-tuned parameter. We have used the correlation value in order to investigate statistically the similarity between the distribution of arrival directions of UHECRs and that of galaxies. We have found that each correlation value for each parameter set begins to converge when the number of the detected events becomes ${\large O}(10^3)$. Since the expected number counts by the experiment of the next generation such as TA, HiRes, Auger, and EUSO are thought to be the order of $10^3$, we will be able to determine the source distribution and values of the parameters in this study in the very near future. Compared with the AGASA data, the significant anisotropy on the arrival directions of UHECRs are found in the analysis of first and second harmonics. This may originate from the incompleteness of the ORS database. This problem may be also solved if the source distribution is slightly changed. For example, this problem may be solved if we assume that UHECRs come from some of the galaxies such as AGNs and radio galaxies.

The Earth effects on the energy spectra of supernova neutrinos are studied. We analyze numerically the time-integrated energy spectra of neutrino in a mantle-core-mantle step function model of the Earth's matter density profile. We consider a realistic frame-work in which there are three active neutrinos whose mass squared differences and mixings are constrained by the present understanding of solar and atmospheric neutrinos. We find that the energy spectra change for some allowed mixing parameters. Especially, the expected number of events at SNO shows characteristic behavior with respect to energy, i.e. a great dip and peak. We show that observations of the Earth effect allow us to identify the solar neutrino solution and to probe the mixing angle $\theta_{2}$.

From the observation of a millisecond pulsar in SN 1987A, the following implications are obtained. 1) The pulsar spindown in SN 1987A is caused by radiating gravitational waves rather than by magnetic dipole radiation and/or relativistic pulsar winds. 2) A mildly deformed shock wave would be formed at the core-collapse and explosion in SN 1987A, which is consistent with the conclusion given in Nagataki (2000). 3) The gravitational waves from the pulsar should be detected in several years using a Fabry-Perot-Michelson interferometer as the gravitational detector, such as LIGO and TAMA. 4) The neutrino oscillation model is not a promising model for the explanation of the kick velocity of the pulsar in SN 1987A. The hydrodynamical instability model is more favored.

Geometrical structure of matter at subnuclear densities is investigated in the presence of a degenerate gas of neutrinos as encountered in stellar collapse. The crystalline phases with spherical, cylindrical and planar nuclei as well as with spherical and cylindrical nuclear bubbles are considered by using a compressible liquid-drop model. This model allows for uncertainties in the lepton fraction $Y_{\rm L}$ in addition to those in the nuclear surface tension $E_{\rm surf}$ and in the proton chemical potential in bulk neutron matter $\mu_{\rm p}^{(0)}$. The phase diagrams obtained at zero temperature show that only the phases with rod-like and slab-like nuclei occur at typical values of $Y_{\rm L}$, $E_{\rm surf}$ and $\mu_{\rm p}^{(0)}$, whereas the bubble phases, especially with spherical bubbles, are at best expected at hypothetically low values of $Y_{\rm L}$ and/or $E_{\rm surf}$. For the rod-like and slab-like nuclei, thermally induced displacements are calculated from their respective elastic constants. It is found that at temperatures appropriate to supernova cores, thermal fluctuations would destroy the layered lattice of slab-like nuclei almost independently of the nuclear models and of the degree of the neutrino degeneracy.

In the framework of the Hartle-Hawking no-boundary proposal, we investigated quantum creation of the multidimensional universe with a cosmological constant ($\Lambda$) but without matter fields. We have found that the classical solutions of the Euclidean Einstein equations in this model have ``quasi-attractors'', i.e., most trajectories on the a-b plane, where a and b are the scale factors of external and internal spaces, go around a point. It is presumed that the wave function of the universe has a hump near this quasi-attractor point. In the case that both the curvatures of external and internal spaces are positive, and $\Lambda>0$, there exist Lorentzian solutions which start near the quasi-attractor, the internal space remains microscopic, and the external space evolves into our macroscopic universe.

Equilibrium phase diagrams for neutron star matter at subnuclear densities are obtained at zero temperature. Spherical, rod-like and slab-like nuclei as well as spherical and rod-like nuclear bubbles are taken into account by using a compressible liquid-drop model. This model is designed to incorporate uncertainties in the nuclear surface tension and in the proton chemical potential in a gas of dripped neutrons. The resultant phase diagrams show that for typical values of these quantities, the phases with rod-like nuclei and with slab-like nuclei occur in the form of Coulomb lattice at densities below a point where the system becomes uniform. Thermal fluctuations leading to displacements of such nuclei from their equilibrium positions are considered through explicit evaluations of their elastic constants; these fluctuations can be effective at destroying the layered lattice of slab-like nuclei in the temperature region typical of matter in the neutron star crust.

The influence of the presence of hyperons in dense hadronic matter on the quantum nucleation of quark matter is examined at low temperatures relevant to neutron star cores. We calculate the equation of state and the composition of matter before and after deconfinement by using a relativistic mean-field theory and an MIT bag model, respectively; the case in which hyperons are present in the hadronic system is considered, together with the case of the system without hyperons. We find that strangeness contained in hyperons acts to reduce a density jump at deconfinement as well as a lepton fraction in the hadronic phase. As a result of these reductions, a quark matter droplet being in a virtual or real state has its effective mass lightened and its electric charge diminished into nearly zero. The Coulomb screening of leptons on the droplet charge, which has significance to the droplet growth after nucleation in the absence of hyperons, is thus shown to be of little consequence. If the effective droplet mass is small enough to become comparable to the height of the potential barrier, the effect of relativity brings about an exponential increase in the rate of droplet formation via quantum tunneling, whereas the role played by energy dissipation in decelerating the droplet formation, dominant for matter without hyperons, becomes of less importance. For matter with and without hyperons, we estimate the overpressure needed to form the first droplet in the star during the compression due to stellar spin-down or mass accretion from a companion star. The temperature at which a crossover from the quantum nucleation to the Arrhenius-type thermal nucleation takes place is shown to be large compared with the temperature of matter in the core.

We have estimated the dispersion of the inferred relative frequencies of Type-Ia and Type- II supernovae ($N_{Ia}/ N_{II}$) in our Galaxy by fitting the numerical results of supernova nucleosynthesis to the solar-system abundances. The ratio $N_{Ia}/ N_{II}$ is estimated to be 0.056--0.14, which is consistent with the observation, if the model of Woosley and Weaver (1995, WW95) is adopted for Type-II supernovae (SNe II). On the other hand, the upper limit of $N_{Ia}/ N_{II}$ becomes too large for the model of Hashimoto (1995, Ha95). However, Ha95 can fit the solar values better than WW95 as far as the abundant nuclei are concerned. These results mean that Ha95 can reproduce the main nuclei of the solar-system well and that WW95 can reproduce the solar values over a wide mass-number range. We also note that $N_{Ia}/ N_{II}$ tends to become smaller if the delayed detonation model is adopted for Type-Ia supernovae (SNe Ia). The dispersion of $N_{Ia}/ N_{II}$ obtained in the present investigation is larger than that concluded by Tsujimoto et al. (1995), which will have an influence on the estimate of the average life time of SNe Ia's progenitors and/or the star-formation rate history in our galaxy.

The large abundance ratio of $^{44}Ti/^{56}Ni$ in Cas A is puzzling. In fact, the ratio seems to be larger than the theoretical constraint derived by Woosley & Hoffman (1991). However, this constraint is obtained on the assumption that the explosion is spherically symmetric, whereas Cas A is famous for the asymmetric form of the remnant. Recently, Nagataki et al. (1997) calculated the explosive nucleosynthesis of axisymmetrically deformed collapse-driven supernova. They reported that the ratio of $^{44}Ti/^{56}Ni$ was enhanced by the stronger alpha-rich freezeout in the polar region. In this paper, we apply these results to Cas A and examine whether this effect can explain the large amount of $^{44}Ti$ and the large ratio of $^{44}Ti/^{56}Ni$. We demonstrate that the conventional spherically symmetric explosion model can not explain the $^{44}$Ti mass produced in Cas A if its lifetime is shorter than $\sim$ 80 years and the intervening space is transparent to the gamma-ray line from the decay of $^{44}$Ti. On the other hand, we show the axisymmetric explosion models can solve the problem. We expect the same effect from a three dimensionally asymmetric explosion, since the stronger alpha-rich freezeout will also occur in that case in the region where the larger energy is deposited.

The origin of the chemical composition of the intracluster medium (ICM) is discussed in this paper. In particular, the contribution from Type Ia supernovae (SNe Ia) to the ICM enrichment is shown to exist by adopting the fitting formulas which have been used in the analysis of the solar system abundances. Our analysis means that we can use the frequency of SNe Ia relative to SNe II as the better measure than $M_{Fe, SN Ia}/M_{Fe, total}$ for estimating the contribution of SNe Ia. Moreover, the chemical compositions of ICMs are shown to be similar to that of the solar system abundances. We can also reproduce the sulfur/iron abundance ratio within a factor of 2, which means that the abundance problem of sulfur needs not to be emphasized too strongly. We need more precise observations to conclude whether ICMs really suffer the shortage problem of sulfur or not.

Comparing Big Bang Nucleosynthesis predictions with the light element abundances inferred from observational data, we can obtain the strong constraints on some neutrino properties, e.g. number of neutrino species, mass, lifetime. Recently the deuterium abundances were measured in high red-shift QSO absorption systems. It is expected that they are close to the primordial values, however, two groups have reported inconsistent values which are different in one order of magnitude. In this paper we show how we can constrain on $\tau$ neutrino mass and its lifetime in each case when we adopt either high or low deuterium data. We find that if $0.01 \sec \lesssim \tau_{\nutau} \lesssim 1 \sec$ and $10\mev \lesssim m_{\nutau} \lesssim 24\mev$, the theoretical predictions agree with the low D/H abundances. On the other hand if we adopt the high D/H abundances, we obtain the upper bound of $\tau$ neutrino mass, $m_{\nutau}\lesssim 20 \mev$.

Future detection of a supernova neutrino burst by large underground detectors would give important information for the explosion mechanism of collapse-driven supernovae. We studied the statistical analysis for the future detection of a nearby supernova by using a numerical supernova model and realistic Monte-Carlo simulations of detection by the Super-Kamiokande detector. We mainly discuss the detectability of the signatures of the delayed explosion mechanism in the time evolution of the $\anue$ luminosity and spectrum. For a supernova at 10 kpc away from the Earth, we find that not only the signature is clearly discernible, but also the deviation of energy spectrum from the Fermi-Dirac (FD) distribution can be observed. The deviation from the FD distribution would, if observed, provide a test for the standard picture of neutrino emission from collapse-driven supernovae. For the $D$ = 50 kpc case, the signature of the delayed explosion is still observable, but statistical fluctuation is too large to detect the deviation from the FD distribution. We also propose a method for statistical reconstruction of the time evolution of $\anue$ luminosity and spectrum from data, by which we can get a smoother time evolution and smaller statistical errors than a simple, time-binning analysis. This method is useful especially when the available number of events is relatively small, e.g., a supernova in the LMC or SMC. Neutronization burst of $\nu_e$'s produces about 5 scattering events when $D$ = 10 kpc and this signal is difficult to distinguish from $\anue p$ events.

Explosive nucleosynthesis under the axisymmetric explosion in Type II supernova has been performed by means of two dimensional hydrodynamical calculations. We have compared the results with the observations of SN 1987A. Our chief findings are as follows: (1) $^{44}Ti$ is synthesized so much as to explain the tail of the bolometric light curve of SN 1987A. We think this is because the alpha-rich freezeout takes place more actively under the axisymmetric explosion. (2) $^{57}Ni$ and $^{58}Ni$ tend to be overproduced compared with the observations. However, this tendency relies strongly on the progenitor's model. We have also compared the abundance of each element in the mass number range $A= 16-73$ with the solar values. We have found three outstanding features. (1) For the nuclei in the range $A=16-40$, their abundances are insensitive to the initial form of the shock wave. This insensitivity is favored since the spherical calculations thus far can explain the solar system abundances in this mass range. (2) There is an enhancement around A=45 in the axisymmetric explosion compared with the spherical explosion fairly well. In particular, $^{44}Ca$, which is underproduced in the present spherical calculations, is enhanced significantly. (3) In addition, there is an enhancement around A=65. This tendency does not rely on the form of the mass cut but of the initial shock wave. This enhancement may be the problem of the overproduction in this mass range, although this effect would be relatively small since Type I supernovae are chiefly responsible for this mass number range.

Growth of Rayleigh-Taylor (R-T) instabilities under the axisymmetric explosion are investigated by two-dimensional hydrodynamical calculations. The degree of the axisymmetric explosion and amplitude of the initial perturbation are varied parametrically to find the most favorable parameter for reproducing the observed line profile of heavy elements. It is found that spherical explosion can not produce $^{56}Ni$ travelling at high velocity ($\sim 3000$km/sec), the presence of which is affirmed by the observation, even if the amplitude of initial perturbation is as large as 30%. On the other hand, strong axisymmetric explosion model produce high velocity $^{56}Ni$ too much. Weak axisymmetric explosion are favored for the reproduction of the observed line profile. We believe this result shows upper limit of the degree of the axisymmetric explosion. This fact will be important for the simulation of the collapse-driven supernova including rotation, magnetic field, and axisymmetric neutrino radiation, which have a possibility to cause axisymmetric supernova explosion. In addition, the origin of such a large perturbation does not seem to be the structure of the progenitor but the dynamics of the core collapse explosion itself since small perturbation can not produce the high velocity element even if the axisymmetric explosion models are adopted.

Recently it is reported that there is the discrepancy between big bang nucleosynthesis theory and observations (BBN crisis). We show that BBN predictions agree with the primordial abundances of light elements, He4, D, He3 and Li7 inferred from the observational data if an electron neutrino has a net chemical potential xi_{nu_e} due to lepton asymmetry. We estimate that xi_{nu_e} = 0.043^{+0.040}_{-0.040} (95% C.L.) and Omega_bh^2 = 0.015^{+0.006}_{-0.003} (95% C.L.).

Rates for nucleation of two-flavor quark matter in a neutron star core, originally composed of nuclear matter in $\beta$ equilibrium, are calculated at zero temperature by a quantum tunneling analysis incorporating the electrostatic energy. We find that a nucleated droplet would develop into bulk matter due to electron screening effects.

We show that evolution of the luminosity density of galaxies in the universe provides a powerful test for the geometry of the universe. Using reasonable galaxy evolution models of population synthesis which reproduce the colors of local galaxies of various morphological types, we have calculated the luminosity density of galaxies as a function of redshift $z$. Comparison of the result with recent measurements by the Canada-France Redshift Survey in three wavebands of 2800{\AA}, 4400{\AA}, and 1 micron at z<1 indicates that the \Lambda-dominated flat universe with \lambda_0 \sim 0.8 is favored, and the lower limit on \lambda_0 yields 0.37 (99% C.L.) or 0.53 (95% C.L.) if \Omega_0+\lambda_0=1. The Einstein-de Sitter universe with (\Omega_0, \lambda_0)=(1, 0) and the low-density open universe with (0.2, 0) are however ruled out with 99.86% C.L. and 98.6% C.L., respectively. The confidence levels quoted apply unless the standard assumptions on galaxy evolution are drastically violated. We have also calculated a global star formation rate in the universe to be compared with the observed rate beyond z \sim 2. We find from this comparison that spiral galaxies are formed from material accretion over an extended period of a few Gyrs, while elliptical galaxies are formed from initial star burst at z >~ 5 supplying enough amount of metals and ionizing photons in the intergalactic medium.

Transitions of nuclear compositions in the crust of a neutron star induced by stellar spin-down are evaluated at zero temperature. We construct a compressible liquid-drop model for the energy of nuclei immersed in a neutron gas, including pairing and shell correction terms, in reference to the known properties of the ground state of matter above neutron drip density, $4.3 \times 10^{11} g cm^{-3}$. Recent experimental values and extrapolations of nuclear masses are used for a description of matter at densities below neutron drip. Changes in the pressure of matter in the crust due to the stellar spin-down are calculated by taking into account the structure of the crust of a slowly and uniformly rotating relativistic neutron star. If the initial rotation period is $\sim$ ms, these changes cause nuclei, initially being in the ground-state matter above a mass density of about $3 \times 10^{13} g cm^{-3}$, to absorb neutrons in the equatorial region where the matter undergoes compression, and to emit them in the vicinity of the rotation axis where the matter undergoes decompression. Heat generation by these processes is found to have significant effects on the thermal evolution of old neutron stars with low magnetic fields; the surface emission predicted from this heating is compared with the $ROSAT$ observations of X-ray emission from millisecond pulsars and is shown to be insufficient to explain the observed X-ray luminosities.

The neutrino spin-flavor conversion of \bar\nu_e and \nu_\mu which is induced by the interaction of the Majorana neutrino magnetic moment and magnetic fields in the collapse-driven supernova is investigated in detail. We calculate the conversion probability by using the latest precollapse models of Woosley and Weaver (1995), and also those of Nomono and Hashimoto (1988), changing the stellar mass and metallicity in order to estimate the effect of the astrophysical uncertainties. Contour maps of the conversion probability are given for all the models as a function of neutrino mass squared difference and the neutrino magnetic moment times magnetic fields. It is shown that in the solar metallicity models some observational effects are expected with \Delta m^2 = 10^{-5}--10^{-1} [eV^2] and \mu_\nu >~ 10^{-12} (10^9 G / B_0) [\mu_B], where B_0 is the strength of the magnetic fields at the surface of the iron core. We also find that although the dependence on the stellar models or stellar mass is not so large, the metallicity of precollapse stars has considerable effects on this conversion. Such effects may be seen in a supernova in the Large or Small Magellanic Clouds, and should be taken into account when one considers an upper bound on \mu_\nu from the SN1987A data.

We investigate to what extent the oscillation or conversion of neutrinos enhances the expected event rate of the supernova relic neutrino background (SRN) at the Super-Kamiokande detector (SK). The SRN \bar\nu_e's can be almost completely exchanged with \nu_\mu-like neutrinos by the MSW oscillation under the inverse mass hierarchy with \Delta m^2 ~ 10^{-8}--10^5 [eV^2], or by the magnetic moment of Majorana neutrinos with \mu_\nu >~ 10^{-12} \mu_B and \Delta m^2 ~ 10^{-4}--10^0 [eV^2]. In the standard calculation of the SRN flux, the event rate of the SRN \bar\nu_e's at the SK in the observable energy range of 15--40 MeV can be enhanced from 1.2 yr^{-1} to 2.4 yr^{-1} if all \bar\nu_e's are exchanged with \nu_\mu-like neutrinos. The enhancement is prominent especially in the high energy range (\gtilde 25 MeV). In the astrophysically optimistic calculation, the event rate becomes as high as 9.4 yr^{-1}. Because the theoretical upper bound of the SRN events without oscillation is about 5 yr^{-1} taking account of the various astrophysical uncertainties, we might have to resort to the neutrino oscillation if more than 5 events in a year, as well as a significantly harder spectrum, were observed in the SK.

The spectrum of the supernova relic neutrino background (SRN) from collapse-driven supernovae ever occurred in the universe is calculated by using a realistic, time-dependent supernova rate derived from a standard model of galaxy evolution based on the population synthesis method. The SRN spectrum we show here is the most realistic at present, because the largest uncertainty in previous theoretical predictions has come from unrealistic assumptions of the supernova rate so far made. The SRN is one of the targets of the Superkamiokande (SK) detector which will be constructed in a year and the SRN, if at all detected, would provide a new tool to probe the history of supernova explosions in the universe. The expected event rate at the SK is therefore calculated in this paper. Our major results include: (1) the supernova rate is much higher in the early phase of evolution of galaxies and there appears a hump in the SRN spectrum in the low-energy region of $\ltilde 5$ MeV, (2) the SRN flux depends on the Hubble constant ($H_0$) in a way approximately proportional to $H_0^2$ and only weakly on the density parameter of the universe ($\Omega_0$) and a cosmological constant ($\lambda_0$), and (3) the plausible event rate at the SK is 1.2 yr$^{-1}$ in the observable energy range. Such a low event rate is due mainly to a quite low supernova rate at present which is averaged over the morphological types of galaxies. The most optimistic rate in our model is found to be 4.7 yr$^{-1}$, and if more events are detected, we will have to reconsider our current understanding of collapse-driven supernovae and evolution of galaxies.

It is greatly expected that the relic neutrino background from past supernovae is detected by Superkamiokande (SK) which is now under construction. We calculate the spectrum and the event rate at SK systematically by using the results of simulations of a supernova explosion and reasonable supernova rates. We also investigate the effect of a cosmological constant, $\Lambda$, on the spectrum, since some recent cosmological observations strongly suggest the existence of $\Lambda$. We find following results. 1) The spectrum has a peak at about 3 MeV, which is much lower than that of previous estimate ($6 \sim$ 10 MeV). 2) The event rate at SK in the range from 10 MeV to 50 MeV, where the relic neutrinos from past supernovae is dominant, is about $25 \, {h_{50}}^2 \, \left( \frac{R_{SN}}{0.1 {\rm yr^{-1}}} \right) \left(\frac{n_G \, h_{50}^{-3}}{0.02{\rm Mpc}^{-3}}\right)$ events per year, where $R_{SN}$ is the supernova rate in a galaxy, $n_G$ is the number density of galaxies, and $h_{50} = H_0 /$ (50km/s/Mpc), where $H_0$ is the Hubble constant. 3) The event rate is almost insensitive to $\Lambda$. The flux increases in the low energy side ($< 10$ MeV) with increasing $\Lambda$, but decreases in the high energy side (10 MeV $<$) in models in which the integrated number of supernovae in one galaxy is fixed.