The mean, variance and limiting distribution of two statistics sensitive to phylogenetic tree balance

For two decades, the Colless index has been the most frequently used statistic for assessing the balance of phylogenetic trees. In this article, this statistic is studied under the Yule and uniform model of phylogenetic trees. The main tool of analysis is a coupling argument with another well-known index called the Sackin statistic. Asymptotics for the mean, variance and covariance of these two statistics are obtained, as well as their limiting joint distribution for large phylogenies. Under the Yule model, the limiting distribution arises as a solution of a functional fixed point equation. Under the uniform model, the limiting distribution is the Airy distribution. The cornerstone of this study is the fact that the probabilistic models for phylogenetic trees are strongly related to the random permutation and the Catalan models for binary search trees.Comment: Published at http://dx.doi.org/10.1214/105051606000000547 in the Annals of Applied Probability (http://www.imstat.org/aap/) by the Institute of Mathematical Statistics (http://www.imstat.org

Measurement of thermal conductance of silicon nanowires at low temperature

We have performed thermal conductance measurements on individual single crystalline silicon suspended nanowires. The nanowires (130 nm thick and 200 nm wide) are fabricated by e-beam lithography and suspended between two separated pads on Silicon On Insulator (SOI) substrate. We measure the thermal conductance of the phonon wave guide by the 3 method. The cross-section of the nanowire approaches the dominant phonon wavelength in silicon which is of the order of 100 nm at 1K. Above 1.3K the conductance behaves as T3, but a deviation is measured at the lowest temperature which can be attributed to the reduced geometry

The various manifestations of collisionless dissipation in wave propagation

The propagation of an electrostatic wave packet inside a collisionless and initially Maxwellian plasma is always dissipative because of the irreversible acceleration of the electrons by the wave. Then, in the linear regime, the wave packet is Landau damped, so that in the reference frame moving at the group velocity, the wave amplitude decays exponentially with time. In the nonlinear regime, once phase mixing has occurred and when the electron motion is nearly adiabatic, the damping rate is strongly reduced compared to the Landau one, so that the wave amplitude remains nearly constant along the characteristics. Yet, we show here that the electrons are still globally accelerated by the wave packet, and, in one dimension, this leads to a non local amplitude dependence of the group velocity. As a result, a freely propagating wave packet would shrink, and, therefore, so would its total energy. In more than one dimension, not only does the magnitude of the group velocity nonlinearly vary, but also its direction. In the weakly nonlinear regime, when the collisionless damping rate is still significant compared to its linear value, this leads to an effective defocussing effect which we quantify, and which we compare to the self-focussing induced by wave front bowing.Comment: 23 pages, 6 figure

The Initial-Boundary Value Problem in General Relativity

In this article we summarize what is known about the initial-boundary value problem for general relativity and discuss present problems related to it.Comment: 11 pages, 2 figures. Contribution to a special volume for Mario Castagnino's seventy fifth birthda

Merging several separate data acquisition systems at GANIL

International audienceThe Grand Accélérateur National d'Ions (GANIL) as a National User's facility must answer to the problems of synchronization for the coupling of resident and external acquisition systems. For this a solution made up of a C-size VXI electronics called CENTRUM and a software named MERGER was adopted. The technique used for synchronization is based on the distribution of a 32-bit event number or a 48-bit timestamp. The tagged events sent by the various acquisition systems on the network are assembled by MERGER in a single event and sent for the analysis and storag

Predictions of Native American Population Structure Using Linguistic Covariates in a Hidden Regression Framework

International audienceBACKGROUND: The mainland of the Americas is home to a remarkable diversity of languages, and the relationships between genes and languages have attracted considerable attention in the past. Here we investigate to which extent geography and languages can predict the genetic structure of Native American populations. METHODOLOGY/PRINCIPAL FINDINGS: Our approach is based on a Bayesian latent cluster regression model in which cluster membership is explained by geographic and linguistic covariates. After correcting for geographic effects, we find that the inclusion of linguistic information improves the prediction of individual membership to genetic clusters. We further compare the predictive power of Greenberg's and The Ethnologue classifications of Amerindian languages. We report that The Ethnologue classification provides a better genetic proxy than Greenberg's classification at the stock and at the group levels. Although high predictive values can be achieved from The Ethnologue classification, we nevertheless emphasize that Choco, Chibchan and Tupi linguistic families do not exhibit a univocal correspondence with genetic clusters. CONCLUSIONS/SIGNIFICANCE: The Bayesian latent class regression model described here is efficient at predicting population genetic structure using geographic and linguistic information in Native American populations

A 2.5D time-frequency domain model for railway induced soil-building vibration due to railway defects

A new hybrid time-frequency modelling methodology is proposed to simulate the generation of railway vibration caused by singular defects (e.g. joints, switches, crossings), and its propagation through the track, soil and into nearby buildings. To create the full source-to-received model, first the force density due to wheel-rail-defect interaction is calculated using a time domain finite element vehicle-track-soil model. Next, the frequency domain track-soil transfer function is calculated using a 2.5D boundary/finite element approach and coupled with the force densities to recover the free-field response. Finally, the soil-structure interaction of buildings close to the line is computed using a time domain approach. The effect of defect type, train speed and building type (4-storey office block and 8-storey apartment building) on a variety of commonly used international vibration metrics (one-third octaves, PPV, MTVV) is then investigated. It is found that train speed doesn't correlate with building vibration and different defect types have a complex relationship with vibration levels both in the ground and buildings. The 8-storey apartment building has a frequency response dominated by a narrow frequency range, whereas the modal contribution of the 4-storey office building is over a wider frequency band. This results in the 8-storey building having a higher response

Multi-species pair annihilation reactions

We consider diffusion-limited reactions A_i + A_j -> 0 (1 <= i < j <= q) in d space dimensions. For q > 2 and d >= 2 we argue that the asymptotic density decay for such mutual annihilation processes with equal rates and initial densities is the same as for single-species pair annihilation A + A -> 0. In d = 1, however, particle segregation occurs for all q < oo. The total density decays according to a $q$ dependent power law, rho(t) ~ t^{-\alpha(q)}. Within a simplified version of the model \alpha(q) = (q-1) / 2q can be determined exactly. Our findings are supported through Monte Carlo simulations.Comment: 4 pages, revtex; two figures include

The imprints of the Galactic Bar on the Thick Disk with RAVE

We study the kinematics of a local sample of stars, located within a cylinder of 500 pc radius centered on the Sun, in the RAVE data set. We find clear asymmetries in the v R v∞ velocity distributions of thin and thick disk stars: there are more stars moving radially outward for low azimuthal velocities and more radially inward for high azimuthal velocities. Such asymmetries have been previously reported for the thin disk as being due to the Galactic bar, but this is the first time that the same type of structures are seen in the thick disk. Our findings imply that the velocities of thick-disk stars should no longer be described by Schwarzschilds, multivariate Gaussian or purely axisymmetric distributions. Furthermore, the nature of previously reported substructures in the thick disk needs to be revisited as these could be associated with dynamical resonances rather than to accretion events. It is clear that dynamical models of the Galaxy must fit the 3D velocity distributions of the disks, rather than the projected 1D, if we are to understand the Galaxy fully