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.

The presence of dusty debris around main sequence stars denotes the existence of planetary systems. Such debris disks are often identified by the presence of excess continuum emission at infrared and (sub-)millimetre wavelengths, with measurements at longer wavelengths tracing larger and cooler dust grains. The exponent of the slope of the disk emission at sub-millimetre wavelengths, `q', defines the size distribution of dust grains in the disk. This size distribution is a function of the rigid strength of the dust producing parent planetesimals. As part of the survey `PLAnetesimals around TYpical Pre-main seqUence Stars' (PLATYPUS) we observed six debris disks at 9-mm using the Australian Telescope Compact Array. We obtain marginal (~3-\sigma) detections of three targets: HD 105, HD 61005, and HD 131835. Upper limits for the three remaining disks, HD20807, HD109573, and HD109085, provide further constraint of the (sub-)millimetre slope of their spectral energy distributions. The values of q (or their limits) derived from our observations are all smaller than the oft-assumed steady state collisional cascade model (q = 3.5), but lie well within the theoretically expected range for debris disks q ~ 3 to 4. The measured q values for our targets are all < 3.3, consistent with both collisional modelling results and theoretical predictions for parent planetesimal bodies being `rubble piles' held together loosely by their self-gravity.

Radial and azimuthal features (such as disc offsets and eccentric rings) seen in high resolution images of debris discs, provide us with the unique opportunity of finding potential planetary companions which betray their presence by gravitationally sculpting such asymmetric features. The young debris disc around HD 115600, imaged recently by the Gemini Planet Imager, is such a disc with an eccentricity 0.1<e<0.2 and a projected offset from the star of 4 AU. Using our modified N-body code which incorporates radiation forces, we firstly aim to determine the orbit of a hidden planetary companion potentially responsible for shaping the disc. We run a suite of simulations covering a broad range of planetary parameters using a Monte Carlo Markov Chain sampling method and create synthetic images from which we extract the geometric disc parameters to be compared with the observed and model-derived quantities. We then repeat the study using a traditional grid to explore the planetary parameter space and aim secondly to compare the efficiency of both sampling methods. We find a planet of 7.8 Jupiter mass orbiting at 30 AU with an eccentricity of e = 0.2 to be the best fit to the observations of HD 115600. Technically, such planet has a contrast detectable by direct imaging, however the system's orientation does not favour such detection. In this study, at equal number of explored planetary configurations, the Monte Carlo Markov Chain not only converges faster but provides a better fit than a traditional grid.

Resolved debris disc images can exhibit a range of radial and azimuthal structures, including gaps and rings, which can result from planetary companions shaping the disc by their gravitational influence. Currently there are no tools available to determine the architecture of potential companions from disc observations. Recent work by Rodigas et al. (2014) presents how one can estimate the maximum mass and minimum semi major axis of a hidden planet empirically from the width of the disc in scattered light. In this work, we use the predictions of Rodigas et al. applied to two debris discs HD 202628 and HD 207129. We aim to test if the predicted orbits of the planets can explain the features of their debris disc, such as eccentricity and sharp inner edge. We first run dynamical simulations using the predicted planetary parameters of Rodigas et al., and then numerically search for better parameters. Using a modified N-body code including radiation forces, we perform simulations over a broad range of planet parameters and compare synthetics images from our simulations to the observations. We find that the observational features of HD 202628 can be reproduced with a planet five times smaller than expected, located 30 AU beyond the predicted value, while the best match for HD 207129 is for a planet located 5-10 AU beyond the predicted location with a smaller eccentricity. We conclude that the predictions of Rodigas et al. provide a good starting point but should be complemented by numerical simulations.

Numerical simulations are a crucial tool to understand the relationship between debris discs and planetary companions. However, simulations throughout the literature have been conducted with various initial conditions often with little or no justification. In this paper, we aim to study the dependence on the initial conditions of N-body simulations modelling the interaction between a massive and eccentric planet on an exterior debris disc. To achieve this, we first classify three broad approaches used in the literature and provide some physical context for when each category should be used. We then run a series of N-body simulations, that include radiation forces acting on small grains, with varying initial conditions across the three categories. We test the influence of the initial parent body belt width, eccentricity, and alignment with the planet on the resulting debris disc structure and compare the final peak emission location, disc width and offset of synthetic disc images produced with a radiative transfer code. We also track the evolution of the forced eccentricity of the dust grains induced by the planet, as well as resonance dust trapping. We find that an initially broad parent body belt always results in a broader debris disc than an initially narrow parent body belt. While simulations with a parent body belt with low initial eccentricity (e ~ 0) and high initial eccentricity (0 < e < 0.3) resulted in similar broad discs, we find that purely secular forced initial conditions, where the initial disc eccentricity is set to the forced value and the disc is aligned with the planet, always result in a narrower disc. We conclude that broad debris discs can be modelled by using either a dynamically cold or dynamically warm parent belt, while in contrast eccentric narrow debris rings are reproduced using a secularly forced parent body belt.

The debris disk closest to Earth is the one around the star epsilon Eridani at a distance of 3.2 pc. It is the prime target for detailed studies of a belt of planetesimals left from the early phase of planet formation other than the Kuiper Belt. The non-uniform ring-like structure around epsilon Eridani, originally discovered at lambda=850 microns with the bolometer camera SCUBA, could be the signpost of unseen long-period planets interior to the disk that gravitationally interact with it through mean-motion resonances. However, the reliability of the structure at 850 microns, which has been debated, has not been verified with independent observations until now. We present a high signal-to-noise ratio image of this structure at lambda=1.2 mm made with the bolometer camera MAMBO and compare this with the SCUBA image. We have found that three of the four emission clumps (NE, NW, SW) and the two deep hollows to the east and west are at the same positions in the MAMBO and SCUBA images within astrometric uncertainty. The SE clump is at odds, significantly brighter and more extended in the SCUBA than in the MAMBO images, but it is possible that this mismatch is an artifact. We conclude that this degree of positional coincidence provides tentative evidence that the observed structure is robust. In addition, we present the radial brightness profile of our MAMBO image and show that the width of the planetesimal belt around epsilon Eridani is narrower than 22 AU, a more stringent upper limit than determined from previous observations. The corresponding relative width is $0.1 < \Delta R / R < 0.4$, which is lower than for the Kuiper Belt.

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.

Star formation does not occur until the onset of gravitational collapse inside giant molecular clouds. However, the conditions that initiate cloud collapse and regulate the star formation process remain poorly understood. Local processes such as turbulence and magnetic fields can act to promote or prevent collapse. On larger scales, the galactic potential can also influence cloud stability and is traditionally assessed by the tidal and shear effects. In this paper, we examine the stability of giant molecular clouds (GMCs) in the Large Magellanic Cloud (LMC) against shear and the galactic tide using CO data from the Magellanic Mopra Assessment (MAGMA) and rotation curve data from the literature. We calculate the tidal acceleration experienced by individual GMCs and determine the minimum cloud mass required for tidal stability. We also calculate the shear parameter, which is a measure of a clouds susceptibility to disruption via shearing forces in the galactic disk. We examine whether there are correlations between the properties and star forming activity of GMCs and their stability against shear and tidal disruption. We find that the GMCs are in approximate tidal balance in the LMC, and that shear is unlikely to affect their further evolution. GMCs with masses close to the minimal stable mass against tidal disruption are not unusual in terms of their mass, location, or CO brightness, but we note that GMCs with large velocity dispersion tend to be more sensitive to tidal instability. We also note that GMCs with smaller radii, which represent the majority of our sample, tend to more strongly resist tidal and shear disruption. Our results demonstrate that star formation in the LMC is not inhibited by to tidal or shear instability.

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.