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.

AU Microscopii (AU Mic) is the second closest pre main sequence star, at a distance of 9.79 parsecs and with an age of 22 million years. AU Mic possesses a relatively rare and spatially resolved3 edge-on debris disk extending from about 35 to 210 astronomical units from the star, and with clumps exhibiting non-Keplerian motion. Detection of newly formed planets around such a star is challenged by the presence of spots, plage, flares and other manifestations of magnetic activity on the star. Here we report observations of a planet transiting AU Mic. The transiting planet, AU Mic b, has an orbital period of 8.46 days, an orbital distance of 0.07 astronomical units, a radius of 0.4 Jupiter radii, and a mass of less than 0.18 Jupiter masses at 3 sigma confidence. Our observations of a planet co-existing with a debris disk offer the opportunity to test the predictions of current models of planet formation and evolution.

We report measurements of the sky-projected spin-orbit angle for AU\,Mic\,b, a Neptune-size planet orbiting a very young ($\sim20$\,Myr) nearby pre-main sequence M dwarf star which also hosts a bright, edge-on, debris disk. The planet was recently discovered from preliminary analysis of radial velocity observations and confirmed to be transiting its host star from photometric data from the NASA's \textit{TESS} mission. We obtained radial velocity measurements of the Rossiter--McLaughlin effect for AU\,Mic\,b over the course of two partially observable transits and one full transit from high-resolution spectroscopic observations made with the {\textsc{Minerva}}-Australis telescope array. We find that the planet's orbit is aligned with the spin-axis of its host star with a sky-projected spin-orbit angle of $\lambda = 5{^{+16}_{-15}}^{\circ}$. This measurement is in excellent agreement with observations carried out on other instruments around the same time. AU\,Mic is the youngest exoplanetary system for which the projected spin-orbit angle has been measured, making it a key data point in the study of the formation and migration of exoplanets -- particularly given that the system is also host to a bright debris disk. Given that we find this system (star, planet, and disk) to be in near spin-orbit alignment, it seems likely that the planet formed farther out beyond the ice-line and then migrated inwards to its current location via quiescent migration through the disk, rather than through mechanisms that would excite the inclination of the planet's orbit.

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.

Lithium depletion and enrichment in the cosmos is not yet well understood. To help tighten constraints on stellar and Galactic evolution models, we present the largest high-resolution analysis of Li abundances A(Li) to date, with results for over 100 000 GALAH field stars spanning effective temperatures $5900\,\mathrm{K} \lesssim \rm{T_{eff}} \lesssim7000\,\mathrm{K}$ and metallicities $-3 \lesssim \rm[Fe/H] \lesssim +0.5$. We separated these stars into two groups, on the warm and cool side of the so-called Li-dip, a localised region of the Kiel diagram wherein lithium is severely depleted. We discovered that stars in these two groups show similar trends in the A(Li)-[Fe/H] plane, but with a roughly constant offset in A(Li) of 0.4 dex, the warm group having higher Li abundances. At $\rm[Fe/H]\gtrsim-0.5$, a significant increasing in Li abundance with increasing metallicity is evident in both groups, signalling the onset of significant Galactic production. At lower metallicity, stars in the cool group sit on the Spite plateau, showing a reduced lithium of around 0.4 dex relative to the primordial value predicted from Big Bang nucleosynthesis (BBN). However, stars in the warm group between [Fe/H] = -1.0 and -0.5, form an elevated plateau that is largely consistent with the BBN prediction. This may indicate that these stars in fact preserve the primordial Li produced in the early Universe.

We present a neural network autoencoder structure that is able to extract essential latent spectral features from observed spectra and then reconstruct a spectrum from those features. Because of the training with a set of unpeculiar spectra, the network is able to reproduce a spectrum of high signal-to-noise ratio that does not show any spectral peculiarities even if they are present in an observed spectrum. Spectra generated in this manner were used to identify various emission features among spectra acquired by multiple surveys using the HERMES spectrograph at the Anglo-Australian telescope. Emission features were identified by a direct comparison of the observed and generated spectra. Using the described comparison procedure, we discovered 10,364 candidate spectra with a varying degree of H$\alpha$/H$\beta$ emission component produced by different physical mechanisms. A fraction of those spectra belongs to the repeated observation that shows temporal variability in their emission profile. Among emission spectra, we find objects that feature contributions of a nearby rarefied gas (identified through the emission of [NII] and [SII] lines) that was identified in 4004 spectra, which were not all identified as having H$\alpha$ emission. Positions of identified emission-line objects coincide with multiple known regions that harbour young stars. Similarly, detected nebular emission spectra coincide with visually-prominent nebular clouds observable in the red all-sky photographic composites.

We investigate the properties of 1306 red giant stars with high photospheric abundances of lithium observed by the GALAH, K2-HERMES and TESS-HERMES surveys, and discuss them in the context of proposed mechanisms for lithium enrichment in giant stars. We confirm that Li-rich giants are rare, making up only 1.1 per cent of our giant star sample. We use GALAH+ DR3 stellar parameters and a Bayesian isochrone analysis to divide the sample into first-ascent red giant branch and red clump stars, and confirm these classifications using asteroseismic data from K2. We find that red clump stars are more than three times as likely to be lithium-rich as red giant branch stars, and the occurrence rate of lithium richness with metallicity is quite different between the two populations. The probability for a rapidly rotating giant to be lithium-rich is distinctly higher as compared to a non rapidly rotating giant, independent of its evolutionary state. There is also a clear correspondence on the red clump between rapid rotation and extreme levels of lithium enrichment. There are almost no stars on the secondary red clump with lithium abundances above the primordial level, indicating a mass dependence in the red clump lithium enrichment mechanism(s). The complex distribution of lithium-rich giants across evolutionary phase, metallicity, rotation rate and mass implies multiple independent mechanisms for producing lithium enrichment.

Accurate and precise radius estimates of transiting exoplanets are critical for understanding their compositions and formation mechanisms. To know the planet, we must know the host star in as much detail as possible. We present complete results for planet-candidate hosts from the K2-HERMES survey, which uses the HERMES multi-object spectrograph on the Anglo-Australian Telescope to obtain R$\sim$28,000 spectra for more than 30,000 K2 stars. We present complete host-star parameters and planet-candidate radii for 224 K2 candidate planets from C1-C13. Our results cast severe doubt on 30 K2 candidates, as we derive unphysically large radii, larger than 2 $R_{Jup}$. This work highlights the importance of obtaining accurate, precise, and self-consistent stellar parameters for ongoing large planet search programs - something that will only become more important in the coming years, as TESS begins to deliver its own harvest of exoplanets.

The Jovian Trojans are two swarms of objects located around the L$_4$ and L$_5$ Lagrange points. The population is thought to have been captured by Jupiter during the Solar system's youth. Within the swarms, six collisional families have been identified in previous work, with four in the L$_4$ swarm, and two in the L$_5$. Our aim is to investigate the stability of the two Trojan swarms, with a particular focus on these collisional families. We find that the members of Trojan swarms escape the population at a linear rate, with the primordial L$_4$ (23.35% escape) and L$_5$ (24.89% escape) population sizes likely 1.31 and 1.35 times larger than today. Given that the escape rates were approximately equal between the two Trojan swarms, our results do not explain the observed asymmetry between the two groups, suggesting that the numerical differences are primordial in nature, supporting previous studies. Upon leaving the Trojan population, the escaped objects move onto orbits that resemble those of the Centaur and short-period comet populations. Within the Trojan collisional families, the 1996 RJ and 2001 UV$_{209}$ families are found to be dynamically stable over the lifetime of the Solar system, whilst the Hektor, Arkesilos and Ennomos families exhibit various degrees of instability. The larger Eurybates family shows 18.81\% of simulated members escaping the Trojan population. Unlike the L4 swarm, the escape rate from the Eurybates family is found to increase as a function of time, allowing an age estimation of approximately $1.045\pm 0.364 \times 10^9$ years.

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.

We explore the fundamental relations governing the radial and vertical velocity dispersions of stars in the Milky Way, from combined studies of complementary surveys including GALAH, LAMOST, APOGEE, the NASA $Kepler$ and K2 missions, and $Gaia$ DR2. We find that different stellar samples, even though they target different tracer populations and employ a variety of age estimation techniques, follow the same set of fundamental relations. We provide the clearest evidence to date that, in addition to the well-known dependence on stellar age, the velocity dispersions of stars depend on orbital angular momentum $L_z$, metallicity and height above the plane $|z|$, and are well described by a multiplicatively separable functional form. The dispersions have a power-law dependence on age with exponents of 0.441$\pm 0.007$ and 0.251$\pm 0.006$ for $\sigma_z$ and $\sigma_R$ respectively, and the power law is valid even for the oldest stars. For the solar neighborhood stars, the apparent break in the power law for older stars, as seen in previous studies, is due to the anti-correlation of $L_z$ with age. The dispersions decrease with increasing $L_z$ until we reach the Sun's orbital angular momentum, after which $\sigma_z$ increases (implying flaring in the outer disc) while $\sigma_R$ flattens. The dispersions increase with decreasing metallicity, suggesting that the dispersions increase with birth radius. The dispersions also increase linearly with $|z|$. The same set of relations that work in the solar neighborhood also work for stars between $3<R/{\rm kpc}<20$. Finally, the high-[$\alpha$/Fe] stars follow the same relations as the low-[$\alpha$/Fe] stars.

We report the discovery of a warm sub-Saturn, TOI-257b (HD 19916b), based on data from NASA's Transiting Exoplanet Survey Satellite (TESS). The transit signal was detected by TESS and confirmed to be of planetary origin based on radial-velocity observations with the Minerva-Australis telescope array. An analysis of the TESS photometry, the Minerva-Australis, FEROS, and HARPS radial velocities, and the asteroseismic data of the stellar oscillations reveals that TOI-257b has a mass of $M_P=0.134^{+0.023}_{-0.022}$$\rm{M_J}$ ($42.6^{+7.3}_{-7.0}$$\rm{M_{\oplus}}$), a radius of $R_P=0.626^{+0.013}_{-0.012}$$\rm{R_J}$ ($7.02^{+0.15}_{-0.13}$$\rm{R_{\oplus}}$), and an orbit with eccentricity $0.242^{+0.040}_{-0.065}$ and period $18.38827\pm0.00072$$\rm{days}$. TOI-257b orbits a bright ($\mathrm{V}=7.570$mag) somewhat evolved late F-type star with $M_*=1.390\pm0.046$$\rm{M_{\odot}}$, $R_*=1.888\pm0.033$$\rm{R_{\odot}}$, $T_{\rm eff}=6075\pm90$$\rm{K}$, and $v\sin{i}=11.3\pm0.5$km/s. Additionally, we statistically validate a second non-transiting sub-Saturn mass planet on a $\sim71$day orbit using the radial velocity data. This system joins the ranks of a small number of exoplanet host stars that have been characterized with asteroseismology. Warm sub-Saturns are rare in the known sample of exoplanets, and thus the discovery of TOI-257b is important in the context of future work studying the formation and migration history of similar planetary systems.

A wealth of Earth-sized exoplanets will be discovered in the coming years, proving a large pool of candidates from which the targets for the search for life beyond the Solar system will be chosen. The target selection process will require the leveraging of all available information in order to maximise the robustness of the target list and make the most productive use of follow-up resources. Here, we present the results of a suite of $n$-body simulations that demonstrate the degree to which the orbital architecture of the Solar system impacts the variability of Earth's orbital elements. By varying the orbit of Jupiter and keeping the initial orbits of the other planets constant, we demonstrate how subtle changes in Solar system architecture could alter the Earth's orbital evolution -- a key factor in the Milankovitch cycles that alter the amount and distribution of solar insolation, thereby driving periodic climate change on our planet. The amplitudes and frequencies of Earth's modern orbital cycles fall in the middle of the range seen in our runs for all parameters considered -- neither unusually fast nor slow, nor large nor small. This finding runs counter to the `Rare Earth' hypothesis, which suggests that conditions on Earth are so unusual that life elsewhere is essentially impossible. Our results highlight how dynamical simulations of newly discovered exoplanetary systems could be used as an additional means to assess the potential targets of biosignature searches, and thereby help focus the search for life to the most promising targets.

Our understanding of planetary systems different to our own has grown dramatically in the past 30 years. However, our efforts to ascertain the degree to which the Solar system is abnormal or unique have been hindered by the observational biases inherent to the methods that have yielded the greatest exoplanet hauls. On the basis of such surveys, one might consider our planetary system highly unusual - but the reality is that we are only now beginning to uncover the true picture. In this work, we use the full eighteen-year archive of data from the Anglo-Australian Planet Search to examine the abundance of 'Cool Jupiters' - analogs to the Solar system's giant planets, Jupiter and Saturn. We find that such planets are intrinsically far more common through the cosmos than their siblings, the hot Jupiters. We find that the occurrence rate of such 'Cool Jupiters' is $6.73^{+2.09}_{-1.13}$\%, almost an order of magnitude higher than the occurrence of hot Jupiters (at $0.84^{+0.70}_{-0.20}$\%). We also find that the occurrence rate of giant planets is essentially constant beyond orbital distances of $\sim$1\,au. Our results reinforce the importance of legacy radial velocity surveys for the understanding of the Solar system's place in the cosmos.

We present the discoveries of KELT-25b (TIC 65412605, TOI-626.01) and KELT-26b (TIC 160708862, TOI-1337.01), two transiting companions orbiting relatively bright, early A-stars. The transit signals were initially detected by the KELT survey, and subsequently confirmed by \textit{TESS} photometry. KELT-25b is on a 4.40-day orbit around the V = 9.66 star CD-24 5016 ($T_{\rm eff} = 8280^{+440}_{-180}$ K, $M_{\star}$ = $2.18^{+0.12}_{-0.11}$ $M_{\odot}$), while KELT-26b is on a 3.34-day orbit around the V = 9.95 star HD 134004 ($T_{\rm eff}$ =$8640^{+500}_{-240}$ K, $M_{\star}$ = $1.93^{+0.14}_{-0.16}$ $M_{\odot}$), which is likely an Am star. We have confirmed the sub-stellar nature of both companions through detailed characterization of each system using ground-based and \textit{TESS} photometry, radial velocity measurements, Doppler Tomography, and high-resolution imaging. For KELT-25, we determine a companion radius of $R_{\rm P}$ = $1.64^{+0.039}_{-0.043}$ $R_{\rm J}$, and a 3-sigma upper limit on the companion's mass of $\sim64~M_{\rm J}$. For KELT-26b, we infer a planetary mass and radius of $M_{\rm P}$ = $1.41^{+0.43}_{-0.51}$ $M_{\rm J}$ and $R_{\rm P}$ = $1.940^{+0.060}_{-0.058}$ $R_{\rm J}$. From Doppler Tomographic observations, we find KELT-26b to reside in a highly misaligned orbit. This conclusion is weakly corroborated by a subtle asymmetry in the transit light curve from the \textit{TESS} data. KELT-25b appears to be in a well-aligned, prograde orbit, and the system is likely a member of a cluster or moving group.

Our knowledge of the populations and occurrence rates of planets orbiting evolved intermediate-mass stars lags behind that for solar-type stars by at least a decade. Some radial velocity surveys have targeted these low-luminosity giant stars, providing some insights into the properties of their planetary systems. Here we present the final data release of the Pan-Pacific Planet Search, a 5-year radial velocity survey using the 3.9m Anglo-Australian Telescope. We present 1293 precise radial velocity measurements for 129 stars, and highlight six potential substellar-mass companions which require additional observations to confirm. Correcting for the substantial incompleteness in the sample, we estimate the occurrence rate of giant planets orbiting low-luminosity giant stars to be approximately 7.8$^{+9.1}_{-3.3}$\%. This result is consistent with the frequency of such planets found to orbit main-sequence A-type stars, from which the PPPS stars have evolved.

We report the discovery of TOI-677 b, first identified as a candidate in light curves obtained within Sectors 9 and 10 of the Transiting Exoplanet Survey Satellite (TESS) mission and confirmed with radial velocities. TOI-677 b has a mass of M_p = 1.236$^{+0.069}_{-0.067}$ M_J, a radius of R_p = 1.170 +- 0.03 R_J,and orbits its bright host star (V=9.8 mag) with an orbital period of 11.23660 +- 0.00011 d, on an eccentric orbit with e = 0.435 +- 0.024. The host star has a mass of M_* = 1.181 +- 0.058 M_sun, a radius of R_* = 1.28 +- 0.03 R_sun, an age of 2.92$^{+0.80}_{-0.73}$ Gyr and solar metallicity, properties consistent with a main sequence late F star with T_eff = 6295 +- 77 K. We find evidence in the radial velocity measurements of a secondary long term signal which could be due to an outer companion. The TOI-677 b system is a well suited target for Rossiter-Mclaughlin observations that can constrain migration mechanisms of close-in giant planets.

We present a computationally tractable implementation of spectro-perfectionism, a method which minimizes error imparted by spectral extraction. We develop our method in conjunction with a full raw reduction pipeline for the MINiature Exoplanet Radial Velocity Array (MINERVA), capable of performing both optimal extraction and spectro-perfectionism. Although spectro-perfectionism remains computationally expensive, our implementation can extract a MINERVA exposure in approximately $30\,\text{min}$. We describe our localized extraction procedure and our approach to point spread function fitting. We compare the performance of both extraction methods on a set of 119 exposures on HD122064, an RV standard star. Both the optimal extraction and spectro-perfectionism pipelines achieve nearly identical RV precision under a six-exposure chronological binning. We discuss the importance of reliable calibration data for point spread function fitting and the potential of spectro-perfectionism for future precise radial velocity exoplanet studies.

We present isochrone ages and initial bulk metallicities ($\rm [Fe/H]_{bulk}$, by accounting for diffusion) of 163,722 stars from the GALAH Data Release 2, mainly composed of main sequence turn-off stars and subgiants ($\rm 7000 K>T_{eff}>4000 K$ and $\rm log g>3$ dex). The local age-metallicity relationship (AMR) is nearly flat but with significant scatter at all ages; the scatter is even higher when considering the observed surface abundances. After correcting for selection effects, the AMR appear to have intrinsic structures indicative of two star formation events, which we speculate are connected to the thin and thick disks in the solar neighborhood. We also present abundance ratio trends for 16 elements as a function of age, across different $\rm [Fe/H]_{bulk}$ bins. In general, we find the trends in terms of [X/Fe] vs age from our far larger sample to be compatible with studies based on small ($\sim$ 100 stars) samples of solar twins but we now extend it to both sub- and super-solar metallicities. The $\alpha$-elements show differing behaviour: the hydrostatic $\alpha$-elements O and Mg show a steady decline with time for all metallicities while the explosive $\alpha$-elements Si, Ca and Ti are nearly constant during the thin disk epoch (ages $\lessapprox $ 12 Gyr). The s-process elements Y and Ba show increasing [X/Fe] with time while the r-process element Eu have the opposite trend, thus favouring a primary production from sources with a short time-delay such as core-collapse supernovae over long-delay events such as neutron star mergers.

Using kinematics from Gaia and the large elemental abundance space of the second data release of the GALAH survey, we identify two new members of the Fimbulthul stellar stream, and chemically tag them to massive, multi-metallic globular cluster $\omega$ Centauri. Recent analysis of the second data release of Gaia had revealed the Fimbulthul stellar stream in the halo of the Milky Way. It had been proposed that the stream is associated with the $\omega$ Centauri, but this proposition relied exclusively upon the kinematics and metallicities of the stars to make the association. In this work, we find our two new members of the stream to be metal-poor stars that are enhanced in sodium and aluminium, typical of second population globular cluster stars, but not otherwise seen in field stars. Furthermore, the stars share the s-process abundance pattern seen in $\omega$ Centauri, which is rare in field stars. Apart from one star within 1.5 deg of $\omega$ Centauri, we find no other stars observed by GALAH spatially near $\omega$ Centauri or the Fimbulthul stream that could be kinematically and chemically linked to the cluster. Chemically tagging stars in the Fimbulthul stream to $\omega$ Centauri confirms the earlier work, and further links this tidal feature in the Milky Way halo to $\omega$ Centauri.

The Minerva-Australis telescope array is a facility dedicated to the follow-up, confirmation, characterisation, and mass measurement of bright transiting planets discovered by the Transiting Exoplanet Survey Satellite (TESS) -- a category in which it is almost unique in the southern hemisphere. It is located at the University of Southern Queensland's Mount Kent Observatory near Toowoomba, Australia. Its flexible design enables multiple 0.7m robotic telescopes to be used both in combination, and independently, for high-resolution spectroscopy and precision photometry of TESS transit planet candidates. Minerva-Australis also enables complementary studies of exoplanet spin-orbit alignments via Doppler observations of the Rossiter-McLaughlin effect, radial velocity searches for non-transiting planets, planet searches using transit timing variations, and ephemeris refinement for TESS planets. In this first paper, we describe the design, photometric instrumentation, software, and science goals of Minerva-Australis, and note key differences from its Northern hemisphere counterpart -- the Minerva array. We use recent transit observations of four planets--WASP-2b, WASP-44b, WASP-45b, and HD 189733b to demonstrate the photometric capabilities of Minerva-Australis.

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.

A significant gender disparity is widely known to occur within the fields of science, technology, engineering and mathematics (STEM). Women are both under-represented in their participation in Australian STEM education (compared to their fraction within the population at large), and face a much higher attrition rate from STEM subjects throughout secondary and tertiary education. Even within STEM careers the proportion of women frequently generally decreases with increasing position of seniority. In 2015 organisers and members of the Australian Space Research Community began to analyse the data from participants of the Australian Space Research Conference, in order to derive statistics relating to the gender balance of the conference. In this paper, the data from the most recent conference (held on the Gold Coast, in September 2018) is analysed, considering the gender demographics of delegates, presenters (within numerous categories), and awards. This year, we also present the dimensions of career type and academic level of conference attendees. The resulting trends are compared to those of other national space research conferences - the Scientific Meeting of the Astronomical Society of Australia. The results show that the representation of females dropped to the lowest levels since formal data collection began in 2015, with females comprising 23% of conference attendees, and 21% of all session presentations (posters and talks combined). Gender parity was achieved however amongst the invited Plenary presentations. Recommendations for removing some of the barriers to female participation in the conference and encouraging greater gender balance and diversity in future are discussed.

Transiting Exoplanet Survey Satellite (TESS) observations have revealed a compact multi-planet system around the sixth-magnitude star HR 858 (TIC 178155732, TOI 396), located 32 parsecs away. Three planets, each about twice the size of Earth, transit this slightly-evolved, late F-type star, which is also a member of a visual binary. Two of the planets may be in mean motion resonance. We analyze the TESS observations, using novel methods to model and remove instrumental systematic errors, and combine these data with follow-up observations taken from a suite of ground-based telescopes to characterize the planetary system. The HR 858 planets are enticing targets for precise radial velocity observations, secondary eclipse spectroscopy, and measurements of the Rossiter-McLaughlin effect.

We present a detailed analysis of the orbital stability of the HD 181433 planetary system, finding it to exhibit strong dynamical instability across a wide range of orbital eccentricities, semi-major axes, and mutual inclinations. We also analyse the behaviour of an alternative system architecture, proposed by Campanella (2011), and find that it offers greater stability than the original solution, as a result of the planets being trapped in strong mutual resonance. We take advantage of more recent observations to perform a full refit of the system, producing a new planetary solution. The best-fit orbit for HD 181433 d now places the planet at a semi-major axis of 6.60$\pm$0.22 au, with an eccentricity of 0.469$\pm$0.013. Extensive simulations of this new system architecture reveal it to be dynamically stable across a broad range of potential orbital parameter space, increasing our confidence that the new solution represents the ground truth of the system. Our work highlights the advantage of performing dynamical simulations of candidate planetary systems in concert with the orbital fitting process, as well as supporting the continuing monitoring of radial velocity planet search targets.

The sensitivity of radial velocity (RV) surveys for exoplanet detection are extending to increasingly long orbital periods, where companions with periods of several years are now being regularly discovered. Companions with orbital periods that exceed the duration of the survey manifest in the data as an incomplete orbit or linear trend, a feature that can either present as the sole detectable companion to the host star, or as an additional signal overlain on the signatures of previously discovered companion(s). A diagnostic that can confirm or constrain scenarios in which the trend is caused by an unseen stellar, rather than planetary, companion is the use of high-contrast imaging observations. Here, we present RV data from the Anglo-Australian Planet Search (AAPS) for twenty stars that show evidence of orbiting companions. Of these, six companions have resolved orbits, with three that lie in the planetary regime. Two of these (HD~92987b and HD~221420b) are new discoveries. Follow-up observations using the Differential Speckle Survey Instrument (DSSI) on the Gemini South telescope revealed that five of the twenty monitored companions are likely stellar in nature. We use the sensitivity of the AAPS and DSSI data to place constraints on the mass of the companions for the remaining systems. Our analysis shows that a planetary-mass companion provides the most likely self-consistent explanation of the data for many of the remaining systems.

The latest Gaia data release enables us to accurately identify stars that are more luminous than would be expected on the basis of their spectral type and distance. During an investigation of the 329 best Solar twin candidates uncovered among the spectra acquired by the GALAH survey, we identified 64 such over-luminous stars. In order to investigate their exact composition, we developed a data-driven methodology that can generate a synthetic photometric signature and spectrum of a single star. By combining multiple such synthetic stars into an unresolved binary or triple system and comparing the results to the actual photometric and spectroscopic observations, we uncovered 6 definitive triple stellar system candidates and an additional 14 potential candidates whose combined spectrum mimics the Solar spectrum. Considering the volume correction factor for a magnitude limited survey, the fraction of probable unresolved triple stars with long orbital periods is ~2 %. Possible orbital configurations of the candidates were investigated using the selection and observational limits. To validate the discovered multiplicity fraction, the same procedure was used to evaluate the multiplicity fraction of other stellar types.

Centaurs--icy bodies orbiting beyond Jupiter and interior to Neptune--are believed to be dynamically related to Jupiter Family Comets (JFCs), which have aphelia near Jupiter's orbit and perihelia in the inner Solar System. Previous dynamical simulations have recreated the Centaur/JFC conversion, but the mechanism behind that process remains poorly described. We have performed a numerical simulation of Centaur analogues that recreates this process, generating a dataset detailing over 2.6 million close planet/planetesimal interactions. We explore scenarios stored within that database and, from those, describe the mechanism by which Centaur objects are converted into JFCs. Because many JFCs have perihelia in the terrestrial planet region, and since Centaurs are constantly resupplied from the Scattered Disk and other reservoirs, the JFCs are an ever-present impact threat.

Asteroseismology is a promising tool to study Galactic structure and evolution because it can probe the ages of stars. Earlier attempts comparing seismic data from the {\it Kepler} satellite with predictions from Galaxy models found that the models predicted more low-mass stars compared to the observed distribution of masses. It was unclear if the mismatch was due to inaccuracies in the Galactic models, or the unknown aspects of the selection function of the stars. Using new data from the K2 mission, which has a well-defined selection function, we find that an old metal-poor thick disc, as used in previous Galactic models, is incompatible with the asteroseismic information. We show that spectroscopic measurements of [Fe/H] and [$\alpha$/Fe] elemental abundances from the GALAH survey indicate a mean metallicity of $\log (Z/Z_{\odot})=-0.16$ for the thick disc. Here $Z$ is the effective solar-scaled metallicity, which is a function of [Fe/H] and [$\alpha$/Fe]. With the revised disc metallicities, for the first time, the theoretically predicted distribution of seismic masses show excellent agreement with the observed distribution of masses. This provides an indirect verification of the asteroseismic mass scaling relation is good to within five percent. Using an importance-sampling framework that takes the selection function into account, we fit a population synthesis model of the Galaxy to the observed seismic and spectroscopic data. Assuming the asteroseismic scaling relations are correct, we estimate the mean age of the thick disc to be about 10 Gyr, in agreement with the traditional idea of an old $\alpha$-enhanced thick disc.

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.

Radial velocity (RV) searches for exoplanets have surveyed many of the nearest and brightest stars for long-term velocity variations indicative of a companion body. Such surveys often detect high-amplitude velocity signatures of objects that lie outside the planetary mass regime, most commonly those of a low-mass star. Such stellar companions are frequently discarded as false-alarms to the main science goals of the survey, but high-resolution imaging techniques can be employed to either directly detect or place significant constraints on the nature of the companion object. Here, we present the discovery of a compact companion to the nearby star HD~118475. Our Anglo-Australian Telescope (AAT) RV data allow the extraction of the full Keplerian orbit of the companion, found to have a minimum mass of 0.445~$M_\odot$. Follow-up speckle imaging observations at the predicted time of maximum angular separation rule out a main sequence star as the source of the RV signature at the 3.3$\sigma$ significance level, implying that the companion must be a low-luminosity compact object, most likely a white dwarf. We provide an isochrone analysis combined with our data that constrain the possible inclinations of the binary orbit. We discuss the eccentric orbit of the companion in the context of tidal circularization timescales and show that non-circular orbit was likely inherited from the progenitor. Finally, we emphasize the need for utilizing such an observation method to further understand the demographics of white dwarf companions around nearby stars.

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.

We present an approach that is able to both rapidly assess the dynamical stability of multiple planet systems, and determine whether an exoplanet system would be capable of hosting a dynamically stable Earth-mass companion in its habitable zone. We conduct a suite of numerical simulations using a swarm of massless test particles in the vicinity of the orbit of a massive planet, in order to develop a predictive tool which can be used to achieve these desired outcomes. In this work, we outline both the numerical methods we used to develop the tool, and demonstrate its use. We find that the test particles survive in systems either because they are unperturbed due to being so far removed from the massive planet, or due to being trapped in stable mean motion resonant orbits with the massive planet. The resulting unexcited test particle swarm produces a unique signature in (a,e) space that represents the stable regions within the system. We are able to scale and translate this stability signature, and combine several together in order to conservatively assess the dynamical stability of newly discovered multiple planet systems. We also assess the stability of a system's habitable zone and determine whether an Earth-mass companion could remain on a stable orbit, without the need for exhaustive numerical simulations.

We report the Transiting Exoplanet Survey Satellite ($TESS$) detection of a multi-planet system orbiting the $V=10.9$ K0 dwarf TOI 125. We find evidence for up to five planets, with varying confidence. Three high signal-to-noise transit signals correspond to sub-Neptune-sized planets ($2.76$, $2.79$, and $2.94\ R_{\oplus}$), and we statistically validate the planetary nature of the two inner planets ($P_b = 4.65$ days, $P_c = 9.15$ days). With only two transits observed, we report the outer object ($P_{.03} = 19.98$ days) as a high signal-to-noise ratio planet candidate. We also detect a candidate transiting super-Earth ($1.4\ R_{\oplus}$) with an orbital period of only $12.7$ hours and a candidate Neptune-sized planet ($4.2\ R_{\oplus}$) with a period of $13.28$ days, both at low signal-to-noise. This system is amenable to mass determination via radial velocities and transit timing variations, and provides an opportunity to study planets of similar size while controlling for age and environment. The ratio of orbital periods between TOI 125 b and c ($P_c/P_b = 1.97$) is slightly smaller than an exact 2:1 commensurability and is atypical of multiple planet systems from $Kepler$, which show a preference for period ratios just $wide$ of first-order period ratios. A dynamical analysis refines the allowed parameter space through stability arguments and suggests that, despite the nearly commensurate periods, the system is unlikely to be in resonance.

When, in the course of searching for exoplanets, sparse sampling and noisy data make it necessary to disentangle possible solutions to the observations, one must consider the possibility that what appears to be a single eccentric Keplerian signal may in reality be attributed to two planets in near-circular orbits. There is precedent in the literature for such outcomes, whereby further data or new analysis techniques reveal hitherto occulted signals. Here, we perform suites of simulations to explore the range of possible two-planet configurations that can result in such confusion. We find that a single Keplerian orbit with $e>$0.5 can virtually never be mimicked by such deceptive system architectures. This result adds credibility to the most eccentric planets that have been found to date, and suggests that it could well be worth revisiting the catalogue of moderately eccentric 'confirmed' exoplanets in the coming years, as more data become available, to determine whether any such deceptive couplets are hidden in the observational data.

We examine eight known single-eccentric planetary systems in light of recently released large data archives and new analysis techniques. For four of these systems (HD 7449, HD 65216, HD 89744, HD 92788) we find evidence for additional long-period companions. HD 65216c is a Jupiter analog, with a period of 14.7 yr, $e=0.18$, and m sin $i$ of 2M_Jup, whilst the remaining candidate companions move on as-yet-incomplete orbits. Our results highlight the importance of revisiting the analysis of known exoplanetary systems when new data become available, particularly given the possibility that poorly-sampled data might previously have led to the detection of a 'false-positive' single eccentric planet, when the system in question actually contains two (or more) planets on near-circular orbits.

We present the chemodynamic structure of the solar neighbourhood using 62 814 stars within a 500 pc sphere of the Sun observed by GALAH and with astrometric parameters from Gaia DR2. We measure the velocity dispersion for all three components (vertical, radial, and tangential) and find that it varies smoothly with [Fe/H] and [$\alpha$/Fe] for each component. The vertical component is especially clean, with $\sigma_{v_z}$ increasing from a low of $8$ km s$^{-1}$ at solar-[$\alpha$/Fe] and [Fe/H] to a high of more than 50 km s$^{-1}$ for more metal-poor and [$\alpha$/Fe] enhanced populations. We find no evidence of a large decrease in the velocity dispersion of the higher-[$\alpha$/Fe] populations as claimed in analysis prior to Gaia DR2, although the trend of increasing velocity dispersion with [$\alpha$/Fe] for the same metallicity does significantly flatten at high-[$\alpha$/Fe]. The eccentricity distribution for local stars varies most strongly as a function of [$\alpha$/Fe], where stars with [$\alpha$/Fe]$<0.1$ dex having generally circular orbits ($e<0.15$), while the median eccentricity increases rapidly for more [$\alpha$/Fe]-enhanced stellar populations up to $e\sim0.35$. These [$\alpha$/Fe]-enhanced populations have guiding radii consistent with origins in the inner Galaxy. Of the stars with metallicities much higher than the local ISM ([Fe/H]>0.1 dex), we find that more than 70\% have $e<0.2$ and are likely observed in the solar neighbourhood through churning/migration rather than blurring effects, as the epicyclic motion for these stars is not large enough to reach the radii at which they were likely born based on their metallicity.

Swan bands - characteristic molecular absorption features of the C$_2$ molecule - are a spectroscopic signature of carbon-enhanced stars. They can also be used to identify carbon-enhanced metal-poor (CEMP) stars. The GALAH (GALactic Archaeology with Hermes) is a magnitude-limited survey of stars producing high-resolution, high signal-to-noise spectra. We used 627,708 GALAH spectra to search for carbon-enhanced stars with a supervised and unsupervised classification algorithm, relying on the imprint of the Swan bands. We identified 918 carbon-enhanced stars, including 12 already described in the literature. An unbiased selection function of the GALAH survey allows us to perform a population study of carbon-enhanced stars. Most of them are giants, out of which we find 28 CEMP candidates. A large fraction of our carbon-enhanced stars with repeated observations show variation in radial velocity, hinting that there is a large fraction of variables among them. 32 of the detected stars also show strong Lithium enhancement in their spectra.

Accurate stellar properties are crucial for determining exoplanet characteristics. Gaia DR2 presents revised distances, luminosities, and radii for 1.6 billion stars. Here, we report the calculation of revised radii and densities for 320 exoplanets using this data and present updated calculations of the incident flux received by 690 known exoplanets. This allows the likelihood that those planets orbit in the habitable zone of their host stars to be reassessed. As a result of this analysis, three planets can be added to the catalogue of potentially habitable worlds: HIP~67537~b, HD~148156~b, and HD~106720~b. In addition, the changed parameterisation of BD~+49~898 means that its planet, BD~+49~898~b, now receives an incident flux that places it outside the optimistic habitable zone region, as defined by \citep{Kopparapu2013,Kopparapu2014}. We find that use of the new \textit{Gaia} data results in a mean increase in calculated exoplanet radius of 3.76\%. Previously, CoRoT-3 b had been reported as having the highest density of all known exoplanets. Here, we use updated information to revise the calculated density of CoRoT-3~b from 26.4$g\:cm^{-3}$ to 17.3$\pm2.9g\:cm^{-3}$. We also report the densest exoplanet in our dataset, KELT-1~b, with a density of 23.7$\pm4.0g\:cm^{-3}$. Overall, our results highlight the importance of ensuring the the parameterisation of known exoplanets be revisited whenever significant improvements are made to the precision of the stellar parameters upon which they are based.

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.

If the Galaxy is axisymmetric and in dynamical equilibrium, we expect negligible fluctuations in the residual line-of-sight velocity field. Recent results using the \apg{} survey find significant fluctuations in velocity for stars in the midplane ($|z|<$0.25 kpc) out to 5 kpc, suggesting that the dynamical influence of non-axisymmetric features i.e., the Milky Way's bar, spiral arms and merger events extends out to the Solar neighborhood. Their measured power spectrum has a characteristic amplitude of 11 \kms{} on a scale of 2.5 kpc. The existence of such large-scale streaming motions has important implications for determining the Sun's motion about the Galactic Centre. Using Red Clump stars from \glh{} and \apg{}, we map the line-of-sight velocities around the Sun (d$<$5 kpc), and $|z|<$1.25 kpc from the midplane. By subtracting a smooth axisymmetric model for the velocity field, we study the residual fluctuations and compare our findings with mock survey generated by \glx{}. We find negligible large-scale fluctuations away from the plane. In the mid-plane, we reproduce the earlier \apg{} power spectrum but with 20\% smaller amplitude (9.3 \kms{}) after taking into account a few systematics (e.g., volume completeness). Using a flexible axisymmetric model the power-amplitude is further reduced to 6.3 \kms{}. Additionally, our simulations show that, in the plane, distances are underestimated for high-mass Red Clump stars which can lead to spurious power-amplitude of about 5.2 \kms{}. Taking this into account, we estimate the amplitude of real fluctuations to be $<$4.6 \kms{}, about a factor of three less than the \apg{} result.

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.

In this work, we hunt for the best places to find exo-Earths in the currently known exoplanet population. While it is still unclear whether Jupiter had a beneficial or detrimental effect on the creation of the right environment for a habitable Earth to develop, we focus on the 51 multiple planet systems that have at least one Jupiter-like planet and aim to identify which would be good candidates to host an exo-Earth. We conduct a series of numerical simulations to identify dynamically stable regions of the habitable zone of the multiple exoplanet systems capable of hosting an Earth-mass planet. We produce a candidate list of 16 systems that could host such a stable exo-Earth in their habitable zone, and for which the induced radial velocity signal of a hypothetical one, two or four Earth-mass planet on the host star would be detectable with the ESPRESSO spectrograph. We find that whilst the gravitational interactions with the massive planet nearest the habitable zone are critical in determining stability, the secular resonant interactions between multiple planets can also have a dramatic influence on the overall stability of the habitable zone.

GALAH is a large-scale magnitude-limited southern stellar spectroscopic survey. Its second data release (GALAH DR2) provides values of stellar parameters and abundances of 23 elements for 342,682 stars (Buder et al.). Here we add a description of the public release of radial velocities with a typical accuracy of 0.1 km/s for 336,215 of these stars, achievable due to the large wavelength coverage, high resolving power and good signal to noise ratio of the observed spectra, but also because convective motions in stellar atmosphere and gravitational redshift from the star to the observer are taken into account. In the process we derive medians of observed spectra which are nearly noiseless, as they are obtained from between 100 and 1116 observed spectra belonging to the same bin with a width of 50 K in temperature, 0.2 dex in gravity, and 0.1 dex in metallicity. Publicly released 1181 median spectra have a resolving power of 28,000 and trace the well-populated stellar types with metallicities between -0.6 and +0.3. Note that radial velocities from GALAH are an excellent match to the accuracy of velocity components along the sky plane derived by Gaia for the same stars. The level of accuracy achieved here is adequate for studies of dynamics within stellar clusters, associations and streams in the Galaxy. So it may be relevant for studies of the distribution of dark matter.

Open cluster members are coeval and share the same initial bulk chemical compositions. Consequently, differences in surface abundances between members of a cluster that are at different evolutionary stages can be used to study the effects of mixing and internal chemical processing. We carry out an abundance analysis of seven elements (Li, O, Na, Mg, Al, Si, Fe) in 66 stars belonging to the open cluster M67, based on high resolution GALAH spectra, 1D MARCS model atmospheres, and, for the first time for a large sample of M67 stars, non-local thermodynamic equilibrium (non-LTE) radiative transfer. From the non-LTE analysis, we find a typical star-to-star scatter in the abundance ratios of around 0.05 dex; this scatter is slightly but systematically larger when LTE is assumed instead. We find trends in the abundance ratios with effective temperature, indicating systematic differences in the surface abundances between turn-off and giant stars; these trends are more pronounced when LTE is assumed. However, in the non-LTE analysis, most of the element trends have been flattened. Two species are exceptions to this behaviour, namely Al and Si, which both clearly display remaining trends in the non-LTE analysis. We comment on the possible origin of these trends, by comparing them with recent stellar models that include atomic diffusion.

Sparse open clusters can be found at high galactic latitudes where loosely populated clusters are more easily detected against the lower stellar background. Because most star formation takes place in the thin disk, the observed population of clusters far from the Galactic plane is hard to explain. We combined spectral parameters from the GALAH survey with the Gaia DR2 catalogue to study the dynamics and chemistry of five old sparse high-latitude clusters in more detail. We find that four of them (NGC 1252, NGC 6994, NGC 7772, NGC 7826) - originally classified in 1888 - are not clusters but are instead chance projections on the sky. Member stars quoted in the literature for these four clusters are unrelated in our multi-dimensional physical parameter space; the published cluster properties are therefore irrelevant. We confirm the existence of morphologically similar NGC 1901 for which we provide a probabilistic membership analysis. An overdensity in three spatial dimensions proves to be enough to reliably detect sparse clusters, but the whole 6-dimensional space must be used to identify members with high confidence, as demonstrated in the case of NGC 1901.

Super Earths, planets between 5-10 Earth masses, are the most common type of exoplanet known, yet are completely absent from our Solar system. As a result, their detailed properties, compositions, and formation mechanisms are poorly understood. NASA's Transiting Exoplanet Survey Satellite (TESS) will identify hundreds of Super-Earths orbiting bright stars, for the first time allowing in-depth characterisation of these planets. At the University of Southern Queensland, we are host to the MINERVA-Australis project, dedicated wholly to the follow-up characterisation and mass measurement of TESS planets. We give an update on the status of MINERVA-Australis and our expected performance.

We present a study using the second data release of the GALAH survey of stellar parameters and elemental abundances of 15 pairs of stars identified by Oh et al 2017. They identified these pairs as potentially co-moving pairs using proper motions and parallaxes from Gaia DR1. We find that 11 very wide (>1.7 pc) pairs of stars do in fact have similar Galactic orbits, while a further four claimed co-moving pairs are not truly co-orbiting. Eight of the 11 co-orbiting pairs have reliable stellar parameters and abundances, and we find that three of those are quite similar in their abundance patterns, while five have significant [Fe/H] differences. For the latter, this indicates that they could be co-orbiting because of the general dynamical coldness of the thin disc, or perhaps resonances induced by the Galaxy, rather than a shared formation site. Stars such as these, wide binaries, debris of past star formation episodes, and coincidental co-orbiters, are crucial for exploring the limits of chemical tagging in the Milky Way.

Jupiter and Saturn each have complex systems of satellites and rings. These satellites can be classified into dynamical groups, implying similar formation scenarios. Recently, a larger number of additional irregular satellites have been discovered around both gas giants that have yet to be classified. The aim of this paper is to examine the relationships between the satellites and rings of the gas giants, using an analytical technique called cladistics. Cladistics is traditionally used to examine relationships between living organisms, the `tree of life'. In this work, we perform the first cladistical study of objects in a planetary science context. Our method uses the orbital, physical and compositional characteristics of satellites to classify the objects in the Jovian and Saturnian systems. We find that the major relationships between the satellites in the two systems, such as families, as presented in previous studies, are broadly preserved. In addition, based on our analysis of the Jovian system, we identify a new retrograde irregular family, the Iocaste family, and suggest that the Phoebe family of the Saturnian system can be further divided into two subfamilies. We also propose that the Saturnian irregular families be renamed, to be consistent with the convention used in Jovian families. Using cladistics, we are also able to assign the new unclassified irregular satellites into families. Taken together, the results of this study demonstrate the potential use of the cladistical technique in the investigation of relationships between orbital bodies.

We report the serendipitous observations of 571 luminous supergiants in the Magellanic Clouds by the spectroscopic GALAH and TESS-HERMES surveys: 434 stars in the Large Magellanic Cloud and 137 in the Small Magellanic Cloud. We also find one star that appears associated with structured star formation in the Magellanic Bridge. Both of these surveys are aimed at the local volume of the Galaxy but have simple, magnitude-limited selection functions that mean they include some observations of luminous extra-Galactic stars. The surveys determine stellar parameter and abundances using The Cannon, a data-driven generative modelling approach. In this work, we explore the results from The Cannon when it is fed the spectra of these intrinsically luminous supergiants in the Magellanic Clouds, which are well outside the normal bounds of The Cannon's training set. We find that, although the parameters are astrophysically incorrect, the $v\sin i$ and the abundances of lithium, barium, and magnesium are excellent discriminants of these stars. It shows that in the future, with an expanded training set, it should be possible to determine accurate values for these types of stars.

While the number of exoplanets discovered continues to increase at a rapid rate, we are still to discover any system that truly resembles the Solar system. Existing and near future surveys will likely continue this trend of rapid discovery. To see if these systems are Solar system analogues, we will need to efficiently allocate resources to carry out intensive follow-up observations. We seek to uncover the properties and trends across systems that indicate how much of the habitable zone is stable in each system to provide focus for planet hunters. We study the dynamics of all known single Jovian planetary systems, to assess the dynamical stability of the habitable zone around their host stars. We perform a suite of simulations of all systems where the Jovian planet will interact gravitationally with the habitable zone, and broadly classify these systems. Besides the system's mass ratio (Mpl/Mstar), and the Jovian planet's semi-major axis (apl) and eccentricity (epl), we find that there are no underlying system properties which are observable that indicate the potential for planets to survive within the system's habitable zone. We use Mpl/Mstar, apl and epl to generate a parameter space over which the unstable systems cluster, thus allowing us to predict which systems to exclude from future observational or numerical searches for habitable exoplanets. We also provide a candidate list of 20 systems that have completely stable habitable zones and Jovian planets orbiting beyond the habitable zone as potential first order Solar system analogues.

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.

Using GALAH survey data of nearby stars, we look at how structure in the planar (u,v) velocity distribution depends on metallicity and on viewing direction within the Galaxy. In nearby stars, with distance d < 1 kpc, the Hercules stream is most strongly seen in higher metallicity stars [Fe/H] > 0.2. The Hercules stream peak v value depends on viewed galactic longitude, which we interpret as due to the gap between the stellar stream and more circular orbits being associated with a specific angular momentum value of about 1640 km/s kpc. The association of the gap with a particular angular momentum value supports a bar resonant model for the Hercules stream. Moving groups previously identified in Hipparcos observations are easiest to see in stars nearer than 250 pc, and their visibility and peak velocities in the velocity distributions depends on both viewing direction (galactic longitude and hemisphere) and metallicity. We infer that there is fine structure in local velocity distributions that varies over distances of a few hundred pc in the Galaxy.

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.

Accurate and precise radius estimates of transiting exoplanets are critical for understanding their compositions and formation mechanisms. To know the planet, we must know the host star in as much detail as possible. We present first results from the K2-HERMES project, which uses the HERMES multi-object spectrograph on the Anglo-Australian Telescope to obtain R$\sim$28,000 spectra of up to 360 stars in one exposure. This ongoing project aims to derive self-consistent spectroscopic parameters for about half of K2 target stars. We present complete stellar parameters and isochrone-derived masses and radii for 46 stars hosting 57 K2 candidate planets in Campaigns 1-3. Our revised host-star radii cast severe doubt on three candidate planets: EPIC\,201407812.01, EPIC\,203070421.01, and EPIC\,202843107.01, all of which now have inferred radii well in excess of the largest known inflated Jovian planets.

The surprising discovery of a ring system around the Centaur 10199 Chariklo in 2013 led to a reanalysis of archival stellar occultation data for the Centaur 2060 Chiron by Ortiz et al. One possible interpretation of that data is that a system of rings exists around Chiron. In this work, we study the dynamical history of the proposed Chiron ring system by integrating nearly 36,000 clones of the Centaur backward in time for 100 Myr under the influence of the Sun and the four giant planets. The severity of all close encounters between the clones and planets while the clones are in the Centaur region is recorded, along with the mean time between close encounters. We find that severe and extreme close encounters are very rare, making it possible that the Chiron ring system has remained intact since its injection into the Centaur region, which we find likely occurred within the past 8.5 Myr. Our simulations yield a backward dynamical half-life for Chiron of 0.7 Myr. The dynamical classes of a sample of clones are found. It is found that, on average, the Centaur lifetimes of resonance hopping clones are twice those of random-walk clones because of resonance sticking in mean motion resonances. In addition, we present MEGNO and chaotic lifetime maps of the region bound by 13 au $\le a \le$14 au and e $\le$0.5. We confirm that the current mean orbital parameters of Chiron are located in a highly chaotic region of a-e phase space.

Radial velocity observations from three instruments reveal the presence of a 4 M_jup planet candidate orbiting the K giant HD 76920. HD 76920b has an orbital eccentricity of 0.856$\pm$0.009, making it the most eccentric planet known to orbit an evolved star. There is no indication that HD 76920 has an unseen binary companion, suggesting a scattering event rather than Kozai oscillations as a probable culprit for the observed eccentricity. The candidate planet currently approaches to about four stellar radii from its host star, and is predicted to be engulfed on a $\sim$100 Myr timescale due to the combined effects of stellar evolution and tidal interactions.

The technique of chemical tagging uses the elemental abundances of stellar atmospheres to `reconstruct' chemically homogeneous star clusters that have long since dispersed. The GALAH spectroscopic survey --which aims to observe one million stars using the Anglo-Australian Telescope -- allows us to measure up to 30 elements or dimensions in the stellar chemical abundance space, many of which are not independent. How to find clustering reliably in a noisy high-dimensional space is a difficult problem that remains largely unsolved. Here we explore t-distributed stochastic neighbour embedding (t-SNE) -- which identifies an optimal mapping of a high-dimensional space into fewer dimensions -- whilst conserving the original clustering information. Typically, the projection is made to a 2D space to aid recognition of clusters by eye. We show that this method is a reliable tool for chemical tagging because it can: (i) resolve clustering in chemical space alone, (ii) recover known open and globular clusters with high efficiency and low contamination, and (iii) relate field stars to known clusters. t-SNE also provides a useful visualization of a high-dimensional space. We demonstrate the method on a dataset of 13 abundances measured in the spectra of 187,000 stars by the GALAH survey. We recover 7 of the 9 observed clusters (6 globular and 3 open clusters) in chemical space with minimal contamination from field stars and low numbers of outliers. With chemical tagging, we also identify two Pleiades supercluster members (which we confirm kinematically), one as far as 6$^\circ$ -- one tidal radius away from the cluster centre.

In recent years, there has been significant debate and discussion about the glaring gender disparity in the physical sciences. To better understand and address this within the Australian Space Research Community, in 2015 we began the process of keeping a record of the gender balance at the annual Australian Space Research Conference. In addition, we began holding an annual 'Women in Space Research' lunch at that meeting, to discuss the situation, and search for routes by which issues of equity can be addressed, and the situation improved. We present an update based on the 16th Australian Space Research Conference, held at RMIT, Melbourne, in September 2016. As in 2015, male attendees outnumbered female attendees approximately 3:1. However, there was a small shift (~2.3%) in the balance, with female delegates now making up 26.4% of the total, up from 24.1% in 2015. This shift was mirrored in the gender distribution of talks, with 28.9% of all oral presentations being given by women (up from 25.2%). More striking, however, were the changes in the distribution of plenary presentations (44.4% female, up from 22.2%), poster presentations (31.8% female, up from 7.7%), and student awards (33.3% female, up from 12.5%). These changes are encouraging, and will hopefully continue in years to come. The conference organising committee again mirrored the gender balance of the delegates as a whole (27.3 % female vs. 26.4% of delegates), though the program committee was markedly more male-dominated this year than last (82.4% male, against last year's 72.2%). At this year's 'Women in Space Research' lunchtime event, a number of suggestions were made that could help to make future conferences, and the wider community, a more equitable place... --- abstract abridged ---

In the coming years, it is likely that the first potentially Earth-like planets will be discovered orbiting other stars. Once found, the characterisation of those planets will play a vital role in determining which will be chosen as the first targets for the search for life beyond the Solar System. We must thus be able to gauge the relative importance of the various factors proposed to influence potential planetary habitability, in order to best focus that search. One of the plethora of factors to be considered in that process is the climatic variability of the exo-Earths in question. In the 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 the Solar System adds 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 preliminary results of three suites of integrations that together examine the influence of Solar System architecture on the Earth's Milankovitch cycles. We consider separately the influence of the planets Jupiter, Mars and Venus, each of which contributes to the forcing of Earth's orbital evolution. Our results illustrate how small changes to the architecture of a given planetary system can result in marked changes in the potential habitability of the planets therein, and are an important first step in developing a means by which the nature of climate variability on planets beyond our Solar System can be characterised.

With continued improvement in telescope sensitivity and observational techniques, the search for rocky planets in stellar habitable zones is entering an exciting era. With so many exoplanetary systems available for follow-up observations to find potentially habitable planets, one needs to prioritise the ever-growing list of candidates. We aim to determine which of the known planetary systems are dynamically capable of hosting rocky planets in their habitable zones, with the goal of helping to focus future planet search programs. We perform an extensive suite of numerical simulations to identify regions in the habitable zones of single Jovian planet systems where Earth mass planets could maintain stable orbits, specifically focusing on the systems in the Catalog of Earth-like Exoplanet Survey Targets (CELESTA). We find that small, Earth-mass planets can maintain stable orbits in cases where the habitable zone is largely, or partially, unperturbed by a nearby Jovian, and that mutual gravitational interactions and resonant mechanisms are capable of producing stable orbits even in habitable zones that are significantly or completely disrupted by a Jovian. Our results yield a list of 13 single Jovian planet systems in CELESTA that are not only capable of supporting an Earth-mass planet on stable orbits in their habitable zone, but for which we are also able to constrain the orbits of the Earth-mass planet such that the induced radial velocity signals would be detectable with next generation instruments.

Chariklo is the only small Solar system body confirmed to have rings. Given the instability of its orbit, the presence of rings is surprising, and their origin remains poorly understood. In this work, we study the dynamical history of the Chariklo system by integrating almost 36,000 Chariklo clones backwards in time for one Gyr under the influence of the Sun and the four giant planets. By recording all close encounters between the clones and planets, we investigate the likelihood that Chariklo's rings could have survived since its capture to the Centaur population. Our results reveal that Chariklo's orbit occupies a region of stable chaos, resulting in its orbit being marginally more stable than those of the other Centaurs. Despite this, we find that it was most likely captured to the Centaur population within the last 20 Myr, and that its orbital evolution has been continually punctuated by regular close encounters with the giant planets. The great majority (> 99%) of those encounters within one Hill radius of the planet have only a small effect on the rings. We conclude that close encounters with giant planets have not had a significant effect on the ring structure. Encounters within the Roche limit of the giant planets are rare, making ring creation through tidal disruption unlikely.

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.

We present the discovery of two extended $\sim$0.12 mag dimming events of the weak-lined T-Tauri star V1334. The start of the first event was missed but came to an end in late 2003, and the second began in February 2009, and continues as of November 2016. Since the egress of the current event has not yet been observed, it suggests a period of $>$13 years if this event is periodic. Spectroscopic observations suggest the presence of a small inner disk, although the spectral energy distribution shows no infrared excess. We explore the possibility that the dimming events are caused by an orbiting body (e.g. a disk warp or dust trap), enhanced disk winds, hydrodynamical fluctuations of the inner disk, or a significant increase in the magnetic field flux at the surface of the star. We also find a $\sim$0.32 day periodic photometric signal that persists throughout the 2009 dimming which appears to not be due to ellipsoidal variations from a close stellar companion. High precision photometric observations of V1334 Tau during K2 campaign 13, combined with simultaneous photometric and spectroscopic observations from the ground, will provide crucial information about the photometric variability and its origin.

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 report the discovery of a second long-period giant planet orbiting HD 30177, a star previously known to host a massive Jupiter analog (HD 30177b: a=3.8$\pm$0.1 au, m sin $i=9.7\pm$0.5 Mjup). HD 30177c can be regarded as a massive Saturn analog in this system, with a=9.9$\pm$1.0 au and m sin $i=7.6\pm$3.1 Mjup. The formal best fit solution slightly favours a closer-in planet at $a\sim$7 au, but detailed n-body dynamical simulations show that configuration to be unstable. A shallow local minimum of longer-period, lower-eccentricity solutions was found to be dynamically stable, and hence we adopt the longer period in this work. The proposed $\sim$32 year orbit remains incomplete; further monitoring of this and other stars is necessary to reveal the population of distant gas giant planets with orbital separations $a\sim$10 au, analogous to that of Saturn.

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.

The Galactic Archaeology with HERMES (GALAH) Survey is a massive observational project to trace the Milky Way's history of star formation, chemical enrichment, stellar migration and minor mergers. Using high-resolution (R$\simeq$28,000) spectra taken with the High Efficiency and Resolution Multi-Element Spectrograph (HERMES) instrument at the Anglo-Australian Telescope (AAT), GALAH will determine stellar parameters and abundances of up to 29 elements for up to one million stars. Selecting targets from a colour-unbiased catalogue built from 2MASS, APASS and UCAC4 data, we expect to observe dwarfs at 0.3 to 3 kpc and giants at 1 to 10 kpc. This enables a thorough local chemical inventory of the Galactic thin and thick disks, and also captures smaller samples of the bulge and halo. In this paper we present the plan, process and progress as of early 2016 for GALAH survey observations. In our first two years of survey observing we have accumulated the largest high-quality spectroscopic data set at this resolution, over 200,000 stars. We also present the first public GALAH data catalogue: stellar parameters (Teff, log(g), [Fe/H], [alpha/Fe]), radial velocity, distance modulus and reddening for 10680 observations of 9860 Tycho-2 stars that may be included in the first Gaia data release.

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.

The Kilodegree Extremely Little Telescope (KELT) survey is a ground-based program designed to search for transiting exoplanets orbiting relatively bright stars. To achieve this, the KELT Science Team operates two planets facilities - KELT-North, at Winer Observatory, Arizona, and KELT-South, at the South African Astronomical Observatory. The telescopes used at these observatories have particularly wide fields of view, allowing KELT to study a large number of potential exoplanet host stars. One of the major advantages of targeting bright stars is that the exoplanet candidates detected can be easily followed up by small, ground-based observatories distributed around the world. This paper will provide a brief overview of the KELT-North and the KELT-South surveys, the follow-up observations preformed by the KELT Follow-up Collaboration, and exoplanet discoveries confirmed thus far, before concluding with a brief discussion of the future for the KELT program.

In recent years, the striking gender imbalance in the physical sciences has been a topic for much debate. National bodies and professional societies in the astronomical and space sciences are now taking active steps to understand and address this imbalance. In order to begin this process in the Australian Space Research community, we must first understand the current state of play. In this work, we therefore present a short 'snapshot' of the current gender balance in our community, as observed at the 15th Australian Space Research Conference. We find that, at this year's conference, male attendees outnumbered female attendees by a ratio of 3:1 (24% female). This gender balance was repeated in the distribution of conference talks and plenary presentations (25 and 22% female, respectively). Of the thirteen posters presented at the conference, twelve were presented by men (92%), a pattern repeated in the awards for the best student presentations (seven male recipients vs one female). The program and organising committees for the meeting fairly represented the gender balance of the conference attendees (28% and 30% female, respectively). These figures provide a baseline for monitoring future progress in increasing the participation of women in the field. They also suggest that the real barrier is not speaking, but in enabling conference attendance and retaining female scientists through their careers - in other words, addressing and repairing the 'leaky pipeline'.

We present updated simulations of the detectability of Jupiter analogs by the 17-year Anglo-Australian Planet Search. The occurrence rate of Jupiter-like planets that have remained near their formation locations beyond the ice line is a critical datum necessary to constrain the details of planet formation. It is also vital in our quest to fully understand how common (or rare) planetary systems like our own are in the Galaxy. From a sample of 202 solar-type stars, and correcting for imperfect detectability on a star-by-star basis, we derive a frequency of $6.2^{+2.8}_{-1.6}$% for giant planets in orbits from 3-7 AU. When a consistent definition of "Jupiter analog" is used, our results are in agreement with those from other legacy radial velocity surveys.

The Neptune Trojans are the most recent addition to the panoply of Solar system small body populations. The orbit of the first discovered member, 2001 QR$_{322}$, was investigated shortly after its discovery, based on early observations of the object, and it was found to be dynamically stable on timescales comparable to the age of the Solar system. As further observations were obtained of the object over the following years, the best-fit solution for its orbit changed. We therefore carried out a new study of 2001 QR$_{322}$'s orbit in 2010, finding that it lay on the boundary between dynamically stable and unstable regions in Neptune's Trojan cloud, and concluding that further observations were needed to determine the true stability of the object's orbit. Here we follow up on that earlier work, and present the preliminary results of a dynamical study using an updated fit to 2001 QR$_{322}$'s orbit. Despite the improved precision with which the orbit of 2001 QR$_{322}$ is known, we find that the best-fit solution remains balanced on a knife-edge, lying between the same regions of stability and instability noted in our earlier work. In the future, we intend to carry out new observations that should hopefully refine the orbit to an extent that its true nature can finally be disentangled.

We report the discovery of two giant planets orbiting the K giant HD 33844 based on radial velocity data from three independent campaigns. The planets move on nearly circular orbits with semimajor axes $a_b=1.60\pm$0.02 AU and $a_c=2.24\pm$0.05 AU, and have minimum masses (m sin $i$) of $M_b=1.96\pm$0.12 Mjup and $M_c=1.76\pm$0.18 Mjup. Detailed N-body dynamical simulations show that the two planets remain on stable orbits for more than $10^6$ years for low eccentricities, and are most likely trapped in a mutual 3:5 mean-motion resonance.

We report the detection of two new long-period giant planets orbiting the stars HD 95872 and HD 162004 (psi1 Draconis B) by the McDonald Observatory planet search. The planet HD 95872b has a minimum mass of 4.6 M_Jup and an orbital semi-major axis of 5.2 AU. The giant planet psi1 Dra Bb has a minimum mass of 1.5 M_Jup and an orbital semi-major axis of 4.4 AU. Both of these planets qualify as Jupiter analogs. These results are based on over one and a half decades of precise radial velocity measurements collected by our program using the McDonald Observatory Tull Coude spectrograph at the 2.7 m Harlan J. Smith telescope. In the case of psi1 Draconis B we also detect a long-term non-linear trend in our data that indicates the presence of an additional giant planet, similar to the Jupiter-Saturn pair. The primary of the binary star system, psi1 Dra A, exhibits a very large amplitude radial velocity variation due to another stellar companion. We detect this additional member using speckle imaging. We also report two cases - HD 10086 and HD 102870 (beta Virginis) - of significant radial velocity variation consistent with the presence of a planet, but that are probably caused by stellar activity, rather than reflexive Keplerian motion. These two cases stress the importance of monitoring the magnetic activity level of a target star, as long-term activity cycles can mimic the presence of a Jupiter-analog planet.

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.

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.

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.

The bispectrum of single-field inflationary trajectories in which the speed of sound of the inflationary trajectories $c_s$ is constant but not equal to the speed of light $c=1$ is explored. The trajectories are generated as random realisations of the Hubble Slow-Roll (HSR) hierarchy and the bispectra are calculated using numerical techniques that extends previous work. This method allows for out-of-slow-roll models with non-trivial time dependence and arbitrarily low $c_s$. The ensembles obtained using this method yield distributions for the shape and scale-dependence of the bispectrum and their relations with the standard inflationary parameters such as scalar spectral tilt $n_s$ and tensor-to-scalar ratio $r$. The distributions demonstrate the squeezed-limit consistency relations for arbitrary single-field inflationary models.

In this work we present the results of a suite of dynamical simulations following the orbital evolution of 8,022 hypothetical Centaur objects. These Centaurs begin our integrations on orbits in the vicinity of the 2:1 mean motion resonance with Neptune, and we follow their dynamical evolution for a period of 3 Myr under the gravitational influence of a motionless Sun and the four Jovian planets. The great majority of the test particles studied rapidly escaped from the vicinity of the 2:1 mean motion resonance of Neptune and diffused throughout the Solar System. The average libration time of Centaurs in the vicinity of 2:1 mean motion resonance of Neptune was found to be just 27 kyr. Although two particles did remain near the resonance for more than 1 Myr. Upon leaving the vicinity of the 2:1 resonance, the majority of test particles evolved by a process of random walk in semi-major axis, due to repeated close encounters with the giant planets.

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.

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.

The HR 8799 system, with its four giant planets and two debris belts, has an architecture closely mirroring that of our Solar system where the inner, warm asteroid belt and outer, cool Edgeworth-Kuiper belt bracket the giant planets. As such, it is a valuable laboratory for examining exoplanetary dynamics and debris disk-exoplanet interactions. Whilst the outer debris belt of HR 8799 has been well resolved by previous observations, the spatial extent of the inner disk remains unknown. This leaves a significant question mark over both the location of the planetesimals responsible for producing the belt's visible dust and the physical properties of those grains. We have performed the most extensive simulations to date of the inner, unresolved debris belt around HR 8799, using UNSW Australia's Katana supercomputing facility to follow the dynamical evolution of a model inner disk comprising 300,298 particles for a period of 60 million years. These simulations have enabled the characterisation of the extent and structure of the inner disk in detail, and will in future allow us to provide a first estimate of the small-body impact rate and water delivery prospects for possible (as-yet undetected) terrestrial planet(s) in the inner system.

Precise radial velocities from the Anglo-Australian Telescope confirm the presence of a rare short-period planet around the K0 giant HD 121056. An independent two-planet solution using the AAT data shows that the inner planet has P=89.1+/-0.1 days, and m sin i=1.35+/-0.17 Mjup. These data also confirm the planetary nature of the outer companion, with m sin i=3.9+/-0.6 Mjup and a=2.96+/-0.16 AU. HD 121056 is the most-evolved star to host a confirmed multiple-planet system, and is a valuable example of a giant star hosting both a short-period and a long-period planet.

We have investigated the dynamical stability of the proposed companions orbiting the Algol type short-period eclipsing binary SW Lyncis (Kim et al. 2010). The two candidate companions are of stellar to sub-stellar nature, and were inferred from timing measurements of the system's primary and secondary eclipses. We applied well-tested numerical techniques to accurately integrate the orbits of the two companions and to test for chaotic dynamical behaviour. We carried out the stability analysis within a systematic parameter survey varying both the geometries and orientation of the orbits of the companions, as well as their masses. In all our numerical integrations we found that the proposed SW Lyn multi-body system is highly unstable on time-scales on the order of 1000 years. Our results cast doubt on the interpretation that the timing variations are caused by two companions. This work demonstrates that a straightforward dynamical analysis can help to test whether a best-fit companion-based model is a physically viable explanation for measured eclipse timing variations. We conclude that dynamical considerations reveal that the propsed SW Lyncis multi-body system most likely does not exist or the companions have significantly different orbital properties as conjectured in Kim et al. (2010).

The simplest interpretation of the Bicep2 result is that the scalar primordial power spectrum is slightly suppressed at large scales. These models result in a large tensor-to-scalar ratio $r$. In this work we show that the type of inflationary trajectory favoured by Bicep2 also leads to a larger non-Gaussian signal at large scales, roughly an order of magnitude larger than a standard slow-roll trajectory.

We report the detection of GJ 832c, a super-Earth orbiting near the inner edge of the habitable zone of GJ 832, an M dwarf previously known to host a Jupiter analog in a nearly-circular 9.4-year orbit. The combination of precise radial-velocity measurements from three telescopes reveals the presence of a planet with a period of 35.68+/-0.03 days and minimum mass (m sin i) of 5.4+/-1.0 Earth masses. GJ 832c moves on a low-eccentricity orbit (e=0.18+/-0.13) towards the inner edge of the habitable zone. However, given the large mass of the planet, it seems likely that it would possess a massive atmosphere, which may well render the planet inhospitable. Indeed, it is perhaps more likely that GJ 832c is a "super-Venus," featuring significant greenhouse forcing. With an outer giant planet and an interior, potentially rocky planet, the GJ 832 planetary system can be thought of as a miniature version of our own Solar system.

In this work we revisit the proposed multi-circumbinary system RZ Dra. We find the proposed system to be highly unstable. We attempt to find a best-fit light-travel time model rendering the orbits to follow stable orbits. We found a best-fit solution, but the resulting orbits also exhibit short-term orbital instabilities. We therefore conclude that the observed timing variation must be due to other astrophysical effects causing a change in the binary orbital period. Follow-up observations of this system is encouraged.

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)

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.

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.

In coming years, the first truly Earth-like planets will be discovered orbiting other stars, and the search for signs of life on these worlds will begin. However, such observations will be hugely time-consuming and costly, and so it will be important to determine which of those planets represent the best prospects for life elsewhere. One of the key factors in such a decision will be the climate variability of the planet in question - too chaotic a climate might render a planet less promising as a target for our initial search for life elsewhere. On the Earth, the climate of the last few million years has been dominated by a series of glacial and interglacial periods, driven by periodic variations in the Earth's orbital elements and axial tilt. These Milankovitch cycles are driven by the gravitational influence of the other planets, and as such are strongly dependent on the architecture of the Solar system. Here, we present the first results of a study investigating the influence of the orbit of Jupiter on the Milankovitch cycles at Earth - a first step in developing a means to characterise the nature of periodic climate change on planets beyond our Solar system.

In the last few years, a number of planets have been proposed to orbit several post main-sequence binary star systems on the basis of observed variations in the timing of eclipses between the binary components. A common feature of these planet candidates is that the best-fit orbits are often highly eccentric, such that the multiple planet systems proposed regularly feature mutually crossing orbits - a scenario that almost always leads to unstable planetary systems. In this work, we present the results of dynamical studies of all multiple-planet systems proposed to orbit these highly evolved binary stars, finding that most do not stand up to dynamical scrutiny. In one of the potentially stable cases (the NN Serpentis 2-planet system), we consider the evolution of the binary star system, and show that it is highly unlikely that planets could survive from the main sequence to obtain their current orbits - again casting doubt on the proposed planets. We conclude by considering alternative explanations for the observed variation in eclipse timings for these systems.

We report the discovery of two long-period giant planets from the Anglo-Australian Planet Search. HD 154857c is in a multiple-planet system, while HD 114613b appears to be solitary. HD 114613b has an orbital period P=10.5 years, and a minimum mass m sin i of 0.48 Jupiter masses; HD 154857c has P=9.5 years and m sin i=2.6 Jupiter masses. These new data confirm the planetary nature of the previously unconstrained long-period object in the HD 154857 system. We have performed detailed dynamical stability simulations which show that the HD 154857 two-planet system is stable on timescales of at least 100 million years. These results highlight the continued importance of "legacy" surveys with long observational baselines; these ongoing campaigns are critical for determining the population of Jupiter analogs, and hence of those planetary systems with architectures most like our own Solar system.

We present a detailed dynamical analysis of the orbital stability of the BD +20 2457 system, which features planets or brown dwarfs moving on relatively eccentric orbits. We find that the system exhibits strong dynamical instability on astronomically short timescales across a wide range of plausible orbital eccentricities, semi-major axes, and inclinations. If the system truly hosts massive planets or brown dwarfs, our results suggest that they must move on orbits significantly different to those proposed in the discovery work. If that is indeed the case, then it is likely that the best-fit orbital solutions for the proposed companions will change markedly as future observations are made. Such observations may result in the solution shifting to a more dynamically-stable regime, potentially one where stability is ensured by mutually resonant motion.

We use the Hamilton--Jacobi formalism to constrain the space of possible single field, inflationary Hubble flow trajectories when compared to the WMAP and Planck satellites Cosmic Microwave Background (CMB) results. This method yields posteriors on the space of Hubble Slow Roll (HSR) parameters that uniquely determine the history of the Hubble parameter during the inflating epoch. The trajectories are used to numerically determine the observable primordial power spectrum and bispectra that can then be compared to observations. Our analysis is used to infer the most likely shape of the inflaton potential $V(\phi)$ and also yields a prediction for, $f_{\rm NL}$, the dimensionless amplitude of the non-Gaussian bispectrum.

We present six years of new radial-velocity data from the Anglo-Australian and Magellan Telescopes on the HD 73526 2:1 resonant planetary system. We investigate both Keplerian and dynamical (interacting) fits to these data, yielding four possible configurations for the system. The new data now show that both resonance angles are librating, with amplitudes of 40 degrees and 60 degrees, respectively. We then perform long-term dynamical stability tests to differentiate these solutions, which only differ significantly in the masses of the planets. We show that while there is no clearly preferred system inclination, the dynamical fit with i=90 degrees provides the best combination of goodness-of-fit and long-term dynamical stability.

In this work we carry out an analysis of the observed times of primary and secondary eclipses of the post-common envelope binary NSVS14256825. Recently, \cite{Almeida2013} proposed that two circumbinary companions orbit this short-period eclipsing binary, in order to explain observed variations in the timing of mutual eclipses between the two binary components. Using a standard weighted least-squares minimisation technique, we have extensively explored the topology of $\chi^2$ parameter space of a single planet model. We find the data set to be insufficient to reliably constrain a one-companion model. Various models, each with similar statistical significance, result in substantially different orbital architectures for the additional companion. No evidence is seen for a second companion of planetary nature. We suspect insufficient coverage (baseline) of timing data causing the best-fit parameters to be unconstrained.

We carry out a numerical calculation of the bispectrum in generalised trajectories of canonical, single--field inflation. The trajectories are generated in the Hamilton-Jacobi (HJ) formalism based on Hubble Slow Roll (HSR) parameters. The calculation allows generally shape and scale dependent bispectra, or dimensionless $f_{NL}$, in the out-of-slow-roll regime. The distributions of $f_{NL}$ for various shapes and HSR proposals are shown as an example of how this procedure can be used within the context of Monte Carlo exploration of inflationary trajectories. We also show how allowing out-of-slow-roll behaviour can lead to a bispectrum that is relatively large for equilateral shapes.

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.

We use the most up to date Milky Way model and solar orbit data in order to test the hypothesis that the Sun's galactic spiral arm crossings cause mass extinction events on Earth. To do this, we created a new model of the Milky Way's spiral arms by combining a large quantity of data from several surveys. We then combined this model with a recently derived solution for the solar orbit to determine the timing of the Sun's historical passages through the Galaxy's spiral arms. Our new model was designed with a symmetrical appearance, with the major alteration being the addition of a spur at the far side of the Galaxy. A correlation was found between the times at which the Sun crosses the spiral arms and six known mass extinction events. Furthermore, we identify five additional historical mass extinction events that might be explained by the motion of the Sun around our Galaxy. These five additional significant drops in marine genera that we find include significant reductions in diversity at 415, 322, 300, 145 and 33 Myr ago. Our simulations indicate that the Sun has spent ~60% of its time passing through our Galaxy's various spiral arms. Also, we briefly discuss and combine previous work on the Galactic Habitable Zone with the new Milky Way model.

Recent observations of the NN Serpentis post-common envelope binary system have revealed eclipse timing variations that have been attributed to the presence of two Jovian-mass exoplanets. Under the assumption that these planets are real and survived from the binary's Main Sequence state, we reconstruct initial binaries that give rise to the present NN Ser configuration and test the dynamical stability of the original system. Under standard assumptions about binary evolution, we find that survival of the planets through the entire Main Sequence life-time is very unlikely. Hence, we conclude that the planets are not survivors from before the Common Envelope phase, implying that either they formed recently out of material ejected from the primary, or that the observed signals are of non-planetary origin.

Dynamical studies of new exoplanet systems are a critical component of the discovery and characterisation process. Such studies can provide firmer constraints on the parameters of the newly discovered planets, and may even reveal that the proposed planets do not stand up to dynamical scrutiny. Here, we demonstrate how dynamical studies can assist the characterisation of such systems through two examples: QS Virginis and HD 73526.

In recent years, a number of planetary systems have been proposed to orbit evolved binary star systems. The presence of planets is invoked to explain observed variations in the timing of mutual eclipses between the primary and secondary components of the binary star system. The planets recently proposed orbiting the cataclysmic variable system QS Virginis are the latest in this on-going series of "extreme planets". The two planets proposed to orbit QS Virginis would move on mutually crossing orbits - a situation that is almost invariably unstable on very short timescales. In this work, we present the results of a detailed dynamical study of the orbital evolution of the two proposed planets, revealing that they are dynamically unstable on timescales of less than one thousand years across the entire range of orbital elements that provide a plausible fit to the observational data, and regardless of their mutual orbital inclination. We conclude that the proposed planets around the cataclysmic variable QS Virginis simply cannot exist.

Determining the orbital eccentricity of an extrasolar planet is critically important for understanding the system's dynamical environment and history. However, eccentricity is often poorly determined or entirely mischaracterized due to poor observational sampling, low signal-to-noise, and/or degeneracies with other planetary signals. Some systems previously thought to contain a single, moderate-eccentricity planet have been shown, after further monitoring, to host two planets on nearly-circular orbits. We investigate published apparent single-planet systems to see if the available data can be better fit by two lower-eccentricity planets. We identify nine promising candidate systems and perform detailed dynamical tests to confirm the stability of the potential new multiple-planet systems. Finally, we compare the expected orbits of the single- and double-planet scenarios to better inform future observations of these interesting systems.

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.

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.

The Hamilton-Jacobi (HJ) approach for exploring inflationary trajectories is employed in the generation of generalised inflationary non-Gaussian signals arising from single field inflation. Scale dependent solutions for $f_{NL}$ are determined via the numerical integration of the three--point function in the curvature perturbation. This allows the full exploration of single field inflationary dynamics in the out-of-slow-roll regime and opens up the possibility of using future observations of non-Gaussianity to constraint the inflationary potential using model-independent methods. The distribution of `equilateral' $f_{NL}$ arising from single field inflation with both canonical and non-canonical kinetic terms are show as an example of the application of this procedure.

As a direct result of ongoing efforts to detect more exoplanetary systems, an ever-increasing number of multiple-planet systems are being announced. But how many of these systems are truly what they seem? In many cases, such systems are announced solely on the basis of orbital fits to observational data, and no attempt is made to see whether the proposed orbits are actually dynamically feasible. As a result, it is certain that planetary systems are being announced that involve planets moving on orbits that would be dynamically unstable on timescales of just a few hundred years. Here, we present the results of dynamical simulations that investigate the orbital stability and evolution of a number of recently discovered exoplanetary systems. These simulations have enabled us to create highly detailed dynamical maps of those systems, allowing us to better constrain the orbits of the planets contained therein. In some cases, our results have even led to the very existence of the planets themselves being called into question.

We present a detailed dynamical analysis of the orbital stability of the two circumbinary planets recently proposed to orbit the evolved eclipsing binary star system NSVS 14256825. As is the case for other recently proposed circumbinary planetary systems detected through the timing of mutual eclipses between the central binary stars, the proposed planets do not stand up to dynamical scrutiny. The proposed orbits for the two planets are extremely unstable on timescales of less than a thousand years, regardless of the mutual inclination between the planetary orbits. For the scenario where the planetary orbits are coplanar, a small region of moderate stability was observed, featuring orbits that were somewhat protected from destabilisation by the influence of mutual 2:1 mean-motion resonance between the orbits of the planets. Even in this stable region, however, the systems tested typically only survived on timescales of order 1 million years, far shorter than the age of the system. Our results suggest that, if there are planets in the NSVS 14256825 system, they must move on orbits dramatically different to those proposed in the discovery work. More observations are clearly critically required in order to constrain the nature of the suggested orbital bodies.

The details of the Solar system's formation are still heavily debated. Questions remain about the formation locations of the giant planets, and the degree to which volatile material was mixed throughout the proto-planetary system. One diagnostic which offers great promise in helping to unravel the history of planet formation is the study of the level of deuteration in various Solar system bodies. For example, the D/H ratio of methane in the atmosphere of Titan can be used as a diagnostic of the initial conditions of the solar nebula within the region of giant planet formation, and can help us to determine where Titan (and, by extension, the Saturnian system) accreted its volatile material. The level of Titanian deuteration also has implications for both the sources and long term evolution of Titan's atmospheric composition. We present the results of observations taken in the 1.58 microns window using the NIFS spectrometer on the Gemini telescope, and model our data using the VSTAR line--by--line transfer model, which yields a D/H ratio for Titan's atmosphere of 143+/-16) x 10^{-6} [1]. We are currently in the process of modeling the Gemini high resolution GNIRS spectra using new sets of line parameters derived for methane in the region between 1.2-1.7 microns [2].

It is now accepted that the Solar system's youth was a dynamic and chaotic time. The giant planets migrated significant distances to reach their current locations, and evidence of that migration's influence on the Solar system abounds. That migration's pace, and the distance over which it occurred, is still heavily debated. Some models feature systems in which the giant planets were initially in an extremely compact configuration, in which Uranus and Neptune are chaotically scattered into the outer Solar system. Others feature architectures that were initially more relaxed, and smoother, more sedate migration. To determine which of these scenarios best represents the formation of our Solar system, we must turn to the structure of the system's small body populations, in which the scars of that migration are still clearly visible. We present the first results of a program investigating the effect of giant planet migration on the reservoirs of small bodies that existed at that time. As the planets migrate, they stir these reservoirs, scattering vast numbers of small bodies onto dynamically unstable orbits in the outer Solar system. The great majority of those bodies are rapidly removed from the system, through collisions and ejections, but a small number become captured as planetary Trojans or irregular satellites. Others are driven by the migration, leading to a significant sculpting of the asteroid belt and trans-Neptunian region. The capture and retention efficiencies to these stable reservoirs depend on the particular migration scenario used. Advocates of chaotic migration from an initially compact scenario argue that smoother, more sedate migration cannot explain the observed populations of Trojans and irregular satellites. Our results draw a strikingly different picture, revealing that such smooth migration is perfectly capable of reproducing the observed populations.

The abilities of radial velocity exoplanet surveys to detect the lowest-mass extra-solar planets are currently limited by a combination of instrument precision, lack of data, and "jitter". Jitter is a general term for any unknown features in the noise, and reflects a lack of detailed knowledge of stellar physics (asteroseismology, starspots, magnetic cycles, granulation, and other stellar surface phenomena), as well as the possible underestimation of instrument noise. We study an extensive set of radial velocities for the star HD 10700 ($\tau$ Ceti) to determine the properties of the jitter arising from stellar surface inhomogeneities, activity, and telescope-instrument systems, and perform a comprehensive search for planetary signals in the radial velocities. We perform Bayesian comparisons of statistical models describing the radial velocity data to quantify the number of significant signals and the magnitude and properties of the excess noise in the data. We reach our goal by adding artificial signals to the "flat" radial velocity data of HD 10700 and by seeing which one of our statistical noise models receives the greatest posterior probabilities while still being able to extract the artificial signals correctly from the data. We utilise various noise components to assess properties of the noise in the data and analyse the HARPS, AAPS, and HIRES data for HD 10700 to quantify these properties and search for previously unknown low-amplitude Keplerian signals. ...

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.

We perform several suites of highly detailed dynamical simulations to investigate the architectures of the 24 Sextantis and HD 200964 planetary systems. The best fit orbital solution for the two planets in the 24 Sex system places them on orbits with periods that lie very close to 2:1 commensurability, while that for the HD 200964 system places the two planets therein in orbits whose periods lie close to a 4:3 commensurability. In both cases, the proposed best-fit orbits are mutually crossing - a scenario that is only dynamically feasible if the planets are protected from close encounters by the effects of mutual mean motion resonance. Our simulations reveal that the best fit orbits for both systems lie within narrow islands of dynamical stability, and are surrounded by much larger regions of extreme instability. As such, we show that the planets are only feasible if they are currently trapped in mutual mean-motion resonance - the 2:1 resonance in the case of 24 Sex b and c, and the 4:3 resonance in the case of HD 200964 b and c. In both cases, the region of stability is strongest and most pronounced when the planetary orbits are mutually coplanar. As the inclination of planet c with respect to planet b is increased, the stability of both systems rapidly collapses.

The observed properties of transiting exoplanets are an exceptionally rich source of information that allows us to understand and characterize their physical properties. Unfortunately, only a relatively small fraction of the known exoplanets discovered using the radial velocity technique are known to transit their host, due to the stringent orbital geometry requirements. For each target, the transit probability and predicted transit time can be calculated to great accuracy with refinement of the orbital parameters. However, the transit probability of short period and eccentric orbits can have a reasonable time dependence due to the effects of apsidal and nodal precession, thus altering their transit potential and predicted transit time. Here we investigate the magnitude of these precession effects on transit probabilities and apply this to the known radial velocity exoplanets. We assess the refinement of orbital parameters as a path to measuring these precessions and cyclic transit probabilities.

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.

We have performed a detailed dynamical study of the recently identified Neptunian Trojan 2004 KV18, only the second object to be discovered librating around Neptune's trailing Lagrange point, L5. We find that 2004 KV18 is moving on a highly unstable orbit, and was most likely captured from the Centaur population at some point in the last ~1 Myr, having originated in the Scattered Disk, beyond the orbit of Neptune. The instability of 2004 KV18 is so great that many of the test particles studied leave the Neptunian Trojan cloud within just ~0.1 - 0.3 Myr, and it takes just 37 million years for half of the 91125 test particles created to study its dynamical behaviour to be removed from the Solar system entirely. Unlike the other Neptunian Trojans previously found to display dynamical instability on hundred million year timescales (2001 QR322 and 2008 LC18), 2004 KV18 displays such extreme instability that it must be a temporarily captured Trojan, rather than a primordial member of the Neptunian Trojan population. As such, it offers a fascinating insight into the processes through which small bodies are transferred around the outer Solar system, and represents an exciting addition to the menagerie of the Solar system's small bodies.

The post-main sequence eclipsing binary NN Serpentis was recently announced as the potential host of at least two massive planetary companions. In that work, the authors put forward two potential architectures that fit the observations of the eclipsing binary with almost identical precision. In this work, we present the results of a dynamical investigation of the orbital stability of both proposed system architectures, finding that they are only stable for scenarios in which the planets are locked in mutual mean motion resonance. In the discovery work, the authors artificially fixed the orbital eccentricity of the more massive planet, NN Ser(AB) c, at 0. Here, we reanalyse the observational data on NN Serpentis without this artificial constraint, and derive a new orbital solution for the two proposed planets. We detail the results of further dynamical simulations investigating the stability of our new orbital solution, and find that allowing a small non-zero eccentricity for the outer planet renders the system unstable. We conclude that, although the original orbits proposed for the NN Serpentis planetary system prove dynamically feasible, further observations of the system are vital in order to better constrain the system's true architecture.

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.

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.

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.

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.

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.

The Neptune Trojans are the most recently discovered population of small bodies in the Solar System. To date, only eight have been discovered, though it is thought likely that the total population at least rivals that of the asteroid belt. Their origin is still the subject of some debate. Here, we detail the results of dynamical studies of two Neptune Trojans, 2001 QR322 and 2008 LC18. We find that both objects lie very close to boundaries between dynamically stable and unstable regions, with a significant probability that either or both of the objects are actually unstable on timescales of a few hundred million years. Such instability supports the idea that at least these two Neptune Trojans are dynamically captured objects, rather than objects that formed in situ. This that does not, however, rule out the possibility that these two objects were captured during Neptune's proposed post-formation migration, and have remained as Trojans ever since.

It has long been thought that the presence of a giant planet is a pre-requisite for the development of life on potentially habitable planets. Without Jupiter, it was argued, the Earth would have been subject to a punishing impact regime, which would have significantly retarded or outright prevented the development of life on our planet. Although this idea is widely embraced, little research has previously been carried out to support it. Here, we present the results of several suites of dynamical integrations used to model the influence of Jupiter's mass and orbit on the impact rate that would be experienced by the Earth. We find that, far from being a simple shield, Jupiter's role in determining the terrestrial impact flux is significantly more complicated than previously thought. Far from being a simple friend, such giant planets are perhaps more likely to imperil the development of life on otherwise habitable planets.

We present high-precision radial velocity (RV) observations of four solar-type (F7-G5) stars - HD 79498, HD 155358, HD 197037, and HD 220773 - taken as part of the McDonald Observatory Planet Search Program. For each of these stars, we see evidence of Keplerian motion caused by the presence of one or more gas giant planets in long-period orbits. We derive orbital parameters for each system, and note the properties (composition, activity, etc.) of the host stars. While we have previously announced the two-gas-giant HD 155358 system, we now report a shorter period for planet c. This new period is consistent with the planets being trapped in mutual 2:1 mean-motion resonance. We therefore perform an in-depth stability analysis, placing additional constraints on the orbital parameters of the planets. These results demonstrate the excellent long-term RV stability of the spectrometers on both the Harlan J. Smith 2.7 m telescope and the Hobby-Eberly telescope.

In early 2011, the discovery of two planets moving on surprisingly extreme orbits around the eclipsing polar cataclysmic variable system HU Aquraii was announced based on variations in the timing of mutual eclipses between the two central stars. We perform a detailed dynamical analysis of the stability of the exoplanet system as proposed in that work, revealing that it is simply dynamically unfeasible. We then apply the latest rigorous methods used by the Anglo-Australian Planet Search to analyse radial velocity data to re-examine the data used to make the initial claim. Using that data, we arrive at a significantly different orbital solution for the proposed planets, which we then show through dynamical analysis to be equally unfeasible. Finally, we discuss the need for caution in linking eclipse-timing data for cataclysmic variables to the presence of planets, and suggest a more likely explanation for the observed signal.

Recently, the first collisional family was identified in the trans-Neptunian belt. The family consists of Haumea and at least ten other ~100km-sized trans-Neptunian objects (TNOs) located in the region a = 42 - 44.5 AU. In this work, we model the long-term orbital evolution of an ensemble of fragments representing hypothetical post-collision distributions at the time of the family's birth. We consider three distinct scenarios, in which the kinetic energy of dispersed particles were varied such that their mean ejection velocities (veje) were of order 200 m/s, 300 m/s and 400 m/s, respectively. Each simulation considered resulted in collisional families that reproduced that currently observed. The results suggest that 60-75% of the fragments created in the collision will remain in the trans-Neptunian belt, even after 4 Gyr of dynamical evolution. The surviving particles were typically concentrated in wide regions of orbital element space centred on the initial impact location, with their orbits spread across a region spanning {\Delta}a ~ 6-12 AU, {\Delta}e ~ 0.1-0.15 and {\Delta}i ~ 7-10{\deg}. Most of the survivors populated the so-called Classical and Detached regions of the trans-Neptunian belt, whilst a minor fraction entered the Scattered Disk reservoir (<1%), or were captured in Neptunian mean motion resonances (<10%). In addition, except for those fragments located near strong resonances, the great majority displayed negligible long-term orbital variation. This implies that the orbital distribution of the intrinsic Haumean family can be used to constrain the orbital conditions and physical nature of the collision that created the family, billions of years ago. Indeed, our results suggest that the formation of the Haumean collisional family most likely occurred after the bulk of Neptune's migration was complete, or even some time after the migration had completely ceased.

For many years, it was assumed that Jupiter prevented the Earth from being subject to a punishing impact regime that would greatly hinder the development of life. Here, we present the 4th in a series of studies investigating this hypothesis. Previously, we examined the effect of Jupiter's mass on the impact rate experienced by Earth. Here, we extend that approach to consider the influence of Jupiter's orbital eccentricity and inclination on the impact rate. We first consider scenarios in which Jupiter's orbital eccentricity was somewhat higher and somewhat lower than that in our Solar System. We find that Jupiter's orbital eccentricity plays a moderate role in determining the impact flux at Earth, with more eccentric orbits resulting in a higher impact rate of asteroids than for more circular orbits. This is particularly pronounced at high "Jupiter" masses. For short-period comets, the same effect is clearly apparent, albeit to a lesser degree. The flux of short-period comets impacting the Earth is slightly higher for more eccentric Jovian orbits. We also consider scenarios in which Jupiter's orbital inclination was greater than that in our Solar System. Increasing Jupiter's orbital inclination greatly increased the flux of asteroidal impactors. However, at the highest tested inclination, the disruption to the Asteroid belt was so great that the belt would be entirely depleted after an astronomically short period of time. In such a system, the impact flux from asteroid bodies would therefore be very low, after an initial period of intense bombardment. By contrast, the influence of Jovian inclination on impacts from short-period comets was very small. A slight reduction in the impact flux was noted for the moderate and high inclination scenarios considered in this work - the results for inclinations of five and twenty-five degrees were essentially identical.

It has recently been proposed, on the basis of eclipse-timing data, that the eclipsing polar cataclysmic variable HU Aquarii is host to at least two giant planets. However, that result has been called into question based upon the dynamical stability of the proposed planets. In this work, we present a detailed re-analysis of all eclipse timing data available for the HU Aquarii system, making use of standard techniques used to fit orbits to radial-velocity data. We find that the eclipse timings can be used to obtain a two-planet solution that does not require the presence of additional bodies within the system. We then perform a highly detailed dynamical analysis of the proposed planetary system. We show that the improved orbital parameters we have derived correspond to planets that are dynamically unstable on unfeasibly short timescales (of order 10^4 years or less). Given these results, we discuss briefly how the observed signal might in fact be the result of the intrinsic properties of the eclipsing polar, rather than being evidence of dynamically improbable planets. Taken in concert, our results highlight the need for caution in interpreting such timing variations as being planetary in nature.

It has recently been suggested that the eclipsing polar HU Aquarii is host to at least two giant planets. We have performed highly detailed dynamical analysis of the orbits of those planets and show that the proposed system is highly unstable on timescales of < 5\times 10$^3$ years. For the coplanar orbits suggested in the discovery letter, we find stable orbital solutions for the planetary system only if the outer body moves on an orbit that brings it no closer to the host star than ~ 6 AU. The required periastron distance for the outer planet lies approximately 5 Hill radii beyond the orbit of the inner planet, and well beyond the 1-{\sigma} error bars placed on the orbit of the outer planet in the discovery letter. If the orbits of the proposed planets are significantly inclined with respect to one another, the median stability increases slightly, but such systems still become destabilised on astronomically minute timescales (typically within a few 10$^4$ years). Only in the highly improbable scenario where the outer planet follows a retrograde but coplanar orbit (i.e. inclined by 180 degrees to the orbit of the inner planet) is there any significant region of stability within the original 1-{\sigma} orbital uncertainties. Our results suggest that, if there is a second (and potentially, a third planet) in the HU Aquarii system, its orbit is dramatically different to that suggested in the discovery paper, and that more observations are critically required in order to constrain the nature of the suggested orbital bodies.

Determining the occurrence rate of terrestrial-mass planets (m_p < 10M_earth) is a critically important step on the path towards determining the frequency of Earth-like planets (eta-Earth), and hence the uniqueness of our Solar system. Current radial-velocity surveys, achieving precisions of 1 m/s, are able to detect terrestrial-mass planets and provide meaningful estimates of their occurrence rate. We present an analysis of 67 solar-type stars from the Anglo-Australian Planet Search specifically targeted for very high-precision observations. When corrected for incompleteness, we find that the planet occurrence rate increases sharply with decreasing planetary mass. Our results are consistent with those from other surveys: in periods shorter than 50 days, we find that 1.5% of stars host a giant (m_p > 100M_earth) planet, and that 18.5% of stars host a planet with m_p < 10M_earth. The preponderance of low-mass planets in short-period orbits is in conflict with formation simulations in which the majority of terrestrial-mass planets reside at larger orbital distances. This work gives a hint as to the size of eta-Earth, but to make meaningful predictions on the frequency of terrestrial planets in longer, potentially habitable orbits, low-mass terrestrial planet searches at periods of 100-200 days must be made an urgent priority for ground-based Doppler planet searches in the years ahead.

We present results examining the fate of the Trojan clouds produced in our previous work. We find that the stability of Neptunian Trojans seems to be strongly correlated to their initial post-migration orbital elements, with those objects that survive as Trojans for billions of years displaying negligible orbital evolution. The great majority of these survivors began the integrations with small eccentricities (e < 0.2) and small libration amplitudes (A < 30 - 40{\deg}). The survival rate of "pre-formed" Neptunian Trojans (which in general survived on dynamically cold orbits (e < 0.1, i < 5 - 10{\deg})) varied between ~5 and 70%. By contrast, the survival rate of "captured" Trojans (on final orbits spread across a larger region of e-i element space) were markedly lower, ranging between 1 and 10% after 4 Gyr. Taken in concert with our earlier work, we note that planetary formation scenarios which involve the slow migration (a few tens of millions of years) of Neptune from an initial planetary architecture that is both resonant and compact (aN < 18 AU) provide the most promising fit of those we considered to the observed Trojan population. In such scenarios, we find that the current day Trojan population would number ~1% of that which was present at the end of the planet's migration, with the bulk being sourced from captured, rather than pre-formed objects. We note, however, that even those scenarios still fail to reproduce the currently observed portion of the Neptune Trojan population moving on orbits with e < 0.1 but i > 20{\deg}. Dynamical integrations of the currently observed Trojans show that five out of the seven are dynamically stable on 4 Gyr timescales, while 2001 QR322, exhibits significant dynamical instability. The seventh Trojan object, 2008 LC18, has such large orbital uncertainties that only future studies will be able to determine its stability.

We present a short review of the impact regime experienced by the terrestrial planets within our own Solar system, describing the three populations of potentially hazardous objects which move on orbits that take them through the inner Solar system. Of these populations, the origins of two (the Near-Earth Asteroids and the Long-Period Comets) are well understood, with members originating in the Asteroid belt and Oort cloud, respectively. By contrast, the source of the third population, the Short-Period Comets, is still under debate. The proximate source of these objects is the Centaurs, a population of dynamically unstable objects that pass perihelion between the orbits of Jupiter and Neptune. However, a variety of different origins have been suggested for the Centaur population. Here, we present evidence that at least a significant fraction of the Centaur population can be sourced from the planetary Trojan clouds, stable reservoirs of objects moving in 1:1 mean-motion resonance with the giant planets (primarily Jupiter and Neptune). Focusing on simulations of the Neptunian Trojan population, we show that an ongoing flux of objects should be leaving that region to move on orbits within the Centaur population. With conservative estimates of the flux from the Neptunian Trojan clouds, we show that their contribution to that population could be of order ~3%, while more realistic estimates suggest that the Neptune Trojans could even be the main source of fresh Centaurs. We suggest that further observational work is needed to constrain the contribution made by the Neptune Trojans to the ongoing flux of material to the inner Solar system, and believe that future studies of the habitability of exoplanetary systems should take care not to neglect the contribution of resonant objects (such as planetary Trojans) to the impact flux that could be experienced by potentially habitable worlds.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

Within the next few years, the first Earth-mass planets will be discovered around other stars. Some of those worlds will certainly lie within the classical "habitable zone" of their parent stars, and we will quickly move from knowing of no exoEarths to knowing many. For the first time, we will be in a position to carry out a detailed search for the first evidence of life beyond our Solar System. However, such observations will be hugely taxing and time consuming to perform, and it is almost certain that far more potentially habitable worlds will be known than it is possible to study. It is therefore important to catalogue and consider the various effects which make a promising planet more or less suitable for the development of life. In this work, we review the various planetary, dynamical and stellar influences that could influence the habitability of exoEarths. The various influences must be taken in concert when we attempt to decide where to focus our first detailed search for life. While there is no guarantee that any given planet will be inhabited, it is vitally important to ensure that we focus our time and effort on those planets most likely to yield a positive result.

Thermal emission from Kuiper Belt object (136108) Haumea was measured with Herschel-PACS at 100 and 160 micrometers for almost a full rotation period. Observations clearly indicate a 100-micrometer thermal lightcurve with an amplitude of a factor of ~ 2, which is positively correlated with the optical lightcurve. This confirms that both are primarily due to shape effects. A 160-micrometer lightcurve is marginally detected. Radiometric fits of the mean Herschel- and Spitzer- fluxes indicate an equivalent diameter D ~ 1300 km and a geometric albedo p_v ~ 0.70-0.75. These values agree with inferences from the optical lightcurve, supporting the hydrostatic equilibrium hypothesis. The large amplitude of the 100-micrometer lightcurve suggests that the object has a high projected a/b axis ratio (~ 1.3) and a low thermal inertia as well as possible variable infrared beaming. This may point to fine regolith on the surface, with a lunar-type photometric behavior. The quality of the thermal data is not sufficient to clearly detect the effects of a surface dark spot.

We present results from the earliest observations of DEBRIS, a Herschel Key Programme to conduct a volume- and flux-limited survey for debris discs in A-type through M-type stars. PACS images (from chop/nod or scan-mode observations) at 100 and 160 micron are presented toward two A-type stars and one F-type star: beta Leo, beta UMa and eta Corvi. All three stars are known disc hosts. Herschel spatially resolves the dust emission around all three stars (marginally, in the case of beta UMa), providing new information about discs as close as 11 pc with sizes comparable to that of the Solar System. We have combined these data with existing flux density measurements of the discs to refine the SEDs and derive estimates of the fractional luminosities, temperatures and radii of the discs.

The goal of the Herschel Open Time Key programme "TNOs are Cool!" is to derive the physical and thermal properties for a large sample of Centaurs and trans-Neptunian objects (TNOs), including resonant, classical, detached and scattered disk objects. We present results for seven targets either observed in PACS point-source, or in mini scan-map mode. Spitzer-MIPS observations were included for three objects. The sizes of these targets range from 100 km to almost 1000 km, five have low geometric albedos below 10%, (145480) 2005 TB190 has a higher albedo above 15%. Classical thermal models driven by an intermediate beaming factor of $\eta$=1.2 or $\eta$-values adjusted to the observed colour temperature fit the multi-band observations well in most cases. More sophisticated thermophysical models give very similar diameter and albedo values for thermal inertias in the range 0-25 Jm-2s-0.5K-1, consistent with very low heat conductivities at temperatures far away from the Sun. The early experience with observing and model strategies will allow us to derive physical and thermal properties for our complete Herschel TNO sample of 140 targets as a benchmark for understanding the solar system debris disk, and extra-solar ones as well.

Of the four giant planets in the Solar system, only Jupiter and Neptune are currently known to possess swarms of Trojan asteroids - small objects that experience a 1:1 mean motion resonance with their host planet. In Lykawka et al. (2009), we performed extensive dynamical simulations, including planetary migration, to investigate the origin of the Neptunian Trojan population. Utilising the vast amount of simulation data obtained for that work, together with fresh results from new simulations, we here investigate the dynamical capture of Trojans by all four giant planets from a primordial trans-Neptunian disk. We find the likelihood of a given planetesimal from this region being captured onto an orbit within Jupiter's Trojan cloud lies between several times 10^-6 and 10^-5. For Saturn, the probability is found to be in the range <10^-6 to 10^-5, whilst for Uranus the probabilities range between 10^-5 and 10^-4. Finally, Neptune displays the greatest probability of Trojan capture, with values ranging between 10^-4 and 10^-3. Our results suggest that all four giant planets are able to capture and retain a significant population of Trojan objects from the disk by the end of planetary migration. As a result of encounters with the giant planets prior to Trojan capture, these objects tend to be captured on orbits that are spread over a wide range of orbital eccentricities and inclinations. The bulk of captured objects are to some extent dynamically unstable, and therefore the populations of these objects tend to decay over the age of the Solar System, providing an important ongoing source of new objects moving on dynamically unstable orbits among the giant planets. Given that a huge population of objects would be displaced by Neptune's outward migration (with a potential cumulative mass a number of times that of the Earth), we conclude that the surviving remnant of the Trojans captured during the migration of the outer planets might be sufficient to explain the currently known Trojan populations in the outer Solar system.

The fact that the Centaurs are the primary source of the Short Period Comets is well established. However, the origin of the Centaurs themselves is still under some debate, with a variety of different source reservoirs being proposed in the last decade. In this work, we suggest that the Neptune Trojans (together with the Jovian Trojans) could represent an additional significant source of Centaurs. Using dynamical simulations of the first Neptune Trojan discovered (2001 QR322), together with integrations following the evolution of clouds of theoretical Neptune Trojans obtained during simulations of planetary migration, we show that the Neptune Trojan population contains a great number of objects which are unstable on both Myr and Gyr timescales. Using individual examples, we show how objects that leave the Neptunian Trojan cloud evolve onto orbits indistinguishable from those of the known Centaurs, before providing a range of estimates of the flux from this region to the Centaur population. With only moderate assumptions, it is shown that the Trojans can contribute a significant proportion of the Centaur population, and may even be the dominant source reservoir. This result is supported by past work on the colours of the Trojans and the Centaurs, but it will take future observations to determine the full scale of the contribution of the escaped Trojans to the Centaur population.

We present the results of detailed dynamical simulations of the effect of the migration of the four giant planets on both the transport of pre-formed Neptune Trojans, and the capture of new Trojans from a trans-Neptunian disk. We find that scenarios involving the slow migration of Neptune over a large distance (50Myr to migrate from 18.1AU to its current location) provide the best match to the properties of the known Trojans. Scenarios with faster migration (5Myr), and those in which Neptune migrates from 23.1AU to its current location, fail to adequately reproduce the current day Trojan population. Scenarios which avoid disruptive perturbation events between Uranus and Neptune fail to yield any significant excitation of pre-formed Trojans (transported with efficiencies between 30 and 98% whilst maintaining the dynamically cold nature of these objects). Conversely, scenarios with periods of strong Uranus-Neptune perturbation lead to the almost complete loss of such pre-formed objects. In these cases, a small fraction (~0.15%) of these escaped objects are later recaptured as Trojans prior to the end of migration, with a wide range of eccentricities (<0.35) and inclinations (<40 deg). In all scenarios (including those with such disruptive interaction between Uranus and Neptune) the capture of objects from the trans-Neptunian disk (through which Neptune migrates) is achieved with efficiencies between ~0.1 and ~1%. The captured Trojans display a wide range of inclinations (<40 deg for slow migration, and <20 deg for rapid migration) and eccentricities (<0.35), and we conclude that, given the vast amount of material which undoubtedly formed beyond the orbit of Neptune, such captured objects may be sufficient to explain the entire Neptune Trojan population. (Shortened version)

In this work, we investigate the dynamical stability of pre-formed Neptune Trojans under the gravitational influence of the four giant planets in compact planetary architectures, over 10 Myr. In our modelling, the initial orbital locations of Uranus and Neptune (aN) were varied to produce systems in which those planets moved on non-resonant orbits, or in which they lay in their mutual 1:2, 2:3 and 3:4 mean-motion resonances (MMRs). In total, 420 simulations were carried out, examining 42 different architectures, with a total of 840000 particles across all runs. In the non-resonant cases, the Trojans suffered only moderate levels of dynamical erosion, with the most compact systems (those with aN less than or equal 18 AU) losing around 50% of their Trojans by the end of the integrations. In the 2:3 and 3:4 MMR scenarios, however, dynamical erosion was much higher with depletion rates typically greater than 66% and total depletion in the most compact systems. The 1:2 resonant scenarios featured disruption on levels intermediate between the non-resonant cases and other resonant scenarios, with depletion rates of the order of tens of percent. Overall, the great majority of plausible pre-migration planetary architectures resulted in severe levels of depletion of the Neptunian Trojan clouds. In particular, if Uranus and Neptune formed near their mutual 2:3 or 3:4 MMR and at heliocentric distances within 18 AU (as favoured by recent studies), we found that the great majority of pre-formed Trojans would have been lost prior to Neptune's migration. This strengthens the case for the great bulk of the current Neptunian Trojan population having been captured during that migration.

Since early work on the stability of the first Neptunian Trojan, 2001 QR322, suggested that it was a dynamically stable, primordial body, it has been assumed this applies to both that object, and its more recently discovered brethren. However, it seems that things are no longer so clear cut. In this work, we present the results of detailed dynamical simulations of the orbital behaviour of 2001 QR322. Using an ephemeris for the object that has significantly improved since earlier works, we follow the evolution of 19683 test particles, placed on orbits within the observational error ellipse of 2001 QR322's orbit, for a period of 1 Gyr. We find that majority of these "clones" of 2001 QR322 are dynamically unstable, exhibiting a near-exponential decay from both the Neptunian Trojan cloud (decay halflife ~550 Myr) and the Solar system (decay halflife ~590 Myr). The stability of the object within Neptune's Trojan cloud is found to be strongly dependent on the initial semi-major axis used, with those objects located at semimajor axis equal or greater than 30.30 AU being significantly less stable than those interior to this value, as a result of their having initial libration amplitudes very close to a critical threshold dividing regular and irregular motion, located at ~70-75 deg (full extent of angular motion). This result suggests that, if 2001 QR322 is a primordial Neptunian Trojan, it must be a representative of a population that was once significantly larger than that we see today, and adds weight to the idea that the Neptune Trojans may represent a significant source of objects moving on unstable orbits between the giant planets (the Centaurs).

The difference between the measured atmospheric abundances of neon, argon, krypton and xenon for Venus, the Earth and Mars is striking. Because these abundances drop by at least two orders of magnitude as one moves outward from Venus to Mars, the study of the origin of this discrepancy is a key issue that must be explained if we are to fully understand the different delivery mechanisms of the volatiles accreted by the terrestrial planets. In this work, we aim to investigate whether it is possible to quantitatively explain the variation of the heavy noble gas abundances measured on Venus, the Earth and Mars, assuming that cometary bombardment was the main delivery mechanism of these noble gases to the terrestrial planets. To do so, we use recent dynamical simulations that allow the study of the impact fluxes of comets upon the terrestrial planets during the course of their formation and evolution. Assuming that the mass of noble gases delivered by comets is proportional to rate at which they collide with the terrestrial planets, we show that the krypton and xenon abundances in Venus and the Earth can be explained in a manner consistent with the hypothesis of cometary bombardment. In order to explain the krypton and xenon abundance differences between the Earth and Mars, we need to invoke the presence of large amounts of CO2-dominated clathrates in the Martian soil that would have efficiently sequestered these noble gases.

It has long been assumed that the planet Jupiter acts as a giant shield, significantly lowering the impact rate of small bodies on the Earth. However, until recently, very little work had been carried out examining the role played by Jupiter in determining the frequency of such collisions. In this work, the third of a series of papers, we examine the degree to which the impact rate on Earth resulting from the Oort cloud comets is enhanced or lessened by the presence of a giant planet in a Jupiter-like orbit, in an attempt to more fully understand the impact regime under which life on Earth has developed. Our results show that the presence of a giant planet in a Jupiter-like orbit significantly alters the impact rate of Oort cloud comets on the Earth, decreasing the rate as the mass of the giant increases. The greatest bombardment flus is observed when no giant planet is present.

Over 300 extrasolar planets (exoplanets) have been detected orbiting nearby stars. We now hope to conduct a census of all planets around nearby stars and to characterize their atmospheres and surfaces with spectroscopy. Rocky planets within their star's habitable zones have the highest priority, as these have the potential to harbor life. Our science goal is to find and characterize all nearby exoplanets; this requires that we measure the mass, orbit, and spectroscopic signature of each one at visible and infrared wavelengths. The techniques for doing this are at hand today. Within the decade we could answer long-standing questions about the evolution and nature of other planetary systems, and we could search for clues as to whether life exists elsewhere in our galactic neighborhood.

It has long been assumed that the planet Jupiter acts as a giant shield, significantly lowering the impact rate of minor bodies upon the Earth, and thus enabling the development and evolution of life in a collisional environment which is not overly hostile. However, in the past, little work has been carried out to examine the validity of this idea. In the second of a series of papers, we examine the degree to which the impact risk resulting from objects on Centaur-like orbits is affected by the presence of a giant planet, in a continuing attempt to fully understand the impact regime under which life on Earth has developed. The Centaurs, which occupy orbits beyond Jupiter, have their origins in the Edgeworth-Kuiper belt that extends beyond Neptune. The giant planets peturb the Centaurs, sending a significanr fraction into the inner Solar System where they become visible as short-period comets. In this work we present results which show that the presence of a giant planet can act to significantly change the impact rate of short-period comets on the Earth, and that a giant planet often actually increases the impact flux greatly over that which would be expected were it not present. (Shortened version of abstract.)

The asteroids are the major source of potential impactors on the Earth today. It has long been assumed that the giant planet Jupiter acts as a shield, significantly lowering the impact rate on the Earth from both cometary and asteroidal bodies. Such shielding, it is claimed, enabled the development and evolution of life in a collisional environment which is not overly hostile. The reduced frequency of impacts, and of related mass extinctions, would have allowed life the time to thrive, where it would otherwise have been suppressed. However, in the past, little work has been carried out to examine the validity of this idea. In the first of several papers, we examine the degree to which the impact risk resulting from a population representative of the asteroids is enhanced or lessened by the presence of a giant planet, in an attempt to fully understand the impact regime under which life on Earth has developed. Our results show that the situation is far less clear cut that has previously been assumed - for example, the presence of a giant planet can act to enhance the impact rate of asteroids on the Earth significantly.

We consider the role of deuterium as a potential marker of location and ambient conditions during the formation of small bodies in our Solar system. We concentrate in particular on the formation of the regular icy satellites of Jupiter and the other giant planets, but include a discussion of the implications for the Trojan asteroids and the irregular satellites. We examine in detail the formation of regular planetary satellites within the paradigm of a circum-Jovian subnebula. Particular attention is paid to the two extreme potential subnebulae - "hot" and "cold". In particular, we show that, for the case of the "hot" subnebula model, the D:H ratio in water ice measured from the regular satellites would be expected to be near-Solar. In contrast, satellites which formed in a "cold" subnebula would be expected to display a D:H ratio that is distinctly over-Solar. We then compare the results obtained with the enrichment regimes which could be expected for other families of icy small bodies in the outer Solar system - the Trojan asteroids and the irregular satellites. In doing so, we demonstrate how measurements by Laplace, the James Webb Space Telescope, HERSCHEL and ALMA will play an important role in determining the true formation locations and mechanisms of these objects.

We propose a scenario for the formation of the Main Belt in which asteroids incorporated icy particles formed in the outer Solar Nebula. We calculate the composition of icy planetesimals formed beyond a heliocentric distance of 5 AU in the nebula by assuming that the abundances of all elements, in particular that of oxygen, are solar. As a result, we show that ices formed in the outer Solar Nebula are composed of a mix of clathrate hydrates, hydrates formed above 50 K and pure condensates produced at lower temperatures. We then consider the inward migration of solids initially produced in the outer Solar Nebula and show that a significant fraction may have drifted to the current position of the Main Belt without encountering temperature and pressure conditions high enough to vaporize the ices they contain. We propose that, through the detection and identification of initially buried ices revealed by recent impacts on the surfaces of asteroids, it could be possible to infer the thermodynamic conditions that were present within the Solar Nebula during the accretion of these bodies, and during the inward migration of icy planetesimals. We also investigate the potential influence that the incorporation of ices in asteroids may have on their porosities and densities. In particular, we show how the presence of ices reduces the value of the bulk density of a given body, and consequently modifies its macro-porosity from that which would be expected from a given taxonomic type.

Protoplanetary disks start their lives with a dust free inner region where the temperatures are higher than the sublimation temperature of solids. As the star illuminates the innermost particles, which are immersed in gas at the sublimation edge, these particles are subject to a photophoretic force. We examine the motion of dust particles at the inner edge of protoplanetary disks due to photophoretic drag. We find that, in a laminar disk, photophoretic drag increases the size of the inner hole after accretion onto the central body has become subdued. This region within the hole becomes an optically transparent zone containing gas and large dusty particles (>>10 cm), but devoid of, or strongly depleted in, smaller dust aggregates. Photophoresis can clear the inner disk of dust out to 10 AU in less than 1 Myr. At late times, the edge reaches a stable equilibrium between inward drift and photophoretic outward drift, at a distance of some tens of AU. Eventually, the edge will move inwards again as the disk disperses, shifting the equilibrium position back from about 40 AU to below 30 AU in 1-2 Myr in the disk model. In a turbulent disk, diffusion can delay the clearing of a disk by photophoresis. Smaller and/or age-independent holes of radii of a few AU are also possible outcomes of turbulent diffusion counteracting photophoresis. This outward and then inward moving edge marks a region of high dust concentration. This density enhancement, and the efficient transport of particles from close to the star to large distances away, can explain features of comets such as high measured ratios of crystalline to amorphous silicates, and has a large number of other applications.

Large scale simulations of Centaurs have yielded vast amounts of data, the analysis of which allows interesting but uncommon scenarios to be studied. One such rare phenomenon is the temporary capture of Centaurs as Trojans of the giant planets. Such captures are generally short (10 kyr to 100 kyr), but occur with sufficient frequency (about 40 objects larger than 1 km in diameter every Myr) that they may well contribute to the present-day populations. Uranus and Neptune seem to have great difficulty capturing Centaurs into the 1:1 resonance, while Jupiter captures some, and Saturn the most (80 %). We conjecture that such temporary capture from the Centaur population may be the dominant delivery route into the Saturnian Trojans. Photometric studies of the Jovian Trojans may reveal outliers with Centaur-like as opposed to asteroidal characteristics, and these would be prime candidates for captured Centaurs.

Detailed orbit integrations of clones of five Centaurs -- namely, 1996 AR20, 2060 Chiron, 1995 SN55, 2000 FZ53 and 2002 FY36 -- for durations of 3 Myr are presented. One of our Centaur sample starts with perihelion initially under the control of Jupiter (1996 AR20), two start under the control of Saturn (Chiron and 1995 SN55) and one each starts under the control of Uranus (2000 FZ53) and Neptune (2002 FY36) respectively. A variety of interesting pathways are illustrated with detailed examples including: capture into the Jovian Trojans, repeated bursts of short-period comet behaviour, capture into mean-motion resonances with the giant planets and into Kozai resonances, as well as traversals of the entire Solar system. For each of the Centaurs, we provide statistics on the numbers (i) ejected, (ii) showing short-period comet behaviour and (iii) becoming Earth and Mars crossing. For example, Chiron has over 60 % of its clones becoming short-period objects, whilst 1995 SN55 has over 35 %. Clones of these two Centaurs typically make numerous close approaches to Jupiter. At the other extreme, 2000 FZ53 has roughly 2 % of its clones becoming short-period objects. In our simulations, typically 20 % of the clones which become short-period comets subsequently evolve into Earth-crossers.

Large-scale simulations of the Centaur population are carried out. The evolution of 23328 particles based on the orbits of 32 well-known Centaurs is followed for up to 3 Myr in the forward and backward direction under the influence of the 4 massive planets. The objects exhibit a rich variety of dynamical behaviour with half-lives ranging from 540 kyr (1996 AR20) to 32 Myr (2000 FZ53). The mean half-life of the entire sample of Centaurs is 2.7 Myr. The data are analyzed using a classification scheme based on the controlling planets at perihelion and aphelion, previously given in Horner et al (2003). Transfer probabilities are computed and show the main dynamical pathways of the Centaur population. The total number of Centaurs with diameters larger than 1 km is estimated as roughly 44300, assuming an inward flux of one new short-period comet every 200 yrs. The flux into the Centaur region from the Edgeworth-Kuiper belt is estimated to be 1 new object every 125 yrs. Finally, the flux from the Centaur region to Earth-crossing orbits is 1 new Earth-crosser every 880 yrs

A new classification scheme is introduced for comet-like bodies in the Solar system. It covers the traditional comets as well as the Centaurs and Edgeworth-Kuiper belt objects. At low inclinations, close encounters with planets often result in near-constant perihelion or aphelion distances, or in perihelion-aphelion interchanges, so the minor bodies can be labelled according to the planets predominantly controlling them at perihelion and aphelion. For example, a JN object has a perihelion under the control of Jupiter and aphelion under the control of Neptune, and so on. This provides 20 dynamically distinct categories of outer Solar system objects in the Jovian and trans-Jovian regions. The Tisserand parameter with respect to the planet controlling perihelion is also often roughly constant under orbital evolution. So, each category can be further sub-divided according to the Tisserand parameter. The dynamical evolution of comets, however, is dominated not by the planets nearest at perihelion or aphelion, but by the more massive Jupiter. The comets are separated into four categories -- Encke-type, short-period, intermediate and long-period -- according to aphelion distance. The Tisserand parameter categories now roughly correspond to the well-known Jupiter-family comets, transition-types and Halley-types. In this way, the nomenclature for the Centaurs and Edgeworth-Kuiper belt objects is based on, and consistent with, that for comets.

Two sets of investigators -- Murray (1999) and Matese, Whitman & Whitmire (1999) -- have recently claimed evidence for an undiscovered Solar System planet from possible great circle alignments in the aphelia directions of the long period comets. However, comet discoveries are bedevilled by selection effects. These include anomalies caused by the excess of observers in the northern as against the southern hemisphere, seasonal and diurnal biases, directional effects which make it harder to discover comets in certain regions of the sky, as well as sociological biases. The stream proposed by Murray is shown on an equal area Hammer-Aitoff projection. The addition of newer data weakens the case for the alignment. There is also evidence that the subsample in the stream is affected by seasonal and north-south biases. The stream proposed by Matese et al. is most obvious in the sample of dynamically new comets, and especially in those whose orbits are best known. The most recent data continues to maintain the overpopulation in this great circle. This pattern in the data occurs about once in a thousand times by chance. Numerical integrations are used to demonstrate that a planet by itself can reduce the perihelia of comets in its orbital plane to sufficiently small values so that they could be discovered from the Earth. There is a need for a sample of long period comets that is free from unknown or hard-to-model selection effects. Such will be provided by the European Space Agency satellite GAIA.