-
Predicting unobserved climate time series data at distant areas via spatial correlation using reservoir computing
Authors:
Shihori Koyama,
Daisuke Inoue,
Hiroaki Yoshida,
Kazuyuki Aihara,
Gouhei Tanaka
Abstract:
Collecting time series data spatially distributed in many locations is often important for analyzing climate change and its impacts on ecosystems. However, comprehensive spatial data collection is not always feasible, requiring us to predict climate variables at some locations. This study focuses on a prediction of climatic elements, specifically near-surface temperature and pressure, at a target…
▽ More
Collecting time series data spatially distributed in many locations is often important for analyzing climate change and its impacts on ecosystems. However, comprehensive spatial data collection is not always feasible, requiring us to predict climate variables at some locations. This study focuses on a prediction of climatic elements, specifically near-surface temperature and pressure, at a target location apart from a data observation point. Our approach uses two prediction methods: reservoir computing (RC), known as a machine learning framework with low computational requirements, and vector autoregression models (VAR), recognized as a statistical method for analyzing time series data. Our results show that the accuracy of the predictions degrades with the distance between the observation and target locations. We quantitatively estimate the distance in which effective predictions are possible. We also find that in the context of climate data, a geographical distance is associated with data correlation, and a strong data correlation significantly improves the prediction accuracy with RC. In particular, RC outperforms VAR in predicting highly correlated data within the predictive range. These findings suggest that machine learning-based methods can be used more effectively to predict climatic elements in remote locations by assessing the distance to them from the data observation point in advance. Our study on low-cost and accurate prediction of climate variables has significant value for climate change strategies.
△ Less
Submitted 5 June, 2024;
originally announced June 2024.
-
The Architecture of Sponge Choanocyte Chambers Maximizes Mechanical Pumping Efficiency
Authors:
Takumi Ogawa,
Shuji Koyama,
Toshihiro Omori,
Kenji Kikuchi,
Helene de Maleprade,
Raymond E. Goldstein,
Takuji Ishikawa
Abstract:
Sponges, the basalmost members of the animal kingdom, exhibit a range of complex architectures in which microfluidic channels connect multitudes of spherical chambers lined with choanocytes, flagellated filter-feeding cells. Choanocyte chambers can possess scores or even hundreds of such cells, which drive complex flows entering through porous walls and exiting into the sponge channels. One of the…
▽ More
Sponges, the basalmost members of the animal kingdom, exhibit a range of complex architectures in which microfluidic channels connect multitudes of spherical chambers lined with choanocytes, flagellated filter-feeding cells. Choanocyte chambers can possess scores or even hundreds of such cells, which drive complex flows entering through porous walls and exiting into the sponge channels. One of the mysteries of the choanocyte chamber is its spherical shape, as it seems inappropriate for inducing directional transport since many choanocyte flagella beat in opposition to such a flow. Here we combine direct imaging of choanocyte chambers in living sponges with computational studies of many-flagella models to understand the connection between chamber architecture and directional flow. We find that those flagella that beat against the flow play a key role in raising the pressure inside the choanocyte chamber, with the result that the mechanical pumping efficiency, calculated from the pressure rise and flow rate, reaches a maximum at a small outlet opening angle. Comparison between experimental observations and the results of numerical simulations reveal that the chamber diameter, flagellar wave number and the outlet opening angle of the freshwater sponge $E. muelleri$, as well as several other species, are related in a manner that maximizes the mechanical pumping efficiency. These results indicate the subtle balances at play during morphogenesis of choanocyte chambers, and give insights into the physiology and body design of sponges.
△ Less
Submitted 22 February, 2024;
originally announced February 2024.
-
Mechanical Optimization of Skateboard Pumping
Authors:
Florian Kogelbauer,
Shinsuke Koyama,
Daniel E. Callan,
Shigeru Shinomoto
Abstract:
Skateboarders perform a reciprocating motion on a curved ramp, called pumping, by moving their bodies up and down perpendicular to the ramp's surface. We propose a simple mechanical model for this pumping motion and solve the equation of motion explicitly in angular coordinates. This allows us to derive an optimal control strategy to maximize amplitude by dynamically adjusting the center of mass o…
▽ More
Skateboarders perform a reciprocating motion on a curved ramp, called pumping, by moving their bodies up and down perpendicular to the ramp's surface. We propose a simple mechanical model for this pumping motion and solve the equation of motion explicitly in angular coordinates. This allows us to derive an optimal control strategy to maximize amplitude by dynamically adjusting the center of mass of the skateboarder. This optimal strategy is compared to experimental results for the motion of a skilled and an unskilled skateboarder in a half-pipe, validating that a skilled skateboarder follows the optimal control strategy more closely.
△ Less
Submitted 16 March, 2024; v1 submitted 29 December, 2023;
originally announced December 2023.
-
Dynamical cooperation model for mitigating the segregation phase in Schelling's model
Authors:
Akihisa Okada,
Daisuke Inoue,
Shihori Koyama,
Tadayoshi Matsumori,
Hiroaki Yoshida
Abstract:
We consider a Schelling-like segregation model, in which the behavior of individual agents is determined by a mixed individual and global utility. With a high ratio of global utility being incorporated, the agents are cooperative in order to realize a homogenized state, otherwise the agents are less cooperative, leading to an undesired Nash equilibrium with low utility. In the present study, we in…
▽ More
We consider a Schelling-like segregation model, in which the behavior of individual agents is determined by a mixed individual and global utility. With a high ratio of global utility being incorporated, the agents are cooperative in order to realize a homogenized state, otherwise the agents are less cooperative, leading to an undesired Nash equilibrium with low utility. In the present study, we introduce a dynamically varying cooperation degree parameter to prevent the agents from falling into such a low-utility equilibrium state. More precisely, a large cooperation degree is assigned when the agents are in high-utility regions, whereas agents having low utility behave more individually. Simulation results show that homogenized phases with globally high utility are achieved with the present dynamical control, even for the case of a low mean value of cooperation degree. Since the cooperation degree represents the magnitude with which Pigouvian tax is enforced in the model of residential movement within a city, this result suggests the possibility of tax intervention to circumvent the undesired segregation of residents.
△ Less
Submitted 17 October, 2022;
originally announced October 2022.
-
Anomalous Diffusion in a Monolayer of Lightweight Spheres Fluidized in Airflow
Authors:
Shihori Koyama,
Tomoki Matsuno,
Takashi Noguchi
Abstract:
This paper presents statistical analyses of random motions in a single layer of fluidized lightweight spherical particles. Foam polystyrene spheres were driven by an upward airflow through the sieve mesh, and their two-dimensional motion was acquired using image analysis. In the bulk region, the particle velocity distributions changed from Gaussian to heavy-tailed distribution as the bulk packing…
▽ More
This paper presents statistical analyses of random motions in a single layer of fluidized lightweight spherical particles. Foam polystyrene spheres were driven by an upward airflow through the sieve mesh, and their two-dimensional motion was acquired using image analysis. In the bulk region, the particle velocity distributions changed from Gaussian to heavy-tailed distribution as the bulk packing fraction $φ_b$ was increased. The mean square displacement of the particles exhibited transition to subdiffusion at much lower $φ_b$ than observed in previous studies using similar setup but with heavier particles. A slight superdiffusion and significant growth of the correlation length in the two-body velocity correlation was observed at further large $φ_b$. The effect of the wall on the dynamics of the particles was also investigated and the anisotropy of the granular temperature was found to be a useful index to discriminate between the wall region and the bulk. The turbulence statistics in the wake of a particle indicated a strong wall-normal asymmetry of aerodynamic forcing as the ``thermal'' agitation in the wall region.
△ Less
Submitted 13 October, 2021;
originally announced October 2021.
-
Electro-osmotic diode based on colloidal nano-valves between double membranes
Authors:
Shihori Koyama,
Daisuke Inoue,
Akihisa Okada,
Hiroaki Yoshida
Abstract:
The rectification of electro-osmotic flows is important in micro/nano fluidics applications such as micro-pumps and energy conversion devices. Here, we propose a simple electro-osmotic diode in which colloidal particles are contained between two parallel membranes with different pore densities. While the flow in the forward direction just pushes the colloidal particles toward the high-pore-density…
▽ More
The rectification of electro-osmotic flows is important in micro/nano fluidics applications such as micro-pumps and energy conversion devices. Here, we propose a simple electro-osmotic diode in which colloidal particles are contained between two parallel membranes with different pore densities. While the flow in the forward direction just pushes the colloidal particles toward the high-pore-density membrane, the backward flow is blocked by the particles near the low-pore-density membrane, which clog the pores. Nonequilibrium molecular dynamics simulations show a strong nonlinear dependence on the electric field for both the electric current and electro-osmotic flow, indicating diode characteristics. A mathematical model to reproduce the electro-osmotic diode behavior is constructed, introducing an effective pore diameter as a model for pores clogged by the colloidal particles. Good agreement is obtained between the proposed model with estimated parameter values and the results of direct molecular dynamics simulations. The proposed electro-osmotic diode has potential application in downsized microfluidic pumps, e.g., the pump induced under AC electric fields.
△ Less
Submitted 13 September, 2021; v1 submitted 9 September, 2021;
originally announced September 2021.
-
Membranes for spontaneous separation of pedestrian counter flows
Authors:
Shihori Koyama,
Daisuke Inoue,
Akihisa Okada,
Hiroaki Yoshida
Abstract:
Designing efficient traffic lanes for pedestrians is a critical aspect of urban planning as walking remains the most common form of mobility among the increasingly diverse methods of transportation. Herein, we investigate pedestrian counter flows in a straight corridor, in which two groups of people are walking in opposite directions. We demonstrate, using a molecular dynamics approach applying th…
▽ More
Designing efficient traffic lanes for pedestrians is a critical aspect of urban planning as walking remains the most common form of mobility among the increasingly diverse methods of transportation. Herein, we investigate pedestrian counter flows in a straight corridor, in which two groups of people are walking in opposite directions. We demonstrate, using a molecular dynamics approach applying the social force model, that a simple array of obstacles improves flow rates by producing flow separations even in crowded situations. We also report on a developed model describing the separation behavior that regards an array of obstacles as a membrane and induces spontaneous separation of pedestrians groups. When appropriately designed, those obstacles are fully capable of controlling the filtering direction so that pedestrians tend to keep moving to their left (or right) spontaneously. These results have the potential to provide useful guidelines for industrial designs aimed at improving ubiquitous human mobility.
△ Less
Submitted 30 March, 2020; v1 submitted 5 March, 2020;
originally announced March 2020.
-
The statistical physics of discovering exogenous and endogenous factors in a chain of events
Authors:
Shinsuke Koyama,
Shigeru Shinomoto
Abstract:
Event occurrence is not only subject to the environmental changes, but is also facilitated by the events that have occurred in a system. Here, we develop a method for estimating such extrinsic and intrinsic factors from a single series of event-occurrence times. The analysis is performed using a model that combines the inhomogeneous Poisson process and the Hawkes process, which represent exogenous…
▽ More
Event occurrence is not only subject to the environmental changes, but is also facilitated by the events that have occurred in a system. Here, we develop a method for estimating such extrinsic and intrinsic factors from a single series of event-occurrence times. The analysis is performed using a model that combines the inhomogeneous Poisson process and the Hawkes process, which represent exogenous fluctuations and endogenous chain-reaction mechanisms, respectively. The model is fit to a given dataset by minimizing the free energy, for which statistical physics and a path-integral method are utilized. Because the process of event occurrence is stochastic, parameter estimation is inevitably accompanied by errors, and it can ultimately occur that exogenous and endogenous factors cannot be captured even with the best estimator. We obtained four regimes categorized according to whether respective factors are detected. By applying the analytical method to real time series of debate in a social-networking service, we have observed that the estimated exogenous and endogenous factors are close to the first comments and the follow-up comments, respectively. This method is general and applicable to a variety of data, and we have provided an application program, by which anyone can analyze any series of event times.
△ Less
Submitted 1 March, 2020;
originally announced March 2020.
-
Parabolic Jets from the Spinning Black Hole in M87
Authors:
Masanori Nakamura,
Keiichi Asada,
Kazuhiro Hada,
Hung-Yi Pu,
Scott Noble,
Chihyin Tseng,
Kenji Toma,
Motoki Kino,
Hiroshi Nagai,
Kazuya Takahashi,
Juan-Carlos Algaba,
Monica Orienti,
Kazunori Akiyama,
Akihiro Doi,
Gabriele Giovannini,
Marcello Giroletti,
Mareki Honma,
Shoko Koyama,
Rocco Lico,
Kotaro Niinuma,
Fumie Tazaki
Abstract:
The M87 jet is extensively examined by utilizing general relativistic magnetohydrodynamic (GRMHD) simulations as well as the steady axisymmetric force-free electrodynamic (FFE) solution. Quasi-steady funnel jets are obtained in GRMHD simulations up to the scale of $\sim 100$ gravitational radius ($r_{\rm g}$) for various black hole (BH) spins. As is known, the funnel edge is approximately determin…
▽ More
The M87 jet is extensively examined by utilizing general relativistic magnetohydrodynamic (GRMHD) simulations as well as the steady axisymmetric force-free electrodynamic (FFE) solution. Quasi-steady funnel jets are obtained in GRMHD simulations up to the scale of $\sim 100$ gravitational radius ($r_{\rm g}$) for various black hole (BH) spins. As is known, the funnel edge is approximately determined by the following equipartitions; i) the magnetic and rest-mass energy densities and ii) the gas and magnetic pressures. Our numerical results give an additional factor that they follow the outermost parabolic streamline of the FFE solution, which is anchored to the event horizon on the equatorial plane. We also identify the matter dominated, non-relativistic corona/wind play a dynamical role in shaping the funnel jet into the parabolic geometry. We confirm a quantitative overlap between the outermost parabolic streamline of the FFE jet and the edge of jet sheath in VLBI observations at $\sim 10^{1}$-$10^{5} \, r_{\rm g}$, suggesting that the M87 jet is likely powered by the spinning BH. Our GRMHD simulations also indicate a lateral stratification of the bulk acceleration (i.e., the spine-sheath structure) as well as an emergence of knotty superluminal features. The spin characterizes the location of the jet stagnation surface inside the funnel. We suggest that the limb-brightened feature could be associated with the nature of the BH-driven jet, if the Doppler beaming is a dominant factor. Our findings can be examined with (sub-)mm VLBI observations, giving a clue for the origin of the M87 jet.
△ Less
Submitted 23 October, 2018;
originally announced October 2018.
-
Identifying exogenous and endogenous activity in social media
Authors:
Kazuki Fujita,
Alexey Medvedev,
Shinsuke Koyama,
Renaud Lambiotte,
Shigeru Shinomoto
Abstract:
The occurrence of new events in a system is typically driven by external causes and by previous events taking place inside the system. This is a general statement, applying to a range of situations including, more recently, to the activity of users in Online social networks (OSNs). Here we develop a method for extracting from a series of posting times the relative contributions of exogenous, e.g.…
▽ More
The occurrence of new events in a system is typically driven by external causes and by previous events taking place inside the system. This is a general statement, applying to a range of situations including, more recently, to the activity of users in Online social networks (OSNs). Here we develop a method for extracting from a series of posting times the relative contributions of exogenous, e.g. news media, and endogenous, e.g. information cascade. The method is based on the fitting of a generalized linear model (GLM) equipped with a self-excitation mechanism. We test the method with synthetic data generated by a nonlinear Hawkes process, and apply it to a real time series of tweets with a given hashtag. In the empirical dataset, the estimated contributions of exogenous and endogenous volumes are close to the amounts of original tweets and retweets respectively. We conclude by discussing the possible applications of the method, for instance in online marketing.
△ Less
Submitted 2 August, 2018;
originally announced August 2018.
-
Fluctuation scaling in neural spike trains
Authors:
Shinsuke Koyama,
Ryota Kobayashi
Abstract:
Fluctuation scaling has been observed universally in a wide variety of phenomena. In time series that describe sequences of events, fluctuation scaling is expressed as power function relationships between the mean and variance of either inter-event intervals or counting statistics, depending on measurement variables. In this article, fluctuation scaling has been formulated for a series of events i…
▽ More
Fluctuation scaling has been observed universally in a wide variety of phenomena. In time series that describe sequences of events, fluctuation scaling is expressed as power function relationships between the mean and variance of either inter-event intervals or counting statistics, depending on measurement variables. In this article, fluctuation scaling has been formulated for a series of events in which scaling laws in the inter-event intervals and counting statistics were related. We have considered the first-passage time of an Ornstein-Uhlenbeck process and used a conductance-based neuron model with excitatory and inhibitory synaptic inputs to demonstrate the emergence of fluctuation scaling with various exponents, depending on the input regimes and the ratio between excitation and inhibition. Furthermore, we have discussed the possible implication of these results in the context of neural coding.
△ Less
Submitted 2 November, 2015; v1 submitted 23 September, 2014;
originally announced September 2014.
-
Development of the Photomultiplier-Tube Readout System for the CTA Large Size Telescope
Authors:
H. Kubo,
R. Paoletti,
Y. Awane,
A. Bamba,
M. Barcelo,
J. A. Barrio,
O. Blanch,
J. Boix,
C. Delgado,
D. Fink,
D. Gascon,
S. Gunji,
R. Hagiwara,
Y. Hanabata,
K. Hatanaka,
M. Hayashida,
M. Ikeno,
S. Kabuki,
H. Katagiri,
J. Kataoka,
Y. Konno,
S. Koyama,
T. Kishimoto,
J. Kushida,
G. Martinez
, et al. (29 additional authors not shown)
Abstract:
We have developed a prototype of the photomultiplier tube (PMT) readout system for the Cherenkov Telescope Array (CTA) Large Size Telescope (LST). Two thousand PMTs along with their readout systems are arranged on the focal plane of each telescope, with one readout system per 7-PMT cluster. The Cherenkov light pulses generated by the air showers are detected by the PMTs and amplified in a compact,…
▽ More
We have developed a prototype of the photomultiplier tube (PMT) readout system for the Cherenkov Telescope Array (CTA) Large Size Telescope (LST). Two thousand PMTs along with their readout systems are arranged on the focal plane of each telescope, with one readout system per 7-PMT cluster. The Cherenkov light pulses generated by the air showers are detected by the PMTs and amplified in a compact, low noise and wide dynamic range gain block. The output of this block is then digitized at a sampling rate of the order of GHz using the Domino Ring Sampler DRS4, an analog memory ASIC developed at Paul Scherrer Institute. The sampler has 1,024 capacitors per channel and four channels are cascaded for increased depth. After a trigger is generated in the system, the charges stored in the capacitors are digitized by an external slow sampling ADC and then transmitted via Gigabit Ethernet. An onboard FPGA controls the DRS4, trigger threshold, and Ethernet transfer. In addition, the control and monitoring of the Cockcroft-Walton circuit that provides high voltage for the 7-PMT cluster are performed by the same FPGA. A prototype named Dragon has been developed that has successfully sampled PMT signals at a rate of 2 GHz, and generated single photoelectron spectra.
△ Less
Submitted 12 July, 2013;
originally announced July 2013.
-
Temporal fluctuation scaling in nonstationary counting processes
Authors:
Shinsuke Koyama
Abstract:
The fluctuation scaling law has universally been observed in a wide variety of phenomena. For counting processes describing the number of events occurred during time intervals, it is expressed as a power function relationship between the variance and the mean of the event count per unit time, the characteristic exponent of which is obtained theoretically in the limit of long duration of counting w…
▽ More
The fluctuation scaling law has universally been observed in a wide variety of phenomena. For counting processes describing the number of events occurred during time intervals, it is expressed as a power function relationship between the variance and the mean of the event count per unit time, the characteristic exponent of which is obtained theoretically in the limit of long duration of counting windows. Here I show that the scaling law effectively appears even in a short timescale in which only a few events occur. Consequently, the counting statistics of nonstationary event sequences are shown to exhibit the scaling law as well as the dynamics at temporal resolution of this timescale. I also propose a method to extract in a systematic manner the characteristic scaling exponent from nonstationary data.
△ Less
Submitted 27 June, 2013;
originally announced June 2013.
-
Approximate Methods for State-Space Models
Authors:
Shinsuke Koyama,
Lucia Castellanos Pérez-Bolde,
Cosma Rohilla Shalizi,
Robert E. Kass
Abstract:
State-space models provide an important body of techniques for analyzing time-series, but their use requires estimating unobserved states. The optimal estimate of the state is its conditional expectation given the observation histories, and computing this expectation is hard when there are nonlinearities. Existing filtering methods, including sequential Monte Carlo, tend to be either inaccurate…
▽ More
State-space models provide an important body of techniques for analyzing time-series, but their use requires estimating unobserved states. The optimal estimate of the state is its conditional expectation given the observation histories, and computing this expectation is hard when there are nonlinearities. Existing filtering methods, including sequential Monte Carlo, tend to be either inaccurate or slow. In this paper, we study a nonlinear filter for nonlinear/non-Gaussian state-space models, which uses Laplace's method, an asymptotic series expansion, to approximate the state's conditional mean and variance, together with a Gaussian conditional distribution. This {\em Laplace-Gaussian filter} (LGF) gives fast, recursive, deterministic state estimates, with an error which is set by the stochastic characteristics of the model and is, we show, stable over time. We illustrate the estimation ability of the LGF by applying it to the problem of neural decoding and compare it to sequential Monte Carlo both in simulations and with real data. We find that the LGF can deliver superior results in a small fraction of the computing time.
△ Less
Submitted 20 April, 2010;
originally announced April 2010.