A dynamic generative model can extract interpretable oscillatory components from multichannel neurophysiological recordings.

Modern neurophysiological recordings are performed using multichannel sensor arrays that are able to record activity in an increasingly high number of channels numbering in the 100's to 1000's. Often, underlying lower-dimensional patterns of activity are responsible for the observed dynamics, but these representations are difficult to reliably identify using existing methods that attempt to summarize multivariate relationships in a post-hoc manner from univariate analyses, or using current blind source separation methods. While such methods can reveal appealing patterns of activity, determining the number of components to include, assessing their statistical significance, and interpreting them requires extensive manual intervention and subjective judgement in practice. These difficulties with component selection and interpretation occur in large part because these methods lack a generative model for the underlying spatio-temporal dynamics. Here we describe a novel component analysis method anchored by a generative model where each source is described by a bio-physically inspired state space representation. The parameters governing this representation readily capture the oscillatory temporal dynamics of the components, so we refer to it as Oscillation Component Analysis (OCA). These parameters - the oscillatory properties, the component mixing weights at the sensors, and the number of oscillations - all are inferred in a data-driven fashion within a Bayesian framework employing an instance of the expectation maximization algorithm. We analyze high-dimensional electroencephalography and magnetoencephalography recordings from human studies to illustrate the potential utility of this method for neuroscience data.

Development and prospective validation of postoperative pain prediction from preoperative EHR data using attention-based set embeddings.

Preoperative knowledge of expected postoperative pain can help guide perioperative pain management and focus interventions on patients with the greatest risk of acute pain. However, current methods for predicting postoperative pain require patient and clinician input or laborious manual chart review and often do not achieve sufficient performance. We use routinely collected electronic health record data from a multicenter dataset of 234,274 adult non-cardiac surgical patients to develop a machine learning method which predicts maximum pain scores on the day of surgery and four subsequent days and validate this method in a prospective cohort. Our method, POPS, is fully automated and relies only on data available prior to surgery, allowing application in all patients scheduled for or considering surgery. Here we report that POPS achieves state-of-the-art performance and outperforms clinician predictions on all postoperative days when predicting maximum pain on the 0-10 NRS in prospective validation, though with degraded calibration. POPS is interpretable, identifying comorbidities that significantly contribute to postoperative pain based on patient-specific context, which can assist clinicians in mitigating cases of acute pain.

Eelbrain, a Python toolkit for time-continuous analysis with temporal response functions.

Even though human experience unfolds continuously in time, it is not strictly linear; instead, it entails cascading processes building hierarchical cognitive structures. For instance, during speech perception, humans transform a continuously varying acoustic signal into phonemes, words, and meaning, and these levels all have distinct but interdependent temporal structures. Time-lagged regression using temporal response functions (TRFs) has recently emerged as a promising tool for disentangling electrophysiological brain responses related to such complex models of perception. Here, we introduce the Eelbrain Python toolkit, which makes this kind of analysis easy and accessible. We demonstrate its use, using continuous speech as a sample paradigm, with a freely available EEG dataset of audiobook listening. A companion GitHub repository provides the complete source code for the analysis, from raw data to group-level statistics. More generally, we advocate a hypothesis-driven approach in which the experimenter specifies a hierarchy of time-continuous representations that are hypothesized to have contributed to brain responses, and uses those as predictor variables for the electrophysiological signal. This is analogous to a multiple regression problem, but with the addition of a time dimension. TRF analysis decomposes the brain signal into distinct responses associated with the different predictor variables by estimating a multivariate TRF (mTRF), quantifying the influence of each predictor on brain responses as a function of time(-lags). This allows asking two questions about the predictor variables: (1) Is there a significant neural representation corresponding to this predictor variable? And if so, (2) what are the temporal characteristics of the neural response associated with it? Thus, different predictor variables can be systematically combined and evaluated to jointly model neural processing at multiple hierarchical levels. We discuss applications of this approach, including the potential for linking algorithmic/representational theories at different cognitive levels to brain responses through computational models with appropriate linking hypotheses.

Switching state-space modeling of neural signal dynamics.

Linear parametric state-space models are a ubiquitous tool for analyzing neural time series data, providing a way to characterize the underlying brain dynamics with much greater statistical efficiency than non-parametric data analysis approaches. However, neural time series data are frequently time-varying, exhibiting rapid changes in dynamics, with transient activity that is often the key feature of interest in the data. Stationary methods can be adapted to time-varying scenarios by employing fixed-duration windows under an assumption of quasi-stationarity. But time-varying dynamics can be explicitly modeled by switching state-space models, i.e., by using a pool of state-space models with different dynamics selected by a probabilistic switching process. Unfortunately, exact solutions for state inference and parameter learning with switching state-space models are intractable. Here we revisit a switching state-space model inference approach first proposed by Ghahramani and Hinton. We provide explicit derivations for solving the inference problem iteratively after applying a variational approximation on the joint posterior of the hidden states and the switching process. We introduce a novel initialization procedure using an efficient leave-one-out strategy to compare among candidate models, which significantly improves performance compared to the existing method that relies on deterministic annealing. We then utilize this state inference solution within a generalized expectation-maximization algorithm to estimate model parameters of the switching process and the linear state-space models with dynamics potentially shared among candidate models. We perform extensive simulations under different settings to benchmark performance against existing switching inference methods and further validate the robustness of our switching inference solution outside the generative switching model class. Finally, we demonstrate the utility of our method for sleep spindle detection in real recordings, showing how switching state-space models can be used to detect and extract transient spindles from human sleep electroencephalograms in an unsupervised manner.

Altered directional functional connectivity underlies post-stroke cognitive recovery.

Cortical ischaemic strokes result in cognitive deficits depending on the area of the affected brain. However, we have demonstrated that difficulties with attention and processing speed can occur even with small subcortical infarcts. Symptoms appear independent of lesion location, suggesting they arise from generalized disruption of cognitive networks. Longitudinal studies evaluating directional measures of functional connectivity in this population are lacking. We evaluated six patients with minor stroke exhibiting cognitive impairment 6-8 weeks post-infarct and four age-similar controls. Resting-state magnetoencephalography data were collected. Clinical and imaging evaluations of both groups were repeated 6- and 12 months later. Network Localized Granger Causality was used to determine differences in directional connectivity between groups and across visits, which were correlated with clinical performance. Directional connectivity patterns remained stable across visits for controls. After the stroke, inter-hemispheric connectivity between the frontoparietal cortex and the non-frontoparietal cortex significantly increased between visits 1 and 2, corresponding to uniform improvement in reaction times and cognitive scores. Initially, the majority of functional links originated from non-frontal areas contralateral to the lesion, connecting to ipsilesional brain regions. By visit 2, inter-hemispheric connections, directed from the ipsilesional to the contralesional cortex significantly increased. At visit 3, patients demonstrating continued favourable cognitive recovery showed less reliance on these inter-hemispheric connections. These changes were not observed in those without continued improvement. Our findings provide supporting evidence that the neural basis of early post-stroke cognitive dysfunction occurs at the network level, and continued recovery correlates with the evolution of inter-hemispheric connectivity.

Association of Intraoperative Opioid Administration With Postoperative Pain and Opioid Use.

Importance: Opioids administered to treat postsurgical pain are a major contributor to the opioid crisis, leading to chronic use in a considerable proportion of patients. Initiatives promoting opioid-free or opioid-sparing modalities of perioperative pain management have led to reduced opioid administration in the operating room, but this reduction could have unforeseen detrimental effects in terms of postoperative pain outcomes, as the relationship between intraoperative opioid usage and later opioid requirements is not well understood. Objective: To characterize the association between intraoperative opioid usage and postoperative pain and opioid requirements. Design, Setting, and Participants: This retrospective cohort study evaluated electronic health record data from a quaternary care academic medical center (Massachusetts General Hospital) for adult patients who underwent noncardiac surgery with general anesthesia from April 2016 to March 2020. Patients who underwent cesarean surgery, received regional anesthesia, received opioids other than fentanyl or hydromorphone, were admitted to the intensive care unit, or who died intraoperatively were excluded. Statistical models were fitted on the propensity weighted data set to characterize the effect of intraoperative opioid exposures on primary and secondary outcomes. Data were analyzed from December 2021 to October 2022. Exposures: Intraoperative fentanyl and intraoperative hydromorphone average effect site concentration estimated using pharmacokinetic/pharmacodynamic models. Main Outcomes and Measures: The primary study outcomes were the maximal pain score during the postanesthesia care unit (PACU) stay and the cumulative opioid dose, quantified in morphine milligram equivalents (MME), administered during the PACU stay. Medium- and long-term outcomes associated with pain and opioid dependence were also evaluated. Results: The study cohort included a total of 61 249 individuals undergoing surgery (mean [SD] age, 55.44 [17.08] years; 32 778 [53.5%] female). Increased intraoperative fentanyl and intraoperative hydromorphone were both associated with reduced maximum pain scores in the PACU. Both exposures were also associated with a reduced probability and reduced total dosage of opioid administration in the PACU. In particular, increased fentanyl administration was associated with lower frequency of uncontrolled pain; a decrease in new chronic pain diagnoses reported at 3 months; fewer opioid prescriptions at 30, 90, and 180 days; and decreased new persistent opioid use, without significant increases in adverse effects. Conclusions and Relevance: Contrary to prevailing trends, reduced opioid administration during surgery may have the unintended outcome of increasing postoperative pain and opioid consumption. Conversely, improvements in long-term outcomes might be achieved by optimizing opioid administration during surgery.

Non-Asymptotic Guarantees for Reliable Identification of Granger Causality via the LASSO.

Granger causality is among the widely used data-driven approaches for causal analysis of time series data with applications in various areas including economics, molecular biology, and neuroscience. Two of the main challenges of this methodology are: 1) over-fitting as a result of limited data duration, and 2) correlated process noise as a confounding factor, both leading to errors in identifying the causal influences. Sparse estimation via the LASSO has successfully addressed these challenges for parameter estimation. However, the classical statistical tests for Granger causality resort to asymptotic analysis of ordinary least squares, which require long data duration to be useful and are not immune to confounding effects. In this work, we address this disconnect by introducing a LASSO-based statistic and studying its non-asymptotic properties under the assumption that the true models admit sparse autoregressive representations. We establish fundamental limits for reliable identification of Granger causal influences using the proposed LASSO-based statistic. We further characterize the false positive error probability and test power of a simple thresholding rule for identifying Granger causal effects and provide two methods to set the threshold in a data-driven fashion. We present simulation studies and application to real data to compare the performance of our proposed method to ordinary least squares and existing LASSO-based methods in detecting Granger causal influences, which corroborate our theoretical results.

NLGC: Network localized Granger causality with application to MEG directional functional connectivity analysis.

Identifying the directed connectivity that underlie networked activity between different cortical areas is critical for understanding the neural mechanisms behind sensory processing. Granger causality (GC) is widely used for this purpose in functional magnetic resonance imaging analysis, but there the temporal resolution is low, making it difficult to capture the millisecond-scale interactions underlying sensory processing. Magnetoencephalography (MEG) has millisecond resolution, but only provides low-dimensional sensor-level linear mixtures of neural sources, which makes GC inference challenging. Conventional methods proceed in two stages: First, cortical sources are estimated from MEG using a source localization technique, followed by GC inference among the estimated sources. However, the spatiotemporal biases in estimating sources propagate into the subsequent GC analysis stage, may result in both false alarms and missing true GC links. Here, we introduce the Network Localized Granger Causality (NLGC) inference paradigm, which models the source dynamics as latent sparse multivariate autoregressive processes and estimates their parameters directly from the MEG measurements, integrated with source localization, and employs the resulting parameter estimates to produce a precise statistical characterization of the detected GC links. We offer several theoretical and algorithmic innovations within NLGC and further examine its utility via comprehensive simulations and application to MEG data from an auditory task involving tone processing from both younger and older participants. Our simulation studies reveal that NLGC is markedly robust with respect to model mismatch, network size, and low signal-to-noise ratio, whereas the conventional two-stage methods result in high false alarms and mis-detections. We also demonstrate the advantages of NLGC in revealing the cortical network-level characterization of neural activity during tone processing and resting state by delineating task- and age-related connectivity changes.

Operational Parameters for Sub-Nano Tesla Field Resolution of PHMR Sensors in Harsh Environments.

The resolution of planar-Hall magnetoresistive (PHMR) sensors was investigated in the frequency range from 0.5 Hz to 200 Hz in terms of its sensitivity, average noise level, and detectivity. Analysis of the sensor sensitivity and voltage noise response was performed by varying operational parameters such as sensor geometrical architectures, sensor configurations, sensing currents, and temperature. All the measurements of PHMR sensors were carried out under both constant current (CC) and constant voltage (CV) modes. In the present study, Barkhausen noise was revealed in 1/f noise component and found less significant in the PHMR sensor configuration. Under measured noise spectral density at optimized conditions, the best magnetic field detectivity was achieved better than 550 pT/√Hz at 100 Hz and close to 1.1 nT/√Hz at 10 Hz for a tri-layer multi-ring PHMR sensor in an unshielded environment. Furthermore, the promising feasibility and possible routes for further improvement of the sensor resolution are discussed.

Bridge Resistance Compensation for Noise Reduction in a Self-Balanced PHMR Sensor.

Advanced microelectromechanical system (MEMS) magnetic field sensor applications demand ultra-high detectivity down to the low magnetic fields. To enhance the detection limit of the magnetic sensor, a resistance compensator integrated self-balanced bridge type sensor was devised for low-frequency noise reduction in the frequency range of 0.5 Hz to 200 Hz. The self-balanced bridge sensor was a NiFe (10 nm)/IrMn (10 nm) bilayer structure in the framework of planar Hall magnetoresistance (PHMR) technology. The proposed resistance compensator integrated with a self-bridge sensor architecture presented a compact and cheaper alternative to marketable MEMS MR sensors, adjusting the offset voltage compensation at the wafer level, and led to substantial improvement in the sensor noise level. Moreover, the sensor noise components of electronic and magnetic origin were identified by measuring the sensor noise spectral density as a function of temperature and operating power. The lowest achievable noise in this device architecture was estimated at ~3.34 nV/Hz at 100 Hz.

Neuro-current response functions: A unified approach to MEG source analysis under the continuous stimuli paradigm.

Characterizing the neural dynamics underlying sensory processing is one of the central areas of investigation in systems and cognitive neuroscience. Neuroimaging techniques such as magnetoencephalography (MEG) and Electroencephalography (EEG) have provided significant insights into the neural processing of continuous stimuli, such as speech, thanks to their high temporal resolution. Existing work in the context of auditory processing suggests that certain features of speech, such as the acoustic envelope, can be used as reliable linear predictors of the neural response manifested in M/EEG. The corresponding linear filters are referred to as temporal response functions (TRFs). While the functional roles of specific components of the TRF are well-studied and linked to behavioral attributes such as attention, the cortical origins of the underlying neural processes are not as well understood. In this work, we address this issue by estimating a linear filter representation of cortical sources directly from neuroimaging data in the context of continuous speech processing. To this end, we introduce Neuro-Current Response Functions (NCRFs), a set of linear filters, spatially distributed throughout the cortex, that predict the cortical currents giving rise to the observed ongoing MEG (or EEG) data in response to continuous speech. NCRF estimation is cast within a Bayesian framework, which allows unification of the TRF and source estimation problems, and also facilitates the incorporation of prior information on the structural properties of the NCRFs. To generalize this analysis to M/EEG recordings which lack individual structural magnetic resonance (MR) scans, NCRFs are extended to free-orientation dipoles and a novel regularizing scheme is put forward to lessen reliance on fine-tuned coordinate co-registration. We present a fast estimation algorithm, which we refer to as the Champ-Lasso algorithm, by leveraging recent advances in optimization, and demonstrate its utility through application to simulated and experimentally recorded MEG data under auditory experiments. Our simulation studies reveal significant improvements over existing methods that typically operate in a two-stage fashion, in terms of spatial resolution, response function reconstruction, and recovering dipole orientations. The analysis of experimentally-recorded MEG data without MR scans corroborates existing findings, but also delineates the distinct cortical distribution of the underlying neural processes at high spatiotemporal resolution. In summary, we provide a principled modeling and estimation paradigm for MEG source analysis tailored to extracting the cortical origin of electrophysiological responses to continuous stimuli.

Impact of Gd Doping on the Microstructural-, Electronic-, Magnetotransport, and Magnetic Properties of Pr0.6Gd0.1Sr0.3MnO3 Nanomanganite.

Gadolinium doping effects on structural-, electronic transport and magnetic properties of Pr0.6Gd0.1Sr0.3MnO3 (PGSMO) nanomanganite (40-65 nm) have been investigated. We have considered the core-shell structure of our PGSMO nanoparticles, which can explain its intriguing magnetic and transport properties. It is reported 65 nm sample exhibits higher metal to insulator transition temperature, TP (∼215 K) and less resistive behavior compare to its nano counterparts. Enhanced low-field magnetoresistance (LFMR) is observed with a reduction of particle size; however maximum change in magnetoresistance (MR), ∼72% is recorded for 65 nm sample around TP at 60 kOe field. Low-temperature magnetization data reveals the existence of glassy phase in the compound. It is revealed that substitution of optimal Gd doping at A site on nanometric scale promotes exotic magnetic and transport properties on nanoscale without adding any external physical perturbation or introducing any significant lattice distortion.

Correlation of dielectric, electrical and magnetic properties near the magnetic phase transition temperature of cobalt zinc ferrite.

Multiferroic composite structures, i.e., composites of magnetostrictive and piezoelectric materials, can be envisioned to achieve the goal of strong room-temperature ME coupling for real practical device applications. Magnetic materials with high magnetostriction, high Néel temperature (TN), high resistivity and large magnetization are required to observe high ME coupling in composite structures. In continuation of our investigations on suitable magnetic candidates for multiferroic composite structures, we have studied the crystal structure, dielectric, transport, and magnetic properties of Co0.65Zn0.35Fe2O4 (CZFO). Rietveld refinement of X-ray diffraction patterns confirms the phase purity with a cubic crystal structure with the (Fd3[combining macron]m) space group; however, we have found a surprisingly large magneto-dielectric anomaly at the Néel temperature, unexpected for a cubic structure. The presence of mixed valences of Fe2+/Fe3+ cations is probed by X-ray photoelectron spectroscopy (XPS), which supports the catonic ordering-mediated large dielectric response. Large dielectric permittivity dispersion with a broad anomaly is observed in the vicinity of the magnetic phase transition temperature (TN) of CZFO suggesting a strong correlation between dielectric and magnetic properties. The evidence of strong spin-polaron coupling has been established from temperature dependent dielectric, ac conductivity and magnetization studies. The ferrimagnetic-paramagnetic phase transition of CZFO has been found at ∼640 K, which is well above room temperature. CZFO exhibits low loss tangent, a high dielectric constant, large magnetization with soft magnetic behavior and magnetodielectric coupling above room temperature, elucidating the possible potential candidates for multiferroic composite structures as well as for multifunctional and spintronics device applications.