Heat conduction in a chain of harmonic and anharmonic oscillators under the presence of an energy conserving noise

Reproducing Fourier's law of heat conduction from a microscopic stochastic model is a long standing challenge in statistical physics. As was shown by Rieder, Lebowitz and Lieb many years ago, a chain of harmonically coupled oscillators connected to two heat baths at different temperatures does not reproduce the diffusive behaviour of Fourier's law, but instead a ballistic one with an infinite thermal conductivity. Since then, there has been a substantial effort from the scientific community in identifying the key mechanism necessary to reproduce such diffusivity, which usually revolved around anharmonicity and the effect of impurities. Recently, it was shown by Dhar, Venkateshan and Lebowitz that Fourier's law can be recovered by introducing an energy conserving noise, whose role is to simulate the elastic collisions between the atoms and other microscopic degrees of freedom, which one would expect to be present in a real solid. For a one-dimensional chain this is accomplished numerically by randomly flipping - under the framework of a Poisson process with a variable “rate of collisions" - the sign of the velocity of an oscillator. In this poster we present Langevin simulations of a one-dimensional chain of oscillators coupled to two heat baths at different temperatures. We consider both harmonic and anharmonic (quartic) interactions, which are studied with and without the energy conserving noise. With these results we are able to map in detail how the heat conductivity k is influenced by both anharmonicity and the energy conserving noise. We also present a detailed analysis of the behaviour of k as a function of the size of the system and the rate of collisions, which includes a finite-size scaling method that enables us to extract the relevant critical exponents. Finally, we show that for harmonic chains, k is independent of temperature, both with and without the noise. Conversely, for anharmonic chains we find that k increases roughly linearly with the temperature of a given reservoir, while keeping the temperature difference fixed.

On the Limit of Frequency of Electrochemical Oscillators and Its Relationship to Kinetic Parameters

The class of electrochemical oscillators characterized by a partially hidden negative differential resistance in an N-shaped current potential curve encompasses a myriad of experimental examples. We present a comprehensive methodological analysis of the oscillation frequency of this class of systems and discuss its dependence on electrical and kinetic parameters. The analysis is developed from a skeleton ordinary differential equation model, and an equation for the oscillation frequency is obtained. Simulations are carried out for a model system, namely, the nickel electrodissolution, and the numerical results are confirmed by experimental data on this system. In addition, the treatment is further applied to the electro-oxidation of ethylene glycol where unusually large oscillation frequencies have been reported. Despite the distinct chemistry underlying the oscillatory dynamics of these systems, a very good agreement between experiments and theoretical predictions is observed. The application of the developed theory is suggested as an important step for primary kinetic characterization.

Identification of Regions Involved in Substrate Binding and Dimer Stabilization within the Central Domains of Yeast Hsp40 Sis1

Protein folding, refolding and degradation are essential for cellular life and are regulated by protein homeostatic processes such those that involve the molecular chaperone DnaK/Hsp70 and its co-chaperone DnaJ. Hsp70 action is initiated when proteins from the DnaJ family bind an unfolded protein for delivery purposes. In eukaryotes, the DnaJ family can be divided into two main groups, Type I and Type II, represented by yeast cytosolic Ydj1 and Sis1, respectively. Although sharing some unique features both members of the DnaJ family, Ydj1 and Sis1 are structurally and functionally distinct as deemed by previous studies, including the observation that their central domains carry the structural and functional information even in switched chimeras. In this study, we combined several biophysical tools for evaluating the stability of Sis1 and mutants that had the central domains (named Gly/Met rich domain and C-terminal Domain I) deleted or switched to those of Ydj1 to gain insight into the role of these regions in the structure and function of Sis1. The mutants retained some functions similar to full length wild-type Sis1, however they were defective in others. We found that: 1) Sis1 unfolds in at least two steps as follows: folded dimer to partially folded monomer and then to an unfolded monomer. 2) The Gly/Met rich domain had intrinsically disordered characteristics and its deletion had no effect on the conformational stability of the protein. 3) The deletion of the C-terminal Domain I perturbed the stability of the dimer. 4) Exchanging the central domains perturbed the conformational stability of the protein. Altogether, our results suggest the existence of two similar subdomains in the C-terminal domain of DnaJ that could be important for stabilizing each other in order to maintain a folded substrate-binding site as well as the dimeric state of the protein.

Modeling the Emergence of Circadian Rhythms in a Clock Neuron Network

Circadian rhythms in pacemaker cells persist for weeks in constant darkness, while in other types of cells the molecular oscillations that underlie circadian rhythms damp rapidly under the same conditions. Although much progress has been made in understanding the biochemical and cellular basis of circadian rhythms, the mechanisms leading to damped or self-sustained oscillations remain largely unknown. There exist many mathematical models that reproduce the circadian rhythms in the case of a single cell of the Drosophila fly. However, not much is known about the mechanisms leading to coherent circadian oscillation in clock neuron networks. In this work we have implemented a model for a network of interacting clock neurons to describe the emergence (or damping) of circadian rhythms in Drosophila fly, in the absence of zeitgebers. Our model consists of an array of pacemakers that interact through the modulation of some parameters by a network feedback. The individual pacemakers are described by a well-known biochemical model for circadian oscillation, to which we have added degradation of PER protein by light and multiplicative noise. The network feedback is the PER protein level averaged over the whole network. In particular, we have investigated the effect of modulation of the parameters associated with (i) the control of net entrance of PER into the nucleus and (ii) the non-photic degradation of PER. Our results indicate that the modulation of PER entrance into the nucleus allows the synchronization of clock neurons, leading to coherent circadian oscillations under constant dark condition. On the other hand, the modulation of non-photic degradation cannot reset the phases of individual clocks subjected to intrinsic biochemical noise.

Mutual Information Rate and Bounds for It

The amount of information exchanged per unit of time between two nodes in a dynamical network or between two data sets is a powerful concept for analysing complex systems. This quantity, known as the mutual information rate (MIR), is calculated from the mutual information, which is rigorously defined only for random systems. Moreover, the definition of mutual information is based on probabilities of significant events. This work offers a simple alternative way to calculate the MIR in dynamical (deterministic) networks or between two time series (not fully deterministic), and to calculate its upper and lower bounds without having to calculate probabilities, but rather in terms of well known and well defined quantities in dynamical systems. As possible applications of our bounds, we study the relationship between synchronisation and the exchange of information in a system of two coupled maps and in experimental networks of coupled oscillators.

SYNCHRONIZATION OF CHAOS AND THE TRANSITION TO WAVE TURBULENCE

We investigated the transition to wave turbulence in a spatially extended three-wave interacting model, where a spatially homogeneous state undergoing chaotic dynamics undergoes spatial mode excitation. The transition to this weakly turbulent state can be regarded as the loss of synchronization of chaos of mode oscillators describing the spatial dynamics.

Period determination in the food-entrainable and methamphetamine-sensitive circadian oscillator(s)

Daily rhythmic processes are coordinated by circadian clocks, which are present in numerous central and peripheral tissues. In mammals, two circadian clocks, the food-entrainable oscillator (FEO) and methamphetamine-sensitive circadian oscillator (MASCO), are "black box" mysteries because their anatomical loci are unknown and their outputs are not expressed under normal physiological conditions. In the current study, the investigation of the timekeeping mechanisms of the FEO and MASCO in mice with disruption of all three paralogs of the canonical clock gene, Period, revealed unique and convergent findings. We found that both the MASCO and FEO in Per1(-/-)/Per2(-/-)/Per3(-/-) mice are circadian oscillators with unusually short (similar to 21 h) periods. These data demonstrate that the canonical Period genes are involved in period determination in the FEO and MASCO, and computational modeling supports the hypothesis that the FEO and MASCO use the same timekeeping mechanism or are the same circadian oscillator. Finally, these studies identify Per1(-/-)/Per2(-/-)/Per3(-/-) mice as a unique tool critical to the search for the elusive anatomical location(s) of the FEO and MASCO.

SINGULARLY PERTURBED DISCOUNTED MARKOV CONTROL PROCESSES IN A GENERAL STATE SPACE

This paper studies the asymptotic optimality of discrete-time Markov decision processes (MDPs) with general state space and action space and having weak and strong interactions. By using a similar approach as developed by Liu, Zhang, and Yin [Appl. Math. Optim., 44 (2001), pp. 105-129], the idea in this paper is to consider an MDP with general state and action spaces and to reduce the dimension of the state space by considering an averaged model. This formulation is often described by introducing a small parameter epsilon > 0 in the definition of the transition kernel, leading to a singularly perturbed Markov model with two time scales. Our objective is twofold. First it is shown that the value function of the control problem for the perturbed system converges to the value function of a limit averaged control problem as epsilon goes to zero. In the second part of the paper, it is proved that a feedback control policy for the original control problem defined by using an optimal feedback policy for the limit problem is asymptotically optimal. Our work extends existing results of the literature in the following two directions: the underlying MDP is defined on general state and action spaces and we do not impose strong conditions on the recurrence structure of the MDP such as Doeblin's condition.

The interaction of postural and voluntary strategies for stability in Parkinson's disease

de Lima-Pardini AC, Papegaaij S, Cohen RG, Teixeira LA, Smith BA, Horak FB. The interaction of postural and voluntary strategies for stability in Parkinson's disease. J Neurophysiol 108: 1244-1252, 2012. First published June 6, 2012; doi:10.1152/jn.00118.2012.-This study assessed the effects of stability constraints of a voluntary task on postural responses to an external perturbation in subjects with Parkinson's disease (PD) and healthy elderly participants. Eleven PD subjects and twelve control subjects were perturbed with backward surface translations while standing and performing two versions of a voluntary task: holding a tray with a cylinder placed with the flat side down [low constraint (LC)] or with the rolling, round side down [high constraint (HC)]. Participants performed alternating blocks of LC and HC trials. PD participants accomplished the voluntary task as well as control subjects, showing slower tray velocity in the HC condition compared with the LC condition. However, the latency of postural responses was longer in the HC condition only for control subjects. Control subjects presented different patterns of hip-shoulder coordination as a function of task constraint, whereas PD subjects had a relatively invariant pattern. Initiating the experiment with the HC task led to 1) decreased postural stability in PD subjects only and 2) reduced peak hip flexion in control subjects only. These results suggest that PD impairs the capacity to adapt postural responses to constraints imposed by a voluntary task.

Keeping your balance while balancing a cylinder: interaction between postural and voluntary goals

The present study investigated whether postural responses are influenced by the stability constraint of a voluntary, manual task. We also examined how task constraint and first experience (the condition with which the participants started the experiment) influence the kinematic strategies used to simultaneously accomplish a postural response and a voluntary task. Twelve healthy, older adults were perturbed during standing, while holding a tray with a cylinder placed with the flat side down (low constraint, LC) or with the rolling, round side down (high constraint, HC). Central set changed according to the task constraint, as shown by a higher magnitude of both the gastrocnemius and tibialis anterior muscle activation bursts in the HC than in the LC condition. This increase in muscle activation was not reflected, however, in changes in the center of pressure or center of mass displacement. Task constraint influenced the peak shoulder flexion for the voluntary tray task but not the peak hip flexion for the postural task. In contrast, first experience influenced the peak hip flexion but not the peak shoulder flexion. These results suggest an interaction between two separate control mechanisms for automatic postural responses and voluntary stabilization tasks.

Entanglement entropy in long-range harmonic oscillators

We study the Von Neumann and Renyi entanglement entropy of long-range harmonic oscillators (LRHO) by both theoretical and numerical means. We show that the entanglement entropy in massless harmonic oscillators increases logarithmically with the sub-system size as S - c(eff)/3 log l. Although the entanglement entropy of LRHO's shares some similarities with the entanglement entropy at conformal critical points we show that the Renyi entanglement entropy presents some deviations from the expected conformal behaviour. In the massive case we demonstrate that the behaviour of the entanglement entropy with respect to the correlation length is also logarithmic as the short-range case. Copyright (c) EPLA, 2012

Fourier´s law from a chain of coupled anharmonic oscillators under energy conserving shot noise.

We analyze the transport of heat along a chain of particles interacting through anharmonic potentials consisting of quartic terms in addition to harmonic quadratic terms and subject to heat reservoirs at its ends. Each particle is also subject to an impulsive shot noise with exponentially distributed waiting times whose effect is to change the sign of its velocity, thus conserving the energy of the chain. We show that the introduction of this energy conserving stochastic noise leads to Fourier's law. That is for large system size L the heat current J behaves as J ‘approximately’ 1/L, which amounts to say that the conductivity k is constant. The conductivity is related to the current by J = kΔT/L, where ΔT is the difference in the temperatures of the reservoirs. The behavior of heat conductivity k for small intensities¸ of the shot noise and large system sizes L are obtained by assuming a scaling behavior of the type k = ‘L POT a Psi’(L’lambda POT a/b’) where a and b are scaling exponents. For the pure harmonic case a = b = 1, characterizing a ballistic conduction of heat when the shot noise is absent. For the anharmonic case we found values for the exponents a and b smaller then 1 and thus consistent with a superdiffusive conduction of heat without the shot noise. We also show that the heat conductivity is not constant but is an increasing function of temperature.

Pullback exponential attractors for evolution processes in Banach spaces: properties and applications

This article is a continuation of our previous work [5], where we formulated general existence theorems for pullback exponential attractors for asymptotically compact evolution processes in Banach spaces and discussed its implications in the autonomous case. We now study properties of the attractors and use our theoretical results to prove the existence of pullback exponential attractors in two examples, where previous results do not apply.