Source author record

Barbara Drossel

Barbara Drossel 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

28works
18topics
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

28 published item(s)

preprint2022arXiv

Explicit expressions for stationary states of the Lindblad equation for a finite state space

The Lindblad equation describes the time evolution of a density matrix of a quantum mechanical system. Stationary solutions are obtained by time-averaging the solution, which will in general depend on the initial state. We provide an analytical expression for the steady states of the Lindblad equation using the quantum jump unraveling, a version of an ergodic theorem, and the stationary probabilities of the corresponding discret-time Markov chains. Our result is valid when the number of states appearing the in quantum trajectory is finite. The classical case of a Markov jump-process is recovered as a special case, and differences between the two are discussed.

preprint2022arXiv

Master Stability Functions for metacommunities with two types of habitats

Current questions in ecology revolve around instabilities in the dynamics on spatial networks and particularly the effect of node heterogeneity. We extend the Master Stability Function formalism to inhomogeneous biregular networks having two types of spatial nodes. Notably, this class of systems also allows the investigation of certain types of dynamics on higher-order networks. Combined with the Generalized Modelling approach to study the linear stability of steady states, this is a powerful tool to numerically asses the stability of large ensembles of systems. We analyze the stability of ecological metacommunities with two distinct types of habitats analytically and numerically in order to identify several sets of conditions under which the dynamics can become stabilized by dispersal. Our analytical approach allows general insights into stabilizing and destabilizing effects in metapopulations. Specifically, we show that maladaptive dispersal may be stable under certain conditions.

preprint2021arXiv

Obtaining the long-term behavior of Master equations with finite state space from the structure of the associated state transition network

The Master equation describes the time evolution of the probabilities of a system with a discrete state space. This time evolution approaches for long times a stationary state that will in general depend on the initial probability distribution. Conditions under which the stationary state is unique are usually given as remarks appended to more comprehensive theories in the mathematical literature. We provide a direct and complete derivation of a necessary and sufficient criterion for when this steady state is unique. We translate this problem into the language of graph theory and show that there is a one-to-one correspondence between minimal absorbing sets within the state-transition network and linearly independent stationary states of the Master equation.

preprint2020arXiv

Bifurcation and chaos in nonlinear Lindblad equations

The Lindblad equation describes the dissipative time evolution of a density matrix that characterizes an open quantum system in contact with its environment. The widespread ensemble interpretation of a density matrix requires its time evolution to be linear. However, when the density matrix is obtained from a mean field theory of interacting quantum systems or from a top-down control by a changing classical environment, the ensemble interpretation is inappropriate and nonlinear dynamics arise naturally. We therefore study the dynamical behavior of nonlinear Lindblad equations using the example of a two-level system. By using techniques developed for classical dynamical systems we show that various types of bifurcations and even chaotic dynamics can occur. We also discuss experimental situations for which our results could be relevant.

preprint2020arXiv

Predicting properties of the stationary probability currents for two-species reaction systems without solving the Fokker-Planck equation

We derive methods for estimating the topology of the stationary probability current $\vec{j}_s$ of the two-species Fokker-Planck equation (FPE) without the need to solve the FPE. These methods are chosen such that they become exact in certain limits, such as infinite system size or vanishing coupling between species in the diffusion matrix. The methods make predictions about the fixed points of $\vec{j}_s$ and their relation to extrema of the stationary probability distribution and to fixed points of the convective field, which is related to the deterministic drift of the system. Furthermore, they predict the rotation sense of $\vec{j}_s$ around extrema of the stationary probability distribution. Even though these methods cannot be proven to be valid away from extrema, the boundary lines between regions with different rotation sense are obtained with surprising accuracy. We illustrate and test these method using simple reaction systems with only one coupling term between the two species as well as a few generic reaction networks taken from literature. We use it also to investigate the shape of non-physical probability currents occurring in reaction systems with detailed balance due to the approximations involved in deriving the Fokker-Planck equation.

preprint2020arXiv

Relation between the convective field and the stationary probability distribution of chemical reaction networks

We investigate the relation between the stationary probability distribution of chemical reaction systems and the convective field derived from the chemical Fokker-Planck equation (CFPE) by comparing predictions of the convective field to the results of stochastic simulations based on Gillespie's algorithm. The convective field takes into account the drift term of the CFPE and the reaction bias introduced by the diffusion term. For one-dimensional systems, fixed points and bifurcations of the convective field correspond to extrema and phenomenological bifurcations of the stationary probability distribution whenever the CFPE is a good approximation to the stochastic dynamics. This provides an efficient way to calculate the effect of system size on the number and location of probability maxima and their phenomenological bifurcations in parameter space. For two-dimensional systems, we study models that have saddle-node and Hopf bifurcations in the macroscopic limit. Here, the existence of two stable fixed points of the convective field correlates either with two peaks of the stationary probability distribution, or with a peak and a shoulder. In contrast, a Hopf bifurcation that occurs in the convective field for decreasing system size is not accompanied by the onset of a crater-shaped probability distribution; decreasing system size rather destroys craters and replaces them by local maxima.

preprint2019arXiv

Emergence of time

Microphysical laws are time reversible, but macrophysics, chemistry and biology are not. This chapter explores how this asymmetry (a classic example of a broken symmetry) arises due to the cosmological context, where a non-local Direction of Time is imposed by the expansion of the universe. This situation is best represented by an Evolving Block Universe, where local arrows of time (thermodynamic, electrodynamic, gravitational, wave, quantum, biological) emerge in concordance with the Direction of Time because a global Past Condition results in the Second Law of Thermodynamics pointing to the future. At the quantum level, the indefinite future changes to the definite past due to quantum wave function collapse events.

preprint2019arXiv

What condensed matter physics and statistical physics teach us about the limits of unitary time evolution

The Schrodinger equation for a macroscopic number of particles is linear in the wave function, deterministic, and invariant under time reversal. In contrast, the concepts used and calculations done in statistical physics and condensed matter physics involve stochasticity, nonlinearities, irreversibility, top-down effects, and elements from classical physics. This paper analyzes several methods used in condensed matter physics and statistical physics and explains how they are in fundamental ways incompatible with the above properties of the Schrodinger equation. The problems posed by reconciling these approaches to unitary quantum mechanics are of a similar type as the quantum measurement problem. This paper therefore argues that rather than aiming at reconciling these contrasts one should use them to identify the limits of quantum mechanics. The thermal wave length and thermal time indicate where these limits are for (quasi-)particles that constitute the thermal degrees of freedom.

preprint2016arXiv

Fragile-to-strong transition in liquid silica

We investigate anomalies in liquid silica with molecular dynamics simulations and present evidence for a fragile-to-strong transition at around 3100 K-3300 K. To this purpose, we studied the structure and dynamical properties of silica over a wide temperature range, finding four indicators of a fragile-to-strong transition. First, there is a density minimum at around 3000 K and a density maximum at 4700 K. The turning point is at 3400 K. Second, the local structure characterized by the tetrahedral order parameter changes dramatically around 3000 K from a higher-ordered, lower-density phase to a less ordered, higher-density phase. Third, the correlation time τ changes from an Arrhenius behavior below 3300 K to a Vogel-Fulcher-Tammann behavior at higher temperatures. Fourth, the Stokes-Einstein relation holds for temperatures below 3000 K, but is replaced by a fractional relation above this temperature. Furthermore, our data indicate that dynamics become again simple above 5000 K, with Arrhenius behavior and a classical Stokes-Einstein relation.

preprint2016arXiv

How many bits specify a quantum state?

Quantum mechanics suggests that nature is discrete, with one state per phase space volume $\hbar^{3N}$. This appears to contradict the idea that the state of an N-particle system can have infinite precision and is described by a set of exponentially many complex numbers. Using a finite-temperature gas confined in a box as an example, this short paper argues that there are indeed limits to the precision of wave functions, and that this may help at understanding the quantum-to-classical transition.

preprint2016arXiv

The influence of dispersal on a predator-prey system with two habitats

Dispersal between different habitats influences the dynamics and stability of populations considerably. Furthermore, these effects depend on the local interactions of a population with other species. Here, we perform a general and comprehensive study of the simplest possible system that includes dispersal and local interactions, namely a 2-patch 2-species system. We evaluate the impact of dispersal on stability and on the occurrence of bifurcations, including pattern forming bifurcations that lead to spatial heterogeneity, in 19 different classes of models with the help of the generalized modelling approach. We find that dispersal often destabilizes equilibria, but it can stabilize them if it increases population losses. If dispersal is nonrandom, i.e. if emigration or immigration rates depend on population densities, the correlation of stability with migration rates is positive in part of the models. We also find that many systems show all four types of bifurcations and that antisynchronous oscillations occur mostly with nonrandom dispersal.

preprint2015arXiv

Evolutionary food web model based on body masses gives realistic networks with permanent species turnover

The networks of predator-prey interactions in ecological systems are remarkably complex, but nevertheless surprisingly stable in terms of long term persistence of the system as a whole. In order to understand the mechanism driving the complexity and stability of such food webs, we developed an eco-evolutionary model in which new species emerge as modifications of existing ones and dynamic ecological interactions determine which species are viable. The food-web structure thereby emerges from the dynamical interplay between speciation and trophic interactions. The proposed model is less abstract than earlier evolutionary food web models in the sense that all three evolving traits have a clear biological meaning, namely the average body mass of the individuals, the preferred prey body mass, and the width of their potential prey body mass spectrum. We observed networks with a wide range of sizes and structures and high similarity to natural food webs. The model networks exhibit a continuous species turnover, but massive extinction waves that affect more than $50 \%$ of the network are not observed.

preprint2015arXiv

The influence of the migration network topology on the stability of a small food web

The stability of ecosystems as well as the relation between topology and dynamics on multilayer networks are important questions that are usually discussed in separate communities. Here, we combine these two topics by investigating the influence of the topology of the migration network on the stability of a four-species foodweb module on six patches. The parameters are chosen such that the dynamics on an isolated patch have a periodic attractor with all four species present as well as an attractor where the prey that is preferred by the top predator dies out. The stability measure used here is robustness, which is the average proportion of surviving species in the system, and which shows a complex dependence on the migration rate. We use principal component analysis to quantify the migration network structure in terms of the most relevant network measures, and we evaluate correlations between these measures and characteristics of the robustness curves. Our most important findings are that higher connectivity of the migration network leads to a larger maximum robustness, that a broad distribution of connectivities favors extinction of the preferred prey at intermediate migration rates, and that migration topologies with a larger betweenness centrality are more prone to extinction of the preferred prey at the onset of synchronization. Our study thus demonstrates a strong correlation between dynamical robustness and spatial topology and can serve as an example for similar studies in other types of multilayer networks.

preprint2014arXiv

On the relation between the second law of thermodynamics and classical and quantum mechanics

In textbooks on statistical mechanics, one finds often arguments based on classical mechanics, phase space and ergodicity in order to justify the second law of thermodynamics. However, the basic equations of motion of classical mechanics are deterministic and reversible, while the second law of thermodynamics is irreversible and not deterministic, because it states that a system forgets its past when approaching equilibrium. I argue that all "derivations" of the second law of thermodynamics from classical mechanics include additional assumptions that are not part of classical mechanics. The same holds for Boltzmann's H-theorem. Furthermore, I argue that the coarse-graining of phase-space that is used when deriving the second law cannot be viewed as an expression of our ignorance of the details of the microscopic state of the system, but reflects the fact that the state of a system is fully specified by using only a finite number of bits, as implied by the concept of entropy, which is related to the number of different microstates that a closed system can have. While quantum mechanics, as described by the Schroedinger equation, puts this latter statement on a firm ground, it cannot explain the irreversibility and stochasticity inherent in the second law.

preprint2013arXiv

Formation of the frozen core in critical Boolean Networks

We investigate numerically and analytically the formation of the frozen core in critical random Boolean networks with biased functions. We demonstrate that a previously used efficient algorithm for obtaining the frozen core, which starts from the nodes with constant functions, fails when the number of inputs per node exceeds 4. We present computer simulation data for the process of formation of the frozen core and its robustness, and we show that several important features of the data can be derived by using a mean-field calculation.

preprint2013arXiv

Scaling laws in critical random Boolean networks with general in- and out-degree distributions

We evaluate analytically and numerically the size of the frozen core and various scaling laws for critical Boolean networks that have a power-law in- and/or out-degree distribution. To this purpose, we generalize an efficient method that has previously been used for conventional random Boolean networks and for networks with power-law in-degree distributions. With this generalization, we can also deal with power-law out-degree distributions. When the power-law exponent is between 2 and 3, the second moment of the distribution diverges with network size, and the scaling exponent of the nonfrozen nodes depends on the degree distribution exponent. Furthermore, the exponent depends also on the dependence of the cutoff of the degree distribution on the system size. Altogether, we obtain an impressive number of different scaling laws depending on the type of cutoff as well as on the exponents of the in- and out-degree distributions. We confirm our scaling arguments and analytical considerations by numerical investigations.

preprint2012arXiv

A one-dimensional model with water-like anomalies and two phase transitions

We investigate a one-dimensional model that shows several properties of water. The model combines the long-range attraction of the van der Waals model with the nearest-neighbor interaction potential by Ben-Naim, which is a step potential that includes a hard core and a potential well. Starting from the analytical expression for the partition function, we determine numerically the Gibbs energy and other thermodynamic quantities. The model shows two phase transitions, which can be interpreted as the liquid-gas transition and a transition between a high-density and a low-density liquid. At zero temperature, the low-density liquid goes into the crystalline phase. Furthermore, we find several anomalies that are considered characteristic for water. We explore a wide range of pressure and temperature values and the dependence of the results on the depth and width of the potential well.

preprint2010arXiv

A three-species model explaining cyclic dominance of pacific salmon

The four-year oscillations of the number of spawning sockeye salmon (Oncorhynchus nerka) that return to their native stream within the Fraser River basin in Canada are a striking example of population oscillations. The period of the oscillation corresponds to the dominant generation time of these fish. Various - not fully convincing - explanations for these oscillations have been proposed, including stochastic influences, depensatory fishing, or genetic effects. Here, we show that the oscillations can be explained as a stable dynamical attractor of the population dynamics, resulting from a strong resonance near a Neimark Sacker bifurcation. This explains not only the long-term persistence of these oscillations, but also reproduces correctly the empirical sequence of salmon abundance within one period of the oscillations. Furthermore, it explains the observation that these oscillations occur only in sockeye stocks originating from large oligotrophic lakes, and that they are usually not observed in salmon species that have a longer generation time.

preprint2010arXiv

Boolean networks with robust and reliable trajectories

We construct and investigate Boolean networks that follow a given reliable trajectory in state space, which is insensitive to fluctuations in the updating schedule, and which is also robust against noise. Robustness is quantified as the probability that the dynamics return to the reliable trajectory after a perturbation of the state of a single node. In order to achieve high robustness, we navigate through the space of possible update functions by using an evolutionary algorithm. We constrain the networks to having the minimum number of connections required to obtain the reliable trajectory. Surprisingly, we find that robustness always reaches values close to 100 percent during the evolutionary optimization process. The set of update functions can be evolved such that it differs only slightly from that of networks that were not optimized with respect to robustness. The state space of the optimized networks is dominated by the basin of attraction of the reliable trajectory.

preprint2010arXiv

Boolean versus continuous dynamics on simple two-gene modules

We investigate the dynamical behavior of simple modules composed of two genes with two or three regulating connections. Continuous dynamics for mRNA and protein concentrations is compared to a Boolean model for gene activity. Using a generalized method, we study within a single framework different continuous models and different types of regulatory functions, and establish conditions under which the system can display stable oscillations. These conditions concern the time scales, the degree of cooperativity of the regulating interactions, and the signs of the interactions. Not all models that show oscillations under Boolean dynamics can have oscillations under continuous dynamics, and vice versa.

preprint2010arXiv

Coexisting patterns of population oscillations: the degenerate Neimark Sacker bifurcation as a generic mechanism

We investigate a population dynamics model that exhibits a Neimark Sacker bifurcation with a period that is naturally close to 4. Beyond the bifurcation, the period becomes soon locked at 4 due to a strong resonance, and a second attractor of period 2 emerges, which coexists with the first attractor over a considerable parameter range. A linear stability analysis and a numerical investigation of the second attractor reveal that the bifurcations producing the second attractor occur naturally in this type of system.

preprint2009arXiv

Boolean networks with reliable dynamics

We investigated the properties of Boolean networks that follow a given reliable trajectory in state space. A reliable trajectory is defined as a sequence of states which is independent of the order in which the nodes are updated. We explored numerically the topology, the update functions, and the state space structure of these networks, which we constructed using a minimum number of links and the simplest update functions. We found that the clustering coefficient is larger than in random networks, and that the probability distribution of three-node motifs is similar to that found in gene regulation networks. Among the update functions, only a subset of all possible functions occur, and they can be classified according to their probability. More homogeneous functions occur more often, leading to a dominance of canalyzing functions. Finally, we studied the entire state space of the networks. We observed that with increasing systems size, fixed points become more dominant, moving the networks close to the frozen phase.

preprint2007arXiv

Evolution of Canalizing Boolean Networks

Boolean networks with canalizing functions are used to model gene regulatory networks. In order to learn how such networks may behave under evolutionary forces, we simulate the evolution of a single Boolean network by means of an adaptive walk, which allows us to explore the fitness landscape. Mutations change the connections and the functions of the nodes. Our fitness criterion is the robustness of the dynamical attractors against small perturbations. We find that with this fitness criterion the global maximum is always reached and that there is a huge neutral space of 100% fitness. Furthermore, in spite of having such a high degree of robustness, the evolved networks still share many features with "chaotic" networks.

preprint2001arXiv

Biological Evolution and Statistical Physics

This review is an introduction to theoretical models and mathematical calculations for biological evolution, aimed at physicists. The methods in the field are naturally very similar to those used in statistical physics, although the majority of publications appeared in biology journals. The review has three parts, which can be read independently. The first part deals with evolution in fitness landscapes and includes Fisher's theorem, adaptive walks, quasispecies models, effects of finite population sizes, and neutral evolution. The second part studies models of coevolution, including evolutionary game theory, kin selection, group selection, sexual selection, speciation, and coevolution of hosts and parasites. The third part discusses models for networks of interacting species and their extinction avalanches. Throughout the review, attention is paid to giving the necessary biological information, and to pointing out the assumptions underlying the models, and their limits of validity.

preprint1996arXiv

Forest fires and other examples of self-organized criticality

We review the properties of the self-organized critical (SOC) forest-fire model. The paradigm of self-organized criticality refers to the tendency of certain large dissipative systems to drive themselves into a critical state independent of the initial conditions and without fine-tuning of the parameters. After an introduction, we define the rules of the model and discuss various large-scale structures which may appear in this system. The origin of the critical behavior is explained, critical exponents are introduced, and scaling relations between the exponents are derived. Results of computer simulations and analytical calculations are summarized. The existence of an upper critical dimension and the universality of the critical behavior under changes of lattice symmetry or the introduction of immunity are discussed. A survey of interesting modifications of the forest-fire model is given. Finally, several other important SOC models are briefly described.

preprint1996arXiv

Self--organized criticality due to a separation of energy scales

Certain systems with slow driving and avalanche-like dissipation events are naturally close to a critical point when the ratio of two energy scales is large. The first energy scale is the threshold above which an avalanche is triggered, the second scale is the threshold above which a site is affected by an avalanche. I present results of computer simulations, and a mean-field theory.

preprint1995arXiv

Energy Barriers to Motion of Flux Lines in Random Media

We propose algorithms for determining both lower and upper bounds for the energy barriers encountered by a flux line in moving through a two-dimensional random potential. Analytical arguments, supported by numerical simulations, suggest that these bounds scale with the length $t$ of the line as $t^{1/3}$ and $t^{1/3}\sqrt{\ln t}$, respectively. This provides the first confirmation of the hypothesis that barriers have the same scaling as the fluctuation in the free energy. \pacs{PACS numbers: 74.60.Ge, 05.70.Ln, 05.40.+j}

preprint1994arXiv

Universality in the One-Dimensional Self-Organized Critical Forest-Fire Model

We modify the rules of the self-organized critical forest-fire model in one dimension by allowing the fire to jump over holes of $\le k$ sites. An analytic calculation shows that not only the size distribution of forest clusters but also the size distribution of fires is characterized by the same critical exponent as in the nearest-neighbor model, i.e. the critical behavior of the model is universal. Computer simulations confirm the analytic results.