title: Pirouette: Higher-Order Typed Functional Choreographies

abstract: We present Pirouette, a language for typed higher-order functional choreographic programming. Pirouette offers programmers the ability to write a centralized functional program and compile it via endpoint projection into programs for each node in a distributed system. Moreover, Pirouette is defined generically over a (local) language of messages, and lifts guarantees about the message type system to its own. Message type soundness also guarantees deadlock freedom. All of our results are verified in Coq.


title: Cross Modality 3D Navigation Using Reinforcement Learning and Neural Style Transfer

abstract: This paper presents the use of Multi-Agent Reinforcement Learning (MARL) to perform navigation in 3D anatomical volumes from medical imaging. We utilize Neural Style Transfer to create synthetic Computed Tomography (CT) agent gym environments and assess the generalization capabilities of our agents to clinical CT volumes. Our framework does not require any labelled clinical data and integrates easily with several image translation techniques, enabling cross modality applications. Further, we solely condition our agents on 2D slices, breaking grounds for 3D guidance in much more difficult imaging modalities, such as ultrasound imaging. This is an important step towards user guidance during the acquisition of standardised diagnostic view planes, improving diagnostic consistency and facilitating better case comparison.


title: Guaranteed blind deconvolution and demixing via hierarchically sparse reconstruction

abstract: The blind deconvolution problem amounts to reconstructing both a signal and a filter from the convolution of these two. It constitutes a prominent topic in mathematical and engineering literature. In this work, we analyze a sparse version of the problem: The filter $h\in \mathbb{R}^\mu$ is assumed to be $s$-sparse, and the signal $b \in \mathbb{R}^n$ is taken to be $\sigma$-sparse, both supports being unknown. We observe a convolution between the filter and a linear transformation of the signal. Motivated by practically important multi-user communication applications, we derive a recovery guarantee for the simultaneous demixing and deconvolution setting. We achieve efficient recovery by relaxing the problem to a hierarchical sparse recovery for which we can build on a flexible framework. At the same time, for this we pay the price of some sub-optimal guarantees compared to the number of free parameters of the problem. The signal model we consider is sufficiently general to capture many applications in a number of engineering fields. Despite their practical importance, we provide first rigorous performance guarantees for efficient and simple algorithms for the bi-sparse and generalized demixing setting. We complement our analytical results by presenting results of numerical simulations. We find evidence that the sub-optimal scaling $s^2\sigma \log(\mu)\log(n)$ of our derived sufficient condition is likely overly pessimistic and that the observed performance is better described by a scaling proportional to $ s\sigma$ up to log-factors.


title: Quantum phase transition in the one-dimensional Dicke-Hubbard model with coupled qubits

abstract: We study the ground state phase diagram of a one-dimensional two qubits Dicke-Hubbard model with XY qubit-qubit interaction. We use a numerical method combing the cluster mean-field theory and the matrix product state(MPS) to obtain the exact wave function of the ground state. When counter-rotating wave terms(CRTs) in the qubit-cavity coupling are neglected, we observe a rich phase diagram including a quantum phase transition between the Mott-insulating phase and the superfluid phase. This phase transition can be either the first-order or the second-order type depending on whether the total angular momentum changes across the phase diagram. Moreover, we observe two quantum triple points, at which three different phases coexist, with both positive and negative XY interactions. By further considering the effect of CRTs, we find that the main feature in the previous phase diagram, including the existence of quantum triple points, is retained. We also show that CRTs extremely demolish the non-local correlations in the coherent phase.


title: SLA-Driven Load Scheduling in Multi-Tier Cloud Computing: Financial Impact Considerations

abstract: A cloud service provider strives to provide a high Quality of Service (QoS) to client jobs. Such jobs vary in computational and Service-Level-Agreement (SLA) obligations, as well as differ with respect to tolerating delays and SLA violations. The job scheduling plays a critical role in servicing cloud demands by allocating appropriate resources to execute client jobs. The response to such jobs is optimized by the cloud provider on a multi-tier cloud computing environment. Typically, the complex and dynamic nature of multi-tier environments incurs difficulties in meeting such demands, because tiers are dependent on each other which in turn makes bottlenecks of a tier shift to escalate in subsequent tiers. However, the optimization process of existing approaches produces single-tier-driven schedules that do not employ the differential impact of SLA violations in executing client jobs. Furthermore, the impact of schedules optimized at the tier level on the performance of schedules formulated in subsequent tiers tends to be ignored, resulting in a less than optimal performance when measured at the multi-tier level. Thus, failing in committing job obligations incurs SLA penalties that often take the form of either financial compensations, or losing future interests and motivations of unsatisfied clients in the service provided. In this paper, a scheduling and allocation approach is proposed to formulate schedules that account for differential impacts of SLA violation penalties and, thus, produce schedules that are optimal in financial performance. A queue virtualization scheme is designed to facilitate the formulation of optimal schedules at the tier and multi-tier levels of the cloud environment. Because the scheduling problem is NP-hard, a biologically inspired approach is proposed to mitigate the complexity of finding optimal schedules.


title: Sound Absorption and Transmission Loss Properties of Open-Celled Aluminum Foams with Stepwise Relative Density Gradients

abstract: We investigate the acoustical properties of uncompressed and compressed open-celled aluminum metal foams fabricated using a directional solidification foaming process. We compressed the fabricated foams using a hydraulic press to different compression ratios and characterized the effect of compression on the cellular microstructure using microtomography and scanning electron microscopy. The static airflow resistances of the samples are measured and related to the observed microstructural changes. We measured the normal incidence acoustical properties using two- and four-microphone impedance tube methods and show that the compression substantially improves their sound absorption and transmission loss performance. We then stack individual disks with different compression ratios to create various stepwise relative density gradient configurations and show that stepwise gradients provide a significant improvement in properties as compared to the uncompressed sample. The effect of increasing and decreasing relative density gradients on the overall absorption and transmission loss behavior is characterized. Finally, we use an experimentally informed and validated transfer matrix method to predict the effect of various layer thicknesses and stacking sequences on the global acoustical properties. Our results show that open-celled metal foams with stepwise relative density gradients can be designed to provide tailored acoustic absorption performance while reducing the overall weight of the noise reduction package.


title: Positive Definite Multi-Kernels for Scattered Data Interpolations

abstract: In this article, we use the knowledge of positive definite tensors to develop a concept of positive definite multi-kernels to construct the kernel-based interpolants of scattered data. By the techniques of reproducing kernel Banach spaces, the optimal recoveries and error analysis of the kernel-based interpolants are shown for a special class of strictly positive definite multi-kernels.


title: Sampling Methods for Bayesian Inference Involving Convergent Noisy Approximations of Forward Maps

abstract: We present Bayesian techniques for solving inverse problems which involve mean-square convergent random approximations of the forward map. Noisy approximations of the forward map arise in several fields, such as multiscale problems and probabilistic numerical methods. In these fields, a random approximation can enhance the quality or the efficiency of the inference procedure, but entails additional theoretical and computational difficulties due to the randomness of the forward map. A standard technique to address this issue is to combine Monte Carlo averaging with Markov chain Monte Carlo samplers, as for example in the pseudo-marginal Metropolis--Hastings methods. In this paper, we consider mean-square convergent random approximations, and quantify how Monte Carlo errors propagate from the forward map to the solution of the inverse problems. Moreover, we review and describe simple techniques to solve such inverse problems, and compare performances with a series of numerical experiments.


title: Fast FPT-Approximation of Branchwidth

abstract: Branchwidth determines how graphs, and more generally, arbitrary connectivity (basically symmetric and submodular) functions could be decomposed into a tree-like structure by specific cuts. We develop a general framework for designing fixed-parameter tractable (FPT) 2-approximation algorithms for branchwidth of connectivity functions. The first ingredient of our framework is combinatorial. We prove a structural theorem establishing that either a sequence of particular refinement operations could decrease the width of a branch decomposition or that the width of the decomposition is already within a factor of 2 from the optimum. The second ingredient is an efficient implementation of the refinement operations for branch decompositions that support efficient dynamic programming. We present two concrete applications of our general framework. $\bullet$ An algorithm that for a given $n$-vertex graph $G$ and integer $k$ in time $2^{2^{O(k)}} n^2$ either constructs a rank decomposition of $G$ of width at most $2k$ or concludes that the rankwidth of $G$ is more than $k$. It also yields a $(2^{2k+1}-1)$-approximation algorithm for cliquewidth within the same time complexity, which in turn, improves to $f(k)n^2$ the running times of various algorithms on graphs of cliquewidth $k$. Breaking the "cubic barrier" for rankwidth and cliquewidth was an open problem in the area. $\bullet$ An algorithm that for a given $n$-vertex graph $G$ and integer $k$ in time $2^{O(k)} n$ either constructs a branch decomposition of $G$ of width at most $2k$ or concludes that the branchwidth of $G$ is more than $k$. This improves over the 3-approximation that follows from the recent treewidth 2-approximation of Korhonen [FOCS 2021].


title: Self-Regularization in turbulence from the Kolmogorov 4/5-Law and Alignment

abstract: A defining feature of 3D hydrodynamic turbulence is that the rate of energy dissipation is bounded away from zero as viscosity is decreased (Reynolds number increased). This phenomenon - anomalous dissipation - is sometimes called the `zeroth law of turbulence' as it underpins many celebrated theoretical predictions. Another robust feature observed in turbulence is that velocity structure functions $S_p(\ell) :=\langle |\delta_\ell u|^p\rangle$ exhibit persistent power-law scaling in the inertial range, namely $S_p(\ell) \sim |\ell|^{\zeta_p}$ for exponents $\zeta_p>0$ over an ever-increasing (with Reynolds) range of scales. This behavior indicates that the velocity field retains some fractional differentiability uniformly in the Reynolds number. The Kolmogorov 1941 theory of turbulence predicts that $\zeta_p=p/3$ for all $p$ and Onsager's 1949 theory establishes the requirement that $\zeta_p\leq p/3$ for $p\geq 3$ for consistency with the zeroth law. Empirically, $\zeta_2 \gtrapprox 2/3$ and $\zeta_3 \lessapprox 1$, suggesting that turbulent Navier-Stokes solutions approximate dissipative weak solutions of the Euler equations possessing (nearly) the minimal degree of singularity required to sustain anomalous dissipation. In this note, we adopt an experimentally supported hypothesis on the anti-alignment of velocity increments with their separation vectors and demonstrate that the inertial dissipation provides a regularization mechanism via the Kolmogorov 4/5-law.


title: Exponential stability of Timoshenko-Gurtin-Pipkin systems with full thermal coupling

abstract: We analyze the stability properties of a linear thermoelastic Timoshenko-Gurtin-Pipkin system with thermal coupling acting on both the shear force and the bending moment. Under either the mixed Dirichlet-Neumann or else the full Dirichlet boundary conditions, we show that the associated solution semigroup in the history space framework of Dafermos is exponentially stable independently of the values of the structural parameters of the model.


title: Automated Supervised Feature Selection for Differentiated Patterns of Care

abstract: An automated feature selection pipeline was developed using several state-of-the-art feature selection techniques to select optimal features for Differentiating Patterns of Care (DPOC). The pipeline included three types of feature selection techniques; Filters, Wrappers and Embedded methods to select the top K features. Five different datasets with binary dependent variables were used and their different top K optimal features selected. The selected features were tested in the existing multi-dimensional subset scanning (MDSS) where the most anomalous subpopulations, most anomalous subsets, propensity scores, and effect of measures were recorded to test their performance. This performance was compared with four similar metrics gained after using all covariates in the dataset in the MDSS pipeline. We found out that despite the different feature selection techniques used, the data distribution is key to note when determining the technique to use.


title: Monitoring geometrical properties of word embeddings for detecting the emergence of new topics

abstract: Slow emerging topic detection is a task between event detection, where we aggregate behaviors of different words on short period of time, and language evolution, where we monitor their long term evolution. In this work, we tackle the problem of early detection of slowly emerging new topics. To this end, we gather evidence of weak signals at the word level. We propose to monitor the behavior of words representation in an embedding space and use one of its geometrical properties to characterize the emergence of topics. As evaluation is typically hard for this kind of task, we present a framework for quantitative evaluation. We show positive results that outperform state-of-the-art methods on two public datasets of press and scientific articles.


title: Markov Chain Approximations to Stochastic Differential Equations by Recombination on Lattice Trees

abstract: We revisit the classical problem of approximating a stochastic differential equation by a discrete-time and discrete-space Markov chain. Our construction iterates Caratheodory's theorem over time to match the moments of the increments locally. This allows to construct a Markov chain with a sparse transition matrix where the number of attainable states grows at most polynomially as time increases. Moreover, the MC evolves on a tree whose nodes lie on a "universal lattice" in the sense that an arbitrary number of different SDEs can be approximated on the same tree. The construction is not tailored to specific models, we discuss both the case of uni-variate and multi-variate case SDEs, and provide an implementation and numerical experiments.


title: Bridging the Gap: Categorizing Gravitational-Wave Events at the Transition Between Neutron Stars and Black Holes

abstract: We search for features in the mass distribution of detected compact binary coalescences which signify the transition between neutron stars and black holes. We analyze all gravitational wave detections by LIGO-Virgo-KAGRA made through the end of the first half of the third observing run, and find clear evidence for two different populations of compact objects based solely on gravitational wave data. We confidently (99.3%) find a steepening relative to a single power law describing neutron stars and low-mass black holes below $2.4^{+0.5}_{-0.5}\text{ M}_{\odot}$, which is consistent with many predictions for the maximum neutron star mass. We find suggestions of the purported lower mass gap between the most massive neutron stars and the least massive black holes, but are unable to conclusively resolve it with current data. If it exists, we find the lower mass gap's edges to lie at $2.2^{+0.7}_{-0.5}\text{ M}_{\odot}$ and $6.0^{+2.4}_{-1.4}\text{ M}_{\odot}$. We re-examine events that have been deemed "exceptional" by the LIGO-Virgo-KAGRA collaborations in the context of these features. We analyze GW190814 self-consistently in the context of the full population of compact binaries, finding support for its secondary to be either a neutron star or a lower mass gap object, consistent with previous claims. Our models are the first to accommodate this event, which is an outlier with respect to the binary black hole population. We find that GW200105 and GW200115 probe the edges of, and may have components within, the lower mass gap. As future data improve global population models, the classification of these events will also improve.


title: Justification of the KP-II approximation in dynamics of two-dimensional FPU systems

abstract: Dynamics of the Fermi-Pasta-Ulam (FPU) system on a two-dimensional square lattice is considered in the limit of small-amplitude long-scale waves with slow transverse modulations. In the absence of transverse modulations, dynamics of such waves, even at an oblique angle with respect to the square lattice, is known to be described by the Korteweg-de Vries (KdV) equation. For the three basic directions (horizontal, vertical, and diagonal), we prove that the modulated waves are well described by the Kadomtsev-Petviashvili (KP-II) equation. The result was expected long ago but proving rigorous bounds on the approximation error turns out to be complicated due to the nonlocal terms of the KP-II equation and the vector structure of the FPU systems on two-dimensional lattices. We have obtained these error bounds by extending the local well-posedness result for the KP-II equation in Sobolev spaces and by controlling the error terms with energy estimates. The bounds are useful in the analysis of transverse stability of solitary and periodic waves in two-dimensional FPU systems due to many results available for the KP-II equation.

categories: math.AP math-ph math.DS math.MP nlin.PS nlin.SI


title: Entanglement Transitions from Stochastic Resetting of Non-Hermitian Quasiparticles

abstract: We put forward a phenomenological theory for entanglement dynamics in monitored quantum many-body systems with well-defined quasiparticles. Within this theory entanglement is carried by ballistically propagating non-Hermitian quasiparticles which are stochastically reset by the measurement protocol with rate given by their finite inverse lifetime. We write down a renewal equation for the statistics of the entanglement entropy and show that depending on the spectrum of quasiparticle decay rates different entanglement scaling can arise and even sharp entanglement phase transitions. When applied to a Quantum Ising chain where the transverse magnetization is measured by quantum jumps, our theory predicts a critical phase with logarithmic scaling of the entanglement, an area law phase and a continuous phase transition between them, with an effective central charge vanishing as a square root at the transition point. We compare these predictions with with exact numerical calculations on the same model and find an excellent agreement.


title: Model Checking Temporal Properties of Recursive Probabilistic Programs

abstract: Probabilistic pushdown automata (pPDA) are a standard operational model for programming languages involving discrete random choices and recursive procedures. Temporal properties are useful for specifying the chronological order of events during program execution. Existing approaches for model checking pPDA against temporal properties have focused mostly on $\omega$-regular and LTL properties. In this paper, we give decidability and complexity results for the model checking problem of pPDA against $\omega$-visibly pushdown languages that can be described by specification logics such as CaRet. These logical formulae allow specifying properties that explicitly take the structured computations arising from procedural programs into account. For example, CaRet is able to match procedure calls with their corresponding future returns, and thus allows to express fundamental program properties such as total and partial correctness.


title: Analytic Hochschild-Kostant-Rosenberg Theorem

abstract: Let $R$ be a Banach ring. We prove that the category of chain complexes of complete bornological $R$-modules (and several related categories) is a derived algebraic context in the sense of Raksit. We then use the framework of derived algebra to prove a version of the Hochschild-Kostant-Rosenberg Theorem, which relates the circle action on the Hochschild algebra to the de Rham-differential-enriched-de Rham algebra of a simplicial, commutative, complete bornological algebra. This has a geometric interpretation in the language of derived analytic geometry, namely, the derived loop stack of a derived analytic stack is equivalent to the shifted tangent stack. Using this geometric interpretation we extend our results to derived schemes.


title: Period Change Rates in Large Magellanic Cloud Cepheids Revisited

abstract: The period-change rate (PCR) of pulsating variable stars is a useful probe of changes in their interior structure, and thus of their evolutionary stages. So far, the PCRs of Classical Cepheids in the Large Magellanic Cloud (LMC) have been explored in a limited sample of the total population of these variables. Here we use a template-based method to build observed minus computed (O-C) period diagrams, from which we can derive PCRs for these stars by taking advantage of the long time baseline afforded by the Digital Access to a Sky Century @ Harvard (DASCH) light curves, combined with additional data from the Optical Gravitational Lensing Experiment (OGLE), the MAssive Compact Halo Object (MACHO) project, Gaia's Data Release 2, and in some cases the All-Sky Automated Survey (ASAS). From an initial sample of 2315 sources, our method provides an unprecedented sample of 1303 LMC Classical Cepheids with accurate PCRs, the largest for any single galaxy, including the Milky Way. The derived PCRs are largely compatible with theoretically expected values, as computed by our team using the Modules for Experiments in Stellar Astrophysics (MESA) code, as well as with similar previous computations available in the literature. Additionally, five long-period (P>50 d) sources display a cyclic behavior in their O-C diagrams, which is clearly incompatible with evolutionary changes. Finally, on the basis of their large positive PCR values, two first-crossing Cepheid candidates are identified.


title: Deep-Learning Based Linear Precoding for MIMO Channels with Finite-Alphabet Signaling

abstract: This paper studies the problem of linear precoding for multiple-input multiple-output (MIMO) communication channels employing finite-alphabet signaling. Existing solutions typically suffer from high computational complexity due to costly computations of the constellation-constrained mutual information. In contrast to existing works, this paper takes a different path of tackling the MIMO precoding problem. Namely, a data-driven approach, based on deep learning, is proposed. In the offline training phase, a deep neural network learns the optimal solution on a set of MIMO channel matrices. This allows the reduction of the computational complexity of the precoder optimization in the online inference phase. Numerical results demonstrate the efficiency of the proposed solution vis-\`a-vis existing precoding algorithms in terms of significantly reduced complexity and close-to-optimal performance.


title: Visualizing the Emergence of Intermediate Visual Patterns in DNNs

abstract: This paper proposes a method to visualize the discrimination power of intermediate-layer visual patterns encoded by a DNN. Specifically, we visualize (1) how the DNN gradually learns regional visual patterns in each intermediate layer during the training process, and (2) the effects of the DNN using non-discriminative patterns in low layers to construct disciminative patterns in middle/high layers through the forward propagation. Based on our visualization method, we can quantify knowledge points (i.e., the number of discriminative visual patterns) learned by the DNN to evaluate the representation capacity of the DNN. Furthermore, this method also provides new insights into signal-processing behaviors of existing deep-learning techniques, such as adversarial attacks and knowledge distillation.


title: Detecting topological phase transitions through entanglement between disconnected partitions in a Kitaev chain with long-range interactions

abstract: We explore the behaviour of the disconnected entanglement entropy (DEE) across the topological phases of a long range interacting Kitaev chain where the long range interactions decay as a power law with an exponent $\alpha$. We show that while the DEE may not remain invariant deep within the topologically non-trivial phase when $\alpha<1$, it nevertheless shows a quantized discontinuous jump at the quantum critical point and can act as a strong marker for the detection of topological phase transition. We also study the time evolution of the DEE after a sudden quench of the chemical potential within the same phase. In the short range limit of a finite chain, the DEE is expected to remain constant upto a critical time after the quench, which diverges in the thermodynamic limit. However, no such critical time is found to exist when the long range interactions dominate (i.e., $\alpha<1$).


title: Localized states emerging from singular and nonsingular flat bands in a frustrated fractal-like photonic lattice

abstract: We report on singular and nonsingular flat bands in a Sierpinski fractal-like photonic lattice. We demonstrate that the the lowest two bands, being isolated and degenerate due to geometrical frustration, are nonsingular and thus can be spanned by a complete set of compact localized states. We experimentally prove these states to propagate diffractionless in the photonic lattice. Our results reveal the interplay between geometrical frustration, degenerate flat bands and compact localized states in a single photonic lattice, and pave the way to photonic spin liquid ground states.


title: On the convergence of the orthogonal spectral sequence

abstract: We show that the orthogonal spectral sequence introduced by the second author is strongly convergent in Voevodsky's triangulated category of motives DM over a field k. In the context of the Morel-Voevodsky motivic stable homotopy category we provide concrete examples where the spectral sequence is not strongly convergent, and give a criterion under which the strong convergence still holds. This criterion holds for Voevodsky's slices, and as a consequence we obtain a spectral sequence which converges strongly to the E1-term of Voevodsky's slice spectral sequence.


title: Regularization Graphs -- A unified framework for variational regularization of inverse problems

abstract: We introduce and study a mathematical framework for a broad class of regularization functionals for ill-posed inverse problems: Regularization Graphs. Regularization graphs allow to construct functionals using as building blocks linear operators and convex functionals, assembled by means of operators that can be seen as generalizations of classical infimal convolution operators. This class of functionals exhaustively covers existing regularization approaches and it is flexible enough to craft new ones in a simple and constructive way. We provide well-posedness and convergence results with the proposed class of functionals in a general setting. Further, we consider a bilevel optimization approach to learn optimal weights for such regularization graphs from training data. We demonstrate that this approach is capable of optimizing the structure and the complexity of a regularization graph, allowing, for example, to automatically select a combination of regularizers that is optimal for given training data.


title: Wright functions of the second kind and Whittaker functions

abstract: In the framework of higher transcendental functions the Wright functions of the second kind have increased their relevance resulting from their applications in probability theory and, in particular, in fractional diffusion processes. Here, these functions are compared with the well-known Whittaker functions in some special cases of fractional order. In addition, we point out two erroneous representations in the literature.


title: Disengagement Cause-and-Effect Relationships Extraction Using an NLP Pipeline

abstract: The advancement in machine learning and artificial intelligence is promoting the testing and deployment of autonomous vehicles (AVs) on public roads. The California Department of Motor Vehicles (CA DMV) has launched the Autonomous Vehicle Tester Program, which collects and releases reports related to Autonomous Vehicle Disengagement (AVD) from autonomous driving. Understanding the causes of AVD is critical to improving the safety and stability of the AV system and provide guidance for AV testing and deployment. In this work, a scalable end-to-end pipeline is constructed to collect, process, model, and analyze the disengagement reports released from 2014 to 2020 using natural language processing deep transfer learning. The analysis of disengagement data using taxonomy, visualization and statistical tests revealed the trends of AV testing, categorized cause frequency, and significant relationships between causes and effects of AVD. We found that (1) manufacturers tested AVs intensively during the Spring and/or Winter, (2) test drivers initiated more than 80% of the disengagement while more than 75% of the disengagement were led by errors in perception, localization & mapping, planning and control of the AV system itself, and (3) there was a significant relationship between the initiator of AVD and the cause category. This study serves as a successful practice of deep transfer learning using pre-trained models and generates a consolidated disengagement database allowing further investigation for other researchers.


title: Data Selection for Efficient Model Update in Federated Learning

abstract: The Federated Learning (FL) workflow of training a centralized model with distributed data is growing in popularity. However, until recently, this was the realm of contributing clients with similar computing capability. The fast expanding IoT space and data being generated and processed at the edge are encouraging more effort into expanding federated learning to include heterogeneous systems. Previous approaches distribute light-weight models to clients are rely on knowledge transfer to distil the characteristic of local data in partitioned updates. However, their additional knowledge exchange transmitted through the network degrades the communication efficiency of FL. We propose to reduce the size of knowledge exchanged in these FL setups by clustering and selecting only the most representative bits of information from the clients. The partitioned global update adopted in our work splits the global deep neural network into a lower part for generic feature extraction and an upper part that is more sensitive to this selected client knowledge. Our experiments show that only 1.6% of the initially exchanged data can effectively transfer the characteristic of the client data to the global model in our FL approach, using split networks. These preliminary results evolve our understanding of federated learning by demonstrating efficient training using strategically selected training samples.


title: Upper and lower bounds for Dunkl heat kernel

abstract: On $\mathbb R^N$ equipped with a normalized root system $R$, a multiplicity function $k(\alpha) > 0$, and the associated measure $$ dw(\mathbf x)=\prod_{\alpha\in R}|\langle \mathbf x,\alpha\rangle|^{k(\alpha)}\, d\mathbf x, $$ let $h_t(\mathbf x,\mathbf y)$ denote the heat kernel of the semigroup generated by the Dunkl Laplace operator $\Delta_k$. Let $d(\mathbf x,\mathbf y)=\min_{\sigma\in G} \| \mathbf x-\sigma(\mathbf y)\|$, where $G$ is the reflection group associated with $R$. We derive the following upper and lower bounds for $h_t(\mathbf x,\mathbf y)$: for all $c_l>1/4$ and $0<c_u<1/4$ there are constants $C_l,C_u>0$ such that $$ C_{l}w(B(\mathbf{x},\sqrt{t}))^{-1}e^{-c_{l}\frac{d(\mathbf{x},\mathbf{y})^2}{t}} \Lambda(\mathbf x,\mathbf y,t) \leq h_t(\mathbf{x},\mathbf{y}) \leq C_{u}w(B(\mathbf{x},\sqrt{t}))^{-1}e^{-c_{u}\frac{d(\mathbf{x},\mathbf{y})^2}{t}} \Lambda(\mathbf x,\mathbf y,t), $$ where $\Lambda(\mathbf x,\mathbf y,t)$ can be expressed by means of some rational functions of $\| \mathbf x-\sigma(\mathbf y)\|/\sqrt{t}$. An exact formula for $\Lambda(\mathbf x,\mathbf y,t)$ is provided.


title: SocialVec: Social Entity Embeddings

abstract: This paper introduces SocialVec, a general framework for eliciting social world knowledge from social networks, and applies this framework to Twitter. SocialVec learns low-dimensional embeddings of popular accounts, which represent entities of general interest, based on their co-occurrences patterns within the accounts followed by individual users, thus modeling entity similarity in socio-demographic terms. Similar to word embeddings, which facilitate tasks that involve text processing, we expect social entity embeddings to benefit tasks of social flavor. We have learned social embeddings for roughly 200,000 popular accounts from a sample of the Twitter network that includes more than 1.3 million users and the accounts that they follow, and evaluate the resulting embeddings on two different tasks. The first task involves the automatic inference of personal traits of users from their social media profiles. In another study, we exploit SocialVec embeddings for gauging the political bias of news sources in Twitter. In both cases, we prove SocialVec embeddings to be advantageous compared with existing entity embedding schemes. We will make the SocialVec entity embeddings publicly available to support further exploration of social world knowledge as reflected in Twitter.


title: Rigorous theory of thin-vapor-layer linear optical properties: The case of quenching of atomic polarization upon collisions of atoms with dielectric walls

abstract: Hot alkali metal vapors enclosed in sub-micron spectroscopic cells provide an ideal system for fundamental studies of the atom-wall and atom-light interactions at nanoscale. Here, we propose a novel approach for calculating the eigenmodes of a thin-vapor-layer beyond the limitations of the first-order perturbation theory in optical density for the case of quenching of atomic polarization upon collisions of atoms with dielectric walls. We show that higher-order optical density corrections lead to a remarkable density-dependent blue shift and deformation of the spectral line shapes of reflection, transmission, and absorption. We also demonstrate that the eigenmodes of the thin-vapor-layer can be calculated independently of the choice of optical boundary conditions. This greatly extends the applicability of the constructed theory for the development of miniature atomic sensors.


title: Solving the Class Imbalance Problem Using a Counterfactual Method for Data Augmentation

abstract: Learning from class imbalanced datasets poses challenges for many machine learning algorithms. Many real-world domains are, by definition, class imbalanced by virtue of having a majority class that naturally has many more instances than its minority class (e.g. genuine bank transactions occur much more often than fraudulent ones). Many methods have been proposed to solve the class imbalance problem, among the most popular being oversampling techniques (such as SMOTE). These methods generate synthetic instances in the minority class, to balance the dataset, performing data augmentations that improve the performance of predictive machine learning (ML) models. In this paper we advance a novel data augmentation method (adapted from eXplainable AI), that generates synthetic, counterfactual instances in the minority class. Unlike other oversampling techniques, this method adaptively combines exist-ing instances from the dataset, using actual feature-values rather than interpolating values between instances. Several experiments using four different classifiers and 25 datasets are reported, which show that this Counterfactual Augmentation method (CFA) generates useful synthetic data points in the minority class. The experiments also show that CFA is competitive with many other oversampling methods many of which are variants of SMOTE. The basis for CFAs performance is discussed, along with the conditions under which it is likely to perform better or worse in future tests.


title: Snapshot Ptychography on Array cameras

abstract: We use convolutional neural networks to recover images optically down-sampled by $6.7\times$ using coherent aperture synthesis over a 16 camera array. Where conventional ptychography relies on scanning and oversampling, here we apply decompressive neural estimation to recover full resolution image from a single snapshot, although as shown in simulation multiple snapshots can be used to improve SNR. In place training on experimental measurements eliminates the need to directly calibrate the measurement system. We also present simulations of diverse array camera sampling strategies to explore how snapshot compressive systems might be optimized.


title: Phonon-assisted luminescence in defect centers from many-body perturbation theory

abstract: Phonon-assisted luminescence is a key property of defect centers in semiconductors, and can be measured to perform the readout of the information stored in a quantum bit, or to detect temperature variations. The investigation of phonon-assisted luminescence usually employs phenomenological models, such as that of Huang and Rhys, with restrictive assumptions that can fail to be predictive. In this work, we predict luminescence and study exciton-phonon couplings within a rigorous many-body perturbation theory framework, an analysis that has never been performed for defect centers. In particular, we study the optical emission of the negatively-charged boron vacancy in 2D hexagonal boron nitride, which currently stands out among defect centers in 2D materials thanks to its promise for applications in quantum information and quantum sensing. We show that phonons are responsible for the observed luminescence, which otherwise would be dark due to symmetry. We also show that the symmetry breaking induced by the static Jahn-Teller effect is not able to describe the presence of the experimentally observed peak at 1.5 eV.


title: S-multi-SNE: Semi-Supervised Classification and Visualisation of Multi-View Data

abstract: An increasing number of multi-view data are being published by studies in several fields. This type of data corresponds to multiple data-views, each representing a different aspect of the same set of samples. We have recently proposed multi-SNE, an extension of t-SNE, that produces a single visualisation of multi-view data. The multi-SNE approach provides low-dimensional embeddings of the samples, produced by being updated iteratively through the different data-views. Here, we further extend multi-SNE to a semi-supervised approach, that classifies unlabelled samples by regarding the labelling information as an extra data-view. We look deeper into the performance, limitations and strengths of multi-SNE and its extension, S-multi-SNE, by applying the two methods on various multi-view datasets with different challenges. We show that by including the labelling information, the projection of the samples improves drastically and it is accompanied by a strong classification performance.


title: Properties of Navier-Stokes mild solutions with initial data in subcritical Lorentz spaces

abstract: For initial data $f$ in a subcritical Lorentz space $L^{p,q}(\mathbb{R}^{n}) \hookrightarrow \dot B^{-\frac np}_{\infty,\infty}(\mathbb{R}^n)$ ($n<p<\infty$, $1\leq q \leq \infty$), we prove results which imply in particular that a local in time mild Navier-Stokes solution cannot become unbounded in the $L^{p,q}(\mathbb{R}^{n})$-norm before it becomes unbounded in the norm of the larger subcritical Besov space $\dot B^{-\frac np}_{\infty,\infty}(\mathbb{R}^n)$. In view of the known local theory in such large Besov spaces, this can be thought of as a `propagation of regularity' type of result; here, we provide a self-contained local theory (including scaling-appropriate `blow-up estimates', similar to those established by J. Leray in the Lebesgue setting) in the Lorentz setting, along with uniqueness results which imply such propagation of regularity. Our existence results are based on the method of T. Kato (1984) with Lorentz spaces replacing Lebesgue spaces throughout, and are given without any reference to the Besov framework. The uniqueness results are similarly self-contained, extending certain Lebesgue space methods of Fabes-Jones-Riviere (1972) to the Lorentz setting. We also establish certain continuity properties of the solutions which are constructed. (A more detailed abstract is provided in the article itself.)


title: Multiple $SO(5)$ isovector pairing and seniority $Sp(2\Omega)$ multi-$j$ algebras with isospin

abstract: With nucleons occupying several shell model $j$ orbits, the isovector pair creation operator $A^1_\mu$ (creates a two particle state with angular momentum $J=0$ and isospin $T=1$) is no longer unique. Choosing it to be a sum of single-$j$ isovector pair creation operators each with a phase, there will be multiple pair $SO(5)$ algebras with isospin; with $r$ number of $j$ orbits, there will be $2^{r-1}$ $SO(5)$ algebras each with a corresponding complementary $Sp(2\Omega)$ algebra [$2\Omega = \sum_j (2j+1)$] that gives seniority and reduced isospin quantum numbers. Three applications of multiple $SO(5)$ algebras are presented demonstrating the usefulness of considering $SO(5)$ pairing algebras with general sign factors.


title: Semantically Consistent Image-to-Image Translation for Unsupervised Domain Adaptation

abstract: Unsupervised Domain Adaptation (UDA) aims to adapt models trained on a source domain to a new target domain where no labelled data is available. In this work, we investigate the problem of UDA from a synthetic computer-generated domain to a similar but real-world domain for learning semantic segmentation. We propose a semantically consistent image-to-image translation method in combination with a consistency regularisation method for UDA. We overcome previous limitations on transferring synthetic images to real looking images. We leverage pseudo-labels in order to learn a generative image-to-image translation model that receives additional feedback from semantic labels on both domains. Our method outperforms state-of-the-art methods that combine image-to-image translation and semi-supervised learning on relevant domain adaptation benchmarks, i.e., on GTA5 to Cityscapes and SYNTHIA to Cityscapes.


title: On the relation between Stratonovich and Ito integrals with functional integrands of conditional measure flows

abstract: In this small note we explicit the relation between Ito and Stratonovich integrals when conditional measure flow components are present in the integrands. The `correction' term involves Lions-type measure derivatives and clarifies which cross-correlations need to be taken into account. We cast the framework in relation to SDEs of mean-field type depending on conditional flows of measure. The result being trivial under full flows of measure.


title: A Data-driven Approach to Neural Architecture Search Initialization

abstract: Algorithmic design in neural architecture search (NAS) has received a lot of attention, aiming to improve performance and reduce computational cost. Despite the great advances made, few authors have proposed to tailor initialization techniques for NAS. However, literature shows that a good initial set of solutions facilitate finding the optima. Therefore, in this study, we propose a data-driven technique to initialize a population-based NAS algorithm. Particularly, we proposed a two-step methodology. First, we perform a calibrated clustering analysis of the search space, and second, we extract the centroids and use them to initialize a NAS algorithm. We benchmark our proposed approach against random and Latin hypercube sampling initialization using three population-based algorithms, namely a genetic algorithm, evolutionary algorithm, and aging evolution, on CIFAR-10. More specifically, we use NAS-Bench-101 to leverage the availability of NAS benchmarks. The results show that compared to random and Latin hypercube sampling, the proposed initialization technique enables achieving significant long-term improvements for two of the search baselines, and sometimes in various search scenarios (various training budgets). Moreover, we analyze the distributions of solutions obtained and find that that the population provided by the data-driven initialization technique enables retrieving local optima (maxima) of high fitness and similar configurations.


title: Image of the phonon spectrum in the 1/f noise of topological insulators

abstract: As reported recently, the 1/f noise intensity in both (Bi,Sb)2Te3 [Islam et al., Appl. Phys. Lett. 111, 062107 (2017)] and BiSbTeSe1.6 [Biswas et al., Appl. Phys. Lett. 115, 131601 (2019)] features noise peaks which develop at some specific temperatures. We compared this noise structure with either phonon density of states or Raman spectrum of each topological insulator (TI), respectively. In (BiSb)2Te3, the comparison revealed that the noise peaks track the van Hove singularities in the phonon density of states. It resulted that bulk atomic oscillators are responsible for the noise peaks. The most intense noise peak observed in (BiSb)2Te3 at 50 K is attributed to the thermal motion of the Bi atoms. Other less intense noise peaks are assigned to either a single phonon mode or multi-phonon combinations. We found that thermal motions of Bi and Te2 atoms in different symmetry directions are involved in most of the phonon combinations, which stand for the signature of the lattice anharmonicity. The noise increase observed in both (Bi,Sb)2Te3 and BiSbTeSe1.6 above a specific temperature threshold is attributed to the anharmonicity-induced strengthening of the carrier-phonon coupling. In the case of BiSbTeSe1.6, we show that all noise singularities are mirrored in the Raman spectrum of a structurally close TI (BiSbTeSe2) in the whole temperature range. This indicates that although transport can be at the surface or in the bulk or both of them, the carrier-phonon interaction is the only source of 1/f fluctuations in TIs. Inherently, these results imply that microscopic origin of 1/f noise in solid is in the perpetual thermal motion of the atoms.


title: Normal limiting distributions for systems of linear equations in random sets

abstract: We consider the binomial random set model $[n]_p$ where each element in $\{1,\dots,n\}$ is chosen independently with probability $p:=p(n)$. We show that for essentially all regimes of $p$ and very general conditions for a matrix $A$ and a column vector $\mathbf{b}$, the count of specific integer solutions to the system of linear equations $A\mathbf{x} = \mathbf{b}$ with the entries of $\mathbf{x}$ in $[n]_p$ follows a (conveniently rescaled) normal limiting distribution. This applies among others to the number of solutions with every variable having a different value, as well as to a broader class of so-called non-trivial solutions in homogeneous strictly balanced systems. Our proof relies on the delicate linear algebraic study both of the subjacent matrices and the corresponding ranks of certain submatrices, together with the application of the method of moments in probability theory.


title: A Fresh Look at AGB stars in Galactic Open Clusters with Gaia: Impact on Stellar Models and the Initial-Final Mass Relation

abstract: Benefiting from the GAIA second and early third releases of photometric and astrometric data we examine the population of asymptotic giant branch (AGB) stars that appear in the fields of intermediate-age and young open star clusters. We identify 49 AGB star candidates, brighter than the tip of the red giant branch, with a good-to-high cluster membership probability. Among them we find 19 TP-AGB stars with known spectral type: 4 M stars, 3 MS/S stars and 12 C stars. By combining observations, stellar models, and radiative transfer calculations that include the effect of circumstellar dust, we characterize each star in terms of initial mass, luminosity, mass-loss rate, core mass, period and mode of pulsation. The information collected helps us shed light on the TP-AGB evolution at solar-like metallicity, placing constraints on the third dredge-up process, the initial masses of carbon stars, stellar winds, and the initial-final mass relation (IFMR). In particular, we find that two bright carbon stars, MSB 75 and BM IV 90, members of the clusters NGC 7789 and NGC 2660 (with similar ages of $\simeq 1.2-1.6$ Gyr and initial masses $ 2.1 \ge M_{\rm i}/M_{\odot} \ge 1.9$), have unusually high core masses, $M_{\rm c} \approx 0.67-0.7\,M_{\odot}$. These results support the findings of a recent work (Marigo et al. 2020) that identified a kink in the IFMR, which interrupts its monotonic trend just at the same initial masses. Finally, we investigate two competing scenarios to explain the $M_{\rm c}$ data: the role of stellar winds in single-star evolution, and binary interactions through the blue-straggler channel.


title: New Streaming Algorithms for High Dimensional EMD and MST

abstract: We study streaming algorithms for two fundamental geometric problems: computing the cost of a Minimum Spanning Tree (MST) of an $n$-point set $X \subset \{1,2,\dots,\Delta\}^d$, and computing the Earth Mover Distance (EMD) between two multi-sets $A,B \subset \{1,2,\dots,\Delta\}^d$ of size $n$. We consider the turnstile model, where points can be added and removed. We give a one-pass streaming algorithm for MST and a two-pass streaming algorithm for EMD, both achieving an approximation factor of $\tilde{O}(\log n)$ and using polylog$(n,d,\Delta)$-space only. Furthermore, our algorithm for EMD can be compressed to a single pass with a small additive error. Previously, the best known sublinear-space streaming algorithms for either problem achieved an approximation of $O(\min\{ \log n , \log (\Delta d)\} \log n)$ [Andoni-Indyk-Krauthgamer '08, Backurs-Dong-Indyk-Razenshteyn-Wagner '20]. For MST, we also prove that any constant space streaming algorithm can only achieve an approximation of $\Omega(\log n)$, analogous to the $\Omega(\log n)$ lower bound for EMD of [Andoni-Indyk-Krauthgamer '08]. Our algorithms are based on an improved analysis of a recursive space partitioning method known generically as the Quadtree. Specifically, we show that the Quadtree achieves an $\tilde{O}(\log n)$ approximation for both EMD and MST, improving on the $O(\min\{ \log n , \log (\Delta d)\} \log n)$ approximation of [Andoni-Indyk-Krauthgamer '08, Backurs-Dong-Indyk-Razenshteyn-Wagner '20].


title: Traveling waves near Couette flow for the 2D Euler equation

abstract: In this paper we reveal the existence of a large family of new, nontrivial and smooth traveling waves for the 2D Euler equation at an arbitrarily small distance from the Couette flow in $H^s$, with $s<3/2$, at the level of the vorticity. The speed of these waves is of order 1 with respect to this distance. This result strongly contrasts with the setting of very high regularity in Gevrey spaces (see arXiv:1306.5028), where the problem exhibits an inviscid damping mechanism that leads to relaxation of perturbations back to nearby shear flows. It also complements the fact that there not exist nontrivial traveling waves in the $H^{\frac{3}{2}+}$ neighborhoods of Couette flow (see arXiv:1004.5149).


title: Second-generation dust in planetary systems: The case of HD 163296

abstract: Observations indicate that large, dust-laden protoplanetary discs are common. Some features, like gaps, rings and spirals, suggest they may host young planets, which can excite the orbits of nearby leftover planetesimals. Energetic collisions among these bodies can lead to the production of second-generation dust. Grains produced by collisions may have a dynamical behaviour different from that of first-generation, primordial dust out of which planetesimals and planets formed. We aim to study these differences for the HD 163296 system and determine whether dynamical signatures in the mixture of the two dust populations can help separate their contributions. We use three-dimensional (3-D) hydrodynamic models to describe the gaseous disc with three, Saturn- to Jupiter-mass, embedded planets. Dust grains, of sizes 1um-1mm, are treated as Lagrangean particles with resolved thermodynamics and mass loss. Initial disc and planet configurations are derived from observation-based work, which indicates low gas viscosity. The 3-D approach also allows us to detect the formation of vortices induced by Rossby waves, where dust becomes concentrated and may contribute to planetesimal formation. We find that the main differences in the dynamical behaviour of first- and second-generation dust occur in the vertical distribution. The two populations have similar distributions around the disc mid-plane, although second-generation dust shows longer residence times close to the radial locations of the planets' gas gaps. Sedimentation rates of um-size grains are comparable to or lower than the production rates by planetesimals' collisions, making this population potentially observable. These outcomes can be extended to similar systems harbouring giant planets.


title: Multigraded Castelnuovo-Mumford regularity via Klyachko filtrations

abstract: In this paper, we consider $\mathbb{Z}^{r}-$graded modules on the $\mathrm{Cl}(X)$ $-$graded Cox ring $\mathbb{C}[x_{1},\dotsc,x_{r}]$ of a smooth complete toric variety $X$. Using the theory of Klyachko filtrations in the reflexive case, we construct a collection of lattice polytopes codifying the multigraded Hilbert function of the module. We apply this approach to reflexive $\mathbb{Z}^{s+r+2}$-graded modules over non-standard bigraded polynomial rings $\mathbb{C}[x_{0},\dotsc,x_{s},y_{0},\dotsc,y_{r}]$. In this case, we give sharp bounds for the multigraded regularity index of their multigraded Hilbert function, and a method to compute their Hilbert polynomial.


title: A Retrospective Analysis using Deep-Learning Models for Prediction of Survival Outcome and Benefit of Adjuvant Chemotherapy in Stage II/III Colorectal Cancer

abstract: Most early-stage colorectal cancer (CRC) patients can be cured by surgery alone, and only certain high-risk early-stage CRC patients benefit from adjuvant chemotherapies. However, very few validated biomarkers are available to accurately predict survival benefit from postoperative chemotherapy. We developed a novel deep-learning algorithm (CRCNet) using whole-slide images from Molecular and Cellular Oncology (MCO) to predict survival benefit of adjuvant chemotherapy in stage II/III CRC. We validated CRCNet both internally through cross-validation and externally using an independent cohort from The Cancer Genome Atlas (TCGA). We showed that CRCNet can accurately predict not only survival prognosis but also the treatment effect of adjuvant chemotherapy. The CRCNet identified high-risk subgroup benefits from adjuvant chemotherapy most and significant longer survival is observed among chemo-treated patients. Conversely, minimal chemotherapy benefit is observed in the CRCNet low- and medium-risk subgroups. Therefore, CRCNet can potentially be of great use in guiding treatments for Stage II/III CRC.


title: An Analysis of Elephants' Movement Data in Sub-Saharan Africa Using Clustering

abstract: Understanding the movement of animals is crucial to conservation efforts. Past research often focuses on factors affecting movement, rather than locations of interest that animals return to or habitat. We explore the use of clustering to identify locations of interest to African Elephants in regions of Sub-Saharan Africa. Our analysis was performed using publicly available datasets for tracking African elephants at Kruger National Park (KNP), South Africa; Etosha National Park, Namibia; as well as areas in Burkina Faso and the Congo. Using the DBSCAN and KMeans clustering algorithms, we calculate clusters and centroids to simplify elephant movement data and highlight important locations of interest. Through a comparison of feature spaces with and without temperature, we show that temperature is an important feature to explain movement clustering. Recognizing the importance of temperature, we develop a technique to add external temperature data from an API to other geospatial datasets that would otherwise not have temperature data. After addressing the hurdles of using external data with marginally different timestamps, we consider the quality of this data, and the quality of the centroids of the clusters calculated based on this external temperature data. Finally, we overlay these centroids onto satellite imagery and locations of human settlements to validate the real-life applications of the calculated centroids to identify locations of interest for elephants. As expected, we confirmed that elephants tend to cluster their movement around sources of water as well as some human settlements, especially those with water holes. Identifying key locations of interest for elephants is beneficial in predicting the movement of elephants and preventing poaching. These methods may in the future be applied to other animals beyond elephants to identify locations of interests for them.


title: Learning Formulas in Finite Variable Logics

abstract: We consider grammar-restricted exact learning of formulas and terms in finite variable logics. We propose a novel and versatile automata-theoretic technique for solving such problems. We first show results for learning formulas that classify a set of positively- and negatively-labeled structures. We give algorithms for realizability and synthesis of such formulas along with upper and lower bounds. We also establish positive results using our technique for other logics and variants of the learning problem, including first-order logic with least fixed point definitions, higher-order logics, and synthesis of queries and terms with recursively-defined functions.


title: Multivariable Super-Twisting Algorithm for Systems with Uncertain Input Matrix and Perturbations

abstract: This paper proposes a Lyapunov approach to the design of a multivariable generalized Super-Twisting algorithm (MGSTA), which is able to control a system with perturbations and uncertain control matrix, both depending on time and the system states. The presented procedure shows that, under reasonable assumptions for the uncertainties, it is always possible to find a set of constant gains for the MGSTA in order to ensure global and robust finite-time stability of the system's outputs. Simulation results on an omnidirectional mobile robot illustrate the performance of the MGSTA.


title: A Unified Game-Theoretic Interpretation of Adversarial Robustness

abstract: This paper provides a unified view to explain different adversarial attacks and defense methods, \emph{i.e.} the view of multi-order interactions between input variables of DNNs. Based on the multi-order interaction, we discover that adversarial attacks mainly affect high-order interactions to fool the DNN. Furthermore, we find that the robustness of adversarially trained DNNs comes from category-specific low-order interactions. Our findings provide a potential method to unify adversarial perturbations and robustness, which can explain the existing defense methods in a principle way. Besides, our findings also make a revision of previous inaccurate understanding of the shape bias of adversarially learned features.


title: A practical analysis of ROP attacks

abstract: Control Flow Hijacking attacks have posed a serious threat to the security of applications for a long time where an attacker can damage the control Flow Integrity of the program and execute arbitrary code. These attacks can be performed by injecting code in the program's memory or reusing already existing code in the program (also known as Code-Reuse Attacks). Code-Reuse Attacks in the form of Return-into-libc Attacks or Return-Oriented Programming Attacks are said to be Turing Complete, providing a guarantee that there will always exist code segments (also called ROP gadgets) within a binary allowing an attacker to perform any kind of function by building a suitable ROP chain (chain of ROP gadgets). Our goal is to study different techniques of performing ROP Attacks and find the difficulties encountered to perform such attacks. For this purpose, we have designed an automated tool which works on 64-bit systems and generates a ROP chain from ROP gadgets to execute arbitrary system calls.


title: Structural analysis of disordered dimer packings

abstract: Jammed disordered packings of non-spherical particles show significant variation in the packing density as a function of particle shape for a given packing protocol. Rotationally symmetric elongated shapes such as ellipsoids, spherocylinders, and dimers, e.g., pack significantly denser than spheres over a narrow range of aspect ratios, exhibiting a characteristic peak at aspect ratios of $\alpha_{\rm max}\approx 1.4-1.5$. However, the structural features that underlie this non-monotonic behaviour in the packing density are unknown. Here, we study disordered packings of frictionless dimers in three dimensions generated by a gravitational pouring protocol in LAMMPS. Focusing on the characteristics of contacts as well as orientational and translational order metrics, we identify a number of structural features that accompany the formation of maximally dense packings as the dimer aspect ratio $\alpha$ is varied from the spherical limit. Our results highlight that dimer packings undergo significant structural changes as $\alpha$ increases up to $\alpha_{\rm max}$ manifest in the reorganisation of the contact configurations between neighbouring dimers, increasing nematic order, and decreasing local translational order. Remarkably, for $\alpha>\alpha_{\rm max}$ our metrics remain largely unchanged, indicating that the peak in the packing density is related to the interplay of structural rearrangements for $\alpha<\alpha_{\rm max}$ and subsequent excluded volume effects with unchanged structure for $\alpha>\alpha_{\rm max}$.

categories: cond-mat.soft cond-mat.dis-nn cond-mat.mtrl-sci cond-mat.stat-mech


title: Optimal Inverted Landing in a Small Aerial Robot with Varied Approach Velocities and Landing Gear Designs

abstract: Inverted landing is a challenging feat to perform in aerial robots, especially without external positioning. However, it is routinely performed by biological fliers such as bees, flies, and bats. Our previous observations of landing behaviors in flies suggest an open-loop causal relationship between their putative visual cues and the kinematics of the aerial maneuvers executed. For example, the degree of rotational maneuver (the amount of body inversion prior to touchdown) and the amount of leg-assisted body swing both depend on the flies' initial body states while approaching the ceiling. In this work, inspired by the inverted landing behavior of flies, we used a physics-based simulation with experimental validation to systematically investigate how optimized inverted landing maneuvers depend on the initial approach velocities with varied magnitude and direction. This was done by analyzing the putative visual cues (that can be derived from onboard measurements) during optimal maneuvering trajectories. We identified a three-dimensional policy region, from which a mapping to a global inverted landing policy can be developed without the use of external positioning data. Through simulation, we also investigated the effects of an array of landing gear designs on the optimized landing performance and identified their advantages and disadvantages. The above results have been partially validated using limited experimental testing and will continue to inform and guide our future experiments, for example by applying the calculated global policy.


title: Sharp ill-posedness for the generalized Camassa-Holm equation in Besov spaces

abstract: In this paper, we consider the Cauchy problem for the generalized Camassa-Holm equation that includes the Camassa-Holm as well as the Novikov equation on the line. We present a new and unified method to prove the sharp ill-posedness for the generalized Camassa-Holm equation in $B^s_{p,\infty}$ with $s>\max\{1+1/p, 3/2\}$ and $1\leq p\leq\infty$ in the sense that the solution map to this equation starting from $u_0$ is discontinuous at $t = 0$ in the metric of $B^s_{p,\infty}$. Our results cover and improve the previous work given in [J. Li, Y. Yu, W. Zhu, Ill-posedness for the Camassa-Holm and related equations in Besov spaces, J. Differential Equations, 306 (2022), 403--417], solving an open problem left in [J. Li, Y. Yu, W. Zhu, Ill-posedness for the Camassa-Holm and related equations in Besov spaces, J. Differential Equations, 306 (2022), 403--417].


title: MetaFEM: A Generic FEM Solver By Meta-expressions

abstract: Current multi-physics Finite Element Method (FEM) solvers are complex systems in terms of both their mathematical complexity and lines of code. This paper proposes a skeleton generic FEM solver, named MetaFEM, in total about 5,000 lines of Julia code, which translates generic input Partial Differential Equation (PDE) weak forms into corresponding GPU-accelerated simulations with a grammar similar to FEniCS or FreeFEM. Two novel approaches differentiate MetaFEM from the common solvers: (1) the FEM kernel is based on an original theory/algorithm which explicitly processes meta-expressions, as the name suggests, and (2) the symbolic engine is a rule-based Computer Algebra System (CAS), i.e., the equations are rewritten/derived according to a set of rewriting rules instead of going through completely fixed routines, supporting easy customization by developers. Example cases in thermal conduction, linear elasticity and incompressible flow are presented to demonstrate utility.


title: Time-resolved velocity and pressure field quantification in a flow-focusing device for ultrafast microbubble production

abstract: Flow-focusing devices have gained great interest in the past decade, due to their capability to produce monodisperse microbubbles for diagnostic and therapeutic medical ultrasound applications. However, up-scaling production to industrial scale requires a paradigm shift from single chip operation to highly parallelized systems. Parallelization gives rise to fluidic interactions between nozzles that, in turn, may lead to a decreased monodispersity. Here, we study the velocity and pressure field fluctuations in a single flow-focusing nozzle during bubble production. We experimentally quantify the velocity field inside the nozzle at 100 ns time resolution, and a numerical model provides insight into both the oscillatory velocity and pressure fields. Our results demonstrate that, at the length scale of the flow focusing channel, the velocity oscillations propagate at fluid dynamical time scale (order of microseconds) whereas the dominant pressure oscillations are linked to the bubble pinch-off and propagate at a much faster time scale (order of nanoseconds).


title: Empirical analysis of representation learning and exploration in neural kernel bandits

abstract: Neural bandits have been shown to provide an efficient solution to practical sequential decision tasks that have nonlinear reward functions. The main contributor to that success is approximate Bayesian inference, which enables neural network (NN) training with uncertainty estimates. However, Bayesian NNs often suffer from a prohibitive computational overhead or operate on a subset of parameters. Alternatively, certain classes of infinite neural networks were shown to directly correspond to Gaussian processes (GP) with neural kernels (NK). NK-GPs provide accurate uncertainty estimates and can be trained faster than most Bayesian NNs. We propose to guide common bandit policies with NK distributions and show that NK bandits achieve state-of-the-art performance on nonlinear structured data. Moreover, we propose a framework for measuring independently the ability of a bandit algorithm to learn representations and explore, and use it to analyze the impact of NK distributions w.r.t.~those two aspects. We consider policies based on a GP and a Student's t-process (TP). Furthermore, we study practical considerations, such as training frequency and model partitioning. We believe our work will help better understand the impact of utilizing NKs in applied settings.


title: Morphogenesis of street networks. A reaction-diffusion system for self-organized cities

abstract: Urban morphogenesis is the process of formation of its elements and the specialization of its suburbs. Street networks are the structural part of the system. Understand their formation reveals crucial information about urban transformation and dynamics behind its functioning. In this work we focus on the morphogenesis of street networks and we study it through a spatial network generator model. This latter is composed by three layers (a reaction-diffusion layer, a dynamic vector field and a spatial network) surrounded by an environment. The emerging network feeds back to its morphogenetic elements, driving the model to an unexpected behaviour. We applied the model in a real urban context (F\'ecamp town, Normandy, France) and we measure properties of obtained networks.


title: ActFloor-GAN: Activity-Guided Adversarial Networks for Human-Centric Floorplan Design

abstract: We present a novel two-stage approach for automated floorplan design in residential buildings with a given exterior wall boundary. Our approach has the unique advantage of being human-centric, that is, the generated floorplans can be geometrically plausible, as well as topologically reasonable to enhance resident interaction with the environment. From the input boundary, we first synthesize a human-activity map that reflects both the spatial configuration and human-environment interaction in an architectural space. We propose to produce the human-activity map either automatically by a pre-trained generative adversarial network (GAN) model, or semi-automatically by synthesizing it with user manipulation of the furniture. Second, we feed the human-activity map into our deep framework ActFloor-GAN to guide a pixel-wise prediction of room types. We adopt a re-formulated cycle-consistency constraint in ActFloor-GAN to maximize the overall prediction performance, so that we can produce high-quality room layouts that are readily convertible to vectorized floorplans. Experimental results show several benefits of our approach. First, a quantitative comparison with prior methods shows superior performance of leveraging the human-activity map in predicting piecewise room types. Second, a subjective evaluation by architects shows that our results have compelling quality as professionally-designed floorplans and much better than those generated by existing methods in terms of the room layout topology. Last, our approach allows manipulating the furniture placement, considers the human activities in the environment, and enables the incorporation of user-design preferences.


title: Persistence of the Brauer-Manin obstruction on cubic surfaces

abstract: Let $X$ be a cubic surface over a global field $k$. We prove that a Brauer-Manin obstruction to the existence of $k$-points on $X$ will persist over every extension $L/k$ with degree relatively prime to $3$. In other words, a cubic surface has nonempty Brauer set over $k$ if and only if it has nonempty Brauer set over some extension $L/k$ with $3\nmid[L:k]$. Therefore, the conjecture of Colliot-Th\'el\`ene and Sansuc on the sufficiency of the Brauer-Manin obstruction for cubic surfaces implies that $X$ has a $k$-rational point if and only if $X$ has a $0$-cycle of degree $1$. This latter statement is a special case of a conjecture of Cassels and Swinnerton-Dyer.


title: POSHAN: Cardinal POS Pattern Guided Attention for News Headline Incongruence

abstract: Automatic detection of click-bait and incongruent news headlines is crucial to maintaining the reliability of the Web and has raised much research attention. However, most existing methods perform poorly when news headlines contain contextually important cardinal values, such as a quantity or an amount. In this work, we focus on this particular case and propose a neural attention based solution, which uses a novel cardinal Part of Speech (POS) tag pattern based hierarchical attention network, namely POSHAN, to learn effective representations of sentences in a news article. In addition, we investigate a novel cardinal phrase guided attention, which uses word embeddings of the contextually-important cardinal value and neighbouring words. In the experiments conducted on two publicly available datasets, we observe that the proposed methodgives appropriate significance to cardinal values and outperforms all the baselines. An ablation study of POSHAN shows that the cardinal POS-tag pattern-based hierarchical attention is very effective for the cases in which headlines contain cardinal values.


title: Spectrum of p-adic linear differential equations I: The shape of the spectrum

abstract: This paper extends our previous works arXiv:1802.07306 [math.NT], arXiv:1808.02382 [math.NT] on determining the spectrum, in the Berkovich sense, of ultrametric linear differential equations. Our previous works focused on equations with constant coefficients or over a field of formal power series. In this paper, we investigate the spectrum of $p$-adic differential equations at a generic point on a quasi-smooth curve. This analysis allows us to establish a significant connection between the spectrum and the spectral radii of convergence of a differential equation when considering the affine line. Furthermore, the spectrum offers a more detailed decomposition compared to Robba's decomposition based on spectral radii.


title: Interpreting Representation Quality of DNNs for 3D Point Cloud Processing

abstract: In this paper, we evaluate the quality of knowledge representations encoded in deep neural networks (DNNs) for 3D point cloud processing. We propose a method to disentangle the overall model vulnerability into the sensitivity to the rotation, the translation, the scale, and local 3D structures. Besides, we also propose metrics to evaluate the spatial smoothness of encoding 3D structures, and the representation complexity of the DNN. Based on such analysis, experiments expose representation problems with classic DNNs, and explain the utility of the adversarial training.


title: On Slice Isolation Options in the Transport Network and Associated Feasibility Indicators

abstract: Isolation is one of the more relevant attributes associated to the idea of network slicing, introduced by 5G services. Through isolation it is expected that slices from different customers could gracefully coexist without interfering each other, in the sense that whatever misbehavior or unforeseen demand from one slice customer could not affect the communication service received by any other slice customer supported atop the same physical transport infrastructure. This paper surveys and compare different technical approaches that can be taken for providing distinct isolation levels in the transport network, as a major component of end-to-end network slices. Furthermore, a number of isolation feasibility indicators are defined and proposed. These indicators are based on the approaches referred before, as a mean of guiding orchestration decisions at the time of provisioning or reconfiguring the transport slices in the network.


title: On the application of incomplete FWH surfaces for aeroacoustic predictions

abstract: This work is motivated by CFD simulations from a realistic landing gear performed modeling only the half bottom of the aircraft fuselage [15]. Hence, in this previous analysis, the aeroacoustic predictions had to employ incomplete FWH surfaces. In the present study, a discussion on the closed surface requirement is provided for the frequency domain FWH formulation, although we believe that the approach presented here can also be applied for the time-domain methodology. We show that all sources computed on the closed permeable surface may play an important role if the acoustic prediction is sought at observer positions distributed along the entire circular arc. However, results obtained for 2D and 3D model problems show that, if wisely designed, a permeable FWH surface composed only by a finite patch, i.e., an incomplete surface, can be employed to accurately predict the noise at specific observer positions. In this case, the surface must be placed between the line of sight from the sources to the observers, for example, in flyover or sideline locations. In the current applications, results from finite surfaces agree with the expected values as long as the source magnitudes decay on elements distant from the true incident source, i.e., the airframe. Hence, in order for this approach to be valid, the sources cannot be truncated. The planar setup is beneficial in airframe noise applications since it avoids contamination by quadrupole sources crossing the boundaries and also may simplify the surface design.


title: SmartDepthSync: Open Source Synchronized Video Recording System of Smartphone RGB and Depth Camera Range Image Frames with Sub-millisecond Precision

abstract: Nowadays, smartphones can produce a synchronized (synced) stream of high-quality data, including RGB images, inertial measurements, and other data. Therefore, smartphones are becoming appealing sensor systems in the robotics community. Unfortunately, there is still the need for external supporting sensing hardware, such as a depth camera precisely synced with the smartphone sensors. In this paper, we propose a hardware-software recording system that presents a heterogeneous structure and contains a smartphone and an external depth camera for recording visual, depth, and inertial data that are mutually synchronized. The system is synced at the time and the frame levels: every RGB image frame from the smartphone camera is exposed at the same moment of time with a depth camera frame with sub-millisecond precision. We provide a method and a tool for sync performance evaluation that can be applied to any pair of depth and RGB cameras. Our system could be replicated, modified, or extended by employing our open-sourced materials.


title: Multiplicative-Accumulative matching of NLO calculations with parton showers

abstract: We propose a new approach for combining next-to-leading order (NLO) and parton shower (PS) calculations so as to obtain three core features: (a) applicability to general showers, as with the MCatNLO and POWHEG methods; (b) positive-weight events, as with the KrkNLO and POWHEG methods; and (c) all showering attributed to the parton shower code, as with the MCatNLO and KrkNLO methods. This is achieved by using multiplicative matching in phase space regions where the shower overestimates the matrix element and accumulative (additive) matching in regions where the shower underestimates the matrix element, an approach that can be viewed as a combination of the MCatNLO and KrkNLO methods.


title: Messy or Ordered? Multi-scale Mechanics DictatesShape-Morphing of Hierarchical 2D Fiber-Networks

abstract: Shape-morphing networks of mesoscale filaments are a common hierarchical feature in biology and hold significant potential for a range of technological applications, from micro-muscles to shape-morphing optical devices. Here, we demonstrate both experimentally and computationally that the shape-morphing of highly-ordered 2D networks constructed of thermoresponsive mesoscale polymeric fibers strongly depends on the physical attributes of the single fiber, in particular on its diameter, as well as on the network's density. We show that based on these parameters, such fiber-networks exhibit one of two thermally driven morphing behaviors: (i) the fibers stay straight, and the network preserves its ordered morphology, exhibiting a bulk-like behavior; or (ii) the fibers buckle and the network becomes messy and highly disordered. Notably, in both cases, the networks display memory and regain their original ordered morphology upon shrinking. This hierarchical principle, demonstrated here on a range of networks, offers a new way for controlling the shape morphing of materials with mesoscale resolutions and elucidates that minute changes in the mesoscale structural attributes can translate to a dramatic change in the morphing behaviors at the macroscale.


title: AUTOKD: Automatic Knowledge Distillation Into A Student Architecture Family

abstract: State-of-the-art results in deep learning have been improving steadily, in good part due to the use of larger models. However, widespread use is constrained by device hardware limitations, resulting in a substantial performance gap between state-of-the-art models and those that can be effectively deployed on small devices. While Knowledge Distillation (KD) theoretically enables small student models to emulate larger teacher models, in practice selecting a good student architecture requires considerable human expertise. Neural Architecture Search (NAS) appears as a natural solution to this problem but most approaches can be inefficient, as most of the computation is spent comparing architectures sampled from the same distribution, with negligible differences in performance. In this paper, we propose to instead search for a family of student architectures sharing the property of being good at learning from a given teacher. Our approach AutoKD, powered by Bayesian Optimization, explores a flexible graph-based search space, enabling us to automatically learn the optimal student architecture distribution and KD parameters, while being 20x more sample efficient compared to existing state-of-the-art. We evaluate our method on 3 datasets; on large images specifically, we reach the teacher performance while using 3x less memory and 10x less parameters. Finally, while AutoKD uses the traditional KD loss, it outperforms more advanced KD variants using hand-designed students.


title: Novel Approach for Evaluating Detector-Related Uncertainties in a LArTPC Using MicroBooNE Data

abstract: Primary challenges for current and future precision neutrino experiments using liquid argon time projection chambers (LArTPCs) include understanding detector effects and quantifying the associated systematic uncertainties. This paper presents a novel technique for assessing and propagating LArTPC detector-related systematic uncertainties. The technique makes modifications to simulation waveforms based on a parameterization of observed differences in ionization signals from the TPC between data and simulation, while remaining insensitive to the details of the detector model. The modifications are then used to quantify the systematic differences in low- and high-level reconstructed quantities. This approach could be applied to future LArTPC detectors, such as those used in SBN and DUNE.


title: Reaching for the Edge I: Probing the Outskirts of Massive Galaxies with HSC, DECaLS, SDSS, and Dragonfly

abstract: The outer light (stellar halos) of massive galaxies has recently emerged as a possible low scatter tracer of dark matter halo mass. To test the robustness of outer light measurements across different data sets, we compare the surface brightness profiles of massive galaxies using four independent data sets: the Hyper Suprime-Cam survey (HSC), the Dark Energy Camera Legacy Survey (DECaLS), the Sloan Digital Sky Survey (SDSS), and the Dragonfly Wide Field Survey (Dragonfly). We use customized pipelines for HSC and DECaLS to achieve better sky background subtraction. For galaxies at $z<0.05$, Dragonfly has the best control of systematics, reaching surface brightness levels of $\mu_r \sim 30$ mag/arcsec$^{2}$. At $0.19<z<0.50$, HSC can reliably recover surface brightness profiles to $\mu_{r} \sim 28.5$ mag/arcsec$^{2}$ reaching $R=100 - 150$ kpc. DECaLS surface brightness profiles show good agreement with HSC but are noisier at large radii. The median profiles of galaxy ensembles in both HSC and DECaLS reach $R > 200$ kpc without significant bias. At $0.19<z<0.50$, DECaLS and HSC measurements of the stellar mass contained within 100 kpc agree within 0.05 dex. Finally, we use weak gravitational lensing to show that measurements of outer light with DECaLS at $0.19<z<0.50$ show a similar promise as HSC as a low scatter proxy of halo mass. The tests and results from this paper represent an important step forward for accurate measurements of the outer light of massive galaxies and demonstrate that outer light measurements from DECam imaging will be a promising method for finding galaxy clusters for DES and DESI.


title: The Relation between Morphological Asymmetry and Nuclear Activity in Low-redshift Galaxies

abstract: The morphology of galaxies reflects their assembly history and ongoing dynamical perturbations from the environment. Analyzing i-band images from the Pan-STARRS1 Survey, we study the optical morphological asymmetry of the host galaxies of a large, well-defined sample of nearby active galactic nuclei (AGNs) to investigate the role of mergers and interactions in triggering nuclear activity. The AGNs, comprising 245 type 1 and 4514 type 2 objects, are compared with 4537 star-forming galaxies matched in redshift and stellar mass. We develop a comprehensive masking strategy to isolate the emission of the target from foreground stars and other contaminating sources, all the while retaining projected companions of comparable brightness that may be major mergers. Among three variants of nonparametric indices, both the popular CAS asymmetry parameter and the outer asymmetry parameter (A_outer) yield robust measures of morphological distortion for star-forming galaxies and type 2 AGNs, while only A_outer is effective for type 1 AGNs. The shape asymmetry, by comparison, is affected more adversely by background noise. Asymmetry indices > 0.4 effectively trace systems that are candidate ongoing mergers. Contrary to theoretical expectations, galaxy interactions and mergers are not the main drivers of nuclear activity, at least not in our sample of low-redshift, relatively low-luminosity AGNs, whose host galaxies are significantly less asymmetric than the control sample of star-forming galaxies. Moreover, type 2 AGNs are morphologically indistinguishable from their type 1 counterparts. The level of AGN activity does not correlate with asymmetry, not even among the major merger candidates. As a by-product, we find, consistent with previous studies, that the average asymmetry of star-forming galaxies increases above the main sequence, although not all major mergers exhibit enhanced star formation.


title: Computability and Beltrami fields in Euclidean space

abstract: In this article, we pursue our investigation of the connections between the theory of computation and hydrodynamics. We prove the existence of stationary solutions of the Euler equations in Euclidean space, of Beltrami type, that can simulate a universal Turing machine. In particular, these solutions possess undecidable trajectories. Heretofore, the known Turing complete constructions of steady Euler flows in dimension 3 or higher were not associated to a prescribed metric. Our solutions do not have finite energy, and their construction makes crucial use of the non-compactness of $\mathbb R^3$, however they can be employed to show that an arbitrary tape-bounded Turing machine can be robustly simulated by a Beltrami flow on $\mathbb T^3$ (with the standard flat metric). This shows that there exist steady solutions to the Euler equations on the flat torus exhibiting dynamical phenomena of (robust) computational complexity as high as desired. We also quantify the energetic cost for a Beltrami field on $\mathbb T^3$ to simulate a tape-bounded Turing machine, thus providing additional support for the space-bounded Church-Turing thesis. Another implication of our construction is that a Gaussian random Beltrami field on Euclidean space exhibits arbitrarily high computational complexity with probability $1$. Finally, our proof also yields Turing complete flows and diffeomorphisms on $\mathbb{S}^2$ with zero topological entropy, thus disclosing a certain degree of independence within different hierarchies of complexity.


title: Optimal Approximate Distance Oracle for Planar Graphs

abstract: A $(1+\epsilon)$-approximate distance oracle of an edge-weighted graph is a data structure that returns an approximate shortest path distance between any two query vertices up to a $(1+\epsilon)$ factor. Thorup (FOCS 2001, JACM 2004) and Klein (SODA 2002) independently constructed a $(1+\epsilon)$-approximate distance oracle with $O(n\log n)$ space, measured in number of words, and $O(1)$ query time when $G$ is an undirected planar graph with $n$ vertices and $\epsilon$ is a fixed constant. Many follow-up works gave $(1+\epsilon)$-approximate distance oracles with various trade-offs between space and query time. However, improving $O(n\log n)$ space bound without sacrificing query time remains an open problem for almost two decades. In this work, we resolve this problem affirmatively by constructing a $(1+\epsilon)$-approximate distance oracle with optimal $O(n)$ space and $O(1)$ query time for undirected planar graphs and fixed $\epsilon$. We also make substantial progress for planar digraphs with non-negative edge weights. For fixed $\epsilon > 0$, we give a $(1+\epsilon)$-approximate distance oracle with space $o(n\log(Nn))$ and $O(\log\log(Nn)$ query time; here $N$ is the ratio between the largest and smallest positive edge weight. This improves Thorup's (FOCS 2001, JACM 2004) $O(n\log(Nn)\log n)$ space bound by more than a logarithmic factor while matching the query time of his structure. This is the first improvement for planar digraphs in two decades, both in the weighted and unweighted setting.


title: On the effective dimension and multilevel Monte Carlo

abstract: I consider the problem of integrating a function $f$ over the $d$-dimensional unit cube. I describe a multilevel Monte Carlo method that estimates the integral with variance at most $\epsilon^{2}$ in $O(d+\ln(d)d_{t}\epsilon^{-2})$ time, for $\epsilon>0$, where $d_{t}$ is the truncation dimension of $f$. In contrast, the standard Monte Carlo method typically achieves such variance in $O(d\epsilon^{-2})$ time. A lower bound of order $d+d_{t}\epsilon^{-2}$ is described for a class of multilevel Monte Carlo methods.


title: Quantification of systematic errors in the electron density and temperature measured with Thomson scattering at W7-X

abstract: The electron density and temperature profiles measured with Thomson scattering at the stellarator Wendelstein 7-X show features which seem to be unphysical, but so far could not be associated with any source of error considered in the data processing. A detailed Bayesian analysis reveals that errors in the spectral calibration cannot explain the features observed in the profiles. Rather, it seems that small fluctuations in the laser position are sufficient to affect the profile substantially. The impact of these fluctuations depends on the laser position itself, which, in turn, provides a method to find the optimum laser alignment in the future.


title: Machine Learning Product State Distributions from Initial Reactant States for a Reactive Atom-Diatom Collision System

abstract: A machine learned (ML) model for predicting product state distributions from specific initial states (state-to-distribution or STD) for reactive atom-diatom collisions is presented and quantitatively tested for the N($^4$S)+O$_{2}$(X$^3 \Sigma_{\rm g}^{-}$) $\rightarrow$ NO(X$^2\Pi$) +O($^3$P) reaction. The reference data set for training the neural network (NN) consists of final state distributions determined from explicit quasi-classical trajectory (QCT) simulations for $\sim 2000$ initial conditions. Overall, the prediction accuracy as quantified by the root-mean-squared difference $(\sim 0.003)$ and the $R^2$ $(\sim 0.99)$ between the reference QCT and predictions of the STD model is high for the test set and off-grid state specific initial conditions and for initial conditions drawn from reactant state distributions characterized by translational, rotational and vibrational temperatures. Compared with a more coarse grained distribution-to-distribution (DTD) model evaluated on the same initial state distributions, the STD model shows comparable performance with the additional benefit of the state resolution in the reactant preparation. Starting from specific initial states also leads to a more diverse range of final state distributions which requires a more expressive neural network to be used compared with DTD. Direct comparison between explicit QCT simulations, the STD model, and the widely used Larsen-Borgnakke (LB) model shows that the STD model is quantitative whereas the LB model is qualitative at best for rotational distributions $P(j')$ and fails for vibrational distributions $P(v')$. As such the STD model can be well-suited for simulating nonequilibrium high-speed flows, e.g., using the direct simulation Monte Carlo method.


title: System Level Disturbance Reachable Sets and their Application to Tube-based MPC

abstract: Tube-based model predictive control (MPC) methods leverage tubes to bound deviations from a nominal trajectory due to uncertainties in order to ensure constraint satisfaction. This paper presents a novel tube-based MPC formulation based on system level disturbance reachable sets (SL-DRS), which leverage the affine system level parameterization (SLP). We show that imposing a finite impulse response (FIR) constraint on the affine SLP guarantees containment of all future deviations in a finite sequence of SL-DRS. This allows us to formulate a system level tube-MPC (SLTMPC) method using the SL-DRS as tubes, which enables concurrent optimization of the nominal trajectory and the tubes, while using a positively invariant terminal set. Finally, we show that the SL-DRS tubes can also be computed offline.


title: A nutritional strategy to promote gilthead seabream performance under low temperatures

abstract: Gilthead seabream (Sparus aurata) is vulnerable to low water temperature, which may occur in the Southern Europe and Mediterranean region during Winter. Fish are poikilothermic animals, so feed intake, digestion, metabolism and ultimately growth are affected by water temperature. This study aimed to evaluate growth performance, feed utilisation, nutrient apparent digestibility, and N losses to the environment in seabream juveniles reared under low temperature (13 degrees Celsius). Three isolipid and isoenergetic diets were formulated: a commercial-like diet (COM) with 44% crude protein and 27.5% fishmeal; and 2 diets with 42% CP (ECO and ECOSup), reduced FM inclusion, and 15% poultry meal. ECOSup diet was supplemented with a mix of feed additives intended to promote fish growth and feed intake. The ECO diets presented lower production costs than the COM diet and included more sustainable ingredients. Seabream juveniles (154.5 g) were randomly assigned to triplicate tanks and fed the diets for 84 days. Fish fed the ECOSup and COM diets attained a similar final body weight. ECOSup fed fish presented significantly higher HSI than COM fed fish, probably due to higher hepatic glycogen reserves. The VSI of ECOSup fed fish were significantly lower compared to COM fed fish, which is a positive achievement from a consumer point of view. Nutrient digestibility was similar in ECOSup and COM diets. Feeding fish with the ECO diets resulted in lower faecal N losses when compared to COM fed fish. Feeding seabream with an eco-friendly diet with a mix of feed additives promoted growth, improved fish nutritional status and minimised N losses to the environment whilst lowering production costs. Nutritional strategies that ultimately promote feed intake and diet utilisation are valuable tools that may help conditioning fish to sustain growth even under adverse conditions.


title: Standard Model soft QCD at ATLAS and CMS

abstract: Recent results in soft QCD at LHC by the ATLAS, CMS and TOTEM Collaborations are presented. Special focus is reserved to studies in diffractive and forward physics, and to the underlying event, with comparison with previous results and highlighting novel techniques.


title: Toward twist-2 $T$-odd transverse-momentum-dependent gluon distributions: the $f$-type linearity function

abstract: We present preliminary results on the leading-twist transverse-momentum-dependent (TMD) gluon linearity function, which is directly connected to nucleon transverse-spin asymmetries originating from the density of linearly-polarized gluons. The function is calculated in a spectator-model framework for the parent proton and for the $f$-type gauge-link structure. Our work represents a further step toward the definition of a complete set of gluon TMD distributions at twist-2 that can be used to access the gluon dynamics inside nucleons and nuclei at new-generation colliders.