Source author record

R. Livi

R. Livi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

28works
14topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2022arXiv

Density And Velocity Fluctuations of Alpha Particles in Magnetic Switchbacks

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_{αp}$, inside vs outside, nor do we observe any signature when separating the switchbacks according to $V_{α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_{αp}$, with the Alfvénic motion of the alpha particles vanishing as the difference $|V_{pw} - V_{αp}|$ decreases. We show how the Alfvénic 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.

preprint2022arXiv

Sub-Alfvenic Solar Wind observed by PSP: Characterization of Turbulence, Anisotropy, Intermittency, and Switchback

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.

preprint2022arXiv

The Radial Evolution of the Solar Wind as Organized by Electron Distribution Parameters

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.

preprint2021arXiv

Inferred Linear Stability of Parker Solar Probe Observations using One- and Two-Component Proton Distributions

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.

preprint2021arXiv

Strong perpendicular velocity-space in proton beams observed by Parker Solar Probe

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.

preprint2021arXiv

Suprathermal Ion Energy spectra and Anisotropies near the Heliospheric Current Sheet crossing observed by the Parker Solar Probe during Encounter 7

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.

preprint2020arXiv

Parker Solar Probe In-Situ Observations of Magnetic Reconnection Exhausts During Encounter 1

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énic 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én velocity shear, whereas the median shear across 43 Alfvénic current sheets examined was 71% of the Alfvén velocity shear. This finding could suggest that large, albeit sub-Alfvénic, velocity shears suppress reconnection. An alternative interpretation is that the Alfvénic current sheets are isolated rotational discontinuities which do not undergo local reconnection.

preprint2020arXiv

Parker Solar Probe observations of proton beams simultaneous with ion-scale waves

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.

preprint2020arXiv

Sunward propagating whistler waves collocated with localized magnetic field holes in the solar wind: Parker Solar Probe observations at 35.7 Sun radii

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.

preprint2019arXiv

Clustering of Intermittent Magnetic and Flow Structures near Parker Solar Probe's First Perihelion -- A Partial-Variance-of-Increments Analysis

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.

preprint2019arXiv

Enhanced Energy Transfer Rate in Solar Wind Turbulence Observed near the Sun from Parker Solar Probe

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ármán 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.

preprint2019arXiv

Identification of Magnetic Flux Ropes from Parker Solar Probe Observations during the First Encounter

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.

preprint2019arXiv

Observations of heating along intermittent structures in the inner heliosphere from PSP data

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.

preprint2019arXiv

Observations of the 2019 April 4 Solar Energetic Particle Event at the Parker Solar Probe

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.

preprint2019arXiv

Predicting the Solar Wind at Parker Solar Probe Using an Empirically Driven MHD Model

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.

preprint2019arXiv

Seed Population Pre-Conditioning and Acceleration Observed by Parker Solar Probe

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.

preprint2019arXiv

The Evolution and Role of Solar Wind Turbulence in the Inner Heliosphere

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énic 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én 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én 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.

preprint2016arXiv

Coupled transport in rotor models

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ödinger 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.

preprint2015arXiv

Neural networks with excitatory and inhibitory components: direct and inverse problems by a mean-field approach

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.

preprint2013arXiv

Off-equilibrium Langevin dynamics of the discrete nonlinear Schroedinger chain

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ödinger 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.

preprint2012arXiv

Negative Temperature States in the Discrete Nonlinear Schroedinger Equation

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.

preprint2011arXiv

Cortical phase transitions, non-equilibrium thermodynamics and the time-dependent Ginzburg-Landau equation

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.

preprint2011arXiv

Hydrodynamics and the fluctuation theorem

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.

preprint2009arXiv

Nonequilibrium dynamics of a stochastic model of anomalous heat transport: numerical analysis

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.

preprint2008arXiv

Stochastic dynamics of model proteins on a directed graph

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.

preprint2002arXiv

On the relationship between directed percolation and the synchronization transition in spatially extended systems

We study the nature of the synchronization transition in spatially extended systems by discussing a simple stochastic model. An analytic argument is put forward showing that, in the limit of discontinuous processes, the transition belongs to the directed percolation (DP) universality class. The analysis is complemented by a detailed investigation of the dependence of the first passage time for the amplitude of the difference field on the adopted threshold. We find the existence of a critical threshold separating the regime controlled by linear mechanisms from that controlled by collective phenomena. As a result of this analysis we conclude that the synchronization transition belongs to the DP class also in continuous models. The conclusions are supported by numerical checks on coupled map lattices too.

preprint2001arXiv

Slow energy relaxation and localization in 1D lattices

We investigate the energy relaxation process produced by thermal baths at zero temperature acting on the boundary atoms of chains of classical anharmonic oscillators. Time-dependent perturbation theory allows us to obtain an explicit solution of the harmonic problem: even in such a simple system nontrivial features emerge from the interplay of the different decay rates of Fourier modes. In particular, a crossover from an exponential to an inverse-square-root law occurs on a time scale proportional to the system size $N$. A further crossover back to an exponential law is observed only at much longer times (of the order $N^3$). In the nonlinear chain, the relaxation process is initially equivalent to the harmonic case over a wide time span, as illustrated by simulations of the $β$ Fermi-Pasta-Ulam model. The distinctive feature is that the second crossover is not observed due to the spontaneous appearance of breathers, i.e. space-localized time-periodic solutions, that keep a finite residual energy in the lattice. We discuss the mechanism yielding such solutions and also explain why it crucially depends on the boundary conditions.

preprint1999arXiv

Noise-driven Synchronization in Coupled Map Lattices

Synchronization is shown to occur in spatially extended systems under the effect of additive spatio-temporal noise. In analogy to low dimensional systems, synchronized states are observable only if the maximum Lyapunov exponent $Λ$ is negative. However, a sufficiently high noise level can lead, in map with finite domain of definition, to nonlinear propagation of information, even in non chaotic systems. In this latter case the transition to synchronization is ruled by a new ingredient : the propagation velocity of information $V_F$. As a general statement, we can affirm that if $V_F$ is finite the time needed to achieve a synchronized trajectory grows exponentially with the system size $L$, while it increases logarithmically with $L$ when, for sufficiently large noise amplitude, $V_F = 0$ .