The training strategy involves sequentially addressing the core regulatory domains of 1) cerebral stability, 2) arousal regulation, 3) autonomic regulation, 4) affect regulation, and 5) executive function. A complementary target is the left heteromodal cortex. We review each in turn.
2.1. Cerebral Stability
The fundamental burden of any self-organizing control system is the maintenance of unconditional stability, sufficient to sustain basic functionality. Brain instabilities violating this criterion include seizures, migraines, panic attacks, asthmatic episodes, syncope, vertigo, sleep apnea, narcoleptic events, and suicidal episodes. The ability to address this disparate range of instabilities with ILF neurofeedback supports the view that it constitutes a comprehensive approach to the improvement of cerebral self-regulatory competence. Since the presence of such instabilities interferes with the subtlety that characterizes good regulation, addressing them must constitute the primary objective of a training strategy toward functional normalization.
The management of seizure susceptibility with neurofeedback has more solid literature support than prevails for any other condition [
61,
62]. Excellent results have also been reported for the remediation of migraine risk in a clinical cohort that had been treated to medical standards [
63]. With a quantitative EEG (qEEG )-based protocol, more than 50% of trainees ended up migraine-free for the one-year follow-up, versus none among those who relied on medication alone. Another 40% improved by more than 50% in incidence, versus 9% for the medication-only group. Only 2% failed to improve incidence significantly.
These excellent results notwithstanding, a significant step forward was taken with the realization that the classic brain instabilities tend to respond best to inter-hemispheric placements at homotopic sites, as shown in
Figure 3a. This discovery was made initially in working with migraines. With lateralized placement, a lateralized migraine may well just escape to the other hemisphere; targeting it there just brings it back. With inter-hemispheric placement, no such escape is available, so the migraine vanishes. Typically, this occurs within session, and often within just a few minutes. In yet other cases, the migraine is set on a path to resolution that takes a little longer. With extended training with the same protocol, migraine incidence typically subsides substantially, paralleling Walker’s findings. An outcome study on migraine headache using this method has just been published [
64].
Reliance on inter-hemispheric placement goes back to the early seventies, when Douglas Quirk used SMR reinforcement in bipolar montage at {C3 – C4} with violent offenders at the Ontario Correctional Institute in Canada. In combination with galvanic skin response biofeedback, he was able to reduce recidivism in selected violent offenders by a factor of more than three (from 65% to 20%, as documented in three-year follow-up) with 30 sessions of training. His work was the first application of inter-hemispheric training to a cerebral instability, as well as the first large-scale clinical utilization of NF [
65]. Quirk claimed to have worked successfully with more than 2700 inmates.
Once the impact of inter-hemispheric placements was recognized for migraines, it did not take long to establish that the same held true for brain instabilities in considerable generality. A key cerebral vulnerability appears to lie in the coordination between the two hemispheres. This coordination is managed at such a subtle level that a prescriptive training strategy is likely unavailing. Here the approach of endogenous neuromodulation is obligatory—both for its capacity for refined individualization and for its preferential access to the deep infra-low frequency domain. To re-establish control, the brain merely needs information about its own functioning; it is not in need of instruction.
Consider, in this regard, the attempt by Gruzelier et al to shift hemispheric negativity in schizophrenia that was referred to previously [
21]. If this procedure were adopted for bipolar disorder, the brain would be driven into mania in one case, and into depression in the other. The balance point would remain elusive to any prescriptive type of training. With endogenous neuromodulation, the balance point can be located by fine adjustment of the training frequency. Training too high takes the brain in a manic direction; training too low takes it toward the depressive state. This can take place in a matter of minutes in rapid responders. Training at the ORF serves to reinforce the network status that maintains stability against excursions into either mania or depression. In practice, the search of the optimal response frequency (ORF) might take a few sessions.
The process is illustrated in
Figure 4. Positive attributes can be combined to yield an Index of Functionality; adverse observables can be combined to yield an Index of Dysregulation. An optimum in both parameters is achievable at the optimal response frequency, and the therapeutic journey typically begins with the search for the ORF, particularly in highly symptomatic, and thus highly reactive/responsive clients.
When the brain is exposed to information relative to its state, just acting upon that information accelerates its journey through state space. The clinician adjusts the available parameters—placement and target frequency—to move the process toward its most propitious outcome. In a highly symptomatic individual, the journey to identify the ORF may well encounter adverse attractors in the attractor landscape [
4], which would call for rapid accommodation on the part of the therapist. There are compensations: the highly symptomatic individual is likely to be very sensitive to parameter change, thus facilitating the optimization procedure by way of prompt reporting. Conversely, when the trainee is less parametrically sensitive, then training at the ORF is also less critical.
In the search for a client’s optimal training parameters throughout the training process, we are compelled to operate clinically with what is termed “ipsative trend assessment,” i.e., the discernment of change induced via the training process in our clinical observables (66). In this process, the client’s report on subjective experience of change is heavily relied upon—if the client is capable of reporting, and particularly if the client is symptomatic. Symptom changes within session are the primary drivers of the optimization procedure. In the absence of symptoms that can be tracked within session, the process relies on state shifts with respect to issues on which the trainee can report: alertness, calmness, emotional ambience, etc. The objective is to reach the state at which the client feels maximally alert, calm, and euthymic. That specifies the conditions under which improved self-regulatory competence is achieved most efficiently in terms of both target frequency and placement. The conventional measures used in biofeedback can also be helpful here: galvanic skin response, finger temperature, heart rate and its variability, as well as scalp muscle tension. Over the course of sessions, additional information is recruited regarding the quality of sleep, and status with respect to drives such a hunger, pain, and cravings, depending on the client’s issues.
The above approach was adopted in the late nineties, when we were still working exclusively in the EEG spectral range. Here the ORF may have to be fine-tuned at the 1 – 3 % level. It may migrate modestly—and slowly—over the course of sessions, or when other protocols are introduced, but we regard it as a basically stable characteristic of a particular nervous system, at least over the time scales relevant in therapy.
Given the centrality of brain stability in neuroregulation, the hypothesis that all our clinical objectives might be met with inter-hemispheric training at various homotopic sites was evaluated in the early 2000s, concurrent with the migration down to the lowest EEG frequencies [
67]. Three placements turned out to be of greatest utility: {T3 - T4}, {P3 - P4,} and {Fp1 - Fp2}. Frontal placement {F3 - F4} played a role mainly in the anxiety-depression spectrum. Continuous performance test data, which we had been acquiring systematically on our clients since 1990, continued to yield satisfactory results [
68]. Nevertheless, over the course of two years the addition of lateralized placements was indicated, as shown in
Figure 3b. All our lateralized montages had either T3 or T4 in common. Very quickly, {T4-P4} and {T4-Fp2} came to dominate in the practice, indicating a right-side priority in the training—a complete reversal from the early days of the field, when the focus of cognitive neuroscience was on left-hemisphere function almost exclusively. At that time, right hemisphere function was still terra incognita.
What has survived the test of time in terms of principal training sites is shown in
Figure 3c. Observe that with the addition of lateralized protocols the inter-hemispheric training defaulted largely to {T3 - T4}. This schema turns out to be independent of whether one is training in the EEG band or at infra-low frequencies. The same networks are being targeted in essentially the same manner; we are merely appealing to different regions of the frequency hierarchy. However, the frequency range selected for training turns out to be determinative for outcomes.
2.8. The Frequency Domain
2.8.1. Extension of Training into the Infra-Low Frequency Domain
The evolution of frequency-based training through its critical transition into the infra-low frequency domain has had little representation in the literature to date. A case study on anxiety and fear in a cancer patient demonstrated that training could be done productively even in the delta and theta bands via endogenous neuromodulation [
72]. A study on Complex Regional Pain Syndrome captured the state of the art just prior to entry into the ILF regime [
73]. The larger scope of neurofeedback in application to chronic pain is covered in a book chapter [
74].
With availability of suitable software at our clinic in 2006, the training frequency could be extended to a bandpass of 0.1 Hz, down from its prior limit of 1.5 Hz (0-3Hz bandpass) of the NeuroCybernetics system. Within the first six months, some two-thirds of our clients ended up training at that lowest frequency, as illustrated in
Figure 5. This finding was a compelling argument to extend the frequency range further. That, however, required new instrumentation design (Cygnet, Bee Group, Switzerland).
10 mHz in target frequency was reached in 2008 with the new instrument. However, the pattern of cases piling up at the lowest frequency just repeated. Within just six months, the training band was extended down to 1mHz. Every step down in the available frequency range presented a somewhat greater challenge to signal processing, and a certain novelty to the clinician. Within a few months, however, the pattern repeated, with some two-thirds training at the lowest frequency, as shown in
Figure 6. In consequence, the range was extended to 0.1mHz in 2010. This once again pushed the limits of the available architecture, which called for yet another upgrade of the signal acquisition system, both hardware and software. This became available in 2013. By 2015, a repeat of the above pattern forced another extension of the range to 10 µHz (microHz). By 2017 we were working down to 1 µHz, and in 2019 the range was extended to 0.1 µHz. At the beginning of 2022, training in the 0.01 µHz range became available—and almost immediately became indispensable.
With every downward step, our clinical reach extended to some people we had not been able to help adequately before. Despite the increasing technical challenges as we went lower, the essential experience of the training—in the perspective of the trainee—was very similar across the entire range. In particular, the response to alteration of the training frequency was just as prompt at low frequencies as at higher ones. The brain is of course keying on the dynamics prevailing at the target frequency; it is not aware of the target frequency itself, which disappears entirely from view. Even signal extraction deep in the ILF spectrum must reflect real-time dynamics or it would be irrelevant to our project. The observable dynamics at issue lie in the spectral regime covered in fMRI measurements, 0.1-200 mHz, and that holds for all target frequencies that lie below 0.1mHz.
2.8.2. Optimal Response Frequency: A Resonance Phenomenon
With each highly responsive client, an optimal response frequency can be found at which the training is clearly ‘better’ than at neighboring frequencies, both above and below. We are still far from an understanding the neurophysiological basis of what we are dealing with here, but the clinical realities are difficult to dispute. The narrowness of the frequency range for best behavior argues strongly in favor of the resonance model for the ORF. The observed behavior also exhibits the characteristic of a resonance phenomenon more specifically, as illustrated in
Figure 7.
To explain our observations, resonance is best viewed in the complex plane, as elegantly described by Feynman [
75]. The real axis yields the magnitude of the response, which peaks sharply at the resonance frequency. It is described by the parameter Q, the ratio of center frequency to the width of the curve. Q falls in the range of 10-20, which serves to confirm the resonance model. The imaginary axis yields the phase response. The resulting system response differs above and below the ORF in a phase-sensitive manner. While the phase response is experientially complex and may be difficult for the client to describe, it is generally perceived as adverse in both cases, above and below the ORF, but differing qualitatively [
76]. That further supports the resonance model.
The strong frequency dependence becomes observable by virtue of the dysregulation status of the trainee. The variables at issue here are multiple, consisting of all that are being tracked through the training as indices to the physiological state of each individual. At the resonance frequency, the functional status is optimized, the dysregulation status is minimized (given the prevailing constraints), and the responsivity, the apparent rate of progress in training, appears to be optimized. Clients have sometimes referred to this as the “sweet spot.” The more severe the dysregulation status, the greater is the necessity to train at the ORF. Fortunately, this covaries with the ability to discriminate it, as indicated in
Figure 7c.
For example, it was observed early on that a migraine aura precipitated in a particular client with training at 3.9 Hz could be disrupted by moving to 4.0 Hz. With return to 3.9 Hz the aura would be rekindled. Movement back to 4.0 would cause the aura to subside once again. This phenomenon was reproducible, all in the course of a single session. Such a sharp frequency dependence is the hallmark of the resonance response, in considerable contrast to the gentler dependence on frequency observed elsewhere across the spectrum. A Q of ~20 is indicated in this case. This kind of evidence established the case for the resonance model for the ORF originally in the EEG spectral range some 25 years ago. The stability of the ORF over the course of sessions—it tends to migrate only incrementally and slowly with a given protocol—argues for the hypothesis that the brain organizes the frequency domain around these key frequencies, which are dispersed throughout the ILF and EEG ranges.
Such a frequency hierarchy has already been postulated by Fujimoto and Kaneko in a theoretical treatment that modeled the coupling between high and low frequency dynamics. Under a set of limiting assumptions, the behavior of chaotic elements at high frequency propagates to low frequency in the presence of nonlinear oscillators broadly and uniformly distributed in time scale [
77]. Specifically, the authors assumed a chain of nonlinear oscillators that were distributed as a power series with a coefficient of just less than two. Although the analysis was performed under very specific conditions, the authors nevertheless conjecture that “…the propagation from faster to slower elements as described in this paper is a universal property of systems of coupled chaotic dynamical elements with distributed time scales…”. The present work furnishes empirical support for the existence of such oscillators extending deep into the infra-low frequency sphere. In the clinical domain, we also know that bipolar and depressive cycles can extend over many months in a quasi-periodic manner, and that could be reflective of one of these deeper rhythms.
2.8.3. The Frequency Rules
The frequency rules that apply to ORFs in lateralized placements are shown in
Figure 8 [
78,
79,
80]. They divide into a high-frequency and a low-frequency domain and are observed to be consistent within their respective regions. Both a linear and a logarithmic plot illustrate the frequency relationships, which are arithmetic in the EEG range and geometric in the ILF range. The left hemisphere optimizes at 2 Hz higher than the right and optimizes at a factor of two higher in the ILF range. The Delta band is the transition region in the two domains, as the two criteria coincide at a right-hemisphere frequency of 2 Hz and a left-hemisphere frequency of 4 Hz (i.e., 2 x 2 = 2 + 2).
Different rules apply to inter-hemispheric placements. This is easily said, but it was a long process to tease them out. That was not the case for the lateralized placements, where the legacy protocols of “C3-beta and C4-SMR” differed by 3 Hz, and it was merely a matter of systematic trial to establish that the optimal difference was really 2 Hz rather than 3. For the inter-hemispheric placements, it was noted that the cortical resting rhythms differed by nominally 4 Hz between the posterior sensory regime (i.e., 10 Hz) and the central somatomotor regime (~14 Hz). Thus, it was found that posterior placements optimize at 4 Hz lower than central (and a factor of four lower in the ILF), except for posterior temporal sites (T5/T6), which train identically to T3/T4. Frontal placements, in turn, optimize at 2 Hz lower than central placements (a factor of two lower in the ILF). In actual practice, with most clients optimizing at the lowest available frequency in the ILF so much of the time, and with T3-T4 having such a broad clinical footprint, the opportunity to probe other inter-hemispheric placements at yet lower frequencies has been limited until recently.
It is appropriate at this point for the authors to acknowledge that some of what has been presented no doubt strains credulity on first exposure. We have had the benefit of years to do reality-testing, progressing incrementally by way of Bayesian inference, cumulatively confirming/disconfirming provisional hypotheses. The entire construct rests substantially on ‘soft data,’ namely the subjective reporting of highly symptomatic clients with respect to categories that are not readily quantifiable—and inherently variable. Particular observations are often non-repeatable, as the state vector is itinerant. And yet what has emerged is the hardest of testable hypotheses: the frequency rules. They are constantly under test within the global practitioner network. The existence of the frequency rules validates the ORF principle, and the existence of the ORF in turn validates the approach for elucidating it. There is no alternative to the within-subject design for proving out these hypotheses. Every clinical case is its own A/B design, with an ongoing protocol contingency at every juncture.
2.9. Mechanisms Implications
The most striking feature of this entire development is the long-term trend toward ever lower target frequencies. This pattern is driven by the issue of arousal regulation. Those whose training begins with the {T4 - P4} placement exclusively tend to train at the lowest frequencies. By contrast, those whose training begins with the inter-hemispheric protocol exclusively have their target frequencies much more broadly distributed. Over the course of training, most clients (>90%) are likely to experience both starting protocols. And it is found that in such cases the inter-hemispheric placement is more tightly constrained in terms of ORF, as is to be expected for cerebral instabilities. In generality, all right-lateralized placements train at the same ORF as inter-hemispheric placement at {T3 - T4}. This fact alone serves to establish the primacy of {T3 - T4} among alternative inter-hemispheric placements.
The implication of the above divergence in spectral distributions is that our two principal failure modes arise from different sources. Arousal dysregulation tends to be the result of environmental insult or persistent duress, whereas brain instabilities tend to have a genetic foundation, one that promotes the vulnerability to hyper-excitability. The training of both failure modes is favored at the ORF, as the frequency domain organization imposes its own constraints.
Another striking consistency is the right-side dominance that has emerged with training in the ILF regime, with the combination of T4-P4 and T4-Fp2. This has all been empirically driven, with the brain effectively in a controlling role with respect to protocol priorities. However, it is consistent with the focus on early developmental priorities. Right-side function is the first to mature in infancy and early childhood [
81]. More broadly, we observe that the two principal failure modes that have been identified are not just relevant to those afflicted with early trauma, or are otherwise highly deficited, but to our clinical population in general. Our two starting protocols are essentially universally applicable, which implies that we are all subject to the identified key vulnerabilities. Right-hemisphere function is therefore our primary clinical target across almost the entire range of mental-health related dysfunctions that we encounter, and across our entire clinical population, with only rare exceptions, if any.
2.9. The Neurofeedback Therapist’s ‘Systems Perspective’ and its Foundations
The ‘systems perspective’, from the standpoint of a neurofeedback therapist, is shown in
Figure 9. We assert the primacy of the regulatory arc that begins with interoception, informs autonomic regulation and affect regulation, and ultimately governs the tonic ambient of central arousal. All these core functions reside in a state of intimate co-regulation. They are under the primary management of the right hemisphere, and this coordination is organized in the ILF regime. Personal history stamps the proceedings via hippocampal, afferent vagal and cranial nerve pathways. The history of traditional biofeedback has long testified to the intimate association of affect regulation with the autonomic nervous system. The latter is trained in order to tame the former, as for example in the management of anxiety.
The foundation for the above conception was laid years ago in the animal work of Nina Aladjalova, whose book became available in the English language in 1964 [
82]. Aladjalova studied the infra-slow rhythmic potential oscillations (ISPOs) at great length. “A single stimulation of the reticular formation immediately elicits an arousal reaction in the EEG of the cortex, but has no effect on infraslow activity,” she writes. “This reaction is apparently regulated by the rapid regulatory system. Stimulation of the ventromedial part of the hypothalamus…intensifies infraslow cortical activity within 30-40 minutes. This reaction is presumably regulated by the slow regulatory system.”
“…infraslow activity is intensified by certain actions after a long latency period, 30-100 and 120-200 minutes later. We conjectured that this phenomenon reflects the activity of the slow control system of the brain…not only to automatically adjust the system to keeping internal environment constant but actively to establish a new level of activity” (81).
Even in an optimal functioning context, this remains the priority in training. We have been engaged with the slow control system that Aladjalova first identified, and it has been more a matter of homeodynamics than of homeostasis. This system is centrally regulated by the hypothalamus, which governs our internal milieu—autonomic function, sleep-wake cycle, ultradian rhythms, etc. The brainstem is of course central to arousal regulation via the ascending reticular activating system, and as such constitutes the top of the cerebral regulatory hierarchy. Phasic arousal is subject to input from the thalamocortical networks. The tonic level of arousal is in turn informed by a second regulatory hierarchy, under hypothalamic control.
It is appropriate at this point to make the connection, as best we can, between the principal protocols and the targeted core state variables. These protocols were already well-established before our acquaintance with intrinsic connectivity networks (ICNs). Our working model in the late nineties was that we were engaging the multi-modal association areas, the highest level of integration of our sensorium. These are the last areas to mature in neural development, being associated with the highest level of plasticity. According to a more recent study, these regions also exhibit the highest connectivity gradient [
83]. That is the equivalent of saying that these sites are the most highly integrative in character.
In the perspective of the ICN model, we are training the sites where the Default Mode Network is accessible to us at the cortical surface at lateralized sites. The two criteria are convergent. Further, the multi-modal association areas serve as input to the salience network, so we are training the nexus of the Sensorium, the Default Mode and the Salience network. Our early work in the ILF regime was influenced by Buckner et al. [
84], which examined the connectivity relationships among the hubs of the DMN. We’ve created our own graphical representation of the data presented there, and this is shown in
Figure 10. In addition to the general argument that it is most efficient to train the relationships between the hubs, there is the subsidiary argument that one would like to train those linkages that exhibit the highest connectivity. Those network linkages that the brain keeps under the tightest control make for the most discriminating sources of information back to the brain. These considerations further underpin the primacy of T3, T4, P3, and P4 in our protocols. Buckner’s data also supports the case for right-side priority in the training. Observe that the connectivities are generally larger on the right side, and that T4 is more intimately inter-connected than T3. This biases us toward right hemisphere training.
In order to illuminate the relationship between the two hemispheres, we draw on a seminal paper that yields information flow among the principal hubs of the Default Mode Network. These hubs were originally identified in the study of microstates by the Lehman group in Switzerland [
85]. Each of the microstates is identified with one of four hubs of the DMN, three in the posterior region and one with an anterior locus. The two principal hubs, anterior and posterior, lie along the midline, and two posterior hubs are lateralized.
Information flow among these hubs has been determined by means of a measure of directed coherence in the alpha and low beta bands, leading to the finding that information flow was dominant from the left hemisphere to the right, as well as from the left to the midline hub, relative to the flows the other way [
86]. The imbalance can be substantial. The clear implication is that with respect to the regulatory role of the lower EEG bands, the left hemisphere is in a commanding position with respect to the right hemisphere.
A division of responsibilities is indicated. We may infer from the above that a reciprocal relationship exists in which the right hemisphere bears the primary burden of organizing our resting states, and in that capacity also governs the left hemisphere. The left hemisphere, in turn, supervises our engagement with the outside world, and in that role also governs right hemisphere function. The ILF regime plays a primary role in organizing the resting state configuration, whereas the EEG regime handles the complexity, the coordination, the immediacy, and the temporal precision required for our interface with the outside world. The upper delta band falls in the middle ground, with right hemisphere primacy extending up to 2 Hz, and left hemisphere primacy extending down to 4 Hz, according to the frequency rules.
Vinod Menon proposed the triple network model of psychopathology in 2010 (preprint), postulating that “aberrant organization and functioning of the Central Executive Network, the Salience Network, and the Default Mode Network are prominent features of several major psychiatric and neurological conditions [
87]. With our limited focus on remediation—as opposed to phenomenology—it now appears that psychopathology is much more explicitly rooted in the Default Mode Network and the Salience Network, and particularly in their tonic regulation within the ILF regime. The Central Executive is essentially missing from the conversation. A right-hemisphere bias in early development has also been identified by Schore [
88]. This, then, defines our agenda with respect to the primary regulatory arc in
Figure 9.