Search

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

We select 1,052,469 (754,635) thin disk stars from {\it Gaia} eDR3 and LAMOST DR7 in the range of Galactocentric radius $R$ (guiding center radius $R_\mathrm{g}$) from 8 to 11\,kpc to investigate the asymmetries between the North and South of the disk midplane. More specifically we analyze the vertical velocity dispersion profiles ($\sigma_{v_{z}}(z$)) in different bins of $R$ ($R_\mathrm{g}$) and $[\mathrm{Fe/H}]$. We find troughs in the profiles of $\sigma_{v_{z}}(z)$ located in both the North ($z \sim 0.7$\,kpc) and South ($z \sim -0.5$\,kpc) of the disk at all radial and chemical bins studied. The difference between the Northern and Southern vertical velocity dispersion profiles ($\Delta\sigma_{v_{z}}(|z|)$) shows a shift between curves of different $R$ and $R_\mathrm{g}$. A similar shift exists in these NS asymmetry profiles further divided into different $[\mathrm{Fe/H}]$ ranges. The sample binned with $R_\mathrm{g}$ more clearly displays the features in the velocity dispersion profiles. The shift in the peaks of the $\Delta\sigma_{v_{z}}$ profiles and the variation in the phase spiral shape binned by metallicity indicate the variation of the vertical potential profiles and the radial metallicity gradient. The wave-like signal in NS asymmetry of $\sigma_{v_{z}}(z)$ largely originates from phase spiral; while the NS asymmetry profiles of [Fe/H] only display a weak wave-like feature near solar radius. We perform a test particle simulation to qualitatively reproduce the observed results. A quantitative explanation of the NS asymmetry in the metallicity profile needs careful consideration of the spiral shape and the perturbation model, and we leave this for future work., Comment: 13 pages, 12 figures, accepted for publication in MNRAS

The 5n-vector ensemble method is a multiple test for the targeted search of continuous gravitational waves from an ensemble of known pulsars. This method can improve the detection probability combining the results from individually undetectable pulsars if few signals are near the detection threshold. In this paper, we apply the 5n-vector ensemble method to the O3 data set from the LIGO and Virgo detectors considering an ensemble of 201 known pulsars. We find no evidence for a signal from the ensemble and set a 95% credible upper limit on the mean ellipticity assuming a common exponential distribution for the pulsars' ellipticities. Using two independent hierarchical Bayesian procedures, we find upper limits of $1.2 \times 10^{-9}$ and $2.5 \times 10^{-9}$ on the mean ellipticity for the 201 analyzed pulsars.

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., Comment: 142 pages. Accepted by the Astrophysical Journal Supplement

As the sample size of repeating fast radio bursts (FRBs) has grown, an increasing diversity of phenomenology has emerged. Through long-term multi-epoch studies of repeating FRBs, it is possible to assess which phenomena are common to the population and which are unique to individual sources. We present a multi-epoch monitoring campaign of the repeating FRB source 20180301A using the ultra-wideband low (UWL) receiver observations with Murriyang, the Parkes 64-m radio telescope. The observations covered a wide frequency band spanning approximately 0.7--4 GHz, and yielded the detection of 46 bursts. None of the repeat bursts displayed radio emission in the range of 1.8--4 GHz, while the burst emission peaked at 1.1 GHz. We discover evidence for secular trends in the burst dispersion measure, indicating a decline at a rate of $-2.7\pm0.2\,{\rm pc\,cm^{-3}\,yr^{-1}}$. We also found significant variation in the Faraday rotation measure of the bursts across the follow-up period, including evidence of a sign reversal. While a majority of bursts did not exhibit any polarization, those that did show a decrease in the linear polarization fraction as a function of frequency, consistent with spectral depolarization due to scattering, as observed in other repeating FRB sources. Surprisingly, no significant variation in the polarization position angles was found, which is in contrast with earlier measurements reported for the FRB source. We measure the burst rate and sub-pulse drift rate variation and compare them with the previous results. These novel observations, along with the extreme polarization properties observed in other repeating FRBs, suggest that a sub-sample of FRB progenitors possess highly dynamic magneto-ionic environments., Comment: 21 pages, 14 figures; accepted for publication in MNRAS

The Survey for Distant Solar Twins (SDST) aims to find stars very similar to the Sun at distances 1-4 kpc, several times more distant than any currently known solar twins and analogues. The goal is to identify the best stars with which to test whether the fine-structure constant, alpha, varies with dark matter density in our Galaxy. Here we use EPIC, our line-by-line differential technique, to measure the stellar parameters - effective temperature Teff, surface gravity log g, metallicity [Fe/H] - from moderate resolution (R < 32,000) spectra of 877 solar twin and analogue candidates (547 at 1-4 kpc) observed with the HERMES spectrograph on the Anglo-Australian Telescope. These are consistent with expectations for Teff and log g from photometry, and for [Fe/H] from the Besancon stellar population model. EPIC provides small enough uncertainties (~90 K, 0.08 dex, 0.05 dex, respectively), even at the low signal-to-noise ratios available (S/N >~ 25 per pixel), to identify 299 new solar analogues (> 90% confidence), and 20 solar twins (>50% confidence), 206 and 12 of which are at 1-4 kpc. By extending EPIC to measure line broadening and lithium abundance from HERMES spectra, and with ages derived from isochrone fitting with our stellar parameters, we identify 174 solar analogues at 1-4 kpc which are relatively inactive, slowly rotating, and with no evidence of spectroscopic binarity. These are the preferred targets for follow-up spectroscopy to measure alpha., Comment: 12 pages, 12 figures. Accepted by MNRAS. The full table of stars used in this work can be found at https://doi.org/10.5281/zenodo.7332606

The fine structure constant, $\alpha$, sets the strength of the electromagnetic force. The Standard Model of particle physics provides no explanation for its value, which could potentially vary. The wavelengths of stellar absorption lines depend on $\alpha$, but are subject to systematic effects owing to astrophysical processes in stellar atmospheres. We measured precise line wavelengths using 17 stars, selected to have almost identical atmospheric properties to those of the Sun (solar twins), which reduces those systematic effects. We found that $\alpha$ varies by $\lesssim$50 parts-per-billion (ppb) within 50 parsecs from Earth. Combining the results from all 17 stars provides an empirical, local reference for stellar measurements of $\alpha$ with an ensemble precision of 12 ppb., Comment: 33 pages, 6 figures. Published in Science (11 November 2022). This is the accepted version which includes 20 pages of Supplementary Materials

Studies of solar twins have key impacts on the astronomical community, but only $\sim$100--200 nearby solar twins ($<$ 1 kpc) have been reliably identified over the last few decades. The aim of our survey (SDST) is to identify $\sim$150--200 distant solar twins and analogues (up to $\lesssim$ 4 kpc) closer to the Galactic Centre. We took advantage of the precise Gaia and Skymapper surveys to select Sun-like candidates in a 2-degree field, which were observed with the HERMES spectrograph on the Anglo-Australian Telescope. We successfully built up the required signal-to-noise ratio (25-per-pixel in the HERMES red band) for most targets as faint as Gaia G of 17.4 mag. The stellar photometric/astrometric parameters (e.g., \teff, \logg, mass) of our candidates are derived in this paper, while the spectroscopic parameters will be presented in the third paper in this SDST series. The selection success rate - the fraction of targets which belong to solar twins or analogues - was estimated from simulated survey data and the Besan\c{c}on stellar population model, and compared with the actual success rate of the survey. We find that expected and actual success rates agree well, indicating that the numbers of solar twins and analogues we discover in SDST are consistent with expectations, affirming the survey approach. These distant solar analogues are prime targets for testing for any variation in the strength of electromagnetism in regions of higher dark matter density, and can make additional contributions to our understanding of, e.g., Galactic chemical evolution in the inner Milky Way., Comment: 15 pages, 13 figures. Accepted by MNRAS

The rich absorption spectra of Sun-like stars are enticing probes for variations in the fine-structure constant, $\alpha$, which gauges the strength of electromagnetism. While individual line wavelengths are sensitive to $\alpha$, they are also sensitive to physical processes in the stellar atmospheres, which has precluded their use so far. Here we demonstrate a new, differential approach using solar twins: velocity separations between close pairs of transitions are compared across stars with very similar physical properties, strongly suppressing astrophysical and instrumental systematic errors. We utilise 423 archival exposures of 18 solar twins from the High-Accuracy Radial velocity Planetary Searcher (HARPS), in which calibration errors can be reduced to $\lesssim$3 m/s. For stars with $\approx$10 high signal-to-noise ratio spectra ($\ge$200 per pixel), velocity separations between pairs are measured with $\approx$10 m/s statistical precision. A companion paper assesses a range of systematic error sources using 130 stars, with a greater range of stellar parameters, providing accurate corrections for astrophysical effects and a residual, intrinsic star-to-star scatter of 0-13 m/s. Within these uncertainties, we find no evidence for velocity separation differences in 17 transition pairs between solar twins. In a second companion paper, this is found to limit local ($\lesssim$50 pc) variations in $\alpha$ to $\approx$50 parts per billion, $\sim$2 orders of magnitude less than other Galactic constraints., Comment: 15 pages, 5 figures. Accepted by MNRAS. The VarConLib (Varying Constants Library) software used in this work is available at https://github.com/DBerke/varconlib. The measurements and results in this work are available at https://github.com/DBerke/Berke_et_alia_2022_supplemental_data

Sun-like stars are a new probe of variations in the fine-structure constant, $\alpha$, via the solar twins approach: velocity separations of close pairs of absorption lines are compared between stars with very similar stellar parameters, i.e. effective temperature, metallicity and surface gravity within 100K, 0.1 dex and 0.2 dex of the Sun's values. Here we assess possible systematic errors in this approach by analysing $\gtrsim$10,000 archival exposures from the High-Accuracy Radial velocity Planetary Searcher (HARPS) of 130 stars covering a much broader range of stellar parameters. We find that each transition pair's separation shows broad, low-order variations with stellar parameters which can be accurately modelled, leaving only a small residual, intrinsic star-to-star scatter of 0-33 m/s (average $\approx$7 m/s, $\approx$10$^{-4}$\r{A} at 5000\r{A}). This limits the precision available from a single pair in one star. We consider potential systematic errors from a range of instrumental and astrophysical sources (e.g. wavelength calibration, charge transfer inefficiency, stellar magnetic activity, line blending) and conclude that variations in elemental abundances, isotope ratios and stellar rotational velocities may explain this star-to-star scatter. Finally, we find that the solar twins approach can be extended to solar analogues - within 300K, 0.3 dex and 0.4 dex of the Sun's parameters - without significant additional systematic errors, allowing a much larger number of stars to be used as probes of variation in $\alpha$, including at much larger distances., Comment: 18 pages, 12 figures. Accepted by MNRAS. The VarConLib (Varying Constants Library) software used in this work is available at https://github.com/DBerke/varconlib. The measurements and results in this work are available at https://github.com/DBerke/Berke_et_alia_2022_supplemental_data

We report the detection of FRB20191107B with the UTMOST radio telescope at a dispersion measure (DM) of 714.9 ${\rm pc~cm^{-3}}$. The burst consists of three components, the brightest of which has an intrinsic width of only 11.3 $\mu$s and a scattering tail with an exponentially decaying time-scale of 21.4 $\mu$s measured at 835 MHz. We model the sensitivity of UTMOST and other major FRB surveys to such narrow events. We find that $>60\%$ of FRBs like FRB20191107B are being missed, and that a significant population of very narrow FRBs probably exists and remains underrepresented in these surveys. The high DM and small scattering timescale of FRB20191107B allows us to place an upper limit on the strength of turbulence in the Intergalactic Medium (IGM), quantified as scattering measure (SM), of ${\rm SM_{IGM} < 8.4 \times 10^{-7} ~kpc~m^{-20/3}}$. Almost all UTMOST FRBs have full phase information due to real-time voltage capture which provides us with the largest sample of coherently dedispersed single burst FRBs. Our 10.24 $\mu$s time resolution data yields accurately measured FRB scattering timescales. We combine the UTMOST FRBs with 10 FRBs from the literature and find no obvious evidence for a DM-scattering relation, suggesting that IGM is not the dominant source of scattering in FRBs. We support the results of previous studies and identify the local environment of the source in the host galaxy as the most likely region which dominates the observed scattering of our FRBs.

We measure the enclosed Milky Way mass profile to Galactocentric distances of $\sim70$ and $\sim50$ kpc using the smooth, diffuse stellar halo samples of Bird et al. The samples are LAMOST and SDSS/SEGUE K giants (KG) and SDSS/SEGUE blue horizontal branch (BHB) stars with accurate metallicities. The 3D kinematics are available through LAMOST and SDSS/SEGUE distances and radial velocities and {\it Gaia} DR2 proper motions. Two methods are used to estimate the enclosed mass: 3D spherical Jeans equation and Evans et al. tracer mass estimator (TME). We remove substructure via the Xue et al. method based on integrals of motion. We evaluate the uncertainties on our estimates due to random sampling noise, systematic distance errors, the adopted density profile, and non-virialization and non-spherical effects of the halo. The tracer density profile remains a limiting systematic in our mass estimates, although within these limits we find reasonable agreement across the different samples and the methods applied. Out to $\sim70$ and $\sim50$ kpc, the Jeans method yields total enclosed masses of $4.3\pm0.95$ (random) $\pm0.6$ (systematic) $\times10^{11}$ M$_\odot$ and $4.1\pm1.2$ (random) $\pm0.6$ (systematic) $\times10^{11}$ M$_\odot$ for the KG and BHB stars, respectively. For the KG and BHB samples we find a dark matter virial mass of $M_{200}=0.55^{+0.15}_{-0.11}$ (random) $\pm0.083$ (systematic) $\times10^{12}$ M$_\odot$ and $M_{200}=1.00^{+0.67}_{-0.33}$ (random) $\pm0.15$ (systematic) $\times10^{12}$ M$_\odot$, respectively., Comment: Accepted for publication in MNRAS. Comments welcomed

Solar twins are stars of key importance to the field of astronomy and offer a multitude of science applications. Only a small number ($\lesssim200$) of solar twins are known today, all of which are relatively close to our Sun ($\lesssim800\,pc$). The goal of our Survey for Distant Solar Twins (SDST) is to identify many more solar twin and solar analogue stars out to much larger distances ($\sim4\,kpc$). In this paper, we present a new method to identify solar twins using relatively low $S/N$, medium resolving power ($R\sim 28{,}000$) spectra that will be typical of such distant targets observed with HERMES on the $3.9\,m$ Anglo-Australian Telescope (AAT). We developed a novel approach, namely EPIC, to measure stellar parameters which we use to identify stars similar to our Sun. EPIC determines the stellar atmospheric parameters (effective temperature $T_{\mathrm{eff}}$, surface gravity $\log g$ and metallicity [Fe/H]) using differential equivalent width (EW) measurements of selected spectroscopic absorption features and a simple model, trained on previously analysed spectra, that connects these EWs to the stellar parameters. The reference for the EW measurements is a high $S/N$ solar spectrum which is used to minimise several systematic effects. EPIC is fast, optimised for Sun-like stars and yields stellar parameter measurements with small enough uncertainties to enable spectroscopic identification of solar twin and analogue stars up to $\sim4\,kpc$ away using AAT/HERMES, i.e.\ $\sigma\left(T_{\mathrm{eff}}, \log g, \textrm{[Fe/H]}\right) = \left(50\,K, 0.08\,dex, 0.03\,dex\right)$ on average at $S/N=25$., Comment: Accepted by MNRAS. The current release of the EPIC algorithm is available at https://dx.doi.org/10.5281/zenodo.6029325

Fast Radio Bursts (FRBs) are highly energetic transient events with duration of order of microseconds to milliseconds and of unknown origin. They are known to lie at cosmological distances, through localisation to host galaxies. Recently, an FRB-like event was seen from the Milky Way magnetar SGR 1935+2154 by the CHIME and STARE2 telescopes. This is the only magnetar that has produced FRB events in our galaxy. Finding similar events in the Milky Way is of great interest to understanding FRB progenitors. Such events will be strongly affected by the turbulent interstellar medium in the Milky Way, their intrinsic energy distribution and their spatial locations within the plane of the Milky Way. We examine these effects using models for the distribution of electrons in the ISM to estimate the dispersion measure and pulse scattering of mock events, and a range of models for the spatial distribution and luminosity functions, including models motivated by the spatial distribution of the Milky Way's magnetars. We evaluate the fraction of FRB events in the Milky Way that are detectable by STARE2 for a range of ISM models, spatial distributions and burst luminosity functions. In all the models examined, only a fraction of burst events are detectable, mainly due to the scattering effects of the ISM. We find that GReX, a proposed all-sky experiment, could increase the detection rate of Milky Way FRB events by an order of magnitude, depending on assumptions made about the luminosity function and scale-height of the FRBs., Comment: 11 pages, 14 figures, accepted for publication in MNRAS

Fast Radio Bursts (FRBs) are a promising tool for studying the low-density universe as their dispersion measures (DM) are extremely sensitive probes of electron column density. Active Galactic Nuclei (AGN) inject energy into the intergalactic medium, affecting the DM and its scatter. To determine the effectiveness of FRBs as a probe of AGN feedback, we analysed three different AGN models from the EAGLE simulation series. We measured the mean DM-redshift relation, and the scatter around it, using $2.56 \times 10^8$ sightlines at 131 redshift ($z$) bins between $0 \leq z \leq 3$. While the DM-redshift relation itself is highly robust against different AGN feedback models, significant differences are detected in the scatter around the mean: weaker feedback leads to more scatter. We find $\sim 10^4$ localised FRBs are needed to discriminate between the scatter in standard feedback and stronger, more intermittent feedback models. The number of FRBs required is dependent on the redshift distribution of the detected population. A log-normal redshift distribution at $z=0.5$ requires approximately 50% fewer localised FRBs than a distribution centred at $z=1$. With the Square Kilometre Array expected to detect $>10^3$ FRBs per day, in the future, FRBs will be able to provide constraints on AGN feedback., Comment: 5 pages, 4 figures; Submitted to MNRAS Letters. Comments welcome

The mechanism that produces fast radio burst (FRB) emission is poorly understood. Targeted monitoring of repeating FRB sources provides the opportunity to fully characterize the emission properties in a manner impossible with one-off bursts. Here, we report observations of the source of FRB 20201124A, with the Australian Square Kilometre Array Pathfinder (ASKAP) and the ultra-wideband low (UWL) receiver at the Parkes 64-m radio telescope (Murriyang). The source entered a period of emitting bright bursts during early 2021 April. We have detected 16 bursts from this source. One of the bursts detected with ASKAP is the brightest burst ever observed from a repeating FRB source with an inferred fluence of $640\pm70$ Jy ms. Of the five bursts detected with the Parkes UWL, none display any emission in the range 1.1--4 GHz. All UWL bursts are highly polarized, with their Faraday rotation measures (RMs) showing apparent variations. We obtain an average RM of $-614$ rad m$^{-2}$ for this FRB source with a standard deviation of $16$ rad m$^{-2}$ in the UWL bursts. In one of the UWL bursts, we see evidence of significant circularly polarized emission with a fractional extent of $47\pm1$ per cent. Such a high degree of circular polarization has never been seen before in bursts from repeating FRB sources. We also see evidence for significant variation in the linear polarization position angle in the pulse profile of this UWL repeat burst. Models for repeat burst emission will need to account for the increasing diversity in the burst polarization properties., Comment: 14 pages, 8 figures; accepted for publication in MNRAS

We present the localization and host galaxies of one repeating and two apparently non-repeating Fast Radio Bursts. FRB20180301A was detected and localized with the Karl G. Jansky Very Large Array to a star-forming galaxy at $z=0.3304$. FRB20191228A, and FRB20200906A were detected and localized by the Australian Square Kilometre Array Pathfinder to host galaxies at $z=0.2430$ and $z=0.3688$, respectively. We combine these with 13 other well-localized FRBs in the literature, and analyze the host galaxy properties. We find no significant differences in the host properties of repeating and apparently non-repeating FRBs. FRB hosts are moderately star-forming, with masses slightly offset from the star-forming main-sequence. Star formation and low-ionization nuclear emission-line region (LINER) emission are major sources of ionization in FRB host galaxies, with the former dominant in repeating FRB hosts. FRB hosts do not track stellar mass and star formation as seen in field galaxies (more than 95% confidence). FRBs are rare in massive red galaxies, suggesting that progenitor formation channels are not solely dominated by delayed channels which lag star formation by Gigayears. The global properties of FRB hosts are indistinguishable from core-collapse supernovae (CCSNe) and short gamma-ray bursts (SGRBs) hosts, and the spatial offset (from galaxy centers) of FRBs is mostly inconsistent with that of the Galactic neutron star population (95% confidence). The spatial offsets of FRBs (normalized to the galaxy effective radius) also differ from those of globular clusters (GCs) in late- and early-type galaxies with 95% confidence., Comment: 27 pages, 12 figures, 8 tables. Accepted for publication in the Astrophysical Journal

Galactic electron density distribution models are crucial tools for estimating the impact of the ionised interstellar medium on the impulsive signals from radio pulsars and fast radio bursts. The two prevailing Galactic electron density models are YMW16 (Yao et al., 2017) and NE2001 (Cordes & Lazio, 2002). Here, we introduce a software package PyGEDM which provides a unified application programming interface (API) for these models and the YT20 (Yamasaki & Totani, 2020) model of the Galactic halo. We use PyGEDM to compute all-sky maps of Galactic dispersion measure (DM) for YMW16 and NE2001, and compare the large-scale differences between the two. In general, YMW16 predicts higher DM values toward the Galactic anticentre. YMW16 predicts higher DMs at low Galactic latitudes, but NE2001 predicts higher DMs in most other directions. We identify lines of sight for which the models are most discrepant, using pulsars with independent distance measurements. YMW16 performs better on average than NE2001, but both models show significant outliers. We suggest that future campaigns to determine pulsar distances should focus on targets where the models show large discrepancies, so future models can use those measurements to better estimate distances along those line of sight. We also suggest that the Galactic halo should be considered as a component in future GEDMs, to avoid overestimating the Galactic DM contribution for extragalactic sources such as FRBs., Comment: 13 pages, 7 figures. Accepted to PASA

We present our 500 pc distance-limited study of stellar fares using the Dark Energy Camera as part of the Deeper, Wider, Faster Program. The data was collected via continuous 20-second cadence g band imaging and we identify 19,914 sources with precise distances from Gaia DR2 within twelve, ~3 square-degree, fields over a range of Galactic latitudes. An average of ~74 minutes is spent on each field per visit. All light curves were accessed through a novel unsupervised machine learning technique designed for anomaly detection. We identify 96 flare events occurring across 80 stars, the majority of which are M dwarfs. Integrated are energies range from $\sim 10^{31}-10^{37}$ erg, with a proportional relationship existing between increased are energy with increased distance from the Galactic plane, representative of stellar age leading to declining yet more energetic are events. In agreement with previous studies we observe an increase in flaring fraction from M0 -> M6 spectral types. Furthermore, we find a decrease in the flaring fraction of stars as vertical distance from the galactic plane is increased, with a steep decline present around ~100 pc. We find that ~70% of identified flares occur on short timescales of ~8 minutes. Finally we present our associated are rates, finding a volumetric rate of $2.9 \pm 0.3 \times 10^{-6}$ flares pc$^{-3}$ hr$^{-1}$., Comment: 20 pages, 14 figures, 6 tables

The dispersion measure (DM) of fast radio bursts (FRBs) provides a unique way to probe ionised baryons in the intergalactic medium (IGM). Cosmological models with different parameters lead to different DM-redshift ($\mathrm{DM}-z$) relations. Additionally, the over/under-dense regions in the IGM and the circumgalactic medium of intervening galaxies lead to scatter around the mean $\mathrm{DM}-z$ relations. We have used the Evolution and Assembly of GaLaxies and their Environments (EAGLE) simulations to measure the mean $\mathrm{DM}-z$ relation and the scatter around it using over one billion lines-of-sight between redshifts $0

Dedicated surveys using different detection pipelines are being carried out at multiple observatories to find more Fast Radio Bursts (FRBs). Understanding the efficiency of detection algorithms and the survey completeness function is important to enable unbiased estimation of the underlying FRB population properties. One method to achieve end-to-end testing of the system is by injecting mock FRBs in the live data-stream and searching for them blindly. Mock FRB injection is particularly effective for machine-learning-based classifiers, for which analytic characterisation is impractical. We describe a first-of-its-kind implementation of a real-time mock FRB injection system at the upgraded Molonglo Observatory Synthesis Telescope (UTMOST) and present our results for a set of 20,000 mock FRB injections. The injections have yielded clear insight into the detection efficiencies and have provided a survey completeness function for pulse width, fluence and DM. Mock FRBs are recovered with uniform efficiency over the full range of injected DMs, however the recovery fraction is found to be a strong function of the width and Signal-to-Noise (SNR). For low widths ($\lesssim 20$ ms) and high SNR ($\gtrsim$ 9) the recovery is highly effective with recovery fractions exceeding 90%. We find that the presence of radio frequency interference causes the recovered SNR values to be systematically lower by up to 20% compared to the injected values. We find that wider FRBs become increasingly hard to recover for the machine-learning-based classifier employed at UTMOST. We encourage other observatories to implement live injection set-ups for similar testing of their surveys., Comment: 12 pages, 9 figures. Accepted for publication in MNRAS

The fast radio burst (FRB) population is observationally divided into sources that have been observed to repeat and those that have not. There is tentative evidence that the bursts from repeating sources have different properties than the non-repeating ones. In order to determine the occurrence rate of repeating sources and characterize the nature of repeat emission, we have been conducting sensitive searches for repetitions from bursts detected with the Australian Square Kilometre Array Pathfinder (ASKAP) with the 64-m Parkes radio telescope, using the recently commissioned Ultra-wideband Low (UWL) receiver system, over a band spanning 0.7$-$4.0 GHz. We report the detection of a repeat burst from the source of FRB 20190711A. The detected burst is 1 ms wide and has a bandwidth of just 65 MHz. We find no evidence of any emission in the remaining part of the 3.3 GHz UWL band. While the emission bandwidths of the ASKAP and UWL bursts show $\nu^{-4}$ scaling consistent with a propagation effect, the spectral occupancy is inconsistent with diffractive scintillation. This detection rules out models predicting broad-band emission from the FRB 20190711A source and puts stringent constraints on the emission mechanism. The low spectral occupancy highlights the importance of sub-banded search methods in detecting FRBs., Comment: 7 pages, 5 figures; accepted for publication in MNRAS

We present a new fast radio burst (FRB) at 920 MHz discovered during commensal observations conducted with the Australian Square Kilometre Array Pathfinder (ASKAP) as part of the Commensal Real-time ASKAP Fast Transients (CRAFT) survey. FRB 191001 was detected at a dispersion measure (DM) of 506.92(4) pc cm$^{-3}$ and its measured fluence of 143(15) Jy ms is the highest of the bursts localized to host galaxies by ASKAP to date. The subarcsecond localization of the FRB provided by ASKAP reveals that the burst originated in the outskirts of a highly star-forming spiral in a galaxy pair at redshift $z=0.2340(1)$. Radio observations show no evidence for a compact persistent radio source associated with the FRB 191001 above a flux density of $15 \mu$Jy. However, we detect diffuse synchrotron radio emission from the disk of the host galaxy that we ascribe to ongoing star formation. FRB 191001 was also detected as an image-plane transient in a single 10 s snapshot with a flux density of 19.3 mJy in the low-time-resolution visibilities obtained simultaneously with CRAFT data. The commensal observation facilitated a search for repeating and slowly varying radio emissions 8 hr before and 1 hr after the burst. We found no variable radio emission on timescales ranging from 1 ms to 1.4 hr. We report our upper limits and briefly review FRB progenitor theories in the literature that predict radio afterglows. Our data are still only weakly constraining of any afterglows at the redshift of the FRB. Future commensal observations of more nearby and bright FRBs will potentially provide stronger constraints., Comment: 12 pages, 6 figures, Published in ApJ Letters. We have corrected the R.A. uncertainty in the position of FRB 191001 from 0.006s to 0.02s

Identification of anomalous light curves within time-domain surveys is often challenging. In addition, with the growing number of wide-field surveys and the volume of data produced exceeding astronomers ability for manual evaluation, outlier and anomaly detection is becoming vital for transient science. We present an unsupervised method for transient discovery using a clustering technique and the Astronomaly package. As proof of concept, we evaluate 85553 minute-cadenced light curves collected over two 1.5 hour periods as part of the Deeper, Wider, Faster program, using two different telescope dithering strategies. By combining the clustering technique HDBSCAN with the isolation forest anomaly detection algorithm via the visual interface of Astronomaly, we are able to rapidly isolate anomalous sources for further analysis. We successfully recover the known variable sources, across a range of catalogues from within the fields, and find a further 7 uncatalogued variables and two stellar flare events, including a rarely observed ultra fast flare (5 minute) from a likely M-dwarf., Comment: Accepted 7 Aug 2020, 19 pages, 8 figures

The pulse morphology of fast radio bursts (FRBs) provides key information in both understanding progenitor physics and the plasma medium through which the burst propagates. We present a study of the profiles of 33 bright FRBs detected by the Australian Square Kilometre Array Pathfinder. We identify seven FRBs with measureable intrinsic pulse widths, including two FRBs that have been seen to repeat. In our modest sample we see no evidence for bimodality in the pulse width distribution. We also identify five FRBs with evidence of millisecond timescale pulse broadening caused by scattering in inhomogeneous plasma. We find no evidence for a relationship between pulse broadening and extragalactic dispersion measure. The scattering could be either caused by extreme turbulence in the host galaxy or chance propagation through foreground galaxies. With future high time resolution observations and detailed study of host galaxy properties we may be able to probe line-of-sight turbulence on gigaparsec scales., Comment: 10 pages, 6 figures, Accepted to MNRAS

In the era of advanced electromagnetic and gravitational wave detectors, it has become increasingly important to effectively combine and study the impact of stellar evolution on binaries and dynamical systems of stars. Systematic studies dedicated to exploring uncertain parameters in stellar evolution are required to account for the recent observations of the stellar populations. We present a new approach to the commonly used Single-Star Evolution (SSE) fitting formulae, one that is more adaptable: Method of Interpolation for Single Star Evolution (METISSE). It makes use of interpolation between sets of pre-computed stellar tracks to approximate evolution parameters for a population of stars. We have used METISSE with detailed stellar tracks computed by the Modules for Experiments in Stellar Astrophysics (MESA), Bonn Evolutionary Code (BEC) and Cambridge STARS code. METISSE better reproduces stellar tracks computed using the STARS code compared to SSE, and is on average three times faster. Using stellar tracks computed with MESA and BEC, we apply METISSE to explore the differences in the remnant masses, the maximum radial expansion, and the main-sequence lifetime of massive stars. We find that different physical ingredients used in the evolution of stars, such as the treatment of radiation dominated envelopes, can impact their evolutionary outcome. For stars in the mass range 9 to 100 M$_\odot$, the predictions of remnant masses can vary by up to 20 M$_\odot$, while the maximum radial expansion achieved by a star can differ by an order of magnitude between different stellar models., Comment: Accepted for publication in MNRAS

The Australian SKA Pathfinder (ASKAP) telescope has started to localize Fast Radio Bursts (FRBs) to arcsecond accuracy from the detection of a single pulse, allowing their host galaxies to be reliably identified. We discuss the global properties of the host galaxies of the first four FRBs localized by ASKAP, which lie in the redshift range $0.11

Combining high time and frequency resolution full-polarisation spectra of Fast Radio Bursts (FRBs) with knowledge of their host galaxy properties provides an opportunity to study both the emission mechanism generating them and the impact of their propagation through their local environment, host galaxy, and the intergalactic medium. The Australian Square Kilometre Array Pathfinder (ASKAP) telescope has provided the first ensemble of bursts with this information. In this paper, we present the high time and spectral resolution, full polarisation observations of five localised FRBs to complement the results published for the previously studied ASKAP FRB~181112. We find that every FRB is highly polarised, with polarisation fractions ranging from 80 -- 100\%, and that they are generally dominated by linear polarisation. While some FRBs in our sample exhibit properties associated with an emerging archetype (i.e., repeating or apparently non-repeating), others exhibit characteristic features of both, implying the existence of a continuum of FRB properties. When examined at high time resolution, we find that all FRBs in our sample have evidence for multiple sub-components and for scattering at a level greater than expected from the Milky Way. We find no correlation between the diverse range of FRB properties (e.g., scattering time, intrinsic width, and rotation measure) and any global property of their host galaxy. The most heavily scattered bursts reside in the outskirts of their host galaxies, suggesting that the source-local environment rather than the host interstellar medium is likely the dominant origin of the scattering in our sample., Comment: 18 pages, 4 figures, 5 tables, accepted by MNRAS

We analyze the anisotropy profile of the Milky Way's smooth, diffuse stellar halo using SDSS/SEGUE blue horizontal branch stars and SDSS/SEGUE and LAMOST K giants. These intrinsically luminous stars allow us to probe the halo to approximately 100 kpc from the Galactic center. Line-of-sight velocities, distances, metallicities, and proper motions are available for all stars via SDSS/SEGUE, LAMOST, and Gaia, and we use these data to construct a full 7D set consisting of positions, space motions, and metallicity. We remove substructure from our samples using integrals of motion based on the method of Xue et al. We find radially dominated kinematic profiles with nearly constant anisotropy within 20 kpc, beyond which the anisotropy profile gently declines although remains radially dominated to the furthest extents of our sample. Independent of star type or substructure removal, the anisotropy depends on metallicity, such that the orbits of the stars become less radial with decreasing metallicity. For $-1.7<$ [Fe/H] $<-1$, the smooth, diffuse halo anisotropy profile begins to decline at Galactocentric distances $\sim20$ kpc, from $\beta\sim0.9$ to 0.7 for K giants and from $\beta\sim0.8$ to 0.1 for blue horizontal branch stars. For [Fe/H] $<-1.7$, the smooth, diffuse halo anisotropy remains constant along all distances with $0.2<\beta<0.7$ depending on the metallicity range probed, although independent on star type. These samples are ideal for estimating the total Galactic mass as they represent the virialized stellar halo system., Comment: 35 pages, 20 figures, 2 tables, accepted by ApJ, comments welcome

With the increasing numbers of large stellar survey projects, the quality and quantity of excellent tracers to study the Milky Way is rapidly growing, one of which is the classical Cepheids. Classical Cepheids are high precision standard candles with very low typical uncertainties ($<$ 3%) available via the mid-infrared period-luminosity relation. About 3500 classical Cepheids identified from OGLE, ASAS-SN, Gaia, WISE and ZTF survey data have been analyzed in this work, and their spatial distributions show a clear signature of Galactic warp. Two kinematical methods are adopted to measure the Galactic rotation curve in the Galactocentric distance range of $4\lesssim R_{\rm GC} \lesssim 19$ kpc. Gently declining rotation curves are derived by both the proper motion (PM) method and 3-dimensional velocity vector (3DV) method. The largest sample of classical Cepheids with most accurate 6D phase-space coordinates available to date are modeled in the 3DV method, and the resulting rotation curve is found to decline at the relatively smaller gradient of ($-1.33\pm0.1$) ${\rm km\,s^{-1}\,kpc^{-1}}$. Comparing to results from the PM method, a higher rotation velocity (($232.5\pm0.83$) ${\rm km\,s^{-1}}$) is derived at the position of Sun in the 3DV method. The virial mass and local dark matter density are estimated from the 3DV method which is the more reliable method, $M_{\rm vir} = (0.822\pm0.052)\times 10^{12}\,M_\odot$ and $\rho_{\rm DM,\odot} = 0.33\pm0.03$ GeV ${\rm cm^{-3}}$, respectively., Comment: 21 pages, 5 figures, 1 table. Accepted for publication in Astrophysical Journal Letters. Following Minor Revision

While pulsars possess exceptional rotational stability, large scale timing studies have revealed at least two distinct types of irregularities in their rotation: red timing noise and glitches. Using modern Bayesian techniques, we investigated the timing noise properties of 300 bright southern-sky radio pulsars that have been observed over 1.0-4.8 years by the upgraded Molonglo Observatory Synthesis Telescope (MOST). We reanalysed the spin and spin-down changes associated with nine previously reported pulsar glitches, report the discovery of three new glitches and four unusual glitch-like events in the rotational evolution of PSR J1825$-$0935. We develop a refined Bayesian framework for determining how red noise strength scales with pulsar spin frequency ($\nu$) and spin-down frequency ($\dot{\nu}$), which we apply to a sample of 280 non-recycled pulsars. With this new method and a simple power-law scaling relation, we show that red noise strength scales across the non-recycled pulsar population as $\nu^{a} |\dot{\nu}|^{b}$, where $a = -0.84^{+0.47}_{-0.49}$ and $b = 0.97^{+0.16}_{-0.19}$. This method can be easily adapted to utilise more complex, astrophysically motivated red noise models. Lastly, we highlight our timing of the double neutron star PSR J0737$-$3039, and the rediscovery of a bright radio pulsar originally found during the first Molonglo pulsar surveys with an incorrectly catalogued position., Comment: Accepted by MNRAS. 28 pages, 8 figures, 8 tables

We report the first detection of a glitch in the radio pulsar PSR J0908$-$4913 (PSR B0906$-$49) during regular timing observations by the Molonglo Observatory Synthesis Telescope (MOST) as part of the UTMOST project., Comment: 2 pages, 1 table

We have implemented prescriptions for modelling pulsars in the rapid binary population synthesis code COMPAS. We perform a detailed analysis of the double neutron star (DNS) population, accounting for radio survey selection effects. The surface magnetic field decay timescale (${\sim}1000$\,Myr) and mass scale (${\sim}0.02$\,M$_\odot$) are the dominant uncertainties in our model. Mass accretion during common envelope evolution plays a non-trivial role in recycling pulsars. We find a best-fit model that is in broad agreement with the observed Galactic DNS population. Though the pulsar parameters (period and period derivative) are strongly biased by radio selection effects, the observed orbital parameters (orbital period and eccentricity) closely represent the intrinsic distributions. The number of radio observable DNSs in the Milky Way at present is about 2500 in our model, corresponding to approximately 10\% of the predicted total number of DNSs in the galaxy. Using our model calibrated to the Galactic DNS population, we make predictions for DNS mergers observed in gravitational waves. The DNS chirp mass distribution varies from 1.1M$_\odot$ to 2.1M$_\odot$ and the median is found to be 1.14\,M$_\mathrm{\odot}$. The expected effective spin $\chi_\mathrm{eff}$ for isolated DNSs is $\lesssim$0.03 from our model. We predict that 34\% of the current Galactic isolated DNSs will merge within a Hubble time, and have a median total mass of 2.7\,M$_\mathrm{\odot}$. Finally, we discuss implications for fast radio bursts and post-merger remnant gravitational-waves., Comment: 26 pages, Accepted for publication by MNRAS

Present-day galaxies are surrounded by cool and enriched halo gas extending to hundreds of kiloparsecs. This halo gas is thought to be the dominant reservoir of material available to fuel future star formation, but direct constraints on its mass and physical properties have been difficult to obtain. We report the detection of a fast radio burst (FRB 181112) with arcsecond precision, which passes through the halo of a foreground galaxy. Analysis of the burst shows the halo gas has low net magnetization and turbulence. Our results imply predominantly diffuse gas in massive galactic halos, even those hosting active supermassive black holes, contrary to some previous results., Comment: Published in Science on 2019 September 26; Main (3 figures; 1 Table) + Supp (12 figures; 7 Tables)

The rate of fast radio bursts (FRBs) in the direction of nearby galaxy clusters is expected to be higher than the mean cosmological rate if intrinsically faint FRBs are numerous. In this paper, we describe a targeted search for faint FRBs near the core of the Virgo cluster using the Australian Square Kilometer Array Pathfinder telescope. During 300 hr of observations, we discovered one burst, FRB 180417, with dispersion measure DM $=474.8$ cm$^{-3}$pc. The FRB was promptly followed up by several radio telescopes for 27 h, but no repeat bursts were detected. An optical follow-up of FRB 180417 using the PROMPT5 telescope revealed no new sources down to an $R$-band magnitude of 20.1. We argue that FRB 180417 is likely behind the Virgo cluster as the Galactic and intracluster DM contribution are small compared to the DM of the FRB, and there are no galaxies in the line of sight. The non-detection of FRBs from Virgo constrains the faint-end slope, $\alpha<1.52$ (at 68\% confidence limit), and the minimum luminosity, $L_{\rm min}\gtrsim 2\times 10^{40}$ erg s$^{-1}$ (at 68\% confidence limit), of the FRB luminosity function assuming cosmic FRB rate of $10^4$ FRBs sky$^{-1}$ day$^{-1}$ with flux above 1 Jy located out to redshift of 1. Further FRB surveys of galaxy clusters with high-sensitivity instruments will tighten the constraints on the faint end of the luminosity function and, thus, are strongly encouraged., Comment: 9 pages, 6 figures, Accepted for publication in MNRAS

We report the detection of repeat bursts from the source of FRB 171019, one of the brightest fast radio bursts (FRBs) detected in the Australian Square Kilometre Array Pathfinder (ASKAP) fly's eye survey. Two bursts from the source were detected with the Green Bank Telescope in observations centered at 820 MHz. The repetitions are a factor of $\sim 590$ fainter than the ASKAP-discovered burst. All three bursts from this source show no evidence of scattering and have consistent pulse widths. The pulse spectra show modulation that could be evidence for either steep spectra or patchy emission. The two repetitions were the only ones found in an observing campaign for this FRB totaling 1000 hr, which also included ASKAP and the 64-m Parkes radio telescope, over a range of frequencies (720$-$2000 MHz) at epochs spanning two years. The inferred scaling of repetition rate with fluence of this source agrees with the other repeating source, FRB 121102. The detection of faint pulses from FRB 171019 shows that at least some FRBs selected from bright samples will repeat if follow-up observations are conducted with more sensitive telescopes., Comment: 9 pages, 4 figures; accepted for publication in ApJL; updated plots, revised scattering analysis

We report on a search for the host galaxy of FRB171020, the fast radio burst with the smallest recorded dispersion measure (DM $=114$ pc cm$^{-3}$) of our on-ongoing ASKAP survey. The low DM confines the burst location within a sufficiently small volume to rigorously constrain the identity of the host galaxy. We identify 16 candidate galaxies in the search volume and single out ESO 601-G036, a Sc galaxy at redshift $z=0.00867$, as the most likely host galaxy. UV and optical imaging and spectroscopy reveal this galaxy has a star-formation rate of approximately 0.1 M$_\odot$ yr$^{-1}$ and oxygen abundance $12 + \log({\rm O/H}) = 8.3 \pm 0.2$, properties remarkably consistent with the galaxy hosting the repeating FRB121102. However, in contrast to FRB121102, follow-up radio observations of ESO 601-G036 show no compact radio emission above a 5$\sigma$ limit of $L_{2.1{\rm GHz}}=3.6\times 10^{19}$ W Hz$^{-1}$. Using radio continuum observations of the field, combined with archival optical imaging data, we find no analog to the persistent radio source associated with FRB121102 within the localization region of FRB171020 out to $z=0.06$. These results suggest that FRBs are not necessarily associated with a luminous and compact radio continuum source., Comment: Submitted to ApJL

We report the detection of a glitch event in the pulsar J1709$-$4429 (also known as B1706$-$44) during regular monitoring observations with the Molonglo Observatory Synthesis Telescope (UTMOST). The glitch was found during timing operations, in which we regularly observe over 400 pulsars with up to daily cadence, while commensally searching for Rotating Radio Transients, pulsars, and FRBs. With a fractional size of $\Delta\nu/\nu \approx 52.4 \times10^{-9}$, the glitch reported here is by far the smallest known for this pulsar, attesting to the efficacy of glitch searches with high cadence using UTMOST., Comment: 3 pages, 1 figure

The anisotropy parameter $\beta$ characterizes the extent to which orbits in stellar systems are predominantly radial or tangential, and is likely to constrain, for the stellar halo of the Milky Way, scenarios for its formation and evolution. We have measured the anisotropy $\beta$ as a function of Galactocentric radius from $5-100$ kpc for over 8600 metal poor ([Fe/H] $<-1.3$) halo K giants from the LAMOST catalog with line-of-sight velocities and distances, matched to proper motions from the second $Gaia$ data release. We construct full 6-D positions and velocities for the K giants to directly measure the 3 components of the velocity dispersion $(\sigma_r, \sigma_\theta, \sigma_\phi)$ (in spherical coordinates). We find that the orbits in the halo are radial over our full Galactocentric distance range reaching over 100 kpc. The anisotropy remains remarkably unchanged with Galactocentric radius from approximately 5 to 25 kpc, with an amplitude that depends on the metallicity of the stars, dropping from $\beta \approx 0.9$ for $-1.8 \leq$ [Fe/H] $< -1.3$ (for the bulk of the stars) to $\beta \approx 0.6$ for the lowest metallicities ([Fe/H] $< -1.8$). Considering our sample as a whole, $\beta\approx0.8$ and, beyond 25 kpc, the orbits gradually become less radial and anisotropy decreases to $\beta<0.3$ past 100 kpc. Within 8 kpc, $\beta<0.8$. The measurement of anisotropy is affected by substructure and streams, particularly beyond a Galactocentric distance of approximately 25 kpc, where the Sagittarius stream is prominent in the data. These results are complimentary to recent analysis of simulations by Loebman et al. and of SDSS/$Gaia$ DR1 data by Belokurov et al., Comment: 9 pages, 5 figures, published version

We report the detection of an ultra-bright fast radio burst (FRB) from a modest, 3.4-day pilot survey with the Australian Square Kilometre Array Pathfinder. The survey was conducted in a wide-field fly's-eye configuration using the phased-array-feed technology deployed on the array to instantaneously observe an effective area of $160$ deg$^2$, and achieve an exposure totaling $13200$ deg$^2$ hr. We constrain the position of FRB 170107 to a region $8'\times8'$ in size (90% containment) and its fluence to be $58\pm6$ Jy ms. The spectrum of the burst shows a sharp cutoff above $1400$ MHz, which could be either due to scintillation or an intrinsic feature of the burst. This confirms the existence of an ultra-bright ($>20$ Jy ms) population of FRBs., Comment: 8 pages, 4 figures, accepted in ApJ letters. Version 2: Fixed galactic coordinates in Table 2

We present a hard X-ray NuSTAR observation of PSR J0437-4715, the nearest millisecond pulsar. The known pulsations at the apparent pulse period ~5.76 ms are detected at energies up to 20 keV. We measure a photon index $\Gamma= 1.65\pm0.24$ (90% confidence) for the power law fit to the non-thermal emission. It had been shown that spectral models with two or three thermal components fit the XMM-Newton spectrum of PSR J0437-4715, depending on the slope of the power-law component, and the amount of absorption of soft X-rays. The new constraint on the high-energy emission provided by NuSTAR removes ambiguities regarding the thermal components of the emission below 3 keV. We performed a simultaneous spectral analysis of the XMM-Newton and NuSTAR data to confirm that three thermal components and a power law are required to fit the 0.3-20 keV emission of PSR J0437-4715. Adding a ROSAT-PSPC spectrum further confirmed this result and allowed us to better constrain the temperatures of the three thermal components. A phase-resolved analysis of the NuSTAR data revealed no significant change in the photon index of the high-energy emission. This NuSTAR observation provides further impetus for future observations with the NICER mission (Neutron Star Interior Composition Explorer) whose sensitivity will provide much stricter constraints on the equation of state of nuclear matter by combining model fits to the pulsar's phase-folded lightcurve with the pulsar's well-defined mass and distance from radio timing observations., Comment: Accepted for publication in MNRAS. 12 pages, 8 figures, 3 table

We test the long-term kinematical stability of a Galactic stellar halo model, due to Kafle, et al. (2012), who study the kinematics of approximately 5000 blue horizontal branch (BHB) stars in the Sloan Digital Sky Survey (SDSS). The velocity dispersion $\sigma$ and anisotropy parameter $\beta$ of the stars have been determined as functions of Galactocentric radius, over the range $6

We present accurate models of the gravitational potential produced by a radially exponential disk mass distribution. The models are produced by combining three separate Miyamoto-Nagai disks. Such models have been used previously to model the disk of the Milky Way, but here we extend this framework to allow its application to disks of any mass, scalelength, and a wide range of thickness from infinitely thin to near spherical (ellipticities from 0 to 0.9). The models have the advantage of simplicity of implementation, and we expect faster run speeds over a double exponential disk treatment. The potentials are fully analytical, and differentiable at all points. The mass distribution of our models deviates from the radial mass distribution of a pure exponential disk by <0.4% out to 4 disk scalelengths, and <1.9% out to 10 disk scalelengths. We tabulate fitting parameters which facilitate construction of exponential disks for any scalelength, and a wide range of disk thickness (a user-friendly, web-based interface is also available). Our recipe is well suited for numerical modelling of the tidal effects of a giant disk galaxy on star clusters or dwarf galaxies. We consider three worked examples; the Milky Way thin and thick disk, and a disky dwarf galaxy., Comment: Accepted to MNRAS February 2nd 2015, 7 pages, 7 figure, 2 tables

Context. Many large stellar surveys have been and are still being carried out, providing huge amounts of data, for which stellar physical parameters will be derived. Solar twins and analogues provide a means to test the calibration of these stellar catalogues because the Sun is the best-studied star and provides precise fundamental parameters. Solar twins should be centred on the solar values. Aims. This spectroscopic study of solar analogues selected from the Geneva-Copenhagen Survey (GCS) at a resolution of 48,000 provides effective temperatures and metallicities for these stars. We test whether our spectroscopic parameters, as well as the previous photometric calibrations, are properly centred on the Sun. In addition, we search for more solar twins in our sample. Methods. The methods used in this work are based on literature methods for solar twin searches and on methods we developed in previous work to distinguish the metallicity-temperature degeneracies in the differential comparison of spectra of solar analogues versus a reference solar reflection spectrum. Results. We derive spectroscopic parameters for 148 solar analogues (about 70 are new entries to the literature) and verify with a-posteriori differential tests that our values are well-centred on the solar values. We use our dataset to assess the two alternative calibrations of the GCS parameters; our methods favour the latest revision. We show that the choice of spectral line list or the choice of asteroid or time of observation does not affect the results. We also identify seven solar twins in our sample, three of which are published here for the first time. Conclusions. Our methods provide an independent means to differentially test the calibration of stellar catalogues around the values of a well-known benchmark star, which makes our work interesting for calibration tests of upcoming Galactic surveys., Comment: 12 pages, 11 figures. Published version

The capability to reconstruct dissolved stellar systems in dynamical and chemical space is a key factor in improving our understanding of the evolution of the Milky Way. Here we concentrate on the dynamical aspect and given that a significant portion of the stars in the Milky Way have been born in stellar associations or clusters that have lived a few Myr up to several Gyr, we further restrict our attention to the evolution of star clusters. We have carried out our simulations in two steps: (1) we create a simulation of dissolution and mixing processes which yields a close fit to the present-day Milky Way dynamics and (2) we have evolved three sets of stellar clusters with masses of 400, 1000 and 15000 M$_{\odot}$ to dissolution. The birth location of these sets was 4, 6, 8 and 10 kpc for the 400 and 1000 M$_{\odot}$ clusters and 4, 6, 8, 10 and 12 kpc for the 15000 M$_{\odot}$. We have focused our efforts on studying the state of the escapers from these clusters after 4.5 Gyr of evolution with particular attention to stars that reach the Solar annulus, i.e. 7.5 $\le$ R$_{\rm{gc}}$ $\le$ 8.5 kpc. We give results for Solar twins and siblings over a wide range of radii and cluster masses for two dissolution mechanisms. From kinematics alone, we conclude that the Sun was $\sim$50 per cent more likely to have been born near its current Galactocentric radius, rather than have migrated (radially) $\sim$2 kpc since birth. We conclude our analysis by calculating magnitudes and colours of our single stars for comparison with the samples that the $Gaia$, $Gaia$-ESO and GALAH-AAO surveys will obtain. In terms of reconstructing dissolved star clusters we find that on short time-scales we cannot rely on kinematic evolution alone and thus it will be necessary to extend our study to include information on chemical space., Comment: Accepted for publication in MNRAS