Showing 1–50 of 114 results for author: Sasada, M

X-ray and multiwavelength polarization of Mrk 501 from 2022 to 2023

Authors: Chien-Ting J. Chen, Ioannis Liodakis, Riccardo Middei, Dawoon E. Kim, Laura Di Gesu, Alessandro Di Marco, Steven R. Ehlert, Manel Errando, Michela Negro, Svetlana G. Jorstad, Alan P. Marscher, Kinwah Wu, Iván Agudo, Juri Poutanen, Tsunefumi Mizuno, Pouya M. Kouch, Elina Lindfors, George A. Borman, Tatiana S. Grishina, Evgenia N. Kopatskaya, Elena G. Larionova, Daria A. Morozova, Sergey S. Savchenko, Ivan S. Troitsky, Yulia V. Troitskaya , et al. (121 additional authors not shown)

Abstract: We present multiwavelength polarization measurements of the luminous blazar Mrk~501 over a 14-month period. The 2--8 keV X-ray polarization was measured with the Imaging X-ray Polarimetry Explorer (IXPE) with six 100-ks observations spanning from 2022 March to 2023 April. Each IXPE observation was accompanied by simultaneous X-ray data from NuSTAR, Swift/XRT, and/or XMM-Newton. Complementary optic… ▽ More We present multiwavelength polarization measurements of the luminous blazar Mrk~501 over a 14-month period. The 2--8 keV X-ray polarization was measured with the Imaging X-ray Polarimetry Explorer (IXPE) with six 100-ks observations spanning from 2022 March to 2023 April. Each IXPE observation was accompanied by simultaneous X-ray data from NuSTAR, Swift/XRT, and/or XMM-Newton. Complementary optical-infrared polarization measurements were also available in the B, V, R, I, and J bands, as were radio polarization measurements from 4.85 GHz to 225.5 GHz. Among the first five IXPE observations, we did not find significant variability in the X-ray polarization degree and angle with IXPE. However, the most recent sixth observation found an elevated polarization degree at $>3σ$ above the average of the other five observations. The optical and radio measurements show no apparent correlations with the X-ray polarization properties. Throughout the six IXPE observations, the X-ray polarization degree remained higher than, or similar to, the R-band optical polarization degree, which remained higher than the radio value. This is consistent with the energy-stratified shock scenario proposed to explain the first two IXPE observations, in which the polarized X-ray, optical, and radio emission arises from different regions. △ Less

Submitted 15 July, 2024; originally announced July 2024.

Comments: Accepted for publication in ApJ

Testing particle acceleration in blazar jets with continuous high-cadence optical polarization observations

Authors: Ioannis Liodakis, Sebastian Kiehlmann, Alan P. Marscher, Haocheng Zhang, Dmitry Blinov, Svetlana G. Jorstad, Iván Agudo, Erika Benítez, Andrei Berdyugin, Giacomo Bonnoli, Carolina Casadio, Chien-Ting Chen, Wen-Ping Chen, Steven R. Ehlert, Juan Escudero, Tatiana S. Grishina, David Hiriart, Angela Hsu, Ryo Imazawa, Helen E. Jermak, Jincen Jose, Philip Kaaret, Evgenia N. Kopatskaya, Bhavana Lalchand, Elena G. Larionova , et al. (22 additional authors not shown)

Abstract: Variability can be the pathway to understanding the physical processes in astrophysical jets, however, the high-cadence observations required to test particle acceleration models are still missing. Here we report on the first attempt to produce continuous, >24 hour polarization light curves of blazars using telescopes distributed across the globe and the rotation of the Earth to avoid the rising S… ▽ More Variability can be the pathway to understanding the physical processes in astrophysical jets, however, the high-cadence observations required to test particle acceleration models are still missing. Here we report on the first attempt to produce continuous, >24 hour polarization light curves of blazars using telescopes distributed across the globe and the rotation of the Earth to avoid the rising Sun. Our campaign involved 16 telescopes in Asia, Europe, and North America. We observed BL Lacertae and CGRaBS J0211+1051 for a combined 685 telescope hours. We find large variations in the polarization degree and angle for both sources in sub-hour timescales as well as a ~180 degree rotation of the polarization angle in CGRaBS J0211+1051 in less than two days. We compared our high-cadence observations to Particle-In-Cell magnetic reconnection and turbulent plasma simulations. We find that although the state of the art simulation frameworks can produce a large fraction of the polarization properties, they do not account for the entirety of the observed polarization behavior in blazar jets. △ Less

Submitted 21 June, 2024; originally announced June 2024.

Comments: 20 pages, 15 figures, 2 tables, accepted for publication in A&A. The data used in the paper are available here: https://doi.org/10.7910/DVN/IETSXS

The Japanese Vision for the Black Hole Explorer Mission

Authors: Kazunori Akiyama, Kotaro Niinuma, Kazuhiro Hada, Akihiro Doi, Yoshiaki Hagiwara, Aya E. Higuchi, Mareki Honma, Tomohisa Kawashima, Dimitar Kolev, Shoko Koyama, Sho Masui, Ken Ohsuga, Hidetoshi Sano, Hideki Takami, Yuh Tsunetoe, Yoshinori Uzawa, Takuya Akahori, Yuto Akiyama, Peter Galison, Takayuki J. Hayashi, Tomoya Hirota, Makoto Inoue, Yuhei Iwata, Michael D. Johnson, Motoki Kino , et al. (21 additional authors not shown)

Abstract: The Black Hole Explorer (BHEX) is a next-generation space very long baseline interferometry (VLBI) mission concept that will extend the ground-based millimeter/submillimeter arrays into space. The mission, closely aligned with the science priorities of the Japanese VLBI community, involves an active engagement of this community in the development of the mission, resulting in the formation of the B… ▽ More The Black Hole Explorer (BHEX) is a next-generation space very long baseline interferometry (VLBI) mission concept that will extend the ground-based millimeter/submillimeter arrays into space. The mission, closely aligned with the science priorities of the Japanese VLBI community, involves an active engagement of this community in the development of the mission, resulting in the formation of the Black Hole Explorer Japan Consortium. Here we present the current Japanese vision for the mission, ranging from scientific objectives to instrumentation. The Consortium anticipates a wide range of scientific investigations, from diverse black hole physics and astrophysics studied through the primary VLBI mode, to the molecular universe explored via a potential single-dish observation mode in the previously unexplored 50-70\,GHz band that would make BHEX the highest-sensitivity explorer ever of molecular oxygen. A potential major contribution for the onboard instrument involves supplying essential elements for its high-sensitivity dual-band receiving system, which includes a broadband 300\,GHz SIS mixer and a space-certified multi-stage 4.5K cryocooler akin to those used in the Hitomi and XRISM satellites by the Japan Aerospace Exploration Agency. Additionally, the Consortium explores enhancing and supporting BHEX operations through the use of millimeter/submillimeter facilities developed by the National Astronomical Observatory of Japan, coupled with a network of laser communication stations operated by the National Institute of Information and Communication Technology. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: Submitted to SPIE Conference Series, 17 pages for the main text, 33 pages for the references, 4 figures, SPIE Astronomical Telescopes + Instrumentation 2024

IXPE observation of PKS 2155-304 reveals the most highly polarized blazar

Authors: Pouya M. Kouch, Ioannis Liodakis, Riccardo Middei, Dawoon E. Kim, Fabrizio Tavecchio, Alan P. Marscher, Herman L. Marshall, Steven R. Ehlert, Laura Di Gesu, Svetlana G. Jorstad, Iván Agudo, Grzegorz M. Madejski, Roger W. Romani, Manel Errando, Elina Lindfors, Kari Nilsson, Ella Toppari, Stephen B. Potter, Ryo Imazawa, Mahito Sasada, Yasushi Fukazawa, Koji S. Kawabata, Makoto Uemura, Tsunefumi Mizuno, Tatsuya Nakaoka , et al. (111 additional authors not shown)

Abstract: We report the X-ray polarization properties of the high-synchrotron-peaked (HSP) blazar PKS 2155$-$304 based on observations with the Imaging X-ray Polarimetry Explorer (IXPE). We observed the source between Oct 27 and Nov 7, 2023. We also conducted an extensive contemporaneous multiwavelength (MW) campaign. We find that during the first half ($T_1$) of the IXPE pointing, the source exhibited the… ▽ More We report the X-ray polarization properties of the high-synchrotron-peaked (HSP) blazar PKS 2155$-$304 based on observations with the Imaging X-ray Polarimetry Explorer (IXPE). We observed the source between Oct 27 and Nov 7, 2023. We also conducted an extensive contemporaneous multiwavelength (MW) campaign. We find that during the first half ($T_1$) of the IXPE pointing, the source exhibited the highest X-ray polarization degree detected for an HSP blazar thus far, (30.7$\pm$2.0)%, which dropped to (15.3$\pm$2.1)% during the second half ($T_2$). The X-ray polarization angle remained stable during the IXPE pointing at 129.4$^\circ$$\pm$1.8$^\circ$ and 125.4$^\circ$$\pm$3.9$^\circ$ during $T_1$ and $T_2$, respectively. Meanwhile, the optical polarization degree remained stable during the IXPE pointing, with average host-galaxy-corrected values of (4.3$\pm$0.7)% and (3.8$\pm$0.9)% during the $T_1$ and $T_2$, respectively. During the IXPE pointing, the optical polarization angle changed achromatically from $\sim$140$^\circ$ to $\sim$90$^\circ$ and back to $\sim$130$^\circ$. Despite several attempts, we only detected (99.7% conf.) the radio polarization once (during $T_2$, at 225.5 GHz): with degree (1.7$\pm$0.4)% and angle 112.5$^\circ$$\pm$5.5$^\circ$. The direction of the broad pc-scale jet is rather ambiguous and has been found to point to the east and south at different epochs; however, on larger scales (> 1.5 pc) the jet points toward the southeast ($\sim$135$^\circ$), similar to all of the MW polarization angles. Moreover, the X-ray to optical polarization degree ratios of $\sim$7 and $\sim$4 during $T_1$ and $T_2$, respectively, are similar to previous IXPE results for several HSP blazars. These findings, combined with the lack of correlation of temporal variability between the MW polarization properties, agree with an energy-stratified shock-acceleration scenario in HSP blazars. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: 17 pages, 10 figures, 4 tables, Accepted for publication in A&A

Characterization of Gradient Condition for Asymmetric Partial Exclusion Processes and Their Scaling Limits

Authors: Patrícia Gonçalves, Kohei Hayashi, Makiko Sasada

Abstract: We consider partial exclusion processes~(PEPs) on the one-dimensional square lattice, that is, a system of interacting particles where each particle random walks according to a jump rate satisfying an exclusion rule that allows up to a certain number of particles can exist on each site. Particularly, we assume that the jump rate is given as a product of two functions depending on occupation variab… ▽ More We consider partial exclusion processes~(PEPs) on the one-dimensional square lattice, that is, a system of interacting particles where each particle random walks according to a jump rate satisfying an exclusion rule that allows up to a certain number of particles can exist on each site. Particularly, we assume that the jump rate is given as a product of two functions depending on occupation variables on the original and target sites. Our interest is to study the limiting behavior, especially to derive some macroscopic PDEs by means of (fluctuating) hydrodynamics, of fluctuation fields associated with PEPs, starting from an invariant measure. The so-called gradient condition, meaning that the symmetric part of the instantaneous current is written in a gradient form, and that the invariant measures are given as a product measure is technically crucial. Our first main result is to clarify the relationship between these two conditions, and we show that the gradient condition and the existence of product invariant measures are mutually equivalent, provided the jump rate is given in the above simple form, as it is imposed in most of the literature, and the dynamics is asymmetric. Moreover, when the width of the lattice tends to zero and the process is accelerated in diffusive time-scaling, we show that the family of fluctuation fields converges to the stationary energy solution of the stochastic Burgers equation (SBE), under the setting that the jump rate to the right neighboring site is a bit larger than the one to the left side, of which discrepancy is given as square root of the width of the underlying lattice. This fills the gap at the level of universality of SBE since it has been proved for the exclusion process (a special case of PEP) and for the zero-range process. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 25 pages

MSC Class: 60H15; 60K35; 82B44

Ordered magnetic fields around the 3C 84 central black hole

Authors: G. F. Paraschos, J. -Y. Kim, M. Wielgus, J. Röder, T. P. Krichbaum, E. Ros, I. Agudo, I. Myserlis, M. Moscibrodzka, E. Traianou, J. A. Zensus, L. Blackburn, C. -K. Chan, S. Issaoun, M. Janssen, M. D. Johnson, V. L. Fish, K. Akiyama, A. Alberdi, W. Alef, J. C. Algaba, R. Anantua, K. Asada, R. Azulay, U. Bach , et al. (258 additional authors not shown)

Abstract: 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures a… ▽ More 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures and understand the physical conditions in the compact region of 3C84. We used EHT 228GHz observations and, given the limited (u,v)-coverage, applied geometric model fitting to the data. We also employed quasi-simultaneously observed, multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure. We report the detection of a highly ordered, strong magnetic field around the central, SMBH of 3C84. The brightness temperature analysis suggests that the system is in equipartition. We determined a turnover frequency of $ν_m=(113\pm4)$GHz, a corresponding synchrotron self-absorbed magnetic field of $B_{SSA}=(2.9\pm1.6)$G, and an equipartition magnetic field of $B_{eq}=(5.2\pm0.6)$G. Three components are resolved with the highest fractional polarisation detected for this object ($m_\textrm{net}=(17.0\pm3.9)$%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017-2018. We report a steeply negative slope of the spectrum at 228GHz. We used these findings to test models of jet formation, propagation, and Faraday rotation in 3C84. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C84. However, systematic uncertainties due to the limited (u,v)-coverage, however, cannot be ignored. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 15 pages, 6 figures, published in A&A

Journal ref: Issue: A&A Volume 682, February 2024; Article number: L3; Number of pages: 15

Insights into the broad-band emission of the TeV blazar Mrk 501 during the first X-ray polarization measurements

Authors: S. Abe, J. Abhir, V. A. Acciari, A. Aguasca-Cabot, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, K. Asano, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, A. Bautista, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, J. Bernete, A. Berti, J. Besenrieder , et al. (239 additional authors not shown)

Abstract: We present the first multi-wavelength study of Mrk 501 including very-high-energy (VHE) gamma-ray observations simultaneous to X-ray polarization measurements from the Imaging X-ray Polarimetry Explorer (IXPE). We use radio-to-VHE data from a multi-wavelength campaign organized between 2022-03-01 and 2022-07-19. The observations were performed by MAGIC, Fermi-LAT, NuSTAR, Swift (XRT and UVOT), and… ▽ More We present the first multi-wavelength study of Mrk 501 including very-high-energy (VHE) gamma-ray observations simultaneous to X-ray polarization measurements from the Imaging X-ray Polarimetry Explorer (IXPE). We use radio-to-VHE data from a multi-wavelength campaign organized between 2022-03-01 and 2022-07-19. The observations were performed by MAGIC, Fermi-LAT, NuSTAR, Swift (XRT and UVOT), and several instruments covering the optical and radio bands. During the IXPE pointings, the VHE state is close to the average behavior with a 0.2-1 TeV flux of 20%-50% the emission of the Crab Nebula. Despite the average VHE activity, an extreme X-ray behavior is measured for the first two IXPE pointings in March 2022 with a synchrotron peak frequency >1 keV. For the third IXPE pointing in July 2022, the synchrotron peak shifts towards lower energies and the optical/X-ray polarization degrees drop. The X-ray polarization is systematically higher than at lower energies, suggesting an energy-stratification of the jet. While during the IXPE epochs the polarization angle in the X-ray, optical and radio bands align well, we find a clear discrepancy in the optical and radio polarization angles in the middle of the campaign. We model the broad-band spectra simultaneous to the IXPE pointings assuming a compact zone dominating in the X-rays and VHE, and an extended zone stretching further downstream the jet dominating the emission at lower energies. NuSTAR data allow us to precisely constrain the synchrotron peak and therefore the underlying electron distribution. The change between the different states observed in the three IXPE pointings can be explained by a change of magnetization and/or emission region size, which directly connects the shift of the synchrotron peak to lower energies with the drop in polarization degree. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: Accepted for publication in Astronomy & Astrophysics. 19 pages, 9 figures. Corresponding authors: Lea Heckmann, Axel Arbet Engels, David Paneque

Detection of X-ray Polarization from the Blazar 1ES 1959+650 with the Imaging X-ray Polarimetry Explorer

Authors: Manel Errando, Ioannis Liodakis, Alan P. Marscher, Herman L. Marshall, Riccardo Middei, Michela Negro, Abel Lawrence Peirson, Matteo Perri, Simonetta Puccetti, Pazit L. Rabinowitz, Iván Agudo, Svetlana G. Jorstad, Sergey S. Savchenko, Dmitry Blinov, Ioakeim G. Bourbah, Sebastian Kiehlmann, Evangelos Kontopodis, Nikos Mandarakas, Stylianos Romanopoulos, Raphael Skalidis, Anna Vervelaki, Francisco José Aceituno, Maria I. Bernardos, Giacomo Bonnoli, Víctor Casanova , et al. (121 additional authors not shown)

Abstract: Observations of linear polarization in the 2-8 keV energy range with the Imaging X-ray Polarimetry Explorer (IXPE) explore the magnetic field geometry and dynamics of the regions generating non-thermal radiation in relativistic jets of blazars. These jets, particularly in blazars whose spectral energy distribution peaks at X-ray energies, emit X-rays via synchrotron radiation from high-energy part… ▽ More Observations of linear polarization in the 2-8 keV energy range with the Imaging X-ray Polarimetry Explorer (IXPE) explore the magnetic field geometry and dynamics of the regions generating non-thermal radiation in relativistic jets of blazars. These jets, particularly in blazars whose spectral energy distribution peaks at X-ray energies, emit X-rays via synchrotron radiation from high-energy particles within the jet. IXPE observations of the X-ray selected BL Lac-type blazar 1ES 1959+650 in 2022 May 3-4 showed a significant linear polarization degree of $Π_\mathrm{x} = 8.0\% \pm 2.3\%$ at an electric-vector position angle $ψ_\mathrm{x} = 123^\circ \pm 8^\circ$. However, in 2022 June 9-12, only an upper limit of $Π_\mathrm{x} \leq 5.1\%$ could be derived (at the 99% confidence level). The degree of optical polarization at that time $Π_\mathrm{O} \sim 5\%$ is comparable to the X-ray measurement. We investigate possible scenarios for these findings, including temporal and geometrical depolarization effects. Unlike some other X-ray selected BL Lac objects, there is no significant chromatic dependence of the measured polarization in 1ES 1959+650, and its low X-ray polarization may be attributed to turbulence in the jet flow with dynamical timescales shorter than 1 day. △ Less

Submitted 9 January, 2024; originally announced January 2024.

Comments: 16 pages, 4 figures, accepted for publication in The Astrophysical Journal

Journal ref: The Astrophysical Journal, 963 (2024) 5

First characterization of the emission behavior of Mrk421 from radio to VHE gamma rays with simultaneous X-ray polarization measurements

Authors: S. Abe, J. Abhir, V. A. Acciari, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, J. Bernete, A. Berti, J. Besenrieder, C. Bigongiari, A. Biland , et al. (229 additional authors not shown)

Abstract: We perform the first broadband study of Mrk421 from radio to TeV gamma rays with simultaneous measurements of the X-ray polarization from IXPE. The data were collected within an extensive multiwavelength campaign organized between May and June 2022 using MAGIC, Fermi-LAT, NuSTAR, XMM-Newton, Swift, and several optical and radio telescopes to complement IXPE. During the IXPE exposures, the measured… ▽ More We perform the first broadband study of Mrk421 from radio to TeV gamma rays with simultaneous measurements of the X-ray polarization from IXPE. The data were collected within an extensive multiwavelength campaign organized between May and June 2022 using MAGIC, Fermi-LAT, NuSTAR, XMM-Newton, Swift, and several optical and radio telescopes to complement IXPE. During the IXPE exposures, the measured 0.2-1 TeV flux is close to the quiescent state and ranges from 25% to 50% of the Crab Nebula without intra-night variability. Throughout the campaign, the VHE and X-ray emission are positively correlated at a $4σ$ significance level. The IXPE measurements unveil a X-ray polarization degree that is a factor of 2-5 higher than in the optical/radio bands; that implies an energy-stratified jet in which the VHE photons are emitted co-spatially with the X-rays, in the vicinity of a shock front. The June 2022 observations exhibit a rotation of the X-ray polarization angle. Despite no simultaneous VHE coverage being available during a large fraction of the swing, the Swift-XRT monitoring unveils an X-ray flux increase with a clear spectral hardening. It suggests that flares in high synchrotron peaked blazars can be accompanied by a polarization angle rotation, as observed in some flat spectrum radio quasars. Finally, during the polarization angle rotation, NuSTAR data reveal two contiguous spectral hysteresis loops in opposite directions (clockwise and counter-clockwise), implying important changes in the particle acceleration efficiency on $\sim$hour timescales. △ Less

Submitted 17 December, 2023; originally announced December 2023.

Comments: Accepted for publication in Astronomy & Astrophysics. 29 pages, 22 figures. Corresponding authors: Axel Arbet Engels, Felix Schmuckermaier, David Paneque

An elementary and direct computation of cohomology with and without a group action

Authors: Makiko Sasada

Abstract: Recently, we introduced a configuration space with interaction structure and a uniform local cohomology on it with co-authors in arXiv:2009.04699. The notion is used to understand a common structure of infinite product spaces appeared in the proof of Varadhan's non-gradient method. For this, the cohomology of the configuration space with a group action is the main target to study, but the cohomolo… ▽ More Recently, we introduced a configuration space with interaction structure and a uniform local cohomology on it with co-authors in arXiv:2009.04699. The notion is used to understand a common structure of infinite product spaces appeared in the proof of Varadhan's non-gradient method. For this, the cohomology of the configuration space with a group action is the main target to study, but the cohomology is easily obtained from that of the configuration space without a group action by applying a well-known property on the group cohomology. In fact, the analysis of the cohomology of the configuration space without a group action is the essential part of arXiv:2009.04699. In this article, we give an elementary and direct proof to obtain the cohomology of a space with a group action from that without a group action under a certain condition including the setting of the configuration space with interaction structure. In particular, no knowledge of group cohomology is required. △ Less

Submitted 4 December, 2023; originally announced December 2023.

Comments: The manuscript was already published (without peer review) in RIMS Kokyuroku No.2177 in 2021

Journal ref: RIMS Kokyuroku No.2177 (2021)

Observations of Low and Intermediate Spectral Peak Blazars with the Imaging X-ray Polarimetry Explorer

Authors: Herman L. Marshall, Ioannis Liodakis, Alan P. Marscher, Niccolo Di Lalla, Svetlana G. Jorstad, Dawoon E. Kim, Riccardo Middei, Michela Negro, Nicola Omodei, Abel L. Peirson, Matteo Perri, Simonetta Puccetti, Ivan Agudo, Giacomo Bonnoli, Andrei V. Berdyugin, Elisabetta Cavazzuti, Nicole Rodriguez Cavero, Immacolata Donnarumma, Laura Di Gesu, Jenni Jormanainen, Henric Krawczynski, Elina Lindfors, Frederic Marin, Francesco Massaro, Luigi Pacciani , et al. (133 additional authors not shown)

Abstract: We present X-ray polarimetry observations from the Imaging X-ray Polarimetry Explorer (IXPE) of three low spectral peak and one intermediate spectral peak blazars, namely 3C 273, 3C 279, 3C 454.3, and S5 0716+714. For none of these objects was IXPE able to detect X-ray polarization at the 3$σ$ level. However, we placed upper limits on the polarization degree at $\sim$10-30\%. The undetected polari… ▽ More We present X-ray polarimetry observations from the Imaging X-ray Polarimetry Explorer (IXPE) of three low spectral peak and one intermediate spectral peak blazars, namely 3C 273, 3C 279, 3C 454.3, and S5 0716+714. For none of these objects was IXPE able to detect X-ray polarization at the 3$σ$ level. However, we placed upper limits on the polarization degree at $\sim$10-30\%. The undetected polarizations favor models where the X-ray band is dominated by unpolarized photons upscattered by relativistic electrons in the jets of blazars, although hadronic models are not completely eliminated. We discuss the X-ray polarization upper limits in the context of our contemporaneous multiwavelength polarization campaigns. △ Less

Submitted 19 October, 2023; v1 submitted 17 October, 2023; originally announced October 2023.

Comments: 19 pages, 10 figures, 6 tables; submitted to the Astrophysical Journal

Magnetic Field Properties inside the Jet of Mrk 421: Multiwavelength Polarimetry Including the Imaging X-ray Polarimetry Explorer

Authors: Dawoon E. Kim, Laura Di Gesu, Ioannis Liodakis, Alan P. Marscher, Svetlana G. Jorstad, Riccardo Midde, Herman L. Marshall, Luigi Pacciani, Iván Agudo, Fabrizio Tavecchio, Nicolò Cibrario, Stefano Tugliani, Raffaella Bonino, Michela Negro, Simonetta Puccetti, Francesco Tombesi, Enrico Costa, Immacolata Donnarumma, Paolo Soffitta, Tsunefumi Mizuno, Yasushi Fukazawa, Koji S. Kawabata, Tatsuya Nakaoka, Makoto Uemura, Ryo Imazawa , et al. (111 additional authors not shown)

Abstract: We conducted a polarimetry campaign from radio to X-ray wavelengths of the high-synchrotron-peak (HSP) blazar Mrk 421, including Imaging X-ray Polarimetry Explorer (IXPE) measurements on 2022 December 6-8. We detected X-ray polarization of Mrk 421 with a degree of $Π_{\rm X}$=14$\pm$1$\%$ and an electric-vector position angle $ψ_{\rm X}$=107$\pm$3$^{\circ}$ in the 2-8 keV band. From the time varia… ▽ More We conducted a polarimetry campaign from radio to X-ray wavelengths of the high-synchrotron-peak (HSP) blazar Mrk 421, including Imaging X-ray Polarimetry Explorer (IXPE) measurements on 2022 December 6-8. We detected X-ray polarization of Mrk 421 with a degree of $Π_{\rm X}$=14$\pm$1$\%$ and an electric-vector position angle $ψ_{\rm X}$=107$\pm$3$^{\circ}$ in the 2-8 keV band. From the time variability analysis, we find a significant episodic variation in $ψ_{\rm X}$. During 7 months from the first IXPE pointing of Mrk 421 in 2022 May, $ψ_{\rm X}$ varied across the range of 0$^{\circ}$ to 180$^{\circ}$, while $Π_{\rm X}$ maintained similar values within $\sim$10-15$\%$. Furthermore, a swing in $ψ_{\rm X}$ in 2022 June was accompanied by simultaneous spectral variations. The results of the multiwavelength polarimetry show that the X-ray polarization degree was generally $\sim$2-3 times greater than that at longer wavelengths, while the polarization angle fluctuated. Additionally, based on radio, infrared, and optical polarimetry, we find that rotation of $ψ$ occurred in the opposite direction with respect to the rotation of $ψ_{\rm X}$ over longer timescales at similar epochs. The polarization behavior observed across multiple wavelengths is consistent with previous IXPE findings for HSP blazars. This result favors the energy-stratified shock model developed to explain variable emission in relativistic jets. The accompanying spectral variation during the $ψ_{\rm X}$ rotation can be explained by a fluctuation in the physical conditions, e.g., in the energy distribution of relativistic electrons. The opposite rotation direction of $ψ$ between the X-ray and longer-wavelength polarization accentuates the conclusion that the X-ray emitting region is spatially separated from that at longer wavelengths. △ Less

Submitted 9 October, 2023; originally announced October 2023.

Comments: 17 pages, 13 figures, 4 tables; Accepted for publication in A&A

Journal ref: A&A 681, A12 (2024)

Gamma rays from a reverse shock with turbulent magnetic fields in GRB 180720B

Authors: Makoto Arimoto, Katsuaki Asano, Koji S. Kawabata, Kenji Toma, Ramandeep Gill, Jonathan Granot, Masanori Ohno, Shuta Takahashi, Naoki Ogino, Hatsune Goto, Kengo Nakamura, Tatsuya Nakaoka, Kengo Takagi, Miho Kawabata, Masayuki Yamanaka, Mahito Sasada, Soebur Razzaque

Abstract: Gamma-ray bursts (GRBs) are the most electromagnetically luminous cosmic explosions. They are powered by collimated streams of plasma (jets) ejected by a newborn stellar-mass black hole or neutron star at relativistic velocities (near the speed of light). Their short-lived (typically tens of seconds) prompt $γ$-ray emission from within the ejecta is followed by long-lived multi-wavelength afterglo… ▽ More Gamma-ray bursts (GRBs) are the most electromagnetically luminous cosmic explosions. They are powered by collimated streams of plasma (jets) ejected by a newborn stellar-mass black hole or neutron star at relativistic velocities (near the speed of light). Their short-lived (typically tens of seconds) prompt $γ$-ray emission from within the ejecta is followed by long-lived multi-wavelength afterglow emission from the ultra-relativistic forward shock. This shock is driven into the circumburst medium by the GRB ejecta that are in turn decelerated by a mildly-relativistic reverse shock. Forward shock emission was recently detected up to teraelectronvolt-energy $γ$-rays, and such very-high-energy emission was also predicted from the reverse shock. Here we report the detection of optical and gigaelectronvolt-energy $γ$-ray emission from GRB 180720B during the first few hundred seconds, which is explained by synchrotron and inverse-Compton emission from the reverse shock propagating into the ejecta, implying a low-magnetization ejecta. Our optical measurements show a clear transition from the reverse shock to the forward shock driven into the circumburst medium, accompanied by a 90-degree change in the mean polarization angle and fluctuations in the polarization degree and angle. This indicates turbulence with large-scale toroidal and radially-stretched magnetic field structures in the reverse and forward shocks, respectively, which tightly couple to the physics of relativistic shocks and GRB jets -- launching, composition, dissipation and particle acceleration. △ Less

Submitted 6 October, 2023; originally announced October 2023.

Comments: 5 pages, 4 figures (main) plus Methods and Supplementary Methods, accepted for publication

A search for pulsars around Sgr A* in the first Event Horizon Telescope dataset

Authors: Pablo Torne, Kuo Liu, Ralph P. Eatough, Jompoj Wongphechauxsorn, James M. Cordes, Gregory Desvignes, Mariafelicia De Laurentis, Michael Kramer, Scott M. Ransom, Shami Chatterjee, Robert Wharton, Ramesh Karuppusamy, Lindy Blackburn, Michael Janssen, Chi-kwan Chan, Geoffrey B. Crew, Lynn D. Matthews, Ciriaco Goddi, Helge Rottmann, Jan Wagner, Salvador Sanchez, Ignacio Ruiz, Federico Abbate, Geoffrey C. Bower, Juan J. Salamanca , et al. (261 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission… ▽ More The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission spectra - are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic Center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most-sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the Fast-Folding-Algorithm and single pulse search targeting both pulsars and burst-like transient emission; using the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction ($\lesssim$2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

X-ray Polarization of BL Lacertae in Outburst

Authors: Abel L. Peirson, Michela Negro, Ioannis Liodakis, Riccardo Middei, Dawoon E. Kim, Alan P. Marscher, Herman L. Marshall, Luigi Pacciani, Roger W. Romani, Kinwah Wu, Alessandro Di Marco, Niccolo Di Lalla, Nicola Omodei, Svetlana G. Jorstad, Ivan Agudo, Pouya M. Kouch, Elina Lindfors, Francisco Jose Aceituno, Maria I. Bernardos, Giacomo Bonnoli, Victor Casanova, Maya Garcia-Comas, Beatriz Agis-Gonzalez, Cesar Husillos, Alessandro Marchini , et al. (106 additional authors not shown)

Abstract: We report the first $> 99\%$ confidence detection of X-ray polarization in BL Lacertae. During a recent X-ray/$γ$-ray outburst, a 287 ksec observation (2022 November 27-30) was taken using the Imaging X-ray Polarimetry Explorer ({\it IXPE}), together with contemporaneous multiwavelength observations from the Neil Gehrels {\it Swift} observatory and {\it XMM-Newton} in soft X-rays (0.3--10~keV), {\… ▽ More We report the first $> 99\%$ confidence detection of X-ray polarization in BL Lacertae. During a recent X-ray/$γ$-ray outburst, a 287 ksec observation (2022 November 27-30) was taken using the Imaging X-ray Polarimetry Explorer ({\it IXPE}), together with contemporaneous multiwavelength observations from the Neil Gehrels {\it Swift} observatory and {\it XMM-Newton} in soft X-rays (0.3--10~keV), {\it NuSTAR} in hard X-rays (3--70~keV), and optical polarization from the Calar Alto, and Perkins Telescope observatories. Our contemporaneous X-ray data suggest that the {\it IXPE} energy band is at the crossover between the low- and high-frequency blazar emission humps. The source displays significant variability during the observation, and we measure polarization in three separate time bins. Contemporaneous X-ray spectra allow us to determine the relative contribution from each emission hump. We find $>99\%$ confidence X-ray polarization $Π_{2-4{\rm keV}} = 21.7^{+5.6}_{-7.9}\%$ and electric vector polarization angle $ψ_{2-4{\rm keV}} = -28.7 \pm 8.7^{\circ}$ in the time bin with highest estimated synchrotron flux contribution. We discuss possible implications of our observations, including previous {\it IXPE} BL Lacertae pointings, tentatively concluding that synchrotron self-Compton emission dominates over hadronic emission processes during the observed epochs. △ Less

Submitted 3 July, 2023; v1 submitted 23 May, 2023; originally announced May 2023.

Comments: Accepted to ApJL. 16 pages, 8 figures

Discovery of X-ray polarization angle rotation in active galaxy Mrk 421

Authors: Laura Di Gesu, Herman L. Marshall, Steven R. Ehlert, Dawoon E. Kim, Immacolata Donnarumma, Fabrizio Tavecchio, Ioannis Liodakis, Sebastian Kiehlmann, Iván Agudo, Svetlana G. Jorstad, Fabio Muleri, Alan P. Marscher, Simonetta Puccetti, Riccardo Middei, Matteo Perri, Luigi Pacciani, Michela Negro, Roger W. Romani, Alessandro Di Marco, Dmitry Blinov, Ioakeim G. Bourbah, Evangelos Kontopodis, Nikos Mandarakas, Stylianos Romanopoulos, Raphael Skalidis , et al. (118 additional authors not shown)

Abstract: The magnetic field conditions in astrophysical relativistic jets can be probed by multiwavelength polarimetry, which has been recently extended to X-rays. For example, one can track how the magnetic field changes in the flow of the radiating particles by observing rotations of the electric vector position angle $Ψ$. Here we report the discovery of a $Ψ_{\mathrm x}$ rotation in the X-ray band in th… ▽ More The magnetic field conditions in astrophysical relativistic jets can be probed by multiwavelength polarimetry, which has been recently extended to X-rays. For example, one can track how the magnetic field changes in the flow of the radiating particles by observing rotations of the electric vector position angle $Ψ$. Here we report the discovery of a $Ψ_{\mathrm x}$ rotation in the X-ray band in the blazar Mrk 421 at an average flux state. Across the 5 days of Imaging X-ray Polarimetry Explorer (IXPE) observations of 4-6 and 7-9 June 2022, $Ψ_{\mathrm x}$ rotated in total by $\geq360^\circ$. Over the two respective date ranges, we find constant, within uncertainties, rotation rates ($80 \pm 9$ and $91 \pm 8 ^\circ/\rm day$) and polarization degrees ($Π_{\mathrm x}=10\%\pm1\%$). Simulations of a random walk of the polarization vector indicate that it is unlikely that such rotation(s) are produced by a stochastic process. The X-ray emitting site does not completely overlap the radio/infrared/optical emission sites, as no similar rotation of $Ψ$ was observed in quasi-simultaneous data at longer wavelengths. We propose that the observed rotation was caused by a helical magnetic structure in the jet, illuminated in the X-rays by a localized shock propagating along this helix. The optically emitting region likely lies in a sheath surrounding an inner spine where the X-ray radiation is released. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comparison of Polarized Radiative Transfer Codes used by the EHT Collaboration

Authors: Ben S. Prather, Jason Dexter, Monika Moscibrodzka, Hung-Yi Pu, Thomas Bronzwaer, Jordy Davelaar, Ziri Younsi, Charles F. Gammie, Roman Gold, George N. Wong, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, Anne-Kathrin Baczko, David Ball, Mislav Baloković, John Barrett, Michi Bauböck, Bradford A. Benson, Dan Bintley , et al. (248 additional authors not shown)

Abstract: Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curve… ▽ More Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U , and V respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ

Follow-up survey for the binary black hole merger GW200224_222234 using Subaru/HSC and GTC/OSIRIS

Authors: Takayuki Ohgami, Josefa Becerra Gonzalez, Nozomu Tominaga, Tomoki Morokuma, Yousuke Utsumi, Yuu Niino, Masaomi Tanaka, Smaranika Banerjee, Frederick Poidevin, Jose Antonio Acosta-Pulido, Ismael Perez-Fournon, Teo Munoz-Darias, Hiroshi Akitaya, Kenshi Yanagisawa, Mahito Sasada, Michitoshi Yoshida, Mirko Simunovic, Ryou Ohsawa, Ichi Tanaka, Tsuyoshi Terai, Yuhei Takagi, The J-GEM collaboration

Abstract: The LIGO/Virgo detected a gravitational wave (GW) event, named GW200224_222234 (a.k.a. S200224ca) and classified as a binary-black-hole coalescence, on February 24, 2020. Given its relatively small localization skymap (71 deg$^2$ for a 90% credible region; revised to 50 deg$^2$ in GWTC-3), we performed target-of-opportunity observations using the Subaru/Hyper Suprime-Cam (HSC) in the $r2$- and… ▽ More The LIGO/Virgo detected a gravitational wave (GW) event, named GW200224_222234 (a.k.a. S200224ca) and classified as a binary-black-hole coalescence, on February 24, 2020. Given its relatively small localization skymap (71 deg$^2$ for a 90% credible region; revised to 50 deg$^2$ in GWTC-3), we performed target-of-opportunity observations using the Subaru/Hyper Suprime-Cam (HSC) in the $r2$- and $z$-bands. Observations were conducted on February 25 and 28 and March 23, 2020, with the first epoch beginning 12.3 h after the GW detection. The survey covered the highest probability sky area of 56.6 deg$^2$, corresponding to a 91% probability. This was the first deep follow-up ($m_{r}\gtrsim24, m_{z}\gtrsim23$) for a binary-black-hole merger covering $>$90% of the localization. By performing image subtraction and candidate screening including light curve fitting with transient templates and examples, we found 22 off-nucleus transients that were not ruled out as the counterparts of GW200224_222234 with only our Subaru/HSC data. We also performed GTC/OSIRIS spectroscopy of the probable host galaxies for five candidates; two are likely to be located within the 3D skymap, whereas the others are not. In conclusion, 19 transients remain as possible optical counterparts of GW200224_222234; however, we could not identify a unique promising counterpart. If there are no counterparts in the remaining candidates, the upper limits of optical luminosity are $νL_ν < 5.2^{+2.4}_{-1.9}\times 10^{41}$ erg s$^{-1}$ and $νL_ν < 1.8^{+0.8}_{-0.6}\times 10^{42}$ erg s$^{-1}$ in the $r2$- and $z$-bands, respectively, at $\sim$12 h after GW detection. We also discuss improvements in the strategies of optical follow-ups for future GW events. △ Less

Submitted 18 February, 2023; originally announced February 2023.

Comments: 22 pages, 15 figures. Accepted for publication in the Astrophysical Journal

Relationships between two linearizations of the box-ball system : Kerov-Kirillov-Reschetikhin bijection and slot configuration

Authors: Matteo Mucciconi, Makiko Sasada, Tomohiro Sasamoto, Hayate Suda

Abstract: The box-ball system (BBS), which was introduced by Takahashi and Satsuma in 1990, is a soliton cellular automaton. Its dynamics can be linearized by a few methods, among which the best known is the Kerov-Kirillov-Reschetikhin (KKR) bijection using rigged partitions. Recently a new linearization method in terms of "slot configurations" was introduced by Ferrari-Nguyen-Rolla-Wang, but its relations… ▽ More The box-ball system (BBS), which was introduced by Takahashi and Satsuma in 1990, is a soliton cellular automaton. Its dynamics can be linearized by a few methods, among which the best known is the Kerov-Kirillov-Reschetikhin (KKR) bijection using rigged partitions. Recently a new linearization method in terms of "slot configurations" was introduced by Ferrari-Nguyen-Rolla-Wang, but its relations to existing ones have not been clarified. In this paper we investigate this issue and clarify the relation between the two linearizations. For this we introduce a novel way of describing the BBS dynamics using a carrier with seat numbers. We show that the seat number configuration also linearizes the BBS and reveals explicit relations between the KKR bijection and the slot configuration. In addition, by using these explicit relations, we also show that even in case of finite carrier capacity the BBS can be linearized via the slot configuration. △ Less

Submitted 31 December, 2022; originally announced January 2023.

Comments: 33 pages, 12 figures

MSC Class: 37B15; 82C22; 82C23

Yang-Baxter maps and independence preserving property

Authors: Makiko Sasada, Ryosuke Uozumi

Abstract: We study a surprising relationship between two properties for bijective functions $F : \mathcal{X} \times \mathcal{X} \to \mathcal{X} \times \mathcal{X}$ for a set $\mathcal{X}$ which are introduced from very different backgrounds. One of the property is that $F$ is a Yang-Baxter map, namely it satisfies the "set-theoretical" Yang-Baxter equation, and the other property is the independence preserv… ▽ More We study a surprising relationship between two properties for bijective functions $F : \mathcal{X} \times \mathcal{X} \to \mathcal{X} \times \mathcal{X}$ for a set $\mathcal{X}$ which are introduced from very different backgrounds. One of the property is that $F$ is a Yang-Baxter map, namely it satisfies the "set-theoretical" Yang-Baxter equation, and the other property is the independence preserving property (IP property for short), which means that there exist independent (non-constant) $\mathcal{X}$-valued random variables $X,Y$ such that $U,V$ are also independent with $(U,V)=F(X,Y)$. Recently in the study of invariant measures for a discrete integrable system, a class of functions having these two properties were found. Motivated by this, we analyze a relationship between the Yang-Baxter maps and the IP property, which has never been studied as far as we are aware, focusing on the case $\mathcal{X}=\mathbb{R}_+$. Our first main result is that all quadrirational Yang-Baxter maps $F : \mathbb{R}_+ \times \mathbb{R}_+ \to \mathbb{R}_+ \times \mathbb{R}_+$ in the most interesting subclass have the independence preserving property. In particular, we find new classes of bijections having the IP property. Our second main result is that these newly introduce bijections are fundamental in the class of (known) bijections with the IP property, in the sense that most of known bijections having the IP property are derived from these maps by taking special parameters or performing some limiting procedure. This reveals that the IP property, which has been investigated for specific functions individually, can be understood in a unified manner. △ Less

Submitted 14 December, 2022; v1 submitted 1 December, 2022; originally announced December 2022.

Comments: 19 pages

X-ray Polarization Observations of BL Lacertae

Authors: Riccardo Middei, Ioannis Liodakis, Matteo Perri, Simonetta Puccetti, Elisabetta Cavazzuti, Laura Di Gesu, Steven R. Ehlert, Grzegorz Madejski, Alan P. Marscher, Herman L. Marshall, Fabio Muleri, Michela Negro, Svetlana G. Jorstad, Beatriz Agís-González, Iván Agudo, Giacomo Bonnoli, Maria I. Bernardos, Víctor Casanova, Maya García-Comas, César Husillos, Alessandro Marchini, Alfredo Sota, Pouya M. Kouch, George A. Borman, Evgenia N. Kopatskaya , et al. (121 additional authors not shown)

Abstract: Blazars are a class of jet-dominated active galactic nuclei with a typical double-humped spectral energy distribution. It is of common consensus the Synchrotron emission to be responsible for the low frequency peak, while the origin of the high frequency hump is still debated. The analysis of X-rays and their polarization can provide a valuable tool to understand the physical mechanisms responsibl… ▽ More Blazars are a class of jet-dominated active galactic nuclei with a typical double-humped spectral energy distribution. It is of common consensus the Synchrotron emission to be responsible for the low frequency peak, while the origin of the high frequency hump is still debated. The analysis of X-rays and their polarization can provide a valuable tool to understand the physical mechanisms responsible for the origin of high-energy emission of blazars. We report the first observations of BL Lacertae performed with the Imaging X-ray Polarimetry Explorer ({IXPE}), from which an upper limit to the polarization degree $Π_X<$12.6\% was found in the 2-8 keV band. We contemporaneously measured the polarization in radio, infrared, and optical wavelengths. Our multiwavelength polarization analysis disfavors a significant contribution of proton synchrotron radiation to the X-ray emission at these epochs. Instead, it supports a leptonic origin for the X-ray emission in BL Lac. △ Less

Submitted 24 November, 2022; originally announced November 2022.

Comments: 17 pages, 5 figures, accepted for publication in ApJL

Journal ref: 2023, ApJ Letters, 942, L10

The microvariability and wavelength dependence of polarization degree/angle of BL Lacertae in the outburst 2020 to 2021

Authors: Ryo Imazawa, Mahito Sasada, Natsuko Hazama, Yasushi Fukazawa, Koji S. Kawabata, Tatsuya Nakaoka, Hiroshi Akitaya, Thomas Bohn, Anjasha Gangopadhyay

Abstract: We have obtained simultaneous and continuous photo-polarization observations of the blazar BL Lacertae in optical and near-infrared (NIR) bands during a historical outburst from 2020 to 2021. In total, fourteen nights of observations were performed where ten observations show microvariability on timescales of a few minutes to several hours. This suggests a compact emission region, and the timescal… ▽ More We have obtained simultaneous and continuous photo-polarization observations of the blazar BL Lacertae in optical and near-infrared (NIR) bands during a historical outburst from 2020 to 2021. In total, fourteen nights of observations were performed where ten observations show microvariability on timescales of a few minutes to several hours. This suggests a compact emission region, and the timescales are difficult to explain by a one-zone shock-in-jet model. Moreover, we found significant differences in the polarization degree (PD) and angle between optical and NIR bands. Nine nights showed a PD in the optical band that is greater than or equal to that in the NIR band, which can be explained by either a shock-in-jet model or the Turbulent Extreme Multi-Zone (TEMZ) model. On the other hand, five nights showed higher PD in a NIR band than an optical band, which cannot be explained by simple shock-in-jet models nor the simple TEMZ model. The observed timescales and wavelength dependency of the PD and polarization angle suggest the existence of complicated multiple emission regions such as an irregular TEMZ model. △ Less

Submitted 21 November, 2022; v1 submitted 17 November, 2022; originally announced November 2022.

Comments: 19 pages, 8 figures

Journal ref: published, Publications of the Astronomical Society of Japan, 2022, psac084

Varadhan's decomposition of shift-invariant closed uniform forms for large scale interacting systems on general crystal lattices

Authors: Kenichi Bannai, Makiko Sasada

Abstract: We prove a uniform version of Varadhan decomposition for shift-invariant closed uniform forms associated to large scale interacting systems on general crystal lattices. In particular, this result includes the case of translation invariant processes on Euclidean lattices $\mathbf{Z}^d$ with finite range. Our result generalizes the result of arXiv:2009.04699 which was valid for systems on transferab… ▽ More We prove a uniform version of Varadhan decomposition for shift-invariant closed uniform forms associated to large scale interacting systems on general crystal lattices. In particular, this result includes the case of translation invariant processes on Euclidean lattices $\mathbf{Z}^d$ with finite range. Our result generalizes the result of arXiv:2009.04699 which was valid for systems on transferable graphs. In subsequent research, we will use the result of this article to prove Varadhan's decomposition of closed $L^2$-forms for large scale interacting systems on general crystal lattices. △ Less

Submitted 31 August, 2022; originally announced August 2022.

Comments: 28 pages, 2 figures,

MSC Class: Primary: 82C22; Secondary: 05C63; 55N91; 60J60; 70G40

Varadhan's Decomposition of Shift-Invariant Closed $L^2$-forms for Large Scale Interacting Systems on the Euclidean Lattice

Authors: Kenichi Bannai, Makiko Sasada

Abstract: We rigorously formulate and prove for a relatively general class of interactions Varadhan's Decomposition of shift-invariant closed $L^2$-forms for a large scale interacting system on a Euclidean lattice $(\mathbb{Z}^d,\mathbb{E}^d)$. Such decomposition of closed forms has played an essential role in proving the diffusive scaling limit of nongradient systems. The universal expression in terms of c… ▽ More We rigorously formulate and prove for a relatively general class of interactions Varadhan's Decomposition of shift-invariant closed $L^2$-forms for a large scale interacting system on a Euclidean lattice $(\mathbb{Z}^d,\mathbb{E}^d)$. Such decomposition of closed forms has played an essential role in proving the diffusive scaling limit of nongradient systems. The universal expression in terms of conserved quantities was sought from observations for specific models, but a precise formulation or rigorous proof up until now had been elusive. Our result is based on the universal decomposition of shift-invariant closed uniform forms studied in our previous article (arXiv:2009.04699). In the present article, we show that the same universal structure also appears for $L^2$-forms. The essential assumptions are: (i) the set of states on each vertex is a finite set, (ii) the measure on the configuration space is the product measure, and (iii) there is a certain uniform spectral gap estimate for the mean field version of the interaction. Our result is applicable for generalized exclusion processes, multi-species exclusion processes, and more general lattice gas models. △ Less

Submitted 18 March, 2022; v1 submitted 17 November, 2021; originally announced November 2021.

MSC Class: Primary: 82C22; Secondary: 05C63; 55N91; 60J60; 70G40

Event Horizon Telescope observations of the jet launching and collimation in Centaurus A

Authors: Michael Janssen, Heino Falcke, Matthias Kadler, Eduardo Ros, Maciek Wielgus, Kazunori Akiyama, Mislav Baloković, Lindy Blackburn, Katherine L. Bouman, Andrew Chael, Chi-kwan Chan, Koushik Chatterjee, Jordy Davelaar, Philip G. Edwards, Christian M. Fromm, José L. Gómez, Ciriaco Goddi, Sara Issaoun, Michael D. Johnson, Junhan Kim, Jun Yi Koay, Thomas P. Krichbaum, Jun Liu, Elisabetta Liuzzo, Sera Markoff , et al. (215 additional authors not shown)

Abstract: Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supe… ▽ More Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supermassive black holes (SMBHs) in Messier 87 and our galactic center. A large southern declination of $-43^{\circ}$ has however prevented VLBI imaging of Centaurus A below $λ1$cm thus far. Here, we show the millimeter VLBI image of the source, which we obtained with the Event Horizon Telescope at $228$GHz. Compared to previous observations, we image Centaurus A's jet at a tenfold higher frequency and sixteen times sharper resolution and thereby probe sub-lightday structures. We reveal a highly-collimated, asymmetrically edge-brightened jet as well as the fainter counterjet. We find that Centaurus A's source structure resembles the jet in Messier 87 on ${\sim}500r_g$ scales remarkably well. Furthermore, we identify the location of Centaurus A's SMBH with respect to its resolved jet core at $λ1.3$mm and conclude that the source's event horizon shadow should be visible at THz frequencies. This location further supports the universal scale invariance of black holes over a wide range of masses. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 27 pages, 9 figures. This is a post-peer-review, pre-copyedit version of an article published in Nature Astronomy. The final authenticated version is available online at: https://doi.org/10.1038/s41550-021-01417-w

Journal ref: Nature Astronomy, July 2021, Volume 5, p. 1017-1028

The Variability of the Black-Hole Image in M87 at the Dynamical Time Scale

Authors: Kaushik Satapathy, Dimitrios Psaltis, Feryal Ozel, Lia Medeiros, Sean T. Dougall, Chi-kwan Chan, Maciek Wielgus, Ben S. Prather, George N. Wong, Charles F. Gammie, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David R. Ball, Mislav Baloković, John Barrett, Bradford A. Benson, Dan Bintley, Lindy Blackburn, Raymond Blundell , et al. (213 additional authors not shown)

Abstract: The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expect… ▽ More The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure phase measurements on all six linearly independent non-trivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of $\sim3-5^\circ$. The only triangles that exhibit substantially higher variability ($\sim90-180^\circ$) are the ones with baselines that cross visibility amplitude minima on the $u-v$ plane, as expected from theoretical modeling. We used two sets of General Relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black-hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black-hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: Accepted for Publication in ApJ

J-GEM optical and near-infrared follow-up of gravitational wave events during LIGO's and Virgo's third observing run

Authors: Mahito Sasada, Yousuke Utsumi, Ryosuke Itoh, Nozomu Tominaga, Masaomi Tanaka, Tomoki Morokuma, Kenshi Yanagisawa, Koji S. Kawabata, Takayuki Ohgami, Michitoshi Yoshida, Fumio Abe, Ryo Adachi, Hiroshi Akitaya, Yang Chong, Kazuki Daikuhara, Ryo Hamasaki, Satoshi Honda, Ryohei Hosokawa, Kota Iida, Fumiya Imazato, Chihiro Ishioka, Takumi Iwasaki, Mingjie Jian, Yuhei Kamei, Takahiro Kanai , et al. (46 additional authors not shown)

Abstract: The Laser Interferometer Gravitational-wave Observatory Scientific Collaboration and Virgo Collaboration (LVC) sent out 56 gravitational-wave (GW) notices during the third observing run (O3). Japanese collaboration for Gravitational wave ElectroMagnetic follow-up (J-GEM) performed optical and near-infrared observations to identify and observe an electromagnetic (EM) counterpart. We constructed web… ▽ More The Laser Interferometer Gravitational-wave Observatory Scientific Collaboration and Virgo Collaboration (LVC) sent out 56 gravitational-wave (GW) notices during the third observing run (O3). Japanese collaboration for Gravitational wave ElectroMagnetic follow-up (J-GEM) performed optical and near-infrared observations to identify and observe an electromagnetic (EM) counterpart. We constructed web-based system which enabled us to obtain and share information of candidate host galaxies for the counterpart, and status of our observations. Candidate host galaxies were selected from the GLADE catalog with a weight based on the three-dimensional GW localization map provided by LVC. We conducted galaxy-targeted and wide-field blind surveys, real-time data analysis, and visual inspection of observed galaxies. We performed galaxy-targeted follow-ups to 23 GW events during O3, and the maximum probability covered by our observations reached to 9.8%. Among them, we successfully started observations for 10 GW events within 0.5 days after the detection. This result demonstrates that our follow-up observation has a potential to constrain EM radiation models for a merger of binary neutron stars at a distance of up to $\sim$100~Mpc with a probability area of $\leq$ 500~deg$^2$. △ Less

Submitted 9 June, 2021; originally announced June 2021.

Comments: 25 pages, 10 figures, 3 tables, published at https://academic.oup.com/ptep/article/2021/5/05A104/6122445

Constraints on black-hole charges with the 2017 EHT observations of M87*

Authors: Prashant Kocherlakota, Luciano Rezzolla, Heino Falcke, Christian M. Fromm, Michael Kramer, Yosuke Mizuno, Antonios Nathanail, Hector Olivares, Ziri Younsi, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David Ball, Mislav Balokovic, John Barrett, Bradford A. Benson, Dan Bintley, Lindy Blackburn, Raymond Blundell, Wilfred Boland , et al. (212 additional authors not shown)

Abstract: Our understanding of strong gravity near supermassive compact objects has recently improved thanks to the measurements made by the Event Horizon Telescope (EHT). We use here the M87* shadow size to infer constraints on the physical charges of a large variety of nonrotating or rotating black holes. For example, we show that the quality of the measurements is already sufficient to rule out that M87*… ▽ More Our understanding of strong gravity near supermassive compact objects has recently improved thanks to the measurements made by the Event Horizon Telescope (EHT). We use here the M87* shadow size to infer constraints on the physical charges of a large variety of nonrotating or rotating black holes. For example, we show that the quality of the measurements is already sufficient to rule out that M87* is a highly charged dilaton black hole. Similarly, when considering black holes with two physical and independent charges, we are able to exclude considerable regions of the space of parameters for the doubly-charged dilaton and the Sen black holes. △ Less

Submitted 19 May, 2021; originally announced May 2021.

Comments: 15 pages, 4 figures, published in PRD on May 19

A Decomposition Theorem of Varadhan Type for Co-local Forms on Large Scale Interacting Systems

Authors: Kenichi Bannai, Makiko Sasada

Abstract: In our previous article with Yukio Kametani, we investigated the geometric structure underlying a large scale interacting system on infinite graphs, via constructing a suitable cohomology theory called uniformly local cohomology, which reflects the geometric property of the microscopic model, using a class of functions called the uniformly local functions. In this article, we introduce the co-loca… ▽ More In our previous article with Yukio Kametani, we investigated the geometric structure underlying a large scale interacting system on infinite graphs, via constructing a suitable cohomology theory called uniformly local cohomology, which reflects the geometric property of the microscopic model, using a class of functions called the uniformly local functions. In this article, we introduce the co-local functions on the geometric structure associated to a large scale interacting system. We may similarly define the notion of uniformly local functions for co-local functions. However, contrary to the functions appearing in our previous article, the co-local functions reflect the stochastic property of the model, namely the probability measure on the configuration space. We then prove a decomposition theorem of Varadhan type for closed co-local forms. The space of co-local functions and forms contain the space of $L^2$-functions and forms. In the last section, we state a conjecture concerning the decomposition theorem for the $L^2$-case. △ Less

Submitted 18 March, 2022; v1 submitted 12 May, 2021; originally announced May 2021.

Comments: 23 pages. Included reference to subsequent research

MSC Class: Primary: 82C22; Secondary: 05C63; 55N91; 60J60; 70G40

The Polarized Image of a Synchrotron Emitting Ring of Gas Orbiting a Black Hole

Authors: Ramesh Narayan, Daniel C. M. Palumbo, Michael D. Johnson, Zachary Gelles, Elizabeth Himwich, Dominic O. Chang, Angelo Ricarte, Jason Dexter, Charles F. Gammie, Andrew A. Chael, The Event Horizon Telescope Collaboration, :, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David Ball, Mislav Balokovic, John Barrett, Bradford A. Benson, Dan Bintley , et al. (215 additional authors not shown)

Abstract: Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equ… ▽ More Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov (2002) and conservation of the Walker-Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images. △ Less

Submitted 13 May, 2021; v1 submitted 4 May, 2021; originally announced May 2021.

Comments: 29 pages, 11 figures, published in ApJ on May 3

Journal ref: ApJ 912 35 (2021)

Broadband Multi-wavelength Properties of M87 during the 2017 Event Horizon Telescope Campaign

Authors: J. C. Algaba, J. Anczarski, K. Asada, M. Balokovic, S. Chandra, Y. -Z. Cui, A. D. Falcone, M. Giroletti, C. Goddi, K. Hada, D. Haggard, S. Jorstad, A. Kaur, T. Kawashima, G. Keating, J. -Y. Kim, M. Kino, S. Komossa, E. V. Kravchenko, T. P. Krichbaum, S. -S. Lee, R. -S. Lu, M. Lucchini, S. Markoff, J. Neilsen , et al. (14 additional authors not shown)

Abstract: In 2017, the Event Horizon Telescope (EHT) Collaboration succeeded in capturing the first direct image of the center of the M87 galaxy. The asymmetric ring morphology and size are consistent with theoretical expectations for a weakly accreting supermassive black hole of mass approximately 6.5 x 10^9 M_solar. The EHTC also partnered with several international facilities in space and on the ground,… ▽ More In 2017, the Event Horizon Telescope (EHT) Collaboration succeeded in capturing the first direct image of the center of the M87 galaxy. The asymmetric ring morphology and size are consistent with theoretical expectations for a weakly accreting supermassive black hole of mass approximately 6.5 x 10^9 M_solar. The EHTC also partnered with several international facilities in space and on the ground, to arrange an extensive, quasi-simultaneous multi-wavelength campaign. This Letter presents the results and analysis of this campaign, as well as the multi-wavelength data as a legacy data repository. We captured M87 in a historically low state, and the core flux dominates over HST-1 at high energies, making it possible to combine core flux constraints with the more spatially precise very long baseline interferometry data. We present the most complete simultaneous multi-wavelength spectrum of the active nucleus to date, and discuss the complexity and caveats of combining data from different spatial scales into one broadband spectrum. We apply two heuristic, isotropic leptonic single-zone models to provide insight into the basic source properties, but conclude that a structured jet is necessary to explain M87's spectrum. We can exclude that the simultaneous gamma-ray emission is produced via inverse Compton emission in the same region producing the EHT mm-band emission, and further conclude that the gamma-rays can only be produced in the inner jets (inward of HST-1) if there are strongly particle-dominated regions. Direct synchrotron emission from accelerated protons and secondaries cannot yet be excluded. △ Less

Submitted 14 April, 2021; originally announced April 2021.

Comments: 43 pages, 19 figures, 10 tables; complete author list available in manuscript; The Astrophysical Journal Letters, 2021, 911, L11; publication doi: 3847/2041-8213/abef71, data doi: 10.25739/mhh2-cw46

On the stationary solutions of random polymer models and their zero-temperature limits

Authors: David A. Croydon, Makiko Sasada

Abstract: We derive stationary measures for certain zero-temperature random polymer models, which we believe are new in the case of the zero-temperature limit of the beta random polymer (that has been called the river delta model). To do this, we apply techniques developed for understanding the stationary measures of the corresponding positive-temperature random polymer models (and some deterministic integr… ▽ More We derive stationary measures for certain zero-temperature random polymer models, which we believe are new in the case of the zero-temperature limit of the beta random polymer (that has been called the river delta model). To do this, we apply techniques developed for understanding the stationary measures of the corresponding positive-temperature random polymer models (and some deterministic integrable systems). More specifically, the article starts with a survey of results for the four basic beta-gamma models (i.e. the inverse-gamma, gamma, inverse-beta and beta random polymers), highlighting how the maps underlying the systems in question can each be reduced to one of two basic bijections, and that through an `independence preservation' property, these bijections characterise the associated stationary measures. We then derive similar descriptions for the corresponding zero-temperature maps, whereby each is written in terms of one of two bijections. One issue with this picture is that, unlike in the positive-temperature case, the change of variables required is degenerate in general, and so whilst the argument yields stationary solutions, it does not provide a complete characterisation of them. On the other hand, this degeneracy does allow us to explain the appearance of atoms in the stationary measures of certain zero-temperature models. We also derive from our viewpoint various links between random polymer models, some of which recover known results, some of which are novel, and some of which lead to further questions. △ Less

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

MSC Class: 82D60 (primary); 37K60; 37L40; 60E05; 60K35; 82B23

Gravitational Wave Physics and Astronomy in the nascent era

Authors: Makoto Arimoto, Hideki Asada, Michael L. Cherry, Michiko S. Fujii, Yasushi Fukazawa, Akira Harada, Kazuhiro Hayama, Takashi Hosokawa, Kunihito Ioka, Yoichi Itoh, Nobuyuki Kanda, Koji S. Kawabata, Kyohei Kawaguchi, Nobuyuki Kawai, Tsutomu Kobayashi, Kazunori Kohri, Yusuke Koshio, Kei Kotake, Jun Kumamoto, Masahiro N. Machida, Hideo Matsufuru, Tatehiro Mihara, Masaki Mori, Tomoki Morokuma, Shinji Mukohyama , et al. (28 additional authors not shown)

Abstract: The detections of gravitational waves (GW) by LIGO/Virgo collaborations provide various possibilities to physics and astronomy. We are quite sure that GW observations will develop a lot both in precision and in number owing to the continuous works for the improvement of detectors, including the expectation to the newly joined detector, KAGRA, and the planned detector, LIGO-India. In this occasion,… ▽ More The detections of gravitational waves (GW) by LIGO/Virgo collaborations provide various possibilities to physics and astronomy. We are quite sure that GW observations will develop a lot both in precision and in number owing to the continuous works for the improvement of detectors, including the expectation to the newly joined detector, KAGRA, and the planned detector, LIGO-India. In this occasion, we review the fundamental outcomes and prospects of gravitational wave physics and astronomy. We survey the development focusing on representative sources of gravitational waves: binary black holes, binary neutron stars, and supernovae. We also summarize the role of gravitational wave observations as a probe of new physics. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Comments: 89pages, Prog. Theor. Exp. Phys. accepted

ASASSN-18aan: An Eclipsing SU UMa-type Cataclysmic Variable with a 3.6-hour Orbital Period and a Late G-type Secondary Star

Authors: Yasuyuki Wakamatsu, John R. Thorstensen, Naoto Kojiguchi, Keisuke Isogai, Mariko Kimura, Ryuhei Ohnishi, Taichi Kato, Hiroshi Itoh, Yuki Sugiura, Sho Sumiya, Hanami Matsumoto, Daiki Ito, Kengo Nikai, Hiroshi Akitaya, Chihiro Ishioka, Kohei Oide, Takahiro Kanai, Yoshinori Uzawa, Yumiko Oasa, Tamás Tordai, Tonny Vanmunster, Sergey Yu. Shugarov, Masayuki Yamanaka, Mahito Sasada, Kengo Takagi , et al. (19 additional authors not shown)

Abstract: We report photometric and spectroscopic observations of the eclipsing SU UMa-type dwarf nova ASASSN-18aan. We observed the 2018 superoutburst with 2.3 mag brightening and found the orbital period ($P_{\rm orb}$) to be 0.149454(3) d, or 3.59 hr. This is longward of the period gap, establishing ASASSN-18aan as one of a small number of long-$P_{\rm orb}$ SU UMa-type dwarf novae. The estimated mass ra… ▽ More We report photometric and spectroscopic observations of the eclipsing SU UMa-type dwarf nova ASASSN-18aan. We observed the 2018 superoutburst with 2.3 mag brightening and found the orbital period ($P_{\rm orb}$) to be 0.149454(3) d, or 3.59 hr. This is longward of the period gap, establishing ASASSN-18aan as one of a small number of long-$P_{\rm orb}$ SU UMa-type dwarf novae. The estimated mass ratio, ($q=M_2/M_1 = 0.278(1)$), is almost identical to the upper limit of tidal instability by the 3:1 resonance. From eclipses, we found that the accretion disk at the onset of the superoutburst may reach the 3:1 resonance radius, suggesting that the superoutburst of ASASSN-18aan results from the tidal instability. Considering the case of long-$P_{\rm orb}$ WZ Sge-type dwarf novae, we suggest that the tidal dissipation at the tidal truncation radius is enough to induce SU UMa-like behavior in relatively high-$q$ systems such as SU UMa-type dwarf novae, but that this is no longer effective in low-$q$ systems such as WZ Sge-type dwarf novae. The unusual nature of the system extends to the secondary star, for which we find a spectral type of G9, much earlier than typical for the orbital period, and a secondary mass $M_2$ of around 0.18 M$_{\odot}$, smaller than expected for the orbital period and the secondary's spectral type. We also see indications of enhanced sodium abundance in the secondary's spectrum. Anomalously hot secondaries are seen in a modest number of other CVs and related objects. These systems evidently underwent significant nuclear evolution before the onset of mass transfer. In the case of ASASSN-18aan, this apparently resulted in a mass ratio lower than typically found at the system's $P_{\rm orb}$, which may account for the occurrence of a superoutburst at this relatively long period. △ Less

Submitted 8 February, 2021; v1 submitted 8 February, 2021; originally announced February 2021.

Comments: 17 pages, 14 figures, accepted for publication in PASJ

Optical follow-up observation for GW event S190510g using Subaru/Hyper Suprime-Cam

Authors: Takayuki Ohgami, Nozomu Tominaga, Yousuke Utsumi, Yuu Niino, Masaomi Tanaka, Smaranika Banerjee, Ryo Hamasaki, Michitoshi Yoshida, Tsuyoshi Terai, Yuhei Takagi, Tomoki Morokuma, Mahito Sasada, Hiroshi Akitaya, Naoki Yasuda, Kenshi Yanagisawa, Ryou Ohsawa

Abstract: A gravitational wave event, S190510g, which was classified as a binary-neutron-star coalescence at the time of preliminary alert, was detected by LIGO/Virgo collaboration on May 10, 2019. At 1.7 hours after the issue of its preliminary alert, we started a target-of-opportunity imaging observation in Y-band to search for its optical counterpart using the Hyper Suprime-Cam (HSC) on the Subaru Telesc… ▽ More A gravitational wave event, S190510g, which was classified as a binary-neutron-star coalescence at the time of preliminary alert, was detected by LIGO/Virgo collaboration on May 10, 2019. At 1.7 hours after the issue of its preliminary alert, we started a target-of-opportunity imaging observation in Y-band to search for its optical counterpart using the Hyper Suprime-Cam (HSC) on the Subaru Telescope. The observation covers a 118.8 deg$^2$ sky area corresponding to 11.6% confidence in the localization skymap released in the preliminary alert and 1.2% in the updated skymap. We divided the observed area into two fields based on the availability of HSC reference images. For the fields with the HSC reference images, we applied an image subtraction technique; for the fields without the HSC reference images, we sought individual HSC images by matching a catalog of observed objects with the PS1 catalog. The search depth is 22.28 mag in the former method and the limit of search depth is 21.3 mag in the latter method. Subsequently, we performed visual inspection and obtained 83 candidates using the former method and 50 candidates using the latter method. Since we have only the 1-day photometric data, we evaluated probability to be located inside the 3D skymap by estimating their distances with photometry of associated extended objects. We found three candidates are likely located inside the 3D skymap and concluded they could be an counterpart of S190510g, while most of 133 candidates were likely to be supernovae because the number density of candidates was consistent with the expected number of supernova detections. By comparing our observational depth with a light curve model of such a kilonova reproducing AT2017gfo, we show that early-deep observations with the Subaru/HSC can capture the rising phase of blue component of kilonova at the estimated distance of S190510g (~230 Mpc). △ Less

Submitted 4 January, 2021; originally announced January 2021.

Origin of the UV to X-ray emission of radio galaxy NGC 1275 explored by analyzing its variability

Authors: Fumiya Imazato, Yasushi Fukazawa, Mahito Sasada, Takanori Sakamoto

Abstract: We analyze the ultraviolet (UV) and X-ray data of NGC 1275 obtained with {\it Swift}/UVOT, XRT, BAT and {\it Fermi} Large Area Telescope over about 10 years to investigate the origin of the nuclear emission from NGC 1275. We confirm that the UV and soft/hard X-ray fluxes gradually increased along with the GeV gamma rays. At times, short-term variations in the UV or soft X-ray spectral regions show… ▽ More We analyze the ultraviolet (UV) and X-ray data of NGC 1275 obtained with {\it Swift}/UVOT, XRT, BAT and {\it Fermi} Large Area Telescope over about 10 years to investigate the origin of the nuclear emission from NGC 1275. We confirm that the UV and soft/hard X-ray fluxes gradually increased along with the GeV gamma rays. At times, short-term variations in the UV or soft X-ray spectral regions showed rapid variations correlated with the GeV gamma-rays. However there was no significant correlation between the UV and soft X-rays. The UV spectrum had a narrow spectral shape that could be represented by single-temperature blackbody radiation. These results could possibly indicate that the long-term variability of UV and X-ray emissions is caused by the jet, while the emissions from the accretion disk contribute to the UV and X-ray bands to some extent. △ Less

Submitted 20 November, 2020; originally announced November 2020.

Comments: 38 pages, 14 figures. Accepted in ApJ

Follow-up Observations for IceCube-170922A: Detection of Rapid Near-Infrared Variability and Intensive Monitoring of TXS 0506+056

Authors: Tomoki Morokuma, Yousuke Utsumi, Kouji Ohta, Masayuki Yamanaka, Koji S. Kawabata, Yoshiyuki Inoue, Masaomi Tanaka, Michitoshi Yoshida, Ryosuke Itoh, Mahito Sasada, Nozomu Tominaga, Hiroki Mori, Miho Kawabata, Tatsuya Nakaoka, Maiko Chogi, Taisei Abe, Ruochen Huang, Naoki Kawahara, Hiroki Kimura, Hiroki Nagashima, Kengo Takagi, Yuina Yamazaki, Wei Liu, Ryou Ohsawa, Shigeyuki Sako , et al. (13 additional authors not shown)

Abstract: We present our follow-up observations to search for an electromagnetic counterpart of the IceCube high-energy neutrino, IceCube-170922A. Monitoring observations of a likely counterpart, TXS 0506+056, are also described. First, we quickly took optical and near-infrared images of 7 flat-spectrum radio sources within the IceCube error region right after the neutrino detection and found a rapid flux d… ▽ More We present our follow-up observations to search for an electromagnetic counterpart of the IceCube high-energy neutrino, IceCube-170922A. Monitoring observations of a likely counterpart, TXS 0506+056, are also described. First, we quickly took optical and near-infrared images of 7 flat-spectrum radio sources within the IceCube error region right after the neutrino detection and found a rapid flux decline of TXS 0506+056 in Kanata/HONIR J-band data. Motivated by this discovery, intensive follow-up observations of TXS 0506+056 are continuously done, including our monitoring imaging observations, spectroscopic observations, and polarimetric observations in optical and near-infrared wavelengths. TXS 0506+056 shows a large amplitude (~1.0 mag) variability in a time scale of several days or longer, although no significant variability is detected in a time scale of a day or shorter. TXS 0506+056 also shows a bluer-when-brighter trend in optical and near-infrared wavelengths. Structure functions of variabilities are examined and indicate that TXS 0506+056 is not a special blazar in terms of optical variability. Polarization measurement results of TXS 0506+056 are also discussed. △ Less

Submitted 10 November, 2020; originally announced November 2020.

Comments: 19 pages, 10 figures, 4 tables. Accepted for publication in Publications of the Astronomical Society of Japan (PASJ)

Bi-infinite solutions for KdV- and Toda-type discrete integrable systems based on path encodings

Authors: David A. Croydon, Makiko Sasada, Satoshi Tsujimoto

Abstract: We define bi-infinite versions of four well-studied discrete integrable models, namely the ultra-discrete KdV equation, the discrete KdV equation, the ultra-discrete Toda equation, and the discrete Toda equation. For each equation, we show that there exists a unique solution to the initial value problem when the given data lies within a certain class, which includes the support of many shift ergod… ▽ More We define bi-infinite versions of four well-studied discrete integrable models, namely the ultra-discrete KdV equation, the discrete KdV equation, the ultra-discrete Toda equation, and the discrete Toda equation. For each equation, we show that there exists a unique solution to the initial value problem when the given data lies within a certain class, which includes the support of many shift ergodic measures. Our unified approach, which is also applicable to other integrable systems defined locally via lattice maps, involves the introduction of a path encoding (that is, a certain antiderivative) of the model configuration, for which we are able to describe the dynamics more generally than in previous work on finite size systems, periodic systems and semi-infinite systems. In particular, in each case we show that the behaviour of the system is characterized by a generalization of the classical 'Pitman's transformation' of reflection in the past maximum, which is well-known to probabilists. The picture presented here also provides a means to identify a natural 'carrier process' for configurations within the given class, and is convenient for checking that the systems we discuss are all-time reversible. Finally, we investigate links between the different systems, such as showing that bi-infinite all-time solutions for the ultra-discrete KdV (resp. Toda) equation may appear as ultra-discretizations of corresponding solutions for the discrete KdV (resp. Toda) equation. △ Less

Submitted 17 May, 2022; v1 submitted 1 November, 2020; originally announced November 2020.

MSC Class: 37K10 (primary); 35Q53; 37K60

Gravitational Test Beyond the First Post-Newtonian Order with the Shadow of the M87 Black Hole

Authors: Dimitrios Psaltis, Lia Medeiros, Pierre Christian, Feryal Ozel, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Keiichi Asada, Rebecca Azulay, David Ball, Mislav Balokovic, John Barrett, Dan Bintley, Lindy Blackburn, Wilfred Boland, Geoffrey C. Bower, Michael Bremer, Christiaan D. Brinkerink, Roger Brissenden, Silke Britzen, Dominique Broguiere, Thomas Bronzwaer, Do-Young Byun, John E. Carlstrom, Andrew Chael , et al. (163 additional authors not shown)

Abstract: The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the p… ▽ More The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the predicted black-hole shadows that are inconsistent with even the current EHT measurements. We use numerical calculations of regular, parametric, non-Kerr metrics to identify the common characteristic among these different parametrizations that control the predicted shadow size. We show that the shadow-size measurements place significant constraints on deviation parameters that control the second post-Newtonian and higher orders of each metric and are, therefore, inaccessible to weak-field tests. The new constraints are complementary to those imposed by observations of gravitational waves from stellar-mass sources. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: Physical Review Letters

Time-resolved spectroscopy and photometry of an M dwarf flare star YZ Canis Minoris with OISTER and TESS: Blue asymmetry in H$α$ line during the non-white light flare

Authors: Hiroyuki Maehara, Yuta Notsu, Kousuke Namekata, Satoshi Honda, Adam F. Kowalski, Noriyuki Katoh, Tomohito Ohshima, Kota Iida, Motoki Oeda, Katsuhiro L. Murata, Masayuki Yamanaka, Kengo Takagi, Mahito Sasada, Hiroshi Akitaya, Kai Ikuta, Soshi Okamoto, Daisaku Nogami, Kazunari Shibata

Abstract: In this paper, we present the results from spectroscopic and photometric observations of the M-type flare star YZ CMi in the framework of the Optical and Infrared Synergetic Telescopes for Education and Research (OISTER) collaborations during the Transiting Exoplanet Survey Satellite (TESS) observation period. We detected 145 white-light flares from the TESS light curve and 4 H$α$ flares from the… ▽ More In this paper, we present the results from spectroscopic and photometric observations of the M-type flare star YZ CMi in the framework of the Optical and Infrared Synergetic Telescopes for Education and Research (OISTER) collaborations during the Transiting Exoplanet Survey Satellite (TESS) observation period. We detected 145 white-light flares from the TESS light curve and 4 H$α$ flares from the OISTER observations performed between 2019-01-16 and 2019-01-18. Among them, 3 H$α$ flares were associated with white-light flares. However, one of them did not show clear brightening in continuum; during this flare, the H$α$ line exhibited blue-asymmetry which has lasted for $\sim 60$ min. The line of sight velocity of the blue-shifted component is $-80$ - $-100$ km s$^{-1}$. This suggests that there can be upward flows of chromospheric cool plasma even without detectable red/NIR continuum brightening. By assuming that the blue-asymmetry in H$α$ line was caused by a prominence eruption on YZ CMi, we estimated the mass and kinetic energy of the upward-moving material to be $10^{16}$ - $10^{18}$ g and $10^{29.5}$ - $10^{31.5}$ erg, respectively. The estimated mass is comparable to expectations from the empirical relation between the flare X-ray energy and mass of upward-moving material for stellar flares and solar CMEs. In contrast, the estimated kinetic energy for the non-white-light flare on YZ CMi is roughly $2$ orders of magnitude smaller than that expected from the relation between flare X-ray energy and kinetic energy for solar CMEs. This could be understood by the difference in the velocity between CMEs and prominence eruptions. △ Less

Submitted 1 October, 2020; v1 submitted 29 September, 2020; originally announced September 2020.

Comments: 39 pages, 15 figures, supplementary table (CSV format); accepted for publication in PASJ

Topological Structures of Large Scale Interacting Systems via Uniform Functions and Forms

Authors: Kenichi Bannai, Yukio Kametani, Makiko Sasada

Abstract: In this article, we investigate the topological structure of large scale interacting systems on infinite graphs, by constructing a suitable cohomology which we call the uniform cohomology. This provides a new perspective for the identification of macroscopic observables from the microscopic system. As a straightforward application of our theory when the underlying graph has a free action of a grou… ▽ More In this article, we investigate the topological structure of large scale interacting systems on infinite graphs, by constructing a suitable cohomology which we call the uniform cohomology. This provides a new perspective for the identification of macroscopic observables from the microscopic system. As a straightforward application of our theory when the underlying graph has a free action of a group, we prove a certain decomposition theorem for shift-invariant closed uniform forms. This result is a uniform version in a very general setting of the decomposition result for shift-invariant closed $L^2$-forms originally proposed by Varadhan, which has repeatedly played a key role in the proof of the hydrodynamic limits of nongradient large scale interacting systems. In a subsequent article, we use this result as a key to prove Varadhan's decomposition theorem for a general class of large scale interacting systems. △ Less

Submitted 22 March, 2022; v1 submitted 10 September, 2020; originally announced September 2020.

Comments: 61 pages, 6 figures. Updated title and Included reference to subsequent research concerning Varadhan's decomposition of $L^2$-forms, which is proved using the main result of this article

MSC Class: 82C22; 05C63; 55N91; 60J60; 70G40

Multi-wavelength photometry during the 2018 superoutburst of the WZ Sge-type dwarf nova EG Cancri

Authors: Mariko Kimura, Keisuke Isogai, Taichi Kato, Naoto Kojiguchi, Yasuyuki Wakamatsu, Ryuhei Ohnishi, Yuki Sugiura, Hanami Matsumoto, Sho Sumiya, Daiki Ito, Kengo Nikai, Katsura Matsumoto, Sergey Yu. Shugarov, Natalia Kathysheva, Hiroshi Itoh, Pavol A. Dubovsky, Igor Kudzej, Hiroshi Akitaya, Kohei Oide, Takahiro Kanai, Chihiro Ishioka, Yumiko Oasa, Tonny Vanmunster, Arto Oksanen, Tamás Tordai , et al. (23 additional authors not shown)

Abstract: We report on the multi-wavelength photometry of the 2018 superoutburst in EG Cnc. We have detected stage A superhumps and long-lasting late-stage superhumps via the optical photometry and have constrained the binary mass ratio and its possible range. The median value of the mass ratio is 0.048 and the upper limit is 0.057, which still implies that EG Cnc is one of the possible candidates for the p… ▽ More We report on the multi-wavelength photometry of the 2018 superoutburst in EG Cnc. We have detected stage A superhumps and long-lasting late-stage superhumps via the optical photometry and have constrained the binary mass ratio and its possible range. The median value of the mass ratio is 0.048 and the upper limit is 0.057, which still implies that EG Cnc is one of the possible candidates for the period bouncer. This object also showed multiple rebrightenings in this superoutburst, which are the same as those in its previous superoutburst in 1996--1997 despite the difference in the main superoutburst. This would represent that the rebrightening type is inherent to each object and is independent of the initial disk mass at the beginning of superoutbursts. We also found that $B-I$ and $J-K_{\rm S}$ colors were unusually red just before the rebrightening phase and became bluer during the quiescence between rebrightenings, which would mean that the low-temperature mass reservoir at the outermost disk accreted with time after the main superoutburst. Also, the ultraviolet flux was sensitive to rebrightenings as well as the optical flux, and the $U-B$ color became redder during the rebrightening phase, which would indicate that the inner disk became cooler when this object repeated rebrightenings. Our results thus basically support the idea that the cool mass reservoir in the outermost disk is responsible for rebrightenings. △ Less

Submitted 26 August, 2020; originally announced August 2020.

Comments: Accepted for publication in PASJ. 13 pages and 9 figures

Discrete integrable systems and Pitman's transformation

Authors: David A. Croydon, Makiko Sasada

Abstract: We survey recent work that relates Pitman's transformation to a variety of classical integrable systems, including the box-ball system, the ultra-discrete and discrete KdV equations, and the ultra-discrete and discrete Toda lattice equations. It is explained how this connection enables the dynamics of the integrable systems to be initiated from infinite configurations, which is important in the st… ▽ More We survey recent work that relates Pitman's transformation to a variety of classical integrable systems, including the box-ball system, the ultra-discrete and discrete KdV equations, and the ultra-discrete and discrete Toda lattice equations. It is explained how this connection enables the dynamics of the integrable systems to be initiated from infinite configurations, which is important in the study of invariant measures. In the special case of spatially independent and identically distributed configurations, progress on the latter topic is also reported. △ Less

Submitted 13 July, 2020; originally announced July 2020.

Detailed balance and invariant measures for discrete KdV- and Toda-type systems

Authors: David A. Croydon, Makiko Sasada

Abstract: In order to study the invariant measures of discrete KdV- and Toda-type systems, this article focusses on models, discretely indexed in space and time, whose dynamics are deterministic and defined locally via lattice equations. A detailed balance criterion is presented that, amongst the measures that describe spatially independent and identically/alternately distributed configurations, characteriz… ▽ More In order to study the invariant measures of discrete KdV- and Toda-type systems, this article focusses on models, discretely indexed in space and time, whose dynamics are deterministic and defined locally via lattice equations. A detailed balance criterion is presented that, amongst the measures that describe spatially independent and identically/alternately distributed configurations, characterizes those that are temporally invariant in distribution. A condition for establishing ergodicity of the dynamics is also given. These results are applied to various examples of discrete integrable systems, namely the ultra-discrete and discrete KdV equations, for which it is shown that the relevant invariant measures are of exponential/geometric and generalized inverse Gaussian form, respectively, as well as the ultra-discrete and discrete Toda lattice equations, for which the relevant invariant measures are found to be of exponential/geometric and gamma form. Ergodicity is demonstrated in the case of the KdV-type models. Links between the invariant measures of the different systems are presented, as are connections with stochastic integrable models and iterated random functions. Furthermore, a number of conjectures concerning the characterization of standard distributions are posed. △ Less

Submitted 9 December, 2021; v1 submitted 13 July, 2020; originally announced July 2020.

Calcium-rich Transient SN 2019ehk in A Star-Forming Environment: Yet Another Candidate for An Ultra-Stripped Envelope Supernova

Authors: Tatsuya Nakaoka, Keiichi Maeda, Masayuki Yamanaka, Masaomi Tanaka, Miho Kawabata, Takashi J. Moriya, Koji S. Kawabata, Nozomu Tominaga, Kengo Takagi, Fumiya Imazato, Tomoki Morokuma, Shigeyuki Sako, Ryou Ohsawa, Takashi Nagao, Ji-an Jiang, Umut Burgaz, Kenta Taguchi, Makoto Uemura, Hiroshi Akitaya, Mahito Sasada, Keisuke Isogai, Masaaki Otsuka, Hiroyuki Maehara

Abstract: We present optical and near-infrared observations of SN~Ib~2019ehk. We show that it evolved to a Ca-rich transient according to its spectral properties and evolution in late phases. It, however, shows a few distinguishable properties from the canonical Ca-rich transients: a short-duration first peak in the light curve, high peak luminosity, and association with a star-forming environment. Indeed,… ▽ More We present optical and near-infrared observations of SN~Ib~2019ehk. We show that it evolved to a Ca-rich transient according to its spectral properties and evolution in late phases. It, however, shows a few distinguishable properties from the canonical Ca-rich transients: a short-duration first peak in the light curve, high peak luminosity, and association with a star-forming environment. Indeed, some of these features are shared with iPTF14gqr and iPTF16hgs, which are candidates for a special class of core-collapse SNe (CCSNe): the so-called ultra-stripped envelope SNe, i.e., a relatively low-mass He (or C+O) star explosion in a binary as a precursor of double neutron star binaries. The estimated ejecta mass ($0.43 M_\odot$) and explosion energy ($1.7 \times 10^{50} $~erg) are consistent with this scenario. The analysis of the first peak suggests existence of dense circumstellar material in the vicinity of the progenitor, implying a CCSN origin. Based on these analyses, we suggest SN 2019ehk is another candidate for an ultra-stripped envelope SN. These ultra-stripped envelope SN candidates seem to form a subpopulation among Ca-rich transients, associated with young population. We propose that the key to distinguishing this population is the early first peak in their light curves. △ Less

Submitted 6 May, 2020; originally announced May 2020.

Comments: submitted to ApJ on Mar 27, 2020

SYMBA: An end-to-end VLBI synthetic data generation pipeline

Authors: F. Roelofs, M. Janssen, I. Natarajan, R. Deane, J. Davelaar, H. Olivares, O. Porth, S. N. Paine, K. L. Bouman, R. P. J. Tilanus, I. M. van Bemmel, H. Falcke, K. Akiyama, A. Alberdi, W. Alef, K. Asada, R. Azulay, A. Baczko, D. Ball, M. Baloković, J. Barrett, D. Bintley, L. Blackburn, W. Boland, G. C. Bower , et al. (183 additional authors not shown)

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

Submitted 2 April, 2020; originally announced April 2020.

Comments: 20 pages, 15 figures, accepted for publication in A&A

Generalized hydrodynamic limit for the box-ball system

Authors: David A. Croydon, Makiko Sasada

Abstract: We deduce a generalized hydrodynamic limit for the box-ball system, which explains how the densities of solitons of different sizes evolve asymptotically under Euler space-time scaling. To describe the limiting soliton flow, we introduce a continuous state-space analogue of the soliton decomposition of Ferrari, Nguyen, Rolla and Wang (cf. the original work of Takahashi and Satsuma), namely we rela… ▽ More We deduce a generalized hydrodynamic limit for the box-ball system, which explains how the densities of solitons of different sizes evolve asymptotically under Euler space-time scaling. To describe the limiting soliton flow, we introduce a continuous state-space analogue of the soliton decomposition of Ferrari, Nguyen, Rolla and Wang (cf. the original work of Takahashi and Satsuma), namely we relate the densities of solitons of given sizes in space to corresponding densities on a scale of 'effective distances', where the dynamics are linear. For smooth initial conditions, we further show that the resulting evolution of the soliton densities in space can alternatively be characterised by a partial differential equation, which naturally links the time-derivatives of the soliton densities and the 'effective speeds' of solitons locally. △ Less

Submitted 13 May, 2022; v1 submitted 13 March, 2020; originally announced March 2020.

MSC Class: 37B15 (primary); 82C22; 82C23; 82C70 (secondary)

SN 2018hna: 1987A-like supernova with a signature of shock breakout

Authors: Avinash Singh, D. K. Sahu, G. C. Anupama, Brajesh Kumar, Harsh Kumar, Masayuki Yamanaka, Petr V. Baklanov, Nozomu Tominaga, Sergei I. Blinnikov, Keiichi Maeda, Anirban Dutta, Varun Bhalerao, Ramya M. Anche, Sudhanshu Barway, Hiroshi Akitaya, Tatsuya Nakaoka, Miho Kawabata, Koji S Kawabata, Mahito Sasada, Kengo Takagi, Hiroyuki Maehara, Keisuke Isogai, Masaru Kino, Kenta Taguchi, Takashi Nagao

Abstract: High cadence ultraviolet, optical and near-infrared photometric and low-resolution spectroscopic observations of the peculiar Type II supernova (SN) 2018hna are presented. The early phase multiband light curves exhibit the adiabatic cooling envelope emission following the shock breakout up to ~14 days from the explosion. SN~2018hna has a rise time of $\sim$\,88 days in the V-band, similar to SN 19… ▽ More High cadence ultraviolet, optical and near-infrared photometric and low-resolution spectroscopic observations of the peculiar Type II supernova (SN) 2018hna are presented. The early phase multiband light curves exhibit the adiabatic cooling envelope emission following the shock breakout up to ~14 days from the explosion. SN~2018hna has a rise time of $\sim$\,88 days in the V-band, similar to SN 1987A. A $\rm^{56}Ni$ mass of ~0.087$\pm$0.004 $\rm M_{\odot}$ is inferred for SN 2018hna from its bolometric light curve. Hydrodynamical modelling of the cooling phase suggests a progenitor with a radius ~50 $\rm R_{\odot}$, a mass of ~14-20 $\rm M_{\odot}$ and explosion energy of ~1.7-2.9$\rm \times$ $\rm 10^{51}\ erg$. The smaller inferred radius of the progenitor than a standard red supergiant is indicative of a blue supergiant progenitor of SN 2018hna. A sub-solar metallicity (~0.3 $\rm Z_{\odot}$) is inferred for the host galaxy UGC 07534, concurrent with the low-metallicity environments of 1987A-like events. △ Less

Submitted 21 August, 2019; originally announced August 2019.

Comments: 8 pages, 5 figures, Accepted for Publication in ApJ Letters

Type Ia SN 2019ein: New Insights into the Similarities and diversities among High-Velocity SNe Ia

Authors: Miho Kawabata, Keiichi Maeda, Masayuki Yamanaka, Tatsuya Nakaoka, Koji Kawabata, Ryo Adachi, Hiroshi Akitaya, Umut Burgaz, Hidekazu Hanayama, Takashi Horiuchi, Ryohei Hosokawa, Kota Iida, Fumiya Imazato, Keisuke Isogai, Ji-an Jiang, Noriyuki Katoh, Hiroki Kimura, Masaru Kino, Daisuke Kuroda, Hiroyuki Maehara, Kazuya Matsubayashi, Kumiko Morihana, Katsuhiro Murata, Takashi Nagao, Masafumi Niwano , et al. (16 additional authors not shown)

Abstract: We present optical observations of type Ia supernova (SN) 2019ein, starting at 2 days after the estimated explosion date. The spectra and the light curves show that SN 2019ein belongs to the High-Velocity (HV) and Bload Line groups with relatively rapid decline in the light curves (Delta m15(B) = 1.36 +- 0.02 mag) and the short rise time (15.37 +- 0.55 days). The Si II 6355 velocity, associated wi… ▽ More We present optical observations of type Ia supernova (SN) 2019ein, starting at 2 days after the estimated explosion date. The spectra and the light curves show that SN 2019ein belongs to the High-Velocity (HV) and Bload Line groups with relatively rapid decline in the light curves (Delta m15(B) = 1.36 +- 0.02 mag) and the short rise time (15.37 +- 0.55 days). The Si II 6355 velocity, associated with a photospheric component but not with a detached high-velocity feature, reached ~ 20,000 km s-1 at 12 days before the B-band maximum. The line velocity however decreased very rapidly and smoothly toward the maximum light, where it was ~ 13,000 km s-1 as relatively low among HV SNe. This indicates that the speed of the spectral evolution of HV SNe Ia is correlated not only to the velocity at the maximum light, but also to the light curve decline rate like the case for Normal-Velocity (NV) SNe Ia. Spectral synthesis modeling shows that the outermost layer at > 17,000 km s-1 is well described by the O-Ne-C burning layer extending to at least 25,000 km s-1, and there is no unburnt carbon below 30,000 km s-1; these properties are largely consistent with the delayed detonation scenario, and are shared with the prototypical HV SN 2002bo despite the large difference in Delta m15(B). This structure is strikingly different from that derived for the well-studied NV SN 2011fe. We suggest that the relation between the mass of 56Ni (or Delta m15) and the extent of the O-Ne-C burning layer provides an important constraint on the explosion mechanism(s) of HV and NV SNe. △ Less

Submitted 16 March, 2020; v1 submitted 8 August, 2019; originally announced August 2019.

Comments: 15 pages, 12 figres, Accepted to ApJ

Duality between box-ball systems of finite box and/or carrier capacity

Authors: David A. Croydon, Makiko Sasada

Abstract: We construct the dynamics of the box-ball system with box capacity $J$ and carrier capacity $K$, which we abbreviate to BBS($J$,$K$), in the case of infinite initial configurations, and show that this system is dual to the analogous BBS($K$,$J$) model. Towards this end, we build on previous work for the original box-ball system, that is BBS($1$,$\infty$), to show that when the box capacity $J$ and… ▽ More We construct the dynamics of the box-ball system with box capacity $J$ and carrier capacity $K$, which we abbreviate to BBS($J$,$K$), in the case of infinite initial configurations, and show that this system is dual to the analogous BBS($K$,$J$) model. Towards this end, we build on previous work for the original box-ball system, that is BBS($1$,$\infty$), to show that when the box capacity $J$ and carrier capacity $K$ satisfy $J<K$ the dynamics can be represented by a Pitman-type transformation. These ideas are applied in the case of random initial configurations to show that the distributional properties of spatial stationarity and invariance under the BBS dynamics are dual. Moreover, for independent and identically distributed configurations, we derive a characterisation of invariant measures in terms of a detailed balance equation, which captures the duality of the system locally; this is used to find all invariant measures in this class. Finally, we deduce the speed of a tagged particle, and show that this also satisfies a natural duality relation. △ Less

Submitted 1 May, 2019; originally announced May 2019.

MSC Class: 37B15 (primary); 60G50; 60K35; 82B99 (secondary)