Global stability analysis of birhythmicity in a self-sustained oscillator

We analyze global stability properties of birhythmicity in a self-sustained system with random excitations. The model is a multi-limit cycles variation of the van der Pol oscillatorintroduced to analyze enzymatic substrate reactions in brain waves. We show that the two frequencies are strongly influenced by the nonlinear coefficients $\alpha$ and $\beta$. With a random excitation, such as a Gaussian white noise, the attractor's global stability is measured by the mean escape time $\tau$ from one limit-cycle. An effective activation energy barrier is obtained by the slope of the linear part of the variation of the escape time $\tau$ versus the inverse noise-intensity 1/D. We find that the trapping barriers of the two frequencies can be very different, thus leaving the system on the same attractor for an overwhelming time. However, we also find that the system is nearly symmetric in a narrow range of the parameters.Comment: 17 pages, 8 figures, to appear on Choas, 201

Interface Motion in Random Media at Finite Temperature

We have studied numerically the dynamics of a driven elastic interface in a random medium, focusing on the thermal rounding of the depinning transition and on the behavior in the $T=0$ pinned phase. Thermal effects are quantitatively more important than expected from simple dimensional estimates. For sufficient low temperature the creep velocity at a driving force equal to the $T=0$ depinning force exhibits a power-law dependence on $T$, in agreement with earlier theoretical and numerical predictions for CDW's. We have also examined the dynamics in the $T=0$ pinned phase resulting from slowly increasing the driving force towards threshold. The distribution of avalanche sizes $S_\|$ decays as $S_\|^{-1-\kappa}$, with $\kappa = 0.05\pm 0.05$, in agreement with recent theoretical predictions.Comment: harvmac.tex, 30 pages, including 9 figures, available upon request. SU-rm-94073

Computational Complexity of Determining the Barriers to Interface Motion in Random Systems

The low-temperature driven or thermally activated motion of several condensed matter systems is often modeled by the dynamics of interfaces (co-dimension-1 elastic manifolds) subject to a random potential. Two characteristic quantitative features of the energy landscape of such a many-degree-of-freedom system are the ground-state energy and the magnitude of the energy barriers between given configurations. While the numerical determination of the former can be accomplished in time polynomial in the system size, it is shown here that the problem of determining the latter quantity is NP-complete. Exact computation of barriers is therefore (almost certainly) much more difficult than determining the exact ground states of interfaces.Comment: 8 pages, figures included, to appear in Phys. Rev.

Determining Pair Interactions from Structural Correlations

We examine metastable configurations of a two-dimensional system of interacting particles on a quenched random potential landscape and ask how the configurational pair correlation function is related to the particle interactions and the statistical properties of the potential landscape. Understanding this relation facilitates quantitative studies of magnetic flux line interactions in type II superconductors, using structural information available from Lorentz microscope images or Bitter decorations. Previous work by some of us supported the conjecture that the relationship between pair correlations and interactions in pinned flux line ensembles is analogous to the corresponding relationship in the theory of simple liquids. The present paper aims at a more thorough understanding of this relation. We report the results of numerical simulations and present a theory for the low density behavior of the pair correlation function which agrees well with our simulations and captures features observed in experiments. In particular, we find that the resulting description goes beyond the conjectured classical liquid type relation and we remark on the differences.Comment: 7 pages, 6 figures. See also http://rainbow.uchicago.edu/~grier

The ARIEL Instrument Control Unit design for the M4 Mission Selection Review of the ESA's Cosmic Vision Program

The Atmospheric Remote-sensing Infrared Exoplanet Large-survey mission (ARIEL) is one of the three present candidates for the ESA M4 (the fourth medium mission) launch opportunity. The proposed Payload will perform a large unbiased spectroscopic survey from space concerning the nature of exoplanets atmospheres and their interiors to determine the key factors affecting the formation and evolution of planetary systems. ARIEL will observe a large number (>500) of warm and hot transiting gas giants, Neptunes and super-Earths around a wide range of host star types, targeting planets hotter than 600 K to take advantage of their well-mixed atmospheres. It will exploit primary and secondary transits spectroscopy in the 1.2-8 um spectral range and broad-band photometry in the optical and Near IR (NIR). The main instrument of the ARIEL Payload is the IR Spectrometer (AIRS) providing low-resolution spectroscopy in two IR channels: Channel 0 (CH0) for the 1.95-3.90 um band and Channel 1 (CH1) for the 3.90-7.80 um range. It is located at the intermediate focal plane of the telescope and common optical system and it hosts two IR sensors and two cold front-end electronics (CFEE) for detectors readout, a well defined process calibrated for the selected target brightness and driven by the Payload's Instrument Control Unit (ICU).Comment: Experimental Astronomy, Special Issue on ARIEL, (2017

Simulation of the Zero Temperature Behavior of a 3-Dimensional Elastic Medium

We have performed numerical simulation of a 3-dimensional elastic medium, with scalar displacements, subject to quenched disorder. We applied an efficient combinatorial optimization algorithm to generate exact ground states for an interface representation. Our results indicate that this Bragg glass is characterized by power law divergences in the structure factor $S(k)\sim A k^{-3}$. We have found numerically consistent values of the coefficient $A$ for two lattice discretizations of the medium, supporting universality for $A$ in the isotropic systems considered here. We also examine the response of the ground state to the change in boundary conditions that corresponds to introducing a single dislocation loop encircling the system. Our results indicate that the domain walls formed by this change are highly convoluted, with a fractal dimension $d_f=2.60(5)$. We also discuss the implications of the domain wall energetics for the stability of the Bragg glass phase. As in other disordered systems, perturbations of relative strength $\delta$ introduce a new length scale $L^* \sim \delta^{-1/\zeta}$ beyond which the perturbed ground state becomes uncorrelated with the reference (unperturbed) ground state. We have performed scaling analysis of the response of the ground state to the perturbations and obtain $\zeta = 0.385(40)$. This value is consistent with the scaling relation $\zeta=d_f/2- \theta$, where $\theta$ characterizes the scaling of the energy fluctuations of low energy excitations.Comment: 20 pages, 13 figure

Scaling of interfaces in brittle fracture and perfect plasticity

The roughness properties of two-dimensional fracture surfaces as created by the slow failure of random fuse networks are considered and compared to yield surfaces of perfect plasticity with similar disorder. By studying systems up to a linear size L=350 it is found that in the cases studied the fracture surfaces exhibit self-affine scaling with a roughness exponent close to 2/3, which is asymptotically exactly true for plasticity though finite-size effects are evident for both. The overlap of yield or minimum energy and fracture surfaces with exactly the same disorder configuration is shown to be a decreasing function of the system size and to be of a rather large magnitude for all cases studied. The typical ``overlap cluster'' length between pairs of such interfaces converges to a constant with $L$ increasing.Comment: Accepted for publication in Phys. Rev.

Hysteresis and hierarchies: dynamics of disorder-driven first-order phase transformations

We use the zero-temperature random-field Ising model to study hysteretic behavior at first-order phase transitions. Sweeping the external field through zero, the model exhibits hysteresis, the return-point memory effect, and avalanche fluctuations. There is a critical value of disorder at which a jump in the magnetization (corresponding to an infinite avalanche) first occurs. We study the universal behavior at this critical point using mean-field theory, and also present preliminary results of numerical simulations in three dimensions.Comment: 12 pages plus 2 appended figures, plain TeX, CU-MSC-747