Showing 1–50 of 247 results for author: Chakrabarty, D

The Radius of the High Mass Pulsar PSR J0740+6620 With 3.6 Years of NICER Data

Authors: Tuomo Salmi, Devarshi Choudhury, Yves Kini, Thomas E. Riley, Serena Vinciguerra, Anna L. Watts, Michael T. Wolff, Zaven Arzoumanian, Slavko Bogdanov, Deepto Chakrabarty, Keith Gendreau, Sebastien Guillot, Wynn C. G. Ho, Daniela Huppenkothen, Renee M. Ludlam, Sharon M. Morsink, Paul S. Ray

Abstract: We report an updated analysis of the radius, mass, and heated surface regions of the massive pulsar PSR J0740+6620 using NICER data from 2018 September 21 to 2022 April 21, a substantial increase in data set size compared to previous analyses. Using a tight mass prior from radio timing measurements and jointly modeling the new NICER data with XMM-Newton data, the inferred equatorial radius and gra… ▽ More We report an updated analysis of the radius, mass, and heated surface regions of the massive pulsar PSR J0740+6620 using NICER data from 2018 September 21 to 2022 April 21, a substantial increase in data set size compared to previous analyses. Using a tight mass prior from radio timing measurements and jointly modeling the new NICER data with XMM-Newton data, the inferred equatorial radius and gravitational mass are $12.49_{-0.88}^{+1.28}$ km and $2.073_{-0.069}^{+0.069}$ $M_\odot$ respectively, each reported as the posterior credible interval bounded by the $16\,\%$ and $84\,\%$ quantiles, with an estimated systematic error $\lesssim 0.1$ km. This result was obtained using the best computationally feasible sampler settings providing a strong radius lower limit but a slightly more uncertain radius upper limit. The inferred radius interval is also close to the $R=12.76_{-1.02}^{+1.49}$ km obtained by Dittmann et al. 2024, when they require the radius to be less than $16$ km as we do. The results continue to disfavor very soft equations of state for dense matter, with $R<11.15$ km for this high mass pulsar excluded at the $95\,\%$ probability. The results do not depend significantly on the assumed cross-calibration uncertainty between NICER and XMM-Newton. Using simulated data that resemble the actual observations, we also show that our pipeline is capable of recovering parameters for the inferred models reported in this paper. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: 17 pages, 9 figures (2 of which are figure sets), 2 tables, accepted for publication in ApJ

X-ray and Radio campaign of the Z-source GX 340+0: discovery of X-ray polarization and its implications

Authors: Yash Bhargava, Mason Ng, Liang Zhang, Arvind Balasubramanian, Thomas D. Russell, Aman Kaushik, Vishal Jadoliya, Swati Ravi, Sudip Bhattacharyya, Mayukh Pahari, Jeroen Homan, Herman L. Marshall, Deepto Chakrabarty, Francesco Carotenuto

Abstract: We present the discovery of X-ray polarization from the neutron star low-mass X-ray binary and Z-source, GX~340$+$0, using an Imaging X-ray Polarimetry Explorer (IXPE) observation in March 2024. Along with the IXPE observation, we conducted an extensive X-ray and radio monitoring campaign to ascertain the source properties during and around the IXPE observation. The source was within the horizonta… ▽ More We present the discovery of X-ray polarization from the neutron star low-mass X-ray binary and Z-source, GX~340$+$0, using an Imaging X-ray Polarimetry Explorer (IXPE) observation in March 2024. Along with the IXPE observation, we conducted an extensive X-ray and radio monitoring campaign to ascertain the source properties during and around the IXPE observation. The source was within the horizontal branch throughout the multiwavelength campaign. We measured a significant X-ray polarization in 2--8 keV with polarization degree (PD) = $4.02 \pm 0.35$% and polarization angle (PA) = $37.6 \pm 2.5^\circ$. The energy-dependent polarization indicates that in the 2-2.5 keV energy range, the PA is much lower, $\sim9\pm8^\circ$, while other energy bands are consistent with the PA found over 2.5--8 keV. The simultaneous AstroSat-IXPE spectro-polarimetric observations provide some evidence for independent polarization from various spectral components, hinting at a disparity in the PA from the accretion disk and the Comptonized emission, while suggesting an unpolarized emission from the blackbody component. Radio observations in the 0.7--9 GHz frequency range reveal a non-detection of radio emission in 0.7-1.5 GHz and a significant detection in 5.5--9 GHz, suggesting the presence of a spectral break in 1.5-5.5 GHz. Using ATCA observation we place upper limits on the radio polarization at $<$6% on the linear polarization and $<$4% on the circular polarization at 3$σ$ level. We discuss the origin of the X-ray polarization and its implications on the geometry of the spectral components. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: Submitted in ApJL, 4 figures, 3 tables

NICER Discovery that SRGA J144459.2-604207 is an Accreting Millisecond X-ray Pulsar

Authors: Mason Ng, Paul S. Ray, Andrea Sanna, Tod E. Strohmayer, Alessandro Papitto, Giulia Illiano, Arianna C. Albayati, Diego Altamirano, Tuğba Boztepe, Tolga Güver, Deepto Chakrabarty, Zaven Arzoumanian, D. J. K. Buisson, Elizabeth C. Ferrara, Keith C. Gendreau, Sebastien Guillot, Jeremy Hare, Gaurava K. Jaisawal, Christian Malacaria, Michael T. Wolff

Abstract: We present the discovery, with the Neutron Star Interior Composition Explorer (NICER), that SRGA J144459.2-604207 is a 447.9 Hz accreting millisecond X-ray pulsar (AMXP), which underwent a four-week long outburst starting on 2024 February 15. The AMXP resides in a 5.22 hr binary, orbiting a low-mass companion donor with $M_d>0.1M_\odot$. We report on the temporal and spectral properties from NICER… ▽ More We present the discovery, with the Neutron Star Interior Composition Explorer (NICER), that SRGA J144459.2-604207 is a 447.9 Hz accreting millisecond X-ray pulsar (AMXP), which underwent a four-week long outburst starting on 2024 February 15. The AMXP resides in a 5.22 hr binary, orbiting a low-mass companion donor with $M_d>0.1M_\odot$. We report on the temporal and spectral properties from NICER observations during the early days of the outburst, from 2024 February 21 through 2024 February 23, during which NICER also detected a type-I X-ray burst that exhibited a plateau lasting ~6 s. The spectra of the persistent emission were well described by an absorbed thermal blackbody and power-law model, with blackbody temperature $kT\approx0.9{\rm\,keV}$ and power-law photon index $Γ\approx1.9$. Time-resolved burst spectroscopy confirmed the thermonuclear nature of the burst, where an additional blackbody component reached a maximum temperature of nearly $kT\approx3{\rm\,keV}$ at the peak of the burst. We discuss the nature of the companion as well as the type-I X-ray burst. △ Less

Submitted 14 May, 2024; v1 submitted 30 April, 2024; originally announced May 2024.

Comments: 11 pages, 4 figures, 3 tables. Accepted by ApJL

Distribution and Recovery Phase of Geomagnetic Storms During Solar Cycles 23 and 24

Authors: Wageesh Mishra, Preity Sukla Sahani, Soumyaranjan Khuntia, Dibyendu Chakrabarty

Abstract: Coronal mass ejections (CMEs) and Stream Interaction Regions (SIRs) are the main drivers of intense geomagnetic storms. We study the distribution of geomagnetic storms associated with different drivers during solar cycles 23 and 24 (1996-2019). Although the annual occurrence rate of geomagnetic storms in both cycles tracks the sunspot cycle, the second peak in storm activity lags the second sunspo… ▽ More Coronal mass ejections (CMEs) and Stream Interaction Regions (SIRs) are the main drivers of intense geomagnetic storms. We study the distribution of geomagnetic storms associated with different drivers during solar cycles 23 and 24 (1996-2019). Although the annual occurrence rate of geomagnetic storms in both cycles tracks the sunspot cycle, the second peak in storm activity lags the second sunspot peak. SIRs contribute significantly to the second peak in storm numbers in both cycles, particularly for moderate to stronger-than-moderate storms. We note semiannual peaks in storm numbers much closer to equinoxes for moderate storms, and slightly shifted from equinoxes for intense and stronger-than-intense storms. We note a significant fraction of multiple-peak storms in both cycles due to isolated ICMEs/SIRs, while single-peak storms from multiple interacting drivers, suggesting a complex relationship between storm steps and their drivers. Our study focuses on investigating the recovery phases of geomagnetic storms and examining their dependencies on various storm parameters. Multiple-peak storms in both cycles have recovery phase duration strongly influenced by slow and fast decay phases with no correlation with the main phase buildup rate and Dst peak. However, the recovery phase in single-peak storms for both cycles depends to some extent on the main phase buildup rate and Dst peak, in addition to slow and fast decay phases. Future research should explore recovery phases of single and multiple-peak storms incorporating in-situ solar wind observations for a deeper understanding of storm evolution and decay processes. △ Less

Submitted 14 April, 2024; originally announced April 2024.

Comments: 12 pages, 5 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

The black hole low mass X-ray binary V404 Cygni is part of a wide hierarchical triple, and formed without a kick

Authors: Kevin B. Burdge, Kareem El-Badry, Erin Kara, Claude Canizares, Deepto Chakrabarty, Anna Frebel, Sarah C. Millholland, Saul Rappaport, Rob Simcoe, Andrew Vanderburg

Abstract: Evidence suggests that when compact objects such as black holes and neutron stars form, they may receive a ``natal kick,'' where the stellar remnant gains momentum. Observational evidence for neutron star kicks is substantial, yet limited for black hole natal kicks, and some proposed black hole formation scenarios result in very small kicks. Here, we report the discovery that the canonical black h… ▽ More Evidence suggests that when compact objects such as black holes and neutron stars form, they may receive a ``natal kick,'' where the stellar remnant gains momentum. Observational evidence for neutron star kicks is substantial, yet limited for black hole natal kicks, and some proposed black hole formation scenarios result in very small kicks. Here, we report the discovery that the canonical black hole low-mass X-ray binary V404 Cygni is part of a wide hierarchical triple with a tertiary companion at least 3500 astronomical units away from the inner binary. Given the orbital configuration, the black hole likely received a sub-5 kilometer per second kick to have avoided unbinding the tertiary. This discovery reveals that at least some black holes form with nearly no natal kick. Furthermore, the tertiary in this system lends credence to evolutionary models of low-mass X-ray binaries involving a hierarchical triple structure. Remarkably, the tertiary is evolved, indicating that the system formed 3-5 billion years ago, and that the black hole has removed at least half a solar mass of matter from its evolved secondary companion. During the event in which the black hole formed, it is likely that at least half of the mass of the black hole progenitor collapsed into the black hole; it may even have undergone a complete implosion, enabling the tertiary to remain loosely bound. △ Less

Submitted 4 April, 2024; originally announced April 2024.

Comments: 36 pages, 10 figures, 1 table

Prolonged and Extremely Non-radial Solar Wind Flows

Authors: Susanta Kumar Bisoi, Diptiranjan Rout, P. Janardhan, K. Fujiki, Dibyendu Chakrabarty, Karan Sahu

Abstract: We present a study of three highly non-radial solar wind events when the azimuthal solar wind flow angle exceeds > 6 degrees for one day or more. None of the events are associated with coronal mass ejections and co-rotating interaction regions observed at 1 AU. For all events, the solar wind outflows at 1 AU have low solar wind velocity and solar wind density. Based on the significant increase in… ▽ More We present a study of three highly non-radial solar wind events when the azimuthal solar wind flow angle exceeds > 6 degrees for one day or more. None of the events are associated with coronal mass ejections and co-rotating interaction regions observed at 1 AU. For all events, the solar wind outflows at 1 AU have low solar wind velocity and solar wind density. Based on the significant increase in the Oxygen charge state ratio of O7+/O6+ at 1 AU for all of the events, we have traced them back to the Sun and found that their source regions originated in an active region and coronal hole (AR-CH) pairs mainly located at the central meridian. Further, examining the dynamical evolutions in their source regions using both the Extreme ultra-violet Imaging Telescope and Michelson Doppler Imager, it is found that the changes taking place in AR-CH boundaries eventually disturbed the stable CH configurations, resulting in a reduction of the CH area and finally its disappearance, leaving only with the AR. Our study provides a possible explanation to discuss the origin of the prolonged and highly non-radial solar wind flows. △ Less

Submitted 19 January, 2024; originally announced January 2024.

Comments: Bulletin de la Société Royale des Sciences de Liège (BSRSL), 93 (2), 1-13 (in press)

SWASTi-CME: A physics-based model to study CME evolution and its interaction with Solar Wind

Authors: Prateek Mayank, Bhargav Vaidya, Wageesh Mishra, D. Chakrabarty

Abstract: Coronal mass ejections (CMEs) are primary drivers of space weather and studying their evolution in the inner heliosphere is vital to prepare for a timely response. Solar wind streams, acting as background, influence their propagation in the heliosphere and associated geomagnetic storm activity. This study introduces SWASTi-CME, a newly developed MHD-based CME model integrated into the Space Weathe… ▽ More Coronal mass ejections (CMEs) are primary drivers of space weather and studying their evolution in the inner heliosphere is vital to prepare for a timely response. Solar wind streams, acting as background, influence their propagation in the heliosphere and associated geomagnetic storm activity. This study introduces SWASTi-CME, a newly developed MHD-based CME model integrated into the Space Weather Adaptive SimulaTion (SWASTi) framework. It incorporates a non-magnetized elliptic cone and a magnetized flux rope CME model. To validate the model's performance with in-situ observation at L1, two Carrington rotations were chosen: one during solar maxima with multiple CMEs, and one during solar minima with a single CME. The study also presents a quantitative analysis of CME-solar wind interaction using this model. To account for ambient solar wind effects, two scenarios of different complexity in solar wind conditions were established. The results indicate that ambient conditions can significantly impact some of the CME properties in the inner heliosphere. We found that the drag force on the CME front exhibits a variable nature, resulting in asymmetric deformation of the CME leading edge. Additionally, the study reveals that the impact on the distribution of CME internal pressure primarily occurs during the initial stage, while the CME density distribution is affected throughout its propagation. Moreover, regardless of the ambient conditions, it was observed that after a certain propagation time (t), the CME volume follows a non-fractal power-law expansion ($\propto t^{3.03-3.33}$) due to the attainment of a balanced state with ambient. △ Less

Submitted 27 October, 2023; originally announced October 2023.

Comments: Accepted for publication in ApJS

X-ray and Radio Monitoring of the Neutron Star Low Mass X-ray Binary 1A 1744-361: Quasi Periodic Oscillations, Transient Ejections, and a Disk Atmosphere

Authors: Mason Ng, Andrew K. Hughes, Jeroen Homan, Jon M. Miller, Sean N. Pike, Diego Altamirano, Peter Bult, Deepto Chakrabarty, D. J. K. Buisson, Benjamin M. Coughenour, Rob Fender, Sebastien Guillot, Tolga Güver, Gaurava K. Jaisawal, Amruta D. Jaodand, Christian Malacaria, James C. A. Miller-Jones, Andrea Sanna, Gregory R. Sivakoff, Tod E. Strohmayer, John A. Tomsick, Jakob van den Eijnden

Abstract: We report on X-ray (NICER/NuSTAR/MAXI/Swift) and radio (MeerKAT) timing and spectroscopic analysis from a three-month monitoring campaign in 2022 of a high-intensity outburst of the dipping neutron star low-mass X-ray binary 1A 1744-361. The 0.5-6.8 keV NICER X-ray hardness-intensity and color-color diagrams of the observations throughout the outburst suggests that 1A 1744-361 spent most of its ou… ▽ More We report on X-ray (NICER/NuSTAR/MAXI/Swift) and radio (MeerKAT) timing and spectroscopic analysis from a three-month monitoring campaign in 2022 of a high-intensity outburst of the dipping neutron star low-mass X-ray binary 1A 1744-361. The 0.5-6.8 keV NICER X-ray hardness-intensity and color-color diagrams of the observations throughout the outburst suggests that 1A 1744-361 spent most of its outburst in an atoll-state, but we show that the source exhibited Z-state-like properties at the peak of the outburst, similar to a small sample of other atoll-state sources. A timing analysis with NICER data revealed several instances of an $\approx8$ Hz quasi-periodic oscillation (QPO; fractional rms amplitudes of ~5%) around the peak of the outburst, the first from this source, which we connect to the normal branch QPOs (NBOs) seen in the Z-state. Our observations of 1A 1744-361 are fully consistent with the idea of the mass accretion rate being the main distinguishing parameter between atoll- and Z-states. Radio monitoring data by MeerKAT suggests that the source was at its radio-brightest during the outburst peak, and that the source transitioned from the 'island' spectral state to the 'banana' state within ~3 days of the outburst onset, launching transient jet ejecta. The observations present the strongest evidence for radio flaring, including jet ejecta, during the island-to-banana spectral state transition at low accretion rates (atoll-state). The source also exhibited Fe XXV, Fe XXVI K$α$, and K$β$ X-ray absorption lines, whose origins likely lie in an accretion disk atmosphere. △ Less

Submitted 30 April, 2024; v1 submitted 2 October, 2023; originally announced October 2023.

Comments: 30 pages, 13 figures, and 8 tables. Accepted by ApJ (before proofs)

Suprathermal population associated with stream interaction regions observed by STEREO-A: New insights

Authors: Bijoy Dalal, Dibyendu Chakrabarty, Nandita Srivastava, Aveek Sarkar

Abstract: Stream interaction regions (SIRs) are often thought to be responsible for the generation of suprathermal population in the interplanetary medium. Despite the source being same, wide variations in spectral indices of suprathermal populations are observed at 1 au during SIRs. This poses significant uncertainty in understanding the generation of suprathermal ion populations by SIRs and indicates inte… ▽ More Stream interaction regions (SIRs) are often thought to be responsible for the generation of suprathermal population in the interplanetary medium. Despite the source being same, wide variations in spectral indices of suprathermal populations are observed at 1 au during SIRs. This poses significant uncertainty in understanding the generation of suprathermal ion populations by SIRs and indicates interplay of multiple source mechanisms. In the present work, by analyzing variations in suprathermal 4He, O, and Fe for 20 SIR events recorded by STEREO-A during 2007 - 2014, we find that the spectral indices of these elements vary in the range of 2.06-4.08, 1.85-4.56, and 2.11-4.04 respectively for 19 events. However, in one special case, all the three suprathermal elements show nearly identical (1.5) spectral indices. We offer possible mechanisms, which can cause significant variations in the spectral indices of suprathermal particles. More importantly, we show the possible role of merging and/or contraction of small-scale magnetic islands near 1 au in producing nearly identical spectral indices for three different elements with different first ionization potential and mass-to-charge ratio. The occurrence of these magnetic islands near 1 au also supports the minimal modulation in spectral indices of these particles. The role of a possible solar flare in generating these magnetic islands near the heliospheric current sheet is also suggested. △ Less

Submitted 19 September, 2023; originally announced September 2023.

NICER observations of thermonuclear bursts from 4U 1728-34: Detection of oscillations prior to the onset of two bursts

Authors: Funda Bostanci, Tugba Boztepe, Tolga Guver, Tod E. Strohmayer, Yuri Cavecchi, Ersin Gogus, Diego Altamirano, Peter Bult, Deepto Chakrabarty, Sebastien Guillot, Gaurava K. Jaisawal, Christian Malacaria, Giulio C. Mancuso, Andrea Sanna, Jean H. Swank

Abstract: We present temporal and time-resolved spectral analyses of all the thermonuclear X-ray bursts observed from the neutron star low-mass X-ray binary (LMXB) 4U 1728-34 with NICER from June 2017 to September 2019. In total, we detected 11 X-ray bursts from the source and performed time-resolved spectroscopy. Unlike some of the earlier results for other bursting sources from NICER, our spectral results… ▽ More We present temporal and time-resolved spectral analyses of all the thermonuclear X-ray bursts observed from the neutron star low-mass X-ray binary (LMXB) 4U 1728-34 with NICER from June 2017 to September 2019. In total, we detected 11 X-ray bursts from the source and performed time-resolved spectroscopy. Unlike some of the earlier results for other bursting sources from NICER, our spectral results indicate that the use of a scaling factor for the persistent emission is not statistically necessary. This is primarily a result of the strong interstellar absorption in the line of sight towards 4U 1728-34, which causes the count rates to be significantly lower at low energies. We also searched for burst oscillations and detected modulations in six different bursts at around the previously known burst oscillation frequency of 363 Hz. Finally, we report the detection of oscillations prior to two bursts at 356 and 359 Hz, respectively. This is the first time in the literature where burst oscillations are detected before the rapid rise in X-ray flux, from any known burster. These oscillations disappear as soon as the burst rise starts and occur at a somewhat lower frequency than the oscillations we detect during the bursts. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: Accepted for publication in the Astrophysical Journal

An updated mass-radius analysis of the 2017-2018 NICER data set of PSR J0030+0451

Authors: Serena Vinciguerra, Tuomo Salmi, Anna L. Watts, Devarshi Choudhury, Thomas E. Riley, Paul S. Ray, Slavko Bogdanov, Yves Kini, Sebastien Guillot, Deepto Chakrabarty, Wynn C. G. Ho, Daniela Huppenkothen, Sharon M. Morsink, Zorawar Wadiasingh

Abstract: In 2019 the NICER collaboration published the first mass and radius inferred for PSR J0030+0451, thanks to NICER observations, and consequent constraints on the equation of state characterising dense matter. Two independent analyses found a mass of $\sim 1.3-1.4\,\mathrm{M_\odot}$ and a radius of $\sim 13\,$km. They also both found that the hot spots were all located on the same hemisphere, opposi… ▽ More In 2019 the NICER collaboration published the first mass and radius inferred for PSR J0030+0451, thanks to NICER observations, and consequent constraints on the equation of state characterising dense matter. Two independent analyses found a mass of $\sim 1.3-1.4\,\mathrm{M_\odot}$ and a radius of $\sim 13\,$km. They also both found that the hot spots were all located on the same hemisphere, opposite to the observer, and that at least one of them had a significantly elongated shape. Here we reanalyse, in greater detail, the same NICER data set, incorporating the effects of an updated NICER response matrix and using an upgraded analysis framework. We expand the adopted models and jointly analyse also XMM-Newton data, which enables us to better constrain the fraction of observed counts coming from PSR J0030+0451. Adopting the same models used in previous publications, we find consistent results, although with more stringent inference requirements. We also find a multi-modal structure in the posterior surface. This becomes crucial when XMM-Newton data is accounted for. Including the corresponding constraints disfavors the main solutions found previously, in favor of the new and more complex models. These have inferred masses and radii of $\sim [1.4 \mathrm{M_\odot}, 11.5$ km] and $\sim [1.7 \mathrm{M_\odot}, 14.5$ km], depending on the assumed model. They display configurations that do not require the two hot spots generating the observed X-rays to be on the same hemisphere, nor to show very elongated features, and point instead to the presence of temperature gradients and the need to account for them. △ Less

Submitted 18 August, 2023; originally announced August 2023.

Comments: 35 pages, 11 figures, 7 tables

Journal ref: ApJ, 961, 62 (2024)

New insights on the behaviour of solar wind protons and alphas in the Stream Interaction Region in solar cycle 23 and 24

Authors: Yogesh, D. Chakrabarty, Nandita Srivastava

Abstract: Although the enhancements in the alpha-proton ratio in the solar wind (expressed as $A_{He} = N_a/N_p*100$) in the Interplanetary Coronal Mass Ejections (ICMEs) have been studied in the past, $A_{He}$ enhancements at the stream interface region received very little attention so far. In this letter, by extensively analyzing the stream interaction region (SIR) events observed in solar cycle 23 and 2… ▽ More Although the enhancements in the alpha-proton ratio in the solar wind (expressed as $A_{He} = N_a/N_p*100$) in the Interplanetary Coronal Mass Ejections (ICMEs) have been studied in the past, $A_{He}$ enhancements at the stream interface region received very little attention so far. In this letter, by extensively analyzing the stream interaction region (SIR) events observed in solar cycle 23 and 24, we show that the stream interface of alphas starts separating out from that of protons from the minimum of solar cycle 23. The population of alpha particles are enhanced compared to protons at higher angles between bulk velocity and local magnetic field (henceforth, bulk velocity angle) in the fast wind region of SIRs if the background solar wind is taken as reference. The analysis of differential velocities between alphas and protons also reveals that the faster alpha particles accumulate near the fast wind side of the stream interface region leading to enhancement of $A_{He}$. The investigation brings out, for the first time, the salient changes in $A_{He}$ in SIRs for the two solar cycles and highlight the importances of bulk velocity angle and differential velocity in the fast wind region for the changes in $A_{He}$ in SIRs. △ Less

Submitted 8 August, 2023; v1 submitted 1 April, 2023; originally announced April 2023.

Orbital decay in an accreting and eclipsing 13.7 minute orbital period binary with a luminous donor

Authors: Kevin B. Burdge, Kareem El-Badry, Saul Rappaport, Tin Long Sunny Wong, Evan B. Bauer, Lars Bildsten, Ilaria Caiazzo, Deepto Chakrabarty, Emma Chickles, Matthew J. Graham, Erin Kara, S. R. Kulkarni, Thomas R. Marsh, Melania Nynka, Thomas A. Prince, Robert A. Simcoe, Jan van Roestel, Zach Vanderbosch, Eric C. Bellm, Richard G. Dekany, Andrew J. Drake, George Helou, Frank J. Masci, Jennifer Milburn, Reed Riddle , et al. (2 additional authors not shown)

Abstract: We report the discovery of ZTF J0127+5258, a compact mass-transferring binary with an orbital period of 13.7 minutes. The system contains a white dwarf accretor, which likely originated as a post-common envelope carbon-oxygen (CO) white dwarf, and a warm donor ($T_{\rm eff,\,donor}= 16,400\pm1000\,\rm K$). The donor probably formed during a common envelope phase between the CO white dwarf and an e… ▽ More We report the discovery of ZTF J0127+5258, a compact mass-transferring binary with an orbital period of 13.7 minutes. The system contains a white dwarf accretor, which likely originated as a post-common envelope carbon-oxygen (CO) white dwarf, and a warm donor ($T_{\rm eff,\,donor}= 16,400\pm1000\,\rm K$). The donor probably formed during a common envelope phase between the CO white dwarf and an evolving giant which left behind a helium star or helium white dwarf in a close orbit with the CO white dwarf. We measure gravitational wave-driven orbital inspiral with $\sim 35σ$ significance, which yields a joint constraint on the component masses and mass transfer rate. While the accretion disk in the system is dominated by ionized helium emission, the donor exhibits a mixture of hydrogen and helium absorption lines. Phase-resolved spectroscopy yields a donor radial-velocity semi-amplitude of $771\pm27\,\rm km\, s^{-1}$, and high-speed photometry reveals that the system is eclipsing. We detect a {\it Chandra} X-ray counterpart with $L_{X}\sim 3\times 10^{31}\,\rm erg\,s^{-1}$. Depending on the mass-transfer rate, the system will likely evolve into either a stably mass-transferring helium CV, merge to become an R Crb star, or explode as a Type Ia supernova in the next million years. We predict that the Laser Space Interferometer Antenna (LISA) will detect the source with a signal-to-noise ratio of $24\pm6$ after 4 years of observations. The system is the first \emph{LISA}-loud mass-transferring binary with an intrinsically luminous donor, a class of sources that provide the opportunity to leverage the synergy between optical and infrared time domain surveys, X-ray facilities, and gravitational-wave observatories to probe general relativity, accretion physics, and binary evolution. △ Less

Submitted 23 March, 2023; originally announced March 2023.

Comments: 13 pages, 7 figures, 2 tables, submitted to ApJL

Clumpy wind studies and the non-detection of cyclotron line in OAO 1657-415

Authors: Pragati Pradhan, Carlo Ferrigno, Biswajit Paul, Enrico Bozzo, Ileyk El Mellah, David P. Huenemoerder, James F. Steiner, Victoria Grinberg, Felix Furst, Chandreyee Maitra, Patrizia Romano, Peter Kretschma, Jamie Kennea, Deepto Chakrabarty

Abstract: Winds of massive stars are suspected to be inhomogeneous (or clumpy), which biases the measures of their mass loss rates. In High Mass X-ray Binaries (HMXBs), the compact object can be used as an orbiting X-ray point source to probe the wind and constrain its clumpiness. We perform spectro-timing analysis of the HMXB OAO 1657-415 with non-simultaneous NuSTAR and NICER observations. We compute the… ▽ More Winds of massive stars are suspected to be inhomogeneous (or clumpy), which biases the measures of their mass loss rates. In High Mass X-ray Binaries (HMXBs), the compact object can be used as an orbiting X-ray point source to probe the wind and constrain its clumpiness. We perform spectro-timing analysis of the HMXB OAO 1657-415 with non-simultaneous NuSTAR and NICER observations. We compute the hardness ratio from the energy-resolved light curves, and using an adaptive rebinning technique, we thus select appropriate time segments to search for rapid spectral variations on timescales of a few hundreds to thousands of seconds. Column density and intensity of Iron K$α$ line were strongly correlated, and the recorded spectral variations were consistent with accretion from a clumpy wind. We also illustrate a novel framework to measure clump sizes, masses in HMXBs more accurately based on absorption measurements and orbital parameters of the source. We then discuss the limitations posed by current X-ray spacecrafts in such measurements and present prospects with future X-ray missions. We find that the source pulse profiles show a moderate dependence on energy. We identify a previously undetected dip in the pulse profile visible throughout the NuSTAR observation near spin phase 0.15 possibly caused by intrinsic changes in accretion geometry close to the neutron star. We do not find any evidence for the debated cyclotron line at $\sim$ 36\,keV in the time-averaged or the phase-resolved spectra with NuSTAR. △ Less

Submitted 21 February, 2023; originally announced February 2023.

Comments: 32 pages, 11 figures in main text, 7 figures in Appendix, Accepted for publication in ApJ

Global asymmetry in $Δ$X variations during the 06 April 2000 geomagnetic storm: Relative roles of IMF Bz and By

Authors: Sumanjit Chakraborty, D. Chakrabarty

Abstract: This investigation is directed to understand the asymmetry in $Δ$X variations caused due to the relative roles played by IMF Bz and IMF By in a particular interval (22:22 - 22:55 UT), during the main phase of a strong geomagnetic storm event of April 06, 2000 (Ap = 236). Two pairs of antipodal stations, being part of the SuperMAG network, are considered here. Ionospheric convection maps from Super… ▽ More This investigation is directed to understand the asymmetry in $Δ$X variations caused due to the relative roles played by IMF Bz and IMF By in a particular interval (22:22 - 22:55 UT), during the main phase of a strong geomagnetic storm event of April 06, 2000 (Ap = 236). Two pairs of antipodal stations, being part of the SuperMAG network, are considered here. Ionospheric convection maps from SuperDARN network are used to understand spatio-temporal evolution of the DP2 ionospheric convection patterns over high-latitudes. The two-dimensional maps of equivalent currents are used to show signatures of global DP2 currents associated with the interplay effect between the two IMF components. Observations show increases in the difference in $Δ$X variations between nearly antipodal stations from the Japanese-European/African sector with respect to the same between the nearly antipodal stations from the Pacific/American-Indian sector. This asymmetry is observed during the period when the absolute magnitude of IMF By is larger than that of IMF Bz resulting in a significant and conspicuous enhancement in IMF |By/Bz|. It is suggested that the distortions in DP2 cells and associated rotation of electrodynamic day-night divider, bring one pair of stations under the same DP2 cell and one station of the other pair under a different DP2 cell and throat flow region leading to the asymmetry in $Δ$X variations between the antipodal stations. Therefore, the work highlights the importance of the interplay between IMF Bz and IMF By in determining the ionospheric impact over low latitudes during strong geomagnetic conditions. △ Less

Submitted 24 January, 2023; originally announced January 2023.

Comments: Accepted for publication in JGR Space Physics on 04 January 2023

Journal ref: Journal of Geophysical Research: Space Physics (2023)

The Aditya-L1 mission of ISRO

Authors: Durgesh Tripathi, D. Chakrabarty, A. Nandi, B. Raghvendra Prasad, A. N. Ramaprakash, Nigar Shaji, K. Sankarasubramanian, R. Satheesh Thampi, V. K. Yadav

Abstract: The Aditya-L1 is the first space-based solar observatory of the Indian Space Research Organization (ISRO). The spacecraft will carry seven payloads providing uninterrupted observations of the Sun from the first Lagrangian point. Aditya-L1 comprises four remote sensing instruments, {\it viz.} a coronagraph observing in visible and infrared, a full disk imager in Near Ultra-Violet (NUV), and two ful… ▽ More The Aditya-L1 is the first space-based solar observatory of the Indian Space Research Organization (ISRO). The spacecraft will carry seven payloads providing uninterrupted observations of the Sun from the first Lagrangian point. Aditya-L1 comprises four remote sensing instruments, {\it viz.} a coronagraph observing in visible and infrared, a full disk imager in Near Ultra-Violet (NUV), and two full-sun integrated spectrometers in soft X-ray and hard X-ray. In addition, there are three instruments for in-situ measurements, including a magnetometer, to study the magnetic field variations during energetic events. Aditya-L1 is truly a mission for multi-messenger solar astronomy from space that will provide comprehensive observations of the Sun across the electromagnetic spectrum and in-situ measurements in a broad range of energy, including magnetic field measurements at L1. △ Less

Submitted 30 December, 2022; v1 submitted 26 December, 2022; originally announced December 2022.

Comments: 10 pages, 6 figures, Submitted to the Proceedings of IAU 372: The Era of Multi-Messenger Solar Physics"

Spectral Evolution of Ultraluminous X-ray Pulsar NGC 300 ULX-1

Authors: Mason Ng, Ronald A. Remillard, James F. Steiner, Deepto Chakrabarty, Dheeraj R. Pasham

Abstract: We report on results from a one-year soft X-ray observing campaign of the ultraluminous X-ray pulsar NGC 300 ULX-1 by the Neutron star Interior Composition Explorer (NICER) during 2018--2019. Our analysis also made use of data from Swift/XRT and XMM-Newton in order to model and remove contamination from the nearby eclipsing X-ray binary NGC 300 X-1. We constructed and fitted a series of 5-day aver… ▽ More We report on results from a one-year soft X-ray observing campaign of the ultraluminous X-ray pulsar NGC 300 ULX-1 by the Neutron star Interior Composition Explorer (NICER) during 2018--2019. Our analysis also made use of data from Swift/XRT and XMM-Newton in order to model and remove contamination from the nearby eclipsing X-ray binary NGC 300 X-1. We constructed and fitted a series of 5-day averaged NICER spectra of NGC 300 ULX-1 in the 0.4--4.0 keV range to evaluate the long-term spectral evolution of the source, and found that an absorbed power-law model provided the best fit overall. Over the course of our observations, the source flux (0.4--4.0 keV; absorbed) dimmed from $2\times10^{-12}$ to below $10^{-13}{\rm\,erg\,s^{-1}\,cm^{-2}}$ and the spectrum softened, with the photon index going from $Γ\approx1.6$ to $Γ\approx2.6$. We interpret the spectral softening as reprocessed emission from the accretion disk edge coming into view while the pulsar was obscured by the possibly precessing disk. Some spectral fits were significantly improved by the inclusion of a disk blackbody component, and we surmise that this could be due to the pulsar emerging in between obscuration episodes by partial covering absorbers. We posit that we observed a low-flux state of the system (due to line-of-sight absorption) punctuated by the occasional appearance of the pulsar, indicating short-term source variability nested in longer-term accretion disk precession timescales. △ Less

Submitted 11 October, 2022; originally announced October 2022.

Comments: 13 pages, including 7 figures and 2 tables. Accepted for publication in the Astrophysical Journal

Observations of 4U 1626-67 with the Imaging X-ray Polarimetry Explorer

Authors: Herman L. Marshall, Mason Ng, Daniele Rogantini, Jeremy Heyl, Sergey S. Tsygankov, Juri Poutanen, Enrico Costa, Silvia Zane, Christian Malacaria, Ivan Agudo, Lucio A. Antonelli, Matteo Bachetti, Luca Baldini, Wayne H. Baumgartner, Ronaldo Bellazzini, Stefano Bianchi, Stephen D. Bongiorno, Raffaella Bonino, Alessandro Brez, Niccolo Bucciantini, Fiamma Capitanio, Simone Castellano, Elisabetta Cavazzuti, Stefano Ciprini, Alessandra De Rosa , et al. (69 additional authors not shown)

Abstract: We present measurements of the polarization of X-rays in the 2-8 keV band from the pulsar in the ultracompact low mass X-ray binary 4U1626-67 using data from the Imaging X-ray Polarimetry Explorer (IXPE). The 7.66 s pulsations were clearly detected throughout the IXPE observations as well as in the NICER soft X-ray observations, which we use as the basis for our timing analysis and to constrain th… ▽ More We present measurements of the polarization of X-rays in the 2-8 keV band from the pulsar in the ultracompact low mass X-ray binary 4U1626-67 using data from the Imaging X-ray Polarimetry Explorer (IXPE). The 7.66 s pulsations were clearly detected throughout the IXPE observations as well as in the NICER soft X-ray observations, which we use as the basis for our timing analysis and to constrain the spectral shape over 0.4-10 keV energy band. Chandra HETGS high-resolution X-ray spectra were also obtained near the times of the IXPE observations for firm spectral modeling. We find an upper limit on the pulse-averaged linear polarization of <4% (at 95% confidence). Similarly, there was no significant detection of polarized flux in pulse phase intervals when subdividing the bandpass by energy. However, spectropolarimetric modeling over the full bandpass in pulse phase intervals provide a marginal detection of polarization of the power-law spectral component at the 4.8 +/- 2.3% level (90% confidence). We discuss the implications concerning the accretion geometry onto the pulsar, favoring two-component models of the pulsed emission. △ Less

Submitted 6 October, 2022; originally announced October 2022.

Comments: 19 pages, 7 figures, 7 tables; accepted for publication in the Astrophysical Journal

A dense $\mathbf{0.1 M_{\rm \odot}}$ star in a 51-minute orbital period eclipsing binary

Authors: Kevin B. Burdge, Kareem El-Badry, Thomas R. Marsh, Saul Rappaport, Warren R. Brown, Ilaria Caiazzo, Deepto Chakrabarty, V. S. Dhillon, Jim Fuller, Boris T. Gänsicke, Matthew J. Graham, Erin Kara, S. R. Kulkarni, S. P. Littlefair, Przemek Mróz, Pablo Rodríguez-Gil, Jan van Roestel, Robert A. Simcoe, Eric C. Bellm, Andrew J. Drake, Richard G. Dekany, Steven L. Groom, Russ R. Laher, Frank J. Masci, Reed Riddle , et al. (2 additional authors not shown)

Abstract: In over a thousand known cataclysmic variables (CVs), where a white dwarf is accreting from a hydrogen-rich star, only a dozen have orbital periods below 75 minutes. One way to achieve these short periods requires the donor star to have undergone substantial nuclear evolution prior to interacting with the white dwarf, and it is expected that these objects will transition to helium accretion. These… ▽ More In over a thousand known cataclysmic variables (CVs), where a white dwarf is accreting from a hydrogen-rich star, only a dozen have orbital periods below 75 minutes. One way to achieve these short periods requires the donor star to have undergone substantial nuclear evolution prior to interacting with the white dwarf, and it is expected that these objects will transition to helium accretion. These transitional CVs have been proposed as progenitors of helium CVs. However, no known transitional CV is expected to reach an orbital period short enough to account for most of the helium CV population, leaving the role of this evolutionary pathway unclear. Here we report observations of ZTF J1813+4251, a 51-minute orbital period, fully eclipsing binary system consisting of a star with a temperature comparable to that of the Sun but a density 100 times greater due to its helium-rich composition, accreting onto a white dwarf. Phase-resolved spectra, multi-band light curves and the broadband spectral energy distribution allow us to obtain precise and robust constraints on the masses, radii and temperatures of both components. Evolutionary modeling shows that ZTF J1813+4251 is destined to become a helium CV binary, reaching an orbital period under 20 minutes, rendering ZTF J1813+4251 a previously missing link between helium CV binaries and hydrogen-rich CVs. △ Less

Submitted 4 October, 2022; originally announced October 2022.

Comments: 48 Pages, 12 figures, 2 tables, Published online by Nature on Oct 5, 2022

Journal ref: Nature 610 467-471 (2022)

MAXI J1957+032: a new accreting millisecond X-ray pulsar in an ultra-compact binary

Authors: A. Sanna, P. Bult, M. NG, P. S. Ray, G. K. Jaisawal, L. Burderi, T. Di Salvo, A. Riggio, D. Altamirano, T. E. Strohmayer, A. Manca, K. C. Gendreau, D. Chakrabarty, W. Iwakiri, R. Iaria

Abstract: The detection of coherent X-ray pulsations at ~314 Hz (3.2 ms) classifies MAXI J1957+032 as a fast-rotating, accreting neutron star. We present the temporal and spectral analysis performed using NICER observations collected during the latest outburst of the source. Doppler modulation of the X-ray pulsation revealed the ultra-compact nature of the binary system characterised by an orbital period of… ▽ More The detection of coherent X-ray pulsations at ~314 Hz (3.2 ms) classifies MAXI J1957+032 as a fast-rotating, accreting neutron star. We present the temporal and spectral analysis performed using NICER observations collected during the latest outburst of the source. Doppler modulation of the X-ray pulsation revealed the ultra-compact nature of the binary system characterised by an orbital period of ~1 hour and a projected semi-major axis of 14 lt-ms. The neutron star binary mass function suggests a minimum donor mass of 1.7e-2 Msun, assuming a neutron star mass of 1.4 Msun and a binary inclination angle lower than 60 degrees. This assumption is supported by the lack of eclipses or dips in the X-ray light curve of the source. We characterised the 0.5-10 keV energy spectrum of the source in outburst as the superposition of a relatively cold black-body-like thermal emission compatible with the emission from the neutron star surface and a Comptonisation component with photon index consistent with a typical hard state. We did not find evidence for iron K-alpha lines or reflection components. △ Less

Submitted 11 August, 2022; originally announced August 2022.

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

The discovery of the 528.6 Hz accreting millisecond X-ray pulsar MAXI J1816-195

Authors: Peter Bult, Diego Altamirano, Zaven Arzoumanian, Deepto Chakrabarty, Jérôme Chenevez, Elizabeth C. Ferrara, Keith C. Gendreau, Sebastien Guillot, Tolga Güver, Wataru Iwakiri, Gaurava K. Jaisawal, Giulio C. Mancuso, Christian Malacaria, Mason Ng, Andrea Sanna, Tod E. Strohmayer, Zorawar Wadiasingh, Michael T. Wolff

Abstract: We present the discovery of 528.6 Hz pulsations in the new X-ray transient MAXI J1816-195. Using NICER, we observed the first recorded transient outburst from the neutron star low-mass X-ray binary MAXI J1816-195 over a period of 28 days. From a timing analysis of the 528.6 Hz pulsations, we find that the binary system is well described as a circular orbit with an orbital period of 4.8 hours and a… ▽ More We present the discovery of 528.6 Hz pulsations in the new X-ray transient MAXI J1816-195. Using NICER, we observed the first recorded transient outburst from the neutron star low-mass X-ray binary MAXI J1816-195 over a period of 28 days. From a timing analysis of the 528.6 Hz pulsations, we find that the binary system is well described as a circular orbit with an orbital period of 4.8 hours and a projected semi-major axis of 0.26 light-seconds for the pulsar, which constrains the mass of the donor star to $0.10-0.55 M_\odot$. Additionally, we observed 15 thermonuclear X-ray bursts showing a gradual evolution in morphology over time, and a recurrence time as short as 1.4 hours. We did not detect evidence for photospheric radius expansion, placing an upper limit on the source distance of 8.6 kpc. △ Less

Submitted 9 August, 2022; originally announced August 2022.

Comments: 8 pages, 3 figures, 2 tables. Accepted for publication in ApJ Letters

SWASTi-SW: Space Weather Adaptive SimulaTion framework for Solar Wind and its relevance to ADITYA-L1 mission

Authors: Prateek Mayank, Bhargav Vaidya, D. Chakrabarty

Abstract: Solar wind streams, acting as background, govern the propagation of space weather drivers in the heliosphere, which induce geomagnetic storm activities. Therefore, predictions of the solar wind parameters are the core of space weather forecasts. This work presents an indigenous three-dimensional (3D) Solar Wind model (SWASTi-SW). This numerical framework for forecasting the ambient solar wind is b… ▽ More Solar wind streams, acting as background, govern the propagation of space weather drivers in the heliosphere, which induce geomagnetic storm activities. Therefore, predictions of the solar wind parameters are the core of space weather forecasts. This work presents an indigenous three-dimensional (3D) Solar Wind model (SWASTi-SW). This numerical framework for forecasting the ambient solar wind is based on a well-established scheme that uses a semi-empirical coronal model and a physics-based inner heliospheric model. This study demonstrates a more generalized version of Wang-Sheeley-Arge (WSA) relation, which provides a speed profile input to the heliospheric domain. Line-of-sight observations of GONG and HMI magnetograms are used as inputs for the coronal model, which in turn, provides the solar wind plasma properties at 0.1 AU. These results are then used as an initial boundary condition for the magnetohydrodynamics (MHD) model of the inner heliosphere to compute the solar wind properties up to 2.1 AU. Along with the validation run for multiple Carrington rotations, the effect of variation of specific heat ratio and study of stream interaction region (SIR) is also presented. This work showcases the multi-directional features of SIRs and provides synthetic measurements for potential observations from the Solar Wind Ion Spectrometer (SWIS) subsystem of Aditya Solar wind Particle EXperiment (ASPEX) payload on-board ISRO's upcoming solar mission Aditya-L1. △ Less

Submitted 27 July, 2022; originally announced July 2022.

Comments: Accepted for publication in ApJS

Outflows and spectral evolution in the eclipsing AMXP SWIFT J1749.4-2807 with NICER, XMM-Newton and NuSTAR

Authors: A. Marino, A. Anitra, S. M. Mazzola, T. Di Salvo, A. Sanna, P. Bult, S. Guillot, G. Mancuso, M. Ng, A. Riggio, A. C. Albayati, D. Altamirano, Z. Arzoumanian, L. Burderi, C. Cabras, D. Chakrabarty, N. Deiosso, K. C. Gendreau, R. Iaria, A. Manca, T. E. Strohmayer

Abstract: The neutron star low-mass X-ray binary SWIFT J1749.4-2807 is the only known eclipsing accreting millisecond X-ray pulsar. In this manuscript we perform a spectral characterization of the system throughout its 2021, two-week-long outburst, analyzing 11 NICER observations and quasi-simultaneous XMM-Newton and NuSTAR single observations at the outburst peak. The broadband spectrum is well-modeled wit… ▽ More The neutron star low-mass X-ray binary SWIFT J1749.4-2807 is the only known eclipsing accreting millisecond X-ray pulsar. In this manuscript we perform a spectral characterization of the system throughout its 2021, two-week-long outburst, analyzing 11 NICER observations and quasi-simultaneous XMM-Newton and NuSTAR single observations at the outburst peak. The broadband spectrum is well-modeled with a black body component with a temperature of $\sim$0.6 keV, most likely consistent with a hot spot on the neutron star surface, and a Comptonisation spectrum with power-law index $Γ\sim 1.9$, arising from a hot corona at $\sim$12 keV. No direct emission from the disc was found, possibly due to it being too cool. A high truncation radius for the disc, i.e., at $\sim$20--30 R$_{G}$ , was obtained from the analysis of the broadened profile of the Fe line in the reflection component. The significant detection of a blue-shifted Fe XXVI absorption line at $\sim$7 keV indicates weakly relativistic X-ray disc winds, which are typically absent in the hard state of X-ray binaries. By comparing the low flux observed during the outburst and the one expected in a conservative mass-transfer, we conclude that mass-transfer in the system is highly non-conservative, as also suggested by the wind detection. Finally, using the Nicer spectra alone, we followed the system while it was fading to quiescence. During the outburst decay, as the spectral shape hardened, the hot spot on the neutron star surface cooled down and shrank, a trend which could be consistent with the pure power-law spectrum observed during quiescence. △ Less

Submitted 18 July, 2022; originally announced July 2022.

Comments: 15 pages, 9 figure; accepted for publication on MNRAS

A novel method to identify IMF $B_y$ influence over low-latitude ionosphere from the correlation of magnetic field variations at nearly antipodal stations

Authors: Sumanjit Chakraborty, D. Chakrabarty

Abstract: This work is directed to explore the efficacy of strong and consistent correlations between the $Δ$X component of magnetic field at two equatorial/low-latitude stations at nearly antipodal locations during geomagnetic storms. The antipodal stations considered are Huancayo (HUA: 12.06$^\circ$S, 75.21$^\circ$W geographic; magnetic dip 0.3$^\circ$N) in the Peruvian longitude sector and Alibag (ABG: 1… ▽ More This work is directed to explore the efficacy of strong and consistent correlations between the $Δ$X component of magnetic field at two equatorial/low-latitude stations at nearly antipodal locations during geomagnetic storms. The antipodal stations considered are Huancayo (HUA: 12.06$^\circ$S, 75.21$^\circ$W geographic; magnetic dip 0.3$^\circ$N) in the Peruvian longitude sector and Alibag (ABG: 18.64$^\circ$N, 72.87$^\circ$E geographic; magnetic dip 10$^\circ$N) in the Indian longitude sector. Six strong geomagnetic storm events during the period 2000-2005, falling in the maximum-to-declining phase of solar cycle 23, are analyzed for this study. These stations are part of the SuperMAG network and data from these stations are openly available. It is noted that although $Δ$X variations over these two stations are, in general, uncorrelated, significant correlations are observed on certain occasions. Correlation coefficient of at least 75$\%$ for 10 minutes is defined as a requisite criterion to infer the possible connection between the $Δ$X variations over these two stations. The ionospheric convection maps from the SuperDARN network are also used to understand the Spatio-temporal evolution of the two-cell ionospheric convection patterns over high-latitudes during these periods of observations. This exercise reveals that the $Δ$X variations over the antipodal locations are significantly correlated when the two-cell convection maps show appropriate rotations and both the stations possibly come under the single cell. Therefore, this investigation brings out a novel method to identify the IMF $B_y$ influence over the low/equatorial latitudes based on the openly available data. △ Less

Submitted 24 January, 2023; v1 submitted 13 July, 2022; originally announced July 2022.

Comments: Manuscript will be submitted soon

SRGA J181414.6-225604: A new Galactic symbiotic X-ray binary outburst triggered by an intense mass loss episode of a heavily obscured Mira variable

Authors: Kishalay De, Ilya Mereminskiy, Roberto Soria, Charlie Conroy, Erin Kara, Shreya Anand, Michael C. B. Ashley, Martha L. Boyer, Deepto Chakrabarty, Brian Grefenstette, Matthew J. Hankins, Lynne A. Hillenbrand, Jacob E. Jencson, Viraj Karambelkar, Mansi M. Kasliwal, Ryan M. Lau, Alexander Lutovinov, Anna M. Moore, Mason Ng, Christos Panagiotou, Dheeraj R. Pasham, Andrey Semena, Robert Simcoe, Jamie Soon, Gokul P. Srinivasaragavan , et al. (2 additional authors not shown)

Abstract: We present the discovery and multi-wavelength characterization of SRGA J181414.6-225604, a Galactic hard X-ray transient discovered during the ongoing SRG/ART-XC sky survey. Using data from the Palomar Gattini-IR survey, we identify a spatially and temporally coincident variable infrared (IR) source, IRAS 18111-2257, and classify it as a very late-type (M7-M8), long period ($1502 \pm 24$ days) and… ▽ More We present the discovery and multi-wavelength characterization of SRGA J181414.6-225604, a Galactic hard X-ray transient discovered during the ongoing SRG/ART-XC sky survey. Using data from the Palomar Gattini-IR survey, we identify a spatially and temporally coincident variable infrared (IR) source, IRAS 18111-2257, and classify it as a very late-type (M7-M8), long period ($1502 \pm 24$ days) and luminous ($M_K\approx -9.9 \pm 0.2$) O-rich Mira donor star located at a distance of $\approx 14.6^{+2.9}_{-2.3}$ kpc. Combining multi-color photometric data over the last $\approx 25$ years, we show that the IR counterpart underwent a recent (starting $\approx 800$ days before the X-ray flare) enhanced mass loss (reaching $\approx 2.1 \times 10^{-5}$ M$_\odot$ yr$^{-1}$) episode resulting in an expanding dust shell obscuring the underlying star. Multi-epoch follow-up from Swift, NICER and NuSTAR reveal a $\approx 200$ day long X-ray outburst reaching a peak luminosity of $L_X \approx 2.5 \times 10^{36}$ erg s$^{-1}$, characterized by a heavily absorbed ($N_{\rm H} \approx 6\times 10^{22}$ cm$^{-2}$) X-ray spectrum consistent with an optically thick Comptonized plasma. The X-ray spectral and timing behavior suggest the presence of clumpy wind accretion together with a dense ionized nebula overabundant in silicate material surrounding the compact object. Together, we show that SRGA J181414.6-225604 is a new symbiotic X-ray binary in outburst, triggered by an intense dust formation episode of a highly evolved donor. Our results offer the first direct confirmation for the speculated connection between enhanced late-stage donor mass loss and active lifetimes of the symbiotic X-ray binaries. △ Less

Submitted 18 May, 2022; originally announced May 2022.

Comments: 32 pages, 15 figures, Accepted for publication in ApJ

A 62-minute orbital period black widow binary in a wide hierarchical triple

Authors: Kevin B. Burdge, Thomas R. Marsh, Jim Fuller, Eric C. Bellm, Ilaria Caiazzo, Deepto Chakrabarty, Michael W. Coughlin, Kishalay De, V. S. Dhillon, Matthew J. Graham, Pablo Rodrí guez-Gil, Amruta D. Jaodand, David L. Kaplan, Erin Kara, Albert K. H. Kong, S. R. Kulkarni, Kwan-Lok Li, S. P. Littlefair, Walid A. Majid, Przemek Mróz, Aaron B. Pearlman, E. S. Phinney, Jan van Roestel, Robert A. Simcoe, Igor Andreoni , et al. (8 additional authors not shown)

Abstract: Over a dozen millisecond pulsars are ablating low-mass companions in close binary systems. In the original "black widow", the 8-hour orbital period eclipsing pulsar PSR J1959+2048 (PSR B1957+20), high energy emission originating from the pulsar is irradiating and may eventually destroy a low-mass companion. These systems are not only physical laboratories that reveal the dramatic result of exposin… ▽ More Over a dozen millisecond pulsars are ablating low-mass companions in close binary systems. In the original "black widow", the 8-hour orbital period eclipsing pulsar PSR J1959+2048 (PSR B1957+20), high energy emission originating from the pulsar is irradiating and may eventually destroy a low-mass companion. These systems are not only physical laboratories that reveal the dramatic result of exposing a close companion star to the relativistic energy output of a pulsar, but are also believed to harbour some of the most massive neutron stars, allowing for robust tests of the neutron star equation of state. Here, we report observations of ZTF J1406+1222, a wide hierarchical triple hosting a 62-minute orbital period black widow candidate whose optical flux varies by a factor of more than 10. ZTF J1406+1222 pushes the boundaries of evolutionary models, falling below the 80 minute minimum orbital period of hydrogen-rich systems. The wide tertiary companion is a rare low metallicity cool subdwarf star, and the system has a Galactic halo orbit consistent with passing near the Galactic center, making it a probe of formation channels, neutron star kick physics, and binary evolution. △ Less

Submitted 4 May, 2022; originally announced May 2022.

Comments: 63 pages, 15 figures, 2 tables, published in Nature on May 4, 2022

Journal ref: Nature 605, 41-45 (2022)

A holistic approach to understand Helium enrichment in Interplanetary coronal mass ejections: New insights

Authors: Yogesh, D. Chakrabarty, N. Srivastava

Abstract: Despite helium abundance (AHe = nH/nHe) is ~ 8 % at the solar photospheric/chromospheric heights, AHe can be found to exceed 8% in interplanetary coronal mass ejections (ICMEs) on many occasions. Although various factors like interplanetary shocks, chromospheric evaporation and "sludge removal" have been separately invoked in the past to address the AHe enhancements in ICMEs, none of these process… ▽ More Despite helium abundance (AHe = nH/nHe) is ~ 8 % at the solar photospheric/chromospheric heights, AHe can be found to exceed 8% in interplanetary coronal mass ejections (ICMEs) on many occasions. Although various factors like interplanetary shocks, chromospheric evaporation and "sludge removal" have been separately invoked in the past to address the AHe enhancements in ICMEs, none of these processes could explain the variability of AHe in ICMEs comprehensively. Based on extensive analysis of 275 ICME events, we show that there is a solar activity variation of ICME averaged AHe values. The investigation also reveals that the first ionization potential effect as well as coronal temperature are not the major contributing factors for AHe enhancements in ICMEs. Investigation on concurrent solar flares and ICME events for 63 cases reveals that chromospheric evaporation in tandem with gravitational settling determine the AHe enhancements and variabilities beyond 8% in ICMEs. While chromospheric evaporation releases the helium from chromosphere into the corona, the gravitationally settled helium is thrown out during the ICMEs.We show that the intensity and timing of the preceding flares from the same active region from where the CME erupts are important factors to understand the AHe enhancements in ICMEs. △ Less

Submitted 3 February, 2022; originally announced February 2022.

Radius Constraints from Reflection Modeling of Cygnus X-2 with NuSTAR and NICER

Authors: R. M. Ludlam, E. M. Cackett, J. A. García, J. M. Miller, A. L. Stevens, A. C. Fabian, J. Homan, M. NG, S. Guillot, D. J. K. Buisson, D. Chakrabarty

Abstract: We present a spectral analysis of NuSTAR and NICER observations of the luminous, persistently accreting neutron star (NS) low-mass X-ray binary Cygnus X-2. The data were divided into different branches that the source traces out on the Z-track of the X-ray color-color diagram; namely the horizontal branch, normal branch, and the vertex between the two. The X-ray continuum spectrum was modeled in t… ▽ More We present a spectral analysis of NuSTAR and NICER observations of the luminous, persistently accreting neutron star (NS) low-mass X-ray binary Cygnus X-2. The data were divided into different branches that the source traces out on the Z-track of the X-ray color-color diagram; namely the horizontal branch, normal branch, and the vertex between the two. The X-ray continuum spectrum was modeled in two different ways that produced a comparable quality fit. The spectra showed clear evidence of a reflection component in the form of a broadened Fe K line, as well as a lower energy emission feature near 1 keV likely due to an ionized plasma located far from the innermost accretion disk. We account for the reflection spectrum with two independent models (relxillns and rdblur*rfxconv). The inferred inclination is in agreement with earlier estimates from optical observations of ellipsoidal light curve modeling (relxillns: $i=67^{\circ}\pm4^{\circ}$, rdblur*rfxconv: $i=60^{\circ}\pm10^{\circ}$). The inner disk radius remains close to the NS ($R_{\rm in}\leq1.15\ R_{\mathrm{ISCO}}$) regardless of the source position along the Z-track or how the 1 keV feature is modeled. Given the optically determined NS mass of $1.71\pm0.21\ M_{\odot}$, this corresponds to a conservative upper limit of $R_{\rm in}\leq19.5$ km for $M=1.92\ M_{\odot}$ or $R_{\rm in}\leq15.3$ km for $M=1.5\ M_{\odot}$. We compare these radius constraints to those obtained from NS gravitational wave merger events and recent NICER pulsar light curve modeling measurements. △ Less

Submitted 27 January, 2022; originally announced January 2022.

Comments: Accepted for publication in ApJ, 13 pages, 5 tables, 5 figures

A massive AGB donor in Scutum X-1: Identification of the first Mira variable in an X-ray binary

Authors: Kishalay De, Deepto Chakrabarty, Roberto Soria, Michael C. B. Ashley, Charlie Conroy, Matthew J. Hankins, Mansi M. Kasliwal, Ryan M. Lau, Anna M. Moore, Robert Simcoe, Jamie Soon, Tony Travouillon

Abstract: The symbiotic X-ray binary Sct X-1 was suggested as the first known neutron star accreting from a red supergiant companion. Although known for nearly 50 years, detailed characterization of the donor remains lacking, particularly due to the extremely high reddening towards the source ($A_V\gtrsim25$ mag). Here, we present i) improved localization of the counterpart using Gaia and Chandra observatio… ▽ More The symbiotic X-ray binary Sct X-1 was suggested as the first known neutron star accreting from a red supergiant companion. Although known for nearly 50 years, detailed characterization of the donor remains lacking, particularly due to the extremely high reddening towards the source ($A_V\gtrsim25$ mag). Here, we present i) improved localization of the counterpart using Gaia and Chandra observations, ii) the first broadband infrared spectrum ($\approx1-5\,μ$m; $R\approx 2000$) obtained with SpeX on the NASA Infrared Telescope Facility and iii) $J$-band light curve from the Palomar Gattini-IR survey. The infrared spectrum is characterized by i) deep water absorption features (H$_2$O index $\approx 40$%), ii) strong TiO, VO and CO features, and iii) weak/absent CN lines. We show that these features are inconsistent with known red supergiants, but suggest a M8-9 III type O-rich Mira donor star. We report the discovery of large amplitude ($ΔJ\approx3.5$ mag) periodic photometric variability suggesting a pulsation period of $621\pm36\,{\rm(systematic)}\pm8\,{\rm(statistical)}$ days, which we use to constrain the donor to be a relatively luminous Mira ($M_K=-8.6\pm0.3$ mag) at a distance of $3.6^{+0.8}_{-0.7}$ kpc. Comparing these characteristics to recent models, we find the donor to be consistent with a $\approx 3-5$ M$_\odot$ star at an age of $\approx 0.1-0.3$ Gyr. Together, we show that Sct X-1 was previously mis-classified as an evolved High Mass X-ray Binary; instead it is an intermediate mass system with the first confirmed Mira donor in an X-ray binary. We discuss the implications of Mira donors in symbiotic X-ray binaries, and highlight the potential of wide-field infrared time domain surveys and broadband infrared spectroscopy to unveil their demographics. △ Less

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

Comments: 10 pages, 4 figures. Accepted for publication in ApJ Letters

Differential behaviors of suprathermal $^4$He and Fe populations in the interplanetary medium during solar cycle 24

Authors: Bijoy Dalal, Dibyendu Chakrabarty, Nandita Srivastava

Abstract: Investigations on the solar cycle variation of the properties of suprathermal populations (H and other heavy ions like $^4$He, $^3$He, C, O and Fe) in the solar wind are sparse and hence, poorly understood. In the present investigation, solar cycle variations of "quiet" time suprathermal elements are investigated using $<$ $\sim$ 1 MeV/n particle flux data obtained from Ultra-Low Energy Isotope Sp… ▽ More Investigations on the solar cycle variation of the properties of suprathermal populations (H and other heavy ions like $^4$He, $^3$He, C, O and Fe) in the solar wind are sparse and hence, poorly understood. In the present investigation, solar cycle variations of "quiet" time suprathermal elements are investigated using $<$ $\sim$ 1 MeV/n particle flux data obtained from Ultra-Low Energy Isotope Spectrometer on board Advanced Composition Explorer satellite during the solar cycle 23 and 24. The analysis reveals that helium ($^4$He) shows zero or positive lags with respect to sunspot numbers in solar cycle 23 while it shows zero or negative lag in solar cycle 24. On the contrary, although iron (Fe) shows zero or positive lag in cycle 23 similar to $^4$He, it shows only zero lag in cycle 24 and no negative lag is seen. Further, significant differences in the spectral indices are seen between $^4$He and Fe in cycle 24 compared to the cycle 23. These results suggest that generation mechanisms responsible for suprathermal $^4$He and Fe underwent changes in cycle 24 and these mechanisms are probably dependent on the first ionization potential and mass to charge ratio. This proposition gets credence from the fact that changes in the lag and spectral slopes for C and O are not significantly different in cycle 23 and 24. △ Less

Submitted 11 August, 2022; v1 submitted 25 December, 2021; originally announced December 2021.

Comments: 18 pages, 7 (+12 supplementary) figures, accepted for publication in the Astrophysical Journal (ApJ)

Evidence for a Compact Object in the Aftermath of the Extra-Galactic Transient AT2018cow

Authors: Dheeraj R. Pasham, Wynn C. G. Ho, William Alston, Ronald Remillard, Mason Ng, Keith Gendreau, Brian D. Metzger, Diego Altamirano, Deepto Chakrabarty, Andrew Fabian, Jon Miller, Peter Bult, Zaven Arzoumanian, James F. Steiner, Tod Strohmayer, Francesco Tombesi, Jeroen Homan, Edward M. Cackett, Alice Harding

Abstract: The brightest Fast Blue Optical Transients (FBOTs) are mysterious extragalactic explosions that may represent a new class of astrophysical phenomena. Their fast time to maximum brightness of less than a week and decline over several months and atypical optical spectra and evolution are difficult to explain within the context of core-collapse of massive stars which are powered by radioactive decay… ▽ More The brightest Fast Blue Optical Transients (FBOTs) are mysterious extragalactic explosions that may represent a new class of astrophysical phenomena. Their fast time to maximum brightness of less than a week and decline over several months and atypical optical spectra and evolution are difficult to explain within the context of core-collapse of massive stars which are powered by radioactive decay of Nickel-56 and evolve more slowly. AT2018cow (at redshift of 0.014) is an extreme FBOT in terms of rapid evolution and high luminosities. Here we present evidence for a high-amplitude quasi-periodic oscillation (QPO) of AT2018cow's soft X-rays with a frequency of 224 Hz (at 3.7$σ$ significance level or false alarm probability of 0.02%) and fractional root-mean-squared amplitude of >30%. This signal is found in the average power density spectrum taken over the entire 60-day outburst and suggests a highly persistent signal that lasts for a billion cycles. The high frequency (rapid timescale) of 224 Hz (4.4 ms) argues for a compact object in AT2018cow, which can be a neutron star or black hole with a mass less than 850 solar masses. If the QPO is the spin period of a neutron star, we can set limits on the star's magnetic field strength. Our work highlights a new way of using high time-resolution X-ray observations to study FBOTs. △ Less

Submitted 20 December, 2021; v1 submitted 8 December, 2021; originally announced December 2021.

Comments: Published in Nature astronomy on 13th December 2021

A NICER look at thermonuclear X-ray bursts from Aql X-1

Authors: Tolga Guver, Tugba Boztepe, David R. Ballantyne, Z. Funda Bostanci, Peter Bult, Gaurava K. Jaisawal, Ersin Gogus, Tod E. Strohmayer, Diego Altamirano, Sebastien Guillot, Deepto Chakrabarty

Abstract: We present spectral and temporal properties of all the thermonuclear X-ray bursts observed from Aql X-1 by the Neutron Star Interior and Composition Explorer (NICER) between 2017 July and 2021 April. This is the first systematic investigation of a large sample of type I X-ray bursts from Aql X-1 with improved sensitivity at low energies. We detect 22 X-ray bursts including two short recurrence bur… ▽ More We present spectral and temporal properties of all the thermonuclear X-ray bursts observed from Aql X-1 by the Neutron Star Interior and Composition Explorer (NICER) between 2017 July and 2021 April. This is the first systematic investigation of a large sample of type I X-ray bursts from Aql X-1 with improved sensitivity at low energies. We detect 22 X-ray bursts including two short recurrence burst events in which the separation was only 451 s and 496 s. We perform time resolved spectroscopy of the bursts using the fixed and scaled background (f_a method) approaches. We show that the use of a scaling factor to the pre-burst emission is the statistically preferred model in about 68% of all the spectra compared to the fixed background approach. Typically the f_a values are clustered around 1-3, but can reach up to 11 in a burst where photospheric radius expansion is observed. Such f_a values indicate a very significant increase in the pre-burst emission especially at around the peak flux moments of the bursts. We show that the use of the f_a factor alters the best fit spectral parameters of the burst emission. Finally, we employed a reflection model instead of scaling the pre-burst emission. We show that reflection models also do fit the spectra and improve the goodness of the fits. In all cases we see that the disc is highly ionized by the burst emission and the fraction of the reprocessed emission to the incident burst flux is typically clustered around 20%. △ Less

Submitted 23 November, 2021; originally announced November 2021.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society

Search for continuous gravitational waves from 20 accreting millisecond X-ray pulsars in O3 LIGO data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, N. Adhikari, R. X. Adhikari, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, T. Akutsu, S. Albanesi, A. Allocca, P. A. Altin, A. Amato, C. Anand , et al. (1612 additional authors not shown)

Abstract: Results are presented of searches for continuous gravitational waves from 20 accreting millisecond X-ray pulsars with accurately measured spin frequencies and orbital parameters, using data from the third observing run of the Advanced LIGO and Advanced Virgo detectors. The search algorithm uses a hidden Markov model, where the transition probabilities allow the frequency to wander according to an… ▽ More Results are presented of searches for continuous gravitational waves from 20 accreting millisecond X-ray pulsars with accurately measured spin frequencies and orbital parameters, using data from the third observing run of the Advanced LIGO and Advanced Virgo detectors. The search algorithm uses a hidden Markov model, where the transition probabilities allow the frequency to wander according to an unbiased random walk, while the $\mathcal{J}$-statistic maximum-likelihood matched filter tracks the binary orbital phase. Three narrow sub-bands are searched for each target, centered on harmonics of the measured spin frequency. The search yields 16 candidates, consistent with a false alarm probability of 30% per sub-band and target searched. These candidates, along with one candidate from an additional target-of-opportunity search done for SAX J1808.4$-$3658, which was in outburst during one month of the observing run, cannot be confidently associated with a known noise source. Additional follow-up does not provide convincing evidence that any are a true astrophysical signal. When all candidates are assumed non-astrophysical, upper limits are set on the maximum wave strain detectable at 95% confidence, $h_0^{95\%}$. The strictest constraint is $h_0^{95\%} = 4.7\times 10^{-26}$ from IGR J17062$-$6143. Constraints on the detectable wave strain from each target lead to constraints on neutron star ellipticity and $r$-mode amplitude, the strictest of which are $ε^{95\%} = 3.1\times 10^{-7}$ and $α^{95\%} = 1.8\times 10^{-5}$ respectively. This analysis is the most comprehensive and sensitive search of continuous gravitational waves from accreting millisecond X-ray pulsars to date. △ Less

Submitted 21 January, 2022; v1 submitted 19 September, 2021; originally announced September 2021.

Comments: 40 pages, 6 figures. This version contains minor typographical revisions to match published article

Report number: LIGO-P2100221

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

A month of monitoring the new magnetar Swift J1555.2-5402 during an X-ray outburst

Authors: Teruaki Enoto, Mason Ng, Chin-ping Hu, Tolga Guver, Gaurava K. Jaisawal, Brendan O'Connor, Ersin Gogus, Amy Lien, Shota Kisaka, Zorawar Wadiasingh, Walid A. Majid, Aaron B. Pearlman, Zaven Arzoumanian, Karishma Bansal, Harsha Blumer, Deepto Chakrabarty, Keith Gendreau, Wynn C. G. Ho, Chryssa Kouveliotou, Paul S. Ray, Tod E. Strohmayer, George Younes, David M. Palmer, Takanori Sakamoto, Takuya Akahori , et al. (1 additional authors not shown)

Abstract: The soft gamma-ray repeater Swift J1555.2-5402 was discovered by means of a 12-ms duration short burst detected with Swift BAT on 2021 June 3. Then 1.6 hours after the first burst detection, NICER started daily monitoring of this X-ray source for a month. The absorbed 2-10 keV flux stays nearly constant at around 4e-11 erg/s/cm2 during the monitoring timespan, showing only a slight gradual decline… ▽ More The soft gamma-ray repeater Swift J1555.2-5402 was discovered by means of a 12-ms duration short burst detected with Swift BAT on 2021 June 3. Then 1.6 hours after the first burst detection, NICER started daily monitoring of this X-ray source for a month. The absorbed 2-10 keV flux stays nearly constant at around 4e-11 erg/s/cm2 during the monitoring timespan, showing only a slight gradual decline. A 3.86-s periodicity is detected, and the time derivative of this period is measured to be 3.05(7)e-11 s/s. The soft X-ray pulse shows a single sinusoidal shape with a root-mean-square pulsed fraction that increases as a function of energy from 15% at 1.5 keV to 39% at 7 keV. The equatorial surface magnetic field, characteristic age, and spin-down luminosity are derived under the dipole field approximation to be 3.5e+14 G, 2.0 kyr, and 2.1e+34 erg/s, respectively. An absorbed blackbody with a temperature of 1.1 keV approximates the soft X-ray spectrum. Assuming a source distance of 10 kpc, the peak X-ray luminosity is ~8.5e+35 erg/s in the 2--10 keV band. During the period of observations, we detect 5 and 37 short bursts with Swift/BAT and NICER, respectively. Based on these observational properties, especially the inferred strong magnetic field, this new source is classified as a magnetar. We also coordinated hard X-ray and radio observations with NuSTAR, DSN, and VERA. A hard X-ray power-law component that extends up to at least 40 keV is detected at 3-sigma significance. The 10-60 keV flux, which is dominated by the power-law component, is ~9e-12 erg/s/cm2 with a photon index of ~1.2. The pulsed fraction has a sharp cutoff above 10 keV, down to ~10% in the hard-tail component band. No radio pulsations are detected during the DSN nor VERA observations. We place 7σ upper limits of 0.043mJy and 0.026 mJy on the flux density at S-band and X-band, respectively. △ Less

Submitted 6 August, 2021; originally announced August 2021.

Comments: 26 pages, 7 figures (3 figures in the appendix), 10 tables (9 tables in the appendix). Submitted to ApJL

Chandra/LETGS Studies of the Collisional Plasma in 4U 1626-67

Authors: Paul B. Hemphill, Norbert S. Schulz, Herman M. Marshall, Deepto Chakrabarty

Abstract: We present an analysis of \chandra/LETGS observations of the ultracompact X-ray binary (UCXB) 4U 1626$-$67, continuing our project to analyze the existing Chandra gratings data of this interesting source. The extremely low mass, hydrogen-depleted donor star provides a unique opportunity to study the properties and structure of the metal-rich accreted plasma. There are strong, double-peak… ▽ More We present an analysis of \chandra/LETGS observations of the ultracompact X-ray binary (UCXB) 4U 1626$-$67, continuing our project to analyze the existing Chandra gratings data of this interesting source. The extremely low mass, hydrogen-depleted donor star provides a unique opportunity to study the properties and structure of the metal-rich accreted plasma. There are strong, double-peaked emission features of OVII-VIII and Ne IX-X, but no other identified emission lines are detected. Our spectral fit simultaneously models the emission line profiles and the plasma parameters, using a two-temperature collisionally-ionized plasma. Based on our line profile fitting, we constrain the inclination of the system to 25--60$^{\circ}$ and the inner disk radius to $\sim$1500 gravitational radii, in turn constraining the donor mass to $\lesssim$0.026 M_sun, while our plasma modeling confirms previous reports of high neon abundance in the source, establishing a Ne/O ratio in the system of $0.47 \pm 0.04$, while simultaneously estimating a very low Fe/O ratio of $0.0042 \pm 0.0008$ and limiting the Mg/O ratio to less than 1% by number. We discuss these results in light of previous work. △ Less

Submitted 11 June, 2021; originally announced June 2021.

Comments: accepted by the Astrophysical Main Journal, 16 pages, 13 figures

Evidence for distinctive changes in the solar wind helium abundance in cycle 24

Authors: Yogesh, D. Chakrabarty, N. Srivastava

Abstract: The relative abundance of alpha particles with respect to proton, usually expressed as $A_{He}$ = ($n_α/n_p$)*100, is known to respond to solar activity although changes in its behaviour in the last four solar cycles are not known. In this letter, by systematically analysing inter-calibrated $A_{He}$ data obtained from the first Lagrangian point of the Sun-Earth system, we show that $A_{He}$ varia… ▽ More The relative abundance of alpha particles with respect to proton, usually expressed as $A_{He}$ = ($n_α/n_p$)*100, is known to respond to solar activity although changes in its behaviour in the last four solar cycles are not known. In this letter, by systematically analysing inter-calibrated $A_{He}$ data obtained from the first Lagrangian point of the Sun-Earth system, we show that $A_{He}$ variations are distinctively different in solar cycle 24 as compared to the last three cycles. The frequency of $A_{He}$ = 2-3% events is significantly higher in slow/intermediate solar winds in cycle 24 as opposed to the dominance of the typical $A_{He}$ = 4-5% events in the previous three cycles. Further, the occurrence of $A_{He}$ $\geq$ 10% events is significantly reduced in cycle 24. Not only that, the changes in delay of $A_{He}$ with respect to peak sunspot numbers are less sensitive to changes in solar wind velocity in cycle 24. The investigation suggests that the coronal magnetic field configuration started undergoing systematic changes starting from cycle 23 and this altered magnetic field configuration affected the way helium got processed and depleted in the solar atmosphere. △ Less

Submitted 10 February, 2021; originally announced February 2021.

Thermonuclear X-ray Bursts with late secondary peaks observed from 4U 1608-52

Authors: Tolga Guver, Tugba Boztepe, Ersin Gogus, Manoneeta Chakraborty, Tod E. Strohmayer, Peter Bult, Diego Altamirano, Gaurava K. Jaisawal, Tugce Kocabiyik, Christian Malacaria, Unnati Kashyap, Keith C. Gendreau, Zaven Arzoumanian, Deepto Chakrabarty

Abstract: We report the temporal and spectral analysis of three thermonuclear X-ray bursts from 4U 1608-52, observed by the Neutron Star Interior Composition Explorer (NICER) during and just after the outburst observed from the source in 2020. In two of the X-ray bursts, we detect secondary peaks, 30 and 18 seconds after the initial peaks. The secondary peaks show a fast rise exponential decay-like shape re… ▽ More We report the temporal and spectral analysis of three thermonuclear X-ray bursts from 4U 1608-52, observed by the Neutron Star Interior Composition Explorer (NICER) during and just after the outburst observed from the source in 2020. In two of the X-ray bursts, we detect secondary peaks, 30 and 18 seconds after the initial peaks. The secondary peaks show a fast rise exponential decay-like shape resembling a thermonuclear X-ray burst. Time-resolved X-ray spectral analysis reveals that the peak flux, blackbody temperature, and apparent emitting radius values of the initial peaks are in agreement with X-ray bursts previously observed from 4U 1608-52, while the same values for the secondary peaks tend toward the lower end of the distribution of bursts observed from this source. The third X-ray burst, which happened during much lower accretion rates did not show any evidence for a deviation from an exponential decay and was significantly brighter than the previous bursts. We present the properties of the secondary peaks and discuss the events within the framework of short recurrence time bursts or bursts with secondary peaks. We find that the current observations do not fit in standard scenarios and challenge our understanding of flame spreading. △ Less

Submitted 27 January, 2021; originally announced January 2021.

Comments: Accepted for publication in the Astrophysical Journal

NICER Discovery of Millisecond X-ray Pulsations and an Ultracompact Orbit in IGR J17494-3030

Authors: Mason Ng, Paul S. Ray, Peter Bult, Deepto Chakrabarty, Gaurava K. Jaisawal, Christian Malacaria, Diego Altamirano, Zaven Arzoumanian, Keith C. Gendreau, Tolga Güver, Matthew Kerr, Tod E. Strohmayer, Zorawar Wadiasingh, Michael T. Wolff

Abstract: We report the detection of 376.05 Hz (2.66 ms) coherent X-ray pulsations in NICER observations of a transient outburst of the low-mass X-ray binary IGR J17494-3030 in 2020 October/November. The system is an accreting millisecond X-ray pulsar in a 75 minute ultracompact binary. The mass donor is most likely a $\simeq 0.02 M_\odot$ finite-entropy white dwarf composed of He or C/O. The fractional rms… ▽ More We report the detection of 376.05 Hz (2.66 ms) coherent X-ray pulsations in NICER observations of a transient outburst of the low-mass X-ray binary IGR J17494-3030 in 2020 October/November. The system is an accreting millisecond X-ray pulsar in a 75 minute ultracompact binary. The mass donor is most likely a $\simeq 0.02 M_\odot$ finite-entropy white dwarf composed of He or C/O. The fractional rms pulsed amplitude is 7.4%, and the soft (1-3 keV) X-ray pulse profile contains a significant second harmonic. The pulsed amplitude and pulse phase lag (relative to our mean timing model) are energy-dependent, each having a local maximum at 4 keV and 1.5 keV, respectively. We also recovered the X-ray pulsations in archival 2012 XMM-Newton observations, allowing us to measure a long-term pulsar spin-down rate of $\dotν= -2.1(7)\times10^{-14}$ Hz/s and to infer a pulsar surface dipole magnetic field strength of $\simeq 10^9$ G. We show that the mass transfer in the binary is likely non-conservative, and we discuss various scenarios for mass loss from the system. △ Less

Submitted 1 February, 2021; v1 submitted 15 January, 2021; originally announced January 2021.

Comments: 10 pages, 3 figures, and 1 table; accepted for publication in ApJL

The X-ray bursts of XTE J1739-285: a NICER sample

Authors: Peter Bult, Diego Altamirano, Zaven Arzoumanian, Anna V. Bilous, Deepto Chakrabarty, Keith C. Gendreau, Tolga Güver, Gaurava K. Jaisawal, Erik Kuulkers, Christian Malacaria, Mason Ng, Andrea Sanna, Tod E. Strohmayer

Abstract: In this work we report on observations with the Neutron Star Interior Composition Explorer of the known neutron star X-ray transient XTE J1739-285. We observed the source in 2020 February and March, finding it in a highly active bursting state. Across a 20-day period, we detected 32 thermonuclear X-ray bursts, with an average burst recurrence time of $2.0^{+0.4}_{-0.3}$ hr. A timing and spectral a… ▽ More In this work we report on observations with the Neutron Star Interior Composition Explorer of the known neutron star X-ray transient XTE J1739-285. We observed the source in 2020 February and March, finding it in a highly active bursting state. Across a 20-day period, we detected 32 thermonuclear X-ray bursts, with an average burst recurrence time of $2.0^{+0.4}_{-0.3}$ hr. A timing and spectral analysis of the ensemble of X-ray bursts reveals homogeneous burst properties, evidence for short-recurrence time bursts, and the detection of a 386.5 Hz burst oscillation candidate. The latter is especially notable, given that a previous study of this source claimed a 1122 Hz burst oscillation candidate. We did not find any evidence of variability near 1122 Hz, and instead find that the 386.5 Hz oscillation is the more prominent signal of the two burst oscillation candidates. Hence, we conclude it is unlikely that XTE J1739-285 has a sub-millisecond rotation period. △ Less

Submitted 21 December, 2020; originally announced December 2020.

Comments: 15 pages, 9 figures, 2 tables. Accepted for publication in ApJ

Simultaneous NICER and NuSTAR Observations of the Ultra-compact X-ray Binary 4U 1543-624

Authors: R. M. Ludlam, A. D. Jaodand, J. A. García, N. Degenaar, J. A. Tomsick, E. M. Cackett, A. C. Fabian, P. Gandhi, D. J. K. Buisson, A. W. Shaw, D. Chakrabarty

Abstract: We present the first joint NuSTAR and NICER observations of the ultra-compact X-ray binary (UCXB) 4U 1543$-$624 obtained in 2020 April. The source was at a luminosity of $L_{0.5-50\ \mathrm{keV}} = 4.9 (D/7\ \mathrm{kpc})^{2}\times10^{36}$ ergs s$^{-1}$ and showed evidence of reflected emission in the form of an O VIII line, Fe K line, and Compton hump within the spectrum. We used a full reflectio… ▽ More We present the first joint NuSTAR and NICER observations of the ultra-compact X-ray binary (UCXB) 4U 1543$-$624 obtained in 2020 April. The source was at a luminosity of $L_{0.5-50\ \mathrm{keV}} = 4.9 (D/7\ \mathrm{kpc})^{2}\times10^{36}$ ergs s$^{-1}$ and showed evidence of reflected emission in the form of an O VIII line, Fe K line, and Compton hump within the spectrum. We used a full reflection model, known as xillverCO, that is tailored for the atypical abundances found in UCXBs, to account for the reflected emission. We tested the emission radii of the O and Fe line components and conclude that they originate from a common disk radius in the innermost region of the accretion disk ($R_{\rm in} \leq1.07\ R_{\mathrm{ISCO}}$). Assuming that the compact accretor is a neutron star (NS) and the position of the inner disk is the Alfvén radius, we placed an upper limit on the magnetic field strength to be $B\leq0.7(D/7\ \mathrm {kpc})\times10^{8}$ G at the poles. Given the lack of pulsations detected and position of $R_{\rm in}$, it was likely that a boundary layer region had formed between the NS surface and inner edge of the accretion disk with an extent of 1.2 km. This implies a maximum radius of the neutron star accretor of $R_{\mathrm{NS}}\leq 12.1$ km when assuming a canonical NS mass of 1.4 $M_{\odot}$. △ Less

Submitted 9 March, 2021; v1 submitted 18 December, 2020; originally announced December 2020.

Comments: Accepted for publication in ApJ, 10 pages, 6 figures, 1 table

A Small Satellite Version of a Broad-band Soft X-ray Polarimeter

Authors: Herman L. Marshall, Sarah N. T. Heine, Alan Garner, Eric M. Gullikson, H. Moritz Günther, Christopher Leitz, Rebecca Masterson, Eric D. Miller, William Zhang, Rozenn Boissay-Malaquin, Ilaria Caiazzo, Deepto Chakrabarty, Rosemary Davidson, Luigi C. Gallo, Ralf K. Heilmann, Jeremy Heyl, Erin Kara, Alan Marscher, Norbert S. Schulz

Abstract: We describe a new implementation of a broad-band soft X-ray polarimeter, substantially based on a previous design. This implementation, the Pioneer Soft X-ray Polarimeter (PiSoX) is a SmallSat, designed for NASA's call for Astrophysics Pioneers, small missions that could be CubeSats, balloon experiments, or SmallSats. As in the REDSoX Polarimeter, the grating arrangement is designed optimally for… ▽ More We describe a new implementation of a broad-band soft X-ray polarimeter, substantially based on a previous design. This implementation, the Pioneer Soft X-ray Polarimeter (PiSoX) is a SmallSat, designed for NASA's call for Astrophysics Pioneers, small missions that could be CubeSats, balloon experiments, or SmallSats. As in the REDSoX Polarimeter, the grating arrangement is designed optimally for the purpose of polarimetry with broad-band focussing optics by matching the dispersion of the spectrometer channels to laterally graded multilayers (LGMLs). The system can achieve polarization modulation factors over 90%. For PiSoX, the optics are lightweight Si mirrors in a one-bounce parabolic configuration. High efficiency, blazed gratings from opposite sectors are oriented to disperse to a LGML forming a channel covering the wavelength range from 35 to 75 Angstroms (165 - 350 eV). Upon satellite rotation, the intensities of the dispersed spectra, after reflection and polarizing by the LGMLs, give the three Stokes parameters needed to determine a source's linear polarization fraction and orientation. The design can be extended to higher energies as LGMLs are developed further. We describe examples of the potential scientific return from instruments based on this design. △ Less

Submitted 4 December, 2020; originally announced December 2020.

Comments: 20 pages, 8 figures, to appear in Proceedings SPIE, volume 11444

Discovery of Thermonuclear Type-I X-ray Bursts from the X-ray binary MAXI J1807+132

Authors: A. C. Albayati, D. Altamirano, G. K. Jaisawal, P. Bult, S. Rapisarda, G. C. Mancuso, T. Güver, Z. Arzoumanian, D. Chakrabarty, J. Chenevez, J. M. C. Court, K. C. Gendreau, S. Guillot, L. Keek, C. Malacaria, T. E. Strohmayer

Abstract: MAXI J1807+132 is a low-mass X-ray binary (LMXB) first detected in outburst in 2017. Observations during the 2017 outburst did not allow for an unambiguous identification of the nature of the compact object. MAXI J1807+132 was detected in outburst again in 2019 and was monitored regularly with NICER. In this paper we report on five days of observations during which we detected three thermonuclear… ▽ More MAXI J1807+132 is a low-mass X-ray binary (LMXB) first detected in outburst in 2017. Observations during the 2017 outburst did not allow for an unambiguous identification of the nature of the compact object. MAXI J1807+132 was detected in outburst again in 2019 and was monitored regularly with NICER. In this paper we report on five days of observations during which we detected three thermonuclear (Type-I) X-ray bursts, identifying the system as a neutron star LMXB. Time-resolved spectroscopy of the three Type-I bursts revealed typical characteristics expected for these phenomena. All three Type-I bursts show slow rises and long decays, indicative of mixed H/He fuel. We find no strong evidence that any of the Type-I bursts reached the Eddington Luminosity; however, under the assumption that the brightest X-ray burst underwent photospheric radius expansion, we estimate a <12.4kpc upper limit for the distance. We searched for burst oscillations during the Type-I bursts from MAXI J1807+132 and found none (<10% amplitude upper limit at 95% confidence level). Finally, we found that the brightest Type-I burst shows a ~1.6sec pause during the rise. This pause is similar to one recently found with NICER in a bright Type-I burst from the accreting millisecond X-ray pulsar SAX J1808.4-3658. The fact that Type-I bursts from both sources can show this type of pause suggests that the origin of the pauses is independent of the composition of the burning fuel, the peak luminosity of the Type-I bursts, or whether the NS is an X-ray pulsar. △ Less

Submitted 20 November, 2020; originally announced November 2020.

Comments: 9 pages, 4 figures

Journal ref: MNRAS, 501 (2021) 261-268

Discovery of thermonuclear (Type I) X-ray bursts in the X-ray binary Swift J1858.6-0814 observed with NICER and NuSTAR

Authors: D. J. K. Buisson, D. Altamirano, P. Bult, G. C. Mancuso, T. Güver, G. K. Jaisawal, J. Hare, A. C. Albayati, Z. Arzoumanian, N. Castro Segura, D. Chakrabarty, P. Gandhi, S. Guillot, J. Homan, K. C. Gendreau, J. Jiang, C. Malacaria, J. M. Miller, M. Özbey Arabacı, R. Remillard, T. E. Strohmayer, F. Tombesi, J. A. Tomsick, F. M. Vincentelli, D. J. Walton

Abstract: Swift J1858.6-0814 is a recently discovered X-ray binary notable for extremely strong variability (by factors $>100$ in soft X-rays) in its discovery state. We present the detection of five thermonuclear (Type I) X-ray bursts from Swift J1858.6-0814, implying that the compact object in the system is a neutron star. Some of the bursts show photospheric radius expansion, so their peak flux can be us… ▽ More Swift J1858.6-0814 is a recently discovered X-ray binary notable for extremely strong variability (by factors $>100$ in soft X-rays) in its discovery state. We present the detection of five thermonuclear (Type I) X-ray bursts from Swift J1858.6-0814, implying that the compact object in the system is a neutron star. Some of the bursts show photospheric radius expansion, so their peak flux can be used to estimate the distance to the system. The peak luminosity, and hence distance, can depend on several system parameters; for the most likely values, a high inclination and a helium atmosphere, $D=12.8_{-0.6}^{+0.8}$ kpc, although systematic effects allow a conservative range of $9-18$ kpc. Before one burst, we detect a QPO at $9.6\pm0.5$ mHz with a fractional rms amplitude of $2.2\pm0.2$% ($0.5-10$ keV), likely due to marginally stable burning of helium; similar oscillations may be present before the other bursts but the light curves are not long enough to allow their detection. We also search for burst oscillations but do not detect any, with an upper limit in the best case of 15% fractional amplitude (over $1-8$ keV). Finally, we discuss the implications of the neutron star accretor and this distance on other inferences which have been made about the system. In particular, we find that Swift J1858.6-0814 was observed at super-Eddington luminosities at least during bright flares during the variable stage of its outburst. △ Less

Submitted 7 September, 2020; originally announced September 2020.

Comments: 14 pages, 8 figures, 2 tables, MNRAS accepted

Data Processing Software for Chandrayaan-2 Solar X-ray Monitor

Authors: N. P. S. Mithun, Santosh V. Vadawale, Arpit R. Patel, M. Shanmugam, D. Chakrabarty, Partha Konar, Tejas N. Sarvaiya, Girish D. Padia, Aveek Sarkar, Prashant Kumar, Prashant Jangid, Aaditya Sarda, Manan S. Shah, Anil Bhardwaj

Abstract: Solar X-ray Monitor (XSM) instrument of India's Chandrayaan-2 lunar mission carries out broadband spectroscopy of the Sun in soft X-rays. XSM, with its unique features such as low background, high time cadence, and high spectral resolution, provides the opportunity to characterize transient and quiescent X-ray emission from the Sun even during low activity periods. It records the X-ray spectrum at… ▽ More Solar X-ray Monitor (XSM) instrument of India's Chandrayaan-2 lunar mission carries out broadband spectroscopy of the Sun in soft X-rays. XSM, with its unique features such as low background, high time cadence, and high spectral resolution, provides the opportunity to characterize transient and quiescent X-ray emission from the Sun even during low activity periods. It records the X-ray spectrum at one-second cadence, and the data recorded on-board are downloaded at regular intervals along with that of other payloads. During ground pre-processing, the XSM data is segregated, and the level-0 data is made available for higher levels of processing at the Payload Operations Center (POC). XSM Data Analysis Software (XSMDAS) is developed to carry out the processing of the level-0 data to higher levels and to generate calibrated light curves and spectra for user-defined binning parameters such that it is suitable for further scientific analysis. A front-end for the XSMDAS named XSM Quick Look Display (XSMQLD) is also developed to facilitate a first look at the data without applying calibration. XSM Data Management-Monitoring System (XSMDMS) is designed to carry out automated data processing at the POC and to maintain an SQLite database with relevant information on the data sets and an internal web application for monitoring data quality and instrument health. All XSM raw and calibrated data products are in FITS format, organized into day-wise files, and the data archive follows Planetary Data System-4 (PDS4) standards. The XSM data will be made available after a lock-in period along with the XSM Data Analysis Software from ISRO Science Data Archive (ISDA) at Indian Space Science Data Center(ISSDC). Here we discuss the design and implementation of all components of the software for the XSM data processing and the contents of the XSM data archive. △ Less

Submitted 22 July, 2020; originally announced July 2020.

Comments: Submitted to Astronomy and Computing, 13 pages, 9 figures

Journal ref: Astronomy and Computing 34, 100449 (2021)

A strongly changing accretion morphology during the outburst decay of the neutron star X-ray binary 4U 1608-52

Authors: J. van den Eijnden, N. Degenaar, R. M. Ludlam, A. S. Parikh, J. M. Miller, R. Wijnands, K. C. Gendreau, Z. Arzoumanian, D. Chakrabarty, P. Bult

Abstract: It is commonly assumed that the properties and geometry of the accretion flow in transient low-mass X-ray binaries (LMXBs) significantly change when the X-ray luminosity decays below $\sim 10^{-2}$ of the Eddington limit ($L_{\rm Edd}$). However, there are few observational cases where the evolution of the accretion flow is tracked in a single X-ray binary over a wide dynamic range. In this work,… ▽ More It is commonly assumed that the properties and geometry of the accretion flow in transient low-mass X-ray binaries (LMXBs) significantly change when the X-ray luminosity decays below $\sim 10^{-2}$ of the Eddington limit ($L_{\rm Edd}$). However, there are few observational cases where the evolution of the accretion flow is tracked in a single X-ray binary over a wide dynamic range. In this work, we use NuSTAR and NICER observations obtained during the 2018 accretion outburst of the neutron star LMXB 4U 1608-52, to study changes in the reflection spectrum. We find that the broad Fe-K$α$ line and Compton hump, clearly seen during the peak of the outburst when the X-ray luminosity is $\sim 10^{37}$ erg/s ($\sim 0.05$ $L_{\rm Edd}$), disappear during the decay of the outburst when the source luminosity drops to $\sim 4.5 \times 10^{35}$ erg/s ($\sim 0.002$ $L_{\rm Edd}$). We show that this non-detection of the reflection features cannot be explained by the lower signal-to-noise at lower flux, but is instead caused by physical changes in the accretion flow. Simulating synthetic NuSTAR observations on a grid of inner disk radius, disk ionisation, and reflection fraction, we find that the disappearance of the reflection features can be explained by either increased disk ionisation ($\log ξ\geq 4.1$) or a much decreased reflection fraction. A changing disk truncation alone, however, cannot account for the lack of reprocessed Fe-K$α$ emission. The required increase in ionisation parameter could occur if the inner accretion flow evaporates from a thin disk into a geometrically thicker flow, such as the commonly assumed formation of an radiatively inefficient accretion flow at lower mass accretion rates. △ Less

Submitted 10 February, 2020; originally announced February 2020.

Comments: Accepted for publication in MNRAS

NICER X-ray Observations of Seven Nearby Rotation-Powered Millisecond Pulsars

Authors: S. Guillot, M. Kerr, P. S. Ray, S. Bogdanov, S. Ransom, J. S. Deneva, Z. Arzoumanian, P. Bult, D. Chakrabarty, K. C. Gendreau, W. C. G. Ho, G. K. Jaisawal, C. Malacaria, M. C. Miller, T. E. Strohmayer, M. T. Wolff, K. S. Wood, N. A. Webb, L. Guillemot, I. Cognard, G. Theureau

Abstract: NICER observed several rotation-powered millisecond pulsars to search for or confirm the presence of X-ray pulsations. When broad and sine-like, these pulsations may indicate thermal emission from hot polar caps at the magnetic poles on the neutron star surface. We report confident detections ($\ge4.7σ$ after background filtering) of X-ray pulsations for five of the seven pulsars in our target sam… ▽ More NICER observed several rotation-powered millisecond pulsars to search for or confirm the presence of X-ray pulsations. When broad and sine-like, these pulsations may indicate thermal emission from hot polar caps at the magnetic poles on the neutron star surface. We report confident detections ($\ge4.7σ$ after background filtering) of X-ray pulsations for five of the seven pulsars in our target sample: PSR J0614-3329, PSR J0636+5129, PSR J0751+1807, PSR J1012+5307, and PSR J2241-5236, while PSR J1552+5437 and PSR J1744-1134 remain undetected. Of those, only PSR J0751+1807 and PSR J1012+5307 had pulsations previously detected at the 1.7$σ$ and almost 3$σ$ confidence levels, respectively, in XMM-Newton data. All detected sources exhibit broad sine-like pulses, which are indicative of surface thermal radiation. As such, these MSPs are promising targets for future X-ray observations aimed at constraining the neutron star mass-radius relation and the dense matter equation of state using detailed pulse profile modeling. Furthermore, we find that three of the detected millisecond pulsars exhibit a significant phase offset between their X-ray and radio pulses. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: 25 pages, 11 tables, 4 figures. In press in The Astrophysical Journal Letters

Journal ref: ApJL, 887, L27 (2019)

Constraining the Neutron Star Mass-Radius Relation and Dense Matter Equation of State with NICER. II. Emission from Hot Spots on a Rapidly Rotating Neutron Star

Authors: Slavko Bogdanov, Frederick K. Lamb, Simin Mahmoodifar, M. Coleman Miller, Sharon M. Morsink, Thomas E. Riley, Tod E. Strohmayer, Anna L. Watts, Alexander J. Dittmann, Deepto Chakrabarty, Sebastien Guillot, Zaven Arzoumanian, Keith C. Gendreau

Abstract: We describe the model of surface emission from a rapidly rotating neutron star that is applied to Neutron Star Interior Composition Explorer X-ray data of millisecond pulsars in order to statistically constrain the neutron star mass-radius relation and dense matter equation of state. To ensure that the associated calculations are both accurate and precise, we conduct an extensive suite of verifica… ▽ More We describe the model of surface emission from a rapidly rotating neutron star that is applied to Neutron Star Interior Composition Explorer X-ray data of millisecond pulsars in order to statistically constrain the neutron star mass-radius relation and dense matter equation of state. To ensure that the associated calculations are both accurate and precise, we conduct an extensive suite of verification tests between our numerical codes for both the Schwarzschild + Doppler and Oblate Schwarzschild approximations, and compare both approximations against exact numerical calculations. We find superb agreement between the code outputs, as well as in comparison against a set of analytical and semi-analytical calculations, which combined with their speed, demonstrates that the codes are well-suited for large-scale statistical sampling applications. A set of verified, high-precision reference synthetic pulse profiles is provided to the community to facilitate testing of other independently developed codes. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: 31 pages, 14 figures; accepted for publication in the Astrophysical Journal Letters

Constraining the Neutron Star Mass-Radius Relation and Dense Matter Equation of State with NICER. I. The Millisecond Pulsar X-ray Data Set

Authors: Slavko Bogdanov, Sebastien Guillot, Paul S. Ray, Michael T. Wolff, Deepto Chakrabarty, Wynn C. G. Ho, Matthew Kerr, Frederick K. Lamb, Andrea Lommen, Renee M. Ludlam, Reilly Milburn, Sergio Montano, M. Coleman Miller, Michi Baubock, Feryal Ozel, Dimitrios Psaltis, Ronald A. Remillard, Thomas E. Riley, James F. Steiner, Tod E. Strohmayer, Anna L. Watts, Kent S. Wood, Jesse Zeldes, Teruaki Enoto, Takashi Okajima , et al. (5 additional authors not shown)

Abstract: We present the set of deep Neutron Star Interior Composition Explorer (NICER) X-ray timing observations of the nearby rotation-powered millisecond pulsars PSRs J0437-4715, J0030+0451, J1231-1411, and J2124-3358, selected as targets for constraining the mass-radius relation of neutron stars and the dense matter equation of state via modeling of their pulsed thermal X-ray emission. We describe the i… ▽ More We present the set of deep Neutron Star Interior Composition Explorer (NICER) X-ray timing observations of the nearby rotation-powered millisecond pulsars PSRs J0437-4715, J0030+0451, J1231-1411, and J2124-3358, selected as targets for constraining the mass-radius relation of neutron stars and the dense matter equation of state via modeling of their pulsed thermal X-ray emission. We describe the instrument, observations, and data processing/reduction procedures, as well as the series of investigations conducted to ensure that the properties of the data sets are suitable for parameter estimation analyses to produce reliable constraints on the neutron star mass-radius relation and the dense matter equation of state. We find that the long-term timing and flux behavior and the Fourier-domain properties of the event data do not exhibit any anomalies that could adversely affect the intended measurements. From phase-selected spectroscopy, we find that emission from the individual pulse peaks is well described by a single-temperature hydrogen atmosphere spectrum, with the exception of PSR J0437-4715, for which multiple temperatures are required. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: 28 pages, 9 figures; accepted for publication in the Astrophysical Journal Letters

A NICER view of PSR J0030+0451: evidence for a global-scale multipolar magnetic field

Authors: Anna V. Bilous, Anna L. Watts, Alice K. Harding, Thomas E. Riley, Zaven Arzoumanian, Slavko Bogdanov, Keith C. Gendreau, Paul S. Ray, Sebastien Guillot, Wynn C. G. Ho, Deepto Chakrabarty

Abstract: Recent modeling of NICER observations of thermal X-ray pulsations from the surface of the isolated millisecond pulsar PSR J0030+0451 suggests that the hot emitting regions on the pulsar's surface are far from antipodal, which is at odds with the classical assumption that the magnetic field in the pulsar magnetosphere is predominantly that of a centered dipole. Here, we review these results and exa… ▽ More Recent modeling of NICER observations of thermal X-ray pulsations from the surface of the isolated millisecond pulsar PSR J0030+0451 suggests that the hot emitting regions on the pulsar's surface are far from antipodal, which is at odds with the classical assumption that the magnetic field in the pulsar magnetosphere is predominantly that of a centered dipole. Here, we review these results and examine previous attempts to constrain the magnetospheric configuration of PSR J0030+0451. To the best of our knowledge, there is in fact no direct observational evidence that PSR J0030+0451's magnetic field is a centered dipole. Developing models of physically motivated, non-canonical magnetic field configurations and the currents that they can support poses a challenging task. However, such models may have profound implications for many aspects of pulsar research, including pulsar braking, estimates of birth velocities, and interpretations of multi-wavelength magnetospheric emission. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: Appears in ApJ Letters Focus Issue on NICER Constraints on the Dense Matter Equation of State

Journal ref: Bilous et al. 2019, ApJ Letters, 887, L23

A NICER View of PSR J0030+0451: Millisecond Pulsar Parameter Estimation

Authors: Thomas E. Riley, Anna L. Watts, Slavko Bogdanov, Paul S. Ray, Renee M. Ludlam, Sebastien Guillot, Zaven Arzoumanian, Charles L. Baker, Anna V. Bilous, Deepto Chakrabarty, Keith C. Gendreau, Alice K. Harding, Wynn C. G. Ho, James M. Lattimer, Sharon M. Morsink, Tod E. Strohmayer

Abstract: We report on Bayesian parameter estimation of the mass and equatorial radius of the millisecond pulsar PSR J0030$+$0451, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer (NICER) X-ray spectral-timing event data. We perform relativistic ray-tracing of thermal emission from hot regions of the pulsar's surface. We assume two distinct hot regions based on two clear p… ▽ More We report on Bayesian parameter estimation of the mass and equatorial radius of the millisecond pulsar PSR J0030$+$0451, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer (NICER) X-ray spectral-timing event data. We perform relativistic ray-tracing of thermal emission from hot regions of the pulsar's surface. We assume two distinct hot regions based on two clear pulsed components in the phase-folded pulse-profile data; we explore a number of forms (morphologies and topologies) for each hot region, inferring their parameters in addition to the stellar mass and radius. For the family of models considered, the evidence (prior predictive probability of the data) strongly favors a model that permits both hot regions to be located in the same rotational hemisphere. Models wherein both hot regions are assumed to be simply-connected circular single-temperature spots, in particular those where the spots are assumed to be reflection-symmetric with respect to the stellar origin, are strongly disfavored. For the inferred configuration, one hot region subtends an angular extent of only a few degrees (in spherical coordinates with origin at the stellar center) and we are insensitive to other structural details; the second hot region is far more azimuthally extended in the form of a narrow arc, thus requiring a larger number of parameters to describe. The inferred mass $M$ and equatorial radius $R_\mathrm{eq}$ are, respectively, $1.34_{-0.16}^{+0.15}$ M$_{\odot}$ and $12.71_{-1.19}^{+1.14}$ km, whilst the compactness $GM/R_\mathrm{eq}c^2 = 0.156_{-0.010}^{+0.008}$ is more tightly constrained; the credible interval bounds reported here are approximately the $16\%$ and $84\%$ quantiles in marginal posterior mass. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: Appears in ApJ Letters Focus Issue on NICER Constraints on the Dense Matter Equation of State; 76 pages, 24 figures, 7 tables, 8 figure sets (available in the online journal or from the authors)

Journal ref: ApJL, 887, L21 (2019)