Beyond inverse Ising model: structure of the analytical solution for a class of inverse problems

I consider the problem of deriving couplings of a statistical model from measured correlations, a task which generalizes the well-known inverse Ising problem. After reminding that such problem can be mapped on the one of expressing the entropy of a system as a function of its corresponding observables, I show the conditions under which this can be done without resorting to iterative algorithms. I find that inverse problems are local (the inverse Fisher information is sparse) whenever the corresponding models have a factorized form, and the entropy can be split in a sum of small cluster contributions. I illustrate these ideas through two examples (the Ising model on a tree and the one-dimensional periodic chain with arbitrary order interaction) and support the results with numerical simulations. The extension of these methods to more general scenarios is finally discussed.Comment: 15 pages, 6 figure

First astronomical unit scale image of the GW Ori triple. Direct detection of a new stellar companion

Young and close multiple systems are unique laboratories to probe the initial dynamical interactions between forming stellar systems and their dust and gas environment. Their study is a key building block to understanding the high frequency of main-sequence multiple systems. However, the number of detected spectroscopic young multiple systems that allow dynamical studies is limited. GW Orionis is one such system. It is one of the brightest young T Tauri stars and is surrounded by a massive disk. Our goal is to probe the GW Orionis multiplicity at angular scales at which we can spatially resolve the orbit. We used the IOTA/IONIC3 interferometer to probe the environment of GW Orionis with an astronomical unit resolution in 2003, 2004, and 2005. By measuring squared visibilities and closure phases with a good UV coverage we carry out the first image reconstruction of GW Ori from infrared long-baseline interferometry. We obtain the first infrared image of a T Tauri multiple system with astronomical unit resolution. We show that GW Orionis is a triple system, resolve for the first time the previously known inner pair (separation $\rho\sim$1.4 AU) and reveal a new more distant component (GW Ori C) with a projected separation of $\sim$8 AU with direct evidence of motion. Furthermore, the nearly equal (2:1) H-band flux ratio of the inner components suggests that either GW Ori B is undergoing a preferential accretion event that increases its disk luminosity or that the estimate of the masses has to be revisited in favour of a more equal mass-ratio system that is seen at lower inclination. Accretion disk models of GW Ori will need to be completely reconsidered because of this outer companion C and the unexpected brightness of companion B.Comment: 5 pages, 9 figures, accepted Astronomy and Astrophysics Letters. 201

Inference algorithms for gene networks: a statistical mechanics analysis

The inference of gene regulatory networks from high throughput gene expression data is one of the major challenges in systems biology. This paper aims at analysing and comparing two different algorithmic approaches. The first approach uses pairwise correlations between regulated and regulating genes; the second one uses message-passing techniques for inferring activating and inhibiting regulatory interactions. The performance of these two algorithms can be analysed theoretically on well-defined test sets, using tools from the statistical physics of disordered systems like the replica method. We find that the second algorithm outperforms the first one since it takes into account collective effects of multiple regulators

Thermal diffusivity measurements of metastable austenite during continuous cooling

The thermal diffusivity of the metastable undercooled austenite is relevant for the quantitative analysis of the carbon and low-alloy steel quench. The standard laser-flash method requires prior thermal equilibrium between the sample and the furnace, which may not be possible to achieve without allowing the metastable phase to transform. Nevertheless, depending upon the steel's hardenability, the thermal transient due to a laser pulse may be much shorter than a cooling transient sufficiently steep to prevent the transformation of the austenite. In one such case, flash measurements were performed during continuous sample cooling and the thermal diffusivity of the metastable austenite was determined by using an extension of the standard analytical model. The adopted analytical model and data reduction procedure are described and the limitations and uncertainties of this method are discussed, also with the aid of a non-linear numerical simulation. The measured thermal diffusivity of the under cooled low-alloy austenite decreases linearly from 5.4•10−6 m2 s−1 at 1133 K to 4.3•10−6 m2 s−1 at 755 K; this trend is in broad agreement with one previous set of measurements upon a low-alloy undercooled austenite and with a large number of previous standard measurements upon stable (high-alloy) austenitic stainless steels

Reverse Engineering Gene Networks with ANN: Variability in Network Inference Algorithms

Motivation :Reconstructing the topology of a gene regulatory network is one of the key tasks in systems biology. Despite of the wide variety of proposed methods, very little work has been dedicated to the assessment of their stability properties. Here we present a methodical comparison of the performance of a novel method (RegnANN) for gene network inference based on multilayer perceptrons with three reference algorithms (ARACNE, CLR, KELLER), focussing our analysis on the prediction variability induced by both the network intrinsic structure and the available data. Results: The extensive evaluation on both synthetic data and a selection of gene modules of "Escherichia coli" indicates that all the algorithms suffer of instability and variability issues with regards to the reconstruction of the topology of the network. This instability makes objectively very hard the task of establishing which method performs best. Nevertheless, RegnANN shows MCC scores that compare very favorably with all the other inference methods tested. Availability: The software for the RegnANN inference algorithm is distributed under GPL3 and it is available at the corresponding author home page (http://mpba.fbk.eu/grimaldi/regnann-supmat

Direct constraint on the distance of y2 Velorum from AMBER/VLTI observations

In this work, we present the first AMBER observations, of the Wolf-Rayet and O (WR+O) star binary system y2 Velorum. The AMBER instrument was used with the telescopes UT2, UT3, and UT4 on baselines ranging from 46m to 85m. It delivered spectrally dispersed visibilities, as well as differential and closure phases, with a resolution R = 1500 in the spectral band 1.95-2.17 micron. We interpret these data in the context of a binary system with unresolved components, neglecting in a first approximation the wind-wind collision zone flux contribution. We show that the AMBER observables result primarily from the contribution of the individual components of the WR+O binary system. We discuss several interpretations of the residuals, and speculate on the detection of an additional continuum component, originating from the free-free emission associated with the wind-wind collision zone (WWCZ), and contributing at most to the observed K-band flux at the 5% level. The expected absolute separation and position angle at the time of observations were 5.1±0.9mas and 66±15° respectively. However, we infer a separation of 3.62+0.11-0.30 mas and a position angle of 73+9-11°. Our analysis thus implies that the binary system lies at a distance of 368+38-13 pc, in agreement with recent spectrophotometric estimates, but significantly larger than the Hipparcos value of 258+41-31 pc

Detecting a stochastic gravitational wave background with the Laser Interferometer Space Antenna

The random superposition of many weak sources will produce a stochastic background of gravitational waves that may dominate the response of the LISA (Laser Interferometer Space Antenna) gravitational wave observatory. Unless something can be done to distinguish between a stochastic background and detector noise, the two will combine to form an effective noise floor for the detector. Two methods have been proposed to solve this problem. The first is to cross-correlate the output of two independent interferometers. The second is an ingenious scheme for monitoring the instrument noise by operating LISA as a Sagnac interferometer. Here we derive the optimal orbital alignment for cross-correlating a pair of LISA detectors, and provide the first analytic derivation of the Sagnac sensitivity curve.Comment: 9 pages, 11 figures. Significant changes to the noise estimate

Search for Relativistic Magnetic Monopoles with IceCube

We present the first results in the search for relativistic magnetic monopoles with the IceCube detector, a subsurface neutrino telescope located in the South Polar ice cap containing a volume of 1 km$^{3}$. This analysis searches data taken on the partially completed detector during 2007 when roughly 0.2 km$^{3}$ of ice was instrumented. The lack of candidate events leads to an upper limit on the flux of relativistic magnetic monopoles of \Phi_{\mathrm{90%C.L.}}\sim 3\e{-18}\fluxunits for $\beta\geq0.8$. This is a factor of 4 improvement over the previous best experimental flux limits up to a Lorentz boost $\gamma$ below $10^{7}$. This result is then interpreted for a wide range of mass and kinetic energy values.Comment: 11 pages, 11 figures. v2 is minor text edits, no changes to resul