Nonequilibrium hyperuniform states in active turbulence.

We demonstrate that the complex spatiotemporal structure in active fluids can feature characteristics of hyperuniformity. Using a hydrodynamic model, we show that the transition from hyperuniformity to nonhyperuniformity and antihyperuniformity depends on the strength of active forcing and can be related to features of active turbulence without and with scaling characteristics of inertial turbulence. Combined with identified signatures of Levy walks and nonuniversal diffusion in these systems, this allows for a biological interpretation and the speculation of nonequilibrium hyperuniform states in active fluids as optimal states with respect to robustness and strategies of evasion and foraging.

Engineering and detection of light scattering directionalities in dewetted nanoresonators through dark-field scanning microscopy.

Dewetted, SiGe nanoparticles have been successfully exploited for light management in the visible and near-infrared, although their scattering properties have been so far only qualitatively studied. Here, we demonstrate that the Mie resonances sustained by a SiGe-based nanoantenna under tilted illumination, can generate radiation patterns in different directions. We introduce a novel dark-field microscopy setup that exploits the movement of the nanoantenna under the objective lens to spectrally isolate Mie resonances contribution to the total scattering cross-section during the same measurement. The knowledge of islands' aspect ratio is then benchmarked by 3D, anisotropic phase-field simulations and contributes to a correct interpretation of the experimental data.

Deformable active nematic particles and emerging edge currents in circular confinements.

We consider a microscopic field theoretical approach for interacting active nematic particles. With only steric interactions the self-propulsion strength in such systems can lead to different collective behaviour, e.g. synchronized self-spinning and collective translation. The different behaviour results from the delicate interplay between internal nematic structure, particle shape deformation and particle-particle interaction. For intermediate active strength an asymmetric particle shape emerges and leads to chirality and self-spinning crystals. For larger active strength the shape is symmetric and translational collective motion emerges. Within circular confinements, depending on the packing fraction, the self-spinning regime either stabilizes positional and orientational order or can lead to edge currents and global rotation which destroys the synchronized self-spinning crystalline structure.

Topological fine structure of smectic grain boundaries and tetratic disclination lines within three-dimensional smectic liquid crystals.

Observing and characterizing the complex ordering phenomena of liquid crystals subjected to external constraints constitutes an ongoing challenge for chemists and physicists alike. To elucidate the delicate balance appearing when the intrinsic positional order of smectic liquid crystals comes into play, we perform Monte-Carlo simulations of rod-like particles in a range of cavities with a cylindrical symmetry. Based on recent insights into the topology of smectic orientational grain boundaries in two dimensions, we analyze the emerging three-dimensional defect structures from the perspective of tetratic symmetry. Using an appropriate three-dimensional tetratic order parameter constructed from the Steinhardt order parameters, we show that those grain boundaries can be interpreted as a pair of tetratic disclination lines that are located on the edges of the nematic domain boundary. Thereby, we shed light on the fine structure of grain boundaries in three-dimensional confined smectics.

Mesoscale Defect Motion in Binary Systems: Effects of Compositional Strain and Cottrell Atmospheres.

The velocity of dislocations is derived analytically to incorporate and predict the intriguing effects induced by the preferential solute segregation and Cottrell atmospheres in both two-dimensional and three-dimensional binary systems of various crystalline symmetries. The corresponding mesoscopic description of defect dynamics is constructed through the amplitude formulation of the phase-field crystal model, which has been shown to accurately capture elasticity and plasticity in a wide variety of systems. Modifications of the Peach-Koehler force as a result of solute concentration variations and compositional stresses are presented, leading to interesting new predictions of defect motion due to effects of Cottrell atmospheres. These include the deflection of dislocation glide paths, the variation of climb speed and direction, and the change or prevention of defect annihilation, all of which play an important role in determining the fundamental behaviors of complex defect network and dynamics. The analytic results are verified by numerical simulations.

Hyperuniform Monocrystalline Structures by Spinodal Solid-State Dewetting.

Materials featuring anomalous suppression of density fluctuations over large length scales are emerging systems known as disordered hyperuniform. The underlying hidden order renders them appealing for several applications, such as light management and topologically protected electronic states. These applications require scalable fabrication, which is hard to achieve with available top-down approaches. Theoretically, it is known that spinodal decomposition can lead to disordered hyperuniform architectures. Spontaneous formation of stable patterns could thus be a viable path for the bottom-up fabrication of these materials. Here, we show that monocrystalline semiconductor-based structures, in particular Si_{1-x}Ge_{x} layers deposited on silicon-on-insulator substrates, can undergo spinodal solid-state dewetting featuring correlated disorder with an effective hyperuniform character. Nano- to micrometric sized structures targeting specific morphologies and hyperuniform character can be obtained, proving the generality of the approach and paving the way for technological applications of disordered hyperuniform metamaterials. Phase-field simulations explain the underlying nonlinear dynamics and the physical origin of the emerging patterns.

Properties of surface Landau-de Gennes Q-tensor models.

Uniaxial nematic liquid crystals whose molecular orientation is subjected to tangential anchoring on a curved surface offer a non trivial interplay between the geometry and the topology of the surface and the orientational degree of freedom. We consider a general thin film limit of a Landau-de Gennes Q-tensor model which retains the characteristics of the 3D model. From this, previously proposed surface models follow as special cases. We compare fundamental properties, such as the alignment of the orientational degrees of freedom with principle curvature lines, order parameter symmetry and phase transition type for these models, and suggest experiments to identify suitable model assumptions.

Smectic monolayer confined on a sphere: topology at the particle scale.

The impact of topology on the structure of a smectic monolayer confined to a sphere is explored by particle-resolved computer simulations of hard rods. The orientations of the particles are tangential to the sphere and either free or restricted to a prescribed director field with a latitude or longitude orderings. Depending on the imprinted topology, a wealth of different states are found including equatorial smectic with isotropic poles, equatorial smectic with empty poles, a broken egg-shell like modulated smectic, a capped nematic with equatorial bald patches, equatorial nematic with empty poles, and a situation with 4 or 8 half-strength topological defects. Potentially these states could be verified in experiments with Pickering emulsions of droplets with colloidal rods. The unique nature of dipolar structures consisting of positive and negative half-strength disclinations is revealed. These structures, classified by their density and interaction with other defects in the system, relieve the strain of the poles by separating closely positioned half-strength defects. The proximity of these structures to the half-strength defects might enhance the structural diffusion of the defects across the system.

Stress Induced Branching of Growing Crystals on Curved Surfaces.

If two-dimensional crystals grow on a curved surface, the Gaussian curvature of the surface induces elastic stress and affects the growth pathway. The elastic stress can be alleviated by incorporating defects or, if this is energetically unfavorable, via an elastic instability which leads to anisotropic growth with branched ribbonlike structures. This instability provides a generic route to grow defect-free crystals on curved surfaces. Depending on the elastic properties of the crystal and the geometric properties of the surface, different growth morphologies with two-, four-, and sixfold symmetry develop. Using a phase field crystal type modeling approach, we provide a microscopic understanding of the morphology selection.

A Mechanistic Collective Cell Model for Epithelial Colony Growth and Contact Inhibition.

We present a mechanistic hybrid continuum-discrete model to simulate the dynamics of epithelial cell colonies. Collective cell dynamics are modeled using continuum equations that capture plastic, viscoelastic, and elastic deformations in the clusters while providing single-cell resolution. The continuum equations can be viewed as a coarse-grained version of previously developed discrete models that treat epithelial clusters as a two-dimensional network of vertices or stochastic interacting particles and follow the framework of dynamic density functional theory appropriately modified to account for cell size and shape variability. The discrete component of the model implements cell division and thus influences cell size and shape that couple to the continuum component. The model is validated against recent in vitro studies of epithelial cell colonies using Madin-Darby canine kidney cells. In good agreement with experiments, we find that mechanical interactions and constraints on the local expansion of cell size cause inhibition of cell motion and reductive cell division. This leads to successively smaller cells and a transition from exponential to quadratic growth of the colony that is associated with a constant-thickness rim of growing cells at the cluster edge, as well as the emergence of short-range ordering and solid-like behavior. A detailed analysis of the model reveals a scale invariance of the growth and provides insight into the generation of stresses and their influence on the dynamics of the colonies. Compared to previous models, our approach has several advantages: it is independent of dimension, it can be parameterized using classical elastic properties (Poisson's ratio and Young's modulus), and it can easily be extended to incorporate multiple cell types and general substrate geometries.

Multi-scale simulation of plant stem reinforcement by brachysclereids: A case study in apple fruit peduncles.

Sclereid formation in addition to or in gaps of fragmented fibre rings is common in dicotyledonous plant stems. Whether this sclereid formation is force-triggered remains open so far. In fruit peduncles of several Malus species as modified plant stems, for example, the persistent fibre ring is displaced to the centre by formation of cortex parenchyma during growth. Parenchyma cells subsequently differentiate into an additional layer of brachysclereids, previously interpreted as an adaptation to continuously rising fruit loads. The present study pursues a multi-scale numerical modelling approach, to verify the important effect for different cellular architectures in both sclerenchyma categories on the stiffness of these tissues and the entire peduncle. First, different material properties are simulated analogue to plant tissues on the basis of three cell types. A regular three-dimensional and a random Voronoi microstructure combined with various mechanical cell wall parameters are applied. Using homogenisation simulations based on HILL's principle, numerical calculations predict a lower effective homogenised tissue stiffness of isodiametric brachysclereids compared to those of fibres, confirming experimentally obtained data from Malus fruit peduncles. Furthermore, a curved peduncle model with a complex arrangement of different material layers is generated. Diverse material sets are tested under three representative loadings, using an adaptive diffuse domain approach (AMDiS). The model explains the function of sclereids as considerable contributors to the stiffness against bending and tensile deformations, as well as torsion, especially in consequence of superimposed load conditions in the case of a curved plant stem.

A Navier-Stokes phase-field crystal model for colloidal suspensions.

We develop a fully continuous model for colloidal suspensions with hydrodynamic interactions. The Navier-Stokes Phase-Field Crystal model combines ideas of dynamic density functional theory with particulate flow approaches and is derived in detail and related to other dynamic density functional theory approaches with hydrodynamic interactions. The derived system is numerically solved using adaptive finite elements and is used to analyze colloidal crystallization in flowing environments demonstrating a strong coupling in both directions between the crystal shape and the flow field. We further validate the model against other computational approaches for particulate flow systems for various colloidal sedimentation problems.

Crystalline order and topological charges on capillary bridges.

We numerically investigate crystalline order on negative Gaussian curvature capillary bridges. In agreement with the experimental results in [W. Irvine et al., Nature, Pleats in crystals on curved surfaces, 2010, 468, 947] we observe for decreasing integrated Gaussian curvature, a sequence of transitions, from no defects to isolated dislocations, pleats, scars and isolated sevenfold disclinations. We especially focus on the dependency of topological charge on the integrated Gaussian curvature, for which we observe, again in agreement with the experimental results, no net disclination for an integrated curvature down to -10, and an approximately linear behavior from there on until the disclinations match the integrated curvature of -12. In contrast to previous studies in which ground states for each geometry are searched for, we here show that the experimental results, which are likely to be in a metastable state, can be best resembled by mimicking the experimental settings and continuously changing the geometry. The obtained configurations are only low energy local minima. The results are computed using a phase field crystal approach on catenoid-like surfaces and are highly sensitive to the initialization.

Signaling networks and cell motility: a computational approach using a phase field description.

The processes of protrusion and retraction during cell movement are driven by the turnover and reorganization of the actin cytoskeleton. Within a reaction-diffusion model which combines processes along the cell membrane with processes within the cytoplasm a Turing type instability is used to form the necessary polarity to distinguish between cell front and rear and to initiate the formation of different organizational arrays within the cytoplasm leading to protrusion and retraction. A simplified biochemical network model for the activation of GTPase which accounts for the different dimensionality of the cell membrane and the cytoplasm is used for this purpose and combined with a classical Helfrich type model to account for bending and stiffness effects of the cell membrane. In addition streaming within the cytoplasm and the extracellular matrix is taken into account. Combining these phenomena allows to simulate the dynamics of cells and to reproduce the primary phenomenology of cell motility. The coupled model is formulated within a phase field approach and solved using adaptive finite elements.

Active smectics on a sphere.

The dynamics of active smectic liquid crystals confined on a spherical surface is explored through an active phase field crystal model. Starting from an initially randomly perturbed isotropic phase, several types of topological defects are spontaneously formed, and then annihilate during a coarsening process until a steady state is achieved. The coarsening process is highly complex involving several scaling laws of defect densities as a function of time where different dynamical exponents can be identified. In general the exponent for the final stage towards the steady state is significantly larger than that in the passive and in the planar case, i.e. the coarsening is getting accelerated both by activity and by the topological and geometrical properties of the sphere. A defect type characteristic for this active system is a rotating spiral of evolving smectic layering lines. On a sphere this defect type also determines the steady state. Our results can in principle be confirmed by dense systems of synthetic or biological active particles.

Wetting resistance at its topographical limit: the benefit of mushroom and serif T structures.

Springtails (Collembola) are wingless arthropods adapted to cutaneous respiration in temporarily rain-flooded habitats. They immediately form a plastron, protecting them against suffocation upon immersion into water and even low-surface-tension liquids such as alkanes. Recent experimental studies revealed a high-pressure resistance of such plastrons against collapse. In this work, skin sections of Orthonychiurus stachianus are studied by transmission electron microscopy. The micrographs reveal cavity side-wall profiles with characteristic overhangs. These were fitted by polynomials to allow access for analytical and numerical calculations of the breakthrough pressure, that is, the barrier against plastron collapse. Furthermore, model profiles with well-defined geometries were used to set the obtained results into context and to develop a general design principle for the most robust surface structures. Our results indicate the decisive role of the sectional profile of overhanging structures to form a robust heterogeneous wetting state for low-surface-tension liquids that enables the omniphobicity. Furthermore, the design principles of mushroom and serif T structures pave the way for omniphobic surfaces with a high-pressure resistance irrespective of solid surface chemistry.

Robust statistical properties of T1 transitions in a multi-phase field model of cell monolayers.

Large-scale tissue deformation which is fundamental to tissue development hinges on local cellular rearrangements, such as T1 transitions. In the realm of the multi-phase field model, we analyse the statistical and dynamical properties of T1 transitions in a confluent monolayer. We identify an energy profile that is robust to changes in several model parameters. It is characterized by an asymmetric profile with a fast increase in energy before the T1 transition and a sudden drop after the T1 transition, followed by a slow relaxation. The latter being a signature of the fluidity of the cell monolayer. We show that T1 transitions are sources of localised large deformation of the cells undergoing the neighbour exchange, and they induce other T1 transitions in the nearby cells leading to a chaining of events that propagate local cell deformation to large scale tissue flows.

Axis convergence in C. elegans embryos.

Embryos develop in a surrounding that guides key aspects of their development. For example, the anteroposterior (AP) body axis is always aligned with the geometric long axis of the surrounding eggshell in fruit flies and worms. The mechanisms that ensure convergence of the AP axis with the long axis of the eggshell remain unresolved. We investigate axis convergence in early C. elegans development, where the nascent AP axis, when misaligned, actively re-aligns to converge with the long axis of the egg. We identify two physical mechanisms that underlie axis convergence. First, bulk cytoplasmic flows, driven by actomyosin cortical flows, can directly reposition the AP axis. Second, active forces generated within the pseudocleavage furrow, a transient actomyosin structure similar to a contractile ring, can drive a mechanical re-orientation such that it becomes positioned perpendicular to the long axis of the egg. This in turn ensures AP axis convergence. Numerical simulations, together with experiments that either abolish the pseudocleavage furrow or change the shape of the egg, demonstrate that the pseudocleavage-furrow-dependent mechanism is a major driver of axis convergence. We conclude that active force generation within the actomyosin cortical layer drives axis convergence in the early nematode.

The influence of membrane bound proteins on phase separation and coarsening in cell membranes.

A theoretical explanation of the existence of lipid rafts in cell membranes remains a topic of lively debate. Large, micrometer sized rafts are readily observed in artificial membranes and can be explained using thermodynamic models for phase separation and coarsening. In live cells such domains are not observed and various models are proposed to describe why the systems do not coarsen. We review these attempts critically and show within a phase field approach that membrane bound proteins have the potential to explain the different behaviour observed in vitro and in vivo. Large scale simulations are performed to compute scaling laws and size distribution functions under the influence of membrane bound proteins and to observe a significant slow down of the domain coarsening at longer times and a breakdown of the self-similarity of the size-distribution function.

Faceting of Si and Ge crystals grown on deeply patterned Si substrates in the kinetic regime: phase-field modelling and experiments.

The development of three-dimensional architectures in semiconductor technology is paving the way to new device concepts for various applications, from quantum computing to single photon avalanche detectors. In most cases, such structures are achievable only under far-from-equilibrium growth conditions. Controlling the shape and morphology of the growing structures, to meet the strict requirements for an application, is far more complex than in close-to-equilibrium cases. The development of predictive simulation tools can be essential to guide the experiments. A versatile phase-field model for kinetic crystal growth is presented and applied to the prototypical case of Ge/Si vertical microcrystals grown on deeply patterned Si substrates. These structures, under development for innovative optoelectronic applications, are characterized by a complex three-dimensional set of facets essentially driven by facet competition. First, the parameters describing the kinetics on the surface of Si and Ge are fitted on a small set of experimental results. To this goal, Si vertical microcrystals have been grown, while for Ge the fitting parameters have been obtained from data from the literature. Once calibrated, the predictive capabilities of the model are demonstrated and exploited for investigating new pattern geometries and crystal morphologies, offering a guideline for the design of new 3D heterostructures. The reported methodology is intended to be a general approach for investigating faceted growth under far-from-equilibrium conditions.