Binary-corrected velocity dispersions from single- and multi-epoch radial velocities: massive stars in R136 as a test case

Orbital motions from binary stars can broaden the observed line-of-sight velocity distribution of a stellar system, artificially inflating the measured line-of-sight velocity dispersion, which can in turn lead to erroneous conclusions about the dynamical state of the system. Cottaar et al. (2012b) proposed a maximum likelihood procedure to recover the intrinsic velocity dispersion of a resolved star cluster from a single epoch of radial velocity data of individual stars, which they achieved by simultaneously fitting the intrinsic velocity distribution of the single stars and the centres of mass of the binaries along with the velocity shifts caused by binary orbital motions. Assuming well-characterized binary properties, they showed that this procedure can accurately reproduce intrinsic velocity dispersions below 1 km s$^{-1}$ for solar-type stars. Here we investigate the systematic offsets induced in cases where the binary properties are uncertain, and we show how two epochs of radial velocity data with an appropriate baseline can help to mitigate these systematic effects. We first test the method above using Monte Carlo simulations, taking into account the large uncertainties in the binary properties of OB stars. We then apply it to radial velocity data in the young massive cluster R136, an example for which the intrinsic velocity dispersion of O-type stars is known from an intensive multi-epoch approach. For typical velocity dispersions of young massive clusters ($\gtrsim 4$ km s$^{-1}$) and with a single epoch of data, we demonstrate that the method can just about distinguish between a cluster in virial equilibrium and an unbound cluster. This is due to the higher spectroscopic binary fraction and more loosely constrained distributions of orbital parameters of OB stars compared to solar-type stars. By extending the maximum likelihood method to multi-epoch data, Comment: Accepted by A&A; minor corrections made on November 2

Can low metallicity binaries avoid merging?

Rapid mass transfer in a binary system can drive the accreting star out of thermal equilibrium, causing it to expand. This can lead to a contact system, strong mass loss from the system and possibly merging of the two stars. In low metallicity stars the timescale for heat transport is shorter due to the lower opacity. The accreting star can therefore restore thermal equilibrium more quickly and possibly avoid contact. We investigate the effect of accretion onto main sequence stars with radiative envelopes with different metallicities. We find that a low metallicity (Z<0.001), 4 solar mass star can endure a 10 to 30 times higher accretion rate before it reaches a certain radius than a star at solar metallicity. This could imply that up to two times fewer systems come into contact during rapid mass transfer when we compare low metallicity. This factor is uncertain due to the unknown distribution of binary parameters and the dependence of the mass transfer timescale on metallicity. In a forthcoming paper we will present analytic fits to models of accreting stars at various metallicities intended for the use in population synthesis models.Comment: To appear in the proceedings of "First Stars III", Santa Fe, New Mexico, July 16-20, 2007, 3 pages, 2 figure

Large-scale mantle discontinuity topography beneath Europe: Signature of akimotoite in subducting slabs

The mantle transition zone is delineated by seismic discontinuities around 410 and 660 km, which are generally related to mineral phase transitions. Study of the topography of the discontinuities further constrains which phase transitions play a role and, combined with their Clapeyron slopes, what temperature variations occur. Here we use P to S converted seismic waves or receiver functions to study the topography of the mantle seismic discontinuities beneath Europe and the effect of subducting and ponding slabs beneath southern Europe on these features. We combine roughly 28,000 of the highest quality receiver functions into a common conversion point stack. In the topography of the discontinuity around 660 km, we find broadscale depressions of 30 km beneath central Europe and around the Mediterranean. These depressions do not correlate with any topography on the discontinuity around 410 km. Explaining these strong depressions by purely thermal effects on the dissociation of ringwoodite to bridgmanite and periclase requires unrealistically large temperature reductions. Presence of several wt % water in ringwoodite leads to a deeper phase transition, but complementary observations, such as elevated Vp/Vs ratio, attenuation, and electrical conductivity, are not observed beneath central Europe. Our preferred hypothesis is the dissociation of ringwoodite into akimotoite and periclase in cold downwelling slabs at the bottom of the transition zone. The strongly negative Clapeyron slope predicted for the subsequent transition of akimotoite to bridgmanite explains the depression with a temperature reduction of 200–300 K and provides a mechanism to pond slabs in the first place

Het beeld van Nederland: hoe zien Molukkers, Chinezen, woonwagenbewoners en Turken de Nederlanders en zichzelf?

Contains fulltext : 3768.pdf (publisher's version ) (Open Access

Large-scale mantle discontinuity topography beneath Europe: Signature of akimotoite in subducting slabs

The mantle transition zone is delineated by seismic discontinuities around 410 and 660 km, which are generally related to mineral phase transitions. Study of the topography of the discontinuities further constrains which phase transitions play a role and, combined with their Clapeyron slopes, what temperature variations occur. Here we use P to S converted seismic waves or receiver functions to study the topography of the mantle seismic discontinuities beneath Europe and the effect of subducting and ponding slabs beneath southern Europe on these features. We combine roughly 28,000 of the highest quality receiver functions into a common conversion point stack. In the topography of the discontinuity around 660 km, we find broadscale depressions of 30 km beneath central Europe and around the Mediterranean. These depressions do not correlate with any topography on the discontinuity around 410 km. Explaining these strong depressions by purely thermal effects on the dissociation of ringwoodite to bridgmanite and periclase requires unrealistically large temperature reductions. Presence of several wt % water in ringwoodite leads to a deeper phase transition, but complementary observations, such as elevated Vp/Vs ratio, attenuation, and electrical conductivity, are not observed beneath central Europe. Our preferred hypothesis is the dissociation of ringwoodite into akimotoite and periclase in cold downwelling slabs at the bottom of the transition zone. The strongly negative Clapeyron slope predicted for the subsequent transition of akimotoite to bridgmanite explains the depression with a temperature reduction of 200–300 K and provides a mechanism to pond slabs in the first place.SC is funded by the Drapers’ Company Research Fellowship through Pembroke College, Cambridge, UK. AD was funded by the European Research Council under the European Community’s Seventh Framework Programme (FP7/20072013/ERC grant agreement 204995) and by a Philip Leverhulme Prize.This is the final version of the article. It first appeared from Wiley via http://dx.doi.org/10.1002/2015JB012452 The data used are freely available from the IRIS (www.iris.edu) and ORFEUS (http://www.orfeus-eu.org) databases

Morphology of seismically slow lower-mantle structures

Large low shear velocity provinces (LLSVPs), whose origin and dynamic implication remain enigmatic, dominate the lowermost mantle. For decades, seismologists have created increasingly detailed pictures of the LLSVPs through tomographic models constructed with different modeling methodologies, data sets, parametrizations and regularizations. Here, we extend the cluster analysis methodology of Lekic $\textit{et al.}$, to classify seismic mantle structure in five recent global shear wave speed ($\textit{V}_S$) tomographic models into three groups. By restricting the analysis to moving depth windows of the radial profiles of $\textit{V}_S$, we assess the vertical extent of features. We also show that three clusters are better than two (or four) when representing the entire lower mantle, as the boundaries of the three clusters more closely follow regions of high lateral $\textit{V}_S$ gradients. Qualitatively, we relate the anomalously slow cluster to the LLSVPs, the anomalously fast cluster to slab material entering the lower mantle and the neutral cluster to ‘background’ lower mantle material. We obtain compatible results by repeating the analysis on recent global $\textit{P}$-wave speed ($\textit{V}_P$) models, although we find less agreement across $\textit{V}_P$ models. We systematically show that the clustering results, even in detail, agree remarkably well with a wide range of local waveform studies. This suggests that the two LLSVPs consist of multiple internal anomalies with a wide variety of morphologies, including shallowly to steeply sloping, and even overhanging, boundaries. Additionally, there are indications of previously unrecognized meso-scale features, which, like the Perm anomaly, are separated from the two main LLSVPs beneath the Pacific and Africa. The observed wide variety of structure size and morphology offers a challenge to recreate in geodynamic models; potentially, the variety can result from various degrees of mixing of several compositionally distinct components. Finally, we obtain new, much larger estimates of the volume/mass occupied by LLSVPs— 8.0 per cent ±0.9 ($\mu$ ± 1$\sigma$) of whole mantle volume and 9.1 per cent ±1.0 ($\mu$ ± 1$\sigma$) of whole mantle mass—and discuss implications for associating the LLSVPs with the hidden reservoir enriched in heat producing elements.National Science Foundation (EAR1352214), Packard Foundation, Pembroke College, Cambridge (Drapers’ Company Research Fellowship

Converted phases from sharp 1000 km depth mid-mantle heterogeneity beneath Western Europe

Until recently, most of the lower mantle was generally considered to be well-mixed with strong heterogeneity restricted to the lowermost several hundred kilometres above the core–mantle boundary, known as the $\textbf{D}$″ layer. However several recent studies have started to hint at a potential change in Earth's structure at mid-mantle depths beneath the transition zone. Here we present a continental-wide search of Europe and the North Atlantic for mid-mantle P-to-s wave converted phases. Our data set consists of close to 50,000 high quality receiver functions. These are combined in slowness and depth stacks to identify seismic discontinuities in the range of 800–1400 km depth to determine at which depths and in which tectonic settings these features exist. Receiver functions are computed in different frequency bands to resolve the sharpness of the observed discontinuities. We find most seismic velocity jumps are observed between 975–1050 km depth, localised beneath western Europe and Iceland. The shear wave velocity jumps are roughly 1–2.5% velocity increase with depth occurring over less than 8 km in width. The most robust observations are coincident with areas of active upwelling (under Iceland) and an elongate lateral low velocity anomaly imaged in recent tomographic models which has been interpreted as diverted plume material at depth. The lack of any suggested phase change in a normal pyrolitic mantle composition at around 1000 km depth indicates the presence of regional chemical heterogeneity within the mid-mantle, potentially caused by diverted plume material. We hypothesise that our observations represent either a phase change within chemically distinct plume material itself, or are caused by small scale chemical heterogeneities entrained within the upwelling plume, either in the form of recycled basaltic material or deep sourced chemically distinct material from LLSVPs. Our observations, which cannot be directly linked to an area of either active or ancient subduction, along with observations in other hotspot regions, suggest that such mid-mantle seismic features are not unique to subduction zones despite the large number of observations that have previously been made in such settings.The facilities of IRIS Data Services, and specifically the IRIS Data Management Center, as well as the ORFEUS data centre were used for access to waveforms, related metadata, and/or derived products used in this study. IRIS Data Services are funded through the Seismological Facilities for the Advancement of Geoscience and EarthScope (SAGE) Proposal of the National Science Foundation under Cooperative Agreement EAR-1261681. For full citation list of all FDSN networks please see the Supplementary Material. Seismometers for the Cambridge seismic network in Iceland were borrowed from the Natural Environment Research Council (NERC) SEIS-UK (loans 857, 968 and 1022), and funded by research grants from the NERC and the European Community's Seventh Framework Programme Grant No. 308377 (Project FUTUREVOLC), to Robert S. White. J.J. was funded by a graduate studentship from NERC (LBAG/148 Task 5). S.C. is funded by the Drapers' Company Research Fellowship through Pembroke College, Cambridge. Thanks are also extended to the Icelandic Meteorological office for sharing data that were used in this study. A.D. and J.J. were funded by the European Research Council under the European Community's Seventh Framework Programme (FP7/2007–2013/ERC grant agreement 204995) and by a Philip Leverhulme Prize. Data was downloaded from IRIS DMC and figures made using GMT (Wessel and Smith, 2001)

Characterizing a cluster's dynamic state using a single epoch of radial velocities

Radial velocity measurements can be used to constrain the dynamical state of a stellar cluster. However, for clusters with velocity dispersions smaller than a few km/s the observed radial velocity distribution tends to be dominated by the orbital motions of binaries rather than the stellar motions through the potential well of the cluster. Our goal is to characterize the intrinsic velocity distribution of a cluster from a single epoch of radial velocity data, even for a cluster with a velocity dispersion of a fraction of a km/s, using a maximum likelihood procedure. Assuming a period, mass ratio, and eccentricity distribution for the binaries in the observed cluster this procedure fits a dynamical model describing the velocity distribution for the single stars and center of masses of the binaries, simultaneously with the radial velocities caused by binary orbital motions, using all the information available in the observed velocity distribution. We find that the fits to the intrinsic velocity distribution depend only weakly on the binary properties assumed, so the uncertainty in the fitted parameters tends to be dominated by statistical uncertainties. Based on Monte Carlo simulations we provide an estimate of how these statistical uncertainties vary with the velocity dispersion, binary fraction, and the number of observed stars, which can be used to estimate the sample size needed to reach a specific accuracy. Finally we test the method on the well-studied open cluster NGC 188, showing that it can reproduce a velocity dispersion of only 0.5 km/s using a single epoch of the multi-epoch radial velocity data. If the binary period, mass ratio, and eccentricity distribution of the observed stars are roughly known, this procedure can be used to correct for the effect of binary orbital motions on an observed velocity distribution. [Abridged]Comment: 11 pages, 6 figures, accepted by A&