title: Stack graphs: Name resolution at scale

abstract: We present stack graphs, an extension of Visser et al.'s scope graphs framework. Stack graphs power Precise Code Navigation at GitHub, allowing users to navigate name binding references both within and across repositories. Like scope graphs, stack graphs encode the name binding information about a program in a graph structure, in which paths represent valid name bindings. Resolving a reference to its definition is then implemented with a simple path-finding search. GitHub hosts millions of repositories, containing petabytes of total code, implemented in hundreds of different programming languages, and receiving thousands of pushes per minute. To support this scale, we ensure that the graph construction and path-finding judgments are file-incremental: for each source file, we create an isolated subgraph without any knowledge of, or visibility into, any other file in the program. This lets us eliminate the storage and compute costs of reanalyzing file versions that we have already seen. Since most commits change a small fraction of the files in a repository, this greatly amortizes the operational costs of indexing large, frequently changed repositories over time. To handle type-directed name lookups (which require "pausing" the current lookup to resolve another name), our name resolution algorithm maintains a stack of the currently paused (but still pending) lookups. Stack graphs can be constructed via a purely syntactic analysis of the program's source code, using a new declarative graph construction language. This means that we can extract name binding information for every repository without any per-package configuration, and without having to invoke an arbitrary, untrusted, package-specific build process.


title: Design of cycloidal rays in optical waveguides in analogy to the fastest descent problem

abstract: In this work, we present the design of cycloidal waveguides from a gradient refractive index (GRIN) medium in analogy to the fastest descending problem in classical mechanics. Light rays propagate along cycloids in this medium, of which the refractive index can be determined from relating to the descending speed under gravity force. It can be used as GRIN lenses or waveguides, and the frequency specific focusing and imaging properties have been discussed. The results suggest that the waveguide can be viewed as an optical filter. Its frequency response characteristics change with the refractive index profile and the device geometries.


title: DEArt: Dataset of European Art

abstract: Large datasets that were made publicly available to the research community over the last 20 years have been a key enabling factor for the advances in deep learning algorithms for NLP or computer vision. These datasets are generally pairs of aligned image / manually annotated metadata, where images are photographs of everyday life. Scholarly and historical content, on the other hand, treat subjects that are not necessarily popular to a general audience, they may not always contain a large number of data points, and new data may be difficult or impossible to collect. Some exceptions do exist, for instance, scientific or health data, but this is not the case for cultural heritage (CH). The poor performance of the best models in computer vision - when tested over artworks - coupled with the lack of extensively annotated datasets for CH, and the fact that artwork images depict objects and actions not captured by photographs, indicate that a CH-specific dataset would be highly valuable for this community. We propose DEArt, at this point primarily an object detection and pose classification dataset meant to be a reference for paintings between the XIIth and the XVIIIth centuries. It contains more than 15000 images, about 80% non-iconic, aligned with manual annotations for the bounding boxes identifying all instances of 69 classes as well as 12 possible poses for boxes identifying human-like objects. Of these, more than 50 classes are CH-specific and thus do not appear in other datasets; these reflect imaginary beings, symbolic entities and other categories related to art. Additionally, existing datasets do not include pose annotations. Our results show that object detectors for the cultural heritage domain can achieve a level of precision comparable to state-of-art models for generic images via transfer learning.


title: Conformalized survival analysis with adaptive cutoffs

abstract: This paper introduces an assumption-lean method that constructs valid and efficient lower predictive bounds (LPBs) for survival times with censored data. We build on recent work by Cand\`es et al. (2021), whose approach first subsets the data to discard any data points with early censoring times, and then uses a reweighting technique (namely, weighted conformal inference (Tibshirani et al., 2019)) to correct for the distribution shift introduced by this subsetting procedure. For our new method, instead of constraining to a fixed threshold for the censoring time when subsetting the data, we allow for a covariate-dependent and data-adaptive subsetting step, which is better able to capture the heterogeneity of the censoring mechanism. As a result, our method can lead to LPBs that are less conservative and give more accurate information. We show that in the Type I right-censoring setting, if either of the censoring mechanism or the conditional quantile of survival time is well estimated, our proposed procedure achieves nearly exact marginal coverage, where in the latter case we additionally have approximate conditional coverage. We evaluate the validity and efficiency of our proposed algorithm in numerical experiments, illustrating its advantage when compared with other competing methods. Finally, our method is applied to a real dataset to generate LPBs for users' active times on a mobile app.


title: Global total precipitable water variations and trends during 1958-2021

abstract: Global responses of the hydrological cycle to climate change have been widely studied but uncertainties of temperature responses to lower-tropospheric water vapor still remain. Here, we investigate the trends in global total precipitable water (TPW) and surface temperature from 1958 to 2021 using improved ERA5 and JRA-55 reanalysis datasets and further validate these trends by using radiosonde, Atmospheric Infrared Sounder (AIRS), and Microwave Satellite (SSMI(S)) observations. Our results indicate a global increase in total precipitable water (TPW) of 0.66% per decade according to ERA5 data and 0.88 % per decade in JRA-55 data.These variations in TPW reflect the interactions of global warming feedback mechanisms across different spatial scales. Our results also revealed a significant near-surface temperature (T2m) warming trend at the rate of 0.14 K dec-1 and a strong water vapor response to temperature at a rate of 4-6 % K -1 globally, with land areas warming approximately twice as fast as the oceans. The relationship between TPW and T2m or surface skin temperature Ts showed a variation around 6 - 8 % K -1 in the 15-60 latitude band, aligning with theoretical estimates from the Clausius Clapeyron equation.


title: Fast convergent PML method for scattering with periodic surfaces: the exceptional case

abstract: In the author's previous paper (Zhang et al. 2022), exponential convergence was proved for the perfectly matched layers (PML) approximation of scattering problems with periodic surfaces in 2D. However, due to the overlapping of singularities, an exceptional case, i.e., when the wave number is a half integer, has to be excluded in the proof. However, numerical results for these cases still have fast convergence rate and this motivates us to go deeper into these cases. In this paper, we focus on these cases and prove that the fast convergence result for the discretized form. Numerical examples are also presented to support our theoretical results.


title: Pump--probe x-ray microscopy of photo-induced magnetization dynamics at MHz repetition rates

abstract: We present time-resolved scanning x-ray microscopy measurements with picosecond photo-excitation via a tailored infrared pump laser at a scanning transmission x-ray microscope. Specifically, we image the laser-induced demagnetization and remagnetization of thin ferrimagnetic GdFe films proceeding on a few nanoseconds time scale. Controlling the heat load on the sample via additional reflector and heat-sink layers allows us to conduct destruction-free measurements at a repetition rate of \SI{50}{\mega\hertz}. Near-field enhancement of the photo-excitation and controlled annealing effects lead to laterally heterogeneous magnetization dynamics which we trace with \SI{30}{\nano\meter} spatial resolution. Our work opens new opportunities to study photo-induced dynamics on the nanometer scale, with access to picosecond to nanosecond timescales, which is of technological relevance, especially in the field of magnetism.


title: Interval Markov Decision Processes with Continuous Action-Spaces

abstract: Interval Markov Decision Processes (IMDPs) are finite-state uncertain Markov models, where the transition probabilities belong to intervals. Recently, there has been a surge of research on employing IMDPs as abstractions of stochastic systems for control synthesis. However, due to the absence of algorithms for synthesis over IMDPs with continuous action-spaces, the action-space is assumed discrete a-priori, which is a restrictive assumption for many applications. Motivated by this, we introduce continuous-action IMDPs (caIMDPs), where the bounds on transition probabilities are functions of the action variables, and study value iteration for maximizing expected cumulative rewards. Specifically, we decompose the max-min problem associated to value iteration to $|\mathcal{Q}|$ max problems, where $|\mathcal{Q}|$ is the number of states of the caIMDP. Then, exploiting the simple form of these max problems, we identify cases where value iteration over caIMDPs can be solved efficiently (e.g., with linear or convex programming). We also gain other interesting insights: e.g., in certain cases where the action set $\mathcal{A}$ is a polytope, synthesis over a discrete-action IMDP, where the actions are the vertices of $\mathcal{A}$, is sufficient for optimality. We demonstrate our results on a numerical example. Finally, we include a short discussion on employing caIMDPs as abstractions for control synthesis.


title: The Luminosity Phase Space of Galactic and Extragalactic X-ray Transients Out to Intermediate Redshifts

abstract: We present a detailed compilation and analysis of the X-ray phase space of low- to intermediate-redshift ($ 0\le z \le 1$) transients that consolidates observed light curves (and theory where necessary) for a large variety of classes of transient/variable phenomena in the 0.3--10 keV energy band. We include gamma-ray burst afterglows, supernovae, supernova shock breakouts and shocks interacting with the environment, tidal disruption events and active galactic nuclei, fast blue optical transients, cataclysmic variables, magnetar flares/outbursts and fast radio bursts, cool stellar flares, X-ray binary outbursts, and ultraluminous X-ray sources. Our overarching goal is to offer a comprehensive resource for the examination of these ephemeral events, extending the X-ray duration-luminosity phase space (DLPS) to show luminosity evolution. We use existing observations (both targeted and serendipitous) to characterize the behavior of various transient/variable populations. Contextualizing transient signals in the larger DLPS serves two primary purposes: to identify areas of interest (i.e., regions in the parameter space where one would expect detections, but in which observations have historically been lacking) and to provide initial qualitative guidance in classifying newly discovered transient signals. We find that while the most luminous (largely extragalactic) and least luminous (largely Galactic) part of the phase space is well-populated at $t > 0.1$ days, intermediate luminosity phenomena (L$_x = 10^{34} - 10^{42}$ erg s$^{-1}$) represent a gap in the phase space. We thus identify L$_x = 10^{34} - 10^{42}$ erg s$^{-1}$ and $t = 10^{-4} - 0.1$ days as a key discovery phase space in transient X-ray astronomy.


title: Attention-based Neural Cellular Automata

abstract: Recent extensions of Cellular Automata (CA) have incorporated key ideas from modern deep learning, dramatically extending their capabilities and catalyzing a new family of Neural Cellular Automata (NCA) techniques. Inspired by Transformer-based architectures, our work presents a new class of $\textit{attention-based}$ NCAs formed using a spatially localized$\unicode{x2014}$yet globally organized$\unicode{x2014}$self-attention scheme. We introduce an instance of this class named $\textit{Vision Transformer Cellular Automata}$ (ViTCA). We present quantitative and qualitative results on denoising autoencoding across six benchmark datasets, comparing ViTCA to a U-Net, a U-Net-based CA baseline (UNetCA), and a Vision Transformer (ViT). When comparing across architectures configured to similar parameter complexity, ViTCA architectures yield superior performance across all benchmarks and for nearly every evaluation metric. We present an ablation study on various architectural configurations of ViTCA, an analysis of its effect on cell states, and an investigation on its inductive biases. Finally, we examine its learned representations via linear probes on its converged cell state hidden representations, yielding, on average, superior results when compared to our U-Net, ViT, and UNetCA baselines.


title: Uncertainty-Aware DNN for Multi-Modal Camera Localization

abstract: Camera localization, i.e., camera pose regression, represents an important task in computer vision since it has many practical applications such as in the context of intelligent vehicles and their localization. Having reliable estimates of the regression uncertainties is also important, as it would allow us to catch dangerous localization failures. In the literature, uncertainty estimation in Deep Neural Networks (DNNs) is often performed through sampling methods, such as Monte Carlo Dropout (MCD) and Deep Ensemble (DE), at the expense of undesirable execution time or an increase in hardware resources. In this work, we considered an uncertainty estimation approach named Deep Evidential Regression (DER) that avoids any sampling technique, providing direct uncertainty estimates. Our goal is to provide a systematic approach to intercept localization failures of camera localization systems based on DNNs architectures, by analyzing the generated uncertainties. We propose to exploit CMRNet, a DNN approach for multi-modal image to LiDAR map registration, by modifying its internal configuration to allow for extensive experimental activity on the KITTI dataset. The experimental section highlights CMRNet's major flaws and proves that our proposal does not compromise the original localization performances but also provides, at the same time, the necessary introspection measures that would allow end-users to act accordingly.


title: Capacitary density and removable sets for Newton-Sobolev functions in metric spaces

abstract: In a complete metric space equipped with a doubling measure and supporting a $(1,1)$-Poincar\'e inequality, we show that every set satisfying a suitable capacitary density condition is removable for Newton-Sobolev functions.


title: Isometric Representations in Neural Networks Improve Robustness

abstract: Artificial and biological agents cannon learn given completely random and unstructured data. The structure of data is encoded in the metric relationships between data points. In the context of neural networks, neuronal activity within a layer forms a representation reflecting the transformation that the layer implements on its inputs. In order to utilize the structure in the data in a truthful manner, such representations should reflect the input distances and thus be continuous and isometric. Supporting this statement, recent findings in neuroscience propose that generalization and robustness are tied to neural representations being continuously differentiable. In machine learning, most algorithms lack robustness and are generally thought to rely on aspects of the data that differ from those that humans use, as is commonly seen in adversarial attacks. During cross-entropy classification, the metric and structural properties of network representations are usually broken both between and within classes. This side effect from training can lead to instabilities under perturbations near locations where such structure is not preserved. One of the standard solutions to obtain robustness is to add ad hoc regularization terms, but to our knowledge, forcing representations to preserve the metric structure of the input data as a stabilising mechanism has not yet been studied. In this work, we train neural networks to perform classification while simultaneously maintaining within-class metric structure, leading to isometric within-class representations. Such network representations turn out to be beneficial for accurate and robust inference. By stacking layers with this property we create a network architecture that facilitates hierarchical manipulation of internal neural representations. Finally, we verify that isometric regularization improves the robustness to adversarial attacks on MNIST.


title: On symmetric 2-(70,24,8) designs with an automorphism of order 6

abstract: In this paper we analyze possible actions of an automorphism of order six on a $2$-$(70, 24, 8)$ design, and give a complete classification for the action of the cyclic automorphism group of order six $G= \langle \rho \rangle \cong Z_6 \cong Z_2 \times Z_3$ where $\rho^3$ fixes exactly $14$ points (blocks) and $\rho^2$ fixes $4$ points (blocks). Up to isomorphism, there are $3718$ such designs. This result significantly increases the number of known $2$-$(70,24,8)$ designs.


title: Late lumping of transformation-based feedback laws for boundary control systems

abstract: Late-lumping feedback design for infinite-dimensional linear systems with unbounded input operators is considered. The proposed scheme is suitable for the approximation of backstepping and flatness-based designs and relies on a decomposition of the feedback into a bounded and an unbounded part. Approximation applies to the bounded part only, while the unbounded part is assumed to allow for an exact realization. Based on spectral results, the convergence of the closed-loop dynamics to the desired dynamics is established. By duality, similar results apply to the approximation of the observer output-injection gains for systems with boundary observation. The proposed design and approximation steps are demonstrated and illustrated based on a hyperbolic infinite-dimensional system.


title: Gas clumping and its effect on hydrostatic bias in the MACSIS simulations

abstract: We use the MACSIS hydrodynamical simulations to estimate the extent of gas clumping in the intracluster medium of massive galaxy clusters and how it affects the hydrostatic mass bias. By comparing the clumping to the azimuthal scatter in the emission measure, an observational proxy, we find that they both increase with radius and are larger in higher-mass and dynamically perturbed systems. Similar trends are also seen for the azimuthal temperature scatter and non-thermal pressure fraction, both of which correlate with density fluctuations, with these values also increasing with redshift. However, in agreement with recent work, we find only a weak correlation between the clumping, or its proxies, and the hydrostatic mass bias. To reduce the effect of clumping in the projected profiles, we compute the azimuthal median following recent observational studies, and find this reduces the scatter in the bias. We also attempt to correct the cluster masses by using a non-thermal pressure term and find over-corrected mass estimates ($1-b=0.86$ to $1-b=1.15$) from 3D gas profiles but improved mass estimates ($1-b=0.75$ to $1-b=0.85$) from projected gas profiles, with the caveat of systematically increased scatter. We conclude that the cluster-averaged mass bias is minimised from applying a non-thermal pressure correction ($1-b=0.85$) with more modest reductions from selecting clusters that have low clumping ($1-b=0.79$) or are dynamically relaxed ($1-b=0.80$). However, the latter selection is most effective at minimising the scatter for individual objects. Such results can be tested with next generation X-ray missions equipped with high-resolution spectrometers such as Athena.


title: On The Equivalence Of The Mean Variance Criterion And Stochastic Dominance Criteria

abstract: We study the necessary and sufficient conditions under which the Mean-Variance Criterion (MVC) is equivalent to the Maximum Expected Utility Criterion (MEUC), for two lotteries. Based on Chamberlain (1983), we conclude that the MVC is equivalent to the Second-order Stochastic Dominance Rule (SSDR) under any symmetric Elliptical distribution. We then discuss the work of Schuhmacher et al. (2021). Although their theoretical findings deduce that the Mean-Variance Analysis remains valid under Skew-Elliptical distributions, we argue that this does not entail that the MVC coincides with the SSDR. In fact, generating multiple MV-pairs that follow a Skew-Normal distribution it becomes evident that the MVC fails to coincide with the SSDR for some types of risk-averse investors. In the second part of this work, we examine the premise of Levy and Markowitz (1979) that "the MVC deduces the maximization of the expected utility of an investor, under any approximately quadratic utility function, without making any further assumption on the distribution of the lotteries". Using Monte Carlo Simulations, we find out that the set of approximately quadratic utility functions is too narrow. Specifically, our simulations indicate that $\log{(a+Z)}$ and $(1+Z)^a$ are almost quadratic, while $-e^{-a(1+Z)}$ and $-(1+Z)^{-a}$ fail to approximate a quadratic utility function under either an Extreme Value or a Stable Pareto distribution.


title: WiserVR: Semantic Communication Enabled Wireless Virtual Reality Delivery

abstract: Virtual reality (VR) over wireless is expected to be one of the killer applications in next-generation communication networks. Nevertheless, the huge data volume along with stringent requirements on latency and reliability under limited bandwidth resources makes untethered wireless VR delivery increasingly challenging. Such bottlenecks, therefore, motivate this work to seek the potential of using semantic communication, a new paradigm that promises to significantly ease the resource pressure, for efficient VR delivery. To this end, we propose a novel framework, namely WIreless SEmantic deliveRy for VR (WiserVR), for delivering consecutive 360{\deg} video frames to VR users. Specifically, deep learning-based multiple modules are well-devised for the transceiver in WiserVR to realize high-performance feature extraction and semantic recovery. Among them, we dedicatedly develop a concept of semantic location graph and leverage the joint-semantic-channel-coding method with knowledge sharing to not only substantially reduce communication latency, but also to guarantee adequate transmission reliability and resilience under various channel states. Moreover, implementation of WiserVR is presented, followed by corresponding initial simulations for performance evaluation compared with benchmarks. Finally, we discuss several open issues and offer feasible solutions to unlock the full potential of WiserVR.


title: The dissolution of a miscible drop rising or falling in another liquid at low Reynolds number

abstract: 'A basic and basically unsolved problem in fluid dynamics is to determine the evolution of rising bubbles and falling drops of one miscible liquid in another' [1]. Here, we address this important literature gap and present the first theory predicting the velocity, volume and composition of such drops at low Reynolds numbers. For the case where the diffusion out of the drop is negligible, we obtain a universal scaling law. For the more general case where diffusion occurs into and out of the drop, the full dynamics is governed by a parameter-free first-order ordinary differential equation, whose closed form solution exists, and only depends on the initial condition. Our analysis depends primarily on 'drop-scale' effective parameters for the diffusivity through the interfacial boundary layer. We validate our results against experimental data for water drops suspended in syrup, corresponding to certain regimes of the mass exchange ratio between water and syrup, and by this explicitly identify the drop-scale parameters of the theory. [1] Joseph, D.D. and Renardy, Y.Y., 2013. Fundamentals of two-fluid dynamics: part II: lubricated transport, drops and miscible liquids (Vol. 4). Springer Science & Business Media.


title: Impact of trans-Planckian excitations on black-hole radiation in dipolar condensates

abstract: We consider a quasi-one-dimensional dipolar condensate in an analogue black hole setup. It is shown that the existence of a roton minimum in the condensate dispersion relation leaves deep imprints onto the Hawking radiation spectrum. In particular, the emitted radiation can be either more intense or suppressed, depending on the depth of the roton minimum in the excitation spectrum. In addition, we find that spontaneous particle creation occurs even when the horizon is removed. Our results establish that dipolar condensates offer a richer and more versatile environment for the simulation of particle production from the quantum vacuum in the presence of horizon-interfaces than their contact-interaction counterparts.


title: EquiMod: An Equivariance Module to Improve Self-Supervised Learning

abstract: Self-supervised visual representation methods are closing the gap with supervised learning performance. These methods rely on maximizing the similarity between embeddings of related synthetic inputs created through data augmentations. This can be seen as a task that encourages embeddings to leave out factors modified by these augmentations, i.e. to be invariant to them. However, this only considers one side of the trade-off in the choice of the augmentations: they need to strongly modify the images to avoid simple solution shortcut learning (e.g. using only color histograms), but on the other hand, augmentations-related information may be lacking in the representations for some downstream tasks (e.g. color is important for birds and flower classification). Few recent works proposed to mitigate the problem of using only an invariance task by exploring some form of equivariance to augmentations. This has been performed by learning additional embeddings space(s), where some augmentation(s) cause embeddings to differ, yet in a non-controlled way. In this work, we introduce EquiMod a generic equivariance module that structures the learned latent space, in the sense that our module learns to predict the displacement in the embedding space caused by the augmentations. We show that applying that module to state-of-the-art invariance models, such as SimCLR and BYOL, increases the performances on CIFAR10 and ImageNet datasets. Moreover, while our model could collapse to a trivial equivariance, i.e. invariance, we observe that it instead automatically learns to keep some augmentations-related information beneficial to the representations.


title: An efficient algorithm for the $\ell_{p}$ norm based metric nearness problem

abstract: Given a dissimilarity matrix, the metric nearness problem is to find the nearest matrix of distances that satisfy the triangle inequalities. This problem has wide applications, such as sensor networks, image processing, and so on. But it is of great challenge even to obtain a moderately accurate solution due to the $O(n^{3})$ metric constraints and the nonsmooth objective function which is usually a weighted $\ell_{p}$ norm based distance. In this paper, we propose a delayed constraint generation method with each subproblem solved by the semismooth Newton based proximal augmented Lagrangian method (PALM) for the metric nearness problem. Due to the high memory requirement for the storage of the matrix related to the metric constraints, we take advantage of the special structure of the matrix and do not need to store the corresponding constraint matrix. A pleasing aspect of our algorithm is that we can solve these problems involving up to $10^{8}$ variables and $10^{13}$ constraints. Numerical experiments demonstrate the efficiency of our algorithm. In theory, firstly, under a mild condition, we establish a primal-dual error bound condition which is very essential for the analysis of local convergence rate of PALM. Secondly, we prove the equivalence between the dual nondegeneracy condition and nonsingularity of the generalized Jacobian for the inner subproblem of PALM. Thirdly, when $q(\cdot)=\|\cdot\|_{1}$ or $\|\cdot\|_{\infty}$, without the strict complementarity condition, we also prove the equivalence between the the dual nondegeneracy condition and the uniqueness of the primal solution.

categories: math.OC cs.LG cs.NA math.NA stat.CO stat.ML


title: data2vec-aqc: Search for the right Teaching Assistant in the Teacher-Student training setup

abstract: In this paper, we propose a new Self-Supervised Learning (SSL) algorithm called data2vec-aqc, for speech representation learning from unlabeled speech data. Our goal is to improve SSL for speech in domains where both unlabeled and labeled data are limited. Building on the recently introduced data2vec, we introduce additional modules to the data2vec framework that leverage the benefit of data augmentations, quantized representations, and clustering. The interaction between these modules helps solve the cross-contrastive loss as an additional self-supervised objective. data2vec-aqc achieves up to 14.1% and 20.9% relative WER improvement over the existing state-of-the-art data2vec system over the test-clean and test-other sets, respectively of LibriSpeech, without the use of any language model (LM). Our proposed model also achieves up to 17.8\% relative WER gains over the baseline data2vec when fine-tuned on a subset of the Switchboard dataset. Code: https://github.com/Speech-Lab-IITM/data2vec-aqc.


title: Superposition Formulae for the Geometric B\"acklund Transformations of the Hyperbolic and Elliptic Sine-Gordon and Sinh-Gordon Equations

abstract: We provide superposition formulae for the six cases of B\"acklund transformations corresponding to space-like and time-like surfaces in the 3-dimensional pseudo-Euclidean space. In each case, the surfaces have constant negative or positive Gaussian curvature and they correspond to solutions of one of the following equations: the sine-Gordon, the sinh-Gordon, the elliptic sine-Gordon and the elliptic sinh-Gordon equation. The superposition formulae provide infinitely many solutions algebraically after the first integration of the B\"acklund transformation. Such transformations and the corresponding superposition formulae provide solutions of the same hyperbolic equation, while they show an unusual property for the elliptic equations. The B\"acklund transformation alternates solutions of the elliptic sinh-Gordon equation with those of the elliptic sine-Gordon equation and the superposition formulae provide solutions of the same elliptic equation. Explicit examples and illustrations are given.


title: Exact Completeness of LP Hierarchies for Linear Codes

abstract: Determining the maximum size $A_2(n,d)$ of a binary code of blocklength $n$ and distance $d$ remains an elusive open question even when restricted to the important class of linear codes. Recently, two linear programming hierarchies extending Delsarte's LP were independently proposed to upper bound $A_2^{\text{Lin}}(n,d)$ (the analogue of $A_2(n,d)$ for linear codes). One of these hierarchies, by the authors, was shown to be approximately complete in the sense that the hierarchy converges to $A_2^{\text{Lin}}(n,d)$ as the level grows beyond $n^2$. Despite some structural similarities, not even approximate completeness was known for the other hierarchy by Loyfer and Linial. In this work, we prove that both hierarchies recover the exact value of $A_2^{\text{Lin}}(n,d)$ at level $n$. We also prove that at this level the polytope of Loyfer and Linial is integral.Even though these hierarchies seem less powerful than general hierarchies such as Sum-of-Squares, we show that they have enough structure to yield exact completeness via pseudoprobabilities.


title: How Democracies Polarize: A Multilevel Perspective

abstract: Democracies employ elections at various scales to select officials at the corresponding levels of administration. The geographical distribution of political opinion, the policy issues delegated to each level, and the multilevel interactions between elections can all greatly impact the makeup of these representative bodies. This perspective is not new: the adoption of federal systems has been motivated by the idea that they possess desirable traits not provided by democracies on a single scale. Yet most existing models of polarization do not capture how nested local and national elections interact with heterogeneous political geographies. We begin by developing a framework to describe the multilevel distribution of opinions and analyze the flow of variance among geographic scales, applying it to historical data in the United States from 1912 to 2020. We describe how unstable elections can arise due to the spatial distribution of opinions and how tradeoffs occur between national and local elections. We also examine multi-dimensional spaces of political opinion, for which we show that a decrease in local salience can constrain the dimensions along which elections occur, preventing a federal system from serving as an effective safeguard against polarization. These analyses, based on the interactions between elections and opinion distributions at various scales, offer insights into how democracies can be strengthened to mitigate polarization and increase electoral representation.


title: Phase space geometry and optimal state preparation in quantum metrology with collective spins

abstract: We revisit well-known protocols in quantum metrology using collective spins and propose a unifying picture for optimal state preparation based on a semiclassical description in phase space. We show how this framework allows for quantitative predictions of the timescales required to prepare various metrologically useful states, and that these predictions remain accurate even for moderate system sizes, surprisingly far from the classical limit. Furthermore, this framework allows us to build a geometric picture that relates optimal (exponentially fast) entangled probe preparation to the existence of separatrices connecting saddle points in phase space. We illustrate our results with the paradigmatic examples of the two-axis counter-twisting and twisting-and-turning Hamiltonians, where we provide analytical expressions for all the relevant optimal time scales. Finally, we propose a generalization of these models to include $p$-body collective interaction (or $p$-order twisting), beyond the usual case of $p=2$. Using our geometric framework, we prove a no-go theorem for the local optimality of these models for $p>2$.


title: Spin-orbit coupling-enhanced valley ordering of malleable bands in twisted bilayer graphene on WSe2

abstract: New phases of matter can be stabilized by a combination of diverging electronic density of states, strong interactions, and spin-orbit coupling. Recent experiments in magic-angle twisted bilayer graphene (TBG) have uncovered a wealth of novel phases as a result of interaction-driven spin-valley flavour polarization. In this work, we explore correlated phases appearing due to the combined effect of spin-orbit coupling-enhanced valley polarization and large density of states below half filling ($\nu \lesssim 2$) of the moir\'e band in a TBG coupled to tungsten diselenide. We observe anomalous Hall effect, accompanied by a series of Lifshitz transitions, that are highly tunable with carrier density and magnetic field. Strikingly, the magnetization shows an abrupt sign change in the vicinity of half-filling, confirming its orbital nature. The coercive fields reported are about an order of magnitude higher than previous studies in graphene-based moir\'e systems, presumably aided by a Stoner instability favoured by the van Hove singularities in the malleable bands. While the Hall resistance is not quantized at zero magnetic fields, indicative of a ground state with partial valley polarization, perfect quantization and complete valley polarization are observed at finite fields. Our findings illustrate that singularities in the flat bands in the presence of spin-orbit coupling can stabilize ordered phases even at non-integer moir\'e band fillings.


title: Quantum algorithms for classical Boolean functions via adaptive measurements: Exponential reductions in space-time resources

abstract: The limited computational power of constant-depth quantum circuits can be boosted by adapting future gates according to the outcomes of mid-circuit measurements. We formulate computation of a variety of Boolean functions in the framework of adaptive measurement-based quantum computation using a cluster state resource and a classical side-processor that can add bits modulo 2, so-called $l2$-MBQC. Our adaptive approach overcomes a known challenge that computing these functions in the nonadaptive setting requires a resource state that is exponentially large in the size of the computational input. In particular, we construct adaptive $l2$-MBQC algorithms based on the quantum signal processing technique that compute the mod-$p$ functions with the best known scaling in the space-time resources (i.e., qubit count, quantum circuit depth, classical memory size, and number of calls to the side-processor). As the subject is diverse and has a long history, the paper includes reviews of several previously constructed algorithms and recasts them as adaptive $l2$-MBQCs using cluster state resources. Our results constitute an alternative proof of an old theorem regarding an oracular separation between the power of constant-depth quantum circuits and constant-depth classical circuits with unbounded fan-in NAND and mod-$p$ gates for any prime $p$.


title: Fair Visual Recognition via Intervention with Proxy Features

abstract: Deep learning models often learn to make predictions that rely on sensitive social attributes like gender and race, which poses significant fairness risks, especially in societal applications, e.g., hiring, banking, and criminal justice. Existing work tackles this issue by minimizing information about social attributes in models for debiasing. However, the high correlation between target task and social attributes makes bias mitigation incompatible with target task accuracy. Recalling that model bias arises because the learning of features in regard to bias attributes (i.e., bias features) helps target task optimization, we explore the following research question: \emph{Can we leverage proxy features to replace the role of bias feature in target task optimization for debiasing?} To this end, we propose \emph{Proxy Debiasing}, to first transfer the target task's learning of bias information from bias features to artificial proxy features, and then employ causal intervention to eliminate proxy features in inference. The key idea of \emph{Proxy Debiasing} is to design controllable proxy features to on one hand replace bias features in contributing to target task during the training stage, and on the other hand easily to be removed by intervention during the inference stage. This guarantees the elimination of bias features without affecting the target information, thus addressing the fairness-accuracy paradox in previous debiasing solutions. We apply \emph{Proxy Debiasing} to several benchmark datasets, and achieve significant improvements over the state-of-the-art debiasing methods in both of accuracy and fairness.


title: CircleSnake: Instance Segmentation with Circle Representation

abstract: Circle representation has recently been introduced as a medical imaging optimized representation for more effective instance object detection on ball-shaped medical objects. With its superior performance on instance detection, it is appealing to extend the circle representation to instance medical object segmentation. In this work, we propose CircleSnake, a simple end-to-end circle contour deformation-based segmentation method for ball-shaped medical objects. Compared to the prevalent DeepSnake method, our contribution is three-fold: (1) We replace the complicated bounding box to octagon contour transformation with a computation-free and consistent bounding circle to circle contour adaption for segmenting ball-shaped medical objects; (2) Circle representation has fewer degrees of freedom (DoF=2) as compared with the octagon representation (DoF=8), thus yielding a more robust segmentation performance and better rotation consistency; (3) To the best of our knowledge, the proposed CircleSnake method is the first end-to-end circle representation deep segmentation pipeline method with consistent circle detection, circle contour proposal, and circular convolution. The key innovation is to integrate the circular graph convolution with circle detection into an end-to-end instance segmentation framework, enabled by the proposed simple and consistent circle contour representation. Glomeruli are used to evaluate the performance of the benchmarks. From the results, CircleSnake increases the average precision of glomerular detection from 0.559 to 0.614. The Dice score increased from 0.804 to 0.849. The code has been released: https://github.com/hrlblab/CircleSnake


title: Task-Oriented Over-the-Air Computation for Multi-Device Edge AI

abstract: Departing from the classic paradigm of data-centric designs, the 6G networks for supporting edge AI features task-oriented techniques that focus on effective and efficient execution of AI task. Targeting end-to-end system performance, such techniques are sophisticated as they aim to seamlessly integrate sensing (data acquisition), communication (data transmission), and computation (data processing). Aligned with the paradigm shift, a task-oriented over-the-air computation (AirComp) scheme is proposed in this paper for multi-device split-inference system. In the considered system, local feature vectors, which are extracted from the real-time noisy sensory data on devices, are aggregated over-the-air by exploiting the waveform superposition in a multiuser channel. Then the aggregated features as received at a server are fed into an inference model with the result used for decision making or control of actuators. To design inference-oriented AirComp, the transmit precoders at edge devices and receive beamforming at edge server are jointly optimized to rein in the aggregation error and maximize the inference accuracy. The problem is made tractable by measuring the inference accuracy using a surrogate metric called discriminant gain, which measures the discernibility of two object classes in the application of object/event classification. It is discovered that the conventional AirComp beamforming design for minimizing the mean square error in generic AirComp with respect to the noiseless case may not lead to the optimal classification accuracy. The reason is due to the overlooking of the fact that feature dimensions have different sensitivity towards aggregation errors and are thus of different importance levels for classification. This issue is addressed in this work via a new task-oriented AirComp scheme designed by directly maximizing the derived discriminant gain.


title: An Aggregation of Aggregation Methods in Computational Pathology

abstract: Image analysis and machine learning algorithms operating on multi-gigapixel whole-slide images (WSIs) often process a large number of tiles (sub-images) and require aggregating predictions from the tiles in order to predict WSI-level labels. In this paper, we present a review of existing literature on various types of aggregation methods with a view to help guide future research in the area of computational pathology (CPath). We propose a general CPath workflow with three pathways that consider multiple levels and types of data and the nature of computation to analyse WSIs for predictive modelling. We categorize aggregation methods according to the context and representation of the data, features of computational modules and CPath use cases. We compare and contrast different methods based on the principle of multiple instance learning, perhaps the most commonly used aggregation method, covering a wide range of CPath literature. To provide a fair comparison, we consider a specific WSI-level prediction task and compare various aggregation methods for that task. Finally, we conclude with a list of objectives and desirable attributes of aggregation methods in general, pros and cons of the various approaches, some recommendations and possible future directions.


title: Improving Performance of Higher-Order Codirectional Raman Amplifiers Using Phase-Modulated Signals -- Functional Principle

abstract: Existing unrepeatered submarine links are increasingly upgraded to the most advanced modulation format currently available for commercial applications. Quite often the use of third-order codirectional Raman amplifiers is necessary. Power fluctuations of the involved high-power pump induces phase shifts in phase modulated signals via the nonlinear Kerr effects. An approach for reducing the impact of this effect on system performance is explained in this document.


title: Instance-Dependent Generalization Bounds via Optimal Transport

abstract: Existing generalization bounds fail to explain crucial factors that drive the generalization of modern neural networks. Since such bounds often hold uniformly over all parameters, they suffer from over-parametrization and fail to account for the strong inductive bias of initialization and stochastic gradient descent. As an alternative, we propose a novel optimal transport interpretation of the generalization problem. This allows us to derive instance-dependent generalization bounds that depend on the local Lipschitz regularity of the learned prediction function in the data space. Therefore, our bounds are agnostic to the parametrization of the model and work well when the number of training samples is much smaller than the number of parameters. With small modifications, our approach yields accelerated rates for data on low-dimensional manifolds and guarantees under distribution shifts. We empirically analyze our generalization bounds for neural networks, showing that the bound values are meaningful and capture the effect of popular regularization methods during training.


title: Revealing microcanonical phases and phase transitions of strongly correlated electrons via time-averaged classical shadows

abstract: Quantum computers and simulators promise to enable the study of strongly correlated quantum systems. Yet, surprisingly, it is hard for them to compute ground states. They can, however, efficiently compute the dynamics of closed quantum systems. We propose a method to study the quantum thermodynamics of strongly correlated electrons from quantum dynamics. We define time-averaged classical shadows (TACS) and prove it is a classical shadow(CS) of the von Neumann ensemble, the time-averaged density matrix. We then show that the diffusion maps, an unsupervised machine learning algorithm, can efficiently learn the phase diagram and phase transition of the one-dimensional transverse field Ising model both for ground states using CS \emph{and state trajectories} using TACS. It does so from state trajectories by learning features that appear to be susceptibility and entropy from a total of 90,000 shots taken along a path in the microcanonical phase diagram. Our results suggest a low number of shots from quantum simulators can produce quantum thermodynamic data with a quantum advantage.


title: Executable Models and Instance Tracking for Decentralized Applications -- Towards an Architecture Based on Blockchains and Cloud Platforms

abstract: The execution of decentralized applications on blockchains is limited today by technical and organizational barriers, including scalability and the high complexity to specify execution correctly for developers as well as for domain experts in organizations. Overcoming these limitations could allow for decentralized coordination beyond data, where distributed parties rely on higher-level abstractions for coordinating their actions using decentralized applications, not limited to organizations. Towards this goal, the paper at hand proposes executable models as high-level abstraction that can be observed and tracked by distributed parties. In particular, it is investigated how executable models on cloud platforms can be coupled with smart contracts for tracking their execution, concluding with an architecture as exploratory research result towards supporting scalability and decentralized coordination.


title: Where Do We Go From Here? Guidelines For Offline Recommender Evaluation

abstract: Various studies in recent years have pointed out large issues in the offline evaluation of recommender systems, making it difficult to assess whether true progress has been made. However, there has been little research into what set of practices should serve as a starting point during experimentation. In this paper, we examine four larger issues in recommender system research regarding uncertainty estimation, generalization, hyperparameter optimization and dataset pre-processing in more detail to arrive at a set of guidelines. We present a TrainRec, a lightweight and flexible toolkit for offline training and evaluation of recommender systems that implements these guidelines. Different from other frameworks, TrainRec is a toolkit that focuses on experimentation alone, offering flexible modules that can be can be used together or in isolation. Finally, we demonstrate TrainRec's usefulness by evaluating a diverse set of twelve baselines across ten datasets. Our results show that (i) many results on smaller datasets are likely not statistically significant, (ii) there are at least three baselines that perform well on most datasets and should be considered in future experiments, and (iii) improved uncertainty quantification (via nested CV and statistical testing) rules out some reported differences between linear and neural methods. Given these results, we advocate that future research should standardize evaluation using our suggested guidelines.


title: Thermophotovoltaic conversion efficiency measurement at high view factors

abstract: A standardized method for measuring thermophotovoltaic (TPV) efficiency has not been yet established, which makes the reported results difficult to compare. Besides, most of the TPV efficiencies reported to date have been obtained using small view factors, i.e., large cell-to-emitter distances, so the impact of the series resistance is usually underestimated, and the optical cavity effects, i.e., the multiple reflections taking place between the emitter and the cell, are not accounted for experimentally. In this work, we present an experimental setup able to measure the TPV efficiency under high view factors (up to 0.98), by using small emitter-to-cell distances (< 1 mm). This allows a more accurate direct measurement of the TPV efficiency at higher power densities than previous works. As a result, a TPV efficiency of 26.4+-0.1 % and a power density of 4.3+-0.8 W/cm2 have been obtained for an InGaAs TPV cell with a back surface reflector irradiated by a graphite thermal emitter at 1592 C.


title: A Quantum Kernel Learning Approach to Acoustic Modeling for Spoken Command Recognition

abstract: We propose a quantum kernel learning (QKL) framework to address the inherent data sparsity issues often encountered in training large-scare acoustic models in low-resource scenarios. We project acoustic features based on classical-to-quantum feature encoding. Different from existing quantum convolution techniques, we utilize QKL with features in the quantum space to design kernel-based classifiers. Experimental results on challenging spoken command recognition tasks for a few low-resource languages, such as Arabic, Georgian, Chuvash, and Lithuanian, show that the proposed QKL-based hybrid approach attains good improvements over existing classical and quantum solutions.


title: Towards a Unified Description of Isoscalar Giant Monopole Resonances in a Self-Consistent Quasiparticle-Vibration Coupling Approach

abstract: "Why is the EoS for tin so soft?" is a longstanding question, which prevents us from determining the nuclear incompressibility $K_\infty$ accurately. To solve this puzzle, a fully self-consistent quasiparticle random phase approximation (QRPA) plus quasiparticle-vibration coupling (QPVC) approach based on Skyrme-Hartree-Fock-Bogoliubov is developed. We show that the many-body correlations introduced by QPVC, which shift the ISGMR energy in Sn isotopes by about 0.4 MeV more than the energy in $^{208}$Pb, play a crucial role in providing a unified description of the ISGMR in Sn and Pb isotopes. The best description of the experimental strength functions is given by SV-K226 and KDE0, which are characterized by incompressibility values $K_\infty=$ 226 MeV and 229 MeV, respectively, at mean field level.


title: On the Navier-Stokes Equations with Stochastic Lie Transport

abstract: We prove the existence and uniqueness of maximal solutions to the 3D SALT (Stochastic Advection by Lie Transport, [Holm arXiv:1410.8311]) Navier-Stokes Equation in velocity and vorticity form, on the torus and the bounded domain respectively. The current work partners the paper [Goodair et al, arXiv:2209.09137] as an application of the abstract framework presented there, justifying the results first announced in [Goodair, arXiv:2202.09242v2]. In particular this represents the first well-posedness result for a fluid equation perturbed by a general transport type noise on a bounded domain.


title: Knowing the Past to Predict the Future: Reinforcement Virtual Learning

abstract: Reinforcement Learning (RL)-based control system has received considerable attention in recent decades. However, in many real-world problems, such as Batch Process Control, the environment is uncertain, which requires expensive interaction to acquire the state and reward values. In this paper, we present a cost-efficient framework, such that the RL model can evolve for itself in a Virtual Space using the predictive models with only historical data. The proposed framework enables a step-by-step RL model to predict the future state and select optimal actions for long-sight decisions. The main focuses are summarized as: 1) how to balance the long-sight and short-sight rewards with an optimal strategy; 2) how to make the virtual model interacting with real environment to converge to a final learning policy. Under the experimental settings of Fed-Batch Process, our method consistently outperforms the existing state-of-the-art methods.


title: Multi-Vector Retrieval as Sparse Alignment

abstract: Multi-vector retrieval models improve over single-vector dual encoders on many information retrieval tasks. In this paper, we cast the multi-vector retrieval problem as sparse alignment between query and document tokens. We propose AligneR, a novel multi-vector retrieval model that learns sparsified pairwise alignments between query and document tokens (e.g. `dog' vs. `puppy') and per-token unary saliences reflecting their relative importance for retrieval. We show that controlling the sparsity of pairwise token alignments often brings significant performance gains. While most factoid questions focusing on a specific part of a document require a smaller number of alignments, others requiring a broader understanding of a document favor a larger number of alignments. Unary saliences, on the other hand, decide whether a token ever needs to be aligned with others for retrieval (e.g. `kind' from `kind of currency is used in new zealand}'). With sparsified unary saliences, we are able to prune a large number of query and document token vectors and improve the efficiency of multi-vector retrieval. We learn the sparse unary saliences with entropy-regularized linear programming, which outperforms other methods to achieve sparsity. In a zero-shot setting, AligneR scores 51.1 points nDCG@10, achieving a new retriever-only state-of-the-art on 13 tasks in the BEIR benchmark. In addition, adapting pairwise alignments with a few examples (<= 8) further improves the performance up to 15.7 points nDCG@10 for argument retrieval tasks. The unary saliences of AligneR helps us to keep only 20% of the document token representations with minimal performance loss. We further show that our model often produces interpretable alignments and significantly improves its performance when initialized from larger language models.


title: On the Role of Fiducial Structures in Minisuperspace Reduction and Quantum Fluctuations in LQC

abstract: We study the homogeneous minisuperspace reduction within the canonical framework for a scalar field theory and gravity. Symmetry reduction is implemented via second class constraints for the field modes over a partitioning of the non-compact spatial slice $\Sigma$ into disjoint cells. The canonical structure of the resulting homogeneous theories is obtained via the associated Dirac bracket which can only be defined on a finite number of cells homogeneously patched together and agrees with the full theory Poisson bracket for the averaged fields. This identifies a finite region $V_o$, the fiducial cell, whose size $L$ sets the physical scale over which homogeneity is imposed, namely a wavelength cutoff. The reduced theory results from 1) selecting a subset of $V_o$-averaged observables of the full theory; 2) neglecting inhomogeneous $\vec k\neq\mathbf0$ modes with wavelengths $\lambda\geq L$ and $\lambda< L$; 3) neglecting boundary terms encoding interactions between neighbouring cells. The error made is of order $\mathcal O(1/kL)$. As a result, the off-shell structures of the reduced theory depend on the size of $V_o$ and different $V_o$ identify canonically inequivalent theories whose dynamics though is $V_o$-independent. Their quantisation leads then to a family of $V_o$-labeled quantum representations and the quantum version of an active rescaling of $V_o$ is implemented via a suitable dynamics-preserving isomorphism between the different theories. We discuss the consequences for statistical moments, fluctuations, and semiclassical states in both a standard and polymer quantisation. For a scalar field of mass $m$, we also sketch the quantum reduction and identify a subsector of the QFT where the results of the"first reduced, then quantised" theories can be reproduced with good approximation as long as $m\gg1/L$. Finally, a strategy to include inhomogeneities in cosmology is outlined.


title: Periods, Power Series, and Integrated Algebraic Numbers

abstract: Periods are defined as integrals of semialgebraic functions defined over the rationals. Periods form a countable ring not much is known about. Examples are given by taking the antiderivative of a power series which is algebraic over the polynomial ring over the rationals and evaluate it at a rational number. We follow this path and close these algebraic power series under taking iterated antiderivatives and nearby algebraic and geometric operations. We obtain a system of rings of power series whose coefficients form a countable real closed field. Using techniques from o-minimality we are able to show that every period belongs to this field. In the setting of o-minimality we define exponential integrated algebraic numbers and show that the Euler constant is an exponential integrated algebraic number. Hence they are a good candiate for a natural number system extending the period ring and containing important mathematical constants.


title: Secure and Efficient Privacy-preserving Authentication Scheme using Cuckoo Filter in Remote Patient Monitoring Network

abstract: With the ubiquitous advancement in smart medical devices and systems, the potential of Remote Patient Monitoring (RPM) network is evolving in modern healthcare systems. The medical professionals (doctors, nurses, or medical experts) can access vitals and sensitive physiological information about the patients and provide proper treatment to improve the quality of life through the RPM network. However, the wireless nature of communication in the RPM network makes it challenging to design an efficient mechanism for secure communication. Many authentication schemes have been proposed in recent years to ensure the security of the RPM network. Pseudonym, digital signature, and Authenticated Key Exchange (AKE) protocols are used for the Internet of Medical Things (IoMT) to develop secure authorization and privacy-preserving communication. However, traditional authentication protocols face overhead challenges due to maintaining a large set of key-pairs or pseudonyms results on the hospital cloud server. In this research work, we identify this research gap and propose a novel secure and efficient privacy-preserving authentication scheme using cuckoo filters for the RPM network. The use of cuckoo filters in our proposed scheme provides an efficient way for mutual anonymous authentication and a secret shared key establishment process between medical professionals and patients. Moreover, we identify the misbehaving sensor nodes using a correlation-based anomaly detection model to establish secure communication. The security analysis and formal security validation using SPAN and AVISPA tools show the robustness of our proposed scheme against message modification attacks, replay attacks, and man-in-the-middle attacks.


title: Initiator-free photo-crosslinkable cellulose-based resists for fabricating submicron patterns via direct laser writing

abstract: Novel bifunctional cellulose diacetate derivatives were synthesized in order to achieve bio-based photoresists, which can be structured by two-photon absorption via direct laser writing (DLW) without the need to use a photoinitiator. Therefore, cellulose diacetate is functionalized with thiol moieties and olefinic or methacrylic side groups enabling thiol-conjugated crosslinking. These cellulose derivatives are also photo-crosslinkable via UV irradiation ($\lambda$ = 254 nm and 365 nm) without using an initiator.


title: Measuring the polarization content of gravitational waves with strongly lensed binary black hole mergers

abstract: Alternative theories of gravity predict up to six distinct polarization modes for gravitational waves. Strong gravitational lensing of gravitational waves allows us to probe the polarization content of these signals by effectively increasing the number of observations from the same astrophysical source. The lensing time delays due to the multiple observed lensed images combined with the rotation of the Earth allows for effective non-collocated interferometers to be defined with respect to the source location and hence probe the alternative polarization amplitudes with more observations. To measure these amplitudes, we jointly fit the image observations to a single gravitational wave signal model that takes into account the image magnifications, time delays, and polarization mode amplitudes. We show that for certain systems, we can make a measurement of the relative mode amplitudes for lensed events with two detectable images.


title: The effects of delay on the HKB model of human motor coordination

abstract: In this paper, we analyse the celebrated Haken-Kelso-Bunz (HKB) model, describing the dynamics of bimanual coordination, in the presence of delay. We study the linear dynamics, stability, nonlinear behaviour and bifurcations of this model by both theoretical and numerical analysis. We calculate in-phase and anti-phase limit cycles as well as quasi-periodic solutions via double Hopf bifurcation analysis and centre manifold reduction. Moreover, we uncover further details on the global dynamic behaviour by numerical continuation, including the occurrence of limit cycles in phase quadrature and 1-1 locking of quasi-periodic solutions.


title: Self-attractive semiflexible polymers under an external force field

abstract: The dynamical response of a tethered semiflexible polymer with self-attractive interactions and subjected to an external force field is numerically investigated by varying stiffness and self-interaction strength. The chain is confined in two spatial dimensions and placed in contact with a heat bath described by the Brownian multiparticle collision method. For strong self-attraction the equilibrium conformations range from compact structures to double-stranded chains, and to rods when increasing the stiffness. Under the external field at small rigidities, the initial close-packed chain is continuously unwound by the force before being completely elongated. For double-stranded conformations the transition from the folded state to the open one is sharp being steeper for larger stiffnesses. The discontinuity in the transition appears in the force-extension relation as well as in the probability distribution function of the gyration radius. The relative deformation with respect to the equilibrium case along the direction normal to the force is found to decay as the inverse of the applied force.


title: Continuous Symmetry Breaking in a Trapped-Ion Spin Chain

abstract: One-dimensional systems exhibiting a continuous symmetry can host quantum phases of matter with true long-range order only in the presence of sufficiently long-range interactions. In most physical systems, however, the interactions are short-ranged, hindering the emergence of such phases in one dimension. Here we use a one-dimensional trapped-ion quantum simulator to prepare states with long-range spin order that extends over the system size of up to $23$ spins and is characteristic of the continuous symmetry-breaking phase of matter. Our preparation relies on simultaneous control over an array of tightly focused individual-addressing laser beams, generating long-range spin-spin interactions. We also observe a disordered phase with frustrated correlations. We further study the phases at different ranges of interaction and the out-of-equilibrium response to symmetry-breaking perturbations. This work opens an avenue to study new quantum phases and out-of-equilibrium dynamics in low-dimensional systems.

categories: quant-ph cond-mat.quant-gas cond-mat.stat-mech cond-mat.str-el


title: Estimating the technical wind energy potential of Kansas that incorporates the atmospheric response for policy applications

abstract: Energy scenarios and transition pathways need estimates of technical wind energy potentials. However, the standard policy-side approach uses observed wind speeds, thereby neglecting the effects of kinetic energy (KE) removal by the wind turbines that depletes the regional wind resource, lowers wind speeds, and reduces capacity factors. The standard approach therefore significantly overestimates the wind resource potential relative to estimates using numerical models of the atmosphere with interactive wind farm parameterizations. Here, we test the extent to which these effects of KE removal can be accounted for by our KE Budget of the Atmosphere (KEBA) approach over Kansas in the central US, a region with a high wind energy resource. We find that KEBA reproduces the simulated estimates within 10 - 11%, which are 30 - 50% lower than estimates using the standard approach. We also evaluate important differences in the depletion of the wind resource between daytime and nighttime conditions, which are due to effects of stability. Our results indicate that the KEBA approach is a simple yet adequate approach to evaluating regional-scale wind resource potentials, and that resource depletion effects need to be accounted for at such scales in policy applications.


title: Modern GPR Target Recognition Methods

abstract: Traditional GPR target recognition methods include pre-processing the data by removal of noisy signatures, dewowing (high-pass filtering to remove low-frequency noise), filtering, deconvolution, migration (correction of the effect of survey geometry), and can rely on the simulation of GPR responses. The techniques usually suffer from the loss of information, inability to adapt from prior results, and inefficient performance in the presence of strong clutter and noise. To address these challenges, several advanced processing methods have been developed over the past decade to enhance GPR target recognition. In this chapter, we provide an overview of these modern GPR processing techniques. In particular, we focus on the following methods: adaptive receive processing of range profiles depending on the target environment; adoption of learning-based methods so that the radar utilizes the results from prior measurements; application of methods that exploit the fact that the target scene is sparse in some domain or dictionary; application of advanced classification techniques; and convolutional coding which provides succinct and representatives features of the targets. We describe each of these techniques or their combinations through a representative application of landmine detection.


title: Pseudo-laminar chaos from on-off intermittency

abstract: In finite-dimensional, chaotic, Lorenz-like wave-particle dynamical systems one can find diffusive trajectories, which share their appearance with that of laminar chaotic diffusion [Phys. Rev. Lett. 128, 074101 (2022)] known from delay systems with lag-time modulation. Applying, however, to such systems a test for laminar chaos, as proposed in [Phys. Rev. E 101, 032213 (2020)], these signals fail such test, thus leading to the notion of pseudo-laminar chaos. The latter can be interpreted as integrated periodically driven on-off intermittency. We demonstrate that, on a signal level, true laminar and pseudo-laminar chaos are hardly distinguishable in systems with and without dynamical noise. However, very pronounced differences become apparent when correlations of signals and increments are considered. We compare and contrast these properties of pseudo-laminar chaos with true laminar chaos.


title: How to minimize the environmental contamination caused by hydrocarbon releases by onshore pipelines: The key role of a three-dimensional three-phase fluid flow numerical model

abstract: The contamination impact and the migration of the contaminant into the surrounding environment due to the presence of a spilled oil pipeline will cause significant damage to the natural ecosystem. For this reason, it is decisive to develop a rapid response strategy that might include accurate predictions of oil migration trajectories from numerical simulation modeling. In this paper, a three-dimensional model based on a high-resolution shock-capturing conservative method to resolve the nonlinear governing partial differential equations of the migration of a spilled light nonaqueous liquid oil contaminant in a variably saturated zone is employed to investigate the migration of the oil pipeline leakage with great accuracy. The effects on the oil type density, gasoline and diesel oil, the unsaturated zone depth, its saturation, the hydraulic gradient, and the pressure oil pipeline are investigated through the temporal evolution of the contaminant migration following the saturation profiles of the three-phase fluids flow in the variably saturated zone. The calculation results indicate that the leaking oil's pressure is the parameter that significantly affects the contaminants' arrival time to the groundwater table. Also, the water saturation of the unsaturated zone influences the arrival time as the water saturation increases for a fixed depth. The unsaturated zone depth significantly influences the contaminant migration unsaturated zone. At the same time, the oil density and the hydraulic gradient have limited effects on the contaminant migration in the variably saturated zone.


title: An Exponentially Converging Particle Method for the Mixed Nash Equilibrium of Continuous Games

abstract: We consider the problem of computing mixed Nash equilibria of two-player zero-sum games with continuous sets of pure strategies and with first-order access to the payoff function. This problem arises for example in game-theory-inspired machine learning applications, such as distributionally-robust learning. In those applications, the strategy sets are high-dimensional and thus methods based on discretisation cannot tractably return high-accuracy solutions. In this paper, we introduce and analyze a particle-based method that enjoys guaranteed local convergence for this problem. This method consists in parametrizing the mixed strategies as atomic measures and applying proximal point updates to both the atoms' weights and positions. It can be interpreted as a time-implicit discretization of the "interacting" Wasserstein-Fisher-Rao gradient flow. We prove that, under non-degeneracy assumptions, this method converges at an exponential rate to the exact mixed Nash equilibrium from any initialization satisfying a natural notion of closeness to optimality. We illustrate our results with numerical experiments and discuss applications to max-margin and distributionally-robust classification using two-layer neural networks, where our method has a natural interpretation as a simultaneous training of the network's weights and of the adversarial distribution.


title: Critical scaling through Gini index

abstract: In the systems showing critical behavior, various response functions have a singularity at the critical point. Therefore, as the driving field is tuned towards its critical value, the response functions change drastically, typically diverging with universal critical exponents. In this work, we quantify the inequality of response functions with measures traditionally used in economics, namely by constructing a Lorenz curve and calculating the corresponding Gini index. The scaling of such a response function, when written in terms of the Gini index, shows singularity at a point that is at least as universal as the corresponding critical exponent. The critical scaling, therefore, becomes a single parameter fit, which is a considerable simplification from the usual form where the critical point and critical exponents are independent. We also show that another measure of inequality, the Kolkata index, crosses the Gini index at a point just prior to the critical point. Therefore, monitoring these two inequality indices for a system where the critical point is not known, can produce a precursory signal for the imminent criticality. This could be useful in many systems, including that in condensed matter, bio- and geophysics to atmospheric physics. The generality and numerical validity of the calculations are shown with the Monte Carlo simulations of the two dimensional Ising model, site percolation on square lattice and the fiber bundle model of fracture.


title: Numerical integration of Schr\"odinger maps via the Hasimoto transform

abstract: We introduce a numerical approach to computing the Schr\"odinger map (SM) based on the Hasimoto transform which relates the SM flow to a cubic nonlinear Schr\"odinger (NLS) equation. In exploiting this nonlinear transform we are able to introduce the first fully explicit unconditionally stable symmetric integrators for the SM equation. Our approach consists of two parts: an integration of the NLS equation followed by the numerical evaluation of the Hasimoto transform. Motivated by the desire to study rough solutions to the SM equation, we also introduce a new symmetric low-regularity integrator for the NLS equation. This is combined with our novel fast low-regularity Hasimoto (FLowRH) transform, based on a tailored analysis of the resonance structures in the Magnus expansion and a fast realisation based on block-Toeplitz partitions, to yield an efficient low-regularity integrator for the SM equation. This scheme in particular allows us to obtain approximations to the SM in a more general regime (i.e. under lower regularity assumptions) than previously proposed methods. The favorable properties of our methods are exhibited both in theoretical convergence analysis and in numerical experiments.


title: Experimental superposition of time directions

abstract: In the macroscopic world, time is intrinsically asymmetric, flowing in a specific direction, from past to future. However, the same is not necessarily true for quantum systems, as some quantum processes produce valid quantum evolutions under time reversal. Supposing that such processes can be probed in both time directions, we can also consider quantum processes probed in a coherent superposition of forwards and backwards time directions. This yields a broader class of quantum processes than the ones considered so far in the literature, including those with indefinite causal order. In this work, we demonstrate for the first time an operation belonging to this new class: the quantum time flip. Using a photonic realisation of this operation, we apply it to a game formulated as a discrimination task between two sets of operators. This game not only serves as a witness of an indefinite time direction, but also allows for a computational advantage over strategies using a fixed time direction, and even those with an indefinite causal order.


title: Microscopic theory for the pair correlation function of liquidlike colloidal suspensions under shear flow

abstract: We present a theoretical framework to investigate the microscopic structure of concentrated hard-sphere colloidal suspensions under strong shear flows by fully taking into account the boundary-layer structure of convective diffusion. We solve the pair Smoluchowski equation with shear separately in the compressing and extensional sectors of the solid angle, by means of matched asymptotics. A proper, albeit approximate, treatment of the hydrodynamic interactions in the different sectors allows us to construct a potential of mean force containing the effect of the flow field on pair correlations. We insert the obtained pair potential in the Percus-Yevick relation and use the latter as a closure to solve the Ornstein-Zernike integral equation. For a wide range of either the packing fraction $\eta$ and the P\'eclet ($\textrm{Pe}$) number, we compute the pair correlation function and extract scaling laws for its value at contact. For all the considered value of $\textrm{Pe},$ we observe a very good agreement between theoretical findings and numerical results from literature, up to rather large values of $\eta.$ The theory predicts a consistent enhancement of the structure factor $ S(k)$ at $k \to 0,$ upon increasing the $\textrm{Pe}$ number. We argue this behaviour may signal the onset of a phase transition from the isotropic phase to a non-uniform one, induced by the external shear flow.

categories: cond-mat.soft cond-mat.dis-nn cond-mat.stat-mech physics.chem-ph physics.flu-dyn


title: FeynGKZ: a Mathematica package for solving Feynman integrals using GKZ hypergeometric systems

abstract: In the Lee-Pomeransky representation, Feynman integrals can be identified as a subset of Euler-Mellin integrals, which are known to satisfy Gel'fand-Kapranov-Zelevinsky (GKZ) system of partial differential equations. Here we present an automated package to derive the associated GKZ system for a given Feynman diagram and solve it in terms of hypergeometric functions using two equivalent algorithms, namely the triangulation method and the Gr\"obner deformation method. We present our code in the form of a Mathematica package FeynGKZ.wl which requires the softwares polymake, Macaulay2 and TOPCOM, and the packages AMBRE and Olsson.wl as dependencies. As applications of the package, we find series solutions to the GKZ systems of several one-loop and two-loop Feynman integrals. These are included in the file Examples.nb that can be downloaded along with the package from https://github.com/anant-group/FeynGKZ.


title: Novel results obtained by modeling of dynamic processes in superconductors: phase-slip centers as cooling engines

abstract: Based on a time-dependent Ginzburg-Landau system of equations and finite element modeling, we present novel results related with the physics of phase-slippage in superconducting wires surrounded by a non-superconductive environment. These results are obtained within our previously reported approach related to superconducting rings and superconductive gravitational wave detector transducers. It is shown that the phase-slip centers (PSCs) can be effective in originating not only positive but also negative thermal fluxes. With an appropriate design utilizing thermal diodes, PSCs can serve as cryocooling engines. Operating at $T\sim 1$ K cryostat cold-finger, they can achieve sub-Kelvin temperatures without using $^3$He.


title: Evaluating Impact of Social Media Posts by Executives on Stock Prices

abstract: Predicting stock market movements has always been of great interest to investors and an active area of research. Research has proven that popularity of products is highly influenced by what people talk about. Social media like Twitter, Reddit have become hotspots of such influences. This paper investigates the impact of social media posts on close price prediction of stocks using Twitter and Reddit posts. Our objective is to integrate sentiment of social media data with historical stock data and study its effect on closing prices using time series models. We carried out rigorous experiments and deep analysis using multiple deep learning based models on different datasets to study the influence of posts by executives and general people on the close price. Experimental results on multiple stocks (Apple and Tesla) and decentralised currencies (Bitcoin and Ethereum) consistently show improvements in prediction on including social media data and greater improvements on including executive posts.


title: Characterizing Intrinsic Compositionality in Transformers with Tree Projections

abstract: When trained on language data, do transformers learn some arbitrary computation that utilizes the full capacity of the architecture or do they learn a simpler, tree-like computation, hypothesized to underlie compositional meaning systems like human languages? There is an apparent tension between compositional accounts of human language understanding, which are based on a restricted bottom-up computational process, and the enormous success of neural models like transformers, which can route information arbitrarily between different parts of their input. One possibility is that these models, while extremely flexible in principle, in practice learn to interpret language hierarchically, ultimately building sentence representations close to those predictable by a bottom-up, tree-structured model. To evaluate this possibility, we describe an unsupervised and parameter-free method to \emph{functionally project} the behavior of any transformer into the space of tree-structured networks. Given an input sentence, we produce a binary tree that approximates the transformer's representation-building process and a score that captures how "tree-like" the transformer's behavior is on the input. While calculation of this score does not require training any additional models, it provably upper-bounds the fit between a transformer and any tree-structured approximation. Using this method, we show that transformers for three different tasks become more tree-like over the course of training, in some cases unsupervisedly recovering the same trees as supervised parsers. These trees, in turn, are predictive of model behavior, with more tree-like models generalizing better on tests of compositional generalization.


title: Boosting word frequencies in authorship attribution

abstract: In this paper, I introduce a simple method of computing relative word frequencies for authorship attribution and similar stylometric tasks. Rather than computing relative frequencies as the number of occurrences of a given word divided by the total number of tokens in a text, I argue that a more efficient normalization factor is the total number of relevant tokens only. The notion of relevant words includes synonyms and, usually, a few dozen other words in some ways semantically similar to a word in question. To determine such a semantic background, one of word embedding models can be used. The proposed method outperforms classical most-frequent-word approaches substantially, usually by a few percentage points depending on the input settings.


title: A Categorical Framework for Modeling with Stock and Flow Diagrams

abstract: Stock and flow diagrams are already an important tool in epidemiology, but category theory lets us go further and treat these diagrams as mathematical entities in their own right. In this chapter we use communicable disease models created with our software, StockFlow.jl, to explain the benefits of the categorical approach. We first explain the category of stock-flow diagrams and note the clear separation between the syntax of these diagrams and their semantics, demonstrating three examples of semantics already implemented in the software: ODEs, causal loop diagrams, and system structure diagrams. We then turn to two methods for building large stock-flow diagrams from smaller ones in a modular fashion: composition and stratification. Finally, we introduce the open-source ModelCollab software for diagram-based collaborative modeling. The graphical user interface of this web-based software lets modelers take advantage of the ideas discussed here without any knowledge of their categorical foundations.


title: SoK: A Stratified Approach to Blockchain Decentralization

abstract: Decentralization has been touted as the principal security advantage which propelled blockchain systems at the forefront of developments in the financial technology space. Its exact semantics nevertheless remain highly contested and ambiguous, with proponents and critics disagreeing widely on the level of decentralization offered by existing systems. To address this, we put forth a systematization of the current landscape with respect to decentralization and we derive a methodology that can help direct future research towards defining and measuring decentralization. Our approach dissects blockchain systems into multiple layers, or strata, each possibly encapsulating multiple categories, and it enables a unified method for measuring decentralization in each one. Our layers are (1) hardware, (2) software, (3) network, (4) consensus, (5) economics ("tokenomics"), (6) client API, (7) governance, and (8) geography. Armed with this stratification, we examine for each layer which pertinent properties of distributed ledgers (safety, liveness, privacy, stability) can be at risk due to centralization and in what way. We also introduce a practical test, the "Minimum Decentralization Test" which can provide quick insights about the decentralization state of a blockchain system. To demonstrate how our stratified methodology can be used in practice, we apply it fully (layer by layer) to Bitcoin, and we provide examples of systems which comprise one or more "problematic" layers that cause them to fail the MDT. Our work highlights the challenges in measuring and achieving decentralization, and suggests various potential directions where future research is needed.


title: Learning an Artificial Language for Knowledge-Sharing in Multilingual Translation

abstract: The cornerstone of multilingual neural translation is shared representations across languages. Given the theoretically infinite representation power of neural networks, semantically identical sentences are likely represented differently. While representing sentences in the continuous latent space ensures expressiveness, it introduces the risk of capturing of irrelevant features which hinders the learning of a common representation. In this work, we discretize the encoder output latent space of multilingual models by assigning encoder states to entries in a codebook, which in effect represents source sentences in a new artificial language. This discretization process not only offers a new way to interpret the otherwise black-box model representations, but, more importantly, gives potential for increasing robustness in unseen testing conditions. We validate our approach on large-scale experiments with realistic data volumes and domains. When tested in zero-shot conditions, our approach is competitive with two strong alternatives from the literature. We also use the learned artificial language to analyze model behavior, and discover that using a similar bridge language increases knowledge-sharing among the remaining languages.


title: DC-cycleGAN: Bidirectional CT-to-MR Synthesis from Unpaired Data

abstract: Magnetic resonance (MR) and computer tomography (CT) images are two typical types of medical images that provide mutually-complementary information for accurate clinical diagnosis and treatment. However, obtaining both images may be limited due to some considerations such as cost, radiation dose and modality missing. Recently, medical image synthesis has aroused gaining research interest to cope with this limitation. In this paper, we propose a bidirectional learning model, denoted as dual contrast cycleGAN (DC-cycleGAN), to synthesize medical images from unpaired data. Specifically, a dual contrast loss is introduced into the discriminators to indirectly build constraints between real source and synthetic images by taking advantage of samples from the source domain as negative samples and enforce the synthetic images to fall far away from the source domain. In addition, cross-entropy and structural similarity index (SSIM) are integrated into the DC-cycleGAN in order to consider both the luminance and structure of samples when synthesizing images. The experimental results indicate that DC-cycleGAN is able to produce promising results as compared with other cycleGAN-based medical image synthesis methods such as cycleGAN, RegGAN, DualGAN, and NiceGAN. The code will be available at https://github.com/JiayuanWang-JW/DC-cycleGAN.


title: Driver Digital Twin for Online Prediction of Personalized Lane Change Behavior

abstract: Connected and automated vehicles (CAVs) are supposed to share the road with human-driven vehicles (HDVs) in a foreseeable future. Therefore, considering the mixed traffic environment is more pragmatic, as the well-planned operation of CAVs may be interrupted by HDVs. In the circumstance that human behaviors have significant impacts, CAVs need to understand HDV behaviors to make safe actions. In this study, we develop a Driver Digital Twin (DDT) for the online prediction of personalized lane change behavior, allowing CAVs to predict surrounding vehicles' behaviors with the help of the digital twin technology. DDT is deployed on a vehicle-edge-cloud architecture, where the cloud server models the driver behavior for each HDV based on the historical naturalistic driving data, while the edge server processes the real-time data from each driver with his/her digital twin on the cloud to predict the lane change maneuver. The proposed system is first evaluated on a human-in-the-loop co-simulation platform, and then in a field implementation with three passenger vehicles connected through the 4G/LTE cellular network. The lane change intention can be recognized in 6 seconds on average before the vehicle crosses the lane separation line, and the Mean Euclidean Distance between the predicted trajectory and GPS ground truth is 1.03 meters within a 4-second prediction window. Compared to the general model, using a personalized model can improve prediction accuracy by 27.8%. The demonstration video of the proposed system can be watched at https://youtu.be/5cbsabgIOdM.


title: A Unified Framework for Symmetry Handling

abstract: Handling symmetries in optimization problems is essential for devising efficient solution methods. In this article, we present a general framework that captures many of the already existing symmetry handling methods. While these methods are mostly discussed independently from each other, our framework allows to apply different methods simultaneously and thus outperforming their individual effect. Moreover, most existing symmetry handling methods only apply to binary variables. Our framework allows to easily generalize these methods to general variable types. Numerical experiments confirm that our novel framework is superior to the state-of-the-art symmetry handling methods as implemented in the solver SCIP on a broad set of instances.


title: Machine Learning for Metasurfaces Design and Their Applications

abstract: Metasurfaces (MTSs) are increasingly emerging as enabling technologies to meet the demands for multi-functional, small form-factor, efficient, reconfigurable, tunable, and low-cost radio-frequency (RF) components because of their ability to manipulate waves in a sub-wavelength thickness through modified boundary conditions. They enable the design of reconfigurable intelligent surfaces (RISs) for adaptable wireless channels and smart radio environments, wherein the inherently stochastic nature of the wireless environment is transformed into a programmable propagation channel. In particular, space-limited RF applications, such as communications and radar, that have strict radiation requirements are currently being investigated for potential RIS deployment. The RIS comprises sub-wavelength units or meta-atoms, which are independently controlled and whose geometry and material determine the spectral response of the RIS. Conventionally, designing RIS to yield the desired EM response requires trial and error by iteratively investigating a large possibility of various geometries and materials through thousands of full-wave EM simulations. In this context, machine/deep learning (ML/DL) techniques are proving critical in reducing the computational cost and time of RIS inverse design. Instead of explicitly solving Maxwell's equations, DL models learn physics-based relationships through supervised training data. The ML/DL techniques also aid in RIS deployment for numerous wireless applications, which requires dealing with multiple channel links between the base station (BS) and the users. As a result, the BS and RIS beamformers require a joint design, wherein the RIS elements must be rapidly reconfigured. This chapter provides a synopsis of DL techniques for both inverse RIS design and RIS-assisted wireless systems.


title: Self-Attentive Sequential Recommendation with Cheap Causal Convolutions

abstract: Sequential Recommendation is a prominent topic in current research, which uses user behavior sequence as an input to predict future behavior. By assessing the correlation strength of historical behavior through the dot product, the model based on the self-attention mechanism can capture the long-term preference of the sequence. However, it has two limitations. On the one hand, it does not effectively utilize the items' local context information when determining the attention and creating the sequence representation. On the other hand, the convolution and linear layers often contain redundant information, which limits the ability to encode sequences. In this paper, we propose a self-attentive sequential recommendation model based on cheap causal convolution. It utilizes causal convolutions to capture items' local information for calculating attention and generating sequence embedding. It also uses cheap convolutions to improve the representations by lightweight structure. We evaluate the effectiveness of the proposed model in terms of both accurate and calibrated sequential recommendation. Experiments on benchmark datasets show that the proposed model can perform better in single- and multi-objective recommendation scenarios.


title: Contract Composition for Dynamical Control Systems: Definition and Verification using Linear Programming

abstract: Designing large-scale control systems to satisfy complex specifications is hard in practice, as most formal methods are limited to systems of modest size. Contract theory has been proposed as a modular alternative to formal methods in control, in which specifications are defined by assumptions on the input to a component and guarantees on its output. However, current contract-based methods for control systems either prescribe guarantees on the state of the system, going against the spirit of contract theory, or can only support rudimentary compositions. In this paper, we present a contract-based modular framework for discrete-time dynamical control systems. We extend the definition of contracts by allowing the assumption on the input at a time $k$ to depend on outputs up to time $k-1$, which is essential when considering the feedback connection of an unregulated dynamical system and a controller. We also define contract composition for arbitrary interconnection topologies, under the pretence of well-posedness, and prove that this notion supports modular design, analysis and verification. This is done using graph theory methods, and specifically using the notions of topological ordering and backward-reachable nodes. Lastly, we use $k$-induction to present an algorithm for verifying vertical contracts, which are claims of the form "the conjugation of given component-level contracts is a stronger specification than a given contract on the integrated system". These algorithms are based on linear programming, and scale linearly with the number of components in the interconnected network. A numerical example is provided to demonstrate the scalability of the presented approach, as well as the modularity achieved by using it.