Marys Medicine

Mathematical models of bipolar disorder

Commun Nonlinear Sci Numer Simulat 14 (2009) 2897–2908 Contents lists available at Commun Nonlinear Sci Numer Simulat Mathematical models of bipolar disorder Darryl Daugherty a, Tairi Roque-Urrea a, John Urrea-Roque a, Jessica Troyer b, Stephen Wirkus c,*,Mason A. Porter d a Department of Mathematical Sciences, Binghamton University, Binghamton, NY 13902-6000, USAb L-3 Communications Corporation, Huntsville, AL 35806, USAc Mathematical and Natural Sciences Division, Arizona State University, Glendale, AZ 85306-4908, USAd Oxford Centre for Industrial and Applied Mathematics, Mathematical Institute, University of Oxford, OX1 3LB, UK We use limit cycle oscillators to model bipolar II disorder, which is characterized by alter- Received 17 September 2008 nating hypomanic and depressive episodes and afflicts about 1% of the United States adult Accepted 26 October 2008 population. We consider two non-linear oscillator models of a single bipolar patient. In Available online 14 November 2008 both frameworks, we begin with an untreated individual and examine the mathematicaleffects and resulting biological consequences of treatment. We also briefly consider the dynamics of interacting bipolar II individuals using weakly-coupled, weakly-damped har- monic oscillators. We discuss how the proposed models can be used as a framework for refined models that incorporate additional biological data. We conclude with a discussion of possible generalizations of our work, as there are several biologically-motivated exten- sions that can be readily incorporated into the series of models presented here.
Limit cycle oscillators Ó 2008 Elsevier B.V. All rights reserved.
We posted the original version of this article on the arXiv preprint server in November 2003. There it lay dormant for several years, although some authors have since built on our work . Meanwhile, we have collectively experienced sev-eral changes in our institutional affiliations (with an average of at least one move per author) and personal lives (including atleast one marriage, at least three children, and several heartbreaking losses in baseball playoff games).
Recently, Richard Rand sent us an article by Nana that was published in the present journal and takes a similar per- spective to our own original paper. Two of us have already commented on Nana's article in Ref. We are pleased that ouroriginal manuscript has now been revived, so that we can finally present our research in published form. Our hope is that thetoy models and analysis presented below (constituting the second version of Ref. which we posted on the arXiv server inJune 2004) will ultimately be refined considerably with the incorporation of clinical data.
2. Characterizations of bipolar II disorder About 1% of the United States adult population is afflicted with bipolar disorder (manic depression) which poses myriad difficulties to clinical practitioners. It is difficult to diagnose, as bipolar patients often do not adhere to treatment * Corresponding author. Tel.: +1 602 543 8236.
E-mail address: (S. Wirkus).
1007-5704/$ - see front matter Ó 2008 Elsevier B.V. All rights reserved.
doi:10.1016/j.cnsns.2008.10.027 D. Daugherty et al. / Commun Nonlinear Sci Numer Simulat 14 (2009) 2897–2908 and/or medication, and most drugs can be toxic if taken individually . The available data on bipolar disorder is sparse—due in part to a lack of agreement on well-suited trial design—so it is very difficult to study using clinical trials .
Psychiatrists have established a broad range of criteria for classifying bipolar disorder; see, for example, the Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition (DSM-IV) Its characteristics, which need not all be present in everybipolar individual, include mixed episodes, in which it is possible to simultaneously experience symptoms of both mania anddepression, and rapid cycling, in which a patient experiences at least four cycles per year .
There are two primary types of bipolar disorder. Bipolar I disorder is characterized by a combination of manic and depres- sive episodes with the possibility of mixed episodes, whereas bipolar II disorder is characterized by a combination of hypo-and depressive episodes Bipolar II patients tend to be more prone to rapid cycling—especially if initially treatedonly with antidepressants In this paper, we focus on bipolar II disorder for two reasons. First, individuals who sufferfrom it are more often observed to exhibit approximately periodic mood swings than bipolar I individuals. Second, bipolar IIindividuals do not experience mixed episodes.
Bipolar II disorder is often misdiagnosed as either unipolar depression or a severe personality disorder To be correctly diagnosed, a patient seeking treatment must give an accurate description of his past behavior. Because of the nature of hypo-manic episodes (increased energy, decreased need for sleep, etc.) patients often do not describe these conditions to doctors andare therefore diagnosed with unipolar depression . According to the results of the National Depressive and Manic-DepressiveAssociation 2000 survey of people with bipolar disorder, over one third of respondents sought professional help within a year ofonset of symptoms. However, it took as many as 10 years and four physicians for some patients to be correctly diagnosed .
Treatment for bipolar disorder ideally includes a combination of medication and therapy. Typical drug treatments include mood stabilizers, antipsychotics, antidepressants, and select anticonvulsants. Among the more commonly used drugs areLithium, Valproate (also know as Depakote), Carbamazepine (also known as Tegretol), and Prozac. Mood controlling drugssuch as Lithium take 4–10 days to reach therapeutic levels in the blood stream, so initial treatment is likely to include bothantidepressants and antipsychotics . During the maintenance state, antidepressants and antipsychotics are generally sup-plemented by mood stabilizers. Monotherapy (i.e., single-drug therapy) is generally avoided by clinicians because of thestrong side effects of some of the drugs used for treatment. Other drugs that have been employed include selective serotoninreuptake inhibitors (SSRI) and monoamineoxidase inhibitors (MAOI), both of which are typically used for depression. In somecases, special care must be taken to ensure that a bipolar individual does not fall into a pattern of rapid cycling or becomeaddicted to the medication used for treatment. In fact, substance abuse has been associated with bipolar disorder .
Bipolar II disorder is also known to be highly heritable. It has been reported, for example, that the offspring of people with bipolar II disorder have a 35% chance of being afflicted as well. In particular, there are known cases of family units with mul-tiple bipolar II individuals, so the dynamics of interacting bipolar patients is also of interest to psychologists In this paper, we propose and examine two mathematical models that attempt to represent the qualitative dynamics of individual bipolar patients. We also consider a closely interacting pair of bipolar patients (which can occur, for example, inhouseholds with bipolar sibling pairs The importance of the present work is that it suggests ways of thinking about bipolar disorder. In particular, mathemat- ical models have the potential to provide significant insight into this disorder, provided there is adequate data to constrainthe models. Studies using time-series analysis (of a relatively small number of bipolar and normal individuals) suggest thatrelatively simple mechanisms may be responsible for the complex mood variations in bipolar disorder The models wedevelop provide a first step toward using mathematical modeling to increase scientific understanding of this disorder. Ulti-mately, our models will need refinement to allow greater predictability and ties to candidate biological processes. Our hopeis for the present work to motivate the collection of the additional time-series data necessary to achieve this goal.
3. Scope of the paper Numerous research articles discuss bipolar disorder, focusing primarily on clinical trials and case studies. Unfortunately, very little of this prior research is amenable to mathematical treatment, whose employment could prove vital in advancingour knowledge of this medical condition . Nevertheless, several studies of bipolar disorder—including the recent compre-hensive article the moods of bipolar patients in terms of the fraction of weeks during a year when partic-ular symptoms were (or were not) displayed. While the models we examine in the present paper can also be interpreted inthese terms, a weekly or daily time-series description of the symptoms could inspire the creation or modification of math-ematical models that may shed additional light on the underlying dynamics of bipolar disorder. For example, early work byWehr and Goodwin (which demonstrates oscillatory behavior in the moods of bipolar individuals) could be used, in princi-ple, to help refine our mathematical framework. However, their time-series description included only one bipolar patient, soits practical utility in the development of a mathematical description of bipolar II disorder is limited A comprehensive study (with a large number of clinical trials) like that of Judd et al. that uses the time-series format of Wehr and Goodwin could lead to far more realistic mathematical models than the available data currently allows. Although 1 The main distinction between hypomanic episodes and manic ones is that the former are much less severe than the latter. Additionally, mania can last much longer than a week, whereas hypomania has been shown to have a median duration of about 1–2 days and need only last 4 days to reach DSM-IVcriteria .
D. Daugherty et al. / Commun Nonlinear Sci Numer Simulat 14 (2009) 2897–2908 arduous, such an undertaking would yield important advancements in modeling and understanding bipolar disorder. Ourgoal in this paper is considerably more modest—we use dynamical systems theory to develop minimalist models of bipolarII disorder that we hope will lead to additional qualitative understanding of the behavioral oscillations associated with thismedical condition. The novelty of our work lies in our mathematical approach to the modeling of bipolar disorder.
Although there exist biological models for this medical condition, mathematical equations have never (to our knowledge) previously been employed to study bipolar individuals. In this work, we propose and analyze two mathematical models ofbipolar II disorder. Rather than focusing on the method of diagnosis (which is a difficult medical problem), we instead con-centrate on the dynamics of our models under the proposed treatment strategy, which is understood to involve a combina-tion of drugs and/or therapy. The models we study are not meant to represent a specific treatment but are instead intendedto provide some insight into the complicated dynamics of bipolar disorder. In this respect, we consider our work to be a firststep in developing mathematical models of the mood swings of bipolar individuals. The shift in thinking from a statisticalmodel to a non-linear mathematical model for mood variations has been suggested by some authors, including Totterdellet al. Additionally, Ehlers , Frank et al. and Gottschalk et al. have noted the importance of non-linearand chaotic dynamics and the employment of simple models incorporating such features to achieve better understandingof mood and behavior. In fact, the latter authors attempted to identify and quantify the non-linear behavior in the moodof bipolar patients by using time-series analysis to study power spectra and fractal dimension. However, they did not at-tempt to develop mathematical equations to model mood swings in bipolar individuals (although they did speculate thatutilizing van der Pol oscillators may be useful), which is the goal of this work. In this respect, our theoretical work nicelycomplements the data analysis of Gottschalk et al.
Reinterpreting the time-series results of Totterdell et al. and gathering similar time-series data for bipolar patients has considerable merit. We hope that our work in developing dynamical equations describing the mood swings of bipolar indi-viduals, in conjunction with additional studies like those of Totterdell et al., and Gottschalk et al., leads to the eventual devel-opment of more detailed models that incorporate clinical data from both bipolar and ‘‘normal" individuals.
Numerous studies have proposed possible underlying mechanisms for bipolar disorder, and some have even examined the effect of light stimuli and seasonal variations. Potentially, insightful mathematical models can be motivated from suchwork, but the perspective we take is to employ minimalist mathematical models describing biological and physical oscilla-tors that assume a priori the existence of asymptotic oscillations (i.e., limit cycles) . This allows us to gain insight into thedynamics of the mood swings of bipolar II individuals, despite the fact that the underlying dynamics are not fully understoodat the level of chemicals in the body. This kind of minimalist perspective has been employed successfully to study, for exam-ple, the circadian rhythms of the avian chick eye and the dynamics of coupled microwave oscillators Minimalist models of this sort are not intended to suggest the mechanisms that underly a given phenomenon but rather to gain a qualitative understanding of the dynamics of the existing mechanisms (especially in situations, like this one, inwhich they are not known or poorly understood). We also hope that our work will spur the data collection necessary to de-velop more realistic mathematical models that can be tied closer to the application at hand. Mathematical studies have thepotential to help our understanding of bipolar disorder considerably, but a lot more work must be done to reach that stage.
This paper is intended to be a first step along that path.
4. Limit cycle oscillators for bipolar II individuals For modeling purposes, some simplifying assumptions are necessary. First, although bipolar II disorder can be somewhat erratic, episodes are known to exhibit recurrent patterns. For a group of patients with this disorder, there is a periodicitygoverning the manic and depressive episodes . Moreover, it is commonly assumed that the disorder will progress se-verely if left untreated.
As previously mentioned, we are not attempting to provide an explanation of the underlying mechanisms of bipolar dis- order (although we hope that we can provide a stepping-stone toward achieving this highly desirable objective). Rather, ourimmediate goal is to better understand the dynamics of the mood swings by assuming a priori the existence of an oscillatoror oscillators that might approximate the observed behavior of bipolar individuals.
of early work by Wehr and Goodwin suggests that the behavior of the single bipolar patient they studied can be described qualitatively as a limit cycle oscillator . Moreover, a typical bipolar patient undergoes roughly four symptomchanges per year, which corresponds roughly to observations of Wehr and Goodwin The rapid cycling that occursunder the treatment of desipramine hydrochloride also lends itself to a mathematical interpretation in terms of limit cycleoscillators. Our analysis provides insight into how these limit cycle oscillations are induced with treatment.
As hypomanic and depressive episodes are periodic and—if untreated—increase in severity over time, one can model the mood and mood change of an untreated individual using a negatively damped harmonic oscillator, where x denotes the emotional state of the patient, _x is the rate of mood changes between hypomania and severe depression,and a > 0 and x are parameters. The main drawback of such a simple model is its unbounded oscillations.
Most patients with bipolar II disorder are diagnosed when they are in a depressive episode, as hypomania ordinarily does not prevent the normal function of individuals. (Mania, in contrast, causes significant occupational and social discomfort.) In D. Daugherty et al. / Commun Nonlinear Sci Numer Simulat 14 (2009) 2897–2908 many cases, hypomania can even enhance short-term functionality Because it often takes several years for bipolar dis-order to be properly diagnosed, we assume treatment begins when a patient is in his early 20s. We also suppose that treat-ment can be modeled using an autonomous forcing function (although non-autonomous forcing, involving time-periodicfunctions such as a sequence of delta functions or trigonometric functions, is also worth considering). The inclusion of thisforcing term turns Eq. into the well-known van der Pol oscillator, The forcing function gðx; _xÞ ¼ bx2 _x represents aggregate treatment and includes a combination of antidepressants, mood stabilizers, psychotherapy, and either antipsychotics or tranquilizers. Finding the correct way to treat a given patient maytake as many as 10–15 trials on different medications. The most comprehensive treatment, however, involves both medica-tion and psychotherapy One can suppose that Eq. describes a treated bipolar individual, as it possesses a unique stable limit cycle surrounding the origin It is a specific example of a more general class of equations of the form x þ f ðxÞ_x þ hðxÞ ¼ 0; which is known as the Liénard equation . The presence of a limit cycle indicates that after treatment, the bipolar pa-tient's mood variations (asymptotically) oscillate within a range of x values to be determined by the parameters in Eq.
Although this simplistic model gives some insight into the dynamical properties of periodic mood variations of bipolar individuals, the characterization of untreated individuals is not realistic because of the unbounded oscillations (moodswings) that result from every initial condition. We thus seek a model that not only captures the behavior of treated bipolarII individuals, but also does a better job of capturing the dynamics of untreated individuals.
4.1. Model 1: van der Pol oscillator with autonomous forcing In this model, we use a van der Pol oscillator to represent an untreated bipolar II individual. This characterization is more realistic than that in Eq. , as the mood swings of an untreated patient can still grow large but now approach abounded limit cycle rather than eventually becoming infinitely severe. We again apply treatment in the form of an auton-omous forcing function gðx; _xÞ, x  a _x þ x2x  bx2 _x ¼ gðx; _xÞ: In our analysis, we used gðx; _xÞ ¼ cx4 _x þ dx2 _x: For our simulations, we considered the case d = 0, but one can apply a medication of the more general form to adjust thecoefficient b.
Because normal individuals also have mood swings, one must designate how severe such mood variation must be in order for someone to be diagnosed as bipolar. In other words, only limit cycles with some minimal amplitude (say, 0.1) correspondto bipolar mood swings. In an untreated patient described by with g(x)  0, the parameters a = 0.36, b = 100, and x = 5lead to a limit cycle with an amplitude of approximately 0.12, as shown in . In the context of the present model,such an amplitude threshold dividing functional individuals from bipolar II ones is arbitrary, so the scales in can be ad-justed to account for individuals who suffer from bipolar disorder to varying degree (e.g., individuals who have mania ratherthan hypomania).
Starting with the van der Pol oscillator caricature of the untreated patient, we apply the forcing function gðx; _xÞ to repre- sent medication. We use to represent a treated bipolar II patient, who possesses an unstable, large limit cycle (greaterthan the 0.1 amplitude, which is the threshold that we have chosen for ‘normal' behavior) and a stable, small limit cycle.
Choosing a = 0.1, b = 100, and c = 5000 yields the plots in Using the method of averaging we study the number and amplitude (denoted A) of the limit cycles of Eq. The dynamics of A are described by slow flow equations of motion, whose non-zero equilibria correspond to limit cycles of Eq.
Depending on the values of the parameters, there can be zero, one, or two limit cycles. The slow flow equations forthe amplitude contain both positive and negative equilibria, but only the non-negative solutions are relevant to the originaldynamical system. Limit cycle amplitudes must satisfy Once we find the limit cycle amplitudes, we compute the Jacobian (in this case, the second derivative) to determine theirstability as functions of a, b, and c. There is always a slow-flow equilibrium at A = 0, which corresponds to the equilibriumat (0, 0) of . We obtain the following conditions for the existence and stability of limit cycles: 1. Zero limit cycles occur when either b2  8ac < 0 (so that A2 is not real) or a, b, and c all have the same sign (so that A2 < 0).
In these cases, there are no slow-flow equilibria with positive A.
D. Daugherty et al. / Commun Nonlinear Sci Numer Simulat 14 (2009) 2897–2908 Emotional State x Fig. 1. Time-series representing the behavior of an untreated patient. The parameters are a = 0.36, b = 100, and x = 5. This leads to a stable limit cyclewith an amplitude of about 0.12, which we designate as larger than that describing ‘normal' mood swings.
Rate of Mood Change Emotional State x Fig. 2. Limit cycle corresponding to the time-series in . This describes the relationship between the individual's moods and how fast they change.
2. One limit cycle:  There is one limit cycle when b/c > 0 and b2 = 8ac. There is a bifurcation at this point (corresponding to the coalescence of two limit cycles), and stability cannot be determined by computing a Jacobian.
 There is one limit cycle when b2  8ac >j b j, regardless of the sign of b/c. This occurs exactly when ac < 0. Such a limit cycle is always unstable.
 There is one limit cycle when a = 0 and b/c < 0. In this situation, the Jacobian is b2/2c, so the limit cycle is stable if and only if c > 0.
D. Daugherty et al. / Commun Nonlinear Sci Numer Simulat 14 (2009) 2897–2908 Emotional State x –0.05 Fig. 3. Treated patient with parameters a = 0.1, b = 100, and c = 5000, which yields an unstable limit cycle with amplitude greater than 0.1 and a stablelimit cycle with amplitude less than 0.1. This patient begins to receive treatment at about age 5 or 10. Individuals usually first become bipolar between theages of 18 and 22, but the onset of bipolar disorder can occur in early childhood Rate of Mood Change Emotional State x Fig. 4. Limit cycles corresponding to the time-series in . The depicted trajectory approaches the smaller (stable) limit cycle as the time t increases. Theamplitude of this limit cycle is less than 0.1, so it describes the mood swings of a functional individual.
3. Two limit cycles occur if b/c < 0 and b2 > 8ac > 0. This happens for c, a > 0 when b < 0 and for c, a < 0 when b > 0. In the biologically relevant situation describing a treated individual, the smaller limit cycle is stable. An example of this isdepicted in with the parameter values a = 0.1, b = 100, and c = 5000.

D. Daugherty et al. / Commun Nonlinear Sci Numer Simulat 14 (2009) 2897–2908 4.2. Model 2: Liénard oscillator with autonomous forcing The introduction of a forcing function (representing treatment) to Model I results in smaller mood swings; in so doing, however, a larger unstable limit cycle was introduced. This has the potential drawback of predicting that if an individual goestoo long without being diagnosed (and thus the amplitude of the mood swings is too large), one would need to ensure thatthe initial condition when treatment begins is within the basin of attraction of the smaller stable limit cycle. In other words,if the mood swings are too large, it might be necessary to drastically reduce the mood amplitude before introducing ‘normal'treatment.
Model 2 provides an alternative minimalist framework to study bipolar II disorder that does not have this drawback. In this situation, we will demonstrate the presence of three limit cycles. We will also show that the ones with the largest andsmallest amplitude are stable, whereas the limit cycle between them is unstable. Toward this end, we consider an equationof the form x þ f ðxÞ _x þ hðxÞ ¼ gðx; _xÞ; gðx; _xÞ ¼ q _x3 þ l _x5 þ m_x11: For this model, we consider the constant function f(x) = 0.38 and the linear function h(x) = 180x together with the param-eters l = 0.78 and m = 0.00093. With q = 0.38, we obtain a large stable limit cycle with mood amplitude approximatelyequal to 0.44. As with Model 1, an untreated patient has bounded mood swings.
Without treatment, this model describes an individual with steadily worsening mood swings throughout his childhood and adolescent years. At approximately age 20, the individual's mood variations increase dramatically to the point wherethe individual can be clinically diagnosed with bipolar II disorder. For the given model and any initial conditions, the indi-vidual's mood variations settle asymptotically to the stable limit cycle with mood swing amplitude jxj 6 0.45, as shown in. As this system's limit cycle is globally stable, every trajectory eventually spirals toward it, yielding the time-seriesdepicted in Observe that the mood varies between ±0.4465. Moreover, the bipolar II individual described by thistime-series can diagnosed with both hypomania and severe depression at age 20.
In considering the role of treatment gðx; _xÞ in this model, we note that q is the key parameter that will be adjusted. For a treated patient, the value q = 3.302 (for example) yields qualitatively appropriate dynamics. A treated bipolar II patient isthen modeled by x  0:38 _x þ 180x ¼ 3:302 _x3 þ 0:78_x5  0:00093 _x11: This model encompasses all aspects of treatment as a single function gðx; _xÞ. In , we see that by inputting the new initialconditions into the model that includes treatment, the patient exhibits favorable emotional patterns. There is a large stablelimit cycle and an unstable limit cycle just inside it. A smaller stable limit cycle is the desired emotional pattern of the pa-tient; all initial conditions lying inside the large unstable cycle yield solutions that spiral toward this smaller cycle, whoseamplitude is sufficiently small so that the associated mood swings are deemed normal. One obtains the same qualitativebehavior with smaller jqj (such as jqj = 2), although the separation between the unstable limit cycle and the larger stable Fig. 5. Untreated bipolar patient for Model 2. Only a large stable limit cycle is present.

D. Daugherty et al. / Commun Nonlinear Sci Numer Simulat 14 (2009) 2897–2908 Fig. 6. Time-series corresponding to The initial condition lies inside the limit cycle and thus represents an individual who gradually worsens overtime.
Fig. 7. Treated bipolar II patient for Model 2. There is a large outer limit cycle which is stable and an unstable limit cycle lying just inside it. There is also asmaller stable limit cycle inside the unstable one; this orbit represents the desired mood swings of a treated patient.
limit cycle surrounding it becomes larger as jqj is decreased. The chosen value of q yields limit cycles that are almost on topof each other; increasing jqj changes this qualitative structure.
The larger limit cycle prevents a treated patient from having unbounded mood variations that could otherwise occur, for example, if a perturbation were to make the current mood amplitude too large. The smaller stable limit cycle represents thedesired low-amplitude mood variations. Comparing reveals that after treatment, the hypomanic and severedepressive episodes have both disappeared. At the beginning of treatment, the amplitude of the mood swings stays relativelyclose to 0.2. It then tends asymptotically to the stable limit cycle, where the mood x varies between approximately ±0.03. Animportant feature of the model demonstrated in is that treatment can begin at any time because the small limit cycle

D. Daugherty et al. / Commun Nonlinear Sci Numer Simulat 14 (2009) 2897–2908 Fig. 8. Time-series corresponding to . The initial condition lies inside the unstable limit cycle and is thus in the basin of attraction of the smaller stablelimit cycle representing normal mood swings.
attracts all solutions that lie inside the unstable limit cycle. If the mood swings of the bipolar individual are initially too large,it may again be necessary for a brief drastic treatment to bring them into the basin of attraction of the smaller limit cyclebefore applying normal medication.
5. Bipolar patients as interacting oscillators Examining the possible effects two patients have on each other while undergoing treatment is of significant mathematical interest. Moreover, as mentioned previously, there are known cases of family units with multiple bipolar II individuals. Thus,the interaction of two bipolar patients is also of psychological interest . It is known to occur, for example, in householdswith bipolar sibling pairs . There has also been some work on using group treatment as a prophylaxis with patients inremission in order to prevent further episodes We again stress, however, that clinical data needs to be collected to exam-ine these ideas more extensively.
Consider, for example, two treated bipolar II individuals, each of whom is separately described by Eq. We transform into polar coordinates to provide a natural manner of adding coupling terms to represent the interaction between the twopatients. Denoting by Kj(r1, h1, r2, h2) the interactions between the patients, the dynamical system of interest is r1 ¼ br3ðcos h þ r1 sin h1½a sinðh1Þ þ cosðh1Þ  x2 cosðh1Þ þ K1ðr1; h1; r2; h2Þ; h1 ¼ a sinðh1Þ cosðh1Þ  ðsin h1Þ2  x2ðcos h1Þ2  br2ðcos h 1Þ3 sinðh1Þ þ K 2ðr1; h1; r2; h2Þ; r2 ¼ br3 cos h2ðsin h þ r2 sin h2½a sinðh2Þ þ cosðh2Þ  x2 cosðh2Þ þ K3ðr1; h1; r2; h2Þ; h2 ¼ a sinðh2Þ cosðh2Þ  ðsin h2Þ2  x2ðcos h2Þ2  br2ðcos h 2Þ3 sinðh2Þ þ K 4ðr1; h1; r2; h2Þ: For simplicity, we assume that the phase difference affects only the phase terms and the amplitude difference affects only the amplitude terms. We also assume linear proportionality, so that the coupling terms simplify to K1 ¼ k1ðr2  r1Þ; K2 ¼ k2ðh2  h1Þ; K3 ¼ k3ðr1  r2Þ; K4 ¼ k4ðh1  h2Þ; D. Daugherty et al. / Commun Nonlinear Sci Numer Simulat 14 (2009) 2897–2908 where the ki are the proportionality constants. An important detail to note is that the angular coupling terms in includeno limitation on the phase angle, so that h1  h2 can increase beyond 2p. There is no difficulty, however, when we convertback to rectangular coordinates. The coupling terms remain bounded, as the radial coupling terms are proportional to x2 þ _x2 and the angular ones are proportional to arctanð _x=xÞ. When we numerically integrate these dynamical equations, we utilize arctanð _x=xÞ 2 ðp=2; p=2Þ. Although the arctangents in rectangular coordinates create a discontinuous vector field,the behavior near the discontinuities correctly models the biology observed in situations in which two patients with bipolardisorder are interacting. (That is, we observe fluctuations near in-phase and out-of-phase modes).
Other choices of coupling terms avoid unbounded phase differences. Consider, for example, K1 ¼ k1ðr2  r1Þ; K2 ¼ k2 sinðh2  h1Þ; K3 ¼ k3ðr1  r2Þ; K4 ¼ k4 sinðh1  h2Þ; where the ki are again proportionality constants. We have conducted some mathematical analysis of with both and discuss one of our results briefly. Note, however, that although these results are mathematically interesting, we onlymention these models in passing because insufficient clinical data is available to properly constrain the numerous forms ofcoupling available when constructing mathematical models describing interacting bipolar II individuals. Consequently, weonly state results for two common types of solutions in such systems.
In coupled systems, important types of behavior include in-phase modes, for which x1 = x2 and y1 = y2, and out-of-phase modes, for which x1 = x2 and y1 = y2. For the first form of coupling , both in-phase and out-of-phase modes exist andare stable with relatively small basins of attractions. Larger stable motions exist that attract initial conditions starting out-side those basins of attraction. Both in-phase and out-of-phase modes also exist with the second form of coupling andthe former appears to be the only stable motion in the system. Initial conditions starting far away from the in-phase modeeventually approach this mode.
Interpreted biologically, two coupled bipolar II individuals with similar moods (x1  x2) and mood cycling rates (y1  y2) tend to remain in-phase. They are synchronized in the sense that they enter hypomanic and depressive states (roughly) con-currently. The same argument applies to out-of-phase modes, so if two bipolar II individuals are almost completely out-of-phase, they will remain that way. One individual enters a hypomanic phase when the other enters a depressive phase.
In this work, we proposed two limit cycle oscillator models that attempt to explain qualitatively the moods swings of bipolar II individuals. There is insufficient data at this time to construct a quantitatively detailed model describing bipolardisorder, so this first step employing minimalist models is an essential one. Indeed, similar phenomenological models haveyielded significant insight into the dynamics of several oscillatory biological and physical phenomena .
Even without time-series data to describe the moods of bipolar patients, the models we have developed could prove important in increasing the understanding of the effect of treatment on the cyclic behavior of bipolar individuals. Becausethe two-dimensional dynamical systems we utilize provide simple caricatures of the behavior of bipolar individuals, it isconcomitantly easy to analyze these models. Moreover, despite their simplicity, they successfully capture qualitativelythe known behavior of treated and untreated bipolar II individuals. A quantitative analysis of their dynamics, which isour long-term goal, will require refining our models using time-series data from clinical trials. It is our hope that the presentwork will stimulate further mathematical analyses that incorporate such data as well as the data collection that will permitthese modeling efforts.
In Model 1, we considered treatment as a forcing function that led to a stable limit cycle with smaller amplitude, which we interpreted as a decrease in the severity of the mood swings of the bipolar patient. This model suggests that individualsthat are not diagnosed at a sufficiently young age might need drastic intervention to bring their mood swings to a reasonablelevel. Given that individuals are usually diagnosed in the depressive state and that bipolar disorder is frequently misdiag-nosed initially, this aspect of Model 1 is especially noteworthy Model 2 provides an alternative framework in which a bipolar II individual possesses two stable limit cycles even without treatment; satisfactory treatment then confines the patient to the one with smaller amplitude. The rapid cycling induced bydesiramine hydrochloride that was reported by Wehr and Goodwin can be interpreted as a slight increase in the ampli-tude of the oscillation but—more centrally—as an increase in the oscillation frequency. This phenomenon can be addressedwithin the modeling framework we have developed by allowing the frequency parameter to also be perturbed slightly bytreatment.
At this point, it is also important to briefly consider how our models can be both refined and generalized. First, incorpo- rating biological data is of paramount importance. We used minimal data because little data is available in a form that facil-itates the creation or refinement of mathematical models. This is the largest deficiency in our analysis; employing ourmodeling framework along with data from psychological testing of bipolar II individuals will yield a significantly increasedunderstanding of this medical condition. In fact, one of the primary purposes of this work is to motivate the collection of suchdata so that detailed models of bipolar disorder can ultimately be developed. In addition to time-series data describing thebehavior of a large number of bipolar patients, it is also desirable to have data describing the mood variations of ‘‘normal"individuals to provide a basis for comparison.
D. Daugherty et al. / Commun Nonlinear Sci Numer Simulat 14 (2009) 2897–2908 Numerous generalizations of our models can be considered to further examine mood disorders using limit cycle oscilla- tors. For example, one can incorporate the fact that the medication administered to bipolar II individuals typically takes sev-eral days to reach therapeutic levels in the blood stream by adding a time-delay to the treatment function. One can also gainfurther mathematical and biological insight by directly examining the bifurcations that occur in our models. Moreover, onecan incorporate explicit time-dependence (for example, via delta functions representing instantaneous medication) into thetreatment function by considering a non-autonomous forcing function g(x, t).
Limit cycle oscillators can also be used to study behavioral cycling in mood disorders besides bipolar disorder, as such phenomena are of considerable interest to both psychologists and psychiatrists Pertinent disorders include unipolardepression and cyclothymia. Individuals suffering from the latter disorder rapidly flip back and forth between depressedand euphoric moods. Their high and low periods, which are relatively short, are much less intense than those of bipolar indi-viduals. Although it is not known precisely what leads individuals from one pole to the other, possible causes include psy-chosocial stress, disrupted sleep patterns, or some endogenous or internal biological mechanism that is not connected toenvironmental stimuli It is also not known what happens to cycling over time. It is speculated that an individual's cy-cling may slow down with age. From a mathematical perspective, we wish to highlight that the effects of disrupted sleeppatterns can be studied by coupling a cyclothymic (or bipolar) individual with a limit cycle oscillator representing the humansleep–wake cycle .
Another important issue to consider is the threshold for flipping poles. It has been theorized, for example, that the thresh- old for depression (and for flipping poles) gets lower over time. It has been speculated (by Robert Post , for example) thatthe longer one is ill, the more autonomous or disconnected episodes become from the environment. This is termed ‘‘kindlingtheory" and is believed to be analogous to what happens in epilepsy, in which there is a gradual kindling of biological dis-turbances in the brain that eventually surpasses the threshold for a seizure.
In this paper, we provide a mathematical framework for the modeling of bipolar disorder in terms of low-dimensional limit cycle oscillators. We propose and analyze two phenomenological models of bipolar individuals. Rather than focusingon the method of diagnosis, which is a difficult medical problem, we instead concentrate on the dynamics of our modelsunder the proposed treatment strategy, which involves a combination of drugs and therapy. The purpose of our modelingefforts is not to suggest a specific treatment for bipolar individuals but rather to provide some insight into the complicateddynamics of bipolar disorder. In this respect, we view our work as a first step in developing mathematical models of themood swings of bipolar individuals. Our intent is to provide a mathematical framework that ultimately leads to the devel-opment of more detailed models of bipolar disorder that incorporate clinical data. With this work, we hope to motivate thecollection of time-series data from clinical trials that will lead to refinements of our model that incorporate such data. In ourview, dynamical systems theory and mathematical modeling in general can lead to important advancements in the under-standing of bipolar disorder.
Valuable discussions with Michael Stubna, Carlos Castillo-Chavez, Richard Rand, Abdul-Aziz Yakubu, John Franke, Roxana Lopez-Cruz, Terry Kupers, Dane Quinn, Richard Frenette, Larry Riso, and Iain Macmillan are gratefully acknowledged. We alsothank Richard Rand for advice and discussions in reviving this article in 2008. Our research has been partially supported bygrants given by the National Science Foundation, National Security Agency, the Sloan Foundation (through the Cornell-SloanNational Pipeline Program in the Mathematical Sciences), and the NSF VIGRE program at Georgia Tech. Substantial financialand moral support was also provided by several groups at Cornell University: the Office of the Provost, the College of Agri-culture and Life Science (CALS), and the Department of Biological Statistics and Computational Biology.
[1] American Psychiatric Association. Diagnostic and statistical manual of mental disorders: DSM-IV-TR. Washington, DC: American Psychiatric Association; 2000.
[2] Bennett P, Segurado R, et al. The Wellcome Trust UK–Irish bipolar affective disorder sibling-pair genome screen: first stage report. Mol Psychiat [3] Camacho Erika T. Mathematical models of retinal dynamics. PhD thesis, Center for Applied Mathematics, Cornell University, Ithaca, NY; 2003.
[4] Camacho Erika T, Rand Richard H, Howland Howard H. Dynamics of two van der Pol oscillators coupled via a bath. Spec Boley Issue Int J Solid Struct [5] Colom Francesco, Vieta Eduard, et al. A randomized trial on the efficacy of group psychoeducation in the prophylaxis of recurrences in bipolar patients whose disease is in remission. Arch Gen Psychiat 2003;60(4):402–7.
[6] Daugherty Darryl, Roque-Urrea Tairi, Urrea-Roque John, Snyder Jessica, Wirkus Stephen, Porter Mason A. Mathematical models of bipolar disorder; 2003. Available from: <arXiv: nlin.CD/0311032>.
[7] Ehlers Cindy L. Chaos and complexity: can it help us to understand mood and behavior? Arch Gen Psych 1995;52:960–4.
[8] Fawcett Jan, Golden Bernard, Rosenfeld Nancy. New hope for people with bipolar disorder. Roseville, CA: Prima Publishing; 2000.
[9] Ferrier I Nicol, Macmillan Iain C, Young Allan H. The search for the wandering thymostat: a review of some developments in bipolar disorder research.
Brit J Psychiat 2001;178(41):s103–6.
[10] Ford-Martin Paula. Bipolar disorder. The Gale Encyclopedia of Alternative Medicine. In: Longe Jacqueline L, editor. Florence, KY: Gale Group; 2000.
D. Daugherty et al. / Commun Nonlinear Sci Numer Simulat 14 (2009) 2897–2908 [11] Frank Ellen, Swartz Holly A, Kupfer David J. Interpersonal and social rhythm therapy: managing the chaos of bipolar disorder. Biol Psychiat [12] Giacomini Hector, Neukirch Sébastian. Number of limit cycles of the Liénard equation. Phys Rev E 1997;56(4):3809–13.
[14] Gottschalk A, Bauer MS, Whybrow PC. Evidence of chaotic mood variation in bipolar disorder. Arch Gen Psychiat 1995;52:947–59.
[15] Griswold Kim. Management of bipolar disorder. Am Family Phys; 15 September 2000.
[17] Haslam David, Kennedy Sidney, Kusumakar Vivek, Kutcher Stan, Matte Raymond, Parikh Sagar, et al. A summary of clinical issues and treatment options. Bipolar disorder sub-committee Canadian network for mood and anxiety treatments (CANMAT). [18] Mobaji Hossein, Karahalios Karrie G. HCI applications for aiding children with mental disorders. Crossroads: ACM Stud Mag 2005;12(2):3.
[19] Anand Mohan M. Nonlinear dynamic system model of bipolar disorder. Proc IEEE 2007 [SoutheastCon].
[20] Judd Lewis L, Akiskal Hagop S, et al. A prospective investigation of the natural history of the long-term weekly symptomatic status of bipolar II disorder. Arch Gen Psychiat 2003;60(3):261–9.
[21] Kupers Terry. Private communication. Wright Institute; 29 June 2002.
[22] Lewis L, Vornik LA. Perceptions and impact of bipolar disorder: How far have we really come? Results of the National Depressive and Manic-Depressive Association 2000 survey of individuals with bipolar disorder. J Clin Psychiat 2003;64(2):161–74.
[23] Macmillan Iain C. Private communication. Consultant psychiatrist, early intervention service, Norfolk Mental Health Care NHS Trust; 25 May 2004.
[24] Nana Laurent. Bifurcation analysis of parametrically excited bipolar disorder model. Commun Nonlinear Sci Numer Simulat 2009;14:351–60.
[25] Post Robert M. Transduction of psychosocial stress into the neurobiology of recurrent affective disorder. Am J Psychiat 1992;149:999–1010.
[26] Post Robert M, Luckenbaugh David A. Unique design issues in clinical trials of patients with bipolar affective disorder. J Psychiat Res 2003;37(1):61–73.
[27] Rand Richard H. Topics in nonlinear dynamics with computer algebra. Computation in education: mathematics, science and engineering series. USA: Gordon and Breach Science Publishers; 1994.
[28] Riso Lawrence. Private communication. Department of Psychology, Georgia State University; 9 July 2003.
[29] Strogatz Steven H. Nonlinear Dynamics andChaos: With Applications in Physics, Biology, Chemistry, andEngineering. Reading, Mass: Addison-Wesley; [30] Strogatz Steven H. The mathematical structure of the human sleep–wake cycle. Lect Note Biomath, vol. 69. New York, NY: Springer-Verlag; 1986.
[31] Totterdell Peter, Briner Rob B, Parkinson Brian, Reynolds Shirley. Fingerprinting time-series: dynamic patterns in self-report and performance measures uncovered by a graphical nonlinear method. Brit J Psychol 1996;87:43–60.
[32] Wehr Thomas A, Goodwin Frederick K. Rapid cycling in manic-depressives induced by tricyclic antidepressants. Arch Gen Psychiat 1979;36:555–9.
[33] Weiss Roger D, Griffin Margaret L, et al. Group therapy for patients with bipolar disorder and substance dependence: results of a pilot study. J Clin [34] Wirkus Stephen A, Porter Mason A. Comment on bifurcation analysis of parametrically excited bipolar disorder model. Commun Nonlinear Sci Numer Simulat [in press].
[35] Wirkus Stephen A, Rand Richard H. The dynamics of two coupled van der Pol oscillators with delay coupling. Nonlinear Dyn 2002;30(3):205–21.

Source: http://eprints.maths.ox.ac.uk/1813/1/mason.pdf