brain rhythms
play

Brain Rhythms Sue Ann Campbell Department of Applied Mathematics - PowerPoint PPT Presentation

Brain Rhythms Sue Ann Campbell Department of Applied Mathematics University of Waterloo October 14, 2017 Outline Biological Background 1 Mathematical Background 2 Modelling Rhythms 3 Summary 4 Sue Ann Campbell (Waterloo) FILOMACS 2 /


  1. Brain Rhythms Sue Ann Campbell Department of Applied Mathematics University of Waterloo October 14, 2017

  2. Outline Biological Background 1 Mathematical Background 2 Modelling Rhythms 3 Summary 4 Sue Ann Campbell (Waterloo) FILOMACS 2 / 30

  3. What are brain rhythms? The electrical activity of the brain exhibits characteristic wave forms

  4. What are brain rhythms? The electrical activity of the brain exhibits characteristic wave forms These vary depending on the brain state and are characterized by their frequency

  5. How do brain rhythms arise? Need to understand a bit of physiology

  6. How do brain rhythms arise? A magnified slice of the brain

  7. How do brain rhythms arise? An individual neuron

  8. How do brain rhythms arise? Behaviour of individual neurons

  9. How do brain rhythms arise? Behaviour of brain - network of neurons

  10. What is an oscillator? Anything that varies periodically in time

  11. What is an oscillator? Anything that varies periodically in time Pendulum

  12. What is an oscillator? Anything that varies periodically in time Metronome

  13. What is an oscillator? Anything that varies periodically in time Firefly

  14. What is an oscillator? Anything that varies periodically in time 3 2 θ 1 0 2 4 6 8 t –1 –2 –3

  15. Neurons are oscillators

  16. What happens when oscillators are connected? Movies Brain slice: https://www.youtube.com/watch?v=t3TaMU_qXMc Fireflies: https://www.youtube.com/watch?v=ZGvtnE1Wy6U 32 metronomes: https://www.youtube.com/watch?v=5v5eBf2KwF8 2 metronomes: https://www.youtube.com/watch?v=yysnkY4WHyM

  17. What happens when oscillators are connected? Synchronization - oscillators all reach maximum at same time

  18. What happens when oscillators are connected? Synchronization - oscillators all reach maximum at same time Phase-locking - oscillators have fixed phase difference

  19. Some Mathematics - Recursively Defined Sequences Recall the Fibonacci Sequence: 1 , 1 , 2 , 3 , 5 , 8 , 13 , . . .

  20. Some Mathematics - Recursively Defined Sequences Recall the Fibonacci Sequence: 1 , 1 , 2 , 3 , 5 , 8 , 13 , . . . Can write a general formula/algorithm for generating the terms of this sequence x 0 = 1 = 1 x 1 x n = x n − 1 + x n − 2 , n = 2 , 3 , 4 , . . .

  21. Some Mathematics - Recursively Defined Sequences Recall the Fibonacci Sequence: 1 , 1 , 2 , 3 , 5 , 8 , 13 , . . . Can write a general formula/algorithm for generating the terms of this sequence x 0 = 1 = 1 x 1 x n = x n − 1 + x n − 2 , n = 2 , 3 , 4 , . . . A recursively defined sequence

  22. Modelling With Recursively Defined Sequences Let x n represent the value of some variable at time n Basic Idea: current value = (previous value) + (change)

  23. Modelling With Recursively Defined Sequences Let x n represent the value of some variable at time n Basic Idea: current value = (previous value) + (change) x n + 1 = x n + ∆

  24. Modelling With Recursively Defined Sequences Let x n represent the value of some variable at time n Basic Idea: current value = (previous value) + (change) x n + 1 = x n + ∆ Assumption: change in current value only depends on previous value

  25. Modelling With Recursively Defined Sequences Let x n represent the value of some variable at time n Basic Idea: current value = (previous value) + (change) x n + 1 = x n + ∆ Assumption: change in current value only depends on previous value x n + 1 = x n + g ( x n )

  26. Modelling With Recursively Defined Sequences Let x n represent the value of some variable at time n Basic Idea: current value = (previous value) + (change) x n + 1 = x n + ∆ Assumption: change in current value only depends on previous value x n + 1 = x n + g ( x n ) Recursive formula for an unknown sequence x 0 , x 1 , x 2 , x 3 , ...

  27. Simple Model Assumption: Change is proportional to current value g ( x n ) = ax n where a is a constant.

  28. Simple Model Assumption: Change is proportional to current value g ( x n ) = ax n where a is a constant. Model: x n + 1 = x n + ax n = ( 1 + a ) x n

  29. Simple Model Solving the Model Let starting value (time 0) be arbitrary: x 0

  30. Simple Model Solving the Model Let starting value (time 0) be arbitrary: x 0 Day 1: x 1 = ( 1 + a ) x 0 = ( 1 + a ) x 0

  31. Simple Model Solving the Model Let starting value (time 0) be arbitrary: x 0 Day 1: x 1 = ( 1 + a ) x 0 = ( 1 + a ) x 0 Day 2: x 2 = ( 1 + a ) x 1 ( 1 + a ) 2 x 0 =

  32. Simple Model Solving the Model Let starting value (time 0) be arbitrary: x 0 Day 1: x 1 = ( 1 + a ) x 0 = ( 1 + a ) x 0 Day 2: x 2 = ( 1 + a ) x 1 ( 1 + a ) 2 x 0 = ( 1 + a ) n x 0 Day n : x n =

  33. Simple Model Solving the Model Let starting value (time 0) be arbitrary: x 0 Day 1: x 1 = ( 1 + a ) x 0 = ( 1 + a ) x 0 Day 2: x 2 = ( 1 + a ) x 1 ( 1 + a ) 2 x 0 = ( 1 + a ) n x 0 Day n : x n = A geometric sequence

  34. Simple Model Solving the Model Let starting value (time 0) be arbitrary: x 0 Day 1: x 1 = ( 1 + a ) x 0 = ( 1 + a ) x 0 Day 2: x 2 = ( 1 + a ) x 1 ( 1 + a ) 2 x 0 = ( 1 + a ) n x 0 Day n : x n = A geometric sequence What happens to population if a = 0? a > 0? − 1 < a < 0?

  35. Plotting the Solution a = 0 , x 0 = 5 a = 0 : Constant solution

  36. Plotting the Solution a > 0 , x 0 = 1 a > 0 : Exponential Growth

  37. Plotting the Solution − 1 < a < 0 , x 0 = 10 a < 0 : Exponential Decay

  38. Modelling Oscillators Define the phase of the oscillator Amplitude t (time) t 0 t+ ∆ t t +T 0 0 ∆ t/T 0 1 θ (phase) 2π 2π∆ t/T 0

  39. Modelling Oscillators with Recursive Sequences Think of oscillator as angle of point moving on circle with fixed radius y θ x Change in θ in a small time interval θ n + 1 = θ n + Ω

  40. Two Coupled Oscillators θ 1 , n + 1 = θ 1 n + Ω + A sin( θ 2 n − θ 1 n ) = θ 2 n + Ω + A sin( θ 1 n − θ 2 n ) θ 2 , n + 1 where A is a small number ( | A | << 1 ).

  41. Two Coupled Oscillators θ 1 , n + 1 = θ 1 n + Ω + A sin( θ 2 n − θ 1 n ) = θ 2 n + Ω + A sin( θ 1 n − θ 2 n ) θ 2 , n + 1 where A is a small number ( | A | << 1 ). Define phase difference φ n = θ 2 n − θ 1 n φ n + 1 = φ n − 2 A sin( φ n )

  42. Two Coupled Oscillators Simulations of model with A = 0 . 1 φ 0 = 6 , 5 , 4 , 3 , 2 , 1 , 0 φ n + 1 = φ n − 2 A sin( φ n )

  43. Two Coupled Oscillators Simulations of model with A = 0 . 1 φ 0 = 6 , 5 , 4 , 3 , 2 , 1 , 0 φ n + 1 = φ n − 2 A sin( φ n )

  44. Two Coupled Oscillators Simulations of model with A = 0 . 1 φ 0 = 6 , 5 , 4 , 3 , 2 , 1 , 0 φ n + 1 = φ n − 2 A sin( φ n )

  45. Two Coupled Oscillators Simulations of model with A = 0 . 1 φ 0 = 6 , 5 , 4 , 3 , 2 , 1 , 0 φ n + 1 = φ n − 2 A sin( φ n )

  46. Two Coupled Oscillators Simulations of model with A = 0 . 1 φ 0 = 6 , 5 , 4 , 3 , 2 , 1 , 0 φ n + 1 = φ n − 2 A sin( φ n )

  47. Two Coupled Oscillators Simulations of model with A = 0 . 1 φ 0 = 6 , 5 , 4 , 3 , 2 , 1 , 0 φ n + 1 = φ n − 2 A sin( φ n )

  48. Two Coupled Oscillators Simulations of model with A = 0 . 1 φ 0 = 6 , 5 , 4 , 3 , 2 , 1 , 0 φ n + 1 = φ n − 2 A sin( φ n )

  49. Two Coupled Oscillators Simulations of model with A = 0 . 1 φ 0 = 6 , 5 , 4 , 3 , 2 , 1 , 0 φ n + 1 = φ n − 2 A sin( φ n )

  50. Two Coupled Oscillators φ n + 1 = φ n − 2 A sin( φ n ) Special constant solutions: φ n = φ ∗ , n = 1 , 2 , . . . (equilibrium solutions) Correspond to state the phase difference between the two oscillators is fixed (phase locking).

  51. Two Coupled Oscillators φ n + 1 = φ n − 2 A sin( φ n ) Special constant solutions: φ n = φ ∗ , n = 1 , 2 , . . . (equilibrium solutions) Correspond to state the phase difference between the two oscillators is fixed (phase locking). Occur when sin( φ ∗ ) = 0

  52. Two Coupled Oscillators φ n + 1 = φ n − 2 A sin( φ n ) Two possibilities φ = 0 ( 2 π ) : periodic solutions with two oscillators in-phase φ = π : periodic solutions with two oscillators anti-phase (one half period out of phase).

  53. Two Coupled Oscillators

  54. Two Coupled Oscillators φ n + 1 = φ n − 2 A sin( φ n ) How the equilibrium solutions affect the behaviour depends on the sign of A .

  55. Two Coupled Oscillators φ n + 1 = φ n − 2 A sin( φ n ) How the equilibrium solutions affect the behaviour depends on the sign of A . A > 0

  56. Two Coupled Oscillators φ n + 1 = φ n − 2 A sin( φ n ) How the equilibrium solutions affect the behaviour depends on the sign of A . A < 0

  57. Two Coupled Oscillators φ n + 1 = φ n − 2 A sin( φ n ) How the equilibrium solutions affect the behaviour depends on the sign of A . A < 0

  58. Two Coupled Oscillators - Different Frequencies θ 1 , n + 1 = Ω 1 + A sin( θ 2 n − θ 1 n ) θ 2 , n + 1 = Ω 2 + A sin( θ 1 n − θ 2 n ) ⇓ φ n + 1 = ω + φ n − 2 A sin( φ n )

  59. Two Coupled Oscillators - Different Frequencies θ 1 , n + 1 = Ω 1 + A sin( θ 2 n − θ 1 n ) = Ω 2 + A sin( θ 1 n − θ 2 n ) θ 2 , n + 1 ⇓ = ω + φ n − 2 A sin( φ n ) φ n + 1 Equilibrium solutions: φ ∗ such that sin( φ ∗ ) = ω .

Recommend


More recommend