Supplemental Figures A Nonlinear Dynamical Theory of Cell Injury
|
|
- Jemimah Mathews
- 5 years ago
- Views:
Transcription
1 upplementary Figures, A Nonlinear ynamical Theory of Cell Injury - upplemental Figures A Nonlinear ynamical Theory of Cell Injury onald J. egracia, Zhi-Feng Huang, ui Huang A. f ( X ) = a Θ n n Θ + X n B..0 f ( X ) n = n = 2 n = 4 n = Θ = X X upplemental Figure. Characteristics of the Hill function. (A) Equation of the Hill function. a is an amplitude parameter that can be used as a normalizing factor. The ratio Θ n /(Θ n + X n ) ranges from 0 to, where n is the Hill coefficient. (B) Plots of a Hill function at various values of n. At low n, the function approximates a linear relationship. As n increases, the function assumes its characteristic sigmoidal form. At large n, the Hill function approximates a square wave, clearly illustrating an on/off effect with a threshold Θ. At Θ = X, f(x) = 0.5a; the threshold Θ then indicates a 50% inhibition effect. Conceptually, the use of the Hill function in the nonlinear dynamic model of this paper provides a threshold value of (i.e., Θ ) above which total damage will effectively inhibit total stress responses. And vice versa, there will be a threshold value of (i.e., Θ ) above which stress responses will be able to effectively inhibit damage. A key notion of Eq. (4) (main text) is that Θ and Θ change with injury magnitude, I. At low values of injury, the stress responses will have high Θ and be effective in inhibiting damage. As I increases so does Θ, making it progressively easier for damage to inhibit stress responses. In addition, Θ decreases with increasing I, making stress responses progressively less effective as injury magnitudes increase.
2 upplementary Files, A Nonlinear ynamical Theory of Cell Injury - 2 upplemental Figure 2. Eq. (4) in the main text is a normalized form (called rescaled in the physics literature) such that at each I, the maximum values of and are or 00%. We seek to modify Eq. (4) to output absolute values (referred to as unscaled ) of and. We therefore make one additional modeling assumption over that of Eq. (4): the v and v parameters are also defined as functions of I along with the Θ and Θ parameters [Eq. (2)], giving the following OE system: where the v and Θ parameters are the I-dependant functions: n d ' Θ' = v ' ' n n k dt Θ ' + ' n d ' Θ' = v' k' n n dt Θ ' + ' () v' = v Ie Θ ' = c I e + λ3i 2 + λi v' = v Ie Θ ' = c I e λ4i 2 λ2i (2) The transformation between the normalized/rescaled form Eq. (4) and the unscaled form Eq. () is: Ie Ie + λ3i λ4i ' = ( ), ' = ( ) which implies that λ = λ + λ 4 and λ = λ 2 + λ 3 in Eq. (4), and v, v, c, c, k, k, n, and n are unchanged from those in Eq. (4). Figure 2 shows the output of unscaled Eq. () (panels A, C) as compared to scaled Eq.(4) (panels B, ). (panels A and B) and * (panels C and ) shown in red. (panels A and B) and * (panels C and ) shown in green. There is no qualitative difference in the dynamics of Eq. (4) vs. Eq. () when λ 3 = λ 4 = 0. (3)
3 upplemental Figure 3. A brief primer on graphically solving an ordinary differential equation (OE). (A) Because in general nonlinear OEs cannot be solved exactly by integration, a graphical method invented by the French mathematician Henri Poincaré is now widely used to study OEs. This method is in fact relatively simple conceptually, as illustrated here. In short, a point in the plane, along with the OE parameters, serves as input to the OE (step ). The OE expression is then solved with the input numbers, and the OE outputs two magnitudes (d/dt and d/dt; the numerical values are listed as v and v 2 ). As the ratios d/dt and d/dt indicate, these are slopes and are vector quantities in the (, ) plane. The resultant of v and v 2 (labeled v 2 ) is then calculated (step 2) and the vector v 2 plotted at the point that originally serves as the input to the OE (step 3). One repeats this calculation, looping through the points in the plane, and plotting the A. Graphical construction of a vector field 0 0 (, ) Iλ n d ( cie ) = v I n n k = v λ dt ( cie ) + Iλ n d ( cie ) = v k I = v λ n n 2 dt ( c Ie ) + (2.) [ ij, ij, I, c, λ, c, λ, v, v, k, k, n, n ] = [0.2, 0.2, 8, 0., 0., 2.5, 0.9,,,,, 4, 4 ] vector output of the OE at each point to generate a vector field (also called a slope field). (B) A solution to the OE is any expression that satisfies the differential equation. For linear OEs, the exact solutions can be obtained by integration, but not for most nonlinear OEs. Therefore the vector field provides a graphical means to see all solutions to the OE. The solutions are known as trajectories. A trajectory is simply a flow path along the vector field. There are several well-known methods for numerically approximating trajectories on the vector field (e.g. the Runge-Kutta method). A specific trajectory is determined by the initial conditions, or the starting point of the trajectory. Illustrated are several trajectories from different initial conditions, where the initial conditions are depicted by stars. In general, the space of all trajectories is called a state space. ince our OE is two dimensional, the state space is referred to as the phase plane. In this example, it can be seen by visual inspection that every trajectory ends at the same point on the phase plane. This point is indicated by the yellow circle with red outline. ince all trajectories are attracted to this point, it is called an attractor point and is designated (*, *). More specifically, the attractor is the fixed point, or equilibrium point. Here the system no longer changes in time; the system is fixed at the attractor state. Because there is only one attractor on this example of phase plane, this phase plane is called monostable. A specific trajectory is mathematically equivalent to a matching pair of and time course. The trajectory itself is a curve, where each point on the curve indicates a pair (, ) which are the values of and that occur at the same point on a time course. This is illustrated in the figure (B) by showing three (, ) pairs (in red, green, and blue), and their corresponding locations on the time course. Thus, a trajectory is a more general way to view a time course. Finally, the fixed points or equilibrium points of the phase plane are, in a sense, the steady solutions to the OE because all trajectories in a specific phase plane go to these fixed points. The parameter values used for this example are listed in the figure. (.) v2 = v +v2 B. Phase Planes, trajectories, initial conditions and time courses [I, c, λ, c, λ, v, v, k, k, n, n ] = [3.5, 0., 0., 35, 0.9,,,,, 4, 4 ] Monostable phase plane with attractor (*, *) at (0.06, ) ( 0, 0 ) = (0, 0) and (3.) time, t 8
4 upplementary Files, A Nonlinear ynamical Theory of Cell Injury - 4 upplemental Figure 4: A bistable phase plane. Prior to the widespread availability of computers, it was generally unappreciated how complex OE solutions were: the focus on linear OEs gave a misleadingly simple impression of OE behavior. However, when considering the larger class of nonlinear OEs, one cannot, in general, know a priori the complicated flow patterns of the vector field. An element of surprise can occur, called emergence, when the vector field is plotted and reveals unanticipated features in the state space. Our nonlinear model embodied in Eq. (4) is relatively simple, yet under some parameter conditions, a vector field more complex than the monostable type can emerge, as shown here. When one visually inspects the vector field generated by the above set of parameters, and traces out trajectories, it is seen that there are two possible fixed point attractors to which the trajectories can converge. One of these is in the upper left (green) at (0.03, ), and the other in the lower right (red) at (0.99, 0.06) of the - unit plane. This type of phase plane is called bistable because there are two possible stable fixed points. In addition there is a third point at (0.57, 0.38) (blue) to which the trajectories appear to be drawn toward; but as they get closer to this point, they are repelled from it and towards one of the attractors. This point is thus called a repeller. When there is more than one fixed point, the attractors and repeller fixed points always appear together (trogatz, 994). The repeller point divides the phase plane into those trajectories that approach the attractor at (0.03, ) and those that approach (0.99, 0.06). The yellow dotted line indicates this division and is called the separatrix. Thus, all trajectories to the left of the separatrix form the basin of attraction for the attractor at (0.03, ). Trajectories to the right of the separatrix form the basin of attraction for the attractor at (0.99, 0.06). The emergence of bistable phase planes in the solutions to Eq. (4) constitutes one of the most important features of our nonlinear cell injury model. Because our theory associates attractors of * > * with survival, and those of * > * with death, bistable phase planes represent clear-cut avenues for therapeutic treatments based on the injury dynamics. As discussed in the article text, an injured system within the basin of attraction of the death attractor can be diverted, at least in principle, via some manipulation of the system, to the survival attractor. Initial conditions for the two plotted trajectories are marked by stars and their coordinates listed on the plot. The input parameters for this phase plane are also listed in the figure.
5 upplemental Figure 5: Often an OE system contains various parameters and it is of interest to study the behavior of the OE system as the parameters change. Parameters are not variables. A simple example explains this distinction. Condier the familiar equation of a line: y=mx+b. For this equation, x is the variable, and m and b are parameters. ifferent lines are obained by changing the m and b parameters, but all of them are lines and all are functions of the variable, x. An analogous logic holds for the nonlinear dynamic model of cell injury. Its parameters give different injury systems analogous to how the m and b parameters generate different lines. It is of interest to understand how the OE system behaves when the parameters change; this study is called variously bifurcation analysis, qualitative analysis or fixed point analysis and is illustrated here. Again, the basic notion is relatively simple. In general, the fixed points change as a parameter changes. (A) This is illustrated by showing several vector fields resulting from changing the parameter I (leaving all others constant). As can be seen, the fixed points are different amongst the vector fields. In particular, the number of fixed points changes, in this instance, from one to three. When the number of fixed points changes, a bifurcation occurs. In this example, the system bifurcates twice: it goes from one to three fixed points, and then goes from three to one fixed point. The presence of bifurcations is an example of emergence in nonlinear dynamics (see trogatz, 994, for further information). The numerical values of the fixed points (*, *) are listed below each respective vector field. (B) To characterize a system as the parameters change, one constructs a bifurcation diagram by simply plotting the fixed points of the OE vs. the parameter of interest. Because our system is two dimensional, every fixed point is a (*, *) pair. The system bifurcates at the same parameter values for * and * so only one needs to be studied. Here we plot * vs. I. When the system bifurcates, all three values of *, corresponding to each fixed point, are plotted at the given value of I (arrows). The repellers are here colored blue to distinguish them from the attractors. In bistable regions of OE solutions, there are two possible outcomes. In our model one outcome is interpreted as a survival phenotype (green), the other as a death phenotype (red). The * value from each fixed point is similarly color-coded in the bifurcation diagram of (B). The type of bifurcation diagram shown, which results from varying only the single parameter, I, is the simplest way to analyze Eq. (4). Because this OE has 2 variables and parameters, a full bifurcation analysis would vary all 3 quantities within the model constraints (e.g. c > c and λ > 0 and λ > 0) resulting in 3-dimensional bifurcation diagrams. Parameter values are listed on the figure.
6 upplementary Files, A Nonlinear ynamical Theory of Cell Injury - 6 upplemental Figure 6: The four types of injury courses are cross-sections of surfaces in 3 parameter space and interconvert via the surface topology. (A) Illustration of type A bistable injury course dynamics converting to monostable injury course dynamics. Yellow and orange dashed lines show characteristic 2 bifurcation diagram profiles for * and *, respectively, for type A bistable injury course dynamics. (B) Illustration of type B bistable injury course dynamics converting to monostable injury course dynamics. Yellow and orange dashed lines show characteristic 2 bifurcation diagram profiles for * and *, respectively, for type B bistable injury course dynamics. For 3 plots of *, attractors are indicated as red and repellers blue; for 3 plots of *, attractors are green and repellers are blue. Insets are rotations of surfaces. ecay time vs. I plots shown under 3 plots correspond to the
7 cross sections at c values indicated in the 3 panels; decay times are in red, decay times in green, and dashed magenta lines are for I = I X. Parameters (c, λ, λ, c ) used to generate plots in A and B were [0., 0., 0.9, 2.5 < c < 50] and [0.075, 0.0, 0.0, < c < 0.57], respectively. Figure 6 illustrates the bifurcation diagrams of * and * along families of injury courses resulting from varying c and I. Figure 6A shows how type A bistable injury courses transform into monostable ones. The repellers form a curved wall within the surface that disappears at a discreet value of c : the subsequent injury courses are then monostable. In the (t R, t ) vs. I plots below the 3 figures, at lower c the decay time plots are clearly of the type A bistable variety; however, with increasing c, they gradually morph into the monostable decay time plots with the characteristic spike around I X. Figure 6B shows a family of injury courses that transform from the type B bistable into the monostable dynamical pattern. A different 2 surface emerges with a repeller wall forming an almost conical shape. In the decay time plots, as c increases, the post-i X decay hump, characteristic of the type B bistable dynamical pattern, merges with and transforms into the spike characteristic of the monostable dynamical pattern. These examples illustrate that the dynamics of a specific injury course can be obtained from the cross-sections within the higher dimension parameter space of the solutions to Eq. (4). They also illustrate that the model can output many variations of the four basic injury course dynamic patterns. In fact, multi-dimensional bifurcation analysis suggests an intimate relationship between different injury systems and their dynamics. That is, overtly different forms of injury could be related to each other based on the dynamical patterns of their injury courses. At present, different injury modalities, for example, ischemia, trauma, shock, seizures, etc., are treated as completely different. While there is increasing recognition of common damage mechanisms (e.g. oxidative damage, induction of stress responses such as the heat shock response), our model shows there are natural relationships amongst the dynamics of injury courses.
8 upplementary Files, A Nonlinear ynamical Theory of Cell Injury - 8 upplemental Figure 7. Manipulation of initial conditions to cause a death outcome at I < I X. As discussed in the main text, altering initial conditions can cause the system to flip state in a bistable phase plane. (A) Plots of Θ and Θ vs. I; I X = 3.3. (B) Phase plane at I = 3. A number of survival (green) and death (red) trajectories are plotted and reveal that 0 = is the cutoff below which the system survives and above which it dies. Also indicated are the survival attractor (green), death attractor (red), and repeller (blue). (C) The injury course of the system. ashed magenta line is I = I X, dotted blue line is I = 3, which corresponds to the level of the phase plane in B. This is a type B bistable injury course. The parameter values used for this example are listed in the figure. In this bistable phase plane at I < I X, from initial conditions (0, 0) the system approaches the survival attractor (green circle) and survives. However, if 0 > (interpreted as the system being damaged to 6.9% of its maximum), the system will divert to the death attractor and die from what would be a sublethal insult when induced from the uninjured [e.g. ( 0, 0 ) = (0, 0)] state. This example is reminiscent of the worsened stroke outcome experienced by diabetic patients. It is well-known that diabetics suffer worse stroke outcome than non-diabetics. Research in this area seeks to elucidate specific diabetes-induced damage pathways that worsen stroke outcome. While it has been known that increased cardio-vascular disease in diabetics increases stroke risk, this does not fully explain the worsened outcome. Our model suggests that the worsened outcome is due to the underlying injury dynamics of the system, and that diabetes serves as a role of setting the initial conditions of the stroke injury such that 0 > 0.
Networks in systems biology
Networks in systems biology Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4500, Spring 2017 M. Macauley (Clemson) Networks in systems
More informationPrinciples of Synthetic Biology: Midterm Exam
Principles of Synthetic Biology: Midterm Exam October 28, 2010 1 Conceptual Simple Circuits 1.1 Consider the plots in figure 1. Identify all critical points with an x. Put a circle around the x for each
More information= F ( x; µ) (1) where x is a 2-dimensional vector, µ is a parameter, and F :
1 Bifurcations Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 A bifurcation is a qualitative change
More informationFundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad
Fundamentals of Dynamical Systems / Discrete-Time Models Dr. Dylan McNamara people.uncw.edu/ mcnamarad Dynamical systems theory Considers how systems autonomously change along time Ranges from Newtonian
More information... it may happen that small differences in the initial conditions produce very great ones in the final phenomena. Henri Poincaré
Chapter 2 Dynamical Systems... it may happen that small differences in the initial conditions produce very great ones in the final phenomena. Henri Poincaré One of the exciting new fields to arise out
More informationThe Nonlinear Pendulum
The Nonlinear Pendulum - Pádraig Ó Conbhuí - 08531749 TP Monday 1. Abstract This experiment was performed to examine the effects that linearizing equations has on the accuracy of results and to find ways
More informationDynamical Systems. 1.0 Ordinary Differential Equations. 2.0 Dynamical Systems
. Ordinary Differential Equations. An ordinary differential equation (ODE, for short) is an equation involves a dependent variable and its derivatives with respect to an independent variable. When we speak
More informationProject 1 Modeling of Epidemics
532 Chapter 7 Nonlinear Differential Equations and tability ection 7.5 Nonlinear systems, unlike linear systems, sometimes have periodic solutions, or limit cycles, that attract other nearby solutions.
More informationThe Big Picture. Discuss Examples of unpredictability. Odds, Stanisław Lem, The New Yorker (1974) Chaos, Scientific American (1986)
The Big Picture Discuss Examples of unpredictability Odds, Stanisław Lem, The New Yorker (1974) Chaos, Scientific American (1986) Lecture 2: Natural Computation & Self-Organization, Physics 256A (Winter
More informationPH 120 Project # 2: Pendulum and chaos
PH 120 Project # 2: Pendulum and chaos Due: Friday, January 16, 2004 In PH109, you studied a simple pendulum, which is an effectively massless rod of length l that is fixed at one end with a small mass
More informationCOMPOSITE MATERIAL WITH NEGATIVE STIFFNESS INCLUSION FOR VIBRATION DAMPING: THE EFFECT OF A NONLINEAR BISTABLE ELEMENT
11 th International Conference on Vibration Problems Z. Dimitrovová et.al. (eds.) Lisbon, Portugal, 9 12 September 2013 COMPOSITE MATERIAL WITH NEGATIVE STIFFNESS INCLUSION FOR VIBRATION DAMPING: THE EFFECT
More informationThe construction and use of a phase diagram to investigate the properties of a dynamic system
File: Phase.doc/pdf The construction and use of a phase diagram to investigate the properties of a dynamic system 1. Introduction Many economic or biological-economic models can be represented as dynamic
More informationSystems of Linear ODEs
P a g e 1 Systems of Linear ODEs Systems of ordinary differential equations can be solved in much the same way as discrete dynamical systems if the differential equations are linear. We will focus here
More informationLinearization of Differential Equation Models
Linearization of Differential Equation Models 1 Motivation We cannot solve most nonlinear models, so we often instead try to get an overall feel for the way the model behaves: we sometimes talk about looking
More informationNonlinear dynamics & chaos BECS
Nonlinear dynamics & chaos BECS-114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes
More informationDynamical Systems: Lecture 1 Naima Hammoud
Dynamical Systems: Lecture 1 Naima Hammoud Feb 21, 2017 What is dynamics? Dynamics is the study of systems that evolve in time What is dynamics? Dynamics is the study of systems that evolve in time a system
More informationTransitioning to Chaos in a Simple Mechanical Oscillator
Transitioning to Chaos in a Simple Mechanical Oscillator Hwan Bae Physics Department, The College of Wooster, Wooster, Ohio 69, USA (Dated: May 9, 8) We vary the magnetic damping, driver frequency, and
More informationKey ideas about dynamic models
GSTDMB 202: DYNAMICAL MODELLING FOR BIOLOGY AND MEDICINE Lecture 2 Introduction to modelling with differential equations Markus Owen ey ideas about dynamic models Dynamic models describe the change in
More informationIntroduction to Dynamical Systems Basic Concepts of Dynamics
Introduction to Dynamical Systems Basic Concepts of Dynamics A dynamical system: Has a notion of state, which contains all the information upon which the dynamical system acts. A simple set of deterministic
More information16 Period doubling route to chaos
16 Period doubling route to chaos We now study the routes or scenarios towards chaos. We ask: How does the transition from periodic to strange attractor occur? The question is analogous to the study of
More informationSupporting Information
Supporting Information Visualizing the Effect of Partial Oxide Formation on Single Silver Nanoparticle Electrodissolution Vignesh Sundaresan, Joseph W. Monaghan, and Katherine A. Willets* Department of
More informationNumerical Algorithms as Dynamical Systems
A Study on Numerical Algorithms as Dynamical Systems Moody Chu North Carolina State University What This Study Is About? To recast many numerical algorithms as special dynamical systems, whence to derive
More informationLecture 3. Dynamical Systems in Continuous Time
Lecture 3. Dynamical Systems in Continuous Time University of British Columbia, Vancouver Yue-Xian Li November 2, 2017 1 3.1 Exponential growth and decay A Population With Generation Overlap Consider a
More informationNonlinear Autonomous Systems of Differential
Chapter 4 Nonlinear Autonomous Systems of Differential Equations 4.0 The Phase Plane: Linear Systems 4.0.1 Introduction Consider a system of the form x = A(x), (4.0.1) where A is independent of t. Such
More informationElectric Field and Electric Potential
Electric Field and Electric Potential INTRODUCTION Physicists use the concept of a field 1 to explain the interaction of particles or bodies through space, i.e., the action-at-a-distance 2 force between
More information3 Mathematical modeling of the torsional dynamics of a drill string
3 Mathematical modeling of the torsional dynamics of a drill string 3.1 Introduction Many works about torsional vibrations on drilling systems [1, 12, 18, 24, 41] have been published using different numerical
More informationANALYSIS OF BIOLOGICAL NETWORKS USING HYBRID SYSTEMS THEORY. Nael H. El-Farra, Adiwinata Gani & Panagiotis D. Christofides
ANALYSIS OF BIOLOGICAL NETWORKS USING HYBRID SYSTEMS THEORY Nael H El-Farra, Adiwinata Gani & Panagiotis D Christofides Department of Chemical Engineering University of California, Los Angeles 2003 AIChE
More informationExample of a Blue Sky Catastrophe
PUB:[SXG.TEMP]TRANS2913EL.PS 16-OCT-2001 11:08:53.21 SXG Page: 99 (1) Amer. Math. Soc. Transl. (2) Vol. 200, 2000 Example of a Blue Sky Catastrophe Nikolaĭ Gavrilov and Andrey Shilnikov To the memory of
More informationA Producer-Consumer Model With Stoichiometry
A Producer-Consumer Model With Stoichiometry Plan B project toward the completion of the Master of Science degree in Mathematics at University of Minnesota Duluth Respectfully submitted by Laura Joan Zimmermann
More informationDynamical Embodiments of Computation in Cognitive Processes James P. Crutcheld Physics Department, University of California, Berkeley, CA a
Dynamical Embodiments of Computation in Cognitive Processes James P. Crutcheld Physics Department, University of California, Berkeley, CA 94720-7300 and Santa Fe Institute, 1399 Hyde Park Road, Santa Fe,
More informationDynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0
More informationChapter 24 BIFURCATIONS
Chapter 24 BIFURCATIONS Abstract Keywords: Phase Portrait Fixed Point Saddle-Node Bifurcation Diagram Codimension-1 Hysteresis Hopf Bifurcation SNIC Page 1 24.1 Introduction In linear systems, responses
More informationARTIFICIAL NEURAL NETWORK PART I HANIEH BORHANAZAD
ARTIFICIAL NEURAL NETWORK PART I HANIEH BORHANAZAD WHAT IS A NEURAL NETWORK? The simplest definition of a neural network, more properly referred to as an 'artificial' neural network (ANN), is provided
More informationNeural Excitability in a Subcritical Hopf Oscillator with a Nonlinear Feedback
Neural Excitability in a Subcritical Hopf Oscillator with a Nonlinear Feedback Gautam C Sethia and Abhijit Sen Institute for Plasma Research, Bhat, Gandhinagar 382 428, INDIA Motivation Neural Excitability
More informationMore Details Fixed point of mapping is point that maps into itself, i.e., x n+1 = x n.
More Details Fixed point of mapping is point that maps into itself, i.e., x n+1 = x n. If there are points which, after many iterations of map then fixed point called an attractor. fixed point, If λ
More informationNumerical method for approximating the solution of an IVP. Euler Algorithm (the simplest approximation method)
Section 2.7 Euler s Method (Computer Approximation) Key Terms/ Ideas: Numerical method for approximating the solution of an IVP Linear Approximation; Tangent Line Euler Algorithm (the simplest approximation
More informationSimplest Chaotic Flows with Involutional Symmetries
International Journal of Bifurcation and Chaos, Vol. 24, No. 1 (2014) 1450009 (9 pages) c World Scientific Publishing Company DOI: 10.1142/S0218127414500096 Simplest Chaotic Flows with Involutional Symmetries
More informationVectors a vector is a quantity that has both a magnitude (size) and a direction
Vectors In physics, a vector is a quantity that has both a magnitude (size) and a direction. Familiar examples of vectors include velocity, force, and electric field. For any applications beyond one dimension,
More informationSupplementary Figures
Supplementary Figures a x 1 2 1.5 1 a 1 = 1.5 a 1 = 1.0 0.5 a 1 = 1.0 b 0 0 20 40 60 80 100 120 140 160 2 t 1.5 x 2 1 0.5 a 1 = 1.0 a 1 = 1.5 a 1 = 1.0 0 0 20 40 60 80 100 120 140 160 t Supplementary Figure
More informationSTUDY OF SYNCHRONIZED MOTIONS IN A ONE-DIMENSIONAL ARRAY OF COUPLED CHAOTIC CIRCUITS
International Journal of Bifurcation and Chaos, Vol 9, No 11 (1999) 19 4 c World Scientific Publishing Company STUDY OF SYNCHRONIZED MOTIONS IN A ONE-DIMENSIONAL ARRAY OF COUPLED CHAOTIC CIRCUITS ZBIGNIEW
More informationMagnetic Fields. Goals. Introduction. Mapping magnetic fields with iron filings
Lab 7. Magnetic Fields Goals To visualize the magnetic fields produced by several different configurations of simple bar magnets using iron filings. To use small magnetic compasses to trace out the magnetic
More informationComputer Problems for Taylor Series and Series Convergence
Computer Problems for Taylor Series and Series Convergence The two problems below are a set; the first should be done without a computer and the second is a computer-based follow up. 1. The drawing below
More informationDaba Meshesha Gusu and O.Chandra Sekhara Reddy 1
International Journal of Basic and Applied Sciences Vol. 4. No. 1 2015. Pp.22-27 Copyright by CRDEEP. All Rights Reserved. Full Length Research Paper Solutions of Non Linear Ordinary Differential Equations
More informationSjoerd Verduyn Lunel (Utrecht University)
Universiteit Utrecht, February 9, 016 Wasserstein distances in the analysis of time series and dynamical systems Sjoerd Verduyn Lunel (Utrecht University) S.M.VerduynLunel@uu.nl Abstract A new approach
More informationStability of Dynamical systems
Stability of Dynamical systems Stability Isolated equilibria Classification of Isolated Equilibria Attractor and Repeller Almost linear systems Jacobian Matrix Stability Consider an autonomous system u
More informationAPPPHYS217 Tuesday 25 May 2010
APPPHYS7 Tuesday 5 May Our aim today is to take a brief tour of some topics in nonlinear dynamics. Some good references include: [Perko] Lawrence Perko Differential Equations and Dynamical Systems (Springer-Verlag
More informationJohns Hopkins Math Tournament Proof Round: Automata
Johns Hopkins Math Tournament 2018 Proof Round: Automata February 9, 2019 Problem Points Score 1 10 2 5 3 10 4 20 5 20 6 15 7 20 Total 100 Instructions The exam is worth 100 points; each part s point value
More informationAnalyzing Autonomous DEs: Spotted Owls
Analyzing Autonomous DEs: Spotted Owls A group of biologists are making predictions about the spotted owl population in a forest in the Pacific Northwest. The autonomous differential equation the scientist
More informationDynamical Systems: Ecological Modeling
Dynamical Systems: Ecological Modeling G Söderbacka Abstract Ecological modeling is becoming increasingly more important for modern engineers. The mathematical language of dynamical systems has been applied
More informationMATH 415, WEEK 11: Bifurcations in Multiple Dimensions, Hopf Bifurcation
MATH 415, WEEK 11: Bifurcations in Multiple Dimensions, Hopf Bifurcation 1 Bifurcations in Multiple Dimensions When we were considering one-dimensional systems, we saw that subtle changes in parameter
More information8.1 Bifurcations of Equilibria
1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations
More information. For each initial condition y(0) = y 0, there exists a. unique solution. In fact, given any point (x, y), there is a unique curve through this point,
1.2. Direction Fields: Graphical Representation of the ODE and its Solution Section Objective(s): Constructing Direction Fields. Interpreting Direction Fields. Definition 1.2.1. A first order ODE of the
More informationarxiv: v4 [cond-mat.dis-nn] 30 Apr 2018
Attractor metadynamics in terms of target points in slow fast systems: adiabatic vs. symmetry protected flow in a recurrent neural network arxiv:1611.00174v4 [cond-mat.dis-nn] 30 Apr 2018 Hendrik Wernecke,
More informationPHYS-1050 Hydrogen Atom Energy Levels Solutions Spring 2013
1 Introduction Read through this information before proceeding on with the lab. 1.1 Energy Levels 1.1.1 Hydrogen Atom A Hydrogen atom consists of a proton and an electron which are bound together the proton
More informationLecture 1. Scott Pauls 1 3/28/07. Dartmouth College. Math 23, Spring Scott Pauls. Administrivia. Today s material.
Lecture 1 1 1 Department of Mathematics Dartmouth College 3/28/07 Outline Course Overview http://www.math.dartmouth.edu/~m23s07 Matlab Ordinary differential equations Definition An ordinary differential
More informationComplex Numbers. The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition,
Complex Numbers Complex Algebra The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition, and (ii) complex multiplication, (x 1, y 1 )
More informationJulia Sets and the Mandelbrot Set
Julia Sets and the Mandelbrot Set Julia sets are certain fractal sets in the complex plane that arise from the dynamics of complex polynomials. These notes give a brief introduction to Julia sets and explore
More informationMaps and differential equations
Maps and differential equations Marc R. Roussel November 8, 2005 Maps are algebraic rules for computing the next state of dynamical systems in discrete time. Differential equations and maps have a number
More informationMAS212 Assignment #2: The damped driven pendulum
MAS Assignment #: The damped driven pendulum Sam Dolan (January 8 Introduction In this assignment we study the motion of a rigid pendulum of length l and mass m, shown in Fig., using both analytical and
More informationSolutions of a PT-symmetric Dimer with Constant Gain-loss
Solutions of a PT-symmetric Dimer with Constant Gain-loss G14DIS Mathematics 4th Year Dissertation Spring 2012/2013 School of Mathematical Sciences University of Nottingham John Pickton Supervisor: Dr
More informationK. Pyragas* Semiconductor Physics Institute, LT-2600 Vilnius, Lithuania Received 19 March 1998
PHYSICAL REVIEW E VOLUME 58, NUMBER 3 SEPTEMBER 998 Synchronization of coupled time-delay systems: Analytical estimations K. Pyragas* Semiconductor Physics Institute, LT-26 Vilnius, Lithuania Received
More informationMULTISTABILITY IN A BUTTERFLY FLOW
International Journal of Bifurcation and Chaos, Vol. 23, No. 12 (2013) 1350199 (10 pages) c World Scientific Publishing Company DOI: 10.1142/S021812741350199X MULTISTABILITY IN A BUTTERFLY FLOW CHUNBIAO
More informationPhysics: spring-mass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics
Applications of nonlinear ODE systems: Physics: spring-mass system, planet motion, pendulum Chemistry: mixing problems, chemical reactions Biology: ecology problem, neural conduction, epidemics Economy:
More informationSYMBOLIC AND NUMERICAL COMPUTING FOR CHEMICAL KINETIC REACTION SCHEMES
SYMBOLIC AND NUMERICAL COMPUTING FOR CHEMICAL KINETIC REACTION SCHEMES by Mark H. Holmes Yuklun Au J. W. Stayman Department of Mathematical Sciences Rensselaer Polytechnic Institute, Troy, NY, 12180 Abstract
More informationDYNAMICAL SYNTHESIS OF OSCILLATING SYSTEM MASS-SPRING-MAGNET
DYNAMICAL SYNTHESIS OF OSCILLATING SYSTEM MASS-SPRING-MAGNET Rumen NIKOLOV, Todor TODOROV Technical University of Sofia, Bulgaria Abstract. On the basis of a qualitative analysis of nonlinear differential
More informationLecture 7: Simple genetic circuits I
Lecture 7: Simple genetic circuits I Paul C Bressloff (Fall 2018) 7.1 Transcription and translation In Fig. 20 we show the two main stages in the expression of a single gene according to the central dogma.
More informationBehaviour of simple population models under ecological processes
J. Biosci., Vol. 19, Number 2, June 1994, pp 247 254. Printed in India. Behaviour of simple population models under ecological processes SOMDATTA SINHA* and S PARTHASARATHY Centre for Cellular and Molecular
More informationSupplementary Figure 3
Supplementary Figure 3 a 1 (i) (ii) (iii) (iv) (v) log P gene Q group, % ~ ε nominal 2 1 1 8 6 5 A B C D D' G J L M P R U + + ε~ A C B D D G JL M P R U -1 1 ε~ (vi) Z group 2 1 1 (vii) (viii) Z module
More informationSmith Chart Figure 1 Figure 1.
Smith Chart The Smith chart appeared in 1939 as a graph-based method of simplifying the complex math (that is, calculations involving variables of the form x + jy) needed to describe the characteristics
More informationAdditional on-line material for paper A comprehensive set of simulations of high-velocity collisions between main sequence stars.
Mon. Not. R. Astron. Soc. 000, 1 4 (2004) Printed 10 November 2004 (MN LATEX style file v2.2) Additional on-line material for paper A comprehensive set of simulations of high-velocity collisions between
More informationSlope Fields: Graphing Solutions Without the Solutions
8 Slope Fields: Graphing Solutions Without the Solutions Up to now, our efforts have been directed mainly towards finding formulas or equations describing solutions to given differential equations. Then,
More informationControlling the cortex state transitions by altering the oscillation energy
Controlling the cortex state transitions by altering the oscillation energy Valery Tereshko a and Alexander N. Pisarchik b a School of Computing, University of Paisley, Paisley PA 2BE, Scotland b Centro
More information2 Discrete growth models, logistic map (Murray, Chapter 2)
2 Discrete growth models, logistic map (Murray, Chapter 2) As argued in Lecture 1 the population of non-overlapping generations can be modelled as a discrete dynamical system. This is an example of an
More informationLecture 1: A Preliminary to Nonlinear Dynamics and Chaos
Lecture 1: A Preliminary to Nonlinear Dynamics and Chaos Autonomous Systems A set of coupled autonomous 1st-order ODEs. Here "autonomous" means that the right hand side of the equations does not explicitly
More information1.5 Phase Line and Bifurcation Diagrams
1.5 Phase Line and Bifurcation Diagrams 49 1.5 Phase Line and Bifurcation Diagrams Technical publications may use special diagrams to display qualitative information about the equilibrium points of the
More informationA New Science : Chaos
A New Science : Chaos Li Shi Hai Department of Mathematics National University of Singapore In the new movie Jurassic Park [C], Malcolm, a mathematician specialized in Chaos Theory, explained that Hammond's
More informationIntroduction. First Experiment
Course : Bsc Applied Physical Science(Computer Science) IInd Year (Semester IV) Paper no : 14 Paper title : Electromagnetic Theory Lecture No : 14 Tittle : Faraday s Law of Induction Introduction Hello
More informationNonlinear Control Lecture 2:Phase Plane Analysis
Nonlinear Control Lecture 2:Phase Plane Analysis Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 r. Farzaneh Abdollahi Nonlinear Control Lecture 2 1/53
More informationAnalysis of Mode Transitions in Biological Networks
Analysis of Mode Transitions in Biological Networks Nael H. El-Farra, Adiwinata Gani, and Panagiotis D. Christofides Dept. of Chemical Engineering, University of California, Los Angeles, CA 90095 DOI 10.1002/aic.10499
More informationAnalyzing Nepal earthquake epicenters
Analyzing Nepal earthquake epicenters The Nepal Earthquake Epicenters map shows epicenters of the earthquakes that occurred in and around Nepal. The year of the earthquake, its epicenter, and its magnitude
More informationMathematical Foundations of Neuroscience - Lecture 7. Bifurcations II.
Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland Winter 2009/2010 Filip
More informationChaos in the Hénon-Heiles system
Chaos in the Hénon-Heiles system University of Karlstad Christian Emanuelsson Analytical Mechanics FYGC04 Abstract This paper briefly describes how the Hénon-Helies system exhibits chaos. First some subjects
More informationDesigning Information Devices and Systems I Spring 2019 Lecture Notes Note 2
EECS 6A Designing Information Devices and Systems I Spring 9 Lecture Notes Note Vectors and Matrices In the previous note, we introduced vectors and matrices as a way of writing systems of linear equations
More informationComputational Models of Human Cognition
Computational Models of Human Cognition Models A model is a means of representing the structure or workings of a system or object. e.g., model car e.g., economic model e.g., psychophysics model 2500 loudness
More informationThe Big Picture. Python environment? Discuss Examples of unpredictability. Chaos, Scientific American (1986)
The Big Picture Python environment? Discuss Examples of unpredictability Email homework to me: chaos@cse.ucdavis.edu Chaos, Scientific American (1986) Odds, Stanislaw Lem, The New Yorker (1974) 1 Nonlinear
More informationCompensator Design to Improve Transient Performance Using Root Locus
1 Compensator Design to Improve Transient Performance Using Root Locus Prof. Guy Beale Electrical and Computer Engineering Department George Mason University Fairfax, Virginia Correspondence concerning
More informationChapter 1. Introduction
Chapter 1 Introduction 1.1 What is Phase-Locked Loop? The phase-locked loop (PLL) is an electronic system which has numerous important applications. It consists of three elements forming a feedback loop:
More informationWhat happens when things change. Transient current and voltage relationships in a simple resistive circuit.
Module 4 AC Theory What happens when things change. What you'll learn in Module 4. 4.1 Resistors in DC Circuits Transient events in DC circuits. The difference between Ideal and Practical circuits Transient
More informationTime varying networks and the weakness of strong ties
Supplementary Materials Time varying networks and the weakness of strong ties M. Karsai, N. Perra and A. Vespignani 1 Measures of egocentric network evolutions by directed communications In the main text
More informationProject Management Prof. Raghunandan Sengupta Department of Industrial and Management Engineering Indian Institute of Technology Kanpur
Project Management Prof. Raghunandan Sengupta Department of Industrial and Management Engineering Indian Institute of Technology Kanpur Module No # 07 Lecture No # 35 Introduction to Graphical Evaluation
More information520 Chapter 9. Nonlinear Differential Equations and Stability. dt =
5 Chapter 9. Nonlinear Differential Equations and Stabilit dt L dθ. g cos θ cos α Wh was the negative square root chosen in the last equation? (b) If T is the natural period of oscillation, derive the
More informationIntroduction to bifurcations
Introduction to bifurcations Marc R. Roussel September 6, Introduction Most dynamical systems contain parameters in addition to variables. A general system of ordinary differential equations (ODEs) could
More informationReading Question 24.1
Reading Question 24.1 A compass in a magnetic field will line up A. With the north pole pointing in the direction of the magnetic field. B. With the north pole pointing opposite the direction of the magnetic
More informationNonlinear Dynamics. Moreno Marzolla Dip. di Informatica Scienza e Ingegneria (DISI) Università di Bologna.
Nonlinear Dynamics Moreno Marzolla Dip. di Informatica Scienza e Ingegneria (DISI) Università di Bologna http://www.moreno.marzolla.name/ 2 Introduction: Dynamics of Simple Maps 3 Dynamical systems A dynamical
More information) (d o f. For the previous layer in a neural network (just the rightmost layer if a single neuron), the required update equation is: 2.
1 Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.034 Artificial Intelligence, Fall 2011 Recitation 8, November 3 Corrected Version & (most) solutions
More informationA minimal mathematical model combining several regulatory cycles from the budding yeast cell cycle
A minimal mathematical model combining several regulatory cycles from the budding yeast cell cycle K. Sriram, G. Bernot and F. Képès Abstract: A novel topology of regulatory networks abstracted from the
More informationTWO DIMENSIONAL FLOWS. Lecture 5: Limit Cycles and Bifurcations
TWO DIMENSIONAL FLOWS Lecture 5: Limit Cycles and Bifurcations 5. Limit cycles A limit cycle is an isolated closed trajectory [ isolated means that neighbouring trajectories are not closed] Fig. 5.1.1
More informationElectrophysiology of the neuron
School of Mathematical Sciences G4TNS Theoretical Neuroscience Electrophysiology of the neuron Electrophysiology is the study of ionic currents and electrical activity in cells and tissues. The work of
More informationDiscrete Time Coupled Logistic Equations with Symmetric Dispersal
Discrete Time Coupled Logistic Equations with Symmetric Dispersal Tasia Raymer Department of Mathematics araymer@math.ucdavis.edu Abstract: A simple two patch logistic model with symmetric dispersal between
More informationAP Physics C Mechanics Objectives
AP Physics C Mechanics Objectives I. KINEMATICS A. Motion in One Dimension 1. The relationships among position, velocity and acceleration a. Given a graph of position vs. time, identify or sketch a graph
More information