The Arches is one of the youngest, densest and most massive clusters in the Galaxy. As such it provides a unique insight into the lifecycle of the most massive stars known and the formation and survival of such stellar aggregates in the extreme conditions of the Galactic Centre. In a previous study we presented an initial stellar census for the Arches and in this work we expand upon this, providing new and revised classifications for ~30% of the 105 spectroscopically identified cluster members as well as distinguishing potential massive runaways. The results of this survey emphasise the homogeneity and apparent co-evality of the Arches and confirm the absence of H-free Wolf-Rayets of WC sub-type and predicted luminosities. The increased depth of our complete dataset also provides significantly better constraints on the main sequence population; with the identification of O9.5 V stars for the first time we now spectroscopically sample stars with initial masses ranging from ~16Msun to >120Msun. Indeed, following from our expanded stellar census we might expect >50 stars within the Arches to have been born with masses >60Msun, while all 105 spectroscopically confirmed cluster members are massive enough to leave relativistic remnants upon their demise. Moreover the well defined observational properties of the main sequence cohort will be critical to the construction of an extinction law appropriate for the Galactic Centre and consequently the quantitative analysis of the Arches population and subsequent determination of the cluster initial mass function.

Observational studies suggest that many OB stars are found within binary systems which may be expected to interact during their lifetimes. Significant mass transfer or merger of both components will modify evolutionary pathways, facilitating the production of exceptionally massive stars which will present as blue stragglers. Identification and characterisation of such objects is crucial if the outcomes of binary evolutionary channels are to be quantified. The massive cluster Westerlund 1 hosts a rich population of X-ray bright stars where the emission is thought to derive from the wind collision zones of massive binaries. Selected on this basis, we present the results of a multiwavelength analysis of the X-ray luminous O stars Wd1-27 and -30a. We find both to be early/mid-O hypergiants with luminosities, temperatures and masses significantly in excess of other early stars within Wd1, hence qualifying as massive blue stragglers. The nature of Wd1-27 remains uncertain but the detection of radial velocity changes and the X-ray properties of Wd1-30a suggest that it is a binary. Analysis of Gaia proper motion and parallactic data indicates that both are cluster members; we also provide a membership list for Wd1 based on this analysis. The presence of hypergiants of spectral types O to M within Wd1 cannot be understood via single-star evolution. We suppose that the early-B and mid-O hypergiants formed via binary-induced mass-stripping of the primary and mass-transfer to the secondary, respectively. This implies that for a subset of objects massive star-formation may be regarded as a two-stage process, with binary-driven mass-transfer or merger yielding stars with masses significantly in excess of their initial birth mass (Abridged).

The Quintuplet is one of the most massive galactic clusters known, but appears to host a diverse stellar population that is difficult to reconcile with an instantaneous formation event. We present HST photometry and VLT spectroscopy in order to improve observational constraints, finding the Quintuplet to be far more homogeneous than previously thought. O7-8 Ia and O9-B0 Ia supergiants form a smooth morphological sequence with a cohort of seven early-B hypergiants and six luminous blue variables and WN9-11h stars, which together comprise the richest population of such stars of any known stellar aggregate. No further H-free Wolf-Rayets were found, resulting in a 13:1 ratio for WC/WN stars. However a small population of late-O hypergiants and WN8-9ha stars was identified, while a subset of the supergiants are unexpectedly faint, suggesting they are both less massive and older than the greater cluster population. We find an impressive coincidence between cluster members preceding the H-free WR phase and the evolutionary predictions for a 60Msun star, suggesting an age of ~3.0-3.6Myr for the Quintuplet. Neither the late-O hypergiants nor the low luminosity supergiants are predicted; we suggest that the former either result from rapid rotators or are the products of binary driven mass-stripping, while the latter may be interlopers. The H-free WRs must have had initial masses >60Msun, but it is difficult to reconcile their observational properties with theoretical expectations. Since the WRs represent an evolutionary phase directly preceding core-collapse, they are are crucial to understanding both this process and the nature of the resultant relativistic remnants, providing unique constraints on the evolution and death of the most massive stars forming in the local, high metallicity Universe (Abridged).

As part of a wider investigation of evolved massive stars in Galactic open clusters, we have spectroscopically identified three candidate classical Cepheids in the little-studied clusters Berkeley 51, Berkeley 55 and NGC 6603. Using new multi-epoch photometry, we confirm that Be 51 #162 and Be 55 #107 are bona fide Cepheids, with pulsation periods of 9.83+/-0.01 d and 5.850+/-0.005 d respectively, while NGC 6603 star W2249 does not show significant photometric variability. Using the period-luminosity relationship for Cepheid variables, we determine a distance to Be 51 of 5.3(+1.0,-0.8) kpc and an age of 44(+9,-8) Myr, placing it in a sparsely-attested region of the Perseus arm. For Be 55, we find a distance of 2.2+/-0.3 kpc and age of 63(+12,-11) Myr, locating the cluster in the Local arm. Taken together with our recent discovery of a long-period Cepheid in the starburst cluster VdBH222, these represent an important increase in the number of young, massive Cepheids known in Galactic open clusters. We also consider new Gaia (data release 2) parallaxes and proper motions for members of Be 51 and Be 55; the uncertainties on the parallaxes do not allow us to refine our distance estimates to these clusters, but the well-constrained proper motion measurements furnish further confirmation of cluster membership. However, future final Gaia parallaxes for such objects should provide valuable independent distance measurements, improving the calibration of the period-luminosity relationship, with implications for the distance ladder out to cosmological scales.

We have carried out a spectroscopic variability survey of some of the most massive stars in the Arches cluster, using K-band observations obtained with SINFONI on the VLT. One target, F2, exhibits substantial changes in radial velocity; in combination with new KMOS and archival SINFONI spectra, its primary component is found to undergo radial velocity variation with a period of 10.483+/-0.002 d and an amplitude of ~350 km/s. A secondary radial velocity curve is also marginally detectable. We reanalyse archival NAOS-CONICA photometric survey data in combination with our radial velocity results to confirm this object as an eclipsing SB2 system, and the first binary identified in the Arches. We model it as consisting of an 82+/-12 M_sun WN8-9h primary and a 60+/-8 M_sun O5-6 Ia+ secondary, and as having a slightly eccentric orbit, implying an evolutionary stage prior to strong binary interaction. As one of four X-ray bright Arches sources previously proposed as colliding-wind massive binaries, it may be only the first of several binaries to be discovered in this cluster, presenting potential challenges to recent models for the Arches' age and composition. It also appears to be one of the most massive binaries detected to date; the primary's calculated initial mass of >~120 M_sun would arguably make this the most massive binary known in the Galaxy.

Located within the central region of the Galaxy, the Arches cluster appears to be one of the youngest, densest and most massive stellar aggregates within the Milky Way. As such it has the potential to be a uniquely instructive laboratory for the study of star formation in extreme environments and the physics of very massive stars. In order to determine the fundamental physical properties of both cluster and constituent stars, we provide and analyse new HST+VLT near-IR datasets. Stacking multiple epochs of spectroscopy results in the deepest view of the cluster ever obtained, allowing us to to identify candidate giant and main sequence stars for the first time. All cluster members are found to be WNLh or O stars, with the smooth and continuous progression in spectral morphologies from O super-/hypergiants through to the WNLh cohort implying a direct evolutionary connection. Importantly no H-free Wolf-Rayets are found, while no products of binary interaction/mass-transfer may be unambiguously identified, despite the presence of massive binaries within the Arches. We infer a main sequence turn-off around O4-5V, corresponding to ~30-38Msun, while the eclipsing binary F2 implies current masses of ~80Msun and ~60Msun for the WNLh and O hypergiant cohorts, respectively. A cluster age of ~2-3Myr is suggested by the location of the main-sequence turn-off. While the absence of H-free Wolf-Rayets argues against the prior occurrence of SNe, such an age does accommodate such events for exceptionally massive stars. Future progress requires quantitative analysis of cluster members combined with additional spectroscopic observations to better constrain the binary population; nevertheless it is already abundantly clear that the Arches offers an unprecedented insight into the formation, evolution and death of the most massive stars Nature allows to form in the local universe (Abridged).

We searched the SuperWASP archive for objects displaying multiply periodic photometric variations. Specifically we sought evidence for eclipsing binary stars displaying a further non-harmonically related signal in their power spectra. The object SWASP J050634.16-353648.4 is identified as a relatively bright (V ~ 11.5) semi-detached eclipsing binary with a 5.104 d orbital period that displays coherent pulsations with a semi-amplitude of 65 mmag at a frequency of 13.45 per day. Follow-up radial velocity spectroscopy with the Southern African Large Telescope confirmed the binary nature of the system. Using the phoebe code to model the radial velocity curve with the SuperWASP photometry allowed parameters of both stellar components to be determined. This yielded a primary (pulsating) star with a mass of 1.73 +/- 0.11 solar mass and a radius of 2.41 +/- 0.06 solar radii plus a Roche-lobe filling secondary star with a mass of 0.41 +/- 0.03 solar mass and a radius of 4.21 +/- 0.11 solar radii. 1SWASP J050634.16-353648.4 is therefore a bright delta Sct pulsator in a semi-detached eclipsing binary with one of the largest pulsation amplitudes of any such system known. The pulsation constant indicates that the mode is likely a first overtone radial pulsation.

Galactic starburst clusters play a twin role in astrophysics, serving as laboratories for the study of stellar physics and also delineating the structure and recent star formation history of the Milky Way. In order to exploit these opportunities we have undertaken a multi-epoch spectroscopic survey of the red supergiant dominated young massive clusters thought to be present at both near and far ends of the Galactic Bar. Significant spectroscopic variability suggestive of radial pulsations was found for the yellow supergiant VdBH 222 #505. Follow-up photometric investigations revealed modulation with a period of ~23.325d; both timescale and pulsational profile are consistent with a Cepheid classification. As a consequence #505 may be recognised as one of the longest period Galactic cluster Cepheids identified to date and hence of considerable use in constraining the bright end of the period/luminosity relation at solar metallicities. In conjunction with extant photometry we infer a distance of ~6kpc for VdBH222 and an age of ~20Myr. This results in a moderate reduction in both integrated cluster mass (~2x10^4Msun) and the initial stellar masses of the evolved cluster members (~10Msun). As such, VdBH222 becomes an excellent test-bed for studying the properties of some of the lowest mass stars observed to undergo type-II supernovae. Moreover, the distance is in tension with a location of VdBH 222 at the far end of the Galactic Bar. Instead a birthsite in the near 3kpc arm is suggested; providing compelling evidence of extensive recent star formation in a region of the inner Milky Way which has hitherto been thought to be devoid of such activity.

Since M4.5 - M9 dwarfs exhibit equatorial rotation velocities of order 10 km/s on average, radial velocity surveys targeting this stellar population will likely need to find methods to effectively remove starspot jitter. We present the first high resolution Doppler images of the M4.5 dwarf, GJ 791.2A, and the M9 dwarf, LP 944-20. The time series spectra of both objects reveal numerous line profile distortions over the rotation period of each star which we interpret as starspots. The transient distortions are modelled with spot/photosphere contrast ratios that correspond to model atmosphere temperature differences of Tphot-Tspot = 200 and 300 K. GJ 791.2A is a fully convective star with vsini = 35.1 km/s. Although we find more starspot structure at high latitudes, we reconstruct spots at a range of phases and latitudes with a mean spot filling of ~3%. LP 944-20 is one of the brightest known late-M dwarfs, with spectral type M9V and vsini = 30.8 km/s. Its spectral time series exhibits two dominant transient line distortions that are reconstructed as high latitude spots, while a mean spot filling factor of only 1.5% is found. The occurrence of low-contrast spots at predominantly high latitudes, which we see in both targets here, is in general likely to be responsible for the low amplitude photometric variability seen in late-M dwarfs. For GJ 791.2A, the radial velocities induced by the starspot features yield an rms velocity variability of 138 m/s, which can be reduced by a factor of 1.9 using our reconstructed surface brightness distributions.

Orbital period changes of binary stars may be caused by the presence of a third massive body in the system. Here we have searched the archive of the Wide Angle Search for Planets (SuperWASP) project for evidence of period variations in 13927 eclipsing binary candidates. Sinusoidal period changes, strongly suggestive of third bodies, were detected in 2% of cases; however, linear period changes were observed in a further 22% of systems. We argue on distributional grounds that the majority of these apparently linear changes are likely to reflect longer-term sinusoidal period variations caused by third bodies, and thus estimate a higher-order multiplicity fraction of 24% for SuperWASP binaries, in good agreement with other recent figures for the fraction of triple systems amongst binary stars in general.

The Super Wide Angle Search for Planets (SuperWASP) is a whole-sky high-cadence optical survey which has searched for exoplanetary transit signatures since 2004. Its archive contains long-term light curves for ~30 million 8-15 V magnitude stars, making it a valuable serendipitous resource for variable star research. We have concentrated on the evidence it provides for eclipsing binaries, in particular those exhibiting orbital period variations, and have developed custom tools to measure periods precisely and detect period changes reliably. Amongst our results are: a collection of 143 candidate contact or semi-detached eclipsing binaries near the short-period limit in the main sequence binary period distribution; a probable hierarchical triple exhibiting dramatic sinusoidal period variations; a new doubly-eclipsing quintuple system; and new evidence for period change or stability in 12 post-common-envelope eclipsing binaries, which may support the existence of circumbinary planets in such systems. A large-scale search for period changes in ~14000 SuperWASP eclipsing binary candidates also yields numerous examples of sinusoidal period change, suggestive of tertiary companions, and may allow us to constrain the frequency of triple systems amongst low-mass stars.

Our discovery of 1SWASP J093010.78+533859.5 as a probable doubly eclipsing quadruple system containing a contact binary with P~0.23 d and a detached binary with P~1.31 d was announced in 2013. Subsequently Koo et al. confirmed the detached binary spectroscopically and identified a fifth set of static spectral lines at its location, corresponding to a further non-eclipsing component of the system. Here we present new spectroscopic and photometric observations, allowing confirmation of the contact binary and improved modelling of all four eclipsing components. The detached binary is found to contain components of masses 0.837(8) and 0.674(7) M_sol, with radii of 0.832(18) and 0.669(18) R_sol and effective temperatures of 5185(-20,+25) and 4325(-15,+20) K respectively, the contact system has masses 0.86(2) and 0.341(11) M_sol, radii of 0.79(4) and 0.52(5) R_sol respectively, and a common T_eff of 4700(50) K. The fifth star is of similar temperature and spectral type to the primaries in the two binaries. Long-term photometric observations indicate the presence of a spot on one component of the detached binary, moving at an apparent rate of approximately one rotation every two years. Both binaries have consistent system velocities around -11 to -12 km/s, which match the average radial velocity of the fifth star, consistent distance estimates for both subsystems of d=78(3) and d=73(4) pc are also found, and (with some further assumptions) of d=83(9) pc for the fifth star. These findings strongly support the claim that both binaries (and very probably all five stars) are gravitationally bound in a single system. The consistent angles of inclination found for the two binaries (88.2(3) and 86(4) degrees) may also indicate that they originally formed by fragmentation (~9-10 Gyr ago) from a single protostellar disk and subsequently remained in the same orbital plane.

Period or amplitude variations in eclipsing binaries may reveal the presence of additional massive bodies in the system, such as circumbinary planets. Here, we have studied twelve previously-known eclipsing post-common-envelope binaries for evidence of such light curve variations, on the basis of multi-year observations in the SuperWASP archive. The results for HW Vir provided strong evidence for period changes consistent with those measured by previous studies, and help support a two-planet model for the system. ASAS J102322-3737.0 exhibited plausible evidence for a period increase not previously suggested; while NY Vir, QS Vir and NSVS 14256825 afforded less significant support for period change, providing some confirmation to earlier claims. In other cases, period change was not convincingly observed; for AA Dor and NSVS 07826147, previous findings of constant period were confirmed. This study allows us to present hundreds of new primary eclipse timings for these systems, and further demonstrates the value of wide-field high-cadence surveys like SuperWASP for the investigation of variable stars.

The two objects 1SWASP J150822.80-054236.9 and 1SWASP J160156.04+202821.6 were initially detected from their SuperWASP archived light curves as candidate eclipsing binaries with periods close to the short-period cut-off of the orbital period distribution of main sequence binaries, at ~0.2 d. Here, using INT spectroscopic data, we confirm them as double-lined spectroscopic and eclipsing binaries, in contact configuration. Following modelling of their visual light curves and radial velocity curves, we determine their component and system parameters to precisions between ~2 and 11%. The former system contains 1.07 and 0.55 M_sun components, with radii of 0.90 and 0.68 R_sun respectively; its primary exhibits pulsations with period 1/6 the orbital period of the system. The latter contains 0.86 and 0.57 M_sun components, with radii of 0.75 and 0.63R_sun respectively.

We report our investigation of 1SWASP J234401.81-212229.1, a variable with a 18461.6 s period. After identification in a 2011 search of the SuperWASP archive for main-sequence eclipsing binary candidates near the distribution's short-period limit of approx. 0.20 d, it was measured to be undergoing rapid period decrease in our earlier work, though later observations supported a cyclic variation in period length. Spectroscopic data obtained in 2012 with the Southern African Large Telescope did not, however, support the interpretation of the object as a normal eclipsing binary. Here, we consider three possible explanations consistent with the data: a single-star oblique rotator model in which variability results from stable cool spots on opposite magnetic poles; a two-star model in which the secondary is a brown dwarf; and a three-star model involving a low-mass eclipsing binary in a hierarchical triple system. We conclude that the latter is the most likely model.

Building on previous work, a new search of the SuperWASP archive was carried out to identify eclipsing binary systems near the short-period limit. 143 candidate objects were detected with orbital periods between 16000 and 20000 s, of which 97 are new discoveries. Period changes significant at 1 sigma or more were detected in 74 of these objects, and in 38 the changes were significant at 3 sigma or more. The significant period changes observed followed an approximately normal distribution with a half-width at half-maximum of ~0.1 s/yr. There was no apparent relationship between period length and magnitude or direction of period change. Amongst several interesting individual objects studied, 1SWASP J093010.78+533859.5 is presented as a new doubly eclipsing quadruple system, consisting of a contact binary with a 19674.575 s period and an Algol-type binary with a 112799.109 s period, separated by 66.1 AU, being the sixth known system of this type.

SuperWASP light curves for 53 W UMa-type eclipsing binary (EB) candidates, identified in previous work as being close to the contact binary short-period limit, were studied for evidence of period change. The orbital periods of most of the stars were confirmed, and period decrease, significant at more than 5 sigma, was observed in three objects: 1SWASP J174310.98+432709.6 (-0.055 \pm0.003 s/yr), 1SWASP J133105.91+121538.0 (-0.075 \pm0.013 s/yr) and 1SWASP J234401.81-212229.1 (-0.313 \pm0.019 s/yr). The magnitudes of the observed period changes cannot be explained by magnetic braking or gravitational radiation effects, and are most likely primarily due to unstable mass transfer from primary to secondary components, possibly accompanied by unstable mass and angular momentum loss from the systems. If these period decreases persist, the systems could merge on a relatively short timescale.