Source author record

Rastko Sknepnek

Rastko Sknepnek appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2026arXiv

Machine learning short-ranged many-body interactions in colloidal systems using descriptors based on Voronoi cells

Machine learning (ML) strategies are opening the door to faster computer simulations, allowing us to simulate more realistic colloidal systems. Since the interactions in colloidal systems are often highly many-body, stemming from e.g. depletion and steric interactions, one of the challenges for these algorithms is capturing the many-body nature of these interactions. In this paper, we introduce a new ML-based strategy for fitting many-body interactions in colloidal systems where the many-body interaction is highly local. To this end, we develop Voronoi-based descriptors for capturing the local environment and fit the effective potential using a simple neural network. To test this algorithm, we consider a simple two-dimensional model for a colloid-polymer mixture, where the colloid-colloid interactions and colloid-polymer interactions are hard-disk like, while the polymers themselves interact as ideal gas particles. We find that a Voronoi-based description is sufficient to accurately capture the many-body nature of this system. Moreover, we find that the Pearson correlation function alone is insufficient to determine the predictive power of the network emphasizing the importance of additional metrics when assessing the quality of ML-based potentials.

preprint2022arXiv

Generating active T1 transitions through mechanochemical feedback

Convergence-extension in embryos is controlled by chemical and mechanical signalling. A key cellular process is the exchange of neighbours via T1 transitions. We propose and analyse a model with positive feedback between recruitment of myosin motors and mechanical tension in cell junctions. The model produces active T1 events, which act to elongate the tissue perpendicular to external pulling. Using an idealized tissue patch comprising several active cells embedded in a matrix of passive hexagonal cells we identified an optimal range of pulling forces to trigger an active T1 event. We show that pulling also generates tension chains in a realistic patch made entirely of active cells of random shapes, and leads to convergence-extension over a range of parameters. Our findings show that active intercalations can generate stress that activates T1 events in neighbouring cells resulting in tension dependent tissue re-organisation, in qualitative agreement with experiments on gastrulation in chick embryos.

preprint2020arXiv

Dense active matter model of motion patterns in confluent cell monolayers

Epithelial cell monolayers show remarkable displacement and velocity correlations over distances of ten or more cell sizes that are reminiscent of supercooled liquids and active nematics. We show that many observed features can be described within the framework of dense active matter, and argue that persistent uncoordinated cell motility coupled to the collective elastic modes of the cell sheet is sufficient to produce swirl-like correlations. We obtain this result using both continuum active linear elasticity and a normal modes formalism, and validate analytical predictions with numerical simulations of two agent-based cell models, soft elastic particles and the self-propelled Voronoi model together with in-vitro experiments of confluent corneal epithelial cell sheets. Simulations and normal mode analysis perfectly match when tissue-level reorganisation occurs on times longer than the persistence time of cell motility. Our analytical model quantitatively matches measured velocity correlation functions over more than a decade with a single fitting parameter.

preprint2019arXiv

Wrinkle patterns in active viscoelastic thin sheets

We show that a viscoelastic thin sheet driven out of equilibrium by active structural remodelling develops a rich variety of shapes as a result of a competition between viscous relaxation and activity. In the regime where active processes are faster than viscoelastic relaxation, wrinkles that are formed due to remodelling are unable to relax to a configuration that minimises the elastic energy and the sheet is inherently out of equilibrium. We argue that this non-equilibrium regime is of particular interest in biology as it allows the system to access morphologies that are unavailable if restricted to the adiabatic evolution between configurations that minimise the elastic energy alone. Here, we introduce activity using the formalism of evolving target metric and showcase the diversity of wrinkling morphologies arising from out of equilibrium dynamics.

preprint2016arXiv

Active Vertex Model for Cell-Resolution Description of Epithelial Tissue Mechanics

We introduce an Active Vertex Model (AVM) for cell-resolution studies of the mechanics of confluent epithelial tissues consisting of tens of thousands of cells, with a level of detail inaccessible to similar methods. The AVM combines the Vertex Model for confluent epithelial tissues with active matter dynamics. This introduces a natural description of the cell motion and accounts for motion patterns observed on multiple scales. Furthermore, cell contacts are generated dynamically from positions of cell centres. This not only enables efficient numerical implementation, but provides a natural description of the T1 transition events responsible for local tissue rearrangements. The AVM also includes cell alignment, cell-specific mechanical properties, cell growth, division and apoptosis. In addition, the AVM introduces a flexible, dynamically changing boundary of the epithelial sheet allowing for studies of phenomena such as the fingering instability or wound healing. We illustrate these capabilities with a number of case studies.

preprint2014arXiv

Active swarms on a sphere

Here we show that coupling to curvature has profound effects on collective motion in active systems, leading to patterns not observed in flat space. Biological examples of such active motion in curved environments are numerous: curvature and tissue folding are crucial during gastrulation, epithelial and endothelial cells move on constantly growing, curved crypts and vili in the gut, and the mammalian corneal epithelium grows in a steady-state vortex pattern. On the physics side, droplets coated with actively driven microtubule bundles show active nematic patterns. We study a model of self-propelled particles with polar alignment on a sphere. Hallmarks of these motion patterns are a polar vortex and a circulating band arising due to the incompatibility between spherical topology and uniform motion - a consequence of the hairy ball theorem. We present analytical results showing that frustration due to curvature leads to stable elastic distortions storing energy in the band.

preprint2014arXiv

Defect dynamics in active nematics

Topological defects are distinctive signatures of liquid crystals. They profoundly affect the viscoelastic behavior of the fluid by constraining the orientational structure in a way that inevitably requires global changes not achievable with any set of local deformations. In active nematic liquid crystals topological defects not only dictate the global structure of the director, but also act as local sources of motion, behaving as self-propelled particles. In this article we present a detailed analytical and numerical study of the mechanics of topological defects in active nematic liquid crystals.

preprint2014arXiv

Effects of scars on crystalline shell stability under external pressure

We study how the stability of spherical crystalline shells under external pressure is influenced by the defect structure. In particular, we compare stability for shells with a minimal set of topologically-required defects to shells with extended defect arrays (grain boundary "scars" with non-vanishing net disclination charge). We perform Monte Carlo simulations to compare how shells with and without scars deform quasi-statically under external hydrostatic pressure. We find that the critical pressure at which shells collapse is lowered for scarred configurations that break icosahedral symmetry and raised for scars that preserve icosahedral symmetry. The particular shapes which arise from breaking of an initial icosahedrally-symmetric shell depend on the Föppl-von Kármán number.

preprint2013arXiv

Optimal shapes and stresses of adherent cells on patterned substrates

We investigate a continuum mechanical model for an adherent cell on two dimensional adhesive micropatterned substrates. The cell is modeled as an isotropic and homogeneous elastic material subject to uniform internal contractile stresses. The build-up of tension from cortical actin bundles at the cell periphery is incorporated by introducing an energy cost for bending of the cell boundary, resulting to a resistance to changes in local curvature. Integrin-based adhesions are modeled as harmonic springs, that pin the cell to adhesive patches of a predefined geometry. Using Monte Carlo simulations and analytical techniques we investigate the competing effects of bulk contractility and cortical bending rigidity in regulating cell shapes on non-adherent regions. We show that the crossover from convex to concave cell edges is controlled by the interplay between contractile stresses and boundary bending rigidity. In particular, the cell boundary becomes concave beyond a critical value of the contractile stress that is proportional to the cortical bending rigidity. Furthermore, the intracellular stresses are found largely concentrated at the concave edge of the cell. The model can be used to generate a cell-shape phase diagram for each specific adhesion geometry.

preprint2013arXiv

Pathways to faceting of vesicles

The interplay between geometry, topology and order can lead to geometric frustration that profoundly affects the shape and structure of a curved surface. In this commentary we show how frustration in this context can result in the faceting of elastic vesicles. We show that, under the right conditions, an assortment of regular and irregular polyhedral structures may be the low energy states of elastic membranes with spherical topology. In particular, we show how topological defects, necessarily present in any crystalline lattice confined to spherical topology, naturally lead to the formation of icosahedra in a homogeneous elastic vesicle. Furthermore, we show that introducing heterogeneities in the elastic properties, or allowing for non-linear bending response of a homogeneous system, opens non-trivial pathways to the formation of faceted, yet non-icosahedral, structures.

preprint2013arXiv

Shapes of pored membranes

We study the shapes of pored membranes within the framework of the Helfrich theory under the constraints of fixed area and pore size. We show that the mean curvature term leads to a budding- like structure, while the Gaussian curvature term tends to flatten the membrane near the pore; this is corroborated by simulation. We propose a scheme to deduce the ratio of the Gaussian rigidity to the bending rigidity simply by observing the shape of the pored membrane. This ratio is usually difficult to measure experimentally. In addition, we briefly discuss the stability of a pore by relaxing the constraint of a fixed pore size and adding the line tension. Finally, the flattening effect due to the Gaussian curvature as found in studying pored membranes is extended to two-component membranes. We find that sufficiently high contrast between the components' Gaussian rigidities leads to budding which is distinct from that due to the line tension.

preprint2012arXiv

A curvature-driven effective attraction in multicomponent membranes

We study closed liquid membranes that segregate into three phases due to differences in the chemical and physical properties of its components. The shape and in-plane membrane arrangement of the phases are coupled through phase-specific bending energies and line tensions. We use simulated annealing Monte Carlo simulations to find low-energy structures, allowing both phase arrangement and membrane shape to relax. The three-phase system is the simplest one in which there are multiple interface pairs, allowing us to analyze interfacial preferences and pairwise distinct line tensions. We observe the system's preference for interface pairs that maximize differences in spontaneous curvature. From a pattern selection perspective, this acts as an effective attraction between phases of most disparate spontaneous curvature. We show that this effective attraction is robust enough to persist even when the interface between these phases is the most penalized by line tension. This effect is driven by geometry and not by any explicit component-component interaction.

preprint2012arXiv

Nonlinear elastic model for faceting of vesicles with soft grain boundaries

We use an elastic model to explore faceting of solid-wall vesicles with elastic heterogeneities. We show that faceting occurs in regions where the vesicle wall is softer, such as areas of reduced wall thicknesses or concentrated in crystalline defects. The elastic heterogeneities are modeled as a second component with reduced elastic parameters. Using simulated annealing Monte Carlo simulations we obtain the vesicle shape by optimizing the distributions of facets and boundaries. Our model allows us to reduce the effects of the residual stress generated by crystalline defects, and reveals a robust faceting mechanism into polyhedra other than the icosahedron.

preprint2012arXiv

Planar sheets meet negative curvature liquid interfaces

If an inextensible thin sheet is adhered to a substrate with a negative Gaussian curvature it will experience stress due to geometric frustration. We analyze the consequences of such geometric frustration using analytic arguments and numerical simulations. Both concentric wrinkles and eye-like folds are shown to be compatible with negative curvatures. Which pattern will be realized depends on the curvature of the substrate. We discuss both types of folding patterns and determine the phase diagram governing their appearance.

preprint2010arXiv

Spectral analysis for the iron-based superconductors: Anisotropic spin fluctuations and fully gapped s^{\pm}-wave superconductivity

Spin fluctuations are considered to be one of the candidates that drive a sign-reversed s^{\pm} superconducting state in the iron pnictides. In the magnetic scenario, whether the spin fluctuation spectrum exhibits certain unique fine structures is an interesting aspect for theoretical study in order to understand experimental observations. We investigate the detailed momentum dependence of the short-range spin fluctuations using a 2-orbital model in the self-consistent fluctuation exchange approximation and find that a common feature of those fluctuations that are capable of inducing a fully gapped s^{\pm} state is the momentum anisotropy with lengthened span along the direction transverse to the antiferromagnetic momentum transfer. Performing a qualitative analysis based on the orbital character and the deviation from perfect nesting of the electronic structure for the 2-orbital and a more complete 5-orbital model, we gain the insight that this type of anisotropic spin fluctuations favor superconductivity due to their enhancement of intra-orbital, but inter-band, pair scattering processes. The momentum anisotropy leads to elliptically shaped magnetic responses which have been observed in inelastic neutron scattering measurements. Meanwhile, our detailed study on the magnetic and the electronic spectrum shows that the dispersion of the magnetic resonance mode in the nearly isotropic s^{\pm} superconducting state exhibits anisotropic propagating behavior in an upward pattern and the coupling of the resonance mode to fermions leads to a dip feature in the spectral function.

preprint2009arXiv

Anisotropic pairing in the iron pnictides

We determine the anisotropy of the spin fluctuation induced pairing gap on the Fermi surface of the FeAs based superconductors as function of the exchange and Hund's coupling $J_{H}$. We find that for sufficiently large $J_{H}$, nearly commensurate magnetic fluctuations yield a fully gapped $s^{\pm}$-pairing state with small anisotropy of the gap amplitude on each Fermi surface sheet, but significant variations of the gap amplitude for different sheets of the Fermi surface. In particular, we obtain the large variation of the gap amplitude on different Fermi surface sheets, as seen in ARPES experiments. For smaller values of Hund's coupling incommensurate magnetic fluctuations yield an $s^{\pm}$-pairing state with line nodes. Such a state is also possible once the anisotropy of the material is reduced and three dimensional effects come into play.

preprint2009arXiv

Orbital coupling and superconductivity in the iron pnictides

We demonstrate that strong inter-orbital interaction is very efficient to achieve superconductivity due to magnetic fluctuations in the iron pnictides. Fermi surface states that are coupled by the antiferromagnetic wave vector are often of different orbital nature, causing pair-hopping interactions between distinct Fe-3d orbitals to become important. Performing a self-consistent FLEX calculation below $T_{c}$ we determine the superconducting order parameter as function of intra- and inter-orbital couplings. We find an $s^{\pm}$-pairing state with $T_{c}\simeq80\mathrm{K}$ for realistic parameters.