Feigenbaum graphs: a complex network perspective of chaos

The recently formulated theory of horizontal visibility graphs transforms time series into graphs and allows the possibility of studying dynamical systems through the characterization of their associated networks. This method leads to a natural graph-theoretical description of nonlinear systems with qualities in the spirit of symbolic dynamics. We support our claim via the case study of the period-doubling and band-splitting attractor cascades that characterize unimodal maps. We provide a universal analytical description of this classic scenario in terms of the horizontal visibility graphs associated with the dynamics within the attractors, that we call Feigenbaum graphs, independent of map nonlinearity or other particulars. We derive exact results for their degree distribution and related quantities, recast them in the context of the renormalization group and find that its fixed points coincide with those of network entropy optimization. Furthermore, we show that the network entropy mimics the Lyapunov exponent of the map independently of its sign, hinting at a Pesin-like relation equally valid out of chaos.Comment: Published in PLoS ONE (Sep 2011

The Backward Monte Carlo Method for Semiconductor Device Simulation

© The Author(s) 2018A backward Monte Carlo method for the numerical solution of the semiconductor Boltzmann equation is presented. The method is particularly suited to simulate rare events. The general theory of the backward Monte Carlo method is described, and several estimators for the contact current are derived from that theory. The transition probabilities for the construction of the backward trajectories are chosen so as to satisfy the principle of detailed balance. This property guarantees stability of the numerical method and allows for a clear physical interpretation of the estimators. A symmetric sampling method which generates wave vectors always in pairs symmetric to the origin can be shown to yield zero current exactly as thermal equilibrium is approached. The properties of the different estimators are evaluated by simulation of an n-channel MOSFET. Quantities varying over many orders of magnitude can be resolved with ease. Such quantities are the drain current in the sub-threshold region, the high-energy tail of the carrier distribution function, and the so-called acceleration integral which varies over 30 orders in the example shown.Austrian Research Promotion Agency (FFG)149215041