Regression with strongly correlated data

This paper discusses linear regression of strongly correlated data that arises, for example, in magnetohydrodynamic equilibrium reconstructions. We have proved that, generically, the covariance matrix of the estimated regression parameters for fixed sample size goes to zero as the correlations become unity. That is, in this limit the estimated parameters are known with perfect accuracy. Simple examples are shown to illustrate this effect and the nature of the exceptional cases in which the estimate covariance does not go to zero

Online open neuroimaging mass meta-analysis

We describe a system for meta-analysis where a wiki stores numerical data in a simple format and a web service performs the numerical computation. We initially apply the system on multiple meta-analyses of structural neuroimaging data results. The described system allows for mass meta-analysis, e.g., meta-analysis across multiple brain regions and multiple mental disorders.Comment: 5 pages, 4 figures SePublica 2012, ESWC 2012 Workshop, 28 May 2012, Heraklion, Greec

Modeling Maxwell's demon with a microcanonical Szilard engine

Following recent work by Marathe and Parrondo [PRL, 104, 245704 (2010)], we construct a classical Hamiltonian system whose energy is reduced during the adiabatic cycling of external parameters, when initial conditions are sampled microcanonically. Combining our system with a device that measures its energy, we propose a cyclic procedure during which energy is extracted from a heat bath and converted to work, in apparent violation of the second law of thermodynamics. This paradox is resolved by deriving an explicit relationship between the average work delivered during one cycle of operation, and the average information gained when measuring the system's energy

Turning Characteristics of U.S. Coast Guard Icebreaker M-5

http://deepblue.lib.umich.edu/bitstream/2027.42/96616/1/39015087359009.pd

Binary inspiral, gravitational radiation, and cosmology

Observations of binary inspiral in a single interferometric gravitational wave detector can be cataloged according to signal-to-noise ratio $\rho$ and chirp mass $\cal M$. The distribution of events in a catalog composed of observations with $\rho$ greater than a threshold $\rho_0$ depends on the Hubble expansion, deceleration parameter, and cosmological constant, as well as the distribution of component masses in binary systems and evolutionary effects. In this paper I find general expressions, valid in any homogeneous and isotropic cosmological model, for the distribution with $\rho$ and $\cal M$ of cataloged events; I also evaluate these distributions explicitly for relevant matter-dominated Friedmann-Robertson-Walker models and simple models of the neutron star mass distribution. In matter dominated Friedmann-Robertson-Walker cosmological models advanced LIGO detectors will observe binary neutron star inspiral events with $\rho>8$ from distances not exceeding approximately $2\,\text{Gpc}$, corresponding to redshifts of $0.48$ (0.26) for $h=0.8$ ($0.5$), at an estimated rate of 1 per week. As the binary system mass increases so does the distance it can be seen, up to a limit: in a matter dominated Einstein-deSitter cosmological model with $h=0.8$ ($0.5$) that limit is approximately $z=2.7$ (1.7) for binaries consisting of two $10\,\text{M}_\odot$ black holes. Cosmological tests based on catalogs of the kind discussed here depend on the distribution of cataloged events with $\rho$ and $\cal M$. The distributions found here will play a pivotal role in testing cosmological models against our own universe and in constructing templates for the detection of cosmological inspiraling binary neutron stars and black holes.Comment: REVTeX, 38 pages, 9 (encapsulated) postscript figures, uses epsf.st

Topological Entropy of Braids on the Torus

A fast method is presented for computing the topological entropy of braids on the torus. This work is motivated by the need to analyze large braids when studying two-dimensional flows via the braiding of a large number of particle trajectories. Our approach is a generalization of Moussafir's technique for braids on the sphere. Previous methods for computing topological entropies include the Bestvina--Handel train-track algorithm and matrix representations of the braid group. However, the Bestvina--Handel algorithm quickly becomes computationally intractable for large braid words, and matrix methods give only lower bounds, which are often poor for large braids. Our method is computationally fast and appears to give exponential convergence towards the exact entropy. As an illustration we apply our approach to the braiding of both periodic and aperiodic trajectories in the sine flow. The efficiency of the method allows us to explore how much extra information about flow entropy is encoded in the braid as the number of trajectories becomes large.Comment: 19 pages, 44 figures. SIAM journal styl

Wetting and Minimal Surfaces

We study minimal surfaces which arise in wetting and capillarity phenomena. Using conformal coordinates, we reduce the problem to a set of coupled boundary equations for the contact line of the fluid surface, and then derive simple diagrammatic rules to calculate the non-linear corrections to the Joanny-de Gennes energy. We argue that perturbation theory is quasi-local, i.e. that all geometric length scales of the fluid container decouple from the short-wavelength deformations of the contact line. This is illustrated by a calculation of the linearized interaction between contact lines on two opposite parallel walls. We present a simple algorithm to compute the minimal surface and its energy based on these ideas. We also point out the intriguing singularities that arise in the Legendre transformation from the pure Dirichlet to the mixed Dirichlet-Neumann problem.Comment: 22 page

Improving the efficiency of the detection of gravitational wave signals from inspiraling compact binaries: Chebyshev interpolation

Inspiraling compact binaries are promising sources of gravitational waves for ground and space-based laser interferometric detectors. The time-dependent signature of these sources in the detectors is a well-characterized function of a relatively small number of parameters; thus, the favored analysis technique makes use of matched filtering and maximum likelihood methods. Current analysis methodology samples the matched filter output at parameter values chosen so that the correlation between successive samples is 97% for which the filtered output is closely correlated. Here we describe a straightforward and practical way of using interpolation to take advantage of the correlation between the matched filter output associated with nearby points in the parameter space to significantly reduce the number of matched filter evaluations without sacrificing the efficiency with which real signals are recognized. Because the computational cost of the analysis is driven almost exclusively by the matched filter evaluations, this translates directly into an increase in computational efficiency, which in turn, translates into an increase in the size of the parameter space that can be analyzed and, thus, the science that can be accomplished with the data. As a demonstration we compare the present "dense sampling" analysis methodology with our proposed "interpolation" methodology, restricted to one dimension of the multi-dimensional analysis problem. We find that the interpolated search reduces by 25% the number of filter evaluations required by the dense search with 97% correlation to achieve the same efficiency of detection for an expected false alarm probability. Generalized to higher dimensional space of a generic binary including spins suggests an order of magnitude increase in computational efficiency.Comment: 23 pages, 5 figures, submitted to Phys. Rev.

Gravitational Waves from coalescing binaries: Estimation of parameters

The paper presents a statistical model which reproduces the results of Monte Carlo simulations to estimate the parameters of the gravitational wave signal from a coalesing binary system. The model however is quite general and would be useful in other parameter estimation problems.Comment: LaTeX with RevTeX macros, 4 figure