Source author record

Johannes Zierenberg

Johannes Zierenberg 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

13works
9topics
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

13 published item(s)

preprint2022arXiv

Autocorrelations from emergent bistability in homeostatic spiking neural networks on neuromorphic hardware

A unique feature of neuromorphic computing is that memory is an implicit part of processing through traces of past information in the system's collective dynamics. The extent of memory about past inputs is commonly quantified by the autocorrelation time of collective dynamics. Based on past experimental evidence, a potential explanation for the underlying autocorrelations are close-to-critical fluctuations. Here, we show for self-organized networks of excitatory and inhibitory leaky integrate-and-fire neurons that autocorrelations can originate from emergent bistability upon reducing external input strength. We identify the bistability as a fluctuation-induced stochastic switching between metastable active and quiescent states in the vicinity of a non-equilibrium phase transition. This bistability occurs for networks with fixed heterogeneous weights as a consequence of homeostatic self-organization during development. Specifically, in our experiments on neuromorphic hardware and in computer simulations, the emergent bistability gives rise to autocorrelation times exceeding 500 ms despite single-neuron timescales of only 20 ms. Our results provide the first verification of biologically compatible autocorrelation times in networks of leaky integrate-and-fire neurons, which here are not generated by close-to-critical fluctuations but by emergent bistability in homeostatically regulated networks. Our results thereby constitute a new, complementary mechanism for emergent autocorrelations in networks of spiking neurons, with implications for biological and artificial networks, and introduces the general paradigm of fluctuation-induced bistability for driven systems with absorbing states.

preprint2020arXiv

Characterizing spreading dynamics of subsampled systems with non-stationary external input

Many systems with propagation dynamics, such as spike propagation in neural networks and spreading of infectious diseases, can be approximated by autoregressive models. The estimation of model parameters can be complicated by the experimental limitation that one observes only a fraction of the system (subsampling) and potentially time-dependent parameters, leading to incorrect estimates. We show analytically how to overcome the subsampling bias when estimating the propagation rate for systems with certain non-stationary external input. This approach is readily applicable to trial-based experimental setups and seasonal fluctuations, as demonstrated on spike recordings from monkey prefrontal cortex and spreading of norovirus and measles.

preprint2020arXiv

The challenges of containing SARS-CoV-2 via test-trace-and-isolate

Without a cure, vaccine, or proven long-term immunity against SARS-CoV-2, test-trace-and-isolate (TTI) strategies present a promising tool to contain its spread. For any TTI strategy, however, mitigation is challenged by pre- and asymptomatic transmission, TTI-avoiders, and undetected spreaders, who strongly contribute to hidden infection chains. Here, we studied a semi-analytical model and identified two tipping points between controlled and uncontrolled spread: (1) the behavior-driven reproduction number of the hidden chains becomes too large to be compensated by the TTI capabilities, and (2) the number of new infections exceeds the tracing capacity. Both trigger a self-accelerating spread. We investigated how these tipping points depend on challenges like limited cooperation, missing contacts, and imperfect isolation. Our model results suggest that TTI alone is insufficient to contain an otherwise unhindered spread of SARS-CoV-2, implying that complementary measures like social distancing and improved hygiene remain necessary.

preprint2019arXiv

Description of spreading dynamics by microscopic network models and macroscopic branching processes can differ due to coalescence

Spreading processes are conventionally monitored on a macroscopic level by counting the number of incidences over time. The spreading process can then be modeled either on the microscopic level, assuming an underlying interaction network, or directly on the macroscopic level, assuming that microscopic contributions are negligible. The macroscopic characteristics of both descriptions are commonly assumed to be identical. In this work, we show that these characteristics of microscopic and macroscopic descriptions can be different due to coalescence, i.e., a node being activated at the same time by multiple sources. In particular, we consider a (microscopic) branching network (probabilistic cellular automaton) with annealed connectivity disorder, record the macroscopic activity, and then approximate this activity by a (macroscopic) branching process. In this framework, we analytically calculate the effect of coalescence on the collective dynamics. We show that coalescence leads to a universal non-linear scaling function for the conditional expectation value of successive network activity. This allows us to quantify the difference between the microscopic model parameter and established macroscopic estimates. To overcome this difference, we propose a non-linear estimator that correctly infers the model branching parameter for all system sizes.

preprint2019arXiv

Tailored ensembles of neural networks optimize sensitivity to stimulus statistics

The dynamic range of stimulus processing in living organisms is much larger than a single neural network can explain. For a generic, tunable spiking network we derive that while the dynamic range is maximal at criticality, the interval of discriminable intensities is very similar for any network tuning due to coalescence. Compensating coalescence enables adaptation of discriminable intervals. Thus, we can tailor an ensemble of networks optimized to the distribution of stimulus intensities, e.g., extending the dynamic range arbitrarily. We discuss potential applications in machine learning.

preprint2015arXiv

Exploring different regimes in finite-size scaling of the droplet condensation-evaporation transition

We present a finite-size scaling analysis of the droplet condensation-evaporation transition of a lattice gas (in two and three dimensions) and a Lennard-Jones gas (in three dimensions) at fixed density. Parallel multicanonical simulations allow sampling of the required system sizes with precise equilibrium estimates. In the limit of large systems, we verify the theoretical leading-order scaling prediction for both the transition temperature and the finite-size rounding. In addition, we present an emerging intermediate scaling regime, consistent in all considered cases and with similar recent observations for polymer aggregation. While the intermediate regime locally may show a different effective scaling, we show that it is a gradual crossover to the large-system scaling behavior by including empirical higher-order corrections. This implies that care has to be taken when considering scaling ranges, possibly leading to completely wrong predictions for the thermodynamic limit. In this study, we consider a crossing of the phase boundary orthogonal to the usual fixed temperature studies. We show that this is an equivalent approach and, under certain conditions, may show smaller finite-size corrections.

preprint2015arXiv

From amorphous aggregates to polymer bundles: The role of stiffness on structural phases in polymer aggregation

We study the aggregation transition of a finite theta-polymer system in dependence on the bending stiffness $κ$ with the help of parallel multicanonical simulations. In order to distinguish amorphous aggregates from polymer bundles we introduce an order parameter, measuring the correlation of the end-to-end vectors. With the help of this order parameter, we construct generic $T$-$κ$ phase diagrams for systems with $2$ and $8$ polymers and discuss the occurring phases from amorphous aggregates to bundle structures. For an intermediate stiffness range we find multiple aggregated phases which change with increasing number of polymers and discuss their nature with the help of microcanonical analyses. We show that the stiffness of semiflexible theta polymers is the distinguishing parameter for the emergent structural motifs.

preprint2015arXiv

Parallel multicanonical study of the three-dimensional Blume-Capel model

We study the thermodynamic properties of the three-dimensional Blume-Capel model on the simple cubic lattice by means of computer simulations. In particular, we implement a parallelized variant of the multicanonical approach and perform simulations by keeping a constant temperature and crossing the phase boundary along the crystal-field axis. We obtain numerical data for several temperatures in both the first- and second-order regime of the model. Finite-size scaling analyses provide us with transition points and the dimensional scaling behavior in the numerically demanding first-order regime, as well as a clear verification of the expected Ising universality in the respective second-order regime. Finally, we discuss the scaling behavior in the vicinity of the tricritical point.

preprint2014arXiv

Aggregation of theta-polymers in spherical confinement

We investigate the aggregation transition of theta polymers in spherical confinement with multicanonical simulations. This allows for a systematic study of the effect of density on the aggregation transition temperature for up to 24 monodisperse polymers. Our results for solutions in the dilute regime show that polymers can be considered isolated for all temperatures larger than the aggregation temperature, which is shown to be a function of the density. The resulting competition between single-polymer collapse and aggregation yields the lower temperature bound of the isolated chain approximation. We provide entropic and energetic arguments to describe the density dependence and finite-size effects of the aggregation transition for monodisperse solutions in finite systems. This allows us to estimate the aggregation transition temperature of dilute systems in a spherical cavity, using a few simulations of small, sufficiently dilute polymer systems.

preprint2014arXiv

Application of the parallel multicanonical method to lattice gas condensation

We present the speedup from a novel parallel implementation of the multicanonical method on the example of a lattice gas in two and three dimensions. In this approach, all cores perform independent equilibrium runs with identical weights, collecting their sampled histograms after each iteration in order to estimate consecutive weights. The weights are then redistributed to all cores. These steps are repeated until the weights are converged. This procedure benefits from a minimum of communication while distributing the necessary amount of statistics efficiently. Using this method allows us to study a broad temperature range for a variety of large and complex systems. Here, a gas is modeled as particles on the lattice, which interact only with their nearest neighbors. For a fixed density this model is equivalent to the Ising model with fixed magnetization. We compare our results to an analytic prediction for equilibrium droplet formation, confirming that a single macroscopic droplet forms only above a critical density.

preprint2013arXiv

Scaling Properties of Parallelized Multicanonical Simulations

We implemented a parallel version of the multicanonical algorithm and applied it to a variety of systems with phase transitions of first and second order. The parallelization relies on independent equilibrium simulations that only communicate when the multicanonical weight function is updated. That way, the Markov chains efficiently sample the temporary distributions allowing for good estimations of consecutive weight functions. The systems investigated range from the well known Ising and Potts spin systems to bead-spring polymers. We estimate the speedup with increasing number of parallel processes. Overall, the parallelization is shown to scale quite well. In the case of multicanonical simulations of the $q$-state Potts model ($q\ge6$) and multimagnetic simulations of the Ising model, the optimal performance is limited due to emerging barriers.

preprint2012arXiv

Influence of lattice disorder on the structure of persistent polymer chains

We study the static properties of a semiflexible polymer exposed to a quenched random environment by means of computer simulations. The polymer is modeled as two-dimensional Heisenberg chain. For the random environment we consider hard disks arranged on a square lattice. We apply an off-lattice growth algorithm as well as the multicanonical Monte Carlo method to investigate the influence of both disorder occupation probability and polymer stiffness on the equilibrium properties of the polymer. We show that the additional length scale induced by the stiffness of the polymer extends the well-known phenomenology considerably. The polymer's response to the disorder is either contraction or extension depending on the ratio of polymer stiffness and void space extension. Additionally, the periodic structure of the lattice is reflected in the observables that characterize the polymer.

preprint2012arXiv

Simulating flexible polymers in a potential of randomly distributed hard disks

We perform equilibrium computer simulations of a two-dimensional pinned flexible polymer exposed to a quenched disorder potential consisting of hard disks. We are especially interested in the high-density regime of the disorder, where subtle structures such as cavities and channels play a central role. We apply an off-lattice growth algorithm proposed by Garel and Orland [J. Phys. A 23, L621 (1990)], where a distribution of polymers is constructed in parallel by growing each of them monomer by monomer. In addition we use a multicanonical Monte Carlo method in order to cross-check the results of the growth algorithm. We measure the end-to-end distribution and the tangent-tangent correlations. We also investigate the scaling behavior of the mean square end-to-end distance in dependence of the monomer number. While the influence of the potential in the low-density case is merely marginal, it dominates the configurational properties of the polymer for high densities.