Green-Kubo expressions for transport coefficients from dissipative particle dynamics simulations revisited.

This article addresses the debate about the correct application of Green-Kubo expressions for transport coefficients from dissipative particle dynamics simulations. We demonstrate that the Green-Kubo expressions are valid provided that (i) the dynamic model conserves the physical property, whose transport is studied, and (ii) the fluctuations satisfy detailed balance. As a result, the traditional expressions used in molecular dynamics can also be applied to dissipative particle dynamics simulations. However, taking the calculation of the shear viscosity as a paradigmatic example, a random contribution, whose strength scales as 1/δt1/2, with δt the time-step, can cause difficulties if the stress tensor is not separated into the different contributions. We compare our expression to that of Ernst and Brito (M. H. Ernst and R. Brito, Europhys. Lett., 2006, 73, 183-189), which arises from a diametrically different perspective. We demonstrate that the two expressions are completely equivalent and find exactly the same result both analytically and numerically. We show that the differences are not due to the lack of time-reversibility but instead from a pre-averaging of the random contributions. Despite the overall validity of Green-Kubo expressions, we find that the Einstein-Helfand relations (D. C. Malaspina et al. Phys. Chem. Chem. Phys., 2023, 25, 12025-12040) do not suffer from the need to decompose the stress tensor and can readily be used with a high degree of accuracy. Consequently, Einstein-Helfand relations should be seen as the preferred method to calculate transport coefficients from dissipative particle dynamics simulations.

Transport coefficients from Einstein-Helfand relations using standard and energy-conserving dissipative particle dynamics methods.

In this article we demonstrate that contrary to general belief, the standard Einstein-Helfand (EH) formulas are valid for the evaluation of transport coefficients of systems containing dissipative and random forces provided that for these mesoscopic systems: (i) the corresponding conservation laws are satisfied, and (ii) the transition probabilities satisfy detailed balance. Dissipative particle dynamics (DPD) and energy-conserving DPD methods (DPDE), for instance, are archetypical of such mesoscopic approaches satisfying these properties. To verify this statement, we have derived a mesoscopic heat flux form for the DPDE method, suitable for the calculation of the thermal conductivity from an EH expression. We have compared EH measurements against non-equilibrium simulation values for different scenarios, including many-body potentials, and have found excellent agreement in all cases. The expressions are valid notably for systems with density- and temperature-dependent potentials, such as the recently developed generalised DPDE method (GenDPDE) [Avalos et al., Phys. Chem. Chem. Phys., 2019, 21, 24891]. We thus demonstrate that traditional EH formulas in equilibrium simulations can be widely used to obtain transport coefficients, provided that the appropriate expression for the associated flux is used.

Prompt Response of the Dayside Magnetosphere to Discrete Structures Within the Sheath Region of a Coronal Mass Ejection.

A sequence of discrete solar wind structures within the sheath region of an interplanetary coronal mass ejection on November 6, 2015, caused a series of compressions and releases of the dayside magnetosphere. Each compression resulted in a brief adiabatic enhancement of ions (electrons) driving bursts of electromagnetic ion cyclotron (EMIC; whistler mode chorus) wave growth across the dayside magnetosphere. Fine-structured rising tones were observed in the EMIC wave bursts, resulting in nonlinear scattering of relativistic electrons in the outer radiation belt. Multipoint observations allow us to study the spatial structure and evolution of these sheath structures as they propagate Earthward from L1 as well as the spatio-temporal characteristics of the magnetospheric response. This event highlights the importance of fine-scale solar wind structure, in particular within complex sheath regions, in driving dayside phenomena within the inner magnetosphere.

Electromagnetic power of lightning superbolts from Earth to space.

Lightning superbolts are the most powerful and rare lightning events with intense optical emission, first identified from space. Superbolt events occurred in 2010-2018 could be localized by extracting the high energy tail of the lightning stroke signals measured by the very low frequency ground stations of the World-Wide Lightning Location Network. Here, we report electromagnetic observations of superbolts from space using Van Allen Probes satellite measurements, and ground measurements, and with two events measured both from ground and space. From burst-triggered measurements, we compute electric and magnetic power spectral density for very low frequency waves driven by superbolts, both on Earth and transmitted into space, demonstrating that superbolts transmit 10-1000 times more powerful very low frequency waves into space than typical strokes and revealing that their extreme nature is observed in space. We find several properties of superbolts that notably differ from most lightning flashes; a more symmetric first ground-wave peak due to a longer rise time, larger peak current, weaker decay of electromagnetic power density in space with distance, and a power mostly confined in the very low frequency range. Their signal is absent in space during day times and is received with a long-time delay on the Van Allen Probes. These results have implications for our understanding of lightning and superbolts, for ionosphere-magnetosphere wave transmission, wave propagation in space, and remote sensing of extreme events.

Parker Solar Probe Enters the Magnetically Dominated Solar Corona.

The high temperatures and strong magnetic fields of the solar corona form streams of solar wind that expand through the Solar System into interstellar space. At 09:33 UT on 28 April 2021 Parker Solar Probe entered the magnetized atmosphere of the Sun 13 million km above the photosphere, crossing below the Alfvén critical surface for five hours into plasma in casual contact with the Sun with an Alfvén Mach number of 0.79 and magnetic pressure dominating both ion and electron pressure. The spectrum of turbulence below the Alfvén critical surface is reported. Magnetic mapping suggests the region was a steady flow emerging on rapidly expanding coronal magnetic field lines lying above a pseudostreamer. The sub-Alfvénic nature of the flow may be due to suppressed magnetic reconnection at the base of the pseudostreamer, as evidenced by unusually low densities in this region and the magnetic mapping.

Preliminary Findings Associate Hippocampal 1H-MR Spectroscopic Metabolite Concentrations with Psychotic and Manic Symptoms in Patients with Schizophrenia.

BACKGROUND AND PURPOSE: Previous hippocampal proton MR spectroscopic imaging distinguished patients with schizophrenia from controls by elevated Cr levels and significantly more variable NAA and Cho concentrations. This goal of this study was to ascertain whether this metabolic variability is associated with clinical features of the syndrome, possibly reflecting heterogeneous hippocampal pathologies and perhaps variability in its "positive" (psychotic) and "negative" (social and emotional deficits) symptoms. MATERIALS AND METHODS: In a sample of 15 patients with schizophrenia according to the Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition, we examined the association of NAA and Cho levels with research diagnostic interviews and clinical symptom ratings of the patients. Metabolite concentrations were previously obtained with 3D proton MR spectroscopic imaging at 3T, a technique that facilitates complete coverage of this small, irregularly shaped, bilateral, temporal lobe structure. RESULTS: The patient cohort comprised 8 men and 7 women (mean age, 39.1 [SD, 10.8] years, with a mean disease duration of 17.2 [SD, 10.8] years. Despite the relatively modest cohort size, we found the following: 1) Elevated Cho levels predict the positive (psychotic, r = 0.590, P = .021) and manic (r = 0.686, P = .005) symptom severity; and 2) lower NAA levels trend toward negative symptoms (r = 0.484, P = .08). No clinical symptoms were associated with Cr level or hippocampal volume (all, P ≥ .055). CONCLUSIONS: These preliminary findings suggest that NAA and Cho variations reflect different pathophysiologic processes, consistent with microgliosis/astrogliosis and/or lower vitality (reduced NAA) and demyelination (elevated Cho). In particular, the active state-related symptoms, including psychosis and mania, were associated with demyelination. Consequently, their deviations from the means of healthy controls may be a marker that may benefit precision medicine in selection and monitoring of schizophrenia treatment.

Plasma Double Layers at the Boundary Between Venus and the Solar Wind.

The solar wind is slowed, deflected, and heated as it encounters Venus's induced magnetosphere. The importance of kinetic plasma processes to these interactions has not been examined in detail, due to a lack of constraining observations. In this study, kinetic-scale electric field structures are identified in the Venusian magnetosheath, including plasma double layers. The double layers may be driven by currents or mixing of inhomogeneous plasmas near the edge of the magnetosheath. Estimated double-layer spatial scales are consistent with those reported at Earth. Estimated potential drops are similar to electron temperature gradients across the bow shock. Many double layers are found in few high cadence data captures, suggesting that their amplitudes are high relative to other magnetosheath plasma waves. These are the first direct observations of plasma double layers beyond near-Earth space, supporting the idea that kinetic plasma processes are active in many space plasma environments.

Highly structured slow solar wind emerging from an equatorial coronal hole.

During the solar minimum, when the Sun is at its least active, the solar wind1,2 is observed at high latitudes as a predominantly fast (more than 500 kilometres per second), highly Alfvénic rarefied stream of plasma originating from deep within coronal holes. Closer to the ecliptic plane, the solar wind is interspersed with a more variable slow wind3 of less than 500 kilometres per second. The precise origins of the slow wind streams are less certain4; theories and observations suggest that they may originate at the tips of helmet streamers5,6, from interchange reconnection near coronal hole boundaries7,8, or within coronal holes with highly diverging magnetic fields9,10. The heating mechanism required to drive the solar wind is also unresolved, although candidate mechanisms include Alfvén-wave turbulence11,12, heating by reconnection in nanoflares13, ion cyclotron wave heating14 and acceleration by thermal gradients1. At a distance of one astronomical unit, the wind is mixed and evolved, and therefore much of the diagnostic structure of these sources and processes has been lost. Here we present observations from the Parker Solar Probe15 at 36 to 54 solar radii that show evidence of slow Alfvénic solar wind emerging from a small equatorial coronal hole. The measured magnetic field exhibits patches of large, intermittent reversals that are associated with jets of plasma and enhanced Poynting flux and that are interspersed in a smoother and less turbulent flow with a near-radial magnetic field. Furthermore, plasma-wave measurements suggest the existence of electron and ion velocity-space micro-instabilities10,16 that are associated with plasma heating and thermalization processes. Our measurements suggest that there is an impulsive mechanism associated with solar-wind energization and that micro-instabilities play a part in heating, and we provide evidence that low-latitude coronal holes are a key source of the slow solar wind.

Electron Bulk Acceleration and Thermalization at Earth's Quasiperpendicular Bow Shock.

Electron heating at Earth's quasiperpendicular bow shock has been surmised to be due to the combined effects of a quasistatic electric potential and scattering through wave-particle interaction. Here we report the observation of electron distribution functions indicating a new electron heating process occurring at the leading edge of the shock front. Incident solar wind electrons are accelerated parallel to the magnetic field toward downstream, reaching an electron-ion relative drift speed exceeding the electron thermal speed. The bulk acceleration is associated with an electric field pulse embedded in a whistler-mode wave. The high electron-ion relative drift is relaxed primarily through a nonlinear current-driven instability. The relaxed distributions contain a beam traveling toward the shock as a remnant of the accelerated electrons. Similar distribution functions prevail throughout the shock transition layer, suggesting that the observed acceleration and thermalization is essential to the cross-shock electron heating.

Fast Diffusion of Ultrarelativistic Electrons in the Outer Radiation Belt: 17 March 2015 Storm Event.

Inward radial diffusion driven by ULF waves has long been known to be capable of accelerating radiation belt electrons to very high energies within the heart of the belts, but more recent work has shown that radial diffusion values can be highly event-specific, and mean values or empirical models may not capture the full significance of radial diffusion to acceleration events. Here we present an event of fast inward radial diffusion, occurring during a period following the geomagnetic storm of 17 March 2015. Ultrarelativistic electrons up to ∼8 MeV are accelerated in the absence of intense higher-frequency plasma waves, indicating an acceleration event in the core of the outer belt driven primarily or entirely by ULF wave-driven diffusion. We examine this fast diffusion rate along with derived radial diffusion coefficients using particle and fields instruments on the Van Allen Probes spacecraft mission.

Chemistry-based molecular signature underlying the atypia of clozapine.

The central nervous system is functionally organized as a dynamic network of interacting neural circuits that underlies observable behaviors. At higher resolution, these behaviors, or phenotypes, are defined by the activity of a specific set of biomolecules within those circuits. Identification of molecules that govern psychiatric phenotypes is a major challenge. The only organic molecular entities objectively associated with psychiatric phenotypes in humans are drugs that induce psychiatric phenotypes and drugs used for treatment of specific psychiatric conditions. Here, we identified candidate biomolecules contributing to the organic basis for psychosis by deriving an in vivo biomolecule-tissue signature for the atypical pharmacologic action of the antipsychotic drug clozapine. Our novel in silico approach identifies the ensemble of potential drug targets based on the drug's chemical structure and the region-specific gene expression profile of each target in the central nervous system. We subtracted the signature of the action of clozapine from that of a typical antipsychotic, chlorpromazine. Our results implicate dopamine D4 receptors in the pineal gland and muscarinic acetylcholine M1 (CHRM1) and M3 (CHRM3) receptors in the prefrontal cortex (PFC) as significant and unique to clozapine, whereas serotonin receptors 5-HT2A in the PFC and 5-HT2C in the caudate nucleus were common significant sites of action for both drugs. Our results suggest that D4 and CHRM1 receptor activity in specific tissues may represent underappreciated drug targets to advance the pharmacologic treatment of schizophrenia. These findings may enhance our understanding of the organic basis of psychiatric disorders and help developing effective therapies.

Evidence for genetic heterogeneity between clinical subtypes of bipolar disorder.

We performed a genome-wide association study of 6447 bipolar disorder (BD) cases and 12 639 controls from the International Cohort Collection for Bipolar Disorder (ICCBD). Meta-analysis was performed with prior results from the Psychiatric Genomics Consortium Bipolar Disorder Working Group for a combined sample of 13 902 cases and 19 279 controls. We identified eight genome-wide significant, associated regions, including a novel associated region on chromosome 10 (rs10884920; P=3.28 × 10-8) that includes the brain-enriched cytoskeleton protein adducin 3 (ADD3), a non-coding RNA, and a neuropeptide-specific aminopeptidase P (XPNPEP1). Our large sample size allowed us to test the heritability and genetic correlation of BD subtypes and investigate their genetic overlap with schizophrenia and major depressive disorder. We found a significant difference in heritability of the two most common forms of BD (BD I SNP-h2=0.35; BD II SNP-h2=0.25; P=0.02). The genetic correlation between BD I and BD II was 0.78, whereas the genetic correlation was 0.97 when BD cohorts containing both types were compared. In addition, we demonstrated a significantly greater load of polygenic risk alleles for schizophrenia and BD in patients with BD I compared with patients with BD II, and a greater load of schizophrenia risk alleles in patients with the bipolar type of schizoaffective disorder compared with patients with either BD I or BD II. These results point to a partial difference in the genetic architecture of BD subtypes as currently defined.

Follow-up with primary care providers for elevated glycated haemoglobin identified at the dental visit.

OBJECTIVES: This study examined patient experiences after receiving elevated diabetes screening values using blood collected at a dental clinic. It explores patients' reactions to screening, whether or not they sought recommended medical follow-up, and facilitating factors and barriers to obtaining follow-up care. METHODS: At the comprehensive care clinics at a large, urban College of Dentistry in the United States, haemoglobin A1C (HbA1C) values were obtained from 379 study participants who had not been previously diagnosed with diabetes. In all, 169 (44.6%) had elevated HbA1C values. We analysed quantitative and qualitative data concerning these patients' follow-up with primary care providers (PCPs). RESULTS: We were able to contact 112 (66.3%) of the 169 study participants who had an elevated HbA1C reading. Of that group, 61 (54.5%) received recommended follow-up care from a PCP within 3 months, and an additional 28 (25.0%) said they intended to seek such care. Qualitative themes included the following: the screening letter - opportunity or burden, appreciation for the 3-month follow-up call and barriers to medical follow-up that included the following: lack of knowledge about diabetes, not understanding the importance of follow-up, busyness, financial concerns, fear and denial. CONCLUSIONS: Quantitative and qualitative data demonstrate that dentists, dental hygienists and nurses are well poised to discover and translate new models of patient-centred, comprehensive care to patients with oral and systemic illness.

Live birth sex ratios and father's geographic origins in Jerusalem, 1964-1976.

OBJECTIVE: To examine whether ancestry influenced sex ratios of offspring in a birth cohort before parental antenatal sex selection influenced offspring sex. METHODS: We measured the sex ratio as the percent of males according to countries of birth of paternal and maternal grandfathers in 91,459 live births from 1964 to 1976 in the Jerusalem Perinatal Study. Confidence limits (CI) were computed based on an expected sex ratio of 1.05, which is 51.4% male. RESULTS: Of all live births recorded, 51.4% were male. Relative to Jewish ancestry (51.4% males), significantly more males (1,761) were born to Muslim ancestry (54.5, 95% CI = 52.1-56.8, P = 0.01). Among the former, sex ratios were not significantly associated with paternal or maternal age, education, or offspring's birth order. Consistent with a preference for male offspring, the sex ratio decreased despite increasing numbers of births over the 13-year period. Sex ratios were not affected by maternal or paternal origins in North Africa or Europe. However, the offspring whose paternal grandfathers were born in Western Asia included fewer males than expected (50.7, 50.1-51.3, P = 0.02), whether the father was born abroad (50.7) or in Israel (50.8). This was observed for descendents of paternal grandfathers born in Lebanon (47.6), Turkey (49.9), Yemen & Aden (50.2), Iraq (50.5), Afghanistan (50.5), Syria (50.6), and Cyprus (50.7); but not for those from India (51.5) or Iran (51.9). The West Asian group showed the strongest decline in sex ratios with increasing paternal family size. CONCLUSIONS: A decreased sex ratio associated with ancestry in Western Asia is consistent with reduced ability to bear sons by a subset of Jewish men in the Jerusalem cohort. Lower sex ratios may be because of pregnancy stress, which may be higher in this subgroup. Alternatively, a degrading Y chromosome haplogroup or other genetic or epigenetic differences on male germ lines could affect birth ratios, such as differential exposure to an environmental agent, dietary differences, or stress. Differential stopping behaviors that favor additional pregnancies following the birth of a daughter might exacerbate these lower sex ratios.

Metabolic Abnormalities in the Hippocampus of Patients with Schizophrenia: A 3D Multivoxel MR Spectroscopic Imaging Study at 3T.

BACKGROUND AND PURPOSE: Schizophrenia is well-known to be associated with hippocampal structural abnormalities. We used 1H-MR spectroscopy to test the hypothesis that these abnormalities are accompanied by NAA deficits, reflecting neuronal dysfunction, in patients compared with healthy controls. MATERIALS AND METHODS: Nineteen patients with schizophrenia (11 men; mean age, 40.6 ± 10.1 years; mean disease duration, 19.5 ± 10.5 years) and 11 matched healthy controls (5 men; mean age, 33.7 ± 10.1 years) underwent MR imaging and multivoxel point-resolved spectroscopy (TE/TR, 35/1400 ms) 1H-MRS at 3T to obtain their hippocampal GM absolute NAA, Cr, Cho, and mIns concentrations. Unequal variance t tests and ANCOVA were used to compare patients with controls. Bilateral volumes from manually outlined hippocampal masks were compared by using unequal variance t tests. RESULTS: Patients' average hippocampal GM Cr concentrations were 19% higher than that of controls, 8.7 ± 2.2 versus 7.4 ± 1.2 mmol/L (P < .05); showing no differences, concentrations in NAA were 8.8 ± 1.6 versus 8.7 ± 1.2 mmol/L; in Cho, 2.3 ± 0.7 versus 2.1 ± 0.3 mmol/L; and in mIns, 6.1 ± 1.5 versus 5.2 ± 0.9 (all P > .1). There was a positive correlation between mIns and Cr in patients (r = 0.57, P = .05) but not in controls. The mean bilateral hippocampal volume was ∼10% lower in patients: 7.5 ± 0.9 versus 8.4 ± 0.7 cm3 (P < .05). CONCLUSIONS: These findings suggest that the hippocampal volume deficit in schizophrenia is not due to net loss of neurons, in agreement with histopathology studies but not with prior 1H-MR spectroscopy reports. Elevated Cr is consistent with hippocampal hypermetabolism, and its correlation with mIns may also suggest an inflammatory process affecting some cases; these findings may suggest treatment targets and markers to monitor them.

Magnetospheric Multiscale Observations of the Electron Diffusion Region of Large Guide Field Magnetic Reconnection.

We report observations from the Magnetospheric Multiscale (MMS) satellites of a large guide field magnetic reconnection event. The observations suggest that two of the four MMS spacecraft sampled the electron diffusion region, whereas the other two spacecraft detected the exhaust jet from the event. The guide magnetic field amplitude is approximately 4 times that of the reconnecting field. The event is accompanied by a significant parallel electric field (E_{∥}) that is larger than predicted by simulations. The high-speed (∼300 km/s) crossing of the electron diffusion region limited the data set to one complete electron distribution inside of the electron diffusion region, which shows significant parallel heating. The data suggest that E_{∥} is balanced by a combination of electron inertia and a parallel gradient of the gyrotropic electron pressure.

Magnetospheric Multiscale Satellites Observations of Parallel Electric Fields Associated with Magnetic Reconnection.

We report observations from the Magnetospheric Multiscale satellites of parallel electric fields (E_{∥}) associated with magnetic reconnection in the subsolar region of the Earth's magnetopause. E_{∥} events near the electron diffusion region have amplitudes on the order of 100 mV/m, which are significantly larger than those predicted for an antiparallel reconnection electric field. This Letter addresses specific types of E_{∥} events, which appear as large-amplitude, near unipolar spikes that are associated with tangled, reconnected magnetic fields. These E_{∥} events are primarily in or near a current layer near the separatrix and are interpreted to be double layers that may be responsible for secondary reconnection in tangled magnetic fields or flux ropes. These results are telling of the three-dimensional nature of magnetopause reconnection and indicate that magnetopause reconnection may be often patchy and/or drive turbulence along the separatrix that results in flux ropes and/or tangled magnetic fields.

Dust observations at orbital altitudes surrounding Mars.

Dust is common close to the martian surface, but no known process can lift appreciable concentrations of particles to altitudes above ~150 kilometers. We present observations of dust at altitudes ranging from 150 to above 1000 kilometers by the Langmuir Probe and Wave instrument on the Mars Atmosphere and Volatile Evolution spacecraft. Based on its distribution, we interpret this dust to be interplanetary in origin. A comparison with laboratory measurements indicates that the dust grain size ranges from 1 to 12 micrometers, assuming a typical grain velocity of ~18 kilometers per second. These direct observations of dust entering the martian atmosphere improve our understanding of the sources, sinks, and transport of interplanetary dust throughout the inner solar system and the associated impacts on Mars's atmosphere.

Global-scale coherence modulation of radiation-belt electron loss from plasmaspheric hiss.

Over 40 years ago it was suggested that electron loss in the region of the radiation belts that overlaps with the region of high plasma density called the plasmasphere, within four to five Earth radii, arises largely from interaction with an electromagnetic plasma wave called plasmaspheric hiss. This interaction strongly influences the evolution of the radiation belts during a geomagnetic storm, and over the course of many hours to days helps to return the radiation-belt structure to its 'quiet' pre-storm configuration. Observations have shown that the long-term electron-loss rate is consistent with this theory but the temporal and spatial dynamics of the loss process remain to be directly verified. Here we report simultaneous measurements of structured radiation-belt electron losses and the hiss phenomenon that causes the losses. Losses were observed in the form of bremsstrahlung X-rays generated by hiss-scattered electrons colliding with the Earth's atmosphere after removal from the radiation belts. Our results show that changes of up to an order of magnitude in the dynamics of electron loss arising from hiss occur on timescales as short as one to twenty minutes, in association with modulations in plasma density and magnetic field. Furthermore, these loss dynamics are coherent with hiss dynamics on spatial scales comparable to the size of the plasmasphere. This nearly global-scale coherence was not predicted and may affect the short-term evolution of the radiation belts during active times.

Age-related sperm DNA methylation changes are transmitted to offspring and associated with abnormal behavior and dysregulated gene expression.

Advanced paternal age (APA) has been shown to be a significant risk factor in the offspring for neurodevelopmental psychiatric disorders, such as schizophrenia and autism spectrum disorders. During aging, de novo mutations accumulate in the male germline and are frequently transmitted to the offspring with deleterious effects. In addition, DNA methylation during spermatogenesis is an active process, which is susceptible to errors that can be propagated to subsequent generations. Here we test the hypothesis that the integrity of germline DNA methylation is compromised during the aging process. A genome-wide DNA methylation screen comparing sperm from young and old mice revealed a significant loss of methylation in the older mice in regions associated with transcriptional regulation. The offspring of older fathers had reduced exploratory and startle behaviors and exhibited similar brain DNA methylation abnormalities as observed in the paternal sperm. Offspring from old fathers also had transcriptional dysregulation of developmental genes implicated in autism and schizophrenia. Our findings demonstrate that DNA methylation abnormalities arising in the sperm of old fathers are a plausible mechanism to explain some of the risks that APA poses to resulting offspring.