Generic low-density corrections to the equation of state of chain molecules with repulsive intermolecular forces.

Molecular-based equations of state for describing the thermodynamics of chain molecules are often based on mean-field like arguments that reduce the problem of describing the interactions between chains to a simpler one involving only nonbonded monomers. While for dense liquids such arguments are known to work well, at low density they are typically less appropriate due to an incomplete description of the effect of chain connectivity on the local environment of the chains' monomer segments. To address this issue, we develop three semi-empirical approaches that significantly improve the thermodynamic description of chain molecules at low density. The approaches are developed for chain molecules with repulsive intermolecular forces; therefore, they could be used as reference models for developing equations of the state of real fluids based on perturbation theory. All three approaches are extensions of Wertheim's first-order thermodynamic perturbation theory (TPT1) for polymerization. The first model, referred to as TPT1-v, incorporates a second-virial correction that is scaled to zero at liquid-like densities. The second model, referred to as TPT1-y, introduces a Helmholtz-energy contribution to account for correlations between next-nearest-neighbor segments within chain molecules. The third approach, called TPT-E, directly modifies TPT1 without utilizing an additional Helmholtz energy contribution. By employing TPT1 at the core of these approaches, we ensure an accurate description of mixtures and enable a seamless extension from chains of tangentially bonded hard-sphere segments of equal size to hetero-segmented chains, fused chains, and chains of soft repulsive segments (which are influenced by temperature). The low-density corrections implemented in TPT1 are designed to preserve these good characteristics, as confirmed through comparisons with novel molecular simulation results for the pressure of various chain fluids. TPT1-v exhibits excellent transferability across different chain types, but it relies on knowing the second virial coefficient of the chain molecules, which is non-trivial to obtain and determined here using Monte Carlo simulation. The TPT1-y model, on the other hand, achieves comparable accuracy to TPT1-v while being fully predictive, requiring no input besides the geometry of the chain molecules.

Transferable Anisotropic Mie Potential Force Field for Alkanediols.

The development of force fields for polyfunctional molecules, such as alkanediols, requires a careful account of different average intramolecular conformations for gas states compared to dense liquid states, where intra- and intermolecular hydrogen bonds compete. In the present work, the transferable anisotropic Mie (TAMie) potential is extended to 1,n-alkanediols. Using the convention that intramolecular nonbonded interactions up to and including the third neighbor are excluded, all force field parameters developed previously for 1-alcohols were transferred to 1,5-pentanediol and beyond, with good agreement with experimental phase equilibrium data. To obtain trans-gauche ratios of 1,2-ethanediol and 1,3-propanediol that are consistent with experimental results, the propensities for intra- and intermolecular hydrogen bonds had to be balanced. This was achieved by parameterizing the intramolecular dihedral energy functions governing the O-C-C-O and O-C-C-C angles while intramolecular charge-charge interactions were active. All partial charges belonging to a functional group are collected in a charge group and all interactions among two charge groups are evaluated even if they are separated by less than three bonds. With this approach, it is possible to apply the nonbonded parameters from 1-alcohols to alkanediols without further refinement. The agreement with experimental phase equilibrium and shear viscosity data is of similar quality as for the 1-alcohols and the trans-gauche ratio agrees with literature results from spectroscopic measurements and ab initio calculations.

Brain-body states embody complex temporal dynamics.

We propose a computational framework for high-dimensional brain-body states as transient embodiments of nested internal and external dynamics governed by interoception. Unifying recent theoretical work, we suggest ways to reduce arbitrary state complexity to an observable number of features in order to accurately predict and intervene in pathological trajectories.

Altered cortical synaptic lipid signaling leads to intermediate phenotypes of mental disorders.

Excitation/inhibition (E/I) balance plays important roles in mental disorders. Bioactive phospholipids like lysophosphatidic acid (LPA) are synthesized by the enzyme autotaxin (ATX) at cortical synapses and modulate glutamatergic transmission, and eventually alter E/I balance of cortical networks. Here, we analyzed functional consequences of altered E/I balance in 25 human subjects induced by genetic disruption of the synaptic lipid signaling modifier PRG-1, which were compared to 25 age and sex matched control subjects. Furthermore, we tested therapeutic options targeting ATX in a related mouse line. Using EEG combined with TMS in an instructed fear paradigm, neuropsychological analysis and an fMRI based episodic memory task, we found intermediate phenotypes of mental disorders in human carriers of a loss-of-function single nucleotide polymorphism of PRG-1 (PRG-1R345T/WT). Prg-1R346T/WT animals phenocopied human carriers showing increased anxiety, a depressive phenotype and lower stress resilience. Network analysis revealed that coherence and phase-amplitude coupling were altered by PRG-1 deficiency in memory related circuits in humans and mice alike. Brain oscillation phenotypes were restored by inhibtion of ATX in Prg-1 deficient mice indicating an interventional potential for mental disorders.

Predicting Solvation Free Energies from the Minnesota Solvation Database Using Classical Density Functional Theory Based on the PC-SAFT Equation of State.

We critically assess the capabilities of classical density functional theory (DFT) based on the perturbed-chain statistical associating fluid theory (PC-SAFT) equation of state to predict the solvation free energies of small molecules in various hydrocarbon solvents. We compare DFT results with experimental data from the Minnesota solvation database and utilize statistical methods to analyze the accuracy of our approach, as well as its weaknesses. The mean absolute error of the solvation free energies is 3.7 kJ mol-1 for n-alkane solvents, ranging from pentane to hexadecane, with 473 solute-solvent systems. For solvents consisting of cyclic hydrocarbons (cyclohexane, benzene, toluene, and ethylbenzene) with 245 solute-solvent systems, we report a slightly larger mean absolute error of 4.2 kJ mol-1. We identify three possible sources of errors: (i) the neglect of solute-solvent and solvent-solvent Coulomb interactions, which limits the applicability of PC-SAFT DFT to nonpolar and weakly polar molecules; (ii) the solute's Lennard-Jones parameters supplied by the general AMBER force field, which are not parametrized toward solvation free energies; and (iii) the application of the Lorentz-Berthelot combining rules to the dispersive interactions between a segment of the PC-SAFT solvent and a Lennard-Jones interaction site of the solute. The approach is more accurate than standard implementations of phenomenological models in common chemistry software packages, which exhibit mean absolute errors larger than 9.12 kJ mol-1, even though newer phenomenological models achieve a mean absolute error of about 2 kJ mol-1. PC-SAFT DFT is more computationally efficient than state of the art explicit molecular simulations in combination with free energy perturbation methods. It is predictive with respect to solvation free energies, i.e., the input for the model is the (element-specific) molecular force field, the solute configuration from molecular dynamics simulations, and the (substance-specific) PC-SAFT parameters. The PC-SAFT parametrization uses pure-component data and does not require experimental solvation free energies. The PC-SAFT equation of state, without applying a DFT formalism, can also be used to calculate solvation free energies, provided that the PC-SAFT parameters for the solute are available. A large number of substances was recently parametrized by members of our group (Esper, T.; Bauer, G.; Rehner, P.; Gross, J. Ind. Eng. Chem. Res. 2023, 62), which enables a comparison to the DFT approach for 103 substances. Accurate results are obtained from the PC-SAFT equation of state with an MAE below 2.51 kJ mol-1. The DFT approach does not require PC-SAFT parameters for the solute and can be applied to all solutes that can be represented by the molecular force field.

Comparing the performance of beamformer algorithms in estimating orientations of neural sources.

The efficacy of transcranial electric stimulation (tES) to effectively modulate neuronal activity depends critically on the spatial orientation of the targeted neuronal population. Therefore, precise estimation of target orientation is of utmost importance. Different beamforming algorithms provide orientation estimates; however, a systematic analysis of their performance is still lacking. For fixed brain locations, EEG and MEG data from sources with randomized orientations were simulated. The orientation was then estimated (1) with an EEG and (2) with a combined EEG-MEG approach. Three commonly used beamformer algorithms were evaluated with respect to their abilities to estimate the correct orientation: Unit-Gain (UG), Unit-Noise-Gain (UNG), and Array-Gain (AG) beamformer. Performance depends on the signal-to-noise ratios for the modalities and on the chosen beamformer. Overall, the UNG and AG beamformers appear as the most reliable. With increasing noise, the UG estimate converges to a vector determined by the leadfield, thus leading to insufficient orientation estimates.

A Systematic Evaluation of Machine Learning-Based Biomarkers for Major Depressive Disorder.

Importance: Biological psychiatry aims to understand mental disorders in terms of altered neurobiological pathways. However, for one of the most prevalent and disabling mental disorders, major depressive disorder (MDD), no informative biomarkers have been identified. Objective: To evaluate whether machine learning (ML) can identify a multivariate biomarker for MDD. Design, Setting, and Participants: This study used data from the Marburg-Münster Affective Disorders Cohort Study, a case-control clinical neuroimaging study. Patients with acute or lifetime MDD and healthy controls aged 18 to 65 years were recruited from primary care and the general population in Münster and Marburg, Germany, from September 11, 2014, to September 26, 2018. The Münster Neuroimaging Cohort (MNC) was used as an independent partial replication sample. Data were analyzed from April 2022 to June 2023. Exposure: Patients with MDD and healthy controls. Main Outcome and Measure: Diagnostic classification accuracy was quantified on an individual level using an extensive ML-based multivariate approach across a comprehensive range of neuroimaging modalities, including structural and functional magnetic resonance imaging and diffusion tensor imaging as well as a polygenic risk score for depression. Results: Of 1801 included participants, 1162 (64.5%) were female, and the mean (SD) age was 36.1 (13.1) years. There were a total of 856 patients with MDD (47.5%) and 945 healthy controls (52.5%). The MNC replication sample included 1198 individuals (362 with MDD [30.1%] and 836 healthy controls [69.9%]). Training and testing a total of 4 million ML models, mean (SD) accuracies for diagnostic classification ranged between 48.1% (3.6%) and 62.0% (4.8%). Integrating neuroimaging modalities and stratifying individuals based on age, sex, treatment, or remission status does not enhance model performance. Findings were replicated within study sites and also observed in structural magnetic resonance imaging within MNC. Under simulated conditions of perfect reliability, performance did not significantly improve. Analyzing model errors suggests that symptom severity could be a potential focus for identifying MDD subgroups. Conclusion and Relevance: Despite the improved predictive capability of multivariate compared with univariate neuroimaging markers, no informative individual-level MDD biomarker-even under extensive ML optimization in a large sample of diagnosed patients-could be identified.

Intact Mismatch Negativity Responses in Clinical High Risk for Psychosis and First-Episode Psychosis: Evidence From Source-Reconstructed Event-Related Fields and Time-Frequency Data.

BACKGROUND: This study examined whether mismatch negativity (MMN) responses are impaired in participants at clinical high risk for psychosis (CHR-P) and patients with first-episode psychosis (FEP) and whether MMN deficits predict clinical outcomes in CHR-Ps. METHODS: Magnetoencephalography data were collected during a duration-deviant MMN paradigm for a group of 116 CHR-P participants, 33 FEP patients (15 antipsychotic-naïve), clinical high risk negative group (n = 38) with substance abuse and affective disorder, and 49 healthy control participants. Analysis of group differences of source-reconstructed event-related fields as well as time-frequency and intertrial phase coherence focused on the bilateral Heschl's gyri and bilateral superior temporal gyri. RESULTS: Significant magnetic MMN responses were found across participants in the bilateral Heschl's gyri and bilateral superior temporal gyri. However, MMN amplitude as well as time-frequency and intertrial phase coherence responses were intact in CHR-P participants and FEP patients compared with healthy control participants. Furthermore, MMN deficits were not related to persistent attenuated psychotic symptoms or transitions to psychosis in CHR-P participants. CONCLUSIONS: Our data suggest that magnetic MMN responses in magnetoencephalography data are not impaired in early-stage psychosis and may not predict clinical outcomes in CHR-P participants.

Normative tDCS over V5 and FEF reveals practice-induced modulation of extraretinal smooth pursuit mechanisms, but no specific stimulation effect.

The neural networks subserving smooth pursuit eye movements (SPEM) provide an ideal model for investigating the interaction of sensory processing and motor control during ongoing movements. To better understand core plasticity aspects of sensorimotor processing for SPEM, normative sham, anodal or cathodal transcranial direct current stimulation (tDCS) was applied over visual area V5 and frontal eye fields (FEF) in sixty healthy participants. The identical within-subject paradigm was used to assess SPEM modulations by practice. While no specific tDCS effects were revealed, within- and between-session practice effects indicate plasticity of top-down extraretinal mechanisms that mainly affect SPEM in the absence of visual input and during SPEM initiation. To explore the potential of tDCS effects, individual electric field simulations were computed based on calibrated finite element head models and individual functional localization of V5 and FEF location (using functional MRI) and orientation (using combined EEG/MEG) was conducted. Simulations revealed only limited electric field target intensities induced by the applied normative tDCS montages but indicate the potential efficacy of personalized tDCS for the modulation of SPEM. In sum, results indicate the potential susceptibility of extraretinal SPEM control to targeted external neuromodulation (e.g., personalized tDCS) and intrinsic learning protocols.

Neural responses to natural and enhanced speech edges in children with and without dyslexia.

Sensory-neural studies indicate that children with developmental dyslexia show impairments in processing acoustic speech envelope information. Prior studies suggest that this arises in part from reduced sensory sensitivity to amplitude rise times (ARTs or speech "edges") in the envelope, accompanied by less accurate neural encoding of low-frequency envelope information. Accordingly, enhancing these characteristics of the speech envelope may enhance neural speech processing in children with dyslexia. Here we applied an envelope modulation enhancement (EME) algorithm to a 10-min story read in child-directed speech (CDS), enhancing ARTs and also enhancing low-frequency envelope information. We compared neural speech processing (as measured using MEG) for the EME story with the same story read in natural CDS for 9-year-old children with and without dyslexia. The EME story affected neural processing in the power domain for children with dyslexia, particularly in the delta band (0.5-4 Hz) in the superior temporal gyrus. This may suggest that prolonged experience with EME speech could ameliorate some of the impairments shown in natural speech processing by children with dyslexia.

CutFEM forward modeling for EEG source analysis.

Introduction: Source analysis of Electroencephalography (EEG) data requires the computation of the scalp potential induced by current sources in the brain. This so-called EEG forward problem is based on an accurate estimation of the volume conduction effects in the human head, represented by a partial differential equation which can be solved using the finite element method (FEM). FEM offers flexibility when modeling anisotropic tissue conductivities but requires a volumetric discretization, a mesh, of the head domain. Structured hexahedral meshes are easy to create in an automatic fashion, while tetrahedral meshes are better suited to model curved geometries. Tetrahedral meshes, thus, offer better accuracy but are more difficult to create. Methods: We introduce CutFEM for EEG forward simulations to integrate the strengths of hexahedra and tetrahedra. It belongs to the family of unfitted finite element methods, decoupling mesh and geometry representation. Following a description of the method, we will employ CutFEM in both controlled spherical scenarios and the reconstruction of somatosensory-evoked potentials. Results: CutFEM outperforms competing FEM approaches with regard to numerical accuracy, memory consumption, and computational speed while being able to mesh arbitrarily touching compartments. Discussion: CutFEM balances numerical accuracy, computational efficiency, and a smooth approximation of complex geometries that has previously not been available in FEM-based EEG forward modeling.

DISCOVER-EEG: an open, fully automated EEG pipeline for biomarker discovery in clinical neuroscience.

Biomarker discovery in neurological and psychiatric disorders critically depends on reproducible and transparent methods applied to large-scale datasets. Electroencephalography (EEG) is a promising tool for identifying biomarkers. However, recording, preprocessing, and analysis of EEG data is time-consuming and researcher-dependent. Therefore, we developed DISCOVER-EEG, an open and fully automated pipeline that enables easy and fast preprocessing, analysis, and visualization of resting state EEG data. Data in the Brain Imaging Data Structure (BIDS) standard are automatically preprocessed, and physiologically meaningful features of brain function (including oscillatory power, connectivity, and network characteristics) are extracted and visualized using two open-source and widely used Matlab toolboxes (EEGLAB and FieldTrip). We tested the pipeline in two large, openly available datasets containing EEG recordings of healthy participants and patients with a psychiatric condition. Additionally, we performed an exploratory analysis that could inspire the development of biomarkers for healthy aging. Thus, the DISCOVER-EEG pipeline facilitates the aggregation, reuse, and analysis of large EEG datasets, promoting open and reproducible research on brain function.

Modulatory dynamics of periodic and aperiodic activity in respiration-brain coupling.

Bodily rhythms such as respiration are increasingly acknowledged to modulate neural oscillations underlying human action, perception, and cognition. Conversely, the link between respiration and aperiodic brain activity - a non-oscillatory reflection of excitation-inhibition (E:I) balance - has remained unstudied. Aiming to disentangle potential respiration-related dynamics of periodic and aperiodic activity, we applied recently developed algorithms of time-resolved parameter estimation to resting-state MEG and EEG data from two labs (total N = 78 participants). We provide evidence that fluctuations of aperiodic brain activity (1/f slope) are phase-locked to the respiratory cycle, which suggests that spontaneous state shifts of excitation-inhibition balance are at least partly influenced by peripheral bodily signals. Moreover, differential temporal dynamics in their coupling to non-oscillatory and oscillatory activity raise the possibility of a functional distinction in the way each component is related to respiration. Our findings highlight the role of respiration as a physiological influence on brain signalling.

Predictive coordination of breathing during intra-personal speaking and listening.

It has long been known that human breathing is altered during listening and speaking compared to rest: during speaking, inhalation depth is adjusted to the air volume required for the upcoming utterance. During listening, inhalation is temporally aligned to inhalation of the speaker. While evidence for the former is relatively strong, it is virtually absent for the latter. We address both phenomena using recordings of speech envelope and respiration in 30 participants during 14 min of speaking and listening to one's own speech. First, we show that inhalation depth is positively correlated with the total power of the speech envelope in the following utterance. Second, we provide evidence that inhalation during listening to one's own speech is significantly more likely at time points of inhalation during speaking. These findings are compatible with models that postulate alignment of internal forward models of interlocutors with the aim to facilitate communication.

Spatiotemporal dynamics characterise spectral connectivity profiles of continuous speaking and listening.

Speech production and perception are fundamental processes of human cognition that both rely on intricate processing mechanisms that are still poorly understood. Here, we study these processes by using magnetoencephalography (MEG) to comprehensively map connectivity of regional brain activity within the brain and to the speech envelope during continuous speaking and listening. Our results reveal not only a partly shared neural substrate for both processes but also a dissociation in space, delay, and frequency. Neural activity in motor and frontal areas is coupled to succeeding speech in delta band (1 to 3 Hz), whereas coupling in the theta range follows speech in temporal areas during speaking. Neural connectivity results showed a separation of bottom-up and top-down signalling in distinct frequency bands during speaking. Here, we show that frequency-specific connectivity channels for bottom-up and top-down signalling support continuous speaking and listening. These findings further shed light on the complex interplay between different brain regions involved in speech production and perception.

Breathing in waves: Understanding respiratory-brain coupling as a gradient of predictive oscillations.

Breathing plays a crucial role in shaping perceptual and cognitive processes by regulating the strength and synchronisation of neural oscillations. Numerous studies have demonstrated that respiratory rhythms govern a wide range of behavioural effects across cognitive, affective, and perceptual domains. Additionally, respiratory-modulated brain oscillations have been observed in various mammalian models and across diverse frequency spectra. However, a comprehensive framework to elucidate these disparate phenomena remains elusive. In this review, we synthesise existing findings to propose a neural gradient of respiratory-modulated brain oscillations and examine recent computational models of neural oscillations to map this gradient onto a hierarchical cascade of precision-weighted prediction errors. By deciphering the computational mechanisms underlying respiratory control of these processes, we can potentially uncover new pathways for understanding the link between respiratory-brain coupling and psychiatric disorders.

Tinnitus and Multimodal Cortical Interaction. / Tinnitus und multimodale kortikale Interaktion.

The term of subjective tinnitus is used to describe a perceived noise without an external sound source. Therefore, it seems to be obvious that tinnitus can be understood as purely auditory, sensory problem. From a clinical point of view, however, this is a very inadequate description, as there are significant comorbidities associated with chronic tinnitus. Neurophysiological investigations with different imaging techniques give a very similar picture, because not only the auditory system is affected in chronic tinnitus patients, but also a widely ramified subcortical and cortical network. In addition to auditory processing systems, networks consisting of frontal and parietal regions are particularly disturbed. For this reason, some authors conceptualize tinnitus as a network disorder rather than a disorder of a circumscribed system. These findings and this concept suggest that tinnitus must be diagnosed and treated in a multidisciplinary and multimodal manner.

Thalamo-cortical circuits during sensory attenuation in emerging psychosis: a combined magnetoencephalography and dynamic causal modelling study.

Evidence suggests that schizophrenia (ScZ) involves impairments in sensory attenuation. It is currently unclear, however, whether such deficits are present during early-stage psychosis as well as the underlying network and the potential as a biomarker. To address these questions, Magnetoencephalography (MEG) was used in combination with computational modeling to examine M100 responses that involved a "passive" condition during which tones were binaurally presented, while in an "active" condition participants were asked to generate a tone via a button press. MEG data were obtained from 109 clinical high-risk for psychosis (CHR-P) participants, 23 people with a first-episode psychosis (FEP), and 48 healthy controls (HC). M100 responses at sensor and source level in the left and right thalamus (THA), Heschl's gyrus (HES), superior temporal gyrus (STG) and right inferior parietal cortex (IPL) were examined and dynamic causal modeling (DCM) was performed. Furthermore, the relationship between sensory attenuation and persistence of attenuated psychotic symptoms (APS) and transition to psychosis was investigated in CHR-P participants. Sensory attenuation was impaired in left HES, left STG and left THA in FEP patients, while in the CHR-P group deficits were observed only in right HES. DCM results revealed that CHR-P participants showed reduced top-down modulation from the right IPL to the right HES. Importantly, deficits in sensory attenuation did not predict clinical outcomes in the CHR-P group. Our results show that early-stage psychosis involves impaired sensory attenuation in auditory and thalamic regions but may not predict clinical outcomes in CHR-P participants.

Physically based equation of state for Mie ν-6 fluids.

We develop a physically based equation of state that describes Mie ν-6 fluids with an accuracy comparable to that of state-of-the-art empirical models. The equation of state is developed within the framework of the uv-theory [T. van Westen and J. Gross, J. Chem. Phys. 155, 244501 (2021)], which is modified by incorporating the third virial coefficient B3 in the low-density description of the model. The new model interpolates between a first-order Weeks-Chandler-Andersen (WCA) perturbation theory at high densities and a modified first-order WCA theory that recovers the virial expansion up to B3 at low densities. A new algebraic equation for the third virial coefficient of Mie ν-6 fluids is developed-other inputs are taken from previous work. Predicted thermodynamic properties and phase equilibria are compared to a comprehensive database of molecular simulation results from the literature, including Mie fluids of repulsive exponents 9 ≤ ν ≤ 48. The new equation of state is applicable to states with densities up to ρ*(T*)⪅1.1+0.12T* and temperatures T* > 0.3. For the Lennard-Jones fluid (ν = 12), the performance of the model is comparable to that of the best empirical equations of state available. As compared to empirical models, the physical basis of the new model provides several advantages, however: (1) the new model is applicable to Mie fluids of repulsive exponents 9 ≤ ν ≤ 48 instead of only ν = 12, (2) the model leads to a better description of the meta-stable and unstable region (which is important for describing interfacial properties by classical density functional theory), and (3) being a first-order perturbation theory, the new model (potentially) allows an easier and more rigorous extension to non-spherical (chain) fluids and mixtures.