A Fast-Slow Analysis of the Dynamics of REM Sleep

Waking and sleep states are regulated by the coordinated activity of a number of neuronal population in the brainstem and hypothalamus whose synaptic interactions compose a sleep-wake regulatory network. Physiologically based mathematical models of the sleep-wake regulatory network contain mechanisms operating on multiple time scales including relatively fast synaptic-based interations between neuronal populations, and much slower homeostatic and circadian processes that modulate sleep-wake temporal patterning. In this study, we exploit the naturally arising slow time scale of the homeostatic sleep drive in a reduced sleep-wake regulatory network model to utilize fast-slow analysis to investigate the dynamics of rapid eye movement (REM) sleep regulation. The network model consists of a reduced number of wake-, non-REM (NREM) sleep-, and REM sleep-promoting neuronal populations with the synaptic interactions reflecting the mutually inhibitory flip-flop conceptual model for sleep-wake regulation and the reciprocal interaction model for REM sleep regulation. Network dynamics regularly alternate between wake and sleep states as goverend by the slow homeostatic sleep drive. By varying a parameter associated with the activation of the REM-promoting population, we cause REM dynamics during sleep episodes to vary from supression to single activations to regular REM-NREM cycling, corresponding to changes in REM patterning induced by circadian modulation and observed in different mammalian species. We also utilize fast-slow analysis to explain complex effects on sleep-wake patterning of simulated experiments in which agonists and antagonists of different neurotransmitters are microinjected into specific neuronal populations participating in the sleep-wake regulatory network

Convolution of multifractals and the local magnetization in a random field Ising chain

The local magnetization in the one-dimensional random-field Ising model is essentially the sum of two effective fields with multifractal probability measure. The probability measure of the local magnetization is thus the convolution of two multifractals. In this paper we prove relations between the multifractal properties of two measures and the multifractal properties of their convolution. The pointwise dimension at the boundary of the support of the convolution is the sum of the pointwise dimensions at the boundary of the support of the convoluted measures and the generalized box dimensions of the convolution are bounded from above by the sum of the generalized box dimensions of the convoluted measures. The generalized box dimensions of the convolution of Cantor sets with weights can be calculated analytically for certain parameter ranges and illustrate effects we also encounter in the case of the measure of the local magnetization. Returning to the study of this measure we apply the general inequalities and present numerical approximations of the D_q-spectrum. For the first time we are able to obtain results on multifractal properties of a physical quantity in the one-dimensional random-field Ising model which in principle could be measured experimentally. The numerically generated probability densities for the local magnetization show impressively the gradual transition from a monomodal to a bimodal distribution for growing random field strength h.Comment: An error in figure 1 was corrected, small additions were made to the introduction and the conclusions, some typos were corrected, 24 pages, LaTeX2e, 9 figure

The randomly driven Ising ferromagnet, Part I: General formalism and mean field theory

We consider the behavior of an Ising ferromagnet obeying the Glauber dynamics under the influence of a fast switching, random external field. After introducing a general formalism for describing such systems, we consider here the mean-field theory. A novel type of first order phase transition related to spontaneous symmetry breaking and dynamic freezing is found. The non-equilibrium stationary state has a complex structure, which changes as a function of parameters from a singular-continuous distribution with Euclidean or fractal support to an absolutely continuous one.Comment: 12 pages REVTeX/LaTeX format, 12 eps/ps figures. Submitted to Journal of Physics

Orbits and phase transitions in the multifractal spectrum

We consider the one dimensional classical Ising model in a symmetric dichotomous random field. The problem is reduced to a random iterated function system for an effective field. The D_q-spectrum of the invariant measure of this effective field exhibits a sharp drop of all D_q with q < 0 at some critical strength of the random field. We introduce the concept of orbits which naturally group the points of the support of the invariant measure. We then show that the pointwise dimension at all points of an orbit has the same value and calculate it for a class of periodic orbits and their so-called offshoots as well as for generic orbits in the non-overlapping case. The sharp drop in the D_q-spectrum is analytically explained by a drastic change of the scaling properties of the measure near the points of a certain periodic orbit at a critical strength of the random field which is explicitly given. A similar drastic change near the points of a special family of periodic orbits explains a second, hitherto unnoticed transition in the D_q-spectrum. As it turns out, a decisive role in this mechanism is played by a specific offshoot. We furthermore give rigorous upper and/or lower bounds on all D_q in a wide parameter range. In most cases the numerically obtained D_q coincide with either the upper or the lower bound. The results in this paper are relevant for the understanding of random iterated function systems in the case of moderate overlap in which periodic orbits with weak singularity can play a decisive role.Comment: The article has been completely rewritten; the title has changed; a section about the typical pointwise dimension as well as several references and remarks about more general systems have been added; to appear in J. Phys. A; 25 pages, 11 figures, LaTeX2

Einfluss der Aufschlusstemperatur auf die morphologischen Eigenschaften von TMP aus Kiefernholz

Chips from pine wood were subjected to thermomechanical pulping (TMP) at 140 and 180 degrees C for 5 minutes, whereas the cooked chips were defibrated using a single disk pressurized refiner at the same temperatures (140 and 180 degrees C). The fibres were tested for some of their morphological properties including fibre length, fibre width, cell-wall thickness. Moreover, the fine fibre fraction (zero fibres) and the content of splinters were also estimated. The results reveal, that increasing the temperature during thermomechanical pulping decreases the fibre length, the cell width and the fibre wall thickness. It also increases the amount of fine fibres and increases the curl factor

Phase diagram of the random field Ising model on the Bethe lattice

The phase diagram of the random field Ising model on the Bethe lattice with a symmetric dichotomous random field is closely investigated with respect to the transition between the ferromagnetic and paramagnetic regime. Refining arguments of Bleher, Ruiz and Zagrebnov [J. Stat. Phys. 93, 33 (1998)] an exact upper bound for the existence of a unique paramagnetic phase is found which considerably improves the earlier results. Several numerical estimates of transition lines between a ferromagnetic and a paramagnetic regime are presented. The obtained results do not coincide with a lower bound for the onset of ferromagnetism proposed by Bruinsma [Phys. Rev. B 30, 289 (1984)]. If the latter one proves correct this would hint to a region of coexistence of stable ferromagnetic phases and a stable paramagnetic phase.Comment: Article has been condensed and reorganized; Figs 3,5,6 merged; Fig 4 omitted; Some discussion added at end of Sec. III; 9 pages, 5 figs, RevTeX4, AMSTe

MeltMigrator : a MATLAB-based software for modeling three-dimensional melt migration and crustal thickness variations at mid-ocean ridges following a rules-based approach

Author Posting. © American Geophysical Union, 2017. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry, Geophysics, Geosystems 18 (2017): 445–456, doi:10.1002/2016GC006686.MeltMigrator is a MATLAB®-based melt migration software developed to process three-dimensional mantle temperature and velocity data from user-supplied numerical models of mid-ocean ridges, calculate melt production and melt migration trajectories in the mantle, estimate melt flux along plate boundaries, and predict crustal thickness distribution on the seafloor. MeltMigrator is also capable of calculating compositional evolution depending on the choice of petrologic melting model. Programmed in modules, MeltMigrator is highly customizable and can be expanded to a wide range of applications. We have applied it to complex mid-ocean ridge model settings, including transform faults, oblique segments, ridge migration, asymmetrical spreading, background mantle flow, and ridge-plume interaction. In this technical report, we include an example application to a segmented mid-ocean ridge. MeltMigrator is available as a supplement to this paper, and it is also available from GitHub and the University of Maryland Geodynamics Group website.National Science Foundation Grant Number: OCE-0937277 and OCE-14582012017-07-2

Randomly Evolving Idiotypic Networks: Structural Properties and Architecture

We consider a minimalistic dynamic model of the idiotypic network of B-lymphocytes. A network node represents a population of B-lymphocytes of the same specificity (idiotype), which is encoded by a bitstring. The links of the network connect nodes with complementary and nearly complementary bitstrings, allowing for a few mismatches. A node is occupied if a lymphocyte clone of the corresponding idiotype exists, otherwise it is empty. There is a continuous influx of new B-lymphocytes of random idiotype from the bone marrow. B-lymphocytes are stimulated by cross-linking their receptors with complementary structures. If there are too many complementary structures, steric hindrance prevents cross-linking. Stimulated cells proliferate and secrete antibodies of the same idiotype as their receptors, unstimulated lymphocytes die. Depending on few parameters, the autonomous system evolves randomly towards patterns of highly organized architecture, where the nodes can be classified into groups according to their statistical properties. We observe and describe analytically the building principles of these patterns, which allow to calculate number and size of the node groups and the number of links between them. The architecture of all patterns observed so far in simulations can be explained this way. A tool for real-time pattern identification is proposed.Comment: 19 pages, 15 figures, 4 table