A dozen short-period detached binaries are known to host transiting circumbinary planets. In all circumbinary systems so far, the planetary and binary orbits are aligned within a couple of degrees. However, the obliquity of the primary star, which is an important tracer of their formation, evolution, and tidal history, has only been measured in one circumbinary system until now. EBLM J0608-59/TOI-1338 is a low-mass eclipsing binary system with a recently discovered circumbinary planet identified by TESS. Here, we perform high-resolution spectroscopy during primary eclipse to measure the projected stellar obliquity of the primary component. The obliquity is low, and thus the primary star is aligned with the binary and planetary orbits with a projected spin-orbit angle $\beta = 2.8 \pm 17.1$ deg. The rotation period of $18.1 \pm 1.6$ days implied by our measurement of $v\sin{i_\star}$ suggests that the primary has not yet pseudo-synchronized with the binary orbit, but is consistent with gyrochronology and weak tidal interaction with the binary companion. Our result, combined with the known coplanarity of the binary and planet orbits, is suggestive of formation from a single disc. Finally, we considered whether the spectrum of the faint secondary star could affect our measurements. We show through simulations that the effect is negligible for our system, but can lead to strong biases in $v\sin{i_\star}$ and $\beta$ for higher flux ratios. We encourage future studies in eclipse spectroscopy test the assumption of a dark secondary for flux ratios $\gtrsim 1$ ppt.

The interiors of giant planets remain poorly understood. Even for the planets in the Solar System, difficulties in observation lead to major uncertainties in the properties of planetary cores. Exoplanets that have undergone rare evolutionary pathways provide a new route to understanding planetary interiors. We present the discovery of TOI-849b, the remnant core of a giant planet, with a radius smaller than Neptune but an anomalously high mass $M_p=40.8^{+2.4}_{-2.5}M_{\oplus}$ and density of $5.5\pm0.8$gcm$^{-3}$, similar to the Earth. Interior structure models suggest that any gaseous envelope of pure hydrogen and helium consists of no more than $3.9^{+0.8}_{-0.9}$% of the total mass of the planet. TOI-849b transits a late G type star (T$_{\rm mag}=11.5$) with an orbital period of 18.4 hours, leading to an equilibrium temperature of 1800K. The planet's mass is larger than the theoretical threshold mass for runaway gas accretion. As such, the planet could have been a gas giant before undergoing extreme mass loss via thermal self-disruption or giant planet collisions, or it avoided substantial gas accretion, perhaps through gap opening or late formation. Photoevaporation rates cannot provide the mass loss required to reduce a Jupiter-like gas giant, but can remove a few $M_\oplus$ hydrogen and helium envelope on timescales of several Gyr, implying that any remaining atmosphere is likely to be enriched by water or other volatiles from the planetary interior. TOI-849b represents a unique case where material from the primordial core is left over from formation and available to study.

We present the detection of high energy white-light flares from pre-main sequence stars associated with the Orion complex, observed as part of the Next Generation Transit Survey (NGTS). With energies up to $5.2\times10^{35}$ erg these flares are some of the most energetic white-light flare events seen to date. We have used the NGTS observations of flaring and non-flaring stars to measure the average flare occurrence rate for 4 Myr M0-M3 stars. We have also combined our results with those from previous studies to predict average rates for flares above $1\times10^{35}$ ergs for early M stars in nearby young associations.

We report the discovery of NGTS-11 b (=TOI-1847 b), a transiting Saturn in a 35.46-day orbit around a mid K-type star (Teff=5050 K). We initially identified the system from a single-transit event in a TESS full-frame image light-curve. Following seventy-nine nights of photometric monitoring with an NGTS telescope, we observed a second full transit of NGTS-11 b approximately one year after the TESS single-transit event. The NGTS transit confirmed the parameters of the transit signal and restricted the orbital period to a set of 13 discrete periods. We combined our transit detections with precise radial velocity measurements to determine the true orbital period and measure the mass of the planet. We find NGTS-11 b has a radius of 0.817+0.028-0.032 $R_J$, a mass of 0.344+0.092-0.073 $M_J$, and an equilibrium temperature of just 435+34-32 K, making it one of the coolest known transiting gas giants. NGTS-11 b is the first exoplanet to be discovered after being initially identified as a TESS single-transit event, and its discovery highlights the power of intense photometric monitoring in recovering longer-period transiting exoplanets from single-transit events.

We present the discovery of NGTS J214358.5-380102, an eccentric M-dwarf binary discovered by the Next Generation Transit Survey. The system period of 7.618 days is greater than many known eclipsing M-dwarf binary systems. Its orbital eccentricity of $0.323^{+0.0014}_{-0.0037}$, is large relative to the period and semi-major axis of the binary. Global modelling of photometry and radial velocities indicate stellar masses of $M_A$=$0.426 ^{+0.0056}_{-0.0049}$, $M_B$=$0.455 ^{+0.0058}_{-0.0052}$ and stellar radii $R_A$=$0.461 ^{+0.038}_{-0.025}$ $R_B$=$0.411 ^{+0.027}_{-0.039}$, respectively. Comparisons with stellar models for low mass stars show that one star is consistent with model predictions whereas the other is substantially oversized. Spectral analysis of the system suggests a primary of spectral type M3V, consistent with both modelled masses and radii, and with SED fitting of NGTS photometry. As the most eccentric eclipsing M-dwarf binary known, NGTS J214358.5-380102 provides an interesting insight into the strength of tidal effects in the circularisation of stellar orbits.

The Transiting Exoplanet Survey Satellite (TESS) has produced a large number of single transit event candidates which are being monitored by the Next Generation Transit Survey (NGTS). We observed a second epoch for the TIC-231005575 system (Tmag = 12.06, Teff = 5500 +- 85 K) with NGTS and a third epoch with Las Cumbres Observatory's (LCO) telescope in South Africa to constrain the orbital period (P = 61.777 d). Subsequent radial velocity measurements with CORALIE revealed the transiting object has a mass of M2 = 0.128 +- 0.003 M$_\odot$, indicating the system is a G-M binary. The radius of the secondary is R2 = 0.154 +- 0.008 R$_\odot$ and is consistent with models of stellar evolution to better than 1-$\sigma$.

We report the detection of the first circumbinary planet found by TESS. The target, a known eclipsing binary, was observed in sectors 1 through 12 at 30-minute cadence and in sectors 4 through 12 at two-minute cadence. It consists of two stars with masses of 1.1 MSun and 0.3 MSun on a slightly eccentric (0.16), 14.6-day orbit, producing prominent primary eclipses and shallow secondary eclipses. The planet has a radius of ~6.9 REarth and was observed to make three transits across the primary star of roughly equal depths (~0.2%) but different durations -- a common signature of transiting circumbinary planets. Its orbit is nearly circular (e ~ 0.09) with an orbital period of 95.2 days. The orbital planes of the binary and the planet are aligned to within ~1 degree. To obtain a complete solution for the system, we combined the TESS photometry with existing ground-based radial-velocity observations in a numerical photometric-dynamical model. The system demonstrates the discovery potential of TESS for circumbinary planets, and provides further understanding of the formation and evolution of planets orbiting close binary stars.

We observed a transit of WASP-166 b using nine NGTS telescopes simultaneously with TESS observations of the same transit. We achieved a photometric precision of 152 ppm per 30 minutes with the nine NGTS telescopes combined, matching the precision reached by TESS for the transit event around this bright (T=8.87) star. The individual NGTS light curve noise is found to be dominated by scintillation noise and appears free from any time-correlated noise or any correlation between telescope systems. We fit the NGTS data for $T_C$ and $R_p/R_*$. We find $T_C$ to be consistent to within 0.25$\sigma$ of the result from the TESS data, and the difference between the TESS and NGTS measured $R_p/R_*$ values is 0.9$\sigma$. This experiment shows that multi-telescope NGTS photometry can match the precision of TESS for bright stars, and will be a valuable tool in refining the radii and ephemerides for bright TESS candidates and planets. The transit timing achieved will also enable NGTS to measure significant transit timing variations in multi-planet systems.

We report the period, eccentricity, and mass determination for the TESS single-transit event candidate TOI-222, which displayed a single 3000 ppm transit in the TESS two-minute cadence data from Sector 2. We determine the orbital period via radial velocity measurements (P=33.9,days), which allowed for ground-based photometric detection of two subsequent transits. Our data show that the companion to TOI-222 is a low mass star, with a radius of $0.18_{-0.10}^{+0.39}$ Rsun and a mass of $0.23\pm0.01$ Msun. This discovery showcases the ability to efficiently discover long-period systems from TESS single transit events using a combination of radial velocity monitoring coupled with high precision ground-based photometry.

We report the discovery, by the Next Generation Transit Survey (NGTS), of two hot-Jupiters NGTS-8b and NGTS-9b. These orbit a V = 13.68 K0V star (Teff = 5241 +/- 50 K) with a period of 2.49970 days, and a V = 12.80 F8V star (Teff = 6330 +/- 130 K) in 4.43527 days, respectively. The transits were independently verified by follow-up photometric observations with the SAAO 1.0-m and Euler telescopes, and we report on the planetary parameters using HARPS, FEROS and CORALIE radial velocities. NGTS-8b has a mass, 0.93 +0.04 -0.03 MJ and a radius, 1.09 +/- 0.03 RJ similar to Jupiter, resulting in a density of 0.89 +0.08 -0.07 g cm-3. This is in contrast to NGTS-9b, which has a mass of 2.90 +/- 0.17 MJ and a radius of 1.07 +/- 0.06 RJ , resulting in a much greater density of 2.93 +0.53 -0.49 g cm-3. Statistically, the planetary parameters put both objects in the regime where they would be expected to exhibit larger than predicted radii. However, we find that their radii are in agreement with predictions by theoretical non-inflated models.

The hot-Neptune HD219666b (= TOI-118.01) was an early discovery from Sector 1 of the TESS transit survey. Being a rare "Neptune desert" planet, and transiting a bright, V = 9.9 star, it is a prime target for further study including atmospheric characterisation. Targeting exoplanet transits with major facilities such as HST and the imminent JWST depends on an accurate ephemeris. Here we report a recovery of the transit in WASP-South data, enabling an ephemeris from data spanning 8 years that has a period uncertainty which is smaller by a factor 60 compared to the discovery paper.

The Transiting Exoplanet Survey Satellite (\tess) produces a large number of single-transit event candidates, since the mission monitors most stars for only $\sim$27\,days. Such candidates correspond to long-period planets or eclipsing binaries. Using the \tess\ Sector 1 full-frame images, we identified a 7750\,ppm single-transit event with a duration of 7\,hours around the moderately evolved F-dwarf star \tic\ (Tmag=10.23, \teff=6280$\pm{85}$\,K). Using archival WASP photometry we constrained the true orbital period to one of three possible values. We detected a subsequent transit-event with NGTS, which revealed the orbital period to be 38.20\,d. Radial velocity measurements from the CORALIE Spectrograph show the secondary object has a mass of $M_2$= $0.148\pm{0.003}$\,M$_{\odot}$, indicating this system is an F-M eclipsing binary. The radius of the M-dwarf companion is $R_2$ = $0.171\pm{0.003}$\,R$_{\odot}$, making this one of the most well characterised stars in this mass regime. We find that its radius is 2.3-$\sigma$ lower than expected from stellar evolution models.

Transits in the planetary system WASP-4 were recently found to occur 80s earlier than expected in observations from the TESS satellite. We present 22 new times of mid-transit that confirm the existence of transit timing variations, and are well fitted by a quadratic ephemeris with period decay dP/dt = -9.2 +/- 1.1 ms/yr. We rule out instrumental issues, stellar activity and the Applegate mechanism as possible causes. The light-time effect is also not favoured due to the non-detection of changes in the systemic velocity. Orbital decay and apsidal precession are plausible but unproven. WASP-4b is only the third hot Jupiter known to show transit timing variations to high confidence. We discuss a variety of observations of this and other planetary systems that would be useful in improving our understanding of WASP-4 in particular and orbital decay in general.

We report the discovery of a new ultra-short period hot Jupiter from the Next Generation Transit Survey. NGTS-6b orbits its star with a period of 21.17~h, and has a mass and radius of $1.330^{+0.024}_{-0.028}$\mjup\, and $1.271^{+0.197}_{-0.188}$\rjup\, respectively, returning a planetary bulk density of 0.711$^{+0.214}_{-0.136}$~g~cm$^{-3}$. Conforming to the currently known small population of ultra-short period hot Jupiters, the planet appears to orbit a metal-rich star ([Fe/H]$=+0.11\pm0.09$~dex). Photoevaporation models suggest the planet should have lost 5\% of its gaseous atmosphere over the course of the 9.6~Gyrs of evolution of the system. NGTS-6b adds to the small, but growing list of ultra-short period gas giant planets, and will help us to understand the dominant formation and evolutionary mechanisms that govern this population.

We present the discovery of NGTS-7Ab, a high mass brown dwarf transiting an M dwarf with a period of 16.2 hours, discovered as part of the Next Generation Transit Survey (NGTS). This is the shortest period transiting brown dwarf around a main or pre-main sequence star to date. The M star host (NGTS-7A) has an age of roughly 55 Myr and is in a state of spin-orbit synchronisation, which we attribute to tidal interaction with the brown dwarf acting to spin up the star. The host star is magnetically active and shows multiple flares across the NGTS and follow up lightcurves, which we use to probe the flare-starspot phase relation. The host star also has an M star companion at a separation of 1.13 arcseconds with very similar proper motion and systemic velocity, suggesting the NGTS-7 system is a hierarchical triple. The combination of tidal synchronisation and magnetic braking is expected to drive ongoing decay of the brown dwarf orbit, with a remaining lifetime of only 5-10 Myr.

Vetting of exoplanet candidates in transit surveys is a manual process, which suffers from a large number of false positives and a lack of consistency. Previous work has shown that Convolutional Neural Networks (CNN) provide an efficient solution to these problems. Here, we apply a CNN to classify planet candidates from the Next Generation Transit Survey (NGTS). For training datasets we compare both real data with injected planetary transits and fully-simulated data, as well as how their different compositions affect network performance. We show that fewer hand labelled lightcurves can be utilised, while still achieving competitive results. With our best model, we achieve an AUC (area under the curve) score of $(95.6\pm{0.2})\%$ and an accuracy of $(88.5\pm{0.3})\%$ on our unseen test data, as well as $(76.5\pm{0.4})\%$ and $(74.6\pm{1.1})\%$ in comparison to our existing manual classifications. The neural network recovers 13 out of 14 confirmed planets observed by NGTS, with high probability. We use simulated data to show that the overall network performance is resilient to mislabelling of the training dataset, a problem that might arise due to unidentified, low signal-to-noise transits. Using a CNN, the time required for vetting can be reduced by half, while still recovering the vast majority of manually flagged candidates. In addition, we identify many new candidates with high probabilities which were not flagged by human vetters.

Some M-dwarfs around F-/G-type stars have been measured to be hotter and larger than predicted by stellar evolution models. Inconsistencies between observations and models need addressing with more mass, radius and luminosity measurements of low-mass stars to test and refine evolutionary models. Our aim is to measure the masses, radii and ages of the stars in five low-mass eclipsing binary systems discovered by the WASP survey. We use WASP photometry to establish eclipse-time ephemerides and to obtain initial estimates for the transit depth and width. Radial velocity measurements were simultaneously fitted with follow-up photometry to find the best-fitting orbital solution. This solution was combined with measurements of atmospheric parameters to interpolate evolutionary models and estimate the mass of the primary star, and the mass and radius of the M-dwarf companion. We assess how the best fitting orbital solution changes if an alternative limb-darkening law is used and quantify the systematic effects of unresolved companions. We also gauge how the best-fitting evolutionary model changes if different values are used for the mixing length parameter and helium enhancement. We report the mass and radius of five M-dwarfs and find little evidence of inflation with respect to evolutionary models. The primary stars in two systems are near the ``blue hook'' stage of their post sequence evolution, resulting in two possible solutions for mass and age. We find that choices in helium enhancement and mixing-length parameter can introduce an additional 3-5\,\% uncertainty in measured M-dwarf mass. Unresolved companions can introduce an additional 3-8\% uncertainty in the radius of an M-dwarf, while the choice of limb-darkening law can introduce up to an additional 2\% uncertainty.

Measurements of the physical properties of stars at the lower end of the main sequence are scarce. In this context we report masses, radii and surface gravities of ten very-low-mass stars in eclipsing binary systems, with orbital periods of the order of several days. The objects probe the stellar mass-radius relation in the fully convective regime, $M_\star \lesssim 0.35$ M$_\odot$, down to the hydrogen burning mass-limit, $M_{\mathrm{HB}} \sim 0.07$ M$_\odot$. The stars were detected by the WASP survey for transiting extra-solar planets, as low-mass, eclipsing companions orbiting more massive, F- and G-type host stars. We use eclipse observations of the host stars (TRAPPIST, Leonhard Euler, SPECULOOS telescopes), and radial velocities of the host stars (CORALIE spectrograph), to determine physical properties of the low-mass companions. Companion surface gravities are derived from the eclipse and orbital parameters of each system. Spectroscopic measurements of the host star effective temperature and metallicity are used to infer the host star mass and age from stellar evolution models. Masses and radii of the low-mass companions are then derived from the eclipse and orbital parameters of each system. The objects are compared to stellar evolution models for low-mass stars, to test for an effect of the stellar metallicity and orbital period on the radius of low-mass stars in close binary systems. Measurements are in good agreement with stellar models; an inflation of the radii of low-mass stars with respect to model predictions is limited to 1.6 $\pm$ 1.2% in the fully convective regime. The sample of ten objects indicates a scaling of the radius of low-mass stars with the host star metallicity. No correlation between stellar radii and orbital periods of the binary systems is determined. A combined analysis with comparable objects from the literature is consistent with this result.

The power-2 law, $I_{\lambda}(\mu)=1- c (1 - \mu^{\alpha})$, accurately represents the limb-darkening profile for cool stars. It has been implemented in a few transit models to-date using numerical integration but there is as-yet no implementation of the power-2 law in analytic form that is generally available. An algorithm to implement the power-2 law is derived using a combination of an approximation to the required integral and a Taylor expansion of the power-2 law. The accuracy of stellar and planetary radii derived by fitting transit light curves with this approximation is tested using light curves computed by numerical integration of limb-darkening profiles from 3D stellar model atmospheres. Our algorithm (qpower2) is accurate to about 100 ppm for broad-band optical light curves of systems with a star-planet radius ratio p = 0.1. The implementation requires less than 40 lines of python code so can run extremely fast on graphical processing units (GPUs; $\sim$ 1 million models per second for the analysis of 1000 data points). Least-squares fits to simulated light curves show that the star and planet radius are recovered to better than 1% for p < 0.2. The qpower2 algorithm can be used to efficiently and accurately analyse large numbers of high-precision transit light curves using Monte Carlo methods.

The Extremely Low-Resource Optical Identifier (ELROI) beacon is a concept for a milliwatt optical "license plate" that can provide unique ID numbers for everything that goes into space. Using photon counting to enable extreme background rejection in real time, the ID number can be uniquely identified from the ground in a few minutes, even if the ground station detects only a few photons per second. The ELROI concept has been validated in long-range ground tests, and orbital prototypes are scheduled for launch in 2018 and beyond. We discuss the design and signal characteristics of these prototypes, including a PC-104 form factor unit which was integrated into a CubeSat and is currently scheduled to launch in May 2018, and basic requirements on ground stations for observing them. We encourage others to consider observing our test flights.

Atmospheric properties of F-,G- and K-type stars can be measured by spectral model fitting. These methods require data with good signal-to-noise ratio and reliable continuum normalisation. This is particularly challenging for the spectra we have obtained with the CORALIE \'{e}chelle spectrograph for FGK stars with transiting M-dwarf companions. The spectra tend to have low signal-to-noise ratios, which makes it difficult to analyse them using existing methods. Our aim is to create a reliable automated spectral analysis routine to determine $T_{\rm eff}$, [Fe/H], $V \sin i$ from the CORALIE spectra of FGK stars. We use wavelet decomposition to distinguish between noise, continuum trends, and stellar spectral features in the CORALIE spectra. A subset of wavelet coefficients from the target spectrum are compared to those from a grid of models in a Bayesian framework to determine the posterior probability distributions of the atmospheric parameters. By testing our method using synthetic spectra we found that our method converges on the best fitting atmospheric parameters. We test the wavelet method on 20 FGK exoplanet host stars for which higher quality data have been independently analysed using equivalent width measurements. We find that we can determine $T_{\rm eff}$ to a precision of $85$K, [Fe/H] to a precision of 0.06 dex and $V \sin i$ to a precision of 1.35 kms$^{-1}$ for stars with $V \sin i$ $\geq$ 5 kms$^{-1}$. We find an offset in metallicity $\approx -$0.18 dex relative to the equivelent width fitting method. We can determine $\log g$ to a precision of $0.13$ dex but measurements are only reliable to confirm dwarf-like surface gravity. The wavelet method can be used to determine $T_{\rm eff}$, [Fe/H] and $V \sin i$ for FGK stars from \'{e}chelle spectra. We find that our method is self consistent, and robust for spectra with SNR$\gtrapprox 40$.

We report the discovery of an eclipsing binary system with mass-ratio q $\sim$ 0.07. After identifying a periodic photometric signal received by WASP, we obtained CORALIE spectroscopic radial velocities and follow-up light curves with the Euler and TRAPPIST telescopes. From a joint fit of these data we determine that EBLM J0555-57 consists of a sun-like primary star that is eclipsed by a low-mass companion, on a weakly eccentric 7.8-day orbit. Using a mass estimate for the primary star derived from stellar models, we determine a companion mass of $85 \pm 4 M_{\rm Jup}$ ($0.081M_{\odot}$) and a radius of $0.84^{+0.14}_{-0.04} R_{\rm Jup}$ ($0.084 R_{\odot}$) that is comparable to that of Saturn. EBLM J0555-57Ab has a surface gravity $\log g_\mathrm{2} = 5.50^{+0.03}_{-0.13}$ and is one of the densest non-stellar-remnant objects currently known. These measurements are consistent with models of low-mass stars.

Observational evidence for the radial alignment of satellites with their dark matter host has been accumulating steadily in the past few years. The effect is seen over a wide range of scales, from massive clusters of galaxies down to galaxy-sized systems, yet the underlying physical mechanism has still not been established. To this end, we have carried out a detailed analysis of the shapes and orientations of dark matter substructures in high-resolution N-body cosmological simulations. We find a strong tendency for radial alignment of the substructure with its host halo: the distribution of halo major axes is very anisotropic, with the majority pointing towards the center of mass of the host. The alignment peaks once the sub-halo has passed the virial radius of the host for the first time, but is not subsequently diluted, even after the halos have gone through as many as four pericentric passages. This evidence points to the existence of a very rapid dynamical mechanism acting on these systems and we argue that tidal torquing throughout their orbits is the most likely candidate.

The cold dark matter (CDM) scenario generically predicts the existence of triaxial dark matter haloes which contain notable amounts of substructure. However, analytical halo models with smooth, spherically symmetric density profiles are routinely adopted in the modelling of light propagation effects through such objects. In this paper, we address the biases introduced by this procedure by comparing the surface mass densities of actual N-body haloes against the widely used analytical model suggested by Navarro, Frenk and White (1996) (NFW). We conduct our analysis in the redshift range of 0.0 - 1.5. In cluster sized haloes, we find that triaxiality can cause scatter in the surface mass density of the haloes up to sigma_+ = +60% and sigma_- = -70%, where the 1-sigma limits are relative to the analytical NFW model given value. Subhaloes can increase this scatter to sigma_+ = +70% and sigma_- = -80%. In galaxy sized haloes, the triaxial scatter can be as high as sigma_+ = +80% and sigma_- = -70%, and with subhaloes the values can change to sigma_+ = +40% and sigma_- = -80%. We present an analytical model for the surface mass density scatter as a function of distance to the halo centre, halo redshift and halo mass. The analytical description enables one to investigate the reliability of results obtained with simplified halo models. Additionally, it provides the means to add simulated surface density scatter to analytical density profiles. As an example, we discuss the impact of our results on the calculation of microlensing optical depths for MACHOs in CDM haloes.

Using eight dark matter haloes extracted from fully-self consistent cosmological N-body simulations, we perform microlensing experiments. A hypothetical observer is placed at a distance of 8.5 kpc from the centre of the halo measuring optical depths, event durations and event rates towards the direction of the Large Magellanic Cloud. We simulate 1600 microlensing experiments for each halo. Assuming that the whole halo consists of MACHOs, f = 1.0, and a single MACHO mass is m_M = 1.0 Msun, the simulations yield mean values of tau = 4.7^{+5.0}_{-2.2} 10^{-7} and Gamma = 1.6^{+1.3}_{-0.6} 10^{-6} events/star/year. We find that triaxiality and substructure can have major effects on the measured values so that tau and Gamma values of up to three times the mean can be found. If we fit our values of tau and Gamma to the MACHO collaboration observations (Alcock et al. 2000), we find f = 0.23^{+0.15}_{-0.13} and m_M = 0.44^{+0.24}_{-0.16}. Five out of the eight haloes under investigation produce f and m_M values mainly concentrated within these bounds.

We investigate the importance of interactions between dark matter substructures for the mass loss they suffer whilst orbiting within a sample of high resolution galaxy cluster mass Cold Dark Matter haloes formed in cosmological N-body simulations. We have defined a quantitative measure that gauges the degree to which interactions are responsible for mass loss from substructures. This measure indicates that interactions are more prominent in younger systems when compared to older more relaxed systems. We show that this is due to the increased number of encounters a satellite experiences and a higher mass fraction in satellites. This is in spite of the uniformity in the distributions of relative distances and velocities of encounters between substructures within the different host systems in our sample. Using a simple model to relate the net force felt by a single satellite to the mass loss it suffers, we show that interactions with other satellites account for ~30% of the total mass loss experienced over its lifetime. The relation between the age of the host and the importance of interactions increases the scatter about this mean value from ~25% for the oldest to ~45% for the youngest system we have studied. We conclude that satellite interactions play a vital role in the evolution of substructure in dark matter halos and that a significant fraction of the tidally stripped material can be attributed to these interactions.

Seven cosmological hydrodynamic simulations of disk galaxy formation are analyzed to determine the alignment of the disk within the dark matter halo and the internal structure of the halo. We find that the orientation of the outer halo, beyond ~0.1 r_vir, is unaffected by the presence of the disk. In contrast, the inner halo is aligned such that the halo minor axis aligns with the disk axis. The relative orientation of these two regions of the halo are uncorrelated. The alignment of the disk and inner halo appears to take place simultaneously through their joint evolution. The disconnect between these two regions of the halo should be taken into account when modelling tidal streams in the halos of disk galaxies and when calculating intrinsic alignments of disk galaxies based on the properties of dark matter halos.

We present an investigation of satellite galaxies in the outskirts of galaxy clusters taken from a series of high-resolution N-body simulations. We focus on the so-called "backsplash population", i.e. satellite galaxies that once were inside the virial radius of the host but now reside beyond it. We find that this population is significant in number and needs to be appreciated when interpreting the various galaxy morphology environmental relationships and decoupling the degeneracy between nature and nurture. Specifically, we find that approximately half of the galaxies with current clustercentric distance in the interval 1-2 virial radii of the host are backsplash galaxies which once penetrated deep into the cluster potential, with 90% of these entering to within 50% of the virial radius. These galaxies have undergone significant tidal disruption, loosing on average 40% of their mass. This results in a mass function for the backsplash population different to those galaxies infalling for the first time. We further show that these two populations are kinematically distinct and should be observable within existent spectroscopic surveys.

We present a detailed study of the real and integrals-of-motion space distributions of a satellite obtained from a self-consistent high-resolution simulation of a galaxy cluster and re-simulated using various analytical halo potentials. We found that the disrupted satellite appears as a coherent structure in integrals-of-motion space in all models (``live'' and analytical potential) although the distribution is significantly smeared for the live host halo. Further the primary mechanism for this smearing is the mass growth of the host, which changes both the energy and angular momentum of the satellite debris. Hence, this must be considered when searching for (stellar) streams with future observational experiments such as RAVE and GAIA.

We describe our new "MLAPM-halo-finder" (MHF) which is based on the adaptive grid structure of the N-body code MLAPM. We then extend the MHF code in order to track the orbital evolution of gravitationally bound objects through any given cosmological N-body simulation - our so-called "MLAPM-halo-tracker" (MHT). The mode of operation of MHT is demonstrated using a series of eight high-resolution N-body simulations of galaxy clusters. Each of these halos hosts more than one million particles within their virial radii Rvir. We use MHT as well as MHF to follow the temporal evolution of hundreds of individual satellites, and show that the radial distribution of these substructure satellites follows a "universal" radial distribution irrespective of the host halo's environment and formation history. This in fact might pose another problem for simulations of CDM structure formation as there are recent findings by Taylor et al. (2003) that the Milky Way satellites are found preferentially closer to the galactic centre and simulations underestimate the amount of central substructure, respectively. Further, this universal substructure profile is anti-biased with respect to the underlying dark matter profile. Both the halo finder MHF and the halo tracker MHT will become part of the open source MLAPM distribution.

We present results from a series of high-resolution N-body simulations that focus on the formation and evolution of eight dark matter halos, each of order a million particles within the virial radius. We follow the time evolution of hundreds of satellite galaxies with unprecedented time resolution, relating their physical properties to the differing halo environmental conditions. The self-consistent cosmological framework in which our analysis was undertaken allows us to explore satellite disruption within live host potentials, a natural complement to earlier work conducted within static potentials. Our host halos were chosen to sample a variety of formation histories, ages, and triaxialities; despite their obvious differences, we find striking similarities within the associated substructure populations. Namely, all satellite orbits follow nearly the same eccentricity distribution with a correlation between eccentricity and pericentre. We also find that the destruction rate of the substructure population is nearly independent of the mass, age, and triaxiality of the host halo. There are, however, subtle differences in the velocity anisotropy of the satellite distribution. We find that the local velocity bias at all radii is greater than unity for all halos and this increases as we move closer to the halo centre, where it varies from 1.1 to 1.4. For the global velocity bias we find a small but slightly positive bias, although when we restrict the global velocity bias calculation to satellites that have had at least one orbit, the bias is essentially removed.

We present a statistical analysis of the interactions between satellite galaxies in cosmological dark matter halos taken from fully self-consistent high-resolution simulations of galaxy clusters. We show that the number distribution of satellite encounters has a tail that extends to as many as 3-4 encounters per orbit. On average 30% of the substructure population had at least one encounter (per orbit) with another satellite galaxy. However, this result depends on the age of the dark matter host halo with a clear trend for more interactions in younger systems. We also report a correlation between the number of encounters and the distance of the satellites to the centre of the cluster: satellite galaxies closer to the centre experience more interactions. However, this can be simply explained by the radial distribution of the substructure population and merely reflects the fact that the density of satellites is higher in those regions. In order to find substructure galaxies we applied (and present) a new technique based upon the N-body code MLAPM. This new halo finder MHF (MLAPM's-Halo-Finder) acts with exactly the same accuracy as the N-body code itself and is therefore free of any bias and spurious mismatch between simulation data and halo finding precision related to numerical effects.

We present preliminary results from a series of high-resolution N-body simulations that focus on 8 dark matter halos each of order a million particles within the virial radius. We follow the time evolution of hundreds of individually tracked satellite galaxies and relate their physical properties to the differing halo environmental conditions. Our main science driver is to understand how satellite galaxies lose their mass and react to tidal stripping. Unlike previous work our results are performed in a fully self-consistent cosmological context. The preliminary results demonstrate that while environment may vary considerably with respects to formation time and richness of substructure, the satellites evolve similarly.

Nearby clusters such as Virgo and Coma possess galaxy distributions which tend to be aligned with the principal axis of the cluster itself. This has also been confirmed by a recent statistical analysis of some 300 Abell clusters where the effect has been linked to the dynamical state of the cluster. Moreover, the orbits of satellite galaxies in galactic systems like our own Milky Way also demonstrate a high degree of anisotropy - the so-called Holmberg effect, the origin of which has been the subject of debate for more than 30 years. This study presents the analysis of cosmological simulations focusing on the orbits of satellite galaxies within dark matter halos. The apocentres of the orbits of these satellites are preferentially found within a cone of opening angle ~40 around the major axis of the host halo, in accordance with the observed anisotropy found in galaxy clusters. We do, however, note that a link to the dynamical age of the cluster is not well established as both our oldest dark matter halos do show a clear anisotropy signal. Further analysis connects this distribution to the infall pattern of satellites along the filaments: the orbits are determined rather by the environment of the host halo than some "dynamical selection" during their life within the host's virial radius.

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