Source author record

Chris Flynn

Chris Flynn appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

39works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

39 published item(s)

preprint2023arXiv

Search for gravitational wave signals from known pulsars in LIGO-Virgo O3 data using the 5n-vector ensemble method

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.

preprint2022arXiv

Circularly polarized radio emission from the repeating fast radio burst source FRB 20201124A

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.

preprint2022arXiv

Milky Way Mass with K Giants and BHB Stars Using LAMOST, SDSS/SEGUE, and Gaia: 3D Spherical Jeans Equation and Tracer Mass Estimator

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.

preprint2022arXiv

Probing Galactic variations in the fine-structure constant using solar twin stars: methodology and results

The rich absorption spectra of Sun-like stars are enticing probes for variations in the fine-structure constant, $α$, which gauges the strength of electromagnetism. While individual line wavelengths are sensitive to $α$, 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 $α$ to $\approx$50 parts per billion, $\sim$2 orders of magnitude less than other Galactic constraints.

preprint2022arXiv

Probing Galactic variations in the fine-structure constant using solar twin stars: systematic errors

Sun-like stars are a new probe of variations in the fine-structure constant, $α$, 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}$Å at 5000Å). 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 $α$, including at much larger distances.

preprint2022arXiv

Survey for Distant Solar Twins (SDST) -- I. EPIC method for stellar parameter measurement

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.\ $σ\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$.

preprint2022arXiv

The ultra narrow FRB20191107B, and the origins of FRB scattering

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 $μ$s and a scattering tail with an exponentially decaying time-scale of 21.4 $μ$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 $μ$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.

preprint2021arXiv

Characterizing the FRB host galaxy population and its connection to transients in the local and extragalactic Universe

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.

preprint2021arXiv

Fast radio bursts as probes of feedback from active galactic nuclei

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.

preprint2020arXiv

A population analysis of pulse broadening in ASKAP Fast Radio Bursts

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.

preprint2020arXiv

Anisotropy of the Milky Way's stellar halo using K giants from LAMOST and $Gaia$

The anisotropy parameter $β$ 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 $β$ 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 $(σ_r, σ_θ, σ_ϕ)$ (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 $β\approx 0.9$ for $-1.8 \leq$ [Fe/H] $< -1.3$ (for the bulk of the stars) to $β\approx 0.6$ for the lowest metallicities ([Fe/H] $< -1.8$). Considering our sample as a whole, $β\approx0.8$ and, beyond 25 kpc, the orbits gradually become less radial and anisotropy decreases to $β<0.3$ past 100 kpc. Within 8 kpc, $β<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.

preprint2020arXiv

High time resolution and polarisation properties of ASKAP-localised fast radio bursts

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.

preprint2020arXiv

Modelling Double Neutron Stars: Radio and Gravitational Waves

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 $χ_\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.

preprint2020arXiv

The fates of massive stars: exploring uncertainties in stellar evolution with METISSE

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.

preprint2020arXiv

The host galaxies and progenitors of Fast Radio Bursts localized with the Australian Square Kilometre Array Pathfinder

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<z<0.48$. All four are massive galaxies (log( $M_{*}/ M_{\odot}$) $\sim 9.4 -10.4$) with modest star-formation rates of up to $2M_{\odot}$yr$^{-1}$ -- very different to the host galaxy of the first repeating FRB 121102, which is a dwarf galaxy with a high specific star-formation rate. The FRBs localized by ASKAP typically lie in the outskirts of their host galaxies, which appears to rule out FRB progenitor models that invoke active galactic nuclei (AGN) or free-floating cosmic strings. The stellar population seen in these host galaxies also disfavors models in which all FRBs arise from young magnetars produced by superluminous supernovae (SLSNe), as proposed for the progenitor of FRB 121102. A range of other progenitor models (including compact-object mergers and magnetars arising from normal core-collapse supernovae) remain plausible.

preprint2020arXiv

The Rotation Curve, Mass Distribution and Dark Matter Content of the Milky Way from Classical Cepheids

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 $ρ_{\rm DM,\odot} = 0.33\pm0.03$ GeV ${\rm cm^{-3}}$, respectively.

preprint2020arXiv

The UTMOST pulsar timing programme II: Timing noise across the pulsar population

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 ($ν$) and spin-down frequency ($\dotν$), 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 $ν^{a} |\dotν|^{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.

preprint2019arXiv

The low density and magnetization of a massive galaxy halo exposed by a fast radio burst

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.

preprint2016arXiv

The NuSTAR view of the non-thermal emission from PSR J0437-4715

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 $Γ= 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.

preprint2015arXiv

Fading Features Found in the Kinematics of the Far-Reaching Milky Way Stellar Halo

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 $σ$ and anisotropy parameter $β$ of the stars have been determined as functions of Galactocentric radius, over the range $6<R_\mathrm{GC} < 25$ kpc, and show a strong dip in the anisotropy profile at $R_\mathrm{GC}\sim17$ kpc. By directly integrating orbits of particles in a 3-D model of the Galactic potential with these characteristics, we show that the $σ$ and $β$ profiles quickly evolve on a time scale of a $\mathrm{few}\times10$ Myr whereas the density $ρ$ profile remains largely unaffected. We suggest that the feature is therefore transient. The origin of such features in the Galactic halo remains unclear.

preprint2015arXiv

Simple and accurate modelling of the gravitational potential produced by thick and thin exponential disks

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.

preprint2015arXiv

Spectroscopic study of solar twins and analogues

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.

preprint2015arXiv

Tracking Cluster Debris (TraCD)- I. Dissolution of clusters and searching the solar cradle

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.

preprint2014arXiv

Dissecting simulated disc galaxies I: the structure of mono-age populations

We study seven simulated disc galaxies, three with a quiescent merger history, and four with mergers in their last 9 Gyr of evolution. We compare their structure at z=0 by decomposing them into "mono-age populations" (MAPs) of stars within 500 Myr age bins. All studied galaxies undergo a phase of merging activity at high redshift, so that stars older than 9 Gyr are found in a centrally concentrated component, while younger stars are mostly found in discs. We find that most MAPs have simple exponential radial and vertical density profiles, with a scale-height that typically increases with age. Because a large range of merger histories can create populations with simple structures, this suggests that the simplicity of the structure of mono-abundance populations observed in the Milky Way by Bovy et al. (2012b,c) is not necessarily a direct indicator of a quiescent history for the Milky Way. Similarly, the anti-correlation between scale-length and scale-height does not necessarily imply a merger-free history. However, mergers produce discontinuities between thin and thick disc components, and jumps in the age-velocity relation. The absence of a structural discontinuity between thin and thick disc observed in the Milky Way would seem to be a good indicator that no merger with a mass ratio larger than 1:15-1:10 occurred in the last 9 Gyr. Mergers at higher redshift might nevertheless be necessary to produce the thickest, hottest components of the Milky Way's disc.

preprint2014arXiv

Dissecting simulated disc galaxies II: the age-velocity relation

We study the relation between stellar ages and vertical velocity dispersion (the age-velocity relation, or AVR) in a sample of seven simulated disc galaxies. In our simulations, the shape of the AVR for stars younger than 9 Gyr depends strongly on the merger history at low redshift, with even 1:10 - 1:15 mergers being able to create jumps in the AVR (although these jumps might not be detectable if the errors on stellar ages are on the order of 30%). For galaxies with a quiescent history at low redshift, we find that the vertical velocity dispersion rises smoothly for ages up to 8-9 Gyr, following a power law with a slope of ~0.5, similar to what is observed in the solar neighbourhood by the Geneva-Copenhagen Survey. For these galaxies, we show that the slope of the AVR is not imprinted at birth, but is the result of subsequent heating. By contrast, in all our simulations, the oldest stars form a significantly different population, with a high velocity dispersion. These stars are usually born kinematically hot in a turbulent phase of intense mergers at high redshift, and also include some stars accreted from satellites. This maximum in velocity dispersion is strongly decreased when age errors are included, suggesting that observations can easily miss such a jump with the current accuracy of age measurements.

preprint2014arXiv

On the colour of the dark side of the Moon

Aims. 'Earthshine' is the dim light seen on the otherwise dark side of the Moon, particularly when it is close to new. 'Earthlight', or reflected sunlight from the Earth, is the source of Earthshine. Using B and V band CCD images of both the dark and bright side of the Moon, we aim to estimate the Johnson photometry B-V colour of the Earthshine for the first time since the late 1960s. From these measurements we are also able to quantify the colour of Earthlight. Methods. We present images of the Moon taken with a small refractor in Hawaii, in B and V bands and taken under favourable conditions so that scattered light in both bands almost completely cancels, yielding a map of the surface in B-V colour. Co-addition of 100 such images taken in rapid succession substantially improves the signal to noise ratio, and several sources of photometric bias are eliminated by use of relative methods. Results. The earthlit dark side of the Moon is observed to be 0.150 +/- 0.005 mag bluer in B-V than the sunlit bright side, in good agreement with the only known previous measurement of this quantity from 1967. Arguing on the basis of the change in B - V for sunlight reflected once o? the Moon, we derive a colour for earthlight of B-V = 0.44+/-0.02 mag (without applying a small, uncertain, phase-dependent reddening correction). The absence of a colour-gradient in the B-V image implies that the scattering properties of the atmosphere+optical system are almost exactly matched in the two wavelength bands, the consequences of which are discussed.

preprint2014arXiv

Red giants in the outer halo of the elliptical galaxy NGC 5128 / Centaurus A

We used VIMOS on VLT to perform $V$ and $I$ band imaging of the outermost halo of NGC 5128 / Centaurus A ($(m-M)_0=27.91\pm0.08$), 65 kpc from the galaxy's center and along the major axis. The stellar population has been resolved to $I_0 \approx 27$ with a $50\%$ completeness limit of $I_0 = 24.7$, well below the tip of the red-giant branch (TRGB), which is seen at $I_0 \approx 23.9$. The surface density of NGC 5128 halo stars in our fields was sufficiently low that dim, unresolved background galaxies were a major contaminant in the source counts. We isolated a clean sample of red-giant-branch (RGB) stars extending to $\approx 0.8$ mag below the TRGB through conservative magnitude and color cuts, to remove the (predominantly blue) unresolved background galaxies. We derived stellar metallicities from colors of the stars via isochrones and measured the density falloff of the halo as a function of metallicity by combining our observations with HST imaging taken of NGC 5128 halo fields closer to the galaxy center. We found both metal-rich and metal-poor stellar populations and found that the falloff of the two follows the same de Vaucouleurs' law profiles from $\approx 8$ kpc out to $\approx$ 70 kpc. The metallicity distribution function (MDF) and the density falloff agree with the results of two recent studies of similar outermost halo fields in NGC 5128. We found no evidence of a "transition" in the radial profile of the halo, in which the metal-rich halo density would drop rapidly, leaving the underlying metal-poor halo to dominate by default out to greater radial extent, as has been seen in the outer halo of two other large galaxies. If NGC 5128 has such a transition, it must lie at larger galactocentric distances.

preprint2014arXiv

Solar analogues and Solar twins in the HARPS archive

We present 63 Solar analogues and twins for which high S/N archival data are available for the HARPS high resolution spectrograph at the ESO 3.6m telescope. We perform a differential analysis of these stellar spectra relative to the Solar spectrum, similar to previous work using ESO 2.2m/FEROS data, and expand our analysis by introducing a new method to test the temperature and metallicity calibration of Sun-like stars in the Geneva-Copenhagen-Survey (GCS). The HARPS data are significantly better than the FEROS data, with improvements in S/N, spectral resolution, and number of lines we can analyse. We confirm the offsets to the photometric scale found in our FEROS study. We confirm 3 Solar twins found in the FEROS data as Solar twins in the HARPS data, as well as identify 6 new twins.

preprint2012arXiv

New solar twins and the metallicity and temperature scales of the Geneva Copenhagen Survey

We search for "solar twins" in the Geneva-Copenhagen Survey (GCS) using high resolution optical spectroscopy. We initially select Sun-like stars from the GCS by absolute magnitude, b-y colour and metallicity close to the solar values. Our aim is to find the stars which are spectroscopically very close to the Sun using line depth ratios and the median equivalent widths and depths of selected lines with a range of excitation potentials. We present the ten best stars fulfilling combined photometric and spectroscopic criteria, of which six are new twins. We use our full sample of Sun-like stars to examine the calibration of the metallicity and temperature scale in the GCS. Our results give rise to the conclusion that the GCS may be offset from the solar temperature and metallicity for sun-like stars by 100K and 0.1dex, respectively.

preprint2011arXiv

On the vertical equilibrium of the local Galactic disk and the search for disk dark matter

Estimates of the dynamical surface mass density at the solar Galactocentric distance are commonly derived assuming that the disk is in vertical equilibrium with the Galactic potential. This assumption has recently been called into question, based on the claim that the ratio between the kinetic and the gravitational energy in such solutions is a factor of 3 larger than required if Virial equilibrium is to hold. Here we show that this ratio between energies was overerestimated and that the disk solutions are likely to be in Virial equilibrium after all. We additionally demonstrate, using one-dimensional numerical simulations, that the disks are indeed in equilibrium. Hence, given the uncertainties, we find no reason to cast doubt on the steady-state solutions which are traditionally used to measure the matter density of the disk.

preprint2010arXiv

Probing the Galaxy's bars via the Hercules stream

It has been suggested that a resonance between a rotating bar and stars in the solar neighbourhood can produce the so called 'Hercules stream'. Recently, a second bar may have been identified in the Galactic centre, the so called 'long bar', which is longer and much flatter than the traditional Galactic bar, and has a similar mass. We looked at the dynamical effects of both bars, separately and together, on orbits of stars integrated backwards from local position and velocities, and a model of the Galactic potential which includes the bars directly. Both bars can produce Hercules like features, and allow us to measure the rotation rate of the bar(s). We measure a pattern speed, for both bars, of 1.87 +/- 0.02 times the local circular frequency. This is on par with previous measurements for the Galactic bar, although we do adopt a slightly different Solar motion. Finally, we identify a new kinematic feature in local velocity space, caused by the long bar, which is tempting to identify with the high velocity 'Arcturus' stream.

preprint2010arXiv

Revisiting Delta Y/Delta Z from multiple main sequences in Globular Clusters: insight from nearby stars

For nearby K dwarfs, the broadening of the observed Main Sequence at low metallicities is much narrower than expected from isochrones with the standard helium-to-metal enrichment ratio DY/DZ=2. Though the latter value fits well the Main Sequence around solar metallicity, and agrees with independent measurements from HII regions as well as with theoretical stellar yields and chemical evolution models, a much higher DY/DZ~10 is necessary to reproduce the broadening observed for nearby subdwarfs. This result resembles, on a milder scale, the very high DY/DZ estimated from the multiple Main Sequences in Omega Cen and NGC 2808. Although not "inverted" as in omega Cen, where the metal-rich Main Sequence is bluer than the metal-poor one, the broadening observed for nearby subdwarfs is much narrower than stellar models predict for a standard helium content. We use this empirical evidence to argue that a revision of lower Main Sequence stellar models, suggested from nearby stars, could significantly reduce the helium content inferred for the subpopulations of those globular clusters. A simple formula based on empirically calibrated homology relations is constructed, for an alternative estimate of DY/DZ in multiple main sequences. We find that, under the most favourable assumptions, the estimated helium content for the enriched populations could decrease from Y~0.4 to as low as Y~0.3.

preprint2010arXiv

The Inner Halo of M87: A First Direct View of the Red-Giant Population

An unusually deep (V,I) imaging dataset for the Virgo supergiant M87 with the Hubble Space Telescope ACS successfully resolves its brightest red-giant stars, reaching M_I(lim) = -2.5. After assessing the photometric completeness and biasses, we use this material to estimate the metallicity distribution for the inner halo of M87, finding that the distribution is very broad and likely to peak near [m/H] ~ -0.4 and perhaps higher. The shape of the MDF strongly resembles that of the inner halo for the nearby giant E galaxy NGC 5128. As a byproduct of our study, we also obtain a preliminary measurement of the distance to M87 with the TRGB (red-giant branch tip) method; the result is (m-M)_0 = 31.12 +- 0.14 (d = 16.7 +- 0.9 Mpc). Averaging this result with three other recent techniques give a weighted mean d(M87) = (16.4 +- 0.5) Mpc.

preprint2008arXiv

Dynamical effects of the long bar in the Milky Way

We examine the dynamical effects on disk stars of a "long bar" in the Milky Way by inserting a triaxial rotating bar into an axisymmetric disk+bulge+dark halo potential and integrating 3-D orbits of 104 tracer stars over a period of 2 Gyr. The long bar has been detected via "clump giants" in the IR by Lopez-Corredoira et al. (2007), and is estimated to have semi-major axes of (3.9 : 0.6 : 0.1) kpc and a mass of 6 10^9 Msun. We find such a structure has a slight impact on the inner disk-system, moving tracers near to the bar into the bar-region, as well as into the bulge. These effects are under continuing study.

preprint2005arXiv

Kinematics of the Galactic Halo from Horizontal Branch stars in the Hamburg/ESO Survey

Large samples of Field Horizontal Branch (FHB) stars make excellent tracers of the Galactic halo; by studying their kinematics, one can infer important physical properties of our Galaxy. Here we present the results of a medium-resolution spectroscopic survey of 530 FHB stars selected from the Hamburg/ESO survey. The stars have a mean distance of ~7 kpc and thus probe the inner parts of the Milky Way halo. We measure radial velocities from the spectra in order to test the model of Sommer-Larsen et al., who suggested that the velocity ellipsoid of the halo changes from radially-dominated orbits to tangentially-dominated orbits as one proceeds from the inner to the outer halo. We find that the present data are unable to discriminate between this model and a more simple isothermal ellipsoid; we suggest that additional observations towards the Galactic centre might help to differentiate them.

preprint2003arXiv

M Dwarfs from Hubble Space Telescope Star Counts. V. The I-band Luminosity Function

We derive the disk I-band luminosity function from the Zheng et al. sample of \~1400 disk M dwarfs observed with the Hubble Space Telescope. We adopt a Galactic-height-dependent color-magnitude relation to account for the metallicity gradient above the Galactic plane. The resultant I-band luminosity function peaks at M_I~9.5 and drops sharply toward M_I~10.5.

preprint2002arXiv

Cosmology on a Mesh

An adaptive multi grid approach to simulating the formation of structure from collisionless dark matter is described. MLAPM (Multi-Level Adaptive Particle Mesh) is one of the most efficient serial codes available on the cosmological 'market' today. As part of Swinburne University's role in the development of the Square Kilometer Array, we are implementing hydrodynamics, feedback, and radiative transfer within the MLAPM adaptive mesh, in order to simulate baryonic processes relevant to the interstellar and intergalactic media at high redshift. We will outline our progress to date in applying the existing MLAPM to a study of the decay of satellite galaxies within massive host potentials.

preprint1996arXiv

Kinematics of the Outer Stellar Halo

We have tested whether the simple model for the kinematics of the Galactic stellar halo (in particular the outer halo) proposed by Sommer-Larsen, Flynn and Christensen (SLFC) is physically realizable, by directly integrating particles in a 3-D model of the Galactic potential. We are able to show that the SLFC solution can be realized in terms of a distribution of particles with stationary statistical properties in phase-space. Hence, the SLFC model, which shows a notable change in the anisotropy from markedly radial at the sun to markedly tangential beyond about Galactocentric radius r=20 kpc, seems a tenable description of outer halo kinematics.

preprint1996arXiv

Metallicities and kinematics of G and K dwarfs

We have developed a photometric abundance indicator for G and K dwarfs, based on accurate, spectroscopically determined abundances for G and K dwarfs of Morell, Källander and Butcher (1992) and Morell (1994). Broadband Cousins R-I photometry is used to estimate effective temperature and the Geneva b_1 colour to estimate line blanketing in the blue and hence abundance. Abundances can be derived in the range -2.0<[Fe/H]<0.5 for G0 to K3 dwarfs, with a scatter in [Fe/H] of 0.2 dex. We apply the method to a sample of Gliese catalog G and K dwarfs, and examine the metallicity and kinematic properties of the stars. The stars show the well established observational features of the disk, thick disk and halo in the solar neighbourhood. We find that the distribution of local K dwarf metallicity is quite similar to local G dwarfs, indicating that there is a ``K-dwarf'' as well as a G-dwarf problem.