Spatial correlations among proto-planetary disk orientations carry unique information on physics of multiple star formation processes. We select five nearby star-forming regions that comprise a number of proto-planetary disks with spatially-resolved images with ALMA and HST, and search for the mutual alignment of the disk axes. Specifically, we apply the Kuiper test to examine the statistical uniformity of the position angle (PA: the angle of the major axis of the projected disk ellipse measured counter-clockwise from the north) distribution. The disks located in the star-forming regions, except the Lupus clouds, do not show any signature of the alignment, supporting the random orientation. Rotational axes of 16 disks with spectroscopic measurement of PA in the Lupus III cloud, a sub-region of the Lupus field, however, exhibit a weak and possible departure from the random distribution at a $2\sigma$ level, and the inclination angles of the 16 disks are not uniform as well. Furthermore, the mean direction of the disk PAs in the Lupus III cloud is parallel to the direction of its filament structure, and approximately perpendicular to the magnetic field direction. We also confirm the robustness of the estimated PAs in the Lupus clouds by comparing the different observations and estimators based on three different methods including sparse modeling. The absence of the significant alignment of the disk orientation is consistent with the turbulent origin of the disk angular momentum. Further observations are required to confirm/falsify the possible disk alignment in the Lupus III cloud.

Blazars can be divided into two subtypes, flat spectrum radio quasars (FSRQs) and BL Lac objects, which have been distinguished phenomenologically by the strength of their optical emission lines, while their physical nature and relationship are still not fully understood. In this paper, we focus on the differences in their variability. We characterize the blazar variability using the Ornstein-Uhlenbeck (OU) process, and investigate the features that are discriminative for the two subtypes. We used optical photometric and polarimetric data obtained with the 1.5-m Kanata telescope for 2008-2014. We found that four features, namely the variation amplitude, characteristic timescale, and non-stationarity of the variability obtained from the light curves and the median of the degree of polarization (PD), are essential for distinguishing between FSRQs and BL Lac objects. FSRQs are characterized by rare and large flares, while the variability of BL Lac objects can be reproduced with a stationary OU process with relatively small amplitudes. The characteristics of the variability are governed not by the differences in the jet structure between the subtypes, but by the peak frequency of the synchrotron emission. This implies that the nature of the variation in the jets is common in FSRQs and BL Lac objects. We found that BL Lac objects tend to have high PD medians, which suggests that they have a stable polarization component. FSRQs have no such component, possibly because of a strong Compton cooling effect in sub-pc scale jets.

With an emphasis on improving the fidelity even in super-resolution regimes, new imaging techniques have been intensively developed over the last several years, which may provide substantial improvements to the interferometric observation of protoplanetary disks. In this study, sparse modeling (SpM) is applied for the first time to observational data sets taken by the Atacama Large Millimeter/submillimeter Array (ALMA). The two data sets used in this study were taken independently using different array configurations at Band 7 (330 GHz), targeting the protoplanetary disk around HD 142527; one in the shorter-baseline array configuration (~ 430 m), and the other in the longer-baseline array configuration (~ 1570 m). The image resolutions reconstructed from the two data sets are different by a factor of ~ 3. We confirm that the previously known disk structures appear on the images produced by both SpM and CLEAN at the standard beam size. The image reconstructed from the shorter-baseline data using the SpM matches that obtained with the longer-baseline data using CLEAN, achieving a super-resolution image from which a structure finer than the beam size can be reproduced. Our results demonstrate that on-going intensive development in the SpM imaging technique is beneficial to imaging with ALMA.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

We investigate the cosmological applications of scalar-tensor theories that arise effectively from the Lorentz fiber bundle of a Finsler-like geometry. We first show that the involved nonlinear connection induces a new scalar degree of freedom and eventually a scalar-tensor theory. Using both a holonomic and a nonholonomic basis, we show the appearance of an effective dark-energy sector, which additionally acquires an explicit interaction with the matter sector, arising purely from the internal structure of the theory. Applying the theory at late times we find that we can obtain the thermal history of the Universe, namely the sequence of matter and dark-energy epochs, and moreover the effective dark-energy equation-of-state parameter can be quintessencelike, phantomlike, or experience the phantom-divide crossing during cosmological evolution. Furthermore, applying the scenario at early times, we see that one can acquire an exponential de Sitter solution as well as obtain an inflationary realization with the desired scale-factor evolution. These features arise purely from the intrinsic geometrical structure of Finsler-like geometry and reveal the capabilities of the construction.

Galaxy imaging surveys enable us to map the cosmic matter density field through weak gravitational lensing analysis. The density reconstruction is compromised by a variety of noise originating from observational conditions, galaxy number density fluctuations, and intrinsic galaxy properties. We propose a deep-learning approach based on generative adversarial networks (GANs) to reduce the noise in the weak lensing map under realistic conditions. We perform image-to-image translation using conditional GANs in order to produce noiseless lensing maps using the first-year data of the Subaru Hyper Suprime-Cam (HSC) survey. We train the conditional GANs by using 30000 sets of mock HSC catalogs that directly incorporate observational effects. We show that an ensemble learning method with GANs can reproduce the one-point probability distribution function (PDF) of the lensing convergence map within a $0.5-1\sigma$ level. We use the reconstructed PDFs to estimate a cosmological parameter $S_{8} = \sigma_{8}\sqrt{\Omega_{\rm m0}/0.3}$, where $\Omega_{\rm m0}$ and $\sigma_{8}$ represent the mean and the scatter in the cosmic matter density. The reconstructed PDFs place tighter constraint, with the statistical uncertainty in $S_8$ reduced by a factor of $2$ compared to the noisy PDF. This is equivalent to increasing the survey area by $4$ without denoising by GANs. Finally, we apply our denoising method to the first-year HSC data, to place $2\sigma$-level cosmological constraints of $S_{8} < 0.777 \, ({\rm stat}) + 0.105 \, ({\rm sys})$ and $S_{8} < 0.633 \, ({\rm stat}) + 0.114 \, ({\rm sys})$ for the noisy and denoised data, respectively.

Using a prototype of the Tomo-e Gozen wide-field CMOS mosaic camera, we acquire wide-field optical images at a cadence of 2 Hz and search them for transient sources of duration 1.5 to 11.5 seconds. Over the course of eight nights, our survey encompasses the equivalent of roughly two days on one square degree, to a fluence equivalent to a limiting magnitude about $V = 15.6$ in a 1-second exposure. After examining by eye the candidates identified by a software pipeline, we find no sources which meet all our criteria. We compute upper limits to the rate of optical transients consistent with our survey, and compare those to the rates expected and observed for representative sources of ephemeral optical light.

We report observations of a stellar occultation by the classical Kuiper belt object (50000) Quaoar occurred on 28 June 2019. A single-chord high-cadence (2 Hz) photometry dataset was obtained with the Tomo-e Gozen CMOS camera mounted on the 1.05 m Schmidt telescope at Kiso Observatory. The obtained ingress and egress data do not show any indication of atmospheric refraction and allow to set new $1\sigma$ and $3\sigma$ upper limits of 6 and 16 nbar, respectively, for the surface pressure of a pure methane atmosphere. These upper limits are lower than the saturation vapor pressure of methane at Quaoar's expected mean surface temperature ($T \sim 44$ K) and imply the absence of a $\sim$10 nbar-level global atmosphere formed by methane ice on Quaoar's surface.

We present a concept of the X-ray imaging system with high angular-resolution and moderate sensitivity. In this concept, a two-dimensional detector, i.e., imager, is put at a slightly out-of-focused position of the focusing mirror, rather than just at the mirror focus as in the standard optics, to capture the miniature image of objects. In addition, a set of multi-grid masks (or a modulation collimator) is installed in front of the telescope. We find that the masks work as a coded aperture camera and that they boost the angular resolution of the focusing optics. The major advantage of this concept is that a much better angular resolution an order of 2--3 or more than in the conventional optics is achievable, while a high throughput (large effective area) is maintained, which is crucial in photon-limited high-energy astronomy, because any type of mirrors, including light-weight reflective mirrors, can be employed in our concept. If the signal-to-noise ratio is sufficiently high, we estimate that angular resolutions at the diffraction limit of 4" and 0"4 at about 7 keV can be achieved with a pair of masks at separations of 1~m and 100~m, respectively, at the diffraction limit.

Weak gravitational lensing is a powerful probe of the large-scale cosmic matter distribution. Wide-field galaxy surveys allow us to generate the so-called weak lensing maps, but actual observations suffer from noise due to imperfect measurement of galaxy shape distortions and to the limited number density of the source galaxies. In this paper, we explore a deep-learning approach to reduce the noise. We develop an image-to-image translation method with conditional adversarial networks (CANs), which learn efficient mapping from an input noisy weak lensing map to the underlying noise field. We train the CANs using $30000$ image pairs obtained from $1000$ ray-tracing simulations of weak gravitational lensing. We show that the trained CANs reproduce the true one-point probability distribution function (PDF) of the noiseless lensing map with a bias less than $1\sigma$ on average, where $\sigma$ is the statistical error. We perform a Fisher analysis to make forecast for cosmological parameter inference with the one-point lensing PDF. By our denoising method using CANs, the first derivative of the PDF with respect to the cosmic mean matter density and the amplitude of the primordial curvature perturbations becomes larger by $\sim50\%$. This allows us to improve the cosmological constraints by $\sim30-40\%$ with using observational data from ongoing and upcoming galaxy imaging surveys.

We present an overview of a deep transient survey of the COSMOS field with the Subaru Hyper Suprime-Cam (HSC). The survey was performed for the 1.77 deg$^2$ ultra-deep layer and 5.78 deg$^2$ deep layer in the Subaru Strategic Program over 6- and 4-month periods from 2016 to 2017, respectively. The ultra-deep layer shows a median depth per epoch of 26.4, 26.3, 26.0, 25.6, and 24.6 mag in $g$, $r$, $i$, $z$, and $y$ bands, respectively; the deep layer is $\sim0.6$ mag shallower. In total, 1,824 supernova candidates were identified. Based on light curve fitting and derived light curve shape parameter, we classified 433 objects as Type Ia supernovae (SNe); among these candidates, 129 objects have spectroscopic or COSMOS2015 photometric redshifts and 58 objects are located at $z > 1$. Our unique dataset doubles the number of Type Ia SNe at $z > 1$ and enables various time-domain analyses of Type II SNe, high redshift superluminous SNe, variable stars, and active galactic nuclei.

Faraday tomography (or rotation measure synthesis) is a procedure to convert linear polarization spectra into the Faraday dispersion function, which provides us with unique information of magneto-ionic media along the line of sight. Mathematical formulation of Faraday tomography is similar to polarimetric imaging of radio interferometry, where many new methods have been actively developed and shown to outperform the standard CLEAN approaches. In this paper, we propose a sparse reconstruction technique to Faraday tomography. This technique is being developed for interferometric imaging and utilizes computationally less expensive convex regularization functions such as $\ell_1$-norm and total variation (TV) or total squared variation (TSV). The proposed technique solves a convex optimization, and therefore its solution is determined uniquely regardless of the initial condition for given regularization parameters that can be optimized by data themselves. Using a physically-motivated model of turbulent galactic magnetized plasma, we demonstrate that the proposed technique outperforms RM-CLEAN and provides higher-fidelity reconstruction. The proposed technique would be a powerful tool in broadband polarimetry with the Square Kilometre Array (SKA) and its precursors.

We propose an image reconstruction method for an X-ray telescope system with an angular resolution booster proposed by Maeda et al.(2018). The system consists of double multi-grid masks in front of an X-ray mirror and an off-focused two-dimensional imager. Because the obtained image is off-focused, additional image reconstruction process is assumed to be included. Our image reconstruction method is an extension of the traditional Richardson-Lucy algorithm with two regularization terms, one for sparseness and the other for smoothness. Such a combination is desirable for astronomical imaging because astronomical objects have variety in shape from point sources, diffuse sources to mixtures of them. The performance of the system is demonstrated with simulated data for point sources and diffused X-ray sources such as Cas A and Crab Nebula. The image resolution is improved from a few arcmin of focused image without the booster to a few arcsec with the booster. Through the demonstration, the angular resolution booster with the image reconstruction method is shown to be feasible.

Imaging observations of faint meteors were carried out on April 11 and 14, 2016 with a wide-field CMOS mosaic camera, Tomo-e PM, mounted on the 105-cm Schmidt telescope at Kiso Observatory, the University of Tokyo. Tomo-e PM, which is a prototype model of Tomo-e Gozen, can monitor a sky of ${\sim}1.98\,\mathrm{deg^2}$ at 2\,Hz. The numbers of detected meteors are 1514 and 706 on April 11 and 14, respectively. The detected meteors are attributed to sporadic meteors. Their absolute magnitudes range from $+4$ to $+10\,\mathrm{mag}$ in the $V$-band, corresponding to about $8.3{\times}10^{-2}$ to $3.3{\times}10^{-4}\,\mathrm{g}$ in mass. The present magnitude distributions we obtained are well explained by a single power-law luminosity function with a slope parameter $r = 3.1{\pm}0.4$ and a meteor rate $\log_{10}N_0 = -5.5{\pm}0.5$. The results demonstrate a high performance of telescopic observations with a wide-field video camera to constrain the luminosity function of faint meteors. The performance of Tomo-e Gozen is about two times higher than that of Tomo-e PM. A survey with Tomo-e Gozen will provide a more robust measurement of the luminosity function.

Chopping observations with a tip-tilt secondary mirror have conventionally been used in ground-based mid-infrared observations. However, it is not practical for next generation large telescopes to have a large tip-tilt mirror that moves at a frequency larger than a few Hz. We propose an alternative observing method, a "slow-scanning" observation. Images are continuously captured as movie data, while the field-of-view is slowly moved. The signal from an astronomical object is extracted from the movie data by a low-rank and sparse matrix decomposition. The performance of the "slow-scanning" observation was tested in an experimental observation with Subaru/COMICS. The quality of a resultant image in the "slow-scanning" observation was as good as in a conventional chopping observation with COMICS, at least for a bright point-source object. The observational efficiency in the "slow-scanning" observation was better than that in the chopping observation. The results suggest that the "slow-scanning" observation can be a competitive method for the Subaru telescope and be of potential interest to other ground-based facilities to avoid chopping.

We propose a new superresolution imaging technique for interferometry using sparse modeling, utilizing two regularization terms: the $\ell_1$-norm and a new function named Total Squared Variation (TSV) of the brightness distribution. TSV is an edge-smoothing variant of Total Variation (TV), leading to reducing the sum of squared gradients. First, we demonstrate that our technique may achieve super-resolution of $\sim 30$% compared to the traditional CLEAN beam size using synthetic observations of two point sources. Second, we present simulated observations of three physically motivated static models of Sgr A* with the Event Horizon Telescope (EHT) to show the performance of proposed techniques in greater detail. We find that $\ell_1$+TSV regularization outperforms $\ell_1$+TV regularization with the popular isotropic TV term and the Cotton-Schwab CLEAN algorithm, demonstrating that TSV is well-matched to the expected physical properties of the astronomical images, which are often nebulous. Remarkably, in both the image and gradient domains, the optimal beam size minimizing root-mean-squared errors is $\lesssim 10$% of the traditional CLEAN beam size for $\ell_1$+TSV regularization, and non-convolved reconstructed images have smaller errors than beam-convolved reconstructed images. This indicates that the traditional post-processing technique of Gaussian convolution in interferometric imaging may not be required for the $\ell_1$+TSV regularization. We also propose a feature extraction method to detect circular features from the image of a black hole shadow with the circle Hough transform (CHT) and use it to evaluate the performance of the image reconstruction. With our imaging technique and the CHT, the EHT can constrain the radius of the black hole shadow with an accuracy of $\sim 10-20$% in present simulations for Sgr A*.

We develop an approximation formula for the cross-validation error (CVE) of a sparse linear regression penalized by $\ell_1$-norm and total variation terms, which is based on a perturbative expansion utilizing the largeness of both the data dimensionality and the model. The developed formula allows us to reduce the necessary computational cost of the CVE evaluation significantly. The practicality of the formula is tested through application to simulated black-hole image reconstruction on the event-horizon scale with super resolution. The results demonstrate that our approximation reproduces the CVE values obtained via literally conducted cross-validation with reasonably good precision.

We propose a new technique for radio interferometry to obtain super-resolution full polarization images in all four Stokes parameters using sparse modeling. The proposed technique reconstructs the image in each Stokes parameter from the corresponding full-complex Stokes visibilities by utilizing two regularization functions: the $\ell _1$-norm and total variation (TV) of the brightness distribution. As an application of this technique, we present simulated linear polarization observations of two physically motivated models of M87 with the Event Horizon Telescope (EHT). We confirm that $\ell _1$+TV regularization can achieve an optimal resolution of $\sim 25-30$\% of the diffraction limit $\lambda/D_{\rm max}$, which is the nominal spatial resolution of a radio interferometer for both the total intensity (i.e. Stokes $I$) and linear polarizations (i.e. Stokes $Q$ and $U$). This optimal resolution is better than that obtained from the widely used Cotton-Schwab CLEAN algorithm or from using $\ell _1$ or TV regularizations alone. Furthermore, we find that $\ell _1$+TV regularization can achieve much better image fidelity in linear polarization than other techniques over a wide range of spatial scales, not only in the super-resolution regime, but also on scales larger than the diffraction limit. Our results clearly demonstrate that sparse reconstruction is a useful choice for high-fidelity full-polarimetric interferometric imaging.

We propose a new imaging technique for radio and optical/infrared interferometry. The proposed technique reconstructs the image from the visibility amplitude and closure phase, which are standard data products of short-millimeter very long baseline interferometers such as the Event Horizon Telescope (EHT) and optical/infrared interferometers, by utilizing two regularization functions: the $\ell_1$-norm and total variation (TV) of the brightness distribution. In the proposed method, optimal regularization parameters, which represent the sparseness and effective spatial resolution of the image, are derived from data themselves using cross validation (CV). As an application of this technique, we present simulated observations of M87 with the EHT based on four physically motivated models. We confirm that $\ell_1$+TV regularization can achieve an optimal resolution of $\sim 20-30$% of the diffraction limit $\lambda/D_{\rm max}$, which is the nominal spatial resolution of a radio interferometer. With the proposed technique, the EHT can robustly and reasonably achieve super-resolution sufficient to clearly resolve the black hole shadow. These results make it promising for the EHT to provide an unprecedented view of the event-horizon-scale structure in the vicinity of the super-massive black hole in M87 and also the Galactic center Sgr A*.

We apply three data science techniques, Nonnegative Matrix Factorization (NMF), Principal Component Analysis (PCA) and Independent Component Analysis (ICA), to simulated X-ray energy spectra of a particular class of super-massive black holes. Two competing physical models, one whose variable components are additive and the other whose variable components are multiplicative, are known to successfully describe X-ray spectral variation of these super-massive black holes, within accuracy of the contemporary observation. We hope to utilize these techniques to compare the viability of the models by probing the mathematical structure of the observed spectra, while comparing advantages and disadvantages of each technique. We find that PCA is best to determine the dimensionality of a dataset, while NMF is better suited for interpreting spectral components and comparing them in terms of the physical models in question. ICA is able to reconstruct the parameters responsible for spectral variation. In addition, we find that the results of these techniques are sufficiently different that applying them to observed data may be a useful test in comparing the accuracy of the two spectral models.

Optical wide-field surveys with a high cadence are expected to create a new field of astronomy, so-called "movie astronomy," in the near future. The amount of data of the observations will be huge, and hence efficient data compression will be indispensable. Here we propose a low-rank matrix approximation with sparse matrix decomposition as a promising solution to reduce the data size effectively, while preserving sufficient scientific information. We apply one of the methods to the movie data obtained with the prototype model of the Tomo-e Gozen mounted on the 1.0-m Schmidt telescope of Kiso Observatory. Once the full-scale observation of the Tomo-e Gozen commences, it will generate ~30 TB of data per night. We demonstrate that the data are compressed by a factor of about 10 in size without losing transient events like optical short transient point-sources and meteors. The intensity of point sources can be recovered from the compressed data. The processing runs sufficiently fast, compared with the expected data-acquisition rate in the actual observing runs.

We present an application of machine-learning (ML) techniques to source selection in the optical transient survey data with Hyper Suprime-Cam (HSC) on the Subaru telescope. Our goal is to select real transient events accurately and in a timely manner out of a large number of false candidates, obtained with the standard difference-imaging method. We have developed the transient selector which is based on majority voting of three ML machines of AUC Boosting, Random Forest, and Deep Neural Network. We applied it to our observing runs of Subaru-HSC in 2015 May and August, and proved it to be efficient in selecting optical transients. The false positive rate was 1.0% at the true positive rate of 90% in the magnitude range of 22.0--25.0 mag for the former data. For the latter run, we successfully detected and reported ten candidates of supernovae within the same day as the observation. From these runs, we learned the following lessons: (1) the training using artificial objects is effective in filtering out false candidates, especially for faint objects, and (2) combination of ML by majority voting is advantageous.

For short-wavelength VLBI observations, it is difficult to measure the phase of the visibility function accurately. The closure phases are reliable measurements under this situation, though it is not sufficient to retrieve all of the phase information. We propose a new method, Phase Retrieval from Closure Phase (PRECL). PRECL estimates all the visibility phases only from the closure phases. Combining PRECL with a sparse modeling method we have already proposed, imaging process of VLBI does not rely on dirty image nor self-calibration. The proposed method is tested numerically and the results are promising.

We discuss what is an appropriate set of explanatory variables in order to predict the absolute magnitude at the maximum of Type Ia supernovae. In order to have a good prediction, the error for future data, which is called the "generalization error," should be small. We use cross-validation in order to control the generalization error and LASSO-type estimator in order to choose the set of variables. This approach can be used even in the case that the number of samples is smaller than the number of candidate variables. We studied the Berkeley supernova database with our approach. Candidates of the explanatory variables include normalized spectral data, variables about lines, and previously proposed flux-ratios, as well as the color and light-curve widths. As a result, we confirmed the past understanding about Type Ia supernova: i) The absolute magnitude at maximum depends on the color and light-curve width. ii) The light-curve width depends on the strength of Si II. Recent studies have suggested to add more variables in order to explain the absolute magnitude. However, our analysis does not support to add any other variables in order to have a better generalization error.

We propose a new technique to obtain super-resolution images with radio interferometer using sparse modeling. In standard radio interferometry, sampling of ($u$, $v$) is quite often incomplete and thus obtaining an image from observed visibilities becomes an underdetermined problem, and a technique so-called "zero-padding" is often used to fill up unsampled grids in ($u$, $v$) plane, resulting in image degradation by finite beam size as well as numerous side-lobes. In this paper we show that directly solving such an underdetermined problem based on sparse modeling (in this paper LASSO) avoids the above problems introduced by zero-padding, leading to super-resolution images in which structure finer than the standard beam size (diffraction limit) can be reproduced. We present results of one-dimensional and two-dimensional simulations of interferometric imaging, and discuss its implications to super-resolution imaging, particularly focusing on imaging of black hole shadows with millimeter VLBI.

We study the image reconstruction problem of a Compton camera which consists of semiconductor detectors. The image reconstruction is formulated as a statistical estimation problem. We employ a bin-mode estimation (BME) and extend an existing framework to a Compton camera with multiple scatterers and absorbers. Two estimation algorithms are proposed: an accelerated EM algorithm for the maximum likelihood estimation (MLE) and a modified EM algorithm for the maximum a posteriori (MAP) estimation. Numerical simulations demonstrate the potential of the proposed methods.

We report on a proton radiation damage experiment on P-channel CCD newly developed for an X-ray CCD camera onboard the Astro-H satellite. The device was exposed up to 10^9 protons cm^{-2} at 6.7 MeV. The charge transfer inefficiency (CTI) was measured as a function of radiation dose. In comparison with the CTI currently measured in the CCD camera onboard the Suzaku satellite for 6 years, we confirmed that the new type of P-channel CCD is radiation tolerant enough for space use. We also confirmed that a charge-injection technique and lowering the operating temperature efficiently work to reduce the CTI for our device. A comparison with other P-channel CCD experiments is also discussed.

Monte Carlo simulation is one of the best tools to study the complex spectra of Compton-thick AGNs and to figure out the relation between their nuclear structures and X-ray spectra. We have simulated X-ray spectra of Compton-thick AGNs obscured by an accretion torus whose structure is characterized by a half-opening angle, an inclination angle of the torus relative to the observer, and a column density along the equatorial plane. We divided the simulated spectra into three components: one direct component, an absorbed reflection component and an unabsorbed reflection component. We then deduced the dependencies of these components on the parameters describing the structure of the torus. Our simulation results were applied to fit the wide-band spectrum of the Seyfert 2 galaxy Mrk 3 obtained by $Suzaku$. The spectral analysis indicates that we observe the nucleus along a line of sight intercepting the torus near its edge, and the column density along the equatorial plane was estimated to be ~10^24 cm^-2. Using this model, we can estimate the luminosities of both the direct emission and the emission irradiating the surrounding matter. This is useful to find the time variability and time lag between the direct and reflected light.

We obtained a wide-band spectrum of the Compton-thick Seyfert 2 galaxy Mrk 3 with Suzaku. The observed spectrum was clearly resolved into weak, soft power-law emission, a heavily absorbed power-law component, cold reflection, and many emission lines. The heavily absorbed component, absorbed by gas with a column density of 1.1x10^24 cm^-2, has an intrinsic 2--10 keV luminosity of ~1.6x10^43 erg s^-1, and is considered to be direct emission from the Mrk 3 nucleus. The reflection component was interpreted as reflection of the direct light off cold, thick material; the reflection fraction $R$ was 1.36+/-0.20. The cold material is inferred to be located > 1 pc from the central black hole of Mrk 3 due to the low ionization parameter of iron (xi < 1 erg cm s^-1) and the narrow iron line width (s < 22 eV). A Compton shoulder to the iron line was detected, but the intensity of the shoulder component was less than that expected from spherically distributed Compton-thick material. The weak, soft power-law emission is considered to be scattered light by ionized gas. The existence of many highly-ionized lines from O, Ne, Mg, Si, S, and Fe in the observed spectrum indicates that the ionized gas has a broad ionized structure, with xi=10--1000. The scattering fraction with respect to the direct light was estimated to be 0.9+/-0.2%, which indicates that the column density of the scattering region is about 3.6x10^22 cm^-2. This high-quality spectrum obtained by Suzaku can be considered a template for studies of Seyfert 2 galaxies.