  One of the best examples of traditional analysis of evolutionary stable
strategies (ESS) is provided by the so called Dove-Hawk model. In this paper we
present several enhancements to the model aimed at describing the evolution of
cooperative behavior. In addition to Doves and Hawks we introduce several
groups of Cooperators, who act as Doves within their own group, but as Hawks
outside it. This allows to study how cooperating groups may grow and achieve
stability within a nonrepeating evolutionary games framework. Depending on
initial conditions, the final stable population may have one, all-encompassing
Cooperator population or several competing cliques. After taking into account
that Cooperators bear costs necessary to recognize members of one's own group
it is possible to see populations where Cooperators eventually lose against
Hawks or populations where several cliques of Cooperators coexist with a Hawk
population.


  In a recent experiment, a single DNA double helix is stretched and relaxed in
the presence of spermidine, a short positive polyelectrolyte, and the pulling
force is measured as a function of DNA extension. In a certain range of
spermidine concentration, a force plateau appears whose value shows maximum as
a function of spermidine concentration. We present a quantitative theory of
this plateau force based on the theory of reentrant condensation and derive
almost parabolic behavior of the plateau force as a function of the logarithm
of the spermidine concentration in the range of condensation. Our result is in
good agreement with experimental data.


  We introduce a 2-layer network model for the study of the immunization
dynamics in epidemics. Spreading of an epidemic is modeled as an excitatory
process in a small-world network (body layer) while immunization by prevention
for the disease as a dynamic process in a scale-free network (head layer). It
is shown that prevention indeed turns periodic rages of an epidemic into small
fluctuation. The study also reveals that, in a certain situation, prevention
actually plays an adverse role and helps the disease survive. We argue that the
presence of two different characteristic time scales contributes to the
immunization dynamics observed.


  We introduce a simple model of population dynamics which considers birth and
death rates for every individual that depend on the number of particles in its
neighborhood. The model shows an inhomogeneous quasistationary pattern with
many different clusters of particles.
  We derive the equation for the macroscopic density of particles, perform a
linear stability analysis on it, and show that there is a finite-wavelength
instability leading to pattern formation. This is the responsible for the
approximate periodicity with which the clusters of particles arrange in the
microscopic model.
  In addition, we consider the population when immersed in a fluid medium and
analyze the influence of advection on global properties of the model.


  Competitve hybridization, at the surface and in the bulk, lowers the
sensitivity of DNA chips. Competitive surface hybridization occurs when
different targets can hybridize with the same probe. Competitive bulk
hybridization takes place when the targets can hybridize with free
complementary chains in the solution. The effects of competitive hybridization
on the thermodynamically attainable performance of DNA chips are quantified in
terms of the hybridization isotherms of the spots. These relate the equilibrium
degree of the hybridization to the bulk composition. The hybridization isotherm
emerges as a Langmuir isotherm modified for electrostatic interactions within
the probe layer. The sensitivity of the assay in equilibrium is directly
related to the slope of the isotherm. A simpler description is possible in
terms of $c_{50}$s specifying the bulk composition corresponding to 50%
hybridization at the surface. The effects of competitive hybridization are
important for the quantitative analysis of DNA chip results especially when
used to study point mutations.


  Norms, defined as generally accepted behaviour in societies without central
authority (and thus distinguished from laws), are very powerful mechanism
leading to coherent behaviour of the society members. This paper examines,
within a simple numerical simulation, the various effects that may lead to norm
formation and stability. The approach has been first used by Axelrod, who
proposed two step model of norm and meta-norm enforcement. We present here an
extension and detailed analysis of the original work, as well as several new
ideas that may bear on the norm establishment mechanisms in societies. It turns
out that a relatively simple model for simulated norm enforcement predicts
persistent norm breaking even when it is associated with high punishment
levels. The key factors appear to be the combination of the level of penalty
for breaking the norm and proximity of norm enforcers. We also study a totally
different mechanism of norm establishment, without meta-norms but using instead
the direct bonus mechanism to norm-enforcers.


  The effect of sequence heterogeneity on polynucleotide translocation across a
pore and on simple models of molecular motors such as helicases, DNA
polymerase/exonuclease and RNA polymerase is studied in detail. Pore
translocation of RNA or DNA is biased due to the different chemical
environments on the two sides of the membrane, while the molecular motor motion
is biased through a coupling to chemical energy. An externally applied force
can oppose these biases. For both systems we solve lattice models exactly both
with and without disorder. The models incorporate explicitly the coupling to
the different chemical environments for polymer translocation and the coupling
to the chemical energy (as well as nucleotide pairing energies) for molecular
motors. Using the exact solutions and general arguments we show that the
heterogeneity leads to anomalous dynamics. Most notably, over a range of forces
around the stall force (or stall tension for DNA polymerase/exonuclease
systems) the displacement grows sublinearly as t^\mu with \mu<1. The range over
which this behavior can be observed experimentally is estimated for several
systems and argued to be detectable for appropriate forces and buffers. Similar
sequence heterogeneity effects may arise in the packing of viral DNA.


  We further develop the large $ N $ formalism presented by some of us in
earlier works in order to recursively calculate the partition function of a
singly pseudoknotted RNA. We demonstrate that this calculation takes time
proportional to the sixth power of the length of the RNA. The algorithm itself
is presented in a self-contained form for the convenience of readers interested
in implementing it.


  We consider the passage of long polymers of length N through a hole in a
membrane. If the process is slow, it is in principle possible to focus on the
dynamics of the number of monomers s on one side of the membrane, assuming that
the two segments are in equilibrium. The dynamics of s(t) in such a limit would
be diffusive, with a mean translocation time scaling as N^2 in the absence of a
force, and proportional to N when a force is applied. We demonstrate that the
assumption of equilibrium must break down for sufficiently long polymers (more
easily when forced), and provide lower bounds for the translocation time by
comparison to unimpeded motion of the polymer. These lower bounds exceed the
time scales calculated on the basis of equilibrium, and point to anomalous
(sub-diffusive) character of translocation dynamics. This is explicitly
verified by numerical simulations of the unforced translocation of a
self-avoiding polymer. Forced translocation times are shown to strongly depend
on the method by which the force is applied. In particular, pulling the polymer
by the end leads to much longer times than when a chemical potential difference
is applied across the membrane. The bounds in these cases grow as N^2 and
N^{1+\nu}, respectively, where \nu is the exponent that relates the scaling of
the radius of gyration to N. Our simulations demonstrate that the actual
translocation times scale in the same manner as the bounds, although influenced
by strong finite size effects which persist even for the longest polymers that
we considered (N=512).


  Microtubules are stiff filamentary proteins that constitute an important
component of the cytoskeleton of cells. These are known to exhibit a dynamic
instability. A steadily growing microtubule can suddenly start depolymerizing
very rapidly; this phenomenon is known as ``catastrophe''. However, often a
shrinking microtubule is ``rescued'' and starts polymerizing again. Here we
develope a model for the polymerization-depolymerization dynamics of
microtubules in the presence of {\it catastrophe-suppressing drugs}. Solving
the dynamical equations in the steady-state, we derive exact analytical
expressions for the length distributions of the microtubules tipped with
drug-bound tubulin subunits as well as those of the microtubules, in the
growing and shrinking phases, tipped with drug-free pure tubulin subunits. We
also examine the stability of the steady-state solutions.


  A generalized-ensemble technique, multicanonical sampling, is used to study
the folding of a 34-residue human parathyroid hormone fragment. An all-atom
model of the peptide is employed and the protein-solvent interactions are
approximated by an implicit solvent. Our results demonstrate that
generalized-ensemble simulations are well suited to sample low-energy
structures of such large polypeptides. Configurations with a root-mean-square
deviation (rmsd) to the crystal structure of less than one \AA are found.
Finally, we discuss limitations of our implicit solvent model.


  We study the role of the noise in the dynamics of two competing species. We
consider generalized Lotka-Volterra equations in the presence of a
multiplicative noise, which models the interaction between the species and the
environment. The interaction parameter between the species is a random process
which obeys a stochastic differential equation with a generalized bistable
potential in the presence of a periodic driving term, which accounts for the
environment temperature variation. We find noise-induced periodic oscillations
of the species concentrations and stochastic resonance phenomenon. We find also
a nonmonotonic behavior of the mean extinction time of one of the two competing
species as a function of the additive noise intensity.


  We study the time evolution of two ecosystems in the presence of external
noise and climatic periodical forcing by a generalized Lotka-Volterra (LV)
model. In the first ecosystem, composed by two competing species, we find noise
induced phenomena such as: (i) quasi deterministic oscillations, (ii)
stochastic resonance, (iii) noise delayed extinction and (iv) spatial patterns.
In the second ecosystem, composed by three interacting species (one predator
and two preys), using a discrete model of the LV equations we find that the
time evolution of the spatial patterns is strongly dependent on the initial
conditions of the three species.


  We separate double stranded lambda phage DNA by applying a fixed force at a
constant temperature ranging from 15C to 50C, and measure the minimum force
required to separate the two strands, providing the first experimental
determination of the phase boundary between single stranded DNA and double
stranded DNA in the temperature- force plane. The measurements also offer
information on the free energy of dsDNA at temperatures where dsDNA does not
thermally denature in the absence of force. While parts of the phase diagram
can be explained using existing models and free energy parameters, others
deviate significantly. Possible reasons for the deviations between theory and
experiment are considered.


  We investigate on a possible way to connect the presence of Low-Complexity
Sequences (LCS) in DNA genomes and the nonstationary properties of base
correlations. Under the hypothesis that these variations signal a change in the
DNA function, we use a new technique, called Non-Stationarity Entropic Index
(NSEI) method, and we prove that this technique is an efficient way to detect
functional changes with respect to a random baseline. The remarkable aspect is
that NSEI does not imply any training data or fitting parameter, the only
arbitrarity being the choice of a marker in the sequence. We make this choice
on the basis of biological information about LCS distributions in genomes. We
show that there exists a correlation between changing the amount in LCS and the
ratio of long- to short-range correlation.


  We investigate the sensitivity of production rates (activities) of the
regulatory proteins CI (repressor) and Cro at the right operator (O_R}) of
bacteriophage lambda. The DNA binding energies of CI, Cro, and RNA polymerase
are perturbed to check the uncertainty of the activity, due to the experimental
error, by means of a computational scattering method according to which the
binding energies are simultaneously chosen at random around the literature
values, with a width corresponding to the experimental error. In a grand
canonical ensemble, with the randomly drawn protein-DNA binding energies as
input, we calculate the corresponding activities of the promoters P_RM and P_R.
By repeating this procedure we obtain a mean value of the activity that roughly
corresponds to wild-type (unperturbed) activity. The standard deviation
emerging from this scheme, a measure of the sensitivity due to experimental
error, is significant (typically > 20% relative to wild-type activity), but
still the promoter activities are sufficiently separated to make the switch
feasible. We also suggest a new compact way of presenting repressor and Cro
data.


  We construct an algorithm that generates large, band-diagonal transition
matrices for a totally asymmetric exclusion process (TASEP) with local hopping
rate inhomogeneities. The matrices are diagonalized numerically to find
steady-state currents of TASEPs with local variations in hopping rate. The
results are then used to investigate clustering of slow codons along mRNA.
Ribosome density profiles near neighboring clusters of slow codons interact,
enhancing suppression of ribosome throughput when such bottlenecks are closely
spaced. Increasing the slow codon cluster size, beyond $\approx 3-4$, does not
significantly reduce ribosome current. Our results are verified by extensive
Monte-Carlo simulations and provide a biologically-motivated explanation for
the experimentally-observed clustering of low-usage codons.


  We calculate thermodynamic quantities of HP lattice proteins by means of a
multicanonical chain growth algorithm that connects the new variants of the
Pruned-Enriched Rosenbluth Method (nPERM) and flat histogram sampling of the
entire energy space. Since our method directly simulates the density of states,
we obtain results for thermodynamic quantities of the system for all
temperatures. In particular, this algorithm enables us to accurately simulate
the usually difficult accessible low-temperature region. Therefore, it becomes
possible to perform detailed analyses of the low-temperature transition between
ground states and compact globules.


  We study the dynamics of epidemic spreading processes aimed at spontaneous
dissemination of information updates in populations with complex connectivity
patterns. The influence of the topological structure of the network in these
processes is studied by analyzing the behavior of several global parameters
such as reliability, efficiency and load. Large scale numerical simulations of
update spreading processes show that while networks with homogeneous
connectivity patterns permit a higher reliability, scale-free topologies allow
for a better efficiency.


  Logistic growth process with nonlocal interactions is considered in one
dimension. Spontaneous breakdown of translational invariance is shown to take
place at some parameter region, and the bifurcation regime is identified for
short and long range interactions. Domain walls between regions of different
order parameter are expressed as soliton solutions of the reduced dynamics for
nearest neighbor interactions. The analytic results are confirmed by numerical
simulations.


  Motivated by novel results in the theory of correlated sequences, we analyze
the dynamics of random walks with long-term memory (binary chains with
long-range correlations). In our model, the probability for a unit bit in a
binary string depends on the fraction of unities preceding it. We show that the
system undergoes a dynamical phase-transition from normal diffusion, in which
the variance D_L scales as the string's length L, into a super-diffusion phase
(D_L ~ L^{1+|alpha|}), when the correlation strength exceeds a critical value.
We demonstrate the generality of our results with respect to alternative
models, and discuss their applicability to various data, such as coarse-grained
DNA sequences, written texts, and financial data.


  We study the effect of the connectivity pattern of complex networks on the
propagation dynamics of epidemics. The growth time scale of outbreaks is
inversely proportional to the network degree fluctuations, signaling that
epidemics spread almost instantaneously in networks with scale-free degree
distributions. This feature is associated with an epidemic propagation that
follows a precise hierarchical dynamics. Once the highly connected hubs are
reached, the infection pervades the network in a progressive cascade across
smaller degree classes. The present results are relevant for the development of
adaptive containment strategies.


  Recent studies show that many types of human social activities, from
scientific collaborations to sexual contacts, can be understood in terms of
complex network of interactions. Such networking paradigm allows to model many
aspects of social behaviour with relatively simple computer models. The present
work investigates the influence of single leaders on opinion formation within
simulations of agent based artificial networked societies. Several types of
network systems (among them random networks, highly clustered, small world and
scale-free) are studied. The strength of the social influence of individuals is
assumed to be dictated by distance from an agent to another, as well as
individual strengths of the agents. We study the effect of different topologies
on the conditions of general acceptance of leader's opinion by the society.


  Grand canonical Monte Carlo simulations are used to explore the metastable
fluid-fluid coexistence curve of the modified Lennard-Jones model of globular
proteins of ten Wolde and Frenkel (Science, v277, 1975 (1997)). Using both
mixed-field finite-size scaling and histogram reweighting methods, the joint
distribution of density and energy fluctuations is analyzed at coexistence to
accurately determine the critical-point parameters. The subcritical coexistence
region is explored using the recently developed hyper-parallel tempering Monte
Carlo simulation method along with histogram reweighting to obtain the density
distributions. The phase diagram for the metastable fluid-fluid coexistence
curve is calculated in close proximity to the critical point, a region
previously unattained by simulation.


  The work investigates the influence of leader's strategy on opinion formation
in artificial networked societies. The strength of the social influence is
assumed to be dictated by distance from one agent to another, as well as
individual strengths of the agents. The leader is assumed to have much greater
resources, which allows him to tune the way he influences the other agents. We
study various strategies of using these resources to optimize the conditions
needed to `convince' the whole society to leader's opinion. The flexibility of
the model allows it to be used in studies of political, social and marketing
influence and opinion formation.


  We propose an alternative mechanism for the gating of biological membrane
channels in response to membrane tension that involves a change in the slope of
the membrane near the channel. Under biological membrane tensions we show that
the energy difference between the closed (tilted) and open (untilted) states
can far exceed kBT and is comparable to what is available under simple
ilational gating. Recent experiments demonstrate that membrane leaflet
asymmetries (spontaneous curvature) can strong effect the gating of some
channels. Such a phenomenon would be more easy to explain under gating-by-tilt,
given its novel intrinsic sensitivity to such asymmetry.


  In this work, we study the critical behavior of an epidemic propagation model
that considers individuals that can develop drug resistance. In our lattice
model, each site can be found in one of four states: empty, healthy, normally
infected (not drug resistant) and strain infected (drug resistant) states. The
most relevant parameters in our model are related to the mortality, cure and
mutation rates. This model presents two distinct stationary active phases: a
phase with co-existing normal and drug resistant infected individuals and an
intermediate active phase with only drug resistant individuals. We employ a
finite-size scaling analysis to compute the critical points the critical
exponents ratio $\beta/\nu$ governing the phase-transitions between these
active states and the absorbing inactive state. Our results are consistent with
the hypothesis that these transitions belong to the directed percolation
universality class.


  Random graphs with prescribed degree sequences have been widely used as a
model of complex networks. Comparing an observed network to an ensemble of such
graphs allows one to detect deviations from randomness in network properties.
Here we briefly review two existing methods for the generation of random graphs
with arbitrary degree sequences, which we call the ``switching'' and
``matching'' methods, and present a new method based on the ``go with the
winners'' Monte Carlo method. The matching method may suffer from nonuniform
sampling, while the switching method has no general theoretical bound on its
mixing time. The ``go with the winners'' method has neither of these drawbacks,
but is slow. It can however be used to evaluate the reliability of the other
two methods and, by doing this, we demonstrate that the deviations of the
switching and matching algorithms under realistic conditions are small compared
to the ``go with the winners'' algorithm. Because of its combination of speed
and accuracy we recommend the use of the switching method for most
calculations.


  A model of interacting random walkers is presented and shown to give rise to
patterns consisting in periodic arrangements of fluctuating particle clusters.
The model represents biological individuals that die or reproduce at rates
depending on the number of neighbors within a given distance. We evaluate the
importance of the discrete and fluctuating character of this particle model on
the pattern forming process. To this end, a deterministic mean-field
description, including a linear stability and a weakly nonlinear analysis, is
given and compared with the particle model. The deterministic approach is shown
to reproduce some of the features of the discrete description, in particular,
the existence of a finite-wavelength instability. Stochasticity in the particle
dynamics, however, has strong effects in other important aspects such as the
parameter values at which pattern formation occurs, or the nature of the
homogeneous phase.


  We propose that semiflexible polymers in poor solvent collapse in two stages.
The first stage is the well known formation of a dense toroidal aggregate.
However, if the solvent is sufficiently poor, the condensate will undergo a
second structural transition to a twisted entangled state, in which individual
filaments lower their bending energy by additionally orbiting around the mean
path along which they wind. This ``topological ripening'' is consistent with
known simulations and experimental results. It connects and rationalizes
various experimental observations ranging from strong DNA entanglement in viral
capsids to the unusually short pitch of the cholesteric phase of DNA in
sperm-heads. We propose that topological ripening of DNA toroids could improve
the efficiency and stability of gene delivery.


  Many biological and synthetic materials self-assemble into helical or twisted
aggregates. The shape is determined by a complex interplay between elastic
forces and the orientation and chirality of the constituent molecules. We study
this interplay through Monte Carlo simulations, with an accelerated algorithm
motivated by the growth of an aggregate out of solution. The simulations show
that the curvature changes smoothly from cylindrical to saddle-like as the
elastic moduli are varied. Remarkably, aggregates of either handedness form
from molecules of a single handedness, depending on the molecular orientation.


  We propose a new non-equilibrium model for spatial pattern formation on the
basis of local information transfer. Unlike standard models of pattern
formation it is not based on the Turing instability. Information is transmitted
through the system via particle-like excitations whose collective dynamics
result in pattern formation and control. Here, a simple problem of domain
formation is addressed by this model in an implementation as stochastic
cellular automata. One observes stable pattern formation, even in the presence
of noise and cell flow. Noise stabilizes the system through the production of
quasi-particles that control the position of the domain boundary.
Self-organized boundaries become sharp for large system with fluctuations
vanishing with a power of the system size. Pattern proportions are
scale-independent with system size. Pattern formation is stable over large
parameter ranges with two phase transitions occurring at vanishing noise and
increased cell flow.


  Cationic lipid membranes are known to form compact complexes with DNA and to
be effective as gene delivery agents both in vitro and in vivo. Here we employ
molecular dynamics simulations for a detailed atomistic study of lipid bilayers
consisting of a mixture of cationic dimyristoyltrimethylammonium propane
(DMTAP) and zwitterionic dimyristoylphosphatidylcholine (DMPC). Our main
objective is to examine how the composition of the bilayers affects their
structural and electrostatic properties in the liquid-crystalline phase. By
varying the mole fraction of DMTAP, we have found that the area per lipid has a
pronounced non-monotonic dependence on the DMTAP concentration, with a minimum
around the point of equimolar mixture. We show that this behavior has an
electrostatic origin and is driven by the interplay between positively charged
TAP headgroups and the zwitterionic PC heads. This interplay leads to
considerable re-orientation of PC headgroups for an increasing DMTAP
concentration, and gives rise to major changes in the electrostatic properties
of the lipid bilayer, including a significant increase of total dipole
potential across the bilayer and prominent changes in the ordering of water in
the vicinity of the membrane. Moreover, chloride counter-ions are bound mostly
to PC nitrogens implying stronger screening of PC heads by Cl ions compared to
TAP head groups. The implications of these findings are briefly discussed.


  Escherichia coli and other bacteria use rotating helical filaments to swim.
Each cell typically has about four filaments, which bundle or disperse
depending on the sense of motor rotation. To study the bundling process, we
built a macroscopic scale model consisting of stepper-motor-driven polymer
helices in a tank filled with a high-viscosity silicone oil. The Reynolds
number, the ratio of viscous to elastic stresses, and the helix geometry of our
experimental model approximately match the corresponding quantities of the full
scale E. coli cells. We analyze digital video images of the rotating helices to
show that the initial rate of bundling is proportional to the motor frequency
and is independent of the characteristic relaxation time of the filament. We
also determine which combinations of helix handedness and sense of motor
rotation lead to bundling.


  Certain bacteria form filamentous colonies when the cells fail to separate
after dividing. In Bacillus subtilis, Bacillus thermus, and cyanobacteria, the
filaments can wrap into complex supercoiled structures as the cells grow. The
structures may be solenoids or plectonemes, with or without branches in the
latter case. Any microscopic theory of these morphological instabilities must
address the nature of pattern selection in the presence of growth, for growth
renders the problem nonautonomous and the bifurcations dynamic. To gain insight
into these phenomena, we formulate a general theory for growing elastic
filaments with bending and twisting resistance in a viscous medium, and study
an illustrative model problem: a growing filament with preferred twist, closed
into a loop. Growth depletes the twist, inducing a twist strain. The closure of
the loop prevents the filament from unwinding back to the preferred twist;
instead, twist relaxation is accomplished by the formation of supercoils.
Growth also produces viscous stresses on the filament which even in the absence
of twist produce buckling instabilities. Our linear stability analysis and
numerical studies reveal two dynamic regimes. For small intrinsic twist the
instability is akin to Euler buckling, leading to solenoidal structures, while
for large twist it is like the classic writhing of a twisted filament,
producing plectonemic windings. This model may apply to situations in which
supercoils form only, or more readily, when axial rotation of filaments is
blocked. Applications to specific biological systems are proposed.


  We study the aggregation of peptides using the discrete molecular dynamics
simulations. At temperatures above the alpha-helix melting temperature of a
single peptide, the model peptides aggregate into a multi-layer parallel
beta-sheet structure. This structure has an inter-strand distance of 0.48 nm
and an inter-sheet distance of 1.0 nm, which agree with experimental
observations. In this model, the hydrogen bond interactions give rise to the
inter-strand spacing in beta-sheets, while the Go interactions among side
chains make beta-strands parallel to each other and allow beta-sheets to pack
into layers. The aggregates also contain free edges which may allow for further
aggregation of model peptides to form elongated fibrils.


  We have performed parallel tempering Monte Carlo simulations using a simple
continuum heteropolymer model for proteins. All ten heteropolymer sequences
which we have studied have shown first-order transitions at low temperature to
ordered states dominated by single chain conformations. These results are in
contrast with the theoretical predictions of the random energy model for
heteropolymers, from which we would expect continuous transitions to glassy
behavior at low temperatures.


  We give a theoretical treatment of the interaction of electronic excitations
(excitons) in biomolecules and quantum dots with the surrounding polar solvent.
Significant quantum decoherence occurs due to the interaction of the electric
dipole moment of the solute with the fluctuating electric dipole moments of the
individual molecules in the solvent. We introduce spin boson models which could
be used to describe the effects of decoherence on the quantum dynamics of
biomolecules which undergo light-induced conformational change and on
biomolecules or quantum dots which are coupled by Forster resonant energy
transfer.


  We investigate the orientational order of transverse polarization vectors of
long columns of nucleosomal core particles and their coupling to positional
order in high density mesophases discovered recently. Inhomogeneous polar
ordering of these columns precipitates crystalization of the 2D sections with
different orientations of the transverse polarization vector on each column in
the unit cell. We propose possible scenarios for going from the 2D hexagonal
phase into distorted lamellar and related phases observed experimentally.


  We show, both numerically and analytically, that the consequences of 'wiring
up' a competitive population depend quite dramatically on the interplay between
the local connectivity and the global resources. With modest global resources,
adding small amounts of local connectivity yields an increasingly heterogeneous
population. With substantial global resources, high-performing yet reasonably
homogenous collective states emerge instead.


  We solve the stochastic equations for a cluster of parallel bonds with shared
constant loading, rebinding and the completely dissociated state as an
absorbing boundary. In the small force regime, cluster lifetime grows only
logarithmically with bond number for weak rebinding, but exponentially for
strong rebinding. Therefore rebinding is essential to ensure physiological
lifetimes. The number of bonds decays exponentially with time for most cases,
but in the intermediate force regime, a small increase in loading can lead to
much faster decay. This effect might be used by cell-matrix adhesions to induce
signaling events through cytoskeletal loading.


  The outbreak of severe acute respiratory syndrome (SARS) is still threatening
the world because of a possible resurgence. In the current situation that
effective medical treatments such as antiviral drugs are not discovered yet,
dynamical features of the epidemics should be clarified for establishing
strategies for tracing, quarantine, isolation, and regulating social behavior
of the public at appropriate costs. Here we propose a network model for SARS
epidemics and discuss why superspreaders emerged and why SARS spread especially
in hospitals, which were key factors of the recent outbreak. We suggest that
superspreaders are biologically contagious patients, and they may amplify the
spreads by going to potentially contagious places such as hospitals. To avoid
mass transmission in hospitals, it may be a good measure to treat suspected
cases without hospitalizing them. Finally, we indicate that SARS probably
propagates in small-world networks associated with human contacts and that the
biological nature of individuals and social group properties are factors more
important than the heterogeneous rates of social contacts among individuals.
This is in marked contrast with epidemics of sexually transmitted diseases or
computer viruses to which scale-free network models often apply.


  We have studied the collective short wavelength dynamics in deuterated DMPC
bilayers by inelastic neutron scattering. The corresponding dispersion relation
$\hbar\omega$(Q) is presented for the gel and fluid phase of this model system.
The temperature dependence of the inelastic excitations indicates a phase
coexistence between the two phases over a broad range and leads to a different
assignment of excitations than that reported in a preceding inelastic x-ray
scattering study [Phys. Rev. Lett. {\bf 86}, 740 (2001)]. As a consequence, we
find that the minimum in the dispersion relation is actually deeper in the gel
than in the fluid phase. Finally, we can clearly identify an additional
non-dispersive (optical) mode predicted by Molecular Dynamics (MD) simulations
[Phys. Rev. Lett. {\bf 87}, 238101 (2001)].


  We present a very simple model for the study of charge transport in a
molecule patterned on B-DNA. In this model we use a discrete non-linear
Schr\"{o}dinger equation to describe electrons propagating along the
sugar-phosphate backbone of the DNA molecule. We find that in this model, for a
given nonlinearity, the transport is controlled by $J$, a parameter which
relates to the electronic coupling between different molecules on the backbone.
For smaller values of $J$ we have localised states while at higher values of
$J$ the soliton field is spread out and through its interaction with the
lattice it has stronger effects on the distortion of the lattice.


  We investigate saturation effects in susceptible-infected-susceptible (SIS)
models of the spread of epidemics in heterogeneous populations. The structure
of interactions in the population is represented by networks with connectivity
distribution $P(k)$,including scale-free(SF) networks with power law
distributions $P(k)\sim k^{-\gamma}$. Considering cases where the transmission
of infection between nodes depends on their connectivity, we introduce a
saturation function $C(k)$ which reduces the infection transmission rate
$\lambda$ across an edge going from a node with high connectivity $k$. A mean
field approximation with the neglect of degree-degree correlation then leads to
a finite threshold $\lambda_{c}>0$ for SF networks with $2<\gamma \leq 3$. We
also find, in this approximation, the fraction of infected individuals among
those with degree $k$ for $\lambda$ close to $\lambda_{c}$. We investigate via
computer simulation the contact process on a heterogeneous regular lattice and
compare the results with those obtained from mean field theory with and without
neglect of degree-degree correlations.


  We present the first inelastic neutron scattering study of the short
wavelength dynamics in a phospholipid bilayer. We show that inelastic neutron
scattering using a triple-axis spectrometer at the high flux reactor of the ILL
yields the necessary resolution and signal to determine the dynamics of model
membranes. The results can quantitatively be compared to recent Molecular
Dynamics simulations. Reflectivity, in-plane correlations and the corresponding
dynamics can be measured simultaneously to gain a maximum amount of
information. With this method, dispersion relations can be measured with a high
energy resolution. Structure and dynamics in phospholipid bilayers, and the
relation between them, can be studied on a molecular length scale.


  Systems with two species of active molecular motors moving on (cytoskeletal)
filaments into opposite directions are studied theoretically using driven
lattice gas models. The motors can unbind from and rebind to the filaments. Two
motors are more likely to bind on adjacent filament sites if they belong to the
same species. These systems exhibit (i) Continuous phase transitions towards
states with spontaneously broken symmetry, where one motor species is largely
excluded from the filament, (ii) Hysteresis of the total current upon varying
the relative concentrations of the two motor species, and (iii) Coexistence of
traffic lanes with opposite directionality in multi-filament systems. These
theoretical predictions should be experimentally accessible.


  The growth dynamics of rigid biopolymers, consisting of $N$ parallel
protofilaments, is investigated theoretically using simple approximate models.
In our approach, the structure of a polymer's growing end and lateral
interactions between protofilaments are explicitly taken into account, and it
is argued that only few conformations are important for biopolymer's growth. As
a result, exact analytic expressions for growth velocity and dispersion are
obtained for {\it any} number of protofilaments and arbitrary geometry of the
growing end of the biopolymer. Our theoretical predictions are compared with a
full description of biopolymer growth dynamics for the simplest N=2 model. It
is found that the results from the approximate theory are approaching the exact
ones for large lateral interactions between the protofilaments. Our theory is
also applied to analyze the experimental data on the growth of microtubules.


  The paper studies the emergence and stability of cooperative behavior in
populations of agents who interact among themselves in Prisoner's Dilemma games
and who are allowed to choose their partners. The population is then subject to
evolutionary pressures, based on individual payoffs from multiple encounters. A
simple formula for signaling and recognition is introduced, which allows the
agents to guess prospective partner strategy and to refuse the game if the
predicted outcome is unfavorable. We present both algebraic formulation of the
average payoffs and results of computer simulations of evolution of such
society. The simulations result in surprising variety of behavior. We discuss
possible interpretations of the results as well as relationship between
artificial `computer time' and real time of possible social and biological
systems that the model might be applied to.


  In an effort to understand ``rafts'' in biological membranes, we propose
phenomenological models for saturated and unsaturated lipid mixtures, and
lipid-cholesterol mixtures. We consider simple couplings between the local
composition and internal membrane structure, and their influence on transitions
between liquid and ``gel'' membrane phases. Assuming that the gel transition
temperature of the saturated lipid is shifted by the presence of the
unsaturated lipid, and that cholesterol acts as an external field on the chain
melting transition, a variety of phase diagrams are obtained. The phase
diagrams for binary mixtures of saturated/unsaturated lipids and
lipid/cholesterol are in semi-quantitative agreement with the experiments. Our
results also apply to regions in the ternary phase diagram of
lipid/lipid/cholesterol systems.


  Single molecule pulling experiments provide information about interactions in
biomolecules that cannot be obtained by any other method. However, the
reconstruction of the molecule's free energy profile from the experimental data
is still a challenge, in particular for the unstable barrier regions. We
propose a new method for obtaining the full profile by introducing a periodic
ramp and using Jarzynski's identity for obtaining equilibrium quantities from
non-equilibrium data. Our simulated experiments show that this method delivers
significant more accurate data than previous methods, under the constraint of
equal experimental effort.


  In this paper we present a general method for information extraction that
exploits the features of data compression techniques. We first define and focus
our attention on the so-called "dictionary" of a sequence. Dictionaries are
intrinsically interesting and a study of their features can be of great
usefulness to investigate the properties of the sequences they have been
extracted from (e.g. DNA strings). We then describe a procedure of string
comparison between dictionary-created sequences (or "artificial texts") that
gives very good results in several contexts. We finally present some results on
self-consistent classification problems.


  An ensemble of Brownian particles in a feedback controlled flashing ratchet
is studied. The ratchet potential is switched on and off depending on the
position of the particles, with the aim of maximizing the current. We study in
detail a protocol which maximizes the instant velocity of the center of mass of
the ensemble at any time. This protocol is optimal for one particle and
performs better than any periodic flashing for ensembles of moderate size,
whereas is defeated by a random or periodic switching for large ensembles.


  Motivated by experiments in which single DNA molecules are stretched and
twisted we consider a perturbative approach around very high forces, where we
determine the writhe distribution in a simple, analytically tractable model.
Our results are in agreement with recent simulations and experiments.


  A DNA-protein complex modelled by a semiflexible chain and an attractive
spherical core is studied in the situation when an external stretching force is
acting on one end monomer of the chain while the other end monomer is kept
fixed in space. Without stretching force, the chain is wrapped around the core.
By applying an external stretching force, unwrapping of the complex is induced.
We study the statics and the dynamics of the unwrapping process by computer
simulation and simple phenomenological theory. We find two different scenarios
depending on the chain stiffness: For a flexible chain, the extension of the
complex scales linearly with the external force applied. The sphere-chain
complex is disordered, i.e. there is no clear winding of the chain around the
sphere. For a stiff chain, on the other hand, the complex structure is ordered,
which is reminiscent to nucleosome. There is a clear winding number and the
unwrapping process under external stretching is discontinuous with jumps of the
distance-force curve. This is associated to discrete unwinding processes of the
complex. Our predictions are of relevance for experiments, which measure
force-extension curves of DNA-protein complexes, such as nucleosome, using
optical tweezers.


  Numerical studies of the average size of trivially knotted polymer loops with
no excluded volume are undertaken. Topology is identified by Alexander and
Vassiliev degree 2 invariants. Probability of a trivial knot, average gyration
radius, and probability density distributions as functions of gyration radius
are generated for loops of up to N=3000 segments. Gyration radii of trivially
knotted loops are found to follow a power law similar to that of self avoiding
walks consistent with earlier theoretical predictions.


  We study the cooperative rupture of multiple adhesion bonds under shared
linear loading. Simulations of the appropriate Master equation are compared
with numerical integration of a rate equation for the mean number of bonds and
its scaling analysis. In general, force-accelerated rupture is rather abrupt.
For small clusters and slow loading, large fluctuations occur regarding the
timepoint of final rupture, but not the typical shape of the rupture
trajectory. For vanishing rebinding, our numerical results confirm three
scaling regimes predicted before for cluster lifetime as a function of loading
rate. For finite rebinding, the intermediate loading regime becomes irrelevant,
and a sequence of two new scaling laws can be identified in the slow loading
regime.


  Charge migration along DNA molecules has attracted scientific interest for
over half a century. Reports on possible high rates of charge transfer between
donor and acceptor through the DNA, obtained in the last decade from solution
chemistry experiments on large numbers of molecules, triggered a series of
direct electrical transport measurements through DNA single molecules, bundles
and networks. These measurements are reviewed and presented here. From these
experiments we conclude that electrical transport is feasible in short DNA
molecules, in bundles and networks, but blocked in long single molecules that
are attached to surfaces. The experimental background is complemented by an
account of the theoretical/computational schemes that are applied to study the
electronic and transport properties of DNA-based nanowires. Examples of
selected applications are given, to show the capabilities and limits of current
theoretical approaches to accurately describe the wires, interpret the
transport measurements, and predict suitable strategies to enhance the
conductivity of DNA nanostructures.


  Models of contagion arise broadly both in the biological and social sciences,
with applications ranging from the transmission of infectious diseases to the
diffusion of innovations and the spread of cultural fads. In this Letter, we
introduce a general model of contagion which, by explicitly incorporating
memory of past exposures to, for example, an infectious agent, rumor, or new
product, includes the main features of existing contagion models and
interpolates between them. We obtain exact solutions for a simple version of
the model, finding that under general conditions only three classes of
collective dynamics exist, two of which correspond to familiar epidemic
threshold and critical mass dynamics, while the third is a distinct
intermediate case. We find that for a given length of memory, the class into
which a particular system falls is determined by two parameters, each of which
ought to be measurable empirically. Our model suggests novel measures for
assessing the susceptibility of a population to large contagion events, and
also a possible strategy for inhibiting or facilitating them.


  We propose that the first stage of membrane fusion need not be the formation
of a single stalk. Instead, we consider a scenario for defect-induced membrane
fusion that proceeds cooperatively via multiple stalk formation. The defects
(stalks or pores) attract each other via membrane-mediated capillary
interactions that result in a condensation transition of the defects. The
resulting dense phase of stalks corresponds to the so-called fusion
intermediate.


  We study the unzipping of double stranded DNA (dsDNA) by applying a pulling
force at a fraction $s$ $(0 \le s \le 1)$ from the anchored end. From exact
analytical and numerical results, the complete phase diagram is presented. The
phase diagram shows a strong ensemble dependence for various values of $s$. In
addition, we show the existence of an ``eye'' phase and a triple point.


  We study the phase synchronization of Kuramoto's oscillators in small parts
of networks known as motifs. We first report on the system dynamics for the
case of a scale-free network and show the existence of a non-trivial critical
point. We compute the probability that network motifs synchronize, and find
that the fitness for synchronization correlates well with motif's
interconnectedness and structural complexity. Possible implications for present
debates about network evolution in biological and other systems are discussed.


  The renowned general epidemic process describes the stochastic evolution of a
population of individuals which are either susceptible, infected or dead. A
second order phase transition belonging to the universality class of dynamic
isotropic percolation lies between endemic or pandemic behavior of the process.
We generalize the general epidemic process by introducing a fourth kind of
individuals, viz. individuals which are weakened by the process but not yet
infected. This sensibilization gives rise to a mechanism that introduces a
global instability in the spreading of the process and therefore opens the
possibility of a discontinuous transition in addition to the usual continuous
percolation transition. The tricritical point separating the lines of first and
second order transitions constitutes a new universality class, namely the
universality class of tricritical dynamic isotropic percolation. Using
renormalized field theory we work out a detailed scaling description of this
universality class. We calculate the scaling exponents in an
$\epsilon$-expansion below the upper critical dimension $d_{c}=5$ for various
observables describing tricritical percolation clusters and their spreading
properties. In a remarkable contrast to the usual percolation transition, the
exponents $\beta$ and ${\beta}^{\prime}$ governing the two order parameters,
viz. the mean density and the percolation probability, turn out to be different
at the tricritical point. In addition to the scaling exponents we calculate for
all our static and dynamic observables logarithmic corrections to the
mean-field scaling behavior at $d_c=5$.


  A theory is set up of spherical proteins interacting by screened
electrostatics and constant adhesion, in which the effective adhesion parameter
is optimized by a variational principle for the free energy. An analytical
approach to the second virial coefficient is first outlined by balancing the
repulsive electrostatics against part of the bare adhesion. A theory similar in
spirit is developed at nonzero concentrations by assuming an appropriate Baxter
model as the reference state. The first-order term in a functional expansion of
the free energy is set equal to zero which determines the effective adhesion as
a function of salt and protein concentrations. The resulting theory is shown to
have fairly good predictive power for the ionic-strength dependence of both the
second virial coefficient and the osmotic pressure or compressibility of
lysozyme up to about 0.2 volume fraction.


  We discuss pore dynamics in osmotically stressed vesicles. A set of equations
which govern the liposomal size, internal solute concentration, and pore
diameter is solved numerically. We find that dependent on the internal solute
concentration and vesicle size, liposomes can stay pore-free, nucleate a short
lived pore, or nucleate a long-lived pore. The phase diagram of pore stability
is constructed, and the different scaling regimes are deduced analytically.


  Carbon nanotubes and biological filaments each spontaneously assemble into
kinked helices, rings, and "tennis racket" shapes due to competition between
elastic and interfacial effects. We show that the slender geometry is a more
important determinant of the morphology than any molecular details. Our
mesoscopic continuum theory is capable of quantifying observations of these
structures, and is suggestive of their occurrence in other filamentous
assemblies as well.


  1. The birth of network science. 2. What are random networks? 3. Adjacency
matrix. 4. Degree distribution. 5. What are simple networks? Classical random
graphs. 6. Birth of the giant component. 7. Topology of the Web. 8.Uncorrelated
networks. 9. What are small worlds? 10. Real networks are mesoscopic objects.
11. What are complex networks? 12. The configuration model. 13. The absence of
degree--degree correlations. 14.Networks with correlated degrees.15.Clustering.
16. What are small-world networks? 17. `Small worlds' is not the same as
`small-world networks'. 18. Fat-tailed degree distributions. 19.Reasons for the
fat-tailed degree distributions. 20. Preferential linking. 21. Condensation of
edges. 22. Cut-offs of degree distributions. 23. Reasons for correlations in
networks. 24. Classical random graphs cannot be used for comparison with real
networks. 25. How to measure degree--degree correlations. 26. Assortative and
disassortative mixing. 27. Disassortative mixing does not mean that vertices of
high degrees rarely connect to each other. 28. Reciprocal links in directed
nets. 29. Ultra-small-world effect. 30. Tree ansatz. 31.Ultraresilience against
random failures. 32. When correlated nets are ultraresilient. 33. Vulnerability
of complex networks. 34. The absence of an epidemic threshold. 35. Search based
on local information. 36.Ultraresilience disappears in finite nets. 37.Critical
behavior of cooperative models on networks. 38. Berezinskii-Kosterlitz-Thouless
phase transitions in networks. 39.Cascading failures. 40.Cliques & communities.
41. Betweenness. 42.Extracting communities. 43. Optimal paths. 44.Distributions
of the shortest-path length & of the loop's length are narrow. 45. Diffusion on
networks. 46. What is modularity? 47.Hierarchical organization of networks. 48.
Convincing modelling of real-world networks:Is it possible? 49. The small Web..


  Recent single-molecule pulling experiments have shown how it is possible to
manipulate RNA molecules using optical tweezers force microscopy. We
investigate a minimal model for the experimental setup which includes a RNA
molecule connected to two polymers (handles) and a bead, trapped in the optical
potential, attached to one of the handles. Initially, we focus on small
single-domain RNA molecules which unfold in a cooperative way. The model
qualitatively reproduces the experimental results and allow us to investigate
the influence of the bead and handles on the unfolding reaction. A main
ingredient of our model is to consider the appropriate statistical ensemble and
the corresponding thermodynamic potential describing thermal fluctuations in
the system. We then investigate several questions relevant to extract
thermodynamic information from the experimental data. Next, we study the
kinetics using a dynamical model. Finally, we address the more general problem
of a multidomain RNA molecule with Mg2+-tertiary contacts that unfolds in a
sequential way and propose techniques to analyze the breakage force data in
order to obtain the reliable kinetics parameters that characterize each domain.


  Current all-atom potential based molecular dynamics (MD) allow the
identification of a protein's functional motions on a wide-range of
time-scales, up to few tens of ns. However, functional large scale motions of
proteins may occur on a time-scale currently not accessible by all-atom
potential based molecular dynamics. To avoid the massive computational effort
required by this approach several simplified schemes have been introduced. One
of the most satisfactory is the Gaussian Network approach based on the energy
expansion in terms of the deviation of the protein backbone from its native
configuration. Here we consider an extension of this model which captures in a
more realistic way the distribution of native interactions due to the
introduction of effective sidechain centroids. Since their location is entirely
determined by the protein backbone, the model is amenable to the same exact and
computationally efficient treatment as previous simpler models. The ability of
the model to describe the correlated motion of protein residues in
thermodynamic equilibrium is established through a series of successful
comparisons with an extensive (14 ns) MD simulation based on the AMBER
potential of HIV-1 protease in complex with a peptide substrate. Thus, the
model presented here emerges as a powerful tool to provide preliminary, fast
yet accurate characterizations of proteins near-native motion.


  The problem of finite-dimensional asymptotics of infinite-dimensional dynamic
systems is studied. A non-linear kinetic system with conservation of supports
for distributions has generically finite-dimensional asymptotics. Such systems
are apparent in many areas of biology, physics (the theory of parametric wave
interaction), chemistry and economics. This conservation of support has a
biological interpretation: inheritance. The finite-dimensional asymptotics
demonstrates effects of "natural" selection. Estimations of the asymptotic
dimension are presented. After some initial time, solution of a kinetic
equation with conservation of support becomes a finite set of narrow peaks that
become increasingly narrow over time and move increasingly slowly. It is
possible that these peaks do not tend to fixed positions, and the path covered
tends to infinity as t goes to infinity. The drift equations for peak motion
are obtained. Various types of distribution stability are studied: internal
stability (stability with respect to perturbations that do not extend the
support), external stability or uninvadability (stability with respect to
strongly small perturbations that extend the support), and stable realizability
(stability with respect to small shifts and extensions of the density peaks).
Models of self-synchronization of cell division are studied, as an example of
selection in systems with additional symmetry. Appropriate construction of the
notion of typicalness in infinite-dimensional space is discussed, and the
notion of "completely thin" sets is introduced.
  Key words: Dynamics; Attractor; Evolution; Entropy; Natural selection


  Unlike most synthetic materials, biological materials often stiffen as they
are deformed. This nonlinear elastic response, critical for the physiological
function of some tissues, has been documented since at least the 19th century,
but the molecular structure and the design principles responsible for it are
unknown. Current models for this response require geometrically complex ordered
structures unique to each material. In this Article we show that a much simpler
molecular theory accounts for strain stiffening in a wide range of molecularly
distinct biopolymer gels formed from purified cytoskeletal and extracellular
proteins. This theory shows that systems of semi-flexible chains such as
filamentous proteins arranged in an open crosslinked meshwork invariably
stiffen at low strains without the need for a specific architecture or multiple
elements with different intrinsic stiffnesses.


  Under constant applied force, the separation of double-stranded DNA into two
single strands is known to proceed through a series of pauses and jumps. Given
experimental traces of constant-force unzipping, we present a method whereby
the locations of pause points can be extracted in the form of a pause point
spectrum. A simple theoretical model of DNA constant-force unzipping is
demonstrated to produce good agreement with the experimental pause point
spectrum of lambda phage DNA. The locations of peaks in the experimental and
theoretical pause point spectra are found to be nearly coincident below 6000
bp. The model only requires the sequence, temperature and a set of empirical
base pair binding and stacking energy parameters, and the good agreement with
experiment suggests that pause points are primarily determined by the DNA
sequence. The model is also used to predict pause point spectra for the
BacterioPhage PhiX174 genome. The algorithm for extracting the pause point
spectrum might also be useful for studying related systems which exhibit
pausing behavior such as molecular motors.


  We show that the exponential length distribution that is typical of actin
filaments under physiological conditions dramatically narrows in the presence
of (i) crosslinker proteins (ii) polyvalent counterions or (iii) depletion
mediated attractions. A simple theoretical model shows that in equilibrium,
short-range attractions enhance the tendency of filaments to align parallel to
each other, eventually leading to an increase in the average filament length
and a decrease in the relative width of the distribution of filament lengths.


  In this paper we develop a theory to describe innovation processes in a
network of interacting units. We introduce a stochastic picture that allows for
the clarification of the role of fluctuations for the survival of innovations
in such a non-linear system. We refer to the theory of complex networks and
introduce the notion of sensitive networks. Sensitive networks are networks in
which the introduction or the removal of a node/vertex dramatically changes the
dynamic structure of the system. As an application we consider interaction
networks of firms and technologies and describe technological innovation as a
specific dynamic process. Random graph theory, percolation, master equation
formalism and the theory of birth and death processes are the mathematical
instruments used in this paper.


  Kinesins are processive motor proteins that move along microtubules in a
stepwise manner, and their motion is powered by the hydrolysis of ATP. Recent
experiments have investigated the coupling between the individual steps of
single kinesin molecules and ATP hydrolysis, taking explicitly into account
forward steps, backward steps and detachments. A theoretical study of
mechanochemical coupling in kinesins, which extends the approach used
successfully to describe the dynamics of conventional motor proteins, is
presented. The possibility of irreversible detachments of kinesins from the
microtubules is also explicitly taken into account. Using the method of first-
passage times, experimental data on the mechanochemical coupling in kinesins
are fully described using the simplest two-state model. It is shown that the
dwell times for the kinesin to move one step forward or backward, or to
dissociate irreversibly are the same, although the probabilities of these
events are different. It is concluded that the current theoretical view, that
only the forward motion of the motor protein molecule is coupled to ATP
hydrolysis, is consistent with all available experimental observations for
kinesins.


  Signaling pathways and networks determine the ability to communicate in
systems ranging from living cells to human society. We investigate how the
network structure constrains communication in social-, man-made and biological
networks. We find that human networks of governance and collaboration are
predictable on teat-a-teat level, reflecting well defined pathways, but
globally inefficient. In contrast, the Internet tends to have better overall
communication abilities, more alternative pathways, and is therefore more
robust. Between these extremes the molecular network of Saccharomyces cerevisea
is more similar to the simpler social systems, whereas the pattern of
interactions in the more complex Drosophilia melanogaster, resembles the robust
Internet.


  The use of reduced models for investigating the self-assembly dynamics
underlying protein shell formation in spherical viruses is described. The
spontaneous self-assembly of these polyhedral, supramolecular structures, in
which icosahedral symmetry is a conspicuous feature, is a phenomenon whose
dynamics remain unexplored; studying the growth process by means of computer
simulation provides access to the mechanisms underlying assembly. In order to
capture the more universal aspects of self-assembly, namely the manner in which
component shapes influence structure and assembly pathway, in this exploratory
study low-resolution approximations are used to represent the basic protein
building blocks. Alternative approaches involving both irreversible and
reversible assembly are discussed, models based on both schemes are introduced,
and examples of the resulting behavior described.


  We use the worm-like chain model to analytically calculate the linear
response of a grafted semiflexible polymer to a uniform force field. The result
is a function of the bending stiffness, the temperature, the total contour
length, and the orientation of the field with respect to that of the grafted
end. We also study the linear response of a worm-like chain with a periodic
alternating sequence of positive and negative charges. This can be considered
as a model for a polyampholyte with intrinsic bending siffness and negligible
intramolecular interactions. We show how the finite intrinsic persistence
length affects the linear response to the external field.


  We study a solution of long polyanions (PA) with shorter polycations (PC) and
focus on the role of Coulomb interaction. A good example is solutions of DNA
and PC which are widely studied for gene therapy. In the solution, each PA
attracts many PCs to form a complex. When the ratio of total charges of PA and
PC in the solution, $x$, equals to 1, complexes are neutral and they condense
in a macroscopic drop. When $x$ is far away from 1, complexes are strongly
charged. The Coulomb repulsion is large and free complexes are stable. As $x$
approaches to 1, PCs attached to PA disproportionate themselves in two
competing ways. One way is inter-complex disproportionation, in which PCs make
some complexes neutral and therefore condensed in a macroscopic drop while
other complexes become even stronger charged and stay free. The other way is
intra-complex disproportionation, in which PCs make one end of a complex
neutral and condensed in a small droplet while the rest of the complex forms a
strongly charged tail. Thus each complex becomes a "tadpole". These two ways
can also combine together to give even lower free energy. We get a phase
diagram of PA-PC solution in a plane of $x$ and inverse screening radius of the
monovalent salt, which includes phases or phase coexistence with both kinds of
disproportionation.


  The importance of nonlinearities in material constitutive relations has long
been appreciated in the continuum mechanics of macroscopic rods. Although the
moment (torque) response to bending is almost universally linear for small
deflection angles, many rod systems exhibit a high-curvature softening. The
signature behavior of these rod systems is a kinking transition in which the
bending is localized. Recent DNA cyclization experiments by Cloutier and Widom
have offered evidence that the linear-elastic bending theory fails to describe
the high-curvature mechanics of DNA. Motivated by this recent experimental
work, we develop a simple and exact theory of the statistical mechanics of
linear-elastic polymer chains that can undergo a kinking transition. We
characterize the kinking behavior with a single parameter and show that the
resulting theory reproduces both the low-curvature linear-elastic behavior
which is already well described by the Wormlike Chain model, as well as the
high-curvature softening observed in recent cyclization experiments.


  The dynamics of molecular motors which occasionally detach from a
heterogeneous track like DNA or RNA is considered. Motivated by recent single
molecule experiments, we study a simple model for a motor moving along a
disordered track using chemical energy while an external force opposes its
motion. The motors also have finite processivity, i.e., they can leave the
track with a position dependent rate. We show that the response of the system
to disorder in the hopping off rate depends on the value of the external force.
For most values of the external force, strong disorder causes the motors which
survive for long times on the track to be localized at preferred positions.
However, near the stall force, localization occurs for any amount of disorder.
Existence of localized states near the top of the band implies a stretched
exponential contribution to the decay of the survival probability. To obtain
these results, we study the complex eigenvalue spectrum of the time evolution
operator. A similar spectral analysis also provides a very efficient method for
studying the dynamics of motors with infinite processivity.


  The transient behavior of an ecosystem with N random interacting species in
the presence of a multiplicative noise is analyzed. The multiplicative noise
mimics the interaction with the environment. We investigate different
asymptotic dynamical regimes and the role of the external noise on the
probability distribution of the local field.


  We introduce a model of proteins in which all of the key atoms in the protein
backbone are accounted for, thus extending the Freely Rotating Chain model. We
use average bond lengths and average angles from the Protein Databank as input
parameters, leaving the number of residues as a single variable. The model is
used to study the stretching of proteins in the entropic regime. The results of
our Monte Carlo simulations are found to agree well with experimental data,
suggesting that the force extension plot is universal and does not depend on
the side chains or primary structure of proteins.


  We discuss similarities and differences between systems of interacting
players maximizing their individual payoffs and particles minimizing their
interaction energy. Long-run behavior of stochastic dynamics of spatial games
with multiple Nash equilibria is analyzed. In particular, we construct an
example of a spatial game with three strategies, where stochastic stability of
Nash equilibria depends on the number of players and the kind of dynamics.


  We study a physical system which, while devoid of the complexity one usually
associates with proteins, nevertheless displays a remarkable array of
protein-like properties. The constructive hypothesis that this striking
resemblance is not accidental leads not only to a unified framework for
understanding protein folding, amyloid formation and protein interactions but
also has implications for natural selection.


  Encapsulation of dsDNA fragments (contour length 54 nm) by the cationic
diblock copolymer poly(butadiene-b-N-methyl 4-vinyl pyridinium) [PBd-b-P4VPQ]
has been studied with phase contrast, polarized light, and fluorescence
microscopy, as well as scanning electron microscopy. Encapsulation was achieved
with a single emulsion technique. For this purpose, an aqueous DNA solution is
emulsified in an organic solvent (toluene) and stabilized by the amphiphilic
diblock copolymer. The PBd block forms an interfacial brush, whereas the
cationic P4VPQ block complexes with DNA. A subsequent change of the quality of
the organic solvent results in a collapse of the PBd brush and the formation of
a capsule. Inside the capsules, the DNA is compacted as shown by the appearance
of birefringent textures under crossed polarizers and the increase in
fluorescence intensity of labeled DNA. The capsules can also be dispersed in
aqueous medium to form vesicles, provided they are stabilized with an osmotic
agent (polyethylene glycol) in the external phase. It is shown that the DNA is
released from the vesicles once the osmotic pressure drops below 105 N/m2 or if
the ionic strength of the supporting medium exceeds 0.1 M. The method has also
proven to be efficient to encapsulate pUC18 plasmid in sub-micron sized
vesicles and the general applicability of the method has been demonstrated by
the preparation of the charge inverse system: cationic poly(ethylene imine)
encapsulated by the anionic diblock poly(styrene-b-acrylic acid).


  We study the elastic response of a worm-like polymer chain with reversible
kink-like structural defects. This is a generic model for (a) the
double-stranded DNA with sharp bends induced by binding of certain proteins,
and (b) effects of trans-gauche rotations in the backbone of the
single-stranded DNA. The problem is solved both analytically and numerically by
generalizing the well-known analogy to the Quantum Rotator. In the small
stretching force regime, we find that the persistence length is renormalized
due to the presence of the kinks. In the opposite regime, the response to the
strong stretching is determined solely by the bare persistence length with
exponential corrections due to the ``ideal gas of kinks''. This high-force
behavior changes significantly in the limit of high bending rigidity of the
chain. In that case, the leading corrections to the mechanical response are
likely to be due to the formation of multi-kink structures, such as kink pairs.


  Electronic transport through DNA wires in the presence of a strong
dissipative environment is investigated. We show that new bath-induced
electronic states are formed within the bandgap. These states show up in the
linear conductance spectrum as a temperature dependent background and lead to a
crossover from tunneling to thermal activated behavior with increasing
temperature. Depending on the strength of the electron-bath coupling, the
conductance at the Fermi level can show a weak exponential or even an algebraic
length dependence. Our results suggest a new environmental-induced transport
mechanism. This might be relevant for the understanding of molecular conduction
experiments in liquid solution, like those recently performed on poly(GC)
oligomers in a water buffer (B. Xu et al., Nano Lett 4, 1105 (2004)).


  The polymerization dynamics of double-stranded polymers, such as actin
filaments, is investigated theoretically using simple chemical kinetic models
that explicitly take into account some microscopic details of the polymer
structure and the lateral interactions between the protofilaments. By
considering all possible molecular configurations, the exact analytical
expressions for the growth velocity and dispersion for two-stranded polymers
are obtained in the case of the growing at only one end, and for the growth
from both polymer ends. Exact theoretical calculations are compared with the
predictions of approximate multi-layer models that consider only a finite
number of the most relevant polymer configurations. Our theoretical approach is
applied to analyze the experimental data on the growth and fluctuations
dynamics of individual single actin filaments.


  We propose a statistical mechanics model for DNA melting in which base
stacking and pairing are explicitly introduced as distinct degrees of freedom.
Unlike previous approaches, this model describes thermal denaturation of DNA
secondary structure in the whole experimentally accessible temperature range.
Base pairing is described through a zipper model, base stacking through an
Ising model. We present experimental data on the unstacking transition,
obtained exploiting the observation that at moderately low pH this transition
is moved down to experimentally accessible temperatures. These measurements
confirm that the Ising model approach is indeed a good description of base
stacking. On the other hand, comparison with the experiments points to the
limitations of the simple zipper model description of base pairing.


  We propose a nonlinear model derived from first principles, to describe
bubble dynamics of DNA. Our model equations include a term derived from the
dissipative effect of intermolecular vibrational modes. Such modes are excited
by the propagating bubble, and we term it "curvature dissipation". The
equations we derive allow for stable pinned localized kinks which form the
bubble. We perform the stability analysis and specify the energy requirements
for the motion of the localized solutions. Our findings are consistent with
properties of DNA dynamics, and can be used in models for denaturation bubbles,
RNA and DNA transcription, nucleotide excision repair and meiotic
recombination.


  We obtain, using transfer matrix methods, the distribution function $P(R)$ of
the end-to-end distance, the loop formation probability and force-extension
relations in a model for short double-stranded DNA molecules. Accounting for
the appearance of ``bubbles'', localized regions of enhanced flexibility
associated with the opening of a few base pairs of double-stranded DNA in
thermal equilibrium, leads to dramatic changes in $P(R)$ and unusual
force-extension curves. An analytic formula for the loop formation probability
in the presence of bubbles is proposed. For short {\em heterogeneous} chains,
we demonstrate a strong dependence of loop formation probabilities on sequence,
as seen in recent experiments.


  We study a lattice model of a multi-species prey-predator system. Numerical
results show that for a small mutation rate the model develops irregular
long-period oscillatory behavior with sizeable changes in a number of species.
The periodicity of extinctions on Earth was suggested by Raup and Sepkoski but
so far is lacking a satisfactory explanation. Our model indicates that this is
a natural consequence of the ecosystem dynamics, not the result of any
extraterrestrial cause.


  The folding mechanism of the Villin headpiece (HP36) is studied by means of a
novel approach which entails an initial coarse-grained Monte Carlo (MC) scheme
followed by all-atom molecular dynamics (MD) simulations in explicit solvent.
The MC evolution occurs in a simplified free-energy landscape and allows an
efficient selection of marginally-compact structures which are taken as viable
initial conformations for the MD. The coarse-grained MC structural
representation is connected to the one with atomic resolution through a
``fine--graining'' reconstruction algorithm. This two-stage strategy is used to
select and follow the dynamics of seven different unrelated conformations of
HP36. In a notable case the MD trajectory rapidly evolves towards the folded
state, yielding a typical RMS deviation of the core region of only 2.4 A from
the closest NMR model (the typical RMSD over the whole structure being 4.0A).
The analysis of the various MC-MD trajectories provides valuable insight into
the details of the folding and mis-folding mechanisms and particularly about
the delicate influence of local and non-local interactions in steering the
folding process.


  We investigate the time evolution and stationary states of a stochastic,
spatially discrete, population model (contact process) with spatial
heterogeneity and imposed drift (wind) in one- and two-dimensions. We consider
in particular a situation in which space is divided into two regions: an oasis
and a desert (low and high death rates). Carrying out computer simulations we
find that the population in the (quasi) stationary state will be zero,
localized, or delocalized, depending on the values of the drift and other
parameters. The phase diagram is similar to that obtained by Nelson and
coworkers from a deterministic, spatially continuous model of a bacterial
population undergoing convection in a heterogeneous medium.


  We propose that supercoiling energizes double-stranded DNA (dsDNA) so as to
facilitate thermal fluctuations to an unzipped state. We support this with a
model of two elastic rods coupled via forces that represent base pair
interactions. Supercoiling is shown to lead to a spatially localized higher
energy state in a small region of dsDNA consisting of a few base pairs. This
causes the distance between specific base pairs to be extended, enhancing the
thermal probability for their disruption. Our theory permits the development of
an analogy between this unzipping transition and a second order phase
transition, for which the possibility of a new set of critical exponents is
identified.


  Molecular dynamics simulations are used to study the equilibrium distribution
of monovalent ions in a nanopore connecting two water reservoirs separated by a
membrane, both for the empty pore and that with a single stranded DNA molecule
inside. In the presence of DNA, the counterions condense on the stretched
macromolecule effectively neutralizing it, and nearly complete depletion of
coions from the pore is observed. The implications of our results for
experiments on DNA translocation through alpha-hemolysin nanopores are
discussed.


  The determination and classification of fixed points of large Boolean
networks is addressed in terms of constraint satisfaction problem. We develop a
general simplification scheme that, removing all those variables and functions
belonging to trivial logical cascades, returns the computational core of the
network. The onset of an easy-to-complex regulatory phase is introduced as a
function of the parameters of the model, identifying both theoretically and
algorithmically the relevant regulatory variables.


  We show that active transport processes in biological systems can be
understood through a local equilibrium description formulated at the mesoscale,
the scale to describe stochastic processes. This new approach uses the method
established by nonequilibrium thermodynamics to account for the irreversible
processes occurring at this scale and provides nonlinear kinetic equations for
the rates in terms of the driving forces. The results show that the application
domain of nonequilibrium thermodynamics method to biological systems goes
beyond the linear domain. A model for transport of Ca$^{2+}$ by the
Ca$^{2+}$-ATPase, coupled to the hydrolysis of adenosine-triphosphate is
analyzed in detail showing that it depends on the reaction Gibbs energy in a
non-linear way. Our results unify thermodynamic and kinetic descriptions,
thereby opening new perspectives in the study of different transport phenomena
in biological systems.


  We study the role of fluctuations in single molecule experimental
measurements of force-extension curves. We use the Worm Like Chain (WLC) model
to bring out the connection between the Helmholtz ensemble characterized by the
Free Energy and the Gibbs ensemble characterized by the Free Energy . We
consider the rigid rod limit of the WLC model as an instructive special case to
bring out the issue of ensemble inequivalence. We point out the need for taking
into account the free energy of transition when one goes from one ensemble to
another. We also comment on the ``phase transition'' noticed in an isometric
setup for semiflexible polymers and propose a realization of its thermodynamic
limit. We present general arguments which rule out non-monotonic
force-extension curves in some ensembles and note that these do not apply to
the isometric ensemble.


  Spider capture silk is a biomaterial with both high strength and high
elasticity, but the structural design principle underlying these remarkable
properties is still unknown. It was revealed recently by atomic force
microscopy that, an exponential force--extension relationship holds both for
capture silk mesostructures and for intact capture silk fibers [N. Becker et
al., Nature Materials 2, 278 (2003)]. In this Letter a simple hierarchical
chain model was proposed to understand and reproduce this striking observation.
In the hierarchical chain model, a polymer is composed of many structural
motifs which organize into structural modules and supra-modules in a
hierarchical manner. Each module in this hierarchy has its own characteristic
force. The repetitive patterns in the amino acid sequence of the major
flagelliform protein of spider capture silk is in support of this model.


  The onset of intermediate states (denaturation bubbles) and their role during
the melting transition of DNA are studied using the Peyrard-Bishop-Daxuois
model by Monte Carlo simulations with no adjustable parameters. Comparison is
made with previously published experimental results finding excellent
agreement. Melting curves, critical DNA segment length for stability of bubbles
and the possibility of a two states transition are studied.


  We present a study on the near equilibrium dynamics of two small proteins in
the family of truncated hemoglobins, developed under the framework of a
Gaussian network approach. Effective beta carbon atoms are taken into account
besides C-alphas for all residues but glycines in the coarse-graining
procedure, without leading to an increase in the degrees of freedom
(betaGaussian Model). Normalized covariance matrix and deformation along
slowest modes with collective character are analyzed, pointing out
anti-correlations between functionally relevant sites for the proteins under
study. In particular we underline the functional motions of an extended
tunnel-cavity system running inside the protein matrix, which provide a pathway
for small ligands binding with the iron in the heme group. We give a rough
estimate of the order of magnitude of the relaxation times of the slowest two
overdamped modes and compare results with previous studies on globins.


  Although the ``scale-free'' literature is large and growing, it gives neither
a precise definition of scale-free graphs nor rigorous proofs of many of their
claimed properties. In fact, it is easily shown that the existing theory has
many inherent contradictions and verifiably false claims. In this paper, we
propose a new, mathematically precise, and structural definition of the extent
to which a graph is scale-free, and prove a series of results that recover many
of the claimed properties while suggesting the potential for a rich and
interesting theory. With this definition, scale-free (or its opposite,
scale-rich) is closely related to other structural graph properties such as
various notions of self-similarity (or respectively, self-dissimilarity).
Scale-free graphs are also shown to be the likely outcome of random
construction processes, consistent with the heuristic definitions implicit in
existing random graph approaches. Our approach clarifies much of the confusion
surrounding the sensational qualitative claims in the scale-free literature,
and offers rigorous and quantitative alternatives.


  The aim of the present paper is the study of the entropic elasticity of the
dsDNA molecule, having a cristallographic length L of the order of 10 to 30
persistence lengths A, when it is subject to spatial obstructions. We have not
tried to obtain the single molecule partition function by solving a
Schodringer-like equation. We prefer to stay within a discretized version of
the WLC model with an added one-monomer potential, simulating the spatial
constraints. We derived directly from the discretized Boltzmann formula the
transfer matrix connecting the partition functions relative to adjacent
"effective monomers". We have plugged adequate Dirac delta-functions in the
functional integral to ensure that the monomer coordinate and the tangent
vector are independent variables. The partition function is, then, given by an
iterative process which is both numerically efficient and physically
transparent. As a test of our discretized approach, we have studied two
configurations involving a dsDNA molecule confined between a pair of parallel
plates.


  We present two coarse-grained models of different levels of detail for the
description of beta-sheet tapes obtained from equilibrium self-assembly of
short rationally designed oligopeptides in solution. Here we only consider the
case of the homopolymer oligopeptides with the identical sidegroups attached,
in which the tapes have a helicoid surface with two equivalent sides. The
influence of the chirality parameter on the geometrical characteristics, namely
the diameter, inter-strand distance and pitch, of the tapes have been
investigated. The two models are found to produceequivalent results suggesting
a considerable degree of universality in conformations of the tapes.


  Many exotic species combine low probability of establishment at each
introduction with rapid population growth once introduction does succeed. To
analyze this phenomenon, we note that invaders often cluster spatially when
rare, and consequently an introduced exotic's population dynamics should depend
on locally structured interactions. Ecological theory for spatially structured
invasion relies on deterministic approximations, and determinism does not
address the observed uncertainty of the exotic-introduction process. We take a
new approach to the population dynamics of invasion and, by extension, to the
general question of invasibility in any spatial ecology. We apply the physical
theory for nucleation of spatial systems to a lattice-based model of
competition between plant species, a resident and an invader, and the analysis
reaches conclusions that differ qualitatively from the standard ecological
theories. Nucleation theory distinguishes between dynamics of single-cluster
and multi-cluster invasion. Low introduction rates and small system size
produce single-cluster dynamics, where success or failure of introduction is
inherently stochastic. Single-cluster invasion occurs only if the cluster
reaches a critical size, typically preceded by a number of failed attempts. For
this case, we identify the functional form of the probability distribution of
time elapsing until invasion succeeds. Although multi-cluster invasion for
sufficiently large systems exhibits spatial averaging and almost-deterministic
dynamics of the global densities, an analytical approximation from nucleation
theory, known as Avrami's law, describes our simulation results far better than
standard ecological approximations.


  Recent Monte Carlo simulations of a grafted semiflexible polymer in 1+1
dimensions have revealed a pronounced bimodal structure in the probability
distribution of the transverse (bending) fluctuations of the free end, when the
total contour length is of the order of the persistence length [G. Lattanzi et
al., Phys. Rev E 69, 021801 (2004)]. In this paper, we show that the emergence
of bimodality is related to a similar behavior observed when a random walker is
driven in the transverse direction by a certain type of shear flow. We adapt an
effective-medium argument, which was first introduced in the context of the
sheared random-walk problem [E. Ben-Naim et al., Phys. Rev. A 45, 7207 (1992)],
in order to obtain a simple analytic approximation of the probability
distribution of the free-end fluctuations. We show that this approximation
captures the bimodality and most of the qualitative features of the free-end
fluctuations. We also predict that relaxing the local inextensibility
constraint of the wormlike chain could lead to the disappearence of bimodality.


  Stationary solutions to a Fokker-Planck equation corresponding to a noisy
logistic equation with correlated Gaussian white noises are constructed.
Stationary distributions exist even if the corresponding deterministic system
displays an unlimited growth. Positive correlations between the noises can lead
to a minimum of the variance of the process and to the stochastic resonance if
the system is additionally driven by a periodic signal.


  Processive molecular motors which drive the traffic of organelles in cells
move in a directed way along cytoskeletal filaments. On large time scales, they
perform motor walks, i.e., peculiar random walks which arise from the repeated
unbinding from and rebinding to filaments. Unbound motors perform Brownian
motion in the surrounding fluid. In addition, the traffic of molecular motors
exhibits many cooperative phenomena. In particular, it faces similar problems
as the traffic on streets such as the occurrence of traffic jams and the
coordination of (two-way) traffic. These issues are studied here theoretically
using lattice models.


  A motile cell, when stimulated, shows a dramatic increase in the activity of
its membrane, manifested by the appearance of dynamic membrane structures such
as lamellipodia, filopodia and membrane ruffles. The external stimulus turns on
membrane bound activators, like Cdc42 and PIP2, which cause increased branching
and polymerization of the actin cytoskeleton in their vicinity leading to a
local protrusive force on the membrane. The emergence of the complex membrane
structures is a result of the coupling between the dynamics of the membrane,
the activators and the protrusive forces. We present a simple model that treats
the dynamics of a membrane under the action of actin polymerization forces that
depend on the local density of freely diffusing activators on the membrane. We
show that, depending on the spontaneous membrane curvature associated with the
activators, the resulting membrane motion can be wave-like, corresponding to
membrane ruffling and actin-waves, or unstable, indicating the tendency of
filopodia to form. Our model also quantitatively explains a variety of related
experimental observations and makes several testable predictions.


  In the present work we extend the Bak-Sneppen model for biological evolution
by introducing local interactions between species. This ``environmental''
perturbation can modify the intrinsic fitness of each element of the ecology,
leading to higher survival probability, even for the less fit. While the system
still self-organizes toward a critical state, the distribution of fitness
broadens, losing the classical step-function shape. A possible application in
economics is discussed, where firms are represented like evolving species
linked by mutual interests.


  Many systems in the realm of nanophysics from both the living and inorganic
world display slow relaxation kinetics of energy fluctuations. In this paper we
propose a general explanation for such phenomenon, based on the effects of
interactions with the solvent. Within a simple harmonic model of the system
fluctuations, we demonstrate that the inhomogeneity of coupling to the solvent
of the bulk and surface atoms suffices to generate a complex spectrum of decay
rates. We show for Myoglobin and for a metal nano-cluster that the result is a
complex, non-exponential relaxation dynamics.


  We present a model for the kinetics of spontaneous membrane domain (raft)
assembly that includes the effect of membrane recycling ubiquitous in living
cells. We show that the domains have a broad power-law distribution with an
average radius that scales with the 1/4 power of the domain lifetime when the
line tension at the domain edges is large. For biologically reasonable
recycling and diffusion rates the average domain radius is in the tens of nm
range, consistent with observations. This represents one possible link between
signaling (involving rafts) and traffic (recycling) in cells. Finally, we
present evidence that suggests that the average raft size may be the same for
all scale-free recycling schemes.


  What is the fastest way of finding a randomly hidden target? This question of
general relevance is of vital importance for foraging animals. Experimental
observations reveal that the search behaviour of foragers is generally
intermittent: active search phases randomly alternate with phases of fast
ballistic motion. In this letter, we study the efficiency of this type of two
states search strategies, by calculating analytically the mean first passage
time at the target. We model the perception mecanism involved in the active
search phase by a diffusive process. In this framework, we show that the search
strategy is optimal when the average duration of "motion phases" varies like
the power either 3/5 or 2/3 of the average duration of "search phases",
depending on the regime. This scaling accounts for experimental data over a
wide range of species, which suggests that the kinetics of search trajectories
is a determining factor optimized by foragers and that the perception activity
is adequately described by a diffusion process.


  Infectious diseases are practically represented by models with multiple
states and complex transition rules corresponding to, for example, birth,
death, infection, recovery, disease progression, and quarantine. In addition,
networks underlying infection events are often much more complex than described
by meanfield equations or regular lattices. In models with simple transition
rules such as the SIS and SIR models, heterogeneous contact rates are known to
decrease epidemic thresholds. We analyze steady states of various multi-state
disease propagation models with heterogeneous contact rates. In many models,
heterogeneity simply decreases epidemic thresholds. However, in models with
competing pathogens and mutation, coexistence of different pathogens for small
infection rates requires network-independent conditions in addition to
heterogeneity in contact rates. Furthermore, models without spontaneous
neighbor-independent state transitions, such as cyclically competing species,
do not show heterogeneity effects.


  We investigate the stretching response of a thick polymer model by means of
extensive stochastic simulations. The computational results are synthesized in
an analytic expression that characterizes how the force versus elongation curve
depends on the polymer structural parameters: its thickness and granularity
(spacing of the monomers). The expression is used to analyze experimental data
for the stretching of various different types of biopolymers: polypeptides,
polysaccharides and nucleic acids. Besides recovering elastic parameters (such
as the persistence length) that are consistent with those obtained from
standard entropic models, the approach allows to extract viable estimates for
the polymers diameter and granularity. This shows that the basic structural
polymer features have such a profound impact on the elastic behaviour that they
can be recovered with the sole input of stretching measurements.