{"id":{"repo_id":"duke","oai_identifier":"oai:dukespace.lib.duke.edu:10161/1320"},"canonical_url":"https://search.dev.ndltd.org/etd/duke/oai:dukespace.lib.duke.edu:10161/1320","repository":{"repo_id":"duke","name":"Duke University","base_url":"https://dukespace.lib.duke.edu/server/oai/request"},"display":{"title":"Bifurcations in the Echebarria-Karma Modulation Equation for Cardiac Alternans in One Dimension","abstract":"<p>While alternans in a single cardiac cell appears through a simple</p><p>period-doubling bifurcation, in extended tissue the exact nature</p><p>of the bifurcation is unclear. In particular, the phase of</p><p>alternans can exhibit wave-like spatial dependence, either</p><p>stationary or traveling, which is known as <italic>discordant</italic></p><p>alternans. We study these phenomena in simple cardiac models</p><p>through a modulation equation proposed by Echebarria-Karma. In</p><p>this dissertation, we perform bifurcation analysis for their</p><p>modulation equation.</p><p>Suppose we have a cardiac fiber of length l, which is</p><p>stimulated periodically at its x=0 end. When the pacing period</p><p>(basic cycle length) B is below some critical value B<sub>c</sub>,</p><p>alternans emerges along the cable. Let a(x,n) be the amplitude</p><p>of the alternans along the fiber corresponding to the n-th</p><p>stimulus. Echebarria and Karma suppose that a(x,n) varies</p><p>slowly in time and it can be regarded as a time-continuous</p><p>function a(x,t). They derive a weakly nonlinear modulation</p><p>equation for the evolution of a(x,t) under some approximation,</p><p>which after nondimensionization is as follows: </p><p> &partial<sub>t</sub> a = σ a + <bold>L</bold> a - g a <super>3</super>,</p><p>where the linear operator</p><p> <bold>L</bold> a = &partial<sub>xx</sub>a - &partial<sub>x</sub> a -Λ<super>-1</super> ∫ <super>0</super> <sub>x</sub> a(x',t)dx'.</p><p>In the equation, σ is dimensionless and proportional to</p><p>B<sub>c</sub> - B, i.e. σ indicates how rapid the pacing is,</p><p>Λ<super>-1</super> is related to the conduction velocity (CV) of the</p><p>propagation and the nonlinear term -ga<super>3</super> limits growth after the</p><p>onset of linear instability. No flux boundary conditions are</p><p>imposed on both ends.</p><p>The zero solution of their equation may lose stability, as the</p><p>pacing rate is increased. To study the bifurcation, we calculate</p><p>the spectrum of operator <bold>L</bold>. We find that the</p><p>bifurcation may be Hopf or steady-state. Which bifurcation occurs</p><p>first depends on parameters in the equation, and for one critical</p><p>case both modes bifurcate together at a degenerate (codimension 2)</p><p>bifurcation.</p><p>For parameters close to the degenerate case, we investigate the</p><p>competition between modes, both numerically and analytically. We</p><p>find that at sufficiently rapid pacing (but assuming a 1:1</p><p>response is maintained), steady patterns always emerge as the only</p><p>stable solution. However, in the parameter range where Hopf</p><p>bifurcation occurs first, the evolution from periodic solution</p><p>(just after the bifurcation) to the eventual standing wave</p><p>solution occurs through an interesting series of secondary</p><p>bifurcations.</p><p>We also find that for some extreme range of parameters, the</p><p>modulation equation also includes chaotic solutions. Chaotic waves</p><p>in recent years has been regarded to be closely related with</p><p>dreadful cardiac arrhythmia. Proceeding work illustrated some</p><p>chaotic phenomena in two- or three-dimensional space, for instance</p><p>spiral and scroll waves. We show the existence of chaotic waves in</p><p>one dimension by the Echebarria-Karma modulation equation for</p><p>cardiac alternans. This new discovery may provide a different</p><p>mechanism accounting for the instabilities in cardiac dynamics.</p>","abstract_html":"&lt;p&gt;While alternans in a single cardiac cell appears through a simple&lt;/p&gt;&lt;p&gt;period-doubling bifurcation, in extended tissue the exact nature&lt;/p&gt;&lt;p&gt;of the bifurcation is unclear. In particular, the phase of&lt;/p&gt;&lt;p&gt;alternans can exhibit wave-like spatial dependence, either&lt;/p&gt;&lt;p&gt;stationary or traveling, which is known as &lt;italic&gt;discordant&lt;/italic&gt;&lt;/p&gt;&lt;p&gt;alternans. We study these phenomena in simple cardiac models&lt;/p&gt;&lt;p&gt;through a modulation equation proposed by Echebarria-Karma. In&lt;/p&gt;&lt;p&gt;this dissertation, we perform bifurcation analysis for their&lt;/p&gt;&lt;p&gt;modulation equation.&lt;/p&gt;&lt;p&gt;Suppose we have a cardiac fiber of length l, which is&lt;/p&gt;&lt;p&gt;stimulated periodically at its x=0 end. When the pacing period&lt;/p&gt;&lt;p&gt;(basic cycle length) B is below some critical value B&lt;sub&gt;c&lt;/sub&gt;,&lt;/p&gt;&lt;p&gt;alternans emerges along the cable. Let a(x,n) be the amplitude&lt;/p&gt;&lt;p&gt;of the alternans along the fiber corresponding to the n-th&lt;/p&gt;&lt;p&gt;stimulus. Echebarria and Karma suppose that a(x,n) varies&lt;/p&gt;&lt;p&gt;slowly in time and it can be regarded as a time-continuous&lt;/p&gt;&lt;p&gt;function a(x,t). They derive a weakly nonlinear modulation&lt;/p&gt;&lt;p&gt;equation for the evolution of a(x,t) under some approximation,&lt;/p&gt;&lt;p&gt;which after nondimensionization is as follows: &lt;/p&gt;&lt;p&gt; &amp;partial&lt;sub&gt;t&lt;/sub&gt; a = σ a + &lt;bold&gt;L&lt;/bold&gt; a - g a &lt;super&gt;3&lt;/super&gt;,&lt;/p&gt;&lt;p&gt;where the linear operator&lt;/p&gt;&lt;p&gt; &lt;bold&gt;L&lt;/bold&gt; a = &amp;partial&lt;sub&gt;xx&lt;/sub&gt;a - &amp;partial&lt;sub&gt;x&lt;/sub&gt; a -Λ&lt;super&gt;-1&lt;/super&gt; ∫ &lt;super&gt;0&lt;/super&gt; &lt;sub&gt;x&lt;/sub&gt; a(x&#x27;,t)dx&#x27;.&lt;/p&gt;&lt;p&gt;In the equation, σ is dimensionless and proportional to&lt;/p&gt;&lt;p&gt;B&lt;sub&gt;c&lt;/sub&gt; - B, i.e. σ indicates how rapid the pacing is,&lt;/p&gt;&lt;p&gt;Λ&lt;super&gt;-1&lt;/super&gt; is related to the conduction velocity (CV) of the&lt;/p&gt;&lt;p&gt;propagation and the nonlinear term -ga&lt;super&gt;3&lt;/super&gt; limits growth after the&lt;/p&gt;&lt;p&gt;onset of linear instability. No flux boundary conditions are&lt;/p&gt;&lt;p&gt;imposed on both ends.&lt;/p&gt;&lt;p&gt;The zero solution of their equation may lose stability, as the&lt;/p&gt;&lt;p&gt;pacing rate is increased. To study the bifurcation, we calculate&lt;/p&gt;&lt;p&gt;the spectrum of operator &lt;bold&gt;L&lt;/bold&gt;. We find that the&lt;/p&gt;&lt;p&gt;bifurcation may be Hopf or steady-state. Which bifurcation occurs&lt;/p&gt;&lt;p&gt;first depends on parameters in the equation, and for one critical&lt;/p&gt;&lt;p&gt;case both modes bifurcate together at a degenerate (codimension 2)&lt;/p&gt;&lt;p&gt;bifurcation.&lt;/p&gt;&lt;p&gt;For parameters close to the degenerate case, we investigate the&lt;/p&gt;&lt;p&gt;competition between modes, both numerically and analytically. We&lt;/p&gt;&lt;p&gt;find that at sufficiently rapid pacing (but assuming a 1:1&lt;/p&gt;&lt;p&gt;response is maintained), steady patterns always emerge as the only&lt;/p&gt;&lt;p&gt;stable solution. However, in the parameter range where Hopf&lt;/p&gt;&lt;p&gt;bifurcation occurs first, the evolution from periodic solution&lt;/p&gt;&lt;p&gt;(just after the bifurcation) to the eventual standing wave&lt;/p&gt;&lt;p&gt;solution occurs through an interesting series of secondary&lt;/p&gt;&lt;p&gt;bifurcations.&lt;/p&gt;&lt;p&gt;We also find that for some extreme range of parameters, the&lt;/p&gt;&lt;p&gt;modulation equation also includes chaotic solutions. Chaotic waves&lt;/p&gt;&lt;p&gt;in recent years has been regarded to be closely related with&lt;/p&gt;&lt;p&gt;dreadful cardiac arrhythmia. Proceeding work illustrated some&lt;/p&gt;&lt;p&gt;chaotic phenomena in two- or three-dimensional space, for instance&lt;/p&gt;&lt;p&gt;spiral and scroll waves. We show the existence of chaotic waves in&lt;/p&gt;&lt;p&gt;one dimension by the Echebarria-Karma modulation equation for&lt;/p&gt;&lt;p&gt;cardiac alternans. This new discovery may provide a different&lt;/p&gt;&lt;p&gt;mechanism accounting for the instabilities in cardiac dynamics.&lt;/p&gt;","abstract_has_math":false,"creators":["Dai, Shu"],"institution":null,"degree_name":null,"degree_level":null,"degree_discipline":null,"degree_department":null,"school":null,"contributors":[],"advisors":["Schaeffer, David G."],"committee_chairs":[],"committee_members":[],"year":2009,"date_issued":"2009","date_published":"2009","updated_at":"2026-07-24T02:07:07Z","subjects":["Mathematics","Biology, Physiology","bifurcation","cardiac alternans","modulation equation"],"languages":["en_US"],"rights":[],"rights_urls":[],"identifier_entries":[]},"links":{"outbound_url":"https://hdl.handle.net/10161/1320","outbound_label":"Handle","outbound_source":"dc:identifier.uri"},"metadata_groups":[{"id":"people","label":"People","entries":[{"key":"dc:contributor.advisor","label":"Advisor","values":["Schaeffer, David G."]},{"key":"dc:creator","label":"Author","values":["Dai, Shu"]}]},{"id":"academic_context","label":"Academic Context","entries":[{"key":"dc:date.accessioned","label":"Dc Date Accessioned","values":["2009-08-27T18:33:54Z"]},{"key":"dc:date.available","label":"Dc Date Available","values":["2009-08-27T18:33:54Z"]},{"key":"dc:date.issued","label":"Date","values":["2009"]},{"key":"dc:type","label":"Dc Type","values":["Dissertation"]}]},{"id":"subjects_keywords","label":"Subjects and Keywords","entries":[{"key":"dc:subject","label":"Dc Subject","values":["Mathematics","Biology, Physiology","bifurcation","cardiac alternans","modulation equation"]}]},{"id":"language_rights","label":"Language and Rights","entries":[{"key":"dc:language.iso","label":"Language (ISO)","values":["en_US"]}]},{"id":"identifiers","label":"Identifiers","entries":[{"key":"dc:identifier.uri","label":"Identifier URI","values":["https://hdl.handle.net/10161/1320"]}]},{"id":"additional","label":"Additional Metadata","entries":[{"key":"dc:description.abstract","label":"Abstract","values":["<p>While alternans in a single cardiac cell appears through a simple</p><p>period-doubling bifurcation, in extended tissue the exact nature</p><p>of the bifurcation is unclear. In particular, the phase of</p><p>alternans can exhibit wave-like spatial dependence, either</p><p>stationary or traveling, which is known as <italic>discordant</italic></p><p>alternans. We study these phenomena in simple cardiac models</p><p>through a modulation equation proposed by Echebarria-Karma. In</p><p>this dissertation, we perform bifurcation analysis for their</p><p>modulation equation.</p><p>Suppose we have a cardiac fiber of length l, which is</p><p>stimulated periodically at its x=0 end. When the pacing period</p><p>(basic cycle length) B is below some critical value B<sub>c</sub>,</p><p>alternans emerges along the cable. Let a(x,n) be the amplitude</p><p>of the alternans along the fiber corresponding to the n-th</p><p>stimulus. Echebarria and Karma suppose that a(x,n) varies</p><p>slowly in time and it can be regarded as a time-continuous</p><p>function a(x,t). They derive a weakly nonlinear modulation</p><p>equation for the evolution of a(x,t) under some approximation,</p><p>which after nondimensionization is as follows: </p><p> &partial<sub>t</sub> a = σ a + <bold>L</bold> a - g a <super>3</super>,</p><p>where the linear operator</p><p> <bold>L</bold> a = &partial<sub>xx</sub>a - &partial<sub>x</sub> a -Λ<super>-1</super> ∫ <super>0</super> <sub>x</sub> a(x',t)dx'.</p><p>In the equation, σ is dimensionless and proportional to</p><p>B<sub>c</sub> - B, i.e. σ indicates how rapid the pacing is,</p><p>Λ<super>-1</super> is related to the conduction velocity (CV) of the</p><p>propagation and the nonlinear term -ga<super>3</super> limits growth after the</p><p>onset of linear instability. No flux boundary conditions are</p><p>imposed on both ends.</p><p>The zero solution of their equation may lose stability, as the</p><p>pacing rate is increased. To study the bifurcation, we calculate</p><p>the spectrum of operator <bold>L</bold>. We find that the</p><p>bifurcation may be Hopf or steady-state. Which bifurcation occurs</p><p>first depends on parameters in the equation, and for one critical</p><p>case both modes bifurcate together at a degenerate (codimension 2)</p><p>bifurcation.</p><p>For parameters close to the degenerate case, we investigate the</p><p>competition between modes, both numerically and analytically. We</p><p>find that at sufficiently rapid pacing (but assuming a 1:1</p><p>response is maintained), steady patterns always emerge as the only</p><p>stable solution. However, in the parameter range where Hopf</p><p>bifurcation occurs first, the evolution from periodic solution</p><p>(just after the bifurcation) to the eventual standing wave</p><p>solution occurs through an interesting series of secondary</p><p>bifurcations.</p><p>We also find that for some extreme range of parameters, the</p><p>modulation equation also includes chaotic solutions. Chaotic waves</p><p>in recent years has been regarded to be closely related with</p><p>dreadful cardiac arrhythmia. Proceeding work illustrated some</p><p>chaotic phenomena in two- or three-dimensional space, for instance</p><p>spiral and scroll waves. We show the existence of chaotic waves in</p><p>one dimension by the Echebarria-Karma modulation equation for</p><p>cardiac alternans. This new discovery may provide a different</p><p>mechanism accounting for the instabilities in cardiac dynamics.</p>"]},{"key":"dc:format.mimetype","label":"Dc Format Mimetype","values":["application/pdf"]},{"key":"dc:title","label":"Title","values":["Bifurcations in the Echebarria-Karma Modulation Equation for Cardiac Alternans in One Dimension"]}]}],"canonical_facts":{"dc:contributor.advisor":["Schaeffer, David G."],"dc:creator":["Dai, Shu"],"dc:date.accessioned":["2009-08-27T18:33:54Z"],"dc:date.available":["2009-08-27T18:33:54Z"],"dc:date.issued":["2009"],"dc:description.abstract":["<p>While alternans in a single cardiac cell appears through a simple</p><p>period-doubling bifurcation, in extended tissue the exact nature</p><p>of the bifurcation is unclear. In particular, the phase of</p><p>alternans can exhibit wave-like spatial dependence, either</p><p>stationary or traveling, which is known as <italic>discordant</italic></p><p>alternans. We study these phenomena in simple cardiac models</p><p>through a modulation equation proposed by Echebarria-Karma. In</p><p>this dissertation, we perform bifurcation analysis for their</p><p>modulation equation.</p><p>Suppose we have a cardiac fiber of length l, which is</p><p>stimulated periodically at its x=0 end. When the pacing period</p><p>(basic cycle length) B is below some critical value B<sub>c</sub>,</p><p>alternans emerges along the cable. Let a(x,n) be the amplitude</p><p>of the alternans along the fiber corresponding to the n-th</p><p>stimulus. Echebarria and Karma suppose that a(x,n) varies</p><p>slowly in time and it can be regarded as a time-continuous</p><p>function a(x,t). They derive a weakly nonlinear modulation</p><p>equation for the evolution of a(x,t) under some approximation,</p><p>which after nondimensionization is as follows: </p><p> &partial<sub>t</sub> a = σ a + <bold>L</bold> a - g a <super>3</super>,</p><p>where the linear operator</p><p> <bold>L</bold> a = &partial<sub>xx</sub>a - &partial<sub>x</sub> a -Λ<super>-1</super> ∫ <super>0</super> <sub>x</sub> a(x',t)dx'.</p><p>In the equation, σ is dimensionless and proportional to</p><p>B<sub>c</sub> - B, i.e. σ indicates how rapid the pacing is,</p><p>Λ<super>-1</super> is related to the conduction velocity (CV) of the</p><p>propagation and the nonlinear term -ga<super>3</super> limits growth after the</p><p>onset of linear instability. No flux boundary conditions are</p><p>imposed on both ends.</p><p>The zero solution of their equation may lose stability, as the</p><p>pacing rate is increased. To study the bifurcation, we calculate</p><p>the spectrum of operator <bold>L</bold>. We find that the</p><p>bifurcation may be Hopf or steady-state. Which bifurcation occurs</p><p>first depends on parameters in the equation, and for one critical</p><p>case both modes bifurcate together at a degenerate (codimension 2)</p><p>bifurcation.</p><p>For parameters close to the degenerate case, we investigate the</p><p>competition between modes, both numerically and analytically. We</p><p>find that at sufficiently rapid pacing (but assuming a 1:1</p><p>response is maintained), steady patterns always emerge as the only</p><p>stable solution. However, in the parameter range where Hopf</p><p>bifurcation occurs first, the evolution from periodic solution</p><p>(just after the bifurcation) to the eventual standing wave</p><p>solution occurs through an interesting series of secondary</p><p>bifurcations.</p><p>We also find that for some extreme range of parameters, the</p><p>modulation equation also includes chaotic solutions. Chaotic waves</p><p>in recent years has been regarded to be closely related with</p><p>dreadful cardiac arrhythmia. Proceeding work illustrated some</p><p>chaotic phenomena in two- or three-dimensional space, for instance</p><p>spiral and scroll waves. We show the existence of chaotic waves in</p><p>one dimension by the Echebarria-Karma modulation equation for</p><p>cardiac alternans. This new discovery may provide a different</p><p>mechanism accounting for the instabilities in cardiac dynamics.</p>"],"dc:format.mimetype":["application/pdf"],"dc:identifier.uri":["https://hdl.handle.net/10161/1320"],"dc:language.iso":["en_US"],"dc:subject":["Mathematics","Biology, Physiology","bifurcation","cardiac alternans","modulation equation"],"dc:title":["Bifurcations in the Echebarria-Karma Modulation Equation for Cardiac Alternans in One Dimension"],"dc:type":["Dissertation"]},"updated_at":"2026-07-24T02:07:07Z"}