Kinetic Transition Networks for the Thomson Problem and Smale's Seventh Problem

Conformal order and Poincaré-Klein mapping underlying electrostatics-driven inhomogeneity in tethered membranes

Understanding the organization of matter under the long-range electrostatic force is a fundamental problem in multiple fields. In this work, based on the electrically charged tethered membrane model, we reveal regular structures underlying the lowest-energy states of inhomogeneously stretched planar lattices by a combination of numerical simulation and analytical geometric analysis. Specifically we show the conformal order characterized by the preserved bond angle in the lattice deformation and reveal the Poincaré-Klein mapping underlying the electrostatics-driven inhomogeneity. The discovery of the Poincaré-Klein mapping, which connects the Poincaré disk and the Klein disk for the hyperbolic plane, implies the connection of long-range electrostatic force and hyperbolic geometry. We also discuss lattices with patterned charges of opposite signs for modulating in-plane inhomogeneity and even creating 3D shapes, which may have a connection to metamaterials design. This work suggests the geometric analysis as a promising approach for elucidating the organization of matter under the long-range force.

Charged elastic rings: deformation and dynamics

We report the counter-intuitive instability of charged elastic rings, and the persistence of sinusoidal deformations in the lowest-energy configurations by the combination of high-precision numerical simulations and analytical perturbation calculation. We also study the dynamical evolution of the charged ring under random disturbance, and reveal the modulation of the dominant frequencies by the electrostatic force. The purely mechanical analysis of the classical ring system presented in this work yields insights into the subtlety of long-range forces in the organization and dynamics of matter.

Configurational entropy of colloidal particles in a confined space

We calculate the configurational entropy of colloidal particles in a confined geometry interacting as hard disks using Monte Carlo integration. In particular, we consider systems with three kinds of boundary conditions: hard, periodic, and spherical. For small to moderate packing fraction ϕ values, we find the entropies per particle for systems with the periodic and the spherical boundary conditions tend to reach the same value with the increase of the particle number N, while that for the system with the hard boundary conditions still has obvious differences compared to them within the studied N range. Surprisingly, despite the small system sizes, the estimated entropies per particle at infinite system size from extrapolations in the periodic and spherical systems are in reasonable agreement with that calculated using thermodynamic integration. Besides, as N increases we find the pair correlation function begins to exhibit similar features as that of a thermally equilibrated hard-disk fluid at the same packing fraction. Our findings may contribute to a better understanding of how the configurational entropy changes with the system size and the influence of boundary conditions, and provide insights relevant to engineering particles in confined spaces.

Tunable order of helically confined charges

We investigate a system of equally charged Coulomb-interacting particles confined to a toroidal helix in the presence of an external electric field. Due to the confinement, the particles experience an effective interaction that oscillates with the particle distance and allows for the existence of stable bound states, despite the purely repulsive character of the Coulomb interaction. We design an order parameter to classify these bound states and use it to identify a structural crossover of the particle order, occurring when the electric field strength is varied. Amorphous particle configurations for a vanishing electric field and crystalline order in the regime of a strong electric field are observed. We study the impact of parameter variations on the particle order and conclude that the crossover occurs for a wide range of parameter values and even holds for different helical systems.

Crystal-like order and defects in metazoan epithelia with spherical geometry

Since Robert Hooke studied cork cell patterns in 1665, scientists have been puzzled by why cells form such ordered structures. The laws underlying this type of organization are universal, and we study them comparing the living and non-living two-dimensional systems self-organizing at the spherical surface. Such-type physical systems often possess trigonal order with specific elongated defects, scars and pleats, where the 5-valence and 7-valence vertices alternate. In spite of the fact that the same physical and topological rules are involved in the structural organization of biological systems, such topological defects were never reported in epithelia. We have discovered them in the follicular spherical epithelium of ascidians that are emerging models in developmental biology. Surprisingly, the considered defects appear in the epithelium even when the number of cells in it is significantly less than the previously known threshold value. We explain this result by differences in the cell sizes and check our hypothesis considering the self-assembly of different random size particles on the spherical surface. Scars, pleats and other complex defects found in ascidian samples can play an unexpected and decisive role in the permanent renewal and reorganization of epithelia, which forms or lines many tissues and organs in metazoans.

Locally and globally optimal configurations of N particles on the sphere with applications in the narrow escape and narrow capture problems

Determination of optimal arrangements of N particles on a sphere is a well-known problem in physics. A famous example of such is the Thomson problem of finding equilibrium configurations of electrical charges on a sphere. More recently, however, similar problems involving other potentials and nonspherical domains have arisen in biophysical systems. Many optimal configurations have previously been computed, especially for the Thomson problem; however, few results exist for potentials that correspond to more applied problems. Here we numerically compute optimal configurations corresponding to the narrow escape and narrow capture problems in biophysics. We provide comprehensive tables of global energy minima for N≤120 and local energy minima for N≤65, and we exclude all saddle points. Local minima up to N=120 are available online.

Command of Collective Dynamics by Topological Defects in Spherical Crystals

Directing individual motions of many constituents to a coherent dynamical state is a fundamental challenge in multiple fields. Here, based on the spherical crystal model, we show that topological defects in particle arrays can be a crucial element in regulating collective dynamics. Specifically, we highlight the defect-driven synchronized breathing modes around disclinations and collective oscillations with strong connection to disruption of crystalline order. This work opens the promising possibility of an organizational principle based on topological defects, and may inspire new strategies for harnessing intriguing collective dynamics in extensive nonequilibrium systems.

Globally optimal volume-trap arrangements for the narrow-capture problem inside a unit sphere

The determination of statistical characteristics for particles undergoing Brownian motion in constrained domains has multiple applications in various areas of research. This work presents an attempt to systematically compute globally optimal configurations of traps inside a three-dimensional domain that minimize the average of the mean first passage time (MFPT) for the narrow capture problem, the average time it takes a particle to be captured by any trap. For a given domain, the mean first passage time satisfies a linear Poisson problem with Dirichlet-Neumann boundary conditions. While no closed-form general solution of such problems is known, approximate asymptotic MFPT expressions for small traps in a unit sphere have been found. These solutions explicitly depend on trap parameters, including locations, through a pairwise potential function. After probing the applicability limits of asymptotic formulas through comparisons with numerical and available exact solutions of the narrow capture problem, full three-dimensional global optimization was performed to find optimal trap positions in the unit sphere for 2≤N≤100 identical traps. The interaction energy values and geometrical features of the putative optimal trap arrangements are presented.

Depletion zones and crystallography on pinched spheres

Understanding the interplay between ordered structures and substrate curvature is an interesting problem with versatile applications, including functionalization of charged supramolecular surfaces and modern microfluidic technologies. In this work, we investigate the two-dimensional packing structures of charged particles confined on a pinched sphere. By continuously pinching the sphere, we observe cleavage of elongated scars into pleats, proliferation of disclinations, and subsequently, emergence of a depletion zone at the negatively curved waist that is completely void of particles. We systematically study the geometrics and energetics of the depletion zone, and reveal its physical origin as a finite size effect, due to the interplay between Coulomb repulsion and concave geometry of the pinched sphere. These results further our understanding of crystallography on curved surfaces, and have implications in design and manipulation of charged, deformable interfaces in various applications.

Energy landscapes for machine learning

Machine learning techniques are being increasingly used as flexible non-linear fitting and prediction tools in the physical sciences. Fitting functions that exhibit multiple solutions as local minima can be analysed in terms of the corresponding machine learning landscape. Methods to explore and visualise molecular potential energy landscapes can be applied to these machine learning landscapes to gain new insight into the solution space involved in training and the nature of the corresponding predictions. In particular, we can define quantities analogous to molecular structure, thermodynamics, and kinetics, and relate these emergent properties to the structure of the underlying landscape. This Perspective aims to describe these analogies with examples from recent applications, and suggest avenues for new interdisciplinary research.

Topological vacancies in spherical crystals

Understanding geometric frustration of ordered phases in two-dimensional condensed matter on curved surfaces is closely related to a host of scientific problems in condensed matter physics and materials science. Here, we show how two-dimensional Lennard-Jones crystal clusters confined on a sphere resolve geometric frustration and form pentagonal vacancy structures. These vacancies, originating from the combination of curvature and physical interaction, are found to be topological defects and they can be further classified into dislocational and disclinational types. We analyze the dual role of these crystallographic defects as both vacancies and topological defects, illustrate their formation mechanism, and present the phase diagram. The revealed dual role of the topological vacancies may find applications in the fabrication of robust nanopores. This work also shows the promising potential of exploiting richness in both physical interactions and substrate geometries to create new types of crystallographic defects, which have strong connections with the design of crystalline materials.

How geometric frustration shapes twisted fibres, inside and out: competing morphologies of chiral filament assembly

Chirality frustrates and shapes the assembly of flexible filaments in rope-like, twisted bundles and fibres by introducing gradients of both filament shape (i.e. curvature) and packing throughout the structure. Previous models of chiral filament bundle formation have shown that this frustration gives rise to several distinct morphological responses, including self-limiting bundle widths, anisotropic domain (tape-like) formation and topological defects in the lateral inter-filament order. In this paper, we employ a combination of continuum elasticity theory and discrete filament bundle simulations to explore how these distinct morphological responses compete in the broader phase diagram of chiral filament assembly. We show that the most generic model of bundle formation exhibits at least four classes of equilibrium structure-finite-width, twisted bundles with isotropic and anisotropic shapes, with and without topological defects, as well as bulk phases of untwisted, columnar assembly (i.e. 'frustration escape'). These competing equilibrium morphologies are selected by only a relatively small number of parameters describing filament assembly: bundle surface energy, preferred chiral twist and stiffness of chiral filament interactions, and mechanical stiffness of filaments and their lateral interactions. Discrete filament bundle simulations test and verify continuum theory predictions for dependence of bundle structure (shape, size and packing defects of two-dimensional cross section) on these key parameters.

Properties of kinetic transition networks for atomic clusters and glassy solids

Small-world and scale-free properties are analysed for kinetic transition networks of clusters and glassy systems.