Search

This paper studies the distributed mode consensus problem in a multi-agent system, in which the agents each possess a certain attribute and they aim to agree upon the mode (the most frequent attribute owned by the agents) via distributed computation. Three algorithms are proposed. The first one directly calculates the frequency of all attributes at every agent, with protocols based on blended dynamics, and then returns the most frequent attribute as the mode. Assuming knowledge at each agent of a lower bound of the mode frequency as a priori information, the second algorithm is able to reduce the number of frequencies to be computed at every agent if the lower bound is large. The third algorithm further eliminates the need for this information by introducing an adaptive updating mechanism. The algorithms find the mode in finite time, and estimates of convergence time are provided. The proposed first and second algorithms enjoy the plug-and-play property with a dwell time.

This study introduces an approach to obtain a neighboring extremal optimal control (NEOC) solution for a closed-loop optimal control problem, applicable to a wide array of nonlinear systems and not necessarily quadratic performance indices. The approach involves investigating the variation incurred in the functional form of a known closed-loop optimal control law due to small, known parameter variations in the system equations or the performance index. The NEOC solution can formally be obtained by solving a linear partial differential equation, akin to those encountered in the iterative solution of a nonlinear Hamilton-Jacobi equation. Motivated by numerical procedures for solving these latter equations, we also propose a numerical algorithm based on the Galerkin algorithm, leveraging the use of basis functions to solve the underlying Hamilton-Jacobi equation of the original optimal control problem. The proposed approach simplifies the NEOC problem by reducing it to the solution of a simple set of linear equations, thereby eliminating the need for a full re-solution of the adjusted optimal control problem. Furthermore, the variation to the optimal performance index can be obtained as a function of both the system state and small changes in parameters, allowing the determination of the adjustment to an optimal control law given a small adjustment of parameters in the system or the performance index. Moreover, in order to handle large known parameter perturbations, we propose a homotopic approach that breaks down the single calculation of NEOC into a finite set of multiple steps. Finally, the validity of the claims and theory is supported by theoretical analysis and numerical simulations.

The paper deals with the spread of two competing viruses over a network of population nodes, accounting for pairwise interactions and higher-order interactions (HOI) within and between the population nodes. We study the competitive networked bivirus susceptible-infected-susceptible (SIS) model on a hypergraph introduced in Cui et al. [1]. We show that the system has, in a generic sense, a finite number of equilibria, and the Jacobian associated with each equilibrium point is nonsingular; the key tool is the Parametric Transversality Theorem of differential topology. Since the system is also monotone, it turns out that the typical behavior of the system is convergence to some equilibrium point. Thereafter, we exhibit a tri-stable domain with three locally exponentially stable equilibria. For different parameter regimes, we establish conditions for the existence of a coexistence equilibrium (both viruses infect separate fractions of each population node).

This paper considers the classical Susceptible--Infected--Susceptible (SIS) network epidemic model, which describes a disease spreading through $n$ nodes, with the network links governing the possible transmission pathways of the disease between nodes. We consider feedback control to eliminate the disease in scenarios where the disease would otherwise persist in an uncontrolled network. We propose a family of decentralised adaptive-gain control algorithms, in which each node has a control gain that adaptively evolves according to a differential equation, independent of the gains of other nodes. The adaptive gain is applied multiplicatively to either decrease the infection rate or increase the recovery rate. To begin, we assume all nodes are controlled, and prove that both infection rate control and recovery rate control algorithms eliminate the disease with the limiting gains being positive and finite. Then, we consider the possibility of controlling a subset of the nodes, for both the infection rate control and recovery rate control. We first identify a necessary and sufficient condition for the existence of a subset of nodes, which if controlled would result in the elimination of the disease. For a given network, there may exist several such viable subsets, and we propose an iterative algorithm to identify such a subset. Simulations are provided to demonstrate the effectiveness of the various proposed controllers., Comment: Extended version of a journal paper submission

Neuromorphic computing shows promise for advancing computing efficiency and capabilities of AI applications using brain-inspired principles. However, the neuromorphic research field currently lacks standardized benchmarks, making it difficult to accurately measure technological advancements, compare performance with conventional methods, and identify promising future research directions. Prior neuromorphic computing benchmark efforts have not seen widespread adoption due to a lack of inclusive, actionable, and iterative benchmark design and guidelines. To address these shortcomings, we present NeuroBench: a benchmark framework for neuromorphic computing algorithms and systems. NeuroBench is a collaboratively-designed effort from an open community of researchers across industry and academia, aiming to provide a representative structure for standardizing the evaluation of neuromorphic approaches. The NeuroBench framework introduces a common set of tools and systematic methodology for inclusive benchmark measurement, delivering an objective reference framework for quantifying neuromorphic approaches in both hardware-independent (algorithm track) and hardware-dependent (system track) settings. In this article, we outline tasks and guidelines for benchmarks across multiple application domains, and present initial performance baselines across neuromorphic and conventional approaches for both benchmark tracks. NeuroBench is intended to continually expand its benchmarks and features to foster and track the progress made by the research community., Comment: System track baselines added

This paper studies the endemic behavior of a multi-competitive networked susceptible-infected-susceptible (SIS) model. Specifically, the paper deals with three competing virus systems (i.e., tri-virus systems). First, we show that a tri-virus system, unlike a bi-virus system, is not a monotone dynamical system. Using the Parametric Transversality Theorem, we show that, generically, a tri-virus system has a finite number of equilibria and that the Jacobian matrices associated with each equilibrium are nonsingular. The endemic equilibria of this system can be classified as follows: a) single-virus endemic equilibria (also referred to as the boundary equilibria), where precisely one of the three viruses is alive; b) 2-coexistence equilibria, where exactly two of the three viruses are alive; and c) 3-coexistence equilibria, where all three viruses survive in the network. We provide a necessary and sufficient condition that guarantees local exponential convergence to a boundary equilibrium. Further, we secure conditions for the nonexistence of 3-coexistence equilibria (resp. for various forms of 2-coexistence equilibria). We also identify sufficient conditions for the existence of a 2-coexistence (resp. 3-coexistence) equilibrium. We identify conditions on the model parameters that give rise to a continuum of coexistence equilibria. More specifically, we establish i) a scenario that admits the existence and local exponential attractivity of a line of coexistence equilibria; and ii) scenarios that admit the existence of, and, in the case of one such scenario, global convergence to, a plane of 3-coexistence equilibria., Comment: arXiv admin note: substantial text overlap with arXiv:2209.11826

This paper investigates the problem of securing exponentially fast consensus (exponential consensus for short) for identical agents with finite-dimensional linear system dynamics over dynamic network topology. Our aim is to find the weakest possible conditions that guarantee exponentially fast consensus using a Lyapunov function consisting of a sum of terms of the same functional form. We first investigate necessary conditions, starting by examining the system (both agent and network) parameters. It is found that controllability of the linear agents is necessary for reaching consensus. Then, to work out necessary conditions incorporating the network topology, we construct a set of Laplacian matrix-valued functions. The precompactness of this set of functions is shown to be a significant generalization of existing assumptions on network topology, including the common assumption that the edge weights are bounded piecewise constant functions or continuous functions. With the aid of such a precompactness assumption and restricting the Lyapunov function to one consisting of a sum of terms of the same functional form, we prove that a joint $(\delta, T)$-connectivity condition on the network topology is necessary for exponential consensus. Finally, we investigate how the above two ``necessities'' work together to guarantee exponential consensus. To partially address this problem, we define a synchronization index to characterize the interplay between agent parameters and network topology. Based on this notion, it is shown that by designing a proper feedback matrix and under the precompactness assumption, exponential consensus can be reached globally and uniformly if the joint $(\delta,T)$-connectivity and controllability conditions are satisfied, and the synchronization index is not less than one.

This note presents several conditions to characterize real matrix similarity between a Hurwitz matrix (and then more generally, a real square matrix) and a diagonal dominant matrix.

In the vector-field guided path-following problem, a sufficiently smooth vector field is designed such that its integral curves converge to and move along a one-dimensional geometric desired path. The existence of singular points where the vector field vanishes creates a topological obstruction to global convergence to the desired path and some associated topological analysis has been conducted in our previous work. In this paper, we strengthen the result in our previous work by showing that the domain of attraction of the desired path, which is a compact asymptotically stable one-dimensional embedded submanifold of an $n$-dimensional ambient manifold $\mathcal{M}$, is homeomorphic to $\mathbb{R}^{n-1} \times \mathbb{S}^1$, and not just homotopy equivalent to $\mathbb{S}^1$. This result is extended for a $k$-dimensional compact manifold for $k \ge 2$.

Future missions to Enceladus would benefit from multi-moon tours that leverage V-infinity on resonant orbits to progressively transfer between moons. Such "resonance family hopping" trajectories present a vast search space for global optimization due to the different combinations of available resonances and flyby speeds. The proposed multi-objective tour design algorithm optimizes entire moon tours from Titan to Enceladus via grid-based dynamic programming, in which the computation time is significantly reduced by utilizing a database of V-infinity-leveraging transfers. The result unveils a complete trade space of the moon tour design to Enceladus in a tractable computation time and global optimality., Comment: This article was presented at the 33rd AAS/AIAA Space Flight Mechanics Meeting, Austin, TX 2023 (paper ID:23-313). The results in this article are obsolete, and readers are encouraged to peruse the alternate version of the article published in Celestial Mechanics and Dynamical Astronomy (https://doi.org/10.1007/s10569-023-10179-8). 19 pages, 11 figures, 8 tables

This paper considers a deterministic Susceptible-Infected-Susceptible (SIS) networked bivirus epidemic model (termed the bivirus model for short), in which two competing viruses spread through a set of populations (nodes) connected by two graphs, which may be different if the two viruses have different transmission pathways. The networked dynamics can give rise to complex equilibria patterns, and most current results identify conditions on the model parameters for convergence to the healthy equilibrium (where both viruses are extinct) or a boundary equilibrium (where one virus is endemic and the other is extinct). However, there are only limited results on coexistence equilibria (where both viruses are endemic). This paper establishes a set of ``counting'' results which provide lower bounds on the number of coexistence equilibria, and perhaps more importantly, establish properties on the local stability/instability properties of these equilibria. In order to do this, we employ the Poincar\'e-Hopf Theorem but with significant modifications to overcome several challenges arising from the bivirus system model, such as the fact that the system dynamics do not evolve on a manifold in the typical sense required to apply Poincar\'e-Hopf Theory. Subsequently, Morse inequalities are used to tighten the counting results, under the reasonable assumption that the bivirus system is a Morse-Smale dynamical system. Numerical examples are provided which demonstrate the presence of multiple attractor equilibria, and multiple coexistence equilibria., Comment: Accepted version of paper to appear in SIAM Journal of Applied Dynamical Systems

Although generative adversarial networks (GANs) have shown promise in medical imaging, they have four main limitations that impeded their utility: computational cost, data requirements, reliable evaluation measures, and training complexity. Our work investigates each of these obstacles in a novel application of StyleGAN2-ADA to high-resolution medical imaging datasets. Our dataset is comprised of liver-containing axial slices from non-contrast and contrast-enhanced computed tomography (CT) scans. Additionally, we utilized four public datasets composed of various imaging modalities. We trained a StyleGAN2 network with transfer learning (from the Flickr-Faces-HQ dataset) and data augmentation (horizontal flipping and adaptive discriminator augmentation). The network's generative quality was measured quantitatively with the Fr\'echet Inception Distance (FID) and qualitatively with a visual Turing test given to seven radiologists and radiation oncologists. The StyleGAN2-ADA network achieved a FID of 5.22 ($\pm$ 0.17) on our liver CT dataset. It also set new record FIDs of 10.78, 3.52, 21.17, and 5.39 on the publicly available SLIVER07, ChestX-ray14, ACDC, and Medical Segmentation Decathlon (brain tumors) datasets. In the visual Turing test, the clinicians rated generated images as real 42% of the time, approaching random guessing. Our computational ablation study revealed that transfer learning and data augmentation stabilize training and improve the perceptual quality of the generated images. We observed the FID to be consistent with human perceptual evaluation of medical images. Finally, our work found that StyleGAN2-ADA consistently produces high-quality results without hyperparameter searches or retraining., Comment: This preprint has not undergone post-submission improvements or corrections. The Version of Record of this contribution is published in LNCS, volume 13570, and is available online at https://doi.org/10.1007/978-3-031-16980-9_14

This paper studies a distributed state estimation problem for both continuous- and discrete-time linear systems. A simply structured distributed estimator (comprising interconnected local estimators) is first described for estimating the state of a continuous and multi-channel linear system whose sensed outputs are distributed across a fixed multi-agent network. The estimator is then extended to non-stationary networks whose graphs switch according to a switching signal. The estimator is guaranteed to solve the problem, provided a network-widely shared high gain condition achieving a form of spectrum separation is satisfied. As an alternative to sharing a common gain across the network, a fully distributed version of the estimator is also studied in which each agent adaptively adjusts a local gain, though the practicality of this approach is subject to a robustness issue common to adaptive control. A discrete-time version of the distributed state estimation problem is also studied, and a corresponding estimator based again on spectrum separation, but not high gain, is proposed for time-varying networks. For each scenario, it is explained how to construct the estimator so that the state estimation errors in the local estimators all converge to zero exponentially fast at a fixed but arbitrarily chosen rate, provided the network's graph is strongly connected for all time. The proposed estimators are inherently resilient to abrupt changes in the number of agents and communication links in the inter-agent communication graph upon which the algorithms depend, provided the network is redundantly strongly connected and redundantly jointly observable., Comment: 17 pages, 6 figures. arXiv admin note: substantial text overlap with arXiv:1903.05486

This paper investigates the problem of cooperative tuning of multi-agent optimal control systems, where a network of agents (i.e. multiple coupled optimal control systems) adjusts parameters in their dynamics, objective functions, or controllers in a coordinated way to minimize the sum of their loss functions. Different from classical techniques for tuning parameters in a controller, we allow tunable parameters appearing in both the system dynamics and the objective functions of each agent. A framework is developed to allow all agents to reach a consensus on the tunable parameter, which minimizes team loss. The key idea of the proposed algorithm rests on the integration of consensus-based distributed optimization for a multi-agent system and a gradient generator capturing the optimal performance as a function of the parameter in the feedback loop tuning the parameter for each agent. Both theoretical results and simulations for a synchronous multi-agent rendezvous problem are provided to validate the proposed method for cooperative tuning of multi-agent optimal control.

This paper studies the endemic behavior of a multi-competitive networked susceptible-infected-susceptible (SIS) model. In particular, we focus on the case where there are three competing viruses (i.e., the tri-virus system). First, we show that the tri-virus system is not a monotone system. Thereafter, we provide a condition that guarantees local exponential convergence to a boundary equilibrium (exactly one virus is endemic, the other two are dead), and identify a special case that admits the existence and local exponential attractivity of a line of coexistence equilibria (at least two viruses are active). Finally, we identify a particular case (subsumed by the aforementioned special case) such that, for all nonzero initial infection levels, the dynamics of the tri-virus system converge to a plane of coexistence equilibria.

Accurately following a geometric desired path in a two-dimensional space is a fundamental task for many engineering systems, in particular mobile robots. When the desired path is occluded by obstacles, it is necessary and crucial to temporarily deviate from the path for obstacle/collision avoidance. In this paper, we develop a composite guiding vector field via the use of smooth bump functions, and provide theoretical guarantees that the integral curves of the vector field can follow an arbitrary sufficiently smooth desired path and avoid collision with obstacles of arbitrary shapes. These two behaviors are reactive since path (re)-planning and global map construction are not involved. To deal with the common deadlock problem, we introduce a switching vector field, and the Zeno behavior is excluded. Simulations are conducted to support the theoretical results.

High-Dimensional Dynamic Factor Models are presented in detail: The main assumptions and their motivation, main results, illustrations by means of elementary examples. In particular, the role of singular ARMA models in the theory and applications of High-Dimensional Dynamic Factor Models is discussed.The emphasis of the paper is on model classes and their structure theory, rather than on estimation in the narrow sense. Our aim is not a comprehensive survey. Rather we try to point out promising lines of research and applications that have not yet been sufficiently developed., Comment: 28 pages

A path-following control algorithm enables a system's trajectories under its guidance to converge to and evolve along a given geometric desired path. There exist various such algorithms, but many of them can only guarantee local convergence to the desired path in its neighborhood. In contrast, the control algorithms using a well-designed guiding vector field can ensure almost global convergence of trajectories to the desired path; here, "almost" means that in some cases, a measure-zero set of trajectories converge to the singular set where the vector field becomes zero (with all other trajectories converging to the desired path). In this paper, we first generalize the guiding vector field from the Euclidean space to a general smooth Riemannian manifold. This generalization can deal with path-following in some abstract configuration space (such as robot arm joint space). Then we show several theoretical results from a topological viewpoint. Specifically, we are motivated by the observation that singular points of the guiding vector field exist in many examples where the desired path is homeomorphic to the unit circle, but it is unknown whether the existence of singular points always holds in general (i.e., is inherent in the topology of the desired path). In the $n$-dimensional Euclidean space, we provide an affirmative answer, and conclude that it is not possible to guarantee global convergence to desired paths that are homeomorphic to the unit circle. Furthermore, we show that there always exist \emph{non-path-converging trajectories} (i.e., trajectories that do not converge to the desired path) starting from the boundary of a ball containing the desired path in an $n$-dimensional Euclidean space where $n \ge 3$. Examples are provided to illustrate the theoretical results., Comment: 16 pages, 8 figures, TAC

We consider the problem of cooperative motion coordination for multiple heterogeneous mobile vehicles subject to various constraints. These include nonholonomic motion constraints, constant speed constraints, holonomic coordination constraints, and equality/inequality geometric constraints. We develop a general framework involving differential-algebraic equations and viability theory to determine coordination feasibility for a coordinated motion control under heterogeneous vehicle dynamics and different types of coordination task constraints. If a coordinated motion solution exists for the derived differential-algebraic equations and/or inequalities, a constructive algorithm is proposed to derive an equivalent dynamical system that generates a set of feasible coordinated motions for each individual vehicle. In case studies on coordinating two vehicles, we derive analytical solutions to motion generation for two-vehicle groups consisting of car-like vehicles, unicycle vehicles, or vehicles with constant speeds, which serve as benchmark coordination tasks for more complex vehicle groups. The motion generation algorithm is well-backed by simulation data for a wide variety of coordination situations involving heterogeneous vehicles. We then extend the vehicle control framework to deal with the cooperative coordination problem with time-varying coordination tasks and leader-follower structure. We show several simulation experiments on multi-vehicle coordination under various constraints to validate the theory and the effectiveness of the proposed schemes., Comment: 23 pages, 4 figures. Extended version of the paper at IEEE ICRA. Text overlap with arXiv:1809.05509. Submitted to an IEEE journal for publication

The large-scale expansion dynamics of quantum gases is a central tool for ultracold gas experiments and poses a significant challenge for theory. In this work we provide an exact reformulation of the Gross-Pitaevskii equation for the ultracold Bose gas in a coordinate frame that adaptively scales with the system size during evolution, enabling simulations of long evolution times during expansion or similar large-scale manipulation. Our approach makes no hydrodynamic approximations, is not restricted to a scaling ansatz, harmonic potentials, or energy eigenstates, and can be generalized readily to non-contact interactions via the appropriate stress tensor of the quantum fluid. As applications, we simulate the expansion of the ideal gas, a cigar-shaped condensate in the Thomas-Fermi regime, and a linear superposition of counter propagating Gaussian wavepackets. We recover known scaling for the ideal gas and Thomas-Fermi regimes, and identify a linear regime of aspect-ratio preserving free expansion; analysis of the scaling dynamics equations shows that an exact, aspect-ratio invariant, free expansion does not exist for nonlinear evolution. Our treatment enables exploration of nonlinear effects in matter-wave dynamics over large scale-changing evolution., Comment: 12 pages, 3 figures, 2 appendices

Inspired and underpinned by the idea of integral feedback, a distributed constant gain algorithm is proposed for multi-agent networks to solve convex optimization problems with local linear constraints. Assuming agent interactions are modeled by an undirected graph, the algorithm is capable of achieving the optimum solution with an exponential convergence rate. Furthermore, inherited from the beneficial integral feedback, the proposed algorithm has attractive requirements on communication bandwidth and good robustness against disturbance. Both analytical proof and numerical simulations are provided to validate the effectiveness of the proposed distributed algorithms in solving constrained optimization problems., Comment: 12 pages

This paper studies a networked bivirus model, in which two competing viruses spread across a network of interconnected populations; each node represents a population with a large number of individuals. The viruses may spread through possibly different network structures, and an individual cannot be simultaneously infected with both viruses. Focusing on convergence and equilibria analysis, a number of new results are provided. First, we show that for networks with generic system parameters, there exist a finite number of equilibria. Exploiting monotone systems theory, we further prove that for bivirus networks with generic system parameters, then convergence to an equilibrium occurs for all initial conditions, except possibly for a set of measure zero. Given the network structure of one virus, a method is presented to construct an infinite family of network structures for the other virus that results in an infinite number of equilibria in which both viruses coexist. Necessary and sufficient conditions are derived for the local stability/instability of boundary equilibria, in which one virus is present and the other is extinct. A sufficient condition for a boundary equilibrium to be almost globally stable is presented. Then, we show how to use monotone systems theory to generate conclusions on the ordering of stable and unstable equilibria, and in some instances identify the number of equilibria via rapid simulation testing. Last, we provide an analytical method for computing equilibria in networks with only two nodes, and show that it is possible for a bivirus network to have an unstable coexistence equilibrium and two locally stable boundary equilibria., Comment: This is a revised manuscript currently under peer review in SIAM Journal on Control and Optimization (SICON). First version submitted to SICON on 14-Oct-2020

We use a deterministic model to study two competing viruses spreading over a two-layer network in the Susceptible--Infected--Susceptible (SIS) framework, and address a central problem of identifying the winning virus in a "survival-of-the-fittest" battle. Existing sufficient conditions ensure that the same virus always wins regardless of initial states. For networks with an arbitrary but finite number of nodes, there exists a necessary and sufficient condition that guarantees local exponential stability of the two equilibria corresponding to each virus winning the battle, meaning that either of the viruses can win, depending on the initial states. However, establishing existence and finding examples of networks with more than three nodes that satisfy such a condition has remained unaddressed. In this paper, we prove that, for any arbitrary number of nodes, such networks exist. We do this by proving that given almost any network layer of one virus, there exists a network layer for the other virus such that the resulting two-layer network satisfies the aforementioned condition. To operationalize our findings, a four-step procedure is developed to reliably and consistently design one of the network layers, when given the other layer. Conclusions from numerical case studies, including a real-world mobility network that captures the commuting patterns for people between $107$ provinces in Italy, extend on the theoretical result and its consequences., Comment: Submitted to Physical Review E on 2021-Nov-02. Revised version submitted 2022-Sept-22. New version submitted to Systems and Control Letters 2023-Jan-19

We experimentally and numerically demonstrate a method to deterministically generate multiply-quantized superfluid circulation about an obstacle in highly oblate Bose-Einstein condensates (BECs). Our method involves spiraling a blue-detuned laser beam, which acts as a repulsive optical potential, around and towards the center of the BEC. This optical potential serves first as a repulsive stirrer to initiate superflow within the BEC, and then as a pinning potential that transports the center of the superfluid circulation to the center of the condensate. By changing the rate at which the beam moves along the spiral trajectory, we selectively control the net circulation introduced into the BEC. We experimentally achieve pinned superflow with winding numbers as high as 11, which persists for at least 4 s. At the end of the spiral trajectory, with the pinning beam on at full power, the BEC has a toroidal geometry with a high winding-number persistent current. Alternatively, the beam power can be ramped off, allowing controlled placement of a cluster of singly-quantized vortices of the same circulation. This technique can serve as a building block in future experimental architectures to create on-demand vortex distributions and superfluid circulation in BECs., Comment: 10 pages, 10 figures

This paper is concerned with a formation shaping problem for point agents in a two-dimensional space, where control avoids the possibility of reflection ambiguities. One solution for this type of problems was given first for three or four agents by considering a potential function which consists of both the distance error and the signed area terms. Then, by exploiting a hierarchical control strategy with such potential functions, the method was extended to any number of agents recently. However, a specific gain on the signed area term must be employed there, and it does not guarantee the global convergence. To overcome this issue, this paper provides a necessary and sufficient condition for the global convergence, subject to the constraint that the desired formation consists of isosceles triangles only. This clarifies the admissible range of the gain on the signed area for this case. In addition, as for formations consisting of arbitrary triangles, it is shown when high gain on the signed area is admissible for global convergence., Comment: Accepted in the 59th IEEE Conference on Decision and Control (CDC 2020). arXiv admin note: text overlap with arXiv:1808.00312

Ionospheric conductance is a crucial factor in regulating the closure of magnetospheric field-aligned currents through the ionosphere as Hall and Pedersen currents. Despite its importance in predictive investigations of the magnetosphere - ionosphere coupling, the estimation of ionospheric conductance in the auroral region is precarious in most global first-principles based models. This impreciseness in estimating the auroral conductance impedes both our understanding and predictive capabilities of the magnetosphere-ionosphere system during extreme space weather events. In this article, we address this concern, with the development of an advanced Conductance Model for Extreme Events (CMEE) that estimates the auroral conductance from field aligned current values. CMEE has been developed using nonlinear regression over a year's worth of one-minute resolution output from assimilative maps, specifically including times of extreme driving of the solar wind-magnetosphere-ionosphere system. The model also includes provisions to enhance the conductance in the aurora using additional adjustments to refine the auroral oval. CMEE has been incorporated within the Ridley Ionosphere Model (RIM) of the Space Weather Modeling Framework (SWMF) for usage in space weather simulations. This paper compares performance of CMEE against the existing conductance model in RIM, through a validation process for six space weather events. The performance analysis indicates overall improvement in the ionospheric feedback to ground-based space weather forecasts. Specifically, the model is able to improve the prediction of ionospheric currents which impact the simulated dB/dt and {\Delta}B, resulting in substantial improvements in dB/dt predictive skill.

We discuss the design of quantum hybrid inertial sensor that combines an optomechanical inertial sensor with the retro-reflector of a cold atom interferometer. This sensor fusion approach provides absolute and high accuracy measurements with cold atom interferometers, while utilizing the optomechanical inertial sensor at frequencies above the repetition rate of the atom interferometer. This improves the overall measurement bandwidth as well as the robustness and field deployment capabilities of these systems. We evaluate which parameters yield an optimal acceleration sensitivity, from which we anticipate a noise floor at nano-g levels from DC to 1 kHz.

Cortical regions without direct neuronal connections have been observed to exhibit synchronized dynamics. A recent empirical study has further revealed that such regions that share more common neighbors are more likely to behave coherently. To analytically investigate the underlying mechanisms, we consider that a set of n oscillators, which have no direct connections, are linked through m intermediate oscillators (called mediators), forming a complete bipartite network structure. Modeling the oscillators by the Kuramoto-Sakaguchi model, we rigorously prove that mediated remote synchronization, i.e., synchronization between those n oscillators that are not directly connected, becomes more robust as the number of mediators increases. Simulations are also carried out to show that our theoretical findings can be applied to other general and complex networks.

This work considers the distributed computation of the one-to-one vertex correspondences between two undirected and connected graphs, which is called \textit{graph matching}, over multi-agent networks. Given two \textit{isomorphic} and \textit{asymmetric} graphs, there is a unique permutation matrix that maps the vertices in one graph to the vertices in the other. Based on a convex relaxation of graph matching in Aflalo et al. (2015), we propose a distributed computation of graph matching as a distributed convex optimization problem subject to equality constraints and a global set constraint, using a network of multiple agents whose interaction graph is connected. Each agent in the network only knows one column of each of the adjacency matrices of the two graphs, and all agents collaboratively learn the graph matching by exchanging information with their neighbors. The proposed algorithm employs a projected primal-dual gradient method to handle equality constraints and a set constraint. Under the proposed algorithm, the agents' estimates of the permutation matrix converge to the optimal permutation globally and exponentially fast. Finally, simulation results are given to illustrate the effectiveness of the method., Comment: 10 pages, 2 figures, an extended version of a paper submitted to CDC20

In this paper, we develop a hybrid distance-angle rigidity theory that involves heterogeneous distances (or unsigned angles) and signed constraints for a framework in the 2-D and 3-D space. The new rigidity theory determines a (locally) unique formation shape up to a translation and a rotation by a set of distance and signed constraints, or up to a translation, a rotation and additionally a scaling factor by a set of unsigned angle and signed constraints. Under this new rigidity theory, we have a clue to resolve the flip (or reflection) and flex ambiguity for a target formation with hybrid distance-angle constraints. In particular, we can completely eliminate the ambiguity issues if formations are under a specific construction which is called \myemph{signed Henneberg construction} in this paper. We then apply the rigidity theory to formation shape control and develop a gradient-based control system that guarantees an exponential convergence close to a desired formation by inter-neighbor measurements. Several numerical simulations on formation shape control with hybrid distance-angle constraints are provided to validate the theoretical results.

This paper focuses on the equilibria and their regions of attraction for continuous-time nonlinear dynamical systems. The classical Poincar\'e--Hopf Theorem from differential topology is used to derive a result on a sufficient condition on the Jacobian for the existence of a unique equilibrium for the system, which is in fact locally exponentially stable. We show how to apply this result to the deterministic SIS networked model, and a nonlinear Lotka--Volterra system. We apply the result further to establish an impossibility conclusion for a class of distributed feedback controllers whose goal is to drive the SIS network to the zero equilibrium. Specifically, we show that if the uncontrolled system has a unique nonzero equilibrium (a diseased steady-state), then the controlled system also has a unique nonzero equilibrium. Applying results from monotone dynamical systems theory, we further show that both the uncontrolled and controlled system will converge to their nonzero equilibrium from all nonzero initial conditions exponentially fast. A counterpart sufficient condition for the existence of a unique equilibrium for a nonlinear discrete-time dynamical system is also presented., Comment: Submitted journal paper

In this paper, we study distributed methods for solving a Sylvester equation in the form of AX+XB=C for matrices A, B, C$\in R^{n\times n}$ with X being the unknown variable. The entries of A, B and C (called data) are partitioned into a number of pieces (or sometimes we permit these pieces to overlap). Then a network with a given structure is assigned, whose number of nodes is consistent with the partition. Each node has access to the corresponding set of data and holds a dynamic state. Nodes share their states among their neighbors defined from the network structure, and we aim to design flows that can asymptotically converge to a solution of this equation. The decentralized data partitions may be resulted directly from networks consisting of physically isolated subsystems, or indirectly from artificial and strategic design for processing large data sets. Natural partial row/column partitions, full row/column partitions and clustering block partitions of the data A, B and C are assisted by the use of the vectorized matrix equation. We show that the existing "consensus + projection" flow and the "local conservation + global consensus" flow for distributed linear algebraic equations can be used to drive distributed flows that solve this kind of equations. A "consensus + projection + symmetrization" flow is also developed for equations with symmetry constraints on the solution matrices. We reveal some fundamental convergence rate limitations for such flows regardless of the choices of node interaction strengths and network structures. For a special case with B=A$^T$, where the equation mentioned is reduced to a classical Lyapunov equation, we demonstrate that by exploiting the symmetry of data, we can obtain flows with lower complexity for certain partitions.