We report the radiation hardness of a p-channel CCD developed for the X-ray CCD camera onboard the XRISM satellite. This CCD has basically the same characteristics as the one used in the previous Hitomi satellite, but newly employs a notch structure of potential for signal charges by increasing the implant concentration in the channel. The new device was exposed up to approximately $7.9 \times 10^{10} \mathrm{~protons~cm^{-2}}$ at 100 MeV. The charge transfer inefficiency was estimated as a function of proton fluence with an ${}^{55} \mathrm{Fe}$ source. A device without the notch structure was also examined for comparison. The result shows that the notch device has a significantly higher radiation hardness than those without the notch structure including the device adopted for Hitomi. This proves that the new CCD is radiation tolerant for space applications with a sufficient margin.

We present a new analysis of the interstellar protons toward the TeV $\gamma$-ray SNR RX J0852.0$-$4622 (G266.2$-$1.2, Vela Jr.). We used the NANTEN2 $^{12}$CO($J$ = 1-0) and ATCA & Parkes HI datasets in order to derive the molecular and atomic gas associated with the TeV $\gamma$-ray shell of the SNR. We find that atomic gas over a velocity range from $V_\mathrm{LSR}$ = $-4$ km s$^{-1}$ to 50 km s$^{-1}$ or 60 km s$^{-1}$ is associated with the entire SNR, while molecular gas is associated with a limited portion of the SNR. The large velocity dispersion of the HI is ascribed to the expanding motion of a few HI shells overlapping toward the SNR but is not due to the Galactic rotation. The total masses of the associated HI and molecular gases are estimated to be $\sim2.5 \times 10^4 $ $M_{\odot}$ and $\sim10^3$ $M_{\odot}$, respectively. A comparison with the H.E.S.S. TeV $\gamma$-rays indicates that the interstellar protons have an average density around 100 cm$^{-3}$ and shows a good spatial correspondence with the TeV $\gamma$-rays. The total cosmic ray proton energy is estimated to be $\sim10^{48}$ erg for the hadronic $\gamma$-ray production, which may still be an underestimate by a factor of a few due to a small filling factor of the SNR volume by the interstellar protons. This result presents a third case, after RX J1713.7$-$3946 and HESS J1731$-$347, of the good spatial correspondence between the TeV $\gamma$-rays and the interstellar protons, lending further support for a hadronic component in the $\gamma$-rays from young TeV $\gamma$-ray SNRs.

We review a non-standard Big-Bang nucleosynthesis (BBN) scenario within the minimal supersymmetric standard model, and propose an idea to solve both ${}^{7}$Li and ${}^{6}$Li problems. Each problem is a discrepancy between the predicted abundance in the standard BBN and observed one. We focus on the stau, a supersymmetric partner of tau lepton, which is a long-lived charged particle when it is the next lightest supersymmetric particle and is degenerate in mass with the lightest supersymmetric particle. The long-lived stau forms a bound state with a nucleus, and provide non-standard nuclear reactions. One of those, the internal conversion process, accelerates the destruction of ${}^{7}$Be and ${}^{7}$Li, and leads to a solution to the ${}^{7}$Li problem. On the other hand, the bound state of the stau and ${}^{4}$He enhances productions of n, d, t, and ${}^{6}$Li. The over-production of ${}^{6}$Li could solve the ${}^{6}$Li problem. While, the over-productions of d and t could conflict with observations, and the relevant parameter space of the stau is strictly constrained. We therefore need to carefully investigate the stau-${}^{4}$He bound state to find a condition of solving the ${}^{6}$Li problem. The scenario of the long-lived stau simultaneously and successfully fit the abundances of light elements (d, t, ${}^{3}$He, ${}^{4}$He, ${}^{6}$Li, and ${}^{7}$Li) and the neutralino dark matter to the observed ones. Consequently parameter space both of the stau and the neutralino is determined with excellent accuracy.

We have carried out a spectral analysis of the Suzaku X-ray data in the 0.4-12 keV range toward the shell-type very-high-energy {\gamma}-ray supernova remnant RX J1713.7-3946. The aims of this analysis are to estimate detailed X-rays spectral properties at a high angular resolution up to 2 arcmin, and to compare them with the interstellar gas. The X-ray spectrum is non-thermal and used to calculate absorbing column density, photon index, and absorption-corrected X-ray flux. The photon index varies significantly from 2.1 to 2.9. It is shown that the X-ray intensity is well correlated with the photon index, especially in the west region, with a correlation coefficient of 0.81. The X-ray intensity tends to increase with the averaged interstellar gas density while the dispersion is relatively large. The hardest spectra having the photon index less than 2.4 are found outside of the central 10 arcmin of the SNR, from the north to the southeast (~430 arcmin^2) and from the southwest to the northwest (~150 arcmin^2). The former region shows low interstellar gas density, while the latter high interstellar gas density. We present discussion for possible scenarios which explain the distribution of the photon index and its relationship with the interstellar gas.

We show that the Li problems can be solved in the next-to-minimal supersymmetric standard model where the slepton as the next-to-lightest SUSY particle is very long-lived. Such a long-lived slepton induces exotic nuclear reactions in big-bang nucleosynthesis, and destroys and produces the $^7$Li and $^6$Li nuclei via bound state formation. We study cases where the lightest SUSY particle is singlino-like neutralino and bino-like neutralino to present allowed regions in the parameter space which is consistent with the observations on the dark matter and the Higgs mass.

We have discovered two molecular features at radial velocities of -35 km/s and 0 km/s toward the infrared Double Helix Nebula (DHN) in the Galactic center with NANTEN2. The two features show good spatial correspondence with the DHN. We have also found two elongated molecular ridges at these two velocities distributed vertically to the Galactic plane over 0.8 degree. The two ridges are linked by broad features in velocity and are likely connected physically with each other. The ratio between the 12CO J=2-1 and J=1-0 transitions is 0.8 in the ridges which is larger than the average value 0.5 in the foreground gas, suggesting the two ridges are in the Galactic center. An examination of the K band extinction reveals a good coincidence with the CO 0 km/s ridge and is consistent with a distance of 8 +/-2 kpc. We discuss the possibility that the DHN was created by a magnetic phenomenon incorporating torsional Alfv\'en waves launched from the circumnuclear disk (Morris, Uchida & Do 2006) and present a first estimate of the mass and energy involved in the DHN.

RX J1713.7-3946 is the most remarkable very-high-energy \gamma-ray supernova remnant which emits synchrotron X-rays without thermal features. We made a comparative study of CO, HI and X-rays in order to better understand the relationship between the X-rays, and the molecular and atomic gas. The results indicate that the X-rays are enhanced around the CO and HI clumps on a pc scale but are decreased inside the clumps on a 0.1 pc scale. Magnetohydrodynamic numerical simulations of the shock interaction with molecular and atomic gas indicate that the interaction between the shock waves and the clumps excite turbulence which amplifies the magnetic field around the clumps (Inoue et al. 2012). We suggest that the amplified magnetic field around the CO and HI clumps enhances the synchrotron X-rays and possibly the acceleration of cosmic-ray electrons.

We studied allowed region on the intergenerational mixing parameters of sleptons from a viewpoint of big-bang nucleosynthesis in a slepton-neutralino coannihilation scenario. In this scenario, $^7$Li and $^6$Li problems can be solved by considering exotic reactions caused by bound-state effects with a long-lived slepton. Light element abundances are calculated as functions of the relic density and lifetime of the slepton which considerably depend on the intergenerational mixing parameters. Compared with observational light element abundances, we obtain allowed regions on the intergenerational mixing. Ratio of selectron component to stau component, $c_e$, is allowed in $2\times 10^{-11} \lesssim c_e \lesssim 2\times 10^{-9}$ with solving both the $^7$Li and $^6$Li problems. Similarly, the ratio for smuon, $c_{\mu}$, is allowed in $%2\times 10^{-7} \lesssim c_{\mu} \lesssim 5\times 10^{-5}$ for mass difference between slepton and neutralino, which is smaller than muon mass, and $ %10^{-11}\lesssim c_{\mu} \lesssim 2\times 10^{-10}$ for the mass difference in range between muon mass and 125 MeV. We also discuss collider signatures of the slepton decays. We find characteristic double peaks in momentum distribution of event number of the slepton decays with allowed mixing parameters. Discoveries of the double peaks at future collider experiments should confirm our scenario.

RX J1713.7$-$3946 is the most remarkable TeV $\gamma$-ray SNR which emits $\gamma$-rays in the highest energy range. We made a new combined analysis of CO and \ion{H}{1} in the SNR and derived the total protons in the interstellar medium (ISM). We have found that the inclusion of the \ion{H}{1} gas provides a significantly better spatial match between the TeV $\gamma$-rays and ISM protons than the H$_2$ gas alone. In particular, the southeastern rim of the $\gamma$-ray shell has a counterpart only in the \ion{H}{1}. The finding shows that the ISM proton distribution is consistent with the hadronic scenario that comic ray (CR) protons react with ISM protons to produce the $\gamma$-rays. This provides another step forward for the hadronic origin of the $\gamma$-rays by offering one of the necessary conditions missing in the previous hadronic interpretations. We argue that the highly inhomogeneous distribution of the ISM protons is crucial in the origin of the $\gamma$-rays. Most of the neutral gas was likely swept up by the stellar wind of an OB star prior to the SNe to form a low-density cavity and a swept-up dense wall. The cavity explains the low-density site where the diffusive shock acceleration of charged particles takes place with suppressed thermal X-rays, whereas the CR protons can reach the target protons in the wall to produce the $\gamma$-rays. The present finding allows us to estimate the total CR proton energy to be $\sim 10^{48}$ ergs, 0.1% of the total energy of a SNe.

We propose helium-4 spallation processes induced by long-lived stau in supersymmetric standard models, and investigate an impact of the processes on light elements abundances. We show that, as long as the phase space of helium-4 spallation processes is open, they are more important than stau-catalyzed fusion and hence constrain the stau property.

RX J1713.7-3946 is one of the TeV {\gamma}-ray supernova remnants (SNRs) emitting synchrotron X rays. The SNR is associated with molecular gas located at ~1 kpc. We made new molecular observations toward the dense cloud cores, peaks A, C and D, in the SNR in the 12CO(J=2-1) and 13CO(J=2-1) transitions at angular resolution of 90". The most intense core in 13CO, peak C, was also mapped in the 12CO(J=4-3) transition at angular resolution of 38". Peak C shows strong signs of active star formation including bipolar outflow and a far-infrared protostellar source and has a steep gradient with a r^{-2.2$\pm$0.4} variation in the average density within radius r. Peak C and the other dense cloud cores are rim-brightened in synchrotron X rays, suggesting that the dense cloud cores are embedded within or on the outer boundary of the SNR shell. This confirms the earlier suggestion that the X rays are physically associated with the molecular gas (Fukui et al. 2003). We present a scenario where the densest molecular core, peak C, survived against the blast wave and is now embedded within the SNR. Numerical simulations of the shock-cloud interaction indicate that a dense clump can indeed survive shock erosion, since shock propagation speed is stalled in the dense clump. Additionally, the shock-cloud interaction induces turbulence and magnetic field amplification around the dense clump that may facilitate particle acceleration in the lower-density inter-clump space leading to the enhanced synchrotron X rays around dense cores.

The Standing Accretion Shock Instability (SASI) is commonly believed to be responsible for large amplitude dipolar oscillations of the stalled shock during core collapse, potentially leading to an asymmetric supernovae explosion. The degree of asymmetry depends on the amplitude of SASI, which nonlinear saturation mechanism has never been elucidated. We investigate the role of parasitic instabilities as a possible cause of nonlinear SASI saturation. As the shock oscillations create both vorticity and entropy gradients, we show that both Kelvin-Helmholtz and Rayleigh-Taylor types of instabilities are able to grow on a SASI mode if its amplitude is large enough. We obtain simple estimates of their growth rates, taking into account the effects of advection and entropy stratification. In the context of the advective-acoustic cycle, we use numerical simulations to demonstrate how the acoustic feedback can be decreased if a parasitic instability distorts the advected structure. The amplitude of the shock deformation is estimated analytically in this scenario. When applied to the set up of Fernandez & Thompson (2009a), this saturation mechanism is able to explain the dramatic decrease of the SASI power when both the nuclear dissociation energy and the cooling rate are varied. Our results open new perspectives for anticipating the effect, on the SASI amplitude, of the physical ingredients involved in the modeling of the collapsing star.

We calculate the relic density of stau at the beginning of the Big-Bang Nucleosynthesis (BBN) era in the coannihilation scenario of minimal supersymmetric standard model (MSSM). In this scenario, stau can be long-lived and form bound states with nuclei. We put constraints on the parameter space of MSSM by connecting the calculation of the relic density of stau to the observation of the light elements abundance, which strongly depends on the relic density of stau. Consistency between the theoretical prediction and the observational result, both of the dark matter abundance and the light elements abundance, requires the mass difference between the lighter stau and the lightest neutralino to be around 100MeV, the stau mass to be 300 -- 400 GeV, and the mixing angle of the left and right-handed staus to be $\sin\theta_{\tau} = (0.65 \textrm{--} 1)$.

The physical processes involved in the advective-acoustic instability are investigated with 2D numerical simulations. Simple toy models, developped in a companion paper, are used to describe the coupling between acoustic and entropy/vorticity waves, produced either by a stationary shock or by the deceleration of the flow. Using two Eulerian codes based on different second order upwind schemes, we confirm the results of the perturbative analysis. The numerical convergence with respect to the computation mesh size is studied with 1D simulations. We demonstrate that the numerical accuracy of the quantities which depend on the physics of the shock is limited to a linear convergence. We argue that this property is likely to be true for most current numerical schemes dealing with SASI in the core-collapse problem, and could be solved by the use of advanced techniques for the numerical treatment of the shock. We propose a strategy to choose the mesh size for an accurate treatment of the advective-acoustic coupling in future numerical simulations.

A scenario of the Big-Bang Nucleosynthesis is analyzed within the Minimal Supersymmetric Standard Model which is consistent with a stau-neutralino coannihilation scenario to explain the relic abundance of dark matter. We find that we can account for the possible descrepancy of the abundance of $\mathrm{^{7}Li}$ between the observation and the prediction of the Big-Bang Nucleosynthesis by taking the mass of the neutralino as $300 \mathrm{GeV}$ and the mass difference between the stau and the neutralino as $(100 -- 120) MeV$. We can therefore simultaneously explain the abundance of the dark matter and that of $\mathrm{^{7}Li}$ by these values of parameters. The lifetime of staus in this scenario is predicted to be $O(100 -- 1000) sec$.

We consider the effect of a supernova (SN) explosion in a very massive binary that is expected to form in a portion of Population III stars with the mass higher than 100$M_\odot$. In a Population III binary system, a more massive star can result in the formation of a BH and a surrounding accretion disc. Such BH accretion could be a significant source of the cosmic reionization in the early universe. However, a less massive companion star evolves belatedly and eventually undergoes a SN explosion, so that the accretion disc around a BH might be blown off in a lifetime of companion star. In this paper, we explore the dynamical impact of a SN explosion on an accretion disc around a massive BH, and elucidate whether the BH accretion disc is totally demolished or not. For the purpose, we perform three-dimensional hydrodynamic simulations of a very massive binary system, where we assume a BH of $10^3 M_{\odot}$ that results from a direct collapse of a very massive star and a companion star of $100 M_{\odot}$ that undergoes a SN explosion. We calculate the remaining mass of a BH accretion disc as a function of time. As a result, it is found that a significant portion of gas disc can survive through three-dimensional geometrical effects even after the SN explosion of a companion star. Even if the SN explosion energy is higher by two orders of magnitude than the binding energy of gas disc, about a half of disc can be left over. The results imply that the Population III BH accretion disc can be a long-lived luminous source, and therefore could be an important ionizing source in the early universe.

Modification of standard big-bang nucleosynthesis is considered in the minimal supersymmetric standard model to resolve the excessive theoretical prediction of the abundance of primordial lithium 7. We focus on the stau as a next-lightest superparticle, which is long lived due to its small mass difference with the lightest superparticle. It provides a number of additional decay processes of $\mathrm{^{7}Li}$ and $\mathrm{^{7}Be}$. A particularly important process is the internal conversion in the stau-nucleus bound state, which destroys the $\mathrm{^{7}Li}$ and $\mathrm{^{7}Be}$ effectively. We show that the modification can lead to a prediction consistent with the observed abundance of $\mathrm{^{7}Li}$.

We propose a novel method to extract the imprint of gravitational lensing by Pop III stars in the light curves of Gamma Ray Bursts (GRBs). Significant portions of GRBs can originate in hypernovae of Pop III stars and be gravitationally lensed by foreground Pop III stars or their remnants. If the lens mass is on the order of $10^2-10^3M_\odot$ and the lens redshift is greater than 10, the time delay between two lensed images of a GRB is $\approx 1$s and the image separation is $\approx 10 \mu$as. Although it is difficult to resolve the two lensed images spatially with current facilities, the light curves of two images are superimposed with a delay of $\approx 1$ s. GRB light curves usually exhibit noticeable variability, where each spike is less than 1s. If a GRB is lensed, all spikes are superimposed with the same time delay. Hence, if the autocorrelation of light curve with changing time interval is calculated, it should show the resonance at the time delay of lensed images. Applying this autocorrelation method to GRB light curves which are archived as the {\it BATSE} catalogue, we demonstrate that more than half light curves can show the recognizable resonance, if they are lensed. Furthermore, in 1821 GRBs we actually find one candidate of GRB lensed by a Pop III star, which may be located at redshift 20-200. The present method is quite straightforward and therefore provides an effective tool to search for Pop III stars at redshift greater than 10. Using this method, we may find more candidates of GRBs lensed by Pop III stars in the data by the {\it Swift} satellite.

We study the stau lifetime in a scenario with the LSP taken to be a neutralino and the NLSP being a stau, based on the minimal supersymmetric Standard Model. The mass difference between the LSP and NLSP, $\delta m$, must satisfy $\delta m/m_{\tilde{\chi}} \sim$ a few % or less for coannihilation to occur, where $m_{\tilde{\chi}}$ is the neutralino mass. We calculate the stau lifetime from the decay modes $\tilde{\tau}\to \tilde{\chi}\tau$, $\tilde{\chi}\nu_\tau\pi$, and $\tilde{\chi}\nu_\tau\mu(e)\nu_{\mu(e)}$ and discuss its dependence on various parameters. We find that the lifetime is in the range $10^{-22}$--$10^{16}$ sec for $10^{-2} \le \delta m \le 10$ GeV. We also discuss the connection with lepton flavor violation if there is mixing between sleptons.

We perform two-dimensional hydrodynamic simulation on the gas stripping by radiation drag from an interstellar cloud moving in uniform radiation fields. To properly include relativistic radiation drag, the radiation hydrodynamic equation is solved with taking into account the dilution of radiation fields by optical depth of the cloud. As a result, it is found that the optically-thin surface layers are effectively stripped by radiation drag from an optically-thick gas cloud, and simultaneously stripped gas loses momentum. The momentum loss time-scale is found to be on the order of $10^{8}$ years under intensive radiation fields which are expected in the early phase of galaxy evolution. The present results show that the radiation drag is an effective mechanism to extract angular momentum from interstellar medium and allows it to accrete onto the galactic center. The mass accretion driven by radiation drag may lead to the formation of a central supermassive black hole.

We compare the photometric accuracy of spectra and images in flares observed with the Ramaty High Energy Solar Spectroscopic Imager (RHESSI)}spacecraft. We test the accuracy of the photometry by comparing the photon fluxes obtained in different energy ranges from the spectral-fitting software SPEX with those fluxes contained in the images reconstructed with the Clean, MEM, MEM-Vis, Pixon, and Forward-fit algorithms. We quantify also the background fluxes, the fidelity of source geometries, and spatial spectra reconstructed with the five image reconstruction algorithms. We investigate the effects of grid selection, pixel size, field-of-view, and time intervals on the quality of image reconstruction. The detailed parameters and statistics are provided in an accompanying CD-ROM and web page. We find that Forward-fit, Pixon, and Clean have a robust convergence behavior and a photometric accuracy in the order of a few percents, while MEM does not converge optimally for large degrees of freedom (for large field-of-views and/or small pixel sizes), and MEM-Vis suffers in the case of time-variable sources. This comparative study documents the current status of the RHESSI spectral and imaging software, one year after launch.

Two-dimensional numerical simulations of an accretion flow in a close binary system are performed by solving the Euler equations with radiative transfer. In the present study, the specific heat ratio is assumed to be constant while radiative cooling effect is included as a non-adiabatic process. The cooling effect of the disc is considered by discharging energy in the vertical directions from the top and bottom surfaces of the disc. We use the flux-limited diffusion approximation to calculate the radiative heat flux values. Our calculations show that a disc structure appears and the spiral shocks are formed on the disc. These features are similar to that observed in the case of an adiabatic gas with a lower specific heat ratio, $\gamma=1.01$. It is found that when radiative cooling effect is accounted for, the mass of the disc becomes larger than that assuming $\gamma=5/3$, and smaller than that assuming $\gamma=1.01$. It is concluded that employing an adiabatic gas with a lower specific heat ratio is almost a valid assumption for simulating accretion disc with radiative cooling effect.

We report the first measurement of genus curves for the two-dimensional mass map in the neighborhood of rich, X-ray luminous galaxy cluster MS1054-03 at z=0.83, reconstructed from weak lensing data obtained by Suprime-Cam on the prime focus of 8.2m Subaru telescope. We find that the genus curve measured in the whole survey field deviates from that expected from a random Gaussian field. We show that this non-Gaussianity is induced by the rich cluster in this region, and that the genus curve for the region without the cluster is consistent with the prediction for a random Gaussian field. The results demonstrate the feasibility of the genus statistics to examine the non-Gaussianity due to the large scale structures and to probe the statistical properties of the large scale clustering.

We apply the morphological descriptions of two-dimensional contour map, the so-called Minkowski functionals (the area fraction, circumference, and Euler characteristics), to the convergence field $\kappa(\bm{\theta})$ of the large-scale structure reconstructed from the shear map produced by the ray-tracing simulations. The perturbation theory of structure formation has suggested that the non-Gaussian features on the Minkowski functionals with respect to the threshold in the weakly nonlinear regime are induced by the three skewness parameters of $\kappa$ that are sensitive to the density parameter of matter, $\Omega_{\rm m}$. We show that, in the absence of noise due to the intrinsic ellipticities of source galaxies with which the perturbation theory results can be recovered, the accuracy of $\Omega_{\rm m}$ determination is improved by $\sim 20%$ using the Minkowski functionals compared to the conventional method of using the direct measure of skewness.