Search

With Hubble Space Telescope Fine Guidance Sensor astrometry and published and previously unpublished radial velocity measures we explore the exoplanetary system mu Arae. Our modeling of the radial velocities results in improved orbital elements for the four previously known components. Our astrometry contains no evidence for any known companion, but provides upper limits for three companion masses. A final summary of all past Fine Guidance Sensor exoplanet astrometry results uncover a bias towards small inclinations (more face-on than edge-on). This bias remains unexplained by either small number statistics, modeling technique, Fine Guidance Sensor mechanical issues, or orbit modeling of noise-dominated data. A numerical analysis using our refined orbital elements suggests that planet d renders the mu Arae system dynamically unstable on a timescale of 10^5 years, in broad agreement with previous work., Comment: Accepted to Astronomical Journal

The Hilda asteroids are among the least studied populations in the asteroid belt, despite their potential importance as markers of Jupiter's migration in the early Solar system. We present new mid-infrared observations of two notable Hildas, (1162) Larissa and (1911) Schubart, obtained using the Faint Object infraRed CAmera for the SOFIA Telescope (FORCAST), and use these to characterise their thermal inertia and physical properties. For (1162) Larissa, we obtain an effective diameter of \textcolor{black}{46.5$^{+2.3}_{-1.7}$~km, an albedo of 0.12~$\pm$~0.02, and a thermal inertia of 15$^{+10}_{-8}$ Jm$^{-2}$s$^{1/2}$K$^{-1}$. In addition, our Larissa thermal measurements are well matched with an ellipsoidal shape with an axis ratio a/b=1.2 for the most-likely spin properties. Our modelling of (1911) Schubart is not as refined, but the thermal data point towards a high-obliquity spin-pole, with a best-fit a/b=1.3 ellipsoidal shape. This spin-shape solution is yielding a diameter of 72$^{+3}_{-4}$ km, an albedo of 0.039$\pm$~0.02, and a thermal inertia below 30 Jm$^{-2}$s$^{1/2}$K$^{-1}$ (or 10$^{+20}_{-5}$Jm$^{-2}$s$^{1/2}$K$^{-1}$).} As with (1162) Larissa, our results suggest that (1911) Schubart is aspherical, and likely elongated in shape. Detailed dynamical simulations of the two Hildas reveal that both exhibit strong dynamical stability, behaviour that suggests that they are primordial, rather than captured objects. The differences in their albedos, along with their divergent taxonomical classification, suggests that despite their common origin, the two have experienced markedly different histories., Comment: Accepted to appear in MNRAS; 13 pages, 5 figures, 6 tables

The orbital solutions of published multi-planet systems are not necessarily dynamically stable on timescales comparable to the lifetime of the system as a whole. For this reason, dynamical tests of the architectures of proposed exoplanetary systems are a critical tool to probe the stability and feasibility of the candidate planetary systems, with the potential to point the way towards refined orbital parameters of those planets. Such studies can even help in the identification of additional companions in such systems. Here we examine the dynamical stability of three planetary systems, orbiting HD 67087, HD 110014, and HD 133131A. We use the published radial velocity measurements of the target stars to determine the best-fit orbital solutions for these planetary systems using the Systemic console. We then employ the n-body integrator Mercury to test the stability of a range of orbital solutions lying within 3-$\sigma$ of the nominal best-fit for a duration of 100 Myr. From the results of the n-body integrations, we infer the best-fit orbital parameters using the Bayesian package Astroemperor. We find that both HD 110014 and HD 133131A have long-term stable architectures that lie within the 1-$\sigma$ uncertainties of the nominal best-fit to their previously determined orbital solutions. However, the HD 67087 system exhibits a strong tendency toward instability on short timescales. We compare these results to the predictions made from consideration of the angular momentum deficit criterion, and find that its predictions are consistent with our findings., Comment: 9 pages, 5 figures, 3 tables. Accepted for publication in MNRAS

In order to accurately determine stellar properties, knowledge of the effective temperature of stars is vital. We implement Gaia and 2MASS photometry in the InfraRed Flux Method and apply it to over 360,000 stars across different evolutionary stages in the GALAH DR3 survey. We derive colour-effective temperature relations that take into account the effect of metallicity and surface gravity over the range 4000 to 8000 kelvin, from very metal-poor stars to super solar metallicities. The internal uncertainty of these calibrations is of order 40-80 kelvin depending on the colour combination used. Comparison against solar-twins, Gaia benchmark stars and the latest interferometric measurements validates the precision and accuracy of these calibrations from F to early M spectral types. We assess the impact of various sources of uncertainties, including the assumed extinction law, and provide guidelines to use our relations. Robust solar colours are also derived., Comment: MNRAS accepted. Include Gaia solar colours. Colour-effective temperature routines for Gaia DR2 and DR3 system available at https://github.com/casaluca/colte

Massive sets of stellar spectroscopic observations are rapidly becoming available and these can be used to determine the chemical composition and evolution of the Galaxy with unprecedented precision. One of the major challenges in this endeavour involves constructing realistic models of stellar spectra with which to reliably determine stellar abundances. At present, large stellar surveys commonly use simplified models that assume that the stellar atmospheres are approximately in local thermodynamic equilibrium (LTE). To test and ultimately relax this assumption, we have performed non-LTE calculations for $13$ different elements (H, Li, C, N, O, Na, Mg, Al, Si, K, Ca, Mn, and Ba), using recent model atoms that have physically-motivated descriptions for the inelastic collisions with neutral hydrogen, across a grid of $3756$ 1D MARCS model atmospheres that spans $3000\leq T_{\mathrm{eff}}/\mathrm{K}\leq8000$, $-0.5\leq\log{g/\mathrm{cm\,s^{-2}}}\leq5.5$, and $-5\leq\mathrm{[Fe/H]}\leq1$. We present the grids of departure coefficients that have been implemented into the GALAH DR3 analysis pipeline in order to complement the extant non-LTE grid for iron. We also present a detailed line-by-line re-analysis of $50126$ stars from GALAH DR3. We found that relaxing LTE can change the abundances by between $-0.7\,\mathrm{dex}$ and $+0.2\,\mathrm{dex}$ for different lines and stars. Taking departures from LTE into account can reduce the dispersion in the $\mathrm{[A/Fe]}$ versus $\mathrm{[Fe/H]}$ plane by up to $0.1\,\mathrm{dex}$, and it can remove spurious differences between the dwarfs and giants by up to $0.2\,\mathrm{dex}$. The resulting abundance slopes can thus be qualitatively different in non-LTE, possibly with important implications for the chemical evolution of our Galaxy., Comment: 19 pages, 25 figures, 2 tables, arXiv abstract abridged; accepted for publication in A&A

Asteroids and comets (planetesimals) are created in gas- and dust-rich protoplanetary discs. The presence of these planetesimals around main-sequence stars is usually inferred from the detection of excess continuum emission at infrared wavelengths from dust grains produced by destructive processes within these discs. Modelling of the disc structure and dust grain properties for those discs is often hindered by the absence of any meaningful constraint on the location and spatial extent of the disc. Multi-wavelength, spatially resolved imaging is thus invaluable in refining the interpretation of these systems. Observations of HD 48682 at far-infrared (Spitzer,Herschel) and sub-millimetre (JCMT, SMA) wavelengths indicated the presence of an extended, cold debris disc with a blackbody temperature of 57.9 +\- 0.7 K. Here, we combined these data to perform a comprehensive study of the disc architecture and its implications for the dust grain properties. The deconvolved images revealed a cold debris belt, verified by combining a 3D radiative transfer dust continuum model with image analysis to replicate thestructure using a single, axisymmetric annulus. A Markov chain Monte Carlo analysis calculated the maximum likelihood of HD 48682's disc radius (Rdisc = 89 +17 -20 au), fractional width(DeltaRdisc = 0.41 +0.27 -0.20), position angle (theta = 66.3 +4.5 -4.9 degrees), and inclination (phi = 112.5 +4.2 -4.2 degrees). HD 48682 has been revealed to host a collisionally active, broad disc whose emission is dominated by small dust grains, smin approx. 0.6 microns, and a size distribution exponent of 3.60 +\- 0.02., Comment: 13 Pages, 6 figures, to be published in MNRAS

Over the past three decades, we have witnessed one of the great revolutions in our understanding of the cosmos - the dawn of the Exoplanet Era. Where once we knew of just one planetary system (the Solar system), we now know of thousands, with new systems being announced on a weekly basis. Of the thousands of planetary systems we have found to date, however, there is only one that we can study up-close and personal - the Solar system. In this review, we describe our current understanding of the Solar system for the exoplanetary science community - with a focus on the processes thought to have shaped the system we see today. In section one, we introduce the Solar system as a single well studied example of the many planetary systems now observed. In section two, we describe the Solar system's small body populations as we know them today - from the two hundred and five known planetary satellites to the various populations of small bodies that serve as a reminder of the system's formation and early evolution. In section three, we consider our current knowledge of the Solar system's planets, as physical bodies. In section four, we discuss the research that has been carried out into the Solar system's formation and evolution, with a focus on the information gleaned as a result of detailed studies of the system's small body populations. In section five, we discuss our current knowledge of planetary systems beyond our own - both in terms of the planets they host, and in terms of the debris that we observe orbiting their host stars. As we learn ever more about the diversity and ubiquity of other planetary systems, our Solar system will remain the key touchstone that facilitates our understanding and modelling of those newly found systems, and we finish section five with a discussion of the future surveys that will further expand that knowledge., Comment: Invited Review, Accepted for publication in Publications of the Astronomical Society of the Pacific; 23 figures, plus a further 18 in the appendix; 4 tables

The world's meteorite collections contain a very rich picture of what the early Solar System would have been made of, however the lack of spatial context with respect to their parent population for these samples is an issue. The asteroid population is equally as rich in surface mineralogies, and mapping these two populations (meteorites and asteroids) together is a major challenge for planetary science. Directly probing asteroids achieves this at a high cost. Observing meteorite falls and calculating their pre-atmospheric orbit on the other hand, is a cheaper way to approach the problem. The Global Fireball Observatory (GFO) collaboration was established in 2017 and brings together multiple institutions (from Australia, USA, Canada, Morocco, Saudi Arabia, the UK, and Argentina) to maximise the area for fireball observation time and therefore meteorite recoveries. The members have a choice to operate independently, but they can also choose to work in a fully collaborative manner with other GFO partners. This efficient approach leverages the experience gained from the Desert Fireball Network (DFN) pathfinder project in Australia. The state-of-the art technology (DFN camera systems and data reduction) and experience of the support teams is shared between all partners, freeing up time for science investigations and meteorite searching. With all networks combined together, the GFO collaboration already covers 0.6% of the Earth's surface for meteorite recovery as of mid-2019, and aims to reach 2% in the early 2020s. We estimate that after 5 years of operation, the GFO will have observed a fireball from virtually every meteorite type. This combined effort will bring new, fresh, extra-terrestrial material to the labs, yielding new insights about the formation of the Solar System., Comment: Accepted in PSS. 19 pages, 9 figures

The MINiature Exoplanet Radial Velocity Array (MINERVA) is a dedicated observatory of four 0.7m robotic telescopes fiber-fed to a KiwiSpec spectrograph. The MINERVA mission is to discover super-Earths in the habitable zones of nearby stars. This can be accomplished with MINERVA's unique combination of high precision and high cadence over long time periods. In this work, we detail changes to the MINERVA facility that have occurred since our previous paper. We then describe MINERVA's robotic control software, the process by which we perform 1D spectral extraction, and our forward modeling Doppler pipeline. In the process of improving our forward modeling procedure, we found that our spectrograph's intrinsic instrumental profile is stable for at least nine months. Because of that, we characterized our instrumental profile with a time-independent, cubic spline function based on the profile in the cross dispersion direction, with which we achieved a radial velocity precision similar to using a conventional "sum-of-Gaussians" instrumental profile: 1.8 m s$^{-1}$ over 1.5 months on the RV standard star HD 122064. Therefore, we conclude that the instrumental profile need not be perfectly accurate as long as it is stable. In addition, we observed 51 Peg and our results are consistent with the literature, confirming our spectrograph and Doppler pipeline are producing accurate and precise radial velocities., Comment: 22 pages, 9 figures, submitted to PASP, Peer-Reviewed and Accepted

A group of trans-Neptunian objects (TNO) are dynamically related to the dwarf planet 136108 Haumea. Ten of them show strong indications of water ice on their surfaces, are assumed to have resulted from a collision, and are accepted as the only known TNO collisional family. Nineteen other dynamically similar objects lack water ice absorptions and are hypothesized to be dynamical interlopers. We have made observations to determine sizes and geometric albedos of six of the accepted Haumea family members and one dynamical interloper. Ten other dynamical interlopers have been measured by previous works. We compare the individual and statistical properties of the family members and interlopers, examining the size and albedo distributions of both groups. We also examine implications for the total mass of the family and their ejection velocities. We use far-infrared space-based telescopes to observe the target TNOs near their thermal peak and combine these data with optical magnitudes to derive sizes and albedos using radiometric techniques. We determine the power-law slope of ejection velocity as a function of effective diameter. The detected Haumea family members have a diversity of geometric albedos $\sim$ 0.3-0.8, which are higher than geometric albedos of dynamically similar objects without water ice. The median geometric albedo for accepted family members is $p_V=0.48_{-0.18}^{+0.28}$, compared to 0.08$_{-0.05}^{+0.07}$ for the dynamical interlopers. In the size range $D=175-300$ km, the slope of the cumulative size distribution is $q$=3.2$_{-0.4}^{+0.7}$ for accepted family members, steeper than the $q$=2.0$\pm$0.6 slope for the dynamical interlopers with D$< $500 km. The total mass of Haumea's moons and family members is 2.4% of Haumea's mass. The ejection velocities required to emplace them on their current orbits show a dependence on diameter, with a power-law slope of 0.21-0.50., Comment: 15 pages, 4 figures

Circumbinary disks are common around post-asymptotic giant branch (post-AGB) stars with a stellar companion on orbital time scales of a few 100 to few 1000 days. The presence of a disk is usually inferred from the system's spectral energy distribution (SED), and confirmed, for a sub-sample, by interferometric observations. We used the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument on the Very Large Telescope to obtain extreme adaptive optics assisted scattered light images of the post-AGB binary system AR Puppis. Data have been obtained in the V, I, and H bands. Our observations have produced the first resolved images of AR Puppis' circumbinary disk and confirm its edge-on orientation. In our high angular-resolution and high dynamic-range images we identify several structural components such as a dark mid-plane, the disk surface, and arc-like features. We discuss the nature of these components and use complementary photometric monitoring to relate them to the orbital phase of the binary system. Because the star is completely obscured by the disk at visible wavelengths, we conclude that the long-term photometric variability of the system must be caused by variable scattering, not extinction, of star light by the disk over the binary orbit. Finally, we discuss how the short disk life times and fast evolution of the host stars compared to the ages at which protoplanetary disks are typically observed make systems like AR Puppis valuable extreme laboratories to study circumstellar disk evolution and constrain the time scale of dust grain growth during the planet formation process., Comment: 17 pages, 9 figures, accepted for publication in AJ

Multi-planet systems around evolved stars are of interest to trace the evolution of planetary systems into the post-main sequence phase. HD 47366, an evolved intermediate mass star, hosts two giant planets on moderately eccentric orbits. Previous analysis of the planetary system has revealed that it is dynamically unstable on timescales much shorter than the stellar age unless the planets are trapped in mutual 2:1 mean motion resonance, inconsistent with the orbital solution presented in \cite{2016Sato} (hereafter: S16), or are moving on mutually retrograde orbits. Here we examine the orbital stability of the system presented in S16 using the $n$-body code {\sc Mercury} over a broad range of $a$--$e$ parameter space consistent with the observed radial velocities, assuming they are on co-planar orbits. Our analysis confirms that the system as proposed in S16 is not dynamically stable. We therefore undertake a thorough re-analysis of the available observational data for the HD 47366 system, through the Levenberg-Marquardt technique and confirmed by MCMC Bayesian methodology. Our re-analysis reveals an alternative, lower eccentricity fit that is vastly preferred over the highly eccentric orbital solution obtained from the nominal best-fit presented in S16. The new, improved dynamical simulation solution reveals the reduced eccentricity of the planetary orbits, shifting the HD 47366 system into the edge of a broad stability region, increasing our confidence that the planets are all that they seem to be. Our rigorous examination of the dynamical stability of HD 47366 stands as a cautionary tale in finding the global best-fit model., Comment: 12 pages, 12 figures, 5 tables, accepted for publication in AJ

Identifying and characterizing young stars in the Solar neighbourhood is essential to find and describe planets in the early stages of their evolution. This work seeks to identify nearby young stars showing a Lithium 6707.78$\,\AA$ absorption line in the GALAH survey. A robust, data-driven approach is used to search for corresponding templates in the pool of 434,215 measured dwarf spectra in the survey. It enables a model-free search for best-matching spectral templates for all stars, including M dwarfs with strong molecular absorption bands. 3147 stars have been found to have measurable Lithium: 1408 G and 892 K0-K5 dwarfs (EW(Li)$>$0.1$\,$\AA), 335 K5-K9 ($>$0.07$\,$\AA) and 512 M0$\sim$M4 dwarfs ($>$0.05$\,$\AA). Stars with such Lithium features are used to investigate the possibility of searching for young stars above the main sequence based merely on their parallaxes and broad-band photometry. Selection of young stars above the main sequence is highly effective for M dwarfs, moderately effective for K dwarfs and ineffective for G dwarfs. Using a combination of the Lithium information and the complete 6D kinematics from Gaia and GALAH, 305 new candidate moving group members have been found, 123 of which belong to the Scorpius-Centaurus association, 36 to the Pleiades and 25 to the Hyades clusters., Comment: Submitted to MNRAS. 12 pages, 8 figures, 3 tables

The Galactic Archaeology with HERMES (GALAH) survey is a large-scale stellar spectroscopic survey of the Milky Way and designed to deliver chemical information complementary to a large number of stars covered by the $Gaia$ mission. We present the GALAH second public data release (GALAH DR2) containing 342,682 stars. For these stars, the GALAH collaboration provides stellar parameters and abundances for up to 23 elements to the community. Here we present the target selection, observation, data reduction and detailed explanation of how the spectra were analysed to estimate stellar parameters and element abundances. For the stellar analysis, we have used a multi-step approach. We use the physics-driven spectrum synthesis of Spectroscopy Made Easy (SME) to derive stellar labels ($T_\mathrm{eff}$, $\log g$, $\mathrm{[Fe/H]}$, $\mathrm{[X/Fe]}$, $v_\mathrm{mic}$, $v \sin i$, $A_{K_S}$) for a representative training set of stars. This information is then propagated to the whole survey with the data-driven method of $The~Cannon$. Special care has been exercised in the spectral synthesis to only consider spectral lines that have reliable atomic input data and are little affected by blending lines. Departures from local thermodynamic equilibrium (LTE) are considered for several key elements, including Li, O, Na, Mg, Al, Si, and Fe, using 1D MARCS stellar atmosphere models. Validation tests including repeat observations, Gaia benchmark stars, open and globular clusters, and K2 asteroseismic targets lend confidence in our methods and results. Combining the GALAH DR2 catalogue with the kinematic information from $Gaia$ will enable a wide range of Galactic Archaeology studies, with unprecedented detail, dimensionality, and scope., Comment: 34+7 pages, 31 Figures, 5 tables, 1 catalog, submitted to MNRAS

Using data from the GALAH pilot survey, we determine properties of the Galactic thin and thick disks near the solar neighbourhood. The data cover a small range of Galactocentric radius ($7.9 \leq R_\mathrm{GC} \leq 9.5$ kpc), but extend up to 4 kpc in height from the Galactic plane, and several kpc in the direction of Galactic anti-rotation (at longitude $260 ^\circ \leq \ell \leq 280^\circ$). This allows us to reliably measure the vertical density and abundance profiles of the chemically and kinematically defined `thick' and `thin' disks of the Galaxy. The thin disk (low-$\alpha$ population) exhibits a steep negative vertical metallicity gradient, at d[M/H]/d$z=-0.18 \pm 0.01$ dex kpc$^{-1}$, which is broadly consistent with previous studies. In contrast, its vertical $\alpha$-abundance profile is almost flat, with a gradient of d[$\alpha$/M]/d$z$ = $0.008 \pm 0.002$ dex kpc$^{-1}$. The steep vertical metallicity gradient of the low-$\alpha$ population is in agreement with models where radial migration has a major role in the evolution of the thin disk. The thick disk (high-$\alpha$ population) has a weaker vertical metallicity gradient d[M/H]/d$z = -0.058 \pm 0.003$ dex kpc$^{-1}$. The $\alpha$-abundance of the thick disk is nearly constant with height, d[$\alpha$/M]/d$z$ = $0.007 \pm 0.002$ dex kpc$^{-1}$. The negative gradient in metallicity and the small gradient in [$\alpha$/M] indicate that the high-$\alpha$ population experienced a settling phase, but also formed prior to the onset of major SNIa enrichment. We explore the implications of the distinct $\alpha$-enrichments and narrow [$\alpha$/M] range of the sub-populations in the context of thick disk formation., Comment: 18 pages, 15 figures, 3 tables. Accepted for publication in MNRAS

Debris discs are the dusty aftermath of planet formation processes around main-sequence stars. Analysis of these discs is often hampered by the absence of any meaningful constraint on the location and spatial extent of the disc around its host star. Multi-wavelength, resolved imaging ameliorates the degeneracies inherent in the modelling process, making such data indispensable in the interpretation of these systems. The Herschel Space Observatory observed HD 105211 (Eta Cru, HIP 59072) with its PACS instrument in three far-infrared wavebands (70, 100 and 160 um). Here we combine these data with ancillary photometry spanning optical to far-infrared wavelengths in order to determine the extent of the circumstellar disc. The spectral energy distribution and multi-wavelength resolved emission of the disc are simultaneously modelled using a radiative transfer and imaging codes. Analysis of the Herschel/PACS images reveals the presence of extended structure in all three PACS images. From a radiative transfer model we derive a disc extent of 87.0 +/- 2.5 au, with an inclination of 70.7 +/- 2.2 degrees to the line of sight and a position angle of 30.1 +/- 0.5 degrees. Deconvolution of the Herschel images reveal a potential asymmetry but this remains uncertain as a combined radiative transfer and image analysis replicate both the structure and the emission of the disc using a single axisymmetric annulus., Comment: 11 pages, 4 figures, 6 tables, Accepted for publication in MNRAS on 25th March, 2017

Galah is an ongoing high-resolution spectroscopic survey with the goal of disentangling the formation history of the Milky Way, using the fossil remnants of disrupted star formation sites which are now dispersed around the Galaxy. It is targeting a randomly selected, magnitude limited ($V \leq 14$) sample of stars, with the goal of observing one million objects. To date, 300,000 spectra have been obtained. Not all of them are correctly processed by parameter estimation pipelines and we need to know about them. We present a semi-automated classification scheme which identifies different types of peculiar spectral morphologies, in an effort to discover and flag potentially problematic spectra and thus help to preserve the integrity of the survey's results. To this end we employ a recently developed dimensionality reduction technique t-SNE (t-distributed Stochastic Neighbour Embedding), which enables us to represent the complex spectral morphology in a two-dimensional projection map while still preserving the properties of the local neighbourhoods of spectra. We find that the majority (178,483) of the 209,533 Galah spectra considered in this study represents normal single stars, whereas 31,050 peculiar and problematic spectra with very diverse spectral features pertaining to 28,579 stars are distributed into 10 classification categories: Hot stars, Cool metal-poor giants, Molecular absorption bands, Binary stars, H$\alpha$/H$\beta$ emission, H$\alpha$/H$\beta$ emission superimposed on absorption, H$\alpha$/H$\beta$ P-Cygni, H$\alpha$/H$\beta$ inverted P-Cygni, Lithium absorption, and Problematic. Classified spectra with supplementary information are presented in the catalogue, indicating candidates for follow-up observations and population studies of the short-lived phases of stellar evolution.

We present a spectropolarimetric snapshot survey of solar-type planet hosting stars. In addition to 14 planet-hosting stars observed as part of the BCool magnetic snapshot survey, we obtained magnetic observations of a further 19 planet-hosting solar-type stars in order to see if the presence of close-in planets had an effect on the measured surface magnetic field (|B$_{\ell}$|). Our results indicate that the magnetic activity of this sample is congruent with that of the overall BCool sample. The effects of the planetary systems on the magnetic activity of the parent star, if any, are too subtle to detect compared to the intrinsic dispersion and correlations with rotation, age and stellar activity proxies in our sample. Four of the 19 newly observed stars, two of which are subgiants, have unambiguously detected magnetic fields and are future targets for Zeeman Doppler Mapping., Comment: 16 pages, 13 figures, Accepted for publication in MNRAS

In many ways, the HR 8799 planetary system strongly resembles our own. It features four giant planets and two debris belts, analogues to the Asteroid and Edgeworth-Kuiper belts. Here, we present the results of dynamical simulations of HR 8799's inner debris belt, to study its structure and collisional environment. Our results suggest that HR 8799's inner belt is highly structured, with gaps between regions of dynamical stability. The belt is likely constrained between sharp inner and outer edges, located at ~6 and ~8 au, respectively. Its inner edge coincides with a broad gap cleared by the 4:1 mean-motion resonance with HR 8799e. Within the belt, planetesimals are undergoing a process of collisional attrition like that observed in the Asteroid belt. However, whilst the mean collision velocity in the Asteroid belt exceeds 5 km/s, the majority of collisions within HR 8799's inner belt occur with velocities of order 1.2 km/s, or less. Despite this, they remain sufficiently energetic to be destructive - giving a source for the warm dust detected in the system. Interior to the inner belt, test particles remain dynamically unstirred, aside from narrow bands excited by distant high-order resonances with HR 8799e. This lack of stirring is consistent with earlier thermal modelling of HR 8799's infrared excess, which predicted little dust inside 6 au. The inner system is sufficiently stable and unstirred that the formation of telluric planets is feasible, although such planets would doubtless be subject to a punitive impact regime, given the intense collisional grinding required in the inner belt to generate the observed infrared excess., Comment: 25 pages, 11 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society on 1st August 2016

In the coming decades, the discovery of the first truly Earth-like exoplanets is anticipated. The characterisation of those planets will play a vital role in determining which are chosen as targets for the search for life beyond the Solar system. One of the many variables that will be considered in that characterisation and selection process is the nature of the potential climatic variability of the exoEarths in question. In our own Solar system, the Earth's long-term climate is driven by several factors - including the modifying influence of life on our atmosphere, and the temporal evolution of Solar luminosity. The gravitational influence of the other planets in our Solar system add an extra complication - driving the Milankovitch cycles that are thought to have caused the on-going series of glacial and interglacial periods that have dominated Earth's climate for the past few million years. Here, we present the results of a large suite of dynamical simulations that investigate the influence of the giant planet Jupiter on the Earth's Milankovitch cycles. If Jupiter was located on a different orbit, we find that the long-term variability of Earth's orbit would be significantly different. Our results illustrate how small differences in the architecture of planetary systems can result in marked changes in the potential habitability of the planets therein, and are an important first step in developing a means to characterise the nature of climate variability on planets beyond our Solar system., Comment: Published in the peer-reviewed proceedings of the 14th Australian Space Research Conference, with ISBN: 13: 978-0-9775740-8-7; Editors: Wayne Short and Iver Cairns. Full proceedings of that conference available online at: http://www.nssa.com.au/14asrc/14ASRC-proceedings.zip Paper contains 14 page, 5 figures

We present estimates of the basic physical properties -- including size and albedo -- of the extreme Centaur 2013 AZ60. These properties have been derived from optical and thermal infrared measurements. Our optical measurements revealed a likely full period of ~9.4 h with a shallow amplitude of 4.5%. By combining optical brightness information and thermal emission data, we are able to derive a diameter of 62.3 +/- 5.3 km and a geometric albedo of 2.9% -- corresponding to an extremely dark surface. Additionally, our finding of ~> 50 Jm^{-2}K^{-1}s^{-1/2} for the thermal inertia is also noticeably for objects in such a distance. The results of dynamical simulations yield an unstable orbit, with a 50% probability that the target will be ejected from the Solar System within 700,000 years. The current orbit of this object as well as its instability could imply a pristine cometary surface. This possibility is in agreement with the observed low geometric albedo and red photometric colour indices for the object, which are a good match for the surface of a dormant comet -- as would be expected for a long-period cometary body approaching perihelion. Despite the fact it was approaching ever closer to the Sun, however, the object exhibited star-like profiles in each of our observations, lacking any sign of cometary activity. By the albedo, 2013 AZ60 is a candidate for the darkest body among the known TNOs., Comment: Accepted for publication in A&A

In many ways, the HR8799 system resembles our Solar system more closely than any other discovered to date - albeit on a larger, younger, and more dramatic scale - featuring four giant planets and two debris belts. The first belt lies beyond the orbit of the outer planet, and mirrors our Solar system's Edgeworth-Kuiper belt. The second belt lies interior to the orbit of the inner planet, HR8799e, and is analogous to our Asteroid Belt. With such a similar architecture, the system is a valuable laboratory for examining exoplanet dynamics, and the interaction between debris disks and planets. In recent years, HR8799's outer disk has been relatively well characterised, primarily using the Herschel Space Observatory. In contrast, the inner disk, too close to HR8799 to be spatially resolved by Herschel, remains poorly understood. This leaves significant questions over both the location of the planetesimals responsible for producing the observed dust, and the physical properties of those grains. We have performed extensive simulations of HR8799's inner, unresolved debris belt, using UNSW Australia's supercomputing facility, Katana. Here, we present the results of integrations following the evolution of a belt of dynamically hot debris interior to the orbit of HR8799e, for a period of 60 Myr, using an initial population of 500,000 massless test particles. These simulations have enable the characterisation of the extent and structure of the inner belt, revealing that its outer edge must lie interior to the 3:1 mean-motion resonance with HR8799, at approximately 7.5au, and highlighting the presence of fine structure analogous to the Solar system's Kirkwood gaps. In the future, out results will allow us to calculate a first estimate of the small-body impact rate and water delivery prospects for any potential terrestrial planet(s) that might lurk, undetected, in the inner system., Comment: Published in the peer-reviewed proceedings of the 14th Australian Space Research Conference, held at the University of South Australia, Adelaide, 29th September - 1st October 2014. Proceedings ISBN: 13: 978-0-9775740-8-7; Editors: Wayne Short & Iver Cairns Note: arXiv abstract slightly abridged to meet the maximum allowed length

The field of exoplanetary science is one of the most rapidly growing areas of astrophysical research. As more planets are discovered around other stars, new techniques have been developed that have allowed astronomers to begin to characterise them. Two of the most important factors in understanding the evolution of these planets, and potentially determining whether they are habitable, are the behaviour of the winds of the host star and the way in which they interact with the planet. The purpose of this project is to reconstruct the magnetic fields of planet hosting stars from spectropolarimetric observations, and to use these magnetic field maps to inform simulations of the stellar winds in those systems using the Block Adaptive Tree Solar-wind Roe Upwind Scheme (BATS-R-US) code. The BATS-R-US code was originally written to investigate the behaviour of the Solar wind, and so has been altered to be used in the context of other stellar systems. These simulations will give information about the velocity, pressure and density of the wind outward from the host star. They will also allow us to determine what influence the winds will have on the space weather environment of the planet. This paper presents the preliminary results of these simulations for the star $\tau$ Bo\"otis, using a newly reconstructed magnetic field map based on previously published observations. These simulations show interesting structures in the wind velocity around the star, consistent with the complex topology of its magnetic field., Comment: 8 pages, 2 figures, accepted for publication in the peer-reviewed proceedings of the 14th Australian Space Research Conference, held at the University of South Australia, 29th September - 1st October 2014

The study of the planet-debris disk connection can shed light on the formation and evolution of planetary systems, and may help predict the presence of planets around stars with certain disk characteristics. In preliminary analyses of the Herschel DEBRIS and DUNES surveys, Wyatt et al. (2012) and Marshall et al. (2014) identified a tentative correlation between debris and low-mass planets. Here we use the cleanest possible sample out these surveys to assess the presence of such a correlation, discarding stars without known ages, with ages < 1 Gyr and with binary companions <100 AU, to rule out possible correlations due to effects other than planet presence. In our sample of 204 FGK stars, we do not find evidence that debris disks are more common or more dusty around stars harboring high-mass or low-mass planets compared to a control sample without identified planets, nor that debris disks are more or less common (or more or less dusty) around stars harboring multiple planets compared to single-planet systems. Diverse dynamical histories may account for the lack of correlations. The data show the correlation between the presence of high-mass planets and stellar metallicity, but no correlation between the presence of low-mass planets or debris and stellar metallicity. Comparing the observed cumulative distribution of fractional luminosity to those expected from a Gaussian distribution, we find that a distribution centered on the Solar system's value fits well the data, while one centered at 10 times this value can be rejected. This is of interest in the context of future terrestrial planet characterization because it indicates that there are good prospects for finding a large number of debris disk systems (i.e. with evidence of harboring the building blocks of planets) with exozodiacal emission low enough to be appropriate targets for an ATLAST-type mission to search for biosignatures., Comment: Accepted for publication in ApJ. 31 pages, 13 figures, 9 tables

The $\textit{Herschel}$ DEBRIS, DUNES and GT programmes observed 37 exoplanet host stars within 25 pc at 70, 100 and 160 $\mu$m with the sensitivity to detect far-infrared excess emission at flux density levels only an order of magnitude greater than that of the Solar system's Edgeworth-Kuiper belt. Here we present an analysis of that sample, using it to more accurately determine the (possible) level of dust emission from these exoplanet host stars and thereafter determine the links between the various components of these exoplanetary systems through statistical analysis. We have fitted the flux densities measured from recent \textit{Herschel} observations with a simple two parameter ($T_{d}$, $L_{\rm IR}/L_{\star}$) black body model (or to the 3-$\sigma$ upper limits at 100 $\mu$m). From this uniform approach we calculate the fractional luminosity, radial extent, dust temperature and disc mass. We then plotted the calculated dust luminosity or upper limits against the stellar properties, e.g. effective temperature, metallicity, age, and identified correlations between these parameters. A total of eleven debris discs are identified around the 37 stars in the sample. An incidence of ten cool debris discs around the Sun-like exoplanet host stars (29 $\pm$ 9 %) is consistent with the detection rate found by DUNES (20.2 $\pm$ 2.0 %). For the debris disc systems, the dust temperatures range from 20 to 80 K, and fractional luminosities ($L_{\rm IR}/L_{\star}$) between 2.4 $\times$10$^{-6}$ and 4.1 $\times$10$^{-4}$. In the case of non-detections, we calculated typical 3-$\sigma$ upper limits to the dust fractional luminosities of a few $\times10^{-6}$. We recover the previously identified correlation between stellar metallicity and hot Jupiter planets in our data set. We find a correlation between the increased presence of dust, lower planet masses and lower stellar metallicities. (abridged), Comment: 14 pages, 2 figures, 4 tables, A&A accepted

With improvements in exoplanet detection techniques, the number of multiple planet systems discovered is increasing, while the detection of potentially habitable Earth-mass planets remains complicated and thus requires new search strategies. Dynamical studies of known multiple planet systems are therefore a vital tool in the search for stable and habitable planet candidates. Here, we present a dynamical study of the three-planet system HD 204313 to determine whether it could harbour an Earth-like planet within its habitable zone for a sufficient time to develop life. We found two semi-stable regions in the system, but neither prove stable for long enough for a terrestrial planet to develop life. Our investigations suggest that overlapping weak and high order resonances may be responsible for these semi-stable regions. This study established a framework for a larger project that will study the dynamical stability of the habitable zone of multiple planet systems, providing a list of interesting targets for future habitable low-mass planet searches., Comment: 13 pages, 6 figures. Accepted for publication in the peer-reviewed proceedings of the 13th annual Australian Space Science Conference

We present far-infrared and sub-millimeter images of the eta Crv debris disk system obtained with Herschel and SCUBA-2, as well as Hubble Space Telescope visible and near-infrared coronagraphic images. In the 70 micron Herschel image, we clearly separate the thermal emission from the warm and cold belts in the system, find no evidence for a putative dust population located between them, and precisely determine the geometry of the outer belt. We also find marginal evidence for azimuthal asymmetries and a global offset of the outer debris ring relative to the central star. Finally, we place stringent upper limits on the scattered light surface brightness of the outer ring. Using radiative transfer modeling, we find that it is impossible to account for all observed properties of the system under the assumption that both rings contain dust populations with the same properties. While the outer belt is in reasonable agreement with the expectations of steady-state collisional cascade models, albeit with a minimum grain size that is four times larger than the blow-out size, the inner belt appears to contain copious amounts of small dust grains, possibly below the blow-out size. This suggests that the inner belt cannot result from a simple transport of grains from the outer belt and rather supports a more violent phenomenon as its origin. We also find that the emission from the inner belt has not declined over three decades, a much longer timescale than its dynamical timescale, which indicates that the belt is efficiently replenished., Comment: 50 pages, 12 figures, 5 tables, accepted for publication in the Astrophysical Journal

Many nearby main-sequence stars have been searched for debris using the far-infrared Herschel satellite, within the DEBRIS, DUNES and Guaranteed-Time Key Projects. We discuss here 11 stars of spectral types A to M where the stellar inclination is known and can be compared to that of the spatially-resolved dust belts. The discs are found to be well aligned with the stellar equators, as in the case of the Sun's Kuiper belt, and unlike many close-in planets seen in transit surveys. The ensemble of stars here can be fitted with a star-disc tilt of ~<10 degrees. These results suggest that proposed mechanisms for tilting the star or disc in fact operate rarely. A few systems also host imaged planets, whose orbits at tens of AU are aligned with the debris discs, contrary to what might be expected in models where external perturbers induce tilts., Comment: accepted by MNRAS Letters

2012 DR30 is a recently discovered Solar System object on a unique orbit, with a high eccentricity of 0.9867, a perihelion distance of 14.54 AU and a semi-major axis of 1109 AU, in this respect outscoring the vast majority of trans-Neptunian objects. We performed Herschel/PACS and optical photometry to uncover the size and albedo of 2012 DR30, together with its thermal and surface properties. The body is 185 km in diameter and has a relatively low V-band geometric albedo of ~8%. Although the colours of the object indicate that 2012 DR30 is an RI taxonomy class TNO or Centaur, we detected an absorption feature in the Z-band that is uncommon among these bodies. A dynamical analysis of the target's orbit shows that 2012 DR30 moves on a relatively unstable orbit and was most likely only recently placed on its current orbit from the most distant and still highly unexplored regions of the Solar System. If categorised on dynamical grounds 2012 DR30 is the largest Damocloid and/or high inclination Centaur observed so far., Comment: Astronomy and Astrophysics, accepted

To understand the frequency, and thus the formation and evolution, of planetary systems like our own solar system, it is critical to detect Jupiter-like planets in Jupiter-like orbits. For long-term radial-velocity monitoring, it is useful to estimate the observational effort required to reliably detect such objects, particularly in light of severe competition for limited telescope time. We perform detailed simulations of observational campaigns, maximizing the realism of the sampling of a set of simulated observations. We then compute the detection limits for each campaign to quantify the effect of increasing the number of observational epochs and varying their time coverage. We show that once there is sufficient time baseline to detect a given orbital period, it becomes less effective to add further time coverage -- rather, the detectability of a planet scales roughly as the square root of the number of observations, independently of the number of orbital cycles included in the data string. We also show that no noise floor is reached, with a continuing improvement in detectability at the maximum number of observations N=500 tested here., Comment: Accepted for publication in PASP

Debris disks have been found primarily around intermediate and solar mass stars (spectral types A-K) but rarely around low mass M-type stars. We have spatially resolved a debris disk around the remarkable M3-type star GJ581 hosting multiple planets using deep PACS images at 70, 100 and 160 microns as part of the DEBRIS Program on the Herschel Space Observatory. This is the second spatially resolved debris disk found around an M-type star, after the one surrounding the young star AU Mic (12 Myr). However, GJ 581 is much older (2-8 Gyr), and is X-ray quiet in the ROSAT data. We fit an axisymmetric model of the disk to the three PACS images and found that the best fit model is for a disk extending radially from 25+/-12 AU to more than 60 AU. Such a cold disk is reminiscent of the Kuiper Belt but it surrounds a low mass star (0.3 M_sol) and its fractional dust luminosity L_dust/L_* of \sim 10^-4 is much higher. The inclination limits of the disk found in our analysis make the masses of the planets small enough to ensure the long-term stability of the system according to some dynamical simulations. The disk is collisionally dominated down to submicron-sized grains and the dust cannot be expelled from the system by radiation or wind pressures because of the low luminosity and low X-ray luminosity of GJ581. We suggest that the correlation between low-mass planets and debris disks recently found for G-type stars also applies to M-type stars.Finally, the known planets, of low masses and orbiting within 0.3 AU from the star, cannot dynamically perturb the disk over the age of the star, suggesting that an additional planet exists at larger distance that is stirring the disk to replenish the dust., Comment: To appear in Astronomy and Astrophysics, 15 pages, 10 figures

In 2009, the discovery of two planets orbiting the evolved binary star system HW Virginis was announced, based on systematic variations in the timing of eclipses between the two stars. The planets invoked in that work were significantly more massive than Jupiter, and moved on orbits that were mutually crossing - an architecture which suggests that mutual encounters and strong gravitational interactions are almost guaranteed. In this work, we perform a highly detailed analysis of the proposed HW Vir planetary system. First, we consider the dynamical stability of the system as proposed in the discovery work. Through a mapping process involving 91,125 individual simulations, we find that the system is so unstable that the planets proposed simply cannot exist, due to mean lifetimes of less than a thousand years across the whole parameter space. We then present a detailed re-analysis of the observational data on HW Vir, deriving a new orbital solution that provides a very good fit to the observational data. Our new analysis yields a system with planets more widely spaced, and of lower mass, than that proposed in the discovery work, and yields a significantly greater (and more realistic) estimate of the uncertainty in the orbit of the outermost body. Despite this, a detailed dynamical analysis of this new solution similarly reveals that it also requires the planets to move on orbits that are simply not dynamically feasible. Our results imply that some mechanism other than the influence of planetary companions must be the principal cause of the observed eclipse timing variations for HW Vir. If the sys- tem does host exoplanets, they must move on orbits differing greatly from those previously proposed. Our results illustrate the critical importance of performing dynamical analyses as a part of the discovery process for multiple-planet exoplanetary systems., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society

This paper describes Herschel observations of the nearby (8.5pc) G5V multi-exoplanet host star 61 Vir at 70-500micron carried out by the DEBRIS survey. These reveal emission that is extended out to >15arcsec with a morphology that can be fitted by a nearly edge-on (77deg inclination) radially broad (from 30AU to >100AU) debris disk of fractional luminosity 2.7x10^-5, with two unrelated sources nearby that are more prominent at longer wavelengths. Chance alignment with a background object seen at 1.4GHz provides potential for confusion, but the star's 1.4"/yr proper motion allows Spitzer 70micron images to confirm that what we are interpreting as disk emission really is circumstellar. Although the exact shape of the disk's inner edge is not well constrained, the region inside 30AU must be significantly depleted in planetesimals. This is readily explained if there are additional planets in the 0.5-30AU region, but is also consistent with collisional erosion. We also find tentative evidence that the presence of detectable debris around nearby stars correlates with the presence of the lowest mass planets that are detectable in current radial velocity surveys. Out of an unbiased sample of the nearest 60 G stars, 11 are known to have planets, of which 6 (including 61 Vir) have planets that are all less massive than Saturn, and 4 of these have evidence for debris. The debris toward one of these planet-hosts (HD20794) is reported here for the first time. This fraction (4/6) is higher than that expected for nearby field stars (15%), and implies that systems that form low-mass planets are also able to retain bright debris disks. We suggest that this correlation could arise because such planetary systems are dynamically stable and include regions that are populated with planetesimals in the formation process where the planetesimals can remain unperturbed over Gyr timescales., Comment: 20 pages. MNRAS, in press

We present 8 years of high-precision radial velocity (RV) data for HD 204313 from the 2.7 m Harlan J. Smith Telescope at McDonald Observatory. The star is known to have a giant planet (M sin i = 3.5 M_J) on a ~1900-day orbit, and a Neptune-mass planet at 0.2 AU. Using our own data in combination with the published CORALIE RVs of Segransan et al. (2010), we discover an outer Jovian (M sin i = 1.6 M_J) planet with P ~ 2800 days. Our orbital fit suggests the planets are in a 3:2 mean motion resonance, which would potentially affect their stability. We perform a detailed stability analysis, and verify the planets must be in resonance., Comment: Accepted for publication in ApJ

We report the detection of two new planets from the Anglo-Australian Planet Search. These planets orbit two stars each previously known to host one planet. The new planet orbiting HD 142 has a period of 6005\pm427 days, and a minimum mass of 5.3M_Jup. HD142c is thus a new Jupiter analog: a gas-giant planet with a long period and low eccentricity (e = 0.21 \pm 0.07). The second planet in the HD 159868 system has a period of 352.3\pm1.3 days, and m sin i=0.73\pm0.05 M_Jup. In both of these systems, including the additional planets in the fitting process significantly reduced the eccentricity of the original planet. These systems are thus examples of how multiple-planet systems can masquerade as moderately eccentric single-planet systems., Comment: 31 pages, 8 figures, accepted for publication in ApJ

We have performed detailed thermophysical and dynamical modelling of Jovian Trojan (1173) Anchises. Our results reveal a most unusual object. By examining observational data taken by IRAS, Akari and WISE between 11.5 and 60 microns, along with variations in its optical lightcurve, we find Anchises is most likely an elongated body, with an axes-ratio of ~1.4. This yields calculated best-fit dimensions of 170x121x121km (an equivalent diameter of 136+18/-11km). We find the observations are best fit by Anchises having a retrograde sense of rotation, and an unusually high thermal inertia (25 to 100 Jm-2s-0.5K-1). The geometric albedo is found to be 0.027 (+0.006/-0.007). Anchises therefore has one of the highest published thermal inertias of any object larger than 100km in diameter, at such large heliocentric distances, and is one of the lowest albedo objects ever observed. More observations are needed to see if there is a link between the very shallow phase curve, with almost no opposition effect, and the derived thermal properties for this large Trojan asteroid. Our dynamical investigation of Anchises' orbit has revealed it to be dynamically unstable on timescales of hundreds of Myr, similar to the unstable Neptunian Trojans 2001 QR322 and 2008 LC18. Unlike those objects, we find that Anchises' dynamical stability is not a function of its initial orbital elements, the result of the exceptional precision with which its orbit is known. This is the first time that a Jovian Trojan has been shown to be dynamically unstable, and adds weight to the idea that planetary Trojans represent a significant ongoing contribution to the Centaur population, the parents of the short-period comets. The observed instability does not rule out a primordial origin for Anchises, but when taken in concert with the result of our thermophysical analysis, suggest that it would be a fascinating target for future study., Comment: 5 figures, 3 tables, accepted for publication in Monthly Notices of the Royal Astronomical Society

The recent discovery of the first Neptune Trojan at the planet's trailing (L5) Lagrange point, 2008 LC18, offers an opportunity to confirm the formation mechanism of a member of this important tracer population for the Solar system's dynamical history. We tested the stability of 2008 LC18's orbit through a detailed dynamical study, using test particles spread across the orbital uncertainties in a, e, i and {\Omega}. This showed that the wide uncertainties of the published orbit span regions of both extreme dynamical instability, with lifetimes < 100 Myr, and with significant stability (> 1 Gyr lifetimes). The stability of 2008 LC18's clones is greatly dependent on their semi-major axis and only weakly correlated with their eccentricity. Test particles on orbits with an initial semi-major axis less than 29.91 AU have dynamical half-lives shorter than 100 Myr; in contrast, particles with an initial semi-major axis greater than 29.91 AU exhibit such strong dynamical stability that almost all are retained over the 1 Gyr of our simulations. More observations of this object are necessary to improve the orbit. If 2008 LC18 is in the unstable region, then our simulations imply that it is either a temporary Trojan capture, or a representative of a slowly decaying Trojan population (like its sibling the L4 Neptunian Trojan 2001 QR322), and that it may not be primordial. Alternatively, if the orbit falls into the larger, stable region, then 2008 LC18 is a primordial member of the highly stable and highly inclined component of the Neptune Trojan population, joining 2005 TN53 and 2007 VL305. We attempted to recover 2008 LC18 using the 2.3m telescope at Siding Spring Observatory to provide this astrometry, but were unsuccessful due to the high stellar density of its current sky location near the galactic centre. The recovery of this object will require a telescope in the 8m class., Comment: 11 pages, 3 figures, 1 table, accepted for publication in Monthly Notices of the Royal Astronomical Society