Researcher profile

Rastko Sknepnek

Rastko Sknepnek contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

11 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.

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

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.