
We present a systematic study of dynamical heterogeneity in a model for
permanent gels, upon approaching the gelation threshold. We find that the
fluctuations of the self intermediate scattering function are increasing
functions of time, reaching a plateau whose value, at large length scales,
coincides with the mean cluster size and diverges at the percolation threshold.
Another measure of dynamical heterogeneities, i.e. the fluctuations of the
selfoverlap, displays instead a peak and decays to zero at long times. The
peak, however, also scales as the mean cluster size. Arguments are given for
this difference in the long time behavior. We also find that nonGaussian
parameter reaches a plateau in the long time limit. The value of the plateau of
the nonGaussian parameter, which is connected to the fluctuations of
diffusivity of clusters, increases with the volume fraction and remains finite
at percolation threshold.

Experimental results have shown that the kinetics of bond formation in
chemical crosslinking of gelatin solutions is strongly affected not only by
gelatin and reactant concentrations but also by the solution pH. We present an
extended numerical investigation of the phase diagram and of the kinetics of
bond formation as a function of the pH, via Monte Carlo simulations of a
lattice model for gelatin chains and reactant agent in solution. We find a
reentrant phase diagram, namely gelation can be hindered either by loop
formation, at low reactant concentrations, or by saturation of active sites of
the chains via formation of single bonds with crosslinkers, at high reactant
concentrations. The ratio of the characteristic times for the formation of the
first and of the second bond between the crosslinker and an active site of a
chain is found to depend on the reactant reactivity, in good agreement with
experimental data.

We investigate the gel formation from the equilibrium sol phase in a simple
model that has the characteristics of (colloidal) gelforming systems at a
finite temperature. At low volume fraction and low temperatures, particles are
linked by longliving bonds and form an open percolating network. By means of
molecular dynamics simulations, we study the lifetime of bonds and nodes of the
gel network in order to relate these quantities to the complex relaxation
dynamics observed.

We study the formation of a colloidal gel by means of Molecular Dynamics
simulations of a model for colloidal suspensions. A slowing down with gellike
features is observed at low temperatures and low volume fractions, due to the
formation of persistent structures. We show that at low volume fraction the
dynamic susceptibility, which describes dynamic heterogeneities, exhibits a
large plateau, dominated by clusters of long living bonds. At higher volume
fraction, where the effect of the crowding of the particles starts to be
present, it crosses over towards a regime characterized by a peak. We introduce
a suitable mean cluster size of clusters of monomers connected by "persistent"
bonds which well describes the dynamic susceptibility.

We compare the slow dynamics of irreversible gels, colloidal gels, glasses
and spin glasses by analyzing the behavior of the so called nonlinear
dynamical susceptibility, a quantity usually introduced to quantitatively
characterize the dynamical heterogeneities. In glasses this quantity typically
grows with the time, reaches a maximum and then decreases at large time, due to
the transient nature of dynamical heterogeneities and to the absence of a
diverging static correlation length. We have recently shown that in
irreversible gels the dynamical susceptibility is instead an increasing
function of the time, as in the case of spin glasses, and tends asymptotically
to the mean cluster size. On the basis of molecular dynamics simulations, we
here show that in colloidal gelation where clusters are not permanent, at very
low temperature and volume fractions, i.e. when the lifetime of the bonds is
much larger than the structural relaxation time, the nonlinear susceptibility
has a behavior similar to the one of the irreversible gel, followed, at higher
volume fractions, by a crossover towards the behavior of glass forming liquids.

In chemical crosslinking of gelatin solutions, two different time scales
affect the kinetics of the gel formation in the experiments. We complement the
experimental study with Monte Carlo numerical simulations of a lattice model.
This approach shows that the two characteristic time scales are related to the
formation of single bonds crosslinkerchain and of bridges between chains. In
particular their ratio turns out to control the kinetics of the gel formation.
We discuss the effect of the concentration of chains. Finally our results
suggest that, by varying the probability of forming bridges as an independent
parameter, one can finely tune the kinetics of the gelation via the ratio of
the two characteristic times.

We use molecular dynamics computer simulations to investigate the relaxation
dynamics of a simple model for a colloidal gel at a low volume fraction. We
find that due to the presence of the open spanning network this dynamics shows
at low temperature a nontrivial dependence on the wavevector which is very
different from the one observed in dense glassforming liquids. At high wave
vectors the relaxation is due to the fast cooperative motion of the branches of
the gel network, whereas at low wave vectors the overall rearrangements of the
heterogeneous structure produce the relaxation process.

In colloidal suspensions, the competition between attractive and repulsive
interactions gives rise to a rich and complex phenomenology. Here, we study the
equilibrium phase diagram of a model system using a DLVO interaction potential
by means of molecular dynamics simulations and a thermodynamical approach. As a
result, we find tubular and lamellar phases at low volume fraction. Such
phases, extremely relevant for designing new materials, may be not easily
observed in the experiments because of the long relaxation times and the
presence of defects.

We study the structure and the dynamics in the formation of irreversible gels
by means of molecular dynamics simulation of a model system where the gelation
transition is due to the random percolation of permanent bonds between
neighboring particles. We analyze the heterogeneities of the dynamics in terms
of the fluctuations of the intermediate scattering functions: In the sol phase
close to the percolation threshold, we find that this dynamical susceptibility
increases with the time until it reaches a plateau. At the gelation threshold
this plateau scales as a function of the wave vector $k$ as $k^{\eta 2}$, with
$\eta$ being related to the decay of the percolation pair connectedness
function. At the lowest wave vector, approaching the gelation threshold it
diverges with the same exponent $\gamma$ as the mean cluster size. These
findings suggest an alternative way of measuring critical exponents in a system
undergoing chemical gelation.

By means of molecular dynamics, we study a model system for colloidal
suspensions where the interaction is based on a competition between attraction
and repulsion. At low temperatures the relaxation time $\tau$ first increases
as a power law as a function of the volume fraction $\phi$ and then, due to the
finite lifetime of the bonded structures, it deviates from this critical
behavior. We show that colloidal gelation at low temperatures and low volume
fractions is crucially related to the formation of spanning long living
cluster. Besides agreeing with experimental findings in different colloidal
systems, our results shed new light on the different role played by the
formation of long living bonds and the crowding of the particles in colloidal
structural arrest.

The increase of the viscosity, which is observed in attractive colloidal
systems by varying the temperature or the volume fraction, can be related to
the formation of structures due to particle aggregation. In particular we have
studied the non trivial dependence of the viscosity from the temperature and
the volume fraction in the copolymermicellar system L64. The comparison of the
experimental data with the results of numerical simulations in a simple model
for gelation phenomena suggests that this intriguing behavior can be explained
in terms of cluster formation and that this picture can be quite generally
extended to other attractive colloidal systems.

We here discuss the results of 3d MonteCarlo simulations of a minimal lattice
model for gelling systems. We focus on the dynamics, investigated by means of
the time autocorrelation function of the density fluctuations and the particle
mean square displacement. We start from the case of chemical gelation, i.e.
with permanent bonds, we characterize the critical dynamics as determined by
the formation of the percolating cluster, as actually observed in polymer gels.
By opportunely introducing a finite bond life time $\tau_b$ the dynamics
displays relevant changes and eventually the onset of a glassy regime. This has
been interpreted in terms of a crossover to dynamics more typical of colloidal
systems and a novel connection between classical gelation and recent results on
colloidal systems is suggested. By systematically comparing the results in the
case of permanent bonds to finite bond lifetime, the crossover and the glassy
regime can be understood in terms of effective clusters.

We investigate the slow dynamics in gelling systems by means of MonteCarlo
simulations on the cubic lattice of a minimal statistical mechanics model. By
opportunely varying some model parameter we are able to describe a crossover
from the chemical gelation behaviour to dynamics more typical of colloidal
systems. The results suggest a novel connection linking classical gelation, as
originally described by Flory, to more recent results on colloidal systems.

Within a recently introduced model based on the bondfluctuation dynamics we
study the viscoelastic behaviour of a polymer solution at the gelation
threshold. We here present the results of the numerical simulation of the model
on a cubic lattice: the percolation transition, the diffusion properties and
the time autocorrelation functions have been studied. From both the diffusion
coefficients and the relaxation times critical behaviour a critical exponent k
for the viscosity coefficient has been extracted: the two results are
comparable within the errors and are in close agreement with the Rouse model
prediction and with some experimental results. In the critical region below the
transition threshold the time autocorrelation functions show a long time tail
which is well fitted by a stretched exponential decay.