Showing 1–50 of 115 results for author: Russell, C

The neutron star mass, distance, and inclination from precision timing of the brilliant millisecond pulsar J0437$-$4715

Authors: Daniel J. Reardon, Matthew Bailes, Ryan M. Shannon, Chris Flynn, Jacob Askew, N. D. Ramesh Bhat, Zu-Cheng Chen, Małgorzata Curyło, Yi Feng, George B. Hobbs, Agastya Kapur, Matthew Kerr, Xiaojin Liu, Richard N. Manchester, Rami Mandow, Saurav Mishra, Christopher J. Russell, Mohsen Shamohammadi, Lei Zhang, Andrew Zic

Abstract: The observation of neutron stars enables the otherwise impossible study of fundamental physical processes. Timing of binary radio pulsars is particularly powerful, as it enables precise characterization of their (three-dimensional) positions and orbits. PSR J0437$-$4715 is an important millisecond pulsar for timing array experiments and is also a primary target for the Neutron Star Interior Compos… ▽ More The observation of neutron stars enables the otherwise impossible study of fundamental physical processes. Timing of binary radio pulsars is particularly powerful, as it enables precise characterization of their (three-dimensional) positions and orbits. PSR J0437$-$4715 is an important millisecond pulsar for timing array experiments and is also a primary target for the Neutron Star Interior Composition ExploreR (NICER). The main aim of the NICER mission is to constrain the neutron star equation of state by inferring the compactness ($M_p/R$) of the star. Direct measurements of the mass $M_p$ from pulsar timing therefore substantially improve constraints on the radius $R$, and the equation of state. Here we use observations spanning 26 years from Murriyang, the 64-m Parkes radio telescope, to improve the timing model for this pulsar. Among the new precise measurements are the pulsar mass $M_p=1.418\pm 0.044$ M$_{\odot}$, distance $D=156.96 \pm 0.11$ pc, and orbital inclination angle $i=137.506 \pm 0.016^\circ$, which can be used to inform the X-ray pulse profile models inferred from NICER observations. We demonstrate that these results are consistent between multiple data sets from the Parkes Pulsar Timing Array (PPTA), each modelled with different noise assumptions. Using the longest available PPTA data set, we measure an apparent second derivative of the pulsar spin frequency and discuss how this can be explained either by kinematic effects due to the proper motion and radial velocity of the pulsar, or excess low-frequency noise such as a gravitational-wave background. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 13 pages, 3 figures, accepted for publication in Astrophysical Journal Letters

Multistructured accretion flow of Sgr A* II: Signatures of a Cool Accretion Disk in Hydrodynamic Simulations of Stellar Winds

Authors: Mayura Balakrishnan, Christopher M. P. Russell, Lia Corrales, Diego Calderón, Jorge Cuadra, Daryl Haggard, Sera Markoff, Joey Neilsen, Michael Nowak, Q. Daniel Wang, Fred Baganoff

Abstract: Hydrodynamic simulations of the stellar winds from Wolf-Rayet stars within the Galactic Center can provide predictions for the X-ray spectrum of supermassive black hole Sgr A*. Herein, we present results from updated smooth particle hydrodynamics simulations, building on the architecture of Cuadra et al. (2015); Russell et al. (2017), finding that a cold gas disk forms around Sgr A* with a simulat… ▽ More Hydrodynamic simulations of the stellar winds from Wolf-Rayet stars within the Galactic Center can provide predictions for the X-ray spectrum of supermassive black hole Sgr A*. Herein, we present results from updated smooth particle hydrodynamics simulations, building on the architecture of Cuadra et al. (2015); Russell et al. (2017), finding that a cold gas disk forms around Sgr A* with a simulation runtime of 3500 years. This result is consistent with previous grid-based simulations, demonstrating that a cold disk can form regardless of numerical method. We examine the plasma scenarios arising from an environment with and without this cold disk, by generating synthetic spectra for comparison to the quiescent Fe K alpha Sgr A* spectrum from Chandra HETG-S, taken through the Chandra X-ray Visionary Program. We find that current and future X-ray missions are unlikely to distinguish between the kinematic signatures in the plasma in these two scenarios. Nonetheless, the stellar wind plasma model presents a good fit to the dispersed Chandra spectra within 1.5" of Sgr A*. We compare our results to the Radiatively Inefficient Accretion Flow (RIAF) model fit to the HETG-S spectrum presented in Paper I and find that the Bayesian model evidence does not strongly favor either model. With 9" angular resolution and high spectral resolution of the X-IFU, NewAthena will offer a clearer differentiation between the RIAF plasma model and hydrodynamic simulations, but only a future X-ray mission with arcsecond resolution will significantly advance our understanding of Sgr A*'s accretion flow in X-rays. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: 14 pages, 6 figures. Accepted for publication in ApJ

Multistructured accretion flow of Sgr A* I: Examination of a RIAF model

Authors: Mayura Balakrishnan, Lia Corrales, Sera Markoff, Michael Nowak, Daryl Haggard, Q. Daniel Wang, Joey Neilsen, Christopher M. P. Russell, Diego Calderón, Jorge Cuadra, Fred Baganoff

Abstract: The extreme low-luminosity supermassive black hole Sgr A* provides a unique laboratory in which to test radiatively inefficient accretion flow (RIAF) models. Previous fits to the quiescent Chandra ACIS-S spectrum found a RIAF model with an equal inflow-outflow balance works well. In this work, we apply the RIAF model to the Chandra HETG-S spectrum obtained through the Chandra X-ray Visionary Progr… ▽ More The extreme low-luminosity supermassive black hole Sgr A* provides a unique laboratory in which to test radiatively inefficient accretion flow (RIAF) models. Previous fits to the quiescent Chandra ACIS-S spectrum found a RIAF model with an equal inflow-outflow balance works well. In this work, we apply the RIAF model to the Chandra HETG-S spectrum obtained through the Chandra X-ray Visionary Program, which displays features suggestive of temperature and velocity structures within the plasma. A comprehensive forward model analysis accounting for the accretion flow geometry and HETG-S instrumental effects is required for a full interpretation of the quiescent Chandra HETG-S spectrum. We present a RIAF model that takes these effects into account. Our fits to the high-resolution gratings spectrum indicate an inflow balanced by an outflow ($s \sim 1$) alongside a temperature profile that appears shallower than what would be expected from a gravitational potential following $1/r$. The data require that the abundance of Iron relative to solar is $Z_{Fe} < 0.32 Z_\odot$ (90\% credible interval), much lower than the $2~Z_\odot$ metallicity measured in nearby late-type giants. While future missions like NewAthena will provide higher spectral resolution, source separation will continue to be a problem. Leveraging Chandra's unparalleled spatial resolution, which is not expected to be surpassed for decades, remains essential for detailed investigations of the densely populated Galactic Center in X-rays. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: 11 pages, 5 figures, 1 table. Accepted for publication in ApJ

Earth's Alfvén wings driven by the April 2023 Coronal Mass Ejection

Authors: Li-Jen Chen, Daniel Gershman, Brandon Burkholder, Yuxi Chen, Menelaos Sarantos, Lan Jian, James Drake, Chuanfei Dong, Harsha Gurram, Jason Shuster, Daniel Graham, Olivier Le Contel, Steven Schwartz, Stephen Fuselier, Hadi Madanian, Craig Pollock, Haoming Liang, Matthew Argall, Richard Denton, Rachel Rice, Jason Beedle, Kevin Genestreti, Akhtar Ardakani, Adam Stanier, Ari Le , et al. (11 additional authors not shown)

Abstract: We report a rare regime of Earth's magnetosphere interaction with sub-Alfvénic solar wind in which the windsock-like magnetosphere transforms into one with Alfvén wings. In the magnetic cloud of a Coronal Mass Ejection (CME) on April 24, 2023, NASA's Magnetospheric Multiscale mission distinguishes the following features: (1) unshocked and accelerated cold CME plasma coming directly against Earth's… ▽ More We report a rare regime of Earth's magnetosphere interaction with sub-Alfvénic solar wind in which the windsock-like magnetosphere transforms into one with Alfvén wings. In the magnetic cloud of a Coronal Mass Ejection (CME) on April 24, 2023, NASA's Magnetospheric Multiscale mission distinguishes the following features: (1) unshocked and accelerated cold CME plasma coming directly against Earth's dayside magnetosphere; (2) dynamical wing filaments representing new channels of magnetic connection between the magnetosphere and foot points of the Sun's erupted flux rope; (3) cold CME ions observed with energized counter-streaming electrons, evidence of CME plasma captured due to reconnection between magnetic-cloud and Alfvén-wing field lines. The reported measurements advance our knowledge of CME interaction with planetary magnetospheres, and open new opportunities to understand how sub-Alfvénic plasma flows impact astrophysical bodies such as Mercury, moons of Jupiter, and exoplanets close to their host stars. △ Less

Submitted 3 May, 2024; v1 submitted 12 February, 2024; originally announced February 2024.

Comments: 14 pages, including 4 figures, Under review in Geophys. Res. Lett

X-ray plasma flow and turbulence in the colliding winds of WR140

Authors: Asca Miyamoto, Yasuharu Sugawara, Yoshitomo Maeda, Manabu Ishida, Kenji Hamaguchi, Christopher M. P. Russell, Anthony F. J. Moffat

Abstract: We analyse $\textit{XMM-Newton}$ RGS spectra of Wolf-Rayet (WR) 140, an archetype long-period eccentric WR+O colliding wind binary. We evaluate the spectra of O and Fe emission lines and find that the plasmas emitting these lines have the largest approaching velocities with the largest velocity dispersions between phases 0.935 and 0.968 where the inferior conjunction of the O star occurs. This beh… ▽ More We analyse $\textit{XMM-Newton}$ RGS spectra of Wolf-Rayet (WR) 140, an archetype long-period eccentric WR+O colliding wind binary. We evaluate the spectra of O and Fe emission lines and find that the plasmas emitting these lines have the largest approaching velocities with the largest velocity dispersions between phases 0.935 and 0.968 where the inferior conjunction of the O star occurs. This behaviour is the same as that of the Ne line-emission plasma presented in our previous paper. We perform diagnosis of electron number density $n_{\rm e}$ using He-like triplet lines of O and Ne-like Fe-L lines. The former results in a conservative upper limit of $n_{\rm e} \lesssim 10^{10}$-10$^{12}$ cm$^{-3}$ on the O line-emission site, while the latter can not impose any constraint on the Fe line-emission site because of statistical limitations. We calculate the line-of-sight velocity and its dispersion separately along the shock cone. By comparing the observed and calculated line-of-sight velocities, we update the distance of the Ne line-emission site from the stagnation point. By assuming radiative cooling of the Ne line-emission plasma using the observed temperature and the local stellar wind density, we estimate the line-emission site extends along the shock cone by at most $\pm$58 per cent (phase 0.816) of the distance from the stagnation point. In this framework, excess of the observed velocity dispersion over the calculated one is ascribed to turbulence in the hot-shocked plasma at earlier orbital phases of 0.816, 0.912, and 0.935, with the largest velocity dispersion of 340-630 km s$^{-1}$ at phase 0.912. △ Less

Submitted 1 January, 2024; originally announced January 2024.

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

A First Look with JWST Aperture Masking Interferometry (AMI): Resolving Circumstellar Dust around the Wolf-Rayet Binary WR 137 beyond the Rayleigh Limit

Authors: Ryan M. Lau, Matthew J. Hankins, Joel Sanchez-Bermudez, Deepashri Thatte, Anthony Soulain, Rachel A. Cooper, Anand Sivaramakrishnan, Michael F. Corcoran, Alexandra Z. Greenbaum, Theodore R. Gull, Yinuo Han, Olivia C. Jones, Thomas Madura, Anthony F. J. Moffat, Mark R. Morris, Takashi Onaka, Christopher M. P. Russell, Noel D. Richardson, Nathan Smith, Peter Tuthill, Kevin Volk, Gerd Weigelt, Peredur M. Williams

Abstract: We present infrared aperture masking interferometry (AMI) observations of newly formed dust from the colliding winds of the massive binary system Wolf-Rayet (WR) 137 with JWST using the Near Infrared Imager and Slitless Spectrograph (NIRISS). NIRISS AMI observations of WR 137 and a point-spread-function calibrator star, HD~228337, were taken using the F380M and F480M filters in 2022 July and Augus… ▽ More We present infrared aperture masking interferometry (AMI) observations of newly formed dust from the colliding winds of the massive binary system Wolf-Rayet (WR) 137 with JWST using the Near Infrared Imager and Slitless Spectrograph (NIRISS). NIRISS AMI observations of WR 137 and a point-spread-function calibrator star, HD~228337, were taken using the F380M and F480M filters in 2022 July and August as part of the Director's Discretionary Early Release Science (DD-ERS) program 1349. Interferometric observables (squared visibilities and closure phases) from the WR 137 "interferogram" were extracted and calibrated using three independent software tools: ImPlaneIA, AMICAL, and SAMpip. The analysis of the calibrated observables yielded consistent values except for slightly discrepant closure phases measured by ImPlaneIA. Based on all three sets of calibrated observables, images were reconstructed using three independent software tools: BSMEM, IRBis, and SQUEEZE. All reconstructed image combinations generated consistent images in both F380M and F480M filters. The reconstructed images of WR 137 reveal a bright central core with a $\sim300$ mas linear filament extending to the northwest. A geometric colliding-wind model with dust production constrained to the orbital plane of the binary system and enhanced as the system approaches periapsis provided a general agreement with the interferometric observables and reconstructed images. Based on a colliding-wind dust condensation analysis, we suggest that dust formation within the orbital plane of WR 137 is induced by enhanced equatorial mass-loss from the rapidly rotating O9 companion star, whose axis of rotation is aligned with that of the orbit. △ Less

Submitted 22 December, 2023; v1 submitted 27 November, 2023; originally announced November 2023.

Comments: 18 pages, 8 figures, Accepted for publication in ApJ. Updated plotting error in Fig. 2

Comparing recent PTA results on the nanohertz stochastic gravitational wave background

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

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

Submitted 1 September, 2023; originally announced September 2023.

Comments: 21 pages, 9 figures, submitted to ApJ

The Parkes Pulsar Timing Array Third Data Release

Authors: Andrew Zic, Daniel J. Reardon, Agastya Kapur, George Hobbs, Rami Mandow, Małgorzata Curyło, Ryan M. Shannon, Jacob Askew, Matthew Bailes, N. D. Ramesh Bhat, Andrew Cameron, Zu-Cheng Chen, Shi Dai, Valentina Di Marco, Yi Feng, Matthew Kerr, Atharva Kulkarni, Marcus E. Lower, Rui Luo, Richard N. Manchester, Matthew T. Miles, Rowina S. Nathan, Stefan Osłowski, Axl F. Rogers, Christopher J. Russell , et al. (9 additional authors not shown)

Abstract: We present the third data release from the Parkes Pulsar Timing Array (PPTA) project. The release contains observations of 32 pulsars obtained using the 64-m Parkes "Murriyang" radio telescope. The data span is up to 18 years with a typical cadence of 3 weeks. This data release is formed by combining an updated version of our second data release with $\sim 3$ years of more recent data primarily ob… ▽ More We present the third data release from the Parkes Pulsar Timing Array (PPTA) project. The release contains observations of 32 pulsars obtained using the 64-m Parkes "Murriyang" radio telescope. The data span is up to 18 years with a typical cadence of 3 weeks. This data release is formed by combining an updated version of our second data release with $\sim 3$ years of more recent data primarily obtained using an ultra-wide-bandwidth receiver system that operates between 704 and 4032 MHz. We provide calibrated pulse profiles, flux-density dynamic spectra, pulse times of arrival, and initial pulsar timing models. We describe methods for processing such wide-bandwidth observations, and compare this data release with our previous release. △ Less

Submitted 17 October, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 15 pages, 6 figures. Accepted for publication in PASA

The gravitational-wave background null hypothesis: Characterizing noise in millisecond pulsar arrival times with the Parkes Pulsar Timing Array

Authors: Daniel J. Reardon, Andrew Zic, Ryan M. Shannon, Valentina Di Marco, George B. Hobbs, Agastya Kapur, Marcus E. Lower, Rami Mandow, Hannah Middleton, Matthew T. Miles, Axl F. Rogers, Jacob Askew, Matthew Bailes, N. D. Ramesh Bhat, Andrew Cameron, Matthew Kerr, Atharva Kulkarni, Richard N. Manchester, Rowina S. Nathan, Christopher J. Russell, Stefan Osłowski, Xing-Jiang Zhu

Abstract: The noise in millisecond pulsar (MSP) timing data can include contributions from observing instruments, the interstellar medium, the solar wind, solar system ephemeris errors, and the pulsars themselves. The noise environment must be accurately characterized in order to form the null hypothesis from which signal models can be compared, including the signature induced by nanohertz-frequency gravita… ▽ More The noise in millisecond pulsar (MSP) timing data can include contributions from observing instruments, the interstellar medium, the solar wind, solar system ephemeris errors, and the pulsars themselves. The noise environment must be accurately characterized in order to form the null hypothesis from which signal models can be compared, including the signature induced by nanohertz-frequency gravitational waves (GWs). Here we describe the noise models developed for each of the MSPs in the Parkes Pulsar Timing Array (PPTA) third data release, which have been used as the basis of a search for the isotropic stochastic GW background. We model pulsar spin noise, dispersion measure variations, scattering variations, events in the pulsar magnetospheres, solar wind variability, and instrumental effects. We also search for new timing model parameters and detected Shapiro delays in PSR~J0614$-$3329 and PSR~J1902$-$5105. The noise and timing models are validated by testing the normalized and whitened timing residuals for Gaussianity and residual correlations with time. We demonstrate that the choice of noise models significantly affects the inferred properties of a common-spectrum process. Using our detailed models, the recovered common-spectrum noise in the PPTA is consistent with a power law with a spectral index of $γ=13/3$, the value predicted for a stochastic GW background from a population of supermassive black hole binaries driven solely by GW emission. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 18 pages, 10 figures. Accepted for publication in ApJL

Search for an isotropic gravitational-wave background with the Parkes Pulsar Timing Array

Authors: Daniel J. Reardon, Andrew Zic, Ryan M. Shannon, George B. Hobbs, Matthew Bailes, Valentina Di Marco, Agastya Kapur, Axl F. Rogers, Eric Thrane, Jacob Askew, N. D. Ramesh Bhat, Andrew Cameron, Małgorzata Curyło, William A. Coles, Shi Dai, Boris Goncharov, Matthew Kerr, Atharva Kulkarni, Yuri Levin, Marcus E. Lower, Richard N. Manchester, Rami Mandow, Matthew T. Miles, Rowina S. Nathan, Stefan Osłowski , et al. (4 additional authors not shown)

Abstract: Pulsar timing arrays aim to detect nanohertz-frequency gravitational waves (GWs). A background of GWs modulates pulsar arrival times and manifests as a stochastic process, common to all pulsars, with a signature spatial correlation. Here we describe a search for an isotropic stochastic gravitational-wave background (GWB) using observations of 30 millisecond pulsars from the third data release of t… ▽ More Pulsar timing arrays aim to detect nanohertz-frequency gravitational waves (GWs). A background of GWs modulates pulsar arrival times and manifests as a stochastic process, common to all pulsars, with a signature spatial correlation. Here we describe a search for an isotropic stochastic gravitational-wave background (GWB) using observations of 30 millisecond pulsars from the third data release of the Parkes Pulsar Timing Array (PPTA), which spans 18 years. Using current Bayesian inference techniques we recover and characterize a common-spectrum noise process. Represented as a strain spectrum $h_c = A(f/1 {\rm yr}^{-1})^α$, we measure $A=3.1^{+1.3}_{-0.9} \times 10^{-15}$ and $α=-0.45 \pm 0.20$ respectively (median and 68% credible interval). For a spectral index of $α=-2/3$, corresponding to an isotropic background of GWs radiated by inspiraling supermassive black hole binaries, we recover an amplitude of $A=2.04^{+0.25}_{-0.22} \times 10^{-15}$. However, we demonstrate that the apparent signal strength is time-dependent, as the first half of our data set can be used to place an upper limit on $A$ that is in tension with the inferred common-spectrum amplitude using the complete data set. We search for spatial correlations in the observations by hierarchically analyzing individual pulsar pairs, which also allows for significance validation through randomizing pulsar positions on the sky. For a process with $α=-2/3$, we measure spatial correlations consistent with a GWB, with an estimated false-alarm probability of $p \lesssim 0.02$ (approx. $2σ$). The long timing baselines of the PPTA and the access to southern pulsars will continue to play an important role in the International Pulsar Timing Array. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 19 pages, 10 figures, Accepted for publication in ApJL

Searching for continuous Gravitational Waves in the second data release of the International Pulsar Timing Array

Authors: M. Falxa, S. Babak, P. T. Baker, B. Bécsy, A. Chalumeau, S. Chen, Z. Chen, N. J. Cornish, L. Guillemot, J. S. Hazboun, C. M. F. Mingarelli, A. Parthasarathy, A. Petiteau, N. S. Pol, A. Sesana, S. B. Spolaor, S. R. Taylor, G. Theureau, M. Vallisneri, S. J. Vigeland, C. A. Witt, X. Zhu, J. Antoniadis, Z. Arzoumanian, M. Bailes , et al. (102 additional authors not shown)

Abstract: The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evi… ▽ More The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evidence for such signals and set sky averaged 95% upper limits on their amplitude h 95 . The most sensitive frequency is 10nHz with h 95 = 9.1 10-15 . We achieved the best upper limit to date at low and high frequencies of the PTA band thanks to improved effective cadence of observations. In our analysis, we have taken into account the recently discovered common red noise process, which has an impact at low frequencies. We also find that the peculiar noise features present in some pulsars data must be taken into account to reduce the false alarm. We show that using custom noise models is essential in searching for continuous gravitational wave signals and setting the upper limit. △ Less

Submitted 19 March, 2023; originally announced March 2023.

3D MHD models of the centrifugal magnetosphere from a massive star with an oblique dipole field

Authors: Asif ud-Doula, Stanley P. Owocki, Christopher Russell, Marc Gagne, Simon Daley-Yates

Abstract: We present results from new self-consistent 3D MHD simulations of the magnetospheres from massive stars with a dipole magnetic axis that has a non-zero obliquity angle ($β$) to the star's rotation axis. As an initial direct application, we compare the global structure of co-rotating disks for nearly aligned ($β=5^o$) versus half-oblique ($β=45^o$) models, both with moderately rapid rotation (… ▽ More We present results from new self-consistent 3D MHD simulations of the magnetospheres from massive stars with a dipole magnetic axis that has a non-zero obliquity angle ($β$) to the star's rotation axis. As an initial direct application, we compare the global structure of co-rotating disks for nearly aligned ($β=5^o$) versus half-oblique ($β=45^o$) models, both with moderately rapid rotation ($\sim$ 0.5 critical). We find that accumulation surfaces broadly resemble the forms predicted by the analytic Rigidly Rotating Magnetosphere (RRM) model, but the mass buildup to near the critical level for centrifugal breakout against magnetic confinement distorts the field from the imposed initial dipole. This leads to an associated warping of the accumulation surface toward the rotational equator, with the highest density concentrated in {\em wings} centered on the intersection between the magnetic and rotational equators. These MHD models can be used to synthesize rotational modulation of photometric absorption and H$α$ emission for a direct comparison with observations. △ Less

Submitted 27 January, 2023; originally announced January 2023.

Comments: Accepted for publication in MNRAS

Wideband timing of the Parkes Pulsar Timing Array UWL data

Authors: Małgorzata Curyło, Timothy T. Pennucci, Matthew Bailes, N. D. Ramesh Bhat, Andrew D. Cameron, Shi Dai, George Hobbs, Agastya Kapur, Richard N. Manchester, Rami Mandow, Matthew T. Miles, Christopher J. Russell, Daniel J. Reardon, Ryan M. Shannon, Renée Spiewak, Andrew Zic, Xing-Jiang Zhu

Abstract: In 2018 an ultra-wide-bandwidth low-frequency (UWL) receiver was installed on the 64-m Parkes Radio Telescope enabling observations with an instantaneous frequency coverage from 704 to 4032 MHz. Here, we present the analysis of a three-year data set of 35 millisecond pulsars observed with the UWL by the Parkes Pulsar Timing Array (PPTA), using wideband timing methods. The two key differences compa… ▽ More In 2018 an ultra-wide-bandwidth low-frequency (UWL) receiver was installed on the 64-m Parkes Radio Telescope enabling observations with an instantaneous frequency coverage from 704 to 4032 MHz. Here, we present the analysis of a three-year data set of 35 millisecond pulsars observed with the UWL by the Parkes Pulsar Timing Array (PPTA), using wideband timing methods. The two key differences compared to typical narrow-band methods are, firstly, generation of two-dimensional templates accounting for pulse shape evolution with frequency and, secondly, simultaneous measurements of the pulse time-of-arrival (ToA) and dispersion measure (DM). This is the first time that wideband timing has been applied to a uniform data set collected with a single large-fractional bandwidth receiver, for which such techniques were originally developed. As a result of our study, we present a set of profile evolution models and new timing solutions including initial noise analysis. Precision of our ToA and DM measurements is in the range of 0.005 $-$ 2.08 $μ$s and (0.043$-$14.24)$\times10^{-4}$ cm$^{-3}$ pc, respectively, with 94% of the pulsars achieving a median ToA uncertainty of less than 1 $μ$s. △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: 29 pages, 12 figures, Accepted for publication in ApJ

The long-term spectral changes of eta Carinae: are they caused by a dissipating occulter as indicated by CMFGEN models?

Authors: A. Damineli, D. J. Hillier, F. Navarete, A. F. J. Moffat, G. Weigelt, M. F. Corcoran, T. R. Gull, N. D. Richardson, T. P. Ho, T. I. Madura, D. Espinoza-Galeas, H. Hartman, P. Morris, C. S. Pickett, I. R. Stevens, C. M. P. Russell, K. Hamaguchi, F. J. Jablonski, M. Teodoro, P. McGee, P. Cacella, B. Heathcote, K. Harrison, M. Johnston, T. Bohlsen , et al. (1 additional authors not shown)

Abstract: Eta Carinae ($η$\,Car) exhibits a unique set of P Cygni profiles with both broad and narrow components. Over many decades, the spectrum has changed -- there has been an increase in observed continuum fluxes and a decrease in FeII and HI emission line equivalent widths. The spectrum is evolving towards that of a P Cygni star such as P~Cygni itself and HDE~316285. The spectral evolution has been att… ▽ More Eta Carinae ($η$\,Car) exhibits a unique set of P Cygni profiles with both broad and narrow components. Over many decades, the spectrum has changed -- there has been an increase in observed continuum fluxes and a decrease in FeII and HI emission line equivalent widths. The spectrum is evolving towards that of a P Cygni star such as P~Cygni itself and HDE~316285. The spectral evolution has been attributed to intrinsic variations such as a decrease in the mass-loss rate of the primary star or differential evolution in a latitudinal-dependent stellar wind. However intrinsic wind changes conflict with three observational results: the steady long-term bolometric luminosity; the repeating X-ray light curve over the binary period; and the constancy of the dust-scattered spectrum from the Homunculus. We extend previous work that showed a secular strengthening of P~Cygni absorptions by adding more orbital cycles to overcome temporary instabilities and by examining more atomic transitions. {\sc cmfgen} modeling of the primary wind shows that a time-decreasing mass-loss rate is not the best explanation for the observations. However, models with a `small' dissipating absorber in our line-of-site can explain both the increase in brightness and changes in the emission and P Cygni absorption profiles. If the spectral evolution is caused by the dissipating circumstellar medium, and not by intrinsic changes in the binary, the dynamical timescale to recover from the Great Eruption is much less than a century, different from previous suggestions. △ Less

Submitted 7 July, 2023; v1 submitted 2 November, 2022; originally announced November 2022.

Comments: 71 pages, 8 figures, 4 long tables, to appear on ApJ

Nested Dust Shells around the Wolf-Rayet Binary WR 140 observed with JWST

Authors: Ryan M. Lau, Matthew J. Hankins, Yinuo Han, Ioannis Argyriou, Michael F. Corcoran, Jan J. Eldridge, Izumi Endo, Ori D. Fox, Macarena Garcia Marin, Theodore R. Gull, Olivia C. Jones, Kenji Hamaguchi, Astrid Lamberts, David R. Law, Thomas Madura, Sergey V. Marchenko, Hideo Matsuhara, Anthony F. J. Moffat, Mark R. Morris, Patrick W. Morris, Takashi Onaka, Michael E. Ressler, Noel D. Richardson, Christopher M. P. Russell, Joel Sanchez-Bermudez , et al. (7 additional authors not shown)

Abstract: Massive colliding-wind binaries that host a Wolf-Rayet (WR) star present a potentially important source of dust and chemical enrichment in the interstellar medium (ISM). However, the chemical composition and survival of dust formed from such systems is not well understood. The carbon-rich WR (WC) binary WR~140 presents an ideal astrophysical laboratory for investigating these questions given its w… ▽ More Massive colliding-wind binaries that host a Wolf-Rayet (WR) star present a potentially important source of dust and chemical enrichment in the interstellar medium (ISM). However, the chemical composition and survival of dust formed from such systems is not well understood. The carbon-rich WR (WC) binary WR~140 presents an ideal astrophysical laboratory for investigating these questions given its well-defined orbital period and predictable dust-formation episodes every 7.93 years around periastron passage. We present observations from our Early Release Science program (ERS1349) with the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) Medium-Resolution Spectrometer (MRS) and Imager that reveal the spectral and spatial signatures of nested circumstellar dust shells around WR~140. MIRI MRS spectroscopy of the second dust shell and Imager detections of over 17 shells formed throughout the past $\gtrsim130$ years confirm the survival of carbonaceous dust grains from WR~140 that are likely carriers of "unidentified infrared" (UIR)-band features at 6.4 and 7.7 $μ$m. The observations indicate that dust-forming WC binaries can enrich the ISM with organic compounds and carbonaceous dust. △ Less

Submitted 12 October, 2022; originally announced October 2022.

Comments: Published in Nature Astronomy on Oct 12, 2022; 21 pages, 5 figures, 2 tables

Journal ref: Lau, R.M., Hankins, M.J., Han, Y. et al. Nested dust shells around the Wolf-Rayet binary WR 140 observed with JWST. Nat Astron (2022)

Constraining ultralight vector dark matter with the Parkes Pulsar Timing Array second data release

Authors: Yu-Mei Wu, Zu-Cheng Chen, Qing-Guo Huang, Xingjiang Zhu, N. D. Ramesh Bhat, Yi Feng, George Hobbs, Richard N. Manchester, Christopher J. Russell, R. M. Shannon

Abstract: Composed of ultralight bosons, fuzzy dark matter provides an intriguing solution to challenges that the standard cold dark matter model encounters on sub-galactic scales. The ultralight dark matter with mass $m\sim10^{-23} \rm{eV}$ will induce a periodic oscillation in gravitational potentials with a frequency in the nanohertz band, leading to observable effects in the arrival times of radio pulse… ▽ More Composed of ultralight bosons, fuzzy dark matter provides an intriguing solution to challenges that the standard cold dark matter model encounters on sub-galactic scales. The ultralight dark matter with mass $m\sim10^{-23} \rm{eV}$ will induce a periodic oscillation in gravitational potentials with a frequency in the nanohertz band, leading to observable effects in the arrival times of radio pulses from pulsars. Unlike scalar dark matter, pulsar timing signals induced by the vector dark matter are dependent on the oscillation direction of the vector fields. In this work, we search for ultralight vector dark matter in the mass range of $[2\times 10^{-24}, 2\times 10^{-22}]{\rm{eV}}$ through its gravitational effect in the Parkes Pulsar Timing Array (PPTA) second data release. Since no statistically significant detection is made, we place $95\%$ upper limits on the local dark matter density as $ρ_{\rm{\tiny{VF}}} \lesssim 5{\rm{GeV/cm^{3}}}$ for $m\lesssim 10^{-23}{\rm{eV}}$. As no preferred direction is found for the vector dark matter, these constraints are comparable to those given by the scalar dark matter search with an earlier 12-year data set of PPTA. △ Less

Submitted 7 October, 2022; originally announced October 2022.

Comments: 9 pages, 3 figures, 2 tables; accepted for publication as a Letter in Phys. Rev. D

Thin current sheet behind the dipolarization front

Authors: Nakamura, R., Baumjohann, W., Nakamura, T. K. M., Panov, E., V., Schmid, D., Varsani, A., S. Apatenkov, V. A. Sergeev, J. Birn, T. Nagai, C. Gabrielse, M. Andre, J. L. Burch, C. Carr, I. S Dandouras, C. P. Escoubet, A, N. Fazakerley , et al. (4 additional authors not shown)

Abstract: We report a unique conjugate observation of fast flows and associated current sheet disturbances in the near-Earth magnetotail by MMS (Magnetospheric Multiscale) and Cluster preceding a positive bay onset of a small substorm at ~14:10 UT, Sep. 8, 2018. MMS and Cluster were located both at X ~-14 RE. A dipolarization front (DF) of a localized fast flow was detected by Cluster and MMS, separated in… ▽ More We report a unique conjugate observation of fast flows and associated current sheet disturbances in the near-Earth magnetotail by MMS (Magnetospheric Multiscale) and Cluster preceding a positive bay onset of a small substorm at ~14:10 UT, Sep. 8, 2018. MMS and Cluster were located both at X ~-14 RE. A dipolarization front (DF) of a localized fast flow was detected by Cluster and MMS, separated in the dawn-dusk direction by ~4 RE, almost simultaneously. Adiabatic electron acceleration signatures revealed from comparison of the energy spectra confirm that both spacecraft encounter the same DF. We analyzed the change in the current sheet structure based on multi-scale multi-point data analysis. The current sheet thickened during the passage of DF, yet, temporally thinned subsequently associated with another flow enhancement centered more on the dawnward side of the initial flow. MMS and Cluster observed intense perpendicular and parallel current in the off-equatorial region mainly during this interval of the current sheet thinning. Maximum field-aligned currents both at MMS and Cluster are directed tailward. Detailed analysis of MMS data showed that the intense field-aligned currents consisted of multiple small-scale intense current layers accompanied by enhanced Hall-currents in the dawn-dusk flow-shear region. We suggest that the current sheet thinning is related to the flow bouncing process and/or to the expansion/activation of reconnection. Based on these mesoscale and small-scale multipoint observations, 3D evolution of the flow and current-sheet disturbances was inferred preceding the development of a substorm current wedge. △ Less

Submitted 26 August, 2022; originally announced August 2022.

Journal ref: Journal of Geophysical Research: Space Physics, 126, e2021JA029518, 2021

The Renaissance of Main Belt Asteroid Science

Authors: Simone Marchi, Carol A. Raymond, Christopher T. Russell

Abstract: The NASA Dawn spacecraft took off from Cape Canaveral in September 2007 atop a Delta II rocket starting an ambitious journey to Vesta and Ceres, the two most massive worlds in the largest reservoir of asteroids in the Solar System, the Main Belt. Prior to the Dawn launch, Earth-bound observations of Vesta and Ceres revealed intriguing features--from Vesta's rugged shape to Ceres' tenuous water exo… ▽ More The NASA Dawn spacecraft took off from Cape Canaveral in September 2007 atop a Delta II rocket starting an ambitious journey to Vesta and Ceres, the two most massive worlds in the largest reservoir of asteroids in the Solar System, the Main Belt. Prior to the Dawn launch, Earth-bound observations of Vesta and Ceres revealed intriguing features--from Vesta's rugged shape to Ceres' tenuous water exosphere--, but these objects remained fuzzy speckles of light even through the lenses of the most powerful telescopes. With Dawn's exploration of Vesta (2011-2012) and Ceres (2015-2018) these two worlds came into focus. Breath-taking details emerged of how large collisions sculpted Vesta liberating massive amounts of material in the inner Main Belt, providing the source of an important family of meteorites recovered on Earth. Ceres' complex geology, which may rival that of the Earth and Mars, unveiled recent cryovolcanic activity. This book is dedicated to these highlights, and many more discoveries of the Dawn mission. By the time Dawn completed its mission in 2018, our understanding of the formation of the Solar System had greatly evolved thanks to new theoretical models and to a new trove of meteorite geochemical data, and Dawn observations of Vesta and Ceres provide new, vital constraints to synergistically interpret models and data. The broader implications of the Dawn legacy are presented in a series of dedicated chapters. The editors hope this book will serve as a solid reference for the younger generations as well as for more seasoned researchers to successfully pursue future exploration of the Main Belt. We certainly have learned a lot thanks to Dawn, and yet we know that we have barely scratched the surface of what Main Belt asteroids can tell us about the dawn of our Solar System. △ Less

Submitted 16 August, 2022; originally announced August 2022.

Comments: This is the Epilogue of the book

Evaluating the prevalence of spurious correlations in pulsar timing array datasets

Authors: Andrew Zic, George Hobbs, R. M. Shannon, Daniel Reardon, Boris Goncharov, N. D. Ramesh Bhat, Andrew Cameron, Shi Dai, J. R. Dawson, Matthew Kerr, R. N. Manchester, Rami Mandow, Tommy Marshman, Christopher J. Russell, Nithyanandan Thyagarajan, X. -J. Zhu

Abstract: Pulsar timing array collaborations have recently reported evidence for a noise process with a common spectrum among the millisecond pulsars in the arrays. The spectral properties of this common-noise process are consistent with expectations for an isotropic gravitational-wave background (GWB) from inspiralling supermassive black-hole binaries. However, recent simulation analyses based on Parkes Pu… ▽ More Pulsar timing array collaborations have recently reported evidence for a noise process with a common spectrum among the millisecond pulsars in the arrays. The spectral properties of this common-noise process are consistent with expectations for an isotropic gravitational-wave background (GWB) from inspiralling supermassive black-hole binaries. However, recent simulation analyses based on Parkes Pulsar Timing Array data indicate that such a detection may arise spuriously. In this paper, we use simulated pulsar timing array datasets to further test the robustness of the inference methods for spectral and spatial correlations from a GWB. Expanding on our previous results, we find strong support (Bayes factors exceeding $10^5$) for the presence of a common-spectrum noise process in datasets where no common process is present, under a wide range of timing noise prescriptions per pulsar. We show that these results are highly sensitive to the choice of Bayesian priors on timing noise parameters, with priors that more closely match the injected distributions of timing noise parameters resulting in diminished support for a common-spectrum noise process. These results emphasize shortcomings in current methods for inferring the presence of a common-spectrum process, and imply that the detection of a common process is not a reliable precursor to detection of the GWB. Future searches for the nanohertz GWB should remain focussed on detecting spatial correlations, and make use of more tailored specifications for a common-spectrum noise process. △ Less

Submitted 25 July, 2022; originally announced July 2022.

Comments: 12 pages, 12 figures. Accepted for publication in MNRAS

NICER X-ray Observations of Eta Carinae During its Most Recent Periastron Passage

Authors: David Espinoza-Galeas, Michael Francis Corcoran, Kenji Hamaguchi, Christopher M. P. Russell, Theodore R. Gull, Anthony Moffat, Noel D. Richardson, Gerd Weigelt, D. John Hillier, Augusto Damineli, Ian R. Stevens, Thomas Madura, Keith Gendreau, Zaven Arzoumanian, Felipe Navarete

Abstract: We report high-precision X-ray monitoring observations in the 0.4-10 keV band of the luminous, long-period colliding-wind binary Eta Carinae up to and through its most recent X-ray minimum/periastron passage in February 2020. Eta Carinae reached its observed maximum X-ray flux on 7 January 2020, at a flux level of $3.30 \times 10^{-10}$ ergs s$^{-1}$ cm$^{-2}$, followed by a rapid plunge to its ob… ▽ More We report high-precision X-ray monitoring observations in the 0.4-10 keV band of the luminous, long-period colliding-wind binary Eta Carinae up to and through its most recent X-ray minimum/periastron passage in February 2020. Eta Carinae reached its observed maximum X-ray flux on 7 January 2020, at a flux level of $3.30 \times 10^{-10}$ ergs s$^{-1}$ cm$^{-2}$, followed by a rapid plunge to its observed minimum flux, $0.03 \times 10^{-10}$ ergs s$^{-1}$ cm$^{-2}$ near 17 February 2020. The NICER observations show an X-ray recovery from minimum of only $\sim$16 days, the shortest X-ray minimum observed so far. We provide new constraints of the "deep" and "shallow" minimum intervals. Variations in the characteristic X-ray temperature of the hottest observed X-ray emission indicate that the apex of the wind-wind "bow shock" enters the companion's wind acceleration zone about 81 days before the start of the X-ray minimum. There is a step-like increase in column density just before the X-ray minimum, probably associated with the presence of dense clumps near the shock apex. During recovery and after, the column density shows a smooth decline, which agrees with previous $N_{H}$ measurements made by SWIFT at the same orbital phase, indicating that changes in mass-loss rate are only a few percent over the two cycles. Finally, we use the variations in the X-ray flux of the outer ejecta seen by NICER to derive a kinetic X-ray luminosity of the ejecta of $\sim 10^{41}$ ergs s$^{-1}$ near the time of the "Great Eruption'. △ Less

Submitted 7 July, 2022; originally announced July 2022.

Consistency of the Parkes Pulsar Timing Array Signal with a Nanohertz Gravitational-wave Background

Authors: Boris Goncharov, Eric Thrane, Ryan M. Shannon, Jan Harms, N. D. Ramesh Bhat, George Hobbs, Matthew Kerr, Richard N. Manchester, Daniel J. Reardon, Christopher J. Russell, Xing-Jiang Zhu, Andrew Zic

Abstract: Pulsar timing array experiments have recently reported strong evidence for a common-spectrum stochastic process with a strain spectral index consistent with that expected of a nanohertz-frequency gravitational-wave background, but with negligible yet non-zero evidence for spatial correlations required for a definitive detection. However, it was pointed out by the Parkes Pulsar Timing Array (PPTA)… ▽ More Pulsar timing array experiments have recently reported strong evidence for a common-spectrum stochastic process with a strain spectral index consistent with that expected of a nanohertz-frequency gravitational-wave background, but with negligible yet non-zero evidence for spatial correlations required for a definitive detection. However, it was pointed out by the Parkes Pulsar Timing Array (PPTA) collaboration that the same models used in recent analyses resulted in strong evidence for a common-spectrum process in simulations where none is present. In this work, we introduce a methodology to distinguish pulsar power spectra with the same amplitude from noise power spectra of similar but distinct amplitudes. The former is the signature of a spatially uncorrelated pulsar term of a nanohertz gravitational-wave background, whereas the latter could represent ensemble pulsar noise properties. We test the methodology on simulated data sets. We find that the reported common process in PPTA pulsars is indeed consistent with the spectral feature of a pulsar term. We recommend this methodology as one of the validity tests that the real astrophysical and cosmological backgrounds should pass, as well as for inferences about the spatially uncorrelated component of the background. △ Less

Submitted 21 June, 2022; v1 submitted 8 June, 2022; originally announced June 2022.

Comments: 9 pages, 2 figures

Journal ref: The Astrophysical Journal Letters, Volume 932, Number 2 (2022)

Eta Carinae: an evolving view of the central binary, its interacting winds and its foreground ejecta

Authors: Theodore R. Gull, D. John Hillier, Henrik Hartman, Michael F. Corcoran, Augusto Damineli, David Espinoza-Galeas, Kenji Hamaguchi, Felipe Navarete, Krister Nielsen, Thomas Madura, Anthony F. J. Moffat, Patrick Morris, Noel D. Richardson, Christopher M. P. Russell, Ian R. Stevens, Gerd Weigelt

Abstract: FUV spectra of Eta Car, recorded across two decades with HST/STIS, document multiple changes in resonant lines caused by dissipating extinction in our line of sight. The FUV flux has increased nearly ten-fold which has led to increased ionization of the multiple shells within the Homunculus and photo-destruction of molecular hydrogen. Comparison of observed resonant line profiles with CMFGEN model… ▽ More FUV spectra of Eta Car, recorded across two decades with HST/STIS, document multiple changes in resonant lines caused by dissipating extinction in our line of sight. The FUV flux has increased nearly ten-fold which has led to increased ionization of the multiple shells within the Homunculus and photo-destruction of molecular hydrogen. Comparison of observed resonant line profiles with CMFGEN model profiles allows separation of wind-wind collision and shell absorptions from the primary wind, P Cygni profiles.The dissipating occulter preferentially obscured the central binary and interacting winds relative to the very extended primary wind. We are now able to monitor changes in the colliding winds with orbital phase. High velocity transient absorptions occurred across the most recent periastron passage, indicating acceleration of the primary wind by the secondary wind which leads to a downstream, high velocity bowshock that is newly generated every orbital period. There is no evidence of changes in the properties of the binary winds. △ Less

Submitted 30 May, 2022; originally announced May 2022.

Comments: 36 pages, 22 figures, accepted Astrophysical Journal

Magnetic Flux Transport Identification of Active Reconnection: MMS Observations in Earth's Magnetosphere

Authors: Yi Qi, Tak Chu Li, Christopher T. Russell, Robert E. Ergun, Ying-Dong Jia, Mark Hubbert

Abstract: Magnetic reconnection plays an important role in converting energy while modifying field topology. This process takes place in varied plasma environments in which the transport of magnetic flux is intrinsic. Identifying active magnetic reconnection sites in in-situ observations is challenging. A new technique, Magnetic Flux Transport (MFT) analysis, has been developed recently and proven in numeri… ▽ More Magnetic reconnection plays an important role in converting energy while modifying field topology. This process takes place in varied plasma environments in which the transport of magnetic flux is intrinsic. Identifying active magnetic reconnection sites in in-situ observations is challenging. A new technique, Magnetic Flux Transport (MFT) analysis, has been developed recently and proven in numerical simulation for identifying active reconnection efficiently and accurately. In this study, we examine the MFT process in 37 previously reported electron diffusion region (EDR)/reconnection-line crossing events at the dayside magnetopause and in the magnetotail and turbulent magnetosheath using Magnetospheric Multiscale measurements. The coexisting inward and outward MFT flows at an X-point provides a signature that magnetic field lines become disconnected and reconnected. The application of MFT analysis to in-situ observations demonstrates that MFT can successfully identify active reconnection sites under complex varied conditions, including asymmetric and turbulent upstream conditions. It also provides a higher rate of identification than plasma outflow jets alone. MFT can be applied to in situ measurements from both single- and multi-spacecraft missions and laboratory experiments. △ Less

Submitted 12 March, 2022; originally announced March 2022.

Comments: 6 pages, 1 table, 3 figures

Journal ref: The Astrophysical Journal Letters, 926, L34, 2022

The International Pulsar Timing Array second data release: Search for an isotropic Gravitational Wave Background

Authors: J. Antoniadis, Z. Arzoumanian, S. Babak, M. Bailes, A. -S. Bak Nielsen, P. T. Baker, C. G. Bassa, B. Becsy, A. Berthereau, M. Bonetti, A. Brazier, P. R. Brook, M. Burgay, S. Burke-Spolaor, R. N. Caballero, J. A. Casey-Clyde, A. Chalumeau, D. J. Champion, M. Charisi, S. Chatterjee, S. Chen, I. Cognard, J. M. Cordes, N. J. Cornish, F. Crawford , et al. (101 additional authors not shown)

Abstract: We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^α$, we found strong evidence for a spectrally… ▽ More We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^α$, we found strong evidence for a spectrally-similar low-frequency stochastic process of amplitude $A = 3.8^{+6.3}_{-2.5}\times10^{-15}$ and spectral index $α= -0.5 \pm 0.5$, where the uncertainties represent 95\% credible regions, using information from the auto- and cross-correlation terms between the pulsars in the array. For a spectral index of $α= -2/3$, as expected from a population of inspiralling supermassive black hole binaries, the recovered amplitude is $A = 2.8^{+1.2}_{-0.8}\times10^{-15}$. Nonetheless, no significant evidence of the Hellings-Downs correlations that would indicate a gravitational-wave origin was found. We also analyzed the constituent data from the individual pulsar timing arrays in a consistent way, and clearly demonstrate that the combined international data set is more sensitive. Furthermore, we demonstrate that this combined data set produces comparable constraints to recent single-array data sets which have more data than the constituent parts of the combination. Future international data releases will deliver increased sensitivity to gravitational wave radiation, and significantly increase the detection probability. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 17 pages, 12 figures, accepted in MNRAS

High-precision search for dark photon dark matter with the Parkes Pulsar Timing Array

Authors: Xiao Xue, Zi-Qing Xia, Xingjiang Zhu, Yue Zhao, Jing Shu, Qiang Yuan, N. D. Ramesh Bhat, Andrew D. Cameron, Shi Dai, Yi Feng, Boris Goncharov, George Hobbs, Eric Howard, Richard N. Manchester, Aditya Parthasarathy, Daniel J. Reardon, Christopher J. Russell, Ryan M. Shannon, Renée Spiewak, Nithyanandan Thyagarajan, Jingbo Wang, Lei Zhang, Songbo Zhang

Abstract: The nature of dark matter remains obscure in spite of decades of experimental efforts. The mass of dark matter candidates can span a wide range, and its coupling with the Standard Model sector remains uncertain. All these unknowns make the etection of dark matter extremely challenging. Ultralight dark matter, with $m \sim10^{-22}$ eV, is proposed to reconcile the disagreements between observations… ▽ More The nature of dark matter remains obscure in spite of decades of experimental efforts. The mass of dark matter candidates can span a wide range, and its coupling with the Standard Model sector remains uncertain. All these unknowns make the etection of dark matter extremely challenging. Ultralight dark matter, with $m \sim10^{-22}$ eV, is proposed to reconcile the disagreements between observations and predictions from simulations of small-scale structures in the cold dark matter paradigm, while remaining consistent with other observations. Because of its large de Broglie wavelength and large local occupation number within galaxies, ultralight dark matter behaves like a coherently oscillating background field with an oscillating frequency dependent on its mass. If the dark matter particle is a spin-1 dark photon, such as the $U(1)_B$ or $U(1)_{B-L}$ gauge boson, it can induce an external oscillating force and lead to displacements of test masses. Such an effect would be observable in the form of periodic variations in the arrival times of radio pulses from highly stable millisecond pulsars. In this study, we search for evidence of ultralight dark photon dark matter (DPDM) using 14-year high-precision observations of 26 pulsars collected with the Parkes Pulsar Timing Array. While no statistically significant signal is found, we place constraints on coupling constants for the $U(1)_B$ and $U(1)_{B-L}$ DPDM. Compared with other experiments, the limits on the dimensionless coupling constant $ε$ achieved in our study are improved by up to two orders of magnitude when the dark photon mass is smaller than $3\times10^{-22}$~eV ($10^{-22}$~eV) for the $U(1)_{B}$ ($U(1)_{B-L}$) scenario. △ Less

Submitted 26 December, 2021; v1 submitted 14 December, 2021; originally announced December 2021.

Comments: 12 pages, 4 figures

Journal ref: Phys. Rev. Research 4, L012022 (2022)

Constraining cosmological phase transitions with the Parkes Pulsar Timing Array

Authors: Xiao Xue, Ligong Bian, Jing Shu, Qiang Yuan, Xingjiang Zhu, N. D. Ramesh Bhat, Shi Dai, Yi Feng, Boris Goncharov, George Hobbs, Eric Howard, Richard N. Manchester, Christopher J. Russell, Daniel J. Reardon, Ryan M. Shannon, Renée Spiewak, Nithyanandan Thyagarajan, Jingbo Wang

Abstract: A cosmological first-order phase transition is expected to produce a stochastic gravitational wave background. If the phase transition temperature is on the MeV scale, the power spectrum of the induced stochastic gravitational waves peaks around nanohertz frequencies, and can thus be probed with high-precision pulsar timing observations. We search for such a stochastic gravitational wave backgroun… ▽ More A cosmological first-order phase transition is expected to produce a stochastic gravitational wave background. If the phase transition temperature is on the MeV scale, the power spectrum of the induced stochastic gravitational waves peaks around nanohertz frequencies, and can thus be probed with high-precision pulsar timing observations. We search for such a stochastic gravitational wave background with the latest data set of the Parkes Pulsar Timing Array. We find no evidence for a Hellings-Downs spatial correlation as expected for a stochastic gravitational wave background. Therefore, we present constraints on first-order phase transition model parameters. Our analysis shows that pulsar timing is particularly sensitive to the low-temperature ($T \sim 1 - 100$ MeV) phase transition with a duration $(β/H_*)^{-1}\sim 10^{-2}-10^{-1}$ and therefore can be used to constrain the dark and QCD phase transitions. △ Less

Submitted 25 November, 2022; v1 submitted 6 October, 2021; originally announced October 2021.

Comments: 9 pages, 4 figures, 2 tables

Journal ref: Phys. Rev. Lett. 127, 251303, Published 15 December 2021

Competitive X-ray and Optical Cooling in the Collisionless Shocks of WR 140

Authors: A. M. T. Pollock, M. F. Corcoran, I. R. Stevens, C. M. P. Russell, K. Hamaguchi, P. M. Williams, A. F. J. Moffat, G. Weigelt, V. Shenavrin, N. D. Richardson, D. Espinoza, S. A. Drake

Abstract: WR 140 is a long-period, highly eccentric Wolf-Rayet star binary system with exceptionally well-determined orbital and stellar parameters. Bright, variable X-ray emission is generated in shocks produced by the collision of the winds of the WC7pd+O5.5fc component stars. We discuss the variations in the context of the colliding-wind model using broad-band spectrometry from the RXTE, SWIFT, and NICER… ▽ More WR 140 is a long-period, highly eccentric Wolf-Rayet star binary system with exceptionally well-determined orbital and stellar parameters. Bright, variable X-ray emission is generated in shocks produced by the collision of the winds of the WC7pd+O5.5fc component stars. We discuss the variations in the context of the colliding-wind model using broad-band spectrometry from the RXTE, SWIFT, and NICER observatories obtained over 20 years and nearly 1000 observations through 3 consecutive 7.94-year orbits including 3 periastron passages. The X-ray luminosity varies as expected with the inverse of the stellar separation over most of the orbit: departures near periastron are produced when cooling shifts to excess optical emission in CIII $\lambda5696$ in particular. We use X-ray absorption to estimate mass-loss rates for both stars and to constrain the system morphology. The absorption maximum coincides closely with inferior conjunction of the WC star and provides evidence of the ion-reflection mechanism that underlies the formation of collisionless shocks governed by magnetic fields probably generated by the Weibel instability. Comparisons with K-band emission and HeI $λ$10830 absorption show that both are correlated after periastron with the asymmetric X-ray absorption. Dust appears within a few days of periastron suggesting formation within shocked gas near the stagnation point. X-ray flares seen in $η$ Carinae have not occurred in WR 140, suggesting the absence of large-scale wind inhomogeneities. Relatively constant soft emission revealed during the X-ray minimum is probably not from recombining plasma entrained in outflowing shocked gas. △ Less

Submitted 20 September, 2021; originally announced September 2021.

Comments: 52 pages, 17+1 figures, Accepted for publication in The Astrophysical Journal 23 August 2021

81 New Candidate Fast Radio Bursts in Parkes Archive

Authors: X. Yang, S. -B. Zhang, J. -S. Wang, G. Hobbs, T. -R. Sun, R. N. Manchester, J. -J. Geng, C. J. Russell, R. Luo, Z. -F. Tang, C. Wang, J. -J. Wei, L. Staveley-Smith, S. Dai, Y. Li, Y. -Y. Yang, X. -F. Wu

Abstract: We have searched for weak fast radio burst (FRB) events using a database containing 568,736,756 transient events detected using the Parkes radio telescope between 1997 and 2001. In order to classify these pulses, and to identify likely FRB candidates, we used a machine learning algorithm based on ResNet. We identified 81 new candidate FRBs and provide details of their positions, event times, and d… ▽ More We have searched for weak fast radio burst (FRB) events using a database containing 568,736,756 transient events detected using the Parkes radio telescope between 1997 and 2001. In order to classify these pulses, and to identify likely FRB candidates, we used a machine learning algorithm based on ResNet. We identified 81 new candidate FRBs and provide details of their positions, event times, and dispersion measures. These events were detected in only one beam of the Parkes multibeam receiver. We used a relatively low S/N cutoff threshold when selecting these bursts and some have dispersion measures only slightly exceeding the expected Galactic contribution. We therefore present these candidate FRBs as a guide for follow-up observations in the search for repeating FRBs. △ Less

Submitted 31 August, 2021; v1 submitted 1 August, 2021; originally announced August 2021.

Comments: 8 pages, 5 figures, 2 tables, accepted by MNRAS

On the evidence for a common-spectrum process in the search for the nanohertz gravitational-wave background with the Parkes Pulsar Timing Array

Authors: Boris Goncharov, R. M. Shannon, D. J. Reardon, G. Hobbs, A. Zic, M. Bailes, M. Curylo, S. Dai, M. Kerr, M. E. Lower, R. N. Manchester, R. Mandow, H. Middleton, M. T. Miles, A. Parthasarathy, E. Thrane, N. Thyagarajan, X. Xue, X. J. Zhu, A. D. Cameron, Y. Feng, R. Luo, C. J. Russell, J. Sarkissian, R. Spiewak , et al. (4 additional authors not shown)

Abstract: A nanohertz-frequency stochastic gravitational-wave background can potentially be detected through the precise timing of an array of millisecond pulsars. This background produces low-frequency noise in the pulse arrival times that would have a characteristic spectrum common to all pulsars and a well-defined spatial correlation. Recently the North American Nanohertz Observatory for Gravitational Wa… ▽ More A nanohertz-frequency stochastic gravitational-wave background can potentially be detected through the precise timing of an array of millisecond pulsars. This background produces low-frequency noise in the pulse arrival times that would have a characteristic spectrum common to all pulsars and a well-defined spatial correlation. Recently the North American Nanohertz Observatory for Gravitational Waves collaboration (NANOGrav) found evidence for the common-spectrum component in their 12.5-year data set. Here we report on a search for the background using the second data release of the Parkes Pulsar Timing Array. If we are forced to choose between the two NANOGrav models $\unicode{x2014}$ one with a common-spectrum process and one without $\unicode{x2014}$ we find strong support for the common-spectrum process. However, in this paper, we consider the possibility that the analysis suffers from model misspecification. In particular, we present simulated data sets that contain noise with distinctive spectra but show strong evidence for a common-spectrum process under the standard assumptions. The Parkes data show no significant evidence for, or against, the spatially correlated Hellings-Downs signature of the gravitational-wave background. Assuming we did observe the process underlying the spatially uncorrelated component of the background, we infer its amplitude to be $A = 2.2^{+0.4}_{-0.3} \times 10^{-15}$ in units of gravitational-wave strain at a frequency of $1\, \text{yr}^{-1}$. Extensions and combinations of existing and new data sets will improve the prospects of identifying spatial correlations that are necessary to claim a detection of the gravitational-wave background. △ Less

Submitted 11 August, 2021; v1 submitted 26 July, 2021; originally announced July 2021.

Comments: 10 pages, 5 figures

Journal ref: The Astrophysical Journal Letters, Volume 917, Number 2 (2021)

The Parkes pulsar timing array second data release: Timing analysis

Authors: D. J. Reardon, R. M. Shannon, A. D. Cameron, B. Goncharov, G. B. Hobbs, H. Middleton, M. Shamohammadi, N. Thyagarajan, M. Bailes, N. D. R. Bhat, S. Dai, M. Kerr, R. N. Manchester, C. J. Russell, R. Spiewak, J. B. Wang, X. J. Zhu

Abstract: The main goal of pulsar timing array experiments is to detect correlated signals such as nanohertz-frequency gravitational waves. Pulsar timing data collected in dense monitoring campaigns can also be used to study the stars themselves, their binary companions, and the intervening ionised interstellar medium. Timing observations are extraordinarily sensitive to changes in path length between the p… ▽ More The main goal of pulsar timing array experiments is to detect correlated signals such as nanohertz-frequency gravitational waves. Pulsar timing data collected in dense monitoring campaigns can also be used to study the stars themselves, their binary companions, and the intervening ionised interstellar medium. Timing observations are extraordinarily sensitive to changes in path length between the pulsar and the Earth, enabling precise measurements of the pulsar positions, distances and velocities, and the shapes of their orbits. Here we present a timing analysis of 25 pulsars observed as part of the Parkes Pulsar Timing Array (PPTA) project over time spans of up to 24 yr. The data are from the second data release of the PPTA, which we have extended by including legacy data. We make the first detection of Shapiro delay in four Southern pulsars (PSRs J1017$-$7156, J1125$-$6014, J1545$-$4550, and J1732$-$5049), and of parallax in six pulsars. The prominent Shapiro delay of PSR J1125$-$6014 implies a neutron star mass of $M_p = 1.5 \pm 0.2 M_\odot$ (68% credibility interval). Measurements of both Shapiro delay and relativistic periastron advance in PSR J1600$-$3053 yield a large but uncertain pulsar mass of $M_p = 2.06^{+0.44}_{-0.41}$ M$_\odot$ (68% credibility interval). We measure the distance to PSR J1909$-$3744 to a precision of 10 lyr, indicating that for gravitational wave periods over a decade, the pulsar provides a coherent baseline for pulsar timing array experiments. △ Less

Submitted 16 July, 2021; v1 submitted 9 July, 2021; originally announced July 2021.

Comments: 44 pages, 32 figures, 7 tables. Accepted for publication in MNRAS; v2 includes corrected typos in parameter tables and reference added

The brittle boulders of dwarf planet Ceres

Authors: Stefan Schröder, Uri Carsenty, Ernst Hauber, Carol Raymond, Christopher Russell

Abstract: We mapped all boulders larger than 105 m on the surface of dwarf planet Ceres using images of the Dawn framing camera acquired in the Low Altitude Mapping Orbit (LAMO). We find that boulders on Ceres are more numerous towards high latitudes and have a maximum lifetime of $150 \pm 50$ Ma, based on crater counts. These characteristics are distinctly different from those of boulders on asteroid (4) V… ▽ More We mapped all boulders larger than 105 m on the surface of dwarf planet Ceres using images of the Dawn framing camera acquired in the Low Altitude Mapping Orbit (LAMO). We find that boulders on Ceres are more numerous towards high latitudes and have a maximum lifetime of $150 \pm 50$ Ma, based on crater counts. These characteristics are distinctly different from those of boulders on asteroid (4) Vesta, an earlier target of Dawn, which implies that Ceres boulders are mechanically weaker. Clues to their properties can be found in the composition of Ceres' complex crust, which is rich in phyllosilicates and salts. As water ice is though to be present only meters below the surface, we suggest that boulders also harbor ice. Furthermore, the boulder size-frequency distribution is best fit by a Weibull distribution rather than the customary power law, just like for Vesta boulders. This finding is robust in light of possible types of size measurement error. △ Less

Submitted 25 May, 2021; originally announced May 2021.

Comments: Accepted by the Planetary Science Journal, 31 pages, 14 figures

Comparative study of electric currents and energetic particle fluxes in a solar flare and Earth magnetospheric substorm

Authors: Anton Artemyev, Ivan Zimovets, Ivan Sharykin, Yukitoshi Nishimura, Cooper Downs, James Weygand, Robyn Fiori, Xiao-Jia Zhang, Andrei Runov, Marco Velli, Vassilis Angelopoulos, Olga Panasenco, Christopher Russell, Yoshizumi Miyoshi, Satoshi Kasahara, Ayako Matsuoka, Shoichiro Yokota, Kunihiro Keika, Tomoaki Hori, Yoichi Kazama, Shiang-Yu Wang, Iku Shinohara, Yasunobu Ogawa

Abstract: Magnetic field-line reconnection is a universal plasma process responsible for the conversion of magnetic field energy to the plasma heating and charged particle acceleration. Solar flares and Earth's magnetospheric substorms are two most investigated dynamical systems where magnetic reconnection is believed to be responsible for global magnetic field reconfiguration and energization of plasma pop… ▽ More Magnetic field-line reconnection is a universal plasma process responsible for the conversion of magnetic field energy to the plasma heating and charged particle acceleration. Solar flares and Earth's magnetospheric substorms are two most investigated dynamical systems where magnetic reconnection is believed to be responsible for global magnetic field reconfiguration and energization of plasma populations. Such a reconfiguration includes formation of a long-living current systems connecting the primary energy release region and cold dense conductive plasma of photosphere/ionosphere. In both flares and substorms the evolution of this current system correlates with formation and dynamics of energetic particle fluxes. Our study is focused on this similarity between flares and substorms. Using a wide range of datasets available for flare and substorm investigations, we compare qualitatively dynamics of currents and energetic particle fluxes for one flare and one substorm. We showed that there is a clear correlation between energetic particle bursts (associated with energy release due to magnetic reconnection) and magnetic field reconfiguration/formation of current system. We then discuss how datasets of in-situ measurements in the magnetospheric substorm can help in interpretation of datasets gathered for the solar flare. △ Less

Submitted 8 May, 2021; originally announced May 2021.

Conditions in the WR 140 wind-collision region revealed by the 1.083-micron He I line profile

Authors: Peredur M. Williams, Watson P. Varricatt, André-Nicolas Chené, Michael F. Corcoran, Ted R. Gull, Kenji Hamaguchi, Anthony F. J. Moffat, Andrew M. T. Pollock, Noel D. Richardson, Christopher M. P. Russell, Andreas A. C. Sander, Ian R. Stevens, Gerd Weigelt

Abstract: We present spectroscopy of the P~Cygni profile of the 1.083-micron He I line in the WC7 + O5 colliding-wind binary (CWB) WR 140 (HD 193793), observed in 2008, before its periastron passage in 2009, and in 2016-17, spanning the subsequent periastron passage. Both absorption and emission components showed strong variations. The variation of the absorption component as the O5 star was occulted by the… ▽ More We present spectroscopy of the P~Cygni profile of the 1.083-micron He I line in the WC7 + O5 colliding-wind binary (CWB) WR 140 (HD 193793), observed in 2008, before its periastron passage in 2009, and in 2016-17, spanning the subsequent periastron passage. Both absorption and emission components showed strong variations. The variation of the absorption component as the O5 star was occulted by the wind-collision region (WCR) sets a tight constraint on its geometry. While the sightline to the O5 star traversed the WCR, the strength and breadth of the absorption component varied significantly on time-scales of days. An emission sub-peak was observed on all our profiles. The variation of its radial velocity with orbital phase was shown to be consistent with formation in the WCR as it swung round the stars in their orbit. Modelling the profile gave a measure of the extent of the sub-peak forming region. In the phase range 0.93-0.99, the flux in the sub-peak increased steadily, approximately inversely proportionally to the stellar separation, indicating that the shocked gas in the WCR where the line was formed was adiabatic. After periastron, the sub-peak flux was anomalously strong and varied rapidly, suggesting formation in clumps down-stream in the WCR. For most of the time, its flux exceeded the 2-10-keV X-ray emission, showing it to be a significant coolant of the shocked wind. △ Less

Submitted 18 February, 2021; originally announced February 2021.

Comments: 19 pages, 14 figures, accepted for publication in the MNRAS

In situ evidence of ion acceleration between consecutive reconnection jet fronts

Authors: Filomena Catapano, Alessandro Retino, Gaetano Zimbardo, Alexandra Alexandrova, Ian J. Cohen, Drew L. Turner, Olivier Le Contel, Giulia Cozzani, Silvia Perri, Antonella Greco, Hugo Breuillard, Dominique Delcourt, Laurent Mirioni, Yuri Khotyaintsev, Andris Vaivads, Barbara L. Giles, Barry H. Mauk, Stephen A. Fuselier, Roy B. Torbert, Christopher T. Russell, Per A. Lindqvist, Robert E. Ergun, Thomas Moore, James L. Burch

Abstract: Processes driven by unsteady reconnection can efficiently accelerate particles in many astrophysical plasmas. An example are the reconnection jet fronts in an outflow region. We present evidence of suprathermal ion acceleration between two consecutive reconnection jet fronts observed by the Magnetospheric Multiscale mission in the terrestrial magnetotail. An earthward propagating jet is approached… ▽ More Processes driven by unsteady reconnection can efficiently accelerate particles in many astrophysical plasmas. An example are the reconnection jet fronts in an outflow region. We present evidence of suprathermal ion acceleration between two consecutive reconnection jet fronts observed by the Magnetospheric Multiscale mission in the terrestrial magnetotail. An earthward propagating jet is approached by a second faster jet. Between the jets, the thermal ions are mostly perpendicular to magnetic field, are trapped and are gradually accelerated in the parallel direction up to 150 keV. Observations suggest that ions are predominantly accelerated by a Fermi-like mechanism in the contracting magnetic bottle formed between the two jet fronts. The ion acceleration mechanism is presumably efficient in other environments where jet fronts produced by variable rates of reconnection are common and where the interaction of multiple jet fronts can also develop a turbulent environment, e.g. in stellar and solar eruptions. △ Less

Submitted 30 November, 2020; originally announced December 2020.

Searching for gravitational wave bursts from cosmic string cusps with the Parkes Pulsar Timing Array

Authors: N. Yonemaru, S. Kuroyanagi, G. Hobbs, K. Takahashi, X. -J. Zhu, W. A. Coles, S. Dai, E. Howard, R. Manchester, D. Reardon, C. Russell, R. Shannon, N. Thyagarajan, R. Spiewak, J. -B. Wang

Abstract: Cosmic strings are potential gravitational wave (GW) sources that can be probed by pulsar timing arrays (PTAs). In this work we develop a detection algorithm for a GW burst from a cusp on a cosmic string, and apply it to Parkes PTA data. We find four events with a false alarm probability less than 1%. However further investigation shows that all of these are likely to be spurious. As there are no… ▽ More Cosmic strings are potential gravitational wave (GW) sources that can be probed by pulsar timing arrays (PTAs). In this work we develop a detection algorithm for a GW burst from a cusp on a cosmic string, and apply it to Parkes PTA data. We find four events with a false alarm probability less than 1%. However further investigation shows that all of these are likely to be spurious. As there are no convincing detections we place upper limits on the GW amplitude for different event durations. From these bounds we place limits on the cosmic string tension of G mu ~ 10^{-5}, and highlight that this bound is independent from those obtained using other techniques. We discuss the physical implications of our results and the prospect of probing cosmic strings in the era of Square Kilometre Array (SKA). △ Less

Submitted 26 November, 2020; originally announced November 2020.

Comments: Accepted for publication by MNRAS

Identifying and mitigating noise sources in precision pulsar timing data sets

Authors: Boris Goncharov, D. J. Reardon, R. M. Shannon, Xing-Jiang Zhu, Eric Thrane, M. Bailes, N. D. R. Bhat, S. Dai, G. Hobbs, M. Kerr, R. N. Manchester, S. Osłowski, A. Parthasarathy, C. J. Russell, R. Spiewak, N. Thyagarajan, J. B. Wang

Abstract: Pulsar timing array projects measure the pulse arrival times of millisecond pulsars for the primary purpose of detecting nanohertz-frequency gravitational waves. The measurements include contributions from a number of astrophysical and instrumental processes, which can either be deterministic or stochastic. It is necessary to develop robust statistical and physical models for these noise processes… ▽ More Pulsar timing array projects measure the pulse arrival times of millisecond pulsars for the primary purpose of detecting nanohertz-frequency gravitational waves. The measurements include contributions from a number of astrophysical and instrumental processes, which can either be deterministic or stochastic. It is necessary to develop robust statistical and physical models for these noise processes because incorrect models diminish sensitivity and may cause a spurious gravitational wave detection. Here we characterise noise processes for the 26 pulsars in the second data release of the Parkes Pulsar Timing Array using Bayesian inference. In addition to well-studied noise sources found previously in pulsar timing array data sets such as achromatic timing noise and dispersion measure variations, we identify new noise sources including time-correlated chromatic noise that we attribute to variations in pulse scattering. We also identify "exponential dip" events in four pulsars, which we attribute to magnetospheric effects as evidenced by pulse profile shape changes observed for three of the pulsars. This includes an event in PSR J1713$+$0747, which had previously been attributed to interstellar propagation. We present noise models to be used in searches for gravitational waves. We outline a robust methodology to evaluate the performance of noise models and identify unknown signals in the data. The detection of variations in pulse profiles highlights the need to develop efficient profile domain timing methods. △ Less

Submitted 1 November, 2020; v1 submitted 12 October, 2020; originally announced October 2020.

Comments: 18 pages, 7 figures, 5 tables

The surface of (1) Ceres in visible light as seen by Dawn/VIR

Authors: B. Rousseau, M. C. De Sanctis, A. Raponi, M. Ciarniello, E. Ammannito, A. Frigeri, M. Ferrari, S. De Angelis, F. Tosi, S. E. Schröder, C. A. Raymond, C. T. Russell

Abstract: We study the surface of Ceres at visible wavelengths, as observed by the Visible and InfraRed mapping spectrometer (VIR) onboard the Dawn spacecraft, and analyze the variations of various spectral parameters across the whole surface. We also focus on several noteworthy areas of the surface of this dwarf planet. We made use of the newly corrected VIR visible data to build global maps of a calibrate… ▽ More We study the surface of Ceres at visible wavelengths, as observed by the Visible and InfraRed mapping spectrometer (VIR) onboard the Dawn spacecraft, and analyze the variations of various spectral parameters across the whole surface. We also focus on several noteworthy areas of the surface of this dwarf planet. We made use of the newly corrected VIR visible data to build global maps of a calibrated radiance factor at 550 nm, with two color composites and three spectral slopes between 400 and 950 nm. We have made these maps available for the community via the Aladin Desktop software. Ceres surface shows diverse spectral behaviors in the visible range. The color composite and the spectral slope between 480 and 800 nm highlight fresh impact craters and young geologic formations of endogenous origin, which appear bluer than the rest of the surface. The steep slope before 465 nm displays very distinct variations and may be a proxy for the absorptions caused by the $O_2^{-}$ -> $Fe^{3+}$ or the $2Fe^{3+}$ -> $Fe^{2+} + Fe^{4+}$ charge transfer, if the latter are found to be responsible for the drop in this spectral range. We notice several similarities between the spectral slopes and the abundance of phyllosilicates detected in the infrared by the VIR, whereas no correlation can be clearly established with carbonate species. The region of the Dantu impact crater presents a peculiar spectral behavior (especially through the color and the spectral slope before 465 nm) suggesting a change in composition or in the surface physical properties that is not observed elsewhere on Ceres. △ Less

Submitted 7 October, 2020; originally announced October 2020.

Comments: 22 pages, 14 + 5 figures. Published online on Astronomy and Astrophysics on 07-10-2020

Journal ref: A&A 642, A74 (2020)

Precision orbital dynamics from interstellar scintillation arcs for PSR J0437-4715

Authors: Daniel J. Reardon, William A. Coles, Matthew Bailes, N. D. Ramesh Bhat, Shi Dai, George B. Hobbs, Matthew Kerr, Richard N. Manchester, Stefan Oslowski, Aditya Parthasarathy, Christopher J. Russell, Ryan M. Shannon, Renee Spiewak, Lawrence Toomey, Artem V. Tuntsov, Willem van Straten, Mark A. Walker, Jingbo Wang, Lei Zhang, Xing-Jiang Zhu

Abstract: Intensity scintillations of radio pulsars are known to originate from interference between waves scattered by the electron density irregularities of interstellar plasma, often leading to parabolic arcs in the two-dimensional power spectrum of the recorded dynamic spectrum. The degree of arc curvature depends on the distance to the scattering plasma and its transverse velocity with respect to the l… ▽ More Intensity scintillations of radio pulsars are known to originate from interference between waves scattered by the electron density irregularities of interstellar plasma, often leading to parabolic arcs in the two-dimensional power spectrum of the recorded dynamic spectrum. The degree of arc curvature depends on the distance to the scattering plasma and its transverse velocity with respect to the line-of-sight. We report the observation of annual and orbital variations in the curvature of scintillation arcs over a period of 16 years for the bright millisecond pulsar, PSR J0437-4715. These variations are the signature of the relative transverse motions of the Earth, pulsar, and scattering medium, which we model to obtain precise measurements of parameters of the pulsar's binary orbit and the scattering medium itself. We observe two clear scintillation arcs in most of our $>$5000 observations and we show that they originate from scattering by thin screens located at distances $D_1 = 89.8 \pm 0.4$ pc and $D_2 = 124 \pm 3$ pc from Earth. The best-fit scattering model we derive for the brightest arc yields the pulsar's orbital inclination angle $i = 137.1 \pm 0.3^\circ$, and longitude of ascending node, $Ω=206.3\pm0.4^\circ$. Using scintillation arcs for precise astrometry and orbital dynamics can be superior to modelling variations in the diffractive scintillation timescale, because the arc curvature is independent of variations in the level of turbulence of interstellar plasma. This technique can be used in combination with pulsar timing to determine the full three-dimensional orbital geometries of binary pulsars, and provides parameters essential for testing theories of gravity and constraining neutron star masses. △ Less

Submitted 27 September, 2020; originally announced September 2020.

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

The boulder population of asteroid 4 Vesta: Size-frequency distribution and survival time

Authors: Stefan E. Schröder, Uri Carsenty, Ernst Hauber, Franziska Schulzeck, Carol A. Raymond, Christopher T. Russell

Abstract: Dawn's framing camera observed boulders on the surface of Vesta when the spacecraft was in its lowest orbit (LAMO). We identified, measured, and mapped boulders in LAMO images, which have a scale of 20 m per pixel. We estimate that our sample is virtually complete down to a boulder size of 4 pixels (80 m). The largest boulder is a 400 m-sized block on the Marcia crater floor. Relatively few boulde… ▽ More Dawn's framing camera observed boulders on the surface of Vesta when the spacecraft was in its lowest orbit (LAMO). We identified, measured, and mapped boulders in LAMO images, which have a scale of 20 m per pixel. We estimate that our sample is virtually complete down to a boulder size of 4 pixels (80 m). The largest boulder is a 400 m-sized block on the Marcia crater floor. Relatively few boulders reside in a large area of relatively low albedo, surmised to be the carbon-rich ejecta of the Veneneia basin, either because boulders form less easily here or live shorter. By comparing the density of boulders around craters with a known age, we find that the maximum boulder lifetime is about 300 Ma. The boulder size-frequency distribution (SFD) is generally assumed to follow a power law. We fit power laws to the Vesta SFD by means of the maximum likelihood method, but they do not fit well. Our analysis of power law exponents for boulders on other small Solar System bodies suggests that the derived exponent is primarily a function of boulder size range. The Weibull distribution mimics this behavior and fits the Vesta boulder SFD well. The Weibull distribution is often encountered in rock grinding experiments, and may result from the fractal nature of cracks propagating in the rock interior. We propose that, in general, the SFD of particles (including boulders) on the surface of small bodies follows a Weibull distribution rather than a power law. △ Less

Submitted 2 September, 2020; originally announced September 2020.

Direct Measurement of the Solar-Wind Taylor Microscale using MMS Turbulence Campaign Data

Authors: Riddhi Bandyopadhyay, William H. Matthaeus, Alexandros Chasapis, Christopher T. Russell, Robert J. Strangeway, Roy B. Torbert, Barbara L. Giles, Daniel J. Gershman, Craig J. Pollock, James L. Burch

Abstract: Using the novel Magnetospheric Multiscale (MMS) mission data accumulated during the 2019 MMS Solar Wind Turbulence Campaign, we calculate the Taylor microscale $(λ_{\mathrm{T}})$ of the turbulent magnetic field in the solar wind. The Taylor microscale represents the onset of dissipative processes in classical turbulence theory. An accurate estimation of Taylor scale from spacecraft data is, howeve… ▽ More Using the novel Magnetospheric Multiscale (MMS) mission data accumulated during the 2019 MMS Solar Wind Turbulence Campaign, we calculate the Taylor microscale $(λ_{\mathrm{T}})$ of the turbulent magnetic field in the solar wind. The Taylor microscale represents the onset of dissipative processes in classical turbulence theory. An accurate estimation of Taylor scale from spacecraft data is, however, usually difficult due to low time cadence, the effect of time decorrelation, and other factors. Previous reports were based either entirely on the Taylor frozen-in approximation, which conflates time dependence, or that were obtained using multiple datasets, which introduces sample-to-sample variation of plasma parameters, or where inter-spacecraft distance were larger than the present study. The unique configuration of linear formation with logarithmic spacing of the 4 MMS spacecraft, during the campaign, enables a direct evaluation of the $λ_{\mathrm{T}}$ from a single dataset, independent of the Taylor frozen-in approximation. A value of $λ_{\mathrm{T}} \approx 7000 \, \mathrm{km}$ is obtained, which is about 3 times larger than the previous estimates. △ Less

Submitted 19 June, 2020; originally announced June 2020.

Comments: Accepted for publication in the Astrophysical Journal

Interplay of Turbulence and Proton-Microinstability Growth in Space Plasmas

Authors: Riddhi Bandyopadhyay, Ramiz A. Qudsi, William H. Matthaeus, Tulasi N. Parashar, Bennett A. Maruca, S. Peter Gary, Vadim Roytershteyn, Alexandros Chasapis, Barbara L. Giles, Daniel J. Gershman, Craig J. Pollock, Christopher T. Russell, Robert J. Strangeway, Roy B. Torbert, Thomas E. Moore, James L. Burch

Abstract: Numerous prior studies have shown that as proton beta increases, a narrower range of proton temperature anisotropy values is observed. This effect has often been ascribed to the actions of kinetic microinstabilities because the distribution of observational data aligns with contours of constant instability growth rates in the beta-anisotropy plane. However, the linear Vlasov theory of instabilitie… ▽ More Numerous prior studies have shown that as proton beta increases, a narrower range of proton temperature anisotropy values is observed. This effect has often been ascribed to the actions of kinetic microinstabilities because the distribution of observational data aligns with contours of constant instability growth rates in the beta-anisotropy plane. However, the linear Vlasov theory of instabilities assumes a uniform background in which perturbations grow. The established success of linear-microinstability theories suggests that the conditions in regions of extreme temperature anisotropy may remain uniform for a long enough time so that the instabilities have the chance to grow to sufficient amplitude. Turbulence, on the other hand, is intrinsically non-uniform and non-linear. Thin current sheets and other coherent structures generated in a turbulent plasma, may destroy the uniformity fast enough. It is therefore not a-priori obvious whether the presence of intermittency and coherent structures favors or disfavors instabilities. To address this question, we examined the statistical distribution of growth rates associated with proton temperature-anisotropy driven microinstabilities and local nonlinear time scales in turbulent plasmas. Linear growth rates are, on average, substantially less than the local nonlinear rates. However, at the regions of extreme values of temperature anisotropy, near the "edges" of the populated part of the proton temperature anisotropy-parallel beta plane, the instability growth rates are comparable or faster than the turbulence time scales. These results provide a possible answer to the question as to why the linear theory appears to work in limiting plasma excursions in anisotropy and plasma beta. △ Less

Submitted 21 September, 2022; v1 submitted 18 June, 2020; originally announced June 2020.

Comments: Accepted for publication in Physics of Plasmas

Applying saliency-map analysis in searches for pulsars and fast radio bursts

Authors: C. Zhang, C. Wang, G. Hobbs, C. J. Russell, D. Li, S. -B. Zhang, S. Dai, J. -W. Wu, Z. -C. Pan, W. -W. Zhu, L. Toomey, Z. -Y. Ren

Abstract: To investigate the use of saliency-map analysis to aid in searches for transient signals, such as fast radio bursts and individual pulses from radio pulsars. We aim to demonstrate that saliency maps provide the means to understand predictions from machine learning algorithms and can be implemented in piplines used to search for transient events. We have implemented a new deep learning methodology… ▽ More To investigate the use of saliency-map analysis to aid in searches for transient signals, such as fast radio bursts and individual pulses from radio pulsars. We aim to demonstrate that saliency maps provide the means to understand predictions from machine learning algorithms and can be implemented in piplines used to search for transient events. We have implemented a new deep learning methodology to predict whether or not any segment of the data contains a transient event. The algorithm has been trained using real and simulated data sets. We demonstrate that the algorithm is able to identify such events. The output results are visually analysed via the use of saliency maps. We find that saliency maps can produce an enhanced image of any transient feature without the need for de-dispersion or removal of radio frequency interference. Such maps can be used to understand which features in the image were used in making the machine learning decision and to visualise the transient event. Even though the algorithm reported here was developed to demonstrate saliency-map analysis, we have detected, in archival data, a single burst event with dispersion measure of $41$\,cm$^{-3}$pc that is not associated with any currently known pulsar. △ Less

Submitted 22 May, 2020; originally announced May 2020.

Comments: 16 pages, 6 figures, accepted for publication in A&A

Journal ref: A&A 642, A26 (2020)

Statistics of Kinetic Dissipation in Earth's Magnetosheath -- MMS Observations

Authors: Riddhi Bandyopadhyay, William H. Matthaeus, Tulasi N. Parashar, Yan Yang, Alexandros Chasapis, Barbara L. Giles, Daniel J. Gershman, Craig J. Pollock, Christopher T. Russell, Robert J. Strangeway, Roy B. Torbert, Thomas E. Moore, James L. Burch

Abstract: A familiar problem in space and astrophysical plasmas is to understand how dissipation and heating occurs. These effects are often attributed to the cascade of broadband turbulence which transports energy from large scale reservoirs to small scale kinetic degrees of freedom. When collisions are infrequent, local thermodynamic equilibrium is not established. In this case the final stage of energy c… ▽ More A familiar problem in space and astrophysical plasmas is to understand how dissipation and heating occurs. These effects are often attributed to the cascade of broadband turbulence which transports energy from large scale reservoirs to small scale kinetic degrees of freedom. When collisions are infrequent, local thermodynamic equilibrium is not established. In this case the final stage of energy conversion becomes more complex than in the fluid case, and both pressure-dilatation and pressure strain interactions (Pi-D $\equiv -Π_{ij} D_{ij}$) become relevant and potentially important. Pi-D in plasma turbulence has been studied so far primarily using simulations. The present study provides a statistical analysis of Pi-D in the Earth's magnetosheath using the unique measurement capabilities of the Magnetospheric Multiscale (MMS) mission. We find that the statistics of Pi-D in this naturally occurring plasma environment exhibit strong resemblance to previously established fully kinetic simulations results. The conversion of energy is concentrated in space and occurs near intense current sheets, but not within them. This supports recent suggestions that the chain of energy transfer channels involves regional, rather than pointwise, correlations. △ Less

Submitted 19 May, 2020; originally announced May 2020.

Comments: Accepted for publication in Physical Review Letters

The VVV Infrared Variability Catalog (VIVA-I)

Authors: C. E. Ferreira Lopes, N. J. G. Cross, M. Catelan, D. Minniti, M. Hempel, P. W. Lucas, R. Angeloni, F. Jablonsky, V. F. Braga, I. C. Leao, F. R. Herpich, J. Alonso-Garcia, A. Papageorgiou, K. Pichara, R. K. Saito, A. Bradley, J. C. Beamin, C. Cortes, J. R. De Medeiros, Christopher. M. P. Russell

Abstract: Thanks to the VISTA Variables in the Via Lactea (VVV) ESO Public Survey it is now possible to explore a large number of objects in those regions. This paper addresses the variability analysis of all VVV point sources having more than 10 observations in VVVDR4 using a novel approach. In total, the near-IR light curves of 288,378,769 sources were analysed using methods developed in the New Insight I… ▽ More Thanks to the VISTA Variables in the Via Lactea (VVV) ESO Public Survey it is now possible to explore a large number of objects in those regions. This paper addresses the variability analysis of all VVV point sources having more than 10 observations in VVVDR4 using a novel approach. In total, the near-IR light curves of 288,378,769 sources were analysed using methods developed in the New Insight Into Time Series Analysis project. As a result, we present a complete sample having 44, 998, 752 variable star candidates (VVV-CVSC), which include accurate individual coordinates, near-IR magnitudes (ZYJHKs), extinctions A(Ks), variability indices, periods, amplitudes, among other parameters to assess the science. Unfortunately, a side effect of having a highly complete sample, is also having a high level of contamination by non-variable (contamination ratio of non-variables to variables is slightly over 10:1). To deal with this, we also provide some flags and parameters that can be used by the community to de-crease the number of variable candidates without heavily decreasing the completeness of the sample. In particular, we cross-identified 339,601 of our sources with Simbad and AAVSO databases, which provide us with information for these objects at other wavelegths. This sub-sample constitutes a unique resource to study the corresponding near-IR variability of known sources as well as to assess the IR variability related with X-ray and Gamma-Ray sources. On the other hand, the other 99.5% sources in our sample constitutes a number of potentially new objects with variability information for the heavily crowded and reddened regions of the Galactic Plane and Bulge. The present results also provide an important queryable resource to perform variability analysis and to characterize ongoing and future surveys like TESS and LSST. △ Less

Submitted 11 May, 2020; originally announced May 2020.

Comments: 27 pages, 14 figures

Parkes transient events: I. Database of single pulses, initial results and missing FRBs

Authors: S. -B. Zhang, G. Hobbs, C. J. Russell, L. Toomey, S. Dai, J. Dempsey, R. N. Manchester, S. Johnston, L. Staveley-Smith, X. -F. Wu, D. Li, Y. -Y. Yang, S. -Q. Wang, H. Qiu, R. Luo, C. Wang, C. Zhang, L. Zhang, R. Mandow

Abstract: A large number of observations from the Parkes 64\,m-diameter radio telescope, recorded with high time resolution, are publicly available. We have re-processed all of the observations obtained during the first four years (from 1997 to 2001) of the Parkes Multibeam receiver system in order to identify transient events and have built a database that records the 568,736,756 pulse candidates generated… ▽ More A large number of observations from the Parkes 64\,m-diameter radio telescope, recorded with high time resolution, are publicly available. We have re-processed all of the observations obtained during the first four years (from 1997 to 2001) of the Parkes Multibeam receiver system in order to identify transient events and have built a database that records the 568,736,756 pulse candidates generated during this search. We have discovered a new fast radio burst (FRB), FRB~010305, with a dispersion measure of 350$\pm$5\,\,cm$^{-3}\,$pc and explored why so few FRBs have been discovered in data prior to 2001. After accounting for the dispersion smearing across the channel bandwidth and the sky regions surveyed, the number of FRBs is found to be consistent with model predictions. We also present five single pulse candidates from unknown sources, but with Galactic dispersion measures. We extract a diverse range of sources from the database, which can be used, for example, as a training set of data for new software being developed to search for FRBs in the presence of radio frequency interference. △ Less

Submitted 3 April, 2020; originally announced April 2020.

Comments: 15 pages, 4 Figures, 6 Tables, 2 Appendix, accepted for publication in ApJS

The Parkes Pulsar Timing Array Project: Second data release

Authors: M. Kerr, D. J. Reardon, G. Hobbs, R. M. Shannon, R. N. Manchester, S. Dai, C. J. Russell, S. -B. Zhang, W. van Straten, S. Osłowski, A. Parthasarathy, R. Spiewak, M. Bailes, N. D. R. Bhat, A. D. Cameron, W. A. Coles, J. Dempsey, X. Deng, B. Goncharov, J. F Kaczmarek, M. J. Keith, P. D. Lasky, M. E. Lower, B. Preisig, J. M. Sarkissian , et al. (5 additional authors not shown)

Abstract: We describe 14 years of public data from the Parkes Pulsar Timing Array (PPTA), an ongoing project that is producing precise measurements of pulse times of arrival from 26 millisecond pulsars using the 64-m Parkes radio telescope with a cadence of approximately three weeks in three observing bands. A comprehensive description of the pulsar observing systems employed at the telescope since 2004 is… ▽ More We describe 14 years of public data from the Parkes Pulsar Timing Array (PPTA), an ongoing project that is producing precise measurements of pulse times of arrival from 26 millisecond pulsars using the 64-m Parkes radio telescope with a cadence of approximately three weeks in three observing bands. A comprehensive description of the pulsar observing systems employed at the telescope since 2004 is provided, including the calibration methodology and an analysis of the stability of system components. We attempt to provide full accounting of the reduction from the raw measured Stokes parameters to pulse times of arrival to aid third parties in reproducing our results. This conversion is encapsulated in a processing pipeline designed to track provenance. Our data products include pulse times of arrival for each of the pulsars along with an initial set of pulsar parameters and noise models. The calibrated pulse profiles and timing template profiles are also available. These data represent almost 21,000 hrs of recorded data spanning over 14 years. After accounting for processes that induce time-correlated noise, 22 of the pulsars have weighted root-mean-square timing residuals of < 1 $μ$s in at least one radio band. The data should allow end users to quickly undertake their own gravitational-wave analyses (for example) without having to understand the intricacies of pulsar polarisation calibration or attain a mastery of radio-frequency interference mitigation as is required when analysing raw data files. △ Less

Submitted 21 March, 2020; originally announced March 2020.

Comments: 26 pages, 6 figures. Accepted for publication in PASA. The data release can be downloaded from https://doi.org/10.25919/5db90a8bdeb59

Observational Evidence for Stochastic Shock Drift Acceleration of Electrons at the Earth's Bow Shock

Authors: T. Amano, T. Katou, N. Kitamura, M. Oka, Y. Matsumoto, M. Hoshino, Y. Saito, S. Yokota, B. L. Giles, W. R. Paterson, C. T. Russell, O. Le Contel, R. E. Ergun, P. -A. Lindqvist, D. L. Turner, J. F. Fennell, J. B. Blake

Abstract: The first-order Fermi acceleration of electrons requires an injection of electrons into a mildly relativistic energy range. However, the mechanism of injection has remained a puzzle both in theory and observation. We present direct evidence for a novel stochastic shock drift acceleration theory for the injection obtained with Magnetospheric Multiscale (MMS) observations at Earth's bow shock. The t… ▽ More The first-order Fermi acceleration of electrons requires an injection of electrons into a mildly relativistic energy range. However, the mechanism of injection has remained a puzzle both in theory and observation. We present direct evidence for a novel stochastic shock drift acceleration theory for the injection obtained with Magnetospheric Multiscale (MMS) observations at Earth's bow shock. The theoretical model can explain electron acceleration to mildly relativistic energies at high-speed astrophysical shocks, which may provide a solution to the long-standing issue of electron injection. △ Less

Submitted 17 February, 2020; originally announced February 2020.

Comments: 7 pages, 4 figures. Published in PRL

Journal ref: Phys. Rev. Lett., 124, 065101, 2020

Colliding Winds in and around the Stellar Group IRS 13E at the Galactic Center

Authors: Q. Daniel Wang, Jun Li, Christopher M. P. Russell, Jorge Cuadra

Abstract: IRS~13E is an enigmatic compact group of massive stars located in projection only 3.6 arcseconds away from Sgr A*. This group has been suggested to be bounded by an intermediate-mass black hole (IMBH). We present a multi-wavelength study of the group and its interplay with the environment. Based on Chandra observations, we find the X-ray spectrum of IRS~13E can be well characterized by an opticall… ▽ More IRS~13E is an enigmatic compact group of massive stars located in projection only 3.6 arcseconds away from Sgr A*. This group has been suggested to be bounded by an intermediate-mass black hole (IMBH). We present a multi-wavelength study of the group and its interplay with the environment. Based on Chandra observations, we find the X-ray spectrum of IRS~13E can be well characterized by an optically thin thermal plasma. The emission peaks between two strongly mass-losing Wolf-Rayet stars of the group. These properties can be reasonably well reproduced by simulated colliding winds of these two stars. However, this scenario under-predicts the X-ray intensity in outer regions. The residual emission likely results from the ram-pressure confinement of the IRS~13E group wind by the ambient medium and is apparently associated with a shell-like warm gas structure seen in Pa-alpha and in ALMA observations. These latter observations also show strongly peaked thermal emission with unusually large velocity spread between the two stars. These results indicate that the group is colliding with the bar of the dense cool gas mini-spiral around Sgr A*. The extended X-ray morphology of IRS~13E and its association with the bar further suggest that the group is physically much farther away than the projected distance from Sgr A*. The presence of an IMBH, while favorable to keep the stars bound together, is not necessary to explain the observed stellar and gas properties of IRS~13E. △ Less

Submitted 27 December, 2019; originally announced December 2019.

Comments: 17 pages, accepted for publication in MNRAS

In situ Measurement of Curvature of Magnetic Field in Turbulent Space Plasmas: A Statistical Study

Authors: Riddhi Bandyopadhyay, Yan Yang, William H. Matthaeus, Alexandros Chasapis, Tulasi N. Parashar, Christopher T. Russell, Robert J. Strangeway, Roy B. Torbert, Barbara L. Giles, Daniel J. Gershman, Craig J. Pollock, Thomas E. Moore, James L. Burch

Abstract: Using in situ data, accumulated in the turbulent magnetosheath by the Magnetospheric Multiscale (MMS) Mission, we report a statistical study of magnetic field curvature and discuss its role in the turbulent space plasmas. Consistent with previous simulation results, the Probability Distribution Function (PDF) of the curvature is shown to have distinct power-law tails for both high and low value li… ▽ More Using in situ data, accumulated in the turbulent magnetosheath by the Magnetospheric Multiscale (MMS) Mission, we report a statistical study of magnetic field curvature and discuss its role in the turbulent space plasmas. Consistent with previous simulation results, the Probability Distribution Function (PDF) of the curvature is shown to have distinct power-law tails for both high and low value limits. We find that the magnetic-field-line curvature is intermittently distributed in space. High curvature values reside near weak magnetic-field regions, while low curvature values are correlated with small magnitude of the force acting normal to the field lines. A simple statistical treatment provides an explanation for the observed curvature distribution. This novel statistical characterization of magnetic curvature in space plasma provides a starting point for assessing, in a turbulence context, the applicability and impact of particle energization processes, such as curvature drift, that rely on this fundamental quantity. △ Less

Submitted 29 March, 2020; v1 submitted 19 December, 2019; originally announced December 2019.

Comments: Accepted for Publication in the Astrophysical Journal Letters

Electrostatic turbulence and Debye-scale structures in collisionless shocks

Authors: R. Wang, I. Y. Vasko, F. S. Mozer, S. D. Bale, A. V. Artemyev, J. W. Bonnell, R. Ergun, B. Giles, P. -A. Lindqvist, C. T. Russell, R. Strangeway

Abstract: We present analysis of more than one hundred large-amplitude bipolar electrostatic structures in a quasi-perpendicular supercritical Earth's bow shock crossing, measured by the Magnetospheric Multiscale spacecraft. The occurrence of the bipolar structures is shown to be tightly correlated with magnetic field gradients in the shock transition region. The bipolar structures have negative electrostat… ▽ More We present analysis of more than one hundred large-amplitude bipolar electrostatic structures in a quasi-perpendicular supercritical Earth's bow shock crossing, measured by the Magnetospheric Multiscale spacecraft. The occurrence of the bipolar structures is shown to be tightly correlated with magnetic field gradients in the shock transition region. The bipolar structures have negative electrostatic potentials and spatial scales of a few Debye lengths. The bipolar structures propagate highly oblique to the shock normal with velocities (in the plasma rest frame) of the order of the ion-acoustic velocity. We argue that the bipolar structures are ion phase space holes produced by the two-stream instability between incoming and reflected ions. This is the first identification of the ion two-stream instability in collisionless shocks. The implications for electron acceleration are discussed. △ Less

Submitted 3 December, 2019; originally announced December 2019.