Researcher profile

Hirokazu Odaka

Hirokazu Odaka contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

Event reconstruction of Compton telescopes using a multi-task neural network

We have developed a neural network model to perform event reconstruction of Compton telescopes. This model reconstructs events that consist of three or more interactions in a detector. It is essential for Compton telescopes to determine the time order of the gamma-ray interactions and whether the incident photon deposits all energy in a detector or it escapes from the detector. Our model simultaneously predicts these two essential factors using a multi-task neural network with three hidden layers of fully connected nodes. For verification, we have conducted numerical experiments using Monte Carlo simulation, assuming a large-area Compton telescope using liquid argon to measure gamma rays with energies up to $3.0\,\mathrm{MeV}$. The reconstruction model shows excellent performance of event reconstruction for multiple scattering events that consist of up to eight hits. The accuracies of hit order prediction are around $60\%$ while those of escape flags are higher than $70\%$ for up to eight-hit events of $4π$ isotropic photons. Compared with two other algorithms, a classical model and a physics-based probabilistic one, the present neural network method shows high performance in estimation accuracy particularly when the number of scattering is small, 3 or 4. Since simulation data easily optimize the network model, the model can be flexibly applied to a wide variety of Compton telescopes.

preprint2022arXiv

NuSTAR Observations of 52 Compton-thick Active Galactic Nuclei Selected by the Swift/BAT All-sky Hard X-Ray Survey

We present the systematic broadband X-ray spectral analysis of 52 Compton-thick ($24 \leq \log N_{\mathrm{H}}^{\mathrm{LOS}}/\mathrm{cm}^{-2}$) active galactic nucleus (CTAGN) candidates selected by the Swift/BAT all-sky hard X-ray survey observed with Chandra, XMM--Newton, Swift/XRT, Suzaku, and NuSTAR. The XMM--Newton data of 10 objects and the NuSTAR data of 15 objects are published for the first time. We use an X-ray spectral model from a clumpy torus (XClumpy) to determine the torus properties. As a result, the hydrogen column density along the line of sight $N_{\mathrm{H}}^{\mathrm{LOS}}$ obtained from the XClumpy model indicate that 24 objects are Compton-thin AGNs and 28 objects are Compton-thick AGNs in the 90\% confidence interval. The main reason is the difference in the torus model applied. The hydrogen column density along the equatorial direction $N_{\mathrm{H}}^{\mathrm{Equ}}$ of CTAGNs inferred from the XClumpy model is larger than that of less obscured AGNs. The Compton-thin torus covering factor $C_{22}$ obtained from the XClumpy model is consistent with that of Ricci et al. (2017) in the low Eddington ratio ($\log R_{\mathrm{Edd}} \leq -1.0$), whereas $C_{22}$ inferred from the XClumpy model is larger than that of Ricci et al. (2017) in the high Eddington ratio ($-1.0 \leq \log R_{\mathrm{Edd}}$). The average value of the Compton-thick torus covering factor $C_{24}$ obtained from the XClumpy model is $36_{-4}^{+4}$\%. This value is larger than that of Ricci et al. (2015) ($C_{24} \simeq 27_{-4}^{+4}$\%) based on the assumption that all AGNs have intrinsically the same torus structure. These results suggest that the structure of CTAGN may be intrinsically different from that of less obscured AGN.

preprint2022arXiv

Reconstruction of multiple Compton scattering events in MeV gamma-ray Compton telescopes towards GRAMS: the physics-based probabilistic model

Aimed at progress in mega-electron volt (MeV) gamma-ray astronomy, which has not yet been well-explored, Compton telescope missions with a variety of detector concepts have been proposed so far. One of the key techniques for these future missions is an event reconstruction algorithm that is able to determine the scattering orders of multiple Compton scattering events and to identify events in which gamma rays escape from the detectors before they deposit all of their energies. We revisit previous event reconstruction methods and propose a modified algorithm based on a probabilistic method. First, we present a general formalism of the probabilistic model of Compton scattering describing physical interactions inside the detector and measurement processes. Then, we also introduce several approximations in the calculation of the probability functions for efficient computation. For validation, the developed algorithm has been applied to simulation data of a Compton telescope using a liquid argon time projection chamber, which is a new type of Compton telescope proposed for the GRAMS project. We have confirmed that it works successfully for up to 8-hit events, including correction of incoming gamma-ray energies for escape events. The proposed algorithm can be used for next-generation MeV gamma-ray missions featured by large-volume detectors, e.g., GRAMS.

preprint2021arXiv

A NuSTAR confirmation of the 36 ks hard X-ray pulse-phase modulation in the magnetar 1E 1547.0$-$5408

This paper describes an analysis of the NuSTAR data of the fastest-rotating magnetar 1E 1547$-$5408, acquired in 2016 April for a time lapse of 151 ks. The source was detected with a 1-60 keV flux of $1.7 \times 10^{-11}$ ergs s$^{-1}$ cm$^{-2}$, and its pulsation at a period of $2.086710(5)$ sec. In 8-25 keV, the pulses were phase-modulated with a period of $T=36.0 \pm 2.3$ ks, and an amplitude of $\sim 0.2$ sec. This reconfirms the Suzaku discovery of the same effect at $T=36.0 ^{+4.5}_{-2.5} $ ks, made in the 2009 outburst. These results strengthen the view derived from the Suzaku data, that this magnetar performs free precession as a result of its axial deformation by $\sim 0.6 \times 10^{-4}$, possibly caused by internal toroidal magnetic fields reaching $\sim 10^{16}$ G. Like in the Suzaku case, the modulation was not detected in energies below $\sim 8$ keV. Above 10 keV, the pulse-phase behaviour, including the 36 ks modulation parameters, exhibited complex energy dependences: at $\sim 22$ keV, the modulation amplitude increased to $\sim 0.5$ sec, and the modulation phase changed by $\sim 65^\circ$ over 10--27 keV, followed by a phase reversal. Although the pulse significance and pulsed fraction were originally very low in $>10$ keV, they both increased noticeably, when the arrival times of individual photons were corrected for these systematic pulse-phase variations. Possible origins of these complex phenomena are discussed, in terms of several physical processes that are specific to ultra-strong magnetic fields.

preprint2021arXiv

Simulation-based spectral analysis of X-ray CCD data affected by photon pile-up

We have developed a simulation-based method of spectral analysis for pile-up affected data of X-ray CCDs without any loss of photon statistics. As effects of the photon pile-up appear as complicated nonlinear detector responses, we employ a detailed simulation to calculate the important processes in an X-ray observation including physical interactions, detector signal generation, detector readout, and a series of data reduction processes. This simulation naturally reproduces X-ray-like and background-like events as results of X-ray photon merging in a single pixel or in a chunk of adjacent pixels, allowing us to construct a nonlinear spectral analysis framework that can treat pile-up affected observation data. For validation, we have performed data analysis of Suzaku XIS observations by using this framework with various parameters of the detector simulation all of which are optimized for that instrument. We present three cases of different pile-up degrees: PKS~2155-304 (negligible pile-up), Aquila~X-1 (moderate pile-up), and the Crab Nebula (strong pile-up); we show that the nonlinear analysis method produces results consistent with a conventional linear analysis for the negligible pile-up condition, and accurately corrects well-known pile-up effects such as spectral hardening and flux decrease for the pile-up cases. These corrected results are consistent with those obtained by a widely used core-exclusion method or by other observatories with much higher timing resolutions (without pile-up). Our framework is applicable to any types of CCDs used for X-ray astronomy including a future mission such as XRISM by appropriate optimization of the simulation parameters.

preprint2021arXiv

The simulation framework of the timing-based localization for future all-sky gamma-ray observations with a fleet of Cubesats

The timing-based localization, which utilize the triangulation principle with the different arrival time of gamma-ray photons, with a fleet of Cubesats is a unique and powerful solution for the future all-sky gamma-ray observation, which is a key for identification of the electromagnetic counterpart of the gravitational wave sources. The Cubesats Applied for MEasuring and Localising Transients (CAMELOT) mission is now being promoted by the Hungarian and Japanese collaboration with a basic concept of the nine Cubesats constellations in low earth orbit. The simulation framework for estimation of the localization capability has been developed including orbital parameters, an algorithm to estimate the expected observed profile of gamma-ray photons, finding the peak of the cross-correlation function, and a statistical method to find a best-fit position and its uncertainty. It is revealed that a degree-scale localization uncertainty can be achieved by the CAMELOT mission concept for bright short gamma-ray bursts, which could be covered by future large field of view ground-based telescopes. The new approach utilizing machine-learning approach is also investigated to make the procedure automated for the future large scale constellations. The trained neural network with 10$^6$ simulated light curves generated by the artificial short burst templates successfully predicts the time-delay of the real light curve and achieves a comparable performance to the cross-correlation algorithm with full automated procedures.

preprint2020arXiv

Application of X-Ray Clumpy Torus Model (XCLUMPY) to 10 Obscured Active Galactic Nuclei Observed with Suzaku and NuSTAR

We apply XCLUMPY, an X-ray spectral model from a clumpy torus in an active galactic nucleus (AGN), to the broadband X-ray spectra of 10 obscured AGNs observed with both Suzaku and NuSTAR. The infrared spectra of these AGNs were analyzed with the CLUMPY code. Since XCLUMPY adopts the same clump distribution as that in the CLUMPY, we can directly compare the torus parameters obtained from the X-ray spectra and those from the infrared ones. The torus angular widths determined from the infrared spectra ($σ_{\mathrm{IR}}$) are systematically larger than those from the X-ray data ($σ_{\mathrm{X}}$); the difference ($σ_{\mathrm{IR}}-σ_{\mathrm{X}}$) correlates with the inclination angle determined from the X-ray spectrum. These results can be explained by the contribution from dusty polar outflows to the observed infrared flux, which becomes more significant at higher inclinations (more edge-on views). The ratio of the hydrogen column density and V-band extinction in the line of sight absorber shows large scatter ($\simeq$1 dex) around the Galactic value, suggesting that a significant fraction of AGNs have dust-rich circumnuclear environments.

preprint2020arXiv

Artifact-less Coded Aperture Imaging in the X-ray Band with Multiple Different Random Patterns

The coded aperture imaging technique is a useful method of X-ray imaging in observational astrophysics. However, the presence of imaging noise or so-called artifacts in a decoded image is a drawback of this method. We propose a new coded aperture imaging method using multiple different random patterns for significantly reducing the image artifacts. This aperture mask contains multiple different patterns each of which generates a different artifact distribution in its decoded image. By summing all decoded images of the different patterns, the artifact distributions are cancelled out, and we obtain a remarkably accurate image. We demonstrate this concept with imaging experiments of a monochromatic 16 keV hard X-ray beam at the synchrotron photon facility SPring-8, using the combination of a CMOS image sensor and an aperture mask that has four different random patterns composed of holes with a diameter of 27 um and a separation of 39 um. The entire imaging system is installed in a 25 cm-long compact size, and achieves an angular resolution of < 30&#39;&#39; (full width at half maximum). In addition, we show by Monte Carlo simulation that the artifacts can be reduced more effectively if the number of different patterns increases to 8 or 16.

preprint2020arXiv

Development of the detector simulation framework for the Wideband Hybrid X-ray Imager onboard FORCE

FORCE is a Japan-US space-based astronomy mission for an X-ray imaging spectroscopy in an energy range of 1--80 keV. The Wideband Hybrid X-ray Imager (WHXI), which is the main focal plane detector, will use a hybrid semiconductor imager stack composed of silicon and cadmium telluride (CdTe). The silicon imager will be a certain type of the silicon-on-insulator (SOI) pixel sensor, named the X-ray pixel (XRPIX) series. Since the sensor has a small pixel size (30--36 $μ$m) and a thick sensitive region (300--500 $μ$m), understanding the detector response is not trivial and is important in order to optimize the camera design and to evaluate the scientific capabilities. We have developed a framework to simulate observations of celestial sources with semiconductor sensors. Our simulation framework was tested and validated by comparing our simulation results to laboratory measurements using the XRPIX 6H sensor. The simulator well reproduced the measurement results with reasonable physical parameters of the sensor including an electric field structure, a Coulomb repulsion effect on the carrier diffusion, and arrangement of the degraded regions. This framework is also applicable to future XRPIX updates including the one which will be part of the WHXI, as well as various types of semiconductor sensors.

preprint2020arXiv

Snowmass 2021 Letter of Interest: The GRAMS Project: MeV Gamma-Ray Observations and Antimatter-Based Dark Matter Searches

The Gamma-Ray and AntiMatter Survey (GRAMS) project is a proposed next-generation balloon/satellite mission targeting both MeV gamma-ray observations and antimatter-based dark matter searches. A cost-effective, large-scale Liquid Argon Time Projection Chamber (LArTPC) detector technology will allow GRAMS to have a significantly improved sensitivity to MeV gamma rays while extensively probing dark matter parameter space via antimatter measurements.

preprint2020arXiv

The thermal-radiative wind in low mass X-ray binary H 1743-322: II. iron line predictions from Monte Carlo radiation transfer

We show the best current simulations of the absorption and emission features predicted from thermal-radiative winds produced from X-ray illumination of the outer accretion disc in binary systems. We use the density and velocity structure derived from a radiation hydrodynamic code as input to a Monte-Carlo radiation transport calculation. The initial conditions are matched to those of the black hole binary system H1743-322 in its soft, disc dominated state, where wind features are seen in Chandra grating data. Our simulation fits well to the observed line profile, showing that these physical wind models can be the origin of the absorption features seen, rather than requiring a magnetically driven wind. We show how the velocity structure is the key observable discriminator between magnetic and thermal winds. Magnetic winds are faster at smaller radii, whereas thermal winds transition to a static atmosphere at smaller radii. New data from XRISM (due for launch Jan 2022) will give an unprecedented view of the physics of the wind launch and acceleration processes, but the existence of static atmospheres in small disc systems already rules out magnetic winds which assume self-similar magnetic fields from the entire disc as the origin of the absorption features seen.

preprint2020arXiv

The thermal-radiative wind in the neutron star low mass X-ray binary GX 13+1

We fit the observed high ionisation X-ray absorption lines in the neutron star binary GX13+1 with a full simulation of a thermal-radiative wind. This uses a radiation hydrodynamic code coupled to Monte Carlo radiation transfer to compute the observed line profiles from Hydrogen and Helium-like iron and Nickel, including all strong Kα and Kβ transitions. The wind is very strong as this object has a very large disc and is very luminous. The absorption lines from Fe Kα are strongly saturated as the ion columns are large, so the line equivalent widths (EWs) depend sensitively on the velocity structure. We additionally simulate the lines including isotropic turbulence at the level of the azimuthal and radial velocities. We fit these models to the Fe xxv and xxvi absorption lines seen in the highest resolution Chandra third order HETGS data. These data already rule out the addition of turbulence at the level of the radial velocity of ~500 km/s. The velocity structure predicted by the thermal-radiative wind alone is a fairly good match to the observed profile, with an upper limit to additional turbulence at the level of the azimuthal velocity of ~100 km/s. This gives stringent constraints on any remaining contribution from magnetic acceleration.

preprint2020arXiv

Uniform distribution of the extremely overionized plasma associated with the supernova remnant G359.1-0.5

We report on the results of our detailed analyses on the peculiar recombining plasma of the supernova remnant (SNR) G359.1$-$0.5, and the interacting CO clouds. Combining {\it Chandra} and {\it Suzaku} data, we estimated the ionization state of the plasma with a careful treatment of the background spectrum. The average spectrum showed a remarkably large deviation of the electron temperature ($\sim$0.17 keV) from the initial temperature ($>$ 16 keV), indicating that the plasma is in a highly recombination-dominant state. On the other hand, the recombination timescale $({\it n_{\rm e} t})$ is comparable to those of the other recombining SNRs ($\sim4.2 \times 10^{11}$ cm$^{-3}$ s). We also searched for spatial variation of the plasma parameters, but found no significant differences. Using $^{12}$CO($J$=2--1) data obtained with NANTEN2, we found a new, plausible candidate for the interacting CO cloud, which has a line-of-sight velocity of $\sim -$20 km s$^{-1}$. This indicates that the SNR is located at a distance of $\sim$4 kpc, which is the foreground of the Galactic center, as previously reported. The associated CO cloud does not show clear spatial coincidence with the nearby GeV/TeV emission, indicating that the origins of the GeV/TeV emission are likely unrelated to G359.1$-$0.5.

preprint2019arXiv

X-Ray Spectral Model from Clumpy Torus and Its Application to Circinus Galaxy

We construct an X-ray spectral model from the clumpy torus in an active galactic nucleus (AGN), designated as &#34;XCLUMPY&#34;, utilizing the Monte Carlo simulation for Astrophysics and Cosmology framework (MONACO: Odaka et al. 2011, 2016). The adopted geometry of the torus is the same as that in Nenkova et al. (2008), who assume a power law distribution of clumps in the radial direction and a normal distribution in the elevation direction. We investigate the dependence of the X-ray continuum and Fe K$α$ fluorescence line profile on the torus parameters. Our model is compared with other torus models: MYTorus model (Murphy & Yaqoob 2009), Ikeda model (Ikeda et al. 2009), and CTorus model (Liu & Li 2014). As an example, we also present the results applied to the broadband X-ray spectra of the Circinus galaxy observed with XMM-Newton, Suzaku, and NuSTAR. Our model can well reproduce the data, yielding a hydrogen column density along the equatorial plane $N_{\mathrm{H}}^{\mathrm{Equ}} = 9.08_{-0.08}^{+0.14} \times 10^{24}$ cm$^{-2}$, a torus angular width $σ= 14.7_{-0.39}^{+0.44}$ degree, and a 2--10 keV luminosity $\log L_{2-10}/\mathrm{erg \ s^{-1}} = 42.8$. These results are discussed in comparison with the observations in other wavelengths.