The generalized Robinson-Foulds metric

The Robinson-Foulds (RF) metric is arguably the most widely used measure of phylogenetic tree similarity, despite its well-known shortcomings: For example, moving a single taxon in a tree can result in a tree that has maximum distance to the original one; but the two trees are identical if we remove the single taxon. To this end, we propose a natural extension of the RF metric that does not simply count identical clades but instead, also takes similar clades into consideration. In contrast to previous approaches, our model requires the matching between clades to respect the structure of the two trees, a property that the classical RF metric exhibits, too. We show that computing this generalized RF metric is, unfortunately, NP-hard. We then present a simple Integer Linear Program for its computation, and evaluate it by an all-against-all comparison of 100 trees from a benchmark data set. We find that matchings that respect the tree structure differ significantly from those that do not, underlining the importance of this natural condition.Comment: Peer-reviewed and presented as part of the 13th Workshop on Algorithms in Bioinformatics (WABI2013

Vlasov-Maxwell, self-consistent electromagnetic wave emission simulations in the solar corona

1.5D Vlasov-Maxwell simulations are employed to model electromagnetic emission generation in a fully self-consistent plasma kinetic model for the first time in the solar physics context. The simulations mimic the plasma emission mechanism and Larmor drift instability in a plasma thread that connects the Sun to Earth with the spatial scales compressed appropriately. The effects of spatial density gradients on the generation of electromagnetic radiation are investigated. It is shown that 1.5D inhomogeneous plasma with a uniform background magnetic field directed transverse to the density gradient is aperiodically unstable to Larmor-drift instability. The latter results in a novel effect of generation of electromagnetic emission at plasma frequency. When density gradient is removed (i.e. when plasma becomes stable to Larmor-drift instability) and a $low$ density, super-thermal, hot beam is injected along the domain, in the direction perpendicular to the magnetic field, plasma emission mechanism generates non-escaping Langmuir type oscillations which in turn generate escaping electromagnetic radiation. It is found that in the spatial location where the beam is injected, the standing waves, oscillating at the plasma frequency, are excited. These can be used to interpret the horizontal strips observed in some dynamical spectra. Quasilinear theory predictions: (i) the electron free streaming and (ii) the beam long relaxation time, in accord with the analytic expressions, are corroborated via direct, fully-kinetic simulation. Finally, the interplay of Larmor-drift instability and plasma emission mechanism is studied by considering $dense$ electron beam in the Larmor-drift unstable (inhomogeneous) plasma. http://www.maths.qmul.ac.uk/~tsiklauri/movie1.mpg * http://www.maths.qmul.ac.uk/~tsiklauri/movie2.mpg * http://www.maths.qmul.ac.uk/~tsiklauri/movie3.mpgComment: Solar Physics (in press, the final, accepted version

Statistical Survey of Type III Radio Bursts at Long Wavelengths Observed by the Solar TErrestrial RElations Observatory (STEREO)/Waves Instruments: Radio Flux Density Variations with Frequency

We have performed a statistical study of $152$ Type III radio bursts observed by Solar TErrestrial RElations Observatory (STEREO)/Waves between May 2007 and February 2013. We have investigated the flux density between $125$kHz and $16$MHz. Both high- and low-frequency cutoffs have been observed in $60\,$ of events suggesting an important role of propagation. As already reported by previous authors, we observed that the maximum flux density occurs at $1$MHz on both spacecraft. We have developed a simplified analytical model of the flux density as a function of radial distance and compared it to the STEREO/Waves data.Comment: published in Solar Physic

A distance for partially labeled trees

In a number of practical situations, data have structure and the relations among its component parts need to be coded with suitable data models. Trees are usually utilized for representing data for which hierarchical relations can be defined. This is the case in a number of fields like image analysis, natural language processing, protein structure, or music retrieval, to name a few. In those cases, procedures for comparing trees are very relevant. An approximate tree edit distance algorithm has been introduced for working with trees labeled only at the leaves. In this paper, it has been applied to handwritten character recognition, providing accuracies comparable to those by the most comprehensive search method, being as efficient as the fastest.This work is supported by the Spanish Ministry projects DRIMS (TIN2009-14247-C02), and Consolider Ingenio 2010 (MIPRCV, CSD2007-00018), partially supported by EU ERDF and the Pascal Network of Excellence

Atmospheric Heating and Wind Acceleration: Results for Cool Evolved Stars based on Proposed Processes

A chromosphere is a universal attribute of stars of spectral type later than ~F5. Evolved (K and M) giants and supergiants (including the zeta Aurigae binaries) show extended and highly turbulent chromospheres, which develop into slow massive winds. The associated continuous mass loss has a significant impact on stellar evolution, and thence on the chemical evolution of galaxies. Yet despite the fundamental importance of those winds in astrophysics, the question of their origin(s) remains unsolved. What sources heat a chromosphere? What is the role of the chromosphere in the formation of stellar winds? This chapter provides a review of the observational requirements and theoretical approaches for modeling chromospheric heating and the acceleration of winds in single cool, evolved stars and in eclipsing binary stars, including physical models that have recently been proposed. It describes the successes that have been achieved so far by invoking acoustic and MHD waves to provide a physical description of plasma heating and wind acceleration, and discusses the challenges that still remain.Comment: 46 pages, 9 figures, 1 table; modified and unedited manuscript; accepted version to appear in: Giants of Eclipse, eds. E. Griffin and T. Ake (Berlin: Springer

New Insights into White-Light Flare Emission from Radiative-Hydrodynamic Modeling of a Chromospheric Condensation

(abridged) The heating mechanism at high densities during M dwarf flares is poorly understood. Spectra of M dwarf flares in the optical and near-ultraviolet wavelength regimes have revealed three continuum components during the impulsive phase: 1) an energetically dominant blackbody component with a color temperature of T $\sim$ 10,000 K in the blue-optical, 2) a smaller amount of Balmer continuum emission in the near-ultraviolet at lambda $<$ 3646 Angstroms and 3) an apparent pseudo-continuum of blended high-order Balmer lines. These properties are not reproduced by models that employ a typical "solar-type" flare heating level in nonthermal electrons, and therefore our understanding of these spectra is limited to a phenomenological interpretation. We present a new 1D radiative-hydrodynamic model of an M dwarf flare from precipitating nonthermal electrons with a large energy flux of $10^{13}$ erg cm$^{-2}$ s$^{-1}$. The simulation produces bright continuum emission from a dense, hot chromospheric condensation. For the first time, the observed color temperature and Balmer jump ratio are produced self-consistently in a radiative-hydrodynamic flare model. We find that a T $\sim$ 10,000 K blackbody-like continuum component and a small Balmer jump ratio result from optically thick Balmer and Paschen recombination radiation, and thus the properties of the flux spectrum are caused by blue light escaping over a larger physical depth range compared to red and near-ultraviolet light. To model the near-ultraviolet pseudo-continuum previously attributed to overlapping Balmer lines, we include the extra Balmer continuum opacity from Landau-Zener transitions that result from merged, high order energy levels of hydrogen in a dense, partially ionized atmosphere. This reveals a new diagnostic of ambient charge density in the densest regions of the atmosphere that are heated during dMe and solar flares.Comment: 50 pages, 2 tables, 13 figures. Accepted for publication in the Solar Physics Topical Issue, "Solar and Stellar Flares". Version 2 (June 22, 2015): updated to include comments by Guest Editor. The final publication is available at Springer via http://dx.doi.org/10.1007/s11207-015-0708-

Measurement of D*+/- meson production in jets from pp collisions at sqrt(s) = 7 TeV with the ATLAS detector

This paper reports a measurement of D*+/- meson production in jets from proton-proton collisions at a center-of-mass energy of sqrt(s) = 7 TeV at the CERN Large Hadron Collider. The measurement is based on a data sample recorded with the ATLAS detector with an integrated luminosity of 0.30 pb^-1 for jets with transverse momentum between 25 and 70 GeV in the pseudorapidity range |eta| < 2.5. D*+/- mesons found in jets are fully reconstructed in the decay chain: D*+ -> D0pi+, D0 -> K-pi+, and its charge conjugate. The production rate is found to be N(D*+/-)/N(jet) = 0.025 +/- 0.001(stat.) +/- 0.004(syst.) for D*+/- mesons that carry a fraction z of the jet momentum in the range 0.3 < z < 1. Monte Carlo predictions fail to describe the data at small values of z, and this is most marked at low jet transverse momentum.Comment: 10 pages plus author list (22 pages total), 5 figures, 1 table, matches published version in Physical Review

Search for supersymmetry in final states with jets, missing transverse momentum and one isolated lepton in sqrt{s} = 7 TeV pp collisions using 1 fb-1 of ATLAS data

We present an update of a search for supersymmetry in final states containing jets, missing transverse momentum, and one isolated electron or muon, using 1.04 fb^-1 of proton-proton collision data at sqrt{s} = 7 TeV recorded by the ATLAS experiment at the LHC in the first half of 2011. The analysis is carried out in four distinct signal regions with either three or four jets and variations on the (missing) transverse momentum cuts, resulting in optimized limits for various supersymmetry models. No excess above the standard model background expectation is observed. Limits are set on the visible cross-section of new physics within the kinematic requirements of the search. The results are interpreted as limits on the parameters of the minimal supergravity framework, limits on cross-sections of simplified models with specific squark and gluino decay modes, and limits on parameters of a model with bilinear R-parity violation.Comment: 18 pages plus author list (30 pages total), 9 figures, 4 tables, final version to appear in Physical Review

Reducing heterotic M-theory to five dimensional supergravity on a manifold with boundary

This paper constructs the reduction of heterotic $M$-theory in eleven dimensions to a supergravity model on a manifold with boundary in five dimensions using a Calabi-Yau three-fold. New results are presented for the boundary terms in the action and for the boundary conditions on the bulk fields. Some general features of dualisation on a manifold with boundary are used to explain the origin of some topological terms in the action. The effect of gaugino condensation on the fermion boundary conditions leads to a `twist' in the chirality of the gravitino which can provide an uplifting mechanism in the vacuum energy to cancel the cosmological constant after moduli stabilisation.Comment: 16 pages, RevTe