Objective
Brain–computer interfaces (BCIs) could potentially be used to interact with pathological brain signals to intervene and ameliorate their effects in disease states. Here, we provide ...proof‐of‐principle of this approach by using a BCI to interpret pathological brain activity in patients with advanced Parkinson disease (PD) and to use this feedback to control when therapeutic deep brain stimulation (DBS) is delivered. Our goal was to demonstrate that by personalizing and optimizing stimulation in real time, we could improve on both the efficacy and efficiency of conventional continuous DBS.
Methods
We tested BCI‐controlled adaptive DBS (aDBS) of the subthalamic nucleus in 8 PD patients. Feedback was provided by processing of the local field potentials recorded directly from the stimulation electrodes. The results were compared to no stimulation, conventional continuous stimulation (cDBS), and random intermittent stimulation. Both unblinded and blinded clinical assessments of motor effect were performed using the Unified Parkinson's Disease Rating Scale.
Results
Motor scores improved by 66% (unblinded) and 50% (blinded) during aDBS, which were 29% (p = 0.03) and 27% (p = 0.005) better than cDBS, respectively. These improvements were achieved with a 56% reduction in stimulation time compared to cDBS, and a corresponding reduction in energy requirements (p < 0.001). aDBS was also more effective than no stimulation and random intermittent stimulation.
Interpretation
BCI‐controlled DBS is tractable and can be more efficient and efficacious than conventional continuous neuromodulation for PD. Ann Neurol 2013;74:449–457
Adaptive deep brain stimulation (aDBS) uses feedback from brain signals to guide stimulation. A recent acute trial of unilateral aDBS showed that aDBS can lead to substantial improvements in ...contralateral hemibody Unified Parkinson's Disease Rating Scale (UPDRS) motor scores and may be superior to conventional continuous DBS in Parkinson's disease (PD). We test whether potential benefits are retained with bilateral aDBS and in the face of concurrent medication.
We applied bilateral aDBS in 4 patients with PD undergoing DBS of the subthalamic nucleus. aDBS was delivered bilaterally with independent triggering of stimulation according to the amplitude of β activity at the corresponding electrode. Mean stimulation voltage was 3.0±0.1 volts. Motor assessments consisted of double-blinded video-taped motor UPDRS scores that included both limb and axial features.
UPDRS scores were 43% (p=0.04; Cohen's d=1.62) better with aDBS than without stimulation. Motor improvement with aDBS occurred despite an average time on stimulation (ToS) of only 45%. Levodopa was well tolerated during aDBS and led to further reductions in ToS.
Bilateral aDBS can improve both axial and limb symptoms and can track the need for stimulation across drug states.
The gradient technique is a promising tool with theoretical foundations based on the fundamental properties of MHD turbulence and turbulent reconnection. Its various incarnations use spectroscopic, ...synchrotron, and intensity data to trace the magnetic field and measure the media magnetization in terms of Alfven Mach number. We provide an analytical theory of gradient measurements and quantify the effects of averaging gradients along the line of sight and over the plane of the sky. We derive analytical expressions that relate the properties of gradient distribution with the Alfven Mach number \(M_A\). We show that these measurements can be combined with measures of sonic Mach number or line broadening to obtain the magnetic field strength. The corresponding technique has advantages to Davis-Chandrasekhar-Fermi way of obtaining the magnetic field strength.
We consider turbulent synchrotron emitting media that also exhibits Faraday rotation and provide a statistical description of synchrotron polarization fluctuations. In particular, we consider these ...fluctuations as a function of the spatial separation of the direction of measurements and as a function of wavelength for the same line-of-sight. On the basis of our general analytical approach, we introduce several measures that can be used to obtain the spectral slopes and correlation scales of both the underlying magnetic turbulence responsible for emission and the spectrum of the Faraday rotation fluctuations. We show the synergetic nature of these measures and discuss how the study can be performed using sparsely sampled interferometric data. We also discuss how additional characteristics of turbulence can be obtained, including the turbulence anisotropy, the three dimensional direction of the mean magnetic field. We consider both cases when the synchrotron emission and Faraday rotation regions coincide and when they are spatially separated. Appealing to our earlier study in Lazarian and Pogosyan (2012) we explain that our new results are applicable to a wide range of spectral indexes of relativistic electrons responsible for synchrotron emission. We expect wide application of our techniques both with existing synchrotron data sets as well as with big forthcoming data sets from LOFAR and SKA.
We introduce two new ways of obtaining the strength of plane-of-sky (POS) magnetic field by simultaneous use of spectroscopic Doppler-shifted lines and the information on magnetic field direction. ...The latter can be obtained either through polarization measurements or using the velocity gradient technique. We show the advantages that our techniques have compared to the traditional Davis-Chandrasekhar-Fermi (DCF) technique of estimating magnetic field strength from observations. The first technique that we describe in detail employs structure functions of velocity centroids and structure functions of Stokes parameters. We provide analytical expressions for obtaining magnetic field strength from observational data. We successfully test our results using synthetic observations obtained with results of MHD turbulence simulations. We measure velocity and magnetic field fluctuations at small scales using two, three and four point structure functions and compare the performance of these tools. We show that, unlike the DCF, our technique is capable of providing the detailed distribution of POS magnetic field and it can measure magnetic field strength in the presence of both velocity and magnetic field distortions arising from external shear and self-gravity. The second technique applies the velocity gradient technique to velocity channel maps in order to obtain the Alfven Mach number and uses the amplitudes of the gradients to obtain the sonic Mach number. The ratio of these two Mach numbers provides the intensity of magnetic field in the region contributing to the emission in the channel map. We test the technique and discuss obtaining the 3D distribution of POS galactic Magnetic field with it. We discuss the application of the second technique to synchrotron data.
We extend the velocity channel analysis (VCA), introduced by Lazarian & Pogosyan, of the intensity fluctuations in the velocity slices of position-position-velocity (PPV) spectroscopic data from ...Doppler broadened lines to study statistical anisotropy of the underlying velocity and density that arises in a turbulent medium from the presence of magnetic field. In particular, we study analytically how the anisotropy of the intensity correlation in the channel maps changes with the thickness of velocity channels. In agreement with the earlier VCA studies we find that the anisotropy in the thick channels reflects the anisotropy of the density field, while the relative contribution of density and velocity fluctuations to the thin velocity channels depends on the density spectral slope. We show that the anisotropies arising from Alfven, slow and fast magnetohydrodynamical modes are different, in particular, the anisotropy in PPV created by fast modes is opposite to that created by Alfven and slow modes, and this can be used to separate their contributions. We successfully compare our results with the recent numerical study of the PPV anisotropies measured with synthetic observations. We also extend our study to the medium with self-absorption as well as to the case of absorption lines. In addition, we demonstrate how the studies of anisotropy can be performed using interferometers.
We provide a theoretical description of synchrotron fluctuations arising from magnetic turbulence. We derive an expression that relates the correlation of synchrotron fluctuations for an arbitrary ...index of relativistic electrons with the correlations of squared fluctuations of magnetic field component perpendicular to the line of sight. The latter correlations we study assuming that the turbulence is axisymmetric. We obtain general relations valid for an arbitrary model of magnetic turbulence and analyse the relations for particular example of magnetic turbulence that is supported by numerical simulations. We predict that the synchrotron intensity fluctuations are anisotropic with larger correlation present along the direction of magnetic field. This anisotropy is dominated by the quadrupole component with the ratio between quadrupole and monopole parts being sensitive to the compressibility of underlying turbulence. Our work opens avenues for quantitative studies of magnetic turbulence in our galaxy and beyond using synchrotron emission. It also outlines the directions of how synchrotron foreground emission can be separated from cosmological signal, i.e. from CMB or highly redshifted HI emission. For the sake of completeness we also provide the expressions for the synchrotron polarization (Stocks parameters and their combinations) for the model of axisymmetric magnetic turbulence.
In this work we investigate the Principal Component Analysis (PCA) sensitivity to the velocity power spectrum in high opacity regimes of the interstellar medium (ISM). For our analysis we use ...synthetic Position-Position-Velocity (PPV) cubes of fractional Brownian motion (fBm) and magnetohydrodynamics (MHD) simulations, post processed to include radiative transfer effects from CO. We find that PCA analysis is very different from the tools based on the traditional power spectrum of PPV data cubes. Our major finding is that PCA is also sensitive to the phase information of PPV cubes and this allows PCA to detect the changes of the underlying velocity and density spectra at high opacities, where the spectral analysis of the maps provides the universal -3 spectrum in accordance with the predictions of Lazarian \& Pogosyan (2004) theory. This makes PCA potentially a valuable tool for studies of turbulence at high opacities provided that the proper gauging of the PCA index is made. The later, however, we found to be not easy, as the PCA results change in an irregular way for data with high sonic Mach numbers. This is in contrast to synthetic Brownian noise data used for velocity and density fields that show monotonic PCA behavior. We attribute this difference to the PCA's sensitivity to Fourier phase information.