-
N-body linear force law allowing analytic solutions
Authors:
Joseph West,
Sean P. Bartz
Abstract:
We present a pair-wise force law in a system of N particles that produces analytic solutions for arbitrary number of particles, masses, and initial conditions. Each pair of particles interacts via a force that is proportional to the product of their masses and their separation distance, with the force directed radially. We show that, despite the N-body interaction, each particle behaves as if it i…
▽ More
We present a pair-wise force law in a system of N particles that produces analytic solutions for arbitrary number of particles, masses, and initial conditions. Each pair of particles interacts via a force that is proportional to the product of their masses and their separation distance, with the force directed radially. We show that, despite the N-body interaction, each particle behaves as if it interacts only with the center of mass of the system. This effective two-body interaction behaves as Hooke's Law with a common frequency for all particles, with the familiar analytic solutions for the trajectories. With these analytic solutions, it is possible to efficiently simulate a collection of these particles and incorporate other external forces. As an example, we simulate the particles within an adiabatically expanding container and calculate pressure and temperature in both the attractive and repulsive cases.
△ Less
Submitted 2 July, 2024;
originally announced July 2024.
-
MR sequence design using digital twins of non-idealized hardware
Authors:
Daniel J West,
Felix Glang,
Jonathan Endres,
David Leitão,
Moritz Zaiss,
Joseph V Hajnal,
Shaihan J Malik
Abstract:
MRI systems are traditionally engineered to produce close to idealized performance, enabling a simplified pulse sequence design philosophy. An example of this is control of eddy currents produced by gradient fields; usually these are compensated by pre-emphasizing demanded waveforms. This process typically happens invisibly to the pulse designer, allowing them to assume the achieved gradient wavef…
▽ More
MRI systems are traditionally engineered to produce close to idealized performance, enabling a simplified pulse sequence design philosophy. An example of this is control of eddy currents produced by gradient fields; usually these are compensated by pre-emphasizing demanded waveforms. This process typically happens invisibly to the pulse designer, allowing them to assume the achieved gradient waveform will be as desired. Whilst convenient, this imposes stricter limits on the sequence design than the hardware can handle (for example, pre-emphasis adds an additional overhead to amplifiers). This strategy can be undesirable particularly for lower performance or resource-limited hardware. Instead we explore the use of a 'digital twin' (i.e. an end-to-end model of the scanner system) to optimize control inputs, resulting in sequences that inherently compensate for known imperfections. We explore digital twin optimization specifically for gradient system imperfections as an exemplar. This is first explored in simulations using a simple exponential eddy current model, then experimentally using an empirical gradient impulse response function on a 7T MRI system. When unconstrained, digital twin optimization reproduces classic pre-emphasis. When strict hardware constraints are imposed (simulating lower performance hardware), it identifies novel sequences for scenarios where classic pre-emphasis would be unachievable. Experimentally, the optimization approach was demonstrated to substantially reduce ghosting effects in echo planar images on a 7T system. Digital twin optimization may allow more efficient use of hardware by taking a whole system approach. Ultimately this could enable the use of cheaper hardware without loss of performance.
△ Less
Submitted 26 March, 2024;
originally announced March 2024.
-
The ZK Spiral Family of Trajectories and Forces
Authors:
Joseph West
Abstract:
Particle trajectories in the form of a logarithmic spiral with specified angular time dependence, "ZK spirals," are shown to be analytic solutions for motion in non-central, but simple force power-laws. Each ZK spiral is a particular solution to a single associated force law. The position and velocity are determined as analytic functions of time in terms of the growth parameter, specified power la…
▽ More
Particle trajectories in the form of a logarithmic spiral with specified angular time dependence, "ZK spirals," are shown to be analytic solutions for motion in non-central, but simple force power-laws. Each ZK spiral is a particular solution to a single associated force law. The position and velocity are determined as analytic functions of time in terms of the growth parameter, specified power law, and the initial conditions. Four examples, each for a different force law and relevant to the advanced classroom, are presented: the well-known attractive central inverse cube force; a bead on a rigid, horizontal, frictionless wire in the shape of the spiral trajectory; a car moving in a changing radius turn; and a known solution for a powered rocket with variable thrust in a Newtonian gravitational field. In this last case, general expressions for position and velocity of powered transfer orbits within the solar system are presented along with a new analytic expression for the rocket mass as a function of time or distance from the Sun. These results can be extended to include powered flight solutions to a generalized Lambert's problem for circular orbits in attractive central force laws of the form F(r) = - Fo/(r^q), for q greater than or equal to 1.
△ Less
Submitted 7 June, 2024; v1 submitted 10 January, 2024;
originally announced January 2024.
-
Neural network models for preferential concentration of particles in two-dimensional turbulence
Authors:
Thibault Maurel-Oujia,
Suhas S. Jain,
Keigo Matsuda,
Kai Schneider,
Jacob R. West,
Kazuki Maeda
Abstract:
Cluster and void formations are key processes in the dynamics of particle-laden turbulence. In this work, we assess the performance of various neural network models for synthesizing preferential concentration fields of particles in turbulence. A database of direct numerical simulations of homogeneous isotropic two-dimensional turbulence with one-way coupled inertial point particles, is used to tra…
▽ More
Cluster and void formations are key processes in the dynamics of particle-laden turbulence. In this work, we assess the performance of various neural network models for synthesizing preferential concentration fields of particles in turbulence. A database of direct numerical simulations of homogeneous isotropic two-dimensional turbulence with one-way coupled inertial point particles, is used to train the models using vorticity as the input to predict the particle number density fields. We compare autoencoder, U--Net, generative adversarial network (GAN), and diffusion model approaches, and assess the statistical properties of the generated particle number density fields. We find that the GANs are superior in predicting clusters and voids, and therefore result in the best performance. Additionally, we explore a concept of ``supersampling", where neural networks can be used to predict full particle data using only the information of few particles, which yields promising perspectives for reducing the computational cost of expensive DNS computations by avoiding the tracking of millions of particles. We also explore the inverse problem of synthesizing the enstrophy fields using the particle number density distribution as the input at different Stokes numbers. Hence, our study also indicates the potential use of neural networks to predict turbulent flow statistics using experimental measurements of inertial particles.
△ Less
Submitted 22 December, 2023;
originally announced December 2023.
-
A Chromatic Treatment of Linear Polarization in the Solar Corona at the 2023 Total Solar Eclipse
Authors:
Ritesh Patel,
Daniel B. Seaton,
Amir Caspi,
Sarah A. Kovac,
Sarah J. Davis,
John P. Carini,
Charles H. Gardner,
Sanjay Gosain,
Viliam Klein,
Shawn A. Laatsch,
Patricia H. Reiff,
Nikita Saini,
Rachael Weir,
Daniel W. Zietlow,
David F. Elmore,
Andrei E. Ursache,
Craig E. DeForest,
Matthew J. West,
Fred Bruenjes,
Jen Winter
Abstract:
The broadband solar K-corona is linearly polarized due to Thomson scattering. Various strategies have been used to represent coronal polarization. Here, we present a new way to visualize the polarized corona, using observations from the 2023 April 20 total solar eclipse in Australia in support of the Citizen CATE 2024 project. We convert observations in the common four-polarizer orthogonal basis (…
▽ More
The broadband solar K-corona is linearly polarized due to Thomson scattering. Various strategies have been used to represent coronal polarization. Here, we present a new way to visualize the polarized corona, using observations from the 2023 April 20 total solar eclipse in Australia in support of the Citizen CATE 2024 project. We convert observations in the common four-polarizer orthogonal basis (0°, 45°, 90°, & 135°) to -60°, 0°, and +60° (MZP) polarization, which is homologous to R, G, B color channels. The unique image generated provides some sense of how humans might visualize polarization if we could perceive it in the same way we perceive color.
△ Less
Submitted 14 November, 2023;
originally announced December 2023.
-
The Critical Coronal Transition Region: A Physics-framed Strategy to Uncover the Genesis of the Solar Wind and Solar Eruptions
Authors:
Angelos Vourlidas,
Amir Caspi,
Yuan-Kuen Ko,
J. Martin Laming,
James P. Mason,
Mari Paz Miralles,
Nour-Eddine Raouafi,
John C. Raymond,
Daniel B. Seaton,
Leonard Strachan,
Nicholeen Viall,
Juliana Vievering,
Matthew J. West
Abstract:
Our current theoretical and observational understanding suggests that critical properties of the solar wind and Coronal Mass Ejections (CMEs) are imparted within 10 Rs, particularly below 4 Rs. This seemingly narrow spatial region encompasses the transition of coronal plasma processes through the entire range of physical regimes from fluid to kinetic, and from primarily closed to open magnetic fie…
▽ More
Our current theoretical and observational understanding suggests that critical properties of the solar wind and Coronal Mass Ejections (CMEs) are imparted within 10 Rs, particularly below 4 Rs. This seemingly narrow spatial region encompasses the transition of coronal plasma processes through the entire range of physical regimes from fluid to kinetic, and from primarily closed to open magnetic field structures. From a physics perspective, therefore, it is more appropriate to refer to this region as the Critical Coronal Transition Region (CCTR) to emphasize its physical, rather than spatial, importance to key Heliophysics science.
This white paper argues that the comprehensive exploration of the CCTR will answer two of the most central Heliophysics questions, "How and where does the solar wind form?" and "How do eruptions form?", by unifying hardware/software/modeling development and seemingly disparate research communities and frameworks. We describe the outlines of decadal-scale plan to achieve that by 2050.
△ Less
Submitted 25 July, 2023;
originally announced July 2023.
-
Magnetic Energy Powers the Corona: How We Can Understand its 3D Storage & Release
Authors:
Amir Caspi,
Daniel B. Seaton,
Roberto Casini,
Cooper Downs,
Sarah E. Gibson,
Holly Gilbert,
Lindsay Glesener,
Silvina E. Guidoni,
J. Marcus Hughes,
David McKenzie,
Joseph Plowman,
Katharine K. Reeves,
Pascal Saint-Hilaire,
Albert Y. Shih,
Matthew J. West
Abstract:
The coronal magnetic field is the prime driver behind many as-yet unsolved mysteries: solar eruptions, coronal heating, and the solar wind, to name a few. It is, however, still poorly observed and understood. We highlight key questions related to magnetic energy storage, release, and transport in the solar corona, and their relationship to these important problems. We advocate for new and multi-po…
▽ More
The coronal magnetic field is the prime driver behind many as-yet unsolved mysteries: solar eruptions, coronal heating, and the solar wind, to name a few. It is, however, still poorly observed and understood. We highlight key questions related to magnetic energy storage, release, and transport in the solar corona, and their relationship to these important problems. We advocate for new and multi-point co-optimized measurements, sensitive to magnetic field and other plasma parameters, spanning from optical to $γ$-ray wavelengths, to bring closure to these long-standing and fundamental questions. We discuss how our approach can fully describe the 3D magnetic field, embedded plasma, particle energization, and their joint evolution to achieve these objectives.
△ Less
Submitted 25 May, 2023;
originally announced May 2023.
-
Improving Multi-Dimensional Data Formats, Access, and Assimilation Tools for the Twenty-First Century
Authors:
Daniel B. Seaton,
Amir Caspi,
Roberto Casini,
Cooper Downs,
Sarah E. Gibson,
Holly Gilbert,
Lindsay Glesener,
Silvina E. Guidoni,
J. Marcus Hughes,
David McKenzie,
Joseph Plowman,
Katharine K. Reeves,
Pascal Saint-Hilaire,
Albert Y. Shih,
Matthew J. West
Abstract:
Heliophysics image data largely relies on a forty-year-old ecosystem built on the venerable Flexible Image Transport System (FITS) data standard. While many in situ measurements use newer standards, they are difficult to integrate with multiple data streams required to develop global understanding. Additionally, most data users still engage with data in much the same way as they did decades ago. H…
▽ More
Heliophysics image data largely relies on a forty-year-old ecosystem built on the venerable Flexible Image Transport System (FITS) data standard. While many in situ measurements use newer standards, they are difficult to integrate with multiple data streams required to develop global understanding. Additionally, most data users still engage with data in much the same way as they did decades ago. However, contemporary missions and models require much more complex support for 3D multi-parameter data, robust data assimilation strategies, and integration of multiple individual data streams required to derive complete physical characterizations of the Sun and Heliospheric plasma environment. In this white paper we highlight some of the 21$^\mathsf{st}$ century challenges for data frameworks in heliophysics, consider an illustrative case study, and make recommendations for important steps the field can take to modernize its data products and data usage models. Our specific recommendations include: (1) Investing in data assimilation capability to drive advanced data-constrained models, (2) Investing in new strategies for integrating data across multiple instruments to realize measurements that cannot be produced from single observations, (3) Rethinking old data use paradigms to improve user access, develop deep understanding, and decrease barrier to entry for new datasets, and (4) Investing in research on data formats better suited for multi-dimensional data and cloud-based computing.
△ Less
Submitted 25 May, 2023;
originally announced May 2023.
-
COMPLETE: A flagship mission for complete understanding of 3D coronal magnetic energy release
Authors:
Amir Caspi,
Daniel B. Seaton,
Roberto Casini,
Cooper Downs,
Sarah E. Gibson,
Holly Gilbert,
Lindsay Glesener,
Silvina E. Guidoni,
J. Marcus Hughes,
David McKenzie,
Joseph Plowman,
Katharine K. Reeves,
Pascal Saint-Hilaire,
Albert Y. Shih,
Matthew J. West
Abstract:
COMPLETE is a flagship mission concept combining broadband spectroscopic imaging and comprehensive magnetography from multiple viewpoints around the Sun to enable tomographic reconstruction of 3D coronal magnetic fields and associated dynamic plasma properties, which provide direct diagnostics of energy release. COMPLETE re-imagines the paradigm for solar remote-sensing observations through purpos…
▽ More
COMPLETE is a flagship mission concept combining broadband spectroscopic imaging and comprehensive magnetography from multiple viewpoints around the Sun to enable tomographic reconstruction of 3D coronal magnetic fields and associated dynamic plasma properties, which provide direct diagnostics of energy release. COMPLETE re-imagines the paradigm for solar remote-sensing observations through purposefully co-optimized detectors distributed on multiple spacecraft that operate as a single observatory, linked by a comprehensive data/model assimilation strategy to unify individual observations into a single physical framework. We describe COMPLETE's science goals, instruments, and mission implementation. With targeted investment by NASA, COMPLETE is feasible for launch in 2032 to observe around the maximum of Solar Cycle 26.
△ Less
Submitted 25 May, 2023;
originally announced May 2023.
-
Unveiling Pseudo-Crucial Events in Noise-Induced Phase Transitions
Authors:
Jacob D. Baxley,
David R. Lambert,
Mauro Bologna,
Bruce J. West,
Paolo Grigolini
Abstract:
Noise-induced phase transitions are common in various complex systems, from physics to biology. In this article, we investigate the emergence of crucial events in noise-induced phase transition processes and their potential significance for understanding complexity in such systems. We utilize the first-passage time technique and coordinate transformations to study the dynamics of the system and id…
▽ More
Noise-induced phase transitions are common in various complex systems, from physics to biology. In this article, we investigate the emergence of crucial events in noise-induced phase transition processes and their potential significance for understanding complexity in such systems. We utilize the first-passage time technique and coordinate transformations to study the dynamics of the system and identify crucial events. Furthermore, we employ Diffusion Entropy Analysis, a powerful statistical tool, to characterize the complexity of the system and quantify the information content of the identified events. Our results show that the emergence of crucial events is closely related to the complexity of the system and can provide insight into its behavior. This approach may have applications in diverse fields, such as climate modeling, financial markets, and biological systems, where understanding the emergence of crucial events is of great importance.
△ Less
Submitted 16 April, 2023; v1 submitted 11 April, 2023;
originally announced April 2023.
-
Radio Studies of the Middle Corona: Current State and New Prospects in the Next Decade
Authors:
Bin Chen,
Jason E. Kooi,
David B. Wexler,
Dale E. Gary,
Sijie Yu,
Surajit Mondal,
Adam R. Kobelski,
Daniel B. Seaton,
Matthew J. West,
Stephen M. White,
Gregory D. Fleishman,
Pascal Saint-Hilaire,
Peijin Zhang,
Chris R. Gilly,
James P. Mason,
Hamish Reid
Abstract:
The "middle corona," defined by West et al. (2022) as the region between ~1.5-6 solar radii, is a critical transition region that connects the highly structured lower corona to the outer corona where the magnetic field becomes predominantly radial. At radio wavelengths, remote-sensing of the middle corona falls in the meter-decameter wavelength range where a critical transition of radio emission m…
▽ More
The "middle corona," defined by West et al. (2022) as the region between ~1.5-6 solar radii, is a critical transition region that connects the highly structured lower corona to the outer corona where the magnetic field becomes predominantly radial. At radio wavelengths, remote-sensing of the middle corona falls in the meter-decameter wavelength range where a critical transition of radio emission mechanisms occurs. In addition, plasma properties of the middle corona can be probed by trans-coronal radio propagation methods including radio scintillation and Faraday rotation techniques. Together they offer a wealth of diagnostic tools for the middle corona, complementing current and planned missions at other wavelengths. These diagnostics include unique means for detecting and measuring the magnetic field and energetic electrons associated with coronal mass ejections, mapping coronal shocks and electron beam trajectories, as well as constraining the plasma density, magnetic field, and turbulence of the "young" solar wind. Following a brief overview of pertinent radio diagnostic methods, this white paper will discuss the current state of radio studies on the middle corona, challenges to obtaining a more comprehensive picture, and recommend an outlook in the next decade. Our specific recommendations for advancing the middle coronal sciences from the radio perspective are: (1) Prioritizing solar-dedicated radio facilities in the ~0.1-1 GHz range with broadband, high-dynamic-range imaging spectropolarimetry capabilities. (2) Developing facilities and techniques to perform multi-perspective, multiple lines-of-sight trans-coronal radio Faraday Rotation measurements.
△ Less
Submitted 28 January, 2023;
originally announced January 2023.
-
First Perihelion of EUI on the Solar Orbiter mission
Authors:
D. Berghmans,
P. Antolin,
F. Auchère,
R. Aznar Cuadrado,
K. Barczynski,
L. P. Chitta,
S. Gissot,
L. Harra,
Z. Huang,
M. Janvier,
E. Kraaikamp,
D. M. Long,
S. Mandal,
M. Mierla,
S. Parenti,
H. Peter,
L. Rodriguez,
U. Schühle,
P. J. Smith,
S. K. Solanki,
K. Stegen,
L. Teriaca,
C. Verbeeck,
M. J. West,
A. N. Zhukov
, et al. (12 additional authors not shown)
Abstract:
Context. The Extreme Ultraviolet Imager (EUI), onboard Solar Orbiter consists of three telescopes: the two High Resolution Imagers in EUV (HRIEUV) and in Lyman-α (HRILya), and the Full Sun Imager (FSI). Solar Orbiter/EUI started its Nominal Mission Phase on 2021 November 27. Aims. EUI images from the largest scales in the extended corona off limb, down to the smallest features at the base of the c…
▽ More
Context. The Extreme Ultraviolet Imager (EUI), onboard Solar Orbiter consists of three telescopes: the two High Resolution Imagers in EUV (HRIEUV) and in Lyman-α (HRILya), and the Full Sun Imager (FSI). Solar Orbiter/EUI started its Nominal Mission Phase on 2021 November 27. Aims. EUI images from the largest scales in the extended corona off limb, down to the smallest features at the base of the corona and chromosphere. EUI is therefore a key instrument for the connection science that is at the heart of the Solar Orbiter mission science goals. Methods. The highest resolution on the Sun is achieved when Solar Orbiter passes through the perihelion part of its orbit. On 2022 March 26, Solar Orbiter reached for the first time a distance to the Sun close to 0.3 au. No other coronal EUV imager has been this close to the Sun. Results. We review the EUI data sets obtained during the period 2022 March-April, when Solar Orbiter quickly moved from alignment with the Earth (2022 March 6), to perihelion (2022 March 26), to quadrature with the Earth (2022 March 29). We highlight the first observational results in these unique data sets and we report on the in-flight instrument performance. Conclusions. EUI has obtained the highest resolution images ever of the solar corona in the quiet Sun and polar coronal holes. Several active regions were imaged at unprecedented cadences and sequence durations. We identify in this paper a broad range of features that require deeper studies. Both FSI and HRIEUV operate at design specifications but HRILya suffered from performance issues near perihelion. We conclude emphasising the EUI open data policy and encouraging further detailed analysis of the events highlighted in this paper.
△ Less
Submitted 13 January, 2023;
originally announced January 2023.
-
Defining the Middle Corona
Authors:
Matthew J. West,
Daniel B. Seaton,
David B. Wexler,
John C. Raymond,
Giulio Del Zanna,
Yeimy J. Rivera,
Adam R. Kobelski,
Craig DeForest,
Leon Golub,
Amir Caspi,
Chris R. Gilly,
Jason E. Kooi,
Benjamin L. Alterman,
Nathalia Alzate,
Dipankar Banerjee,
David Berghmans,
Bin Chen,
Lakshmi Pradeep Chitta,
Cooper Downs,
Silvio Giordano,
Aleida Higginson,
Russel A. Howard,
Emily Mason,
James P. Mason,
Karen A. Meyer
, et al. (9 additional authors not shown)
Abstract:
The middle corona, the region roughly spanning heliocentric altitudes from $1.5$ to $6\,R_\odot$, encompasses almost all of the influential physical transitions and processes that govern the behavior of coronal outflow into the heliosphere. Eruptions that could disrupt the near-Earth environment propagate through it. Importantly, it modulates inflow from above that can drive dynamic changes at low…
▽ More
The middle corona, the region roughly spanning heliocentric altitudes from $1.5$ to $6\,R_\odot$, encompasses almost all of the influential physical transitions and processes that govern the behavior of coronal outflow into the heliosphere. Eruptions that could disrupt the near-Earth environment propagate through it. Importantly, it modulates inflow from above that can drive dynamic changes at lower heights in the inner corona. Consequently, this region is essential for comprehensively connecting the corona to the heliosphere and for developing corresponding global models. Nonetheless, because it is challenging to observe, the middle corona has been poorly studied by major solar remote sensing missions and instruments, extending back to the Solar and Heliospheric Observatory (SoHO) era. Thanks to recent advances in instrumentation, observational processing techniques, and a realization of the importance of the region, interest in the middle corona has increased. Although the region cannot be intrinsically separated from other regions of the solar atmosphere, there has emerged a need to define the region in terms of its location and extension in the solar atmosphere, its composition, the physical transitions it covers, and the underlying physics believed to be encapsulated by the region. This paper aims to define the middle corona and give an overview of the processes that occur there.
△ Less
Submitted 9 March, 2023; v1 submitted 8 August, 2022;
originally announced August 2022.
-
The Sabatier principle for Battery Anodes: Chemical Kinetics and Reversible Electrodeposition at Heterointerfaces
Authors:
Jingxu Zheng,
Yue Deng,
Wenzao Li,
Jiefu Yin,
Patrick J. West,
Tian Tang,
Xiao Tong,
David C. Bock,
Shuo Jin,
Qing Zhao,
Regina Garcia-Mendez,
Kenneth J. Takeuchi,
Esther S. Takeuchi,
Amy C. Marschilok,
Lynden A. Archer
Abstract:
How surface chemistry influences reactions occurring thereupon has been a long-standing question of broad scientific and technological interest for centuries. Recently, it has re-emerged as a critical question in a subdiscipline of chemistry - electrochemistry at heterointerphases, where the answers have implications for both how, and in what forms, humanity stores the rising quantities of renewab…
▽ More
How surface chemistry influences reactions occurring thereupon has been a long-standing question of broad scientific and technological interest for centuries. Recently, it has re-emerged as a critical question in a subdiscipline of chemistry - electrochemistry at heterointerphases, where the answers have implications for both how, and in what forms, humanity stores the rising quantities of renewable electric power generated from solar and wind installations world-wide. Here we consider the relation between the surface chemistry at such interphases and the reversibility of electrochemical transformations at a rechargeable battery electrode. Conventional wisdom holds that stronger chemical interaction between the metal deposits and electrode promotes reversibility. We report instead that a moderate strength of chemical interaction between the deposit and the substrate, neither too weak nor too strong, enables highest reversibility and stability of the plating/stripping redox processes at a battery anode. Analogous to the empirical Sabatier principle for chemical heterogeneous catalysis, our finding arises from the confluence of competing processes - one driven by electrochemistry and the other by chemical alloying. Based on experimental evaluation of metal plating/stripping systems in battery anodes of contemporary interest, we show that such knowledge provides a powerful tool for designing key materials in highly reversible electrochemical energy storage technologies based on earth-abundant, low-cost metals.
△ Less
Submitted 25 September, 2022; v1 submitted 14 March, 2022;
originally announced March 2022.
-
Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids
Authors:
Suhas S. Jain,
Michael C. Adler,
Jacob R. West,
Ali Mani,
Parviz Moin,
Sanjiva K. Lele
Abstract:
This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook (2007), Subramaniam et al., (2018), and Adler and Lele (2019), in which artificial diffusion terms are added to the conservation equations. The second method is…
▽ More
This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook (2007), Subramaniam et al., (2018), and Adler and Lele (2019), in which artificial diffusion terms are added to the conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., (2010) and Tiwari et al., (2013), in which the diffusion and sharpening terms (together called regularization terms) are added to the conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., (2020), in which the regularization terms are added to the conservation equations.
The primary objective of this work is to compare these three methods in terms of their ability to maintain: constant interface thickness; the conservation property; and accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer-Meshkov instability of a copper-aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. This is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.
△ Less
Submitted 20 September, 2021;
originally announced September 2021.
-
A Direct Detection Search for Hidden Sector New Particles in the 3-60 MeV Mass Range
Authors:
A. Ahmidouch,
S. Davis,
A. Gasparian,
T. J. Hague,
S. Mtingwa,
R. Pedroni,
C. Ayerbe-Gayoso,
H. Bhatt,
B. Devkota,
J. Dunne,
D. Dutta,
L. El Fassi,
A. Karki,
P. Mohanmurthy,
C. Peng,
S. Ali,
X. Bai,
J. Boyd,
B. Dharmasena,
V. Gamage,
K. Gnanvo,
S. Jeffas,
S. Jian,
N. Liyanage,
H. Nguyen
, et al. (36 additional authors not shown)
Abstract:
In our quest to understand the nature of dark matter and discover its non-gravitational interactions with ordinary matter, we propose an experiment using a \pbo ~calorimeter to search for or set new limits on the production rate of i) hidden sector particles in the $3 - 60$ MeV mass range via their $e^+e^-$ decay (or $γγ$ decay with limited tracking), and ii) the hypothetical X17 particle, claimed…
▽ More
In our quest to understand the nature of dark matter and discover its non-gravitational interactions with ordinary matter, we propose an experiment using a \pbo ~calorimeter to search for or set new limits on the production rate of i) hidden sector particles in the $3 - 60$ MeV mass range via their $e^+e^-$ decay (or $γγ$ decay with limited tracking), and ii) the hypothetical X17 particle, claimed in multiple recent experiments. The search for these particles is motivated by new hidden sector models and dark matter candidates introduced to account for a variety of experimental and observational puzzles: the small-scale structure puzzle in cosmological simulations, anomalies such as the 4.2$σ$ disagreement between experiments and the standard model prediction for the muon anomalous magnetic moment, and the excess of $e^+e^-$ pairs from the $^8$Be M1 and $^4$He nuclear transitions to their ground states observed by the ATOMKI group. In these models, the $1 - 100$ MeV mass range is particularly well-motivated and the lower part of this range still remains unexplored. Our proposed direct detection experiment will use a magnetic-spectrometer-free setup (the PRad apparatus) to detect all three final state particles in the visible decay of a hidden sector particle allowing for an effective control of the background and will cover the proposed mass range in a single setting. The use of the well-demonstrated PRad setup allows for an essentially ready-to-run and uniquely cost-effective search for hidden sector particles in the $3 - 60$ MeV mass range with a sensitivity of 8.9$\times$10$^{-8}$ - 5.8$\times$10$^{-9}$ to $ε^2$, the square of the kinetic mixing interaction constant between hidden and visible sectors. This updated proposal includes our response to the PAC49 comments.
△ Less
Submitted 4 August, 2022; v1 submitted 30 August, 2021;
originally announced August 2021.
-
Magnetization Transfer-Mediated MR Fingerprinting
Authors:
Daniel J. West,
Gastao Cruz,
Rui P. A. G. Teixeira,
Torben Schneider,
Jacques-Donald Tournier,
Joseph V. Hajnal,
Claudia Prieto,
Shaihan J. Malik
Abstract:
Purpose: Magnetization transfer (MT) and inhomogeneous MT (ihMT) contrasts are used in MRI to provide information about macromolecular tissue content. In particular, MT is sensitive to macromolecules and ihMT appears to be specific to myelinated tissue. This study proposes a technique to characterize MT and ihMT properties from a single acquisition, producing both semiquantitative contrast ratios,…
▽ More
Purpose: Magnetization transfer (MT) and inhomogeneous MT (ihMT) contrasts are used in MRI to provide information about macromolecular tissue content. In particular, MT is sensitive to macromolecules and ihMT appears to be specific to myelinated tissue. This study proposes a technique to characterize MT and ihMT properties from a single acquisition, producing both semiquantitative contrast ratios, and quantitative parameter maps.
Theory and Methods: Building upon previous work that uses multiband radiofrequency (RF) pulses to efficiently generate ihMT contrast, we propose a cyclic-steady-state approach that cycles between multiband and single-band pulses to boost the achieved contrast. Resultant time-variable signals are reminiscent of a magnetic resonance fingerprinting (MRF) acquisition, except that the signal fluctuations are entirely mediated by magnetization transfer effects. A dictionary-based low-rank inversion method is used to reconstruct the resulting images and to produce both semiquantitative MT ratio (MTR) and ihMT ratio (ihMTR) maps, as well as quantitative parameter estimates corresponding to an ihMT tissue model.
Results: Phantom and in vivo brain data acquired at 1.5T demonstrate the expected contrast trends, with ihMTR maps showing contrast more specific to white matter (WM), as has been reported by others. Quantitative estimation of semisolid fraction and dipolar T1 was also possible and yielded measurements consistent with literature values in the brain.
Conclusions: By cycling between multiband and single-band pulses, an entirely magnetization transfer mediated 'fingerprinting' method was demonstrated. This proof-of-concept approach can be used to generate semiquantitative maps and quantitatively estimate some macromolecular specific tissue parameters.
△ Less
Submitted 23 August, 2021; v1 submitted 6 April, 2021;
originally announced April 2021.
-
Three-Dimensional Reconstructions of Coronal Wave Surfaces Using a New Mask-Fitting Method
Authors:
Li Feng,
Lei Lu,
Bernd Inhester,
Joseph Plowman,
Beili Ying,
Marilena Mierla,
Matthew J. West,
Weiqun Gan
Abstract:
Coronal waves are large-scale disturbances often driven by coronal mass ejections (CMEs). We investigate a spectacular wave event on 7 March 2012, which is associated with an X5.4 flare (SOL2012-03-07). By using a running center-median (RCM) filtering method for the detection of temporal variations in extreme ultraviolet (EUV) images, we enhance the EUV disturbance observed by the Atmospheric Imag…
▽ More
Coronal waves are large-scale disturbances often driven by coronal mass ejections (CMEs). We investigate a spectacular wave event on 7 March 2012, which is associated with an X5.4 flare (SOL2012-03-07). By using a running center-median (RCM) filtering method for the detection of temporal variations in extreme ultraviolet (EUV) images, we enhance the EUV disturbance observed by the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) and the Sun Watcher using Active Pixel System detector and Image Processing (SWAP) onboard the PRoject for Onboard Autonomy 2 (PROBA2). In coronagraph images, a halo front is observed to be the upper counterpart of the EUV disturbance. Based on the EUV and coronagraph images observed from three different perspectives, we have made three-dimensional (3D) reconstructions of the wave surfaces using a new mask-fitting method. The reconstructions are compared with those obtained from forward-fitting methods. We show that the mask fitting method can reflect the inhomogeneous coronal medium by capturing the concave shape of the shock wave front. Subsequently, we trace the developing concave structure and derive the deprojected wave kinematics. The speed of the 3D-wave nose increases from a low value below a few hundred $\mathrm{km\,s^{-1}}$ to a maximum value of about 3800 $\mathrm{km\,s^{-1}}$, and then slowly decreases afterwards. The concave structure starts to decelerate earlier and has significantly lower speeds than those of the wave nose. We also find that the 3D-wave in the extended corona has a much higher speed than the speed of EUV disturbances across the solar disk.
△ Less
Submitted 22 September, 2020;
originally announced September 2020.
-
Rolling Polygons with Granular Material down an Incline
Authors:
Sidney Eckert,
Phillip Ingalls,
J. West
Abstract:
The present paper reports a novel behavior involving regular polygons with n sides and filled to varying degrees with granular materials. These are comprised of a set of hollow polygons produced on a 3D printer, and a single larger hollow hexagon fabricated with wooden sides and clear Plexiglas faces. Empty or full polygons stop immediately at shallow ramp angles, and roll the full length of the r…
▽ More
The present paper reports a novel behavior involving regular polygons with n sides and filled to varying degrees with granular materials. These are comprised of a set of hollow polygons produced on a 3D printer, and a single larger hollow hexagon fabricated with wooden sides and clear Plexiglas faces. Empty or full polygons stop immediately at shallow ramp angles, and roll the full length of the ramp for steep ramp angles, approaching a terminal velocity. This is consistent with results previously reported by other investigators of rolling solid polygons and partially filled cylinders. In contrast, partially filled polygons released at shallow ramp angles accelerate to a "terminal velocity," but then come to an abrupt stop. The distance traveled is reproducible and dependent on the ramp angle, number of sides and the volume fill ratio and is minimized when filled to 0.4 of its total internal volume. For larger ramp angles, the partially filled polygons again approach a "terminal velocity," but the terminal velocity is minimized, again near the fill ratio of 0.4. a simple model is introduced, but it is successful only in replicating the overall trend in velocity as a function of time for the large angle cases.
△ Less
Submitted 4 May, 2020;
originally announced May 2020.
-
A new facility for airborne solar astronomy: NASA's WB-57 at the 2017 total solar eclipse
Authors:
Amir Caspi,
Daniel B. Seaton,
Constantine C. C. Tsang,
Craig E. DeForest,
Paul Bryans,
Edward E. DeLuca,
Steven Tomczyk,
Joan T. Burkepile,
Thomas "Tony" Casey,
John Collier,
Donald "DD" Darrow,
Dominic Del Rosso,
Daniel D. Durda,
Peter T. Gallagher,
Leon Golub,
Matthew Jacyna,
David "DJ" Johnson,
Philip G. Judge,
Cary "Diddle" Klemm,
Glenn T. Laurent,
Johanna Lewis,
Charles J. Mallini,
Thomas "Duster" Parent,
Timothy Propp,
Andrew J. Steffl
, et al. (6 additional authors not shown)
Abstract:
NASA's WB-57 High Altitude Research Program provides a deployable, mobile, stratospheric platform for scientific research. Airborne platforms are of particular value for making coronal observations during total solar eclipses because of their ability both to follow the Moon's shadow and to get above most of the atmospheric airmass that can interfere with astronomical observations. We used the 2017…
▽ More
NASA's WB-57 High Altitude Research Program provides a deployable, mobile, stratospheric platform for scientific research. Airborne platforms are of particular value for making coronal observations during total solar eclipses because of their ability both to follow the Moon's shadow and to get above most of the atmospheric airmass that can interfere with astronomical observations. We used the 2017 Aug 21 eclipse as a pathfinding mission for high-altitude airborne solar astronomy, using the existing high-speed visible-light and near-/mid-wave infrared imaging suite mounted in the WB-57 nose cone. In this paper, we describe the aircraft, the instrument, and the 2017 mission; operations and data acquisition; and preliminary analysis of data quality from the existing instrument suite. We describe benefits and technical limitations of this platform for solar and other astronomical observations. We present a preliminary analysis of the visible-light data quality and discuss the limiting factors that must be overcome with future instrumentation. We conclude with a discussion of lessons learned from this pathfinding mission and prospects for future research at upcoming eclipses, as well as an evaluation of the capabilities of the WB-57 platform for future solar astronomy and general astronomical observation.
△ Less
Submitted 20 April, 2020;
originally announced April 2020.
-
The demography of the peripatetic researcher: Evidence on highly mobile scholars from the Web of Science
Authors:
Samin Aref,
Emilio Zagheni,
Jevin West
Abstract:
The policy debate around researchers' geographic mobility has been moving away from a theorized zero-sum game in which countries can be winners (brain gain) or losers (brain drain), and toward the concept of brain circulation, which implies that researchers move in and out of countries and everyone benefits. Quantifying trends in researchers' movements is key to understanding the drivers of the mo…
▽ More
The policy debate around researchers' geographic mobility has been moving away from a theorized zero-sum game in which countries can be winners (brain gain) or losers (brain drain), and toward the concept of brain circulation, which implies that researchers move in and out of countries and everyone benefits. Quantifying trends in researchers' movements is key to understanding the drivers of the mobility of talent, as well as the implications of these patterns for the global system of science, and for the competitive advantages of individual countries. Existing studies have investigated bilateral flows of researchers. However, in order to understand migration systems, determining the extent to which researchers have worked in more than two countries is essential. This study focuses on the subgroup of highly mobile researchers whom we refer to as peripatetic researchers or super-movers.
More specifically, our aim is to track the international movements of researchers who have published in more than two countries through changes in the main affiliation addresses of researchers in over 62 million publications indexed in the Web of Science database over the 1956-2016 period. Using this approach, we have established a longitudinal dataset on the international movements of highly mobile researchers across all subject categories, and in all disciplines of scholarship. This article contributes to the literature by offering for the first time a snapshot of the key features of highly mobile researchers, including their patterns of migration and return migration by academic age, the relative frequency of their disciplines, and the relative frequency of their countries of origin and destination. Among other findings, the results point to the emergence of a global system that includes the USA and China as two large hubs, and England and Germany as two smaller hubs for highly mobile researchers.
△ Less
Submitted 31 July, 2019;
originally announced July 2019.
-
Multipoint study of successive coronal mass ejections driving moderate disturbances at 1 AU
Authors:
Erika Palmerio,
Camilla Scolini,
David Barnes,
Jasmina Magdalenić,
Matthew J. West,
Andrei N. Zhukov,
Luciano Rodriguez,
Marilena Mierla,
Simon W. Good,
Diana E. Morosan,
Emilia K. J. Kilpua,
Jens Pomoell,
Stefaan Poedts
Abstract:
We analyse in this work the propagation and geoeffectiveness of four successive coronal mass ejections (CMEs) that erupted from the Sun during 21--23 May 2013 and that were detected in interplanetary space by the Wind and/or STEREO-A spacecraft. All these CMEs featured critical aspects for understanding so-called "problem space weather storms" at Earth. In the first three events a limb CMEs result…
▽ More
We analyse in this work the propagation and geoeffectiveness of four successive coronal mass ejections (CMEs) that erupted from the Sun during 21--23 May 2013 and that were detected in interplanetary space by the Wind and/or STEREO-A spacecraft. All these CMEs featured critical aspects for understanding so-called "problem space weather storms" at Earth. In the first three events a limb CMEs resulted in moderately geoeffective in-situ structures at their target location in terms of the disturbance storm time (Dst) index (either measured or estimated). The fourth CME, which also caused a moderate geomagnetic response, erupted from close to the disc centre as seen from Earth, but it was not visible in coronagraph images from the spacecraft along the Sun--Earth line and appeared narrow and faint from off-angle viewpoints. Making the correct connection between CMEs at the Sun and their in-situ counterparts is often difficult for problem storms. We investigate these four CMEs using multiwavelength and multipoint remote-sensing observations (extreme ultraviolet, white light, and radio), aided by 3D heliospheric modelling, in order to follow their propagation in the corona and in interplanetary space and to assess their impact at 1 AU. Finally, we emphasise the difficulties in forecasting moderate space weather effects provoked by problematic and ambiguous events and the importance of multispacecraft data for observing and modelling problem storms.
△ Less
Submitted 4 June, 2019;
originally announced June 2019.
-
Forecasting the Arrival Time of Coronal Mass Ejections: Analysis of the CCMC CME Scoreboard
Authors:
Pete Riley,
Leila Mays,
Jesse Andries,
Tanja Amerstorfer,
Douglas Biesecker,
Veronique Delouille,
Mateja Dumbovic,
Xueshang Feng,
Edmund Henley,
Jon A. Linker,
Christian Mostl,
Marlon Nunez,
Vic Pizzo,
Manuela Temmer,
W. K. Tobiska,
C. Verbeke,
Matthew J West,
Xinhua Zhao
Abstract:
Accurate forecasting of the properties of coronal mass ejections as they approach Earth is now recognized as an important strategic objective for both NOAA and NASA. The time of arrival of such events is a key parameter, one that had been anticipated to be relatively straightforward to constrain. In this study, we analyze forecasts submitted to the Community Coordinated Modeling Center (CCMC) at N…
▽ More
Accurate forecasting of the properties of coronal mass ejections as they approach Earth is now recognized as an important strategic objective for both NOAA and NASA. The time of arrival of such events is a key parameter, one that had been anticipated to be relatively straightforward to constrain. In this study, we analyze forecasts submitted to the Community Coordinated Modeling Center (CCMC) at NASA's Goddard Space Flight Center over the last six years to answer the following questions: (1) How well do these models forecast the arrival time of CME-driven shocks? (2) What are the uncertainties associated with these forecasts? (3) Which model(s) perform best? (4) Have the models become more accurate during the past six years? We analyze all forecasts made by 32 models from 2013 through mid 2018, and additionally focus on 28 events all of which were forecasted by six models. We find that the models are generally able to predict CME-shock arrival times -- in an average sense -- to within 10 hours, but with standard deviations often exceeding 20 hours. The best performers, on the other hand, maintained a mean error (bias) of -1 hour, a mean absolute error of 13 hours, and a precision (s.d.) of 15 hours. Finally, there is no evidence that the forecasts have become more accurate during this interval. We discuss the intrinsic simplifications of the various models analyzed, the limitations of this investigation, and suggest possible paths to improve these forecasts in the future.
△ Less
Submitted 16 October, 2018;
originally announced October 2018.
-
Signatures of Parkinson's disease in complexity of low-frequency fluctuations of postural sway velocity
Authors:
Miroslaw Latka,
Slawomir Budrewicz,
Klaudia Kozlowska,
Magdalena Koszewicz,
Bruce J. West
Abstract:
Posturography is routinely used to qualitatively assess one of the cardinal symptoms of Parkinson's disease -- postural instability. While most measures of balance control are derived from displacement of the center of pressure there is evidence that such control is more likely to be velocity-based. We performed static posturographic tests (eyes open and eyes closed) during quiet standing in narro…
▽ More
Posturography is routinely used to qualitatively assess one of the cardinal symptoms of Parkinson's disease -- postural instability. While most measures of balance control are derived from displacement of the center of pressure there is evidence that such control is more likely to be velocity-based. We performed static posturographic tests (eyes open and eyes closed) during quiet standing in narrow stance for n=30 patients with Parkinson's disease (PDP) in the ON state and compared the results with those of n=30 age-matched senior controls (HSC) and n=60 young controls (HYC). We used differentiator filters to generate time series of low-frequency fluctuations of sway velocity and calculated their Lempel-Ziv complexity (LZC). With eyes closed, the mediolateral LZC of HSC 0.21 (0.02) was significantly higher than those of HYC 0.19 (0.02) and PDP 0.18 (002). Thus aging and PD have opposite effects on mediolateral LZC which strongly differentiates between HSC and PDP (92% sensitivity and 87% specificity).
△ Less
Submitted 20 October, 2017;
originally announced October 2017.
-
Stem-ming the Tide: Predicting STEM attrition using student transcript data
Authors:
Lovenoor Aulck,
Rohan Aras,
Lysia Li,
Coulter L'Heureux,
Peter Lu,
Jevin West
Abstract:
Science, technology, engineering, and math (STEM) fields play growing roles in national and international economies by driving innovation and generating high salary jobs. Yet, the US is lagging behind other highly industrialized nations in terms of STEM education and training. Furthermore, many economic forecasts predict a rising shortage of domestic STEM-trained professions in the US for years to…
▽ More
Science, technology, engineering, and math (STEM) fields play growing roles in national and international economies by driving innovation and generating high salary jobs. Yet, the US is lagging behind other highly industrialized nations in terms of STEM education and training. Furthermore, many economic forecasts predict a rising shortage of domestic STEM-trained professions in the US for years to come. One potential solution to this deficit is to decrease the rates at which students leave STEM-related fields in higher education, as currently over half of all students intending to graduate with a STEM degree eventually attrite. However, little quantitative research at scale has looked at causes of STEM attrition, let alone the use of machine learning to examine how well this phenomenon can be predicted. In this paper, we detail our efforts to model and predict dropout from STEM fields using one of the largest known datasets used for research on students at a traditional campus setting. Our results suggest that attrition from STEM fields can be accurately predicted with data that is routinely collected at universities using only information on students' first academic year. We also propose a method to model student STEM intentions for each academic term to better understand the timing of STEM attrition events. We believe these results show great promise in using machine learning to improve STEM retention in traditional and non-traditional campus settings.
△ Less
Submitted 28 August, 2017;
originally announced August 2017.
-
Complexity Measures of Music
Authors:
April Pease,
Korosh Mahmoodi,
Bruce J. West
Abstract:
We present a technique to search for the presence of crucial events in music, based on the analysis of the music volume. Earlier work on this issue was based on the assumption that crucial events correspond to the change of music notes, with the interesting result that the complexity index of the crucial events is mu ~ 2, which is the same inverse power-law index of the dynamics of the brain. The…
▽ More
We present a technique to search for the presence of crucial events in music, based on the analysis of the music volume. Earlier work on this issue was based on the assumption that crucial events correspond to the change of music notes, with the interesting result that the complexity index of the crucial events is mu ~ 2, which is the same inverse power-law index of the dynamics of the brain. The search technique analyzes music volume and confirms the results of the earlier work, thereby contributing to the explanation as to why the brain is sensitive to music, through the phenomenon of complexity matching. Complexity matching has recently been interpreted as the transfer of multifractality from one complex network to another. For this reason we also examine the mulifractality of music, with the observation that the multifractal spectrum of a computer performance is significantly narrower than the multifractal spectrum of a human performance of the same musical score. We conjecture that although crucial events are demonstrably important for information transmission, they alone are not suficient to define musicality, which is more adequately measured by the multifractality spectrum.
△ Less
Submitted 26 August, 2017;
originally announced August 2017.
-
Alpha Wavelet Power as a Biomarker of Antidepressant Treatment Response in Bipolar Depression
Authors:
Wojciech Jernajczyk,
Pawel Gosek,
Miroslaw Latka,
Klaudia Kozlowska,
Lukasz Swiecicki,
Bruce J. West
Abstract:
There is mounting evidence of a link between the properties of electroencephalograms (EEGs) of depressive patients and the outcome of pharmacotherapy. The goal of this study was to develop an EEG biomarker of antidepressant treatment response which would require only a single EEG measurement. We recorded resting, 21-channel EEG in 17 inpatients suffering from bipolar depression in eyes closed and…
▽ More
There is mounting evidence of a link between the properties of electroencephalograms (EEGs) of depressive patients and the outcome of pharmacotherapy. The goal of this study was to develop an EEG biomarker of antidepressant treatment response which would require only a single EEG measurement. We recorded resting, 21-channel EEG in 17 inpatients suffering from bipolar depression in eyes closed and eyes open conditions. The EEG measurement was performed at the end of the short washout period which followed previously unsuccessful pharmacotherapy. We calculated the normalized wavelet power of alpha rhythm using two referential montages and an average reference montage. In particular, in the occipital (O1, O2, Oz) channels the wavelet power of responders was up to 84% higher than that of nonresponders. Using a novel classification algorithm we were able to correctly predict the outcome of treatment with 90% sensitivity and 100% specificity. The proposed biomarker requires only a single EEG measurement and consequently is intrinsically different from biomarkers which exploit the changes in prefrontal EEG induced by pharmacotherapy over a given time.
△ Less
Submitted 16 February, 2017;
originally announced February 2017.
-
Asymmetry of short-term control of spatio-temporal gait parameters during treadmill walking
Authors:
Klaudia Kozlowska,
Miroslaw Latka,
Bruce J. West
Abstract:
Optimization of energy cost determines average values of spatio-temporal gait parameters such as step duration, step length or step speed. However, during walking, humans need to adapt these parameters at every step to respond to exogenous and/or endogenic perturbations. While some neurological mechanisms that trigger these responses are known, our understanding of the fundamental principles gover…
▽ More
Optimization of energy cost determines average values of spatio-temporal gait parameters such as step duration, step length or step speed. However, during walking, humans need to adapt these parameters at every step to respond to exogenous and/or endogenic perturbations. While some neurological mechanisms that trigger these responses are known, our understanding of the fundamental principles governing step-by-step adaptation remains elusive. We determined the gait parameters of 20 healthy subjects with right-foot preference during treadmill walking at speeds of 1.1, 1.4 and 1.7 m/s. We found that when the value of the gait parameter was conspicuously greater (smaller) than the mean value, it was either followed immediately by a smaller (greater) value of the contralateral leg (interleg control), or the deviation from the mean value decreased during the next movement of ipsilateral leg (intraleg control). The selection of step duration and the selection of step length during such transient control events were performed in unique ways. We quantified the symmetry of short-term control of gait parameters and observed the significant dominance of the right leg in short-term control of all three parameters at higher speeds (1.4 and 1.7 m/s).
△ Less
Submitted 31 January, 2017;
originally announced January 2017.
-
Men Set Their Own Cites High: Gender and Self-citation across Fields and over Time
Authors:
Molly M. King,
Carl T. Bergstrom,
Shelley J. Correll,
Jennifer Jacquet,
Jevin D. West
Abstract:
How common is self-citation in scholarly publication, and does the practice vary by gender? Using novel methods and a data set of 1.5 million research papers in the scholarly database JSTOR published between 1779 and 2011, the authors find that nearly 10 percent of references are self-citations by a paper's authors. The findings also show that between 1779 and 2011, men cited their own papers 56 p…
▽ More
How common is self-citation in scholarly publication, and does the practice vary by gender? Using novel methods and a data set of 1.5 million research papers in the scholarly database JSTOR published between 1779 and 2011, the authors find that nearly 10 percent of references are self-citations by a paper's authors. The findings also show that between 1779 and 2011, men cited their own papers 56 percent more than did women. In the last two decades of data, men self-cited 70 percent more than women. Women are also more than 10 percentage points more likely than men to not cite their own previous work at all. While these patterns could result from differences in the number of papers that men and women authors have published rather than gender-specific patterns of self-citation behavior, this gender gap in self-citation rates has remained stable over the last 50 years, despite increased representation of women in academia. The authors break down self-citation patterns by academic field and number of authors and comment on potential mechanisms behind these observations. These findings have important implications for scholarly visibility and cumulative advantage in academic careers.
△ Less
Submitted 12 December, 2017; v1 submitted 30 June, 2016;
originally announced July 2016.
-
The Value of Conflict in Stable Social Networks
Authors:
Pensri Pramukkul,
Adam Svenkeson,
Bruce J. West,
Paolo Grigolini
Abstract:
A cooperative network model of sociological interest is examined to determine the sensitivity of the global dynamics to having a fraction of the members behaving uncooperatively, that is, being in conflict with the majority. We study a condition where in the absence of these uncooperative individuals, the contrarians, the control parameter exceeds a critical value and the network is frozen in a st…
▽ More
A cooperative network model of sociological interest is examined to determine the sensitivity of the global dynamics to having a fraction of the members behaving uncooperatively, that is, being in conflict with the majority. We study a condition where in the absence of these uncooperative individuals, the contrarians, the control parameter exceeds a critical value and the network is frozen in a state of consensus. The network dynamics change with variations in the percentage of contrarians, resulting in a balance between the value of the control parameter and the percentage of those in conflict with the majority. We show that the transmission of information from a network $B$ to a network $A$, with a small fraction of lookout members in $A$ who adopt the behavior of $B$, becomes maximal when both networks are assigned the same critical percentage of contrarians.
△ Less
Submitted 13 June, 2015;
originally announced June 2015.
-
Building the next pyramid
Authors:
Joseph West,
Greg Gallagher,
Kevin Waters,
Stephen Ward,
Tia Ward
Abstract:
The results of experimental tests of a novel method for moving large (pyramid construction size) stone blocks by rolling them are presented. The method is implemented by tying 12 identical rods of appropriately chosen radius to the faces of the block forming a rough dodecagon prism. Experiments using a 1,000 kg block show that it can be moved across level open ground with a dynamic coefficient of…
▽ More
The results of experimental tests of a novel method for moving large (pyramid construction size) stone blocks by rolling them are presented. The method is implemented by tying 12 identical rods of appropriately chosen radius to the faces of the block forming a rough dodecagon prism. Experiments using a 1,000 kg block show that it can be moved across level open ground with a dynamic coefficient of friction of less than 0.06. This value is a factor of five lower than that obtained for dragging the block, and the best values reported for dragging by others, at 0.3. the results are more dramatic than those obtained on smaller scale experiments on a 29.6 kg block, also reported here. For full scale pyramid blocks, the wooden "rods" woudl need to be posts of order 30 cm in diameter, similar in size to those used as masts on ships in the Nile.
△ Less
Submitted 25 February, 2015;
originally announced February 2015.
-
Fractional calculus ties the microscopic and macroscopic scales of complex network dynamics
Authors:
Bruce J. West,
Malgorzata Turalska,
Paolo Grigolini
Abstract:
A two-state master equation based decision making model has been shown to generate phase transitions, to be topologically complex and to manifest temporal complexity through an inverse power-law probability distribution function in the switching times between the two critical states of consensus. These properties are entailed by the fundamental assumption that the network elements in the decision…
▽ More
A two-state master equation based decision making model has been shown to generate phase transitions, to be topologically complex and to manifest temporal complexity through an inverse power-law probability distribution function in the switching times between the two critical states of consensus. These properties are entailed by the fundamental assumption that the network elements in the decision making model imperfectly imitate one another. The process of subordination establishes that a single network element can be described by a fractional master equation whose analytic solution yields the observed inverse power-law probability distribution obtained by numerical integration of the two-state master equation to a high degree of accuracy.
△ Less
Submitted 10 September, 2014;
originally announced September 2014.
-
Imitation versus payoff - duality of the decision-making process demonstrates criticality and consensus formation
Authors:
Malgorzata Turalska,
Bruce J. West
Abstract:
We consider a dual model of decision making, in which an individual forms its opinion based on contrasting mechanisms of imitation and rational calculation. The decision making model (DMM) implements imitating behavior by means of a network of coupled two-state master equations that undergoes a phase transition at a critical value of a control parameter. The evolutionary spatial game (EGM), being…
▽ More
We consider a dual model of decision making, in which an individual forms its opinion based on contrasting mechanisms of imitation and rational calculation. The decision making model (DMM) implements imitating behavior by means of a network of coupled two-state master equations that undergoes a phase transition at a critical value of a control parameter. The evolutionary spatial game (EGM), being a generalization of the Prisoner's dilemma game, is used to determine in objective fashion the cooperative or anti-cooperative strategy adopted by individuals. Interactions between two sources of dynamics increases the domain of initial states attracted to phase transition dynamics beyond that of the DMM network in isolation. Additionally, on average the influence of the DMM on the game increases the final observed fraction of cooperators in the system.
△ Less
Submitted 10 September, 2014;
originally announced September 2014.
-
How they (should have) built the pyramids
Authors:
J. West,
G. Gallagher,
K. Waters
Abstract:
A novel method is proposed for moving large, pyramid construction size, stone blocks. The method is inspired by a well known introductory physics homework problem, and is implemented by tying 12 identical rods of appropriately chosen radius to the faces of the block. The rods form the corners and new faces that transform the square prism into a dodecagon which can then be moved more easily by roll…
▽ More
A novel method is proposed for moving large, pyramid construction size, stone blocks. The method is inspired by a well known introductory physics homework problem, and is implemented by tying 12 identical rods of appropriately chosen radius to the faces of the block. The rods form the corners and new faces that transform the square prism into a dodecagon which can then be moved more easily by rolling than by dragging. Experimental results are presented and compared to independent work by another group which utilized wooden attachments providing a cylindrical shape. It is found that a small scale stone block converted to dodecagons can be moved across level open ground with a dynamic coefficient of friction of the order 0.2. For full scale pyramid blocks, the wooden rods would need to be posts of order 30 cm in diameter, similar in size to those used as masts on ships in the Nile.
△ Less
Submitted 14 August, 2014;
originally announced August 2014.
-
Spyhoppers and Stargazers: Can Whales See the Stars?
Authors:
Michael J. West
Abstract:
In Moby Dick, Herman Melville wondered how - or what - whales see with eyes on opposite sides of their heads. "It is plain that he can never see an object which is exactly ahead... Is his brain so much more comprehensive, combining and subtle than man's that he can at the same moment of time attentively examine two distinct prospects, one on one side of him, and the other in an exactly opposite di…
▽ More
In Moby Dick, Herman Melville wondered how - or what - whales see with eyes on opposite sides of their heads. "It is plain that he can never see an object which is exactly ahead... Is his brain so much more comprehensive, combining and subtle than man's that he can at the same moment of time attentively examine two distinct prospects, one on one side of him, and the other in an exactly opposite direction?" he asked. It's a good question. But if Melville were alive today he might have pondered something perhaps even more intriguing: Can whales see the stars?
△ Less
Submitted 10 September, 2013;
originally announced September 2013.
-
Memory in network flows and its effects on spreading dynamics and community detection
Authors:
Martin Rosvall,
Alcides V. Esquivel,
Andrea Lancichinetti,
Jevin D. West,
Renaud Lambiotte
Abstract:
Random walks on networks is the standard tool for modelling spreading processes in social and biological systems. This first-order Markov approach is used in conventional community detection, ranking, and spreading analysis although it ignores a potentially important feature of the dynamics: where flow moves to may depend on where it comes from. Here we analyse pathways from different systems, and…
▽ More
Random walks on networks is the standard tool for modelling spreading processes in social and biological systems. This first-order Markov approach is used in conventional community detection, ranking, and spreading analysis although it ignores a potentially important feature of the dynamics: where flow moves to may depend on where it comes from. Here we analyse pathways from different systems, and while we only observe marginal consequences for disease spreading, we show that ignoring the effects of second-order Markov dynamics has important consequences for community detection, ranking, and information spreading. For example, capturing dynamics with a second-order Markov model allows us to reveal actual travel patterns in air traffic and to uncover multidisciplinary journals in scientific communication. These findings were achieved only by using more available data and making no additional assumptions, and therefore suggest that accounting for higher-order memory in network flows can help us better understand how real systems are organized and function.
△ Less
Submitted 12 August, 2014; v1 submitted 21 May, 2013;
originally announced May 2013.
-
Blackbody Friction: Analytic expressions for velocity and position
Authors:
Joseph West
Abstract:
The equations of motion are solved analytically for speed and position for the case of blackbody friction on objects traveling at relativistic speeds with respect to observers fixed in the frame of the blackbody radiation. Two model cases are considered, both assume an object which is a perfect absorber of light, and one of which also assumes that the object maintains a constant rest mass by emitt…
▽ More
The equations of motion are solved analytically for speed and position for the case of blackbody friction on objects traveling at relativistic speeds with respect to observers fixed in the frame of the blackbody radiation. Two model cases are considered, both assume an object which is a perfect absorber of light, and one of which also assumes that the object maintains a constant rest mass by emitting photons in a momentum neutral process. The maximum distances attained for a light sail and a micron-size sand grain are determined under the assumption that the relevant blackbody source is that of cosmic background radiation. At a temperature of 3000K, the temperature at the time the universe became transparent, cosmic background radiation would cause a significant decay of speed fluctuations in dust clouds.
△ Less
Submitted 8 November, 2012;
originally announced November 2012.
-
The role of gender in scholarly authorship
Authors:
Jevin D. West,
Jennifer Jacquet,
Molly M. King,
Shelley J. Correll,
Carl T. Bergstrom
Abstract:
Gender disparities appear to be decreasing in academia according to a number of metrics, such as grant funding, hiring, acceptance at scholarly journals, and productivity, and it might be tempting to think that gender inequity will soon be a problem of the past. However, a large-scale analysis based on over eight million papers across the natural sciences, social sciences, and humanities re- revea…
▽ More
Gender disparities appear to be decreasing in academia according to a number of metrics, such as grant funding, hiring, acceptance at scholarly journals, and productivity, and it might be tempting to think that gender inequity will soon be a problem of the past. However, a large-scale analysis based on over eight million papers across the natural sciences, social sciences, and humanities re- reveals a number of understated and persistent ways in which gender inequities remain. For instance, even where raw publication counts seem to be equal between genders, close inspection reveals that, in certain fields, men predominate in the prestigious first and last author positions. Moreover, women are significantly underrepresented as authors of single-authored papers. Academics should be aware of the subtle ways that gender disparities can appear in scholarly authorship.
△ Less
Submitted 7 November, 2012;
originally announced November 2012.
-
Role of Committed Minorities in Times of Crisis
Authors:
Malgorzata Turalska,
Bruce J. West,
Paolo Grigolini
Abstract:
We use a Cooperative Decision Making (CDM) model to study the effect of committed minorities on group behavior in time of crisis. The CDM model has been shown to generate consensus through a phase-transition process that at criticality establishes long-range correlations among the individuals within a model society. In a condition of high consensus, the correlation function vanishes, thereby makin…
▽ More
We use a Cooperative Decision Making (CDM) model to study the effect of committed minorities on group behavior in time of crisis. The CDM model has been shown to generate consensus through a phase-transition process that at criticality establishes long-range correlations among the individuals within a model society. In a condition of high consensus, the correlation function vanishes, thereby making the network recover the ordinary locality condition. However, this state is not permanent and times of crisis occur when there is an ambiguity concerning a given social issue. The correlation function within the cooperative system becomes similarly extended as it is observed at criticality. This combination of independence (free will) and long-range correlation makes it possible for very small but committed minorities to produce substantial changes in social consensus.
△ Less
Submitted 6 July, 2012;
originally announced July 2012.
-
Beyond the Death of Linear Response: 1/f optimal information transport
Authors:
Gerardo Aquino,
Mauro Bologna,
Paolo Grigolini,
Bruce J. West
Abstract:
Non-ergodic renewal processes have recently been shown by several authors to be insensitive to periodic perturbations, thereby apparently sanctioning the death of linear response, a building block of nonequilibrium statistical physics. We show that it is possible to go beyond the ``death of linear response" and establish a permanent correlation between an external stimulus and the response of a co…
▽ More
Non-ergodic renewal processes have recently been shown by several authors to be insensitive to periodic perturbations, thereby apparently sanctioning the death of linear response, a building block of nonequilibrium statistical physics. We show that it is possible to go beyond the ``death of linear response" and establish a permanent correlation between an external stimulus and the response of a complex network generating non-ergodic renewal processes, by taking as stimulus a similar non-ergodic process. The ideal condition of 1/f-noise corresponds to a singularity that is expected to be relevant in several experimental conditions.
△ Less
Submitted 17 July, 2010;
originally announced July 2010.
-
Rolling Release Siege Engines: Teaching an Old Machine a New Trick
Authors:
Joseph West,
Seth Ross,
James Flesher
Abstract:
The analysis of a new "rolling release" mechanism is presented for three different "siege engine" designs. The range obtained with the rolling release is compared to the ranges obtained using the simpler "cup" and the proven "sling" release mechanisms. It is found that the rolling release is a significant improvement over the cup release, but the range of the rolling release is still well short…
▽ More
The analysis of a new "rolling release" mechanism is presented for three different "siege engine" designs. The range obtained with the rolling release is compared to the ranges obtained using the simpler "cup" and the proven "sling" release mechanisms. It is found that the rolling release is a significant improvement over the cup release, but the range of the rolling release is still well short of that attained with the sling release. It is also noted that the familiar "spoon" projectile holder used on many siege engines operates in a rolling release manner, contrary to appearances and to some previous published studies. Instead, the release mechanism used in ancient engines and in modern kit-purchased replicas is the rolling release.
△ Less
Submitted 27 April, 2010;
originally announced May 2010.
-
Time and ensemble averaging in time series analysis
Authors:
Miroslaw Latka,
Massimiliano Ignaccolo,
Wojciech Jernajczyk,
Bruce J. West
Abstract:
In many applications expectation values are calculated by partitioning a single experimental time series into an ensemble of data segments of equal length. Such single trajectory ensemble (STE) is a counterpart to a multiple trajectory ensemble (MTE) used whenever independent measurements or realizations of a stochastic process are available. The equivalence of STE and MTE for stationary systems w…
▽ More
In many applications expectation values are calculated by partitioning a single experimental time series into an ensemble of data segments of equal length. Such single trajectory ensemble (STE) is a counterpart to a multiple trajectory ensemble (MTE) used whenever independent measurements or realizations of a stochastic process are available. The equivalence of STE and MTE for stationary systems was postulated by Wang and Uhlenbeck in their classic paper on Brownian motion (Rev. Mod. Phys. 17, 323 (1945)) but surprisingly has not yet been proved. Using the stationary and ergodic paradigm of statistical physics -- the Ornstein-Uhlenbeck (OU) Langevin equation, we revisit Wang and Uhlenbeck's postulate. In particular, we find that the variance of the solution of this equation is different for these two ensembles. While the variance calculated using the MTE quantifies the spreading of independent trajectories originating from the same initial point, the variance for STE measures the spreading of two correlated random walkers. Thus, STE and MTE refer to two completely different dynamical processes. Guided by this interpretation, we introduce a novel algorithm of partitioning a single trajectory into a phenomenological ensemble, which we name a threshold trajectory ensemble (TTE), that for an ergodic system is equivalent to MTE. We find that in the cohort of healthy volunteers, the ratio of STE and TTE asymptotic variances of stage 4 sleep electroencephalogram is equal to 1.96 \pm 0.04 which is in agreement with the theoretically predicted value of 2.
△ Less
Submitted 13 April, 2010;
originally announced April 2010.
-
Probability flux as a method for detecting scaling
Authors:
M. Ignaccolo,
P. Grigolini,
B. J. West
Abstract:
We introduce a new method for detecting scaling in time series. The method uses the properties of the probability flux for stochastic self-affine processes and is called the probability flux analysis (PFA). The advantages of this method are: 1) it is independent of the finiteness of the moments of the self-affine process; 2) it does not require a binning procedure for numerical evaluation of the t…
▽ More
We introduce a new method for detecting scaling in time series. The method uses the properties of the probability flux for stochastic self-affine processes and is called the probability flux analysis (PFA). The advantages of this method are: 1) it is independent of the finiteness of the moments of the self-affine process; 2) it does not require a binning procedure for numerical evaluation of the the probability density function. These properties make the method particularly efficient for heavy tailed distributions in which the variance is not finite, for example, in Levy alpha-stable processes. This utility is established using a comparison with the diffusion entropy (DE) method.
△ Less
Submitted 2 April, 2010;
originally announced April 2010.
-
The Failure of the Ergodic Assumption
Authors:
M. Ignaccolo,
M. Latka,
B. J. West
Abstract:
The well established procedure of constructing phenomenological ensemble from a single long time series is investigated. It is determined that a time series generated by a simple Uhlenbeck-Ornstein Langevin equation is mean ergodic. However the probability ensemble average yields a variance that is different from that determined using the phenomenological ensemble (time average). We conclude tha…
▽ More
The well established procedure of constructing phenomenological ensemble from a single long time series is investigated. It is determined that a time series generated by a simple Uhlenbeck-Ornstein Langevin equation is mean ergodic. However the probability ensemble average yields a variance that is different from that determined using the phenomenological ensemble (time average). We conclude that the latter ensemble is often neither stationary nor ergodic and consequently the probability ensemble averages can misrepresent the underlying dynamic process.
△ Less
Submitted 22 January, 2010;
originally announced January 2010.
-
Public Perception of Astronomers: Revered, Reviled and Ridiculed
Authors:
Michael J. West
Abstract:
Society's view of astronomers has changed over time and from culture to culture. This review discusses some of the many ways that astronomers have been perceived by their societies and suggests ways that astronomers can influence public perception of ourselves and our profession in the future.
Society's view of astronomers has changed over time and from culture to culture. This review discusses some of the many ways that astronomers have been perceived by their societies and suggests ways that astronomers can influence public perception of ourselves and our profession in the future.
△ Less
Submitted 28 May, 2009; v1 submitted 25 May, 2009;
originally announced May 2009.
-
The differential aging of inertial and non-inertial observers: The Eyewitness observations of a relativistic polygon traveler
Authors:
J. West
Abstract:
The time dilation of non-inertial travelers in circular and polygonal closed paths are well known. In both cases observers completing a round trip will age less than an observer at rest with respect to the circle / polygon. This rapid aging is contrary to the slow aging that would be documented by temporarily comoving inertial observers who might also move along each edge of the polygon. A detai…
▽ More
The time dilation of non-inertial travelers in circular and polygonal closed paths are well known. In both cases observers completing a round trip will age less than an observer at rest with respect to the circle / polygon. This rapid aging is contrary to the slow aging that would be documented by temporarily comoving inertial observers who might also move along each edge of the polygon. A detailed description of what the non-inertial observer might actually see that would explain their own slow aging is presented. It is also argued that the non-inertial observer will disagree with the comoving inertial observers about the size of the diameter of the polygon. In particular, the non-inertial traveler will observer a Lorentz-like contraction of distances measured perpendicular to the direction of relative motion. As a result, the traveler will obtain a value of for the ratio of circumference to diameter, in contradiction to recently published work that predicts a measured ratio less than and a non-Euclidean geometry for the polygon traveler.
△ Less
Submitted 12 May, 2009;
originally announced May 2009.
-
The Dynamics of EEG Entropy
Authors:
M. Ignaccolo,
M. Latka,
W. Jernajczyk,
P. Grigolini,
B. J. West
Abstract:
EEG time series are analyzed using the diffusion entropy method. The resulting EEG entropy manifests short-time scaling, asymptotic saturation and an attenuated alpha-rhythm modulation. These properties are faithfully modeled by a phenomenological Langevin equation interpreted within a neural network context.
EEG time series are analyzed using the diffusion entropy method. The resulting EEG entropy manifests short-time scaling, asymptotic saturation and an attenuated alpha-rhythm modulation. These properties are faithfully modeled by a phenomenological Langevin equation interpreted within a neural network context.
△ Less
Submitted 5 March, 2009;
originally announced March 2009.
-
Dynamics of EEG Entropy: beyond signal plus noise
Authors:
M. Ignaccolo,
M. Latka,
W. Jernajczyk,
P. Grigolini,
B. J. West
Abstract:
EEG time series are analyzed using the diffusion entropy method. The resulting EEG entropy manifests short-time scaling, asymptotic saturation and an attenuated alpha-rhythm modulation. These properties are faithfully modeled by a phenomenological Langevin equation interpreted within a neural network context. Detrended fluctuation analysis of the EEG data is compared with diffusion entropy analy…
▽ More
EEG time series are analyzed using the diffusion entropy method. The resulting EEG entropy manifests short-time scaling, asymptotic saturation and an attenuated alpha-rhythm modulation. These properties are faithfully modeled by a phenomenological Langevin equation interpreted within a neural network context. Detrended fluctuation analysis of the EEG data is compared with diffusion entropy analysis and is found to suppress certain important properties of the EEG time series.
△ Less
Submitted 6 February, 2009;
originally announced February 2009.
-
Differences in Impact Factor Across Fields and Over Time
Authors:
Benjamin M. Althouse,
Jevin D. West,
Theodore Bergstrom,
Carl T. Bergstrom
Abstract:
The bibliometric measure impact factor is a leading indicator of journal influence, and impact factors are routinely used in making decisions ranging from selecting journal subscriptions to allocating research funding to deciding tenure cases. Yet journal impact factors have increased gradually over time, and moreover impact factors vary widely across academic disciplines. Here we quantify infla…
▽ More
The bibliometric measure impact factor is a leading indicator of journal influence, and impact factors are routinely used in making decisions ranging from selecting journal subscriptions to allocating research funding to deciding tenure cases. Yet journal impact factors have increased gradually over time, and moreover impact factors vary widely across academic disciplines. Here we quantify inflation over time and differences across fields in impact factor scores and determine the sources of these differences. We find that the average number of citations in reference lists has increased gradually, and this is the predominant factor responsible for the inflation of impact factor scores over time. Field-specific variation in the fraction of citations to literature indexed by Thomson Scientific's Journal Citation Reports is the single greatest contributor to differences among the impact factors of journals in different fields. The growth rate of the scientific literature as a whole, and cross-field differences in net size and growth rate of individual fields, have had very little influence on impact factor inflation or on cross-field differences in impact factor.
△ Less
Submitted 18 April, 2008;
originally announced April 2008.
-
Multiscaling comparative analysis of time series and geophysical phenomena
Authors:
Nicola Scafetta,
Bruce J. West
Abstract:
Different methods are used to determine the scaling exponents associated with a time series describing a complex dynamical process, such as those observed in geophysical systems. Many of these methods are based on the numerical evaluation of the variance of a diffusion process whose step increments are generated by the data. An alternative method focuses on the direct evaluation of the scaling c…
▽ More
Different methods are used to determine the scaling exponents associated with a time series describing a complex dynamical process, such as those observed in geophysical systems. Many of these methods are based on the numerical evaluation of the variance of a diffusion process whose step increments are generated by the data. An alternative method focuses on the direct evaluation of the scaling coefficient of the Shannon entropy of the same diffusion distribution. The combined use of these methods can efficiently distinguish between fractal Gaussian and Lévy-walk time series and help to discern between alternative underling complex dynamics.
△ Less
Submitted 29 September, 2005; v1 submitted 29 September, 2005;
originally announced September 2005.