Showing 1–50 of 285 results for author: Takahashi, K

Low [O/Fe] Ratio in a Luminous Galaxy at the Early Cosmic Epoch ($z>10$): Signature of Short Delay Time or Bright Hypernovae/Pair-Instability Supernovae?

Authors: Minami Nakane, Masami Ouchi, Kimihiko Nakajima, Yuichi Harikane, Nozomu Tominaga, Koh Takahashi, Hiroto Yanagisawa, Kuria Watanabe, Ken'ichi Nomoto, Yuki Isobe, Daichi Kashino, Moka Nishigaki, Miho N. Ishigaki, Yoshiaki Ono, Yui Takeda

Abstract: We present an [O/Fe] ratio of a luminous galaxy GN-z11 at $z=10.60$ derived with the deep public JWST/NIRSpec data. We fit the medium-resolution grating data with the model spectra consisting of BPASS-stellar and CLOUDY-nebular spectra in the rest-frame UV wavelength ranges with Fe absorption lines, carefully masking out the other emission and absorption lines in the same manner as previous studie… ▽ More We present an [O/Fe] ratio of a luminous galaxy GN-z11 at $z=10.60$ derived with the deep public JWST/NIRSpec data. We fit the medium-resolution grating data with the model spectra consisting of BPASS-stellar and CLOUDY-nebular spectra in the rest-frame UV wavelength ranges with Fe absorption lines, carefully masking out the other emission and absorption lines in the same manner as previous studies conducted for lower redshift ($z\sim 2-6$) galaxies with oxygen abundance measurements. We obtain an Fe-rich abundance ratio $\mathrm{[O/Fe]}=-0.39^{+0.59}_{-0.36}$ that is confirmed with the independent deep prism data as well as by the classic 1978 index method. This [O/Fe] measurement is lower than that at $z\sim 2-3$ and comparable with those at $z\sim 6$. Because $z=10.60$ is an early epoch after the Big Bang ($\sim 430$ Myr) and the first star formation (likely $\sim 200$ Myr), there is a difficulty in efficiently producing Fe by Type Ia supernovae (SNeIa) requiring sufficient delay time for white-dwarf formation and gas accretion. The Fe-rich abundance ratio in GN-z11 suggests that the delay time is short, or that the major Fe enrichment is not accomplished by SNeIa but bright hypernovae (BrHNe) and/or pair-instability supernovae (PISNe), where the yield models of BrHNe and PISNe explain Fe, Ne, and O abundance ratios of GN-z11. The [O/Fe] measurement is not too low to rule out the connection between GN-z11 and globular clusters (GCs) previously suggested by the nitrogen abundance, but rather confirms the connection with a GC population at high [N/O] if a metal dilution process exists. △ Less

Submitted 19 July, 2024; originally announced July 2024.

Comments: 12 pages, 11 figures, submitted to ApJ

Isotropy of cosmic rays beyond $10^{20}$ eV favors their heavy mass composition

Authors: Telescope Array Collaboration, R. U. Abbasi, Y. Abe, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, B. G. Cheon, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon, N. Hayashida, H. He , et al. (118 additional authors not shown)

Abstract: We report an estimation of the injected mass composition of ultra-high energy cosmic rays (UHECRs) at energies higher than 10 EeV. The composition is inferred from an energy-dependent sky distribution of UHECR events observed by the Telescope Array surface detector by comparing it to the Large Scale Structure of the local Universe. In the case of negligible extra-galactic magnetic fields the resul… ▽ More We report an estimation of the injected mass composition of ultra-high energy cosmic rays (UHECRs) at energies higher than 10 EeV. The composition is inferred from an energy-dependent sky distribution of UHECR events observed by the Telescope Array surface detector by comparing it to the Large Scale Structure of the local Universe. In the case of negligible extra-galactic magnetic fields the results are consistent with a relatively heavy injected composition at E ~ 10 EeV that becomes lighter up to E ~ 100 EeV, while the composition at E > 100 EeV is very heavy. The latter is true even in the presence of highest experimentally allowed extra-galactic magnetic fields, while the composition at lower energies can be light if a strong EGMF is present. The effect of the uncertainty in the galactic magnetic field on these results is subdominant. △ Less

Submitted 3 July, 2024; v1 submitted 27 June, 2024; originally announced June 2024.

Comments: 8 pages, 3 figures, accepted for publication in PRL

Mass composition of ultra-high energy cosmic rays from distribution of their arrival directions with the Telescope Array

Authors: Telescope Array Collaboration, R. U. Abbasi, Y. Abe, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, B. G. Cheon, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon, N. Hayashida, H. He , et al. (118 additional authors not shown)

Abstract: We use a new method to estimate the injected mass composition of ultrahigh cosmic rays (UHECRs) at energies higher than 10 EeV. The method is based on comparison of the energy-dependent distribution of cosmic ray arrival directions as measured by the Telescope Array experiment (TA) with that calculated in a given putative model of UHECR under the assumption that sources trace the large-scale struc… ▽ More We use a new method to estimate the injected mass composition of ultrahigh cosmic rays (UHECRs) at energies higher than 10 EeV. The method is based on comparison of the energy-dependent distribution of cosmic ray arrival directions as measured by the Telescope Array experiment (TA) with that calculated in a given putative model of UHECR under the assumption that sources trace the large-scale structure (LSS) of the Universe. As we report in the companion letter, the TA data show large deflections with respect to the LSS which can be explained, assuming small extra-galactic magnetic fields (EGMF), by an intermediate composition changing to a heavy one (iron) in the highest energy bin. Here we show that these results are robust to uncertainties in UHECR injection spectra, the energy scale of the experiment and galactic magnetic fields (GMF). The assumption of weak EGMF, however, strongly affects this interpretation at all but the highest energies E > 100 EeV, where the remarkable isotropy of the data implies a heavy injected composition even in the case of strong EGMF. This result also holds if UHECR sources are as rare as $2 \times 10^{-5}$ Mpc$^{-3}$, that is the conservative lower limit for the source number density. △ Less

Submitted 3 July, 2024; v1 submitted 27 June, 2024; originally announced June 2024.

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

Observation of Declination Dependence in the Cosmic Ray Energy Spectrum

Authors: The Telescope Array Collaboration, R. U. Abbasi, T. Abu-Zayyad, M. Allen, J. W. Belz, D. R. Bergman, I. Buckland, W. Campbell, B. G. Cheon, K. Endo, A. Fedynitch, T. Fujii, K. Fujisue, K. Fujita, M. Fukushima, G. Furlich, Z. Gerber, N. Globus, W. Hanlon, N. Hayashida, H. He, K. Hibino, R. Higuchi, D. Ikeda, T. Ishii , et al. (101 additional authors not shown)

Abstract: We report on an observation of the difference between northern and southern skies of the ultrahigh energy cosmic ray energy spectrum with a significance of ${\sim}8σ$. We use measurements from the two largest experiments$\unicode{x2014}$the Telescope Array observing the northern hemisphere and the Pierre Auger Observatory viewing the southern hemisphere. Since the comparison of two measurements fr… ▽ More We report on an observation of the difference between northern and southern skies of the ultrahigh energy cosmic ray energy spectrum with a significance of ${\sim}8σ$. We use measurements from the two largest experiments$\unicode{x2014}$the Telescope Array observing the northern hemisphere and the Pierre Auger Observatory viewing the southern hemisphere. Since the comparison of two measurements from different observatories introduces the issue of possible systematic differences between detectors and analyses, we validate the methodology of the comparison by examining the region of the sky where the apertures of the two observatories overlap. Although the spectra differ in this region, we find that there is only a $1.8σ$ difference between the spectrum measurements when anisotropic regions are removed and a fiducial cut in the aperture is applied. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 8 pages, 6 figures

(In)stability of the black hole greybody factors and ringdowns against a small-bump correction

Authors: Naritaka Oshita, Kazufumi Takahashi, Shinji Mukohyama

Abstract: Recently, it has been proposed that the black hole greybody factors can be important to model ringdown spectral amplitudes. We study the stability of greybody factors against a small-bump correction in the perturbation equation. We find (I) that the greybody factor is stable in the frequency region relevant to ringdown and (II) that it is destabilized at higher frequencies, especially for a sharpe… ▽ More Recently, it has been proposed that the black hole greybody factors can be important to model ringdown spectral amplitudes. We study the stability of greybody factors against a small-bump correction in the perturbation equation. We find (I) that the greybody factor is stable in the frequency region relevant to ringdown and (II) that it is destabilized at higher frequencies, especially for a sharper bump correction. This behavior is similar to the case of higher overtones, which is also very sensitive to a small correction. We clarify this (in)stability with the WKB analysis. As the greybody factor is stable at the frequency region relevant to the main part of ringdown, we conclude that the greybody factor is suitable to model ringdown amplitude. In order to investigate a bump correction in a self-consistent manner, we consider the small-bump correction that can be realized in the general framework of effective field theory of black hole perturbations. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 14 pages, 4 figures

Report number: YITP-24-69, IPMU24-0025, RIKEN-iTHEMS-Report-24

CMB spectrum in unified EFT of dark energy: scalar-tensor and vector-tensor theories

Authors: Katsuki Aoki, Mohammad Ali Gorji, Takashi Hiramatsu, Shinji Mukohyama, Masroor C. Pookkillath, Kazufumi Takahashi

Abstract: We study the cosmic microwave background (CMB) radiation in the unified description of the effective field theory (EFT) of dark energy that accommodates both scalar-tensor and vector-tensor theories. The boundaries of different classes of theories are universally parameterised by a new EFT parameter $α_V$ characterising the vectorial nature of dark energy and a set of consistency relations associa… ▽ More We study the cosmic microwave background (CMB) radiation in the unified description of the effective field theory (EFT) of dark energy that accommodates both scalar-tensor and vector-tensor theories. The boundaries of different classes of theories are universally parameterised by a new EFT parameter $α_V$ characterising the vectorial nature of dark energy and a set of consistency relations associated with the global/local shift symmetry. After implementing the equations of motion in a Boltzmann code, as a demonstration, we compute the CMB power spectrum based on the $w$CDM background with the EFT parameterisation of perturbations and a concrete Horndeski/generalised Proca theory. We show that the vectorial nature generically prevents modifications of gravity in the CMB spectrum. On the other hand, while the shift symmetry is less significant in the perturbation equations unless the background is close to the $Λ$CDM, it requires that the effective equation of state of dark energy is in the phantom region $w_{\rm DE}<-1$. The latter is particularly interesting in light of the latest result of the DESI+CMB combination as the observational verification of $w_{\rm DE}>-1$ can rule out shift-symmetric theories including vector-tensor theories in one shot. △ Less

Submitted 7 May, 2024; originally announced May 2024.

Comments: 21 pages, 4 figures, 1 table

Report number: YITP-24-56, RUP-24-8, IPMU24-0016

Strong He I Emission Lines in High N/O Galaxies at $z \sim 6$ Identified in JWST Spectra: High He/H Abundance Ratios or High Electron Densities?

Authors: Hiroto Yanagisawa, Masami Ouchi, Kuria Watanabe, Akinori Matsumoto, Kimihiko Nakajima, Hidenobu Yajima, Kentaro Nagamine, Koh Takahashi, Minami Nakane, Nozomu Tominaga, Hiroya Umeda, Hajime Fukushima, Yuichi Harikane, Yuki Isobe, Yoshiaki Ono, Yi Xu, Yechi Zhang

Abstract: We present HeI/H$β$-flux and He/H-abundance ratios in three JWST galaxies with significant constraints on N/O-abundance ratios, GS-NDG-9422, RXCJ2248-ID, and GLASS150008 at $z\sim 6$ mostly with the spectroscopic coverage from HeI$λ$4471 and HeII$λ$4686 to HeI$λ$7065, comparing with 68 local-dwarf galaxies. We find that these high-$z$ galaxies present strong HeI emission with HeI/H$β$ flux ratios… ▽ More We present HeI/H$β$-flux and He/H-abundance ratios in three JWST galaxies with significant constraints on N/O-abundance ratios, GS-NDG-9422, RXCJ2248-ID, and GLASS150008 at $z\sim 6$ mostly with the spectroscopic coverage from HeI$λ$4471 and HeII$λ$4686 to HeI$λ$7065, comparing with 68 local-dwarf galaxies. We find that these high-$z$ galaxies present strong HeI emission with HeI/H$β$ flux ratios generally larger than those of local-dwarf galaxies. We derive He/H with all of the detected HeI, HeII, and $2-3$ hydrogen Balmer lines in the same manner as the local He/H determination conducted for cosmology studies. These high-$z$ galaxies show He overabundance He/H$\gtrsim 0.10$ or high electron density $n_\mathrm{e}\sim 10^{3-4}$ cm$^{-3}$ much larger than local values at low O/H, $12+\log \mathrm{(O/H)}=7-8$. In contrast, we obtain low He/H and $n_\mathrm{e}$ values for our local-dwarf galaxies by the same technique with the same helium and hydrogen lines, and confirm that the difference between the high-$z$ and local-dwarf galaxies are not mimicked by systematics. While two scenarios of 1) He overabundance and 2) high electron density are not clearly concluded, we find that there is a positive correlation on the He/H-N/O or $n_\mathrm{e}$-N/O plane by the comparison of the high-$z$ and local-dwarf galaxies. The scenario 1) suggests that the overabundant helium and nitrogen are not explained by the standard chemical enrichment of core-collapse supernovae, but the CNO-cycle products and equilibrium ratios, respectively. The scenario 2) indicates that the strong helium lines are originated from the central dense clouds of the high-$z$ galaxies by excessive collisional excitation. △ Less

Submitted 2 May, 2024; originally announced May 2024.

Comments: 17 pages, 10 figures, submitted to ApJ

Black holes, multiple propagation speeds, and energy extraction

Authors: Vitor Cardoso, Shinji Mukohyama, Naritaka Oshita, Kazufumi Takahashi

Abstract: The Standard Model of particle physics predicts the speed of light to be a universal speed of propagation of massless carriers. However, other possibilities exist -- including Lorentz-violating theories -- where different fundamental fields travel at different speeds. Black holes are interesting probes of such physics, as distinct fields would probe different horizons. Here, we build an exact spac… ▽ More The Standard Model of particle physics predicts the speed of light to be a universal speed of propagation of massless carriers. However, other possibilities exist -- including Lorentz-violating theories -- where different fundamental fields travel at different speeds. Black holes are interesting probes of such physics, as distinct fields would probe different horizons. Here, we build an exact spacetime for two interacting scalar fields which have different propagation speeds. One of these fields is able to probe the black hole interior of the other, giving rise to energy extraction from the black hole and a characteristic late-time relaxation. Our results provide further stimulus to the search for extra degrees of freedom, black hole instability, and extra ringdown modes in gravitational-wave events. △ Less

Submitted 23 June, 2024; v1 submitted 8 April, 2024; originally announced April 2024.

Comments: 19 pages, 19 figures, matches published version

Report number: YITP-24-44, IPMU24-0014, RIKEN-iTHEMS-Report-24

Journal ref: Phys. Rev. D 109, 124036 (2024)

Destratification in the Progenitor Interior of the Mg-rich Supernova Remnant N49B

Authors: Toshiki Sato, Kai Matsunaga, Ryo Sawada, Koh Takahashi, Yudai Suwa, John P. Hughes, Hiroyuki Uchida, Takuto Narita, Hideyuki Umeda

Abstract: Simulations of pre-supernova evolution suggest that some intense shell burning can be so active that, in extreme cases, it can merge with the outer shell, changing the initial conditions for the supernova explosion. However, such violent activity in the interior of stars has been difficult to confirm from observations of stars. Here we propose that the elemental composition of O-rich ejecta in sup… ▽ More Simulations of pre-supernova evolution suggest that some intense shell burning can be so active that, in extreme cases, it can merge with the outer shell, changing the initial conditions for the supernova explosion. However, such violent activity in the interior of stars has been difficult to confirm from observations of stars. Here we propose that the elemental composition of O-rich ejecta in supernova remnants can be a tool to test for this kind of intense shell burning activity in the final stages of progenitor evolution. As an example, we discuss the origin of "Mg-rich" ejecta in the supernova remnant N49B. A high Mg/Ne mass ratio $\gtrsim 1$ suggests that the Ne- or O-burning shell has broken into or merged with the outer shell before the collapse. Such Mg-rich (or Ne-poor) ejecta has been identified in some other supernova remnants, supporting the idea that some destratification process, such as a shell merger, does indeed occur in the interiors of some massive stars, although they may not be the majority. Our results suggest that X-ray observations of O-rich ejecta in core-collapse supernova remnants will be a unique tool to probe the shell burning activity during the final stage of a massive star's interior. △ Less

Submitted 6 March, 2024; originally announced March 2024.

Comments: Submitted to ApJ and received positive peer review. This is the version in the first submission. 11 pages, 7 figures, 2 tables

A search for auroral radio emission from βPictoris b

Authors: Yuta Shiohira, Yuka Fujii, Hajime Kita, Tomoki Kimura, Yuka Terada, Keitaro Takahashi

Abstract: Magnetized exoplanets can serve as the source of auroral radio emissions, allowing us to characterize the magnetospheric properties of these planets. Successful detections of auroral radio emissions from brown dwarfs, as well as from Jupiter, suggest that Jupiter-like planets in distant orbits may also generate radio emissions through a similar mechanism. In this study, we present our search for 2… ▽ More Magnetized exoplanets can serve as the source of auroral radio emissions, allowing us to characterize the magnetospheric properties of these planets. Successful detections of auroral radio emissions from brown dwarfs, as well as from Jupiter, suggest that Jupiter-like planets in distant orbits may also generate radio emissions through a similar mechanism. In this study, we present our search for 250-500 MHz emissions from $β$ Pictoris b, one of the most extensively studied young Jupiter-like planets. We conducted the search using the upgraded Giant Metrewave Radio Telescope (uGMRT). Despite the favourable orbital inclination, no signal was detected, putting 3$σ$ upper limits on the radiation at 0.18 mJy. We translate this limit into constraints on the ionospheric and magnetospheric parameters, assuming that the emission is powered by the Hill current system. While the upper limit is larger by a factor of a few than the nominal estimate of radio intensity, we put constraints on the magnetospheric and ionospheric parameters. △ Less

Submitted 23 December, 2023; originally announced December 2023.

Comments: 10 pages, 5 figures. Accepted for publication in MNRAS

Exact solution for rotating black holes in parity-violating gravity

Authors: Hiroaki W. H. Tahara, Kazufumi Takahashi, Masato Minamitsuji, Hayato Motohashi

Abstract: It has recently been pointed out that one can construct invertible conformal transformations with a parity-violating conformal factor, which can be employed to generate a novel class of parity-violating ghost-free metric theories from general relativity. We obtain exact solutions for rotating black holes in such theories by performing the conformal transformation on the Kerr solution in general re… ▽ More It has recently been pointed out that one can construct invertible conformal transformations with a parity-violating conformal factor, which can be employed to generate a novel class of parity-violating ghost-free metric theories from general relativity. We obtain exact solutions for rotating black holes in such theories by performing the conformal transformation on the Kerr solution in general relativity, which we dub conformal Kerr solutions. We explore the geodesic motion of a test particle in the conformal Kerr spacetime. While null geodesics remain the same as those in the Kerr spacetime, timelike geodesics exhibit interesting differences due to an effective external force caused by the parity-violating conformal factor. △ Less

Submitted 20 May, 2024; v1 submitted 19 December, 2023; originally announced December 2023.

Comments: 19 pages, 8 figures

Report number: RUP-23-26, YITP-23-147

Journal ref: Progress of Theoretical and Experimental Physics, Volume 2024, Issue 5, May 2024, 053E02

Low-frequency pulse-jitter measurement with the uGMRT I : PSR J0437$-$4715

Authors: Tomonosuke Kikunaga, Shinnosuke Hisano, Neelam Dhanda Batra, Shantanu Desai, Bhal Chandra Joshi, Manjari Bagchi, T. Prabu, Keitaro Takahashi, Swetha Arumugam, Adarsh Bathula, Subhajit Dandapat, Debabrata Deb, Churchil Dwivedi, Yashwant Gupta, Shebin Jose Jacob, Fazal Kareem, Nobleson K, Pragna Mamidipaka, Avinash Kumar Paladi, Arul Pandian B, Prerna Rana, Jaikhomba Singha, Aman Srivastava, Mayuresh Surnis, Pratik Tarafdar

Abstract: High-precision pulsar timing observations are limited in their accuracy by the jitter noise that appears in the arrival time of pulses. Therefore, it is important to systematically characterise the amplitude of the jitter noise and its variation with frequency. In this paper, we provide jitter measurements from low-frequency wideband observations of PSR J0437$-$4715 using data obtained as part of… ▽ More High-precision pulsar timing observations are limited in their accuracy by the jitter noise that appears in the arrival time of pulses. Therefore, it is important to systematically characterise the amplitude of the jitter noise and its variation with frequency. In this paper, we provide jitter measurements from low-frequency wideband observations of PSR J0437$-$4715 using data obtained as part of the Indian Pulsar Timing Array experiment. We were able to detect jitter in both the 300 - 500 MHz and 1260 - 1460 MHz observations of the upgraded Giant Metrewave Radio Telescope (uGMRT). The former is the first jitter measurement for this pulsar below 700 MHz, and the latter is in good agreement with results from previous studies. In addition, at 300 - 500 MHz, we investigated the frequency dependence of the jitter by calculating the jitter for each sub-banded arrival time of pulses. We found that the jitter amplitude increases with frequency. This trend is opposite as compared to previous studies, indicating that there is a turnover at intermediate frequencies. It will be possible to investigate this in more detail with uGMRT observations at 550 - 750 MHz and future high sensitive wideband observations from next generation telescopes, such as the Square Kilometre Array. We also explored the effect of jitter on the high precision dispersion measure (DM) measurements derived from short duration observations. We find that even though the DM precision will be better at lower frequencies due to the smaller amplitude of jitter noise, it will limit the DM precision for high signal-to-noise observations, which are of short durations. This limitation can be overcome by integrating for a long enough duration optimised for a given pulsar. △ Less

Submitted 18 April, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: 13 pages, 12 figures, 3 tables, accepted for Publication of the Astronomical Society of Australia

Iterative extraction of overtones from black hole ringdown

Authors: Kazuto Takahashi, Hayato Motohashi

Abstract: Extraction of multiple quasinormal modes from ringdown gravitational waves emitted from a binary black hole coalescence is a touchstone to test whether a remnant black hole is described by the Kerr spacetime in general relativity. However, it is not straightforward to check the consistency between the ringdown signal and the quasinormal mode frequencies predicted by the linear perturbation theory.… ▽ More Extraction of multiple quasinormal modes from ringdown gravitational waves emitted from a binary black hole coalescence is a touchstone to test whether a remnant black hole is described by the Kerr spacetime in general relativity. However, it is not straightforward to check the consistency between the ringdown signal and the quasinormal mode frequencies predicted by the linear perturbation theory. While the longest-lived mode can be extracted in a stable manner, the higher overtones damp more quickly and hence the fitting of overtones tends to end up with the overfit. To improve the extraction of overtones, we propose an iterative procedure consisting of fitting and subtraction of the longest-lived mode of the ringdown waveform in the time domain. Through the analyses of the mock waveform and numerical relativity waveform, we clarify that the iterative procedure allows us to extract the overtones in a more stable manner. △ Less

Submitted 21 November, 2023; originally announced November 2023.

Comments: 38 pages, 38 figures

Variety of disk wind-driven explosions in massive rotating stars

Authors: Ludovica Crosato Menegazzi, Sho Fujibayashi, Koh Takahashi, Ayako Ishii

Abstract: We perform a set of two-dimensional, non-relativistic, hydrodynamics simulations for supernova-like explosion associated with stellar core collapse of rotating massive stars to a system of a black hole and a disk connected by the transfer of matter and angular momentum. Our model of the central engine also includes the contribution of the disk wind. In this work, we specifically investigate the wi… ▽ More We perform a set of two-dimensional, non-relativistic, hydrodynamics simulations for supernova-like explosion associated with stellar core collapse of rotating massive stars to a system of a black hole and a disk connected by the transfer of matter and angular momentum. Our model of the central engine also includes the contribution of the disk wind. In this work, we specifically investigate the wind-driven explosion of rotating, large-mass progenitor stars with the zero-age main-sequence mass of $M_\mathrm{ZAMS}=20\,M_\odot$ from arXiv:2008.09132 . This study is carried out using the open-source hydrodynamic code Athena++, for which we implement a method to calculate self-gravity for axially symmetric density distributions. We, then, investigate the explosion properties and the $^{56}$Ni production as a function of (varying) some features of the wind injection. We find a large variety of explosion energy with $E_\mathrm{expl}$ ranging from $\sim 0.049\times10^{51}$~erg to $\sim 34\times10^{51}$~erg and ejecta mass $M_\mathrm{ej}$ from 0.58 to 6 $M_\odot$, which shows a bimodal distribution in high- and low-energy branches. We demonstrate that the resulting outcome of a highly- or sub-energetic explosion for a certain stellar structure is mainly determined by the competition between the ram pressure of the injected matter and that of the infalling envelope. In the nucleosynthesis analysis the $^{56}$Ni mass produced in our models goes from $< 0.2~M_\odot$ in the sub-energetic explosions to $2.1~M_\odot$ in the highly-energetic ones. These results are consistent with the observational data of stripped-envelope and high-energy SNe such as broad-lined type Ic SNe. However, we find a tighter correlation between the explosion energy and the ejecta mass than that observationally measured. △ Less

Submitted 28 February, 2024; v1 submitted 7 November, 2023; originally announced November 2023.

Comments: published by MNRAS, 18 pages, 15 figures, 2 tables, comments welcome

Journal ref: Monthly Notices of the Royal Astronomical Society Volume 529 Issue 1 March 2024 Pages 178 - 195

Using low-frequency scatter-broadening measurements for precision estimates of dispersion measures

Authors: Jaikhomba Singha, Bhal Chandra Joshi, M. A. Krishnakumar, Fazal Kareem, Adarsh Bathula, Churchil Dwivedi, Shebin Jose Jacob, Shantanu Desai, Pratik Tarafdar, P. Arumugam, Swetha Arumugam, Manjari Bagchi, Neelam Dhanda Batra, Subhajit Dandapat, Debabrata Deb, Jyotijwal Debnath, A Gopakumar, Yashwant Gupta, Shinnosuke Hisano, Ryo Kato, Tomonosuke Kikunaga, Piyush Marmat, K. Nobleson, Avinash K. Paladi, Arul Pandian B. , et al. (6 additional authors not shown)

Abstract: A pulsar's pulse profile gets broadened at low frequencies due to dispersion along the line of sight or due to multi-path propagation. The dynamic nature of the interstellar medium makes both of these effects time-dependent and introduces slowly varying time delays in the measured times-of-arrival similar to those introduced by passing gravitational waves. In this article, we present a new method… ▽ More A pulsar's pulse profile gets broadened at low frequencies due to dispersion along the line of sight or due to multi-path propagation. The dynamic nature of the interstellar medium makes both of these effects time-dependent and introduces slowly varying time delays in the measured times-of-arrival similar to those introduced by passing gravitational waves. In this article, we present a new method to correct for such delays by obtaining unbiased dispersion measure (DM) measurements by using low-frequency estimates of the scattering parameters. We evaluate this method by comparing the obtained DM estimates with those, where scatter-broadening is ignored using simulated data. A bias is seen in the estimated DMs for simulated data with pulse-broadening with a larger variability for a data set with a variable frequency scaling index, $α$, as compared to that assuming a Kolmogorov turbulence. Application of the proposed method removes this bias robustly for data with band averaged signal-to-noise ratio larger than 100. We report, for the first time, the measurements of the scatter-broadening time and $α$ from analysis of PSR J1643$-$1224, observed with upgraded Giant Metrewave Radio Telescope as part of the Indian Pulsar Timing Array experiment. These scattering parameters were found to vary with epoch and $α$ was different from that expected for Kolmogorov turbulence. Finally, we present the DM time-series after application of the new technique to PSR J1643$-$1224. △ Less

Submitted 28 September, 2023; originally announced September 2023.

Comments: 9 pages, 8 figures, Submitted to MNRAS

Diagnosis of Circumstellar Matter Structure in Interaction-powered Supernovae with Hydrogen Line Feature

Authors: Ayako T. Ishii, Yuki Takei, Daichi Tsuna, Toshikazu Shigeyama, Koh Takahashi

Abstract: Some supernovae (SNe) are powered by collision of the SN ejecta with a dense circumstellar matter (CSM). Their emission spectra show characteristic line shapes of combined broad emission and narrow P-Cyg lines, which should closely relate to the CSM structure and the mass-loss mechanism that creates the dense CSM. We quantitatively investigate the relationship between the line shape and the CSM st… ▽ More Some supernovae (SNe) are powered by collision of the SN ejecta with a dense circumstellar matter (CSM). Their emission spectra show characteristic line shapes of combined broad emission and narrow P-Cyg lines, which should closely relate to the CSM structure and the mass-loss mechanism that creates the dense CSM. We quantitatively investigate the relationship between the line shape and the CSM structure by Monte Carlo radiative transfer simulations, considering two representative cases of dense CSM formed by steady and eruptive mass loss. Comparing the H$α$ emission between the two cases, we find that a narrow P-Cyg line appears in the eruptive case while it does not appear in the steady case, due to the difference in the velocity gradient in the dense CSM. We also reproduce the blue-shifted photon excess observed in some Type IIn SNe, which is formed by photon transport across the shock wave and find the relationship between the velocity of the shocked matter and the amount of the blue shift of the photon excess. We conclude that the presence or absence of narrow P-Cyg lines can distinguish the mass loss mechanism, and suggest high-resolution spectroscopic observations with $λ/ Δλ\gtrsim 10^4$ after the light curve peak for applying this diagnostic method. △ Less

Submitted 23 January, 2024; v1 submitted 11 September, 2023; originally announced September 2023.

Comments: 25 pages, 14 figures, Accepted for publication in The Astrophysical Journal

Comparing recent PTA results on the nanohertz stochastic gravitational wave background

Authors: The International Pulsar Timing Array Collaboration, G. Agazie, J. Antoniadis, A. Anumarlapudi, A. M. Archibald, P. Arumugam, S. Arumugam, Z. Arzoumanian, J. Askew, S. Babak, M. Bagchi, M. Bailes, A. -S. Bak Nielsen, P. T. Baker, C. G. Bassa, A. Bathula, B. Bécsy, A. Berthereau, N. D. R. Bhat, L. Blecha, M. Bonetti, E. Bortolas, A. Brazier, P. R. Brook, M. Burgay , et al. (220 additional authors not shown)

Abstract: The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTA… ▽ More The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within $1σ$. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings and Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which will involve not just adding in additional pulsars, but also including data from all three PTAs where any given pulsar is timed by more than as single PTA. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 21 pages, 9 figures, submitted to ApJ

Precision of localization of single gravitational-wave source with pulsar timing array

Authors: Ryo Kato, Keitaro Takahashi

Abstract: Pulsar Timing Arrays (PTAs) are expected to be able to detect gravitational waves (GWs) from individual supermassive black hole binaries in the near future. In order to identify the host galaxy of a gravitational wave source, the angular resolution of PTAs should be much better than that expected from the conventional methodology of PTAs. We study the potential usefulness of precise pulsar-distanc… ▽ More Pulsar Timing Arrays (PTAs) are expected to be able to detect gravitational waves (GWs) from individual supermassive black hole binaries in the near future. In order to identify the host galaxy of a gravitational wave source, the angular resolution of PTAs should be much better than that expected from the conventional methodology of PTAs. We study the potential usefulness of precise pulsar-distance measurements in the determination of the sky location of a single GW source. Precise distance information from external observations such as astrometry by Very Long Baseline Interferometry is incorporated as priors in the PTA analysis and we evaluate the precision of the sky location of a GW source by simulating PTA data of 12 milli-second pulsars with only the GW signal and the Gaussian white noise in the timing residuals. We show that only a few pulsars with a distance precision of 1 pc will improve the precision of the source location by more than 1 order in the presence of white noise of 10 ns. △ Less

Submitted 19 December, 2023; v1 submitted 20 August, 2023; originally announced August 2023.

Comments: 13 pages, 7 figures

Journal ref: Phys. Rev. D 108, 123535, 2023

Invertible disformal transformations with arbitrary higher-order derivatives

Authors: Kazufumi Takahashi

Abstract: Invertible disformal transformations serve as a useful tool to explore ghost-free scalar-tensor theories. In this paper, we construct a generalization of invertible disformal transformations that involves arbitrary higher-order covariant derivatives of the scalar field. As a result, we obtain a more general class of ghost-free scalar-tensor theories than ever. Notably, our generalization is such t… ▽ More Invertible disformal transformations serve as a useful tool to explore ghost-free scalar-tensor theories. In this paper, we construct a generalization of invertible disformal transformations that involves arbitrary higher-order covariant derivatives of the scalar field. As a result, we obtain a more general class of ghost-free scalar-tensor theories than ever. Notably, our generalization is such that matter fields can be consistently coupled to these theories without introducing an unwanted extra degree of freedom in the unitary gauge. △ Less

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

Comments: 9 pages

Report number: YITP-23-91

Journal ref: Phys. Rev. D 108, 084031 (2023)

The second data release from the European Pulsar Timing Array: IV. Implications for massive black holes, dark matter and the early Universe

Authors: J. Antoniadis, P. Arumugam, S. Arumugam, P. Auclair, S. Babak, M. Bagchi, A. -S. Bak Nielsen, E. Barausse, C. G. Bassa, A. Bathula, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, C. Caprini, A. Chalumeau, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, M. Crisostomi, S. Dandapat, D. Deb , et al. (89 additional authors not shown)

Abstract: The European Pulsar Timing Array (EPTA) and Indian Pulsar Timing Array (InPTA) collaborations have measured a low-frequency common signal in the combination of their second and first data releases respectively, with the correlation properties of a gravitational wave background (GWB). Such signal may have its origin in a number of physical processes including a cosmic population of inspiralling sup… ▽ More The European Pulsar Timing Array (EPTA) and Indian Pulsar Timing Array (InPTA) collaborations have measured a low-frequency common signal in the combination of their second and first data releases respectively, with the correlation properties of a gravitational wave background (GWB). Such signal may have its origin in a number of physical processes including a cosmic population of inspiralling supermassive black hole binaries (SMBHBs); inflation, phase transitions, cosmic strings and tensor mode generation by non-linear evolution of scalar perturbations in the early Universe; oscillations of the Galactic potential in the presence of ultra-light dark matter (ULDM). At the current stage of emerging evidence, it is impossible to discriminate among the different origins. Therefore, in this paper, we consider each process separately, and investigate the implications of the signal under the hypothesis that it is generated by that specific process. We find that the signal is consistent with a cosmic population of inspiralling SMBHBs, and its relatively high amplitude can be used to place constraints on binary merger timescales and the SMBH-host galaxy scaling relations. If this origin is confirmed, this is the first direct evidence that SMBHBs merge in nature, adding an important observational piece to the puzzle of structure formation and galaxy evolution. As for early Universe processes, the measurement would place tight constraints on the cosmic string tension and on the level of turbulence developed by first-order phase transitions. Other processes would require non-standard scenarios, such as a blue-tilted inflationary spectrum or an excess in the primordial spectrum of scalar perturbations at large wavenumbers. Finally, a ULDM origin of the detected signal is disfavoured, which leads to direct constraints on the abundance of ULDM in our Galaxy. △ Less

Submitted 15 May, 2024; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 30 pages, 23 figures, replaced to match the version published in Astronomy & Astrophysics, note the change in the numbering order in the series (now paper IV)

The second data release from the European Pulsar Timing Array V. Search for continuous gravitational wave signals

Authors: J. Antoniadis, P. Arumugam, S. Arumugam, S. Babak, M. Bagchi, A. S. Bak Nielsen, C. G. Bassa, A. Bathula, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, A. Chalumeau, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, S. Dandapat, D. Deb, S. Desai, G. Desvignes, N. Dhanda-Batra, C. Dwivedi , et al. (75 additional authors not shown)

Abstract: We present the results of a search for continuous gravitational wave signals (CGWs) in the second data release (DR2) of the European Pulsar Timing Array (EPTA) collaboration. The most significant candidate event from this search has a gravitational wave frequency of 4-5 nHz. Such a signal could be generated by a supermassive black hole binary (SMBHB) in the local Universe. We present the results o… ▽ More We present the results of a search for continuous gravitational wave signals (CGWs) in the second data release (DR2) of the European Pulsar Timing Array (EPTA) collaboration. The most significant candidate event from this search has a gravitational wave frequency of 4-5 nHz. Such a signal could be generated by a supermassive black hole binary (SMBHB) in the local Universe. We present the results of a follow-up analysis of this candidate using both Bayesian and frequentist methods. The Bayesian analysis gives a Bayes factor of 4 in favor of the presence of the CGW over a common uncorrelated noise process, while the frequentist analysis estimates the p-value of the candidate to be 1%, also assuming the presence of common uncorrelated red noise. However, comparing a model that includes both a CGW and a gravitational wave background (GWB) to a GWB only, the Bayes factor in favour of the CGW model is only 0.7. Therefore, we cannot conclusively determine the origin of the observed feature, but we cannot rule it out as a CGW source. We present results of simulations that demonstrate that data containing a weak gravitational wave background can be misinterpreted as data including a CGW and vice versa, providing two plausible explanations of the EPTA DR2 data. Further investigations combining data from all PTA collaborations will be needed to reveal the true origin of this feature. △ Less

Submitted 25 June, 2024; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 13 figures, 15 pages, accepted

The second data release from the European Pulsar Timing Array II. Customised pulsar noise models for spatially correlated gravitational waves

Authors: J. Antoniadis, P. Arumugam, S. Arumugam, S. Babak, M. Bagchi, A. S. Bak Nielsen, C. G. Bassa, A. Bathula, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, A. Chalumeau, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, S. Dandapat, D. Deb, S. Desai, G. Desvignes, N. Dhanda-Batra, C. Dwivedi , et al. (73 additional authors not shown)

Abstract: The nanohertz gravitational wave background (GWB) is expected to be an aggregate signal of an ensemble of gravitational waves emitted predominantly by a large population of coalescing supermassive black hole binaries in the centres of merging galaxies. Pulsar timing arrays, ensembles of extremely stable pulsars, are the most precise experiments capable of detecting this background. However, the su… ▽ More The nanohertz gravitational wave background (GWB) is expected to be an aggregate signal of an ensemble of gravitational waves emitted predominantly by a large population of coalescing supermassive black hole binaries in the centres of merging galaxies. Pulsar timing arrays, ensembles of extremely stable pulsars, are the most precise experiments capable of detecting this background. However, the subtle imprints that the GWB induces on pulsar timing data are obscured by many sources of noise. These must be carefully characterized to increase the sensitivity to the GWB. In this paper, we present a novel technique to estimate the optimal number of frequency coefficients for modelling achromatic and chromatic noise and perform model selection. We also incorporate a new model to fit for scattering variations in the pulsar timing package temponest and created realistic simulations of the European Pulsar Timing Array (EPTA) datasets that allowed us to test the efficacy of our noise modelling algorithms. We present an in-depth analysis of the noise properties of 25 millisecond pulsars (MSPs) that form the second data release (DR2) of the EPTA and investigate the effect of incorporating low-frequency data from the Indian PTA collaboration. We use enterprise and temponest packages to compare noise models with those reported with the EPTA DR1. We find that, while in some pulsars we can successfully disentangle chromatic from achromatic noise owing to the wider frequency coverage in DR2, in others the noise models evolve in a more complicated way. We also find evidence of long-term scattering variations in PSR J1600$-$3053. Through our simulations, we identify intrinsic biases in our current noise analysis techniques and discuss their effect on GWB searches. The results presented here directly help improve sensitivity to the GWB and are already being used as part of global PTA efforts. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 20 pages, 6 figures, 9 tables

Journal ref: A&A 678, A49 (2023)

The second data release from the European Pulsar Timing Array III. Search for gravitational wave signals

Authors: J. Antoniadis, P. Arumugam, S. Arumugam, S. Babak, M. Bagchi, A. -S. Bak Nielsen, C. G. Bassa, A. Bathula, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, A. Chalumeau, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, S. Dandapat, D. Deb, S. Desai, G. Desvignes, N. Dhanda-Batra, C. Dwivedi , et al. (73 additional authors not shown)

Abstract: We present the results of the search for an isotropic stochastic gravitational wave background (GWB) at nanohertz frequencies using the second data release of the European Pulsar Timing Array (EPTA) for 25 millisecond pulsars and a combination with the first data release of the Indian Pulsar Timing Array (InPTA). We analysed (i) the full 24.7-year EPTA data set, (ii) its 10.3-year subset based on… ▽ More We present the results of the search for an isotropic stochastic gravitational wave background (GWB) at nanohertz frequencies using the second data release of the European Pulsar Timing Array (EPTA) for 25 millisecond pulsars and a combination with the first data release of the Indian Pulsar Timing Array (InPTA). We analysed (i) the full 24.7-year EPTA data set, (ii) its 10.3-year subset based on modern observing systems, (iii) the combination of the full data set with the first data release of the InPTA for ten commonly timed millisecond pulsars, and (iv) the combination of the 10.3-year subset with the InPTA data. These combinations allowed us to probe the contributions of instrumental noise and interstellar propagation effects. With the full data set, we find marginal evidence for a GWB, with a Bayes factor of four and a false alarm probability of $4\%$. With the 10.3-year subset, we report evidence for a GWB, with a Bayes factor of $60$ and a false alarm probability of about $0.1\%$ ($\gtrsim 3σ$ significance). The addition of the InPTA data yields results that are broadly consistent with the EPTA-only data sets, with the benefit of better noise modelling. Analyses were performed with different data processing pipelines to test the consistency of the results from independent software packages. The inferred spectrum from the latest EPTA data from new generation observing systems is rather uncertain and in mild tension with the common signal measured in the full data set. However, if the spectral index is fixed at 13/3, the two data sets give a similar amplitude of ($2.5\pm0.7)\times10^{-15}$ at a reference frequency of $1\,{\rm yr}^{-1}$. By continuing our detection efforts as part of the International Pulsar Timing Array (IPTA), we expect to be able to improve the measurement of spatial correlations and better characterise this signal in the coming years. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 21 pages, 14 figures, 4 appendix figures, accepted for publication in A&A

Journal ref: A&A 678, A50 (2023)

EMPRESS. XIII. Chemical Enrichments of Young Galaxies Near and Far at z ~ 0 and 4-10: Fe/O, Ar/O, S/O, and N/O Measurements with Chemical Evolution Model Comparisons

Authors: Kuria Watanabe, Masami Ouchi, Kimihiko Nakajima, Yuki Isobe, Nozomu Tominaga, Akihiro Suzuki, Miho N. Ishigaki, Ken'ichi Nomoto, Koh Takahashi, Yuichi Harikane, Shun Hatano, Haruka Kusakabe, Takashi J. Moriya, Moka Nishigaki, Yoshiaki Ono, Masato Onodera, Yuma Sugahara

Abstract: We present gas-phase elemental abundance ratios of 7 local extremely metal-poor galaxies (EMPGs) including our new Keck/LRIS spectroscopy determinations together with 33 JWST $z\sim 4-10$ star-forming galaxies in the literature, and compare chemical evolution models. We develop chemical evolution models with the yields of core-collapse supernovae (CCSNe), Type Ia supernovae, hypernovae (HNe), and… ▽ More We present gas-phase elemental abundance ratios of 7 local extremely metal-poor galaxies (EMPGs) including our new Keck/LRIS spectroscopy determinations together with 33 JWST $z\sim 4-10$ star-forming galaxies in the literature, and compare chemical evolution models. We develop chemical evolution models with the yields of core-collapse supernovae (CCSNe), Type Ia supernovae, hypernovae (HNe), and pair-instability supernovae (PISNe), and compare the EMPGs and high-$z$ galaxies in conjunction with dust depletion contributions. We find that high Fe/O values of EMPGs can (cannot) be explained by PISN metal enrichments (CCSN/HN enrichments even with the mixing-and-fallback mechanism enhancing iron abundance), while that the observed Ar/O and S/O values are much smaller than the predictions of the PISN models. The abundance ratios of the EMPGs can be explained by the combination of Type Ia SNe and CCSNe/HNe whose inner layers of argon and sulfur mostly fallback, which are comparable with Sculptor stellar chemical abundance distribution, suggesting that early chemical enrichment is taken place in the EMPGs. Comparing our chemical evolution models with the star-forming galaxies at $z\sim 4-10$, we find that the Ar/O and S/O ratios of the high-$z$ galaxies are comparable with those of the CCSNe/HNe models, while majority of the high-$z$ galaxies do not have constraints good enough to rule out contributions from PISNe. The high N/O ratio recently reported in GN-z11 cannot be explained even by rotating PISNe, but could be reproduced by the winds of rotating Wolf Rayet stars that end up as a direct collapse. △ Less

Submitted 11 March, 2024; v1 submitted 3 May, 2023; originally announced May 2023.

Multi-band Extension of the Wideband Timing Technique

Authors: Avinash Kumar Paladi, Churchil Dwivedi, Prerna Rana, Nobleson K, Abhimanyu Susobhanan, Bhal Chandra Joshi, Pratik Tarafdar, Debabrata Deb, Swetha Arumugam, A Gopakumar, M A Krishnakumar, Neelam Dhanda Batra, Jyotijwal Debnath, Fazal Kareem, Paramasivan Arumugam, Manjari Bagchi, Adarsh Bathula, Subhajit Dandapat, Shantanu Desai, Yashwant Gupta, Shinnosuke Hisano, Divyansh Kharbanda, Tomonosuke Kikunaga, Neel Kolhe, Yogesh Maan , et al. (5 additional authors not shown)

Abstract: The wideband timing technique enables the high-precision simultaneous estimation of pulsar Times of Arrival (ToAs) and Dispersion Measures (DMs) while effectively modeling frequency-dependent profile evolution. We present two novel independent methods that extend the standard wideband technique to handle simultaneous multi-band pulsar data incorporating profile evolution over a larger frequency sp… ▽ More The wideband timing technique enables the high-precision simultaneous estimation of pulsar Times of Arrival (ToAs) and Dispersion Measures (DMs) while effectively modeling frequency-dependent profile evolution. We present two novel independent methods that extend the standard wideband technique to handle simultaneous multi-band pulsar data incorporating profile evolution over a larger frequency span to estimate DMs and ToAs with enhanced precision. We implement the wideband likelihood using the libstempo python interface to perform wideband timing in the tempo2 framework. We present the application of these techniques to the dataset of fourteen millisecond pulsars observed simultaneously in Band 3 (300 - 500 MHz) and Band 5 (1260 - 1460 MHz) of the upgraded Giant Metrewave Radio Telescope (uGMRT) with a large band gap of 760 MHz as a part of the Indian Pulsar Timing Array (InPTA) campaign. We achieve increased ToA and DM precision and sub-microsecond root mean square post-fit timing residuals by combining simultaneous multi-band pulsar observations done in non-contiguous bands for the first time using our novel techniques. △ Less

Submitted 8 November, 2023; v1 submitted 25 April, 2023; originally announced April 2023.

Comments: Published in MNRAS

Effective description of generalized disformal theories

Authors: Kazufumi Takahashi, Masato Minamitsuji, Hayato Motohashi

Abstract: Generalized disformal transformations enable us to construct the generalized disformal Horndeski theories, which form the most general class of ghost-free scalar-tensor theories to this date. We extend the effective field theory (EFT) of cosmological perturbations to incorporate these generalized disformal Horndeski theories. The main difference from the conventional EFT is that our extended EFT i… ▽ More Generalized disformal transformations enable us to construct the generalized disformal Horndeski theories, which form the most general class of ghost-free scalar-tensor theories to this date. We extend the effective field theory (EFT) of cosmological perturbations to incorporate these generalized disformal Horndeski theories. The main difference from the conventional EFT is that our extended EFT involves operators with higher spatial derivatives of the lapse function. Our EFT also accommodates the generalized disformal transformation of U-DHOST theories. △ Less

Submitted 4 July, 2023; v1 submitted 17 April, 2023; originally announced April 2023.

Comments: 16 pages

Report number: YITP-23-46

Journal ref: JCAP 07 (2023) 009

Noise analysis of the Indian Pulsar Timing Array data release I

Authors: Aman Srivastava, Shantanu Desai, Neel Kolhe, Mayuresh Surnis, Bhal Chandra Joshi, Abhimanyu Susobhanan, Aurélien Chalumeau, Shinnosuke Hisano, Nobleson K., Swetha Arumugam, Divyansh Kharbanda, Jaikhomba Singha, Pratik Tarafdar, P Arumugam, Manjari Bagchi, Adarsh Bathula, Subhajit Dandapat, Lankeswar Dey, Churchil Dwivedi, Raghav Girgaonkar, A. Gopakumar, Yashwant Gupta, Tomonosuke Kikunaga, M. A. Krishnakumar, Kuo Liu , et al. (6 additional authors not shown)

Abstract: The Indian Pulsar Timing Array (InPTA) collaboration has recently made its first official data release (DR1) for a sample of 14 pulsars using 3.5 years of uGMRT observations. We present the results of single-pulsar noise analysis for each of these 14 pulsars using the InPTA DR1. For this purpose, we consider white noise, achromatic red noise, dispersion measure (DM) variations, and scattering vari… ▽ More The Indian Pulsar Timing Array (InPTA) collaboration has recently made its first official data release (DR1) for a sample of 14 pulsars using 3.5 years of uGMRT observations. We present the results of single-pulsar noise analysis for each of these 14 pulsars using the InPTA DR1. For this purpose, we consider white noise, achromatic red noise, dispersion measure (DM) variations, and scattering variations in our analysis. We apply Bayesian model selection to obtain the preferred noise models among these for each pulsar. For PSR J1600$-$3053, we find no evidence of DM and scattering variations, while for PSR J1909$-$3744, we find no significant scattering variations. Properties vary dramatically among pulsars. For example, we find a strong chromatic noise with chromatic index $\sim$ 2.9 for PSR J1939+2134, indicating the possibility of a scattering index that doesn't agree with that expected for a Kolmogorov scattering medium consistent with similar results for millisecond pulsars in past studies. Despite the relatively short time baseline, the noise models broadly agree with the other PTAs and provide, at the same time, well-constrained DM and scattering variations. △ Less

Submitted 16 June, 2023; v1 submitted 21 March, 2023; originally announced March 2023.

Comments: Accepted for publication in PRD, 30 pages, 17 figures, 4 tables

Introduction to Faraday tomography and its future prospects

Authors: Keitaro Takahashi

Abstract: Faraday tomography is a new method of the study of cosmic magnetic fields enabled by broadband low-frequency radio observations. By Faraday tomography, it is possible to obtain the Faraday dispersion function which contains information on the line-of-sight distributions of magnetic fields, thermal electron density, and cosmic-ray electron density by measuring the polarization spectrum from a sourc… ▽ More Faraday tomography is a new method of the study of cosmic magnetic fields enabled by broadband low-frequency radio observations. By Faraday tomography, it is possible to obtain the Faraday dispersion function which contains information on the line-of-sight distributions of magnetic fields, thermal electron density, and cosmic-ray electron density by measuring the polarization spectrum from a source of synchrotron radiation over a wide band. Furthermore, by combining it with 2-dimensional imaging, Faraday tomography allows us to explore the 3-dimensional structure of polarization sources. The application of Faraday tomography has been active in the last 20 years, when broadband observation has become technically feasible. However, the Faraday dispersion function is mathematically the Fourier transform of the polarization spectrum, and since the observable band is finite, it is impossible to obtain a complete Faraday dispersion function by performing Fourier transform. In addition, the Faraday dispersion function does not directly reflect the distribution of magnetic field, thermal-electron density, and cosmic-ray electron density in the physical space, and its physical interpretation is not straightforward. Despite these two difficult problems, Faraday tomography is attracting much attention because it has great potential as a new method for studying cosmic magnetic fields and magnetized plasmas. In particular, the next-generation radio telescope SKA (Square Kilometre Array) is capable of polarization observation with unprecedented sensitivity and broad bands, and the application of Faraday tomography is expected to make dramatic progress in the field of cosmic magnetic fields. In this review, we explain the basics of Faraday tomography with simple and instructive examples. Then representative algorithms to realize Faraday tomography are introduced and finally some applications are shown. △ Less

Submitted 14 March, 2023; originally announced March 2023.

Comments: Review article, 35 pages, 32 figures

Journal ref: Publications of the Astronomical Society of Japan, Volume 75, Issue Supplement 1, February 2023, Pages S50-S84

Monotonicity of the cores of massive stars

Authors: Koh Takahashi, Tomoya Takiwaki, Takashi Yoshida

Abstract: Massive stars are linked with diverse astronomical processes and objects including star formation, supernovae and their remnants, cosmic rays, interstellar media, and galaxy evolution. Understanding their properties is of primary importance for modern astronomy, and finding simple rules that characterize them is especially useful. However, theoretical simulations have not yet realized such relatio… ▽ More Massive stars are linked with diverse astronomical processes and objects including star formation, supernovae and their remnants, cosmic rays, interstellar media, and galaxy evolution. Understanding their properties is of primary importance for modern astronomy, and finding simple rules that characterize them is especially useful. However, theoretical simulations have not yet realized such relations, instead finding that the late evolutionary phases are significantly affected by a complicated interplay between nuclear reactions, chemical mixing, and neutrino radiation, leading to non-monotonic initial mass dependencies of the iron core mass and the compactness parameter. We conduct a set of stellar evolution simulations, in which evolutions of He star models are followed until their central densities uniformly reach 10$^{10}$ g cm$^{-3}$, and analyze their final structures as well as their evolutionary properties including the lifetime, surface radius change, and presumable fates after core collapse. Based on the homogeneous data set, we have found that monotonicity is inherent in the cores of massive stars. We show that not only the density, entropy, and chemical distributions, but also their lifetimes and explosion properties such as the proto-neutron-star mass and the explosion energy can be simultaneously ordered into a monotonic sequence. This monotonicity can be regarded as an empirical principle that characterizes the cores of massive stars. △ Less

Submitted 6 February, 2023; originally announced February 2023.

Comments: 31 pages. Accepted for publication in ApJ

Evolution and explosions of metal enriched supermassive stars: proton rich general relativistic instability supernovae

Authors: Chris Nagele, Hideyuki Umeda, Koh Takahashi

Abstract: The assembly of supermassive black holes poses a challenge primarily because of observed quasars at high redshift, but additionally because of the current lack of observations of intermediate mass black holes. One plausible scenario for creating supermassive black holes is direct collapse triggered by the merger of two gas rich galaxies. This scenario allows the creation of supermassive stars with… ▽ More The assembly of supermassive black holes poses a challenge primarily because of observed quasars at high redshift, but additionally because of the current lack of observations of intermediate mass black holes. One plausible scenario for creating supermassive black holes is direct collapse triggered by the merger of two gas rich galaxies. This scenario allows the creation of supermassive stars with solar metallicity, where the enhanced metallicity is enabled by extremely rapid accretion. We investigate the behavior of metal enriched supermassive stars which collapse due to the general relativistic radial instability during hydrogen burning. These stars contain both hydrogen and metals and thus may explode due to the CNO cycle (carbon-nitrogen-oxygen) and the rp process (rapid proton capture). We perform a suite of stellar evolution simulations for a range of masses and metallicities, both including and neglecting mass loss. We evaluate the stability of these supermassive stars by solving the pulsation equation in general relativity. When the stars becomes unstable, we perform 1D general relativistic hydrodynamical simulations coupled to a 153 isotope nuclear network with cooling from neutrino reactions, in order to determine if the stars explode. If the stars do explode, we post process the nucleosynthesis using a 514 isotope network which includes additional proton rich isotopes. These explosions are characterized by enhanced nitrogen and intermediate mass elements ($16\geq\rm{A}\geq25$), and suppressed light elements ($8\geq\rm{A}\geq14$), and we comment on recent observations of super-solar nitrogen in GN-z11. △ Less

Submitted 20 May, 2023; v1 submitted 5 January, 2023; originally announced January 2023.

Cosmological gravity probes: connecting recent theoretical developments to forthcoming observations

Authors: Shun Arai, Katsuki Aoki, Yuji Chinone, Rampei Kimura, Tsutomu Kobayashi, Hironao Miyatake, Daisuke Yamauchi, Shuichiro Yokoyama, Kazuyuki Akitsu, Takashi Hiramatsu, Shin'ichi Hirano, Ryotaro Kase, Taishi Katsuragawa, Yosuke Kobayashi, Toshiya Namikawa, Takahiro Nishimichi, Teppei Okumura, Maresuke Shiraishi, Masato Shirasaki, Tomomi Sunayama, Kazufumi Takahashi, Atsushi Taruya, Junsei Tokuda

Abstract: Since the discovery of the accelerated expansion of the present Universe, significant theoretical developments have been made in the area of modified gravity. In the meantime, cosmological observations have been providing more high-quality data, allowing us to explore gravity on cosmological scales. To bridge the recent theoretical developments and observations, we present an overview of a variety… ▽ More Since the discovery of the accelerated expansion of the present Universe, significant theoretical developments have been made in the area of modified gravity. In the meantime, cosmological observations have been providing more high-quality data, allowing us to explore gravity on cosmological scales. To bridge the recent theoretical developments and observations, we present an overview of a variety of modified theories of gravity and the cosmological observables in the cosmic microwave background and large-scale structure, supplemented with a summary of predictions for cosmological observables derived from cosmological perturbations and sophisticated numerical studies. We specifically consider scalar-tensor theories in the Horndeski and DHOST family, massive gravity/bigravity, vector-tensor theories, metric-affine gravity, and cuscuton/minimally-modified gravity, and discuss the current status of those theories with emphasis on their physical motivations, validity, appealing features, the level of maturity, and calculability. We conclude that the Horndeski theory is one of the most well-developed theories of modified gravity, although several remaining issues are left for future observations. The paper aims to help to develop strategies for testing gravity with ongoing and forthcoming cosmological observations. △ Less

Submitted 18 December, 2022; originally announced December 2022.

Comments: 115 pages, 8 figures

Pulsations of primordial supermassive stars induced by a general relativistic instability; visible to JWST at z$>$12

Authors: Chris Nagele, Hideyuki Umeda, Koh Takahashi, Keiichi Maeda

Abstract: The origin of high-redshift quasars and their supermassive black hole engines is unclear. One promising solution is the collapse of a primordial supermassive star. Observational confirmation of this scenario may be challenging, but a general relativistic instability supernova provides one avenue for such. Previous studies have found that a general relativistic instability supernova has a potential… ▽ More The origin of high-redshift quasars and their supermassive black hole engines is unclear. One promising solution is the collapse of a primordial supermassive star. Observational confirmation of this scenario may be challenging, but a general relativistic instability supernova provides one avenue for such. Previous studies have found that a general relativistic instability supernova has a potentially decades-long plateau phase visible to JWST at high redshift. In this work, we examine stars with mass just below the general relativistic instability supernova mass range. These stars pulsate, ejecting a portion of their envelopes. They then contract quasi-statically back to an equilibrium temperature, at which point they again become unstable and pulsate once more. Because each pulse consumes a small amount of the available nuclear fuel, there exists the possibility of multiple pulsations. We present simulations of the contracting phase, the pulsation, and the light-curve phase. We find that the lower mass pulsating models are even brighter than the higher mass supernovae because the pulsations occur in the late helium burning phase when the stars have extremely large radii. The fact that the pulsations are more luminous and occur in a wider mass range than the supernovae bodes well for observation. △ Less

Submitted 30 January, 2023; v1 submitted 16 October, 2022; originally announced October 2022.

Eleven-year, 22-year and ~90-year solar cycles discovered in nitrate concentrations in a Dome Fuji (Antarctica) ice core

Authors: Yuko Motizuki, Yoichi Nakai, Kazuya Takahashi, Takashi Imamura, Hideaki Motoyama

Abstract: Ice cores are known to yield information about astronomical phenomena as well as information about past climate. We report time series analyses of annually resolved nitrate variations in an ice core, drilled at the Dome Fuji station in East Antarctica, corresponding to the period from CE 1610 to 1904. Our analyses revealed clear evidence of ~11, ~22, and ~90 year periodicities, comparable to the r… ▽ More Ice cores are known to yield information about astronomical phenomena as well as information about past climate. We report time series analyses of annually resolved nitrate variations in an ice core, drilled at the Dome Fuji station in East Antarctica, corresponding to the period from CE 1610 to 1904. Our analyses revealed clear evidence of ~11, ~22, and ~90 year periodicities, comparable to the respective periodicities of the well-known Schwabe, Hale, and Gleissberg solar cycles. Our results show for the first time that nitrate concentrations in an ice core can be used as a proxy for past solar activity on decadal to multidecadal time scales. Furthermore, 11-year and 22-year periodicities were detected in nitrate variations even during the Maunder Minimum (1645-1715), when sunspots were almost absent. This discovery may support cyclic behavior of the solar dynamo during the grand solar minimum. △ Less

Submitted 22 September, 2022; originally announced September 2022.

Comments: Submitted to Proceedings of the Japan Academy, Series B

Generalized disformal Horndeski theories: cosmological perturbations and consistent matter coupling

Authors: Kazufumi Takahashi, Masato Minamitsuji, Hayato Motohashi

Abstract: Invertible disformal transformations are a useful tool to investigate ghost-free scalar-tensor theories. By performing a higher-derivative generalization of the invertible disformal transformation on Horndeski theories, we construct a novel class of ghost-free scalar-tensor theories, which we dub generalized disformal Horndeski theories. Specifically, these theories lie beyond the quadratic/cubic… ▽ More Invertible disformal transformations are a useful tool to investigate ghost-free scalar-tensor theories. By performing a higher-derivative generalization of the invertible disformal transformation on Horndeski theories, we construct a novel class of ghost-free scalar-tensor theories, which we dub generalized disformal Horndeski theories. Specifically, these theories lie beyond the quadratic/cubic DHOST class. We explore cosmological perturbations to identify a subclass where gravitational waves propagate at the speed of light and clarify the conditions for the absence of ghost/gradient instabilities for tensor and scalar perturbations. We also investigate the conditions under which a matter field can be consistently coupled to these theories without introducing unwanted extra degrees of freedom. △ Less

Submitted 13 January, 2023; v1 submitted 5 September, 2022; originally announced September 2022.

Comments: 23 pages, 2 figures

Report number: YITP-22-89

Journal ref: PTEP 2023, 013E01 (2023)

On the Potential of Faraday Tomography to Identify Shock Structures in Supernova Remnants

Authors: Shinsuke Ideguchi, Tsuyoshi Inoue, Takuya Akahori, Keitaro Takahashi

Abstract: Knowledge about the magnetic fields in supernova remnants (SNRs) is of paramount importance for constraining Galactic cosmic ray acceleration models. It could also indirectly provide information on the interstellar magnetic fields. In this paper, we predict the Faraday dispersion functions (FDFs) of SNRs for the first time. For this study, we use the results of three dimensional (3D) ideal magneto… ▽ More Knowledge about the magnetic fields in supernova remnants (SNRs) is of paramount importance for constraining Galactic cosmic ray acceleration models. It could also indirectly provide information on the interstellar magnetic fields. In this paper, we predict the Faraday dispersion functions (FDFs) of SNRs for the first time. For this study, we use the results of three dimensional (3D) ideal magnetohydrodynamic (MHD) simulations of SNRs expanding into a weak, regular magnetic field. We present the intrinsic FDFs of the shocked region of SNRs for different viewing angles. We find that the FDFs are generally Faraday complex, which implies that conventional rotation measure study is not sufficient to obtain the information on the magnetic fields in the shocked region and Faraday tomography is necessary. We also show that the FDF allows to derive the physical-depth distribution of polarization intensity when the line of sight is parallel to the initial magnetic field orientation. Furthermore, we demonstrate that the location of contact discontinuity can be identified from the radial profile of the width of the FDF with the accuracy of 0.1-0.2 pc. △ Less

Submitted 22 August, 2022; originally announced August 2022.

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 513, Issue 3, July 2022, Pages 3289-3301

Probing particle acceleration at trans-relativistic shocks with off-axis gamma-ray burst afterglows

Authors: Kazuya Takahashi, Kunihito Ioka, Yutaka Ohira, Hendrik J. van Eerten

Abstract: Particle acceleration is expected to be different between relativistic and non-relativistic collisionless shocks. We show that electromagnetic counterparts to gravitational waves (GWs), gamma-ray burst (GRB) afterglows, are ideal targets for observing trans-relativistic evolution of accelerated electron distribution because the GWs spot nearby GRBs with off-axis jets, otherwise missed in gamma-ray… ▽ More Particle acceleration is expected to be different between relativistic and non-relativistic collisionless shocks. We show that electromagnetic counterparts to gravitational waves (GWs), gamma-ray burst (GRB) afterglows, are ideal targets for observing trans-relativistic evolution of accelerated electron distribution because the GWs spot nearby GRBs with off-axis jets, otherwise missed in gamma-ray observations. We find that the relativistic spectral slope begins to change steeply near the peak time of the light curve and approaches the non-relativistic limit in about 10 times the peak time. The trans-relativistic evolution of the afterglow synchrotron spectrum is consistent with GRB 170817A observations within errors, and will be measurable in similar but more distant events at a GW horizon $\sim 200$ Mpc in a denser environment. We roughly estimate that such events represent a fraction of 10-50 per cent of the GRB 170817A-like off-axis short GRBs. We also find that the spectral evolution does not depend on the jet structure if their light curves are similar to each other. △ Less

Submitted 18 October, 2022; v1 submitted 12 August, 2022; originally announced August 2022.

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

Nanohertz Gravitational Wave Astronomy during the SKA Era: An InPTA perspective

Authors: Bhal Chandra Joshi, Achamveedu Gopakumar, Arul Pandian, Thiagaraj Prabu, Lankeswar Dey, Manjari Bagchi, Shantanu Desai, Pratik Tarafdar, Prerna Rana, Yogesh Maan, Neelam Dhanda Batra, Raghav Girgaonkar, Nikita Agarwal, Paramasivan Arumugam, Sarmistha Banik, Avishek Basu, Adarsh Bathula, Subhajit Dandapat, Yashwant Gupta, Shinnosuke Hisano, Ryo Kato, Divyansh Kharbanda, Tomonosuke Kikunaga, Neel Kolhe, M. A. Krishnakumar , et al. (12 additional authors not shown)

Abstract: Decades long monitoring of millisecond pulsars, which exhibit highly stable rotational periods, in pulsar timing array experiments is on the threshold of discovering nanohertz stochastic gravitational wave background. This paper describes the Indian Pulsar timing array (InPTA) experiment, which employs the upgraded Giant Metrewave Radio Telescope (uGMRT) for timing an ensemble of millisecond pulsa… ▽ More Decades long monitoring of millisecond pulsars, which exhibit highly stable rotational periods, in pulsar timing array experiments is on the threshold of discovering nanohertz stochastic gravitational wave background. This paper describes the Indian Pulsar timing array (InPTA) experiment, which employs the upgraded Giant Metrewave Radio Telescope (uGMRT) for timing an ensemble of millisecond pulsars for this purpose. We highlight InPTA's observation strategies and analysis methods, which are relevant for a future PTA experiment with the more sensitive Square Kilometer Array (SKA) telescope. We show that the unique multi-sub-array multi-band wide-bandwidth frequency coverage of the InPTA provides Dispersion Measure estimates with unprecedented precision for PTA pulsars, e.g., ~ 2 x 10{-5} pc-cm{-3} for PSR J1909-3744. Configuring the SKA-low and SKA-mid as two and four sub-arrays respectively, it is shown that comparable precision is achievable, using observation strategies similar to those pursued by the InPTA, for a larger sample of 62 pulsars requiring about 26 and 7 hours per epoch for the SKA-mid and the SKA-low telescopes respectively. We also review the ongoing efforts to develop PTA-relevant general relativistic constructs that will be required to search for nanohertz gravitational waves from isolated super-massive black hole binary systems like blazar OJ 287. These efforts should be relevant to pursue persistent multi-messenger gravitational wave astronomy during the forthcoming era of the SKA telescope, the Thirty Meter Telescope, and the next-generation Event Horizon Telescope. △ Less

Submitted 13 July, 2022; originally announced July 2022.

Comments: Accepted for publication in Journal of Astronomy and Astrophysics for Special Issue on Indian Participation in the SKA (Editors : Abhirup Datta, Nirupam Roy, Preeti Kharb and Tirthankar Roy Choudhury)

The Indian Pulsar Timing Array: First data release

Authors: Pratik Tarafdar, Nobleson K., Prerna Rana, Jaikhomba Singha, M. A. Krishnakumar, Bhal Chandra Joshi, Avinash Kumar Paladi, Neel Kolhe, Neelam Dhanda Batra, Nikita Agarwal, Adarsh Bathula, Subhajit Dandapat, Shantanu Desai, Lankeswar Dey, Shinnosuke Hisano, Prathamesh Ingale, Ryo Kato, Divyansh Kharbanda, Tomonosuke Kikunaga, Piyush Marmat, B. Arul Pandian, T. Prabu, Aman Srivastava, Mayuresh Surnis, Sai Chaitanya Susarla , et al. (13 additional authors not shown)

Abstract: We present the pulse arrival times and high-precision dispersion measure estimates for 14 millisecond pulsars observed simultaneously in the 300-500 MHz and 1260-1460 MHz frequency bands using the upgraded Giant Metrewave Radio Telescope (uGMRT). The data spans over a baseline of 3.5 years (2018-2021), and is the first official data release made available by the Indian Pulsar Timing Array collabor… ▽ More We present the pulse arrival times and high-precision dispersion measure estimates for 14 millisecond pulsars observed simultaneously in the 300-500 MHz and 1260-1460 MHz frequency bands using the upgraded Giant Metrewave Radio Telescope (uGMRT). The data spans over a baseline of 3.5 years (2018-2021), and is the first official data release made available by the Indian Pulsar Timing Array collaboration. This data release presents a unique opportunity for investigating the interstellar medium effects at low radio frequencies and their impact on the timing precision of pulsar timing array experiments. In addition to the dispersion measure time series and pulse arrival times obtained using both narrowband and wideband timing techniques, we also present the dispersion measure structure function analysis for selected pulsars. Our ongoing investigations regarding the frequency dependence of dispersion measures have been discussed. Based on the preliminary analysis for five millisecond pulsars, we do not find any conclusive evidence of chromaticity in dispersion measures. Data from regular simultaneous two-frequency observations are presented for the first time in this work. This distinctive feature leads us to the highest precision dispersion measure estimates obtained so far for a subset of our sample. Simultaneous multi-band uGMRT observations in Band 3 and Band 5 are crucial for high-precision dispersion measure estimation and for the prospect of expanding the overall frequency coverage upon the combination of data from the various Pulsar Timing Array consortia in the near future. Parts of the data presented in this work are expected to be incorporated into the upcoming third data release of the International Pulsar Timing Array. △ Less

Submitted 25 October, 2022; v1 submitted 18 June, 2022; originally announced June 2022.

Comments: 23 pages, 21 figures, 3 tables. Published in PASA

Journal ref: Publications of the Astronomical Society of Australia, Volume 39, 2022, e053

Stability analysis of supermassive primordial stars: a new mass range for general relativistic instability supernovae

Authors: Chris Nagele, Hideyuki Umeda, Koh Takahashi, Takashi Yoshida, Kohsuke Sumiyoshi

Abstract: Observed supermassive black holes in the early universe have several proposed formation channels, in part because most of these channels are difficult to probe. One of the more promising channels, the direct collapse of a supermassive star, has several possible probes including the explosion of a helium-core supermassive star triggered by a general relativistic instability. We develop a straightfo… ▽ More Observed supermassive black holes in the early universe have several proposed formation channels, in part because most of these channels are difficult to probe. One of the more promising channels, the direct collapse of a supermassive star, has several possible probes including the explosion of a helium-core supermassive star triggered by a general relativistic instability. We develop a straightforward method for evaluating the general relativistic radial instability without simplifying assumptions and apply it to population III supermassive stars taken from a post Newtonian stellar evolution code. This method is more accurate than previous determinations and it finds that the instability occurs earlier in the evolutionary life of the star. Using the results of the stability analysis, we perform 1D general relativistic hydrodynamical simulations and we find two general relativistic instability supernovae fueled by alpha capture reactions as well as several lower mass pulsations, analogous to the puslational pair instability process. The mass range for the events (2.6-3.0 $\times 10^4$ ${\rm M}_\odot$) is lower than had been suggested by previous works (5.5 $\times 10^4$ ${\rm M}_\odot$) because the instability occurs earlier in the star's evolution. The explosion may be visible to, among others, JWST, while the discovery of the pulsations opens up additional possibilities for observation. △ Less

Submitted 4 September, 2022; v1 submitted 20 May, 2022; originally announced May 2022.

First High-speed Video Camera Observations of a Lightning Flash Associated with a Downward Terrestrial Gamma-ray Flash

Authors: R. U. Abbasi, M. M. F. Saba, J. W. Belz, P. R. Krehbiel, W. Rison, N. Kieu, D. R. da Silva, Dan Rodeheffer, M. A. Stanley, J. Remington, J. Mazich, R. LeVon, K. Smout, A. Petrizze, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, D. R. Bergman, S. A. Blake, I. Buckland, B. G. Cheon, M. Chikawa, T. Fujii , et al. (127 additional authors not shown)

Abstract: In this paper, we present the first high-speed video observation of a cloud-to-ground lightning flash and its associated downward-directed Terrestrial Gamma-ray Flash (TGF). The optical emission of the event was observed by a high-speed video camera running at 40,000 frames per second in conjunction with the Telescope Array Surface Detector, Lightning Mapping Array, interferometer, electric-field… ▽ More In this paper, we present the first high-speed video observation of a cloud-to-ground lightning flash and its associated downward-directed Terrestrial Gamma-ray Flash (TGF). The optical emission of the event was observed by a high-speed video camera running at 40,000 frames per second in conjunction with the Telescope Array Surface Detector, Lightning Mapping Array, interferometer, electric-field fast antenna, and the National Lightning Detection Network. The cloud-to-ground flash associated with the observed TGF was formed by a fast downward leader followed by a very intense return stroke peak current of -154 kA. The TGF occurred while the downward leader was below cloud base, and even when it was halfway in its propagation to ground. The suite of gamma-ray and lightning instruments, timing resolution, and source proximity offer us detailed information and therefore a unique look at the TGF phenomena. △ Less

Submitted 9 August, 2023; v1 submitted 10 May, 2022; originally announced May 2022.

Journal ref: Geophysical Research Letters, 50, e2023GL102958 (2023)

Linear stability of black holes with static scalar hair in full Horndeski theories: generic instabilities and surviving models

Authors: Masato Minamitsuji, Kazufumi Takahashi, Shinji Tsujikawa

Abstract: In full Horndeski theories, we show that the static and spherically symmetric black hole (BH) solutions with a static scalar field~$φ$ whose kinetic term~$X$ is nonvanishing on the BH horizon are generically prone to ghost/Laplacian instabilities. We then search for asymptotically Minkowski hairy BH solutions with a vanishing $X$ on the horizon free from ghost/Laplacian instabilities. We show that… ▽ More In full Horndeski theories, we show that the static and spherically symmetric black hole (BH) solutions with a static scalar field~$φ$ whose kinetic term~$X$ is nonvanishing on the BH horizon are generically prone to ghost/Laplacian instabilities. We then search for asymptotically Minkowski hairy BH solutions with a vanishing $X$ on the horizon free from ghost/Laplacian instabilities. We show that models with regular coupling functions of $φ$ and $X$ result in no-hair Schwarzschild BHs in general. On the other hand, the presence of a coupling between the scalar field and the Gauss-Bonnet (GB) term $R_{\rm GB}^2$, even with the coexistence of other regular coupling functions, leads to the realization of asymptotically Minkowski hairy BH solutions without ghost/Laplacian instabilities. Finally, we find that hairy BH solutions in power-law $F(R_{\rm GB}^2)$ gravity are plagued by ghost instabilities. These results imply that the GB coupling of the form $ξ(φ)R_{\rm GB}^2$ plays a prominent role for the existence of asymptotically Minkowski hairy BH solutions free from ghost/Laplacian instabilities. △ Less

Submitted 2 August, 2022; v1 submitted 28 April, 2022; originally announced April 2022.

Comments: 29 pages, no figures

Report number: YITP-22-41, WUAP-22-03

Journal ref: Phys. Rev. D 106, 044003 (2022)

Particle acceleration by counter-propagating circularly polarized Alfvén waves

Authors: S. Isayama, K. Takahashi, S. Matsukiyo, T. Sano

Abstract: Counterpropagating Alfvén waves are ubiquitously observed in many astrophysical environments, such as a star surface and a planetary foreshock. We discuss an efficient particle acceleration mechanism in two counterpropagating circularly polarized Alfvén waves. Phase transitions of particle behavior occur when wave amplitudes exceed two critical values. Above the critical amplitudes, the numerical… ▽ More Counterpropagating Alfvén waves are ubiquitously observed in many astrophysical environments, such as a star surface and a planetary foreshock. We discuss an efficient particle acceleration mechanism in two counterpropagating circularly polarized Alfvén waves. Phase transitions of particle behavior occur when wave amplitudes exceed two critical values. Above the critical amplitudes, the numerical simulation shows that any particles irreversibly gain relativistic energy within a short time regardless of their initial position and energy once the coherent waveform is formed. The accelerated particles have spatial coherence. Higher wave phase velocity requires smaller critical amplitudes, while the maximum attainable energy increases as the wavenumber and the frequency decrease. The results may be applicable in some astrophysical phenomena, as well as a future laboratory experiment using high-power lasers. △ Less

Submitted 5 April, 2023; v1 submitted 13 April, 2022; originally announced April 2022.

Effective Field Theory of Gravitating Continuum: Solids, Fluids, and Aether Unified

Authors: Katsuki Aoki, Mohammad Ali Gorji, Shinji Mukohyama, Kazufumi Takahashi

Abstract: We investigate the relativistic effective field theory (EFT) describing a non-dissipative gravitating continuum. In addition to ordinary continua, namely solids and fluids, we find an extraordinary more symmetric continuum, aether. In particular, the symmetry of the aether concludes that a homogeneous and isotropic state behaves like a cosmological constant. We formulate the EFT in the unitary/com… ▽ More We investigate the relativistic effective field theory (EFT) describing a non-dissipative gravitating continuum. In addition to ordinary continua, namely solids and fluids, we find an extraordinary more symmetric continuum, aether. In particular, the symmetry of the aether concludes that a homogeneous and isotropic state behaves like a cosmological constant. We formulate the EFT in the unitary/comoving gauge in which the dynamical degrees of freedom of the continuum (phonons) are eaten by the spacetime metric. This gauge choice, which is interpreted as the Lagrangian description in hydrodynamics, offers a neat geometrical understanding of continua. We examine a thread-based spacetime decomposition with respect to the four-velocity of the continuum which is different from the foliation-based Arnowitt-Deser-Misner one. Our thread-based decomposition respects the symmetries of the continua and, therefore, makes it possible to systematically find invariant building blocks of the EFT for each continuum even at higher orders in the derivative expansion. We also discuss the linear dynamics of the system and show that both gravitons and phonons acquire "masses" in a gravitating background. △ Less

Submitted 31 August, 2022; v1 submitted 13 April, 2022; originally announced April 2022.

Comments: 41 pages + appendices, 1 figure, 3 tables, added explanations, published version

Report number: YITP-22-38, IPMU22-0017

A Parkes "Murriyang" Search for Pulsars and Transients in the Large Magellanic Cloud

Authors: Shinnosuke Hisano, Fronefield Crawford, Victoria Bonidie, Md F. Alam, Keitaro Takahashi, Duncan R. Lorimer, Josh P. Ridley, Maura M. McLaughlin, Benetge B. P. Perera

Abstract: We report on the latest results of a Parkes multibeam survey for pulsars and dispersed radio bursts in the Large Magellanic Cloud (LMC). We conducted both periodicity and single-pulse searches at a much larger range of trial dispersion measures (DMs) than previously searched. We detected 229 single pulses with signal-to-noise ratio ($\rm{S/N}) > 7$ that were classified by the deep learning network… ▽ More We report on the latest results of a Parkes multibeam survey for pulsars and dispersed radio bursts in the Large Magellanic Cloud (LMC). We conducted both periodicity and single-pulse searches at a much larger range of trial dispersion measures (DMs) than previously searched. We detected 229 single pulses with signal-to-noise ratio ($\rm{S/N}) > 7$ that were classified by the deep learning network FETCH as being real (with $> 90$% likelihood), of which 9 were from the known giant-pulse-emitting pulsar PSR B0540$-$69. Two possibly repeating sources were detected with DMs suggesting that they lie within the LMC, but these require confirmation. Only 3 of the 220 unknown pulses had S/N greater than 8, and the DM distribution for these detected pulses follows an exponential falloff with increasing DM and does not show any excess of signals at DM values expected for the LMC. These features suggest that the detected pulses are not likely to be real, although they are visually compelling. We also report the discovery of a new pulsar (PSR J0556$-$67) in our periodicity search. This pulsar has a spin period of 791 ms, a DM of 71 cm$^{-3}$ pc, an estimated 1400 MHz flux density of $\sim 0.12$ mJy, and shows no evidence of binary motion. Future observations may be able to confirm whether any of the weak but promising candidates in our single pulse and periodicity searches, including our two possible repeaters, are real or not. △ Less

Submitted 23 February, 2022; v1 submitted 22 February, 2022; originally announced February 2022.

Comments: 17 pages, 11 figures

Stellar models of Betelgeuse constrained using observed surface conditions

Authors: Tianyin Luo, Hideyuki Umeda, Takashi Yoshida, Koh Takahashi

Abstract: We study stellar models for Betelgeuse using the HR diagram and surface abundances as observational constraints. Previous studies on Betelgeuse have not systematically investigated the surface abundances, but we believe they can be impacted by, and thus be used as an observational constraint for various parameters such as initial mass, rotation, and overshoot scheme. We investigate stellar models… ▽ More We study stellar models for Betelgeuse using the HR diagram and surface abundances as observational constraints. Previous studies on Betelgeuse have not systematically investigated the surface abundances, but we believe they can be impacted by, and thus be used as an observational constraint for various parameters such as initial mass, rotation, and overshoot scheme. We investigate stellar models with varying initial mass as they evolve past the main sequence, and we examine the red supergiant (RSG) properties in detail. For each mass, we vary the initial rotation up to $\sim$ 300 $km\:s^{-1}$, and test two different overshoot parameters. Overall, the acceptable initial mass range is 12 to 25 $M_\odot$, but for non-rotating models only, the range is decreased to 15 to 24 $M_\odot$. Also for rotating models, we find that $v/v_{\rm K} = 0.3$ is the upper limit for initial rotation, as more rapidly rotating models are unable to fit to Betelgeuse's surface abundances as an RSG. In addition, we report two possibilities for the current stage of evolution, core helium burning or core carbon burning and beyond. We find that certain 17 $M_\odot$ models could fit to both stages. Finally, we discuss the implications of our results in the context of merger scenarios which have been suggested as a mechanism to attain the observed surface velocity of Betelgeuse. △ Less

Submitted 7 February, 2022; originally announced February 2022.

ULF wave transmission across collisionless shocks: 2.5D local hybrid simulations

Authors: Primoz Kajdic, Yann Pfau-Kempf, Lucile Turc, Andrew P Dimmock, Minna Palmroth, Kazue Takahashi, Eemilia Kilpua, Jan Soucek, Naoko Takahashi, Luis Preisser, Xochitl Blanco-Cano, Domenico Trotta, David Burgess

Abstract: We study the interaction of upstream ultra-low frequency (ULF) waves with collisionless shocks by analyzing the outputs of eleven 2D local hybrid simulation runs. Our simulated shocks have Alfvénic Mach numbers between 4.29-7.42 and their $θ_{BN}$ angles are 15$^\circ$, 30$^\circ$, 45$^\circ$ and 50$^\circ$. The ULF wave foreshocks develop upstream of all of them. The wavelength and the amplitude… ▽ More We study the interaction of upstream ultra-low frequency (ULF) waves with collisionless shocks by analyzing the outputs of eleven 2D local hybrid simulation runs. Our simulated shocks have Alfvénic Mach numbers between 4.29-7.42 and their $θ_{BN}$ angles are 15$^\circ$, 30$^\circ$, 45$^\circ$ and 50$^\circ$. The ULF wave foreshocks develop upstream of all of them. The wavelength and the amplitude of the upstream waves exhibit a complex dependence on the shock's M$_A$ and $θ_{BN}$. The wavelength positively correlates with both parameters, with the dependence on $θ_{BN}$ being much stronger. The amplitude of the ULF waves is proportional to the product of the reflected beam velocity and density, which also depend on M$_A$ and $θ_{BN}$. The interaction of the ULF waves with the shock causes large-scale (several tens of upstream ion inertial lengths) shock rippling. The properties of the shock ripples are related to the ULF wave properties, namely thier wavelength and amplitude. In turn, the ripples have a large impact on the ULF wave transmission across the shock because they change local shock properties ($θ_{BN}$, strength), so that different sections of the same ULF wave front encounter shock with different characteristics. Downstream fluctuations do not resemble the upstream waves in terms the wavefront extension, orientation or their wavelength. However some features are conserved in the Fourier spectra of downstream compressive waves that present a bump or flattening at wavelengths approximately corresponding to those of the upstream ULF waves. In the transverse downstream spectra these features are weaker. △ Less

Submitted 25 January, 2022; originally announced January 2022.

Journal ref: Journal of Geophysical Research: Space Physics, 126, e2021JA029283

Linear stability of black holes in shift-symmetric Horndeski theories with a time-independent scalar field

Authors: Masato Minamitsuji, Kazufumi Takahashi, Shinji Tsujikawa

Abstract: We study linear perturbations about static and spherically symmetric black holes with a time-independent background scalar field in shift-symmetric Horndeski theories, whose Lagrangian is characterized by coupling functions depending only on the kinetic term of the scalar field $X$. We clarify conditions for the absence of ghosts and Laplacian instabilities along the radial and angular directions… ▽ More We study linear perturbations about static and spherically symmetric black holes with a time-independent background scalar field in shift-symmetric Horndeski theories, whose Lagrangian is characterized by coupling functions depending only on the kinetic term of the scalar field $X$. We clarify conditions for the absence of ghosts and Laplacian instabilities along the radial and angular directions in both odd- and even-parity perturbations. For reflection-symmetric theories described by a k-essence Lagrangian and a nonminimal derivative coupling with the Ricci scalar, we show that black holes endowed with nontrivial scalar hair are unstable around the horizon in general. This includes nonasymptotically flat black holes known to exist when the nonminimal derivative coupling to the Ricci scalar is a linear function of $X$. We also investigate several black hole solutions in nonreflection-symmetric theories. For cubic Galileons with the Einstein-Hilbert term, there exists a nonasymptotically flat hairy black hole with no ghosts/Laplacian instabilities. Also, for the scalar field linearly coupled to the Gauss-Bonnet term, asymptotically flat black hole solutions constructed perturbatively with respect to a small coupling are free of ghosts/Laplacian instabilities. △ Less

Submitted 2 May, 2022; v1 submitted 24 January, 2022; originally announced January 2022.

Comments: 23 pages, 1 figure, 1 table, published version

Report number: YITP-22-06, WUCG-22-01

Journal ref: Phys. Rev. D 105,104001(2022)

Low-frequency wideband timing of InPTA pulsars observed with the uGMRT

Authors: K Nobleson, Nikita Agarwal, Raghav Girgaonkar, Arul Pandian, Bhal Chandra Joshi, M A Krishnakumar, Abhimanyu Susobhanan, Shantanu Desai, T Prabu, Adarsh Bathula, Timothy T Pennucci, Sarmistha Banik, Manjari Bagchi, Neelam Dhanda Batra, Arpita Choudhary, Subhajit Dandapat, Lankeswar Dey, Yashwant Gupta, Shinnosuke Hisano, Ryo Kato, Divyansh Kharbanda, Tomonosuke Kikunaga, Neel Kolhe, Yogesh Maan, Piyush Marmat , et al. (7 additional authors not shown)

Abstract: High-precision measurements of the pulsar dispersion measure (DM) are possible using telescopes with low-frequency wideband receivers. We present an initial study of the application of the wideband timing technique, which can simultaneously measure the pulsar times of arrival (ToAs) and DMs, for a set of five pulsars observed with the upgraded Giant Metrewave Radio Telescope (uGMRT) as part of the… ▽ More High-precision measurements of the pulsar dispersion measure (DM) are possible using telescopes with low-frequency wideband receivers. We present an initial study of the application of the wideband timing technique, which can simultaneously measure the pulsar times of arrival (ToAs) and DMs, for a set of five pulsars observed with the upgraded Giant Metrewave Radio Telescope (uGMRT) as part of the Indian Pulsar Timing Array (InPTA) campaign. We have used the observations with the 300-500 MHz band of the uGMRT for this purpose. We obtain high precision in DM measurements with precisions of the order 10^{-6}cm^{-3}pc. The ToAs obtained have sub-μs precision and the root-mean-square of the post-fit ToA residuals are in the sub-μs range. We find that the uncertainties in the DMs and ToAs obtained with this wideband technique, applied to low-frequency data, are consistent with the results obtained with traditional pulsar timing techniques and comparable to high-frequency results from other PTAs. This work opens up an interesting possibility of using low-frequency wideband observations for precision pulsar timing and gravitational wave detection with similar precision as high-frequency observations used conventionally. △ Less

Submitted 23 February, 2022; v1 submitted 13 December, 2021; originally announced December 2021.

Comments: Accepted for publication in MNRAS

Wavelets and sparsity for Faraday tomography

Authors: Suchetha Cooray, Tsutomu T. Takeuchi, Shinsuke Ideguchi, Takuya Akahori, Yoshimitsu Miyashita, Keitaro Takahashi

Abstract: Faraday tomography through broadband polarimetry can provide crucial information on magnetized astronomical objects, such as quasars, galaxies, or galaxy clusters. However, the limited wavelength coverage of the instruments requires that we solve an ill-posed inverse problem when we want to obtain the Faraday dispersion function (FDF), a tomographic distribution of the magnetoionic media along the… ▽ More Faraday tomography through broadband polarimetry can provide crucial information on magnetized astronomical objects, such as quasars, galaxies, or galaxy clusters. However, the limited wavelength coverage of the instruments requires that we solve an ill-posed inverse problem when we want to obtain the Faraday dispersion function (FDF), a tomographic distribution of the magnetoionic media along the line of sight. This paper explores the use of wavelet transforms and the sparsity of the transformed FDFs in the form of wavelet shrinkage (WS) for finding better solutions to the inverse problem. We recently proposed the Constraining and Restoring iterative Algorithm for Faraday Tomography (CRAFT; Cooray et al. 2021), a new flexible algorithm that showed significant improvements over the popular methods such as Rotation Measure Synthesis. In this work, we introduce CRAFT+WS, a new version of CRAFT incorporating the ideas of wavelets and sparsity. CRAFT+WS exhibit significant improvements over the original CRAFT when tested for a complex FDF of realistic Galactic model. Reconstructions of FDFs demonstrate super-resolution in Faraday depth, uncovering previously unseen Faraday complexities in observations. The proposed approach will be necessary for effective cosmic magnetism studies using the Square Kilometre Array and its precursors. △ Less

Submitted 6 May, 2022; v1 submitted 2 December, 2021; originally announced December 2021.

Comments: 10 Pages, 3 Figures, Submitted to PASJ. Code will be made available upon publication

Invertible disformal transformations with higher derivatives

Authors: Kazufumi Takahashi, Hayato Motohashi, Masato Minamitsuji

Abstract: We consider a higher-derivative generalization of disformal transformations in $D$-dimensional spacetime and clarify the conditions under which they form a group with respect to the matrix product and the functional composition. These conditions allow us to systematically construct the inverse transformation in a fully covariant manner. Applying the invertible generalized disformal transformation… ▽ More We consider a higher-derivative generalization of disformal transformations in $D$-dimensional spacetime and clarify the conditions under which they form a group with respect to the matrix product and the functional composition. These conditions allow us to systematically construct the inverse transformation in a fully covariant manner. Applying the invertible generalized disformal transformation to known ghost-free scalar-tensor theories, we obtain a novel class of ghost-free scalar-tensor theories, whose action contains the third- or higher-order derivatives of the scalar field as well as nontrivial higher-order derivative couplings to the curvature tensor. △ Less

Submitted 4 January, 2022; v1 submitted 22 November, 2021; originally announced November 2021.

Comments: 14 pages

Report number: YITP-21-137

Journal ref: Phys. Rev. D 105, 024015 (2022)