skip to main content


Title: Local Stochastic Algorithms for Alignment in Self-Organizing Particle Systems
We present local distributed, stochastic algorithms for alignment in self-organizing particle systems (SOPS) on two-dimensional lattices, where particles occupy unique sites on the lattice, and particles can make spatial moves to neighboring sites if they are unoccupied. Such models are abstractions of programmable matter, composed of individual computational particles with limited memory, strictly local communication abilities, and modest computational capabilities. We consider oriented particle systems, where particles are assigned a vector pointing in one of q directions, and each particle can compute the angle between its direction and the direction of any neighboring particle, although without knowledge of global orientation with respect to a fixed underlying coordinate system. Particles move stochastically, with each particle able to either modify its direction or make a local spatial move along a lattice edge during a move. We consider two settings: (a) where particle configurations must remain simply connected at all times and (b) where spatial moves are unconstrained and configurations can disconnect. Our algorithms are inspired by the Potts model and its planar oriented variant known as the planar Potts model or clock model from statistical physics. We prove that for any q ≥ 2, by adjusting a single parameter, these self-organizing particle systems can be made to collectively align along a single dominant direction (analogous to a solid or ordered state) or remain non-aligned, in which case the fraction of particles oriented along any direction is nearly equal (analogous to a gaseous or disordered state). In the connected SOPS setting, we allow for two distinct parameters, one controlling the ferromagnetic attraction between neighboring particles (regardless of orientation) and the other controlling the preference of neighboring particles to align. We show that with appropriate settings of the input parameters, we can achieve compression and expansion, controlling how tightly gathered the particles are, as well as alignment or nonalignment, producing a single dominant orientation or not. While alignment is known for the Potts and clock models at sufficiently low temperatures, our proof in the SOPS setting are significantly more challenging because the particles make spatial moves, not all sites are occupied, and the total number of particles is fixed.  more » « less
Award ID(s):
2106687
NSF-PAR ID:
10425738
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Approximation, Randomization, and Combinatorial Optimization. Algorithms and Techniques (APPROX/RANDOM 2022)
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. A ring with a cross-section that has a blunt inner and sharper outer edge can attain an equilibrium orientation in a Newtonian fluid subject to a low Reynolds number simple shear flow. This may be contrasted with the continuous rotation exhibited by most rigid bodies. Such rings align along an orientation when the rotation due to fluid vorticity balances the counter-rotation due to the extensional component of the simple shear flow. While the viscous stress on the particle tries to rotate it, the pressure can generate a counter-vorticity torque that aligns the particle. Using boundary integral computations, we demonstrate ways to effectively control this pressure by altering the geometry of the ring cross-section, thus leading to alignment at moderate particle aspect ratios. Aligning rings that lack fore–aft symmetry can migrate indefinitely along the gradient direction. This differs from the periodic spatial trajectories of fore–aft asymmetric axisymmetric particles that rotate in periodic orbits. The mechanism for migration of aligned rings along the gradient direction is elucidated in this work. The migration speed can be controlled by varying the cross-sectional shape and size of the ring. Our results provide new insights into controlling motion of individual particles and thereby open new pathways towards manipulating macroscopic properties of a suspension. 
    more » « less
  2. Using a numerical model, we analyse the effects of shape on both the orientation and transport of anisotropic particles in wavy flows. The particles are idealized as prolate and oblate spheroids, and we consider the regime of small Stokes and particle Reynolds numbers. We find that the particles preferentially align into the shear plane with a mean orientation that is solely a function of their aspect ratio. This alignment, however, differs from the Jeffery orbits that occur in the residual shear flow (that is, the Stokes drift velocity field) in the absence of waves. Since the drag on an anisotropic particle depends on its alignment with the flow, this preferred orientation determines the effective drag on the particles, which in turn impacts their net downstream transport. We also find that the rate of alignment of the particles is not constant and depends strongly on their initial orientation; thus, variations in initial particle orientation result in dispersion of anisotropic-particle plumes. We show that this dispersion is a function of the particle’s eccentricity and the ratio of the settling and wave time scales. Due to this preferential alignment, we find that a plume of anisotropic particles in waves is on average transported farther but dispersed less than it would be if the particles were randomly oriented. Our results demonstrate that accurate prediction of the transport of anisotropic particles in wavy environments, such as microplastic particles in the ocean, requires the consideration of these preferential alignment effects. 
    more » « less
  3. We investigate algorithmic control of a large swarm of mobile particles (such as robots, sensors, or building material) that move in a 2D workspace using a global input signal (such as gravity or a magnetic field). Upon activation of the field, each particle moves maximally in the same direction until forward progress is blocked by a stationary obstacle or another stationary particle. In an open workspace, this system model is of limited use because it has only two controllable degrees of freedom—all particles receive the same inputs and move uniformly. We show that adding a maze of obstacles to the environment can make the system drastically more complex but also more useful. We provide a wide range of results for a wide range of questions. These can be subdivided into external algorithmic problems, in which particle configurations serve as input for computations that are performed elsewhere, and internal logic problems, in which the particle configurations themselves are used for carrying out computations. For external algorithms, we give both negative and positive results. If we are given a set of stationary obstacles, we prove that it is NP-hard to decide whether a given initial configuration of unit-sized particles can be transformed into a desired target configuration. Moreover, we show that finding a control sequence of minimum length is PSPACE-complete. We also work on the inverse problem, providing constructive algorithms to design workspaces that efficiently implement arbitrary permutations between different configurations. For internal logic, we investigate how arbitrary computations can be implemented. We demonstrate how to encode dual-rail logic to build a universal logic gate that concurrently evaluates AND, NAND, NOR, and OR operations. Using many of these gates and appropriate interconnects, we can evaluate any logical expression. However, we establish that simulating the full range of complex interactions present in arbitrary digital circuits encounters a fundamental difficulty: a FAN-OUT gate cannot be generated. We resolve this missing component with the help of 2 9 1 particles, which can create FAN-OUT gates that produce multiple copies of the inputs. Using these gates we provide rules for replicating arbitrary digital circuits. 
    more » « less
  4. null (Ed.)
    Abstract Background Cryo-EM data generated by electron tomography (ET) contains images for individual protein particles in different orientations and tilted angles. Individual cryo-EM particles can be aligned to reconstruct a 3D density map of a protein structure. However, low contrast and high noise in particle images make it challenging to build 3D density maps at intermediate to high resolution (1–3 Å). To overcome this problem, we propose a fully automated cryo-EM 3D density map reconstruction approach based on deep learning particle picking. Results A perfect 2D particle mask is fully automatically generated for every single particle. Then, it uses a computer vision image alignment algorithm (image registration) to fully automatically align the particle masks. It calculates the difference of the particle image orientation angles to align the original particle image. Finally, it reconstructs a localized 3D density map between every two single-particle images that have the largest number of corresponding features. The localized 3D density maps are then averaged to reconstruct a final 3D density map. The constructed 3D density map results illustrate the potential to determine the structures of the molecules using a few samples of good particles. Also, using the localized particle samples (with no background) to generate the localized 3D density maps can improve the process of the resolution evaluation in experimental maps of cryo-EM. Tested on two widely used datasets, Auto3DCryoMap is able to reconstruct good 3D density maps using only a few thousand protein particle images, which is much smaller than hundreds of thousands of particles required by the existing methods. Conclusions We design a fully automated approach for cryo-EM 3D density maps reconstruction (Auto3DCryoMap). Instead of increasing the signal-to-noise ratio by using 2D class averaging, our approach uses 2D particle masks to produce locally aligned particle images. Auto3DCryoMap is able to accurately align structural particle shapes. Also, it is able to construct a decent 3D density map from only a few thousand aligned particle images while the existing tools require hundreds of thousands of particle images. Finally, by using the pre-processed particle images,Auto3DCryoMap reconstructs a better 3D density map than using the original particle images. 
    more » « less
  5. Polymer-networked nanoparticles are the basis for advanced materials useful wearable electron- ics, drug delivery, autonomous computing and other applications. To characterize and predict the physics and underlying mechanisms of the network connections in 2D and 3D engineered nanopar- ticle (ENP) arrays, we developed an analogous Potts model of 3-state sites. Together with dissipa- tive particle dynamics (DPD) simulations, we found that the network structures in polymer-linked nanoparticle assemblies are generally dominated by the number of nearest neighbors and not the topology of the lattice. When the E-field regulates the network connections, the links along the E-field direction always dominate the overall network structure. 
    more » « less