Showing 1–11 of 11 results for author: Hartmann, A

Coexistence of asynchronous and clustered dynamics in noisy inhibitory neural networks

Authors: Yannick Feld, Alexander K. Hartmann, Alessandro Torcini

Abstract: A regime of coexistence of asynchronous and clustered dynamics is analyzed for globally coupled homogeneous and heterogeneous inhibitory networks of quadratic integrate-and-fire (QIF) neurons subject to Gaussian noise. The analysis is based on accurate extensive simulations and complemented by a mean-field description in terms of low-dimensional next generation neural mass models for heterogeneous… ▽ More A regime of coexistence of asynchronous and clustered dynamics is analyzed for globally coupled homogeneous and heterogeneous inhibitory networks of quadratic integrate-and-fire (QIF) neurons subject to Gaussian noise. The analysis is based on accurate extensive simulations and complemented by a mean-field description in terms of low-dimensional next generation neural mass models for heterogeneously distributed synaptic couplings. The asynchronous regime is observable at low noise and becomes unstable via a sub-critical Hopf bifurcation at sufficiently large noise. This gives rise to a coexistence region between the asynchronous and the clustered regime. The clustered phase is characterized by population bursts in the γ-range (30-120 Hz), where neurons are split in two equally populated clusters firing in alternation. This clustering behaviour is quite peculiar: despite the global activity being essentially periodic, single neurons display switching between the two clusters due to heterogeneity and/or noise. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Comments: 36 pages - 22 figures

Vertebrate pollinators: phase transition in a time-dependent generalized traveling-salesperson problem

Authors: M. Jungsbluth, J. Thiele, Y. Winter, H. Schawe, A. K. Hartmann

Abstract: We introduce a model for the global optimization problem of nectar harvesting by flower visitors, e.g., nectar-feeding bats, as a generalization of the (multiple) traveling-salesperson problem (TSP). The model includes multiple independent animals and many flowers with time-dependent content. This provides an ensemble of realistic combinatorial optimization problems, in contrast to previously stud… ▽ More We introduce a model for the global optimization problem of nectar harvesting by flower visitors, e.g., nectar-feeding bats, as a generalization of the (multiple) traveling-salesperson problem (TSP). The model includes multiple independent animals and many flowers with time-dependent content. This provides an ensemble of realistic combinatorial optimization problems, in contrast to previously studied models like random Satisfiability or standard TSP. We numerically studied the optimum harvesting of these foragers, with parameters obtained from experiments, by using genetic algorithms. For the distribution of travel distances, we find a power-law (or Lévy) distribution, as often found for natural foragers. Note, in contrast to many models, we make no assumption about the nature of the flight-distance distribution, the power law just emerges. This is in contrast to the TSP, where we find in the present study an exponential tail. Furthermore, the optimization problem exhibits a {phase transition}, similar to the TSP, at a critical value for the amount of nectar which can be harvested. This phase transition coincides with a dramatic increase in the typical running time of the optimization algorithm. For the value of the critical exponent nu, describing the divergence of the correlation length, we find nu=1.7(4), which is on the other hand compatible with the value found for the TSP. Finally, we also present data from field experiments in Costa Rica for the resource use for freely visiting flower bats. We found that the temporal patterns in experiments and model agree remarkably, confirming our model. Also the data show that the bats are able to memorize the positions of food sources and optimize, at least partially, their routes. △ Less

Submitted 21 March, 2018; originally announced March 2018.

Comments: 12 pages, 9 figures

Large-deviation properties of the extended Moran model

Authors: Alexander K. Hartmann, Thierry E. Huillet

Abstract: The distributions of the times to the first common ancestor t_mrca is numerically studied for an ecological population model, the extended Moran model. This model has a fixed population size N. The number of descendants is drawn from a beta distribution Beta(alpha, 2-alpha) for various choices of alpha. This includes also the classical Moran model (alpha->0) as well as the uniform distribution (al… ▽ More The distributions of the times to the first common ancestor t_mrca is numerically studied for an ecological population model, the extended Moran model. This model has a fixed population size N. The number of descendants is drawn from a beta distribution Beta(alpha, 2-alpha) for various choices of alpha. This includes also the classical Moran model (alpha->0) as well as the uniform distribution (alpha=1). Using a statistical mechanics-based large-deviation approach, the distributions can be studied over an extended range of the support, down to probabilities like 10^{-70}, which allowed us to study the change of the tails of the distribution when varying the value of alpha in [0,2]. We find exponential distributions p(t_mrca)~ delta^{t_mrca} in all cases, with systematically varying values for the base delta. Only for the cases alpha=0 and alpha=1, analytical results are known, i.e., delta=\exp(-2/N^2) and delta=2/3, respectively. We recover these values, confirming the validity of our approach. Finally, we also study the correlations between t_mrca and the number of descendants. △ Less

Submitted 27 October, 2018; v1 submitted 20 October, 2017; originally announced October 2017.

Comments: 8 pages, 8 figures

Journal ref: Phys. Rev. E 98, 042416 (2018)

Precision and accuracy of single-molecule FRET measurements - a worldwide benchmark study

Authors: Björn Hellenkamp, Sonja Schmid, Olga Doroshenko, Oleg Opanasyuk, Ralf Kühnemuth, Soheila Rezaei Adariani, Anders Barth, Victoria Birkedal, Mark E. Bowen, Hongtao Chen, Thorben Cordes, Tobias Eilert, Carel Fijen, Markus Götz, Giorgos Gouridis, Enrico Gratton, Taekjip Ha, Christian A. Hanke, Andreas Hartmann, Jelle Hendrix, Lasse L. Hildebrandt, Johannes Hohlbein, Christian G. Hübner, Eleni Kallis, Achillefs N. Kapanidis , et al. (28 additional authors not shown)

Abstract: Single-molecule Förster resonance energy transfer (smFRET) is increasingly being used to determine distances, structures, and dynamics of biomolecules in vitro and in vivo. However, generalized protocols and FRET standards ensuring both the reproducibility and accuracy of measuring FRET efficiencies are currently lacking. Here we report the results of a worldwide, comparative, blind study, in whic… ▽ More Single-molecule Förster resonance energy transfer (smFRET) is increasingly being used to determine distances, structures, and dynamics of biomolecules in vitro and in vivo. However, generalized protocols and FRET standards ensuring both the reproducibility and accuracy of measuring FRET efficiencies are currently lacking. Here we report the results of a worldwide, comparative, blind study, in which 20 labs determined the FRET efficiencies of several dye-labeled DNA duplexes. Using a unified and straightforward method, we show that FRET efficiencies can be obtained with a standard deviation between $Δ$E = +-0.02 and +-0.05. We further suggest an experimental and computational procedure for converting FRET efficiencies into accurate distances. We discuss potential uncertainties in the experiment and the modelling. Our extensive quantitative assessment of intensity-based smFRET measurements and correction procedures serve as an essential step towards validation of distance networks with the ultimate aim to archive reliable structural models of biomolecular systems obtained by smFRET-based hybrid methods. △ Less

Submitted 29 December, 2017; v1 submitted 10 October, 2017; originally announced October 2017.

Comments: 27 pages, 5 figures, supplementary material, corrections: p20 last sentence; eqn. 13 & 14

Journal ref: Natue Methods, vol. 15, pages 669-676, 2018

Score distributions of gapped multiple sequence alignments down to the low-probability tail

Authors: Pascal Fieth, Alexander K. Hartmann

Abstract: Assessing the significance of alignment scores of optimally aligned DNA or amino acid sequences can be achieved via the knowledge of the score distribution of random sequences. But this requires obtaining the distribution in the biologically relevant high-scoring region, where the probabilities are exponentially small. For gapless local alignments of infinitely long sequences this distribution is… ▽ More Assessing the significance of alignment scores of optimally aligned DNA or amino acid sequences can be achieved via the knowledge of the score distribution of random sequences. But this requires obtaining the distribution in the biologically relevant high-scoring region, where the probabilities are exponentially small. For gapless local alignments of infinitely long sequences this distribution is known analytically to follow a Gumbel distribution. Distributions for gapped local alignments and global alignments of finite lengths can only be obtained numerically. To obtain result for the small-probability region, specific statistical mechanics-based rare-event algorithms can be applied. In previous studies, this was achieved for pairwise alignments. They showed that, contrary to results from previous simple sampling studies, strong deviations from the Gumbel distribution occur in case of finite sequence lengths. Here we extend the studies to the for practical applications in Molecular Biology much more relevant case of multiple sequence alignments with gaps. We study the distributions of scores over a large range of the support, reaching probabilities as small as 10^-160, for global and local (sum-of-pair scores) multiple alignments. We find that even after suitable rescaling, eliminating the sequence-length dependence, the distributions for multiple alignment differ from the pairwise alignment case. Furthermore, we also show that the previously discussed Gaussian correction to the Gumbel distribution needs to be refined, also for the case of pairwise alignments. △ Less

Submitted 7 December, 2015; originally announced December 2015.

Journal ref: Phys. Rev. E 94, 022127 (2016)

Finite-temperature local protein sequence alignment: percolation and free-energy distribution

Authors: S Wolfsheimer, O Melchert, AK Hartmann

Abstract: Sequence alignment is a tool in bioinformatics that is used to find homological relationships in large molecular databases. It can be mapped on the physical model of directed polymers in random media. We consider the finite-temperature version of local sequence alignment for proteins and study the transition between the linear phase and the biologically relevant logarithmic phase, where the free… ▽ More Sequence alignment is a tool in bioinformatics that is used to find homological relationships in large molecular databases. It can be mapped on the physical model of directed polymers in random media. We consider the finite-temperature version of local sequence alignment for proteins and study the transition between the linear phase and the biologically relevant logarithmic phase, where the free-energy grows linearly or logarithmically with the sequence length. By means of numerical simulations and finite-size scaling analysis we determine the phase diagram in the plane that is spanned by the gap costs and the temperature. We use the most frequently used parameter set for protein alignment. The critical exponents that describe the parameter driven transition are found to be explicitly temperature dependent. Furthermore, we study the shape of the (free-) energy distribution close to the transition by rare-event simulations down to probabilities of the order $10^{-64}$. It is well known that, in the logarithmic region, the optimal score distribution (T=0) is described by a modified Gumbel distribution. We confirm that this also applies for the free-energy distribution ($T>0$). However, in the linear phase, the distribution crosses over to a modified Gaussian distribution. △ Less

Submitted 13 July, 2009; originally announced July 2009.

RNA Secondary Structures: Complex Statics and Glassy Dynamics

Authors: S. Wolfsheimer, B. Burghardt, A. Mann, A. K. Hartmann

Abstract: Models for RNA secondary structures (the topology of folded RNA) without pseudo knots are disordered systems with a complex state-space below a critical temperature. Hence, a complex dynamical (glassy) behavior can be expected, when performing Monte Carlo simulation. Interestingly, in contrast to most other complex systems, the ground states and the density of states can be computed in polynomia… ▽ More Models for RNA secondary structures (the topology of folded RNA) without pseudo knots are disordered systems with a complex state-space below a critical temperature. Hence, a complex dynamical (glassy) behavior can be expected, when performing Monte Carlo simulation. Interestingly, in contrast to most other complex systems, the ground states and the density of states can be computed in polynomial time exactly using transfer matrix methods. Hence, RNA secondary structure is an ideal model to study the relation between static/thermodynamic properties and dynamics of algorithms. Also they constitute an ideal benchmark system for new Monte Carlo methods. Here we considered three different recent Monte Carlo approaches: entropic sampling using flat histograms, optimized-weights ensembles, and ParQ, which estimates the density of states from transition matrices. These methods were examined by comparing the obtained density of states with the exact results. We relate the complexity seen in the dynamics of the Monte Carlo algorithms to static properties of the phase space by studying the correlations between tunneling times, sampling errors, amount of meta-stable states and degree of ultrametricity at finite temperature. △ Less

Submitted 2 February, 2008; v1 submitted 30 October, 2007; originally announced October 2007.

RNA secondary structure design

Authors: Bernd Burghardt, Alexander K. Hartmann

Abstract: We consider the inverse-folding problem for RNA secondary structures: for a given (pseudo-knot-free) secondary structure find a sequence that has that structure as its ground state. If such a sequence exists, the structure is called designable. We implemented a branch-and-bound algorithm that is able to do an exhaustive search within the sequence space, i.e., gives an exact answer whether such a… ▽ More We consider the inverse-folding problem for RNA secondary structures: for a given (pseudo-knot-free) secondary structure find a sequence that has that structure as its ground state. If such a sequence exists, the structure is called designable. We implemented a branch-and-bound algorithm that is able to do an exhaustive search within the sequence space, i.e., gives an exact answer whether such a sequence exists. The bound required by the branch-and-bound algorithm are calculated by a dynamic programming algorithm. We consider different alphabet sizes and an ensemble of random structures, which we want to design. We find that for two letters almost none of these structures are designable. The designability improves for the three-letter case, but still a significant fraction of structures is undesignable. This changes when we look at the natural four-letter case with two pairs of complementary bases: undesignable structures are the exception, although they still exist. Finally, we also study the relation between designability and the algorithmic complexity of the branch-and-bound algorithm. Within the ensemble of structures, a high average degree of undesignability is correlated to a long time to prove that a given structure is (un-)designable. In the four-letter case, where the designability is high everywhere, the algorithmic complexity is highest in the region of naturally occurring RNA. △ Less

Submitted 15 September, 2006; originally announced September 2006.

Comments: 11 pages, 10 figures

Dependence of RNA secondary structure on the energy model

Authors: Bernd Burghardt, Alexander K. Hartmann

Abstract: We analyze a microscopic RNA model, which includes two widely used models as limiting cases, namely it contains terms for bond as well as for stacking energies. We numerically investigate possible changes in the qualitative and quantitative behaviour while going from one model to the other; in particular we test, whether a transition occurs, when continuously moving from one model to the other.… ▽ More We analyze a microscopic RNA model, which includes two widely used models as limiting cases, namely it contains terms for bond as well as for stacking energies. We numerically investigate possible changes in the qualitative and quantitative behaviour while going from one model to the other; in particular we test, whether a transition occurs, when continuously moving from one model to the other. For this we calculate various thermodynamic quantities, both at zero temperature as well as at finite temperatures. All calculations can be done efficiently in polynomial time by a dynamic programming algorithm. We do not find a sign for transition between the models, but the critical exponent $ν$ of the correlation length, describing the phase transition in all models to an ordered low-temperature phase, seems to depend continuously on the model. Finally, we apply the epsilon-Coupling method, to study low excitations. The exponent $θ$ describing the energy-scaling of the excitations seems to depend not much on the energy model. △ Less

Submitted 7 September, 2004; originally announced September 2004.

Comments: 8 pages, 9 figures

Sampling rare events: statistics of local sequence alignments

Authors: Alexander K. Hartmann

Abstract: A new method to simulate probability distributions in regions where the events are VERY unlikely (e.g. p ~ 10^{-40}) is presented. The basic idea is to represent the underlying probability space by the phase space of a physical system. The system is held at a temperature T, which is chosen such that the system preferably generates configurations which originally have low probabilities. Since the… ▽ More A new method to simulate probability distributions in regions where the events are VERY unlikely (e.g. p ~ 10^{-40}) is presented. The basic idea is to represent the underlying probability space by the phase space of a physical system. The system is held at a temperature T, which is chosen such that the system preferably generates configurations which originally have low probabilities. Since the distribution of such a physical system is know from statistical physics, the original unbiased distribution can be obtained. As an application, local alignment of protein sequences based on BLOSUM62 substitution scores with (12,1) affine gap costs are considered The distribution of optimum sequence-alignment scores S is studied numerically over a large range of scores. The deviation of p(S) from the extreme-value (or Gumbel) distribution is quantified. This deviation decreases with growing sequence length. △ Less

Submitted 13 August, 2001; originally announced August 2001.

Comments: 5 pages, 4 figures, revtex

Journal ref: Phys. Rev. E 65, 056102 (2002)

Comment on `Glassy Transition in a Disordered Model for the RNA Secondary Structure'

Authors: Alexander K. Hartmann

Abstract: In cond-mat/9907125 the low-temperature behavior of a model for RNA secondary structure was studied. It is claimed that the model exhibits a breaking of the replica symmetry, since the width of the distribution P(q) of overlaps may converge to a finite value at T=0. The authors used an exact enumeration method to obtain all ground states for a given RNA sequence. Because of the exponential growi… ▽ More In cond-mat/9907125 the low-temperature behavior of a model for RNA secondary structure was studied. It is claimed that the model exhibits a breaking of the replica symmetry, since the width of the distribution P(q) of overlaps may converge to a finite value at T=0. The authors used an exact enumeration method to obtain all ground states for a given RNA sequence. Because of the exponential growing degeneracy, only sequences up to length L=256 could be studied. Here it is shown that, in contrast to the previous results, by going to much larger sizes as L=2000 the variance coverges towards zero, i.e. P(q) is a delta-function in the thermodynamic limit. △ Less

Submitted 7 April, 2000; v1 submitted 9 August, 1999; originally announced August 1999.

Comments: completely rewritten, comment to cond-mat/9907125 (PRL 84, 2026)

Journal ref: Phys. Rev. Lett. 86, 1382 (2001)