Source author record

Yueheng Lan

Yueheng Lan 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

15works
19topics
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

15 published item(s)

preprint2022arXiv

A reduced variational approach for searching cycles in high-dimensional systems

Searching recurrent patterns in complex systems with high-dimensional phase spaces is an important task in diverse fields. In the current work, an improved scheme is proposed to accelerate the recently designed variational approach for finding periodic orbits in systems with chaotic dynamics based on the existence of inertial manifold widely observed in various spatially extended systems, especially those with high dimensions. On the premise of keeping exponential convergence of the variational method, an effective loop evolution equation is derived to greatly reduce the storage and computing time. With repeated modification of local coordinates and evolution of the guess loop being carried out alternately, the rapid convergence and the stability of the reduction scheme are effectively achieved. The dimension of local coordinate subspaces is generally larger than the number of nonnegative Lyapunov exponents to ensure the exponential convergence. The proposed scheme is successfully demonstrated on several well-known examples and expected to supply a powerful tool in the exploration of high-dimensional nonlinear systems.

preprint2022arXiv

Wielding Intermittency with Cycle Expansions

As periodic orbit theory works badly on computing the observable averages of dynamical systems with intermittency, we propose a scheme to cooperate with cycle expansion and perturbation theory so that we can deal with intermittent systems and compute the averages more precisely. Periodic orbit theory assumes that the shortest unstable periodic orbits build the framework of the system and provides cycles expansion to compute dynamical quantities based on them, while the perturbation theory can locally analyze the structure of dynamical systems. The dynamical averages may be obtained more precisely by combining the two techniques together. Based on the integrability near the marginal orbits and the hyperbolicity in the part away from the singularities in intermittent systems, the chief idea of this paper is to revise intermittent maps and maintain the natural measure produced by the original maps. We get the natural measure near the singularity through the Taylor expansions and periodic orbit theory captures the natural measure in the other parts of the phase space. We try this method on 1-dimensional intermittent maps with single singularity, and more precise results are achieved.

preprint2021arXiv

Criticality in Reservoir Computer of Coupled Phase Oscillators

Accumulating evidences show that the cerebral cortex is operating near a critical state featured by power-law size distribution of neural avalanche activities, yet evidence of this critical state in artificial neural networks mimicking the cerebral cortex is lacking. Here we design an artificial neural network of coupled phase oscillators and, by the technique of reservoir computing in machine learning, train it for predicting chaos. It is found that when the machine is properly trained, oscillators in the reservoir are synchronized into clusters whose sizes follow a power-law distribution. This feature, however, is absent when the machine is poorly trained. Additionally, it is found that despite the synchronization degree of the original network, once properly trained, the reservoir network is always developed to the same critical state, exemplifying the "attractor" nature of this state in machine learning. The generality of the results is verified in different reservoir models and by different target systems, and it is found that the scaling exponent of the distribution is independent on the reservoir details and the bifurcation parameter of the target system, but is modified when the dynamics of the target system is changed to a different type. The findings shed lights on the nature of machine learning, and are helpful to the design of high-performance machine in physical systems.

preprint2020arXiv

Koopman analysis in oscillator synchronization

Synchronization is an important dynamical phenomenon in coupled nonlinear systems, which has been studied extensively in recent years. However, analysis focused on individual orbits seems hard to extend to complex systems while a global statistical approach is overly cursory. Koopman operator technique seems to well balance the two approaches. In this paper, we extend Koopman analysis to the study of synchronization of coupled oscillators by extracting important eigenvalues and eigenfunctions from the observed time series. A renormalization group analysis is designed to derive an analytic approximation of the eigenfunction in case of weak coupling that dominates the oscillation. For moderate or strong couplings, numerical computation further confirms the importance of the average frequencies and the associated eigenfunctions. The synchronization transition points could be located with quite high accuracy by checking the correlation of neighbouring eigenfunctions at different coupling strengths, which is readily applied to other nonlinear systems.

preprint2020arXiv

Koopman Operator and Phase Space Partition of Chaotic Maps

Koopman operator describes evolution of observables in the phase space, which could be used to extract characteristic dynamical features of a nonlinear system. Here, we show that it is possible to carry out interesting symbolic partitions based on properly constructed eigenfunctions of the operator for chaotic maps. The partition boundaries are the extrema of these eigenfunctions, the accuracy of which is improved by including more basis functions in the numerical computation. The validity of this scheme is demonstrated in well-known 1-d and 2-d maps. It seems no obstacle to extend the computation to nonlinear systems of high dimensions, which provides a possible way of dissecting complex dynamics.

preprint2016arXiv

Nonlinear waves on circle networks with excitable nodes

Nonlinear wave formation and propagation on a complex network with excitable node dynamics is of fundamental interest in diverse fields in science and engineering. Here, we propose a new model of the Kuramoto type to study nonlinear wave generation and propagation on circular subgraphs of a complex network. On circle networks, in the continuum limit, this model is equivalent to the over-damped Frenkel-Kontorova model. The new model is shown to keep the essential features of those well-known models such as the diffusively coupled Bär-Eiswirth model but with much simplified expression such that analytic analysis becomes possible. We classify traveling wave solutions on circle networks and show the universality of its features with perturbation analysis and numerical computation.

preprint2016arXiv

Stochastic Thermodynamics of a Particle in a Box

The piston system (particles in a box) is the simplest and paradigmatic model in traditional thermodynamics. However, the recently established framework of stochastic thermodynamics (ST) fails to apply to this model system due to the embedded singularity in the potential. In this Letter we study the stochastic thermodynamics of a particle in a box by adopting a novel coordinate transformation technique. Through comparing with the exact solution of a breathing harmonic oscillator, we obtain analytical results of work distribution for an arbitrary protocol in the linear response regime, and verify various predictions of the Fluctuation-Dissipation Relation. When applying to the Brownian Szilard's engine model, we obtain the optimal protocol $λ_t = λ_0 2^{t/τ}$ for a given sufficiently long total time $τ$. Our study not only establishes a paradigm for studying ST of a particle in a box, but also bridges the long-standing gap in the development of ST.

preprint2015arXiv

Energy dissipation in an adaptive molecular circuit

The ability to monitor nutrient and other environmental conditions with high sensitivity is crucial for cell growth and survival. Sensory adaptation allows a cell to recover its sensitivity after a transient response to a shift in the strength of extracellular stimulus. The working principles of adaptation have been established previously based on rate equations which do not consider fluctuations in a thermal environment. Recently, G. Lan et al. (Nature Phys., 8:422-8, 2012) performed a detailed analysis of a stochastic model for the E. coli sensory network. They showed that accurate adaptation is possible only when the system operates in a nonequilibrium steady-state (NESS). They further proposed an energy-speed-accuracy (ESA) trade-off relation. We present here analytic results on the NESS of the model through a mapping to a one-dimensional birth-death process. An exact expression for the entropy production rate is also derived. Based on these results, we are able to discuss the ESA relation in a more general setting. Our study suggests that the adaptation error can be reduced exponentially as the methylation range increases. Finally, we show that a nonequilibrium phase transition exists in the infinite methylation range limit, despite the fact that the model contains only two discrete variables.

preprint2014arXiv

SEK: Sparsity exploiting $k$-mer-based estimation of bacterial community composition

Motivation: Estimation of bacterial community composition from a high-throughput sequenced sample is an important task in metagenomics applications. Since the sample sequence data typically harbors reads of variable lengths and different levels of biological and technical noise, accurate statistical analysis of such data is challenging. Currently popular estimation methods are typically very time consuming in a desktop computing environment. Results: Using sparsity enforcing methods from the general sparse signal processing field (such as compressed sensing), we derive a solution to the community composition estimation problem by a simultaneous assignment of all sample reads to a pre-processed reference database. A general statistical model based on kernel density estimation techniques is introduced for the assignment task and the model solution is obtained using convex optimization tools. Further, we design a greedy algorithm solution for a fast solution. Our approach offers a reasonably fast community composition estimation method which is shown to be more robust to input data variation than a recently introduced related method. Availability: A platform-independent Matlab implementation of the method is freely available at http://www.ee.kth.se/ctsoftware; source code that does not require access to Matlab is currently being tested and will be made available later through the above website.

preprint2014arXiv

The stochastic thermodynamics of a rotating Brownian particle in a gradient flow

We compute the entropy production engendered in the environment from a single Brownian particle which moves in a mean flow, and show that it corresponds in expectation to classical near-equilibrium entropy production in the surrounding fluid with specific mesoscopic transport coefficients. With temperature gradient, extra terms are found which results from the nonlinear interaction between the particle and the non-equilibrated environment. The calculations are carried out in the multi-time formalism and in the advection-dissipation limit where the Stokes number (St) of the flow tends to zero and the Peclet number (Pe) diverges but the combination St times Pe remains constant.

preprint2013arXiv

Improved contact prediction in proteins: Using pseudolikelihoods to infer Potts models

Spatially proximate amino acids in a protein tend to coevolve. A protein's three-dimensional (3D) structure hence leaves an echo of correlations in the evolutionary record. Reverse engineering 3D structures from such correlations is an open problem in structural biology, pursued with increasing vigor as more and more protein sequences continue to fill the data banks. Within this task lies a statistical inference problem, rooted in the following: correlation between two sites in a protein sequence can arise from firsthand interaction but can also be network-propagated via intermediate sites; observed correlation is not enough to guarantee proximity. To separate direct from indirect interactions is an instance of the general problem of inverse statistical mechanics, where the task is to learn model parameters (fields, couplings) from observables (magnetizations, correlations, samples) in large systems. In the context of protein sequences, the approach has been referred to as direct-coupling analysis. Here we show that the pseudolikelihood method, applied to 21-state Potts models describing the statistical properties of families of evolutionarily related proteins, significantly outperforms existing approaches to the direct-coupling analysis, the latter being based on standard mean-field techniques. This improved performance also relies on a modified score for the coupling strength. The results are verified using known crystal structures of specific sequence instances of various protein families. Code implementing the new method can be found at http://plmdca.csc.kth.se/.

preprint2011arXiv

A resolution of the turbulence paradox: numerical implementation

Sommerfeld paradox (turbulence paradox) roughly says that mathematically the Couette linear shear flow is linearly stable for all values of the Reynolds number, but experimentally transition from the linear shear to turbulence occurs under perturbations of any size when the Reynolds number is large enough. In [Li, Lin 2011], we offered a resolution of this paradox. The aim of this paper is to provide a numerical implementation of the resolution. The main idea of the resolution is that even though the linear shear is linearly stable, slow orbits (also called quasi-steady states) in arbitrarily small neighborhoods of the linear shear can be linearly unstable. The key is that in infinite dimensions, smallness in one norm does not mean smallness in all norms. Our study here focuses upon a sequence of 2D oscillatory shears which are the Couette linear shear plus small amplitude and high frequency sinusoidal shear perturbations. In the fluid velocity variable, the sequence approaches the Couette linear shear (e.g. in $L^2$ norm of velocity), thus it can be viewed as Couette linear shear plus small noises; while in the fluid vorticity variable, the sequence does not approaches the Couette linear shear (e.g. in $H^1$ norm of velocity). Unlike the Couette linear shear, the sequence of oscillatory shears has inviscid linear instability; furthermore, with the sequence of oscillatory shears as potentials, the Orr-Sommerfeld operator has unstable eigenvalues when the Reynolds number is large enough, this should lead to transient nonlinear growth which manifests as transient turbulence as observed in experiments. The main result of this paper verifies this transient growth.

preprint2011arXiv

Accelerating cycle expansions by dynamical conjugacy

Periodic orbit theory provides two important functions---the dynamical zeta function and the spectral determinant for the calculation of dynamical averages in a nonlinear system. Their cycle expansions converge rapidly when the system is uniformly hyperbolic but greatly slowed down in the presence of non-hyperbolicity. We find that the slow convergence can be associated with singularities in the natural measure. A properly designed coordinate transformation may remove these singularities and results in a dynamically conjugate system where fast convergence is restored. The technique is successfully demonstrated on several examples of one-dimensional maps and some remaining challenges are discussed.

preprint2005arXiv

Experimental evidence for the Volume Conjecture for the simplest hyperbolic non-2-bridge knot

Loosely speaking, the Volume Conjecture states that the limit of the n-th colored Jones polynomial of a hyperbolic knot, evaluated at the primitive complex n-th root of unity is a sequence of complex numbers that grows exponentially. Moreover, the exponential growth rate is proportional to the hyperbolic volume of the knot. We provide an efficient formula for the colored Jones function of the simplest hyperbolic non-2-bridge knot, and using this formula, we provide numerical evidence for the Hyperbolic Volume Conjecture for the simplest hyperbolic non-2-bridge knot.