Astrochemistry

In the last 40 years a wide range of molecules, including neutrals, cations and anions, containing up to 13 atoms—in addition to detections of C₆₀ and C₇₀ — have been found in the harsh environment of the interstellar medium. The exquisite sensitivity and very high spectral and, more recently, spatial resolution, of modern telescopes has enabled the physics of star formation to be probed through rotational line emission. In this article, I review the basic properties of interstellar clouds and the processes that initiate the chemistry and generate chemical complexity, particularly in regions of star and planet formation. Our understanding of astrochemistry has evolved over the years. Before 1990, the general consensus was that molecules were formed in binary, gas-phase, or volume, reactions, most importantly ion-neutral reactions despite the very low ionization in clouds. Since then, observations have indicated unambiguously that there is also a contribution from surface processes, particularly on the icy mantles that form around refractory grain cores in cold, dense gas. The balance between these two processes depends on particular physical conditions and can vary during the life cycle of a particular volume of interstellar cloud.The complex chemistry that occurs in space is driven mostly through interaction of thegas with cosmic ray protons, a source of ionization that enables a rich ion-neutral chemistry. In addition, I show that the interaction between the gas and the dust in cold, dense regionsalso leads to additional chemical complexity through reactions that take place in ices at onlya few tens of degrees above absolute zero. Although densities are low compared to those in terrestrial environments, the extremely long life times of interstellar clouds and their enormous sizes, enable complex molecules to be synthesised and detected. I show that in some instances, particularly in reactions involving deuterium, the rotational populations of reactants, together with spin-selection rules, can determine the detailed abundances. Although the review is mainly focused on regions associated with star formation, I also consider chemistry in other interesting astronomical regions — in the early Universe and in the envelopes formed by mass loss during the final stages of stellar evolution.

The Ultraviolet-Bright, Slowly Declining Transient PS1-11af as a Partial Tidal Disruption Event

We present the Pan-STARRS1 discovery of the long-lived and blue transient PS1-11af, which was also detected by Galaxy Evolution Explorer with coordinated observations in the near-ultraviolet (NUV) band. PS1-11af is associated with the nucleus of an early type galaxy at redshift z = 0.4046 that exhibits no evidence for star formation or active galactic nucleus activity. Four epochs of spectroscopy reveal a pair of transient broad absorption features in the UV on otherwise featureless spectra. Despite the superficial similarity of these features to P-Cygni absorptions of supernovae (SNe), we conclude that PS1-11af is not consistent with the properties of known types of SNe. Blackbody fits to the spectral energy distribution are inconsistent with the cooling, expanding ejecta of a SN, and the velocities of the absorption features are too high to represent material in homologous expansion near a SN photosphere. However, the constant blue colors and slow evolution of the luminosity are similar to previous optically selected tidal disruption events (TDEs). The shape of the optical light curve is consistent with models for TDEs, but the minimum accreted mass necessary to power the observed luminosity is only 0.002 M, which points to a partial disruption model. A full disruption model predicts higher bolometric luminosities, which would require most of the radiation to be emitted in a separate component at high energies where we lack observations. In addition, the observed temperature is lower than that predicted by pure accretion disk models for TDEs and requires reprocessing to a constant, lower temperature. Three deep non-detections in the radio with the Very Large Array over the first two years after the event set strict limits on the production of any relativistic outflow comparable to Swift J1644+57, even if off-axis.

Hydrogen-Poor Superluminous Supernovae and Long-Duration Gamma-Ray Bursts Have Similar Host Galaxies

We present optical spectroscopy and optical/near-IR photometry of 31 host galaxies of hydrogen-poor superluminous supernovae (SLSNe), including 15 events from the Pan-STARRS1 Medium Deep Survey. Our sample spans the redshift range 0.1 ≲ z ≲ 1.6, and is the first comprehensive host galaxy study of this specific subclass of cosmic explosions. Combining the multi-band photometry and emission-line measurements, we determine the luminosities, stellar masses, star formation rates, and metallicities. We find that, as a whole, the hosts of SLSNe are a low-luminosity (〈M_B 〉 ≈ -17.3 mag), low stellar mass (〈M〉 ≈ 2 × 10⁸ M_⊙) population, with a high median specific star formation rate (〈sSFR〉 ≈ 2 Gyr^-1). The median metallicity of our spectroscopic sample is low, 12 + log (O/H) ≈ 8.35 ≈ 0.45 Z_⊙, although at least one host galaxy has solar metallicity. The host galaxies of H-poor SLSNe are statistically distinct from the hosts of GOODS core-collapse SNe (which cover a similar redshift range), but resemble the host galaxies of long-duration gamma-ray bursts (LGRBs) in terms of stellar mass, SFR, sSFR, and metallicity. This result indicates that the environmental causes leading to massive stars forming either SLSNe or LGRBs are similar, and in particular that SLSNe are more effectively formed in low metallicity environments. We speculate that the key ingredient is large core angular momentum, leading to a rapidly spinning magnetar in SLSNe and an accreting black hole in LGRBs.

PS1-10afx at z = 1.388: Pan-STARRS1 Discovery of a New Type of Superluminous Supernova

We present the Pan-STARRS1 discovery of PS1-10afx, a unique hydrogen-deficient superluminous supernova (SLSN) at redshift z = 1.388. The light curve peaked at z P1 = 21.7 mag, making PS1-10afx comparable to the most luminous known SNe, with Mu = -22.3 mag. Our extensive optical and near-infrared observations indicate that the bolometric light curve of PS1-10afx rose on the unusually fast timescale of ~12 days to the extraordinary peak luminosity of 4.1 × 1044 erg s-1 (M bol = -22.8 mag) and subsequently faded rapidly. Equally important, the spectral energy distribution is unusually red for an SLSN, with a color temperature of ~6800 K near maximum light, in contrast to previous hydrogen-poor SLSNe, which are bright in the ultraviolet (UV). The spectra more closely resemble those of a normal SN Ic than any known SLSN, with a photospheric velocity of ~11, 000 km s-1 and evidence for line blanketing in the rest-frame UV. Despite the fast rise, these parameters imply a very large emitting radius (gsim 5 × 1015 cm). We demonstrate that no existing theoretical model can satisfactorily explain this combination of properties: (1) a nickel-powered light curve cannot match the combination of high peak luminosity with the fast timescale; (2) models powered by the spindown energy of a rapidly rotating magnetar predict significantly hotter and faster ejecta; and (3) models invoking shock breakout through a dense circumstellar medium cannot explain the observed spectra or color evolution. The host galaxy is well detected in pre-explosion imaging with a luminosity near L*, a star formation rate of ~15 M ⊙ yr-1, and is fairly massive (~2 × 1010 M ⊙), with a stellar population age of ~108 yr, also in contrast to the young dwarf hosts of known hydrogen-poor SLSNe. PS1-10afx is distinct from known examples of SLSNe in its spectra, colors, light-curve shape, and host galaxy properties, suggesting that it resulted from a different channel than other hydrogen-poor SLSNe.

Ultraluminous Supernovae as a New Probe of the Interstellar Medium in Distant Galaxies

We present the Pan-STARRS1 discovery and light curves, and follow-up MMT and Gemini spectroscopy of an ultraluminous supernova (ULSN; dubbed PS1-11bam) at a redshift of z = 1.566 with a peak brightness of M UV ≈ -22.3 mag. PS1-11bam is one of the highest redshift spectroscopically confirmed SNe known to date. The spectrum exhibits broad absorption features typical of previous ULSNe (e.g., C II, Si III), and strong and narrow Mg II and Fe II absorption lines from the interstellar medium (ISM) of the host galaxy, confirmed by an [O II]λ3727 emission line at the same redshift. The equivalent widths of the Fe II λ2600 and Mg II λ2803 lines are in the top quartile of the quasar intervening absorption system distribution, but are weaker than those of gamma-ray burst intrinsic absorbers (i.e., GRB host galaxies). We also detect the host galaxy in pre-explosion Pan-STARRS1 data and find that its UV spectral energy distribution is best fit with a young stellar population age of τ* ≈ 15-45 Myr and a stellar mass of M * ≈ (1.1-2.6) × 109 M ⊙ (for Z = 0.05-1 Z ⊙). The star formation rate inferred from the UV continuum and [O II]λ3727 emission line is ≈10 M ⊙ yr-1, higher than in previous ULSN hosts. PS1-11bam provides the first direct demonstration that ULSNe can serve as probes of the ISM in distant galaxies. The depth and red sensitivity of PS1 are uniquely suited to finding such events at cosmologically interesting redshifts (z ~ 1-2); the future combination of LSST and 30 m class telescopes promises to extend this technique to z ~ 4.

The host galaxy and late-time evolution of the superluminous supernova PTF12dam

Superluminous supernovae (SLSNe) of Type Ic have a tendency to occur in faint host galaxies which are likely to have low mass and low metallicity. PTF12dam is one of the closest and best-studied superluminous explosions that has a broad and slowly fading light curve similar to SN 2007bi. Here we present new photometry and spectroscopy for PTF12dam from 200-500 d (rest frame) after peak and a detailed analysis of the host galaxy (SDSS J142446.21+461348.6 at z = 0.107). Using deep templates and image subtraction we show that the light curve can be fit with a magnetar model if escape of high-energy gamma rays is taken into account. The full bolometric light curve from -53 to +399 d (with respect to peak) cannot be fit satisfactorily with the pair-instability models. An alternative model of interaction with a dense circumstellar material (CSM) produces a good fit to the data although this requires a very large mass (˜13 M⊙) of hydrogen-free CSM. The host galaxy is a compact dwarf (physical size ˜1.9 kpc) and with Mg = -19.33 ± 0.10, it is the brightest nearby SLSN Ic host discovered so far. The host is a low-mass system (2.8 × 108 M⊙) with a star formation rate (5.0 M⊙ yr-1), which implies a very high specific star formation rate (17.9 Gyr-1). The remarkably strong nebular emission provide detections of the [O III] λ4363 and [O II] λλ7320, 7330auroral lines and an accurate oxygen abundance of 12 + log (O/H) = 8.05 ± 0.09. We show here that they are at the extreme end of the metallicity distribution of dwarf galaxies and propose that low metallicity is a requirement to produce these rare and peculiar SNe.

Locations of peculiar supernovae as a diagnostic of their origins

We put constraints on the properties of the progenitors of peculiar calcium-rich transients using the distribution of locations within their host galaxies. We confirm that this class of transients do not follow the galaxy stellar mass profile and are more likely to be found in remote locations of their apparent hosts. We test the hypothesis that these transients are from low-metallicity progenitors by comparing their spatial distributions with the predictions of self-consistent cosmological simulations that include star formation and chemical enrichment. We find that while metal-poor stars and our transient sample show a consistent preference for large offsets, metallicity alone cannot explain the extreme cases. Invoking a lower age limit on the progenitor helps to improve the match, indicating these events may result from a very old metal-poor population. We also investigate the radial distribution of globular cluster systems, and show that they too are consistent with the class of calcium-rich transients. Because photometric upper limits exist for globular clusters for some members of the class, a production mechanism related to the dense environment of globular clusters is not favoured for the calcium-rich events. However, the methods developed in this paper may be used in the future to constrain the effects of low metallicity on radially distant core-collapse events or help establish a correlation with globular clusters for other classes of peculiar explosions.

The VLT-FLAMES Tarantula Survey XXI. Stellar spin rates of O-type spectroscopic binaries

Context: The initial distribution of spin rates of massive stars is a fingerprint of their elusive formation process. It also sets a key initial condition for stellar evolution and is thus an important ingredient in stellar population synthesis. So far, most studies have focused on single stars. Most O stars are, however, found in multiple systems. 

Aims: By establishing the spin-rate distribution of a sizeable sample of O-type spectroscopic binaries and by comparing the distributions of binary subpopulations with one another and with that of presumed-single stars in the same region, we aim to constrain the initial spin distribution of O stars in binaries, and to identify signatures of the physical mechanisms that affect the evolution of the spin rates of massive stars. 

Methods: We use ground-based optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS) to establish the projected equatorial rotational velocities (ν_esini) for components of 114 spectroscopic binaries in 30 Doradus. The ν_esini values are derived from the full width at half maximum (FWHM) of a set of spectral lines, using a FWHM vs. ν_esini calibration that we derive based on previous line analysis methods applied to single O-type stars in the VFTS sample. 

Results: The overall ν_esini distribution of the primary stars resembles that of single O-type stars in the VFTS, featuring a low-velocity peak (at ν_esini<200 kms^-1) and a shoulder at intermediate velocities (200 <ν_esini<300 kms^-1). The distributions of binaries and single stars, however, differ in two ways. First, the main peak at ν_esini ~ 100kms^-1 is broader and slightly shifted towards higher spin rates in the binary distribution than that of the presumed-single stars. This shift is mostly due to short-period binaries (P_orb~<10 d). Second, the ν_esini distribution of primaries lacks a significant population of stars spinning faster than 300 kms^-1, while such a population is clearly present in the single-star sample. The ν_esini distribution of binaries with amplitudes of radial velocity variation in the range of 20 to 200 kms^-1 (mostly binaries with P_orb ~ 10-1000 d and/or with q<0.5) is similar to that of single O stars below ν_esini_~<170kms^-1. 

Conclusions: Our results are compatible with the assumption that binary components formed with the same spin distribution as single stars, and that this distribution contains few or no fast-spinning stars. The higher average spin rate of stars in short-period binaries may either be explained by spin-up through tides in such tight binary systems, or by spin-down of a fraction of the presumed-single stars and long-period binaries through magnetic braking (or by a combination of both mechanisms). Most primaries and secondaries of SB2 systems with P_orb~<10 d appear to have similar rotational velocities. This is in agreement with tidal locking in close binaries where the components have similar radii. The lack of very rapidly spinning stars among binary systems supports the idea that most stars with ν_esini_~> 300kms^-1 in the single-star sample are actually spun-up post-binary interaction products. Finally, the overall similarities (low-velocity peak and intermediate-velocity shoulder) of the spin distribution of binary and single stars argue for a massive star formation process in which the initial spin is set independently of whether stars are formed as single stars or as components of a binary system.

The VLT-FLAMES tarantula survey:XII. Rotational velocities of the single O-type stars

Context. The 30 Doradus (30 Dor) region of the Large Magellanic Cloud, also known as the Tarantula nebula, is the nearest starburst region. It contains the richest population of massive stars in the Local Group, and it is thus the best possible laboratory to investigate open questions on the formation and evolution of massive stars. Aims. Using ground-based multi-object optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to establish the (projected) rotational velocity distribution for a sample of 216 presumably single O-type stars in 30 Dor. The sample is large enough to obtain statistically significant information and to search for variations among subpopulations - in terms of spectral type, luminosity class, and spatial location - in the field of view. Methods. We measured projected rotational velocities, 3e sin i, by means of a Fourier transform method and a profile fitting method applied to a set of isolated spectral lines. We also used an iterative deconvolution procedure to infer the probability density, P(3e), of the equatorial rotational velocity, 3e. Results. The distribution of 3e sin i shows a two-component structure: a peak around 80 km s1 and a high-velocity tail extending up to 600 km s^-1 This structure is also present in the inferred distribution P(3e) with around 80% of the sample having 0 <3e ≤ 300 km s^-1 and the other 20% distributed in the high-velocity region. The presence of the low-velocity peak is consistent with what has been found in other studies for late O- and early B-type stars. Conclusions. Most of the stars in our sample rotate with a rate less than 20% of their break-up velocity. For the bulk of the sample, mass loss in a stellar wind and/or envelope expansion is not efficient enough to significantly spin down these stars within the first few Myr of evolution. If massive-star formation results in stars rotating at birth with a large portion of their break-up velocities, an alternative braking mechanism, possibly magnetic fields, is thus required to explain the present-day rotational properties of the O-type stars in 30 Dor. The presence of a sizeable population of fast rotators is compatible with recent population synthesis computations that investigate the influence of binary evolution on the rotation rate of massive stars. Even though we have excluded stars that show significant radial velocity variations, our sample may have remained contaminated by post-interaction binary products. That the highvelocity tail may be populated primarily (and perhaps exclusively) by post-binary interaction products has important implications for the evolutionary origin of systems that produce gamma-ray bursts. © 2013 Author(s).

The VLT-FLAMES tarantula survey:X. Evidence for a bimodal distribution of rotational velocities for the single early B-type stars

Aims. Projected rotational velocities (ve sin i) have been estimated for 334 targets in the VLT-FLAMES Tarantula Survey that do not manifest significant radial velocity variations and are not supergiants. They have spectral types from approximately O9.5 to B3. The estimates have been analysed to infer the underlying rotational velocity distribution, which is critical for understanding the evolution of massive stars. Methods. Projected rotational velocities were deduced from the Fourier transforms of spectral lines, with upper limits also being obtained from profile fitting. For the narrower lined stars, metal and non-diffuse helium lines were adopted, and for the broader lined stars, both non-diffuse and diffuse helium lines; the estimates obtained using the different sets of lines are in good agreement. The uncertainty in the mean estimates is typically 4% for most targets. The iterative deconvolution procedure of Lucy has been used to deduce the probability density distribution of the rotational velocities. Results. Projected rotational velocities range up to approximately 450 kms-1 and show a bi-modal structure. This is also present in the inferred rotational velocity distribution with 25% of the sample having 0 <ve <100 km s^-1 and the high velocity component having v_e ∼ 250 km s^-1. There is no evidence from the spatial and radial velocity distributions of the two components that they represent either field and cluster populations or different episodes of star formation. Be-type stars have also been identified. Conclusions. The bi-modal rotational velocity distribution in our sample resembles that found for late-B and early-A type stars.While magnetic braking appears to be a possible mechanism for producing the low-velocity component, we can not rule out alternative explanations. © ESO 2013.

A simple approach to the supernova progenitor–explosion connection

We present a new approach to understand the landscape of supernova explosion energies, ejected nickel masses, and neutron star birth masses. In contrast to other recent parametric approaches, our model predicts the properties of neutrino-driven explosions based on the pre-collapse stellar structure without the need for hydrodynamic simulations. The model is based on physically motivated scaling laws and simple differential equations describing the shock propagation, the contraction of the neutron star, the neutrino emission, the heating conditions, and the explosion energetics. Using model parameters compatible with multi-D simulations and a fine grid of thousands of supernova progenitors, we obtain a variegated landscape of neutron star and black hole formation similar to other parametrized approaches and find good agreement with semi-empirical measures for the ‘explodability’ of massive stars. Our predicted explosion properties largely conform to observed correlations between the nickel mass and explosion energy. Accounting for the coexistence of outflows and downflows during the explosion phase, we naturally obtain a positive correlation between explosion energy and ejecta mass. These correlations are relatively robust against parameter variations, but our results suggest that there is considerable leeway in parametric models to widen or narrow the mass ranges for black hole and neutron star formation and to scale explosion energies up or down. Our model is currently limited to an all-or-nothing treatment of fallback and there remain some minor discrepancies between model predictions and observational constraints.

High-resolution optical spectroscopy of the sharp-lined B-type star HD83206

Very-high-resolution (R~160000) spectroscopic observations are presented for the early B-type star, HD83206. Because it has very sharp metal lines, this star affords an opportunity to test theories of model atmospheres and line formation. Non-LTE model atmosphere calculations have been used to estimate the atmospheric parameters and absolute metal abundances (C, N, O, Mg and Si); an LTE analysis was also undertaken to investigate the validity of this simpler approach and to estimate an iron abundance. For the non-LTE calculations, there is excellent agreement with observations of the Balmer lines Ha and Hd and the lines of Siii and Siiii for atmospheric parameters of Teff~=21700+/-600K and logg~=4.00+/-0.15dex. The agreement is less convincing for the LTE calculations, and a higher gravity is deduced. Careful comparison of the metal line profiles with non-LTE calculations implies that the projected rotational and microturbulent velocities have maximum values of ~=5 and ~=2kms-1, respectively. The latter value is smaller than has often been adopted in LTE model atmosphere analyses of main-sequence stars. Non-LTE absolute metal abundances are estimated, and a comparison with those for normal B-type stars (deduced using similar non-LTE techniques) shows no significant differences. A comparison of the abundances deduced using non-LTE and LTE calculations implies systematic differences of 0.1-0.2dex, showing the importance of using a non-LTE approach when accurate absolute abundances are required. Its location in the Hertzsprung-Russell diagram and normal metal abundance lead us to conclude that HD83206 is probably a main-sequence B-type star. As such, it is among the sharpest-lined young B-type star discovered to date.

Line-profile tomography of exoplanet transits - II. A gas-giant planet transiting a rapidly rotating A5 star

Most of our knowledge of extrasolar planets rests on precise radial-velocity measurements, either for direct detection or for confirmation of the planetary origin of photometric transit signals. This has limited our exploration of the parameter space of exoplanet hosts to solar- and later-type, sharp-lined stars. Here we extend the realm of stars with known planetary companions to include hot, fast-rotating stars. Planet-like transits have previously been reported in the light curve obtained by the SuperWASP survey of the A5 star HD15082 (WASP-33 V = 8.3, v sini = 86 km s-1). Here we report further photometry and time-series spectroscopy through three separate transits, which we use to confirm the existence of a gas-giant planet with an orbital period of 1.22d in orbit around HD15082. From the photometry and the properties of the planet signal travelling through the spectral line profiles during the transit, we directly derive the size of the planet, the inclination and obliquity of its orbital plane and its retrograde orbital motion relative to the spin of the star. This kind of analysis opens the way to studying the formation of planets around a whole new class of young, early-type stars, hence under different physical conditions and generally in an earlier stage of formation than in sharp-lined late-type stars. The reflex orbital motion of the star caused by the transiting planet is small, yielding an upper mass limit of 4.1MJupiter on the planet. We also find evidence of a third body of substellar mass in the system, which may explain the unusual orbit of the transiting planet. In HD 15082, the stellar line profiles also show evidence of non-radial pulsations, clearly distinct from the planetary transit signal. This raises the intriguing possibility that tides raised by the close-in planet may excite or amplify the pulsations in such stars.

Synthesis and characterization of polymer networks and star polymers containing a novel, hydrolyzable acetal-based dimethacrylate cross-linker

An acid-labile dimethaerylate acetal cross-linker,di(methacryloyloxy-l-ethoxy)methane(DMOEM), was synthesized by the reaction of 2-hydroxyethyl methacrylate and paraformaldehyde using p-toluenesulfonic acid and toluene as catalyst and solvent, respectively. Group transfer polymerization was employed to use this cross-linker in the preparation of nine hydrolyzable polymer structures: one neat cross-linker network, one randomly cross-linked network of methyl methacrylate (MMA), and seven star-shaped polymers of MMA. Gel permeation chromatography (GPC) in tetrahydrofuran (THF) confirmed the narrow molecular weight distributions of the linear polymer precursors to the stars and demonstrated the increase in molecular weight upon the addition of cross-linker for the formation of star-shaped polymers. Characterization of the star polymers in THF using static light scattering and GPC showed that the molecular weights and the number of arms of each star polymer increased with an increase in the molar ratio of cross-linker to initiator and with a decrease in the molar ratio of monomer to initiator. The star polymers with DMOEM cores bore a smaller number of arms than those cross-linked with the non-hydrolyzable commercial cross-linker ethylene glycol dimethacrylate due to the bulkier structure of DMOEM. All DMOEM-containing polymer networks and star polymers were completely hydrolyzed within 48 h using hydrochloric acid in THF.

EL CVn-type binaries - Discovery of 17 helium white dwarf precursors in bright eclipsing binary star systems

The star 1SWASP J024743.37-251549.2 was recently discovered to be a binary star in which an A-type dwarf star eclipses the remnant of a disrupted red giant star (WASP 0247-25 B). The remnant is in a rarely observed state evolving to higher effective temperatures at nearly constant luminosity prior to becoming a very low mass white dwarf composed almost entirely of helium, i.e. it is a pre-helium white dwarf (pre-He-WD). We have used the photometric database from theWide Angle Search for Planets (WASP) to find 17 eclipsing binary stars with orbital periods P = 0.7-2.2 d with similar light curves to 1SWASP J024743.37-251549.2. The only star in this group previously identified as a variable star is the brightest one, EL CVn, which we adopt as the prototype for this class of eclipsing binary star. The characteristic light curves of EL CVn-type stars show a total eclipse by an A-type dwarf star of a smaller, hotter star and a secondary eclipse of comparable depth to the primary eclipse. We have used new spectroscopic observations for six of these systems to confirm that the companions to the A-type stars in these binaries have very low masses (≈0.2M⊙). This includes the companion to EL CVn which was not previously known to be a pre-He-WD. EL CVn-type binary star systems will enable us to study the formation of very low mass white dwarfs in great detail, particularly in those cases where the pre-He-WD star shows non-radial pulsations similar to those recently discovered in WASP0247-25 B. © 2013 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society.