Values of H_0 from Models of the Gravitational Lens 0957+561

The lensed double QSO 0957+561 has a well-measured time delay and hence is useful for a global determination of H0. Uncertainty in the mass distribution of the lens is the largest source of uncertainty in the derived H0. We investigate the range of \hn produced by a set of lens models intended to mimic the full range of astrophysically plausible mass distributions, using as constraints the numerous multiply-imaged sources which have been detected. We obtain the first adequate fit to all the observations, but only if we include effects from the galaxy cluster beyond a constant local magnification and shear. Both the lens galaxy and the surrounding cluster must depart from circular symmetry as well. Lens models which are consistent with observations to 95% CL indicate H0=104^{+31}_{-23}(1-\kthirty) km/s/Mpc. Previous weak lensing measurements constrain the mean mass density within 30" of G1 to be kthirty=0.26+/-0.16 (95% CL), implying H0=77^{+29}_{-24}km/s/Mpc (95% CL). The best-fitting models span the range 65--80 km/s/Mpc. Further observations will shrink the confidence interval for both the mass model and \kthirty. The range of H0 allowed by the full gamut of our lens models is substantially larger than that implied by limiting consideration to simple power law density profiles. We therefore caution against use of simple isothermal or power-law mass models in the derivation of H0 from other time-delay systems. High-S/N imaging of multiple or extended lensed features will greatly reduce the H0 uncertainties when fitting complex models to time-delay lenses.Comment: AASTEX, 48 pages 4 figures, 2 tables. Also available at: http://www.astro.lsa.umich.edu:80/users/philf/www/papers/list.htm

Spherical harmonic decomposition applied to spatial-temporal analysis of human high-density EEG

We demonstrate an application of spherical harmonic decomposition to analysis of the human electroencephalogram (EEG). We implement two methods and discuss issues specific to analysis of hemispherical, irregularly sampled data. Performance of the methods and spatial sampling requirements are quantified using simulated data. The analysis is applied to experimental EEG data, confirming earlier reports of an approximate frequency-wavenumber relationship in some bands.Comment: 12 pages, 8 figures, submitted to Phys. Rev. E, uses APS RevTeX style

Spatially-resolved spectrophotometric analysis and modelling of the Superantennae

We have performed spatially-resolved spectroscopy of the double-nucleated Ultra-Luminous Infrared Galaxy IRAS 19254-7245, ``the Superantennae'', along the line connecting the two nuclei. These data are analysed with a spectral synthesis code, to derive the star formation and extinction properties of the galaxy. The star formation history (SFH) of the two nuclei is similarly characterized by two different main episodes: a recent burst, responsible of the observed emission lines, and an older one, occurred roughly 1 Gyr ago. We tentatively associate this bimodal SFH with a double encounter in the dynamical history of the merger. We have complemented our study with a detailed analysis of the broad band spectral energy distribution of the Superantennae, built from published photometry, providing the separate optical-to-mm SEDs of the two nuclei. Our analysis shows that: a) the southern nucleus is responsible for about 80% of the total infrared luminosity of the system, b) the L-band luminosity in the southern nucleus is dominated by the emission from an obscured AGN, providing about 40 to 50% of the bolometric flux between 8 and 1000 microns; c) the northern nucleus does not show evidence for AGN emission and appears to be in a post-starburst phase. As for the relative strengths of the AGN and starburst components, we find that, while they are comparable at FIR and sub-mm wavelengths, in the radio the Sy2 emission dominates by an order of magnitude the starburst.Comment: 18 pages. Accepted for publication on A&

Simulated Annealing for Topological Solitons

The search for solutions of field theories allowing for topological solitons requires that we find the field configuration with the lowest energy in a given sector of topological charge. The standard approach is based on the numerical solution of the static Euler-Lagrange differential equation following from the field energy. As an alternative, we propose to use a simulated annealing algorithm to minimize the energy functional directly. We have applied simulated annealing to several nonlinear classical field theories: the sine-Gordon model in one dimension, the baby Skyrme model in two dimensions and the nuclear Skyrme model in three dimensions. We describe in detail the implementation of the simulated annealing algorithm, present our results and get independent confirmation of the studies which have used standard minimization techniques.Comment: 31 pages, LaTeX, better quality pics at http://www.phy.umist.ac.uk/~weidig/Simulated_Annealing/, updated for publicatio

Wavelet analysis of epileptic spikes

Interictal spikes and sharp waves in human EEG are characteristic signatures of epilepsy. These potentials originate as a result of synchronous, pathological discharge of many neurons. The reliable detection of such potentials has been the long standing problem in EEG analysis, especially after long-term monitoring became common in investigation of epileptic patients. The traditional definition of a spike is based on its amplitude, duration, sharpness, and emergence from its background. However, spike detection systems built solely around this definition are not reliable due to the presence of numerous transients and artifacts. We use wavelet transform to analyze the properties of EEG manifestations of epilepsy. We demonstrate that the behavior of wavelet transform of epileptic spikes across scales can constitute the foundation of a relatively simple yet effective detection algorithm.Comment: 4 pages, 3 figure

Maximum-Reward Motion in a Stochastic Environment: The Nonequilibrium Statistical Mechanics Perspective

We consider the problem of computing the maximum-reward motion in a reward field in an online setting. We assume that the robot has a limited perception range, and it discovers the reward field on the fly. We analyze the performance of a simple, practical lattice-based algorithm with respect to the perception range. Our main result is that, with very little perception range, the robot can collect as much reward as if it could see the whole reward field, under certain assumptions. Along the way, we establish novel connections between this class of problems and certain fundamental problems of nonequilibrium statistical mechanics . We demonstrate our results in simulation examples

On discretization in time in simulations of particulate flows

We propose a time discretization scheme for a class of ordinary differential equations arising in simulations of fluid/particle flows. The scheme is intended to work robustly in the lubrication regime when the distance between two particles immersed in the fluid or between a particle and the wall tends to zero. The idea consists in introducing a small threshold for the particle-wall distance below which the real trajectory of the particle is replaced by an approximated one where the distance is kept equal to the threshold value. The error of this approximation is estimated both theoretically and by numerical experiments. Our time marching scheme can be easily incorporated into a full simulation method where the velocity of the fluid is obtained by a numerical solution to Stokes or Navier-Stokes equations. We also provide a derivation of the asymptotic expansion for the lubrication force (used in our numerical experiments) acting on a disk immersed in a Newtonian fluid and approaching the wall. The method of this derivation is new and can be easily adapted to other cases

The Effect of the Reynolds Number on Lateral Migration of Nonneutrally-Buoyant Spherical Particles in Poiseuille Flow

The lateral migration of nonneutrally-buoyant spherical particles in Poiseuille flow is investigated numerically using the boundary element method. In particular, the steady, Navier-Stokes equations are solved using a classical domain integration method treating the nonlinear terms as pseudo-body forces. The numerical results for the lateral migration velocity are compared with experimental data. The numerical results indicate that the lateral migration velocity does not scale linearly with the Reynolds number. The methodology is extended to include non-Newtonian power-law fluids. The migration velocity is significantly affected for particles suspended in this class of fluids and can actually change direction for large values of the power-law index

Design and Control of Compliant Tensegrity Robots Through Simulation and Hardware Validation

To better understand the role of tensegrity structures in biological systems and their application to robotics, the Dynamic Tensegrity Robotics Lab at NASA Ames Research Center has developed and validated two different software environments for the analysis, simulation, and design of tensegrity robots. These tools, along with new control methodologies and the modular hardware components developed to validate them, are presented as a system for the design of actuated tensegrity structures. As evidenced from their appearance in many biological systems, tensegrity ("tensile-integrity") structures have unique physical properties which make them ideal for interaction with uncertain environments. Yet these characteristics, such as variable structural compliance, and global multi-path load distribution through the tension network, make design and control of bio-inspired tensegrity robots extremely challenging. This work presents the progress in using these two tools in tackling the design and control challenges. The results of this analysis includes multiple novel control approaches for mobility and terrain interaction of spherical tensegrity structures. The current hardware prototype of a six-bar tensegrity, code-named ReCTeR, is presented in the context of this validation