
Spatially nonuniform strain is important for engineering the pseudomagnetic
field and band structure of graphene. Despite the wide interest in strain
engineering, there is still a lack of control on devicecompatible strain
patterns due to the limited understanding of the structurestrain relationship.
Here, we study the effect of substrate corrugation and curvature on the strain
profiles of graphene via combined experimental and theoretical studies of a
model system: graphene on closely packed SiO2 nanospheres with different
diameters (20200 nm). Experimentally, via quantitative Raman analysis, we
observe partial adhesion and wrinkle features and find that smaller nanospheres
induce larger tensile strain in graphene, theoretically, molecular dynamics
simulations confirm the same microscopic structure and size dependence of
strain and reveal that a larger strain is caused by a stronger, inhomogeneous
interaction force between smaller nanospheres and graphene. This
molecularlevel understanding of the strain mechanism is important for strain
engineering of graphene and other twodimensional materials.

We have created and studied artificial magnetic quasicrystals based on
Penrose tiling patterns of interacting nanomagnets that lack the translational
symmetry of spatially periodic artificial spin ices. Vertexlevel degeneracy
and frustration induced by the network topology of the Penrose pattern leads to
a low energy configuration that we propose as a ground state. Topologically
induced emergent frustration means that this ground state cannot be constructed
from vertices in their ground states. It has two parts, a quasionedimensional
rigid "skeleton" that spans the entire pattern and is capable of longrange
order, and "flippable" clusters of macrospins within it. These lead to
macroscopic degeneracy for the array as a whole. Magnetic force microscopy
imaging of Penrose tiling arrays revealed superdomains that are larger for more
strongly coupled arrays. The superdomain size is larger after
ACdemagnetisation and especially after annealing the array above its blocking
temperature.

When strained beyond the linear regime, soft colloidal glasses yield to
steadystate plastic flow in a way that is similar to the deformation of
conventional amorphous solids. Due to the much larger size of the colloidal
particles with respect to the atoms comprising an amorphous solid, colloidal
glasses allow to obtain microscopic insight into the nature of the yielding
transition, as we illustrate here combining experiments, atomistic simulations,
and mesoscopic modeling. Our results unanimously show growing clusters of
nonaffine deformation percolating at yielding. In agreement with percolation
theory, the spanning cluster is fractal with a fractal dimension d_f~2, and the
correlation length diverges upon approaching the critical yield strain. These
results indicate that percolation of highly nonaffine particles is the
hallmark of the yielding transition in disordered glassy systems.

Plastic yield of amorphous solids occurs by power law distributed slip
avalanches whose universality is still debated. Determination of the power law
exponents from experiments and molecular dynamics simulations is hampered by
limited statistical sampling. On the other hand, while existing elastoplastic
depinning models give precise exponent values, these models to date have been
limited to a scalar approximation of plasticity which is difficult to reconcile
with the statistical isotropy of amorphous materials. Here we introduce for the
first time a fully tensorial mesoscale model for the elastoplasticity of
disordered media that can not only reproduce a wide variety of shear band
patterns observed experimentally for different deformation modes, but also
captures the avalanche dynamics of plastic flow in disordered materials. Slip
avalanches are characterized by universal distributions which are
quantitatively different from mean field predictions, both regarding the
exponents and regarding the form of the scaling functions, and which are
independent of system dimensionality (2D vs 3D), boundary and loading
conditions, and unior biaxiality of the stress state. We also measure average
avalanche shapes, which are equally universal and inconsistent with mean field
predictions. Our results provide strong evidence that the universality class of
plastic yield in amorphous materials is distinct from that of mean field
depinning.

Autocatalytic fibril nucleation has recently been proposed to be a
determining factor for the spread of neurodegenerative diseases, but the same
process could also be exploited to amplify minute quantities of protein
aggregates in a diagnostic context. Recent advances in microfluidic technology
allow analysis of protein aggregation in micronscale samples potentially
enabling such diagnostic approaches, but the theoretical foundations for the
analysis and interpretation of such data are so far lacking. Here we study
computationally the onset of protein aggregation in small volumes and show that
the process is ruled by intrinsic fluctuations whose volume dependent
distribution we also estimate theoretically. Based on these results, we develop
a strategy to quantify in silico the statistical errors associated with the
detection of aggregate containing samples. Our work opens a new perspective on
the forecasting of protein aggregation in asymptomatic subjects.

Faithful segregation of genetic material during cell division requires
alignment of chromosomes between two spindle poles and attachment of their
kinetochores to each of the poles. Failure of these complex dynamical processes
leads to chromosomal instability (CIN), a characteristic feature of several
diseases including cancer. While a multitude of biological factors regulating
chromosome congression and biorientation have been identified, it is still
unclear how they are integrated so that coherent chromosome motion emerges from
a large collection of random and deterministic processes. Here we address this
issue by a three dimensional computational model of motordriven chromosome
congression and biorientation during mitosis. Our model reveals that
successful cell division requires control of the total number of microtubules:
if this number is too small biorientation fails, while if it is too large not
all the chromosomes are able to congress. The optimal number of microtubules
predicted by our model compares well with early observations in mammalian cell
spindles. Our results shed new light on the origin of several pathological
conditions related to chromosomal instability.

Motor proteins display widely different stepping patterns as they move on
microtubule tracks, from the deterministic linear or helical motion performed
by the protein kinesin to the uncoordinated random steps made by dynein. How
these different strategies produce an efficient navigation system needed to
ensure correct cellular functioning is still unclear. Here, we show by
numerical simulations that deterministic and random motor steps yield different
outcomes when random obstacles decorate the microtubule tracks: kinesin moves
faster on clean tracks but its motion is strongly hindered on decorated tracks,
while dynein is slower on clean tracks but more efficient in avoiding
obstacles. Further simulations indicate that dynein's advantage on decorated
tracks is due to its ability to step backwards. Our results explain how
different navigation strategies are employed by the cell to optimize motor
driven cargo transport.

Crystalline plasticity is strongly interlinked with dislocation mechanics and
nowadays is relatively well understood. Concepts and physical models of plastic
deformation in amorphous materials on the other hand  where the concept of
linear lattice defects is not applicable  still are lagging behind. We
introduce an eigenstrainbased finite element lattice model for simulations of
shear band formation and strain avalanches. Our model allows us to study the
influence of surfaces and finite size effects on the statistics of avalanches.
We find that even with relatively complex loading conditions and open boundary
conditions, critical exponents describing avalanche statistics are unchanged,
which validates the use of simpler scalar latticebased models to study these
phenomena.

Delamination of coatings and thin films from substrates generates a
fascinating variety of patterns, from circular blisters to wrinkles and
labyrinth domains, in a way that is not completely understood. We report on
largescale numerical simulations of the universal problem of avoidance and
coalescence of delamination wrinkles, focusing on a case study of graphene
sheets on patterned substrates. By nucleating and growing wrinkles in a
controlled way, we are able to characterize how their interactions, mediated by
longrange stress fields, determine their formation and morphology. We also
study how the interplay between geometry and stresses drives a universal
transition from conformation to delamination when sheets are deposited on
particledecorated substrates. Our results are directly applicable to strain
engineering of graphene and also uncover universal phenomena observed at all
scales, as for example in geomembrane deposition.

We perform large scale simulations of a two dimensional lattice model for
amorphous plasticity with random local yield stresses and longrange
quadrupolar elastic interactions. We show that as the external stress increases
towards the yielding phase transition, the scaling behavior of the avalanches
crosses over from meanfield theory to a different universality class. This
behavior is associated with strain localization, which significantly depends on
the shortrange properties of the interaction kernel.

Typically, the plastic yield stress of a sample is determined from a
stressstrain curve by defining a yield strain and reading off the stress
required to attain it. However, it is not a priori clear that yield strengths
of microscale samples measured this way should display the correct finite size
scaling. Here we study plastic yield as a depinning transition of a 1+1
dimensional interface, and consider how finite size effects depend on the
choice of yield strain, as well as the presence of hardening and the strength
of elastic coupling. Our results indicate that in sufficiently large systems,
the choice of yield strain is unimportant, but in smaller systems one must take
care to avoid spurious effects.

Quenched disorder affects how nonequilibrium systems respond to driving. In
the context of artificial spin ice, an athermal system comprised of
geometrically frustrated classical Ising spins with a twofold degenerate
ground state, we give experimental and numerical evidence of how such disorder
washes out edge effects, and provide an estimate of disorder strength in the
experimental system. We prove analytically that a sequence of applied fields
with fixed amplitude is unable to drive the system to its ground state from a
saturated state. These results should be relevant for other systems where
disorder does not change the nature of the ground state.

We have performed a systematic study of the effects of field strength and
quenched disorder on the driven dynamics of square artificial spin ice. We
construct a network representation of the configurational phase space, where
nodes represent the microscopic configurations and a directed link between node
i and node j means that the field may induce a transition between the
corresponding configurations. In this way, we are able to quantitatively
describe how the field and the disorder affect the connectedness of states and
the reversibility of dynamics. In particular, we have shown that for optimal
field strengths, a substantial fraction of all states can be accessed using
external driving fields, and this fraction is increased by disorder. We discuss
how this relates to control and potential information storage applications for
artificial spin ices.

The thermallydriven formation and evolution of vertex domains is studied for
square artificial spin ice. A self consistent mean field theory is used to show
how domains of ground state ordering form spontaneously, and how these evolve
in the presence of disorder. The role of fluctuations is studied, using Monte
Carlo simulations and analytical modelling. Domain wall dynamics are shown to
be driven by a biasing of random fluctuations towards processes that shrink
closed domains, and fluctuations within domains are shown to generate isolated
small excitations, which may stabilise as the effective temperature is lowered.
Domain dynamics and fluctuations are determined by interaction strengths, which
are controlled by interelement spacing. The role of interaction strength is
studied via experiments and Monte Carlo simulations. Our mean field model is
applicable to ferroelectric `spin' ice, and we show that features similar to
that of magnetic spin ice can be expected, but with different characteristic
temperatures and rates.

We report a novel approach to the question of whether and how the ground
state can be achieved in square artificial spin ices where frustration is
incomplete. We identify two types of disorder: quenched disorder in the island
response to fields and disorder in the sequence of driving fields. Numerical
simulations show that quenched disorder can lead to final states with lower
energy, and disorder in the driving fields always lowers the final energy
attained by the system. We use a network picture to understand these two
effects: disorder in island responses creates new dynamical pathways, and
disorder in driving fields allows more pathways to be followed.

The fieldinduced dynamics of artificial spin ice are determined in part by
interactions between magnetic islands, and the switching characteristics of
each island. Disorder in either of these affects the response to applied
fields. Numerical simulations are used to show that disorder effects are
determined primarily by the strength of disorder relative to interisland
interactions, rather than by the type of disorder. Weak and strong disorder
regimes exist and can be defined in a quantitative way.

Local magnetic ordering in artificial spin ices is discussed from the point
of view of how geometrical frustration controls dynamics and the approach to
steady state. We discuss the possibility of using a particle picture based on
vertex configurations to interpret time evolution of magnetic configurations.
Analysis of possible vertex processes allows us to anticipate different
behaviors for open and closed edges and the existence of different field
regimes. Numerical simulations confirm these results and also demonstrate the
importance of correlations and long range interactions in understanding
particle population evolution. We also show that a mean field model of vertex
dynamics gives important insights into finite size effects.