Showing 1–50 of 89 results for author: Fujisawa, K

Status of Xtend telescope onboard X-Ray Imaging and Spectroscopy Mission (XRISM)

Authors: Koji Mori, Hiroshi Tomida, Hiroshi Nakajima, Takashi Okajima, Hirofumi Noda, Hiroyuki Uchida, Hiromasa Suzuki, Shogo Benjamin Kobayashi, Tomokage Yoneyama, Kouichi Hagino, Kumiko Nobukawa, Takaaki Tanaka, Hiroshi Murakami, Hideki Uchiyama, Masayoshi Nobukawa, Hironori Matsumoto, Takeshi Tsuru, Makoto Yamauchi, Isamu Hatsukade, Hirokazu Odaka, Takayoshi Kohmura, Kazutaka Yamaoka, Manabu Ishida, Yoshitomo Maeda, Takayuki Hayashi , et al. (38 additional authors not shown)

Abstract: Xtend is one of the two telescopes onboard the X-ray imaging and spectroscopy mission (XRISM), which was launched on September 7th, 2023. Xtend comprises the Soft X-ray Imager (SXI), an X-ray CCD camera, and the X-ray Mirror Assembly (XMA), a thin-foil-nested conically approximated Wolter-I optics. A large field of view of $38^{\prime}\times38^{\prime}$ over the energy range from 0.4 to 13 keV is… ▽ More Xtend is one of the two telescopes onboard the X-ray imaging and spectroscopy mission (XRISM), which was launched on September 7th, 2023. Xtend comprises the Soft X-ray Imager (SXI), an X-ray CCD camera, and the X-ray Mirror Assembly (XMA), a thin-foil-nested conically approximated Wolter-I optics. A large field of view of $38^{\prime}\times38^{\prime}$ over the energy range from 0.4 to 13 keV is realized by the combination of the SXI and XMA with a focal length of 5.6 m. The SXI employs four P-channel, back-illuminated type CCDs with a thick depletion layer of 200 $μ$m. The four CCD chips are arranged in a 2$\times$2 grid and cooled down to $-110$ $^{\circ}$C with a single-stage Stirling cooler. Before the launch of XRISM, we conducted a month-long spacecraft thermal vacuum test. The performance verification of the SXI was successfully carried out in a course of multiple thermal cycles of the spacecraft. About a month after the launch of XRISM, the SXI was carefully activated and the soundness of its functionality was checked by a step-by-step process. Commissioning observations followed the initial operation. We here present pre- and post-launch results verifying the Xtend performance. All the in-orbit performances are consistent with those measured on ground and satisfy the mission requirement. Extensive calibration studies are ongoing. △ Less

Submitted 28 June, 2024; originally announced June 2024.

Comments: 10 pages, 8 figures. Proceedings of SPIE Astronomical Telescopes and Instrumentation 2024

Initial operations of the Soft X-ray Imager onboard XRISM

Authors: Hiromasa Suzuki, Tomokage Yoneyama, Shogo B. Kobayashi, Hirofumi Noda, Hiroyuki Uchida, Kumiko K. Nobukawa, Kouichi Hagino, Koji Mori, Hiroshi Tomida, Hiroshi Nakajima, Takaaki Tanaka, Hiroshi Murakami, Hideki Uchiyama, Masayoshi Nobukawa, Yoshiaki Kanemaru, Yoshinori Otsuka, Haruhiko Yokosu, Wakana Yonemaru, Hanako Nakano, Kazuhiro Ichikawa, Reo Takemoto, Tsukasa Matsushima, Marina Yoshimoto, Mio Aoyagi, Kohei Shima , et al. (30 additional authors not shown)

Abstract: XRISM (X-Ray Imaging and Spectroscopy Mission) is an astronomical satellite with the capability of high-resolution spectroscopy with the X-ray microcalorimeter, Resolve, and wide field-of-view imaging with the CCD camera, Xtend. The Xtend consists of the mirror assembly (XMA: X-ray Mirror Assembly) and detector (SXI: Soft X-ray Imager). The components of SXI include CCDs, analog and digital electr… ▽ More XRISM (X-Ray Imaging and Spectroscopy Mission) is an astronomical satellite with the capability of high-resolution spectroscopy with the X-ray microcalorimeter, Resolve, and wide field-of-view imaging with the CCD camera, Xtend. The Xtend consists of the mirror assembly (XMA: X-ray Mirror Assembly) and detector (SXI: Soft X-ray Imager). The components of SXI include CCDs, analog and digital electronics, and a mechanical cooler. After the successful launch on September 6th, 2023 (UT) and subsequent critical operations, the mission instruments were turned on and set up. The CCDs have been kept at the designed operating temperature of $-110^\circ$C ~after the electronics and cooling system were successfully set up. During the initial operation phase, which continued for more than a month after the critical operations, we verified the observation procedure, stability of the cooling system, all the observation options with different imaging areas and/or timing resolutions, and operations for protection against South Atlantic Anomaly. We optimized the operation procedure and observation parameters including the cooler settings, imaging areas for the specific modes with higher timing resolutions, and event selection algorithm. We summarize our policy and procedure of the initial operations for SXI. We also report on a couple of issues we faced during the initial operations and lessons learned from them. △ Less

Submitted 28 June, 2024; originally announced June 2024.

Comments: 10 pages, 8 figures, Proceedings of SPIE Astronomical Telescopes + Instrumentation 2024

The optically thick rotating magnetic wind from a massive white dwarf merger product -- II. axisymmetric magnetohydrodynamic simulations

Authors: Yici Zhong, Kazumi Kashiyama, Shinsuke Takasao, Toshikazu Shigeyama, Kotaro Fujisawa

Abstract: We numerically construct a series of axisymmetric rotating magnetic wind solutions, aiming at exploring the observation properties of massive white dwarf (WD) merger remnants with a strong magnetic field, a fast spin, and an intense mass loss, as inferred for WD J005311. We investigate the magnetospheric structure and the resultant spin-down torque exerted to the merger remnant with respect to the… ▽ More We numerically construct a series of axisymmetric rotating magnetic wind solutions, aiming at exploring the observation properties of massive white dwarf (WD) merger remnants with a strong magnetic field, a fast spin, and an intense mass loss, as inferred for WD J005311. We investigate the magnetospheric structure and the resultant spin-down torque exerted to the merger remnant with respect to the surface magnetic flux $Φ_*$, spin angular frequency $Ω_*$ and the mass loss rate $\dot M$. We confirm that the wind properties for $σ\equiv Φ^2_* Ω_*^2/\dot M v_\mathrm{esc}^3 \gtrsim 1$ significantly deviate from those of the spherical Parker wind, where $v_\mathrm{esc}$ is the escape velocity at stellar surface. For such a rotating magnetic wind sequence, we find: (i) quasi-periodic mass eruption triggered by magnetic reconnection along with the equatorial plane (ii) a scaling relation for the spin-down torque $T \approx (1/2) \times \dot{M} Ω_* R^2_* σ^{1/4}$. We apply our results to discuss the spin-down evolution and wind anisotropy of massive WD merger remnants, the latter of which could be probed by a successive observation of WD J005311 using Chandra. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: 15 pages, 6 figures, comments are welcome

A kinematic study of the disc-outflow system around a high-mass protostar G59.783+0.065 probed by methanol and water masers

Authors: M. Nakamura, K. Motogi, H. Nakamura, Y. Yonekura, K. Fujisawa

Abstract: Class II CH3OH masers are used as a convenient tracer of disc-like structures in high-mass star formation. However, more than half of them show a complex distribution in Very Long Baseline Interferometry (VLBI) maps. The origin of such a complex distribution is still unknown. We conducted VLBI monitoring observations to unveil the origin of a complex class II CH3OH maser in the high-mass star-form… ▽ More Class II CH3OH masers are used as a convenient tracer of disc-like structures in high-mass star formation. However, more than half of them show a complex distribution in Very Long Baseline Interferometry (VLBI) maps. The origin of such a complex distribution is still unknown. We conducted VLBI monitoring observations to unveil the origin of a complex class II CH3OH maser in the high-mass star-forming region G59.783+0.065. We observed the CH3OH maser at 6.7 GHz and the H2O maser at 22 GHz to probe detailed circumstellar kinematics and structures by the Japanese VLBI network and the VLBI Exploration of Radio Astrometry. We found similar bipolar distributions in both masers, specifically two clusters located 2000 au apart along the East-West direction. We detected a linear distribution of CH3OH masers in the Western cluster. A position-velocity diagram shows that the Western CH3OH masers trace a rotating disc-wind or infalling component inside an edge-on disc-like structure. In contrast to the simple bipolar expanding motions of the H2O masers, the CH3OH masers exhibited complex motions despite their spatial coincidence. Some of the Eastern CH3OH masers showed bipolar expansions similar to the H2O masers, while others displayed random or even inward motions. Such complex kinematics and their close association with the H2O maser could occur at the boundary between outflow and inflow. We suggest that the complex distribution of class II CH3OH masers, like G59.783+0.065 arises from several distinct circumstellar structures that simultaneously achieve maser excitation. △ Less

Submitted 4 October, 2023; originally announced October 2023.

Comments: Accepted to MNRAS, 22 pages, 17 figures, 13 tables

Quantum algorithm for the Vlasov simulation of the large-scale structure formation with massive neutrinos

Authors: Koichi Miyamoto, Soichiro Yamazaki, Fumio Uchida, Kotaro Fujisawa, Naoki Yoshida

Abstract: Investigating the cosmological implication of the fact that neutrino has finite mass is of importance for fundamental physics. In particular, massive neutrino affects the formation of the large-scale structure (LSS) of the universe, and, conversely, observations of the LSS can give constraints on the neutrino mass. Numerical simulations of the LSS formation including massive neutrino along with co… ▽ More Investigating the cosmological implication of the fact that neutrino has finite mass is of importance for fundamental physics. In particular, massive neutrino affects the formation of the large-scale structure (LSS) of the universe, and, conversely, observations of the LSS can give constraints on the neutrino mass. Numerical simulations of the LSS formation including massive neutrino along with conventional cold dark matter is thus an important task. For this, calculating the neutrino distribution in the phase space by solving the Vlasov equation is a suitable approach, but it requires solving the PDE in the $(6+1)$-dimensional space and is thus computationally demanding: Configuring $n_\mathrm{gr}$ grid points in each coordinate and $n_t$ time grid points leads to $O(n_\mathrm{gr}^6)$ memory space and $O(n_tn_\mathrm{gr}^6)$ queries to the coefficients in the discretized PDE. We propose a quantum algorithm for this task. Linearizing the Vlasov equation by neglecting the relatively weak self-gravity of the neutrino, we perform the Hamiltonian simulation to produce quantum states that encode the phase space distribution of neutrino. We also propose a way to extract the power spectrum of the neutrino density perturbations as classical data from the quantum state by quantum amplitude estimation with accuracy $ε$ and query complexity of order $\widetilde{O}((n_\mathrm{gr} + n_t)/ε)$. Our method also reduces the space complexity to $O(\mathrm{polylog}(n_\mathrm{gr}/ε))$ in terms of the qubit number, while using quantum random access memories with $O(n_\mathrm{gr}^3)$ entries. As far as we know, this is the first quantum algorithm for the LSS simulation that outputs the quantity of practical interest with guaranteed accuracy. △ Less

Submitted 26 February, 2024; v1 submitted 3 October, 2023; originally announced October 2023.

Report number: RESCEU-22/23

Journal ref: Phys. Rev. Research 6, 013200 (2024)

A dynamical model for IRAS 00500+6713: the remnant of a type Iax supernova SN 1181 hosting a double degenerate merger product WD J005311

Authors: Takatoshi Ko, Hiromasa Suzuki, Kazumi Kashiyama, Hiroyuki Uchida, Takaaki Tanaka, Daichi Tsuna, Kotaro Fujisawa, Aya Bamba, Toshikazu Shigeyama

Abstract: IRAS 00500+6713 is a hypothesized remnant of a type Iax supernova SN 1181. Multi-wavelength observations have revealed its complicated morphology; a dusty infrared ring is sandwiched by the inner and outer X-ray nebulae. We analyze the archival X-ray data taken by XMM-Newton and Chandra to constrain the {angular radius}, mass, and metal abundance of the X-ray nebulae, and construct a theoretical m… ▽ More IRAS 00500+6713 is a hypothesized remnant of a type Iax supernova SN 1181. Multi-wavelength observations have revealed its complicated morphology; a dusty infrared ring is sandwiched by the inner and outer X-ray nebulae. We analyze the archival X-ray data taken by XMM-Newton and Chandra to constrain the {angular radius}, mass, and metal abundance of the X-ray nebulae, and construct a theoretical model describing the dynamical evolution of IRAS 00500+6713, including the effects of the interaction between the SN ejecta and the intense wind enriched with carbon burning ashes from the central white dwarf (WD) J005311. We show that the inner X-ray nebula corresponds to the wind termination shock while the outer X-ray nebula to the shocked interface between the SN ejecta and the interstellar matter. The observed X-ray properties can be explained by our model with an {ejecta kinetic} energy of $E_\mathrm{ej} = (0.77 \mbox{--} 1.1)\times 10^{48}$~erg, an ejecta mass of $M_\mathrm{ej} = 0.18\mbox{--}0.53~M_\odot$, if the currently observed wind from WD J005311 started to blow $t_\mathrm{w} \gtrsim 810$ yr after the explosion, i.e., approximately after A.D. 1990. The inferred SN properties are compatible with those of Type Iax SNe and the timing of the wind launch may correspond to the Kelvin-Helmholtz contraction of the oxygen-neon core of WD J005311 that triggered a surface carbon burning. Our analysis supports that IRAS 00500+6713 is the remnant of SN Iax 1181 produced by a double degenerate merger of oxygen-neon and carbon-oxygen WDs, and WD J005311 is the surviving merger product. △ Less

Submitted 26 May, 2024; v1 submitted 28 April, 2023; originally announced April 2023.

Comments: 24 pages, 13 figures, 4 tables, accepted by ApJ

Report number: RESCEU-10/23

Quantum algorithm for collisionless Boltzmann simulation of self-gravitating systems

Authors: Soichiro Yamazaki, Fumio Uchida, Kotaro Fujisawa, Koichi Miyamoto, Naoki Yoshida

Abstract: The collisionless Boltzmann equation (CBE) is a fundamental equation that governs the dynamics of a broad range of astrophysical systems from space plasma to star clusters and galaxies. It is computationally expensive to integrate the CBE directly in a multi-dimensional phase space, and thus the applications to realistic astrophysical problems have been limited so far. Recently, Todorova & Steijl… ▽ More The collisionless Boltzmann equation (CBE) is a fundamental equation that governs the dynamics of a broad range of astrophysical systems from space plasma to star clusters and galaxies. It is computationally expensive to integrate the CBE directly in a multi-dimensional phase space, and thus the applications to realistic astrophysical problems have been limited so far. Recently, Todorova & Steijl (2020) proposed an efficient quantum algorithm to solve the CBE with significantly reduced computational complexity. We extend the algorithm to perform quantum simulations of self-gravitating systems, incorporating the method to calculate gravity with the major Fourier modes of the density distribution extracted from the solution-encoding quantum state. Our method improves the dependency of time and space complexities on Nv , the number of grid points in each velocity coordinate, compared to the classical simulation methods. We then conduct some numerical demonstrations of our method. We first run a 1+1 dimensional test calculation of free streaming motion on 64*64 grids using 13 simulated qubits and validate our method. We then perform simulations of Jeans collapse, and compare the result with analytic and linear theory calculations. It will thus allow us to perform large-scale CBE simulations on future quantum computers. △ Less

Submitted 23 May, 2024; v1 submitted 29 March, 2023; originally announced March 2023.

Comments: 12 pages, 10 figures

Xtend, the Soft X-ray Imaging Telescope for the X-ray Imaging and Spectroscopy Mission (XRISM)

Authors: Koji Mori, Hiroshi Tomida, Hiroshi Nakajima, Takashi Okajima, Hirofumi Noda, Takaaki Tanaka, Hiroyuki Uchida, Kouichi Hagino, Shogo Benjamin Kobayashi, Hiromasa Suzuki, Tessei Yoshida, Hiroshi Murakami, Hideki Uchiyama, Masayoshi Nobukawa, Kumiko Nobukawa, Tomokage Yoneyama, Hironori Matsumoto, Takeshi Tsuru, Makoto Yamauchi, Isamu Hatsukade, Manabu Ishida, Yoshitomo Maeda, Takayuki Hayashi, Keisuke Tamura, Rozenn Boissay-Malaquin , et al. (30 additional authors not shown)

Abstract: Xtend is a soft X-ray imaging telescope developed for the X-Ray Imaging and Spectroscopy Mission (XRISM). XRISM is scheduled to be launched in the Japanese fiscal year 2022. Xtend consists of the Soft X-ray Imager (SXI), an X-ray CCD camera, and the X-ray Mirror Assembly (XMA), a thin-foil-nested conically approximated Wolter-I optics. The SXI uses the P-channel, back-illuminated type CCD with an… ▽ More Xtend is a soft X-ray imaging telescope developed for the X-Ray Imaging and Spectroscopy Mission (XRISM). XRISM is scheduled to be launched in the Japanese fiscal year 2022. Xtend consists of the Soft X-ray Imager (SXI), an X-ray CCD camera, and the X-ray Mirror Assembly (XMA), a thin-foil-nested conically approximated Wolter-I optics. The SXI uses the P-channel, back-illuminated type CCD with an imaging area size of 31 mm on a side. The four CCD chips are arranged in a 2$\times$2 grid and can be cooled down to $-120$ $^{\circ}$C with a single-stage Stirling cooler. The XMA nests thin aluminum foils coated with gold in a confocal way with an outer diameter of 45~cm. A pre-collimator is installed in front of the X-ray mirror for the reduction of the stray light. Combining the SXI and XMA with a focal length of 5.6m, a field of view of $38^{\prime}\times38^{\prime}$ over the energy range from 0.4 to 13 keV is realized. We have completed the fabrication of the flight model of both SXI and XMA. The performance verification has been successfully conducted in a series of sub-system level tests. We also carried out on-ground calibration measurements and the data analysis is ongoing. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: 11 pages, 13 figures. Proceedings of SPIE Astronomical Telescopes and Instrumentation 2022

Journal ref: Proc. of SPIE vol. 12181, 121811T (2022)

Evolutionary deformation toward the elastic limit by a magnetic field confined in the neutron--star crust

Authors: Yasufumi Kojima, Shota Kisaka, Kotaro Fujisawa

Abstract: Occasional energetic outbursts and anomalous X-ray luminosities are expected to be powered by the strong magnetic field in a neutron star. For a very strong magnetic field, elastic deformation becomes excessively large such that it leads to crustal failure. We studied the evolutionary process driven by the Hall drift for a magnetic field confined inside the crust. Assuming that the elastic force a… ▽ More Occasional energetic outbursts and anomalous X-ray luminosities are expected to be powered by the strong magnetic field in a neutron star. For a very strong magnetic field, elastic deformation becomes excessively large such that it leads to crustal failure. We studied the evolutionary process driven by the Hall drift for a magnetic field confined inside the crust. Assuming that the elastic force acts against the Lorentz force, we examined the duration of the elastic regime and maximum elastic energy stored before the critical state. The breakup time was longer than that required for extending the field to the exterior, because the tangential components of the Lorentz force vanished in the fragile surface region. The conversion of large magnetic energy, confined to the interior, into Joule heat is considered to explain the power for central compact objects. This process can function without reaching its elastic limit, unless the magnetic energy exceeds $2\times 10^{47}$ erg, which requires an average field strength of $2\times10^{15}$ G. Thus, the strong magnetic field hidden in the crust is unlikely to cause outbursts. Furthermore, the magnetic field configuration can discriminate between central compact objects and magnetars. △ Less

Submitted 3 March, 2023; originally announced March 2023.

Comments: 13 pages, 7 figures

Strong toroidal magnetic fields sustained by the elastic crust in a neutron star

Authors: Kotaro Fujisawa, Yasufumi Kojima, Shota Kisaka

Abstract: We investigate new solutions for magnetized neutron stars with a barotropic core in magnetohydrodynamic (MHD) equilibrium and a magneto-elastic crust, which was neglected by previous studies concerning stars in MHD equilibrium. The Lorentz force of the barotropic star is purely irrotational and the structures of magnetic fields are constrained. By contrast, a solenoidal component of the Lorentz fo… ▽ More We investigate new solutions for magnetized neutron stars with a barotropic core in magnetohydrodynamic (MHD) equilibrium and a magneto-elastic crust, which was neglected by previous studies concerning stars in MHD equilibrium. The Lorentz force of the barotropic star is purely irrotational and the structures of magnetic fields are constrained. By contrast, a solenoidal component of the Lorentz force exists in the elastic crust and the structures of the magnetic fields are less restricted. We find that the minor solenoidal component in the elastic crust is important for sustaining the strong magnetic field in the core. Unlike previous studies, the toroidal magnetic field exists in the entire region of the core, and we obtain equilibrium states with large toroidal magnetic fields, where the toroidal magnetic energy is larger than the poloidal magnetic energy. The elastic force of the crust sustains an order of $10^{15}~\mathrm{G}$ toroidal magnetic field in the core, and the maximum strength of the toroidal magnetic field is approximately proportional to the crust thickness. △ Less

Submitted 16 December, 2022; originally announced December 2022.

Comments: 10 pages, 5 figures, accepted for publication in MNRAS

Overview of the Observing System and Initial Scientific Accomplishments of the East Asian VLBI Network (EAVN)

Authors: Kazunori Akiyama, Juan-Carlos Algaba, Tao An, Keiichi Asada, Kitiyanee Asanok, Do-Young Byun, Thanapol Chanapote, Wen Chen, Zhong Chen, Xiaopeng Cheng, James O. Chibueze, Ilje Cho, Se-Hyung Cho, Hyun-Soo Chung, Lang Cui, Yuzhu Cui, Akihiro Doi, Jian Dong, Kenta Fujisawa, Wei Gou, Wen Guo, Kazuhiro Hada, Yoshiaki Hagiwara, Tomoya Hirota, Jeffrey A. Hodgson , et al. (79 additional authors not shown)

Abstract: The East Asian VLBI Network (EAVN) is an international VLBI facility in East Asia and is operated under mutual collaboration between East Asian countries, as well as part of Southeast Asian and European countries. EAVN currently consists of 16 radio telescopes and three correlators located in China, Japan, and Korea, and is operated mainly at three frequency bands, 6.7, 22, and 43 GHz with the lon… ▽ More The East Asian VLBI Network (EAVN) is an international VLBI facility in East Asia and is operated under mutual collaboration between East Asian countries, as well as part of Southeast Asian and European countries. EAVN currently consists of 16 radio telescopes and three correlators located in China, Japan, and Korea, and is operated mainly at three frequency bands, 6.7, 22, and 43 GHz with the longest baseline length of 5078 km, resulting in the highest angular resolution of 0.28 milliarcseconds at 43 GHz. One of distinct capabilities of EAVN is multi-frequency simultaneous data reception at nine telescopes, which enable us to employ the frequency phase transfer technique to obtain better sensitivity at higher observing frequencies. EAVN started its open-use program in the second half of 2018, providing a total observing time of more than 1100 hours in a year. EAVN fills geographical gap in global VLBI array, resulting in enabling us to conduct contiguous high-resolution VLBI observations. EAVN has produced various scientific accomplishments especially in observations toward active galactic nuclei, evolved stars, and star-forming regions. These activities motivate us to initiate launch of the 'Global VLBI Alliance' to provide an opportunity of VLBI observation with the longest baselines on the earth. △ Less

Submitted 14 December, 2022; originally announced December 2022.

Comments: 27 pages, appeared in Galaxies special issue 'Challenges in Understanding Black Hole Powered Jets with VLBI' as an invited review

Journal ref: Galaxies 2022, 10(6), 113

Self-Similar solution of rotating eruptive outflows on its equatorial plane

Authors: Takatoshi Ko, Kotaro Fujisawa, Toshikazu Shigeyama

Abstract: We construct axisymmetric self-similar solutions of transonic outflows emanating from a point source including the effect of the rotation. The solutions are constructed exclusively on the equatorial plane. The features of solutions are determined by three parameters; the adiabatic index $γ$, the dimensionless coordinate of the transonic point, and the dimensionless azimuthal velocity at the transo… ▽ More We construct axisymmetric self-similar solutions of transonic outflows emanating from a point source including the effect of the rotation. The solutions are constructed exclusively on the equatorial plane. The features of solutions are determined by three parameters; the adiabatic index $γ$, the dimensionless coordinate of the transonic point, and the dimensionless azimuthal velocity at the transonic point. We classify the solutions into five groups according to the asymptotic behaviors. We find that the behaviors of the self-similar solutions change at $γ= 11/9$. In addition, some solutions show double-power-law density profiles, which are usually seen in ejecta from a binary merger or nova-like explosion. Thus, our self-similar solutions can be applied not only to the outflow blowing from the central spinning objects, but also to the ejecta erupted from the binary merger or nova-like explosion. △ Less

Submitted 5 November, 2022; originally announced November 2022.

Comments: 14 pages, 6 figures, 1 table. Accepted by ApJ

Report number: RESCEU-21/22 MSC Class: 85-10

A Lagrangian construction of rotating stars

Authors: Misa Ogata, Hirotada Okawa, Kotaro Fujisawa, Nobutoshi Yasutake, Yu Yamamoto, Shoichi Yamada

Abstract: We present a new formulation for numerically obtaining axisymmetric equilibrium structures of rotating stars in two spatial dimensions. With a view to apply it to the secular evolution of rotating stars, we base it on the Lagrangian description, i.e., we solve the force-balance equations to find the spatial positions of fluid elements endowed individually with a mass, specific entropy and angular… ▽ More We present a new formulation for numerically obtaining axisymmetric equilibrium structures of rotating stars in two spatial dimensions. With a view to apply it to the secular evolution of rotating stars, we base it on the Lagrangian description, i.e., we solve the force-balance equations to find the spatial positions of fluid elements endowed individually with a mass, specific entropy and angular momentum. The system of nonlinear equations obtained by finite-differencing the basic equations are solved with the W4 method, which is a new multi-dimensional root-finding scheme of our own devising. We augment it with a remapping scheme to avoid distortions of the Lagrangian coordinates. In this first one of a series of papers, we will give a detailed description of these methods initially. We then present the results of some test calculations, which include the construction of both rapidly rotating barotropic and baroclinic equilibrium states. We gauge their accuracies quantitatively with some diagnostic quantities as well as via comparisons with the counterparts obtained with an Eulerian code. For a demonstrative purpose, we apply the code to a toy-model cooling calculation of a rotating white dwarf. △ Less

Submitted 17 October, 2022; originally announced October 2022.

Comments: 15pages, 14 figures, submitted to MNRAS

Modeling photometric variations due to a global inhomogeneity on an obliquely rotating star: application to lightcurves of white dwarfs

Authors: Yasushi Suto, Shin Sasaki, Masataka Aizawa, Kotaro Fujisawa, Kazumi Kashiyama

Abstract: We develop a general framework to compute photometric variations induced by the oblique rotation of a star with an axisymmetric inhomogeneous surface. We apply the framework to compute lightcurves of white dwarfs adopting two simple models of their surface inhomogeneity. Depending on the surface model and the location of the observer, the resulting lightcurve exhibits a departure from a purely sin… ▽ More We develop a general framework to compute photometric variations induced by the oblique rotation of a star with an axisymmetric inhomogeneous surface. We apply the framework to compute lightcurves of white dwarfs adopting two simple models of their surface inhomogeneity. Depending on the surface model and the location of the observer, the resulting lightcurve exhibits a departure from a purely sinusoidal curve that are observed for a fraction of white dwarfs. As a specific example, we fit our model to the observed phase-folded lightcurve of a fast-spinning white dwarf ZTF J190132.9+145808.7 (with the rotation period of 419s). We find that the size and obliquity angle of the spot responsible for the photometric variation are $\dts \approx 60^\circ$ and $\thetaS \approx 60^\circ$ or $90^\circ$, respectively, implying an interesting constraint on the surface distribution of the magnetic field on white dwarfs. △ Less

Submitted 7 November, 2022; v1 submitted 11 October, 2022; originally announced October 2022.

Comments: 25 pages, 9 figures, accepted for publication in PASJ

Magnetically confined mountains on accreting neutron stars with multipole magnetic fields

Authors: Kotaro Fujisawa, Shota Kisaka, Yasufumi Kojima

Abstract: Magnetically confined mountains on accreting neutron stars are candidates for producing continuous gravitational waves. We formulate a magnetically confined mountain on a neutron star with strong multipole magnetic fields and obtain some sequences of numerical solutions. We find that the mass ellipticity of the mountain increases by one order of magnitude if the neutron star has strong multipole m… ▽ More Magnetically confined mountains on accreting neutron stars are candidates for producing continuous gravitational waves. We formulate a magnetically confined mountain on a neutron star with strong multipole magnetic fields and obtain some sequences of numerical solutions. We find that the mass ellipticity of the mountain increases by one order of magnitude if the neutron star has strong multipole magnetic fields. As matter accretes on to the magnetic pole, the size of the mountain increases and the magnetic fields are buried. If the neutron star has a dipole magnetic field, the dipole magnetic field is buried and transformed into multipole components. By contrast, if the neutron star has both dipole and strong multipole magnetic fields, the multipole magnetic fields are buried and transformed into a negative dipole component. We also calculate magnetically confined mountains with toroidal magnetic fields and find that the ellipticity becomes slightly smaller when the mountain has toroidal magnetic fields. If the multipole magnetic fields are buried, they sustain the intense toroidal magnetic field near the stellar surface, and the ratio of the toroidal magnetic field to the poloidal magnetic field is close to 100. The hidden strong toroidal magnetic fields are sustained by the buried multipole magnetic fields. △ Less

Submitted 9 September, 2022; originally announced September 2022.

Comments: 13 pages, 12 figures, accepted for publication in MNRAS

Proton radiation damage tolerance of wide dynamic range SOI pixel detectors

Authors: Shun Tsunomachi, Takayoshi Kohmura, Kouichi Hagino, Masatoshi Kitajima, Toshiki Doi, Daiki Aoki, Asuka Ohira, Yasuyuki Shimizu, Kaito Fujisawa, Shizusa Yamazaki, Yuusuke Uchida, Makoto Shimizu, Naoki Itoh, Yasuo Arai, Toshinobu Miyoshi, Ryutaro Nishimura, Takeshi Go Tsuru, Ikuo Kurachi

Abstract: We have been developing the SOI pixel detector ``INTPIX'' for space use and general purpose applications such as the residual stress measurement of a rail and high energy physics experiments. INTPIX is a monolithic pixel detector composed of a high-resistivity Si sensor, a SiO2 insulator, and CMOS pixel circuits utilizing Silicon-On-Insulator (SOI) technology. We have considered the possibility of… ▽ More We have been developing the SOI pixel detector ``INTPIX'' for space use and general purpose applications such as the residual stress measurement of a rail and high energy physics experiments. INTPIX is a monolithic pixel detector composed of a high-resistivity Si sensor, a SiO2 insulator, and CMOS pixel circuits utilizing Silicon-On-Insulator (SOI) technology. We have considered the possibility of using INTPIX to observe X-ray polarization in space. When the semiconductor detector is used in space, it is subject to radiation damage resulting from high-energy protons. Therefore, it is necessary to investigate whether INTPIX has high radiation tolerance for use in space. The INTPIX8 was irradiated with 6 MeV protons up to a total dose of 2 krad at HIMAC, National Institute of Quantum Science in Japan, and evaluated the degradation of the performance, such as energy resolution and non-uniformity of gain and readout noise between pixels. After 500 rad irradiation, which is the typical lifetime of an X-ray astronomy satellite, the degradation of energy resolution at 14.4 keV is less than 10%, and the non-uniformity of readout noise and gain between pixels is constant within 0.1%. △ Less

Submitted 8 September, 2022; originally announced September 2022.

Comments: 7 pages, 8 figures, published in proceedings for SPIE Astronomical Telescopes + Instrumentation in 2022

Fast optical flares from M dwarfs detected by a one-second-cadence survey with Tomo-e Gozen

Authors: Masataka Aizawa, Kojiro Kawana, Kazumi Kashiyama, Ryou Ohsawa, Hajime Kawahara, Fumihiro Naokawa, Tomoyuki Tajiri, Noriaki Arima, Hanchun Jiang, Tilman Hartwig, Kotaro Fujisawa, Toshikazu Shigeyama, Ko Arimatsu, Mamoru Doi, Toshihiro Kasuga, Naoto Kobayashi, Sohei Kondo, Yuki Mori, Shin-ichiro Okumura, Satoshi Takita, Shigeyuki Sako

Abstract: We report a one-second-cadence wide-field survey for M-dwarf flares using the Tomo-e Gozen camera mounted on the Kiso Schmidt telescope. We detect 22 flares from M3-M5 dwarfs with rise times and amplitudes ranging from $5\, \mathrm{sec} \lesssim t_\mathrm{rise} \lesssim 100\,\mathrm{sec}$ and $0.5 \lesssim ΔF/F_{\star} \lesssim 20$, respectively. The flare light curves mostly show steeper rises an… ▽ More We report a one-second-cadence wide-field survey for M-dwarf flares using the Tomo-e Gozen camera mounted on the Kiso Schmidt telescope. We detect 22 flares from M3-M5 dwarfs with rise times and amplitudes ranging from $5\, \mathrm{sec} \lesssim t_\mathrm{rise} \lesssim 100\,\mathrm{sec}$ and $0.5 \lesssim ΔF/F_{\star} \lesssim 20$, respectively. The flare light curves mostly show steeper rises and shallower decays than those obtained from the Kepler one-minute cadence data and tend to have flat peak structures. Assuming a blackbody spectrum with temperatures of $9,000-15,000\,\mathrm{K}$, the peak luminosities and bolometric energies are estimated to be $10^{29}\,\mathrm{erg\,sec^{-1}} \lesssim L_\mathrm{peak} \lesssim 10^{31}\,\mathrm{erg\,sec^{-1}}$ and $10^{31}\,\mathrm{erg} \lesssim E_{\rm bol} \lesssim 10^{34}\,\mathrm{erg}$, which constitutes the bright end of fast optical flares for M dwarfs. We confirm that more than 90\% of the host stars of the detected flares are magnetically active based on their H$α$ emission line intensities obtained by LAMOST. The estimated occurrence rate of the detected flares is $\sim 0.7$ per day per an active star, indicating they are common in magnetically active M dwarfs. We argue that the flare light curves can be explained by the chromospheric compression model; the rise time is broadly consistent with the Alfvén transit time of a magnetic loop with a length scale of $l_\mathrm{loop} \sim 10^4\,\mathrm{km}$ and a field strength of $1,000\,\mathrm{G}$, while the decay time is likely determined by the radiative cooling of the compressed chromosphere down to near the photosphere with a temperature of $\gtrsim 10,000\,\mathrm{K}$. These flares from M dwarfs could be a major contamination source for a future search of fast optical transients of unknown types. △ Less

Submitted 26 June, 2022; originally announced June 2022.

Comments: 39 pages, 13 figures, accepted for publication in PASJ

A novel formulation for the evolution of relativistic rotating stars

Authors: Hirotada Okawa, Kotaro Fujisawa, Nobutoshi Yasutake, Misa Ogata, Yu Yamamoto, Shoichi Yamada

Abstract: We present a new formulation to construct numerically equilibrium configurations of rotating stars in general relativity. Having in mind the application to their quasi static evolutions, we adopt a Lagrangian formulation of our own devising, in which we solve force balance equations to seek for the positions of fluid elements assigned to the grid points, instead of the ordinary Eulerian formulatio… ▽ More We present a new formulation to construct numerically equilibrium configurations of rotating stars in general relativity. Having in mind the application to their quasi static evolutions, we adopt a Lagrangian formulation of our own devising, in which we solve force balance equations to seek for the positions of fluid elements assigned to the grid points, instead of the ordinary Eulerian formulation. Unlike previous works in the literature, we do not employ the first integral of the Euler equation, which is not obtained by an analytic integration in general. We assign a mass, specific angular momentum and entropy to each fluid element in contrast to the previous methods, in which the spatial distribution of the angular velocity or angular momentum is specified. Those distributions are determined after the positions of all fluid elements (or grid points) are derived in our formulation. We solve the large system of algebraic nonlinear equations that are obtained by discretizing the time-independent Euler and Einstein equations in the finite-elements method by using our new multi-dimensional root-finding scheme, named the W4 method. To demonstrate the capability of our new formulation, we construct some rotational configurations both barotropic and baroclinic. We also solve three evolutionary sequences that mimic the cooling, mass-loss, and mass-accretion as simple toy models. △ Less

Submitted 21 April, 2022; originally announced April 2022.

Comments: 19 pages, 13 figures

Singularity-Avoiding Multi-Dimensional Root-Finder

Authors: Hirotada Okawa, Kotaro Fujisawa, Yu Yamamoto, Nobutoshi Yasutake, Misa Ogata, Shoichi Yamada

Abstract: We proposed in this paper a new method, which we named the W4 method, to solve nonlinear equation systems. It may be regarded as an extension of the Newton-Raphson~(NR) method to be used when the method fails. Indeed our method can be applied not only to ordinary problems with non-singular Jacobian matrices but also to problems with singular Jacobians, which essentially all previous methods that e… ▽ More We proposed in this paper a new method, which we named the W4 method, to solve nonlinear equation systems. It may be regarded as an extension of the Newton-Raphson~(NR) method to be used when the method fails. Indeed our method can be applied not only to ordinary problems with non-singular Jacobian matrices but also to problems with singular Jacobians, which essentially all previous methods that employ the inversion of the Jacobian matrix have failed to solve. In this article, we demonstrate that (i) our new scheme can define a non-singular iteration map even for those problems by utilizing the singular value decomposition, (ii) a series of vectors in the new iteration map converges to the right solution under a certain condition, (iii) the standard two-dimensional problems in the literature that no single method proposed so far has been able to solve completely are all solved by our new method. △ Less

Submitted 21 April, 2022; originally announced April 2022.

Comments: 13 pages, 4 figures

Magnetic field sustained by the elastic force in neutron star crusts

Authors: Yasufumi Kojima, Shota Kisaka, Kotaro Fujisawa

Abstract: We investigate the magneto--elastic equilibrium of a neutron star crust and magnetic energy stored by the elastic force. The solenoidal motion driven by the Lorentz force can be controlled by the magnetic elastic force, so that conditions for the magnetic field strength and geometry are less restrictive. For equilibrium models, the minor solenoidal part of the magnetic force is balanced by a weak… ▽ More We investigate the magneto--elastic equilibrium of a neutron star crust and magnetic energy stored by the elastic force. The solenoidal motion driven by the Lorentz force can be controlled by the magnetic elastic force, so that conditions for the magnetic field strength and geometry are less restrictive. For equilibrium models, the minor solenoidal part of the magnetic force is balanced by a weak elastic force because the irrotational part is balanced by the dominant gravity and pressure forces. Therefore, a strong magnetic field may be confined in the interior, regardless of poloidal or toroidal components. We numerically calculated axially symmetric models with the maximum shear--strain, and found that a magnetic energy $> 10^{46}$ erg can be stored in the crust, even for a normal surface dipole-field-strength ($<10^{13}$ G). The magnetic energy much exceeds the elastic energy ($ 10^{44} -10^{45}$ erg). The shear--stress spatial distribution revealed that the elastic structure is likely to break down near the surface. In particular, the critical position is highly localized at a depth less than 100 m from the surface. △ Less

Submitted 5 January, 2022; originally announced January 2022.

Comments: 11 pages, 4 figures

Magneto-elastic equilibrium of a neutron-star crust

Authors: Yasufumi Kojima, Shota Kisaka, Kotaro Fujisawa

Abstract: We examine the equilibrium of a magnetized neutron-star-crust. We calculate axially symmetric models in which an elastic force balances solenoidal motion driven by a Lorentz force. A large variety of equilibrium models are allowed by incorporating the elastic shear deformation; in addition, toroidal-magnetic-field dominated models are available. These results remarkably differ from those in barotr… ▽ More We examine the equilibrium of a magnetized neutron-star-crust. We calculate axially symmetric models in which an elastic force balances solenoidal motion driven by a Lorentz force. A large variety of equilibrium models are allowed by incorporating the elastic shear deformation; in addition, toroidal-magnetic-field dominated models are available. These results remarkably differ from those in barotropic fluid stars. We demonstrate some models wherein the magnetic energy exceeds the elastic energy. The excess comes from the fact that a large amount of magnetic energy is associated with the irrotational part of the magnetic force, which is balanced with gravity and pressure. It is sufficient for equilibrium models that the minor solenoidal part is balanced by a weak elastic force. We find that the elasticity in the crust plays an important role on the magnetic-field confinement. Further, we present the spatial distribution of the shear-stress at the elastic limit, by which the crust-fracture location can be identified. The result has useful implications for realistic crust-quake models. △ Less

Submitted 28 June, 2021; v1 submitted 27 June, 2021; originally announced June 2021.

Comments: 13 pages, 9 figures

Evolution of magnetic deformation in neutron star crust

Authors: Yasufumi Kojima, Shota Kisaka, Kotaro Fujisawa

Abstract: In this study, we examine the magnetic field evolution occurring in a neutron star crust. Beyond the elastic limit, the lattice ions are assumed to act as a plastic flow. The Ohmic dissipation, Hall drift, and bulk fluid velocity driven by the Lorentz force are considered in our numerical simulation. A magnetically induced quadrupole deformation is observed in the crust during the evolution. Gener… ▽ More In this study, we examine the magnetic field evolution occurring in a neutron star crust. Beyond the elastic limit, the lattice ions are assumed to act as a plastic flow. The Ohmic dissipation, Hall drift, and bulk fluid velocity driven by the Lorentz force are considered in our numerical simulation. A magnetically induced quadrupole deformation is observed in the crust during the evolution. Generally, the ellipticity decreases as the magnetic energy decreases. In a toroidal-field-dominated model, the sign of the ellipticity changes. Namely, the initial prolate shape tends to become oblate. This occurs because the toroidal component decays rapidly on a smaller timescale than the poloidal dipole component. We find that the magnetic dipole component does not change significantly on the Hall timescale of $\sim 1$Myr for the considered simple initial models. Thus, a more complex initial model is required to study the fast decay of surface dipoles on the abovementioned timescale. △ Less

Submitted 6 November, 2020; originally announced November 2020.

Comments: 12 pages, 7 figures

Near-infrared Monitoring of the Accretion Outburst in the MYSO S255-NIRS3

Authors: Mizuho Uchiyama, Takuya Yamashita, Koichiro Sugiyama, Tatsuya Nakaoka, Miho Kawabata, Ryosuke Itoh, Masayuki Yamanaka, Hiroshi Akitaya, Koji Kawabata, Yoshinori Yonekura, Yu Saito, Kazuhito Motogi, Kenta Fujisawa

Abstract: We followed-up the massive young stellar object (MYSO) S255-NIRS3 (=S255-IRS1b) during its recent accretion outburst event in the Ks band with Kanata/HONIR for four years after its burst and obtained a long-term light curve. This is the most complete NIR light-curve of the S255-NIRS3 burst event that has ever been presented. The light curve showed a steep increase reaching a peak flux that was 3.4… ▽ More We followed-up the massive young stellar object (MYSO) S255-NIRS3 (=S255-IRS1b) during its recent accretion outburst event in the Ks band with Kanata/HONIR for four years after its burst and obtained a long-term light curve. This is the most complete NIR light-curve of the S255-NIRS3 burst event that has ever been presented. The light curve showed a steep increase reaching a peak flux that was 3.4 mag brighter than the quiescent phase and then a relatively moderate year-scale fading until the last observation, similar to that of the accretion burst events such as EXors found in lower-mass young stellar objects. The behavior of the Ks band light curve is similar to that observed in 6.7 GHz class II methanol maser emission, with a sudden increase followed by moderate year-scale fading. However, the maser emission peaks appear 30-50 days earlier than that of the Ks band emission. The similarities confirmed that the origins of the maser emission and the Ks band continuum emission is common as previously shown from another infrared and radio observations by Stecklum et al. (2016); Caratti o Garatti et al. (2017a); Moscadelli et al. (2017). However, the differences in energy transfer paths, such as the exciting/emitting/scattering structures, may cause the delay in the flux-peak dates. △ Less

Submitted 16 October, 2019; originally announced October 2019.

Comments: 8 pages, 4 figures, 2 tables, accepted to PASJ

The optically thick rotating magnetic wind from a massive white dwarf merger product

Authors: Kazumi Kashiyama, Kotaro Fujisawa, Toshikazu Shigeyama

Abstract: WD J005311 is a newly identified white dwarf (WD) in a mid-infrared nebula. The spectroscopic observation indicates the existence of a neon-enriched carbon/oxygen wind with a terminal velocity of $v_{\infty,\rm obs}\sim 16,000\,\rm km\,s^{-1}$ and a mass loss rate of $\dot M_{\rm obs}\sim 3.5\times 10^{-6}\,M_\odot$ yr$^{-1}$. Here we consistently explain the properties of WD J005311 using a newly… ▽ More WD J005311 is a newly identified white dwarf (WD) in a mid-infrared nebula. The spectroscopic observation indicates the existence of a neon-enriched carbon/oxygen wind with a terminal velocity of $v_{\infty,\rm obs}\sim 16,000\,\rm km\,s^{-1}$ and a mass loss rate of $\dot M_{\rm obs}\sim 3.5\times 10^{-6}\,M_\odot$ yr$^{-1}$. Here we consistently explain the properties of WD J005311 using a newly constructed wind solution, where the optically thick outflow is launched from the carbon burning shell on an oxygen-neon core and accelerated by the rotating magnetic field to become supersonic and unbound well below the photosphere. Our model implies that WD J005311 has a mass of $M_* \sim 1.1\mbox{-}1.3\,M_\odot$, a magnetic field of $B_* \sim (2\mbox{-}5)\times 10^7\,\rm G$, and a spin angular frequency of $Ω\sim 0.2\mbox{-}0.5 \,\rm s^{-1}$. The large magnetic field and fast spin support the carbon-oxygen WD merger origin. WD J005311 will neither explode as a type Ia supernova nor collapse into a neutron star. If the wind continues to blow another few kyr, WD J005311 will spin down significantly and join to the known sequence of slowly-rotating magnetic WDs. Otherwise it may appear as a fast-spinning magnetic WD and could be a new high energy source. △ Less

Submitted 16 October, 2019; v1 submitted 29 July, 2019; originally announced July 2019.

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

Stellar Overdensity in the Local Arm in Gaia DR2

Authors: Yusuke Miyachi, Nobuyuki Sakai, Daisuke Kawata, Junichi Baba, Mareki Honma, Noriyuki Matsunaga, Kenta Fujisawa

Abstract: Using the cross-matched data of Gaia DR2 and 2MASS Point Source Catalog, we investigated the surface density distribution of stars aged ~1 Gyr in the thin disk in the range of 90° <= l <= 270°. We selected 4,654 stars above the turnoff corresponding to the age ~1 Gyr, that fall within a small box region in the color-magnitude diagram, (J-Ks)0 versus M(Ks), for which the distance and reddening are… ▽ More Using the cross-matched data of Gaia DR2 and 2MASS Point Source Catalog, we investigated the surface density distribution of stars aged ~1 Gyr in the thin disk in the range of 90° <= l <= 270°. We selected 4,654 stars above the turnoff corresponding to the age ~1 Gyr, that fall within a small box region in the color-magnitude diagram, (J-Ks)0 versus M(Ks), for which the distance and reddening are corrected. The selected sample shows an arm-like overdensity at 90° <= l <= 190°.This overdensity is located close to the Local arm traced by high-mass star forming regions (HMSFRs), but its pitch angle is slightly larger than that of the HMSFR-defined arm. Although the significance of the overdensity we report is marginal, its structure poses questions concerning both of the competing scenarios of spiral arms, the density-wave theory and the dynamic spiral arm model. The offset between the arms traced by stars and HMSFRs, i.e., gas, is difficult to be explained by the dynamic arm scenario. On the other hand, the pitch angle of the stellar Local arm, if confirmed, larger than that of the Perseus arm is difficult to be explained with the classical density-wave scenario. The dynamic arm scenario can explain it if the Local arm is in a growing up phase, while the Perseus arm is in a disrupting phase. Our result provides a new and complex picture of the Galactic spiral arms, and encourages further studies. △ Less

Submitted 8 July, 2019; originally announced July 2019.

Comments: 11 pages, 4 figures, accepted for publication in ApJ

New code for equilibriums and quasiequilibrium initial data of compact objects. IV. Rotating relativistic stars with mixed poloidal and toroidal magnetic fields

Authors: Koji Uryu, Shijun Yoshida, Eric Gourgoulhon, Charalampos Markakis, Kotaro Fujisawa, Antonios Tsokaros, Keisuke Taniguchi, Yoshiharu Eriguchi

Abstract: A new code for computing fully general relativistic solutions of strongly magnetized rapidly rotating compact stars is developed as a part of the COCAL (Compact Object CALculator) code. The full set of Einstein's equations, Maxwell's equations and magnetohydrodynamic equations are consistently solved assuming perfect conductivity, stationarity, and axisymmetry, and strongly magnetized solutions as… ▽ More A new code for computing fully general relativistic solutions of strongly magnetized rapidly rotating compact stars is developed as a part of the COCAL (Compact Object CALculator) code. The full set of Einstein's equations, Maxwell's equations and magnetohydrodynamic equations are consistently solved assuming perfect conductivity, stationarity, and axisymmetry, and strongly magnetized solutions associated with mixed poloidal and toroidal components of magnetic fields are successfully obtained in generic (non-circular) spacetimes. We introduce the formulation of the problem and the numerical method in detail, then present examples of extremely magnetized compact star solutions and their convergence tests. It is found that, in extremely magnetized stars, the stellar matter can be expelled from the region of strongest toroidal fields. Hence we conjecture that a toroidal electro-vacuum region may appear inside of the extremely magnetized compact stars, which may seem like the neutron star becoming the strongest toroidal solenoid coil in the universe. △ Less

Submitted 25 June, 2019; originally announced June 2019.

Comments: 30 pages, 4 figures

Journal ref: Phys. Rev. D 100, 123019 (2019)

The Galactic Center Lobe Filled with Thermal Plasma

Authors: Halca Nagoshi, Yuzo Kubose, Kenta Fujisawa, Kazuo Sorai, Yoshinori Yonekura, Koichiro Sugiyama, Kotaro Niinuma, Kazuhito Motogi, Takahiro Aoki

Abstract: An observational result of a radio continuum and H92$α$ radio recombination line of the Galactic Center Lobe (GCL), using the Yamaguchi 32 m radio telescope, is reported. The obtained spatial intensity distribution of the radio recombination line shows two distinctive ridge-like structures extending from the galactic plane vertically to the north at the eastern and western sides of the galactic ce… ▽ More An observational result of a radio continuum and H92$α$ radio recombination line of the Galactic Center Lobe (GCL), using the Yamaguchi 32 m radio telescope, is reported. The obtained spatial intensity distribution of the radio recombination line shows two distinctive ridge-like structures extending from the galactic plane vertically to the north at the eastern and western sides of the galactic center, which are connected to each other at a latitude of $1.2^{\circ}$ to form a loop-like structure as a whole. This suggests that most of the radio continuum emission of the GCL is free-free emission, and that the GCL is filled with thermal plasma. The east ridge of the GCL observed with the radio recombination line separates 30 pc from the radio arc, which has been considered as a part of the GCL, but coincides with a ridge of the radio continuum at a galactic longitude of $0^{\circ}$. The radial velocity of the radio recombination line is found to be between $-4$ and $+10$ km s$^{-1}$ across the GCL. This velocity is much smaller than the one expected from the galactic rotation, and hence indicates that the GCL exists apart from the galactic center. These characteristics of the GCL suggest that the long-standing hypothesis that the GCL was created by an explosive activity in the galactic center is unlikely, but favor that the GCL is a giant HII region. △ Less

Submitted 5 May, 2019; originally announced May 2019.

A Balloon-Borne Very Long Baseline Interferometry Experiment in the Stratosphere: Systems Design and Developments

Authors: Akihiro Doi, Yusuke Kono, Kimihiro Kimura, Satomi Nakahara, Tomoaki Oyama, Nozomi Okada, Yasutaka Satou, Kazuyoshi Yamashita, Naoko Matsumoto, Mitsuhisa Baba, Daisuke Yasuda, Shunsaku Suzuki, Yutaka Hasegawa, Mareki Honma, Hiroaki Tanaka, Kosei Ishimura, Yasuhiro Murata, Reiho Shimomukai, Tomohiro Tachi, Kazuya Saito, Naohiko Watanabe, Nobutaka Bando, Osamu Kameya, Yoshinori Yonekura, Mamoru Sekido , et al. (11 additional authors not shown)

Abstract: The balloon-borne very long baseline interferometry (VLBI) experiment is a technical feasibility study for performing radio interferometry in the stratosphere. The flight model has been developed. A balloon-borne VLBI station will be launched to establish interferometric fringes with ground-based VLBI stations distributed over the Japanese islands at an observing frequency of approximately 20 GHz… ▽ More The balloon-borne very long baseline interferometry (VLBI) experiment is a technical feasibility study for performing radio interferometry in the stratosphere. The flight model has been developed. A balloon-borne VLBI station will be launched to establish interferometric fringes with ground-based VLBI stations distributed over the Japanese islands at an observing frequency of approximately 20 GHz as the first step. This paper describes the system design and development of a series of observing instruments and bus systems. In addition to the advantages of avoiding the atmospheric effects of absorption and fluctuation in high frequency radio observation, the mobility of a station can improve the sampling coverage (`uv-coverage') by increasing the number of baselines by the number of ground-based counterparts for each observation day. This benefit cannot be obtained with conventional arrays that solely comprise ground-based stations. The balloon-borne VLBI can contribute to a future progress of research fields such as black holes by direct imaging. △ Less

Submitted 11 December, 2018; originally announced December 2018.

Comments: 15 pages, 6 figures, 1 table. Advances in Space Research, in press

The W4 method: a new multi-dimensional root-finding scheme for nonlinear systems of equations

Authors: Hirotada Okawa, Kotaro Fujisawa, Yu Yamamoto, Ryosuke Hirai, Nobutoshi Yasutake, Hiroki Nagakura, Shoichi Yamada

Abstract: We propose a new class of method for solving nonlinear systems of equations, which, among other things,has four nice features: (i) it is inspired by the mathematical property of damped oscillators, (ii) it can be regarded as a simple extention to the Newton-Raphson(NR) method, (iii) it has the same local convergence as the NR method does, (iv) it has a significantly wider convergence region or the… ▽ More We propose a new class of method for solving nonlinear systems of equations, which, among other things,has four nice features: (i) it is inspired by the mathematical property of damped oscillators, (ii) it can be regarded as a simple extention to the Newton-Raphson(NR) method, (iii) it has the same local convergence as the NR method does, (iv) it has a significantly wider convergence region or the global convergence than that of the NR method. In this article, we present the evidence of these properties, applying our new method to some examples and comparing it with the NR method. △ Less

Submitted 12 September, 2018; originally announced September 2018.

Comments: 17 pages, 7 figures

Effects of rotation and magnetic field on the revival of a stalled shock in supernova explosions

Authors: Kotaro Fujisawa, Hirotada Okawa, Yu Yamamoto, Shoichi Yamada

Abstract: We investigate axisymmetric steady solutions of (magneto)hydrodynamics equations that describe approximately accretion flows through a standing shock wave and discuss the effects of rotation and magnetic field on the revival of the stalled shock wave in supernova explosions. We develop a new powerful numerical method to calculate the 2-dimensional (2D) steady accretion flows self-consistently. We… ▽ More We investigate axisymmetric steady solutions of (magneto)hydrodynamics equations that describe approximately accretion flows through a standing shock wave and discuss the effects of rotation and magnetic field on the revival of the stalled shock wave in supernova explosions. We develop a new powerful numerical method to calculate the 2-dimensional (2D) steady accretion flows self-consistently. We first confirm the results of preceding papers that there is a critical luminosity of irradiating neutrinos, above which there exists no steady solution in spherical models. If a collapsing star has rotation and/or magnetic field, the accretion flows are no longer spherical owing to the centrifugal force and/or Lorentz force and the critical luminosity is modified.In fact we find that the critical luminosity is reduced by about 50% - 70% for rapid rotations and about 20% - 50% for strong toroidal magnetic fields, depending on the mass accretion rate. These results may be also interpreted as an existence of the critical specific angular momentum or critical magnetic field, above which there exists no steady solution and the standing shock wave will revive for a given combination of mass accretion rate and neutrino luminosity. △ Less

Submitted 12 September, 2018; originally announced September 2018.

Comments: 22 pages, 9 figures, to be submitted

A comparative study of amplitude calibrations for East-Asia VLBI Network: a-priori and template spectrum methods

Authors: Ilje Cho, Taehyun Jung, Guang-Yao Zhao, Kazunori Akiyama, Satoko Sawada-Satoh, Motoki Kino, Do-Young Byun, Bongwon Sohn, Katsunori M. Shibata, Tomoya Hirota, Kotaro Niinuma, Yoshinori Yonekura, Kenta Fujisawa, Tomoaki Oyama

Abstract: We present the results of comparative study of amplitude calibrations for East-Asia VLBI Network (EAVN) at 22 and 43 GHz using two different methods of an "a-priori" and a "template spectrum", particularly on lower declination sources. Using observational data sets of early EAVN observations, we investigated the elevation-dependence of the gain values at seven stations of the KaVA (KVN and VERA Ar… ▽ More We present the results of comparative study of amplitude calibrations for East-Asia VLBI Network (EAVN) at 22 and 43 GHz using two different methods of an "a-priori" and a "template spectrum", particularly on lower declination sources. Using observational data sets of early EAVN observations, we investigated the elevation-dependence of the gain values at seven stations of the KaVA (KVN and VERA Array) and three additional telescopes in Japan (Takahagi 32m, Yamaguchi 32m and Nobeyama 45m). By comparing the independently obtained gain values based on these two methods, we found that the gain values from each method were consistent within 10% at elevations higher than 10 degree. We also found that the total flux densities of two images produced from the different amplitude calibrations were in agreement within 10% at both 22 and 43 GHz. By using the template spectrum method, furthermore, the additional radio telescopes can participate in the KaVA (i.e. EAVN) so that it can give a notable sensitivity increase. Therefore, our results will constrain the detailed conditions to reliably measure the VLBI amplitude using EAVN and give a potential to extend possible telescopes comprising EAVN. △ Less

Submitted 27 October, 2017; originally announced October 2017.

Comments: in press, 26 pages, 8 figures, 4 tables

Journal ref: Publications of the Astronomical Society of Japan, 2017

Subaru Hyper Suprime-Cam Survey for An Optical Counterpart of GW170817

Authors: Nozomu Tominaga, Masaomi Tanaka, Tomoki Morokuma, Yousuke Utsumi, Masaki S. Yamaguchi, Naoki Yasuda, Masayuki Tanaka, Michitoshi Yoshida, Takuya Fujiyoshi, Hisanori Furusawa, Koji S. Kawabata, Chien-Hsiu Lee, Kentaro Motohara, Ryou Ohsawa, Kouji Ohta, Tsuyoshi Terai, Fumio Abe, Wako Aoki, Yuichiro Asakura, Sudhanshu Barway, Ian A. Bond, Kenta Fujisawa, Satoshi Honda, Kunihito Ioka, Youichi Itoh , et al. (9 additional authors not shown)

Abstract: We perform a $z$-band survey for an optical counterpart of a binary neutron star coalescence GW170817 with Subaru/Hyper Suprime-Cam. Our untargeted transient search covers $23.6$ deg$^2$ corresponding to the $56.6\%$ credible region of GW170817 and reaches the $50\%$ completeness magnitude of $20.6$ mag on average. As a result, we find 60 candidates of extragalactic transients, including J-GEM17bt… ▽ More We perform a $z$-band survey for an optical counterpart of a binary neutron star coalescence GW170817 with Subaru/Hyper Suprime-Cam. Our untargeted transient search covers $23.6$ deg$^2$ corresponding to the $56.6\%$ credible region of GW170817 and reaches the $50\%$ completeness magnitude of $20.6$ mag on average. As a result, we find 60 candidates of extragalactic transients, including J-GEM17btc (a.k.a. SSS17a/DLT17ck). While J-GEM17btc is associated with NGC 4993 that is firmly located inside the 3D skymap of GW170817, the other 59 candidates do not have distance information in the GLADE v2 catalog or NASA/IPAC Extragalactic Database (NED). Among 59 candidates, 58 are located at the center of extended objects in the Pan-STARRS1 catalog, while one candidate has an offset. We present location, $z$-band apparent magnitude, and time variability of the candidates and evaluate the probabilities that they are located inside of the 3D skymap of GW170817. The probability for J-GEM17btc is $64\%$ being much higher than those for the other 59 candidates ($9.3\times10^{-3}-2.1\times10^{-1}\%$). Furthermore, the possibility, that at least one of the other 59 candidates is located within the 3D skymap, is only $3.2\%$. Therefore, we conclude that J-GEM17btc is the most-likely and distinguished candidate as the optical counterpart of GW170817. △ Less

Submitted 7 January, 2018; v1 submitted 16 October, 2017; originally announced October 2017.

Comments: 14 pages, 9 figures. Accepted for publication in PASJ (Publications of the Astronomical Society of Japan)

A Face-on Accretion System in High-Mass Star-Formation: Possible Dusty Infall Streams within 100 AU

Authors: Kazuhito Motogi, Tomoya Hirota, Kazuo Sorai, Yoshinori Yonekura, Koichiro Sugiyama, Mareki Honma, Kotaro Niinuma, Kazuya Hachisuka, Kenta Fujisawa, Andrew J. Walsh

Abstract: We report on interferometric observations of a face-on accretion system around the High-Mass young stellar object, G353.273+0.641. The innermost accretion system of 100 au radius was resolved in a 45 GHz continuum image taken with the Jansky-Very Large Array. Our spectral energy distribution analysis indicated that the continuum could be explained by optically thick dust emission. The total mass o… ▽ More We report on interferometric observations of a face-on accretion system around the High-Mass young stellar object, G353.273+0.641. The innermost accretion system of 100 au radius was resolved in a 45 GHz continuum image taken with the Jansky-Very Large Array. Our spectral energy distribution analysis indicated that the continuum could be explained by optically thick dust emission. The total mass of the dusty system is $\sim$ 0.2 $M_{\sun}$ at minimum and up to a few $M_{\sun}$ depending on the dust parameters. 6.7 GHz CH$_{3}$OH masers associated with the same system were also observed with the Australia Telescope Compact Array. The masers showed a spiral-like, non-axisymmetric distribution with a systematic velocity gradient. The line-of-sight velocity field is explained by an infall motion along a parabolic streamline that falls onto the equatorial plane of the face-on system. The streamline is quasi-radial and reaches the equatorial plane at a radius of 16 au. This is clearly smaller than that of typical accretion disks in High-Mass star formation, indicating that the initial angular momentum was very small, or the CH$_{3}$OH masers selectively trace accreting material that has small angular momentum. In the former case, the initial specific angular momentum is estimated to be 8 $\times$ 10$^{20}$ ($M_{*}$$/$10 $M_{\sun}$)$^{0.5}$ cm$^{2}$ s$^{-1}$, or a significant fraction of the initial angular momentum was removed outside of 100 au. The physical origin of such a streamline is still an open question and will be constrained by the higher-resolution ($\sim$ 10 mas) thermal continuum and line observations with ALMA long baselines. △ Less

Submitted 15 October, 2017; originally announced October 2017.

Comments: 17 pages, 11 figures

Development of the Phase-up Technology of the Radio Telescopes: 6.7 GHz Methanol Maser Observations with Phased Hitachi 32 m and Takahagi 32 m Radio Telescopes

Authors: Kazuhiro Takefuji, Koichiro Sugiyama, Yoshinori Yonekura, Tagiru Saito, Kenta Fujisawa, Tetsuro Kondo

Abstract: For the sake of high-sensitivity 6.7 GHz methanol maser observations, we developed a new technology for coherently combining the two signals from the Hitachi 32 m radio telescope and the Takahagi 32 m radio telescope of the Japanese Very long baseline interferometer Network (JVN), where the two telescopes were separated by about 260 m. After the two telescopes were phased as a twofold larger singl… ▽ More For the sake of high-sensitivity 6.7 GHz methanol maser observations, we developed a new technology for coherently combining the two signals from the Hitachi 32 m radio telescope and the Takahagi 32 m radio telescope of the Japanese Very long baseline interferometer Network (JVN), where the two telescopes were separated by about 260 m. After the two telescopes were phased as a twofold larger single telescope, the mean signal-to-noise ratio (SNR) of the 6.7 GHz methanol masers observed by the phased telescopes was improved to 1.254-fold higher than that of the single dish, through a Very Long Baseline Interferometry (VLBI) experiment on the 50 km baseline of the Kashima 34 m telescope and the 1000 km baseline of the Yamaguchi 32 m telescope. Furthermore, we compared the SNRs of the 6.7 GHz maser spectra for two methods. One is a VLBI method and the other is the newly developed digital position switching, which is a similar technology to that used in noise-cancelling headphones. Finally, we confirmed that the mean SNR of method of the digital position switching (ON-OFF) was 1.597-fold higher than that of the VLBI method. △ Less

Submitted 31 August, 2017; originally announced August 2017.

Comments: 17 pages, 10 figures, accepted for publication in Publications of the Astronomical Society of the Pacific (PASP)

OISTER Optical and Near-Infrared Monitoring Observations of a Peculiar Radio-Loud Active Galactic Nucleus SDSS J110006.07+442144.3

Authors: Tomoki Morokuma, Masaomi Tanaka, Yasuyuki T. Tanaka, Ryosuke Itoh, Nozomu Tominaga, Poshak Gandhi, Elena Pian, Paolo Mazzali, Kouji Ohta, Emiko Matsumoto, Takumi Shibata, Hinako Akimoto, Hiroshi Akitaya, Gamal B. Ali, Tsutomu Aoki, Mamoru Doi, Nana Ebisuda, Ahmed Essam, Kenta Fujisawa, Hideo Fukushima, Shuhei Goda, Yuya Gouda, Hidekazu Hanayama, Yasuhito Hashiba, Osamu Hashimoto , et al. (58 additional authors not shown)

Abstract: We present monitoring campaign observations at optical and near-infrared (NIR) wavelengths for a radio-loud active galactic nucleus (AGN) at z=0.840, SDSS~J110006.07+442144.3 (hereafter, J1100+4421), which was identified during a flare phase in late February, 2014. The campaigns consist of three intensive observing runs from the discovery to March, 2015, mostly within the scheme of the OISTER coll… ▽ More We present monitoring campaign observations at optical and near-infrared (NIR) wavelengths for a radio-loud active galactic nucleus (AGN) at z=0.840, SDSS~J110006.07+442144.3 (hereafter, J1100+4421), which was identified during a flare phase in late February, 2014. The campaigns consist of three intensive observing runs from the discovery to March, 2015, mostly within the scheme of the OISTER collaboration. Optical-NIR light curves and simultaneous spectral energy distributions (SEDs) are obtained. Our measurements show the strongest brightening in March, 2015. We found that the optical-NIR SEDs of J1100+4421 show an almost steady shape despite the large and rapid intranight variability. This constant SED shape is confirmed to extend to $\sim5~μ$m in the observed frame using the archival WISE data. Given the lack of absorption lines and the steep power-law spectrum of $α_ν\sim-1.4$, where $f_ν\proptoν^{α_ν}$, synchrotron radiation by a relativistic jet with no or small contributions from the host galaxy and the accretion disk seems most plausible as an optical-NIR emission mechanism. The steep optical-NIR spectral shape and the large amplitude of variability are consistent with this object being a low $ν_{\rm{peak}}$ jet-dominated AGN. In addition, sub-arcsec resolution optical imaging data taken with Subaru Hyper Suprime-Cam does not show a clear extended component and the spatial scales are significantly smaller than the large extensions detected at radio wavelengths. The optical spectrum of a possible faint companion galaxy does not show any emission lines at the same redshift and hence a merging hypothesis for this AGN-related activity is not supported by our observations. △ Less

Submitted 17 July, 2017; originally announced July 2017.

Comments: 22 pages, 7 figures, 2 tables. Accepted for publication in PASJ

White Paper on East Asian Vision for mm/submm VLBI: Toward Black Hole Astrophysics down to Angular Resolution of 1~R$_{S}$

Authors: K. Asada, M. Kino, M. Honma, T. Hirota, R. -S. Lu, M. Inoue, B. -W. Sohn, Z. -Q. Shen, P. T. P. Ho, K. Akiyama, J-C. Algaba, T. An, G. Bower, D-Y. Byun, R. Dodson, A. Doi, P. G. Edwards, K. Fujisawa, M-F. Gu, K. Hada, Y. Hagiwara, P. Jaroenjittichai, T. Jung, T. Kawashima, S. Koyama , et al. (13 additional authors not shown)

Abstract: This White Paper details the intentions and plans of the East Asian Very Long Baseline Interferometry (VLBI) community for pushing the frontiers of millimeter/submillimeter VLBI. To this end, we shall endeavor to actively promote coordinated efforts in the East Asia region. Our goal is to establish firm collaborations among the East Asia VLBI community in partnership with related institutes in Nor… ▽ More This White Paper details the intentions and plans of the East Asian Very Long Baseline Interferometry (VLBI) community for pushing the frontiers of millimeter/submillimeter VLBI. To this end, we shall endeavor to actively promote coordinated efforts in the East Asia region. Our goal is to establish firm collaborations among the East Asia VLBI community in partnership with related institutes in North America and Europe and to expand existing global mm/submm VLBI arrays for (a) exploring the vicinity of black holes with an ultimate angular resolution down to 1~R$_{S}$ (Schwarzschild radius) and (b) investigating the dynamics of circumstellar gas in star-forming regions and late-type stars, and circumnuclear gas around active galactic nuclei (AGNs). In the first half of this White Paper, we highlight scientific accomplishments of the East Asia (EA) VLBI community. Various VLBI research results on M87, Sgr A*, blazars, narrow-line Seyfert~1 galaxies, and compact symmetric objects are described, and future visions of our VLBI science are briefly presented. Maser science of star formation, stellar evolution, and physics of accretion disks around AGNs are also discussed. A new vision for conducting multi-transition maser studies using mm/submm VLBI together with the Atacama Large Millimeter/Submillimeter Array (ALMA) is described. In the second half of this White Paper, we describe the EA community's vision for using mm/submm VLBI arrays in the framework or extended version of the Event Horizon Telescope (EHT) and the Global Millimeter VLBI Array (GMVA). The accomplishment of the aforementioned goal will maximize the overall scientific outcomes of mm/submm VLBI in the world. △ Less

Submitted 12 May, 2017; originally announced May 2017.

Comments: 93 pages, 54 figures. White Paper from East Asia VLBI community

The shortest periodic and flaring flux variability of a methanol maser emission at 6.7 GHz in G 014.23-00.50

Authors: Koichiro Sugiyama, Katsura Nagase, Yoshinori Yonekura, Munetake Momose, Yasutaka Yasui, Yu Saito, Kazuhito Motogi, Mareki Honma, Kazuya Hachisuka, Naoko Matsumoto, Mizuho Uchiyama, Kenta Fujisawa

Abstract: We detected flaring flux variability that regularly occurred with the period of 23.9 days on a 6.7 GHz methanol maser emission at Vlsr = 25.30 km/s in G 014.23-00.50 through highly frequent monitoring using the Hitachi 32-m radio telescope. By analyzing data from 05 January 2013 to 21 January 2016, the periodic variability has persisted in at least 47 cycles, corresponding to approximately 1,100 d… ▽ More We detected flaring flux variability that regularly occurred with the period of 23.9 days on a 6.7 GHz methanol maser emission at Vlsr = 25.30 km/s in G 014.23-00.50 through highly frequent monitoring using the Hitachi 32-m radio telescope. By analyzing data from 05 January 2013 to 21 January 2016, the periodic variability has persisted in at least 47 cycles, corresponding to approximately 1,100 days. The period of 23.9 days is the shortest one observed in masers at around high-mass young stellar objects so far. The flaring component normally falls below the detection limit (3 sigma) of 0.9 Jy. In the flaring periods, the component rises above the detection limit with the ratio of the peak flux density more than 180 in comparison with a quiescent phase, showing intermittent periodic variability. The time-scale of the flux rise was typically two days or shorter, and both symmetric and asymmetric profiles of flux variability were observed through intraday monitoring. These characteristics might be explained by a change in the flux of seed photons by a colliding-wind binary (CWB) system or a variation of the dust temperature by an extra heating source of a shock formed by the CWB system within a gap region in a circumbinary disk, in which the orbital semi-major axes of the binary are 0.26-0.34 au. △ Less

Submitted 10 April, 2017; originally announced April 2017.

Comments: 12 pages, 9 figures, 2 tables, accepted for publication in PASJ

J-GEM Follow-Up Observations of The Gravitational Wave Source GW151226

Authors: Michitoshi Yoshida, Yousuke Utsumi, Nozomu Tominaga, Tomoki Morokuma, Masaomi Tanaka, Yuichiro Asakura, Kazuya Matsubayashi, Kouji Ohta, Fumio Abe, Sho Chimasu, Hisanori Furusawa, Ryosuke Itoh, Yoichi Itoh, Yuka Kanda, Koji S. Kawabata, Miho Kawabata, Shintaro Koshida, Naoki Koshimoto, Daisuke Kuroda, Yuki Moritani, Kentaro Motohara, Katsuhiro L. Murata, Takahiro Nagayama, Tatsuya Nakaoka, Fumiaki Nakata , et al. (13 additional authors not shown)

Abstract: We report the results of optical--infrared follow-up observations of the gravitational wave (GW) event GW151226 detected by the Advanced LIGO in the framework of J-GEM (Japanese collaboration for Gravitational wave ElectroMagnetic follow-up). We performed wide-field optical imaging surveys with Kiso Wide Field Camera (KWFC), Hyper Suprime-Cam (HSC), and MOA-cam3. The KWFC survey started at 2.26 da… ▽ More We report the results of optical--infrared follow-up observations of the gravitational wave (GW) event GW151226 detected by the Advanced LIGO in the framework of J-GEM (Japanese collaboration for Gravitational wave ElectroMagnetic follow-up). We performed wide-field optical imaging surveys with Kiso Wide Field Camera (KWFC), Hyper Suprime-Cam (HSC), and MOA-cam3. The KWFC survey started at 2.26 days after the GW event and covered 778 deg$^2$ centered at the high Galactic region of the skymap of GW151226. We started the HSC follow-up observations from 12 days after the event and covered an area of 63.5 deg$^2$ of the highest probability region of the northern sky with the limiting magnitudes of 24.6 and 23.8 for i band and z band, respectively. MOA-cam3 covered 145 deg$^2$ of the skymap with MOA-red filter 2.5 months after the GW alert. Total area covered by the wide-field surveys was 986.5 deg$^2$. The integrated detection probability of all the observed area was $\sim$29%. We also performed galaxy-targeted observations with six optical and near-infrared telescopes from 1.61 days after the event. Total of 238 nearby (<100 Mpc) galaxies were observed with the typical I band limiting magnitude of $\sim$19.5. We detected 13 supernova candidates with the KWFC survey, and 60 extragalactic transients with the HSC survey. Two third of the HSC transients were likely supernovae and the remaining one third were possible active galactic nuclei. With our observational campaign, we found no transients that are likely to be associated with GW151226. △ Less

Submitted 4 November, 2016; originally announced November 2016.

Comments: 12 pages, 6 figures, 4 tables, accepted for publication in PASJ

Rotational equilibria by Lagrangian variational principle: toward multi-dimensional stellar evolutions

Authors: Nobutoshi Yasutake, Kotaro Fujisawa, Shoichi Yamada

Abstract: We have developed a new formulation to obtain self-gravitating, axisymmetric configurations in permanent rotation. The formulation is based on the Lagrangian variational principle with a triangulated mesh. It treats not only barotropic but also baroclinic equations of state. We compare the various stellar equilibria obtained by our new scheme with those by Hachisu's self-consistent field scheme fo… ▽ More We have developed a new formulation to obtain self-gravitating, axisymmetric configurations in permanent rotation. The formulation is based on the Lagrangian variational principle with a triangulated mesh. It treats not only barotropic but also baroclinic equations of state. We compare the various stellar equilibria obtained by our new scheme with those by Hachisu's self-consistent field scheme for the barotropic case, and those by Fujisawa's self-consistent field scheme for the baroclinic case. Included in these rotational configurations are those with shellular-type rotations, which are commonly assumed in the evolution calculation of rotating stars. Although radiation processes, convections and meridional flows have not been taken into account in this study, we have in mind the application of this method to the two-dimensional evolution calculations of rotating stars, for which the Lagrangian formulation is best suited. △ Less

Submitted 3 September, 2016; originally announced September 2016.

Comments: 23 pages, 28 figures, accepted for publishing in MNRAS

J-GEM Follow-Up Observations to Search for an Optical Counterpart of The First Gravitational Wave Source GW150914

Authors: Tomoki Morokuma, Masaomi Tanaka, Yuichiro Asakura, Fumio Abe, Paul J. Tristram, Yousuke Utsumi, Mamoru Doi, Kenta Fujisawa, Ryosuke Itoh, Yoichi Itoh, Koji S. Kawabata, Nobuyuki Kawai, Daisuke Kuroda, Kazuya Matsubayashi, Kentaro Motohara, Katsuhiro L. Murata, Takahiro Nagayama, Kouji Ohta, Yoshihiko Saito, Yoichi Tamura, Nozomu Tominaga, Makoto Uemura, Kenshi Yanagisawa, Yoichi Yatsu, Michitoshi Yoshida

Abstract: We present our optical follow-up observations to search for an electromagnetic counterpart of the first gravitational wave source GW150914 in the framework of the Japanese collaboration for Gravitational wave ElectroMagnetic follow-up (J-GEM), which is an observing group utilizing optical and radio telescopes in Japan, as well as those in New Zealand, China, South Africa, Chile, and Hawaii. We car… ▽ More We present our optical follow-up observations to search for an electromagnetic counterpart of the first gravitational wave source GW150914 in the framework of the Japanese collaboration for Gravitational wave ElectroMagnetic follow-up (J-GEM), which is an observing group utilizing optical and radio telescopes in Japan, as well as those in New Zealand, China, South Africa, Chile, and Hawaii. We carried out a wide-field imaging survey with Kiso Wide Field Camera (KWFC) on the 1.05-m Kiso Schmidt telescope in Japan and a galaxy-targeted survey with Tripole5 on the B&C 61-cm telescope in New Zealand. Approximately 24 deg2 regions in total were surveyed in i-band with KWFC and 18 nearby galaxies were observed with Tripole5 in g-, r-, and i-bands 4-12 days after the gravitational wave detection. Median 5-sigma depths are i~18.9 mag for the KWFC data and g~18.9 mag, r~18.7 mag, and i~18.3 mag for the Tripole5 data. Probability for a counterpart to be in the observed area is 1.2% in the initial skymap and 0.1% in the final skymap. We do not find any transient source associated to an external galaxy with spatial offset from its center, which is consistent with the local supernova rate. We summarize future prospects and ongoing efforts to pin down electromagnetic counterparts of binary black hole mergers as well as neutron star mergers. △ Less

Submitted 10 May, 2016; originally announced May 2016.

Comments: submitted to PASJ, 7 pages, 1 figure, 4 tables

Supplement: Localization and broadband follow-up of the gravitational-wave transient GW150914

Authors: B. P. Abbott, R. Abbott, T. D. Abbott, M. R. Abernathy, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, A. Allocca, P. A. Altin, S. B. Anderson, W. G. Anderson, K. Arai , et al. (1522 additional authors not shown)

Abstract: This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the dif… ▽ More This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands. △ Less

Submitted 21 July, 2016; v1 submitted 26 April, 2016; originally announced April 2016.

Comments: For the main Letter, see arXiv:1602.08492

Report number: LIGO-P1600137-v2

Journal ref: The Astrophysical Journal Supplement Series, 225:8 (15pp), 2016 July

Localization and broadband follow-up of the gravitational-wave transient GW150914

Authors: B. P. Abbott, R. Abbott, T. D. Abbott, M. R. Abernathy, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, A. Allocca, P. A. Altin, S. B. Anderson, W. G. Anderson, K. Arai , et al. (1522 additional authors not shown)

Abstract: A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared wit… ▽ More A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams. △ Less

Submitted 21 July, 2016; v1 submitted 26 February, 2016; originally announced February 2016.

Comments: For Supplement, see https://arxiv.org/abs/1604.07864

Report number: LIGO-P1500227-v12

Journal ref: The Astrophysical Journal Letters, 826:L13 (8pp), 2016 July 20

Hyperbolic Self-Gravity Solver for Large Scale Hydrodynamical Simulations

Authors: Ryosuke Hirai, Hiroki Nagakura, Hirotada Okawa, Kotaro Fujisawa

Abstract: A new computationally efficient method has been introduced to treat self-gravity in mesh based hydrodynamical simulations. It is applied simply by slightly modifying the Poisson equation into an inhomogeneous wave equation. This roughly corresponds to the weak field limit of the Einstein equations in general relativity, and as long as the gravitation propagation speed is taken to be larger than th… ▽ More A new computationally efficient method has been introduced to treat self-gravity in mesh based hydrodynamical simulations. It is applied simply by slightly modifying the Poisson equation into an inhomogeneous wave equation. This roughly corresponds to the weak field limit of the Einstein equations in general relativity, and as long as the gravitation propagation speed is taken to be larger than the hydrodynamical characteristic speed, the results agree with solutions for the Poisson equation. The solutions almost perfectly agree if the domain is taken large enough, or appropriate boundary conditions are given. Our new method can not only significantly reduce the computational time compared with existent methods, but is also fully compatible with massive parallel computation, nested grids and adaptive mesh refinement techniques, all of which can accelerate the progress in computational astrophysics and cosmology. △ Less

Submitted 29 March, 2016; v1 submitted 18 February, 2016; originally announced February 2016.

Comments: 7 pages, 6 figures, accepted for publication in PRD

Journal ref: Phys. Rev. D 93, 083006 (2016)

The East-Asian VLBI Network

Authors: Kiyoaki Wajima, Yoshiaki Hagiwara, Tao An, Willem A. Baan, Kenta Fujisawa, Longfei Hao, Wu Jiang, Taehyun Jung, Noriyuki Kawaguchi, Jongsoo Kim, Hideyuki Kobayashi, Se-Jin Oh, Duk-Gyoo Roh, Min Wang, Yuanwei Wu, Bo Xia, Ming Zhang

Abstract: The East-Asian VLBI Network (EAVN) is the international VLBI facility in East Asia and is conducted in collaboration with China, Japan, and Korea. The EAVN consists of VLBI arrays operated in each East Asian country, containing 21 radio telescopes and three correlators. The EAVN will be mainly operated at 6.7 (C-band), 8 (X-band), 22 (K-band), and 43 GHz (Q-band), although the EAVN has an ability… ▽ More The East-Asian VLBI Network (EAVN) is the international VLBI facility in East Asia and is conducted in collaboration with China, Japan, and Korea. The EAVN consists of VLBI arrays operated in each East Asian country, containing 21 radio telescopes and three correlators. The EAVN will be mainly operated at 6.7 (C-band), 8 (X-band), 22 (K-band), and 43 GHz (Q-band), although the EAVN has an ability to conduct observations at 1.6 - 129 GHz. We have conducted fringe test observations eight times to date at 8 and 22 GHz and fringes have been successfully detected at both frequencies. We have also conducted science commissioning observations of 6.7 GHz methanol masers in massive star-forming regions. The EAVN will be operational from the second half of 2017, providing complementary results with the FAST on AGNs, massive star-forming regions, and evolved stars with high angular resolution at cm- to mm-wavelengths. △ Less

Submitted 11 December, 2015; originally announced December 2015.

Comments: 6 pages, 3 figures, 2 tables. To appear in the proceedings of "Frontiers in Radio Astronomy and FAST Early Sciences Symposium 2015" ed. Lei Qian (ASP Conf. Ser.)

A versatile numerical method for obtaining structures of rapidly rotating baroclinic stars: self-consistent and systematic solutions with shellular-type rotation

Authors: Kotaro Fujisawa

Abstract: This paper develops a novel numerical method for obtaining structures of rapidly rotating stars based on a self-consistent field scheme. The solution is obtained iteratively. Both rapidly rotating barotropic and baroclinic equilibrium states are calculated self-consistently using this method. Two types of rotating baroclinic stars are investigated by changing the isentropic surfaces inside the sta… ▽ More This paper develops a novel numerical method for obtaining structures of rapidly rotating stars based on a self-consistent field scheme. The solution is obtained iteratively. Both rapidly rotating barotropic and baroclinic equilibrium states are calculated self-consistently using this method. Two types of rotating baroclinic stars are investigated by changing the isentropic surfaces inside the star. Solution sequences of these are calculated systematically and critical rotation models beyond which no rotating equilibrium state exists are also obtained. All of these rotating baroclinic stars satisfy necessarily the Bjerknes-Rosseland rules. Self-consistent solutions of baro-clinic stars with shellular-type rotation are successfully obtained where the isentropic surfaces are oblate and the surface temperature is hotter at the poles than at the equator if it is assumed that the star is an ideal gas star. These are the first self-consistent and systematic solutions of rapidly rotating baroclinic stars with shellular-type rotations. Since they satisfy the stability criterion due to their rapid rotation, these rotating baroclinic stars would be dynamically stable. This novel numerical method and the solutions of the rapidly rotating baroclinic stars will be useful for investigating stellar evolution with rapid rotations. △ Less

Submitted 19 October, 2015; v1 submitted 9 July, 2015; originally announced July 2015.

Comments: 14 pages, 6 figures; published MNRAS 454, 3060--3072 (2015)

Journal ref: 2015MNRAS.454.3060F

Magnetized stars with differential rotation and a differential toroidal field

Authors: Kotaro Fujisawa

Abstract: We have succeeded in obtaining magnetized equilibrium states with differential rotation and differential toroidal magnetic fields. If an internal toroidal field of a proto-neutron star is wound up from the initial poloidal magnetic field by differential rotation, the distribution of the toroidal magnetic field is determined by the profile of this differential rotation. However, the distributions o… ▽ More We have succeeded in obtaining magnetized equilibrium states with differential rotation and differential toroidal magnetic fields. If an internal toroidal field of a proto-neutron star is wound up from the initial poloidal magnetic field by differential rotation, the distribution of the toroidal magnetic field is determined by the profile of this differential rotation. However, the distributions of the toroidal fields in all previous magnetized equilibrium studies do not represent the magnetic winding by the differential rotation of the star. In this paper, we investigate a formulation of a differential toroidal magnetic field that represents the magnetic field wound up by differential rotation. We have developed two functional forms of differential toroidal fields which correspond to a v-constant and a j-constant field in analogy to differential rotations. As the degree of the differential becomes very high, the toroidal magnetic field becomes highly localized and concentrated near the rotational axis. Such a differential toroidal magnetic field would suppress the low-T/|W| instability more efficiently even if the total magnetic field energy is much smaller than that of a non-differential toroidal magnetic field. △ Less

Submitted 27 August, 2015; v1 submitted 22 April, 2015; originally announced April 2015.

Comments: 10 pages, 4 figures; published in MNRAS

Journal ref: 2015MNRAS.450.4016F

Appearance of the prolate and the toroidal magnetic field dominated stars - Analytic approach

Authors: Kotaro Fujisawa, Yoshiharu Eriguchi

Abstract: We have analyzed magnetized equilibrium states and showed a condition for appearance of the prolate and the toroidal magnetic field dominated stars by analytic approaches. Both observations and numerical stability analysis support that the magnetized star would have the prolate and the large internal toroidal magnetic fields. In this context, many investigations concerning magnetized equilibrium s… ▽ More We have analyzed magnetized equilibrium states and showed a condition for appearance of the prolate and the toroidal magnetic field dominated stars by analytic approaches. Both observations and numerical stability analysis support that the magnetized star would have the prolate and the large internal toroidal magnetic fields. In this context, many investigations concerning magnetized equilibrium states have tried to obtain the prolate and the toroidal dominant solutions, but many of them have failed to obtain such configurations. Since the Lorentz force is a cross product of current density and magnetic field, the prolate shaped configurations and the large toroidal magnetic fields in stars require a special relation between current density and the Lorentz force. We have analyzed simple analytical solutions and found that the prolate and the toroidal dominant configuration require non force-free toroidal current density that flows in the opposite direction with respect to the bulk current within the star. Such current density results in the Lorentz force which makes the stellar shape prolate. Satisfying this special relation between the current density and the Lorentz force is a key for appearance of the prolate and the toroidal magnetic field dominated magnetized star. △ Less

Submitted 2 April, 2015; originally announced April 2015.

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

First Parallax Measurements Towards a 6.7 GHz Methanol Maser with the Australian Long Baseline Array - Distance to G339.884-1.259

Authors: V. Krishnan, S. P. Ellingsen, M. J. Reid, A. Brunthaler, A. Sanna, J. McCallum, C. Reynolds, H. E. Bignall, C. J. Phillips, R. Dodson, M. Rioja, J. L. Caswell, X. Chen, J. R. Dawson, K. Fujisawa, S. Goedhart, J. A. Green, K. Hachisuka, M. Honma, K. Menten, Z. Q. Shen, M. A. Voronkov, A. J. Walsh, Y. Xu, B. Zhang , et al. (1 additional authors not shown)

Abstract: We have conducted the first parallax and proper motion measurements of 6.7 GHz methanol maser emission using the Australian Long Baseline Array (LBA). The parallax of G339.884$-$1.259 measured from five epochs of observations is 0.48$\pm $0.08 mas, corresponding to a distance of $2.1^{+0.4}_{-0.3}$ kpc, placing it in the Scutum spiral arm. This is consistent (within the combined uncertainty) with… ▽ More We have conducted the first parallax and proper motion measurements of 6.7 GHz methanol maser emission using the Australian Long Baseline Array (LBA). The parallax of G339.884$-$1.259 measured from five epochs of observations is 0.48$\pm $0.08 mas, corresponding to a distance of $2.1^{+0.4}_{-0.3}$ kpc, placing it in the Scutum spiral arm. This is consistent (within the combined uncertainty) with the kinematic distance estimate for this source at 2.5$\pm $0.5 kpc using the latest Solar and Galactic rotation parameters. We find from the Lyman continuum photon flux that the embedded core of the young star is of spectral type B1, demonstrating that luminous 6.7 GHz methanol masers can be associated with high-mass stars towards the lower end of the mass range. △ Less

Submitted 19 March, 2015; originally announced March 2015.

Accelerating a water maser face-on jet from a high mass young stellar object

Authors: Kazuhito Motogi, Kazuo Sorai, Mareki Honma, Tomoya Hirota, Kazuya Hachisuka, Kotaro Niinuma, Koichiro Sugiyama, Yoshinori Yonekura, Kenta Fujisawa

Abstract: We report on a long-term single-dish and VLBI monitoring for intermittent flare activities of a Dominant Blue-Shifted H$_{2}$O Maser (DBSM) associated with a southern high mass young stellar object, G353.273+0.641. Bi-weekly single-dish monitoring using Hokkaido University Tomakomai 11-m radio telescope has shown that a systematic acceleration continues over four years beyond a lifetime of individ… ▽ More We report on a long-term single-dish and VLBI monitoring for intermittent flare activities of a Dominant Blue-Shifted H$_{2}$O Maser (DBSM) associated with a southern high mass young stellar object, G353.273+0.641. Bi-weekly single-dish monitoring using Hokkaido University Tomakomai 11-m radio telescope has shown that a systematic acceleration continues over four years beyond a lifetime of individual maser features. This fact suggests that the H$_{2}$O maser traces a region where molecular gas is steadily accelerated. There were five maser flares during five-years monitoring, and maser distributions in four of them were densely monitored by the VLBI Exploration of Radio Astrometry (VERA). The overall distribution of the maser features suggests the presence of a bipolar jet, with the 3D kinematics indicating that it is almost face-on (inclination angle of $\sim$ 8$^{\fdg}$--17$^{\fdg}$ from the line-of-sight). Most of maser features were recurrently excited within a region of 100$\times$100 AU$^{2}$ around the radio continuum peak, while their spatial distributions significantly varied between each flare. This confirms that episodic propagations of outflow shocks recurrently invoke intermittent flare activities. We also measured annual parallax, deriving the source distance of 1.70 $^{+0.19}_{-0.16}$ kpc that is consistent with the commonly-used photometric distance. △ Less

Submitted 10 February, 2015; v1 submitted 2 February, 2015; originally announced February 2015.

Comments: Accepted to PASJ, 17 pages, 12 figure

Equilibrium solutions of relativistic rotating stars with mixed poloidal and toroidal magnetic fields

Authors: Koji Uryu, Eric Gourgoulhon, Charalampos Markakis, Kotaro Fujisawa, Antonios Tsokaros, Yoshiharu Eriguchi

Abstract: Stationary and axisymmetric solutions of relativistic rotating stars with strong mixed poloidal and toroidal magnetic fields are obtained numerically. Because of the mixed components of the magnetic field, the underlying stationary and axisymmetric spacetimes are no longer circular. These configurations are computed from the full set of the Einstein-Maxwell equations, Maxwell's equations and from… ▽ More Stationary and axisymmetric solutions of relativistic rotating stars with strong mixed poloidal and toroidal magnetic fields are obtained numerically. Because of the mixed components of the magnetic field, the underlying stationary and axisymmetric spacetimes are no longer circular. These configurations are computed from the full set of the Einstein-Maxwell equations, Maxwell's equations and from first integrals and integrability conditions of the magnetohydrodynamic equilibrium equations. After a brief introduction of the formulation of the problem, we present the first results for highly deformed magnetized rotating compact stars. △ Less

Submitted 14 October, 2014; originally announced October 2014.

Comments: 7 pages, to appear in PRD rapid communication

Report number: DPUR-TH-42

Journal ref: Phys. Rev. D 90, 101501(R) 2014