Showing 1–50 of 91 results for author: Clark, C J

The Orbit and Companion of PSR J1622-0315: Variable Asymmetry and a Massive Neutron Star

Authors: Bidisha Sen, Manuel Linares, Mark R. Kennedy, Rene P. Breton, Devina Misra, Marco Turchetta, Vikram S. Dhillon, Daniel Mata Sanchez, Colin J. Clark

Abstract: The companion to PSR J1622-0315, one of the most compact known redback millisecond pulsars, shows extremely low irradiation despite its short orbital period. We model this system to determine the binary parameters, combining optical observations from NTT in 2017 and NOT in 2022 with the binary modeling code ICARUS. We find a best-fit neutron star mass of $2.3 \pm 0.4\,\text{M}_\odot $, and a compa… ▽ More The companion to PSR J1622-0315, one of the most compact known redback millisecond pulsars, shows extremely low irradiation despite its short orbital period. We model this system to determine the binary parameters, combining optical observations from NTT in 2017 and NOT in 2022 with the binary modeling code ICARUS. We find a best-fit neutron star mass of $2.3 \pm 0.4\,\text{M}_\odot $, and a companion mass of $0.15 \pm 0.02\,\text{M}_\odot$. We detect for the first time low-level irradiation from asymmetry in the minima as well as a change in the asymmetry of the maxima of its light curves over five years. Using star spot models, we find better fits than those from symmetric direct heating models, with consistent orbital parameters. We discuss an alternative scenario where the changing asymmetry is produced by a variable intrabinary shock. In summary, we find that PSR J1622-0315 combines low irradiation with variable light curve asymmetry, and a relatively high neutron star mass. △ Less

Submitted 15 July, 2024; originally announced July 2024.

Comments: 18 pages, 11 figures

TRAPUM search for pulsars in supernova remnants and pulsar wind nebulae -- I. Survey description and initial discoveries

Authors: J. D. Turner, B. W. Stappers, E. Carli, E. D. Barr, W. Becker, J. Behrend, R. P. Breton, S. Buchner, M. Burgay, D. J. Champion, W. Chen, C. J. Clark, D. M. Horn, E. F. Keane, M. Kramer, L. K ünkel, L. Levin, Y. P. Men, P. V. Padmanabh, A. Ridolfi, V. Venkatraman Krishnan

Abstract: We present the description and initial results of the TRAPUM (TRAnsients And PUlsars with MeerKAT) search for pulsars associated with supernova remnants (SNRs), pulsar wind nebulae and unidentified TeV emission. The list of sources to be targeted includes a large number of well-known candidate pulsar locations but also new candidate SNRs identified using a range of criteria. Using the 64-dish Meer… ▽ More We present the description and initial results of the TRAPUM (TRAnsients And PUlsars with MeerKAT) search for pulsars associated with supernova remnants (SNRs), pulsar wind nebulae and unidentified TeV emission. The list of sources to be targeted includes a large number of well-known candidate pulsar locations but also new candidate SNRs identified using a range of criteria. Using the 64-dish MeerKAT radio telescope, we use an interferometric beamforming technique to tile the potential pulsar locations with coherent beams which we search for radio pulsations, above a signal-to-noise of 9, down to an average flux density upper limit of 30 $μ$Jy. This limit is target-dependent due to the contribution of the sky and nebula to the system temperature. Coherent beams are arranged to overlap at their 50 per cent power radius, so the sensitivity to pulsars is not degraded by more than this amount, though realistically averages around 65 per cent if every location in the beam is considered. We report the discovery of two new pulsars; PSR J1831$-$0941 is an adolescent pulsar likely to be the plerionic engine of the candidate PWN G20.0+0.0, and PSR J1818$-$1502 appears to be an old and faint pulsar that we serendipitously discovered near the centre of a SNR already hosting a compact central object. The survey holds importance for better understanding of neutron star birth rates and the energetics of young pulsars. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: 17 pages, 6 figures, 3 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

Discovery and timing of ten new millisecond pulsars in the globular cluster Terzan 5

Authors: P. V. Padmanabh, S. M. Ransom, P. C. C. Freire, A. Ridolfi, J. D. Taylor, C. Choza, C. J. Clark, F. Abbate, M. Bailes, E. D. Barr, S. Buchner, M. Burgay, M. E. DeCesar, W. Chen, A. Corongiu, D. J. Champion, A. Dutta, M. Geyer, J. W. T. Hessels, M. Kramer, A. Possenti, I. H. Stairs, B. W. Stappers, V. Venkatraman Krishnan, L. Vleeschower , et al. (1 additional authors not shown)

Abstract: We report the discovery of ten new pulsars in the globular cluster Terzan 5 as part of the Transients and Pulsars with MeerKAT (TRAPUM) Large Survey Project. We observed Terzan 5 at L-band (856--1712 MHz) with the MeerKAT radio telescope for four hours on two epochs, and performed acceleration searches of 45 out of 288 tied-array beams covering the core of the cluster. We obtained phase-connected… ▽ More We report the discovery of ten new pulsars in the globular cluster Terzan 5 as part of the Transients and Pulsars with MeerKAT (TRAPUM) Large Survey Project. We observed Terzan 5 at L-band (856--1712 MHz) with the MeerKAT radio telescope for four hours on two epochs, and performed acceleration searches of 45 out of 288 tied-array beams covering the core of the cluster. We obtained phase-connected timing solutions for nine discoveries, covering nearly two decades of archival observations from the Green Bank Telescope for all but one. Highlights include PSR J1748$-$2446ao which is an eccentric ($e = 0.32$) wide-orbit (orbital period $P_{\rm b} = 57.55$ d) system. We were able to measure the rate of advance of periastron ($\dotω$) for this system allowing us to determine a total mass of $3.17 \pm \, 0.02\, \rm M_{\odot}$. With a minimum companion mass ($M_{\rm c}$) of $\sim 0.8\, \rm M_{\odot}$, PSR J1748$-$2446ao is a candidate double neutron star (DNS) system. If confirmed to be a DNS, it would be the fastest spinning pulsar ($P = 2.27$ ms) and the longest orbital period measured for any known DNS system. PSR J1748$-$2446ap has the second highest eccentricity for any recycled pulsar ($e \sim 0.905$) and for this system we can measure the total mass ($1.997 \pm 0.006\, \rm M_{\odot}$) and also estimate the individual pulsar and companion masses. PSR J1748$-$2446ar is an eclipsing redback (minimum $M_{\rm c} \sim 0.34\, \rm M_{\odot}$) system whose properties confirm it to be the counterpart to a previously published source identified in radio and X-ray imaging. With these discoveries, the total number of confirmed pulsars in Terzan 5 is 49, the highest for any globular cluster so far. These discoveries further enhance the rich set of pulsars known in Terzan 5 and provide scope for a deeper understanding of binary stellar evolution, cluster dynamics and ensemble population studies. △ Less

Submitted 19 June, 2024; v1 submitted 26 March, 2024; originally announced March 2024.

Comments: 23 pages, 11 figures, 5 tables, published in A&A

Journal ref: A&A 686, A166 (2024)

A targeted radio pulsar survey of redback candidates with MeerKAT

Authors: T. Thongmeearkom, C. J. Clark, R. P. Breton, M. Burgay, L. Nieder, P. C. C. Freire, E. D. Barr, B. W. Stappers, S. M. Ransom, S. Buchner, F. Calore, D. J. Champion, I. Cognard, J. -M. Grießmeier, M. Kramer, L. Levin, P. V. Padmanabh, A. Possenti, A. Ridolfi, V. Venkatraman Krishnan, L. Vleeschower

Abstract: Redbacks are millisecond pulsar binaries with low mass, irradiated companions. These systems have a rich phenomenology that can be used to probe binary evolution models, pulsar wind physics, and the neutron star mass distribution. A number of high-confidence redback candidates have been identified through searches for variable optical and X-ray sources within the localisation regions of unidentifi… ▽ More Redbacks are millisecond pulsar binaries with low mass, irradiated companions. These systems have a rich phenomenology that can be used to probe binary evolution models, pulsar wind physics, and the neutron star mass distribution. A number of high-confidence redback candidates have been identified through searches for variable optical and X-ray sources within the localisation regions of unidentified but pulsar-like Fermi-LAT gamma-ray sources. However, these candidates remain unconfirmed until pulsations are detected. As part of the TRAPUM project, we searched for radio pulsations from six of these redback candidates with MeerKAT. We discovered three new radio millisecond pulsars, PSRs J0838$-$2527, J0955$-$3947 and J2333$-$5526, confirming their redback nature. PSR J0838$-$2827 remained undetected for two years after our discovery despite repeated observations, likely due to evaporated material absorbing the radio emission for long periods of time. While, to our knowledge, this system has not undergone a transition to an accreting state, the disappearance, likely caused by extreme eclipses, illustrates the transient nature of spider pulsars and the heavy selection bias in uncovering their radio population. Radio timing enabled the detection of gamma-ray pulsations from all three pulsars, from which we obtained 15-year timing solutions. All of these sources exhibit complex orbital period variations consistent with gravitational quadrupole moment variations in the companion stars. These timing solutions also constrain the binary mass ratios, allowing us to narrow down the pulsar masses. We find that PSR J2333$-$5526 may have a neutron star mass in excess of 2 M$_{\odot}$. △ Less

Submitted 14 March, 2024; originally announced March 2024.

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

A 350-MHz Green Bank Telescope Survey of Unassociated Fermi LAT Sources: Discovery and Timing of Ten Millisecond Pulsars

Authors: P. Bangale, B. Bhattacharyya, F. Camilo, C. J. Clark, I. Cognard, M. E. DeCesar, E. C. Ferrara, P. Gentile, L. Guillemot, J. W. T. Hessels, T. J. Johnson, M. Kerr, M. A. McLaughlin, L. Nieder, S. M. Ransom, P. S. Ray, M. S. E. Roberts, J. Roy, S. Sanpa-Arsa, G. Theureau, M. T. Wolff

Abstract: We have searched for radio pulsations towards 49 Fermi Large Area Telescope (LAT) 1FGL Catalog $γ$-ray sources using the Green Bank Telescope at 350 MHz. We detected 18 millisecond pulsars (MSPs) in blind searches of the data; 10 of these were discoveries unique to our survey. Sixteen are binaries, with eight having short orbital periods $P_B < 1$ day. No radio pulsations from young pulsars were d… ▽ More We have searched for radio pulsations towards 49 Fermi Large Area Telescope (LAT) 1FGL Catalog $γ$-ray sources using the Green Bank Telescope at 350 MHz. We detected 18 millisecond pulsars (MSPs) in blind searches of the data; 10 of these were discoveries unique to our survey. Sixteen are binaries, with eight having short orbital periods $P_B < 1$ day. No radio pulsations from young pulsars were detected, although three targets are coincident with apparently radio-quiet $γ$-ray pulsars discovered in LAT data. Here, we give an overview of the survey and present radio and $γ$-ray timing results for the 10 MSPs discovered. These include the only isolated MSP discovered in our survey and six short-$P_B$ binary MSPs. Of these, three have very low-mass companions ($M_c$ $\ll$ 0.1M$_{\odot}$) and hence belong to the class of black widow pulsars. Two have more massive, non-degenerate companions with extensive radio eclipses and orbitally modulated X-ray emission consistent with the redback class. Significant $γ$-ray pulsations have been detected from nine of the discoveries. This survey and similar efforts suggest that the majority of Galactic $γ$-ray sources at high Galactic latitudes are either MSPs or relatively nearby non-recycled pulsars, with the latter having on average a much smaller radio/$γ$-ray beaming ratio as compared to MSPs. It also confirms that past surveys suffered from an observational bias against finding short-$P_B$ MSP systems. △ Less

Submitted 14 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJ (25 pages, 15 figues, 4 tables)

Journal ref: ApJ, Vol 966, 20 pp. (2024)

Mass estimates from optical modelling of the new TRAPUM redback PSR J1910-5320

Authors: O. G. Dodge, R. P. Breton, C. J. Clark, M. Burgay, J. Strader, K. -Y. Au, E. D. Barr, S. Buchner, V. S. Dhillon, E. C. Ferrara, P. C. C. Freire, J. -M. Griessmeier, M. R. Kennedy, M. Kramer, K. -L. Li, P. V. Padmanabh, A. Phosrisom, B. W. Stappers, S. J. Swihart, T. Thongmeearkom

Abstract: Spider pulsars continue to provide promising candidates for neutron star mass measurements. Here we present the discovery of PSR~J1910$-$5320, a new millisecond pulsar discovered in a MeerKAT observation of an unidentified \textit{Fermi}-LAT gamma-ray source. This pulsar is coincident with a recently identified candidate redback binary, independently discovered through its periodic optical flux an… ▽ More Spider pulsars continue to provide promising candidates for neutron star mass measurements. Here we present the discovery of PSR~J1910$-$5320, a new millisecond pulsar discovered in a MeerKAT observation of an unidentified \textit{Fermi}-LAT gamma-ray source. This pulsar is coincident with a recently identified candidate redback binary, independently discovered through its periodic optical flux and radial velocity. New multi-color optical light curves obtained with ULTRACAM/NTT in combination with MeerKAT timing and updated SOAR/Goodman spectroscopic radial velocity measurements allow a mass constraint for PSR~J1910$-$5320. \texttt{Icarus} optical light curve modelling, with streamlined radial velocity fitting, constrains the orbital inclination and companion velocity, unlocking the binary mass function given the precise radio ephemeris. Our modelling aims to unite the photometric and spectroscopic measurements available by fitting each simultaneously to the same underlying physical model, ensuring self-consistency. This targets centre-of-light radial velocity corrections necessitated by the irradiation endemic to spider systems. Depending on the gravity darkening prescription used, we find a moderate neutron star mass of either $1.6\pm0.2$ or $1.4\pm0.2$ $M_\odot$. The companion mass of either $0.45\pm0.04$ or $0.43^{+0.04}_{-0.03}$ $M_\odot$ also further confirms PSR~J1910$-$5320 as an irradiated redback spider pulsar.radiated redback spider pulsar. △ Less

Submitted 18 January, 2024; originally announced January 2024.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 18 pages, 9 figures

Exploring the interstellar medium of NGC 891 at millimeter wavelengths using the NIKA2 camera

Authors: S. Katsioli, R. Adam, P. Ade, H. Ajeddig, P. André, E. Artis, H. Aussel, M. Baes, A. Beelen, A. Benoît, S. Berta, L. Bing, O. Bourrion, M. Calvo, A. Catalano, C. J. R. Clark, I. De Looze, M. De Petris, F. -X. Désert, S. Doyle, E. F. C. Driessen, G. Ejlali, M. Galametz, F. Galliano, A. Gomez , et al. (39 additional authors not shown)

Abstract: In the framework of the IMEGIN Large Program, we used the NIKA2 camera on the IRAM 30-m telescope to observe the edge-on galaxy NGC 891 at 1.15 mm and 2 mm and at a FWHM of 11.1" and 17.6", respectively. Multiwavelength data enriched with the new NIKA2 observations fitted by the HerBIE SED code (coupled with the THEMIS dust model) were used to constrain the physical properties of the ISM. Emission… ▽ More In the framework of the IMEGIN Large Program, we used the NIKA2 camera on the IRAM 30-m telescope to observe the edge-on galaxy NGC 891 at 1.15 mm and 2 mm and at a FWHM of 11.1" and 17.6", respectively. Multiwavelength data enriched with the new NIKA2 observations fitted by the HerBIE SED code (coupled with the THEMIS dust model) were used to constrain the physical properties of the ISM. Emission originating from the diffuse dust disk is detected at all wavelengths from mid-IR to mm, while mid-IR observations reveal warm dust emission from compact HII regions. Indications of mm excess emission have also been found in the outer parts of the galactic disk. Furthermore, our SED fitting analysis constrained the mass fraction of the small (< 15 Angstrom) dust grains. We found that small grains constitute 9.5% of the total dust mass in the galactic plane, but this fraction increases up to ~ 20% at large distances (|z| > 3 kpc) from the galactic plane. △ Less

Submitted 6 October, 2023; originally announced October 2023.

Comments: To appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

Journal ref: EPJ Web of conferences 293 (2024) 00026

The stratification of ISM properties in the edge-on galaxy NGC 891 revealed by NIKA2

Authors: S. Katsioli, E. M. Xilouris, C. Kramer, R. Adam, P. Ade, H. Ajeddig, P. André, E. Artis, H. Aussel, M. Baes, A. Beelen, A. Benoît, S. Berta, L. Bing, O. Bourrion, M. Calvo, A. Catalano, C. J. R. Clark, I. De Looze, M. De Petris, F. -X. Désert, S. Doyle, E. F. C. Driessen, G. Ejlali, M. Galametz , et al. (38 additional authors not shown)

Abstract: As the millimeter wavelength range remains a largely unexplored spectral region for galaxies, the IMEGIN large program aims to map the millimeter continuum emission of 22 nearby galaxies at 1.15 and 2 mm. Using the high-resolution maps produced by the NIKA2 camera, we explore the existence of very cold dust and take possible contamination by free-free and synchrotron emission into account. We stud… ▽ More As the millimeter wavelength range remains a largely unexplored spectral region for galaxies, the IMEGIN large program aims to map the millimeter continuum emission of 22 nearby galaxies at 1.15 and 2 mm. Using the high-resolution maps produced by the NIKA2 camera, we explore the existence of very cold dust and take possible contamination by free-free and synchrotron emission into account. We study the IR-to-radio emission coming from different regions along the galactic plane and at large vertical distances. New observations of NGC 891, using the NIKA2 camera on the IRAM 30m telescope, along with a suite of observations at other wavelengths were used to perform a multiwavelength study of the spectral energy distribution in the interstellar medium in this galaxy. This analysis was performed globally and locally, using the advanced hierarchical Bayesian fitting code, HerBIE, coupled with the THEMIS dust model. Our dust modeling is able to reproduce the near-IR to millimeter emission of NGC 891, with the exception of an excess at a level of 25% obtained by the NIKA2 observations in the outermost parts of the disk. The radio continuum and thermal dust emission are distributed differently in the disk and galaxy halo. Different dusty environments are also revealed by a multiwavelength investigation of the emission features. Our detailed decomposition at millimeter and centimeter wavelengths shows that emission at 1 mm is purely originated by dust. Radio components become progressively important with increasing wavelengths. Finally, we find that emission arising from small dust grains accounts for ~ 9.5% of the total dust mass, reaching up to 20% at large galactic latitudes. Shock waves in the outflows that shatter the dust grains might explain this higher fraction of small grains in the halo. △ Less

Submitted 15 September, 2023; originally announced September 2023.

Comments: 15 pages, 10 figures, Accepted for publication in A&A

Journal ref: A&A 679, A7 (2023)

PulsarX: a new pulsar searching package -I. A high performance folding program for pulsar surveys

Authors: Yunpeng Men, Ewan Barr, C. J. Clark, Emma Carli, Gregory Desvignes

Abstract: Pulsar surveys with modern radio telescopes are becoming increasingly computationally demanding. This is particularly true for wide field-of-view pulsar surveys with radio interferometers, and those conducted in real or quasi-real time. These demands result in data analysis bottlenecks that can limit the parameter space covered by the surveys and diminish their scientific return. In this paper, we… ▽ More Pulsar surveys with modern radio telescopes are becoming increasingly computationally demanding. This is particularly true for wide field-of-view pulsar surveys with radio interferometers, and those conducted in real or quasi-real time. These demands result in data analysis bottlenecks that can limit the parameter space covered by the surveys and diminish their scientific return. In this paper, we address the computational challenge of `candidate folding' in pulsar searching, presenting a novel, efficient approach designed to optimise the simultaneous folding of large numbers of pulsar candidates. We provide a complete folding pipeline appropriate for large-scale pulsar surveys including radio frequency interference (RFI) mitigation, dedispersion, folding and parameter optimization. By leveraging the Fast Discrete Dispersion Measure Transform (FDMT) algorithm proposed by Zackay et al. (2017), we have developed an optimized, and cache-friendly implementation that we term the pruned FDMT (pFDMT). The pFDMT approach efficiently reuses intermediate processing results and prunes the unused computation paths, resulting in a significant reduction in arithmetic operations. In addition, we propose a novel folding algorithm based on the Tikhonov-regularised least squares method (TLSM) that can improve the time resolution of the pulsar profile. We present the performance of its real-world application as an integral part of two major pulsar search projects conducted with the MeerKAT telescope: the MPIfR-MeerKAT Galactic Plane Survey (MMGPS) and the Transients and Pulsars with MeerKAT (TRAPUM) project. In our processing, for approximately 500 candidates, the theoretical number of dedispersion operations can be reduced by a factor of around 50 when compared to brute-force dedispersion, which scales with the number of candidates. △ Less

Submitted 5 September, 2023; originally announced September 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 679, A20 (2023)

The Third Fermi Large Area Telescope Catalog of Gamma-ray Pulsars

Authors: David A. Smith, Philippe Bruel, Colin J. Clark, Lucas Guillemot, Matthew T. Kerr, Paul Ray, Soheila Abdollahi, Marco Ajello, Luca Baldini, Jean Ballet, Matthew Baring, Cees Bassa, Josefa Becerra Gonzalez, Ronaldo Bellazzini, Alessandra Berretta, Bhaswati Bhattacharyya, Elisabetta Bissaldi, Raffaella Bonino, Eugenio Bottacini, Johan Bregeon, Marta Burgay, Toby Burnett, Rob Cameron, Fernando Camilo, Regina Caputo , et al. (134 additional authors not shown)

Abstract: We present 294 pulsars found in GeV data from the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope. Another 33 millisecond pulsars (MSPs) discovered in deep radio searches of LAT sources will likely reveal pulsations once phase-connected rotation ephemerides are achieved. A further dozen optical and/or X-ray binary systems co-located with LAT sources also likely harbor gamma-ray M… ▽ More We present 294 pulsars found in GeV data from the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope. Another 33 millisecond pulsars (MSPs) discovered in deep radio searches of LAT sources will likely reveal pulsations once phase-connected rotation ephemerides are achieved. A further dozen optical and/or X-ray binary systems co-located with LAT sources also likely harbor gamma-ray MSPs. This catalog thus reports roughly 340 gamma-ray pulsars and candidates, 10% of all known pulsars, compared to $\leq 11$ known before Fermi. Half of the gamma-ray pulsars are young. Of these, the half that are undetected in radio have a broader Galactic latitude distribution than the young radio-loud pulsars. The others are MSPs, with 6 undetected in radio. Overall, >235 are bright enough above 50 MeV to fit the pulse profile, the energy spectrum, or both. For the common two-peaked profiles, the gamma-ray peak closest to the magnetic pole crossing generally has a softer spectrum. The spectral energy distributions tend to narrow as the spindown power $\dot E$ decreases to its observed minimum near $10^{33}$ erg s$^{-1}$, approaching the shape for synchrotron radiation from monoenergetic electrons. We calculate gamma-ray luminosities when distances are available. Our all-sky gamma-ray sensitivity map is useful for population syntheses. The electronic catalog version provides gamma-ray pulsar ephemerides, properties and fit results to guide and be compared with modeling results. △ Less

Submitted 20 July, 2023; originally announced July 2023.

Comments: 142 pages. Accepted by the Astrophysical Journal Supplement

The MPIfR-MeerKAT Galactic Plane survey I -- System setup and early results

Authors: P. V. Padmanabh, E. D. Barr, S. S. Sridhar, M. R. Rugel, A. Damas-Segovia, A. M. Jacob, V. Balakrishnan, M. Berezina, M. C. i Bernadich, A. Brunthaler, D. J. Champion, P. C. C. Freire, S. Khan, H. -R. Klöckner, M. Kramer, Y. K. Ma, S. A. Mao, Y. P. Men, K. M. Menten, S. Sengupta, V. Venkatraman Krishnan, O. Wucknitz, F. Wyrowski, M. C. Bezuidenhout, S. Buchner , et al. (8 additional authors not shown)

Abstract: Galactic plane radio surveys play a key role in improving our understanding of a wide range of astrophysical phenomena. Performing such a survey using the latest interferometric telescopes produces large data rates necessitating a shift towards fully or quasi-real-time data analysis with data being stored for only the time required to process them. We present here the overview and setup for the 30… ▽ More Galactic plane radio surveys play a key role in improving our understanding of a wide range of astrophysical phenomena. Performing such a survey using the latest interferometric telescopes produces large data rates necessitating a shift towards fully or quasi-real-time data analysis with data being stored for only the time required to process them. We present here the overview and setup for the 3000 hour Max-Planck-Institut fuer Radioastronomie (MPIfR) MeerKAT Galactic Plane survey (MMGPS). The survey is unique by operating in a commensal mode, addressing key science objectives of the survey including the discovery of new pulsars and transients as well as studies of Galactic magnetism, the interstellar medium and star formation rates. We explain the strategy coupled with the necessary hardware and software infrastructure needed for data reduction in the imaging, spectral and time domains. We have so far discovered 78 new pulsars including 17 confirmed binary systems of which two are potential double neutron star systems. We have also developed an imaging pipeline sensitive to the order of a few tens of micro-Jansky with a spatial resolution of a few arcseconds. Further science operations with an in-house built S-Band receiver operating between 1.7-3.5 GHz are about to commence. Early spectral line commissioning observations conducted at S-Band, targeting transitions of the key molecular gas tracer CH at 3.3 GHz already illustrate the spectroscopic capabilities of this instrument. These results lay a strong foundation for future surveys with telescopes like the Square Kilometre Array (SKA). △ Less

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

Comments: 25 pages, 10 figures, Accepted in MNRAS

Tied-Array Beam Localisation of Radio Transients and Pulsars

Authors: M. C. Bezuidenhout, C. J. Clark, R. P. Breton, B. W. Stappers, E. D. Barr, M. Caleb, W. Chen, F. Jankowski, M. Kramer, K. Rajwade, M. Surnis

Abstract: Multi-element interferometers such as MeerKAT, which observe with high time resolution and have a wide field-of-view, provide an ideal opportunity to perform real-time, untargeted transient and pulsar searches. However, because of data storage limitations, it is not always feasible to store the baseband data required to image the field of a discovered transient or pulsar. This limits the ability o… ▽ More Multi-element interferometers such as MeerKAT, which observe with high time resolution and have a wide field-of-view, provide an ideal opportunity to perform real-time, untargeted transient and pulsar searches. However, because of data storage limitations, it is not always feasible to store the baseband data required to image the field of a discovered transient or pulsar. This limits the ability of surveys to effectively localise their discoveries and may restrict opportunities for follow-up science, especially of one-off events like some Fast Radio Bursts (FRBs). Here we present a novel maximum-likelihood estimation approach to localising transients and pulsars detected in multiple MeerKAT tied-array beams at once, which we call Tied Array Beam Localisation (TABLo), as well as a Python implementation of the method named SeeKAT. We provide real-world examples of SeeKAT's use as well as a Monte Carlo analysis to show that it is capable of localising single pulses detected in beamformed MeerKAT data to (sub-)arcsecond precision. △ Less

Submitted 20 February, 2023; originally announced February 2023.

Comments: 15 pages, 13 figures

The Quest for the Missing Dust: II -- Two Orders of Magnitude of Evolution in the Dust-to-Gas Ratio Resolved Within Local Group Galaxies

Authors: Christopher J. R. Clark, Julia C. Roman-Duval, Karl D. Gordon, Caroline Bot, Matthew W. L. Smith, Lea M. Z. Hagen

Abstract: We explore evolution in the dust-to-gas ratio with density within four well-resolved Local Group galaxies - the LMC, SMC, M31, and M33. We do this using new ${\it Herschel}$ maps, which restore extended emission that was missed by previous ${\it Herschel}$ reductions. This improved data allows us to probe the dust-to-gas ratio across 2.5 orders of magnitude in ISM surface density. We find signific… ▽ More We explore evolution in the dust-to-gas ratio with density within four well-resolved Local Group galaxies - the LMC, SMC, M31, and M33. We do this using new ${\it Herschel}$ maps, which restore extended emission that was missed by previous ${\it Herschel}$ reductions. This improved data allows us to probe the dust-to-gas ratio across 2.5 orders of magnitude in ISM surface density. We find significant evolution in the dust-to-gas ratio, with dust-to-gas varying with density within each galaxy by up to a factor 22.4. We explore several possible reasons for this, and our favored explanation is dust grain growth in denser regions of ISM. We find that the evolution of the dust-to-gas ratio with ISM surface density is very similar between M31 and M33, despite their large differences in mass, metallicity, and star formation rate; conversely, we find M33 and the LMC to have very different dust-to-gas evolution profiles, despite their close similarity in those properties. Our dust-to-gas ratios address previous disagreement between UV- and FIR-based dust-to-gas estimates for the Magellanic Clouds, removing the disagreement for the LMC, and considerably reducing it for the SMC - with our new dust-to-gas measurements being factors of 2.4 and 2.0 greater than the previous far-infrared estimates, respectively. We also observe that the dust-to-gas ratio appears to fall at the highest densities for the LMC, M31, and M33; this is unlikely to be an actual physical phenomenon, and we posit that it may be due to a combined effect of dark gas, and changing dust mass opacity. △ Less

Submitted 14 February, 2023; originally announced February 2023.

Comments: Accepted for publication in the Astrophysical Journal

Neutron star mass estimates from gamma-ray eclipses in spider millisecond pulsar binaries

Authors: C. J. Clark, M. Kerr, E. D. Barr, B. Bhattacharyya, R. P. Breton, P. Bruel, F. Camilo, W. Chen, I. Cognard, H. T. Cromartie, J. Deneva, V. S. Dhillon, L. Guillemot, M. R. Kennedy, M. Kramer, A. G. Lyne, D. Mata Sánchez, L. Nieder, C. Phillips, S. M. Ransom, P. S. Ray, M. S. E. Roberts, J. Roy, D. A. Smith, R. Spiewak , et al. (4 additional authors not shown)

Abstract: Reliable neutron star mass measurements are key to determining the equation-of-state of cold nuclear matter, but these are rare. "Black Widows" and "Redbacks" are compact binaries consisting of millisecond pulsars and semi-degenerate companion stars. Spectroscopy of the optically bright companions can determine their radial velocities, providing inclination-dependent pulsar mass estimates. While i… ▽ More Reliable neutron star mass measurements are key to determining the equation-of-state of cold nuclear matter, but these are rare. "Black Widows" and "Redbacks" are compact binaries consisting of millisecond pulsars and semi-degenerate companion stars. Spectroscopy of the optically bright companions can determine their radial velocities, providing inclination-dependent pulsar mass estimates. While inclinations can be inferred from subtle features in optical light curves, such estimates may be systematically biased due to incomplete heating models and poorly-understood variability. Using data from the Fermi Large Area Telescope, we have searched for gamma-ray eclipses from 49 spider systems, discovering significant eclipses in 7 systems, including the prototypical black widow PSR B1957$+$20. Gamma-ray eclipses require direct occultation of the pulsar by the companion, and so the detection, or significant exclusion, of a gamma-ray eclipse strictly limits the binary inclination angle, providing new robust, model-independent pulsar mass constraints. For PSR B1957$+$20, the eclipse implies a much lighter pulsar ($M_{\rm psr} = 1.81 \pm 0.07\,M_{\odot}$) than inferred from optical light curve modelling. △ Less

Submitted 26 January, 2023; originally announced January 2023.

Comments: 31 pages, 4 figures, includes supplementary tables; published in Nature Astronomy

A black widow population dissection through HiPERCAM multi-band light curve modelling

Authors: D. Mata Sánchez, M. R. Kennedy, C. J. Clark, R. P. Breton, V. S. Dhillon, G. Voisin, F. Camilo, S. Littlefair, T. R. Marsh, J. Stringer

Abstract: Black widows are extreme millisecond pulsar binaries where the pulsar wind ablates their low-mass companion stars. Their optical light curves vary periodically due to the high irradiation and tidal distortion of the companion, which allows us to infer the binary parameters. We present simultaneous multi-band observations obtained with the HIPERCAM instrument at the 10.4-m GTC telescope for six of… ▽ More Black widows are extreme millisecond pulsar binaries where the pulsar wind ablates their low-mass companion stars. Their optical light curves vary periodically due to the high irradiation and tidal distortion of the companion, which allows us to infer the binary parameters. We present simultaneous multi-band observations obtained with the HIPERCAM instrument at the 10.4-m GTC telescope for six of these systems. The combination of this five-band fast photometer with the world's largest optical telescope enables us to inspect the light curve range near minima. We present the first light curve for PSR J1641+8049, as well as attain a significant increase in signal-to-noise and cadence compared with previous publications for the remaining 5 targets: PSR J0023+0923, PSR J0251+2606, PSR J0636+5129, PSR J0952-0607 and PSR J1544+4937. We report on the results of the light curve modelling with the Icarus code for all six systems, which reveals some of the hottest and densest companion stars known. We compare the parameters derived with the limited but steadily growing black widow population for which optical modelling is available. We find some expected correlations, such as that between the companion star mean density and the orbital period of the system, but also a puzzling positive correlation between the orbital inclination and the irradiation temperature of the companion. We propose such a correlation would arise if pulsars with magnetic axis orthogonal to their spin axis are capable of irradiating their companions to a higher degree. △ Less

Submitted 17 January, 2023; originally announced January 2023.

Comments: 18 pages (+12 pages for appendix), 12 figures (+13 in the appendix), 3 tables (1 in the appendix). Accepted for publication in MNRAS on 2023 January 17th

A rapid optical and X-ray timing study of the neutron star X-ray binary Swift J1858.6-0814

Authors: T. Shahbaz, J. A. Paice, K. M. Rajwade, A. Veledina, P. Gandhi., V. S. Dhillon, T. R. Marsh, S. Littlefair, M. R. Kennedy, R. P. Breton, C. J. Clark

Abstract: We present a rapid timing analysis of optical (HiPERCAM and ULTRACAM) and X-ray (NICER) observations of the X-ray transient Swift J1858.6-0814 during 2018 and 2019. The optical light curves show relatively slow, large amplitude (~1 mags in g$_s$) `blue' flares (i.e. stronger at shorter wavelengths) on time-scales of ~minutes as well as fast, small amplitude (~0.1 mag in g$_s$) `red' flares (i.e. s… ▽ More We present a rapid timing analysis of optical (HiPERCAM and ULTRACAM) and X-ray (NICER) observations of the X-ray transient Swift J1858.6-0814 during 2018 and 2019. The optical light curves show relatively slow, large amplitude (~1 mags in g$_s$) `blue' flares (i.e. stronger at shorter wavelengths) on time-scales of ~minutes as well as fast, small amplitude (~0.1 mag in g$_s$) `red' flares (i.e. stronger at longer wavelengths) on time-scales of ~seconds. The `blue' and `red' flares are consistent with X-ray reprocessing and optically thin synchrotron emission, respectively, similar to what is observed in other X-ray binaries. The simultaneous optical versus soft- and hard-band X-ray light curves show time- and energy dependent correlations. The 2019 March 4 and parts of the June data show a nearly symmetric positive cross correlations (CCFs) at positive lags consistent with simple X-ray disc reprocessing. The soft- and hard-band CCFs are similar and can be reproduced if disc reprocessing dominates in the optical and one component (disc or synchrotron Comptonization) dominates both the soft and hard X-rays. A part of the 2019 June data shows a very different CCFs. The observed positive correlation at negative lag in the soft-band can be reproduced if the optical synchrotron emission is correlated with the hot flow X-ray emission. The observed timing properties are in qualitative agreement with the hybrid inner hot accretion flow model, where the relative role of the different X-ray and optical components that vary during the course of the outburst, as well as on shorter time-scales, govern the shape of the optical/X-ray CCFs. △ Less

Submitted 16 January, 2023; originally announced January 2023.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 18 pages, 14 figures

The TRAPUM L-band survey for pulsars in Fermi-LAT gamma-ray sources

Authors: C. J. Clark, R. P. Breton, E. D. Barr, M. Burgay, T. Thongmeearkom, L. Nieder, S. Buchner, B. Stappers, M. Kramer, W. Becker, M. Mayer, A. Phosrisom, A. Ashok, M. C. Bezuidenhout, F. Calore, I. Cognard, P. C. C. Freire, M. Geyer, J. -M. Grießmeier, R. Karuppusamy, L. Levin, P. V. Padmanabh, A. Possenti, S. Ransom, M. Serylak , et al. (13 additional authors not shown)

Abstract: More than 100 millisecond pulsars (MSPs) have been discovered in radio observations of gamma-ray sources detected by the Fermi Large Area Telescope (LAT), but hundreds of pulsar-like sources remain unidentified. Here we present the first results from the targeted survey of Fermi-LAT sources being performed by the Transients and Pulsars with MeerKAT (TRAPUM) Large Survey Project. We observed 79 sou… ▽ More More than 100 millisecond pulsars (MSPs) have been discovered in radio observations of gamma-ray sources detected by the Fermi Large Area Telescope (LAT), but hundreds of pulsar-like sources remain unidentified. Here we present the first results from the targeted survey of Fermi-LAT sources being performed by the Transients and Pulsars with MeerKAT (TRAPUM) Large Survey Project. We observed 79 sources identified as possible gamma-ray pulsar candidates by a Random Forest classification of unassociated sources from the 4FGL catalogue. Each source was observed for 10 minutes on two separate epochs using MeerKAT's L-band receiver (856-1712 MHz), with typical pulsed flux density sensitivities of $\sim$100$\,μ$Jy. Nine new MSPs were discovered, eight of which are in binary systems, including two eclipsing redbacks and one system, PSR J1526$-$2744, that appears to have a white dwarf companion in an unusually compact 5 hr orbit. We obtained phase-connected timing solutions for two of these MSPs, enabling the detection of gamma-ray pulsations in the Fermi-LAT data. A follow-up search for continuous gravitational waves from PSR J1526$-$2744 in Advanced LIGO data using the resulting Fermi-LAT timing ephemeris yielded no detection, but sets an upper limit on the neutron star ellipticity of $2.45\times10^{-8}$. We also detected X-ray emission from the redback PSR J1803$-$6707 in data from the first eROSITA all-sky survey, likely due to emission from an intra-binary shock. △ Less

Submitted 16 December, 2022; originally announced December 2022.

Comments: 17 pages, 8 figures, accepted for publication in MNRAS

Journal ref: MNRAS, 519, 5590-5606 (2023)

The resolved scaling relations in DustPedia: Zooming in on the local Universe

Authors: Viviana Casasola, Simone Bianchi, Laura Magrini, Aleksandr V. Mosenkov, Francesco Salvestrini, Maarten Baes, Francesco Calura, Letizia P. Cassara', Christopher J. R. Clark, Edvige Corbelli, Jacopo Fritz, Frederic Galliano, Elisabetta Liuzzo, Suzanne Madden, Angelos Nersesian, Francesca Pozzi, Sambit Roychowdhury, Ivano Baronchelli, Matteo Bonato, Carlotta Gruppioni, Lara Pantoni

Abstract: We perform a homogeneous analysis of an unprecedented set of spatially resolved scaling relations (SRs) between ISM components and other properties in the range of scales 0.3-3.4 kpc. We also study some ratios: dust-to-stellar, dust-to-gas, and dust-to-metal. We use a sample of 18 large, spiral, face-on DustPedia galaxies. All the SRs are moderate/strong correlations except the dust-HI SR that doe… ▽ More We perform a homogeneous analysis of an unprecedented set of spatially resolved scaling relations (SRs) between ISM components and other properties in the range of scales 0.3-3.4 kpc. We also study some ratios: dust-to-stellar, dust-to-gas, and dust-to-metal. We use a sample of 18 large, spiral, face-on DustPedia galaxies. All the SRs are moderate/strong correlations except the dust-HI SR that does not exist or is weak for most galaxies. The SRs do not have a universal form but each galaxy is characterized by distinct correlations, affected by local processes and galaxy peculiarities. The SRs hold starting from 0.3 kpc, and if a breaking down scale exists it is < 0.3 kpc. By evaluating all galaxies at 3.4 kpc, differences due to peculiarities of individual galaxies are cancelled out and the corresponding SRs are consistent with those of whole galaxies. By comparing subgalactic and global scales, the most striking result emerges from the SRs involving ISM components: the dust-total gas SR is a good correlation at all scales, while the dust-H2 and dust-HI SRs are good correlations at subkpc/kpc and total scales, respectively. For the other explored SRs, there is a good agreement between small and global scales and this may support the picture where the main physical processes regulating the properties and evolution of galaxies occur locally. Our results are consistent with the hypothesis of self-regulation of the SF process. The analysis of subgalactic ratios shows that they are consistent with those derived for whole galaxies, from low to high z, supporting the idea that also these ratios could be set by local processes. Our results highlight the heterogeneity of galaxy properties and the importance of resolved studies on local galaxies in the context of galaxy evolution. They also provide observational constraints to theoretical models and updated references for high-z studies. △ Less

Submitted 28 October, 2022; originally announced October 2022.

Comments: 42 pages, 11 figures and 5 tables in the main text, 2 figures and 1 table in Appendix. Accepted for publication in A&A

Journal ref: A&A 668, A130 (2022)

Multi-colour optical light curves of the companion star to the millisecond pulsar PSR J2051-0827

Authors: V. S. Dhillon, M. R. Kennedy, R. P. Breton, C. J. Clark, D. Mata Sánchez, G. Voisin, E. Breedt, A. J. Brown, M. J. Dyer, M. J. Green, P. Kerry, S. P. Littlefair, T. R. Marsh, S. G. Parsons, I. Pelisoli, D. I. Sahman, J. F. Wild, M. H. van Kerkwijk, B. W. Stappers

Abstract: We present simultaneous, multi-colour optical light curves of the companion star to the black-widow pulsar PSR J2051-0827, obtained approximately 10 years apart using ULTRACAM and HiPERCAM, respectively. The ULTRACAM light curves confirm the previously reported asymmetry in which the leading hemisphere of the companion star appears to be brighter than the trailing hemisphere. The HiPERCAM light cu… ▽ More We present simultaneous, multi-colour optical light curves of the companion star to the black-widow pulsar PSR J2051-0827, obtained approximately 10 years apart using ULTRACAM and HiPERCAM, respectively. The ULTRACAM light curves confirm the previously reported asymmetry in which the leading hemisphere of the companion star appears to be brighter than the trailing hemisphere. The HiPERCAM light curves, however, do not show this asymmetry, demonstrating that whatever mechanism is responsible for it varies on timescales of a decade or less. We fit the symmetrical HiPERCAM light curves with a direct-heating model to derive the system parameters, finding an orbital inclination of $55.9^{+4.8}_{-4.1}$ degrees, in good agreement with radio-eclipse constraints. We find that approximately half of the pulsar's spin-down energy is converted to optical luminosity, resulting in temperatures ranging from approximately $5150^{+190}_{-190}$ K on the day side to $2750^{+130}_{-150}$ K on the night side of the companion star. The companion star is close to filling its Roche lobe ($f_{\rm RL} =0.88^{+0.02}_{-0.02}$) and has a mass of $0.039^{+0.010}_{-0.011}$ M$_{\odot}$, giving a mean density of $20.24^{+0.59}_{-0.44}$ g cm$^{-3}$ and an apsidal motion constant in the range $0.0036 < k_2 < 0.0047$. The companion mass and mean density values are consistent with those of brown dwarfs, but the apsidal motion constant implies a significantly more centrally-condensed internal structure than is typical for such objects. △ Less

Submitted 19 August, 2022; originally announced August 2022.

Comments: 10 pages, 4 figures, 3 tables, accepted for publication in Monthly Notices of the Royal Astronomical Society

Dust emissivity in resolved spiral galaxies

Authors: S. Bianchi, V. Casasola, E. Corbelli, F. Galliano, L. Magrini, A. Nersesian, F. Salvestrini, M. Baes, L. P. Cassara, C. J. R. Clark, I. De Looze, A. P. Jones, S. C. Madden, A. Mosenkov, N. Ysard

Abstract: Context: The far-infrared (FIR) and sub-millimeter (submm) emissivity of the Milky Way (MW) cirrus is an important benchmark for dust grain models. Dust masses in other galaxies are generally derived from the FIR/submm using the emission properties of these MW-calibrated models. Aims: We seek to derive the FIR/submm emissivity in nine nearby spiral galaxies to check its compatibility with MW cirru… ▽ More Context: The far-infrared (FIR) and sub-millimeter (submm) emissivity of the Milky Way (MW) cirrus is an important benchmark for dust grain models. Dust masses in other galaxies are generally derived from the FIR/submm using the emission properties of these MW-calibrated models. Aims: We seek to derive the FIR/submm emissivity in nine nearby spiral galaxies to check its compatibility with MW cirrus measurements. Methods: We obtained values of the emissivity at 70 to 500 um, using maps of dust emission from the Herschel satellite and of gas surface density from the THINGS and HERACLES surveys on a scale generally corresponding to 440 pc. We studied the variation of the emissivity with the surface brightness ratio I(250um)/I(500um), a proxy for the intensity of the interstellar radiation field heating the dust. Results: We find that the average value of the emissivity agrees with MW estimates for pixels sharing the same color as the cirrus, namely, for I(250um)/I(500um) = 4.5. For I(250um)/I(500um) > 5, the measured emissivity is instead up to a factor ~2 lower than predicted from MW dust models heated by stronger radiation fields. Regions with higher I(250um)/I(500um) are preferentially closer to the galactic center and have a higher overall (stellar+gas) surface density and molecular fraction. The results do not depend strongly on the adopted CO-to-molecular conversion factor and do not appear to be affected by the mixing of heating conditions. Conclusions: Our results confirm the validity of MW dust models at low density, but are at odds with predictions for grain evolution in higher density environments. If the lower-than-expected emissivity at high I(250um)/I(500um) is the result of intrinsic variations in the dust properties, it would imply an underestimation of the dust mass surface density of up to a factor ~2 when using current dust models. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: 16 pages, 11 figures, A&A accepted

Journal ref: A&A 664, A187 (2022)

METAL: The Metal Evolution, Transport, and Abundance in the Large Magellanic Cloud Hubble program. IV. Calibration of Dust Depletions vs Abundance Ratios in the Milky Way and Magellanic Clouds and Application to Damped Lyman-alpha Systems

Authors: Julia Roman-Duval, Edward B. Jenkins, Kirill Tchernyshyov, Christopher J. R. Clark, Annalisa De Cia, Karl D. Gordon, Aleksandra Hamanowicz, Vianney Lebouteiller, Marc Rafelski, Karin Sandstrom, Jessica Werk, Petia Yanchulova Merica-Jones

Abstract: The evolution of the metal content of the universe can be tracked through rest-frame UV spectroscopy of damped Ly-$α$ systems (DLAs). Gas-phase abundances in DLAs must be corrected for dust depletion effects, which can be accomplished by calibrating the relation between abundance ratios such as [Zn/Fe] and depletions (the fraction of metals in gas, as opposed to dust). Using samples of gas-phase a… ▽ More The evolution of the metal content of the universe can be tracked through rest-frame UV spectroscopy of damped Ly-$α$ systems (DLAs). Gas-phase abundances in DLAs must be corrected for dust depletion effects, which can be accomplished by calibrating the relation between abundance ratios such as [Zn/Fe] and depletions (the fraction of metals in gas, as opposed to dust). Using samples of gas-phase abundances and depletions in the Milky Way (MW), LMC, and SMC, we demonstrate that the relation between [Zn/Fe] and other abundance ratios does not change significantly between these local galaxies and DLAs, indicating that [Zn/Fe] should trace depletions of heavy elements in those systems. The availability of photospheric abundances in young massive stars, a proxy for the total (gas+dust) metallicity of neutral gas, in the MW LMC, and SMC allows us to calibrate the relation between [Zn/Fe] and depletions in these nearby galaxies. We apply the local calibrations of depletions to DLA systems. We find that the fraction of metals in dust, the dust-to-gas-ratio, and total abundances are 2-5 times lower than inferred from previous depletion calibrations based on MW measurements and a different formalism. However, the trend of dust abundance vs. metallicity remains only slightly sub-linear for all existing depletion calibrations, contrary to what is inferred from FIR, 21 cm, and CO emission in nearby galaxies and predicted by chemical evolution models. Observational constraints on the FIR dust opacity and depletions at metallicities lower than 20\% solar will be needed to resolve this tension. △ Less

Submitted 7 June, 2022; originally announced June 2022.

Comments: 37 pages; 22 figures

A Gamma-ray Pulsar Timing Array Constrains the Nanohertz Gravitational Wave Background

Authors: M. Ajello, W. B. Atwood, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, A. Berretta, B. Bhattacharyya, E. Bissaldi, R. D. Blandford, E. Bloom, R. Bonino, P. Bruel, R. Buehler, E. Burns, S. Buson, R. A. Cameron, P. A. Caraveo, E. Cavazzuti, N. Cibrario, S. Ciprini, C. J. Clark, I. Cognard, J. Coronado-Blázquez , et al. (107 additional authors not shown)

Abstract: After large galaxies merge, their central supermassive black holes are expected to form binary systems whose orbital motion generates a gravitational wave background (GWB) at nanohertz frequencies. Searches for this background utilize pulsar timing arrays, which perform long-term monitoring of millisecond pulsars (MSPs) at radio wavelengths. We use 12.5 years of Fermi Large Area Telescope data to… ▽ More After large galaxies merge, their central supermassive black holes are expected to form binary systems whose orbital motion generates a gravitational wave background (GWB) at nanohertz frequencies. Searches for this background utilize pulsar timing arrays, which perform long-term monitoring of millisecond pulsars (MSPs) at radio wavelengths. We use 12.5 years of Fermi Large Area Telescope data to form a gamma-ray pulsar timing array. Results from 35 bright gamma-ray pulsars place a 95\% credible limit on the GWB characteristic strain of $1.0\times10^{-14}$ at 1 yr$^{-1}$, which scales as the observing time span $t_{\mathrm{obs}}^{-13/6}$. This direct measurement provides an independent probe of the GWB while offering a check on radio noise models. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Comments: 3 figures in the main text. 3 figures and 8 tables are in the supplementary material

Irradiated but not eclipsed, the case of PSR J0610-2100

Authors: E. van der Wateren, C. G. Bassa, C. J. Clark, R. P. Breton, I. Cognard, L. Guillemot, G. H. Janssen, A. G. Lyne, B. W. Stappers, G. Theureau

Abstract: We report on radio timing observations of the black widow binary pulsar J0610$-$2100 and optical observations of its binary companion. The radio timing observations extend the timing baseline to 16yr and reveal a marginal detection of the orbital period derivative, but they show no significant evidence of orbital variations such as those seen in other black widow pulsars. Furthermore, no eclipses… ▽ More We report on radio timing observations of the black widow binary pulsar J0610$-$2100 and optical observations of its binary companion. The radio timing observations extend the timing baseline to 16yr and reveal a marginal detection of the orbital period derivative, but they show no significant evidence of orbital variations such as those seen in other black widow pulsars. Furthermore, no eclipses are seen in the observations at observing frequencies ranging from 310 to 2700MHz. The optical VRI light curves were modulated with the orbital period, reaching maximum brightness of $V=26.8$, $R=25.4$, and $I=23.8$ at superior conjunction of the companion, confirming irradiation of the companion by the pulsar. Modelling the light curves indicates that the companion is likely not filling its Roche lobe, while having a moderate inclination ($i > 54\degr$). We find an unusually low temperature and a low irradiation for the irradiated hemisphere of the companion. We investigate the absence of radio eclipses in PSR J0610-2100 and in other black widow systems in relation to their binary, pulsar, and companion properties. We also discuss the suitability of PSR J0610-2100 for pulsar timing array observations aimed at detecting nano-Hertz gravitational waves. △ Less

Submitted 14 March, 2022; v1 submitted 11 March, 2022; originally announced March 2022.

Comments: 9 pages, 5 figures, accepted for publication by A&A

Journal ref: A&A 661, A57 (2022)

Is the black-widow pulsar PSR J1555-2908 in a hierarchical triple system?

Authors: L. Nieder, M. Kerr, C. J. Clark, P. Bruel, H. T. Cromartie, S. M. Ransom, P. S. Ray

Abstract: The 559 Hz black-widow pulsar PSR J1555-2908, originally discovered in radio, is also a bright gamma-ray pulsar. Timing its pulsations using 12 yr of Fermi-LAT gamma-ray data reveals long-term variations in its spin frequency that are much larger than is observed from other millisecond pulsars. While this variability in the pulsar rotation rate could be intrinsic "timing noise", here we consider a… ▽ More The 559 Hz black-widow pulsar PSR J1555-2908, originally discovered in radio, is also a bright gamma-ray pulsar. Timing its pulsations using 12 yr of Fermi-LAT gamma-ray data reveals long-term variations in its spin frequency that are much larger than is observed from other millisecond pulsars. While this variability in the pulsar rotation rate could be intrinsic "timing noise", here we consider an alternative explanation: the variations arise from the presence of a very-low-mass third object in a wide multi-year orbit around the neutron star and its low-mass companion. With current data, this hierarchical-triple-system model describes the pulsar's rotation slightly more accurately than the best-fitting timing-noise model. Future observations will show if this alternative explanation is correct. △ Less

Submitted 11 February, 2022; originally announced February 2022.

Comments: 6 pages, 2 figures, submitted to ApJL

Measuring the mass of the black widow PSR J1555-2908

Authors: M. R. Kennedy, R. P. Breton, C. J. Clark, D. Mata-Sanchez, G. Voisin, V. S. Dhillon, J. P. Halpern, T. R. Marsh, L. Nieder, P. S. Ray, M. H. van Kerkwijk

Abstract: Accurate measurements of the masses of neutron stars are necessary to test binary evolution models, and to constrain the neutron star equation of state. In pulsar binaries with no measurable post-Keplerian parameters, this requires an accurate estimate of the binary system's inclination and the radial velocity of the companion star by other means than pulsar timing. In this paper, we present the r… ▽ More Accurate measurements of the masses of neutron stars are necessary to test binary evolution models, and to constrain the neutron star equation of state. In pulsar binaries with no measurable post-Keplerian parameters, this requires an accurate estimate of the binary system's inclination and the radial velocity of the companion star by other means than pulsar timing. In this paper, we present the results of a new method for measuring this radial velocity using the binary synthesis code Icarus. This method relies on constructing a model spectrum of a tidally distorted, irradiated star as viewed for a given binary configuration. This method is applied to optical spectra of the newly discovered black widow PSR J1555-2908. By modelling the optical spectroscopy alongside optical photometry, we find that the radial velocity of the companion star is $397\pm4$ km s$^{-1}$ (errors quoted at 95\% confidence interval), as well as a binary inclination of $>75^{\rm o}$. Combined with $γ$-ray pulsation timing information, this gives a neutron star mass of 1.67$^{+0.15}_{-0.09}$ M$_\odot$ and a companion mass of 0.060$^{+0.005}_{-0.003}$ M$_\odot$, placing PSR J1555-2908 at the observed upper limit of what is considered a black widow system. △ Less

Submitted 9 February, 2022; originally announced February 2022.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 15 pages, 7 Figures. Underlying data available at https://zenodo.org/record/5653061

Discovery, Timing, and Multiwavelength Observations of the Black Widow Millisecond Pulsar PSR J1555-2908

Authors: Paul S. Ray, Lars Nieder, Colin J. Clark, Scott M. Ransom, H. Thankful Cromartie, Dale A. Frail, Kunal P. Mooley, Huib Intema, Preshanth Jagannathan, Paul Demorest, Kevin Stovall, Jules P. Halpern, Julia Deneva, Sebastien Guillot, Matthew Kerr, Samuel J. Swihart, Philippe Bruel, Ben W. Stappers, Andrew Lyne, Mitch Mickaliger, Fernando Camilo, Elizabeth C. Ferrara, Michael T. Wolff, P. F. Michelson

Abstract: We report the discovery of PSR J1555-2908, a 1.79 ms radio and gamma-ray pulsar in a 5.6 hr binary system with a minimum companion mass of 0.052 $M_\odot$. This fast and energetic ($\dot E = 3 \times 10^{35}$ erg/s) millisecond pulsar was first detected as a gamma-ray point source in Fermi LAT sky survey observations. Guided by a steep spectrum radio point source in the Fermi error region, we perf… ▽ More We report the discovery of PSR J1555-2908, a 1.79 ms radio and gamma-ray pulsar in a 5.6 hr binary system with a minimum companion mass of 0.052 $M_\odot$. This fast and energetic ($\dot E = 3 \times 10^{35}$ erg/s) millisecond pulsar was first detected as a gamma-ray point source in Fermi LAT sky survey observations. Guided by a steep spectrum radio point source in the Fermi error region, we performed a search at 820 MHz with the Green Bank Telescope that first discovered the pulsations. The initial radio pulse timing observations provided enough information to seed a search for gamma-ray pulsations in the LAT data, from which we derive a timing solution valid for the full Fermi mission. In addition to the radio and gamma-ray pulsation discovery and timing, we searched for X-ray pulsations using NICER but no significant pulsations were detected. We also obtained time-series r-band photometry that indicates strong heating of the companion star by the pulsar wind. Material blown off the heated companion eclipses the 820 MHz radio pulse during inferior conjunction of the companion for ~10% of the orbit, which is twice the angle subtended by its Roche lobe in an edge-on system. △ Less

Submitted 9 February, 2022; originally announced February 2022.

Comments: 15 pages, 6 figures, accepted by ApJ

METAL: The Metal Evolution, Transport, and Abundance in the Large Magellanic Cloud Hubble program. III. Interstellar Depletions, Dust-to-Metal, and Dust-to-Gas Ratios Versus Metallicity

Authors: J. Roman-Duval, E. B. Jenkins, K. Tchernyshyov, C. J. R. Clark, A. De Cia, K. D. Gordon, A. Hamanowicz, V. Lebouteiller, M. Rafelski, K. Sandstrom, J. Werk, P. Yanchulova Merica-Jones

Abstract: The metallicity and gas density dependence of interstellar depletions, the dust-to-gas (D/G), and dust-to-metal (D/M) ratios have important implications for how accurately we can trace the chemical enrichment of the universe; either by using FIR dust emission as a tracer of the ISM; or by using spectroscopy of damped Lyman-alpha systems (DLAs) to measure chemical abundances over a wide range of re… ▽ More The metallicity and gas density dependence of interstellar depletions, the dust-to-gas (D/G), and dust-to-metal (D/M) ratios have important implications for how accurately we can trace the chemical enrichment of the universe; either by using FIR dust emission as a tracer of the ISM; or by using spectroscopy of damped Lyman-alpha systems (DLAs) to measure chemical abundances over a wide range of redshifts. We collect and compare large samples of depletion measurements in the Milky Way (MW), LMC (Z=0.5 Zsun), and SMC (Z=0.2 Zsun). The relation between the depletions of different elements do not strongly vary between the three galaxies, implying that abundance ratios should trace depletions accurately down to 20% solar metallicity. From the depletions, we derive D/G and D/M. The D/G increases with density, consistent with the more efficient accretion of gas-phase metals onto dust grains in the denser ISM. For log N(H) > 21 cm^-2, the depletion of metallicity tracers (S, Zn) exceeds -0.5 dex, even at 20% solar metallicity. The gas fraction of metals increases from the MW to the LMC (factor 3) and SMC (factor 6), compensating the reduction in total heavy element abundances and resulting in those three galaxies having the same neutral gas-phase metallicities. The D/G derived from depletions are a factor of 2 (LMC) and 5 (SMC) higher than the D/G derived from FIR, 21 cm, and CO emission, likely due to the combined uncertainties on the dust FIR opacity and on the depletion of carbon and oxygen. △ Less

Submitted 9 February, 2022; originally announced February 2022.

Comments: 26 pages, 12 figures

The Nearby Evolved Stars Survey II: Constructing a volume-limited sample and first results from the James Clerk Maxwell Telescope

Authors: P. Scicluna, F. Kemper, I. McDonald, S. Srinivasan, A. Trejo, S. H. J. Wallström, J. G. A. Wouterloot, J. Cami, J. Greaves, Jinhua He, D. T. Hoai, Hyosun Kim, O. C. Jones, H. Shinnaga, C. J. R. Clark, T. Dharmawardena, W. Holland, H. Imai, J. Th. van Loon, K. M. Menten, R. Wesson, H. Chawner, S. Feng, S. Goldman, F. C. Liu , et al. (67 additional authors not shown)

Abstract: The Nearby Evolved Stars Survey (NESS) is a volume-complete sample of $\sim$850 Galactic evolved stars within 3\,kpc at (sub-)mm wavelengths, observed in the CO $J = $ (2$-$1) and (3$-$2) rotational lines, and the sub-mm continuum, using the James Clark Maxwell Telescope and Atacama Pathfinder Experiment. NESS consists of five tiers, based on distances and dust-production rate (DPR). We define a n… ▽ More The Nearby Evolved Stars Survey (NESS) is a volume-complete sample of $\sim$850 Galactic evolved stars within 3\,kpc at (sub-)mm wavelengths, observed in the CO $J = $ (2$-$1) and (3$-$2) rotational lines, and the sub-mm continuum, using the James Clark Maxwell Telescope and Atacama Pathfinder Experiment. NESS consists of five tiers, based on distances and dust-production rate (DPR). We define a new metric for estimating the distances to evolved stars and compare its results to \emph{Gaia} EDR3. Replicating other studies, the most-evolved, highly enshrouded objects in the Galactic Plane dominate the dust returned by our sources, and we initially estimate a total DPR of $4.7\times 10^{-5}$ M$_\odot$ yr$^{-1}$ from our sample. Our sub-mm fluxes are systematically higher and spectral indices are typically shallower than dust models typically predict. The 450/850 $μ$m spectral indices are consistent with the blackbody Rayleigh--Jeans regime, suggesting a large fraction of evolved stars have unexpectedly large envelopes of cold dust. △ Less

Submitted 24 October, 2021; originally announced October 2021.

Comments: 20 pages, 13 figures, 5 tables, accepted for publications in MNRAS

The HASHTAG project: The First Submillimeter Images of the Andromeda Galaxy from the Ground

Authors: Matthew W. L. Smith, Stephen A. Eales, Thomas G. Williams, Bumhyun Lee, Zongnan Li, Pauline Barmby, Martin Bureau, Scott Chapman, Brian S. Cho, Aeree Chung, Eun Jung Chung, Hui-Hsuan Chung, Christopher J. R. Clark, David L. Clements, Timothy A. Davis, Ilse De Looze, David J. Eden, Gayathri Athikkat-Eknath, George P. Ford, Yu Gao, Walter Gear, Haley L. Gomez, Richard de Grijs, Jinhua He, Luis C. Ho , et al. (24 additional authors not shown)

Abstract: Observing nearby galaxies with submillimeter telescopes on the ground has two major challenges. First, the brightness is significantly reduced at long submillimeter wavelengths compared to the brightness at the peak of the dust emission. Second, it is necessary to use a high-pass spatial filter to remove atmospheric noise on large angular scales, which has the unwelcome by-product of also removing… ▽ More Observing nearby galaxies with submillimeter telescopes on the ground has two major challenges. First, the brightness is significantly reduced at long submillimeter wavelengths compared to the brightness at the peak of the dust emission. Second, it is necessary to use a high-pass spatial filter to remove atmospheric noise on large angular scales, which has the unwelcome by-product of also removing the galaxy's large-scale structure. We have developed a technique for producing high-resolution submillimeter images of galaxies of large angular size by using the telescope on the ground to determine the small-scale structure (the large Fourier components) and a space telescope (Herschel or Planck) to determine the large-scale structure (the small Fourier components). Using this technique, we are carrying out the HARP and SCUBA-2 High Resolution Terahertz Andromeda Galaxy Survey (HASHTAG), an international Large Program on the James Clerk Maxwell Telescope, with one aim being to produce the first high-fidelity high-resolution submillimeter images of Andromeda. In this paper, we describe the survey, the method we have developed for combining the space-based and ground-based data, and present the first HASHTAG images of Andromeda at 450 and 850um. We also have created a method to predict the CO(J=3-2) line flux across M31, which contaminates the 850um band. We find that while normally the contamination is below our sensitivity limit, the contamination can be significant (up to 28%) in a few of the brightest regions of the 10 kpc ring. We therefore also provide images with the predicted line emission removed. △ Less

Submitted 30 September, 2021; originally announced October 2021.

Comments: 26 pages, 19 figures. Submitted to ApJS June 2021, Accepted September 2021

FAST discovery of an extremely radio-faint millisecond pulsar from the Fermi-LAT unassociated source 3FGL J0318.1+0252

Authors: Pei Wang, Di Li, Colin J. Clark, Pablo Saz Parkinson, Xian Hou, Weiwei Zhu, Lei Qian, Youling Yue, Zhichen Pan, Zhijie Liu, Xuhong Yu, Xiaoyao Xie, Qijun Zhi, Hui Zhang, Jumei Yao, Jun Yan, Chengmin Zhang, Paul S. Ray, Matthew Kerr, David A. Smith, Peter F. Michelson, Elizabeth C. Ferrara, David J. Thompson, Zhiqiang Shen, Na Wang , et al. (1 additional authors not shown)

Abstract: High sensitivity radio searches of unassociated $γ$-ray sources have proven to be an effective way of finding new pulsars. Using the Five-hundred-meter Aperture Spherical radio Telescope (FAST) during its commissioning phase, we have carried out a number of targeted deep searches of \textit{Fermi} Large Area Telescope (LAT) $γ$-ray sources. On Feb. 27$^{th}$, 2018 we discovered an isolated millise… ▽ More High sensitivity radio searches of unassociated $γ$-ray sources have proven to be an effective way of finding new pulsars. Using the Five-hundred-meter Aperture Spherical radio Telescope (FAST) during its commissioning phase, we have carried out a number of targeted deep searches of \textit{Fermi} Large Area Telescope (LAT) $γ$-ray sources. On Feb. 27$^{th}$, 2018 we discovered an isolated millisecond pulsar (MSP), PSR J0318+0253, coincident with the unassociated $γ$-ray source 3FGL J0318.1+0252. PSR J0318+0253 has a spin period of $5.19$ milliseconds, a dispersion measure (DM) of $26$ pc cm$^{-3}$ corresponding to a DM distance of about $1.3$ kpc, and a period-averaged flux density of $\sim$11 $\pm$ 2 $μ$Jy at L-band (1.05-1.45 GHz). Among all high energy MSPs, PSR J0318+0253 is the faintest ever detected in radio bands, by a factor of at least $\sim$4 in terms of L-band fluxes. With the aid of the radio ephemeris, an analysis of 9.6 years of \textit{Fermi}-LAT data revealed that PSR J0318+0253 also displays strong $γ$-ray pulsations. Follow-up observations carried out by both Arecibo and FAST suggest a likely spectral turn-over around 350 MHz. This is the first result from the collaboration between FAST and the \textit{Fermi}-LAT teams as well as the first confirmed new MSP discovery by FAST, raising hopes for the detection of many more MSPs. Such discoveries will make a significant contribution to our understanding of the neutron star zoo while potentially contributing to the future detection of gravitational waves, via pulsar timing array (PTA) experiments. △ Less

Submitted 3 September, 2021; v1 submitted 2 September, 2021; originally announced September 2021.

Comments: 7 pages, 2 figures, 1 table

Journal ref: Sci. China-Phys. Mech. Astron. 64: 129562 (2021)

The Quest for the Missing Dust: I -- Restoring Large Scale Emission in Herschel Maps of Local Group Galaxies

Authors: Christopher J. R. Clark, Julia C. Roman-Duval, Karl D. Gordon, Caroline Bot, Matthew W. L. Smith

Abstract: Because the galaxies of the Local Group have such large angular sizes, much of their diffuse, large-angular-scale emission is filtered out by the Herschel data reduction process. In this work, we restore this previously missed dust in Herschel observations of the Large Magellanic Cloud, Small Magellanic Cloud, M31, and M33. We do this by combining Herschel data (including new reductions for the Ma… ▽ More Because the galaxies of the Local Group have such large angular sizes, much of their diffuse, large-angular-scale emission is filtered out by the Herschel data reduction process. In this work, we restore this previously missed dust in Herschel observations of the Large Magellanic Cloud, Small Magellanic Cloud, M31, and M33. We do this by combining Herschel data (including new reductions for the Magellanic Clouds), in Fourier space, with lower-resolution data from all-sky surveys (Planck, IRAS, and COBE) that did not miss the extended emission. With these new maps, we find that a significant amount of emission was missing from uncorrected Herschel data of these galaxies; over 20% in some bands. Our new photometry also resolves the disagreement between fluxes reported from older HERITAGE Magellanic Cloud Herschel reductions, and fluxes reported from other telescopes. More emission is restored in shorter wavelength bands, especially in the galaxies' peripheries, making these regions 20-40% bluer than before. We also find that the Herschel-PACS instrument response conflicts with the all-sky data, over the 20-90' angular scales to which they are both sensitive, by up to 31%. By binning our new data based on hydrogen column density, we are able to detect emission from dust at low ISM densities (at $Σ_{\rm H} < 1\,{\rm M_{\odot} pc^{-2}}$ in some cases), and are able to detect emission at much lower densities (a factor of 2.2 lower on average, and more than a factor of 7 lower in several cases) than was possible with uncorrected data. △ Less

Submitted 9 September, 2021; v1 submitted 29 July, 2021; originally announced July 2021.

Comments: Accepted for publication in ApJ

Journal ref: ApJ 921 35 (2021)

Optical photometry of two transitional millisecond pulsars in the radio pulsar state

Authors: James G. Stringer, Rene P. Breton, Colin J. Clark, Guillaume Voisin, Mark R. Kennedy, Daniel Mata Sànchez, Tariq Shahbaz, Vik S. Dhillon, Marten van Kerkwijk, Tom R. Marsh

Abstract: We present ULTRACAM multiband optical photometry of two transitional millisecond pulsars, PSR J1023+0038 and PSR J1227$-$4853, taken while both were in their radio pulsar states. The light curves show significant asymmetry about the flux maxima in all observed bands, suggesting an asymmetric source of heating in the system. We model the light curves using the Icarus binary code, using models with… ▽ More We present ULTRACAM multiband optical photometry of two transitional millisecond pulsars, PSR J1023+0038 and PSR J1227$-$4853, taken while both were in their radio pulsar states. The light curves show significant asymmetry about the flux maxima in all observed bands, suggesting an asymmetric source of heating in the system. We model the light curves using the Icarus binary code, using models with an additional "hot spot" heating contribution and an anisotropic heat redistribution model to treat the asymmetry. Our modelling reveals companion stars with under-filled Roche lobes in both PSRs J1023+0038 and J1227$-$4853, with Roche lobe filling factors in the range $f \sim 0.82-0.92$. While the volume-averaged filling factors are closer to unity, significant under-filling is unexpected from tMSPs as they must rapidly over-fill their Roche lobes to start transferring mass, which occurs on timescale of weeks or months. We discuss the motivation and validity of our extensions to the models and the implications of the under-filled Roche lobe, and suggest future work to further investigate the role of the filling factor in the tMSP cycle. △ Less

Submitted 20 July, 2021; originally announced July 2021.

Comments: 19 pages, 13 figures

Probing the spectral shape of dust emission with the DustPedia galaxy sample

Authors: Angelos Nersesian, Wouter Dobbels, Manolis E. Xilouris, Maarten Baes, Simone Bianchi, Viviana Casasola, Christopher J. R. Clark, Ilse De Looze, Frédéric Galliano, Suzanne C. Madden, Aleksandr V. Mosenkov, Evangelos-D. Paspaliaris, Ana Trčka

Abstract: The objective of this paper is to understand the variance of the far-infrared (FIR) spectral energy distribution (SED) of the DustPedia galaxies, and its link with the stellar and dust properties. An interesting aspect of the dust emission is the inferred FIR colours which could inform us about the dust content of galaxies, and how it varies with the physical conditions within galaxies. However, t… ▽ More The objective of this paper is to understand the variance of the far-infrared (FIR) spectral energy distribution (SED) of the DustPedia galaxies, and its link with the stellar and dust properties. An interesting aspect of the dust emission is the inferred FIR colours which could inform us about the dust content of galaxies, and how it varies with the physical conditions within galaxies. However, the inherent complexity of dust grains as well as the variety of physical properties depending on dust, hinder our ability to utilise their maximum potential. We use principal component analysis (PCA) to explore new hidden correlations with many relevant physical properties such as the dust luminosity, dust temperature, dust mass, bolometric luminosity, star-formation rate (SFR), stellar mass, specific SFR, dust-to-stellar mass ratio, the fraction of absorbed stellar luminosity by dust (f_abs), and metallicity. We find that 95% of the variance in our sample can be described by two principal components (PCs). The first component controls the wavelength of the peak of the SED, while the second characterises the width. The physical quantities that correlate better with the coefficients of the first two PCs, and thus control the shape of the FIR SED are: the dust temperature, the dust luminosity, the SFR, and f_abs. Finally, we find a weak tendency for low-metallicity galaxies to have warmer and broader SEDs, while on the other hand high-metallicity galaxies have FIR SEDs that are colder and narrower. △ Less

Submitted 9 July, 2021; originally announced July 2021.

Constraining the origin of the puzzling source HESS J1640-465 and the PeVatron candidate HESS J1641-463 using Fermi-LAT observations

Authors: A. Mares, M. Lemoine-Goumard, F. Acero, C. J. Clark, J. Devin, S. Gabici, J. D. Gelfand, D. A. Green, M. -H. Grondin

Abstract: There are only few very-high-energy sources in our Galaxy which might accelerate particles up to the knee of the cosmic-ray spectrum. To understand the mechanisms of particle acceleration in these PeVatron candidates, \textit{Fermi}-LAT and H.E.S.S. observations are essential to characterize their $γ$-ray emission. HESS J1640$-$465 and the PeVatron candidate HESS J1641$-$463 are two neighboring (\… ▽ More There are only few very-high-energy sources in our Galaxy which might accelerate particles up to the knee of the cosmic-ray spectrum. To understand the mechanisms of particle acceleration in these PeVatron candidates, \textit{Fermi}-LAT and H.E.S.S. observations are essential to characterize their $γ$-ray emission. HESS J1640$-$465 and the PeVatron candidate HESS J1641$-$463 are two neighboring (\ang[astroang]{0.25}) $γ$-ray sources, spatially coincident with the radio supernova remnants (SNRs) G338.3$-$0.0 and G338.5+0.1. Detected both by H.E.S.S. and \textit{Fermi}-LAT, we present here a morphological and spectral analysis of these two sources using 8 years of \textit{Fermi}-LAT data between 200 \si{\mega\electronvolt} and 1 \si{\tera\electronvolt} with multi-wavelength observations to assess their nature. The morphology of HESS J1640$-$465 is described by a 2D Gaussian ($σ=$ \ang[astroang]{0.053} $\pm$ \ang[astroang]{0.011}$_{stat}$ $ \pm$ \ang[astroang]{0.03}$_{syst}$) and its spectrum is modeled by a power-law with a spectral index $Γ= 1.8\pm0.1_{\rm stat}\pm0.2_{\rm syst}$. HESS J1641$-$463 is detected as a point-like source and its GeV emission is described by a logarithmic-parabola spectrum with $α= 2.7 \pm 0.1_ {\rm stat} \pm 0.2_ {\rm syst} $ and significant curvature of $β= 0.11 \pm 0.03_ {\rm stat} \pm 0.05_ {\rm syst} $. Radio and X-ray flux upper limits were derived. We investigated scenarios to explain their emission, namely the emission from accelerated particles within the SNRs spatially coincident with each source, molecular clouds illuminated by cosmic rays from the close-by SNRs, and a pulsar/PWN origin. Our new \emph{Fermi}-LAT results and the radio and flux X-ray upper limits pose severe constraints on some of these models. △ Less

Submitted 5 April, 2021; originally announced April 2021.

Comments: 18 pages 7 figures To be published in ApJ

A search for radio pulsars in five nearby supernova remnants

Authors: S. Sett, R. P. Breton, C. J. Clark, M. H. Kerkwijk, D. L. Kaplan

Abstract: Most neutron stars are expected to be born in supernovae, but only about half of supernova remnants (SNRs) are associated with a compact object. In many cases, a supernova progenitor may have resulted in a black hole. However, there are several possible reasons why true pulsar-SNR associations may have been missed in previous surveys: The pulsar's radio beam may not be oriented towards us; the pul… ▽ More Most neutron stars are expected to be born in supernovae, but only about half of supernova remnants (SNRs) are associated with a compact object. In many cases, a supernova progenitor may have resulted in a black hole. However, there are several possible reasons why true pulsar-SNR associations may have been missed in previous surveys: The pulsar's radio beam may not be oriented towards us; the pulsar may be too faint to be detectable; or there may be an offset in the pulsar position caused by a kick. Our goal is to find new pulsars in SNRs and explore their possible association with the remnant. The search and selection of the remnants presented in this paper was inspired by the non-detection of any X-ray bright compact objects in these remnants when previously studied. Five SNRs were searched for radio pulsars with the Green Bank Telescope at 820 MHz with multiple pointings to cover the full spatial extent of the remnants. A periodicity search plus an acceleration search up to 500 m/s^2 and a single pulse search were performed for each pointing in order to detect potential isolated binary pulsars and single pulses, respectively. No new pulsars were detected in the survey. However, we were able to re-detect a known pulsar, PSR J2047+5029, near SNR G89.0+4.7. We were unable to detect the radio-quiet gamma-ray pulsar PSR J2021+4026, but we do find a flux density limit of 0.08 mJy. Our flux density limits make our survey two to 16 times more sensitive than previous surveys, while also covering the whole spatial extent of the same remnants. We discuss potential explanations for the non-detection of a pulsar in the studied SNRs and conclude that sensitivity is still the most likely factor responsible for the lack of pulsars in some remnants. △ Less

Submitted 29 January, 2021; originally announced January 2021.

Comments: Accepted for publication in A&A on 14/01/2021

Report number: aa36108-19

Journal ref: A&A 647, A183 (2021)

METAL: The Metal Evolution, Transport, and Abundance in the Large Magellanic Cloud Hubble program. II. Variations of interstellar depletions and dust-to-gas ratio within the LMC

Authors: Julia Roman-Duval, Edward B. Jenkins, Kirill Tchernyshyov, Benjamin Williams, Christopher J. R. Clark, Karl D. Gordon, Margaret Meixner, Lea Hagen, Joshua Peek, Karin Sandstrom, Jessica Werk, Petia Yanchulova Merica-Jones

Abstract: A key component of the baryon cycle in galaxies is the depletion of metals from the gas to the dust phase in the neutral ISM. The METAL (Metal Evolution, Transport and Abundance in the Large Magellanic Cloud) program on the Hubble Space Telescope acquired UV spectra toward 32 sightlines in the half-solar metallicity LMC, from which we derive interstellar depletions (gas-phase fractions) of Mg, Si,… ▽ More A key component of the baryon cycle in galaxies is the depletion of metals from the gas to the dust phase in the neutral ISM. The METAL (Metal Evolution, Transport and Abundance in the Large Magellanic Cloud) program on the Hubble Space Telescope acquired UV spectra toward 32 sightlines in the half-solar metallicity LMC, from which we derive interstellar depletions (gas-phase fractions) of Mg, Si, Fe, Ni, S, Zn, Cr, and Cu. The depletions of different elements are tightly correlated, indicating a common origin. Hydrogen column density is the main driver for depletion variations. Correlations are weaker with volume density, probed by CI fine structure lines, and distance to the LMC center. The latter correlation results from an East-West variation of the gas-phase metallicity. Gas in the East, compressed side of the LMC encompassing 30 Doradus and the Southeast HI over-density is enriched by up to +0.3dex, while gas in the West side is metal-deficient by up to -0.5dex. Within the parameter space probed by METAL, no correlation with molecular fraction or radiation field intensity are found. We confirm the factor 3-4 increase in dust-to-metal and dust-to-gas ratios between the diffuse (logN(H)~20 cm-2) and molecular (logN(H)~22 cm-2) ISM observed from far-infrared, 21 cm, and CO observations. The variations of dust-to-metal and dust-to-gas ratios with column density have important implications for the sub-grid physics of chemical evolution, gas and dust mass estimates throughout cosmic times, and for the chemical enrichment of the Universe measured via spectroscopy of damped Lyman-alpha systems. △ Less

Submitted 22 January, 2021; originally announced January 2021.

High-resolution, 3D radiative transfer modelling V. A detailed model of the M51 interacting pair

Authors: Angelos Nersesian, Sebastien Viaene, Ilse De Looze, Maarten Baes, Emmanuel M. Xilouris, Matthew W. L. Smith, Simone Bianchi, Viviana Casasola, Letizia P. Cassara, Christopher J. R. Clark, Wouter Dobbels, Jacopo Fritz, Frederic Galliano, Suzanne C. Madden, Aleksandr V. Mosenkov, Ana Trcka

Abstract: Investigating the dust heating mechanisms in galaxies provides a deeper understanding of how the internal energy balance drives their evolution. Over the last decade, radiative transfer simulations based on the Monte Carlo method have underlined the role of the various stellar populations heating the diffuse dust. Beyond the expected heating through ongoing star formation, both older stellar popul… ▽ More Investigating the dust heating mechanisms in galaxies provides a deeper understanding of how the internal energy balance drives their evolution. Over the last decade, radiative transfer simulations based on the Monte Carlo method have underlined the role of the various stellar populations heating the diffuse dust. Beyond the expected heating through ongoing star formation, both older stellar population (> 8Gyr) and even AGN can contribute energy to the infrared emission of diffuse dust. Here, we examine how the radiation of an external heating source, like the less massive galaxy NGC5195, in the M51 interacting system, could affect the heating of the diffuse dust of its parent galaxy, NGC5194, and vice versa. To quantify the exchange of energy between the two galaxies we use SKIRT, a state-of-the-art Monte Carlo radiative transfer code. In the interest of modelling, the assumed centre-to-centre distance separation between the two galaxies is 10kpc. Our model reproduces the global SED of the system, and it closely matches the observed images. In total, 40.7% of the intrinsic stellar radiation of the combined system is absorbed by dust. Furthermore, we quantify the contribution of the various dust heating sources in the system, and find that the young stellar population of NGC5194 is the predominant dust-heating agent, with a global heating fraction of 71.2%. Another 23% is provided by the older stellar population of the same galaxy, while the remaining 5.8% has its origin in NGC5195. Locally, we find that the regions of NGC5194 closer to NGC5195 are significantly affected by the radiation field of the latter, with the absorbed energy fraction rising up to 38%. The contribution of NGC5195 remains under the percentage level in the outskirts of the disc of NGC5194. This is the first time that the heating of the diffuse dust by a companion galaxy is quantified in a nearby interacting system. △ Less

Submitted 18 September, 2020; v1 submitted 15 September, 2020; originally announced September 2020.

Journal ref: A&A 643, A90 (2020)

Discovery of a Gamma-ray Black Widow Pulsar by GPU-accelerated Einstein@Home

Authors: L. Nieder, C. J. Clark, D. Kandel, R. W. Romani, C. G. Bassa, B. Allen, A. Ashok, I. Cognard, H. Fehrmann, P. Freire, R. Karuppusamy, M. Kramer, D. Li, B. Machenschalk, Z. Pan, M. A. Papa, S. M. Ransom, P. S. Ray, J. Roy, P. Wang, J. Wu, C. Aulbert, E. D. Barr, B. Beheshtipour, O. Behnke , et al. (17 additional authors not shown)

Abstract: We report the discovery of 1.97 ms period gamma-ray pulsations from the 75 minute orbital-period binary pulsar now named PSR J1653-0158. The associated Fermi Large Area Telescope gamma-ray source 4FGL J1653.6-0158 has long been expected to harbor a binary millisecond pulsar. Despite the pulsar-like gamma-ray spectrum and candidate optical/X-ray associations -- whose periodic brightness modulations… ▽ More We report the discovery of 1.97 ms period gamma-ray pulsations from the 75 minute orbital-period binary pulsar now named PSR J1653-0158. The associated Fermi Large Area Telescope gamma-ray source 4FGL J1653.6-0158 has long been expected to harbor a binary millisecond pulsar. Despite the pulsar-like gamma-ray spectrum and candidate optical/X-ray associations -- whose periodic brightness modulations suggested an orbit -- no radio pulsations had been found in many searches. The pulsar was discovered by directly searching the gamma-ray data using the GPU-accelerated Einstein@Home distributed volunteer computing system. The multi-dimensional parameter space was bounded by positional and orbital constraints obtained from the optical counterpart. More sensitive analyses of archival and new radio data using knowledge of the pulsar timing solution yield very stringent upper limits on radio emission. Any radio emission is thus either exceptionally weak, or eclipsed for a large fraction of the time. The pulsar has one of the three lowest inferred surface magnetic-field strengths of any known pulsar with $B_{\rm surf} \approx 4 \times 10^{7}\,$G. The resulting mass function, combined with models of the companion star's optical light curve and spectra, suggests a pulsar mass $\gtrsim 2\,M_{\odot}$. The companion is light-weight with mass $\sim 0.01\,M_{\odot}$, and the orbital period is the shortest known for any rotation-powered binary pulsar. This discovery demonstrates the Fermi Large Area Telescope's potential to discover extreme pulsars that would otherwise remain undetected. △ Less

Submitted 22 October, 2020; v1 submitted 3 September, 2020; originally announced September 2020.

Comments: 12 pages, 3 figures, published in ApJL

Radio pulsations from the $γ$-ray millisecond pulsar PSR J2039-5617

Authors: A. Corongiu, R. P. Mignani, A. S. Seyffert, C. J. Clark, C. Venter, L. Nieder, A. Possenti, M. Burgay, A. Belfiore, A. De Luca, A. Ridolfi, Z. Wadiasingh

Abstract: The predicted nature of the candidate redback pulsar 3FGL\,J2039.6$-$5618 was recently confirmed by the discovery of $γ$-ray millisecond pulsations (Clark et al. 2020, hereafter Paper\,I), which identify this $γ$-ray source as \msp. We observed this object with the Parkes radio telescope in 2016 and 2019. We detect radio pulsations at 1.4\,GHz and 3.1\,GHz, at the 2.6ms period discovered in $γ$-ra… ▽ More The predicted nature of the candidate redback pulsar 3FGL\,J2039.6$-$5618 was recently confirmed by the discovery of $γ$-ray millisecond pulsations (Clark et al. 2020, hereafter Paper\,I), which identify this $γ$-ray source as \msp. We observed this object with the Parkes radio telescope in 2016 and 2019. We detect radio pulsations at 1.4\,GHz and 3.1\,GHz, at the 2.6ms period discovered in $γ$-rays, and also at 0.7\,GHz in one 2015 archival observation. In all bands, the radio pulse profile is characterised by a single relatively broad peak which leads the main $γ$-ray peak. At 1.4\,GHz we found clear evidence of eclipses of the radio signal for about half of the orbit, a characteristic phenomenon in redback systems, which we associate with the presence of intra-binary gas. From the dispersion measure of $24.57\pm0.03$\,pc\,cm$^{-3}$ we derive a pulsar distance of $0.9\pm 0.2$\,kpc or $1.7\pm0.7$\,kpc, depending on the assumed Galactic electron density model. The modelling of the radio and $γ$-ray light curves leads to an independent determination of the orbital inclination, and to a determination of the pulsar mass, qualitatively consistent to the results in Paper\,I. △ Less

Submitted 3 November, 2020; v1 submitted 29 July, 2020; originally announced July 2020.

Comments: 18 pages, accepted for publication on MNRAS

Einstein@Home Discovery of the Gamma-ray Millisecond Pulsar PSR J2039-5617 Confirms Its Predicted Redback Nature

Authors: C. J. Clark, L. Nieder, G. Voisin, B. Allen, C. Aulbert, O. Behnke, R. P. Breton, C. Choquet, A. Corongiu, V. S. Dhillon, H. B. Eggenstein, H. Fehrmann, L. Guillemot, A. K. Harding, M. R. Kennedy, B. Machenschalk, T. R. Marsh, D. Mata Sánchez, R. P. Mignani, J. Stringer, Z. Wadiasingh, J. Wu

Abstract: The Fermi Large Area Telescope gamma-ray source 3FGL J2039.6$-$5618 contains a periodic optical and X-ray source that was predicted to be a "redback" millisecond pulsar (MSP) binary system. However, the conclusive identification required the detection of pulsations from the putative MSP. To better constrain the orbital parameters for a directed search for gamma-ray pulsations, we obtained new opti… ▽ More The Fermi Large Area Telescope gamma-ray source 3FGL J2039.6$-$5618 contains a periodic optical and X-ray source that was predicted to be a "redback" millisecond pulsar (MSP) binary system. However, the conclusive identification required the detection of pulsations from the putative MSP. To better constrain the orbital parameters for a directed search for gamma-ray pulsations, we obtained new optical light curves in 2017 and 2018, which revealed long-term variability from the companion star. The resulting orbital parameter constraints were used to perform a targeted gamma-ray pulsation search using the Einstein@Home distributed volunteer computing system. This search discovered pulsations with a period of 2.65 ms, confirming the source as a binary MSP now known as PSR J2039$-$5617. Optical light curve modelling is complicated, and likely biased, by asymmetric heating on the companion star and long-term variability, but we find an inclination $i > 60°$, for a low pulsar mass between $1.1 M_{\odot} < M_{\rm psr} < 1.6 M_{\odot}$ and a companion mass of 0.15--0.22 $M_{\odot}$, confirming the redback classification. Timing the gamma-ray pulsations also revealed significant variability in the orbital period, which we find to be consistent with quadrupole moment variations in the companion star, suggestive of convective activity. We also find that the pulsed flux is modulated at the orbital period, potentially due to inverse Compton scattering between high-energy leptons in the pulsar wind and the companion star's optical photon field. △ Less

Submitted 2 February, 2021; v1 submitted 29 July, 2020; originally announced July 2020.

Comments: 21 pages, 11 figures, published in MNRAS

Journal ref: MNRAS 502 915 (2021)

Nonparametric galaxy morphology from UV to submm wavelengths

Authors: Maarten Baes, Angelos Nersesian, Viviana Casasola, Simone Bianchi, Letizia P. Cassarà, Christopher J. R. Clark, Ilse De Looze, Wouter Dobbels, Jacopo Fritz, Maud Galametz, Frédéric Galliano, Suzanne C. Madden, Aleksandr V. Mosenkov, Sébastien Viaene, Ana Trčka, Emmanuel M. Xilouris

Abstract: We present the first nonparametric morphological analysis of a set of spiral galaxies from UV to submm wavelengths. Our study is based on high-quality multi-wavelength imaging for nine well-resolved spiral galaxies from the DustPedia database, combined with nonparametric morphology indicators calculated in a consistent way using the {\tt{StatMorph}} package. We measure the half-light radius, the c… ▽ More We present the first nonparametric morphological analysis of a set of spiral galaxies from UV to submm wavelengths. Our study is based on high-quality multi-wavelength imaging for nine well-resolved spiral galaxies from the DustPedia database, combined with nonparametric morphology indicators calculated in a consistent way using the {\tt{StatMorph}} package. We measure the half-light radius, the concentration index, the asymmetry index, the smoothness index, the Gini coefficient and the $M_{20}$ indicator in various wavebands from UV to submm wavelengths, as well as in stellar mass, dust mass and star formation rate maps. We find that the interstellar dust in galaxies is distributed in a more extended, less centrally concentrated, more asymmetric, and more clumpy way than the stars. This is particularly evident when comparing morphological indicators based on the stellar mass and dust mass maps. This should serve as a warning sign against treating the dust in galaxies as a simple smooth component. We argue that the nonparametric galaxy morphology of galaxies from UV to submm wavelengths is an interesting test for cosmological hydrodynamics simulations. △ Less

Submitted 17 July, 2020; v1 submitted 4 July, 2020; originally announced July 2020.

Comments: 9 pages, 5 figures, accepted for publication in A&A

Journal ref: A&A 641, A119 (2020)

A model for redistributing heat over the surface of irradiated spider companions

Authors: G. Voisin, M. R. Kennedy, R. P. Breton, C. J. Clark, D. Mata-Sánchez

Abstract: Spider pulsars are binary systems containing an energetic millisecond pulsar that intensely irradiates a closely orbiting low-mass companion. Modelling their companion's optical light curves is essential to the study of the orbital properties of the binary, including the determination of the pulsar mass, characterising the pulsar wind and the star itself. We aim to generalise the traditional direc… ▽ More Spider pulsars are binary systems containing an energetic millisecond pulsar that intensely irradiates a closely orbiting low-mass companion. Modelling their companion's optical light curves is essential to the study of the orbital properties of the binary, including the determination of the pulsar mass, characterising the pulsar wind and the star itself. We aim to generalise the traditional direct heating model of irradiation, whereby energy deposited by the pulsar wind into the stellar envelope is locally re-emitted, by introducing heat redistribution via diffusion and convection within the outer stellar envelope. We approximate the irradiated stellar envelope as a two-dimensional shell. This allows us to propose an effective equation of energy conservation that can be solved at a reduced computational cost. We then implement this model in the \texttt{Icarus} software and use evidence sampling to determine the most likely convection and diffusion laws for the light curve of the redback companion of PSR J2215+5135. Redistribution effects concentrate near the terminator line of pulsar irradiation, and can create apparent hot and cold spots. Among the models tested for PSR J2215+5135, we find that all models with heat redistribution are more likely than symmetric direct heating. The best-fitting redistribution model involves diffusion together with a uniformly rotating envelope. However, we caution that all models still present serious systematic effects, and that prior knowledge from pulsar timing, spectroscopy and distance are key to determine with certainty the most accurate redistribution law. We propose an extension of the direct heating framework that allows for exploring a variety of heat redistribution effects. Future work is necessary to determine the relevant laws from first principles and empirically using complementary observations. △ Less

Submitted 16 November, 2020; v1 submitted 11 June, 2020; originally announced June 2020.

Comments: 13 pages, 5 figures, published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 499, Issue 2, December 2020, Pages 1758-1768

JINGLE -- IV. Dust, HI gas and metal scaling laws in the local Universe

Authors: I. De Looze, I. Lamperti, A. Saintonge, M. Relano, M. W. L. Smith, C. J. R. Clark, C. D. Wilson, M. Decleir, A. P. Jones, R. C. Kennicutt, G. Accurso, E. Brinks, M. Bureau, P. Cigan, D. L. Clements, P. De Vis, L Fanciullo, Y. Gao, W. K. Gear, L. C. Ho, H. S. Hwang, M. J. Michalowski, J. C. Lee, C. Li, L. Lin , et al. (7 additional authors not shown)

Abstract: Scaling laws of dust, HI gas and metal mass with stellar mass, specific star formation rate and metallicity are crucial to our understanding of the buildup of galaxies through their enrichment with metals and dust. In this work, we analyse how the dust and metal content varies with specific gas mass ($M_{\text{HI}}$/$M_{\star}$) across a diverse sample of 423 nearby galaxies. The observed trends a… ▽ More Scaling laws of dust, HI gas and metal mass with stellar mass, specific star formation rate and metallicity are crucial to our understanding of the buildup of galaxies through their enrichment with metals and dust. In this work, we analyse how the dust and metal content varies with specific gas mass ($M_{\text{HI}}$/$M_{\star}$) across a diverse sample of 423 nearby galaxies. The observed trends are interpreted with a set of Dust and Element evolUtion modelS (DEUS) - incluidng stellar dust production, grain growth, and dust destruction - within a Bayesian framework to enable a rigorous search of the multi-dimensional parameter space. We find that these scaling laws for galaxies with $-1.0\lesssim \log M_{\text{HI}}$/$M_{\star}\lesssim0$ can be reproduced using closed-box models with high fractions (37-89$\%$) of supernova dust surviving a reverse shock, relatively low grain growth efficiencies ($ε$=30-40), and long dus lifetimes (1-2\,Gyr). The models have present-day dust masses with similar contributions from stellar sources (50-80\,$\%$) and grain growth (20-50\,$\%$). Over the entire lifetime of these galaxies, the contribution from stardust ($>$90\,$\%$) outweighs the fraction of dust grown in the interstellar medium ($<$10$\%$). Our results provide an alternative for the chemical evolution models that require extremely low supernova dust production efficiencies and short grain growth timescales to reproduce local scaling laws, and could help solving the conundrum on whether or not grains can grow efficiently in the interstellar medium. △ Less

Submitted 2 June, 2020; originally announced June 2020.

Comments: 56 pages, 30 figures, accepted for publication in MNRAS

Exploiting Orbital Constraints from Optical Data to Detect Binary Gamma-ray Pulsars

Authors: L. Nieder, B. Allen, C. J. Clark, H. J. Pletsch

Abstract: It is difficult to discover pulsars via their gamma-ray emission because current instruments typically detect fewer than one photon per million rotations. This creates a significant computing challenge for isolated pulsars, where the typical parameter search space spans wide ranges in four dimensions. It is even more demanding when the pulsar is in a binary system, where the orbital motion introdu… ▽ More It is difficult to discover pulsars via their gamma-ray emission because current instruments typically detect fewer than one photon per million rotations. This creates a significant computing challenge for isolated pulsars, where the typical parameter search space spans wide ranges in four dimensions. It is even more demanding when the pulsar is in a binary system, where the orbital motion introduces several additional unknown parameters. Building on earlier work by Pletsch & Clark (arXiv:1408.6962), we present optimal methods for such searches. These can also incorporate external constraints on the parameter space to be searched, for example, from optical observations of a presumed binary companion. The solution has two parts. The first is the construction of optimal search grids in parameter space via a parameter-space metric, for initial semicoherent searches and subsequent fully coherent follow-ups. The second is a method to demodulate and detect the periodic pulsations. These methods have different sensitivity properties than traditional radio searches for binary pulsars and might unveil new populations of pulsars. △ Less

Submitted 22 October, 2020; v1 submitted 24 April, 2020; originally announced April 2020.

Comments: 33 pages, 7 figures, published in ApJ

High-resolution, 3D radiative transfer modelling III. The DustPedia barred galaxies

Authors: Angelos Nersesian, Sam Verstocken, Sebastien Viaene, Maarten Baes, Emmanuel M. Xilouris, Simone Bianchi, Viviana Casasola, Christopher J. R. Clark, Jonathan I. Davies, Ilse De Looze, Pieter De Vis, Wouter Dobbels, Jacopo Fritz, Maud Galametz, Frederic Galliano, Anthony P. Jones, Suzanne C. Madden, Aleksandr V. Mosenkov, Ana Trcka, Nathalie Ysard

Abstract: Context: Dust in late-type galaxies in the local Universe is responsible for absorbing approximately one third of the energy emitted by stars. It is often assumed that dust heating is mainly attributable to the absorption of UV and optical photons emitted by the youngest (<= 100 Myr) stars. Consequently, thermal re-emission by dust at FIR wavelengths is often linked to the star-formation activity… ▽ More Context: Dust in late-type galaxies in the local Universe is responsible for absorbing approximately one third of the energy emitted by stars. It is often assumed that dust heating is mainly attributable to the absorption of UV and optical photons emitted by the youngest (<= 100 Myr) stars. Consequently, thermal re-emission by dust at FIR wavelengths is often linked to the star-formation activity of a galaxy. However, several studies argue that the contribution to dust heating by much older stars might be more significant. Advances in radiation transfer (RT) simulations finally allow us to actually quantify the heating mechanisms of diffuse dust by the stellar radiation field. Aims: As one of the main goals in the DustPedia project, we have constructed detailed 3D stellar and dust RT models for nearby galaxies. We analyse the contribution of the different stellar populations to the dust heating in four face-on barred galaxies: NGC1365, M83, M95, and M100. We aim to quantify the fraction directly related to young stars, both globally and on local scales, and to assess the influence of the bar on the heating fraction. Results: We derive global attenuation laws for each galaxy and confirm that galaxies of high sSFR have shallower attenuation curves and weaker UV bumps. On average, 36.5% of the bolometric luminosity is absorbed by dust. We report a clear effect of the bar structure on the radial profiles of the dust-heating fraction by the young stars, and the dust temperature. We find that the young stars are the main contributors to the dust heating, donating, on average ~59% of their luminosity to this purpose throughout the galaxy. This dust-heating fraction drops to ~53% in the bar region and ~38% in the bulge region where the old stars are the dominant contributors to the dust heating. We also find a strong link between the heating fraction by the young stars and the sSFR. △ Less

Submitted 7 April, 2020; originally announced April 2020.

High-resolution, 3D radiative transfer modelling. II. The early-type spiral galaxy M81

Authors: Sam Verstocken, Angelos Nersesian, Maarten Baes, Sébastien Viaene, Simone Bianchi, Viviana Casasola, Christopher J. R. Clark, Jonathan I. Davies, Ilse De Looze, Pieter De Vis, Wouter Dobbels, FrédÉric Galliano, Anthony P. Jones, Suzanne C. Madden, Aleksandr V. Mosenkov, Ana Trčka, Emmanuel M. Xilouris

Abstract: Interstellar dust absorbs stellar light very efficiently and thus shapes the energetic output of galaxies. Studying the impact of different stellar populations on the dust heating remains hard because it requires decoupling the relative geometry of stars and dust, and involves complex processes as scattering and non-local dust heating. We aim to constrain the relative distribution of dust and stel… ▽ More Interstellar dust absorbs stellar light very efficiently and thus shapes the energetic output of galaxies. Studying the impact of different stellar populations on the dust heating remains hard because it requires decoupling the relative geometry of stars and dust, and involves complex processes as scattering and non-local dust heating. We aim to constrain the relative distribution of dust and stellar populations in the spiral galaxy M81 and create a realistic model of the radiation field that describes the observations. Investigating the dust-starlight interaction on local scales, we want to quantify the contribution of young and old stellar populations to the dust heating. We aim to standardise the setup and model selection of such inverse radiative transfer simulations so this can be used for comparable modelling of other nearby galaxies. We present a semi-automated radiative transfer modelling pipeline that implements the necessary steps such as the geometric model construction and the normalisation of the components through an optimisation routine. We use the Monte Carlo radiative transfer code SKIRT to calculate a self-consistent, panchromatic model of the interstellar radiation field. By looking at different stellar populations independently, we can quantify to what extent different stellar age populations contribute to the dust heating. Our method takes into account the effects of non-local heating. We obtain a realistic 3D radiative transfer model of the face-on galaxy M81. We find that only 50.2\% of the dust heating can be attributed to young stellar populations. We confirm a tight correlation between the specific star formation rate and the heating fraction by young stellar populations, both in sky projection and in 3D, also found for radiative transfer models of M31 and M51. We conclude that... (abridged) △ Less

Submitted 7 April, 2020; originally announced April 2020.

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

Journal ref: A&A 637, A24 (2020)

Optical, X-ray, and $γ$-ray observations of the candidate transitional millisecond pulsar 4FGL J0427.8-6704

Authors: M. R. Kennedy, R. P. Breton, C. J. Clark, V. S. Dhillon, M. Kerr, D. A. H. Buckley, S. B. Potter, D. Mata Sánchez, J. G. Stringer, T. R. Marsh

Abstract: We present an optical, X-ray, and $γ$-ray study of 1SXPS J042749.2-670434, an eclipsing X-ray binary which has an associated $γ$-ray counterpart, 4FGL J0427.8-6704. This association has led to the source being classified as a transitional millisecond pulsar (tMSP) in an accreting state. We analyse 10.5 years of Fermi LAT data, and detect a $γ$-ray eclipse at the same phase as optical and X-ray ecl… ▽ More We present an optical, X-ray, and $γ$-ray study of 1SXPS J042749.2-670434, an eclipsing X-ray binary which has an associated $γ$-ray counterpart, 4FGL J0427.8-6704. This association has led to the source being classified as a transitional millisecond pulsar (tMSP) in an accreting state. We analyse 10.5 years of Fermi LAT data, and detect a $γ$-ray eclipse at the same phase as optical and X-ray eclipses at the >5$σ$ level, a significant improvement on the 2.8$σ$level of the previous detection. The confirmation of this eclipse solidifies the association between the X-ray source and the $γ$-ray source, strengthening the tMSP classification. However, analysis of several optical data sets and an X-ray observation do not reveal a change in the source's median brightness over long timescales or a bi-modality on short timescales. Instead, the light curve is dominated by flickering which has a correlation time of 2.6 min alongside a potential quasi-periodic oscillation at $\sim$21 min. The mass of the primary and secondary star are constrained to be $M_1=1.43^{+0.33}_{-0.19}$ M$_{\odot}$ and $M_2=0.3^{+0.17}_{-0.12}$ M$_{\odot}$ through modelling of the optical light curve. While this is still consistent with a white dwarf primary, we favour the transitional millisecond pulsar in a low accretion state classification due to the significance of the $γ$-ray eclipse detection. △ Less

Submitted 30 March, 2020; originally announced March 2020.

Comments: 16 pages, 14 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

The ISM scaling relations in DustPedia late-type galaxies: A benchmark study for the Local Universe

Authors: Viviana Casasola, Simone Bianchi, Pieter De Vis, Laura Magrini, Edvige Corbelli, Christopher J. R. Clark, Jacopo Fritz, Angelos Nersesian, Sebastien Viaene, Maarten Baes, Letizia P. Cassara', Jon Davies, Ilse De Looze, Wouter Dobbels, Maud Galametz, Frederic Galliano, Anthony P. Jones, Suzanne C. Madden, Aleksandr V. Mosenkov, Ana Trcka, Emmanuel Xilouris

Abstract: The purpose of this work is the characterization of the main scaling relations between all the ISM components (dust, atomic/molecular/total gas), gas-phase metallicity, and other galaxy properties, such as Mstar and galaxy morphology, for late-type galaxies in the Local Universe. This study is performed by extracting late-type galaxies from the entire DustPedia sample and by exploiting the large a… ▽ More The purpose of this work is the characterization of the main scaling relations between all the ISM components (dust, atomic/molecular/total gas), gas-phase metallicity, and other galaxy properties, such as Mstar and galaxy morphology, for late-type galaxies in the Local Universe. This study is performed by extracting late-type galaxies from the entire DustPedia sample and by exploiting the large and homogeneous dataset available thanks to the DustPedia project. The sample consists of 436 galaxies with morphological stage from T = 1 to 10, Mstar from 6 x 10^7 to 3 x 10^11 Msun, SFR from 6 x 10^(-4) to 60 Msun/yr, and 12 + log(O/H) from 8 to 9.5. The scaling relations involving the molecular gas are studied by assuming both a constant and a metallicity-dependent XCO. The analysis has been performed by means of the survival analysis technique. We confirm that the dust mass correlates very well with the total gas mass, and find -- for the first time -- that the dust mass correlates better with the atomic gas mass than the molecular one. We characterize important mass ratios such as gas fraction, molecular-to-atomic gas mass ratio, dust-to-total gas mass ratio (DGR), and dust-to-stellar mass ratio. Only the assumption of a metallicity-dependent XCO reproduces the expected decrease of the DGR with increasing morphological stage and decreasing gas-phase metallicity, with a slope of about 1. DGR, gas-phase metallicity, and the dust-to-stellar mass ratio are, for our galaxy sample, directly linked to the galaxy morphology. The molecular-to-atomic gas mass ratio and the DGR show a positive correlation for low molecular gas fractions, but for molecular gas rich galaxies this trend breaks down. This trend has never been found previously, to our knowledge. It provides new constraints for theoretical models of galaxy evolution and a reference for high-redshift studies. △ Less

Submitted 19 December, 2019; v1 submitted 20 November, 2019; originally announced November 2019.

Comments: Accepted for publication in A&A, 21 pages, 14 Figures, 7 Tables, 2 Appendices (with 3 Tables). Published version, unchanged results with respect to the previously posted version

Journal ref: A&A 633, A100 (2020)

Estimating the molecular gas mass of low-redshift galaxies from a combination of mid-infrared luminosity and optical properties

Authors: Yang Gao, Ting Xiao, Cheng Li, Xue-Jian Jiang, Qing-hua Tan, Yu Gao, Christine D. Wilson, Martin Bureau, Amelie Saintonge, Jos'e R. S'anchez-Gallego, Toby Brown, Christopher J. Clark, Ho Seong Hwang, Isabella Lamperti, Lin Lin, Lijie Liu, Dengrong Lu, Hsi-An Pan, Jixian Sun, Thomas G. Williams

Abstract: We present CO(J=1-0) and/or CO(J=2-1) spectroscopy for 31 galaxies selected from the ongoing MaNGA survey, obtained with multiple telescopes. This sample is combined with CO observations from the literature to study the correlation of the CO luminosities ($L_{\rm CO(1-0)}$) with the mid-infrared luminosities at 12 ($L_{12 μm}$) and 22 $μ$m ($L_{\rm 22 μm}$), as well as the dependence of the residu… ▽ More We present CO(J=1-0) and/or CO(J=2-1) spectroscopy for 31 galaxies selected from the ongoing MaNGA survey, obtained with multiple telescopes. This sample is combined with CO observations from the literature to study the correlation of the CO luminosities ($L_{\rm CO(1-0)}$) with the mid-infrared luminosities at 12 ($L_{12 μm}$) and 22 $μ$m ($L_{\rm 22 μm}$), as well as the dependence of the residuals on a variety of galaxy properties. The correlation with $L_{\rm 12 μm}$ is tighter and more linear, but galaxies with relatively low stellar masses and blue colors fall significantly below the mean $L_{\rm CO(1-0)}-L_{\rm 12μm}$ relation. We propose a new estimator of the CO(1-0) luminosity (and thus the total molecular gas mass) that is a linear combination of three parameters: $L_{\rm 12 μm}$, $M_\ast$ and $g-r$. We show that, with a scatter of only 0.18 dex in log $(L_{\rm CO(1-0)})$, this estimator provides unbiased estimates for galaxies of different properties and types. An immediate application of this estimator to a compiled sample of galaxies with only CO(J=2-1) observations yields a distribution of the CO(J=2-1) to CO(J=1-0) luminosity ratios ($R21$) that agrees well with the distribution of real observations, in terms of both the median and the shape. Application of our estimator to the current MaNGA sample reveals a gas-poor population of galaxies that are predominantly early-type and show no correlation between molecular gas-to-stellar mass ratio and star formation rate, in contrast to gas-rich galaxies. We also provide alternative estimators with similar scatters, based on $r$ and/or $z$ band luminosities instead of $M_\ast$. These estimators serve as cheap and convenient $M_{\rm mol}$ proxies to be potentially applied to large samples of galaxies, thus allowing statistical studies of gas-related processes of galaxies. △ Less

Submitted 10 November, 2019; v1 submitted 7 November, 2019; originally announced November 2019.

Comments: 24 pages, 12 figures, accepted for publication in ApJ

Predicting the global far-infrared SED of galaxies via machine learning techniques

Authors: W. Dobbels, M. Baes, S. Viaene, S. Bianchi, J. I. Davies, V. Casasola, C. J. R. Clark, J. Fritz, M. Galametz, F. Galliano, A. Mosenkov, A. Nersesian, A. Trčka

Abstract: Dust plays an important role in shaping a galaxy's spectral energy distribution (SED). It absorbs ultraviolet (UV) to near-infrared (NIR) radiation and re-emits this energy in the far-infrared (FIR). The FIR is essential to understand dust in galaxies. However, deep FIR observations require a space mission, none of which are still active today. We aim to infer the FIR emission across six Herschel… ▽ More Dust plays an important role in shaping a galaxy's spectral energy distribution (SED). It absorbs ultraviolet (UV) to near-infrared (NIR) radiation and re-emits this energy in the far-infrared (FIR). The FIR is essential to understand dust in galaxies. However, deep FIR observations require a space mission, none of which are still active today. We aim to infer the FIR emission across six Herschel bands, along with dust luminosity, mass, and effective temperature, based on the available UV to mid-infrared (MIR) observations. We also want to estimate the uncertainties of these predictions, compare our method to energy balance SED fitting, and determine possible limitations of the model. We propose a machine learning framework to predict the FIR fluxes from 14 UV-MIR broadband fluxes. We used a low redshift sample by combining DustPedia and H-ATLAS, and extracted Bayesian flux posteriors through SED fitting. We trained shallow neural networks to predict the far-infrared fluxes, uncertainties, and dust properties. We evaluated them on a test set using a root mean square error (RMSE) in log-space. Our results (RMSE = 0.19 dex) significantly outperform UV-MIR energy balance SED fitting (RMSE = 0.38 dex), and are inherently unbiased. We can identify when the predictions are off, for example when the input has large uncertainties on WISE 22, or when the input does not resemble the training set. The galaxies for which we have UV-FIR observations can be used as a blueprint for galaxies that lack FIR data. This results in a 'virtual FIR telescope', which can be applied to large optical-MIR galaxy samples. This helps bridge the gap until the next FIR mission. △ Less

Submitted 14 October, 2019; originally announced October 2019.

Comments: 25 pages, 16 figures, accepted by A&A. Interactive figures + jupyter notebooks at https://wdobbels.github.io/FIREnet/

Journal ref: A&A 634, A57 (2020)