Showing 1–50 of 96 results for author: Shibata, S

The centrifugal acceleration and the Y-point of the Pulsar Magnetosphere

Authors: Shinpei Shibata, Shota Kisaka

Abstract: We investigate the centrifugal acceleration in an axisymmetric pulsar magnetosphere under the ideal-MHD approximation. We solved the field-aligned equations of motion for flows inside the current sheet with finite thickness. We find that flows coming into the vicinity of a Y-point become super fast. The centrifugal acceleration takes place efficiently, and most of the Poynting energy is converted… ▽ More We investigate the centrifugal acceleration in an axisymmetric pulsar magnetosphere under the ideal-MHD approximation. We solved the field-aligned equations of motion for flows inside the current sheet with finite thickness. We find that flows coming into the vicinity of a Y-point become super fast. The centrifugal acceleration takes place efficiently, and most of the Poynting energy is converted into kinetic energy. However, the super fast flow does not provide enough centrifugal drift current to open the magnetic field. Opening of the magnetic field is possible by the plasmas that are accelerated in the azimuthal direction with a large Lorentz factor in the closed field region. We find that this acceleration takes place if the field strength increases toward the Y-point from inside. The accelerated plasma is transferred from the closed field region to the open field region by magnetic reconnection with plasmoid emission. We also estimate the Lorentz factor to be reached in the centrifugal wind. △ Less

Submitted 26 July, 2024; originally announced July 2024.

Comments: 10 pages, 13figures, Accepted for pulication in The Astrophysical Journal

Fate of a remnant solid disk around an eccentric giant planet

Authors: S. Shibata, R. Helled

Abstract: The composition of giant planets' atmospheres is an important tracer of their formation history. While many theoretical studies investigate the heavy-element accretion within a gaseous protoplanetary disk, the possibility of solid accretion after disk dissipation has not been explored. Here, we focus on the case of a gas giant planet excited to an eccentric orbit and assess the likelihood of solid… ▽ More The composition of giant planets' atmospheres is an important tracer of their formation history. While many theoretical studies investigate the heavy-element accretion within a gaseous protoplanetary disk, the possibility of solid accretion after disk dissipation has not been explored. Here, we focus on the case of a gas giant planet excited to an eccentric orbit and assess the likelihood of solid accretion after disk dissipation. We perform N-body simulations of planetesimals and embryos around an eccentric giant planet. We consider various sizes and orbits for the eccentric planet and determine the fate of planetesimals and embryos. We find that the orbital evolution of solids, such as planetesimals and embryos, is regulated by weak encounters with the eccentric planet rather than strong close encounters. Even in the region where the Safronov number is smaller than unity, most solid materials fall onto the central star or are ejected from the planetary system. We also develop an analytical model of the solid accretion along the orbital evolution of a giant planet, where the accretion probability is obtained as a function of the planetary mass, radius, semi-major axis, eccentricity, inclination, and solid disk thickness. Our model predicts that $\sim$0.01-0.1 $M_\oplus$ of solids is accreted onto an eccentric planet orbiting in the outer disk ($\sim10$ au). The accreted heavy-element mass increases (decreases) with the eccentricity (inclination) of the planet. We also discuss the possibility of collisions of terrestrial planets and find that $\sim10\%$ of the hot Jupiters formed via high-eccentric migration collide with a planet of $10M_\oplus$. However, we find that solid accretion and collisions with terrestrial planets are minor events for planets in the inner orbit, and a different accretion process is required to enrich eccentric giant planets with heavy elements. △ Less

Submitted 16 July, 2024; originally announced July 2024.

Comments: 15 pages, 16 Figures, accepted to Astronomy and Astrophysics

Small-scale cosmic ray anisotropy observed by the GRAPES-3 experiment at TeV energies

Authors: M. Chakraborty, S. Ahmad, A. Chandra, S. R. Dugad, U. D. Goswami, S. K. Gupta, B. Hariharan, Y. Hayashi, P. Jagadeesan, A. Jain, P. Jain, S. Kawakami, T. Koi, H. Kojima, S. Mahapatra, P. K. Mohanty, R. Moharana, Y. Muraki, P. K. Nayak, T. Nonaka, T. Nakamura, A. Oshima, B. P. Pant, D. Pattanaik, S. Paul , et al. (13 additional authors not shown)

Abstract: GRAPES-3 is a mid-altitude (2200 m) and near equatorial ($11.4^{\circ}$ North) air shower array, overlapping in its field of view for cosmic ray observations with experiments that are located in Northern and Southern hemispheres. We analyze a sample of $3.7\times10^9$ cosmic ray events collected by the GRAPES-3 experiment between 1 January 2013 and 31 December 2016 with a median energy of… ▽ More GRAPES-3 is a mid-altitude (2200 m) and near equatorial ($11.4^{\circ}$ North) air shower array, overlapping in its field of view for cosmic ray observations with experiments that are located in Northern and Southern hemispheres. We analyze a sample of $3.7\times10^9$ cosmic ray events collected by the GRAPES-3 experiment between 1 January 2013 and 31 December 2016 with a median energy of $\sim16$ TeV for study of small-scale ($<60^{\circ}$) angular scale anisotropies. We observed two structures labeled as A and B, deviate from the expected isotropic distribution of cosmic rays in a statistically significant manner. Structure `A' spans $50^{\circ}$ to $80^{\circ}$ in the right ascension and $-15^{\circ}$ to $30^{\circ}$ in the declination coordinate. The relative excess observed in the structure A is at the level of $(6.5\pm1.3)\times10^{-4}$ with a statistical significance of 6.8 standard deviations. Structure `B' is observed in the right ascension range of $110^{\circ}$ to $140^{\circ}$. The relative excess observed in this region is at the level of $(4.9\pm1.4)\times10^{-4}$ with a statistical significance of 4.7 standard deviations. These structures are consistent with those reported by Milagro, ARGO-YBJ, and HAWC. These observations could provide a better understanding of the cosmic ray sources, propagation and the magnetic structures in our Galaxy. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Can Uranus and Neptune form concurrently via pebble, gas and planetesimal accretion?

Authors: Linn E. J. Eriksson, Marit A. S. Mol Lous, Sho Shibata, Ravit Helled

Abstract: The origin of Uranus and Neptune has long been challenging to explain, due to the large orbital distances from the Sun. After a planetary embryo has been formed, the main accretion processes are likely pebble, gas and planetesimal accretion. Previous studies of Uranus and Neptune formation typically don't consider all three processes; and furthermore, do not investigate how the formation of the ou… ▽ More The origin of Uranus and Neptune has long been challenging to explain, due to the large orbital distances from the Sun. After a planetary embryo has been formed, the main accretion processes are likely pebble, gas and planetesimal accretion. Previous studies of Uranus and Neptune formation typically don't consider all three processes; and furthermore, do not investigate how the formation of the outer planet impacts the inner planet. In this paper we study the concurrent formation of Uranus and Neptune via both pebble, gas and planetesimal accretion. We use a dust-evolution model to predict the size and mass flux of pebbles, and derive our own fit for gas accretion. We do not include migration, but consider a wide range of formation locations between 12 and 40au. If the planetary embryos form at the same time and with the same mass, our formation model with an evolving dust population is unable to produce Uranus and Neptune analogues. This is because the mass difference between the planets and the H-He mass fractions become too high. However, if the outer planetary embryo forms earlier and/or more massive than the inner embryo, the two planets do form in a few instances when the disk is metal-rich and dissipates after a few Myr. Furthermore, our study suggests that in-situ formation is rather unlikely. Nethertheless, giant impacts and/or migration could potentially aid in the formation, and future studies including these processes could bring us one step closer to understanding how Uranus and Neptune formed. △ Less

Submitted 29 September, 2023; originally announced October 2023.

Comments: Accepted for publication in MNRAS

Magnetic Field Structure of the Crab Pulsar Wind Nebula Revealed with IXPE

Authors: T. Mizuno, H. Ohno, E. Watanabe, N. Bucciantini, S. Gunji, S. Shibata, P. Slane, W. C. Weisskopf

Abstract: We report a detailed study of the magnetic-field structure of the Crab pulsar wind nebula, using the X-ray polarization data in 2--8~keV obtained with the Imaging X-ray Polarimetry Explorer. Contamination of the pulsar emission to the data of the nebula region was removed through application of a stringent pulsation phase-cut, extracting a phase range of 0.7--1.0 only. We found that the electric f… ▽ More We report a detailed study of the magnetic-field structure of the Crab pulsar wind nebula, using the X-ray polarization data in 2--8~keV obtained with the Imaging X-ray Polarimetry Explorer. Contamination of the pulsar emission to the data of the nebula region was removed through application of a stringent pulsation phase-cut, extracting a phase range of 0.7--1.0 only. We found that the electric field vector polarization angle (PA) was about $130^{\circ}$ from north to east with the polarization degree (PD) of about 25\% at the pulsar position, indicating that the direction of the toroidal magnetic field is perpendicular to the pulsar spin axis in the region close to the termination shock. The PA gradually deviated from the angle as an increasing function of the distance from the pulsar. There was a region of a low PD to the west of the X-ray torus. Although such a region is expected to be located at the torus edge, where geometrical depolarization due to a steep spatial variation of the PA is expected, the observed low-PD region positionally deviated from the edge. We found that the region of low PD positionally coincided with a dense filament seen in the optical band, and conjecture that the low-PD region may be produced through deflection of the pulsar wind. By comparing the values of the PD at the pulsar position between the data and a model, in which toroidal and turbulent magnetic fields were considered, we estimated the fractional energy of the turbulent magnetic field to be about $2/3$ of the total. We also evaluated a potential polarization of the northern jet in the nebula and derived the PD and PA to be about $30\%$ and $120^{\circ}$, respectively. △ Less

Submitted 28 September, 2023; originally announced September 2023.

Comments: This is a pre-copyedited, author-produced version of an article accepted for publication in PASJ following peer review

Heavy-element Accretion by Proto-Jupiter in a Massive Planetesimal Disk, Revisited

Authors: Sho Shibata, Ravit Helled, Hiroshi Kobayashi

Abstract: Planetesimal accretion is a key source for heavy-element enrichment in giant planets. It has been suggested that Jupiter's enriched envelope is a result of planetesimal accretion during its growth assuming it formed in a massive planetesimal disk. In this study, we simulate Jupiter's formation in this scenario. We assume in-situ formation and perform N-body simulations to infer the solid accretion… ▽ More Planetesimal accretion is a key source for heavy-element enrichment in giant planets. It has been suggested that Jupiter's enriched envelope is a result of planetesimal accretion during its growth assuming it formed in a massive planetesimal disk. In this study, we simulate Jupiter's formation in this scenario. We assume in-situ formation and perform N-body simulations to infer the solid accretion rate. We find that tens-Earth masses of planetesimals can be captured by proto-Jupiter during the rapid gas accretion phase. However, if several embryos are formed near Jupiter's core, which is an expected outcome in the case of a massive planetesimal disk, scattering from the embryos increases the eccentricity and inclination of planetesimals and therefore significantly reduces the accretion efficiency. We also compare our results with published semi-analytical models and show that these models cannot reproduce the N-body simulations especially when the planetesimal disk has a large eccentricity and inclination. We show that when the dynamical evolution of planetesimals is carefully modelled, the total mass of captured planetesimals $M_\mathrm{cap,tot}$ is $2 M_\oplus \lesssim M_\mathrm{cap,tot}\lesssim 18 M_\oplus$. The metallicity of Jupiter's envelope can be explained by the planetesimal accretion in our massive disk model despite the low accretion efficiency coming from the high eccentricity and inclination of planetesimals. Our study demonstrates the importance of detailed modelling of planetesimal accretion during the planetary growth and its implications to the heavy-element mass in gaseous planets. △ Less

Submitted 8 December, 2022; originally announced December 2022.

Comments: 20 pages, 17 figures, accepted for publication in MNRAS

Exploring the Ability of HST WFC3 G141 to Uncover Trends in Populations of Exoplanet Atmospheres Through a Homogeneous Transmission Survey of 70 Gaseous Planets

Authors: Billy Edwards, Quentin Changeat, Angelos Tsiaras, Kai Hou Yip, Ahmed F. Al-Refaie, Lara Anisman, Michelle F. Bieger, Amelie Gressier, Sho Shibata, Nour Skaf, Jeroen Bouwman, James Y-K. Cho, Masahiro Ikoma, Olivia Venot, Ingo Waldmann, Pierre-Olivier Lagage, Giovanna Tinetti

Abstract: We present the analysis of the atmospheres of 70 gaseous extrasolar planets via transit spectroscopy with Hubble's Wide Field Camera 3 (WFC3). For over half of these, we statistically detect spectral modulation which our retrievals attribute to molecular species. Among these, we use Bayesian Hierarchical Modelling to search for chemical trends with bulk parameters. We use the extracted water abund… ▽ More We present the analysis of the atmospheres of 70 gaseous extrasolar planets via transit spectroscopy with Hubble's Wide Field Camera 3 (WFC3). For over half of these, we statistically detect spectral modulation which our retrievals attribute to molecular species. Among these, we use Bayesian Hierarchical Modelling to search for chemical trends with bulk parameters. We use the extracted water abundance to infer the atmospheric metallicity and compare it to the planet's mass. We also run chemical equilibrium retrievals, fitting for the atmospheric metallicity directly. However, although previous studies have found evidence of a mass-metallicity trend, we find no such relation within our data. For the hotter planets within our sample, we find evidence for thermal dissociation of dihydrogen and water via the H$^-$ opacity. We suggest that the general lack of trends seen across this population study could be due to i) the insufficient spectral coverage offered by HST WFC3 G141, ii) the lack of a simple trend across the whole population, iii) the essentially random nature of the target selection for this study or iv) a combination of all the above. We set out how we can learn from this vast dataset going forward in an attempt to ensure comparative planetology can be undertaken in the future with facilities such as JWST, Twinkle and Ariel. We conclude that a wider simultaneous spectral coverage is required as well as a more structured approach to target selection. △ Less

Submitted 1 November, 2022; originally announced November 2022.

Comments: Accepted for publication in ApJS

Constraining the origin of giant exoplanets via elemental abundance measurements

Authors: Henrik Knierim, Sho Shibata, Ravit Helled

Abstract: The origin of close-in giant planets is a key open question in planet formation theory. The two leading models are (i) formation at the outer disk followed by migration and (ii) in situ formation. In this work we determine the atmospheric composition of warm Jupiters for both formation scenarios. We perform N-body simulations of planetesimal accretion interior and exterior to the water ice-line fo… ▽ More The origin of close-in giant planets is a key open question in planet formation theory. The two leading models are (i) formation at the outer disk followed by migration and (ii) in situ formation. In this work we determine the atmospheric composition of warm Jupiters for both formation scenarios. We perform N-body simulations of planetesimal accretion interior and exterior to the water ice-line for various planetary formation locations, planetary masses, and planetesimal sizes to estimate the accreted heavy-element mass and final planetary composition. We find that the two models differ significantly: migrating giant planets have 2-14 times higher metallicities than planets that form in situ. The ratio between refractories and volatiles is found to be above one for migrating planets but below 0.4 for planets that form in situ. We also identify very different trends between heavy-element enrichment and planetary mass for these two formation mechanisms. While the metallicity of migrating planets is found to increase with decreasing planetary mass, it is about constant for in situ formation. Our study highlights the importance of measuring the atmospheric composition of warm Jupiters and its connection to their formation and evolutionary paths. △ Less

Submitted 6 October, 2022; v1 submitted 2 September, 2022; originally announced September 2022.

Comments: published in A&A, 5 pages, including 4 figures, plus 2 pages appendix; v2 follows the A&A language guidelines

Journal ref: A&A 665, L5 (2022)

Spectral Break of Energetic Pulsar Wind Nebulae Detected with Wideband X-ray Observations

Authors: Aya Bamba, Shinpei Shibata, Shuta J. Tanaka, Koji Mori, Hiroyuki Uchida, Yukikatsu Terada, Wataru Ishizaki

Abstract: Pulsar wind nebulae (PWNe) are one of the most energetic galactic sources with bright emissions from radio waves to very high-energy gamma-rays. We perform wideband X-ray spectroscopy of four energetic PWNe, N157B, PSR J1813-1749, PSR J1400-6325, and G21.5-0.9, with the Suzaku, Chandra, NuSTAR, and Hitomi observatories. A significant spectral break or cutoff feature is found in the hard X-ray band… ▽ More Pulsar wind nebulae (PWNe) are one of the most energetic galactic sources with bright emissions from radio waves to very high-energy gamma-rays. We perform wideband X-ray spectroscopy of four energetic PWNe, N157B, PSR J1813-1749, PSR J1400-6325, and G21.5-0.9, with the Suzaku, Chandra, NuSTAR, and Hitomi observatories. A significant spectral break or cutoff feature is found in the hard X-ray band for all the samples, except for N157B. The break energies in the broken power-law fitting are in the range of 4--14 keV, whereas the cutoff energies in the cutoff power-law fitting are at 22 keV or higher. The break or cutoff energy does not show a significant correlation with either the spin-down energy or characteristic age of the hosting pulsars. A possible correlation is found between the photon index change in the broken power-law fitting and the X-ray emitting efficiency of the pulsars, although its significance is not high enough to be conclusive. We discuss what determines the break parameters based on simple models. △ Less

Submitted 22 July, 2022; originally announced July 2022.

Comments: 21 pages, 22 figures, PASJ, accepted

Proton Penetration Efficiency over a High Altitude Observatory in Mexico

Authors: S. Miyake, T. Koi, Y. Muraki, Y. Matsubara, S. Masuda, P. Miranda, T. Naito, E. Ortiz, A. Oshima, T. Sakai, T. Sako, S. Shibata, H. Takamaru, M. Tokumaru, J. F. Valdes-Galicia

Abstract: In association with a large solar flare on November 7, 2004, the solar neutron detectors located at Mt. Chacaltaya (5,250m) in Bolivia and Mt. Sierra Negra (4,600m) in Mexico recorded very interesting events. In order to explain these events, we have performed a calculation solving the equation of motion of anti-protons inside the magnetosphere. Based on these results, the Mt. Chacaltaya event may… ▽ More In association with a large solar flare on November 7, 2004, the solar neutron detectors located at Mt. Chacaltaya (5,250m) in Bolivia and Mt. Sierra Negra (4,600m) in Mexico recorded very interesting events. In order to explain these events, we have performed a calculation solving the equation of motion of anti-protons inside the magnetosphere. Based on these results, the Mt. Chacaltaya event may be explained by the detection of solar neutrons, while the Mt. Sierra Negra event may be explained by the first detection of very high energy solar neutron decay protons (SNDPs) around 6 GeV. △ Less

Submitted 8 July, 2022; v1 submitted 5 July, 2022; originally announced July 2022.

Comments: Paper presented in the 21st International Symposium on Very High Energy Cosmic Ray Intercations (ISVHE-CRI 2022) by online

Five key exoplanet questions answered via the analysis of 25 hot Jupiter atmospheres in eclipse

Authors: Quentin Changeat, Billy Edwards, Ahmed F. Al-Refaie, Angelos Tsiaras, Jack W. Skinner, James Y-K Cho, Kai H. Yip, Lara Anisman, Masahiro Ikoma, Michelle F. Bieger, Olivia Venot, Sho Shibata, Ingo P. Waldmann, Giovanna Tinetti

Abstract: Population studies of exoplanets are key to unlocking their statistical properties. So far the inferred properties have been mostly limited to planetary, orbital and stellar parameters extracted from, e.g., Kepler, radial velocity, and GAIA data. More recently an increasing number of exoplanet atmospheres have been observed in detail from space and the ground. Generally, however, these atmospheric… ▽ More Population studies of exoplanets are key to unlocking their statistical properties. So far the inferred properties have been mostly limited to planetary, orbital and stellar parameters extracted from, e.g., Kepler, radial velocity, and GAIA data. More recently an increasing number of exoplanet atmospheres have been observed in detail from space and the ground. Generally, however, these atmospheric studies have focused on individual planets, with the exception of a couple of works which have detected the presence of water vapor and clouds in populations of gaseous planets via transmission spectroscopy. Here, using a suite of retrieval tools, we analyse spectroscopic and photometric data of 25 hot Jupiters, obtained with the Hubble and Spitzer Space Telescopes via the eclipse technique. By applying the tools uniformly across the entire set of 25 planets, we extract robust trends in the thermal structure and chemical properties of hot Jupiters not obtained in past studies. With the recent launch of JWST and the upcoming missions Twinkle, and Ariel, population based studies of exoplanet atmospheres, such as the one presented here, will be a key approach to understanding planet characteristics, formation, and evolution in our galaxy. △ Less

Submitted 13 May, 2022; v1 submitted 25 April, 2022; originally announced April 2022.

Comments: 66 pages, 23 figures, 7 tables. Published in The Astrophysical Journal Supplement Series

Journal ref: ApJS 260 3 (2022)

Enrichment of Jupiter's atmosphere by late planetesimal bombardment

Authors: Sho Shibata, Ravit Helled

Abstract: Jupiter's atmosphere is enriched with heavy elements by a factor of about 3 compared to proto-solar. The origin of this enrichment and whether it represent the bulk composition of the planetary envelope remain unknown. Internal structure models of Jupiter suggest that its envelope is separated from the deep interior and that the planet is not fully mixed. This implies that Jupiter's atmosphere was… ▽ More Jupiter's atmosphere is enriched with heavy elements by a factor of about 3 compared to proto-solar. The origin of this enrichment and whether it represent the bulk composition of the planetary envelope remain unknown. Internal structure models of Jupiter suggest that its envelope is separated from the deep interior and that the planet is not fully mixed. This implies that Jupiter's atmosphere was enriched with heavy elements just before the end of its formation. Such enrichment can be a result of late planetesimal accretion. However, in-situ Jupiter formation models suggest the decreasing accretion rate with increasing planetary mass, which cannot explain Jupiter's atmospheric enrichment. In this study, we model Jupiter's formation and show that an migration of proto-Jupiter from $\sim$ 20 AU to its current location can lead to a late planetesimal accretion and atmospheric enrichment. Late planetesimal accretion does not occur if proto-Jupiter migrates only a few AU. We suggest that if Jupiter's outermost layer is fully-mixed and is relatively thin (up to $\sim$ 20\% of its mass), such late accretion can explain its measured atmospheric composition. It is therefore possible that Jupiter underwent significant orbital migration followed by late planetesimal accretion. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: 18 pages, 8 figures, Accepted for ApJL

The origin of the high metallicity of close-in giant exoplanets II The nature of the sweet spot for accretion

Authors: Sho Shibata, Ravit Helled, Masahiro Ikoma

Abstract: The composition of giant planets reflects their formation history. Planetesimal accretion during the phase of planetary migration could lead to the delivery of heavy elements into giant planets. In our previous paper (Shibata et al. 2020) we showed that planetesimal accretion during planetary migration occurs in a rather narrow region of the protoplanetary disk, which we refer as "the sweet spot f… ▽ More The composition of giant planets reflects their formation history. Planetesimal accretion during the phase of planetary migration could lead to the delivery of heavy elements into giant planets. In our previous paper (Shibata et al. 2020) we showed that planetesimal accretion during planetary migration occurs in a rather narrow region of the protoplanetary disk, which we refer as "the sweet spot for accretion". The goal of this paper is to reveal the nature of the sweet spot and investigate the role of the sweet spot in determining the composition of gas giant planets. We analytically derive the required conditions for the sweet spot. Then, we compare the derived equations with the numerical simulations. We find that the conditions required for the sweet spot can be expressed by the ratio of the gas damping timescale of the planetesimal orbits and the planetary migration timescale. If the planetary migration timescale depends on the surface density of disk gas inversely, the location of the sweet spot does not change with the disk evolution. The mass of planetesimals accreted by the planet depends on the amount of planetesimals that are shepherded by mean motion resonances. Our analysis suggests that tens Earth-mass of planetesimals can be shepherded into the sweet spot without planetesimal collisions. However, as more planetesimals are trapped into mean motion resonances, collisional cascade can lead to fragmentation of planetesimals. This could affect the location of the sweet spot and the population of small objects in planetary systems. We conclude that the composition of gas giant planets depends on whether the planets crossed the sweet spot during their formation. Constraining the metallicity of cold giant planets, that are expected to be outer than the sweet spot, would reveal key information for understanding the origin of heavy elements in giant planets. △ Less

Submitted 23 December, 2021; originally announced December 2021.

Comments: 17 pages, 13 figures, A&A accepted

Journal ref: A&A 659, A28 (2022)

Detectability of southern gamma-ray sources beyond 100 TeV with ALPAQUITA, the prototype experiment of ALPACA

Authors: S. Kato, C. A. H. Condori, E. dela Fuente, A. Gomi, K. Hibino, N. Hotta, I. Toledano-Juarez, Y. Katayose, C. Kato, K. Kawata, W. Kihara, Y. Ko, T. Koi, H. Kojima, D. Kurashige, J. Lozoya, F. Orozco-Luna, R. Mayta, P. Miranda, K. Munakata, H. Nakada, Y. Nakamura, Y. Nakazawa, C. Nina, M. Nishizawa , et al. (25 additional authors not shown)

Abstract: Andes Large-area PArticle detector for Cosmic-ray physics and Astronomy (ALPACA) is an international experiment that applies southern very-high-energy (VHE) gamma-ray astronomy to determine the origin of cosmic rays around the knee energy region ($10^{15}\, {\rm eV} - 10^{16}\, {\rm eV}$). The experiment consists of an air shower (AS) array with a surface of $83, 000\, {\rm m}^2$ and an undergroun… ▽ More Andes Large-area PArticle detector for Cosmic-ray physics and Astronomy (ALPACA) is an international experiment that applies southern very-high-energy (VHE) gamma-ray astronomy to determine the origin of cosmic rays around the knee energy region ($10^{15}\, {\rm eV} - 10^{16}\, {\rm eV}$). The experiment consists of an air shower (AS) array with a surface of $83, 000\, {\rm m}^2$ and an underground water Cherenkov muon detector (MD) array covering $5, 400\, {\rm m}^2$. The experimental site is at the Mt. Chacaltaya plateau in La Paz, Bolivia, with an altitude of $4, 740\, {\rm m}$ corresponding to $572\, {\rm g}/{\rm cm}^2$ atmospheric thickness. As the prototype experiment of ALPACA, the ALPAQUITA experiment aims to begin data acquisition in late 2021. The ALPAQUITA array consists of a smaller AS array ($18, 450\, {\rm m}^2$) and underground MD ($900\, {\rm m}^2$), which are now under construction. ALPAQUITA's sensitivity to gamma-ray sources is evaluated with Monte Carlo simulations. The simulation finds that five gamma-ray sources observed by H.E.S.S. and HAWC experiments will be detected by ALPAQUITA beyond 10 TeV and one out of these five - HESS J1702-420A - above 300 TeV in one calendar year observation. The latter finding means that scientific discussions can be made on the emission mechanism of gamma rays beyond 100 TeV from southern sources on the basis of the observational results of this prototype experiment. △ Less

Submitted 28 September, 2021; originally announced September 2021.

Comments: This is the manuscript accepted by Experimental Astronomy. For the published article, see https://doi.org/10.1007/s10686-021-09796-8

Journal ref: Exp Astron (2021)

Exploring the link between star and planet formation with Ariel

Authors: Diego Turrini, Claudio Codella, Camilla Danielski, Davide Fedele, Sergio Fonte, Antonio Garufi, Mario Giuseppe Guarcello, Ravit Helled, Masahiro Ikoma, Mihkel Kama, Tadahiro Kimura, J. M. Diederik Kruijssen, Jesus Maldonado, Yamila Miguel, Sergio Molinari, Athanasia Nikolaou, Fabrizio Oliva, Olja Panic, Marco Pignatari, Linda Podio, Hans Rickman, Eugenio Schisano, Sho Shibata, Allona Vazan, Paulina Wolkenberg

Abstract: The goal of the Ariel space mission is to observe a large and diversified population of transiting planets around a range of host star types to collect information on their atmospheric composition. The planetary bulk and atmospheric compositions bear the marks of the way the planets formed: Ariel's observations will therefore provide an unprecedented wealth of data to advance our understanding of… ▽ More The goal of the Ariel space mission is to observe a large and diversified population of transiting planets around a range of host star types to collect information on their atmospheric composition. The planetary bulk and atmospheric compositions bear the marks of the way the planets formed: Ariel's observations will therefore provide an unprecedented wealth of data to advance our understanding of planet formation in our Galaxy. A number of environmental and evolutionary factors, however, can affect the final atmospheric composition. Here we provide a concise overview of which factors and effects of the star and planet formation processes can shape the atmospheric compositions that will be observed by Ariel, and highlight how Ariel's characteristics make this mission optimally suited to address this very complex problem. △ Less

Submitted 3 September, 2021; v1 submitted 26 August, 2021; originally announced August 2021.

Comments: 67 pages, 11 figures, 1 table. Postprint of the article to appear as open access in the Ariel Special Issue by the journal Experimental Astronomy

On the Angular Momentum Extraction from the Rotation Powered Pulsars

Authors: Shinpei Shibata, Shota Kisaka

Abstract: The rotation powered pulsar loses angular momentum at a rate of the rotation power divided by the angular velocity $Ω_*$. This means that the length of the lever arm of the angular momentum extracted by the photons, relativistic particles and wind must be on average $c/Ω_*$, which is known as the light cylinder radius. Therefore, any deposition of the rotation power within the light cylinder cause… ▽ More The rotation powered pulsar loses angular momentum at a rate of the rotation power divided by the angular velocity $Ω_*$. This means that the length of the lever arm of the angular momentum extracted by the photons, relativistic particles and wind must be on average $c/Ω_*$, which is known as the light cylinder radius. Therefore, any deposition of the rotation power within the light cylinder causes insufficient loss of angular momentum. In this paper, we investigate two cases of this type of energy release: polar cap acceleration and Ohmic heating in the magnetospheric current inside the star. As for the first case, the outer magnetosphere beyond the light cylinder is found to compensate the insufficient loss of the angular momentum. We argue that the energy flux coming from the sub-rotating magnetic field lines must be larger than the solid-angle average value, and as a result, an enhanced energy flux emanating beyond the light cylinder is observed in different phases in the light curve from those of emission inside the light cylinder. As for the second case, the stellar surface rotates more slowly than the stellar interior. We find that the way the magnetospheric current closes inside the star is linked to how the angular momentum is transferred inside the star. We obtain numerical solutions which show that the magnetospheric current inside the star spreads over the polar cap magnetic flux embedded in the star in such a way that electromotive force is gained efficiently. △ Less

Submitted 29 July, 2021; originally announced July 2021.

Comments: Accepted for publication in MNRAS, 9 pages, 16 figures

Quantitative Age Estimation of Supernova Remnants and Associated Pulsars

Authors: Hiromasa Suzuki, Aya Bamba, Shinpei Shibata

Abstract: The age of a supernova remnant (SNR) is, though undoubtedly one of the most important properties for study of its evolution, difficult to estimate reliably in most cases. In this study, we compare the dynamical and plasma ages of the SNRs and characteristic ages of their associated pulsars with the corresponding SNRs' ages that are generally thought to be reliable ($t_{\rm r}$): historical and lig… ▽ More The age of a supernova remnant (SNR) is, though undoubtedly one of the most important properties for study of its evolution, difficult to estimate reliably in most cases. In this study, we compare the dynamical and plasma ages of the SNRs and characteristic ages of their associated pulsars with the corresponding SNRs' ages that are generally thought to be reliable ($t_{\rm r}$): historical and light-echo ages of the SNRs, kinematic ages of the ejecta knots and kinematic ages of the associated neutron stars (NS). The kinematic age of ejecta knots or a NS is the time that they have taken to reach the current positions from the explosion center. We use all of the available 24 systems for which $t_{\rm r}$ is already available (historical, light-echo, and ejecta kinematic ages) or measurable (NS kinematic age). We estimate the NS kinematic ages for eight SNR-NS systems by determining quantitatively the geometric centers of the SNR shells. The obtained $t_{\rm r}$ ranges from 33 yr to $\approx 400$ kyr. We find that the two SNR ages, dynamical and plasma ages, are consistent with $t_{\rm r}$ within a factor of four, whereas the characteristic ages of the pulsars differ from $t_{\rm r}$ by more than a factor of four in some systems. Using the $t_{\rm r}$ summarized in this work, we present the initial spin periods of the associated pulsars, which are more strictly constrained than the previous works, as well. △ Less

Submitted 21 April, 2021; v1 submitted 20 April, 2021; originally announced April 2021.

Comments: 14 pages, 5 figures, accepted for publication in ApJ

Ariel: Enabling planetary science across light-years

Authors: Giovanna Tinetti, Paul Eccleston, Carole Haswell, Pierre-Olivier Lagage, Jérémy Leconte, Theresa Lüftinger, Giusi Micela, Michel Min, Göran Pilbratt, Ludovic Puig, Mark Swain, Leonardo Testi, Diego Turrini, Bart Vandenbussche, Maria Rosa Zapatero Osorio, Anna Aret, Jean-Philippe Beaulieu, Lars Buchhave, Martin Ferus, Matt Griffin, Manuel Guedel, Paul Hartogh, Pedro Machado, Giuseppe Malaguti, Enric Pallé , et al. (293 additional authors not shown)

Abstract: Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths.… ▽ More Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths. It is the first mission dedicated to measuring the chemical composition and thermal structures of hundreds of transiting exoplanets, enabling planetary science far beyond the boundaries of the Solar System. The payload consists of an off-axis Cassegrain telescope (primary mirror 1100 mm x 730 mm ellipse) and two separate instruments (FGS and AIRS) covering simultaneously 0.5-7.8 micron spectral range. The satellite is best placed into an L2 orbit to maximise the thermal stability and the field of regard. The payload module is passively cooled via a series of V-Groove radiators; the detectors for the AIRS are the only items that require active cooling via an active Ne JT cooler. The Ariel payload is developed by a consortium of more than 50 institutes from 16 ESA countries, which include the UK, France, Italy, Belgium, Poland, Spain, Austria, Denmark, Ireland, Portugal, Czech Republic, Hungary, the Netherlands, Sweden, Norway, Estonia, and a NASA contribution. △ Less

Submitted 10 April, 2021; originally announced April 2021.

Comments: Ariel Definition Study Report, 147 pages. Reviewed by ESA Science Advisory Structure in November 2020. Original document available at: https://www.cosmos.esa.int/documents/1783156/3267291/Ariel_RedBook_Nov2020.pdf/

Report number: ESA/SCI(2020)1

Enhanced X-ray Emission Coinciding with Giant Radio Pulses from the Crab Pulsar

Authors: Teruaki Enoto, Toshio Terasawa, Shota Kisaka, Chin-Ping Hu, Sebastien Guillot, Natalia Lewandowska, Christian Malacaria, Paul S. Ray, Wynn C. G. Ho, Alice K. Harding, Takashi Okajima, Zaven Arzoumanian, Keith C. Gendreau, Zorawar Wadiasingh, Craig B. Markwardt, Yang Soong, Steve Kenyon, Slavko Bogdanov, Walid A. Majid, Tolga Güver, Gaurava K. Jaisawal, Rick Foster, Yasuhiro Murata, Hiroshi Takeuchi, Kazuhiro Takefuji , et al. (12 additional authors not shown)

Abstract: Giant radio pulses (GRPs) are sporadic bursts emitted by some pulsars, lasting a few microseconds. GRPs are hundreds to thousands of times brighter than regular pulses from these sources. The only GRP-associated emission outside radio wavelengths is from the Crab Pulsar, where optical emission is enhanced by a few percent during GRPs. We observed the Crab Pulsar simultaneously at X-ray and radio w… ▽ More Giant radio pulses (GRPs) are sporadic bursts emitted by some pulsars, lasting a few microseconds. GRPs are hundreds to thousands of times brighter than regular pulses from these sources. The only GRP-associated emission outside radio wavelengths is from the Crab Pulsar, where optical emission is enhanced by a few percent during GRPs. We observed the Crab Pulsar simultaneously at X-ray and radio wavelengths, finding enhancement of the X-ray emission by $3.8\pm0.7\%$ (a 5.4$σ$ detection) coinciding with GRPs. This implies that the total emitted energy from GRPs is tens to hundreds of times higher than previously known. We discuss the implications for the pulsar emission mechanism and extragalactic fast radio bursts. △ Less

Submitted 7 April, 2021; originally announced April 2021.

Comments: 63 pages, 23 figures, 7 Tables. Published in Science

Estimating the Effects of Urban Green Regions in terms of Diffusion

Authors: Eric K. Tokuda, Florence A. S. Shibata, Henrique F. de Arruda, Guilherme S. Domingues, Cesar H. Comin, Luciano da F. Costa, Roberto M. Cesar-Jr

Abstract: The interaction between cities and their respective green regions corresponds to an interesting issue that has received growing attention over the last decades. These relationships have multiple natures, ranging from providing habitat for animal life to temperature and humidity dynamics. Several methods based on area, size, shape, and distance have been considered in the literature. Given that sev… ▽ More The interaction between cities and their respective green regions corresponds to an interesting issue that has received growing attention over the last decades. These relationships have multiple natures, ranging from providing habitat for animal life to temperature and humidity dynamics. Several methods based on area, size, shape, and distance have been considered in the literature. Given that several important contributions of green regions to urban areas involve temperature, humidity and gases (e.g.~oxygen) exchanges, which are intrinsically related to physical diffusion, it becomes particularly interesting to simulate the diffusion of green effects over urban areas as a means of better understanding the respective influences. The present work reports a related approach. Once the green regions of a given city are automatically identified by semantic segmentation and have eventual artifacts eliminated, successive convolutions are applied as a means to obtain the unfolding of the diffusion of the green effects along time. As illustrated, the diffusion dynamics is intrinsically interesting because it can be strongly affected by the spatial distribution of the green mass. In particular, we observed that smaller green regions could substantially contribute to the diffusion. The reported approach has been illustrated with respect to the Brazilian city of Ribeirao Preto, whose small and medium-sized green regions were found to complement in an effective manner the diffusion of the green effects as inferred from the performed simulations under specific parameter settings. △ Less

Submitted 22 March, 2021; originally announced March 2021.

Solar Neutron Decay Protons observed in November 7, 2004

Authors: Yasushi Muraki, Jose F. Valde-Galicia, Ernesto Ortiz, Yutaka Matsubara, Shoichi Shibata, Takashi Sako, Satoshi Masuda, Munetoshi Tokumaru, Tatsumi Koi, Akitoshi Ooshima, Takasuke Sakai, Tsuguya Naito, Pedro Miranda

Abstract: We have found an interesting event registered by the solar neutron telescopes installed at high mountains in Bolivia (5250 m a.s.l.) and Mexico (4600 m a.s.l.). The event was observed November 7th of 2004 in association with a large solar flare of magnitude X2.0. Some features in our registers and in two satellites (GOES 11 and SOHO) reveal the presence of electrons and protons as possible product… ▽ More We have found an interesting event registered by the solar neutron telescopes installed at high mountains in Bolivia (5250 m a.s.l.) and Mexico (4600 m a.s.l.). The event was observed November 7th of 2004 in association with a large solar flare of magnitude X2.0. Some features in our registers and in two satellites (GOES 11 and SOHO) reveal the presence of electrons and protons as possible products of neutron decay. Solar neutron decay protons (sndp) were recorded on board ISEE3 satellite in June 3rd, 1982 . On October 19th, 1989, the ground level detectors installed in Goose Bay and Deep River revealed the registration of solar neutron decay protons (sndp). Therefore this is the second example that such an evidence is registered on the Earth's surface. △ Less

Submitted 31 December, 2020; originally announced December 2020.

Comments: total 21 pages, including 13 figures submitted to Earth, Planet and Space

Low X-ray Efficiency of a Young High-B Pulsar PSR J1208-6238 Observed with Chandra

Authors: Aya Bamba, Eri Watanabe, Koji Mori, Shinpei Shibata, Yukikatsu Terada, Hidetoshi Sano, Miroslav D. Filipovic

Abstract: High magnetic field (high-B) pulsars are key sources to bridge magnetars and conventional rotation powered pulsars, and thus to understand the origin of magnetar activities. We have estimated a tight upper-limit on the X-ray flux of one of the youngest high-B pulsars PSR J1208-6238 for the first time; a Chandra 10 ks observation shows no significant source. Depending on the emission models, the 3s… ▽ More High magnetic field (high-B) pulsars are key sources to bridge magnetars and conventional rotation powered pulsars, and thus to understand the origin of magnetar activities. We have estimated a tight upper-limit on the X-ray flux of one of the youngest high-B pulsars PSR J1208-6238 for the first time; a Chandra 10 ks observation shows no significant source. Depending on the emission models, the 3sigma upper-limit on the intrinsic 0.5-7 keV flux to (2.2-10.0)e-14 erg/s/cm2. △ Less

Submitted 6 November, 2020; originally announced November 2020.

Comments: 7 pages, 6 figures, Ap&SS, accepted

Report number: RESCEU-21/20

Journal ref: Ap&SS, 365, 178 (2020)

Development of an optical photon-counting imager with a monolithic Geiger APD array

Authors: Takeshi Nakamori, Yuga Ouchi, Risa Ogihara, Toshio Terasawa, Yuhei Kato, Shinpei Shibata

Abstract: We have developed a sensor system based on an optical photon-counting imager with high timing resolution, aiming for highly time-variable astronomical phenomena. The detector is a monolithic Geiger-mode avalanche photodiode array customized in a Multi-Pixel Photon Counter with a response time on the order of nanoseconds. This paper evaluates the basic performance of the sensor and confirms the gai… ▽ More We have developed a sensor system based on an optical photon-counting imager with high timing resolution, aiming for highly time-variable astronomical phenomena. The detector is a monolithic Geiger-mode avalanche photodiode array customized in a Multi-Pixel Photon Counter with a response time on the order of nanoseconds. This paper evaluates the basic performance of the sensor and confirms the gain linearity, uniformity, and low dark count. We demonstrate the system's ability to detect the period of a flashing LED, using a data acquisition system developed to obtain the light curve with a time bin of 100 microseconds. The Crab pulsar was observed using a 35-cm telescope without cooling, and the equipment detected optical pulses with a period consistent with the data from the radio ephemeris. Although improvements to the system will be necessary for more reliability, the system has been proven to be a promising device for exploring the time-domain optical astronomy. △ Less

Submitted 22 October, 2020; originally announced October 2020.

Comments: 11 pages, 18 figures, accepted for publication in PASJ

Energy sensitivity of the GRAPES-3 EAS array for primary cosmic ray protons

Authors: B. Hariharan, S. Ahmad, M. Chakraborty, A. Chandra, S. R. Dugad, S. K. Gupta, Y. Hayashi, H. Kojima, S. S. R. Inbanathan, P. Jagadeesan, A. Jain, P. Jain, V. B. Jhansi, S. Kawakami, P. K. Mohanty, S. D. Morris, P. K. Nayak, A. Oshima, D. Pattanaik, P. S. Rakshe, K. Ramesh, B. S. Rao, L. V. Reddy, S. Shibata, F. Varsi , et al. (1 additional authors not shown)

Abstract: Low energy ground-based cosmic ray air shower experiments generally have energy threshold in the range of a few tens to a few hundreds of TeV. The shower observables are measured indirectly with an array of detectors. The atmospheric absorption of low energy secondaries limits their detection frequencies at the Earth's surface. However, due to selection effects, a tiny fraction of low energy showe… ▽ More Low energy ground-based cosmic ray air shower experiments generally have energy threshold in the range of a few tens to a few hundreds of TeV. The shower observables are measured indirectly with an array of detectors. The atmospheric absorption of low energy secondaries limits their detection frequencies at the Earth's surface. However, due to selection effects, a tiny fraction of low energy showers, which are produced in the lower atmosphere can reach the observational level. But, due to less information of shower observables, the reconstruction of these showers are arduous. Hence, it is believed that direct measurements by experiments aboard on satellites and balloon flights are more reliable at low energies. Despite having very small efficiency ($\sim$0.1%) at low energies, the large acceptance ($\sim$5 m$^2$sr) of GRAPES-3 experiment allows observing primary cosmic rays down below to $\sim$1 TeV and opens up the possibility to measure primary energy spectrum spanning from a few TeV to beyond cosmic ray knee (up to 10$^{16}$ eV), covering five orders of magnitude. The GRAPES-3 energy threshold for primary protons through Monte Carlo simulations are calculated, which gives reasonably good agreement with data. Furthermore, the total efficiencies and acceptance are also calculated for protons primaries. The ability of GRAPES-3 experiment to cover such a broader energy range may provide a unique handle to bridge the energy spectrum between direct measurements at low energies and indirect measurements at ultra-high energies. △ Less

Submitted 3 May, 2020; originally announced May 2020.

Comments: 15 pages, 6 figures

The angular resolution of GRAPES-3 EAS array after correction for the shower front curvature

Authors: V. B. Jhansi, S. Ahmad, M. Chakraborty, S. R. Dugad, S. K. Gupta, B. Hariharan, Y. Hayashi, P. Jagadeesan, A. Jain, P. Jain, S. Kawakami, H. Kojima, S. Mahapatra, P. K. Mohanty, S. D. Morris, P. K. Nayak, A. Oshima, D. Pattanaik, P. S. Rakshe, K. Ramesh, B. S. Rao, L. V. Reddy, S. Shibata, F. Varsi, M. Zuberi

Abstract: The angular resolution of an extensive air shower (EAS) array plays a critical role in determining its sensitivity for the detection of point $γ$-ray sources in the multi-TeV energy range. The GRAPES-3 an EAS array located at Ooty in India (11.4$^{\circ}$N, 76.7$^{\circ}$E, 2200 m altitude) is designed to study $γ$-rays in the TeV-PeV energy range. It comprises of a dense array of 400 plastic scin… ▽ More The angular resolution of an extensive air shower (EAS) array plays a critical role in determining its sensitivity for the detection of point $γ$-ray sources in the multi-TeV energy range. The GRAPES-3 an EAS array located at Ooty in India (11.4$^{\circ}$N, 76.7$^{\circ}$E, 2200 m altitude) is designed to study $γ$-rays in the TeV-PeV energy range. It comprises of a dense array of 400 plastic scintillators deployed over an area of 25000 m$^2$ and a large area (560 m$^2$) muon telescope. A new statistical method allowed real time determination of the propagation delay of each detector in the GRAPES-3 array. The shape of shower front is known to be curved and here the details of a new method developed for accurate measurement of the shower front curvature is presented. These two developments have led to a sizable improvement in the angular resolution of GRAPES-3 array. It is shown that the curvature depends on the size and age of an EAS. By employing two different techniques, namely, the odd-even and the left-right methods, independent estimates of the angular resolution are obtained. The odd-even method estimates the best achievable resolution of the array. For obtaining the angular resolution, the left-right method is used after implementing the size and age dependent curvature corrections. A comparison of the angular resolution as a function of EAS energy by these two methods shows them be virtually indistinguishable. The angular resolution of GRAPES-3 array is 47$^{\prime}$ for energies E$>$5 TeV and improves to 17$^{\prime}$ at E$>$100 TeV and finally approaching 10$^{\prime}$ at E$>$500 TeV. △ Less

Submitted 12 November, 2019; originally announced November 2019.

The origin of the high metallicity of close-in giant exoplanets: Combined effect of the resonant and aerodynamic shepherding

Authors: Sho Shibata, Ravit Helled, Masahiro Ikoma

Abstract: Context.Recent studies suggest that many giant exoplanets are highly enriched with heavy elements compared to their host star andcontain several tens of Earth masses or more of heavy elements. Such enrichment is considered to have been brought by accretionof planetesimals in late formation stages. Previous dynamical simulations, however, show that planets are unable to collect so muchheavy element… ▽ More Context.Recent studies suggest that many giant exoplanets are highly enriched with heavy elements compared to their host star andcontain several tens of Earth masses or more of heavy elements. Such enrichment is considered to have been brought by accretionof planetesimals in late formation stages. Previous dynamical simulations, however, show that planets are unable to collect so muchheavy elements throughin situplanetesimal accretion. Aims.We investigate whether a giant planet migrating inward can capture planetesimals efficiently to significantly increase its metal-licity. Methods.We performed orbital integrations of a migrating giant planet and planetesimals in a protoplanetary gas disc to infer theplanetesimal mass that is accreted by the planet. Results.We find that the two shepherding processes of mean motion resonances trapping and aerodynamic gas drag inhibit plan-etesimal capture of a migrating planet. However, the amplified libration allows the highly-excited planetesimals in the resonances toescape from the resonance trap and be accreted by the planet. Consequently, we show that a migrating giant planet captures planetes-imals with total mass of several tens of Earth masses, if the planet forms at a few tens of AU in a relatively massive disc. We alsofind that planetesimal capture occurs efficiently in a limited range of semi-major axis, and that the total captured planetesimal massincreases with increasing migration distances. Our results have important implications for understanding the relation between giantplanet metallicity and mass, as we suggest that it reflects the formation location of the planet, or more precisely, the location whererunaway gas accretion occurred. We also suggest the observed metal-rich close-in Jupiters migrated to their present locations fromafar, where they formed. △ Less

Submitted 6 November, 2019; originally announced November 2019.

Comments: 14 pages, 11 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 633, A33 (2020)

CO Multi-line Imaging of Nearby Galaxies (COMING) IV. Overview of the Project

Authors: Kazuo Sorai, Nario Kuno, Kazuyuki Muraoka, Yusuke Miyamoto, Hiroyuki Kaneko, Hiroyuki Nakanishi, Naomasa Nakai, Kazuki Yanagitani, Takahiro Tanaka, Yuya Sato, Dragan Salak, Michiko Umei, Kana Morokuma-Matsui, Naoko Matsumoto, Saeko Ueno, Hsi-An Pan, Yuto Noma, Tsutomu T. Takeuchi, Moe Yoda, Mayu Kuroda, Atsushi Yasuda, Yoshiyuki Yajima, Nagisa Oi, Shugo Shibata, Masumichi Seta , et al. (14 additional authors not shown)

Abstract: Observations of the molecular gas in galaxies are vital to understanding the evolution and star-forming histories of galaxies. However, galaxies with molecular gas maps of their whole discs having sufficient resolution to distinguish galactic structures are severely lacking. Millimeter wavelength studies at a high angular resolution across multiple lines and transitions are particularly needed, se… ▽ More Observations of the molecular gas in galaxies are vital to understanding the evolution and star-forming histories of galaxies. However, galaxies with molecular gas maps of their whole discs having sufficient resolution to distinguish galactic structures are severely lacking. Millimeter wavelength studies at a high angular resolution across multiple lines and transitions are particularly needed, severely limiting our ability to infer the universal properties of molecular gas in galaxies. Hence, we conducted a legacy project with the 45 m telescope of the Nobeyama Radio Observatory, called the CO Multi-line Imaging of Nearby Galaxies (COMING), which simultaneously observed 147 galaxies with high far-infrared flux in $^{12}$CO, $^{13}$CO, and C$^{18}$O $J=1-0$ lines. The total molecular gas mass was derived using the standard CO-to-H$_2$ conversion factor and found to be positively correlated with the total stellar mass derived from the WISE $3.4 μ$m band data. The fraction of the total molecular gas mass to the total stellar mass in galaxies does not depend on their Hubble types nor the existence of a galactic bar, although when galaxies in individual morphological types are investigated separately, the fraction seems to decrease with the total stellar mass in early-type galaxies and vice versa in late-type galaxies. No differences in the distribution of the total molecular gas mass, stellar mass, and the total molecular gas to stellar mass ratio was observed between barred and non-barred galaxies, which is likely the result of our sample selection criteria, in that we prioritized observing FIR bright (and thus molecular gas-rich) galaxies. △ Less

Submitted 9 October, 2019; originally announced October 2019.

Comments: Accepted for publication in PASJ; 47 pages, 5 tables, 29 figures. On-line supplementary images are available at this URL (https://astro3.sci.hokudai.ac.jp/~radio/coming/publications/). CO data is available at the Japanese Virtual Observatory (JVO) website (https://jvo.nao.ac.jp/portal/nobeyama/coming.do) and the project website (https://astro3.sci.hokudai.ac.jp/~radio/coming/data/)

Solar Neutrons Observed from September 4 to 10, 2017 by SEDA-FIB

Authors: K. Kamiya, K. Koga, H. Matsumoto, S. Masuda, Y. Muraki, H. Tajima, S. Shibata

Abstract: The SEDA-FIB is a detector designed to measure solar neutrons. This solar neutron detector was operated onboard the ISS on July 16, 2009 and March 31, 2018. Eighteen large solar flares were later observed by the GOES satellite in solar active region 12673 that appeared on September 4 and lasted until September 10, 2017, with intensity higher than > M2. In nine of those solar flares, the SEDA-FIB d… ▽ More The SEDA-FIB is a detector designed to measure solar neutrons. This solar neutron detector was operated onboard the ISS on July 16, 2009 and March 31, 2018. Eighteen large solar flares were later observed by the GOES satellite in solar active region 12673 that appeared on September 4 and lasted until September 10, 2017, with intensity higher than > M2. In nine of those solar flares, the SEDA-FIB detected clear signals of solar neutrons, along with five minor excesses. Among these events, we focus on two associated with the flares of X2.2 (SOL2017-09-06) and X8.2 (SOL2017-09-10) that share a common feature: a process of accelerating electrons into high energies as clearly recorded by the FERMI-GBM detector. These events may provide us with useful information to elucidate the ion acceleration process. The X8.2 event was a limb flare that proved adequate for fixing the parameters needed to explain the process of particle acceleration into high energies. According to our analysis, the electron acceleration process may possibly be explained by the shock acceleration model. However, we found that it would be difficult to explain the simultaneous acceleration of ions with electrons, unless the ions were preheated prior to their rapid acceleration. △ Less

Submitted 22 July, 2019; originally announced July 2019.

Comments: submitted to International Conference on Cosmic Ras and Astroparticle Physics, held in Madison, 2019

Report number: icrc2019_1150.pdf

Journal ref: Proceeding of Science, 2019

Capture of Solids by Growing Proto-gas Giants: Effects of Gap Formation and Supply-limited Growth

Authors: Sho Shibata, Masahiro Ikoma

Abstract: Studies of internal structure of gas giant planets suggest that their envelopes are enriched with heavier elements than hydrogen and helium relative to their central stars. Such enrichment likely occurred by solid accretion during late formation stages of gas giant planets in which gas accretion dominates protoplanetary growth. Some previous studies performed orbital integration of planetesimals a… ▽ More Studies of internal structure of gas giant planets suggest that their envelopes are enriched with heavier elements than hydrogen and helium relative to their central stars. Such enrichment likely occurred by solid accretion during late formation stages of gas giant planets in which gas accretion dominates protoplanetary growth. Some previous studies performed orbital integration of planetesimals around a growing protoplanet with the assumption of an uniform circumstellar disc to investigate how efficiently the protoplanet captures planetesimals. However, not only planetesimals but also disc gas are gravitationally perturbed by the protoplanet in its late formation stages, resulting in gap opening in the circumstellar disc. In this study, we investigate the effects of such gap formation on the capture of planetesimals by performing dynamical simulations of planetesimals around a growing proto-gas giant planet. Gap formation reduces the surface density of disc gas, makes a steep pressure gradient and limits the growth rate of the protoplanet. We find that the first effect enhances the capture of planetesimals, while the others reduce it. Consequently the amount of planetesimals captured during the gas accretion is estimated to be at most ~3 Mearth. We conclude that the in-situ capture of planetesimals needs the initial solid surface density more than five times higher than that of the minimum mass solar nebula for explaining the inferred large amount of heavy element in Jupiter. For highly dense warm Jupiters, we would need additional processes enhancing the capture and/or supply of planetesimals. △ Less

Submitted 13 June, 2019; originally announced June 2019.

Comments: 15 pages, 13 figures, this is a pre-copyedited, author-produced version of an article accepted for publication in MNRAS following peer review

A high-magnetic-field radio pulsar survey with Swift/XRT

Authors: Eri Watanabe, Shinpei Shibata, Takanori Sakamoto, Aya Bamba

Abstract: We present the X-ray survey results of high-magnetic-field radio pulsars (high-B PSRs) with Swift/XRT. X-ray observations of the rotation-powered pulsars with the dipole magnetic field ${\it B_{\rm d}}$ near the quantum critical field ${\it B_{\rm q}}=4.4\times10^{13}$ G is of great importance for understanding the transition between the rotation-powered pulsars and the magnetars, because there ar… ▽ More We present the X-ray survey results of high-magnetic-field radio pulsars (high-B PSRs) with Swift/XRT. X-ray observations of the rotation-powered pulsars with the dipole magnetic field ${\it B_{\rm d}}$ near the quantum critical field ${\it B_{\rm q}}=4.4\times10^{13}$ G is of great importance for understanding the transition between the rotation-powered pulsars and the magnetars, because there are a few objects that have magnetar-like properties. Out of the 27 high-B PSRs that are in the ATNF pulsar catalogue but have not been reported or have no effective upper-limits in the X-ray bands, we analyze the Swift/XRT data for 21 objects, where 6 objects are newly observed and 15 objects are taken from the archival data. As a result, we have new $3σ$ upper-limits for all the 21 objects. Since the upper-limits are tight, we conclude that we do not find any magnetar-like high-B PSRs such as PSR J1819$-$1458. The probability of the high X-ray efficiency in the high-B PSRs is obtained to be $11\%-29\%$. Combining the previous observations, we discuss which parameter causes magnetar-like properties. It may be suggested that the magnetar-like properties appear only when ${\it B_{\rm d}}\gtrsim10^{13.5}$ G for the radio pulsar population. This is true even if the radio-quiet high-B RPP are included. △ Less

Submitted 30 May, 2019; originally announced May 2019.

Comments: 19 pages, 7 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 486, Issue 4, July 2019, Pages 5323-5334

XMM-Newton Spectrum of the magnetar CXOU J171405.7-381031

Authors: Haruka Watanabe, Aya Bamba, Shinpei Shibata, Eri Watanabe

Abstract: We observe the magnetar CXOU J171405.7-381031 with XMM-Newton and obtain the most reliable X-ray spectral parameters for this magnetar. After removing the flux from the surrounding supernova remnant CTB~37B, the radiation of CXOU J171405.7-381031 is best described by a two-component model, consisting of a blackbody and power law. We obtain a blackbody temperature of 0.58^{+0.03}_{-0.03} keV, photo… ▽ More We observe the magnetar CXOU J171405.7-381031 with XMM-Newton and obtain the most reliable X-ray spectral parameters for this magnetar. After removing the flux from the surrounding supernova remnant CTB~37B, the radiation of CXOU J171405.7-381031 is best described by a two-component model, consisting of a blackbody and power law. We obtain a blackbody temperature of 0.58^{+0.03}_{-0.03} keV, photon index of 2.15^{+0.62}_{-0.68}, and unabsorbed 2-10 keV flux of 2.33^{+0.02}_{-0.02} x 10^{-12} erg cm^{-2} s^{-1}. These new parameters enable us to compare CXOU 171405.7-381031 with other magnetars, and it is found that the luminosity, temperature and the photon index of CXOU J171405.7-381031 are aligned with the known trend among the magnetar population with a slightly higher temperature, which could be caused by its young age. All the magnetars with a spin-down age of less than 1~kyr show time variation or bursts except for CXOU J171405.7-381031. We explore the time variability for ten observations in between 2006 and 2015, but there is no variation larger than \sim 10%. △ Less

Submitted 27 May, 2019; originally announced May 2019.

Comments: 10 pages, 4figures, Accepted for publication in Publication of the Astronomical Society of Japan

Measurement of the Electrical Properties of a Thundercloud Through Muon Imaging by the GRAPES-3 Experiment

Authors: B. Hariharan, A. Chandra, S. R. Dugad, S. K. Gupta, P. Jagadeesan, A. Jain, P. K. Mohanty, S. D. Morris, P. K. Nayak, P. S. Rakshe, K. Ramesh, B. S. Rao, L. V. Reddy, M. Zuberi, Y. Hayashi, S. Kawakami, S. Ahmad, H. Kojima, A. Oshima, S. Shibata, Y. Muraki, K. Tanaka

Abstract: The GRAPES-3 muon telescope located in Ooty, India records rapid ($\sim$10 min) variations in the muon intensity during major thunderstorms. Out of a total of 184 thunderstorms recorded during the interval April 2011-December 2014, the one on 1 December 2014 produced a massive potential of 1.3 GV. The electric field measured by four well-separated (up to 6 km) monitors on the ground was used to he… ▽ More The GRAPES-3 muon telescope located in Ooty, India records rapid ($\sim$10 min) variations in the muon intensity during major thunderstorms. Out of a total of 184 thunderstorms recorded during the interval April 2011-December 2014, the one on 1 December 2014 produced a massive potential of 1.3 GV. The electric field measured by four well-separated (up to 6 km) monitors on the ground was used to help estimate some of the properties of this thundercloud including its altitude and area that were found to be 11.4 km above mean sea level (amsl) and $\geq$380 km$^2$, respectively. A charging time of 6 min to reach 1.3 GV implied the delivery of a power of $\geq$2 GW by this thundercloud that was moving at a speed of $\sim$60 km h$^{-1}$. This work possibly provides the first direct evidence for the generation of GV potentials in thunderclouds that could also possibly explain the production of highest energy (100 MeV) $γ$-rays in the terrestrial $γ$-ray flashes. △ Less

Submitted 23 March, 2019; originally announced March 2019.

Comments: Received 6 January 2019, Revised 21 January 2019, Published 15 March 2019

Journal ref: Phys. Rev. Lett. 122, 105101 (2019)

CO Multi-line Imaging of Nearby Galaxies (COMING). III. Dynamical effect on molecular gas density and star formation in the barred spiral galaxy NGC 4303

Authors: Yoshiyuki Yajima, Kazuo Sorai, Nario Kuno, Kazuyuki Muraoka, Yusuke Miyamoto, Hiroyuki Kaneko, Hiroyuki Nakanishi, Naomasa Nakai, Takahiro Tanaka, Yuya Sato, Dragan Salak, Kana Morokuma-Matsui, Naoko Matsumoto, His-An Pan, Yuto Noma, Tsutomu T. Takeuchi, Moe Yoda, Mayu Kuroda, Atsushi Yasuda, Nagisa Oi, Shugo Shibata, Masumichi Seta, Yoshimasa Watanabe, Shoichiro Kita, Ryusei Komatsuzaki , et al. (2 additional authors not shown)

Abstract: We present the results of $^{12}$CO($J$=1-0) and $^{13}$CO($J$=1-0) simultaneous mappings toward the nearby barred spiral galaxy NGC 4303 as a part of the CO Multi-line Imaging of Nearby Galaxies (COMING) project. Barred spiral galaxies often show lower star-formation efficiency (SFE) in their bar region compared to the spiral arms. In this paper, we examine the relation between the SFEs and the v… ▽ More We present the results of $^{12}$CO($J$=1-0) and $^{13}$CO($J$=1-0) simultaneous mappings toward the nearby barred spiral galaxy NGC 4303 as a part of the CO Multi-line Imaging of Nearby Galaxies (COMING) project. Barred spiral galaxies often show lower star-formation efficiency (SFE) in their bar region compared to the spiral arms. In this paper, we examine the relation between the SFEs and the volume densities of molecular gas $n(\rm{H}_2)$ in the eight different regions within the galactic disk with CO data combined with archival far-ultraviolet and 24 $μ$m data. We confirmed that SFE in the bar region is lower by 39% than that in the spiral arms. Moreover, velocity-alignment stacking analysis was performed for the spectra in the individual regions. The integrated intensity ratios of $^{12}$CO to $^{13}$CO ($R_{12/13}$) range from 10 to 17 as the results of stacking. Fixing a kinetic temperature of molecular gas, $n(\rm{H}_2)$ was derived from $R_{12/13}$ via non-local thermodynamic equilibrium (non-LTE) analysis. The density $n(\rm{H}_2)$ in the bar is lower by 31-37% than that in the arms and there is a rather tight positive correlation between SFEs and $n(\rm{H}_2)$, with a correlation coefficient of $\sim 0.8$. Furthermore, we found a dependence of $n(\rm{H}_2)$ on the velocity dispersion of inter-molecular clouds ($ΔV/ \sin i$). Specifically, $n(\rm{H}_2)$ increases as $ΔV/ \sin i$ increases when $ΔV/ \sin i < 100$ km s$^{-1}$. On the other hand, $n(\rm{H}_2)$ decreases as $ΔV/ \sin i$ increases when $ΔV/ \sin i > 100$ km s$^{-1}$. These relations indicate that the variations of SFE could be caused by the volume densities of molecular gas, and the volume densities could be governed by the dynamical influence such as cloud-cloud collisions, shear and enhanced inner-cloud turbulence. △ Less

Submitted 12 February, 2019; originally announced February 2019.

Comments: 15 pages, 8 figures, accepted for publication in PASJ

Journal ref: 2019, PASJ, 71, S13

Angular and Abundance Distribution of High-energy Gamma Rays and Neutrons Simulated by GEANT4 Code for Solar Flares

Authors: K. Kamiya, K. Koga, H. Matsumoto, Y. Muraki, S. Shibata

Abstract: In the solar flare observed on June 3, 2012, high energy gamma-rays and neutrons were observed. The event includes a remarkable feature of a high neutron/gamma-ratio in the secondary particles. We have examined whether this high n/$γ$-ratio can be explained by simulation. As a result of simulations using the GEANT4 program, the high n/$γ$-ratio may be reproduced for the case that helium and other… ▽ More In the solar flare observed on June 3, 2012, high energy gamma-rays and neutrons were observed. The event includes a remarkable feature of a high neutron/gamma-ratio in the secondary particles. We have examined whether this high n/$γ$-ratio can be explained by simulation. As a result of simulations using the GEANT4 program, the high n/$γ$-ratio may be reproduced for the case that helium and other heavy ions were dominantly accelerated in the flare. △ Less

Submitted 7 February, 2019; originally announced February 2019.

Comments: submitted to the Proceeding of The 20th International Symposium on Very High Energy Cosmic Ray Interaction (ISVHECRI 2018, Nagoya, Japan), Europian Physics Journal

CO Multi-line Imaging of Nearby Galaxies (COMING): VI. Radial variations in star formation efficiency

Authors: Kazuyuki Muraoka, Kazuo Sorai, Yusuke Miyamoto, Moe Yoda, Kana Morokuma-matsui, Masato I. N. Kobayashi, Mayu Kuroda, Hiroyuki Kaneko, Nario Kuno, Tsutomu T. Takeuchi, Hiroyuki Nakanishi, Yoshimasa Watanabe, Takahiro Tanaka, Atsushi Yasuda, Yoshiyuki Yajima, Shugo Shibata, Dragan Salak, Daniel Espada, Naoko Matsumoto, Yuto Noma, Shoichiro Kita, Ryusei Komatsuzaki, Ayumi Kajikawa, Yu Yashima, Hsi-An Pan , et al. (3 additional authors not shown)

Abstract: We examined radial variations in molecular-gas based star formation efficiency (SFE), which is defined as star formation rate per unit molecular gas mass, for 80 galaxies selected from the CO Multi-line Imaging of Nearby Galaxies project (Sorai et al. 2019). The radial variations in SFE for individual galaxies are typically a factor of 2 -- 3, which suggests that SFE is nearly constant along galac… ▽ More We examined radial variations in molecular-gas based star formation efficiency (SFE), which is defined as star formation rate per unit molecular gas mass, for 80 galaxies selected from the CO Multi-line Imaging of Nearby Galaxies project (Sorai et al. 2019). The radial variations in SFE for individual galaxies are typically a factor of 2 -- 3, which suggests that SFE is nearly constant along galactocentric radius. We found the averaged SFE in 80 galaxies of $(1.69 \pm 1.1) \times 10^{-9}$ yr$^{-1}$, which is consistent with Leroy et al. 2008 if we consider the contribution of helium to the molecular gas mass evaluation and the difference in the assumed initial mass function between two studies. We compared SFE among different morphological (i.e., SA, SAB, and SB) types, and found that SFE within the inner radii ($r/r_{25} < 0.3$, where $r_{25}$ is $B$-band isophotal radius at 25 mag arcsec$^{-2}$) of SB galaxies is slightly higher than that of SA and SAB galaxies. This trend can be partly explained by the dependence of SFE on global stellar mass, which probably relates to the CO-to-H$_2$ conversion factor through the metallicity. For two representative SB galaxies in our sample, NGC 3367 and NGC 7479, the ellipse of $r/r_{25}$ = 0.3 seems to cover not only the central region but also the inner part of the disk, mainly the bar. These two galaxies show higher SFE in the bar than in spiral arms. However, we found an opposite trend in NGC 4303; SFE is lower in the bar than in spiral arms, which is consistent with earlier studies (e.g., Momose et al. 2010). These results suggest diversity of star formation activities in the bar. △ Less

Submitted 30 January, 2019; originally announced January 2019.

Comments: 28 pages, 12 figures, accepted for publication in PASJ

CO Multi-line Imaging of Nearby Galaxies (COMING). VII. Fourier decomposition of molecular gas velocity fields and bar pattern speed

Authors: Dragan Salak, Yuto Noma, Kazuo Sorai, Yusuke Miyamoto, Nario Kuno, Alex R. Pettitt, Hiroyuki Kaneko, Takahiro Tanaka, Atsushi Yasuda, Shoichiro Kita, Yoshiyuki Yajima, Shugo Shibata, Naomasa Nakai, Masumichi Seta, Kazuyuki Muraoka, Mayu Kuroda, Hiroyuki Nakanishi, Tsutomu T. Takeuchi, Moe Yoda, Kana Morokuma-Matsui, Yoshimasa Watanabe, Naoko Matsumoto, Nagisa Oi, Hsi-An Pan, Ayumi Kajikawa , et al. (2 additional authors not shown)

Abstract: The $^{12}$CO $(J=1\rightarrow0)$ velocity fields of a sample of 20 nearby spiral galaxies, selected from the CO Multi-line Imaging of Nearby Galaxies (COMING) legacy project of Nobeyama Radio Observatory, have been analyzed by Fourier decomposition to determine their basic kinematic properties, such as circular and noncircular velocities. On average, the investigated barred (SAB and SB) galaxies… ▽ More The $^{12}$CO $(J=1\rightarrow0)$ velocity fields of a sample of 20 nearby spiral galaxies, selected from the CO Multi-line Imaging of Nearby Galaxies (COMING) legacy project of Nobeyama Radio Observatory, have been analyzed by Fourier decomposition to determine their basic kinematic properties, such as circular and noncircular velocities. On average, the investigated barred (SAB and SB) galaxies exhibit a ratio of noncircular to circular velocities of molecular gas larger by a factor of 1.5-2 than non-barred (SA) spiral galaxies at radii within the bar semimajor axis $a_\mathrm{b}$ at 1 kpc resolution, with a maximum at a radius of $R/a_\mathrm{b}\sim0.3$. Residual velocity field images, created by subtracting model velocity fields from the data, reveal that this trend is caused by kpc-scale streaming motions of molecular gas in the bar region. Applying a new method based on radial velocity reversal, we estimated the corotation radius $R_\mathrm{CR}$ and bar pattern speed $Ω_\mathrm{b}$ in seven SAB and SB systems. The ratio of the corotation to bar radius is found to be in a range of $\mathcal{R}\equiv R_\mathrm{CR}/a_\mathrm{b}\sim0.8\mathrm{-}1.6$, suggesting that intermediate (SBb-SBc), luminous barred spiral galaxies host fast and slow rotator bars. Tentative negative correlations are found for $Ω_\mathrm{b}$ vs. $a_\mathrm{b}$ and $Ω_\mathrm{b}$ vs. total stellar mass $M_\ast$, indicating that bars in massive disks are larger and rotate slower, possibly a consequence of angular momentum transfer. The kinematic properties of SAB and SB galaxies, derived from Fourier decomposition, are compared with recent numerical simulations that incorporate various rotation curve models and galaxy interactions. △ Less

Submitted 3 January, 2019; originally announced January 2019.

Comments: Accepted for publication in PASJ; 39 pages, 15 figures

Was the cosmic ray burst detected by the GRAPES-3 on 22 June 2015 caused by transient weakening of geomagnetic field or by an interplanetary anisotropy?

Authors: P. K. Mohanty, K. P. Arunbabu, T. Aziz, S. R. Dugad, S. K. Gupta, B. Hariharan, P. Jagadeesan, A. Jain, S. D. Morris, P. K. Nayak, P. S. Rakshe, K. Ramesh, B. S. Rao, M. Zuberi, Y. Hayashi, S. Kawakami, P. Subramanian, S. Raha, S. Ahmad, A. Oshima, S. Shibata, H. Kojima

Abstract: The GRAPES-3 muon telescope in Ooty, India had claimed detection of a 2 hour (h) high-energy ($\sim$20 GeV) burst of galactic cosmic-rays (GCRs) through a $>$50$σ$ surge in GeV muons, was caused by reconnection of the interplanetary magnetic field (IMF) in the magnetosphere that led to transient weakening of Earth's magnetic shield. This burst had occurred during a G4-class geomagnetic storm (stor… ▽ More The GRAPES-3 muon telescope in Ooty, India had claimed detection of a 2 hour (h) high-energy ($\sim$20 GeV) burst of galactic cosmic-rays (GCRs) through a $>$50$σ$ surge in GeV muons, was caused by reconnection of the interplanetary magnetic field (IMF) in the magnetosphere that led to transient weakening of Earth's magnetic shield. This burst had occurred during a G4-class geomagnetic storm (storm) with a delay of $\frac{1}{2}$h relative to the coronal mass ejection (CME) of 22 June 2015 (Mohanty et al., 2016). However, recently a group interpreted the occurrence of the same burst in a subset of 31 neutron monitors (NMs) to have been the result of an anisotropy in interplanetary space (Evenson et al., 2017) in contrast to the claim in (Mohanty et al., 2016). A new analysis of the GRAPES-3 data with a fine 10.6$^{\circ}$ angular segmentation shows the speculation of interplanetary anisotropy to be incorrect, and offers a possible explanation of the NM observations. The observed 28 minutes (min) delay of the burst relative to the CME can be explained by the movement of the reconnection front from the bow shock to the surface of Earth at an average speed of 35 km/s, much lower than the CME speed of 700 km/s. This measurement may provide a more accurate estimate of the start of the storm. △ Less

Submitted 28 March, 2018; originally announced March 2018.

Comments: Accepted for Publication in Physical Review D

Evidence for GeV Cosmic Rays from White Dwarfs in the Local Cosmic Ray Spectra and in the Gamma-ray Emissivity of the Inner Galaxy

Authors: Tuneyoshi Kamae, Shiu-hang Lee, Kazuo Makishima, Shinpei Shibata, Toshikazu Shigeyama

Abstract: Recent observations found that electrons are accelerated to $\sim$10 GeV and emit synchrotron hard X-rays in two magnetic white dwarfs (WDs), also known as cataclysmic variables (CVs). In nova outbursts of WDs, multi-GeV gamma-rays were detected inferring that protons are accelerated to 100 GeV or higher. In recent optical surveys, the WD density is found to be higher near the Sun than in the Gala… ▽ More Recent observations found that electrons are accelerated to $\sim$10 GeV and emit synchrotron hard X-rays in two magnetic white dwarfs (WDs), also known as cataclysmic variables (CVs). In nova outbursts of WDs, multi-GeV gamma-rays were detected inferring that protons are accelerated to 100 GeV or higher. In recent optical surveys, the WD density is found to be higher near the Sun than in the Galactic disk by a factor $\sim$2.5. The cosmic rays (CR) produced by local CVs and novae will accumulate in the local bubble for $10^6$ - $10^7$ yrs. On these findings, we search for CRs from historic CVs and novae in the observed CR spectra. We model the CR spectra at the heliopause as sums of Galactic and local components based on observational data as much as possible. The initial Galactic CR electron and proton spectra are deduced from the gamma-ray emissivity, the local electron spectrum from the hard X-ray spectra at the CVs, and the local proton spectrum inferred by gamma-ray spectrum at novae. These spectral shapes are then expressed in a simple set of polynomial functions of CR energy and regressively fitted until the high-energy ($>$100 GeV) CR spectra near Earth and the Voyager-1 spectra at the heliopause are reproduced. We then extend the modeling to nuclear CR spectra and find that one spectral shape fits all local nuclear CRs and the apparent hardening of the nuclear CR spectra is caused by the roll-down of local nuclear spectra around 100 - 200 GeV. All local CR spectra populate in a limited energy band below 100 - 200 GeV and enhance gamma-ray emissivity below $\sim$10 GeV. Such an enhancement is observed in the inner Galaxy, suggesting the CR fluxes from CVs and novae are substantially higher there. △ Less

Submitted 14 February, 2018; v1 submitted 5 January, 2018; originally announced January 2018.

Comments: 14 pages, 8 figures, 2 tables. Accepted for publication in PASJ (version after proofreading)

Analysis of solar gamma rays and solar neutrons detected on March 7th and September 25th of 2011 by Ground Level Neutron Telescopes, SEDA-FIB and FERMI-LAT

Authors: Y. Muraki, J. F. Valdes-Galicia, L. X. Gonzalez, K. Kamiya, Y. Katayose, K. Koga, H. Matsumoto, S. Masuda, Y. Matsubara, Y. Nagai, M. Ohnishi, S. Ozawa, T. Sako, S. Shibata, M. Takita, Y. Tanaka, H. Tsuchiya, K. Watanabe, J. L. Zhang

Abstract: At the 33rd ICRC, we reported the possible detection of solar gamma rays by a ground level detector and later re-examined this event. On March 7, 2011, the solar neutron telescope (SNT) located at Mt. Sierra Negra, Mexico (4,600 m) observed enhancements of the counting rate from 19:49 to 20:02 UT and from 20:50 to 21:01 UT. The statistical significance was 9.7sigma and 8.5sigma, respectively. This… ▽ More At the 33rd ICRC, we reported the possible detection of solar gamma rays by a ground level detector and later re-examined this event. On March 7, 2011, the solar neutron telescope (SNT) located at Mt. Sierra Negra, Mexico (4,600 m) observed enhancements of the counting rate from 19:49 to 20:02 UT and from 20:50 to 21:01 UT. The statistical significance was 9.7sigma and 8.5sigma, respectively. This paper discusses the possibility of using this mountain detector to detect solar gamma rays. In association with this event, the solar neutron detector SEDA-FIB onboard the International Space Station has also detected solar neutrons with a statistical significance of 7.5sigma. The FERMI-LAT detector also observed high-energy gamma rays from this flare with a statistical significance of 6.7sigma. We thus attempted to make a unified model to explain this data. In this paper, we report on another candidate for solar gamma rays detected on September 25th, 2011 by the SNT located in Tibet (4,300 m) from 04:37 to 04:47 UT with a statistical significance of 8.0sigma (by the Li-Ma method). △ Less

Submitted 27 June, 2017; originally announced June 2017.

Comments: 15 pages

Magnetic and magnetotransport properties of ThCr2Si2-type Ce2O2Bi composed of conducting Bi2- square net and magnetic Ce-O layer

Authors: Shunsuke Shibata, Ryosuke Sei, Tomoteru Fukumura, Tetsuya Hasegawa

Abstract: ThCr2Si2-type Ce2O2Bi epitaxial thin films were grown by recently developed multilayer solid phase epitaxy. The ionic state of Ce was confirmed to be 3+ by x-ray photoelectron spectroscopy, corresponding to the electronic configuration of [Xe]4f1. Electrical resistivity showed the nonmonotonic temperature dependence with a sharp resistivity maximum, concomitant with a magnetization kink at 10 K, s… ▽ More ThCr2Si2-type Ce2O2Bi epitaxial thin films were grown by recently developed multilayer solid phase epitaxy. The ionic state of Ce was confirmed to be 3+ by x-ray photoelectron spectroscopy, corresponding to the electronic configuration of [Xe]4f1. Electrical resistivity showed the nonmonotonic temperature dependence with a sharp resistivity maximum, concomitant with a magnetization kink at 10 K, suggesting antiferromagnetic transition. In addition, magnetoresistance showed a large angular-dependent magnetoresistance. These results imply that carrier transport in the Bi2- square net could be influenced by magnetic ordering in the Ce-O layer owing to its unique layered structure [Bi2-/(Ce2O2)2+], particularly in the form of epitaxial thin films. △ Less

Submitted 27 May, 2017; originally announced May 2017.

Comments: 7 pages, 4 figures

Journal ref: Appl. Phys. Lett. 110, 192410 (2017)

Magnetar Broadband X-ray Spectra Correlated with Magnetic Fields: Suzaku Archive of SGRs and AXPs Combined with NuSTAR, Swift, and RXTE

Authors: Teruaki Enoto, Shinpei Shibata, Takao Kitaguchi, Yudai Suwa, Takahiko Uchide, Hiroyuki Nishioka, Shota Kisaka, Toshio Nakano, Hiroaki Murakami, Kazuo Makishima

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

Submitted 23 April, 2017; originally announced April 2017.

Comments: 22 pages, 18 figures, 7 tables, accepted for publication in The Astrophysical Journal Supplement

Constraints on Pulsed Emission Model for Repeating FRB 121102

Authors: Shota Kisaka, Teruaki Enoto, Shinpei Shibata

Abstract: Recent localization of the repeating Fast Radio Burst (FRB) 121102 revealed the distance of its host galaxy and luminosities of the bursts. We investigated constraints on the young neutron star (NS) model, that (a) the FRB intrinsic luminosity is supported by the spin-down energy, and (b) the FRB duration is shorter than the NS rotation period. In the case of a circular cone emission geometry, con… ▽ More Recent localization of the repeating Fast Radio Burst (FRB) 121102 revealed the distance of its host galaxy and luminosities of the bursts. We investigated constraints on the young neutron star (NS) model, that (a) the FRB intrinsic luminosity is supported by the spin-down energy, and (b) the FRB duration is shorter than the NS rotation period. In the case of a circular cone emission geometry, conditions (a) and (b) determine the NS parameters within very small ranges, compared with that from only condition (a) discussed in previous works. Anisotropy of the pulsed emission does not affect the area of the allowed parameter region by virtue of condition (b). The determined parameters are consistent with those independently limited by the properties of the possible persistent radio counterpart and the circumburst environments such as surrounding materials. Since the NS in the allowed parameter region is older than the spin-down timescale, the hypothetical GRP-like model expects a rapid radio flux decay of $\lesssim1$ Jy within a few years as the spin-down luminosity decreases. The continuous monitoring will give a hint of discrimination of the models. If no flux evolution will be seen, we need to consider an alternative model, e.g., the magnetically powered flare. △ Less

Submitted 17 August, 2017; v1 submitted 9 February, 2017; originally announced February 2017.

Comments: 12 pages, 1 figure. Accepted for publication in PASJ

A photon observed by the SEDA-FIB detector on the arrival of the gravitational wave

Authors: Y. Muraki, K. Kamiya, K. Koga, H. Matsumoto, S. Shibata

Abstract: We reported a photon detection at the arrival time of the gravitational wave on December 26th, 2015. According to the LIGO-Virgo collaboration, the gravitational wave was produced by the merging process of the two black holes. The merged time was determined at 03:38:53.6 UT on December 26th. At 03:38:54.05 GPS time, one of the detectors of SEDA-FIB on board the International Space Station (ISS) de… ▽ More We reported a photon detection at the arrival time of the gravitational wave on December 26th, 2015. According to the LIGO-Virgo collaboration, the gravitational wave was produced by the merging process of the two black holes. The merged time was determined at 03:38:53.6 UT on December 26th. At 03:38:54.05 GPS time, one of the detectors of SEDA-FIB on board the International Space Station (ISS) detected a photon, arriving from the direction of Corona Australis near the Galactic Bulge (GB). The energy of the photon was about (35+-7) MeV. A 3.0 sigma level detection of the photon is discussed. △ Less

Submitted 29 November, 2016; originally announced November 2016.

Comments: 9 pages, 5 figures, presented talk at JSPS in September 016, Miyazaki University

X-ray and Rotational Luminosity Correlation and Magnetic Heating of the Radio Pulsars

Authors: S. Shibata, E. Watanabe, Y. Yatsu, T. Enoto, A. Bamba

Abstract: Previous works have suggested a correlation between the X-ray luminosity Lx and the rotational luminosity Lrot of the radio pulsars.However, none of the obtained regression lines are statistically acceptable due to large scatters. We construct a statistical model which has an intrinsic Lx-Lrot relation and reproduces the observed Lx distribution about it by using a Monte Carlo simulator, which t… ▽ More Previous works have suggested a correlation between the X-ray luminosity Lx and the rotational luminosity Lrot of the radio pulsars.However, none of the obtained regression lines are statistically acceptable due to large scatters. We construct a statistical model which has an intrinsic Lx-Lrot relation and reproduces the observed Lx distribution about it by using a Monte Carlo simulator, which takes into account the effects obscuring the intrinsic relation,i.e., the anisotropy of radiation, additional heating, uncertainty in distance and detection limit of the instruments. From the ATNF pulsar catalog we collect 57 `ordinary radio pulsars' with significant detection and 42 with upper limits.The sample does not include the high-magnetic field pulsars (>10^{13} G), which are separately analyzed. We obtain a statistically acceptable relation Lx (0.5 - 10 keV)= 10^{31.69} (Lrot / L_0)^{c_1} with c_1 = 1.03 \pm 0.27 and L_0 =10^{35.38}. The distribution about the obtained Lx-Lrot relation is reproduced well by the simulator. Pulsars with abnormally high Lx fall into two types: one is the soft gamma-ray pulsars, and the other is thermally bright pulsars in comparison with the standard cooling curve. On the other hand, pulsars showing low Lx are found to have dim pulsar wind nebulae. We argue that there is an unknown mechanism that governs both the magnetospheric emission and the pulsar wind nebulae, and it might involve the production rate of electron-positron pairs. The high-field pulsars form a distinctive population other than the ordinary pulsars due to their excess luminosities. △ Less

Submitted 5 October, 2016; originally announced October 2016.

Comments: 17 pages, 11 figures, to be published in ApJ

SKA-Japan Pulsar Science with the Square Kilometre Array

Authors: Keitaro Takahashi, Takahiro Aoki, Kengo Iwata, Osamu Kameya, Hiroki Kumamoto, Sachiko Kuroyanagi, Ryo Mikami, Atsushi Naruko, Hiroshi Ohno, Shinpei Shibata, Toshio Terasawa, Naoyuki Yonemaru, Chulmoon Yoo

Abstract: The Square Kilometre Array will revolutionize pulsar studies with its wide field-of-view, wide-band observation and high sensitivity, increasing the number of observable pulsars by more than an order of magnitude. Pulsars are of interest not only for the study of neutron stars themselves but for their usage as tools for probing fundamental physics such as general relativity, gravitational waves an… ▽ More The Square Kilometre Array will revolutionize pulsar studies with its wide field-of-view, wide-band observation and high sensitivity, increasing the number of observable pulsars by more than an order of magnitude. Pulsars are of interest not only for the study of neutron stars themselves but for their usage as tools for probing fundamental physics such as general relativity, gravitational waves and nuclear interaction. In this article, we summarize the activity and interests of SKA-Japan Pulsar Science Working Group, focusing on an investigation of modified gravity theory with the supermassive black hole in the Galactic Centre, gravitational-wave detection from cosmic strings and binary supermassive black holes, a study of the physical state of plasma close to pulsars using giant radio pulses and determination of magnetic field structure of Galaxy with pulsar pairs. △ Less

Submitted 7 March, 2016; originally announced March 2016.

Comments: 13 pages, 5 figures

The Search for >35 MeV Neutrons from the June 3, 2012 Impulsive Flare

Authors: K. Koga, S. Masuda, H. Matsumoto, Y. Muraki, T. Obara O. Okudaira, S. Shibata, T. Yamamoto, T. Goka

Abstract: We analyzed a highly impulsive solar flare observed on June 3, 2012. In association with this flare, emissions of hard X-rays, high-energy gamma rays, and neutrons were detected by the detectors onboard the FERMI, RHESSI satellites and the International Space Station. We compared those results with the pictures taken by the UV telescope onboard the Solar Dynamics Observatory satellite and found th… ▽ More We analyzed a highly impulsive solar flare observed on June 3, 2012. In association with this flare, emissions of hard X-rays, high-energy gamma rays, and neutrons were detected by the detectors onboard the FERMI, RHESSI satellites and the International Space Station. We compared those results with the pictures taken by the UV telescope onboard the Solar Dynamics Observatory satellite and found the crossing structure of two magnetic ropes at two positions on the solar surface almost at the same time. High-energy gamma rays were detected by the Fermi Large Area Telescope satellite, implying that the impulsive flare was one of a major source of proton acceleration processes on the solar surface. At the beginning of research, impulsive solar flares were considered to be the main source of particle acceleration processes; our current observations have confirmed this hypothesis. △ Less

Submitted 20 August, 2015; originally announced August 2015.

Comments: 8 pages, 9 figures, Proceeding of the 34th International Cosmic Ray Conference in Hague in August,2015

Report number: PoS #100

Results of a Measurement of Solar Neutrons Emitted on March 5, 2012 using a Fiber-type neutron monitor onboard the SEDA-AP attached to the ISS

Authors: K. Koga, H. Matsumoto, O. Okudaira, T. Goka, T. Obara, S. Masuda, Y. Muraki, S. Shibata, T. Yamamoto

Abstract: The solar neutron detector SEDA-FIB onboard the International Space Station (ISS) has detected several events from the solar direction associated with three large solar flares observed on March 5th (X1.1), 7th (X5.4), and 9th (M6.3) of 2012. In this study, we present the time profiles of those neutrons and discuss the physics that may be related to a possible acceleration scenario for ions over th… ▽ More The solar neutron detector SEDA-FIB onboard the International Space Station (ISS) has detected several events from the solar direction associated with three large solar flares observed on March 5th (X1.1), 7th (X5.4), and 9th (M6.3) of 2012. In this study, we present the time profiles of those neutrons and discuss the physics that may be related to a possible acceleration scenario for ions over the solar surface. We compare our data with the dynamical pictures of the flares obtained by the ultra-violet telescope of the space-based Solar Dynamics Observatory. △ Less

Submitted 20 August, 2015; originally announced August 2015.

Comments: 8 pages, 8 figures, and 1 table Proceeding of the 34th International Cosmic Ray Conference in Hague in August, 2015

Report number: PoS 068

Journal ref: Solar Physics (2017) 292:115

Simultaneous Observation of Solar Neutrons from the ISS and High Mountain Observatories in association with a flare on July 8, 2014

Authors: Y. Muraki, D. Lopez, K. Koga, F. Kakimoto, T. Goka, L. X. Gonzalez, S. Masuda, Y. Matsubara, H. Matsumoto, P. Miranda, O. Okudaira, T. Obara, J. Salinas, T. Sako, S. Shibata, R. Ticona, Y. Tsunesada, J. F. Valdes-Galicia, K. Watanabe, T. Yamamoto

Abstract: An M6.5-class flare was observed at N12E56 of the solar surface at 16:06 UT on July 8, 2014. In association with this flare, solar neutron detectors located on two high mountains, Mt. Sierra Negra and Chacaltaya and at the space station observed enhancements in the neutral channel. The authors analysed these data and a possible scenario of enhancements produced by high-energy protons and neutrons… ▽ More An M6.5-class flare was observed at N12E56 of the solar surface at 16:06 UT on July 8, 2014. In association with this flare, solar neutron detectors located on two high mountains, Mt. Sierra Negra and Chacaltaya and at the space station observed enhancements in the neutral channel. The authors analysed these data and a possible scenario of enhancements produced by high-energy protons and neutrons is proposed, using the data from continuous observation of a solar surface by the ultraviolet telescope onboard the Solar Dynamical Observatory (SDO). △ Less

Submitted 20 August, 2015; originally announced August 2015.

Comments: 8 pages, 7 figures, and 2 Tables. Proceeding of the 34th International Cosmic Ray Conference in Hague in August, 2015

Report number: PoS 050

Journal ref: Solar Physics 291 (2016) 1241-1265

Time-dependent Multi-group Multidimensional Relativistic Radiative Transfer Code Based On Spherical Harmonic Discrete Ordinate Method

Authors: Nozomu Tominaga, Sanshiro Shibata, Sergei I. Blinnikov

Abstract: We develop a time-dependent multi-group multidimensional relativistic radiative transfer code, which is required to numerically investigate radiation from relativistic fluids involved in, e.g., gamma-ray bursts and active galactic nuclei. The code is based on the spherical harmonic discrete ordinate method (SHDOM) that evaluates a source function including anisotropic scattering in spherical harmo… ▽ More We develop a time-dependent multi-group multidimensional relativistic radiative transfer code, which is required to numerically investigate radiation from relativistic fluids involved in, e.g., gamma-ray bursts and active galactic nuclei. The code is based on the spherical harmonic discrete ordinate method (SHDOM) that evaluates a source function including anisotropic scattering in spherical harmonics and implicitly solves the static radiative transfer equation with a ray tracing in discrete ordinates. We implement treatments of time dependence, multi-frequency bins, Lorentz transformation, and elastic Thomson and inelastic Compton scattering to the publicly available SHDOM code. Our code adopts a mixed frame approach; the source function is evaluated in the comoving frame whereas the radiative transfer equation is solved in the laboratory frame. This implementation is validated with various test problems and comparisons with results of a relativistic Monte Carlo code. These validations confirm that the code correctly calculates intensity and its evolution in the computational domain. The code enables us to obtain an Eddington tensor that relates first and third moments of intensity (energy density and radiation pressure) and is frequently used as a closure relation in radiation hydrodynamics calculations. △ Less

Submitted 17 July, 2015; originally announced July 2015.

Comments: 13 pages, 12 figures. Accepted for publication in the Astrophysical Journal Supplement Series

Nuclear Composition of Magnetized GRB Jets

Authors: Sanshiro Shibata, Nozomu Tominaga

Abstract: We investigate the fraction of metal nuclei in the relativistic jets of gamma-ray bursts associated with core-collapse supernovae. We simulate the fallback in jet-induced explosions with two-dimensional relativistic hydrodynamics calculations and the jet acceleration with steady, radial, relativistic magnetohydrodynamics calculations, and derive detail nuclear composition of the jet by postprocess… ▽ More We investigate the fraction of metal nuclei in the relativistic jets of gamma-ray bursts associated with core-collapse supernovae. We simulate the fallback in jet-induced explosions with two-dimensional relativistic hydrodynamics calculations and the jet acceleration with steady, radial, relativistic magnetohydrodynamics calculations, and derive detail nuclear composition of the jet by postprocessing calculation. We found that if the temperature at the jet launch site is above $4.7\times 10^9$K, quasi-statistical equilibrium (QSE) is established and heavy nuclei are dissociated to light particles such as $^4$He during the acceleration of the jets. The criterion for the survival of metal nuclei is written in terms of the isotropic jet luminosity as $L_{\rm j}^{\rm iso} \lesssim 3.9\times 10^{50}(R_{\rm i}/10^7{\rm cm})^2 (1+σ_{\rm i})~{\rm erg~s^{-1}}$, where $R_{\rm i}$ and $σ_{\rm i}$ are the initial radius of the jets and the initial magnetization parameter, respectively. If the jet is initially dominated by radiation field (i.e., $σ_{\rm i} \ll 1$) and the isotropic luminosity is relatively high ($L_{\rm j}^{\rm iso}\gtrsim 4\times 10^{52}~{\rm erg s^{-1}}$), the metal nuclei cannot survive in the jet. On the other hand, if the jet is mainly accelerated by magnetic field (i.e., $σ_{\rm i} \gg1$), metal nuclei initially contained in the jet can survive without serious dissociation even for the case of high luminosity jet. If the jet contains metal nuclei, the dominant nuclei are $^{28}$Si, $^{16}$O, and $^{32}$S and the mean mass number can be $\langle$A$\rangle\sim25$. △ Less

Submitted 12 March, 2015; originally announced March 2015.

Comments: 9 pages, 6 figures, accepted for publication in the Publications of the Astronomical Society of Japan (PASJ)