Source author record

Kazuyuki Aihara

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
21topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

22 published item(s)

preprint2022arXiv

Reservoir Computing with Diverse Timescales for Prediction of Multiscale Dynamics

Machine learning approaches have recently been leveraged as a substitute or an aid for physical/mathematical modeling approaches to dynamical systems. To develop an efficient machine learning method dedicated to modeling and prediction of multiscale dynamics, we propose a reservoir computing (RC) model with diverse timescales by using a recurrent network of heterogeneous leaky integrator (LI) neurons. We evaluate computational performance of the proposed model in two time series prediction tasks related to four chaotic fast-slow dynamical systems. In a one-step-ahead prediction task where input data are provided only from the fast subsystem, we show that the proposed model yields better performance than the standard RC model with identical LI neurons. Our analysis reveals that the timescale required for producing each component of target multiscale dynamics is appropriately and flexibly selected from the reservoir dynamics by model training. In a long-term prediction task, we demonstrate that a closed-loop version of the proposed model can achieve longer-term predictions compared to the counterpart with identical LI neurons depending on the hyperparameter setting.

preprint2021arXiv

Analytical estimation of maximum fraction of infected individuals with one-shot non-pharmaceutical intervention in a hybrid epidemic model

Facing a global epidemic of new infectious diseases such as COVID-19, non-pharmaceutical interventions (NPIs), which reduce transmission rates without medical actions, are being implemented around the world to mitigate spreads. One of the problems in assessing the effects of NPIs is that different NPIs have been implemented at different times based on the situation of each country; therefore, few assumptions can be shared about how the introduction of policies affects the patient population. Mathematical models can contribute to further understanding these phenomena by obtaining analytical solutions as well as numerical simulations. In this study, an NPI was introduced into the SIR model for a conceptual study of infectious diseases under the condition that the transmission rate was reduced to a fixed value only once within a finite time duration, and its effect was analyzed numerically and theoretically. It was analytically shown that the maximum fraction of infected individuals and the final size could be larger if the intervention starts too early. The analytical results also suggested that more individuals may be infected at the peak of the second wave with a stronger intervention. This study provides quantitative relationship between the strength of a one-shot intervention and the reduction in the number of patients with no approximation. This suggests the importance of the strength and time of NPIs, although detailed studies are necessary for the implementation of NPIs in complicated real-world environments as the model used in this study is based on various simplifications.

preprint2021arXiv

Scaling advantage of nonrelaxational dynamics for high-performance combinatorial optimization

The development of physical simulators, called Ising machines, that sample from low energy states of the Ising Hamiltonian has the potential to drastically transform our ability to understand and control complex systems. However, most of the physical implementations of such machines have been based on a similar concept that is closely related to relaxational dynamics such as in simulated, mean-field, chaotic, and quantum annealing. We show that nonrelaxational dynamics that is associated with broken detailed balance and positive entropy production rate can accelerate the sampling of low energy states compared to that of conventional methods. By implementing such dynamics on field programmable gate array, we show that the nonrelaxational dynamics that we propose, called chaotic amplitude control, exhibits a scaling with problem size of the time to finding optimal solutions and its variance that is significantly smaller than that of relaxational schemes recently implemented on Ising machines.

preprint2021arXiv

Traffic Signal Optimization on a Square Lattice with Quantum Annealing

The spread of intelligent transportation systems in urban cities has caused heavy computational loads, requiring a novel architecture for managing large-scale traffic. In this study, we develop a method for globally controlling traffic signals arranged on a square lattice by means of a quantum annealing machine, namely the D-Wave quantum annealer. We first formulate a signal optimization problem that minimizes the imbalance of traffic flows in two orthogonal directions. Then we reformulate this problem as an Ising Hamiltonian, which is fully compatible with quantum annealers. The new control method is compared with a conventional local control method for a large 50-by-50 city, and the results exhibit the superiority of our global control method in suppressing traffic imbalance over wide parameter ranges. Furthermore, the solutions to the global control method obtained with the quantum annealing machine are better than those obtained with conventional simulated annealing. In addition, we prove analytically that the local and the global control methods converge at the limit where cars have equal probabilities for turning and going straight. These results are verified with numerical experiments.

preprint2020arXiv

A Supervised Learning Algorithm for Multilayer Spiking Neural Networks Based on Temporal Coding Toward Energy-Efficient VLSI Processor Design

Spiking neural networks (SNNs) are brain-inspired mathematical models with the ability to process information in the form of spikes. SNNs are expected to provide not only new machine-learning algorithms, but also energy-efficient computational models when implemented in VLSI circuits. In this paper, we propose a novel supervised learning algorithm for SNNs based on temporal coding. A spiking neuron in this algorithm is designed to facilitate analog VLSI implementations with analog resistive memory, by which ultra-high energy efficiency can be achieved. We also propose several techniques to improve the performance on a recognition task, and show that the classification accuracy of the proposed algorithm is as high as that of the state-of-the-art temporal coding SNN algorithms on the MNIST dataset. Finally, we discuss the robustness of the proposed SNNs against variations that arise from the device manufacturing process and are unavoidable in analog VLSI implementation. We also propose a technique to suppress the effects of variations in the manufacturing process on the recognition performance.

preprint2020arXiv

Fully Neural Network based Model for General Temporal Point Processes

A temporal point process is a mathematical model for a time series of discrete events, which covers various applications. Recently, recurrent neural network (RNN) based models have been developed for point processes and have been found effective. RNN based models usually assume a specific functional form for the time course of the intensity function of a point process (e.g., exponentially decreasing or increasing with the time since the most recent event). However, such an assumption can restrict the expressive power of the model. We herein propose a novel RNN based model in which the time course of the intensity function is represented in a general manner. In our approach, we first model the integral of the intensity function using a feedforward neural network and then obtain the intensity function as its derivative. This approach enables us to both obtain a flexible model of the intensity function and exactly evaluate the log-likelihood function, which contains the integral of the intensity function, without any numerical approximations. Our model achieves competitive or superior performances compared to the previous state-of-the-art methods for both synthetic and real datasets.

preprint2020arXiv

Time-domain digital-to-analog converter for spiking neural network hardware

We propose a new digital-to-analog converter (DAC) for realizing a synapse circuit of mixed-signal spiking neural networks. We named this circuit "time-domain DAC (TDAC)". This produces weights for converting a digital input code into voltage using one current waveform. Therefore, a TDAC is more compact than a conventional DAC that comprises many current sources and resistors. Moreover, a TDAC with leak resistance reproduces biological plausible synaptic responses that are expressed as alpha functions or dual exponential equations. We will show numerical analysis results of a TDAC and circuit simulation results of a circuit designed by the TSMC 40 nm CMOS process.

preprint2019arXiv

Quantum Expectation-Maximization Algorithm

Clustering algorithms are a cornerstone of machine learning applications. Recently, a quantum algorithm for clustering based on the k-means algorithm has been proposed by Kerenidis, Landman, Luongo and Prakash. Based on their work, we propose a quantum expectation-maximization (EM) algorithm for Gaussian mixture models (GMMs). The robustness and quantum speedup of the algorithm is demonstrated. We also show numerically the advantage of GMM over k-means for non-trivial cluster data.

preprint2015arXiv

Bifurcations and strange nonchaotic attractors in a phase oscillator model of glacial-interglacial cycles

Glacial-interglacial cycles are large variations in continental ice mass and greenhouse gases, which have dominated climate variability over the Quaternary. The dominant periodicity of the cycles is $\sim $40 kyr before the so-called middle Pleistocene transition between $\sim$1.2 and $\sim$0.7 Myr ago, and it is $\sim $100 kyr after the transition. In this paper, the dynamics of glacial-interglacial cycles are investigated using a phase oscillator model forced by the time-varying incoming solar radiation (insolation). We analyze the bifurcations of the system and show that strange nonchaotic attractors appear through nonsmooth saddle-node bifurcations of tori. The bifurcation analysis indicates that mode-locking is likely to occur for the 41 kyr glacial cycles but not likely for the 100 kyr glacial cycles. The sequence of mode-locked 41 kyr cycles is robust to small parameter changes. However, the sequence of 100 kyr glacial cycles can be sensitive to parameter changes when the system has a strange nonchaotic attractor.

preprint2015arXiv

Contribution of hidden modes to nonlinear epidemic dynamics in urban human proximity networks

Recently developed techniques to acquire high-quality human mobility data allow large-scale simulations of the spread of infectious diseases with high spatial and temporal resolution.Analysis of such data has revealed the oversimplification of existing theoretical frameworks to infer the final epidemic size or influential nodes from the network topology. Here we propose a spectral decomposition-based framework for the quantitative analysis of epidemic processes on realistic networks of human proximity derived from urban mobility data. Common wisdom suggests that modes with larger eigenvalues contribute more to the epidemic dynamics. However, we show that hidden dominant structures, namely modes with smaller eigenvalues but a greater contribution to the epidemic dynamics, exist in the proximity network. This framework provides a basic understanding of the relationship between urban human motion and epidemic dynamics, and will contribute to strategic mitigation policy decisions.

preprint2015arXiv

Feedback Regulation and its Efficiency in Biochemical Networks

Intracellular biochemical networks fluctuate dynamically due to various internal and external sources of fluctuation. Dissecting the fluctuation into biologically relevant components is important for understanding how a cell controls and harnesses noise and how information is transferred over apparently noisy intracellular networks. While substantial theoretical and experimental advancement on the decomposition of fluctuation was achieved for feedforward networks without any loop, we still lack a theoretical basis that can consistently extend such advancement to feedback networks. The main obstacle that hampers is the circulative propagation of fluctuation by feedback loops. In order to define the relevant quantity for the impact of feedback loops for fluctuation, disentanglement of the causally interlocked influence between the components is required. In addition, we also lack an approach that enables us to infer non-perturbatively the influence of the feedback to fluctuation as the dual reporter system does in the feedforward network. In this work, we resolve these problems by extending the work on the fluctuation decomposition and the dual reporter system. For a single-loop feedback network with two components, we define feedback loop gain as the feedback efficiency that is consistent with the fluctuation decomposition for feedforward networks. Then, we clarify the relation of the feedback efficiency with the fluctuation propagation in an open-looped FF network. Finally, by extending the dual reporter system, we propose a conjugate feedback and feedforward system for estimating the feedback efficiency only from the statistics of the system non-perturbatively.

preprint2015arXiv

Smoothing effect for spatially distributed renewable resources and its impact on power grid robustness

In this paper, we show that spatial correlation of renewable energy outputs greatly influences the robustness of power grids. First, we propose a new index for the spatial correlation among renewable energy outputs. We find that the spatial correlation of renewable energy outputs in a short time-scale is as weak as that caused by independent random variables and that in a long time-scale is as strong as that under perfect synchronization. Then, by employing the topology of the power grid in eastern Japan, we analyze the robustness of the power grid with spatial correlation of renewable energy outputs. The analysis is performed by using a realistic differential-algebraic equations model and the result shows that the spatial correlation of the energy resources strongly degrades the robustness of the power grid. Our result suggests that the spatial correlation of the renewable energy outputs should be taken into account when estimating the stability of power grids.

preprint2015arXiv

Steady state and mean recurrence time for random walks on stochastic temporal networks

Random walks are basic diffusion processes on networks and have applications in, for example, searching, navigation, ranking, and community detection. Recent recognition of the importance of temporal aspects on networks spurred studies of random walks on temporal networks. Here we theoretically study two types of event-driven random walks on a stochastic temporal network model that produces arbitrary distributions of interevent-times. In the so-called active random walk, the interevent-time is reinitialized on all links upon each movement of the walker. In the so-called passive random walk, the interevent-time is only reinitialized on the link that has been used last time, and it is a type of correlated random walk. We find that the steady state is always the uniform density for the passive random walk. In contrast, for the active random walk, it increases or decreases with the node's degree depending on the distribution of interevent-times. The mean recurrence time of a node is inversely proportional to the degree for both active and passive random walks. Furthermore, the mean recurrence time does or does not depend on the distribution of interevent-times for the active and passive random walks, respectively.

preprint2014arXiv

Approximated Infomax Early Stopping: Revisiting Gaussian RBMs on Natural Images

We pursue an early stopping technique that helps Gaussian Restricted Boltzmann Machines (GRBMs) to gain good natural image representations in terms of overcompleteness and data fitting. GRBMs are widely considered as an unsuitable model for natural images because they gain non-overcomplete representations which include uniform filters that do not represent useful image features. We have recently found that GRBMs once gain and subsequently lose useful filters during their training, contrary to this common perspective. We attribute this phenomenon to a tradeoff between overcompleteness of GRBM representations and data fitting. To gain GRBM representations that are overcomplete and fit data well, we propose a measure for GRBM representation quality, approximated mutual information, and an early stopping technique based on this measure. The proposed method boosts performance of classifiers trained on GRBM representations.

preprint2014arXiv

Understanding migraine using dynamical network biomarkers

Background: Mathematical modeling approaches are becoming ever more established in clinical neuroscience. They provide insight that is key to understand complex interactions of network phenomena, in general, and interactions within the migraine generator network, in particular. Purpose: In this study, two recent modeling studies on migraine are set in the context of premonitory symptoms that are easy to confuse for trigger factors. This causality confusion is explained, if migraine attacks are initiated by a transition caused by a tipping point. Conclusion: We need to characterize the involved neuronal and autonomic subnetworks and their connections during all parts of the migraine cycle if we are ever to understand migraine. We predict that mathematical models have the potential to dismantle large and correlated fluctuations in such subnetworks as a dynamical network biomarker of migraine.

preprint2013arXiv

Auto-pooling: Learning to Improve Invariance of Image Features from Image Sequences

Learning invariant representations from images is one of the hardest challenges facing computer vision. Spatial pooling is widely used to create invariance to spatial shifting, but it is restricted to convolutional models. In this paper, we propose a novel pooling method that can learn soft clustering of features from image sequences. It is trained to improve the temporal coherence of features, while keeping the information loss at minimum. Our method does not use spatial information, so it can be used with non-convolutional models too. Experiments on images extracted from natural videos showed that our method can cluster similar features together. When trained by convolutional features, auto-pooling outperformed traditional spatial pooling on an image classification task, even though it does not use the spatial topology of features.

preprint2013arXiv

Dynamics between order and chaos in conceptual models of glacial cycles

The dynamics of glacial cycles is studied in terms of the dynamical systems theory. We explore the dependence of the climate state on the phase of astronomical forcing by examining five conceptual models of glacial cycles proposed in the literature. The models can be expressed as quasiperiodically forced dynamical systems. It is shown that four of them exhibit a strange nonchaotic attractor (SNA), which is an intermediate regime between quasiperiodicity and chaos. Then, the dependence of the climate state on the phase of astronomical forcing is not given by smooth relations, but constitutes a geometrically strange set. Our result suggests that the dynamics of SNA is a candidate for the dynamics of glacial cycles, in addition to the well-known dynamics of quasiperiodicity and chaos.

preprint2013arXiv

Non-standard scaling law of fluctuations in finite-size systems of globally coupled oscillators

Universal scaling laws form one of the central issues in physics. A non-standard scaling law or a breakdown of a standard scaling law, on the other hand, can often lead to the finding of a new universality class in physical systems. Recently, we found that a statistical quantity related to fluctuations follows a non-standard scaling law with respect to system size in a synchronized state of globally coupled non-identical phase oscillators [Nishikawa et al., Chaos $\boldsymbol{22}$, 013133 (2012)]. However, it is still unclear how widely this non-standard scaling law is observed. In the present paper, we discuss the conditions required for the unusual scaling law in globally coupled oscillator systems, and we validate the conditions by numerical simulations of several different models.

preprint2013arXiv

Towards dynamical network biomarkers in neuromodulation of episodic migraine

Computational methods have complemented experimental and clinical neursciences and led to improvements in our understanding of the nervous systems in health and disease. In parallel, neuromodulation in form of electric and magnetic stimulation is gaining increasing acceptance in chronic and intractable diseases. In this paper, we firstly explore the relevant state of the art in fusion of both developments towards translational computational neuroscience. Then, we propose a strategy to employ the new theoretical concept of dynamical network biomarkers (DNB) in episodic manifestations of chronic disorders. In particular, as a first example, we introduce the use of computational models in migraine and illustrate on the basis of this example the potential of DNB as early-warning signals for neuromodulation in episodic migraine.

preprint2012arXiv

Chaotic Phase Synchronization in Bursting-neuron Models Driven by a Weak Periodic Force

We investigate the entrainment of a neuron model exhibiting a chaotic spiking-bursting behavior in response to a weak periodic force. This model exhibits two types of oscillations with different characteristic time scales, namely, long and short time scales. Several types of phase synchronization are observed, such as 1 : 1 phase locking between a single spike and one period of the force and 1 : l phase locking between the period of slow oscillation underlying bursts and l periods of the force. Moreover, spiking-bursting oscillations with chaotic firing patterns can be synchronized with the periodic force. Such a type of phase synchronization is detected from the position of a set of points on a unit circle, which is determined by the phase of the periodic force at each spiking time. We show that this detection method is effective for a system with multiple time scales. Owing to the existence of both the short and the long time scales, two characteristic phenomena are found around the transition point to chaotic phase synchronization. One phenomenon shows that the average time interval between successive phase slips exhibits a power-law scaling against the driving force strength and that the scaling exponent has an unsmooth dependence on the changes in the driving force strength. The other phenomenon shows that Kuramoto's order parameter before the transition exhibits stepwise behavior as a function of the driving force strength, contrary to the smooth transition in a model with a single time scale.

preprint2012arXiv

Long-term fluctuations in globally coupled phase oscillators with general coupling: Finite size effects

We investigate the diffusion coefficient of the time integral of the Kuramoto order parameter in globally coupled nonidentical phase oscillators. This coefficient represents the deviation of the time integral of the order parameter from its mean value on the sample average. In other words, this coefficient characterizes long-term fluctuations of the order parameter. For a system of N coupled oscillators, we introduce a statistical quantity D, which denotes the product of N and the diffusion coefficient. We study the scaling law of D with respect to the system size N. In other well-known models such as the Ising model, the scaling property of D is D \sim O(1) for both coherent and incoherent regimes except for the transition point. In contrast, in the globally coupled phase oscillators, the scaling law of D is different for the coherent and incoherent regimes: D \sim O(1/N^a) with a certain constant a>0 in the coherent regime and D \sim O(1) in the incoherent regime. We demonstrate that these scaling laws hold for several representative coupling schemes.

preprint2006arXiv

Pulse Dynamics in Coupled Excitable FIbers: Soliton-like Collision, Recombination, and Overtaking

We study the dynamics of a reaction-diffusion system composed of two mutually coupled excitable fibers. We focus on the situation in which dynamical properties of the two fibers are not identical because of the parameter difference between the fibers. Using the spatially one-dimensional FitzHugh-Nagumo equations as a model of a single excitable fiber, we show that the system exhibits a rich variety of dynamical behavior, including soliton-like collision between two pulses, recombination of a solitary pulse and synchronized pulses, and overtaking of a slow-moving solitary pulse by fast-moving synchronized pulses.