Solving non-uniqueness in agglomerative hierarchical clustering using multidendrograms

In agglomerative hierarchical clustering, pair-group methods suffer from a problem of non-uniqueness when two or more distances between different clusters coincide during the amalgamation process. The traditional approach for solving this drawback has been to take any arbitrary criterion in order to break ties between distances, which results in different hierarchical classifications depending on the criterion followed. In this article we propose a variable-group algorithm that consists in grouping more than two clusters at the same time when ties occur. We give a tree representation for the results of the algorithm, which we call a multidendrogram, as well as a generalization of the Lance and Williams' formula which enables the implementation of the algorithm in a recursive way.Comment: Free Software for Agglomerative Hierarchical Clustering using Multidendrograms available at http://deim.urv.cat/~sgomez/multidendrograms.ph

Fermat Principle in Finsler Spacetimes

It is shown that, on a manifold with a Finsler metric of Lorentzian signature, the lightlike geodesics satisfy the following variational principle. Among all lightlike curves from a point (emission event) to a timelike curve (worldline of receiver), the lightlike geodesics make the arrival time stationary. Here ``arrival time'' refers to a parametrization of the timelike curve. This variational principle can be applied (i) to the vacuum light rays in an alternative spacetime theory, based on Finsler geometry, and (ii) to light rays in an anisotropic non-dispersive medium with a general-relativistic spacetime as background.Comment: 18 pages, submitted to Gen. Rel. Gra

Formal Verification of Hardware Synthesis

Original manuscript: January 21, 2013We report on the implementation of a certified compiler for a high-level hardware description language (HDL) called Fe-Si (FEatherweight SynthesIs). Fe-Si is a simplified version of Bluespec, an HDL based on a notion of guarded atomic actions. Fe-Si is defined as a dependently typed deep embedding in Coq. The target language of the compiler corresponds to a synthesisable subset of Verilog or VHDL. A key aspect of our approach is that input programs to the compiler can be defined and proved correct inside Coq. Then, we use extraction and a Verilog back-end (written in OCaml) to get a certified version of a hardware design.United States. Defense Advanced Research Projects Agency (Agreement FA8750-12-2-0110

Entanglement between Demand and Supply in Markets with Bandwagon Goods

Whenever customers' choices (e.g. to buy or not a given good) depend on others choices (cases coined 'positive externalities' or 'bandwagon effect' in the economic literature), the demand may be multiply valued: for a same posted price, there is either a small number of buyers, or a large one -- in which case one says that the customers coordinate. This leads to a dilemma for the seller: should he sell at a high price, targeting a small number of buyers, or at low price targeting a large number of buyers? In this paper we show that the interaction between demand and supply is even more complex than expected, leading to what we call the curse of coordination: the pricing strategy for the seller which aimed at maximizing his profit corresponds to posting a price which, not only assumes that the customers will coordinate, but also lies very near the critical price value at which such high demand no more exists. This is obtained by the detailed mathematical analysis of a particular model formally related to the Random Field Ising Model and to a model introduced in social sciences by T C Schelling in the 70's.Comment: Updated version, accepted for publication, Journal of Statistical Physics, online Dec 201

Synergistic warm inflation

We consider an alternative warm inflationary scenario in which $n$ scalar fields coupled to a dissipative matter fluid cooperate to produce power--law inflation. The scalar fields are driven by an exponential potential and the bulk dissipative pressure coefficient is linear in the expansion rate. We find that the entropy of the fluid attains its asymptotic value in a characteristic time proportional to the square of the number of fields. This scenario remains nearly isothermal along the inflationary stage. The perturbations in energy density and entropy are studied in the long--wavelength regime and seen to grow roughly as the square of the scale factor. They are shown to be compatible with COBE measurements of the fluctuations in temperature of the CMB.Comment: 13 pages, Revtex 3 To be published in Physical Review

Distributed phase-covariant cloning with atomic ensembles via quantum Zeno dynamics

We propose an interesting scheme for distributed orbital state quantum cloning with atomic ensembles based on the quantum Zeno dynamics. These atomic ensembles which consist of identical three-level atoms are trapped in distant cavities connected by a single-mode integrated optical star coupler. These qubits can be manipulated through appropriate modulation of the coupling constants between atomic ensemble and classical field, and the cavity decay can be largely suppressed as the number of atoms in the ensemble qubits increases. The fidelity of each cloned qubit can be obtained with analytic result. The present scheme provides a new way to construct the quantum communication network.Comment: 5 pages, 4 figure

Uncovering precision phenotype-biomarker associations in traumatic brain injury using topological data analysis

Background: Traumatic brain injury (TBI) is a complex disorder that is traditionally stratified based on clinical signs and symptoms. Recent imaging and molecular biomarker innovations provide unprecedented opportunities for improved TBI precision medicine, incorporating patho-anatomical and molecular mechanisms. Complete integration of these diverse data for TBI diagnosis and patient stratification remains an unmet challenge. Methods and findings: The Transforming Research and Clinical Knowledge in Traumatic Brain Injury (TRACK-TBI) Pilot multicenter study enrolled 586 acute TBI patients and collected diverse common data elements (TBI-CDEs) across the study population, including imaging, genetics, and clinical outcomes. We then applied topology-based data-driven discovery to identify natural subgroups of patients, based on the TBI-CDEs collected. Our hypothesis was two-fold: 1) A machine learning tool known as topological data analysis (TDA) would reveal data-driven patterns in patient outcomes to identify candidate biomarkers of recovery, and 2) TDA-identified biomarkers would significantly predict patient outcome recovery after TBI using more traditional methods of univariate statistical tests. TDA algorithms organized and mapped the data of TBI patients in multidimensional space, identifying a subset of mild TBI patients with a specific multivariate phenotype associated with unfavorable outcome at 3 and 6 months after injury. Further analyses revealed that this patient subset had high rates of post-traumatic stress disorder (PTSD), and enrichment in several distinct genetic polymorphisms associated with cellular responses to stress and DNA damage (PARP1), and in striatal dopamine processing (ANKK1, COMT, DRD2). Conclusions: TDA identified a unique diagnostic subgroup of patients with unfavorable outcome after mild TBI that were significantly predicted by the presence of specific genetic polymorphisms. Machine learning methods such as TDA may provide a robust method for patient stratification and treatment planning targeting identified biomarkers in future clinical trials in TBI patients

Fitting the integrated Spectral Energy Distributions of Galaxies

Fitting the spectral energy distributions (SEDs) of galaxies is an almost universally used technique that has matured significantly in the last decade. Model predictions and fitting procedures have improved significantly over this time, attempting to keep up with the vastly increased volume and quality of available data. We review here the field of SED fitting, describing the modelling of ultraviolet to infrared galaxy SEDs, the creation of multiwavelength data sets, and the methods used to fit model SEDs to observed galaxy data sets. We touch upon the achievements and challenges in the major ingredients of SED fitting, with a special emphasis on describing the interplay between the quality of the available data, the quality of the available models, and the best fitting technique to use in order to obtain a realistic measurement as well as realistic uncertainties. We conclude that SED fitting can be used effectively to derive a range of physical properties of galaxies, such as redshift, stellar masses, star formation rates, dust masses, and metallicities, with care taken not to over-interpret the available data. Yet there still exist many issues such as estimating the age of the oldest stars in a galaxy, finer details ofdust properties and dust-star geometry, and the influences of poorly understood, luminous stellar types and phases. The challenge for the coming years will be to improve both the models and the observational data sets to resolve these uncertainties. The present review will be made available on an interactive, moderated web page (sedfitting.org), where the community can access and change the text. The intention is to expand the text and keep it up to date over the coming years.Comment: 54 pages, 26 figures, Accepted for publication in Astrophysics & Space Scienc

Longitudinal double-spin asymmetry and cross section for inclusive neutral pion production at midrapidity in polarized proton collisions at sqrt(s) = 200 GeV

We report a measurement of the longitudinal double-spin asymmetry A_LL and the differential cross section for inclusive Pi0 production at midrapidity in polarized proton collisions at sqrt(s) = 200 GeV. The cross section was measured over a transverse momentum range of 1 < p_T < 17 GeV/c and found to be in good agreement with a next-to-leading order perturbative QCD calculation. The longitudinal double-spin asymmetry was measured in the range of 3.7 < p_T < 11 GeV/c and excludes a maximal positive gluon polarization in the proton. The mean transverse momentum fraction of Pi0's in their parent jets was found to be around 0.7 for electromagnetically triggered events.Comment: 6 pages, 3 figures, submitted to Phys. Rev. D (RC

High pTp_{T} non-photonic electron production in pp+pp collisions at s\sqrt{s} = 200 GeV

We present the measurement of non-photonic electron production at high transverse momentum ($p_T >$ 2.5 GeV/$c$) in $p$ + $p$ collisions at $\sqrt{s}$ = 200 GeV using data recorded during 2005 and 2008 by the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). The measured cross-sections from the two runs are consistent with each other despite a large difference in photonic background levels due to different detector configurations. We compare the measured non-photonic electron cross-sections with previously published RHIC data and pQCD calculations. Using the relative contributions of B and D mesons to non-photonic electrons, we determine the integrated cross sections of electrons ($\frac{e^++e^-}{2}$) at 3 GeV/$c < p_T <~$10 GeV/$c$ from bottom and charm meson decays to be ${d\sigma_{(B\to e)+(B\to D \to e)} \over dy_e}|_{y_e=0}$ = 4.0$\pm0.5$({\rm stat.})$\pm1.1$({\rm syst.}) nb and ${d\sigma_{D\to e} \over dy_e}|_{y_e=0}$ = 6.2$\pm0.7$({\rm stat.})$\pm1.5$({\rm syst.}) nb, respectively.Comment: 17 pages, 17 figure