skip to main content


Title: Snaking bifurcations of localized patterns on ring lattices
Abstract We study the structure of stationary patterns in bistable lattice dynamical systems posed on rings with a symmetric coupling structure in the regime of small coupling strength. We show that sparse coupling (for instance, nearest-neighbour or next-nearest-neighbour coupling) and all-to-all coupling lead to significantly different solution branches. In particular, sparse coupling leads to snaking branches with many saddle-node bifurcations, while all-to-all coupling leads to branches with six saddle nodes, regardless of the size of the number of nodes in the graph.  more » « less
Award ID(s):
1714429
NSF-PAR ID:
10299815
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
IMA Journal of Applied Mathematics
Volume:
86
Issue:
5
ISSN:
0272-4960
Page Range / eLocation ID:
1112 to 1140
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Rectilinear crawling locomotion is a primitive and common mode of locomotion in slender soft-bodied animals. It requires coordinated contractions that propagate along a body that interacts frictionally with its environment. We propose a simple approach to understand how this coordination arises in a neuromechanical model of a segmented, soft-bodied crawler via an iterative process that might have both biological antecedents and technological relevance. Using a simple reinforcement learning algorithm, we show that an initial all-to-all neural coupling converges to a simple nearest-neighbour neural wiring that allows the crawler to move forward using a localized wave of contraction that is qualitatively similar to what is observed in Drosophila melanogaster larvae and used in many biomimetic solutions. The resulting solution is a function of how we weight gait regularization in the reward, with a trade-off between speed and robustness to proprioceptive noise. Overall, our results, which embed the brain–body–environment triad in a learning scheme, have relevance for soft robotics while shedding light on the evolution and development of locomotion. 
    more » « less
  2. Abstract Engineering non-linear hybrid light-matter states in tailored lattices is a central research strategy for the simulation of complex Hamiltonians. Excitons in atomically thin crystals are an ideal active medium for such purposes, since they couple strongly with light and bear the potential to harness giant non-linearities and interactions while presenting a simple sample-processing and room temperature operability. We demonstrate lattice polaritons, based on an open, high-quality optical cavity, with an imprinted photonic lattice strongly coupled to excitons in a WS 2 monolayer. We experimentally observe the emergence of the canonical band-structure of particles in a one-dimensional lattice at room temperature, and demonstrate frequency reconfigurability over a spectral window exceeding 85 meV, as well as the systematic variation of the nearest-neighbour coupling, reflected by a tunability in the bandwidth of the p-band polaritons by 7 meV. The technology presented in this work is a critical demonstration towards reconfigurable photonic emulators operated with non-linear photonic fluids, offering a simple experimental implementation and working at ambient conditions. 
    more » « less
  3. Summary

    We consider sparse spatial mixed linear models, particularly those described by Besag and Higdon, and develop an h-likelihood method for their statistical inference. The method proposed allows for singular precision matrices, as it produces estimates that coincide with those from the residual maximum likelihood based on appropriate differencing of the data and has a novel approach to estimating precision parameters by a gamma linear model. Furthermore, we generalize the h-likelihood method to include continuum spatial variations by making explicit use of scaling limit connections between Gaussian intrinsic Markov random fields on regular arrays and the de Wijs process. Keeping various applications of spatial mixed linear models in mind, we devise a novel sparse conjugate gradient algorithm that allows us to achieve fast matrix-free statistical computations. We provide two applications. The first is an extensive analysis of an agricultural variety trial that brings forward various new aspects of nearest neighbour adjustment such as effects on statistical analyses to changes of scale and use of implicit continuum spatial formulation. The second application concerns an analysis of a large cotton field which gives a focus to matrix-free computations. The paper closes with some further considerations, such as applications to irregularly spaced data, use of the parametric bootstrap and some generalizations to the Gaussian Matérn mixed effect models.

     
    more » « less
  4. Learning representations of sets of nodes in a graph is crucial for applications ranging from node-role discovery to link prediction and molecule classification. Graph Neural Networks (GNNs) have achieved great success in graph representation learning. However, expressive power of GNNs is limited by the 1-Weisfeiler-Lehman (WL) test and thus GNNs generate identical representations for graph substructures that may in fact be very different. More powerful GNNs, proposed recently by mimicking higher-order-WL tests, only focus on representing entire graphs and they are computationally inefficient as they cannot utilize sparsity of the underlying graph. Here we propose and mathematically analyze a general class of structure related features, termed Distance Encoding (DE). DE assists GNNs in representing any set of nodes, while providing strictly more expressive power than the 1-WL test. DE captures the distance between the node set whose representation is to be learned and each node in the graph. To capture the distance DE can apply various graph-distance measures such as shortest path distance or generalized PageRank scores. We propose two ways for GNNs to use DEs (1) as extra node features, and (2) as controllers of message aggregation in GNNs. Both approaches can utilize the sparse structure of the underlying graph, which leads to computational efficiency and scalability. We also prove that DE can distinguish node sets embedded in almost all regular graphs where traditional GNNs always fail. We evaluate DE on three tasks over six real networks: structural role prediction, link prediction, and triangle prediction. Results show that our models outperform GNNs without DE by up-to 15% in accuracy and AUROC. Furthermore, our models also significantly outperform other state-of-the-art methods especially designed for the above tasks. 
    more » « less
  5. SUMMARY

    As the high-frequency analogue to field-scale earthquakes, acoustic emissions (AEs) provide a valuable complement to study rock deformation mechanisms. During the load-stepping creep experiments with CO2-saturated water injection into a basaltic sample from Carbfix site in Iceland, 8791 AE events are detected by at least one of the seven piezoelectric sensors. Here, we apply a cross-correlation-based source imaging method, called geometric-mean reverse-time migration (GmRTM) to locate those AE events. Besides the attractive picking-free feature shared with other waveform-based methods (e.g. time-reversal imaging), GmRTM is advantageous in generating high-resolution source images with reduced imaging artefacts, especially for experiments with relatively sparse receivers. In general, the imaged AE locations are found to be scattered across the sample, suggesting a complicated fracture network rather than a well-defined major shear fracture plane, in agreement with X-ray computed tomography imaging results after retrieval of samples from the deformation apparatus. Clustering the events in space and time using the nearest-neighbour approach revealed a group of ‘repeaters’, which are spatially co-located over an elongated period of time and likely indicate crack, or shear band growth. Furthermore, we select 2196 AE events with high signal-to-noise-ratio (SNR) and conduct moment tensor estimation using the adjoint (backpropagated) strain tensor fields at the locations of AE sources. The resulting AE locations and focal mechanisms support our previously assertion that creep of basalt at the experimental conditions is accommodated dominantly by distributed microcracking.

     
    more » « less