State-space intermittent feedback stabilization of a dual balancing task

Balancing the body in upright standing and balancing a stick on the fingertip are two examples of unstable tasks that, in spite of strong motor and sensory differences, appear to share a similar motor control paradigm, namely a state-space intermittent feedback stabilization mechanism. In this study subjects were required to perform the two tasks simultaneously, with the purpose of highlighting both the coordination between the two skills and the underlying interaction between the corresponding controllers. The experimental results reveal, in particular, that upright standing (the less critical task) is modified in an adaptive way, in order to facilitate the more critical task (stick balancing), but keeping the overall spatio-temporal signature well known in regular upright standing. We were then faced with the following question: to which extent the physical/biomechanical interaction between the two independent intermittent controllers is capable to explain the dual task coordination patterns, without the need to introduce an additional, supervisory layer/module? By comparing the experimental data with the output of a simulation study we support the former hypothesis, suggesting that it is made possible by the intrinsic robustness of both state-space intermittent feedback stabilization mechanisms.

as the state vector enters one of the unsafe regions (on-phase). What is important is that, during the on-phase, the purpose of the control action is not to attract the state vector towards the nominal equilibrium state but to allow the state vector to approach or cross the stable manifold, thus turning off the control when this event is detected. The bounded stability that can be achieved with this intermittent feedback approach is quite robust because it can work also with delayed information of the state vector, producing a limit cycle as an alternation of segments of hyperbolic orbits (off-phases) and spiral orbits (on-phases).
In previous studies, it was demonstrated that the state-space intermittent feedback stabilization paradigm can explain in a detailed manner the oscillatory patterns of two very different balancing tasks: quiet upright standing [20][21][22] and stabilization of a CIP 23,24 . In both cases, the intermittent control paradigm matches the rationale of the minimum intervention principle 25 and in any case minimization of energetic costs is a relevant feature at stability's edge 26 . The two tasks differ in a number of ways: the former one has been perfected phylogenetically and ontogenetically during neurodevelopment, whereas the latter is an example of an unstable task that requires learning and adaptation to the degree of difficulty, related to the length of the pendulum and thus to the corresponding falling time constant. These tasks are also characterized by a different number of degrees of freedom, different muscle groups, a different role of muscle stiffness, and a different involvement of sensory modalities. Last but not least, this control paradigm can also explain the multi-joint coordination that underlies the apparent over-simplification of the single inverted pendulum model of quiet standing 27 , focusing on the oscillation of the virtual inverted pendulum that links the ankle to the changing CoM. Moreover, neurodevelopment, for the stabilization of upright standing, or learning, for the stabilization of the CIP, both require a careful tuning of the sensory-motor parameters of the internal model that continuously monitors the evolution of the state, by optimally fusing the appropriate sensory channels. Summing up, although the intermittent control paradigm is the same, the implementation details are quite different and thus it seems unlikely these tasks are served by the same internal control model, suggesting instead that two instantiations of the same state-space intermittent control policy are implemented by the central nervous system as a pair of abstract internal control modules. In a sense, this may be considered an extension of the principle of motor equivalence 28,29 and opens the question about the mechanisms used by the brain for dealing at the same time with a pair of unstable tasks, characterized by similar dynamics but very different sensory motor machinery. In particular, it is natural to formulate the following question: in a dual stabilization task, that requires the coordination of the two skills, the corresponding state-space intermittent feedback controllers require an additional control level or can automatically compensate the interaction determined by the biomechanics of the body and the physics of the CIP?
Thus the purpose of this study is to attempt to answer this question by extending the state-space intermittent feedback framework to the investigation of a dual balancing task. Although it is clear that dual vs. single task effects have been extensively investigated, most studies involved motor-cognitive tasks with a single balancing component, in a variety of situations, e.g. in young adults 30 , in elderly 31 , or in patients 32 . Moreover, the attention was focused more on the quantification of the interaction effects rather than the underlying motor coordination and control problem. In this study we investigated the dual task of maintaining upright balance while stabilizing with the hands an inverted pendulum device. Figure 2 shows the experimental paradigm: a young, healthy subject stands on a force platform and is required to keep her feet fixed during the stabilization of a CIP-like device that consists of a wooden bar (0.4 m long, 0.02 m diameter) kept firmly with the two hands. In the middle of the bar there is a ball-bearing connected to a wooden inverted pendulum (1 m long). The motion of the subject and of the CIP-like device was measured by means of a motion capture system and the ground reaction force by a force platform, thus acquiring the time-course of the CoP (Center of Pressure). The combined oscillations of the body and CIP were also compared with the simulation results of a pair of state-space intermittent controllers applied to the biomechanical model of the body-CIP system sketched in Fig. 2.
A similar experimental paradigm has been investigated recently 33,34 , although with different goals, mainly concerned with the influence of the learning process of stick-balancing on the coupling with postural control. In contrast, this study was focused on relevant features of the dual-task balance control after learning. However, the two studies agree on a similar conclusion about the coupling between the two controllers after training, namely the independent activation of the two controllers.

Experiments.
The experiments involved fourteen healthy young subjects after a short training of stick balance. Consider that there is dynamic interaction between the two tasks. The purpose of the arm movements is to apply a force f t ( ) s via the hand to the CIP-like device, capable to induce an acceleration profile of the stick x t ( ) sä ppropriate for balancing the pendulum. This force, with a minus sign, is reflected at the shoulder together with the force f t ( ) a , necessary for accelerating the mass of the arm while it is carrying the hand, back and forth, with the mentioned acceleration profile. Thus, f t f t ( ( ) ( )) s a − + is a self-generated disturbance to the upright stance stabilization induced by the CIP stabilization process. The preliminary question we needed to answer was then the following one: to which extent such dynamic interaction modifies the posturographic spatio-temporal features that characterize quiet standing? Figure 3 shows typical patterns related to one subject indicating that the anticipated dynamic interactions modify only in a rather minor way the posturographic descriptors (see Table 1), while keeping the same structure: 1) the angular sway size (expressed by the standard deviation of θ com ) is markedly increased, 2) the phase portrait is slightly shifted in the direction of the CIP, as indicated by the small change of the mean value of θ com from the single to the dual balancing task, but keeps its structure, 3) the power spectral density graph is basically unchanged: panel E of Fig. 3 shows that the bi-logarithmic plot is simply shifted upward, because the sway amplitude is increased, but is roughly approximated by a straight line with the same negative slope in the low-frequency range, in agreement with the power law scaling regime 35,36 . These authors found that temporal patterns of postural  63 . Markers for motion captures are attached to the body and to the CIP-like device; the subject stands on a force platform; surface electrodes record the electrical activity of different muscles of the legs/trunk/arms, however the analysis of their activation patterns were not included in this study. Right panel: scheme of the dual balancing task. BIP: Body Inverted Pendulum; VIP: Virtual Inverted Pendulum. In the single balancing task there is no CIP-like device and the two arms are kept extended on the two sides of the body: in this case BIP and VIP coincide as well as the two angles θ b and θ com . sway, i.e. the time course of the CoP (Center of Pressure), for healthy adults exhibit power-law-like behavior in the low-frequency regime: increments of the CoP path behave as those of a random walk with negative correlation, corresponding to movements approaching the upright posture. Table 1 reports, for all the subjects, the amplitude of sway (expressed as standard deviation of θ com ) in the single and dual task, respectively, as well as the amount of the shift of the median value of the sway angle: on average, the amplitude is increased from 0.23 deg in the single task to 1.18 deg in the dual task, with a forward shift of 0.9 deg. Figure 4 shows typical combined CIP/body balancing patterns, during the dual task, in the time and frequency domains and Table 2 reports characteristic parameters of the different subjects. As regards the CIP-like device, generally speaking, as one may expect, such patterns are similar to those recorded in studies where the CIP was stabilized by a sitting subject 23,24 . The range of motion of the stick, measured by the standard deviation of the θ s angle, is close to 2 deg; the spectrum of such oscillations has a well-marked peak, below 1 Hz (0.39 Hz on average); the range of hand motion, measured by the standard deviation of x s , is about 0.1 m. However, what is more interesting is the interaction/coordination between the two balancing motions. In particular, there is no correlation between the oscillation of the body θ com and the oscillation of the stick θ s , whereas there is a positive correlation between the hand forward/backward motion x s and the body angle θ com (0.87 on average) and a negative correlation between x s and θ s (−0.58 on average). While the latter anti-correlation may be attributed to a pure mechanical effect, expressing the bounded stability of the stick inverted pendulum, the former one suggests an indirect synergy between the two intermittent stabilization processes, aimed at the extension of the range of movement of the hand motion that is required for improving the chance of stabilizing the stick while keeping the feet fixed in the starting position. We label indirect the synergistic effect in the sense that it is not coded in a specific modification of the controller; in contrast, it expresses the natural consequence of the dynamic interaction between the two stabilization processes. The acceleration profile of the quick back and forth arm movements learned by the subjects for keeping the stick from falling is, at the same time, a disturbance for the standing body stabilization process and the purposive action generated by the stick stabilization process. As we will show in the following section, this hypothesis is supported by the simulation studies, further emphasizing the robustness as well as the flexibility of the state-space intermittent feedback control paradigm.
Simulations. The dynamical model of the CIP-like device has two degrees of freedom (θ x & s s ) but is an under-actuated system because the only control variable available (the force f t ( ) s applied by the hand to the CIP) is unable to regulate simultaneously the state of the two controlled variables. However, the goal of the control process is more limited: it consists of keeping one variable (θ s ) in a small neighborhood of equilibrium, while allowing the other variable (x s ) to oscillate semi-freely in a large range compatible with arm-length. The model, that can be derived by the Lagrange equations and has been already employed in a previous work 24 , is expressed by the following equation:

On phase
Activation condition:

Control Action:
Off phase Dis activation condition: Control Action x are the controller gain parameters: the function of the first two is to constrain the stick oscillations near a limit cycle in the phase plane around the vertical, thus avoiding the fall, while the function of the other two elements of the control action is simply to limit the CIP motion to a physiological range, compatible with the arm length. The two groups of gains must be tuned considering a trade-off between stick motion and hand motion: increasing the hand gains will reduce the range of hand motion but increase the range of stick motion, with a greater risk of fall; in contrast, reducing the hand gains will also reduce the stick risk of fall but may force the hand beyond arm reachable positions. In the simulation experiments the control force was affected by a white noise: its power could be increased in the simulations to the same value of the control force before driving the system to instability.
The dynamical model of the body inverted pendulum has one degree of freedom (θ b ) and is expressed by the following equation: • I b is the moment of inertia of the body around the ankle.
is the gravity destabilizing torque (m b is the body mass and L b is the distance of the body CoM from the ankle).
a b is the torque due to the viscous-elastic properties of the ankle muscles (K a is the stiffness and B a is the corresponding viscous coefficient: in accordance to experimental evaluations 7,8 < K K a g ); • T int is the state-space intermittent control law that supplements the insufficient stabilizing effect of T stiff for counteraction T grav ; is the disturbing torque for the body stabilization due to the control force transmitted to the CIP device (L s is the distance of the shoulder from the ankle); b is the associated disturbance torque due to the acceleration of the arm mass. www.nature.com/scientificreports www.nature.com/scientificreports/ In particular, the state-space intermittent control action, which is formally quite similar to Eq. 2, is described by the following equation:

On phase
Activation condition: Control Action: Off phase Dis activation condition: Control Action: are the controller gain parameters: their function is to constrain the oscillation of the body inverted pendulum near a limit cycle in the phase plane; α is the slope of the switching function in the phase plane.
The state-space intermittent control models described by Eqs. 1-4 were simulated for both the single and dual balancing tasks. In the former case the two characteristic angles of the body inverted pendulum (θ b and θ com ) coincide, whereas diverge in the latter. Figure 5 shows a typical simulation result of the dual balancing task. Tables 1 and 2 compare the simulation (last line of each table) with the experimental results (penultimate line of each table) for both types of tasks. In particular, Table 1 shows the increase of the amplitude of the body sway from the single to the dual balancing task, associated, in the latter case, with a systematic forward shift of the tilt www.nature.com/scientificreports www.nature.com/scientificreports/ angle; Table 2 shows the main parameters of the dual task (range of motion of the stick angle and the stick position, frequency peak of the stick spectrum, and correlation coefficients of θ .
x vs s s and θ .
x vs s c om ). The robustness of state-space intermittent control paradigm of the body is supported by the experimental results on the amplitude of sway which can accommodate the rather massive self-generated disturbance. The simulation experiments allowed us to evaluate that such disturbance, resulting from the acceleration of the arm and the CIP device, is close to 90% of the total torque acting on the ankle (with a standard deviation of about 3.7 Nm) during the dual balancing task. Of course such disturbance is mostly determined by the arm movements because the weight of the arm is much greater than the CIP device.
Overall, the working hypothesis of testing the capability of the two independent state-space intermittent controllers to successfully stabilize the dual balancing task without any specific modification of the control actions is confirmed. The spatio-temporal features of the experimental and simulation behaviors are quite similar. In particular, the positive correlation between the hand motion and the body sway as well as the anti-correlation between the hand motion and the stick oscillation are found in both cases. Moreover, the Power Spectral Density (PSD) of the body sway angle in both balancing tasks exhibits the same power law scaling regime in the low frequency range (up to 2-3 Hz) for all the subjects as well as the model simulations. This is consistent with the analysis of universal and individual characteristics of postural sway during quiet standing 37 that has shown that the power-law behavior at the low-frequency regime is a universal indicator of the control law whereas the high-frequency regime is sensitive to the anthropometric/biomechanical parameters.
We also analyzed the activation/inactivation sequences generated by the pair of state-space intermittent feedback controllers in the dual balancing task. Figure 5 (panel F) shows a typical pair of sequences of on-off activations of both intermittent controllers. It appears clearly that the average switching rate of the body intermittent controller is lower than that the CIP controller one: 0.68 Hz vs. 1.03 Hz. In both cases, the duration of the on-phases is somehow longer than the off-phases: with a mean duration of the on-phase of 0.880 s (on-phase) vs. 0.585 s (off-phase) for the body controller and 0.537 s (on-phase) vs. 0.436 s (off-phase) for the CIP controller. Moreover, the on/off switching times do not appear to be correlated. Thus, the behavioral coordination and coherence that characterizes the dual balancing task, observed in the experiments and reproduced in the simulations, seems to be consistent with two internal control models, with similar design but independent control actions, that interact through the biomechanics of the body and the mechanics of the external environment.

Discussion
It has been suggested 15 that intermittent control is the natural and rational solution devised by the human brain for answering the fundamental question of motor neuroscience: how the human body may succeed to perform a fast motor task in a real-time fashion with the slow sensory-motor system. This is somehow in contrast with the common wisdom that the brain stabilizes unstable body dynamics using impedance control, by regulating co-activation levels of antagonist muscles 38 . In this framework, the strategy of the CNS would be to learn the optimal impedance for compensating the destabilizing effect of gravity, by selecting the appropriate groups of antagonist muscles and optimally tuning co-activation levels in a preprogrammed manner 39 . Such feedforward control strategy 40,41 has the clear advantage of avoiding the risk of delay-induced instability but, on the other, it forces a trade-off between task-related errors and energetic costs. Moreover, accepting such high energetic costs seems to be incompatible with the minimum intervention principle, which has been suggested as a general, rational strategy for many biological systems 25 . As a matter of fact, the intermittent control strategy for defeating instability and shaping purposive actions is clearly in favor of the minimum intervention principle, by exploiting the natural dynamics arising from the interaction between the human body and the external environment. At the same time, this strategy agrees as well with the general principles of biological autonomy advocated by Francisco Varela 42 in the framework of enactivist theories.
The intermittent control strategy is a general approach that can be articulated in a variety of versions of intermittent controllers: one version assumes anticipatory ballistic bias control, mathematically modeled with a state predictor for compensating feedback delay [43][44][45] ; in a second version, off-loop phases are meant to represent a sensory dead-zone 46 ; a third version has been characterized as act-and-wait control 47 ; another version is the one adopted in this study 20,22,48,49 , namely the state-space intermittent feedback strategy. Such control paradigm exploits the stabilizing effect of one part of the saddle, letting the system evolve by alone when it slides on or near the stable manifold; moreover, the key feature of this control strategy is that, although the off-and on-systems are both unstable, the combination of the two according to the switching mechanism can achieve bounded stability, with a limit cycle oscillation. It is somehow surprising that this simple control paradigm could be applied with success to two widely different skills like quiet upright standing and CIP balancing. In the former case the control action is an ankle torque provided by the ankle muscles and in the latter case is the force applied by the hand to the cart; in both cases the dynamics is characterized by the interaction of the CoM with the CoP, the CoM being the controlled variable and CoP the control variable. What is remarkable is that the CoP is a "real" point in CIP balancing (the position of the lower end of the stick) whereas it is a "virtual" point in upright standing (the origin of the ground reaction force vector): in summary, the state-space intermittent feedback control paradigm can be seen as a CoP strategy, organized in a similar way in both cases. The biological plausibility of this control paradigm, in comparison with the conventional PID-like continuous control, is supported by a number of empirical observations. First of all, the intermittent model is much more robust than the standard model because the size of the region in the parameter space of the feedback control gains that produce stable behavior is significantly larger in the former case than in the latter one, both for upright standing 22 and stick balancing 23 . Furthermore, the intermittent controller can use feedback parameters that are much smaller than the standard model: this is also reflected in the fact that the standard model can only achieve, if successful, a stronger and more limiting form of stability, where the PSD of the sway movements is similar to an over-damped second order system without a www.nature.com/scientificreports www.nature.com/scientificreports/ resonance, a behavior that does not match at all the power law scaling regime, typical of physiological sway 35,36 ; in contrast, the intermittent control model can settle in a weaker stable regime (a limit cycle) by using feedback parameters that are much smaller and producing sway movements whose PSD matches the physiological ones. In the specific case of stick balancing a recent study that investigated the biological plausibility of the state-space intermittent feedback control model 23 found that the produced stick-sway movements are characterized by a non-Gaussian, truncated Lévy distribution 50,51 typical of the experimental data.
The dual balancing task introduces two additional issues to the discussion above on the rationale of intermittent control: the issue of dual task, in general, and the more specific issue of anticipatory postural adjustments (APA's). A dual task can be defined as the concurrent performance of two tasks that can be executed independently and have distinct and separate goals. It is well known that when attempting to perform two tasks at the same time, the tasks often interfere with each other. Typically such interference has been studied by using the refractory period paradigm 52 and has been modeled in different ways in terms of how to manage the increased attentional load, e.g. the central bottleneck model 53 or the central capacity sharing model 54 . The dual balancing task investigated in this work, which is just an example of equilibrium skills exhibited in the circus, in our opinion is a different type of challenge, primarily because for all balancing tasks the issue is not accuracy per se but success/failure of the action, namely avoiding the fall (of the body, the stick etc.). For the specifically investigated dual balancing task, the main result of the study is that the two independent state-space intermittent feedback controllers succeed to achieve and maintain dual bounded stability without any additional supervisory control mechanism, while matching the spatio-temporal features of both sway movements, experimentally recorded. This result characterizes subjects trained enough to exhibit a repeatable performance, namely a behavior at the end of learning. Considering that learning was not addressed in our work, the plausibility of an independent control organization of the dual balancing task is also consistent with a work specifically focused on the learning process of the task 33 , which showed that the coordination patterns between finger movements and CoP shifts appear to evolve, in the course of learning, from collective to independent control. At the same time, we should also consider that although the two intermittent controllers appear to be independent, they also have a strong dynamic interaction through biomechanics itself: thus, the fact that they do succeed to assure dual bounded stability emphasizes their robustness and ability to compensate at the same time the intrinsic and extrinsic destabilizing effects. Such functional synergy, without explicit supervisory control, emerges in spite of the strong differences between the two balancing paradigms: the state-spaces are different, the sensory feedback channels are different (primarily visual for stick balancing and proprioceptive in upright standing), the number of degrees of freedom is different, the recruited groups of muscles are different. Consider that both balance task controllers have an important central/ cognitive component because neither of them can be reduced to a reflex in any sense but requires multisensory integration and decision making; on the other hand, different from the typical attentional tasks addressed in the investigation of dual task interference, such cognitive component is strongly grounded in an embodied cognitive framework: more specifically, it is based on the exploitation of the dynamic affordance provided by the saddle-like instability of an inverted pendulum, whether a stick or the human body. However, additional research is suggested by recent work 55 for integrating a task specific cognitive layer capable to introduce a cognitive dependence on top of the control independence supported by this study.
Modeling the CIP-like device as a single inverted pendulum (SIP) is obviously accurate enough, whereas adopting the same paradigm for upright standing, with the implicit assumption that only ankle rotations are relevant for describing and explaining sway movements, is certainly a less accurate approximation. In particular it ignores hip motion and ankle-hip coordination that have been the focus of recent experimental and modeling studies [56][57][58] . The range of variation of the angular displacement, velocity, and acceleration of the hip is comparable to that of the ankle, thus suggesting that the SIP model should be substituted by a two-link or DIP (Double Inverted Pendulum) model, involving the coordinated control of ankle and hip joints. In particular, it has been found that the acceleration profiles of the two joints are strongly characterized by anti-phase correlation and it was proposed that DIP control could be modeled as an optimal bi-axial active controller with the goal of minimizing the acceleration of the global CoM 59 . However, a recent simulation study involving a DIP mechanical model stabilized by a state-space intermittent controller 27 demonstrated that there is no need to introduce a bi-axial optimization process because the state-space intermittent controller applied to the VIP (Virtual Inverted Pendulum that links the ankle to the variable CoM) can fully explain ankle-hip coordination, while hip motion is simply regulated via muscle stiffness.
As regards the issue of APAs, for understanding dual balancing tasks, we should take into account that different types of APAs, that represent generally feed-forward control processes, can be defined. The traditional group of APAs occurs when a standing person performs an action leading to a postural perturbation or expects an external postural perturbation, inducing changes in the activation levels of postural muscles that can be observed prior to the perturbation time 60,61 , typically about 0.1 s prior to the movement initiation or the perturbation time 62 . A second class of APAs occurs when a person prepares to make a whole-body action that may destabilize the standing posture, such as picking up a load from the floor and placing it on a table: also in this case the main action (reaching/picking/placing) is integrated in a feedforward manner with a rearrangement of body masses in order to avoid a destabilization of the body. At a first look either type of APAs might seem to be irrelevant for the paradigm investigated in this study because APAs are typically feed-forward processes whereas the dual balancing task implies dual feedback control. However, it is an intriguing issue, to be addressed in future studies, to consider hybrid experimental situations, approximating complex tasks or task sequences in real life, where the two control paradigms (feedforward anticipation and intermittent feedback stabilization) may need to be integrated and coordinated by the brain.
A preliminary version of this work, that was limited to three subjects and did not include the formulation and simulation of the dual intermittent control model, was presented at a recent conference 63

Methods
Subjects. Fourteen healthy young subjects with no known neurological impairments took part to the experiments ( Table 3). None of them had previous extensive experience of stick balancing. The research conforms to the ethical standards laid down in the 1964 Declaration of Helsinki, which protects research subjects, and was approved by the local ethical committee of Liguria Region (n. 222REG2015). Accordingly, all the subjects provided an informed consent before performing the experiments. The experiments were carried out at the Motor Learning, Assistive and Rehabilitation Robotics Lab of the Italian Institute of Technology (Genoa, Italy).
Experimental protocol. As a baseline, the subjects were initially asked to stay quietly upright, with parallel feet separated by about 0.2 m, the arms extended on the sides of the body, the eyes open required to fixate a point on a white wall at eye-height, on top of a force platform (AMTI, Watertown, Mass, USA) that recorded the different components of the ground reaction force and the CoP displacements in antero-posterior (AP) and medio-lateral (ML) directions. Sway movements of the body were recorded for 120 s, excluding the initial 20 s from the analysis. This is the basic, single balancing task.
In the dual balancing task the subjects stood with the feet in the same position, the elbows flexed at about 90 deg, holding the CIP-like device, consisting of a cylindrical wooden bar to be kept horizontal and aligned in the ML direction with the two hands. In the middle of the bar was mounted a ball bearing connected to a 1 m wooden stick that could rotate freely in the sagittal plane. The total weight of the CIP-like device is 0.375 kg. The task in this case was learning to balance the stick, namely avoiding the fall by keeping it approximately upright, while maintaining the feet fixed in the initial position, for the longest possible time interval. In order to achieve this goal the subjects learnt to move the CIP-like device back and forth, where the control variable is the force transmitted to the device as a function of the state of the stick (the tilt angle and the corresponding angular velocity). On this purpose, all the subjects spontaneously chose to keep their eyes fixed on the oscillating stick, in order to estimate its variable state. Although all the subjects had some experience in their childhood of balancing long sticks on the hand, none of them was a professional or amateur balancer. Therefore, they needed a suitable training until they mastered the task for at least 30 s continuously over consecutive trials. The practice time was different for the different subjects, typically at least 30 minutes, over one or more days if necessary. Practice was stopped after achieving reliable performance.
Kinematic data of the body inverted pendulum and of the CIP-like device were collected using an optoelectronic system (VICON, Oxford Metrics, Oxford, UK); markers were placed on ankles, knees, hips, shoulders, elbows, wrists and forehead of the participants and on the CIP-like device (two markers on the two extremities of the hand-held bar and one marker on the tip of the stick).
The electrical activation of different muscles of the legs/trunk/arms was also recorded by means of surface electrodes, although it was not specifically used for the present study. It could provide preliminary evidence for future extended studies.
Data Analysis. Markers position and force platform data were collected at a frequency of 100 Hz. Acquired signals were post-processed by a Savitzky-Golay low-pass filter (cut-off frequency of 10 Hz), which was also used to obtain velocities. After checking that movements of the CIP-like device in the vertical and medio-lateral directions were negligible and that the device was kept aligned with the body frame of reference (see Table 3), we focused the analysis only on the sagittal plane, i.e. on the AP direction.
From the CoP recorded through the force platform, we reconstructed the CoM of the human body according to subjects biomechanical features 64 and computed the related tilt angle (θ b ) of the body inverted pendulum (BIP). This angle coincides with the tilt angle of the CoM (θ com ) in the single balancing task, i.e. the VIP coincides with   www.nature.com/scientificreports www.nature.com/scientificreports/ the BIP in this case. In the dual task the VIP is slightly tilted slightly forward with respect to the BIP and we computed it from the combination of the CoMs of the body, the arm, and the CIP device, respectively.
From the markers data we obtained the position of the stick (x s ), the tilt angle and tilt angular velocity of the stick (θ θ  , s s ). We set the reference system in order to map anterior angular displacements as positive and posterior movements as negative angles.
Simulation experiments. The biomechanical parameters of the model sketched in Fig. 2 and the parameters of the two state-space intermittent feedback control models expressed by Eqs. 1, 2, 3, 4 are stored in Table 4.
The simulations were carried out with MATLAB by MathWorks, using the forward Euler method with a time step of 0.001 s.

Data availability
The authors confirm that all data supporting the findings of this study are available within this manuscript. The raw data are also available upon request.  Table 4. Biomechanical and control parameters for the simulation study.