Search

We report observations of direct evidence of energetic protons being accelerated above ~400 keV within the reconnection exhaust of a heliospheric current sheet (HCS) crossing by NASA's Parker Solar Probe (PSP) at a distance of ~16.25 solar radii (Rs) from the Sun. Inside the extended exhaust, both the reconnection-generated plasma jets and the accelerated protons propagated toward the Sun, unambiguously establishing their origin from HCS reconnection sites located beyond PSP. Within the core of the exhaust, PSP detected stably trapped energetic protons up to ~400 keV, which is approximately 1000 times greater than the available magnetic energy per particle. The differential energy spectrum of the accelerated protons behaved as a pure power-law with spectral index of about -5. Supporting simulations using the kglobal model suggest that the trapping and acceleration of protons up to ~400 keV in the reconnection exhaust is likely facilitated by merging magnetic islands with a guide field between ~0.2-0.3 of the reconnecting magnetic field, consistent with the observations. These new results, enabled by PSP's proximity to the Sun, demonstrate that magnetic reconnection in the HCS is a significant new source of energetic particles in the near-Sun solar wind. The discovery of in-situ particle acceleration via magnetic reconnection at the HCS provides valuable insights into this fundamental process which frequently converts the large magnetic field energy density in the near-Sun plasma environment and may be responsible for heating the sun's atmosphere, accelerating the solar wind, and energizing charged particles to extremely high energies in solar flares., Comment: 24 pages, 4 figures

We present an event observed by Parker Solar Probe at $\sim$0.2 au on March 2, 2022 in which imaging and \emph{in situ} measurements coincide. During this event, PSP passed through structures on the flank of a streamer blowout CME including an isolated flux tube in front of the CME, a turbulent sheath, and the CME itself. Imaging observations and \emph{in situ} helicity and principal variance signatures consistently show the presence of flux ropes internal to the CME. In both the sheath, and the CME interval, the distributions are more isotropic, the spectra are softer, and the abundance ratios of Fe/O and He/H are lower than those in the isolated flux tube, and yet elevated relative to typical plasma and SEP abundances. These signatures in the sheath and the CME indicate that both flare populations and those from the plasma are accelerated to form the observed energetic particle enhancements. In contrast, the isolated flux tube shows large streaming, hard spectra and large Fe/O and He/H ratios, indicating flare sources. Energetic particle fluxes are most enhanced within the CME interval from suprathermal through energetic particle energies ($\sim$ keV to $>10$ MeV), indicating particle acceleration, and confinement local to the closed magnetic structure. The flux-rope morphology of the CME helps to enable local modulation and trapping of energetic particles, particularly along helicity channels and other plasma boundaries. Thus, the CME acts to build-up energetic particle populations, allowing them to be fed into subsequent higher energy particle acceleration throughout the inner heliosphere where a compression or shock forms on the CME front., Comment: 41 pages, 19 figures, In Press

A central problem of space plasma physics is how protons and electrons are heated in a turbulent, magnetized plasma. The differential heating of charged species due to dissipation of turbulent fluctuations plays a key role in solar wind evolution. Measurements from previous heliophysics missions have provided estimates of proton and electron heating rates beyond 0.27 au. Using Parker Solar Probe (PSP) data accumulated during the first ten encounters, we extend the evaluation of the individual rates of heat deposition for protons and electrons in to a distance of 0.063 au (13.5 solar radii), in the newly formed solar wind. The PSP data in the near-Sun environment show different behavior of the electron heat conduction flux from what was predicted from previous fits to Helios and Ulysses data. Consequently, the empirically derived proton and electron heating rates exhibit significantly different behavior than previous reports, with the proton heating becoming increasingly dominant over electron heating at decreasing heliocentric distances. We find that the protons receive about 80% of the total plasma heating at ~ 13 solar radii, slightly higher than the near-Earth values. This empirically derived heating partition between protons and electrons will help to constrain theoretical models of solar wind heating., Comment: Accepted for publication in Astrophysical Journal Letters

A variety of energy sources, ranging from dynamic processes like magnetic reconnection and waves to quasi-steady terms like the plasma pressure, may contribute to the acceleration of the solar wind. We utilize a combination of charged particle and magnetic field observations from the Parker Solar Probe (PSP) to attempt to quantify the steady-state contribution of the proton pressure, the electric potential, and the wave energy to the solar wind proton acceleration observed by PSP between 13.3 and ~100 solar radii (RS). The proton pressure provides a natural kinematic driver of the outflow. The ambipolar electric potential acts to couple the electron pressure to the protons, providing another definite proton acceleration term. Fluctuations and waves, while inherently dynamic, can act as an additional effective steady-state pressure term. To analyze the contributions of these terms, we utilize radial binning of single-point PSP measurements, as well as repeated crossings of the same stream at different distances on individual PSP orbits (i.e. "fast radial scans"). In agreement with previous work, we find that the electric potential contains sufficient energy to fully explain the acceleration of the slower wind streams. On the other hand, we find that the wave pressure plays an increasingly important role in the faster wind streams. The combination of these terms can explain the continuing acceleration of both slow and fast wind streams beyond 13.3 RS.

This letter describes a new phenomenon on the Parker Solar Probe of recurring plasma density enhancements that have $\Delta$n/n ~10% and that occur at a repetition rate of ~5 Hz. They were observed sporadically for about five hours between 14 and 15 solar radii on Parker Solar Probe orbit 12 and they were also seen in the same radial range on both the inbound and outbound orbits 11. Their apparently steady-state existence suggests that their pressure gradient was balanced by the electric field. The EX electric field component produced from this requirement is in good agreement with that measured. This provides strong evidence for the measurement accuracy of the density fluctuations and the X- and Y-components of the electric field (the Z-component was not measured). The electrostatic density waves were accompanied by an electromagnetic low frequency wave which occurred with the electrostatic harmonics. The amplitudes of these electrostatic and electromagnetic waves at $\ge$ 1 Hz were greater than the amplitude of the Alfvenic turbulence in their vicinity so they can be important for the heating, scattering, and acceleration of the plasma. The existence of this pair of waves is consistent with the observed plasma distributions and is explained by a magneto-acoustic wave theory that produces a low frequency electromagnetic wave and electrostatic harmonics., Comment: 9 pages including 5 figures

We utilize observations from the Parker Solar Probe (PSP) to study the radial evolution of the solar wind in the inner heliosphere. We analyze electron velocity distribution functions observed by the Solar Wind Electrons, Alphas, and Protons suite to estimate the coronal electron temperature and the local electric potential in the solar wind. From the latter value and the local flow speed, we compute the asymptotic solar wind speed. We group the PSP observations by asymptotic speed, and characterize the radial evolution of the wind speed, electron temperature, and electric potential within each group. In agreement with previous work, we find that the electron temperature (both local and coronal) and the electric potential are anti-correlated with wind speed. This implies that the electron thermal pressure and the associated electric field can provide more net acceleration in the slow wind than in the fast wind. We then utilize the inferred coronal temperature and the extrapolated electric + gravitational potential to show that both electric field driven exospheric models and the equivalent thermally driven hydrodynamic models can explain the entire observed speed of the slowest solar wind streams. On the other hand, neither class of model can explain the observed speed of the faster solar wind streams, which thus require additional acceleration mechanisms., Comment: Submitted to the Astrophysical Journal

Magnetic switchbacks, or sudden reversals in the magnetic field's radial direction, are one of the more striking observations of Parker Solar Probe (PSP) thus far in its mission. While their precise production mechanisms are still unknown, the two main theories are via interchange reconnection events and in-situ generation. In this work density and abundance variations of alpha particles are studied inside and outside individual switchbacks. We find no consistent compositional differences in the alpha particle abundance ratio, $n_{\alpha p}$, inside vs outside, nor do we observe any signature when separating the switchbacks according to $V_{\alpha p}/V_{pw}$, the ratio of alpha-proton differential speed to the wave phase speed (speed the switchback is travelling). We argue these measurements cannot be used to rule in favour of one production mechanism over the other, due to the distance between PSP and the postulated interchange reconnection events. In addition we examine the 3D velocity fluctuations of protons and alpha particles within individual switchbacks. While switchbacks are always associated with increases in proton velocity, alpha velocities may be enhanced, unchanged, or decrease. This is due to the interplay between $V_{pw}$ and $V_{\alpha p}$, with the Alfv\'enic motion of the alpha particles vanishing as the difference $|V_{pw} - V_{\alpha p}|$ decreases. We show how the Alfv\'enic motion of both the alphas and the protons through switchbacks can be understood as approximately rigid arm rotation about the location of the wave frame, and illustrate that the wave frame can therefore be estimated using particle measurements alone, via sphere fitting., Comment: To be published in ApJ

In the lower solar coronal regions where the magnetic field is dominant, the Alfven speed is much higher than the wind speed. In contrast, the near-Earth solar wind is strongly super-Alfvenic, i.e., the wind speed greatly exceeds the Alfven speed. The transition between these regimes is classically described as the "Alfven point" but may in fact occur in a distributed Alfven critical region. NASA's Parker Solar Probe (PSP) mission has entered this region, as it follows a series of orbits that gradually approach more closely to the sun. During its 8th and 9th solar encounters, at a distance of 16 solar radii from the Sun, PSP sampled four extended periods in which the solar wind speed was measured to be smaller than the local Alfven speed. These are the first in-situ detections of sub-Alfvenic solar wind in the inner heliosphere by PSP. Here we explore properties of these samples of sub-Alfvenic solar wind, which may provide important previews of the physical processes operating at lower altitude. Specifically, we characterize the turbulence, anisotropy, intermittency, and directional switchback properties of these sub-Alfvenic winds and contrast these with the neighboring super-Alfvenic periods., Comment: Accepted for publication in the Astrophysical Journal Letter

We summarize a series of numerical experiments of collisional dynamics in dense stellar systems such as globular clusters (GCs) and in weakly collisional plasmas using a novel simulation technique, the so-called Multi-particle collision (MPC) method, alternative to Fokker-Planck and Monte Carlo approaches. MPC is related to particle-mesh approaches for the computation of self consistent long-range fields, ensuring that simulation time scales with $N\log N$ in the number of particles, as opposed to $N^2$ for direct $N$-body. The collisional relaxation effects are modelled by computing particle interactions based on a collision operator approach that ensures rigorous conservation of energy and momenta and depends only on particles velocities and cell-based integrated quantities., Comment: 7 pages, 2 figures. To appear in the proceedings of the 362 IAU symposium "Predictive Power of Computational Astrophysics as a Discovery Tool". Chamonix, France 8-12 Nov. 2021

We present observations of >10-100 keV/nucleon suprathermal (ST) H, He, O, and Fe ions associated with crossings of the heliospheric current sheet (HCS) at radial distances <0.1 au from the Sun. Our key findings are: 1) very few heavy ions are detected during the 1st full crossing, the heavy ion intensities are reduced during the 2nd partial crossing and peak just after the 2nd crossing; 2) ion arrival times exhibit no velocity dispersion; 3) He pitch-angle distributions track the magnetic field polarity reversal and show up to ~10:1 anti-sunward, field-aligned flows and beams closer to the HCS that become nearly isotropic further from the HCS; 4) the He spectrum steepens either side of the HCS and the He, O, and Fe spectra exhibit power-laws of the form ~E^4-6; and 5) maximum energies EX increase with the ion's charge-to-mass (Q/M) ratio as EX/EH proportional to [(QX/MX)]^alpha where alpha~0.65-0.76, assuming that the average Q-states are similar to those measured in gradual and impulsive solar energetic particle events at 1 au. The absence of velocity dispersion in combination with strong field-aligned anisotropies closer to the HCS appears to rule out solar flares and near-sun coronal mass ejection-driven shocks. These new observations present challenges not only for mechanisms that employ direct parallel electric fields and organize maximum energies according to E/Q, but also for local diffusive and magnetic reconnection-driven acceleration models. Re-evaluation of our current understanding of the production and transport of energetic ions is necessary to understand this near-solar, current-sheet-associated population of ST ions., Comment: 4 Figures, 2 Tables

The SWEAP instrument suite on Parker Solar Probe (PSP) has detected numerous proton beams associated with coherent, circularly polarized, ion-scale waves observed by PSP's FIELDS instrument suite. Measurements during PSP Encounters 4-8 revealed pronounced complex shapes in the proton velocity distribution functions (VDFs), in which the tip of the beam undergoes strong perpendicular diffusion, resulting in VDF level contours that resemble a `hammerhead.' We refer to these proton beams, with their attendant `hammerhead' features, as the ion strahl. We present an example of these observations occurring simultaneously with a 7-hour ion-scale wave storm and show results from a preliminary attempt at quantifying the occurrence of ion-strahl broadening through 3-component ion-VDF fitting. We also provide a possible explanation of the ion perpendicular scattering based on quasilinear theory and the resonant scattering of beam ions by parallel-propagating, right circularly polarized, fast-magnetosonic/whistler waves., Comment: 13 pages, 8 figures, Under Review in the Astrophysical Journal

The hot and diffuse nature of the Sun's extended atmosphere allows it to persist in non-equilibrium states for long enough that wave-particle instabilities can arise and modify the evolution of the expanding solar wind. Determining which instabilities arise, and how significant a role they play in governing the dynamics of the solar wind, has been a decades-long process involving in situ observations at a variety of radial distances. With new measurements from Parker Solar Probe (PSP), we can study what wave modes are driven near the Sun, and calculate what instabilities are predicted for different models of the underlying particle populations. We model two hours-long intervals of PSP/SPAN-i measurements of the proton phase-space density during PSP's fourth perihelion with the Sun using two commonly used descriptions for the underlying velocity distribution. The linear stability and growth rates associated with the two models are calculated and compared. We find that both selected intervals are susceptible to resonant instabilities, though the growth rates and kind of modes driven unstable vary depending on if the protons are modeled using one or two components. In some cases, the predicted growth rates are large enough to compete with other dynamic processes, such as the nonlinear turbulent transfer of energy, in contrast with relatively slower instabilities at larger radial distances from the Sun., Comment: 10 pages, 9 figures, accepted in ApJ

Switchbacks are discrete angular deflections in the solar wind magnetic field that have been observed throughout the heliosphere. Recent observations by Parker Solar Probe (PSP) have revealed the presence of patches of switchbacks on the scale of hours to days, separated by 'quieter' radial fields. We aim to further diagnose the origin of these patches using measurements of proton temperature anisotropy that can illuminate possible links to formation processes in the solar corona. We fitted 3D bi-Maxwellian functions to the core of proton velocity distributions measured by the SPAN-Ai instrument onboard PSP to obtain the proton parallel, $T_{p,\|}$, and perpendicular, $T_{p,\perp}$, temperature. We show that the presence of patches is highlighted by a transverse deflection in the flow and magnetic field away from the radial direction. These deflections are correlated with enhancements in $T_{p,\|}$, while $T_{p,\perp}$ remains relatively constant. Patches sometimes exhibit small proton and electron density enhancements. We interpret that patches are not simply a group of switchbacks, but rather switchbacks are embedded within a larger-scale structure identified by enhanced $T_{p,\|}$ that is distinct from the surrounding solar wind. We suggest that these observations are consistent with formation by reconnection-associated mechanisms in the corona., Comment: Letter accepted as part of the PSP Special Issue in A&A

We survey the electron heat flux observed by the Parker Solar Probe (PSP) in the near-Sun environment at heliocentric distances of 0.125-0.25 AU. We utilized measurements from the Solar Wind Electrons Alphas and Protons and FIELDS experiments to compute the solar wind electron heat flux and its components and to place these in context. The PSP observations reveal a number of trends in the electron heat flux signatures near the Sun. The magnitude of the heat flux is anticorrelated with solar wind speed, likely as a result of the lower saturation heat flux in the higher-speed wind. When divided by the saturation heat flux, the resulting normalized net heat flux is anticorrelated with plasma beta on all PSP orbits, which is consistent with the operation of collisionless heat flux regulation mechanisms. The net heat flux also decreases in very high beta regions in the vicinity of the heliospheric current sheet, but in most cases of this type the omnidirectional suprathermal electron flux remains at a comparable level or even increases, seemingly inconsistent with disconnection from the Sun. The measured heat flux values appear inconsistent with regulation primarily by collisional mechanisms near the Sun. Instead, the observed heat flux dependence on plasma beta and the distribution of suprathermal electron parameters are both consistent with theoretical instability thresholds associated with oblique whistler and magnetosonic modes.

We systematically search for magnetic flux rope structures in the solar wind to within the closest distance to the Sun of 0.13 AU, using data from the third and fourth orbits of the Parker Solar Probe. We extend our previous magnetic helicity based technique of identifying magnetic flux rope structures. The method is improved upon to incorporate the azimuthal flow, which becomes larger as the spacecraft approaches the Sun. A total of 21 and 34 magnetic flux ropes are identified during the third (21 days period) and fourth (17 days period) orbits of the Parker Solar Probe, respectively. We provide a statistical analysis of the identified structures, including their relation to the streamer belt and heliospheric current sheet crossing.

Parker Solar Probe (PSP), NASA's latest and closest mission to the Sun, is on a journey to investigate fundamental enigmas of the inner heliosphere. This paper reports initial observations made by the Solar Probe Analyzer for Ions (SPAN-I), one of the instruments in the Solar Wind Electrons Alphas and Protons (SWEAP) instrument suite. We address the presence of secondary proton beams in concert with ion-scale waves observed by FIELDS, the electromagnetic fields instrument suite. We show two events from PSP's 2nd orbit that demonstrate signatures consistent with wave-particle interactions. We showcase 3D velocity distribution functions (VDFs) measured by SPAN-I during times of strong wave power at ion-scales. From an initial instability analysis, we infer that the VDFs departed far enough away from local thermodynamic equilibrium (LTE) to provide sufficient free energy to locally generate waves. These events exemplify the types of instabilities that may be present and, as such, may guide future data analysis characterizing and distinguishing between different wave-particle interactions., Comment: 24 pages, 9 figures, 2 tables

Observations by the Parker Solar Probe mission of the solar wind at about 35.7 solar radii reveal the existence of whistler wave packets with frequencies below 0.1 f/fce (20-80 Hz in the spacecraft frame). These waves often coincide with local minima of the magnetic field magnitude or with sudden deflections of the magnetic field that are called switchbacks. Their sunward propagation leads to a significant Doppler frequency downshift from 200-300 Hz to 20-80 Hz (from 0.2 f/fce to 0.5 f/fce). The polarization of these waves varies from quasi-parallel to significantly oblique with wave normal angles that are close to the resonance cone. Their peak amplitude can be as large as 2 to 4 nT. Such values represent approximately 10% of the background magnetic field, which is considerably more than what is observed at 1 a.u. Recent numerical studies show that such waves may potentially play a key role in breaking the heat flux and scattering the Strahl population of suprathermal electrons into a halo population.

Magnetic reconnection in current sheets converts magnetic energy into particle energy. The process may play an important role in the acceleration and heating of the solar wind close to the Sun. Observations from Parker Solar Probe provide a new opportunity to study this problem, as it measures the solar wind at unprecedented close distances to the Sun. During the 1st orbit, PSP encountered a large number of current sheets in the solar wind through perihelion at 35.7 solar radii. We performed a comprehensive survey of these current sheets and found evidence for 21 reconnection exhausts. These exhausts were observed in heliospheric current sheets, coronal mass ejections, and regular solar wind. However, we find that the majority of current sheets encountered around perihelion, where the magnetic field was strongest and plasma beta was lowest, were Alfv\'enic structures associated with bursty radial jets and these current sheets did not appear to be undergoing local reconnection. We examined conditions around current sheets to address why some current sheets reconnected, while others did not. A key difference appears to be the degree of plasma velocity shear across the current sheets: The median velocity shear for the 21 reconnection exhausts was 24% of the Alfv\'en velocity shear, whereas the median shear across 43 Alfv\'enic current sheets examined was 71% of the Alfv\'en velocity shear. This finding could suggest that large, albeit sub-Alfv\'enic, velocity shears suppress reconnection. An alternative interpretation is that the Alfv\'enic current sheets are isolated rotational discontinuities which do not undergo local reconnection., Comment: In Press (accepted by ApJS on 2019-11-08) This paper is part of the Parker Solar Probe ApJS Special Issue Current citation: Phan, T. D., S. D. Bale, J. P. Eastwood, et al. (2019), Parker Solar Probe In-Situ Observations of Magnetic Reconnection Exhausts During Encounter 1, ApJS., in press, doi: 10.3847/1538-4365/ab55ee

The solar wind shows periods of highly Alfv\'enic activity, where velocity fluctuations and magnetic fluctuations are aligned or anti-aligned with each other. It is generally agreed that solar wind plasma velocity and magnetic field fluctuations observed by Parker Solar Probe (PSP) during the first encounter are mostly highly Alfv\'enic. However, quantitative measures of Alfv\'enicity are needed to understand how the characterization of these fluctuations compares with standard measures from prior missions in the inner and outer heliosphere, in fast wind and slow wind, and at high and low latitudes. To investigate this issue, we employ several measures to quantify the extent of Alfv\'enicity -- the Alfv\'en ratio $r_A$, {normalized} cross helicity $\sigma_c$, {normalized} residual energy $\sigma_r$, and the cosine of angle between velocity and magnetic fluctuations $\cos\theta_{vb}$. We show that despite the overall impression that the Alfv\'enicity is large in the solar wind sampled by PSP during the first encounter, during some intervals the cross helicity starts decreasing at very large scales. These length-scales (often $> 1000 d_i$) are well inside inertial range, and therefore, the suppression of cross helicity at these scales cannot be attributed to kinetic physics. This drop at large scales could potentially be explained by large-scale shears present in the inner heliosphere sampled by PSP. In some cases, despite the cross helicity being constant down to the noise floor, the residual energy decreases with scale in the inertial range. These results suggest that it is important to consider all these measures to quantify Alfv\'enicity., Comment: Submitted to special issue of ApJ for Parker Solar Probe

The solar wind proton temperature at 1-au has been found to be correlated with small-scale intermittent magnetic structures, i.e., regions with enhanced temperature are associated with coherent structures such as current sheets. Using Parker Solar Probe data from the first encounter, we study this association using measurements of radial proton temperature, employing the Partial Variance of Increments (PVI) technique to identify intermittent magnetic structures. We observe that the probability density functions of high-PVI events have higher median temperatures than those with lower PVI, The regions in space where PVI peaks were also locations that had enhanced temperatures when compared with similar regions suggesting a heating mechanism in the young solar wind that is associated with intermittency developed by a nonlinear turbulent cascade.n the immediate vicinity., Comment: 6 pages, 3 figures, part of ApJ special issue for PSP

During the Parker Solar Probe's (PSP) first perihelion pass, the spacecraft reached within a heliocentric distance of \(\sim 37~R_\odot\) and observed numerous magnetic and flow structures characterized by sharp gradients. To better understand these intermittent structures in the young solar wind, an important property to examine is their degree of correlation in time and space. To this end, we use the well-tested Partial Variance of Increments (PVI) technique to identify intermittent events in FIELDS and SWEAP observations of magnetic and proton-velocity fields (respectively) during PSP's first solar encounter, when the spacecraft was within 0.25 au from the Sun. We then examine distributions of waiting times between events with varying separation and PVI thresholds. We find power-law distributions for waiting times shorter than a characteristic scale comparable to the correlation time, suggesting a high degree of correlation that may originate in a clustering process. Waiting times longer than this characteristic time are better described by an exponential, suggesting a random memory-less Poisson process at play. These findings are consistent with near-Earth observations of solar wind turbulence. The present study complements the one by Dudok de Wit et al. (2020, present volume), which focuses on waiting times between observed "switchbacks" in the radial magnetic field.

A solar energetic particle event was detected by the Integrated Science Investigation of the Sun (ISOIS) instrument suite on Parker Solar Probe (PSP) on 2019 April 4 when the spacecraft was inside of 0.17 au and less than 1 day before its second perihelion, providing an opportunity to study solar particle acceleration and transport unprecedentedly close to the source. The event was very small, with peak 1 MeV proton intensities of ~0.3 particles (cm^2 sr s MeV)^-1, and was undetectable above background levels at energies above 10 MeV or in particle detectors at 1 au. It was strongly anisotropic, with intensities flowing outward from the Sun up to 30 times greater than those flowing inward persisting throughout the event. Temporal association between particle increases and small brightness surges in the extreme-ultraviolet observed by the Solar TErrestrial RElations Observatory, which were also accompanied by type III radio emission seen by the Electromagnetic Fields Investigation on PSP, indicates that the source of this event was an active region nearly 80 degrees east of the nominal PSP magnetic footpoint. This suggests that the field lines expanded over a wide longitudinal range between the active region in the photosphere and the corona.

Direct evidence of an inertial-range turbulent energy cascade has been provided by spacecraft observations in heliospheric plasmas. In the solar wind, the average value of the derived heating rate near 1 au is $\sim 10^{3}\, \mathrm{J\,kg^{-1}\,s^{-1}}$, an amount sufficient to account for observed departures from adiabatic expansion. Parker Solar Probe (PSP), even during its first solar encounter, offers the first opportunity to compute, in a similar fashion, a fluid-scale energy decay rate, much closer to the solar corona than any prior in-situ observations. Using the Politano-Pouquet third-order law and the von K\'arm\'an decay law, we estimate the fluid-range energy transfer rate in the inner heliosphere, at heliocentric distance $R$ ranging from $54\,R_{\odot}$ (0.25 au) to $36\,R_{\odot}$ (0.17 au). The energy transfer rate obtained near the first perihelion is about 100 times higher than the average value at 1 au. This dramatic increase in the heating rate is unprecedented in previous solar wind observations, including those from Helios, and the values are close to those obtained in the shocked plasma inside the terrestrial magnetosheath., Comment: Accepted for publication in The Astrophysical Journal Supplement, PSP special issue

A series of solar energetic particle (SEP) events were observed at Parker Solar Probe (PSP) by the Integrated Science Investigation of the Sun (\ISOIS) during the period from April 18, 2019 through April 24, 2019. The PSP spacecraft was located near 0.48 au from the Sun on Parker spiral field lines that projected out to 1 au within $\sim 25^\circ$ of near Earth spacecraft. These SEP events, though small compared to historically large SEP events, were amongst the largest observed thus far in the PSP mission and provide critical information about the space environment inside 1 au during SEP events. During this period the Sun released multiple coronal mass ejections (CMEs). One of these CMEs observed was initiated on April 20, 2019 at 01:25 UTC, and the interplanetary CME (ICME) propagated out and passed over the PSP spacecraft. Observations by the Electromagnetic Fields Investigation (FIELDS) show that the magnetic field structure was mostly radial throughout the passage of the compression region and the plasma that followed, indicating that PSP did not directly observe a flux rope internal to the ICME, consistent with the location of PSP on the ICME flank. Analysis using relativistic electrons observed near Earth by the Electron, Proton and Alpha Monitor (EPAM) on the Advanced Composition Explorer (ACE) demonstrates the presence of electron seed populations (40--300 keV) during the events observed. The energy spectrum of the \ISOIS~ observed proton seed population below 1 MeV is close to the limit of possible stationary state plasma distributions out of equilibrium. \ISOIS~ observations reveal the \revise{enhancement} of seed populations during the passage of the ICME, which \revise{likely indicates a key part} of the pre-acceleration process that occurs close to the Sun.

Since the launch on 2018/08/12, Parker Solar Probe (PSP) has completed its first and second orbits around the Sun, having reached down to 35.7 solar radii at each perihelion. In anticipation of the exciting new data at such unprecedented distances, we have simulated the global 3D heliosphere using an MHD model coupled with a semi-empirical coronal model using the best available photospheric magnetograms as input. We compare our heliospheric MHD simulation results with in situ measurements along the PSP trajectory from its launch to the completion of the second orbit, with particular emphasis on the solar wind structure around the first two solar encounters. Furthermore, we show our model prediction for the third perihelion, which occurred on 2019/09/01. Comparison of the MHD results with PSP observations provides a new insight on the solar wind acceleration. Moreover, PSP observations reveal how accurately the ADAPT-WSA predictions work throughout the inner heliosphere.

The Parker Solar Probe (PSP) observed an interplanetary coronal mass ejection (ICME) event during its first orbit around the sun, among many other events. This event is analyzed by applying a wavelet analysis technique to obtain the reduced magnetic helicity, cross helicity, and residual energy, the first two of which are magnetohydrodynamics (MHD) invariants. Our results show that the ICME, as a large scale magnetic flux rope, possesses high magnetic helicity, very low cross helicity, and highly negative residual energy, thus pointing to a magnetic fluctuation dominated structure. Using the same technique, we also search for small-scale coherent magnetic flux rope structures during the period from 2018/10/22--2018/11/21, which are intrinsic to quasi-2D MHD turbulence in the solar wind. Multiple structures with duration between 8 and 300 minutes are identified from PSP in-situ spacecraft measurements. The location and scales of these structures are characterized by wavelet spectrograms of the normalized reduced magnetic helicity, normalized cross helicity and normalized residual energy. Transport theory suggests that these small-scale magnetic flux ropes may contribute to the acceleration of charged particles through magnetic reconnection processes, and the dissipation of these structures may be important for understanding the coronal heating processes., Comment: 22 pages, 5 figures, accepted in ApJS

The first two orbits of the Parker Solar Probe (PSP) spacecraft have enabled the first in situ measurements of the solar wind down to a heliocentric distance of 0.17 au (or 36 Rs). Here, we present an analysis of this data to study solar wind turbulence at 0.17 au and its evolution out to 1 au. While many features remain similar, key differences at 0.17 au include: increased turbulence energy levels by more than an order of magnitude, a magnetic field spectral index of -3/2 matching that of the velocity and both Elsasser fields, a lower magnetic compressibility consistent with a smaller slow-mode kinetic energy fraction, and a much smaller outer scale that has had time for substantial nonlinear processing. There is also an overall increase in the dominance of outward-propagating Alfv\'enic fluctuations compared to inward-propagating ones, and the radial variation of the inward component is consistent with its generation by reflection from the large-scale gradient in Alfv\'en speed. The energy flux in this turbulence at 0.17 au was found to be ~10% of that in the bulk solar wind kinetic energy, becoming ~40% when extrapolated to the Alfv\'en point, and both the fraction and rate of increase of this flux towards the Sun is consistent with turbulence-driven models in which the solar wind is powered by this flux.

Steady non-equilibrium states are investigated in a one-dimensional setup in the presence of two thermodynamic currents. Two paradigmatic nonlinear oscillators models are investigated: an XY chain and the discrete nonlinear Schr\"odinger equation. Their distinctive feature is that the relevant variable is an angle in both cases. We point out the importance of clearly distinguishing between energy and heat flux. In fact, even in the presence of a vanishing Seebeck coefficient, a coupling between (angular) momentum and energy arises, mediated by the unavoidable presence of a "coherent" energy flux. Such a contribution is the result of the "advection" induced by the position-dependent angular velocity. As a result, in the XY model, the knowledge of the two diagonal elements of the Onsager matrix suffices to reconstruct its transport properties. The analysis of the nonequilibrium steady states finally allows to strengthen the connection between the two models.

We study the dynamics of networks with inhibitory and excitatory leaky-integrate-and-fire neurons with short-term synaptic plasticity in the presence of depressive and facilitating mechanisms. The dynamics is analyzed by a Heterogeneous Mean-Field approximation, that allows to keep track of the effects of structural disorder in the network. We describe the complex behavior of different classes of excitatory and inhibitory components, that give rise to a rich dynamical phase-diagram as a function of the fraction of inhibitory neurons. By the same mean field approach, we study and solve a global inverse problem: reconstructing the degree probability distributions of the inhibitory and excitatory components and the fraction of inhibitory neurons from the knowledge of the average synaptic activity field. This approach unveils new perspectives in the numerical study of neural network dynamics and in the possibility of using these models as testbed for the analysis of experimental data.

Suitable Langevin thermostats are introduced which are able to control both the temperature and the chemical potential of a one-dimensional lattice of nonlinear Schr\"odinger oscillators. The resulting non-equilibrium stationary states are then investigated in two limit cases (low temperatures and large particle densities), where the dynamics can be mapped onto that of a coupled-rotor chain with an external torque. As a result, an effective kinetic definition of temperature can be introduced and compared with the general microcanonical (global) definition., Comment: Accepted for publication in JSTAT

We explore the statistical behavior of the discrete nonlinear Schroedinger equation. We find a parameter region where the system evolves towards a state characterized by a finite density of breathers and a negative temperature. Such a state is metastable but the convergence to equilibrium occurs on astronomical time scales and becomes increasingly slower as a result of a coarsening processes. Stationary negative-temperature states can be experimentally generated via boundary dissipation or from free expansions of wave packets initially at positive temperature equilibrium., Comment: 4 pages, 5 figures

The formation of amplitude modulated and phase modulated assemblies of neurons is observed in the brain functional activity. The study of the formation of such structures requires that the analysis has to be organized in hierarchical levels, microscopic, mesoscopic, macroscopic, each with its characteristic space-time scales and the various forms of energy, electric, chemical, thermal produced and used by the brain. In this paper, we discuss the microscopic dynamics underlying the mesoscopic and the macroscopic levels and focus our attention on the thermodynamics of the non-equilibrium phase transitions. We obtain the time-dependent Ginzburg-Landau equation for the non-stationary regime and consider the formation of topologically non-trivial structures such as the vortex solution. The power laws observed in functional activities of the brain is also discussed and related to coherent states characterizing the many-body dissipative model of brain., Comment: 19 pages, 4 figures, research paper

The fluctuation theorem is a pivotal result of statistical physics. It quantifies the probability of observing fluctuations which are in violation of the second law of thermodynamics. More specifically, it quantifies the ratio of the probabilities of observing entropy-producing and entropy-consuming fluctuations measured over a finite volume and time span in terms of the rate of entropy production in the system, the measurement volume and time. We study the fluctuation theorem in computer simulations of planar shear flow. The simulations are performed employing the method of multiparticle collision dynamics which captures both thermal fluctuations and hydrodynamic interactions. The main outcome of our analysis is that the fluctuation theorem is verified at any averaging time provided that the measurement volume exhibits a specific dependence on a hydrodynamic time scale., Comment: 4 pages, 3 figures, to appear on Physical Review Letters

We study heat transport in a chain of harmonic oscillators with random elastic collisions between nearest-neighbours. The equations of motion of the covariance matrix are numerically solved for free and fixed boundary conditions. In the thermodynamic limit, the shape of the temperature profile and the value of the stationary heat flux depend on the choice of boundary conditions. For free boundary conditions, they also depend on the coupling strength with the heat baths. Moreover, we find a strong violation of local equilibrium at the chain edges that determine two boundary layers of size $\sqrt{N}$ (where $N$ is the chain length), that are characterized by a different scaling behaviour from the bulk. Finally, we investigate the relaxation towards the stationary state, finding two long time scales: the first corresponds to the relaxation of the hydrodynamic modes; the second is a manifestation of the finiteness of the system., Comment: Submitted to Journal of Physics A, Mathematical and Theoretical

A method for reconstructing the energy landscape of simple polypeptidic chains is described. We show that we can construct an equivalent representation of the energy landscape by a suitable directed graph. Its topological and dynamical features are shown to yield an effective estimate of the time scales associated with the folding and with the equilibration processes. This conclusion is drawn by comparing molecular dynamics simulations at constant temperature with the dynamics on the graph, defined by a temperature dependent Markov process. The main advantage of the graph representation is that its dynamics can be naturally renormalized by collecting nodes into "hubs", while redefining their connectivity. We show that both topological and dynamical properties are preserved by the renormalization procedure. Moreover, we obtain clear indications that the heteropolymers exhibit common topological properties, at variance with the homopolymer, whose peculiar graph structure stems from its spatial homogeneity. In order to obtain a clear distinction between a "fast folder" and a "slow folder" in the heteropolymers one has to look at kinetic features of the directed graph. We find that the average time needed to the fast folder for reaching its native configuration is two orders of magnitude smaller than its equilibration time, while for the bad folder these time scales are comparable. Accordingly, we can conclude that the strategy described in this paper can be successfully applied also to more realistic models, by studying their renormalized dynamics on the directed graph, rather than performing lengthy molecular dynamics simulations., Comment: 15 pages, 12 figures

We provide an explicit representation of the nonequilibrium invariant measure for a chain of harmonic oscillators with conservative noise in the presence of stationary heat flow. By first determining the covariance matrix, we are able to express the measure as the product of Gaussian distributions aligned along some collective modes that are spatially localized with power-law tails. Numerical studies show that such a representation applies also to a purely deterministic model, the quartic Fermi-Pasta-Ulam chain.

After reviewing the main features of anomalous energy transport in 1D systems, we report simulations performed with chains of noisy anharmonic oscillators. The stochastic terms are added in such a way to conserve total energy and momentum, thus keeping the basic hydrodynamic features of these models. The addition of this "conservative noise" allows to obtain a more efficient estimate of the power-law divergence of heat conductivity kappa(L) ~ L^alpha in the limit of small noise and large system size L. By comparing the numerical results with rigorous predictions obtained for the harmonic chain, we show how finite--size and --time effects can be effectively controlled. For low noise amplitudes, the alpha values are close to 1/3 for asymmetric potentials and to 0.4 for symmetric ones. These results support the previously conjectured two-universality-classes scenario.

A general method to determine covariant Lyapunov vectors in both discrete- and continuous-time dynamical systems is introduced. This allows to address fundamental questions such as the degree of hyperbolicity, which can be quantified in terms of the transversality of these intrinsic vectors. For spatially extended systems, the covariant Lyapunov vectors have localization properties and spatial Fourier spectra qualitatively different from those composing the orthonormalized basis obtained in the standard procedure used to calculate the Lyapunov exponents., Comment: 4 pages, 3 figures, submitted to Physical Review letters

We investigate the diffusive properties of energy fluctuations in a one-dimensional diatomic chain of hard-point particles interacting through a square--well potential. The evolution of initially localized infinitesimal and finite perturbations is numerically investigated for different density values. All cases belong to the same universality class which can be also interpreted as a Levy walk of the energy with scaling exponent 3/5. The zero-pressure limit is nevertheless exceptional in that normal diffusion is found in tangent space and yet anomalous diffusion with a different rate for perturbations of finite amplitude. The different behaviour of the two classes of perturbations is traced back to the "stable chaos" type of dynamics exhibited by this model. Finally, the effect of an additional internal degree of freedom is investigated, finding that it does not modify the overall scenario, Comment: 16 pages, 15 figures

Thermodynamical properties of the helix-coil transition were successfully described in the past by the model of Lifson, Poland and Sheraga. Here we compute the corresponding structure factor and show that it possesses a universal scaling behavior near the transition point, even when the transition is of first order. Moreover, we introduce a dynamical version of this model, that we solve numerically. A Langevin equation is also proposed to describe the dynamics of the density of hydrogen bonds. Analytical solution of this equation shows dynamical scaling near the critical temperature and predicts a gelation phenomenon above the critical temperature. In the case when comparison of the two dynamical approaches is possible, the predictions of our phenomenological theory agree with the results of the Monte Carlo simulations., Comment: 11 pages, 7 figures

A class of non-local contact processes is introduced and studied using mean-field approximation and numerical simulations. In these processes particles are created at a rate which decays algebraically with the distance from the nearest particle. It is found that the transition into the absorbing state is continuous and is characterized by continuously varying critical exponents. This model differs from the previously studied non-local directed percolation model, where particles are created by unrestricted Levy flights. It is motivated by recent studies of non-equilibrium wetting indicating that this type of non-local processes play a role in the unbinding transition. Other non-local processes which have been suggested to exist within the context of wetting are considered as well., Comment: Accepted with minor revisions by Journal of Statistical Mechanics: Theory and experiments

The dynamical behaviour of a weakly diluted fully-inhibitory network of pulse-coupled spiking neurons is investigated. Upon increasing the coupling strength, a transition from regular to stochastic-like regime is observed. In the weak-coupling phase, a periodic dynamics is rapidly approached, with all neurons firing with the same rate and mutually phase-locked. The strong-coupling phase is characterized by an irregular pattern, even though the maximum Lyapunov exponent is negative. The paradox is solved by drawing an analogy with the phenomenon of ``stable chaos'', i.e. by observing that the stochastic-like behaviour is "limited" to a an exponentially long (with the system size) transient. Remarkably, the transient dynamics turns out to be stationary., Comment: 11 pages, 13 figures, submitted to Phys. Rev. E

In the present Letter we present an analytical and numerical solution of the self-consistent mode-coupling equations for the problem of heat conductivity in one-dimensional systems. Such a solution leads us to propose a different scenario to accomodate the known results obtained so far for this problem. More precisely, we conjecture that the universality class is determined by the leading order of the nonlinear interaction potential. Moreover, our analysis allows us determining the memory kernel, whose expression puts on a more firm basis the previously conjectured connection between anomalous heat conductivity and anomalous diffusion., Comment: Submitted to Physical Review E

The random forced Navier-Stokes equation can be obtained as a variational problem of a proper action. By virtue of incompressibility, the integration over transverse components of the fields allows to cast the action in the form of a large deviation functional. Since the hydrodynamic operator is nonlinear, the functional integral yielding the statistics of fluctuations can be practically computed by linearizing around a physical solution of the hydrodynamic equation. We show that this procedure yields the dimensional scaling predicted by K41 theory at the lowest perturbative order, where the perturbation parameter is the inverse Reynolds number. Moreover, an explicit expression of the prefactor of the scaling law is obtained., Comment: 24 pages

It is argued that some phase--transitions observed in models of non-equilibrium wetting phenomena are related to contact processes with long-range interactions. This is investigated by introducing a model where the activation rate of a site at the edge of an inactive island of length $\ell$ is $1+a\ell^{-\sigma}$. Mean--field analysis and numerical simulations indicate that for $\sigma>1$ the transition is continuous and belongs to the universality class of directed percolation, while for $0<\sigma<1$, the transition becomes first order. This criterion is then applied to discuss critical properties of various models of non--equilibrium wetting., Comment: 12 Figures. V Version resubmitted to Phys. Rev. E after referees report