© 1998 by European Society of Cardiology
Copyright © 1998, European Society of Cardiology
Chaos in the cardiovascular system: an update
Institut für Physiologie der Humboldt-Universität zu Berlin (Charité), Tucholskystraße 2, D-10117 Berlin, Germany
* Corresponding author.
Received 14 May 1998; accepted 31 July 1998
| Abstract |
|---|
|
|
|---|
Rhythmic changes of blood pressure, heart rate, and other cardiovascular measures have drawn the attention of several investigators, since these oscillations can shed light onto the activity of the underlying control network. The overwhelming proportion of circulatory variations, however, are not linear, i.e., they do not consist of perfectly rhythmic components. Thus, these fluctuations are more adequately analysed by non-linear techniques, most of which are adopted from chaos theory. A spotlight issue of Cardiovascular Research (Vol. 31, 1996), focused on chaos in the cardiovascular system. This current review outlines today's understanding of this field by presenting the major discoveries and developments which have taken place since then.
KEYWORDS Heart rate variability; Non-linear dynamics; Risk stratification; Vasomotion
| 1 Introduction |
|---|
|
|
|---|
Two years ago the Spotlight Issue on Chaos in the Cardiovascular System appeared [Cardiovascular Research, volume 31, 1996], in which a selection of topics related to chaos in heart and circulation were highlighted. Since then, new studies have deepened our insight and broadened our understanding regarding the nature of chaotic behaviour.
The purpose of this article is to update those topics which have been affected most by the newest trends in this field. Hence, the recent developments of cardiac chaos were chosen as a major focus. Another emphasis was placed on mechanisms responsible for chaotic behaviour observed in vasomotion and risk stratification. A brief outline of the methods used to identify chaos is found at the end of this review.
| 2 Cardiac chaos |
|---|
|
|
|---|
Heart rate (HR) is not constant but varies considerably, even in the absence of physical or mental stress. A power spectrum of HR reveals at least two frequency ranges [47, 53], in which power can accumulate. It is discussed controversially, whether modification of the power of these ranges occurs by influencing either sympathetic [low frequency range, (LF) 0.04–0.15 Hz] or vagal [high frequency range (HF) 0.15–0.4 Hz] traffic to the heart [10, 70]. Beside these more or less periodic components, HR reveals a broad, noise-like variability over a large frequency span [36]. It seems that this irregular variability, which accounts for the largest proportion of HRV, is due to non-linearities in the control network.
Several authors have quantified non-linear measures in order to test their feasibility to identify changes in autonomic nervous outflow. Yeragani et al. compared fractal dimensions of HR time series with HF and LF powers, and found positive correlations of these measures with sympathetic and vagal activity [76]. Others found that non-linear components of HRV are drastically reduced by cholinergic blockade, but not by β-adrenergic blockade [86, 87], a similar observation has been made after vagal denervation in neonatal swine [40], and in patients suffering from congestive heart failure [6]. In another study, Hagerman and colleagues disrupted the autonomic nervous activity to the heart with propranolol and atropine and found a reduction in the largest Lyapunov exponent [26]. This confirms the potential value of using measures of non-linear dynamics as a tool for evaluating autonomic nervous outflow to the heart. Interestingly, Hagerman et al. could not totally eliminate the non-linear structures in HRV in their experiments using combined blockade. Thus they concluded that other mechanisms like circulating hormones, preload, or afterload contribute to the non-linearity observed. However, respiration does not seem to contribute to this effect, since no evidence was found for a non-linear input from the respiration into the cardiovascular centres [32].
It is intriguing, that marked changes in non-linear properties are observed throughout early human development [43, 51]. This seems to reflect changing control complexity during autonomic nervous system maturation [40]. A distinct non-linearity was not observed in preterm infants [69], which indicates that complex heart rhythms mainly develop during the final weeks of gestation [50]. Attempts have been made to use non-linear measures to identify children at risk, however, with differing results. Largest Lyapunov exponents of infants suffering from bronchopulmonary dysplasia, for instance, showed no differences compared with a control group [44, 51]. On the other hand, a symbolic dynamics approach applied to premature babies at high risk for Sudden Infant Death Syndrome revealed significant differences to controls [45].
In healthy adults, non-linear dynamics in HR seems to represent the normal situation [28, 42, 56]. These are obscured under pathological circumstances, e.g., such in which the cardiovascular–respiratory interplay is impaired [30], which can also be observed in patients after heart transplantation [25]. In line with the view, that the healthy state is characterised by a certain degree of chaos, it has been found that abnormalities in autonomic system functions diminish cardiac chaos. This was recently demonstrated for congestive heart failure (CHF). In CHF patients, Poon and Merrill [57]observed an increase in non-chaotic HR fluctuations, which alternate with chaotic short-term variations. The reduction in HRV complexity in CHF went parallel to a decrease in parasympathetic nervous system activity [6, 48, 72]. This supports the above mentioned hypothesis, that part of HRV complexity arises from cardiac vagal activity. Further support of chaos being inherent to the physiological control of heart rate is provided by Narayanan et al. [46]. They discovered unstable periodic orbits (UPO) in RR time series. Periodic orbits reflect the deterministic dynamics of the underlying system, but the UPOs are unstable, which means that they have a positive Lyapunov exponent and are therefore typical for chaotic systems. (A similar approach, however, was used by Le Pape et al. [38], who did not identify UPOs in mice electrocardiograms).
Many authors have addressed the question, whether ventricular fibrillation is chaotic or just some form of noise [16, 33, 81]. It is known by high-resolution mappings, that re-entrant electrical arrhythmias are a cause for the transition from orderly beating to fibrillation. Ventricular fibrillation consists of propagating wavefronts which traverse the myocardium. As the high-resolution mappings have revealed, there exist local spatio-temporal correlations, that become weaker correlated during fibrillation [82]. This observation may indicate that a deterministic mechanism generates fibrillation [3, 17, 81]. Although the underlying processes are not completely random, it remains unclear whether they represent a form of deterministic chaos. Garfinkel et al. recently studied the spatio-temporal dynamics in thin sheets of human and canine ventricular tissue [15]. The authors induced re-entry by application of the potassium channel opener cromakalim. The resulting period-1 re-entrant tachycardia spontaneously changed into period-2 behaviour. This is an example for a so called period doubling and constitutes a major route to chaos [52](see also Fig. 1). In another experimental setup by this group, this transition was achieved by an external pacing stimulus. In both cases, the period-2 behaviour became amplitude and period modulated, and finally switched to a highly aperiodic fibrillation-like behaviour.
|
The transition mechanisms that these authors observed [15], was from periodicity (tachycardia-flutter) to chaos, (fibrillation) via a quasi-periodic route. To further elucidate the transition pattern, as well as to determine the role of ventricular tissue mass in fibrillation, Kim et al. successively reduced tissue mass by sequential cutting [34]. They found that the number of wave fronts decreased with decreasing tissue mass and that there is a critical tissue mass essential for maintenance of fibrillation. Parallel to the reduction of tissue mass, the Kolmogorov entropy, a measure of dynamical complexity, also decreased [34]. With these experiments, the authors were able to resolve more distinctly the region between periodicity (regular arrhythmias) and chaos (fibrillation). Together with the study of Garfinkel et al. [15], these experiments strongly support the hypothesis that fibrillation is indeed chaos, reached by the perioddoubling route [67]. Taken together, there are studies showing that the healthy cardiovascular system is on the verge to chaos. Circulatory disorders can be manifested by a decrease in chaos (as in CHF), or by a greater magnitude of this form of non-linear dynamics (e.g., fibrillation).
| 3 Vascular chaos |
|---|
|
|
|---|
Rhythmic changes in diameter (or resistance [64]) of blood vessels is a commonly observed physiological phenomenon, which occurs in small arteries and arterioles [14]. Although a great amount of resistance variability is accounted for by almost periodic changes, Yamashiro et al. suggested a relaxation oscillator (pacemaker) behaviour responsible for vasomotion [84]. Indeed, the dynamics of vasomotion may be chaotic. Tsuda and colleagues investigated pulsations in a finger's capillary vessels [71]and found the largest and the second largest Lyapunov exponents to be positive and, hence, the system was classified as chaotic. In order to rule out the possibility that chaos in vasomotion might simply be a secondary effect (i.e., the chaotic motion observed in the capillary vessels originates from RR-interval variability (cardiac chaos), the authors made simultaneous measurements of heart beat. In fact, heart rate did reveal deterministic chaos, but its topology differed from that of the capillaries. Since the peripheral vessels are under control of the autonomic nervous system, Tsuda et al. interpreted their results in terms of a non-linear modulation by the nervous system [71]. Other studies, on isolated mesenteric and femoral rat arteries [66]and hindlimb skeletal muscle of anaesthetized rabbits [59], however, indicate that the local pacemaker cells or paracrine factors are responsible for inducing vasomotion.
Thus, Griffith and Edwards investigated vasomotion in isolated rabbit ear resistance arteries, which were elicited by application of histamine [19]. They described an irregular rhythmicity. The transition from regular (or periodic) behaviour to chaos took place via three possible ways: first, a period-doubling route to chaos; second, a quasi-periodic route to chaos; and third, intermittency, a state, where the system is arrested between almost periodic motion and chaos [24]. These chaotic fluctuations were not abolished by endothelium denudation, suggesting vascular smooth muscle cells as the source of the oscillations (with correlation dimensions between 2 and 4).
This interpretation was supported by a further series of experiments, in which it was shown that vascular smooth muscle Ca2+ fluxes can create aperiodic vasomotion [20, 22]. Two major frequency components in the perfusion pressure spectrum were selectively blocked pharmacologically with calcium antagonists: a fast subsystem (5–20 s) was inhibited by verapamil (inhibiting voltage-dependent Ca2+ influx; inhibition of voltage-dependent K+ channels or the Ca2+–ATPase extrusion pump, however, had no effect on the form and complexity of vasomotion [13]), whereas a slow subsystem (1–5 min) was attenuated by ryanodine, suppressing Ca2+-induced Ca2+ release from the sarcoplasmic reticulum [35]. Thus, a plausible reason for this form of vascular chaos is the coupling between both oscillators [23], i.e., cytosolic [Ca2+]. It was finally shown, that this form of vasomotion (chaos, quasi-periodicity, and mixed-mode responses) can be described in terms of iterative circle maps; a model for the dynamics of coupled non-linear oscillators [9].
The aperiodic oscillations observed in the perfused rabbit ear arteries has successfully been controlled. As mentioned in the initial spotlight issue [65], control of chaos is the external maintenance of the system's dynamics within a physiological range. In this study, control was attempted by a negative feedback loop that regulated pump speed according to the difference between instantaneous perfusion pressure and a reference value [21]. Different absolute levels of perfusion pressure or flow converted irregular pressure fluctuations to periodic dynamics (i.e., stabilization of unstable period-1, period-2, and period-4 orbits), and rhythmic activity was almost completely suppressed (stabilization of an unstable steady-state). However, in approximately 40% of the cases, the dynamics remained highly irregular during external control, providing evidence that absolute levels of flow/perfusion pressure are not the sources of the intrinsic complexity observed in vasomotion. Blood pressure appears rather to play the role of a control, or bifurcation parameter, in vasomotion. Autoregulation, for instance, seems to be a mode, where aperiodic or chaotic behaviour occurs, whereas periodic haemodynamics is activated below a specific local arterial blood pressure and flow threshold (at the lower end of autoregulation) [59, 60]. This characteristic was recently found in a mathematical model of an arteriolar network, where diameter oscillations shifted between periodic, quasi-periodic, and chaotic behaviour, depending on arterial pressure level [7].
How do organisms benefit from the chaotic behaviour of vasomotion? The sensitivity of chaotic systems to initial conditions (i.e., to perturbations) allows large changes in state with minimum expenditure of energy to be brought about [61]. This cost-benefit analysis gains significance in situations, where systems frequently skip from one state to another [49]. On the other hand, perturbations can easily destabilize (unwanted) periodic or steady-state behaviour and thereby induce, rather than suppress, chaos. If the system is allowed to operate near the frontier from periodic behaviour to chaotic behaviour, only minute alterations in a single parameter are necessary to switch the system from regular to irregular behaviour [24, 85], thus providing a high degree of flexibility to the control of microcirculatory blood flow. Thus, non-linear dynamics observed in systemic blood pressure [2, 77–79]may be caused partly by non-linear vasomotion.
| 4 Risk stratification and non-linear dynamics |
|---|
|
|
|---|
Can methods from non-linear dynamics help improve the reliability of risk stratification? This potential of non-linear HRV analysis has been demonstrated in a number of reports, focusing on risk detection and stratification in situations related to myocardial infarction (MI).
A bridge connecting linear methods in time and frequency domain with non-linear dynamics are scatter plots, or Lorenz plots, in which one RR interval is plotted against its predecessor, i.e., RRi vs. RRi–1 (see Fig. 2). The same method is used to reconstruct an attractor in two dimensions (see Fig. 1). Copie and co-workers [8]showed that indices of HRV were strongly correlated with the shape of the scatter plot. The scatter plot length is correlated with long-term variability parameters, scatter plot width was closely related to short-term variability. Since scatter plots provide a geometrical representation of beat-to-beat dynamics, including the non-linear features, they may provide a tool for examining connections between physiological states and non-linear markers [8]. Hnatkova and co-workers [27]used scatter plots to distinguish between postinfarct patients who suffered an arrhythmic event during a follow-up period and those without arrhythmic complications. In this study, scatter plots were shown to provide a larger positive predictive accuracy than linear variability measures (i.e., standard deviation of all RR intervals).
|
Mäkikallio et al. [41]obtained a similar result in a retrospective analysis of 24-h HR data from patients with a history of myocardial infarction: the lengths of the scatter plots were significantly reduced in those patients with a propensity for ventricular tachycardia (VT). In agreement, a previous study from this group clearly demonstrated that analyzing the shape of the scatter plot pattern is a potent tool for the prediction of forthcoming VT [29]. Prior to the onset of VT, the scatter plot broadened and simultaneously, the length tended to decrease. Since the width and length of the scatter plot relates directly to short- and long-term HRV, this technique may provide insight into HR dynamics, which are not readily detected with conventional techniques (e.g., spectrum analysis).
Skinner and colleagues [62]applied the point correlation dimension algorithm (PD2), a modification of the correlation dimension, to HRV of conscious pigs during experimental MI. After coronary artery occlusion, they observed a significant decline in correlation dimension, which further decreased minutes prior to the onset of VT. Intriguingly, the second dip in PD2 was not observed in those animals, which did not experience VT. In a further study, Skinner [63]demonstrated that the mentioned decrease of PD2 prior to life-threatening VT also occurs in patients with pre-existing heart disease. In randomly selected control groups, who had no previous history of VF, but who also had non-sustained ventricular tachycardia, the PD2 values did not fall to such low levels.
Voss et al. [73, 74]introduced two new methods from non-linear dynamics into HRV classification in patients after MI: symbolic dynamics and renormalized entropy. In brief, symbolic dynamics facilitates the analysis of dynamic aspects of HRV. The RR-interval series are transformed into symbol sequences, with symbols from a given symbolic alphabet (here 0, 1, 2, 3) from which words were created with a word length of 3. Further analysis was done using the distributions of these words. The renormalized entropy, on the other hand, is a rather complex method derived from thermodynamics, which can measure the complexity of any state in relation to a reference state (for a derivation of renormalized entropy and symbolic dynamics and the definition of symbols, please refer to [74]). When methods from the three parameter domains (time, frequency, and non-linear domains) were applied individually, none of these lead to a satisfactory discrimination between patients with and without life-threatening arrhythmias. However, the non-linear methods were superior to the other means. The main result of their study was that a combined set of parameters from all three domains significantly improved the precision of high risk stratification.
The methods used for the different domains describe various aspects or characteristics of HRV [31](e.g., a high or low HRV, periodic processes due to an autonomic regulation, or the complexity of the RR-interval series). Clearly, it would be beneficial to have a set of a few parameters covering all domains. In a recent publication, Voss and co-workers [75]introduced a set of four parameters (namely Shannon entropy, relative power in the low and very low frequency range, words from symbolic dynamics consisting of symbols 0 and 2 only, and the mean value of RR intervals from a 30-min stationary segment) and demonstrated that the multi-parametric analysis of HRV significantly increased the prediction of risk in patients who survived a MI. These promising studies indicate the potential of non-linear methods, and their combination with other tools, for identifying patients with a high cardiovascular risk.
| 5 The analysis of chaotic data |
|---|
|
|
|---|
In this section we outline some of the most widely used techniques. The algorithms developed for the analysis of non-linear time series have been designed mainly during the past twenty years, and the majority of the present publications still make use of these approaches.
In several cases of time series analysis, only one quantity, or signal, is obtained during a recording session (e.g., an RR or blood pressure series). The goal is to recognize the major characteristics of the system by this data set. This can be a straightforward task, provided that a valid mathematical model is available. Unfortunately, this is rarely the case in cardiovascular research. When we deal with experiments in which we cannot record all the variables, the total number of degrees of freedom may not be known exactly (mathematically speaking, the number of degrees of freedom, or the dimension of the state space, is the same as the number of first-order differential equations necessary to describe the system.) Thus, prior to the investigation of characteristic quantities of the system, like fractal dimensions, entropies, or Lyapunov exponents, the phase space of the underlying dynamical system has to be reconstructed from the observed variable. This is often done by plotting the time series vs. a time delayed sequence of itself. For the choice of correct exact delay times and sufficiently high embedding dimensions of the phase space, please refer to [1, 11].
The definition of dimension in daily life, i.e., for Euclidian space, is clear: a point has the dimension 0, a line or curve has dimension 1, a surface 2, and a volume 3. Mathematically speaking, a dimension gives a statistical measure of the geometry of a cloud of points and can be assigned to any arbitrary set of data points. There are several concepts of estimating dimensions, but in general they all yield common values, to which we refer to as fractal dimension. Unlike the common definition of dimension, strange attractors originating from chaotic systems may also reveal non-integer (fractal) dimensions. The standard algorithms for estimating the fractal dimension require very large data sets (>10 000) and more importantly, they presuppose stationarity of the time series. To overcome this limitation, which would restrict possible physiological and cardiovascular application, algorithms have been developed, which operate with short and/or instationary data (for example, the pointwise correlation dimension [37, 62, 63]).
Lyapunov exponents measure the average local rate of divergence of neighboring trajectories in phase-space embeddings [11]. A positive Lyapunov exponent can be taken as a definition of chaos, provided the system is known to be deterministic (Fig. 1G). An early algorithm developed by Wolf et al. directly uses the rate of divergence between two nearby trajectories [83]. This algorithm works well with large data sets devoid of noise contamination, but may fail in situations of short and/or noisy data sets. Several algorithms have been developed to overcome this limitation [4, 5, 12, 58]and may be superior for the application to cardiovascular time series.
Entropy refers to system randomness and predictability and allows systems to be classified by rate of information loss or generation. In other words, entropy specifies the amount of information one learns about the source of the signal x by a measurement of x (this idea was formulated by Kolmogorov and Sinai 1958–59 and is now known as KS entropy). Since the KS entropy is the sum of positive Lyapunov exponents, a non-linear system with positive entropy is intrinsically unstable. In case of large data sets, KS entropy can be estimated directly from correlation integrals [18], an alternative approach is approximate entropy [54, 55], suitable for short data sets.
As was mentioned above, the entropy of a signal is related to the concept of forecasting the future evolution of the time series. Non-linear forecasting is an approach developed for making short-term predictions about the trajectories of dynamical systems. It usually compares predictions of the time series' evolution, which were assembled from another part of the time series, with the actual evolution. The prediction error is then the normalized root-mean-square difference between predicted and actual evolution [1, 39, 68, 80], which has been utilized to characterize the nature of the underlying dynamics (periodic, deterministic chaotic, or signals contaminated with externally induced noise).
| 6 Conclusions and perspectives |
|---|
|
|
|---|
Methods derived from the field on non-linear dynamics, or chaos theory, have provided useful tools to characterize systems from life sciences. Applications can be found throughout a wide scope of medicine, such as cardiology, physiology, and the neurosciences.
There are limitations to linear and non-linear methods, and there are several situations in which the one method is superior to the other. Linear methods are appropriate for stable linear systems, e.g., those which are characterized by a few distinct peaks in a power spectrum. Accordingly, it is not surprising, that the application of linear methods can fail in case of non-linear systems. The time series of the latter often reveal broad-band spectra instead of pronounced peaks, indicating that their dynamics stem from multivariate sources. Such dynamics are typical for most biological systems in general, and for the cardiovascular control network in particular.
However, both linear and non-linear approaches can complement each other in the analysis of a system's behaviour. Diameter changes in a small artery, for instance, are interpreted as noise by the broad-band power spectrum. The application of non-linear methods, however, shows that this process obeys deterministic rules. HRV is also, to a certain degree, describable by spectrum analysis. The aperiodic processes in HRV, which do not appear as distinct peaks, are only adequately addressed by non-linear techniques. Regarding HRV, there are indications that a reduction in complexity comes along with a decrease in vagal traffic to the heart, suggesting that a considerable amount of non-linear behaviour is promoted by this branch of the autonomic nervous system.
Finally, studies regarding the pathogenesis of myocardial infarction are promising with regard to the future prospects of non-linear dynamics in clinical use. Especially, high-risk stratification after myocardial infarction has demonstrated the advantage of combining methods from the time and frequency domain.
Time for primary review 42 days.
| References |
|---|
|
|
|---|
- Abarbanel H.D.I, Brown R, Sidorowich J.J, Tsimring L.S. The analysis of observed chaotic data in physical systems. Rev Mod Phys (1993) 65:1331–1392.[CrossRef][Web of Science]
- Almog Y, Eliash S, Oz O, Akselrod S. Non-linear analysis of BP signal. Can it detect malfunctions in BP control? Am J Physiol (1996) 271:H396–H403.[Medline]
- Bayly P.V, Johnson E.E, Wolf P.D, et al. A quantitative measurement of spatial order in ventricular fibrillation. J Cardiovasc Electrophysiol (1993) 4:533–546.[Web of Science][Medline]
- Brown R, Bryant P, Abarbanel H.D.I. Computing the Lyapunov spectrum of a dynamical system from an observed time series. Phys Rev A (1991) 43:2787–2806.[CrossRef][Medline]
- Brown R. Calculating Lyapunov exponents for short and/or noisy data sets. Phys Rev E (1993) 47:3962–3969.[CrossRef]
- Butler G.C, Ando S.-I, Floras J.S. Fractal component of variability of heart rate and systolic blood pressure in congestive heart failure. Clin Sci (1997) 92:543–550.[Web of Science][Medline]
- Cavalcanti S, Ursino M. Chaotic oscillations in microvessel arterial networks. Ann Biomed Eng (1996) 24:37–47.[Web of Science][Medline]
- Copie X, Le Heuzey J.-Y, Iliou M.-C, et al. Correlation between time domain measures of heart rate variability and scatter plots in postinfarction patients. PACE (1996) 19:342–347.[Medline]
- de Brouwer S, Edwards D.H, Griffith T.M. Simplification of the quasiperiodic route to chaos in agonist-induced vasomotion by iterative circle maps. Am J Physiol (1998) 274:H1315–H1326.[Web of Science][Medline]
- Eckberg D.L. Sympathovagal balance: a critical appraisal. Circulation (1997) 96:3224–3232.
[Free Full Text] - Eckmann J.-P, Ruelle D. Ergodic theory of chaos and strange attractors. Rev Mod Phys (1985) 57:617–656.[CrossRef][Web of Science]
- Eckmann J.-P, Kamphorst S.O, Ruelle D, Ciliberto S. Lyapunov exponents from time series. Phys Rev A (1986) 34:4971–4979.[CrossRef][Medline]
- Edwards D.H, Griffith T.M. Entrained ion-transport systems generate the membrane component of chaotic agonist-induced vasomotion. Am J Physiol (1997) 273:H909–H920.[Web of Science][Medline]
- Fujii K, Heistad D.D, Maraci F.M. Vasomotion of basilar arteries in vivo. Am J Physiol (1990) 258:H1829–H1834.[Web of Science][Medline]
- Garfinkel A, Chen P.-S, Walter D.O, et al. Quasiperiodicity and chaos in cardiac fibrillation. J Clin Invest (1997) 99:305–314.[Web of Science][Medline]
- Goldberger A.L, Bhargava V, West B.J, Mandell A.J. Some observations on the question: is ventricular fibrillation chaos? Physica (1986) 19D:282–289.[CrossRef]
- Gray R.A, Pertsov A.M, Jalife J. Spatial and temporal organization during cardiac fibrillation. Nature (1998) 392:75–78.[CrossRef][Medline]
- Grassberger P, Procaccia I. Estimation of the Kolmogorov entropy from a chaotic signal. Phys Rev A (1983) 28:2591–2593.[CrossRef]
- Griffith T.M, Edwards D.H. EDRF suppresses chaotic pressure oscillations in isolated resistance artery without influencing intrinsic complexity. Am J Physiol (1994) 266:H1786–H1800.[Web of Science][Medline]
- Griffith T.M, Edwards D.H. Fractal analysis of the role of smooth muscle Ca2+ fluxes in genesis of chaotic arterial pressure oscillations. Am J Physiol (1994) 266:H1801–H1811.[Web of Science][Medline]
- Griffith T.M, Edwards D.H. Complexity of chaotic vasomotion is insensitive to flow and pressure but can be regulated by external control. Am J Physiol (1995) 269:H656–H668.[Web of Science][Medline]
- Griffith T.M. Temporal chaos in the microcirculation. Cardiovasc Res (1996) 31:342–358.
[Abstract/Free Full Text] - Griffith T.M, Edwards D.H. Ca2+ sequestration as a determinant of chaos and mixed-mode dynamics in agonist-induced vasomotion. Am J Physiol (1997) 272:H1696–H1709.[Web of Science][Medline]
- Griffith T.M, Parthimos P, Crombie J, Edwards D.H. Critical scaling and type-III intermittent chaos in isolated rabbit resistance arteries. Phys Rev E (1997) 56:R6287–R6290.[CrossRef]
- Guzzetti S, Signorini M.G, Cogliati C, et al. Non-linear dynamics and chaotic indices in heart rate variability of normal subjects and heart-transplanted patients. Cardiovasc Res (1996) 31:441–446.
[Abstract/Free Full Text] - Hagerman I, Berglund M, Lorin M, Nowak J, Sylvén C. Chaos-related deterministic regulation of heart rate variability in time- and frequency domains: effects of autonomic blockade and exercise. Cardiovasc Res (1996) 31:410–418.
[Abstract/Free Full Text] - Hnatkova K, Copie X, Staunton A, Malik M. Numeric processing of Lorenz plots of R–R intervals from long-term ECGs: Comparison with time-domain measures of heart rate variability for risk stratification after myocardial infarction. J Electrocardiol (1995) 28(Suppl):74–80.[Web of Science][Medline]
- Hoyer D, Schmidt K, Zwiener U, Bauer R. Characterization of complex heart rate dynamics and their pharmacological disorders by non-linear prediction and special data transformations. Cardiovasc Res (1996) 31:434–440.
[Abstract/Free Full Text] - Huikuri H.V, Seppänen T, Koistinen M.J, et al. Abnormalities in beat-to-beat dynamics of heart rate before the spontaneous onset of life-threatening ventricular tachyarrhythmias in patients with prior myocardial infarction. Circulation (1996) 93:1836–1844.
[Abstract/Free Full Text] - Ivanov P.C, Rosenblum M.G, Peng C.-K, et al. Scaling behaviour of heartbeat intervals obtained by wavelet-based time-series analysis. Nature (1996) 383:323–327.[CrossRef][Medline]
- Kanters J.K, Højgaard M.V, Agner E, Holstein-Rathlou N.H. Short- and long-term variations in non-linear dynamics of heart rate variability. Cardiovasc Res (1996) 31:400–409.
[Abstract/Free Full Text] - Kanters J.K, van Højgaard M.V, Agner E, Holstein-Rathlou N.-H. Influence of forced respiration on non-linear dynamics in heart rate variability. Am J Physiol (1997) 272:R1149–R1154.[Web of Science][Medline]
- Kaplan D.T, Cohen R.J. Is fibrillation chaos? Circ Res (1990) 67:886–892.
[Abstract/Free Full Text] - Kim Y.-H, Garfinkel A, Ikeda T, et al. Spatiotemporal complexity of ventricular fibrillation revealed by tissue mass reduction in isolated swine right ventricle. J Clin Invest (1997) 100:2486–2500.[Web of Science][Medline]
- Knot HJ, Standen NB, Nelson MT. Ryanodine receptors regulate arterial diameter and wall [Ca2+] in cerebral arteries of rat via Ca2+-dependent K+ channels. J Physiol (Lond) 1998;508.1:211–221.
- Kobayashi M, Musha T. 1/f fluctuations of heartbeat period. IEEE Trans Biomed Eng (1982) 29:456–457.[Web of Science][Medline]
- Kroll M.W, Fulton K.W. Slope filtered pointwise correlation dimension algorithm and its evaluation with prefibrillation heart rate data. J Electrocardiol (1991) 24(Suppl):97–101.[CrossRef][Web of Science]
- Le Pape G, Giacomini H, Swynghedauw B, Mansier P. A statistical analysis of sequences of cardiac interbeat intervals does not support the chaos hypothesis. J Theor Biol (1997) 184:123–131.[CrossRef][Medline]
- Linsay P.S. An efficient method of forecasting chaotic time series using linear interpolation. Phys Lett A (1991) 153:353–356.[CrossRef]
- Lipsitz L.A, Pincus S.M, Morin R.J, et al. Preliminary evidence for the evolution in complexity of heart rate dynamics during autonomic maturation in neonatal swine. J Auton Nerv Syst (1997) 65:1–9.[CrossRef][Web of Science][Medline]
- Mäkikallio T.H, Seppänen T, Airaksinen K.E.J, et al. Dynamic analysis of heart rate may predict subsequent ventricular tachycardia after myocardial infarction. Am J Cardiol (1997) 80:783–797.[CrossRef][Web of Science][Medline]
- Mansier P, Clairambault J, Charlotte N, et al. Linear and non-linear analyses of heart rate variability: a minireview. Cardiovasc Res (1996) 31:371–379.
[Abstract/Free Full Text] - Mrowka R, Patzak A, Schubert E, Persson P.B. Linear and non-linear properties of heart rate in postnatal maturation. Cardiovasc Res (1996) 31:447–454.
[Abstract/Free Full Text] - Mrowka R, Schlüter B, Gerhard D, Patzak A. Heart rate control in infants at high risk for sudden infant death syndrome (SIDS). Comp Cardiol (1996) 23:173–176.
- Mrowka R, Schlüter B, Patzak A, Persson P.B. Symbolic dynamics approach in the analysis of heart rate in premature babies at high risk for sudden infant death syndrome (SIDS). Comp Cardiol (1997) 24:37–40.
- Narayanan K, Govindan R.B, Gopinathan M.S. Unstable periodic orbits in human cardiac rhythms. Phys Rev E (1998) 57:4594–4603.[CrossRef]
- Pagani M, Lombardi F, Guzzetti S, et al. Power spectral analysis of heart rate and arterial blood pressure variabilities as a marker of sympatho-vagal interaction in man and conscious dog. Circ Res (1986) 59:178–193.
[Abstract/Free Full Text] - Palazzolo J.A, Estefanous F.G, Murray P.A. Entropy measures of heart rate variation in conscious dogs. Am J Physiol (1998) 274:H1099–H1105.[Web of Science][Medline]
- Parthimos D, Edwards D.H, Griffith T.M. Comparison of chaotic and sinusoidal vasomotion in the regulation of microvascular flow. Cardiovasc Res (1996) 31:388–399.
[Abstract/Free Full Text] - Patzak A. Development of linear and non-linear properties of heart rate control during quiet and active sleep in healthy infants. J Electrocardiol (1995) 28:356–357.[CrossRef][Web of Science][Medline]
- Patzak A, Schlüter B, Orlow W, et al. Linear and non-linear properties of heart rate control in infants at risk. Am J Physiol (1997) 273:R540–R547.[Web of Science][Medline]
- Persson P.B, Wagner C.D. General principles of chaotic dynamics. Cardiovasc Res (1996) 31:332–341.
[Free Full Text] - Persson P.B. Spectrum analysis of cardiovascular time series. Am J Physiol (1997) 273:R1201–R1210.[Web of Science][Medline]
- Pincus S.M. Approximate entropy as a measure of system complexity. Proc Natl Acad Sci USA (1991) 88:2297–2301.
[Abstract/Free Full Text] - Pincus S.M, Goldberger A.L. Physiological time-series analysis: what does regularity quantify? Am J Physiol (1994) 266:H1643–H1656.[Web of Science][Medline]
- Pool R. Is it healthy to be chaotic? Science (1989) 243:604–607.
[Free Full Text] - Poon C.-S, Merrill C.K. Decrease of cardiac chaos in congestive heart failure. Nature (1997) 389:492–495.[CrossRef][Medline]
- Rosenstein M.T, Collins J.J, De Luca C.J. A practical method for calculating largest Lyapunov exponents from small data sets. Physica D (1993) 65:117–134.[CrossRef]
- Schmidt J.A, Intaglietta M, Borgström P. Periodic hemodynamics in skeletal muscle during local arterial pressure reduction. J Appl Physiol (1992) 73:1077–1083.
[Abstract/Free Full Text] - Schmidt JA, Borgström P, Bruttig SP, Fronek A, Intaglietta M. Vasomotion as a flow-dependent phenomenon. In: Allegra C, Intaglietta M, Messmer K, editors. Vasomotion and flow motion. Prog Appl Microcirc. vol. 20. Basel: Karger, 1993. p. 34–51.
- Shinbrot T, Grebogi C, Ott E, Yorke J.A. Using small perturbations to control chaos. Nature (1993) 363:411–417.[CrossRef][Web of Science]
- Skinner J.E, Carpeggiani C, Landisman C.E, Fulton K.W. Correlation dimension of heartbeat intervals is reduced in conscious pigs by myocardial ischemia. Circ Res (1991) 68:966–976.
[Abstract/Free Full Text] - Skinner J.E, Pratt C.M, Vybiral T. A reduction in the correlation dimension of heartbeat intervals precedes imminent ventricular fibrillation in human subjects. Am Heart J (1993) 125:731–743.[CrossRef][Web of Science][Medline]
- Slaaf D.W, Oude Vrielink H.H.E, Tangelder G.-J, Reneman R.S. Effective diameter as a determinant of local vascular resistance in presence of vasomotion. Am J Physiol (1988) 255:H1240–H1243.[Web of Science][Medline]
- Spotlight Issue on Chaos in the Cardiovascular System. Cardiovasc Res 1996; 31.
- Stergiopulos N, Porret C.-A, de Brouwer S, Meister J.-J. Arterial vasomotion: effect of flow and evidence of non-linear dynamics. Am J Physiol (1998) 274:H1858–H1864.[Web of Science][Medline]
- Strogatz SH. Non-linear dynamics and chaos. With applications to physics, biology, chemistry, and engineering. Reading MA: Addison-Wesley Publishing Company, 1994.
- Sugihara G, May R.M. Non-linear forecasting as a way of distinguishing chaos from measurement error in time series. Nature (1990) 344:734–741.[CrossRef][Medline]
- Sugihara G, Allan W, Sobel D, Allan K.D. Non-linear control of heart rate variability in human infants. Proc Natl Acad Sci USA (1996) 93:2608–2613.
[Abstract/Free Full Text] - Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. Heart rate variability: standards of measurement, physiological interpretation, and clinical use. Circulation 1996;93:1043–1065.
- Tsuda I, Tahara T, Iwanaga H. Chaotic pulsation in human capillary vessels and its dependence on mental and physical conditions. Int J Bifurc Chaos (1992) 2:313–324.[CrossRef]
- Turcott R.G, Teich M.C. Fractal character of the electrocardiogram: distinguishing heart-failure and normal patients. Ann Biomed Eng (1996) 24:269–293.[Web of Science][Medline]
- Voss A, Kurths J, Kleiner H.J, Witt A, Wessel N. Improved analysis of heart rate variability by methods of non-linear dynamics. J Electrocardiol (1995) 28(Suppl):81–88.[Web of Science][Medline]
- Voss A, Kurths J, Kleiner H.J, et al. The application of methods of non-linear dynamics for the improved and predictive recognition of patients threatened by sudden cardiac death. Cardiovasc Res (1996) 31:419–433.
[Abstract/Free Full Text] - Voss A, Hnatkova K, Wessel N, et al. Multiparametric analysis of heart rate variability used for risk stratification among survivors of acute myocardial infarction. PACE (1998) 21:186–192.[Medline]
- Yeragani V.K, Srinivasan K, Vempati S, Pohl R, Balon R. Fractal dimension of heart rate time series: an effective measure of autonomic function. J Appl Physiol (1993) 75:2429–2438.
[Abstract/Free Full Text] - Wagner C.D, Persson P.B. Non-linear chaotic dynamics of arterial blood pressure and renal blood flow. Am J Physiol (1995) 268:H621–H627.[Web of Science][Medline]
- Wagner C.D, Mrowka R, Nafz B, Persson P.B. Complexity and chaos in blood pressure after baroreceptor denervation of conscious dogs. Am J Physiol (1995) 269:H1760–H1766.[Web of Science][Medline]
- Wagner C.D, Nafz B, Persson P.B. Chaos in blood pressure control. Cardiovasc Res (1996) 31:380–387.
[Abstract/Free Full Text] - Wales D.J. Calculating the rate of loss of information from chaotic time series by forecasting. Nature (1991) 350:485–488.[CrossRef]
- Witkowski F.X, Kavanagh K.M, Penkoske P.A, et al. Evidence for determinism in ventricular fibrillation. Phys Rev Let (1995) 75:1230–1233.[CrossRef][Web of Science][Medline]
- Witkowski F.X, Leon L.J, Penkoske P.A. Spatiotemporal evolution of ventricular fibrillation. Nature (1998) 392:78–82.[CrossRef][Medline]
- Wolf A, Swift J.B, Swinney H.L, Vastano J.A. Determining Lyapunov exponents from a time series. Physica D (1985) 16:285–317.[CrossRef]
- Yamashiro S.M, Slaaf D.W, Reneman R.S, Tangelder G.J, Bassingthwaighte J.B. Fractal analysis of vasomotion. Ann NY Acad Sci (1990) 591:410–416.[Web of Science][Medline]
- Yip K.P, Holstein-Rathlou N.H. Chaos and non-linear phenomena in renal vascular control. Cardiovasc Res (1996) 31:359–370.
[Abstract/Free Full Text] - Zwiener U, Hoyer D, Lüthke B, Schmidt K, Bauer R. Relations between parameters of spectral power densities and deterministic chaos of heartrate variability. J Auton Nerv Syst (1996) 57:132–135.[CrossRef][Web of Science][Medline]
- Zwiener U, Hoyer D, Bauer R, et al. Cardiovasc Res. (1996) 31:455–465.
This article has been cited by other articles:
![]() |
X. Bai, J. Li, L. Zhou, and X. Li Influence of the menstrual cycle on nonlinear properties of heart rate variability in young women Am J Physiol Heart Circ Physiol, August 1, 2009; 297(2): H765 - H774. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. E Aubert, S. Vandeput, F. Beckers, J. Liu, B. Verheyden, and S. Van Huffel Complexity of cardiovascular regulation in small animals Phil Trans R Soc A, April 13, 2009; 367(1892): 1239 - 1250. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. R. Brown, L. A. Cassis, D. L. Silcox, L. V. Brown, and D. C. Randall Empirical and theoretical analysis of the extremely low frequency arterial blood pressure power spectrum in unanesthetized rat Am J Physiol Heart Circ Physiol, December 1, 2006; 291(6): H2816 - H2824. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. J. Kocica, A. F. Corno, F. Carreras-Costa, M. Ballester-Rodes, M. C. Moghbel, C. N.C. Cueva, V. Lackovic, V. I. Kanjuh, and F. Torrent-Guasp The helical ventricular myocardial band: global, three-dimensional, functional architecture of the ventricular myocardium Eur. J. Cardiothorac. Surg., April 1, 2006; 29(Suppl_1): S21 - S40. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Torrent-Guasp, M. J. Kocica, A. Corno, M. Komeda, J. Cox, A. Flotats, M. Ballester-Rodes, and F. Carreras-Costa Systolic ventricular filling Eur. J. Cardiothorac. Surg., March 1, 2004; 25(3): 376 - 386. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

is plotted as a function of a. The Lyapunov exponent is a marker for a system's sensitivity on initial conditions and is negative for periodic evolution (for a<3.57) and positive for chaotic motion. Note the periodic windows in the chaotic range (3.57, 4) where 


