Lyapunov exponents as a dynamical indicator of a phase transition

We study analytically the behavior of the largest Lyapunov exponent $\lambda_1$ for a one-dimensional chain of coupled nonlinear oscillators, by combining the transfer integral method and a Riemannian geometry approach. We apply the results to a simple model, proposed for the DNA denaturation, which emphasizes a first order-like or second order phase transition depending on the ratio of two length scales: this is an excellent model to characterize $\lambda_1$ as a dynamical indicator close to a phase transition.Comment: 8 Pages, 3 Figure

Proportion of various types of thyroid disorders among newborns with congenital hypothyroidism and normally located gland: A regional cohort study

Objective To determine the proportion of the various types of thyroid disorders among newborns detected by the neonatal TSH screening programme, with a normally located thyroid gland.Patients and methods Of the 882 575 infants screened in our centre between 1981 and 2002, 85 infants with a normally located gland had persistent elevation of serum TSH values (an incidence of 1/10 383). Six of these 85 patients were lost to follow-up and were therefore excluded from the study. During follow-up, patients were classified as having permanent or transient hypothyroidism.Results Among the 79 patients included in the study, transient (n = 30, 38% of cases) and permanent (n = 49, 62% of cases) congenital hypothyroidism (CH) was demonstrated during the follow-up at the age of 0.7 +/- 0.6 years and 2.6 +/- 1.8 years (P < 0.0001), respectively. The proportion of premature births was significantly higher in the group with transient CH (57%) than in the group with permanent CH (2%) (P < 0.0001). A history of iatrogenic iodine overload was identified during the neonatal period in 69% of transient cases. Among permanent CH cases (n = 49), patients were classified as having a goitre (n = 27, 55% of cases), a normal sized and shaped thyroid gland (n = 14, 29% of cases) or a hypoplastic gland (n = 8, 16% of cases). The latter patients demonstrated global thyroid hypoplasia (n = 3), a right hemithyroid (n = 2), hypoplasia of the left lobe (n = 2), or asymmetry in the location of the two lobes (n = 1). Patients with a normal sized and shaped thyroid gland showed a significantly less severe form of hypothyroidism than those with a goitre or a hypoplastic thyroid gland (P < 0.0002). Among permanent CH cases, those with a goitre (n = 27) had an iodine organification defect (n = 10), Pendred syndrome (n = 1), a defect of thyroglobulin synthesis (n = 8), or a defect of sodium iodine symporter (n = 1), and in seven patients no aetiology could be determined. Among permanent cases with a normal sized and shaped thyroid gland (n = 14), a specific aetiology was found in only one patient (pseudohypoparathyroidism) and two patients had Down's syndrome. Among those with a globally hypoplastic gland, a TSH receptor gene mutation was found in two patients.Conclusions A precise description of the phenotype can enhance our understanding of various forms of neonatal hypothyroidism as well as their prevalence and management. It also helps to identify cases of congenital hypothyroidism of unknown aetiology, which will need to be investigated in collaboration with molecular biologists

Statistics of first-passage times in disordered systems using backward master equations and their exact renormalization rules

We consider the non-equilibrium dynamics of disordered systems as defined by a master equation involving transition rates between configurations (detailed balance is not assumed). To compute the important dynamical time scales in finite-size systems without simulating the actual time evolution which can be extremely slow, we propose to focus on first-passage times that satisfy 'backward master equations'. Upon the iterative elimination of configurations, we obtain the exact renormalization rules that can be followed numerically. To test this approach, we study the statistics of some first-passage times for two disordered models : (i) for the random walk in a two-dimensional self-affine random potential of Hurst exponent $H$, we focus on the first exit time from a square of size $L \times L$ if one starts at the square center. (ii) for the dynamics of the ferromagnetic Sherrington-Kirkpatrick model of $N$ spins, we consider the first passage time $t_f$ to zero-magnetization when starting from a fully magnetized configuration. Besides the expected linear growth of the averaged barrier $\bar{\ln t_{f}} \sim N$, we find that the rescaled distribution of the barrier $(\ln t_{f})$ decays as $e^{- u^{\eta}}$ for large $u$ with a tail exponent of order $\eta \simeq 1.72$. This value can be simply interpreted in terms of rare events if the sample-to-sample fluctuation exponent for the barrier is $\psi_{width}=1/3$.Comment: 8 pages, 4 figure

The Theoretical Astrophysical Observatory: Cloud-Based Mock Galaxy Catalogues

We introduce the Theoretical Astrophysical Observatory (TAO), an online virtual laboratory that houses mock observations of galaxy survey data. Such mocks have become an integral part of the modern analysis pipeline. However, building them requires an expert knowledge of galaxy modelling and simulation techniques, significant investment in software development, and access to high performance computing. These requirements make it difficult for a small research team or individual to quickly build a mock catalogue suited to their needs. To address this TAO offers access to multiple cosmological simulations and semi-analytic galaxy formation models from an intuitive and clean web interface. Results can be funnelled through science modules and sent to a dedicated supercomputer for further processing and manipulation. These modules include the ability to (1) construct custom observer light-cones from the simulation data cubes; (2) generate the stellar emission from star formation histories, apply dust extinction, and compute absolute and/or apparent magnitudes; and (3) produce mock images of the sky. All of TAO's features can be accessed without any programming requirements. The modular nature of TAO opens it up for further expansion in the future.Comment: 17 pages, 11 figures, 2 tables; accepted for publication in ApJS. The Theoretical Astrophysical Observatory (TAO) is now open to the public at https://tao.asvo.org.au/. New simulations, models and tools will be added as they become available. Contact [email protected] if you have data you would like to make public through TAO. Feedback and suggestions are very welcom

Semi-Analytic Galaxy Evolution (SAGE): Model Calibration and Basic Results

This paper describes a new publicly available codebase for modelling galaxy formation in a cosmological context, the "Semi-Analytic Galaxy Evolution" model, or SAGE for short. SAGE is a significant update to that used in Croton et al. (2006) and has been rebuilt to be modular and customisable. The model will run on any N-body simulation whose trees are organised in a supported format and contain a minimum set of basic halo properties. In this work we present the baryonic prescriptions implemented in SAGE to describe the formation and evolution of galaxies, and their calibration for three N-body simulations: Millennium, Bolshoi, and GiggleZ. Updated physics include: gas accretion, ejection due to feedback, and reincorporation via the galactic fountain; a new gas cooling--radio mode active galactic nucleus (AGN) heating cycle; AGN feedback in the quasar mode; a new treatment of gas in satellite galaxies; and galaxy mergers, disruption, and the build-up of intra-cluster stars. Throughout, we show the results of a common default parameterization on each simulation, with a focus on the local galaxy population.Comment: 15 pages, 9 figures, accepted for publication in ApJS. SAGE is a publicly available codebase for modelling galaxy formation in a cosmological context, available at https://github.com/darrencroton/sage Questions and comments can be sent to Darren Croton: [email protected]

Driven interfaces in random media at finite temperature : is there an anomalous zero-velocity phase at small external force ?

The motion of driven interfaces in random media at finite temperature $T$ and small external force $F$ is usually described by a linear displacement $h_G(t) \sim V(F,T) t$ at large times, where the velocity vanishes according to the creep formula as $V(F,T) \sim e^{-K(T)/F^{\mu}}$ for $F \to 0$. In this paper, we question this picture on the specific example of the directed polymer in a two dimensional random medium. We have recently shown (C. Monthus and T. Garel, arxiv:0802.2502) that its dynamics for F=0 can be analyzed in terms of a strong disorder renormalization procedure, where the distribution of renormalized barriers flows towards some "infinite disorder fixed point". In the present paper, we obtain that for small $F$, this "infinite disorder fixed point" becomes a "strong disorder fixed point" with an exponential distribution of renormalized barriers. The corresponding distribution of trapping times then only decays as a power-law $P(\tau) \sim 1/\tau^{1+\alpha}$, where the exponent $\alpha(F,T)$ vanishes as $\alpha(F,T) \propto F^{\mu}$ as $F \to 0$. Our conclusion is that in the small force region $\alpha(F,T)<1$, the divergence of the averaged trapping time $\bar{\tau}=+\infty$ induces strong non-self-averaging effects that invalidate the usual creep formula obtained by replacing all trapping times by the typical value. We find instead that the motion is only sub-linearly in time $h_G(t) \sim t^{\alpha(F,T)}$, i.e. the asymptotic velocity vanishes V=0. This analysis is confirmed by numerical simulations of a directed polymer with a metric constraint driven in a traps landscape. We moreover obtain that the roughness exponent, which is governed by the equilibrium value $\zeta_{eq}=2/3$ up to some large scale, becomes equal to $\zeta=1$ at the largest scales.Comment: v3=final versio

Interaction of subducted slabs with the mantle transition-zone: A regime diagram from 2-D thermo-mechanical models with a mobile trench and an overriding plate

Transition zone slab deformation influences Earth's thermal, chemical, and tectonic evolution. However, the mechanisms responsible for the wide range of imaged slab morphologies remain debated. Here we use 2-D thermo-mechanical models with a mobile trench, an overriding plate, a temperature and stress-dependent rheology, and a 10, 30, or 100-fold increase in lower mantle viscosity, to investigate the effect of initial subducting and overriding-plate ages on slab-transition zone interaction. Four subduction styles emerge: (i) a "vertical folding" mode, with a quasi-stationary trench, near-vertical subduction, and buckling/folding at depth (VF); (ii) slabs that induce mild trench retreat, which are flattened/"horizontally deflected" and stagnate at the upper-lower mantle interface (HD); (iii) inclined slabs, which result from rapid sinking and strong trench retreat (ISR); (iv) a two-stage mode, displaying backward-bent and subsequently inclined slabs, with late trench retreat (BIR). Transitions from regime (i) to (iii) occur with increasing subducting plate age (i.e., buoyancy and strength). Regime (iv) develops for old (strong) subducting and overriding plates. We find that the interplay between trench motion and slab deformation at depth dictates the subduction style, both being controlled by slab strength, which is consistent with predictions from previous compositional subduction models. However, due to feedbacks between deformation, sinking rate, temperature, and slab strength, the subducting plate buoyancy, overriding plate strength, and upper-lower mantle viscosity jump are also important controls in thermo-mechanical subduction. For intermediate upper-lower mantle viscosity jumps (×30), our regimes reproduce the diverse range of seismically imaged slab morphologies

THERMODYNAMICS OF A BROWNIAN BRIDGE POLYMER MODEL IN A RANDOM ENVIRONMENT

We consider a directed random walk making either 0 or $+1$ moves and a Brownian bridge, independent of the walk, conditioned to arrive at point $b$ on time $T$. The Hamiltonian is defined as the sum of the square of increments of the bridge between the moments of jump of the random walk and interpreted as an energy function over the bridge connfiguration; the random walk acts as the random environment. This model provides a continuum version of a model with some relevance to protein conformation. The thermodynamic limit of the specific free energy is shown to exist and to be self-averaging, i.e. it is equal to a trivial --- explicitly computed --- random variable. An estimate of the asymptotic behaviour of the ground state energy is also obtained.Comment: 20 pages, uuencoded postscrip

Gas of self-avoiding loops on the brickwork lattice

An exact calculation of the phase diagram for a loop gas model on the brickwork lattice is presented. The model includes a bending energy. In the dense limit, where all the lattice sites are occupied, a phase transition occuring at an asymmetric Lifshitz tricritical point is observed as the temperature associated with the bending energy is varied. Various critical exponents are calculated. At lower densities, two lines of transitions (in the Ising universality class) are observed, terminated by a tricritical point, where there is a change in the modulation of the correlation function. To each tricritical point an associated disorder line is found.Comment: 19 pages, 6 figures. to appear in J. Phys. A : Math. & Ge

Single Chain Force Spectroscopy: Sequence Dependence

We study the elastic properties of a single A/B copolymer chain with a specific sequence. We predict a rich structure in the force extension relations which can be addressed to the sequence. The variational method is introduced to probe local minima on the path of stretching and releasing. At given force, we find multiple configurations which are separated by energy barriers. A collapsed globular configuration consists of several domains which unravel cooperatively. Upon stretching, unfolding path shows stepwise pattern corresponding to the unfolding of each domain. While releasing, several cores can be created simultaneously in the middle of the chain resulting in a different path of collapse.Comment: 6 pages 3 figure