Showing 1–46 of 46 results for author: Nicolaou, G

Observation of a Fully-formed Forward--Reverse Shock Pair Due to the Interaction Between Two Coronal Mass Ejections at 0.5 au

Authors: D. Trotta, A. Dimmock, X. Blanco-Cano, R. Forsyth, H. Hietala, N. Fargette, A. Larosa, N. Lugaz, E. Palmerio, S. W. Good, J. E. Soljento, E. K. J. Kilpua, E. Yordanova, O. Pezzi, G. Nicolaou, T. S. Horbury, R. Vainio, N. Dresing, C. J. Owen, R. Wimmer-Schweingruber

Abstract: We report direct observations of a fast magnetosonic forward--reverse shock pair observed by Solar Orbiter on March 8, 2022 at the short heliocentric distance of 0.5 au. The structure, sharing some features with fully-formed stream interaction regions (SIRs), is due to the interaction between two successive coronal mass ejections (CMEs), never previously observed to give rise to a forward--reverse… ▽ More We report direct observations of a fast magnetosonic forward--reverse shock pair observed by Solar Orbiter on March 8, 2022 at the short heliocentric distance of 0.5 au. The structure, sharing some features with fully-formed stream interaction regions (SIRs), is due to the interaction between two successive coronal mass ejections (CMEs), never previously observed to give rise to a forward--reverse shock pair. The scenario is supported by remote observations from the STEREO-A coronographs, where two candidate eruptions compatible with the in-situ signatures have been found. In the interaction region, we find enhanced energetic particle activity, strong non-radial flow deflections and evidence of magnetic reconnection. At 1~au, well radially-aligned \textit{Wind} observations reveal a complex event, with characteristic observational signatures of both SIR and CME--CME interaction, thus demonstrating the importance of investigating the complex dynamics governing solar eruptive phenomena. △ Less

Submitted 1 August, 2024; v1 submitted 26 April, 2024; originally announced April 2024.

Comments: Accepted in The Astrophysical Journal Letters

The impact of non-equilibrium plasma distributions on solar wind measurements by Vigil's Plasma Analyser

Authors: Hongjie Zhang, Daniel Verscharen, Georgios Nicolaou

Abstract: In order to protect society from space weather impacts, we must monitor space weather and obtain early warnings for extreme events if possible. For this purpose, the European Space Agency is currently preparing to launch the Vigil mission towards the end of this decade as a space-weather monitor at the fifth Lagrange point of the Sun--Earth system. Vigil will carry, amongst other instruments, the… ▽ More In order to protect society from space weather impacts, we must monitor space weather and obtain early warnings for extreme events if possible. For this purpose, the European Space Agency is currently preparing to launch the Vigil mission towards the end of this decade as a space-weather monitor at the fifth Lagrange point of the Sun--Earth system. Vigil will carry, amongst other instruments, the Plasma Analyzer (PLA) to provide quasi-continuous measurements of solar wind ions. We model the performance of the PLA instrument, considering typical solar wind plasma conditions, to compare the expected observations of PLA with the assumed input conditions of the solar wind. We evaluate the instrument performance under realistic, non-equilibrium plasma conditions, accounting for temperature anisotropies, proton beams, and the contributions from drifting $α$-particles. We examine the accuracy of the instrument's performance over a range of input solar wind moments. We identify sources of potential errors due to non-equilibrium plasma conditions and link these to instrument characteristics such as its angular and energy resolution and its field of view. We demonstrate the limitations of the instrument and potential improvements such as applying ground-based fitting techniques to obtain more accurate measurements of the solar wind even under non-equilibrium plasma conditions. The use of ground processing of plasma moments instead of on-board processing is crucial for the extraction of reliable measurements. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 26 pages, 11 figures; accepted for publication in Space Weather

Probing Turbulent Scattering Effects on Suprathermal Electrons in the Solar Wind: Modeling, Observations and Implications

Authors: Arnaud Zaslavsky, Justin C. Kasper, Eduard P. Kontar, Davin E. Larson, Milan Maksimovic, José M. D. C. Marques, Georgios Nicolaou, Christopher J. Owen, Orlando Romeo, Phyllis L. Whittlesey

Abstract: This study explores the impact of a turbulent scattering mechanism, akin to those influencing solar and galactic cosmic rays propagating in the interplanetary medium, on the population of suprathermal electrons in the solar wind. We employ a Fokker-Planck equation to model the radial evolution of electron pitch angle distributions under the action of magnetic focusing, which moves the electrons aw… ▽ More This study explores the impact of a turbulent scattering mechanism, akin to those influencing solar and galactic cosmic rays propagating in the interplanetary medium, on the population of suprathermal electrons in the solar wind. We employ a Fokker-Planck equation to model the radial evolution of electron pitch angle distributions under the action of magnetic focusing, which moves the electrons away from isotropy, and of a diffusion process that tends to bring them back to it. We compare the steady-state solutions of this Fokker-Planck equation with data obtained from the Solar Orbiter and Parker Solar Probe missions and find a remarkable agreement, varying the turbulent mean free path as the sole free parameter in our model. The obtained mean free paths are of the order of the astronomical unit, and display weak dependence on electron energy within the $100$ eV to $1$ keV range. This value is notably lower than Coulomb collision estimates but aligns well with observed mean free paths of low-rigidity solar energetic particles events. The strong agreement between our model and observations leads us to conclude that the hypothesis of turbulent scattering at work on electrons at all heliospheric distances is justified. We discuss several implications, notably the existence of a low Knudsen number region at large distances from the Sun, which offers a natural explanation for the presence of an isotropic ``halo'' component at all distances from the Sun -- electrons being isotropized in this distant region before travelling back into the inner part of the interplanetary medium. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: 23 pages, 13 figures, 1 table

Properties of an interplanetary shock observed at 0.07 and 0.7 Astronomical Units by Parker Solar Probe and Solar Orbiter

Authors: D. Trotta, A. Larosa, G. Nicolaou, T. S. Horbury, L. Matteini, H. Hietala, X. Blanco-Cano, L. Franci, C. H. K. Chen, L. Zhao, G. P. Zank, C. M. S. Cohen, S. D. Bale, R. Laker, N. Fargette, F. Valentini, Y. Khotyaintsev, R. Kieokaew, N. Raouafi, E. Davies, R. Vainio, N. Dresing, E. Kilpua, T. Karlsson, C. J. Owen , et al. (1 additional authors not shown)

Abstract: The Parker Solar Probe (PSP) and Solar Orbiter (SolO) missions opened a new observational window in the inner heliosphere, which is finally accessible to direct measurements. On September 05, 2022, a coronal mass ejection (CME)-driven interplanetary (IP) shock has been observed as close as 0.07 au by PSP. The CME then reached SolO, which was well radially-aligned at 0.7 au, thus providing us with… ▽ More The Parker Solar Probe (PSP) and Solar Orbiter (SolO) missions opened a new observational window in the inner heliosphere, which is finally accessible to direct measurements. On September 05, 2022, a coronal mass ejection (CME)-driven interplanetary (IP) shock has been observed as close as 0.07 au by PSP. The CME then reached SolO, which was well radially-aligned at 0.7 au, thus providing us with the opportunity to study the shock properties at so different heliocentric distances. We characterize the shock, investigate its typical parameters and compare its small-scale features at both locations. Using the PSP observations, we investigate how magnetic switchbacks and ion cyclotron waves are processed upon shock crossing. We find that switchbacks preserve their V--B correlation while compressed upon the shock passage, and that the signature of ion cyclotron waves disappears downstream of the shock. By contrast, the SolO observations reveal a very structured shock transition, with a population of shock-accelerated protons of up to about 2 MeV, showing irregularities in the shock downstream, which we correlate with solar wind structures propagating across the shock. At SolO, we also report the presence of low-energy ($\sim$ 100 eV) electrons scattering due to upstream shocklets. This study elucidates how the local features of IP shocks and their environments can be very different as they propagate through the heliosphere. △ Less

Submitted 10 December, 2023; originally announced December 2023.

Comments: In review in ApJ

Solar Electron Beam -- Langmuir Wave Interactions and How They Modify Solar Electron Beam Spectra: Solar Orbiter Observations of a Match Made in the Heliosphere

Authors: Camille Y. Lorfing, Hamish A. S. Reid, Raul Gomez-Herrero, Milan Maksimovic, Georgios Nicolaou, Christopher J. Owen, Javier Rodriguez-Pacheco, Daniel F. Ryan, Domenico Trotta, Daniel Verscharen

Abstract: Solar Orbiter's four in-situ instruments have recorded numerous energetic electron events at heliocentric distances between 0.5 and 1au. We analyse energetic electron fluxes, spectra, pitch angle distributions, associated Langmuir waves, and type III solar radio bursts for 3 events to understand what causes modifications in the electron flux and identify the origin and characteristics of features… ▽ More Solar Orbiter's four in-situ instruments have recorded numerous energetic electron events at heliocentric distances between 0.5 and 1au. We analyse energetic electron fluxes, spectra, pitch angle distributions, associated Langmuir waves, and type III solar radio bursts for 3 events to understand what causes modifications in the electron flux and identify the origin and characteristics of features observed in the electron spectrum. We investigate what electron beam properties and solar wind conditions are associated with Langmuir wave growth and spectral breaks in the electron peak flux as a function of energy. We observe velocity dispersion and quasilinear relaxation in the electron flux caused by the resonant wave-particle interactions in the deca-keV range, at the energies at which we observe breaks in the electron spectrum, co-temporal with the local generation of Langmuir waves. We show, via the evolution of the electron flux at the time of the event, that these interactions are responsible for the spectral signatures observed around 10 and 50keV, confirming the results of simulations by Kontar & Reid (2009). These signatures are independent of pitch angle scattering. Our findings highlight the importance of using overlapping FOVs when working with data from different sensors. In this work, we exploit observations from all in-situ instruments to address, for the first time, how the energetic electron flux is modified by the beam-plasma interactions, and results into specific features to appear in the local spectrum. Our results, corroborated with numerical simulations, can be extended to a wider range of heliocentric distances. △ Less

Submitted 24 November, 2023; originally announced November 2023.

Complex localization mechanisms in networks of coupled oscillators: two case studies

Authors: Zachary G. Nicolaou, Jason J. Bramburger

Abstract: Localized phenomena abound in nature and throughout the physical sciences. Some universal mechanisms for localization have been characterized, such as in the snaking bifurcations of localized steady states in pattern-forming partial differential equations. While much of this understanding has been targeted at steady states, recent studies have noted complex dynamical localization phenomena in syst… ▽ More Localized phenomena abound in nature and throughout the physical sciences. Some universal mechanisms for localization have been characterized, such as in the snaking bifurcations of localized steady states in pattern-forming partial differential equations. While much of this understanding has been targeted at steady states, recent studies have noted complex dynamical localization phenomena in systems of coupled oscillators. These localized states can come in the form of symmetry-breaking chimera patterns that exhibit a coexistence of coherence and incoherence in symmetric networks of coupled oscillators and gap solitons emerging in the band gap of parametrically driven networks of oscillators. Here, we report detailed numerical continuations of localized time-periodic states in systems of coupled oscillators, while also documenting the numerous bifurcations they give way to. We find novel routes to localization involving bifurcations of heteroclinic cycles in networks of Janus oscillators and strange bifurcation diagrams resembling chaotic tangles in a parametrically driven array of coupled pendula. We highlight the important role of discrete symmetries and the symmetric branch points that emerge in symmetric models. △ Less

Submitted 16 February, 2024; v1 submitted 15 September, 2023; originally announced September 2023.

Comments: 15 pages, 9 figures

Journal ref: Chaos 34, 013131 (2024)

Prevalence of multistability and nonstationarity in driven chemical networks

Authors: Zachary G. Nicolaou, Schuyler B. Nicholson, Adilson E. Motter, Jason R. Green

Abstract: External flows of energy, entropy, and matter can cause sudden transitions in the stability of biological and industrial systems, fundamentally altering their dynamical function. How might we control and design these transitions in chemical reaction networks? Here, we analyze transitions giving rise to complex behavior in random reaction networks subject to external driving forces. In the absence… ▽ More External flows of energy, entropy, and matter can cause sudden transitions in the stability of biological and industrial systems, fundamentally altering their dynamical function. How might we control and design these transitions in chemical reaction networks? Here, we analyze transitions giving rise to complex behavior in random reaction networks subject to external driving forces. In the absence of driving, we characterize the uniqueness of the steady state and identify the percolation of a giant connected component in these networks as the number of reactions increases. When subject to chemical driving (influx and outflux of chemical species), the steady state can undergo bifurcations, leading to multistability or oscillatory dynamics. By quantifying the prevalence of these bifurcations, we show how chemical driving and network sparsity tend to promote the emergence of these complex dynamics and increased rates of entropy production. We show that catalysis also plays an important role in the emergence of complexity, strongly correlating with the prevalence of bifurcations. Our results suggest that coupling a minimal number of chemical signatures with external driving can lead to features present in biochemical processes and abiogenesis. △ Less

Submitted 15 June, 2023; originally announced June 2023.

Comments: 12 pages, 4 figures

Journal ref: J. Chem. Phys. 158, 225101 (2023)

Benchmarking sparse system identification with low-dimensional chaos

Authors: Alan A. Kaptanoglu, Lanyue Zhang, Zachary G. Nicolaou, Urban Fasel, Steven L. Brunton

Abstract: Sparse system identification is the data-driven process of obtaining parsimonious differential equations that describe the evolution of a dynamical system, balancing model complexity and accuracy. There has been rapid innovation in system identification across scientific domains, but there remains a gap in the literature for large-scale methodological comparisons that are evaluated on a variety of… ▽ More Sparse system identification is the data-driven process of obtaining parsimonious differential equations that describe the evolution of a dynamical system, balancing model complexity and accuracy. There has been rapid innovation in system identification across scientific domains, but there remains a gap in the literature for large-scale methodological comparisons that are evaluated on a variety of dynamical systems. In this work, we systematically benchmark sparse regression variants by utilizing the dysts standardized database of chaotic systems. In particular, we demonstrate how this open-source tool can be used to quantitatively compare different methods of system identification. To illustrate how this benchmark can be utilized, we perform a large comparison of four algorithms for solving the sparse identification of nonlinear dynamics (SINDy) optimization problem, finding strong performance of the original algorithm and a recent mixed-integer discrete algorithm. In all cases, we used ensembling to improve the noise robustness of SINDy and provide statistical comparisons. In addition, we show very compelling evidence that the weak SINDy formulation provides significant improvements over the traditional method, even on clean data. Lastly, we investigate how Pareto-optimal models generated from SINDy algorithms depend on the properties of the equations, finding that the performance shows no significant dependence on a set of dynamical properties that quantify the amount of chaos, scale separation, degree of nonlinearity, and the syntactic complexity. △ Less

Submitted 4 February, 2023; originally announced February 2023.

Data-driven discovery and extrapolation of parameterized pattern-forming dynamics

Authors: Zachary G. Nicolaou, Guanyu Huo, Yihui Chen, Steven L. Brunton, J. Nathan Kutz

Abstract: Pattern-forming systems can exhibit a diverse array of complex behaviors as external parameters are varied, enabling a variety of useful functions in biological and engineered systems. First-principles derivations of the underlying transitions can be characterized using bifurcation theory on model systems whose governing equations are known. In contrast, data-driven methods for more complicated an… ▽ More Pattern-forming systems can exhibit a diverse array of complex behaviors as external parameters are varied, enabling a variety of useful functions in biological and engineered systems. First-principles derivations of the underlying transitions can be characterized using bifurcation theory on model systems whose governing equations are known. In contrast, data-driven methods for more complicated and realistic systems whose governing evolution dynamics are unknown have only recently been developed. Here we develop a data-driven approach, the {\em sparse identification for nonlinear dynamics with control parameters} (SINDyCP), to discover dynamics for systems with adjustable control parameters, such as an external driving strength. We demonstrate the method on systems of varying complexity, ranging from discrete maps to systems of partial differential equations. To mitigate the impact of measurement noise, we also develop a weak formulation of SINDyCP and assess its performance on noisy data. We demonstrate applications including the discovery of universal pattern-formation equations, and their bifurcation dependencies, directly from data accessible from experiments and the extrapolation of predictions beyond the weakly nonlinear regime near the onset of an instability. △ Less

Submitted 16 November, 2023; v1 submitted 6 January, 2023; originally announced January 2023.

Comments: 6 pages, 4 figures, plus supplement

Journal ref: Phys. Rev. Research 5, L042017 (2023)

The polytropic behavior of Solar Wind protons as observed by Ulysses spacecraft during Solar minimum

Authors: Georgios Nicolaou, George Livadiotis, David J. McComas

Abstract: We analyze proton bulk parameters derived from Ulysses observations and investigate the polytropic behavior of solar wind protons over a wide range of heliocentric distances and latitudes. The large-scale variations of the proton density and temperature over heliocentric distance, indicate that plasma protons are governed by sub-adiabatic processes (polytropic index $γ$ < 5/3), if we assume proton… ▽ More We analyze proton bulk parameters derived from Ulysses observations and investigate the polytropic behavior of solar wind protons over a wide range of heliocentric distances and latitudes. The large-scale variations of the proton density and temperature over heliocentric distance, indicate that plasma protons are governed by sub-adiabatic processes (polytropic index $γ$ < 5/3), if we assume protons with three effective kinetic degrees of freedom. From the correlation between the small-scale variations of the plasma density and temperature in selected sub-intervals, we derive a polytropic index $γ$ ~ 1.4 on average. Further examination shows that the polytropic index does not depend on the solar wind speed. This agrees with the results of previous analyses of solar wind protons at ~ 1 au. We find that the polytropic index varies slightly over the range of the heliocentric distances and heliographic latitudes explored by Ulysses. We also show that the homogeneity of the plasma and the accuracy of the polytropic model applied to the data-points vary over Ulysses orbit. We compare our results with the results of previous studies which derive the polytropic index of solar wind ions within the heliosphere using observations from various spacecraft. We finally discuss the implications of our findings in terms of heating mechanisms and the effective degrees of freedom of the plasma protons. △ Less

Submitted 28 December, 2022; originally announced December 2022.

Thermal energy budget of electrons in the inner heliosphere: Parker Solar Probe Observations

Authors: Joel B. Abraham, Daniel Verscharen, Robert T. Wicks, Jefferson A. Agudelo Rueda, Christopher J. Owen, Georgios Nicolaou, Seong-Yeop Jeong

Abstract: We present an observational analysis of the electron thermal energy budget using data from Parker Solar Probe. We use the macroscopic moments, obtained from our fits to the measured electron distribution function, to evaluate the thermal energy budget based on the second moment of the Boltzmann equation. We separate contributions to the overall budget from reversible and irreversible processes. We… ▽ More We present an observational analysis of the electron thermal energy budget using data from Parker Solar Probe. We use the macroscopic moments, obtained from our fits to the measured electron distribution function, to evaluate the thermal energy budget based on the second moment of the Boltzmann equation. We separate contributions to the overall budget from reversible and irreversible processes. We find that an irreversible thermal energy source must be present in the inner heliosphere over the heliocentric distance range from 0.15 to 0.47 au. The divergence of the heat flux is positive at heliocentric distances below 0.33 au, while beyond 0.33 au, there is a measurable degradation of the heat flux. Expansion effects dominate the thermal energy budget below 0.3 au. Under our steady-state assumption, the free streaming of the electrons is not sufficient to explain the observed thermal energy density budget. We conjecture that the most likely driver for the required heating process is turbulence. Our results are consistent with the known nonadiabatic polytropic index of the electrons, which we measure as 1.18 in the explored range of heliocentric distances. △ Less

Submitted 20 December, 2022; v1 submitted 3 November, 2022; originally announced November 2022.

Journal ref: ApJ 2022 941 145

Distinct polytropic behavior of plasma during ICME-HSS Interaction

Authors: Kalpesh Ghag, Anil Raghav, Zubair Shaikh, Georgios Nicolaou, Omkar Dhamane, Utsav Panchal

Abstract: Interplanetary Coronal Mass Ejections (ICMEs) and High Speed Streams (HSSs) are noteworthy drivers of disturbance of interplanetary space. Interaction between them can cause several phenomena, such as; generation of waves, enhanced geo-effectiveness, particle acceleration, etc. However, how does thermodynamic properties vary during the ICME-HSS interaction remain an open problem. In this study, we… ▽ More Interplanetary Coronal Mass Ejections (ICMEs) and High Speed Streams (HSSs) are noteworthy drivers of disturbance of interplanetary space. Interaction between them can cause several phenomena, such as; generation of waves, enhanced geo-effectiveness, particle acceleration, etc. However, how does thermodynamic properties vary during the ICME-HSS interaction remain an open problem. In this study, we investigated the polytropic behavior of plasma during an ICME-HSS interaction observed by STEREO and Wind spacecraft. We find that the ICME observed by the STEREO-A has polytropic index $α= 1.0$, i.e., exhibiting isothermal process. Moreover, Wind spacecraft observed the HSS region, non-interacting ICME, and ICME-HSS interaction region. During each regions we found $α$=1.8, $α$=0.7, and $α$=2.5, respectively. It implies that the HSS region exhibits a nearly adiabatic behaviour, ICME region is closely isothermal, and the ICME-HSS interaction region exhibits super-adiabatic behaviour. The insufficient expansion of the ICME due to the interaction with HSS triggers the system for heating and cooling mechanisms which dependent on the degrees of freedom of plasma components. △ Less

Submitted 8 October, 2022; originally announced October 2022.

Energy transport during 3D small-scale reconnection driven by anisotropic plasma turbulence

Authors: Jeffersson A. Agudelo Rueda, Daniel Verscharen, Robert T. Wicks, Christopher J. Owen, Georgios Nicolaou, Kai Germaschewski, Andrew P. Walsh, Ioannis Zouganelis, Santiago Vargas Domínguez

Abstract: Energy dissipation in collisionless plasmas is a longstanding fundamental physics problem. Although it is well known that magnetic reconnection and turbulence are coupled and transport energy from system-size scales to sub-proton scales, the details of the energy distribution and energy dissipation channels remain poorly understood. Especially, the energy transfer and transport associated with thr… ▽ More Energy dissipation in collisionless plasmas is a longstanding fundamental physics problem. Although it is well known that magnetic reconnection and turbulence are coupled and transport energy from system-size scales to sub-proton scales, the details of the energy distribution and energy dissipation channels remain poorly understood. Especially, the energy transfer and transport associated with three dimensional (3D) small-scale reconnection that occurs as a consequence of a turbulent cascade is unknown. We use an explicit fully kinetic particle-in-cell code to simulate 3D small scale magnetic reconnection events forming in anisotropic and Alfvénic decaying turbulence. We identify a highly dynamic and asymmetric reconnection event that involves two reconnecting flux ropes. We use a two-fluid approach based on the Boltzmann equation to study the spatial energy transfer associated with the reconnection event and compare the power density terms in the two-fluid energy equations with standard energy-based damping, heating and dissipation proxies. Our findings suggest that the electron bulk flow transports thermal energy density more efficiently than kinetic energy density. Moreover, in our turbulent reconnection event, the energy-density transfer is dominated by plasma compression. This is consistent with turbulent current sheets and turbulent reconnection events, but not with laminar reconnection. △ Less

Submitted 3 August, 2022; originally announced August 2022.

Comments: Accepted for publication in Apj

Radial evolution of thermal and suprathermal electron populations in the slow solar wind from 0.13 to 0.5 au : Parker Solar Probe Observations

Authors: Joel B. Abraham, Christopher J Owen, Daniel Verscharen, Mayur Bakrania, David Stansby, Robert T. Wicks, Georgios Nicolaou, Phyllis L Whittlesey, Jefferson A. Agudelo Rueda, Seong-Yeop Jeong, Laura Bercic

Abstract: We develop and apply a bespoke fitting routine to a large volume of solar wind electron distribution data measured by Parker Solar Probe (PSP) over its first five orbits, covering radial distances from 0.13 to 0.5 au. We characterise the radial evolution of the electron core, halo and strahl populations in the slow solar wind during these orbits. The fractional densities of these three electron po… ▽ More We develop and apply a bespoke fitting routine to a large volume of solar wind electron distribution data measured by Parker Solar Probe (PSP) over its first five orbits, covering radial distances from 0.13 to 0.5 au. We characterise the radial evolution of the electron core, halo and strahl populations in the slow solar wind during these orbits. The fractional densities of these three electron populations provide evidence for the growth of the combined suprathermal halo and strahl populations from 0.13 to 0.17 au. Moreover, the growth in the halo population is not matched by a decrease of the strahl population at these distances, as has been reported for previous observations at distances greater than 0.3 au. We also find that the halo is negligible at small heliocentric distances. The fractional strahl density remains relatively constant ~1 % below 0.2 au, suggesting that the rise in the relative halo density is not solely due to the transfer of strahl electrons into the halo. △ Less

Submitted 1 June, 2022; v1 submitted 11 April, 2022; originally announced April 2022.

Comments: Published in ApJ

Seasonality, density dependence and spatial population synchrony

Authors: Pedro G Nicolau, Sigrunn H Sørbye, Nigel G Yoccoz, Rolf A Ims

Abstract: Spatial population synchrony has been the focus of theoretical and empirical studies for decades, in the hopes of understanding mechanisms and interactions driving ecological dynamics. In many systems, it is well-known that seasonality plays a critical role in the density dependence structure of the populations, yet this has hardly received any attention in synchrony studies. Here, we propose a pr… ▽ More Spatial population synchrony has been the focus of theoretical and empirical studies for decades, in the hopes of understanding mechanisms and interactions driving ecological dynamics. In many systems, it is well-known that seasonality plays a critical role in the density dependence structure of the populations, yet this has hardly received any attention in synchrony studies. Here, we propose a protocol that allows to elucidate deterministic and stochastic sources of spatial synchrony, while accounting for geographic- and season-specific density dependence. We apply our protocol to seasonally-sampled time series of sub-arctic gray-sided voles, known for marked spatial synchrony. Dissociating seasonal density-dependence contributions to the total observed synchrony reveals differential strength and shape of synchrony patterns by season. Mild winter weather reveals to be an important driver of vole spatial synchrony, with lagged effects in the fall. This has direct implications to the future population dynamics of such species when facing climate change. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: 22 pages, 5 main figures, 2 appendix figures, 3 appendix table

The Stability of the Electron Strahl against the Oblique Fast-magnetosonic/Whistler Instability in the Inner Heliosphere

Authors: Seong-Yeop Jeong, Joel B. Abraham, Daniel Verscharen, Laura Berčič, David Stansby, Georgios Nicolaou, Christopher J. Owen, Robert T. Wicks, Andrew N. Fazakerley, Jeffersson A. Agudelo Rueda, Mayur Bakrania

Abstract: We analyze the micro-kinetic stability of the electron strahl in the solar wind depending on heliocentric distance. The oblique fast-magnetosonic/whistler (FM/W) instability has emerged in the literature as a key candidate mechanism for the effective scattering of the electron strahl into the electron halo population. Using data from Parker Solar Probe (PSP) and Helios, we compare the measured str… ▽ More We analyze the micro-kinetic stability of the electron strahl in the solar wind depending on heliocentric distance. The oblique fast-magnetosonic/whistler (FM/W) instability has emerged in the literature as a key candidate mechanism for the effective scattering of the electron strahl into the electron halo population. Using data from Parker Solar Probe (PSP) and Helios, we compare the measured strahl properties with the analytical thresholds for the oblique FM/W instability in the low- and high-$β_{\parallel c}$ regimes, where $β_{\parallel c}$ is the ratio of the core parallel thermal pressure to the magnetic pressure. Our PSP and Helios data show that the electron strahl is on average stable against the oblique FM/W instability in the inner heliosphere. Our analysis suggests that the instability, if at all, can only be excited sporadically and on short timescales. We discuss the caveats of our analysis and potential alternative explanations for the observed scattering of the electron strahl in the solar wind. Furthermore, we recommend the numerical evaluation of the stability of individual distributions in the future to account for any uncertainties in the validity of the analytical expressions for the instability thresholds. △ Less

Submitted 25 January, 2022; v1 submitted 24 January, 2022; originally announced January 2022.

Comments: Accepted in ApJL

The Kinetic Expansion of Solar-Wind Electrons: Transport Theory and Predictions for the very Inner Heliosphere

Authors: Seong-Yeop Jeong, Daniel Verscharen, Christian Vocks, Joel B. Abraham, Christopher J. Owen, Robert T. Wicks, Andrew N. Fazakerley, David Stansby, Laura Berčič, Georgios Nicolaou, Jeffersson A. Agudelo Rueda, Mayur Bakrania

Abstract: We propose a transport theory for the kinetic evolution of solar-wind electrons in the heliosphere. We derive a gyro-averaged kinetic transport equation that accounts for the spherical expansion of the solar wind and the geometry of the Parker-spiral magnetic field. To solve our three-dimensional kinetic equation, we develop a mathematical approach that combines the Crank--Nicolson scheme in veloc… ▽ More We propose a transport theory for the kinetic evolution of solar-wind electrons in the heliosphere. We derive a gyro-averaged kinetic transport equation that accounts for the spherical expansion of the solar wind and the geometry of the Parker-spiral magnetic field. To solve our three-dimensional kinetic equation, we develop a mathematical approach that combines the Crank--Nicolson scheme in velocity space and a finite-difference Euler scheme in configuration space. We initialize our model with isotropic electron distribution functions and calculate the kinetic expansion at heliocentric distances from 5 to 20 solar radii. In our kinetic model, the electrons evolve mainly through the combination of the ballistic particle streaming, the magnetic mirror force, and the electric field. By applying fits to our numerical results, we quantify the parameters of the electron strahl and core part of the electron velocity distributions. The strahl fit parameters show that the density of the electron strahl is around 7% of the total electron density at a distance of 20 solar radii, the strahl bulk velocity and strahl temperature parallel to the background magnetic field stay approximately constant beyond a distance of 15 solar radii, and $β_{\parallel s}$ (i.e., the ratio between strahl parallel thermal pressure to the magnetic pressure) is approximately constant with heliocentric distance at a value of about 0.02. We compare our results with data measured by Parker Solar Probe. Furthermore, we provide theoretical evidence that the electron strahl is not scattered by the oblique fast-magnetosonic/whistler instability in the near-Sun environment. △ Less

Submitted 3 January, 2022; v1 submitted 31 December, 2021; originally announced January 2022.

Comments: Accepted in ApJ

PySINDy: A comprehensive Python package for robust sparse system identification

Authors: Alan A. Kaptanoglu, Brian M. de Silva, Urban Fasel, Kadierdan Kaheman, Andy J. Goldschmidt, Jared L. Callaham, Charles B. Delahunt, Zachary G. Nicolaou, Kathleen Champion, Jean-Christophe Loiseau, J. Nathan Kutz, Steven L. Brunton

Abstract: Automated data-driven modeling, the process of directly discovering the governing equations of a system from data, is increasingly being used across the scientific community. PySINDy is a Python package that provides tools for applying the sparse identification of nonlinear dynamics (SINDy) approach to data-driven model discovery. In this major update to PySINDy, we implement several advanced feat… ▽ More Automated data-driven modeling, the process of directly discovering the governing equations of a system from data, is increasingly being used across the scientific community. PySINDy is a Python package that provides tools for applying the sparse identification of nonlinear dynamics (SINDy) approach to data-driven model discovery. In this major update to PySINDy, we implement several advanced features that enable the discovery of more general differential equations from noisy and limited data. The library of candidate terms is extended for the identification of actuated systems, partial differential equations (PDEs), and implicit differential equations. Robust formulations, including the integral form of SINDy and ensembling techniques, are also implemented to improve performance for real-world data. Finally, we provide a range of new optimization algorithms, including several sparse regression techniques and algorithms to enforce and promote inequality constraints and stability. Together, these updates enable entirely new SINDy model discovery capabilities that have not been reported in the literature, such as constrained PDE identification and ensembling with different sparse regression optimizers. △ Less

Submitted 25 January, 2022; v1 submitted 12 November, 2021; originally announced November 2021.

Deriving the bulk properties of solar wind electrons observed by Solar Orbiter: A preliminary study of electron plasma thermodynamics

Authors: Georgios Nicolaou, Robert T. Wicks, Christopher J. Owen, Dhiren O. Kataria, Anekallu Chandrasekhar, Gethyn R. Lewis, Daniel Verscharen, Vito Fortunato, Gennaro Mele, Rossana DeMarco, Roberto Bruno

Abstract: We demonstrate the calculation of solar wind electron bulk parameters from recent observations by Solar Wind Analyser Electron Analyser System on board Solar Orbiter. We use our methods to derive the electron bulk parameters in a time interval of a few hours. We attempt a preliminary examination of the polytropic behavior of the electrons by analyzing the derived electron density and temperature.… ▽ More We demonstrate the calculation of solar wind electron bulk parameters from recent observations by Solar Wind Analyser Electron Analyser System on board Solar Orbiter. We use our methods to derive the electron bulk parameters in a time interval of a few hours. We attempt a preliminary examination of the polytropic behavior of the electrons by analyzing the derived electron density and temperature. Moreover, we discuss the challenges in analyzing the observations due to the spacecraft charging and photo-electron contamination in the energy range < 10 eV. Aims: We derive bulk parameters of thermal solar wind electrons by analyzing Solar Orbiter observations and we investigate if there is any typical polytropic model that applies to the electron density and temperature fluctuations. Methods: We use the appropriate transformations to convert the observations to velocity distribution functions in the instrument frame. We then derive the electron bulk parameters by a) calculating the statistical moments of the constructed velocity distribution functions and b) by fitting the constructed distributions with analytical expressions. We firstly test our methods by applying them to an artificial data-set, which we produce by using the forward modeling technique. Results: The forward model validates the analysis techniques which we use to derive the electron bulk parameters. The calculation of the statistical moments and the fitting method determines bulk parameters that are identical within uncertainty to the input parameters we use to simulate the plasma electrons in the first place. An application of our analysis technique to the data reveals a nearly isothermal electron "core". The results are affected by the spacecraft potential and the photo-electron contamination, which we need to characterize in detail in future analyses. △ Less

Submitted 17 September, 2021; originally announced September 2021.

Journal ref: A&A 656, A10 (2021)

Heterogeneity-stabilized homogeneous states in driven media

Authors: Zachary G. Nicolaou, Daniel J. Case, Ernest B. van der Wee, Michelle M. Driscoll, Adilson E. Motter

Abstract: Understanding the relationship between symmetry breaking, system properties, and instabilities has been a problem of longstanding scientific interest. Symmetry-breaking instabilities underlie the formation of important patterns in driven systems, but there are many instances in which such instabilities are undesirable. Using parametric resonance as a model process, here we show that a range of sta… ▽ More Understanding the relationship between symmetry breaking, system properties, and instabilities has been a problem of longstanding scientific interest. Symmetry-breaking instabilities underlie the formation of important patterns in driven systems, but there are many instances in which such instabilities are undesirable. Using parametric resonance as a model process, here we show that a range of states that would be destabilized by symmetry-breaking instabilities can be preserved and stabilized by the introduction of suitable system asymmetry. Because symmetric states are spatially homogeneous and asymmetric systems are spatially heterogeneous, we refer to this effect as heterogeneity-stabilized homogeneity. We illustrate this effect theoretically using driven pendulum array models and demonstrate it experimentally using Faraday wave instabilities. Our results have potential implications for the mitigation of instabilities in engineered systems and the emergence of homogeneous states in natural systems with inherent heterogeneities. △ Less

Submitted 2 August, 2021; originally announced August 2021.

Comments: 9 pages, 4 figures, plus Supplementary Information

Journal ref: Nat. Commun. 12, 4486 (2021)

Whistler instability driven by the sunward electron deficit in the solar wind

Authors: Laura Berčič, Daniel Verscharen, Christopher J. Owen, Lucas Colomban, Matthieu Kretzschmar, Thomas Chust, Milan Maksimović, Dhiren Kataria, Etienne Behar, Matthieu Berthomier, Roberto Bruno, Vito Fortunato, Christopher W. Kelly, Yuri. V. Khotyaintsev, Gethyn R. Lewis, Stefano Livi, Philippe Louarn, Gennaro Mele, Georgios Nicolaou, Gillian Watson, Robert T. Wicks

Abstract: Solar wind electrons play an important role in the energy balance of the solar wind acceleration by carrying energy into interplanetary space in the form of electron heat flux. The heat flux is stored in the complex electron velocity distribution functions (VDFs) shaped by expansion, Coulomb collisions, and field-particle interactions. We investigate how the suprathermal electron deficit in the an… ▽ More Solar wind electrons play an important role in the energy balance of the solar wind acceleration by carrying energy into interplanetary space in the form of electron heat flux. The heat flux is stored in the complex electron velocity distribution functions (VDFs) shaped by expansion, Coulomb collisions, and field-particle interactions. We investigate how the suprathermal electron deficit in the anti-strahl direction, which was recently discovered in the near-Sun solar wind, drives a kinetic instability and creates whistler waves with wave vectors that are quasi-parallel to the direction of the background magnetic field. We combine high-cadence measurements of electron pitch-angle distribution functions and electromagnetic waves provided by Solar Orbiter during its first orbit. Our case study is based on a burst-mode data interval from the Electrostatic Analyser System (SWA-EAS) at a distance of 112 $R_S$ (0.52 au) from the Sun, during which several whistler wave packets were detected by Solar Orbiter's Radio and Plasma Waves (RPW) instrument. The sunward deficit creates kinetic conditions under which the quasi-parallel whistler wave is driven unstable. We directly test our predictions for the existence of these waves through solar wind observations. We find whistler waves that are quasi-parallel and almost circularly polarised, propagating away from the Sun, coinciding with a pronounced sunward deficit in the electron VDF. The cyclotron-resonance condition is fulfilled for electrons moving in the direction opposite to the direction of wave propagation, with energies corresponding to those associated with the sunward deficit. △ Less

Submitted 22 July, 2021; originally announced July 2021.

Journal ref: A&A 656, A31 (2021)

Anharmonic classical time crystals: A coresonance pattern formation mechanism

Authors: Zachary G. Nicolaou, Adilson E. Motter

Abstract: Driven many-body systems have been shown to exhibit discrete time crystal phases characterized by broken discrete time-translational symmetry. This has been achieved generally through a subharmonic response, in which the system undergoes one oscillation every other driving period. Here, we demonstrate that classical time crystals do not need to resonate in a subharmonic fashion but instead can als… ▽ More Driven many-body systems have been shown to exhibit discrete time crystal phases characterized by broken discrete time-translational symmetry. This has been achieved generally through a subharmonic response, in which the system undergoes one oscillation every other driving period. Here, we demonstrate that classical time crystals do not need to resonate in a subharmonic fashion but instead can also exhibit a continuously tunable anharmonic response to driving, which we show can emerge through a coresonance between modes in different branches of the dispersion relation in a parametrically driven medium. This response, characterized by a typically incommensurate ratio between the resonant frequencies and the driving frequency, is demonstrated by introducing a time crystal model consisting of an array of coupled pendula with alternating lengths. Importantly, the coresonance mechanism is the result of a bifurcation involving a fixed point and an invariant torus, with no intermediate limit cycles. This bifurcation thus gives rise to a many-body symmetry-breaking phenomenon directly connecting the symmetry-unbroken phase with a previously uncharacterized phase of matter, which we call an anharmonic time crystal phase. The mechanism is shown to generalize to driven media with any number of coupled fields and is expected to give rise to anharmonic responses in a range of weakly damped pattern-forming systems, with potential applications to the study of nonequilibrium phases, frequency conversion, and acoustic cloaking. △ Less

Submitted 11 May, 2021; originally announced May 2021.

Comments: 11 pages, 7 figures

Journal ref: Phys. Rev. Research 3, 023106 (2021)

First-year ion-acoustic wave observations in the solar wind by the RPW/TDS instrument onboard Solar Orbiter

Authors: D. Píša, J. Souček, O. Santolík, M. Hanzelka, G. Nicolaou, M. Maksimovic, S. D. Bale, T. Chust, Y. Khotyaintsev, V. Krasnoselskikh, M. Kretzschmar, E. Lorfèvre, D. Plettemeier, M. Steller, Š. Štverák, P. Trávníček, A. Vaivads, A. Vecchio, T. Horbury, H. O'Brien, V. Evans, V. Angelini, C. J. Owen, P. Louarn

Abstract: Electric field measurements of the Time Domain Sampler (TDS) receiver, part of the Radio and Plasma Waves (RPW) instrument on board Solar Orbiter, often exhibit very intense broadband wave emissions at frequencies below 20~kHz in the spacecraft frame. In this paper, we present a year-long study of electrostatic fluctuations observed in the solar wind at an interval of heliocentric distances from 0… ▽ More Electric field measurements of the Time Domain Sampler (TDS) receiver, part of the Radio and Plasma Waves (RPW) instrument on board Solar Orbiter, often exhibit very intense broadband wave emissions at frequencies below 20~kHz in the spacecraft frame. In this paper, we present a year-long study of electrostatic fluctuations observed in the solar wind at an interval of heliocentric distances from 0.5 to 1~AU. The RPW/TDS observations provide a nearly continuous data set for a statistical study of intense waves below the local plasma frequency. The on-board and continuously collected and processed properties of waveform snapshots allow for the mapping plasma waves at frequencies between 200~Hz and 20~kHz. We used the triggered waveform snapshots and a Doppler-shifted solution of the dispersion relation for wave mode identification in order to carry out a detailed spectral and polarization analysis. Electrostatic ion-acoustic waves are the common wave emissions observed between the local electron and proton plasma frequency in the soler wind. The occurrence rate of ion-acoustic waves peaks around perihelion at distances of 0.5~AU and decreases with increasing distances, with only a few waves detected per day at 0.9~AU. Waves are more likely to be observed when the local proton moments and magnetic field are highly variable. A more detailed analysis of more than 10000 triggered waveform snapshots shows the mean wave frequency at about 3 kHz and wave amplitude about 2.5 mV/m. The wave amplitude varies as 1/R^(1.38) with the heliocentric distance. The relative phase distribution between two components of the E-field shows a mostly linear wave polarization. Electric field fluctuations are closely aligned with the directions of the ambient field lines. Only a small number (3%) of ion-acoustic waves are observed at larger magnetic discontinuities. △ Less

Submitted 6 September, 2021; v1 submitted 7 April, 2021; originally announced April 2021.

Comments: 9 pages, 8 figures

Journal ref: A&A 656, A14 (2021)

Three-dimensional magnetic reconnection in particle-in-cell simulations of anisotropic plasma turbulence

Authors: Jeffersson A. Agudelo Rueda, Daniel Verscharen, Robert T. Wicks, Christopher J. Owen, Georgios Nicolaou, Andrew P. Walsh, Ioannis Zouganelis, Kai Germaschewski, Santiago Vargas Domínguez

Abstract: We use 3D fully kinetic particle-in-cell simulations to study the occurrence of magnetic reconnection in a simulation of decaying turbulence created by anisotropic counter-propagating low-frequency Alfvén waves consistent with critical-balance theory. We observe the formation of small-scale current-density structures such as current filaments and current sheets as well as the formation of magnetic… ▽ More We use 3D fully kinetic particle-in-cell simulations to study the occurrence of magnetic reconnection in a simulation of decaying turbulence created by anisotropic counter-propagating low-frequency Alfvén waves consistent with critical-balance theory. We observe the formation of small-scale current-density structures such as current filaments and current sheets as well as the formation of magnetic flux ropes as part of the turbulent cascade. The large magnetic structures present in the simulation domain retain the initial anisotropy while the small-scale structures produced by the turbulent cascade are less anisotropic. To quantify the occurrence of reconnection in our simulation domain, we develop a new set of indicators based on intensity thresholds to identify reconnection events in which both ions and electrons are heated and accelerated in 3D particle-in-cell simulations. According to the application of these indicators, we identify the occurrence of reconnection events in the simulation domain and analyse one of these events in detail. The event is related to the reconnection of two flux ropes, and the associated ion and electron exhausts exhibit a complex three-dimensional structure. We study the profiles of plasma and magnetic-field fluctuations recorded along artificial-spacecraft trajectories passing near and through the reconnection region. Our results suggest the presence of particle heating and acceleration related to small-scale reconnection events within magnetic flux ropes produced by the anisotropic Alfvénic turbulent cascade in the solar wind. These events are related to current structures of order a few ion inertial lengths in size. △ Less

Submitted 24 March, 2021; originally announced March 2021.

Comments: Accepted for publication in J. Plasma Phys

Matching temporal signatures of solar features to their corresponding solar wind outflows

Authors: Diego de Pablos, David M. Long, Christopher J. Owen, Gherardo Valori, Georgios Nicolaou, Louise K. Harra

Abstract: The role of small-scale coronal eruptive phenomena in the generation and heating of the solar wind remains an open question. Here, we investigate the role played by coronal jets in forming the solar wind by testing whether temporal variations associated with jetting in EUV intensity can be identified in the outflowing solar wind plasma. This type of comparison is challenging due to inherent differ… ▽ More The role of small-scale coronal eruptive phenomena in the generation and heating of the solar wind remains an open question. Here, we investigate the role played by coronal jets in forming the solar wind by testing whether temporal variations associated with jetting in EUV intensity can be identified in the outflowing solar wind plasma. This type of comparison is challenging due to inherent differences between remote-sensing observations of the source and in situ observations of the outflowing plasma, as well as travel time and evolution of the solar wind throughout the heliosphere. To overcome these, we propose a novel algorithm combining signal filtering, two-step solar wind ballistic backmapping, window shifting, and Empirical Mode Decomposition. We first validate the method using synthetic data, before applying it to measurements from the Solar Dynamics Observatory, and Wind spacecraft. The algorithm enables the direct comparison of remote sensing observations of eruptive phenomena in the corona to in situ measurements of solar wind parameters, among other potential uses. After application to these datasets, we find several time windows where signatures of dynamics found in the corona are embedded in the solar wind stream, at a time significantly earlier than expected from simple ballistic backmapping, with the best performing in situ parameter being the solar wind mass flux. △ Less

Submitted 16 March, 2021; originally announced March 2021.

Relationship between polytropic index and temperature anisotropy in space plasmas

Authors: George Livadiotis, George Nicolaou

Abstract: The paper develops a theoretical relationship between the polytropic index and the temperature anisotropy that may characterize space plasmas. The derivation is based on the correlation among the kinetic energies of particles with velocities described by anisotropic kappa distributions. The correlation coefficient depends on the effective dimensionality of the velocity distribution function, havin… ▽ More The paper develops a theoretical relationship between the polytropic index and the temperature anisotropy that may characterize space plasmas. The derivation is based on the correlation among the kinetic energies of particles with velocities described by anisotropic kappa distributions. The correlation coefficient depends on the effective dimensionality of the velocity distribution function, having its shape determined by the temperature anisotropies caused by the ambient magnetic field; on the other hand, the effective dimensionality is directly dependent on the polytropic index. This analysis leads to the connection between the correlation coefficient, effective dimensionality of velocity space, and the polytropic index, with the ratio of temperature anisotropy. Moreover, a data and statistical analysis is performed to test the application of the theoretically developed model, in the solar wind proton plasma near 1 AU. The derived theoretical relationship is in good agreement with observations, showing that the lowest and classical value of the adiabatic polytropic index occurs in the isotropic case, while higher indices characterize anisotropic plasmas. Finally, possible extensions of the theory considering (i) non-adiabatic polytropic behavior, and (ii) more general distributions, are further discussed. △ Less

Submitted 24 December, 2020; originally announced December 2020.

Comments: 23 pages, 8 figures

Anisotropic kappa distributions I: Formulation based on particle correlations

Authors: George Livadiotis, George Nicolaou, Frederic Allegrini

Abstract: We develop the theoretical basis for the connection of the variety of anisotropic distributions with the statistical correlations among particles velocity components. By examining the most common anisotropic distribution, we derive the correlation coefficient among particle energies, show how this correlation is connected to the effective dimensionality of the velocity distribution, and derive the… ▽ More We develop the theoretical basis for the connection of the variety of anisotropic distributions with the statistical correlations among particles velocity components. By examining the most common anisotropic distribution, we derive the correlation coefficient among particle energies, show how this correlation is connected to the effective dimensionality of the velocity distribution, and derive the connection between anisotropy and adiabatic polytropic index. Having established the importance of correlation among particles in the formulation of anisotropic kappa distributions, we generalize these distributions within the framework of nonextensive statistical mechanics and based on the types of homogeneous or heterogeneous correlations among the particles velocity components. The formulation of the developed generalized distributions mediates the main two types of anisotropic kappa distributions, considering (a) equal correlations, or (b) zero correlations, among different velocity components. Finally, the developed anisotropic kappa distributions are expressed in terms of the energy and pitch angle in arbitrary reference frames. △ Less

Submitted 15 December, 2020; originally announced December 2020.

Comments: 46 pages, 15 figures

Design and Optimization of a High-Time-Resolution Magnetic Plasma Analyzer (MPA)

Authors: Benjamin Criton, Georgios Nicolaou, Daniel Verscharen

Abstract: In-situ measurements of space plasma throughout the solar system require high time resolution to understand the plasma's kinetic fine structure and evolution. In this context, research~is conducted to design instruments with the capability to acquire the plasma velocity distribution and its moments with high cadence. We study a new instrument design, using a constant magnetic field generated by tw… ▽ More In-situ measurements of space plasma throughout the solar system require high time resolution to understand the plasma's kinetic fine structure and evolution. In this context, research~is conducted to design instruments with the capability to acquire the plasma velocity distribution and its moments with high cadence. We study a new instrument design, using a constant magnetic field generated by two permanent magnets, to analyze solar wind protons and $α$-particles with high time resolution. We determine the optimal configuration of the instrument in terms of aperture size, sensor position, pixel size and magnetic field strength. We conduct this analysis based on analytical calculations and SIMION simulations of the particle trajectories in our instrument. We evaluate the velocity resolution of the instrument as well as Poisson errors associated with finite counting statistics. Our instrument is able to resolve Maxwellian and $κ$-distributions for both protons and $α$-particles. This method retrieves measurements of the moments (density, bulk speed and temperature) with a relative error below 1%. Our instrument design achieves these results with an acquisition time of only 5 ms, significantly faster than state-of-the-art electrostatic analyzers. Although the instrument only acquires one-dimensional cuts of the distribution function in velocity space, the simplicity and reliability of the presented instrument concept are two key advantages of our new design. △ Less

Submitted 27 November, 2020; originally announced November 2020.

Journal ref: Appl. Sci. 2020, 10(23), 8483

Enhanced proton parallel temperature inside patches of switchbacks in the inner heliosphere

Authors: L. D. Woodham, T. S. Horbury, L. Matteini, T. Woolley, R. Laker, S. D. Bale, G. Nicolaou, J. E. Stawarz, D. Stansby, H. Hietala, D. E. Larson, R. Livi, J. L. Verniero, M. McManus, J. C. Kasper, K. E. Korreck, N. Raouafi, M. Moncuquet, M. P. Pulupa

Abstract: 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 anis… ▽ More 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. △ Less

Submitted 20 October, 2020; originally announced October 2020.

Comments: Letter accepted as part of the PSP Special Issue in A&A

Journal ref: A&A 650, L1 (2021)

Model-based bias correction for short AR(1) and AR(2) processes

Authors: Sigrunn H. Sørbye, Pedro G. Nicolau, Håvard Rue

Abstract: The class of autoregressive (AR) processes is extensively used to model temporal dependence in observed time series. Such models are easily available and routinely fitted using freely available statistical software like R. A potential caveat in analyzing short time series is that commonly applied estimators for the coefficients of AR processes are severely biased. This paper suggests a model-based… ▽ More The class of autoregressive (AR) processes is extensively used to model temporal dependence in observed time series. Such models are easily available and routinely fitted using freely available statistical software like R. A potential caveat in analyzing short time series is that commonly applied estimators for the coefficients of AR processes are severely biased. This paper suggests a model-based approach for bias correction of well-known estimators for the coefficients of first and second-order stationary AR processes, taking the sampling distribution of the original estimator into account. This is achieved by modeling the relationship between the true and estimated AR coefficients using weighted orthogonal polynomial regression, fitted to a huge number of simulations. The finite-sample distributions of the new estimators are approximated using transformations of skew-normal densities and their properties are demonstrated by simulations and in the analysis of a real ecological data set. The new estimators are easily available in our accompanying R-package ARbiascorrect for time series of length n = 10, 11, ... , 50, where original estimates are found using exact or conditional maximum likelihood, Burg's method or the Yule-Walker equations. △ Less

Submitted 12 October, 2020; originally announced October 2020.

Comments: 21 pages, 11 figures

MSC Class: 62M10; 62F10; 62E17

The Solar Orbiter Science Activity Plan: translating solar and heliospheric physics questions into action

Authors: I. Zouganelis, A. De Groof, A. P. Walsh, D. R. Williams, D. Mueller, O. C. St Cyr, F. Auchere, D. Berghmans, A. Fludra, T. S. Horbury, R. A. Howard, S. Krucker, M. Maksimovic, C. J. Owen, J. Rodriiguez-Pacheco, M. Romoli, S. K. Solanki, C. Watson, L. Sanchez, J. Lefort, P. Osuna, H. R. Gilbert, T. Nieves-Chinchilla, L. Abbo, O. Alexandrova , et al. (160 additional authors not shown)

Abstract: Solar Orbiter is the first space mission observing the solar plasma both in situ and remotely, from a close distance, in and out of the ecliptic. The ultimate goal is to understand how the Sun produces and controls the heliosphere, filling the Solar System and driving the planetary environments. With six remote-sensing and four in-situ instrument suites, the coordination and planning of the operat… ▽ More Solar Orbiter is the first space mission observing the solar plasma both in situ and remotely, from a close distance, in and out of the ecliptic. The ultimate goal is to understand how the Sun produces and controls the heliosphere, filling the Solar System and driving the planetary environments. With six remote-sensing and four in-situ instrument suites, the coordination and planning of the operations are essential to address the following four top-level science questions: (1) What drives the solar wind and where does the coronal magnetic field originate? (2) How do solar transients drive heliospheric variability? (3) How do solar eruptions produce energetic particle radiation that fills the heliosphere? (4) How does the solar dynamo work and drive connections between the Sun and the heliosphere? Maximising the mission's science return requires considering the characteristics of each orbit, including the relative position of the spacecraft to Earth (affecting downlink rates), trajectory events (such as gravitational assist manoeuvres), and the phase of the solar activity cycle. Furthermore, since each orbit's science telemetry will be downloaded over the course of the following orbit, science operations must be planned at mission level, rather than at the level of individual orbits. It is important to explore the way in which those science questions are translated into an actual plan of observations that fits into the mission, thus ensuring that no opportunities are missed. First, the overarching goals are broken down into specific, answerable questions along with the required observations and the so-called Science Activity Plan (SAP) is developed to achieve this. The SAP groups objectives that require similar observations into Solar Orbiter Observing Plans (SOOPs), resulting in a strategic, top-level view of the optimal opportunities for science observations during the mission lifetime. △ Less

Submitted 22 September, 2020; originally announced September 2020.

Comments: 20 pages, 1 figure, accepted by Astronomy & Astrophysics

Journal ref: A&A 642, A3 (2020)

Missing links as a source of seemingly variable constants in complex reaction networks

Authors: Zachary G. Nicolaou, Adilson E. Motter

Abstract: A major challenge in network science is to determine parameters governing complex network dynamics from experimental observations and theoretical models. In complex chemical reaction networks, for example, such as those describing processes in internal combustion engines and power generators, rate constant estimates vary significantly across studies despite substantial experimental efforts. Here,… ▽ More A major challenge in network science is to determine parameters governing complex network dynamics from experimental observations and theoretical models. In complex chemical reaction networks, for example, such as those describing processes in internal combustion engines and power generators, rate constant estimates vary significantly across studies despite substantial experimental efforts. Here, we examine the possibility that variability in measured constants can be largely attributed to the impact of missing network information on parameter estimation. Through the numerical simulation of measurements in incomplete chemical reaction networks, we show that unaccountability of network links presumed unimportant (with local sensitivity amounting to less than two percent of that of a measured link) can create apparent rate constant variations as large as one order of magnitude even if no experimental errors are present in the data. Furthermore, the correlation coefficient between the logarithmic deviation of the rate constant estimate and the cumulative relative sensitivity of the neglected reactions was less than $0.5$ in all cases. Thus, for dynamical processes on complex networks, iteratively expanding a model by determining new parameters from data collected under specific conditions is unlikely to produce reliable results. △ Less

Submitted 28 October, 2020; v1 submitted 16 September, 2020; originally announced September 2020.

Comments: 12 pages, 10 figures. Updated after publication

Journal ref: Phys. Rev. Research 2, 043135 (2020)

Coherent Dynamics Enhanced by Uncorrelated Noise

Authors: Zachary G. Nicolaou, Michael Sebek, Istvan Z. Kiss, Adilson E. Motter

Abstract: Synchronization is a widespread phenomenon observed in physical, biological, and social networks, which persists even under the influence of strong noise. Previous research on oscillators subject to common noise has shown that noise can actually facilitate synchronization, as correlations in the dynamics can be inherited from the noise itself. However, in many spatially distributed networks, such… ▽ More Synchronization is a widespread phenomenon observed in physical, biological, and social networks, which persists even under the influence of strong noise. Previous research on oscillators subject to common noise has shown that noise can actually facilitate synchronization, as correlations in the dynamics can be inherited from the noise itself. However, in many spatially distributed networks, such as the mammalian circadian system, the noise that different oscillators experience can be effectively uncorrelated. Here, we show that uncorrelated noise can in fact enhance synchronization when the oscillators are coupled. Strikingly, our analysis also shows that uncorrelated noise can be more effective than common noise in enhancing synchronization. We first establish these results theoretically for phase and phase-amplitude oscillators subject to either or both additive and multiplicative noise. We then confirm the predictions through experiments on coupled electrochemical oscillators. Our findings suggest that uncorrelated noise can promote rather than inhibit coherence in natural systems and that the same effect can be harnessed in engineered systems. △ Less

Submitted 26 May, 2021; v1 submitted 24 August, 2020; originally announced August 2020.

Comments: 6 pages, 4 figures, plus supplement

Journal ref: Phys. Rev. Lett. 125, 094101 (2020)

Non-normality and non-monotonic dynamics in complex reaction networks

Authors: Zachary G. Nicolaou, Takashi Nishikawa, Schuyler B. Nicholson, Jason R. Green, Adilson E. Motter

Abstract: Complex chemical reaction networks, which underlie many industrial and biological processes, often exhibit non-monotonic changes in chemical species concentrations, typically described using nonlinear models. Such non-monotonic dynamics are in principle possible even in linear models if the matrices defining the models are non-normal, as characterized by a necessarily non-orthogonal set of eigenve… ▽ More Complex chemical reaction networks, which underlie many industrial and biological processes, often exhibit non-monotonic changes in chemical species concentrations, typically described using nonlinear models. Such non-monotonic dynamics are in principle possible even in linear models if the matrices defining the models are non-normal, as characterized by a necessarily non-orthogonal set of eigenvectors. However, the extent to which non-normality is responsible for non-monotonic behavior remains an open question. Here, using a master equation to model the reaction dynamics, we derive a general condition for observing non-monotonic dynamics of individual species, establishing that non-normality promotes non-monotonicity but is not a requirement for it. In contrast, we show that non-normality is a requirement for non-monotonic dynamics to be observed in the Rényi entropy. Using hydrogen combustion as an example application, we demonstrate that non-monotonic dynamics under experimental conditions are supported by a linear chain of connected components, in contrast with the dominance of a single giant component observed in typical random reaction networks. The exact linearity of the master equation enables development of rigorous theory and simulations for dynamical networks of unprecedented size (approaching $10^5$ dynamical variables, even for a network of only 20 reactions and involving less than 100 atoms). Our conclusions are expected to hold for other combustion processes, and the general theory we develop is applicable to all chemical reaction networks, including biological ones. △ Less

Submitted 21 August, 2020; originally announced August 2020.

Comments: Software implementing our methods is available as a Github repository at https://github.com/znicolaou/ratematrix and an animated version of Fig. 1 is available at https://northwestern.box.com/s/otn3m2cov9gi3enht3r8jh5kjo9qnv6d

Journal ref: Phys. Rev. Research 2, 043059 (2020)

Polytropic Behavior of Solar Wind Protons Observed by Parker Solar Probe

Authors: Georgios Nicolaou, George Livadiotis, Robert T. Wicks, Daniel Verscharen, Bennett A. Maruca

Abstract: A polytropic process describes the transition of a fluid from one state to another through a specific relationship between the fluid density and temperature. The value of the polytropic index that governs this relationship determines the heat transfer and the effective degrees of freedom during the process. In this study, we analyze solar wind proton plasma measurements, obtained by the Faraday cu… ▽ More A polytropic process describes the transition of a fluid from one state to another through a specific relationship between the fluid density and temperature. The value of the polytropic index that governs this relationship determines the heat transfer and the effective degrees of freedom during the process. In this study, we analyze solar wind proton plasma measurements, obtained by the Faraday cup instrument on-board Parker Solar Probe. We examine the large-scale variations of the proton plasma density and temperature within the inner heliosphere explored by the spacecraft. We also address a polytropic behavior in the density and temperature fluctuations in short-time intervals, which we analyze in order to derive the effective polytropic index of small time-scale processes. The large-scale variations of the solar wind proton density and temperature which are associated with the plasma expansion through the heliosphere, follow a polytropic model with a polytropic index ~5/3. On the other hand, the short time-scale fluctuations which may be associated with turbulence, follow a model with a larger polytropic index. We investigate possible correlations between the polytropic index of short time-scale fluctuations and the plasma speed, plasma beta, and the magnetic field direction. We discuss the scenario of mechanisms including energy transfer or mechanisms that restrict the particle effective degrees of freedom. △ Less

Submitted 19 May, 2020; originally announced May 2020.

Comments: 20 pages, 9 figures

Critical Switching in Globally Attractive Chimeras

Authors: Yuanzhao Zhang, Zachary G. Nicolaou, Joseph D. Hart, Rajarshi Roy, Adilson E. Motter

Abstract: We report on a new type of chimera state that attracts almost all initial conditions and exhibits power-law switching behavior in networks of coupled oscillators. Such switching chimeras consist of two symmetric configurations, which we refer to as subchimeras, in which one cluster is synchronized and the other is incoherent. Despite each subchimera being linearly stable, switching chimeras are ex… ▽ More We report on a new type of chimera state that attracts almost all initial conditions and exhibits power-law switching behavior in networks of coupled oscillators. Such switching chimeras consist of two symmetric configurations, which we refer to as subchimeras, in which one cluster is synchronized and the other is incoherent. Despite each subchimera being linearly stable, switching chimeras are extremely sensitive to noise: arbitrarily small noise triggers and sustains persistent switching between the two symmetric subchimeras. The average switching frequency scales as a power law with the noise intensity, which is in contrast with the exponential scaling observed in typical stochastic transitions. Rigorous numerical analysis reveals that the power-law switching behavior originates from intermingled basins of attraction associated with the two subchimeras, which in turn are induced by chaos and symmetry in the system. The theoretical results are supported by experiments on coupled optoelectronic oscillators, which demonstrate the generality and robustness of switching chimeras. △ Less

Submitted 24 February, 2020; v1 submitted 18 November, 2019; originally announced November 2019.

Journal ref: Phys. Rev. X 10, 011044 (2020)

The Impact of Turbulent Solar Wind Fluctuations on Solar Orbiter Plasma Proton Measurements

Authors: Georgios Nicolaou, Daniel Verscharen, Robert T. Wicks, Christopher J. Owen

Abstract: Solar Orbiter will observe the Sun and the inner heliosphere to study the connections between solar activity, coronal structure, and the origin of the solar wind. The plasma instruments on board Solar Orbiter will determine the three-dimensional velocity distribution functions of the plasma ions and electrons with high time resolution. The analysis of these distributions will determine the plasma… ▽ More Solar Orbiter will observe the Sun and the inner heliosphere to study the connections between solar activity, coronal structure, and the origin of the solar wind. The plasma instruments on board Solar Orbiter will determine the three-dimensional velocity distribution functions of the plasma ions and electrons with high time resolution. The analysis of these distributions will determine the plasma bulk parameters, such as density, velocity, and temperature. This paper examines the effects of short-time-scale plasma variations on particle measurements and the estimated bulk parameters of the plasma. For the purpose of this study, we simulate the expected observations of solar wind protons, taking into account the performance of the Proton-Alpha Sensor (PAS) on board Solar Orbiter. We particularly examine the effects of Alfvénic and slow-mode-like fluctuations, commonly observed in the solar wind on timescales of milliseconds to hours, on the observations. We do this by constructing distribution functions from modeled observations and calculate their statistical moments in order to derive plasma bulk parameters. The comparison between the derived parameters with the known input allows us to estimate the expected accuracy of Solar Orbiter proton measurements in the solar wind under typical conditions. We find that the plasma fluctuations due to these turbulence effects have only minor effects on future SWA-PAS observations △ Less

Submitted 25 October, 2019; originally announced October 2019.

Minimal scattering entanglement in one-dimensional trapped gases

Authors: Zachary G. Nicolaou, Bohan Xu, Adilson E. Motter

Abstract: The prospect of controlling entanglement in interacting quantum systems offers a myriad of technological and scientific promises, given the progress in experimental studies in systems such as ultracold trapped gases. This control is often challenging because of decoherence, the process by which environmental interactions create spurious entanglements that can destroy the desired entanglement. Cons… ▽ More The prospect of controlling entanglement in interacting quantum systems offers a myriad of technological and scientific promises, given the progress in experimental studies in systems such as ultracold trapped gases. This control is often challenging because of decoherence, the process by which environmental interactions create spurious entanglements that can destroy the desired entanglement. Considering the collisional decoherence that is relevant for quantum measurements utilizing scattering in one-dimensional trapped gases, here we derive a relationship between particle masses and wave packet widths that minimizes the entanglement created during scattering. We assess the relevance of our results by directly observing this relationship in the emergent scales of a master equation for a particle undergoing nonthermal scattering. Our relationship is independent of the details of the particle interactions and sheds light on how to design scattering processes that minimize decoherence. △ Less

Submitted 14 January, 2019; originally announced January 2019.

Comments: 8 pages, 4 figures. Supplemental material available on the publisher website: https://journals.aps.org/pra/abstract/10.1103/PhysRevA.99.012316#supplemental

Journal ref: Phys. Rev. A 99 (1), 012316 (2019)

Multifaceted dynamics of Janus oscillator networks

Authors: Zachary G. Nicolaou, Deniz Eroglu, Adilson E. Motter

Abstract: Recent research has led to the discovery of fundamental new phenomena in network synchronization, including chimera states, explosive synchronization, and asymmetry-induced synchronization. Each of these phenomena has thus far been observed only in systems designed to exhibit that one phenomenon, which raises the questions of whether they are mutually compatible and, if so, under what conditions t… ▽ More Recent research has led to the discovery of fundamental new phenomena in network synchronization, including chimera states, explosive synchronization, and asymmetry-induced synchronization. Each of these phenomena has thus far been observed only in systems designed to exhibit that one phenomenon, which raises the questions of whether they are mutually compatible and, if so, under what conditions they co-occur. Here, we introduce a class of remarkably simple oscillator networks that concurrently exhibit all of these phenomena. The dynamical units consist of pairs of nonidentical phase oscillators, which we refer to as Janus oscillators by analogy with Janus particles and the mythological figure from which their name is derived. In contrast to previous studies, these networks exhibit (i) explosive synchronization with identical oscillators; (ii) extreme multistability of chimera states, including traveling, intermittent, and bouncing chimeras; and (iii) asymmetry-induced synchronization in which synchronization is promoted by random oscillator heterogeneity. These networks also exhibit the previously unobserved possibility of inverted synchronization transitions, in which a transition to a more synchronous state is induced by a reduction rather than an increase in the coupling strength. These various phenomena are shown to emerge under rather parsimonious conditions, and even in locally connected ring topologies, which has the potential to facilitate their use to control and manipulate synchronization in experiments. △ Less

Submitted 30 January, 2019; v1 submitted 15 October, 2018; originally announced October 2018.

Comments: 12 pages, 11 figures, Supplemental animation at article webpage: http://link.aps.org/supplemental/10.1103/PhysRevX.9.011017

Journal ref: Phys. Rev. X 9, 011017 (2019)

The Fluid-like Behavior of Kinetic Alfvén Turbulence in Space Plasma

Authors: Honghong Wu, Daniel Verscharen, Robert T. Wicks, Christopher H. Chen, Jiansen He, Georgios Nicolaou

Abstract: Kinetic Alfvén waves (KAWs) are the short-wavelength extension of the MHD Alfvén-wave branch in the case of highly-oblique propagation with respect to the background magnetic field. Observations of space plasma show that small-scale turbulence is mainly KAW-like. We apply two theoretical approaches, collisional two-fluid theory and collisionless kinetic theory, to obtain predictions for the KAW po… ▽ More Kinetic Alfvén waves (KAWs) are the short-wavelength extension of the MHD Alfvén-wave branch in the case of highly-oblique propagation with respect to the background magnetic field. Observations of space plasma show that small-scale turbulence is mainly KAW-like. We apply two theoretical approaches, collisional two-fluid theory and collisionless kinetic theory, to obtain predictions for the KAW polarizations depending on $β_\mathrm{p}$ (the ratio of the proton thermal pressure to the magnetic pressure) at the ion gyroscale in terms of fluctuations in density, bulk velocity, and pressure. We perform a wavelet analysis of MMS magnetosheath measurements and compare the observations with both theories. We find that the two-fluid theory predicts the observations better than kinetic theory, suggesting that the small-scale KAW-like fluctuations exhibit a fluid-like behavior in the magnetosheath although the plasma is weakly collisional. We also present predictions for the KAW polarizations in the inner heliosphere that are testable with Parker Solar Probe and Solar Orbiter. △ Less

Submitted 15 January, 2019; v1 submitted 29 August, 2018; originally announced August 2018.

Comments: 9 pages, 6 figures, submitted

Journal ref: The Astrophysical Journal, 870:106 (7pp), 2019 January 10

The root of a comet tail - Rosetta ion observations at comet 67P/Churyumov--Gerasimenko

Authors: E. Behar, H. Nilsson, P. Henri, L. Bercic, G. Nicolaou, G. Stenberg Wieser, M. Wieser, B. Tabone, M. Saillenfest, C. Goetz

Abstract: The first 1000 km of the ion tail of comet 67P/Churyumov-Gerasimenko were explored by the European Rosetta spacecraft, 2.7 au away from the Sun. We characterised the dynamics of both the solar wind and the cometary ions on the night-side of the comet's atmosphere. We analysed in situ ion and magnetic field measurements and compared the data to a semi-analytical model. The cometary ions are observe… ▽ More The first 1000 km of the ion tail of comet 67P/Churyumov-Gerasimenko were explored by the European Rosetta spacecraft, 2.7 au away from the Sun. We characterised the dynamics of both the solar wind and the cometary ions on the night-side of the comet's atmosphere. We analysed in situ ion and magnetic field measurements and compared the data to a semi-analytical model. The cometary ions are observed flowing close to radially away from the nucleus during the entire excursion. The solar wind is deflected by its interaction with the new-born cometary ions. Two concentric regions appear, an inner region dominated by the expanding cometary ions and an outer region dominated by the solar wind particles. The single night-side excursion operated by Rosetta revealed that the near radial flow of the cometary ions can be explained by the combined action of three different electric field components, resulting from the ion motion, the electron pressure gradients, and the magnetic field draping. The observed solar wind deflection is governed mostly by the motional electric field. △ Less

Submitted 9 May, 2018; originally announced May 2018.

Comments: 11 pages, 7 figures

Journal ref: A&A 616, A21 (2018)

Chimera States in Continuous Media: Existence and Distinctness

Authors: Zachary G. Nicolaou, Hermann Riecke, Adilson E. Motter

Abstract: The defining property of chimera states is the coexistence of coherent and incoherent domains in systems that are structurally and spatially homogeneous. The recent realization that such states might be common in oscillator networks raises the question of whether an analogous phenomenon can occur in continuous media. Here, we show that chimera states can exist in continuous systems even when the c… ▽ More The defining property of chimera states is the coexistence of coherent and incoherent domains in systems that are structurally and spatially homogeneous. The recent realization that such states might be common in oscillator networks raises the question of whether an analogous phenomenon can occur in continuous media. Here, we show that chimera states can exist in continuous systems even when the coupling is strictly local, as in many fluid and pattern forming media. Using the complex Ginzburg-Landau equation as a model system, we characterize chimera states consisting of a coherent domain of a frozen spiral structure and an incoherent domain of amplitude turbulence. We show that in this case, in contrast with discrete network systems, fluctuations in the local coupling field play a crucial role in limiting the coherent regions. We suggest these findings shed light on new possible forms of coexisting order and disorder in fluid systems. △ Less

Submitted 18 December, 2017; v1 submitted 1 December, 2017; originally announced December 2017.

Comments: Supplemental animation at article webpage: http://link.aps.org/supplemental/10.1103/PhysRevLett.119.244101

Journal ref: Phys. Rev. Lett. 119, 244101 (2017)

Misestimation of temperature when applying Maxwellian distributions to space plasmas described by kappa distributions

Authors: Georgios Nicolaou, George Livadiotis

Abstract: This paper presents the misestimation of temperature when observations from a kappa distributed plasma are analyzed as a Maxwellian. One common method to calculate the space plasma parameters is by fitting the observed distributions using known analytical forms. More often, the distribution function is included in a forward model of the instrument's response, which is used to reproduce the observe… ▽ More This paper presents the misestimation of temperature when observations from a kappa distributed plasma are analyzed as a Maxwellian. One common method to calculate the space plasma parameters is by fitting the observed distributions using known analytical forms. More often, the distribution function is included in a forward model of the instrument's response, which is used to reproduce the observed energy spectrograms for a given set of plasma parameters. In both cases, the modeled plasma distribution fits the measurements to estimate the plasma parameters. The distribution function is often considered to be Maxwellian even though in many cases the plasma is better described by a kappa distribution. In this work we show that if the plasma is described by a kappa distribution, the derived temperature assuming Maxwell distribution can be significantly off. More specifically, we derive the plasma temperature by fitting a Maxwell distribution to pseudo-data produced by a kappa distribution, and then examine the difference of the derived temperature as a function of the kappa index. We further consider the concept of using a forward model of a typical plasma instrument to fit its observations. We find that the relative error of the derived temperature is highly depended on the kappa index and occasionally on the instrument's field of view and response. △ Less

Submitted 6 October, 2016; originally announced October 2016.

Comments: 7 Figures and 19 pages The article is currently in press by "Astrophysics and Space Science" journal The articles number is ASTR-D-16-00131R4

New Horizons Solar Wind Around Pluto (SWAP) Observations of the Solar Wind From 11-33 AU

Authors: H. A. Elliott, D. J. McComas, P. Valek, G. Nicolaou, S. Weidner, G. Livadiotis

Abstract: The Solar Wind Around Pluto (SWAP) instrument on NASA's New Horizon Pluto mission has collected solar wind observations en route from Earth to Pluto, and these observations continue beyond Pluto. Few missions have explored the solar wind in the outer heliosphere making this dataset a critical addition to the field. We created a forward model of SWAP count rates, which includes a comprehensive inst… ▽ More The Solar Wind Around Pluto (SWAP) instrument on NASA's New Horizon Pluto mission has collected solar wind observations en route from Earth to Pluto, and these observations continue beyond Pluto. Few missions have explored the solar wind in the outer heliosphere making this dataset a critical addition to the field. We created a forward model of SWAP count rates, which includes a comprehensive instrument response function based on laboratory and flight calibrations. By fitting the count rates with this model, the proton density (n), speed (V), and temperature (T) parameters are determined. Comparisons between SWAP parameters and both propagated 1 AU observations and prior Voyager 2 observations indicate consistency in both the range and mean wind values. These comparisons as well as our additional findings confirm that small and midsized solar wind structures are worn down with increasing distance due to dynamic interaction of parcels of wind with different speed. For instance, the T-V relationship steepens, as the range in V is limited more than the range in T with distance. At times the T-V correlation clearly breaks down beyond 20 AU, which may indicate wind currently expanding and cooling may have an elevated T reflecting prior heating and compression in the inner heliosphere. The power of wind parameters at shorter periodicities decreases with distance as the longer periodicities strengthen. The solar rotation periodicity is present in temperature beyond 20 AU indicating the observed parcel temperature may reflect not only current heating or cooling, but also heating occurring closer to the Sun. △ Less

Submitted 26 January, 2016; originally announced January 2016.

Comments: 55 pages, 29 Figures, accepted for publication in The Astrophysical Journal Supplements (ApJS)

Longitudinal inverted compressibility in super-strained metamaterials

Authors: Zachary G. Nicolaou, Adilson E. Motter

Abstract: We develop a statistical physics theory for solid-solid phase transitions in which a metamaterial undergoes longitudinal contraction in response to increase in external tension. Such transitions, which are forbidden in thermodynamic equilibrium, have recently been shown to be possible during the decay of metastable, super-strained states. We present a first-principles model to predict these transi… ▽ More We develop a statistical physics theory for solid-solid phase transitions in which a metamaterial undergoes longitudinal contraction in response to increase in external tension. Such transitions, which are forbidden in thermodynamic equilibrium, have recently been shown to be possible during the decay of metastable, super-strained states. We present a first-principles model to predict these transitions and validate it using molecular dynamics simulations. Aside from its immediate mechanical implications, our theory points to a wealth of analogous inverted responses, such as inverted susceptibility or heat-capacity transitions, allowed when considering realistic scales. △ Less

Submitted 2 April, 2013; originally announced April 2013.

Journal ref: J. Stat. Phys. 151, 1162-1174 (2013)

Mechanical Metamaterials with Negative Compressibility Transitions

Authors: Zachary G. Nicolaou, Adilson E. Motter

Abstract: When tensioned, ordinary materials expand along the direction of the applied force. Here, we explore network concepts to design metamaterials exhibiting negative compressibility transitions, during which a material undergoes contraction when tensioned (or expansion when pressured). Continuous contraction of a material in the same direction of an applied tension, and in response to this tension, is… ▽ More When tensioned, ordinary materials expand along the direction of the applied force. Here, we explore network concepts to design metamaterials exhibiting negative compressibility transitions, during which a material undergoes contraction when tensioned (or expansion when pressured). Continuous contraction of a material in the same direction of an applied tension, and in response to this tension, is inherently unstable. The conceptually similar effect we demonstrate can be achieved, however, through destabilisations of (meta)stable equilibria of the constituents. These destabilisations give rise to a stress-induced solid-solid phase transition associated with a twisted hysteresis curve for the stress-strain relationship. The strain-driven counterpart of negative compressibility transitions is a force amplification phenomenon, where an increase in deformation induces a discontinuous increase in response force. We suggest that the proposed materials could be useful for the design of actuators, force amplifiers, micro-mechanical controls, and protective devices. △ Less

Submitted 9 July, 2012; originally announced July 2012.

Comments: Supplementary information available at http://www.nature.com/nmat/journal/v11/n7/abs/nmat3331.html

Journal ref: Nature Materials 11, 608-613 (2012)