Closed-loop estimation of retinal network sensitivity reveals signature of efficient coding

According to the theory of efficient coding, sensory systems are adapted to represent natural scenes with high fidelity and at minimal metabolic cost. Testing this hypothesis for sensory structures performing non-linear computations on high dimensional stimuli is still an open challenge. Here we develop a method to characterize the sensitivity of the retinal network to perturbations of a stimulus. Using closed-loop experiments, we explore selectively the space of possible perturbations around a given stimulus. We then show that the response of the retinal population to these small perturbations can be described by a local linear model. Using this model, we computed the sensitivity of the neural response to arbitrary temporal perturbations of the stimulus, and found a peak in the sensitivity as a function of the frequency of the perturbations. Based on a minimal theory of sensory processing, we argue that this peak is set to maximize information transmission. Our approach is relevant to testing the efficient coding hypothesis locally in any context where no reliable encoding model is known

Generating Log-normal Mock Catalog of Galaxies in Redshift Space

We present a public code to generate a mock galaxy catalog in redshift space assuming a log-normal probability density function (PDF) of galaxy and matter density fields. We draw galaxies by Poisson-sampling the log-normal field, and calculate the velocity field from the linearised continuity equation of matter fields, assuming zero vorticity. This procedure yields a PDF of the pairwise velocity fields that is qualitatively similar to that of N-body simulations. We check fidelity of the catalog, showing that the measured two-point correlation function and power spectrum in real space agree with the input precisely. We find that a linear bias relation in the power spectrum does not guarantee a linear bias relation in the density contrasts, leading to a cross-correlation coefficient of matter and galaxies deviating from unity on small scales. We also find that linearising the Jacobian of the real-to-redshift space mapping provides a poor model for the two-point statistics in redshift space. That is, non-linear redshift-space distortion is dominated by non-linearity in the Jacobian. The power spectrum in redshift space shows a damping on small scales that is qualitatively similar to that of the well-known Fingers-of-God (FoG) effect due to random velocities, except that the log-normal mock does not include random velocities. This damping is a consequence of non-linearity in the Jacobian, and thus attributing the damping of the power spectrum solely to FoG, as commonly done in the literature, is misleading.Comment: 38 pages, 16 figures, code publicly available as "lognormal_galaxies" at http://wwwmpa.mpa-garching.mpg.de/~komatsu/codes.html Matches published version : added figures and explanatory comment

On Validating an Astrophysical Simulation Code

We present a case study of validating an astrophysical simulation code. Our study focuses on validating FLASH, a parallel, adaptive-mesh hydrodynamics code for studying the compressible, reactive flows found in many astrophysical environments. We describe the astrophysics problems of interest and the challenges associated with simulating these problems. We describe methodology and discuss solutions to difficulties encountered in verification and validation. We describe verification tests regularly administered to the code, present the results of new verification tests, and outline a method for testing general equations of state. We present the results of two validation tests in which we compared simulations to experimental data. The first is of a laser-driven shock propagating through a multi-layer target, a configuration subject to both Rayleigh-Taylor and Richtmyer-Meshkov instabilities. The second test is a classic Rayleigh-Taylor instability, where a heavy fluid is supported against the force of gravity by a light fluid. Our simulations of the multi-layer target experiments showed good agreement with the experimental results, but our simulations of the Rayleigh-Taylor instability did not agree well with the experimental results. We discuss our findings and present results of additional simulations undertaken to further investigate the Rayleigh-Taylor instability.Comment: 76 pages, 26 figures (3 color), Accepted for publication in the ApJ

Self-Consistent Magnetic Stellar Evolution Models of the Detached, Solar-Type Eclipsing Binary EF Aquarii

We introduce a new one-dimensional stellar evolution code, based on the existing Dartmouth code, that self-consistently accounts for the presence of a globally pervasive magnetic field. The methods involved in perturbing the equations of stellar structure, the equation of state, and the mixing-length theory of convection are presented and discussed. As a first test of the code's viability, stellar evolution models are computed for the components of a solar-type, detached eclipsing binary (DEB) system, EF Aquarii, shown to exhibit large disagreements with stellar models. The addition of the magnetic perturbation corrects the radius and effective temperature discrepancies observed in EF Aquarii. Furthermore, the required magnetic field strength at the model photosphere is within a factor of two of the magnetic field strengths estimated from the stellar X-ray luminosities measured by ROSAT and those predicted from Ca II K line core emission. These models provide firm evidence that the suppression of thermal convection arising from the presence of a magnetic field is sufficient to significantly alter the structure of solar-type stars, producing noticeably inflated radii and cooler effective temperatures. The inclusion of magnetic effects within a stellar evolution model has a wide range of applications, from DEBs and exoplanet host stars to the donor stars of cataclysmic variables.Comment: Accepted for publication in ApJ, 15 pages, 3 figures; Misprints are corrected in version

Cyclic Variability of the Circumstellar Disc of the Be Star ζ\zeta Tau. II. Testing the 2D Global Disc Oscillation Model

Aims. In this paper we model, in a self-consistent way, polarimetric, photometric, spectrophotometric and interferometric observations of the classical Be star $\zeta$ Tauri. Our primary goal is to conduct a critical quantitative test of the global oscillation scenario. Methods. We have carried out detailed three-dimensional, NLTE radiative transfer calculations using the radiative transfer code HDUST. For the input for the code we have used the most up-to-date research on Be stars to include a physically realistic description for the central star and the circumstellar disc. We adopt a rotationally deformed, gravity darkened central star, surrounded by a disc whose unperturbed state is given by a steady-state viscous decretion disc model. We further assume that disc is in vertical hydrostatic equilibrium. Results. By adopting a viscous decretion disc model for $\zeta$ Tauri and a rigorous solution of the radiative transfer, we have obtained a very good fit of the time-average properties of the disc. This provides strong theoretical evidence that the viscous decretion disc model is the mechanism responsible for disc formation. With the global oscillation model we have successfully fitted spatially resolved VLTI/AMBER observations and the temporal V/R variations of the H$\alpha$ and Br$\gamma$ lines. This result convincingly demonstrates that the oscillation pattern in the disc is a one-armed spiral. Possible model shortcomings, as well as suggestions for future improvements, are also discussed.Comment: 14 pages, 9 figures, accepted to A&

Oscillations of rapidly rotating relativistic stars

Non-axisymmetric oscillations of rapidly rotating relativistic stars are studied using the Cowling approximation. The oscillation spectra have been estimated by Fourier transforming the evolution equations describing the perturbations. This is the first study of its kind and provides information on the effect of fast rotation on the oscillation spectra while it offers the possibility in studying the complete problem by including spacetime perturbations. Our study includes both axisymmetric and non-axisymmetric perturbations and provides limits for the onset of the secular bar mode rotational instability. We also present approximate formulae for the dependence of the oscillation spectrum from rotation. The results suggest that it is possible to extract the relativistic star's parameters from the observed gravitational wave spectrum.Comment: this article will be published in Physical Review

Realistic evaluation of hull performance for rowing shells, canoes, and kayaks in unsteady flow

In this study, we investigated the effect of hull dynamics in shallow water on the hydrodynamic performance of rowing shells as well as canoes and kayaks. An approach was developed to generate data in a towing tank using a test rig capable of reproducing realistic speed profiles. The impact of unsteady shallow-water effects on wave-making resistance was examined via experimental measurements on a benchmark hull. The data generated were used to explore the validity of a computational approach developed to predict unsteady shallow-water wave resistance. Comparison of measured and predicted results showed that the computational approach correctly predicted complex unsteady wave-resistance phenomena at low oscillation frequency and speed, but that total resistance was substantially under-predicted at moderate oscillation frequency and speed. It was postulated that this discrepancy arose from unsteady viscous effects. This was investigated via hot-film measurements for a full-scale single scull in unsteady flow in both towing-tank and field-trial conditions. Results suggested a strong link between acceleration and turbulence and demonstrated that the measured real-world viscous-flow behaviour could be successfully reproduced in the tank. Thus a suitable tank-test approach could provide a reliable guide to hull performance characterization in unsteady flow