title: MetaFormer: A Unified Meta Framework for Fine-Grained Recognition

abstract: Fine-Grained Visual Classification(FGVC) is the task that requires recognizing the objects belonging to multiple subordinate categories of a super-category. Recent state-of-the-art methods usually design sophisticated learning pipelines to tackle this task. However, visual information alone is often not sufficient to accurately differentiate between fine-grained visual categories. Nowadays, the meta-information (e.g., spatio-temporal prior, attribute, and text description) usually appears along with the images. This inspires us to ask the question: Is it possible to use a unified and simple framework to utilize various meta-information to assist in fine-grained identification? To answer this problem, we explore a unified and strong meta-framework(MetaFormer) for fine-grained visual classification. In practice, MetaFormer provides a simple yet effective approach to address the joint learning of vision and various meta-information. Moreover, MetaFormer also provides a strong baseline for FGVC without bells and whistles. Extensive experiments demonstrate that MetaFormer can effectively use various meta-information to improve the performance of fine-grained recognition. In a fair comparison, MetaFormer can outperform the current SotA approaches with only vision information on the iNaturalist2017 and iNaturalist2018 datasets. Adding meta-information, MetaFormer can exceed the current SotA approaches by 5.9% and 5.3%, respectively. Moreover, MetaFormer can achieve 92.3% and 92.7% on CUB-200-2011 and NABirds, which significantly outperforms the SotA approaches. The source code and pre-trained models are released athttps://github.com/dqshuai/MetaFormer.


title: Efficient Distinction between Quantum Direct and Common Causes and its Experimental Verification

abstract: Identifying the causal structures between two statistically correlated events has been widely investigated in many fields of science. While some of the well-studied classical methods are carefully generalized to quantum version of causal inference for certain cases, an effective and efficient way to detect the more general quantum causal structures is still lacking. Here, we introduce a quantity named `Causal Determinant' to efficiently identify the quantum causal structures between two quantum systems and experimentally verify the validity of the method. According to the causal determinant, the quantum direct cause imposed by an arbitrary unitary operator can be perfectly discriminated with the quantum common cause, in which the two quantum systems share a joint quantum state. In addition, the causal determinant has the capability to discriminate between more general causal structures and predict the range of their parameters. The ability to detect more general quantum causal structures of our method can shed new light on the field of quantum causal inference.


title: Feeding What You Need by Understanding What You Learned

abstract: Machine Reading Comprehension (MRC) reveals the ability to understand a given text passage and answer questions based on it. Existing research works in MRC rely heavily on large-size models and corpus to improve the performance evaluated by metrics such as Exact Match ($EM$) and $F_1$. However, such a paradigm lacks sufficient interpretation to model capability and can not efficiently train a model with a large corpus. In this paper, we argue that a deep understanding of model capabilities and data properties can help us feed a model with appropriate training data based on its learning status. Specifically, we design an MRC capability assessment framework that assesses model capabilities in an explainable and multi-dimensional manner. Based on it, we further uncover and disentangle the connections between various data properties and model performance. Finally, to verify the effectiveness of the proposed MRC capability assessment framework, we incorporate it into a curriculum learning pipeline and devise a Capability Boundary Breakthrough Curriculum (CBBC) strategy, which performs a model capability-based training to maximize the data value and improve training efficiency. Extensive experiments demonstrate that our approach significantly improves performance, achieving up to an 11.22% / 8.71% improvement of $EM$ / $F_1$ on MRC tasks.


title: Selecting Sub-tables for Data Exploration

abstract: We present a framework for creating small, informative sub-tables of large data tables to facilitate the first step of data science: data exploration. Given a large data table table T, the goal is to create a sub-table of small, fixed dimensions, by selecting a subset of T's rows and projecting them over a subset of T's columns. The question is: which rows and columns should be selected to yield an informative sub-table? We formalize the notion of "informativeness" based on two complementary metrics: cell coverage, which measures how well the sub-table captures prominent association rules in T, and diversity. Since computing optimal sub-tables using these metrics is shown to be infeasible, we give an efficient algorithm which indirectly accounts for association rules using table embedding. The resulting framework can be used for visualizing the complete sub-table, as well as for displaying the results of queries over the sub-table, enabling the user to quickly understand the results and determine subsequent queries. Experimental results show that we can efficiently compute high-quality sub-tables as measured by our metrics, as well as by feedback from user-studies.


title: Multi-level images around Kerr-Newman black holes

abstract: We study the lensing effects of Kerr-Newman (KN) black holes and the multi-level images of luminous sources around the black holes. We obtain the analytic forms of the lens equations in KN spacetime and derive approximately their limiting forms for photon emissions near the bounded photon orbits. After deriving the expressions of three key parameters $\gamma$, $\delta$ and $\tau$, we investigate the primary, the secondary and the higher-order images of spherical sources and equatorial disk-like sources. In particular, we find that the charge of KN black hole has a significant effect on the higher-order images.


title: Data Science vs Putin: How much does each of us pay for Putin's war?

abstract: Putin's Ukraine war has caused gas prices to skyrocket. Because of Europe's dependence on Russian gas supplies, we all pay significantly more for heating, involuntarily helping to fund Russia's war against Ukraine. Based on an analysis of real-time gas price data, we present a calculation that estimates every household's financial contribution for heating paid to Russian gas suppliers daily at current prices - six euros per household and day. We show ways everyone can save energy and help reduce the dependency on Russian gas supply.


title: The M/G/1 retrial queue with event-dependent arrivals

abstract: We introduce a novel single-server queue with general retrial times and event-dependent arrivals. This is a versatile model for the study of service systems, in which the server needs a non-negligible time to retrieve waiting customers upon a service completion, while future arrivals depend on the last realized event. Such a model is motivated by the customers' behaviour in service systems where they decide to join based on the last realized event. We investigate the necessary and sufficient stability condition and derive the stationary distribution both at service completion epochs, and at an arbitrary epoch using the supplementary variable technique. We also study the asymptotic behaviour under high rate of retrials. Performance measures are explicitly derived and extensive numerical examples are performed to investigate the impact of event-dependency. Moreover, constrained optimisation problems are formulated and solved with ultimate goal to investigate the admission control problem.


title: Consequences of neutron decay inside neutron stars

abstract: The hypothesis that neutrons might decay into dark matter is explored using neutron stars as a testing ground. It is found that in order to obtain stars with masses at the upper end of those observed, the dark matter must experience a relatively strong self-interaction. Conservation of baryon number and energy then require that the star must undergo some heating, with a decrease in radius, leading to an increase in speed of rotation over a period of days.


title: Performance of Random Template Banks

abstract: When searching for new gravitational-wave or electromagnetic sources, the $n$ signal parameters (masses, sky location, frequencies,...) are unknown. In practice, one hunts for signals at a discrete set of points in parameter space, called a template bank. These may be constructed systematically as a lattice, or alternatively, by placing templates at randomly selected points in parameter space. Here, we calculate the fraction of signals lost by an $n$-dimensional random template bank (compared to a very finely spaced bank). This fraction is compared to the corresponding loss fraction for the best possible lattice-based template banks containing the same number of grid points. For dimensions $n<4$ the lattice-based template banks significantly outperform the random ones. However, remarkably, for dimensions $n>8$, the difference is negligible. In high dimensions, random template banks outperform the best known lattices.


title: Azimuthal angular correlation as a new boosted top jet substructure

abstract: When a top quark is highly boosted, the $W$ boson from its decay has a substantial linear polarization that results in a $\cos 2\phi$ azimuthal angular correlation among the top decay products. We show that this correlation can be measured for hadronically decayed boosted tops, and its magnitude provides a way to measure the longitudinal polarization of top quark, which is an important probe of new physics that couples to top sector.


title: A randomized singular value decomposition for third-order oriented tensors

abstract: The oriented singular value decomposition (O-SVD) proposed by Zeng and Ng provides a hybrid approach to the t-product based third-order tensor singular value decomposition with the transform matrix being a factor matrix of the higher order singular value decomposition. Continuing along this vein, this paper explores realizing the O-SVD more efficiently by the tensor-train rank-1 decomposition and gives a truncated O-SVD. Motivated by the success of probabilistic algorithms, we develop a randomized version of the O-SVD and present its detailed error analysis. The new algorithm has advantages in efficiency while keeping good accuracy compared with the current tensor decompositions. Our claims are supported by numerical experiments on several oriented tensors from real applications.


title: DrawingInStyles: Portrait Image Generation and Editing with Spatially Conditioned StyleGAN

abstract: The research topic of sketch-to-portrait generation has witnessed a boost of progress with deep learning techniques. The recently proposed StyleGAN architectures achieve state-of-the-art generation ability but the original StyleGAN is not friendly for sketch-based creation due to its unconditional generation nature. To address this issue, we propose a direct conditioning strategy to better preserve the spatial information under the StyleGAN framework. Specifically, we introduce Spatially Conditioned StyleGAN (SC-StyleGAN for short), which explicitly injects spatial constraints to the original StyleGAN generation process. We explore two input modalities, sketches and semantic maps, which together allow users to express desired generation results more precisely and easily. Based on SC-StyleGAN, we present DrawingInStyles, a novel drawing interface for non-professional users to easily produce high-quality, photo-realistic face images with precise control, either from scratch or editing existing ones. Qualitative and quantitative evaluations show the superior generation ability of our method to existing and alternative solutions. The usability and expressiveness of our system are confirmed by a user study.


title: Online List Labeling: Breaking the $\log^2n$ Barrier

abstract: The online list labeling problem is an algorithmic primitive with a large literature of upper bounds, lower bounds, and applications. The goal is to store a dynamically-changing set of $n$ items in an array of $m$ slots, while maintaining the invariant that the items appear in sorted order, and while minimizing the relabeling cost, defined to be the number of items that are moved per insertion/deletion. For the linear regime, where $m = (1 + \Theta(1)) n$, an upper bound of $O(\log^2 n)$ on the relabeling cost has been known since 1981. A lower bound of $\Omega(\log^2 n)$ is known for deterministic algorithms and for so-called smooth algorithms, but the best general lower bound remains $\Omega(\log n)$. The central open question in the field is whether $O(\log^2 n)$ is optimal for all algorithms. In this paper, we give a randomized data structure that achieves an expected relabeling cost of $O(\log^{3/2} n)$ per operation. More generally, if $m = (1 + \varepsilon) n$ for $\varepsilon = O(1)$, the expected relabeling cost becomes $O(\varepsilon^{-1} \log^{3/2} n)$. Our solution is history independent, meaning that the state of the data structure is independent of the order in which items are inserted/deleted. For history-independent data structures, we also prove a matching lower bound: for all $\epsilon$ between $1 / n^{1/3}$ and some sufficiently small positive constant, the optimal expected cost for history-independent list-labeling solutions is $\Theta(\varepsilon^{-1}\log^{3/2} n)$.


title: Bridging the Gap Between Learning in Discrete and Continuous Environments for Vision-and-Language Navigation

abstract: Most existing works in vision-and-language navigation (VLN) focus on either discrete or continuous environments, training agents that cannot generalize across the two. The fundamental difference between the two setups is that discrete navigation assumes prior knowledge of the connectivity graph of the environment, so that the agent can effectively transfer the problem of navigation with low-level controls to jumping from node to node with high-level actions by grounding to an image of a navigable direction. To bridge the discrete-to-continuous gap, we propose a predictor to generate a set of candidate waypoints during navigation, so that agents designed with high-level actions can be transferred to and trained in continuous environments. We refine the connectivity graph of Matterport3D to fit the continuous Habitat-Matterport3D, and train the waypoints predictor with the refined graphs to produce accessible waypoints at each time step. Moreover, we demonstrate that the predicted waypoints can be augmented during training to diversify the views and paths, and therefore enhance agent's generalization ability. Through extensive experiments we show that agents navigating in continuous environments with predicted waypoints perform significantly better than agents using low-level actions, which reduces the absolute discrete-to-continuous gap by 11.76% Success Weighted by Path Length (SPL) for the Cross-Modal Matching Agent and 18.24% SPL for the Recurrent VLN-BERT. Our agents, trained with a simple imitation learning objective, outperform previous methods by a large margin, achieving new state-of-the-art results on the testing environments of the R2R-CE and the RxR-CE datasets.


title: Flare echos from relaxation shocks in perturbed relativistic jets

abstract: One of the main scenarios to account for the multi-wavelength flux variability observed in relativistic jets of AGNs is based on diffusive shock acceleration of a population of relativistic electrons on internal shocks of various origins. Any complete AGN emission scenario has to be able to explain the wide range of observed variability time scales between the radio and gamma-ray band. In addition constraints are also provided by very-long-baseline interferometry (VLBI), which shows a large variety of moving and standing emission zones with distinct behaviors. We aim to characterize the evolution of stationary and moving emission zones in the jet and to study their multi-wavelength signatures through emission maps and light curves. We focus our study on flare events that occur during strong interactions between moving ejecta and stationary recollimation shocks. We simulate relativistic jets with the magneto-hydrodynamic code MPI-AMRVAC and inject non-thermal particle distributions of electrons in shock regions. We follow the propagation of a moving shock and its interactions with a structure of standing re-collimation shocks in the jet. Synchrotron emission and radiative transfer are calculated in the post-processing code RIPTIDE for given observation angles and frequencies. We demonstrate the appearance of trailing components behind the leading moving shock. The latter destabilizes the jet, causing the emergence of oscillating standing shocks and relaxation shocks. Emissions from these regions can dominate the overall flux or lead to flare echos in the light curve. Another observational marker of relaxation shocks appears in time-distance plots of bright VLBI components of the jet. Our scenario provides a plausible explanation for radio VLBI observations of the radio-galaxy 3C 111 where trailing components have been observed during a radio outburst event in 1997.


title: Improved bound for improper colorings of graphs with no odd clique minor

abstract: Strengthening Hadwiger's conjecture, Gerards and Seymour conjectured in 1995 that every graph with no odd $K_t$-minor is properly $(t-1)$-colorable, this is known as the Odd Hadwiger's conjecture. We prove a relaxation of the above conjecture, namely we show that every graph with no odd $K_t$-minor admits a vertex $(2t-2)$-coloring such that all monochromatic components have size at most $\lceil \frac{1}{2}(t-2) \rceil$. The bound on the number of colors is optimal up to a factor of $2$, improves previous bounds for the same problem by Kawarabayashi (2008), Kang and Oum (2019), Liu and Wood (2021), and strengthens a result by van den Heuvel and Wood (2018), who showed that the above conclusion holds under the more restrictive assumption that the graph is $K_t$-minor free. In addition, the bound on the component-size in our result is much smaller than those of previous results, in which the dependency on $t$ was non-explicit. Our short proof combines the method by van den Heuvel and Wood for $K_t$-minor free graphs with some additional ideas, which make the extension to odd $K_t$-minor free graphs possible.


title: Towards Robust Part-aware Instance Segmentation for Industrial Bin Picking

abstract: Industrial bin picking is a challenging task that requires accurate and robust segmentation of individual object instances. Particularly, industrial objects can have irregular shapes, that is, thin and concave, whereas in bin-picking scenarios, objects are often closely packed with strong occlusion. To address these challenges, we formulate a novel part-aware instance segmentation pipeline. The key idea is to decompose industrial objects into correlated approximate convex parts and enhance the object-level segmentation with part-level segmentation. We design a part-aware network to predict part masks and part-to-part offsets, followed by a part aggregation module to assemble the recognized parts into instances. To guide the network learning, we also propose an automatic label decoupling scheme to generate ground-truth part-level labels from instance-level labels. Finally, we contribute the first instance segmentation dataset, which contains a variety of industrial objects that are thin and have non-trivial shapes. Extensive experimental results on various industrial objects demonstrate that our method can achieve the best segmentation results compared with the state-of-the-art approaches.


title: Mediated interactions between ions in quantum degenerate gases

abstract: We explore the interaction between two trapped ions mediated by a surrounding quantum degenerate Bose or Fermi gas. Using perturbation theory valid for weak atom-ion interaction, we show analytically that the interaction mediated by a Bose gas has a power-law behaviour for large distances whereas it has a Yukawa form for intermediate distances. For a Fermi gas, the mediated interaction is given by a power-law for large density and by a Ruderman-Kittel-Kasuya-Yosida form for low density. For strong atom-ion interactions, we use a diagrammatic theory to demonstrate that the mediated interaction can be a significant addition to the bare Coulomb interaction between the ions, when an atom-ion bound state is close to threshold. Finally, we show that the induced interaction leads to substantial and observable shifts in the ion phonon frequencies.


title: Dimits transition in three-dimensional ion-temperature-gradient turbulence

abstract: We extend our previous work on the 2D Dimits transition in ion-scale turbulence (Ivanov et al. 2020) to include variations along the magnetic field. We consider a three-field fluid model for the perturbations of electrostatic potential, ion temperature, and ion parallel flow in a constant-magnetic-curvature geometry without magnetic shear. It is derived in the cold-ion, long-wavelength asymptotic limit of the gyrokinetic theory. Just as in the 2D model, a low-transport (Dimits) regime exists and is found to be dominated by a quasi-static staircase-like arrangement of strong zonal flows and zonal temperature. This zonal staircase is formed and maintained by a negative turbulent viscosity for the zonal flows. Unlike the 2D model, the 3D one does not suffer from an unphysical blow up beyond the Dimits threshold where the staircase becomes nonlinearly unstable. Instead, a well-defined finite-amplitude saturated state is established. This qualitative difference between 2D and 3D is due to the appearance of small-scale `parasitic' modes that exist only if we allow perturbations to vary along the magnetic field lines. These modes extract energy from the large-scale perturbations and provide an effective enhancement of large-scale thermal diffusion, thus aiding the energy transfer from large injection scales to small dissipative ones. We show that in our model, the parasitic modes always favour a zonal-flow-dominated state. In fact, a Dimits state with a zonal staircase is achieved regardless of the strength of the linear drive provided the system is sufficiently extended along the magnetic field and sufficient parallel resolution is provided.


title: Don't Be So Dense: Sparse-to-Sparse GAN Training Without Sacrificing Performance

abstract: Generative adversarial networks (GANs) have received an upsurging interest since being proposed due to the high quality of the generated data. While achieving increasingly impressive results, the resource demands associated with the large model size hinders the usage of GANs in resource-limited scenarios. For inference, the existing model compression techniques can reduce the model complexity with comparable performance. However, the training efficiency of GANs has less been explored due to the fragile training process of GANs. In this paper, we, for the first time, explore the possibility of directly training sparse GAN from scratch without involving any dense or pre-training steps. Even more unconventionally, our proposed method enables directly training sparse unbalanced GANs with an extremely sparse generator from scratch. Instead of training full GANs, we start with sparse GANs and dynamically explore the parameter space spanned over the generator throughout training. Such a sparse-to-sparse training procedure enhances the capacity of the highly sparse generator progressively while sticking to a fixed small parameter budget with appealing training and inference efficiency gains. Extensive experiments with modern GAN architectures validate the effectiveness of our method. Our sparsified GANs, trained from scratch in one single run, are able to outperform the ones learned by expensive iterative pruning and re-training. Perhaps most importantly, we find instead of inheriting parameters from expensive pre-trained GANs, directly training sparse GANs from scratch can be a much more efficient solution. For example, only training with a 80% sparse generator and a 70% sparse discriminator, our method can achieve even better performance than the dense BigGAN.


title: Remiod: Reference-based Controlled Multiple Imputation of Longitudinal Binary and Ordinal Outcomes with non-ignorable missingness

abstract: Missing data on response variables are common in clinical studies. Corresponding to the uncertainty of missing mechanism, theoretical frameworks on controlled imputation have been developed. In practice, it is recommended to conduct a statistically valid analysis under the primary assumptions on missing data, followed by sensitivity analysis under alternative assumptions to assess the robustness of results. Due to the availability of software, controlled multiple imputation (MI) procedures, including delta-based and reference-based approaches, have become popular for analyzing continuous variables under missing-not-at-random assumptions. Similar tools, however, still limit application of these methods to categorical data. In this paper, we introduce the R package \textbf{remiod}, which utilizes the Bayesian framework to perform imputation in regression models on binary and ordinal outcomes. Following outlining theoretical backgrounds, usage and features of \textbf{remiod} are described and illustrated using examples.


title: Rib Suppression in Digital Chest Tomosynthesis

abstract: Digital chest tomosynthesis (DCT) is a technique to produce sectional 3D images of a human chest for pulmonary disease screening, with 2D X-ray projections taken within an extremely limited range of angles. However, under the limited angle scenario, DCT contains strong artifacts caused by the presence of ribs, jamming the imaging quality of the lung area. Recently, great progress has been achieved for rib suppression in a single X-ray image, to reveal a clearer lung texture. We firstly extend the rib suppression problem to the 3D case at the software level. We propose a $\textbf{T}$omosynthesis $\textbf{RI}$b Su$\textbf{P}$pression and $\textbf{L}$ung $\textbf{E}$nhancement $\textbf{Net}$work (TRIPLE-Net) to model the 3D rib component and provide a rib-free DCT. TRIPLE-Net takes the advantages from both 2D and 3D domains, which model the ribs in DCT with the exact FBP procedure and 3D depth information, respectively. The experiments on simulated datasets and clinical data have shown the effectiveness of TRIPLE-Net to preserve lung details as well as improve the imaging quality of pulmonary diseases. Finally, an expert user study confirms our findings.


title: The Physics and Metaphysics of Pure Shape Dynamics

abstract: The goal of this essay is twofold. First, it provides a quick look at the foundations of modern relational mechanics by tracing its development from Julian Barbour and Bruno Bertotti's original ideas until present-day's pure shape dynamics. Secondly, it discusses the most appropriate metaphysics for pure shape dynamics, showing that relationalism is more of a nuanced thesis rather than an elusive one. The chapter ends with a brief assessment of the prospects of pure shape dynamics in light of quantum physics.


title: Algebraic theory of phase retrieval

abstract: The purpose of this article is to discuss recent advances in the growing field of phase retrieval, and to publicize open problems that we believe will be of interest to mathematicians in general, and algebraists in particular.

categories: eess.SP cs.IT math-ph math.AG math.IT math.MP


title: On sheaf cohomology for supergroups arising from simple classical Lie superalgebras

abstract: In this paper the authors study the behavior of the sheaf cohomology functors $R^{\bullet}\text{ind}_{B}^{G}(-)$ where $G$ is an algebraic group scheme corresponding to a simple classical Lie superalgebra and $B$ is a BBW parabolic subgroup as defined by D. Grantcharov, N. Grantcharov, Nakano and Wu. We provide a systematic treatment that allows us to study the behavior of these cohomology groups $R^{\bullet}\text{ind}_{B}^{G}L_{\mathfrak f}(\lambda)$ where $L_{\mathfrak f}(\lambda)$ is an irreducible representation for the detecting subalgebra ${\mathfrak f}$. In particular, we prove an analog of Kempf's vanishing theorem and the Bott-Borel-Weil theorem for large weights.


title: Boosting human decision-making with AI-generated decision aids

abstract: Human decision-making is plagued by many systematic errors. Many of these errors can be avoided by providing decision aids that guide decision-makers to attend to the important information and integrate it according to a rational decision strategy. Designing such decision aids used to be a tedious manual process. Advances in cognitive science might make it possible to automate this process in the future. We recently introduced machine learning methods for discovering optimal strategies for human decision-making automatically and an automatic method for explaining those strategies to people. Decision aids constructed by this method were able to improve human decision-making. However, following the descriptions generated by this method is very tedious. We hypothesized that this problem can be overcome by conveying the automatically discovered decision strategy as a series of natural language instructions for how to reach a decision. Experiment 1 showed that people do indeed understand such procedural instructions more easily than the decision aids generated by our previous method. Encouraged by this finding, we developed an algorithm for translating the output of our previous method into procedural instructions. We applied the improved method to automatically generate decision aids for a naturalistic planning task (i.e., planning a road trip) and a naturalistic decision task (i.e., choosing a mortgage). Experiment 2 showed that these automatically generated decision-aids significantly improved people's performance in planning a road trip and choosing a mortgage. These findings suggest that AI-powered boosting might have potential for improving human decision-making in the real world.


title: Cosine Model Watermarking Against Ensemble Distillation

abstract: Many model watermarking methods have been developed to prevent valuable deployed commercial models from being stealthily stolen by model distillations. However, watermarks produced by most existing model watermarking methods can be easily evaded by ensemble distillation, because averaging the outputs of multiple ensembled models can significantly reduce or even erase the watermarks. In this paper, we focus on tackling the challenging task of defending against ensemble distillation. We propose a novel watermarking technique named CosWM to achieve outstanding model watermarking performance against ensemble distillation. CosWM is not only elegant in design, but also comes with desirable theoretical guarantees. Our extensive experiments on public data sets demonstrate the excellent performance of CosWM and its advantages over the state-of-the-art baselines.


title: A Modular Approach to the Embodiment of Hand Motions from Human Demonstrations

abstract: Manipulating objects with robotic hands is a complicated task. Not only the fingers of the hand, but also the pose of the robot's end effector need to be coordinated. Using human demonstrations of movements is an intuitive and data-efficient way of guiding the robot's behavior. We propose a modular framework with an automatic embodiment mapping to transfer recorded human hand motions to robotic systems. In this work, we use motion capture to record human motion. We evaluate our approach on eight challenging tasks, in which a robotic hand needs to grasp and manipulate either deformable or small and fragile objects. We test a subset of trajectories in simulation and on a real robot and the overall success rates are aligned.


title: Atomic Size Pearls being Dark Matter giving Electron Signal

abstract: We seek to explain both the seeming observation of dark matter by the seasonal variation of the DAMA-LIBRA data and the observation of "electron recoil" events at Xenon1T in which the liquid Xe scintillator was excited by electrons - in excess to the expected background - by the same dark matter model. In our model the dark matter consists of bubbles of a new type of vacuum containing ordinary atomic matter, say diamond, under high pressure ensured by the surface tension of the separation domain wall. We also seek to explain the self interactions of dark matter suggested by astronomical studies of dwarf galaxies and the central structure of galaxy clusters. The interaction with matter in the shielding is responsible for slowing the dark matter down to a low terminal velocity, so underground detectors have insufficient energy for detection. Further we explain the "mysterious" X-ray line of 3.5 keV, even the 3.5 keV X-ray radiation from the Tycho supernova remnant. The DAMA-LIBRA and Xenon1T experiments see decaying dark matter pearls. Both these underground experiments see events with about 3.5 keV energy, just the energy of the X-ray line. We fit numerically the cross section over mass ratio for the self interaction of the dark matter as observed in the study of dwarf galaxies and find the size of the pearls to be rather close to the smallest possible in our model. However DAMA cause some problem to get heavier pearls. Also the total energy of the dark matter pearls stopped in the shield is reasonably matching. Accepting the different phases of the vacuum the model could be realized inside the Standard Model, our whole scheme could be realized inside the Standard Model.


title: The Completeness of 2D Rubik's Shapes

abstract: The Rubik's cube was invented in 1974 by Erno Rubik, who had no idea of the incredible popularity and mathematical fascinations his toy would bring. Through the years of study on the mathematical properties of the cube, the Rubik's Cube group was introduced to represent all possible moves one could perform on the cube. In this paper, we define a planar analogue to the Rubik's cube, which we dub the Rubik's Square, and prove that the Rubik's square is complete in the sense that given any two configurations there is a sequence of moves which changes one to the other. The Rubik's cube does not have this property. We then abstract the concept of the Rubik's Square to a Rubik's Shape and analyse the completeness in this more general setting.


title: Varying fine-structure constant cosmography

abstract: Cosmography is a phenomenological and relatively model-independent approach to cosmology, where physical quantities are expanded as a Taylor series in the cosmological redshift, or in related variables. Here we apply this methodology to constrain possible cosmological variations of the fine-structure constant, $\alpha$. Two peculiarities of this case are the existence of high-redshift data, and the fact that one term in the series is directly and tightly constraint by local laboratory tests with atomic clocks. We use this atomic clock data, together with direct model-independent high-resolution astrophysical spectroscopy measurements of $\alpha$ up to redshift $z\sim7$ and additional model-dependent constraints on $\alpha$ from the cosmic microwave background and big bang nucleosynthesis, to place stringent (parts per million level) constraints on the first two terms in the $\alpha$ cosmographic series.


title: Laplace and Dirac Operators on Graphs

abstract: Discrete versions of the Laplace and Dirac operators haven been studied in the context of combinatorial models of statistical mechanics and quantum field theory. In this paper we introduce several variations of the Laplace and Dirac operators on graphs, and we investigate graph-theoretic versions of the Schr\"odinger and Dirac equation. We provide a combinatorial interpretation for solutions of the equations and we prove gluing identities for the Dirac operator on lattice graphs, as well as for graph Clifford algebras.


title: Entanglement and correlations in fast collective neutrino flavor oscillations

abstract: Collective neutrino oscillations play a crucial role in transporting lepton flavor in astrophysical settings like supernovae and neutron star binary merger remnants, which are characterized by large neutrino densities. In these settings, simulations in the mean-field approximation show that neutrino-neutrino interactions can overtake vacuum oscillations and give rise to fast collective flavor evolution on time-scales $t\propto\mu^{-1}$, with $\mu$ proportional to the local neutrino density. In this work, we study the full out-of-equilibrium flavor dynamics in simple multi-angle geometries displaying fast oscillations in the mean field linear stability analysis. Focusing on simple initial conditions, we analyze the production of pair correlations and entanglement in the complete many-body-dynamics as a function of the number $N$ of neutrinos in the system, for up to thousands of neutrinos. Similarly to simpler geometries with only two neutrino beams, we identify three regimes: stable configurations with vanishing flavor oscillations, marginally unstable configurations with evolution occurring on long time scales $\tau\approx\mu^{-1}\sqrt{N}$, and unstable configurations showing flavor evolution on short time scales $\tau\approx\mu^{-1}\log(N)$. We present evidence that these fast collective modes are generated by the same dynamical phase transition which leads to the slow bipolar oscillations, establishing a connection between these two phenomena and explaining the difference in their time scales. We conclude by discussing a semi-classical approximation which reproduces the entanglement entropy at short to medium time scales and can be potentially useful in situations with more complicated geometries where classical simulation methods starts to become inefficient.


title: Deterministic Identification for Molecular Communications over the Poisson Channel

abstract: Various applications of molecular communications (MC) are event-triggered, and, as a consequence, the prevalent Shannon capacity may not be the right measure for performance assessment. Thus, in this paper, we motivate and establish the identification capacity as an alternative metric. In particular, we study deterministic identification (DI) for the discrete-time Poisson channel (DTPC), subject to an average and a peak power constraint, which serves as a model for MC systems employing molecule counting receivers. It is established that the codebook size for this channel scales as $2^{(n\log n)R}$, where $n$ and $R$ are the codeword length and coding rate, respectively. Lower and upper bounds on the DI capacity of the DTPC are developed. The obtained large capacity of the DI channel sheds light on the performance of natural DI systems such as natural olfaction, which are known for their extremely large chemical discriminatory power in biology. Furthermore, numerical simulations for the empirical miss-identification and false identification error rates are provided for finite length codes. This allows us to quantify the scale of error reduction in terms of the codeword length.


title: The effect of convection on the position of the free liquid surface under zero and terrestrial gravity

abstract: Theresults of numerical modelinghas based on the solution of Navier-Stokes equations for convection oftwo-layer system air-water are presented.The influence of thermocapillary and gravitational convection on the deformation and time averaged position of the interface of air-water system under sudden lateral heating of a square-section cavity are shown.


title: Stratification and segregation under laminar convection

abstract: The paper is devoted to the study of the formation of stratification in an incompressible fluid due to convective laminar flows in horizontal layers heated from the side. Medium and intensive modes of stationary laminar thermal, concentrational and thermo-concentrational (in particular thermohaline) convection are considered, in which nonlinear flow features are manifested, which can radically change the flow structure and the characteristics of heat and mass transfer. The solutions of the problems of laminar problems of convection show the features of the formation of layered structures, vertical temperature and concentration stratification depending on the determining dimensionless parameters. The metastable instability of the direction of the stratification vector (the location of the free surface) in weightlessness in the presence of capillary convection was shown.


title: Heat and mass transfer during crystal growing by the Czochralski method with a submerged vibrator

abstract: It was shown that the vibrations can intensify of heat-mass transfer and mixing of the melt, and also to reduce dynamic and temperature boundary layers and to increase a temperature gradient on melt/crystal interface, that can to increase rate of crystal growth by the Czochralski method with a submerged vibrator. It is can be very impotent for manufacturing and quality of crystals.


title: Constructing Artificial Traffic Fluids by Designing Cruise Controllers

abstract: In this paper, we apply a Control Lyapunov Function methodology to design two families of cruise controllers for the two-dimensional movement of autonomous vehicles on lane-free roads using the bicycle kinematic model. The control Lyapunov functions are based on measures of the energy of the system with the kinetic energy expressed in ways similar to Newtonian or relativistic mechanics. The derived feedback laws (cruise controllers) are decentralized, as each vehicle determines its control input based on its own speed and on the relative speeds and distances from adjacent vehicles and from the boundary of the road. Moreover, the corresponding macroscopic models are derived, obtaining fluid-like models that consist of a conservation equation and a momentum equation with pressure and viscous terms. Finally, we show that, by selecting appropriately the parameters of the feedback laws, we can determine the physical properties of the "traffic fluid", i.e. we get free hand to create an artificial fluid that approximates the emerging traffic flow.


title: A monotonicity version of a concavity theorem of Lieb

abstract: We give a simple proof of a strengthened version of a theorem of Lieb that played a key role in the proof of strong subadditivity of the quantum entropy.


title: A Robust Spectral Algorithm for Overcomplete Tensor Decomposition

abstract: We give a spectral algorithm for decomposing overcomplete order-4 tensors, so long as their components satisfy an algebraic non-degeneracy condition that holds for nearly all (all but an algebraic set of measure $0$) tensors over $(\mathbb{R}^d)^{\otimes 4}$ with rank $n \le d^2$. Our algorithm is robust to adversarial perturbations of bounded spectral norm. Our algorithm is inspired by one which uses the sum-of-squares semidefinite programming hierarchy (Ma, Shi, and Steurer STOC'16, arXiv:1610.01980), and we achieve comparable robustness and overcompleteness guarantees under similar algebraic assumptions. However, our algorithm avoids semidefinite programming and may be implemented as a series of basic linear-algebraic operations. We consequently obtain a much faster running time than semidefinite programming methods: our algorithm runs in time $\tilde O(n^2d^3) \le \tilde O(d^7)$, which is subquadratic in the input size $d^4$ (where we have suppressed factors related to the condition number of the input tensor).


title: Deep Q-Learning Based Resource Allocation in Interference Systems With Outage Constraint

abstract: This correspondence considers the resource allocation problem in wireless interference channel (IC) under link outage constraints. Since the optimization problem is non-convex in nature, existing approaches to find the optimal power allocation are computationally intensive and thus practically infeasible. Recently, deep reinforcement learning has shown promising outcome in solving non-convex optimization problems with reduced complexity. In this correspondence, we utilize a deep Q-learning (DQL) approach which interacts with the wireless environment and learns the optimal power allocation of a wireless IC while maximizing overall sum-rate of the system and maintaining reliability requirement of each link. We have used two separate deep Q-networks to remove the inherent instability in learning process. Simulation results demonstrate that the proposed DQL approach outperforms existing geometric programming based solution.


title: Adversarial Dual-Student with Differentiable Spatial Warping for Semi-Supervised Semantic Segmentation

abstract: A common challenge posed to robust semantic segmentation is the expensive data annotation cost. Existing semi-supervised solutions show great potential for solving this problem. Their key idea is constructing consistency regularization with unsupervised data augmentation from unlabeled data for model training. The perturbations for unlabeled data enable the consistency training loss, which benefits semi-supervised semantic segmentation. However, these perturbations destroy image context and introduce unnatural boundaries, which is harmful for semantic segmentation. Besides, the widely adopted semi-supervised learning framework, i.e. mean-teacher, suffers performance limitation since the student model finally converges to the teacher model. In this paper, first of all, we propose a context friendly differentiable geometric warping to conduct unsupervised data augmentation; secondly, a novel adversarial dual-student framework is proposed to improve the Mean-Teacher from the following two aspects: (1) dual student models are learned independently except for a stabilization constraint to encourage exploiting model diversities; (2) adversarial training scheme is applied to both students and the discriminators are resorted to distinguish reliable pseudo-label of unlabeled data for self-training. Effectiveness is validated via extensive experiments on PASCAL VOC2012 and Cityscapes. Our solution significantly improves the performance and state-of-the-art results are achieved on both datasets. Remarkably, compared with fully supervision, our solution achieves comparable mIoU of 73.4% using only 12.5% annotated data on PASCAL VOC2012. Our codes and models are available at https://github.com/cao-cong/ADS-SemiSeg.


title: Hundreds of new satellites of figure-eight orbit computed with high precision

abstract: Satellites (topological powers) of the famous figure-eight orbit are special periodic solutions of the planar three-body problem. In this paper we use a modified Newton's method based on the Continuous analog of Newton's method and high precision arithmetic for a purposeful numerical search of new satellites of the figure-eight orbit. Over 700 new satellites are found, including 76 new linearly stable ones. 7 of the newly found linearly stable satellites are choreographies. The linear stability is checked by a high precision computing of the eigenvalues of the monodromy matrices. The initial conditions of all found solutions are given with 150 correct decimal digits.

categories: math.NA astro-ph.EP cs.NA nlin.CD physics.comp-ph


title: Machine Learning Applications in Lung Cancer Diagnosis, Treatment and Prognosis

abstract: The recent development of imaging and sequencing technologies enables systematic advances in the clinical study of lung cancer. Meanwhile, the human mind is limited in effectively handling and fully utilizing the accumulation of such enormous amounts of data. Machine learning-based approaches play a critical role in integrating and analyzing these large and complex datasets, which have extensively characterized lung cancer through the use of different perspectives from these accrued data. In this article, we provide an overview of machine learning-based approaches that strengthen the varying aspects of lung cancer diagnosis and therapy, including early detection, auxiliary diagnosis, prognosis prediction and immunotherapy practice. Moreover, we highlight the challenges and opportunities for future applications of machine learning in lung cancer.


title: Revisiting Degeneracy, Strict Feasibility, Stability, in Linear Programming

abstract: Currently, the simplex method and the interior point method are indisputably the most popular algorithms for solving linear programs, LPs. Unlike general conic programs, LPs with a finite optimal value do not require strict feasibility in order to establish strong duality. Hence strict feasibility is seldom a concern, even though strict feasibility is equivalent to stability and a compact dual optimal set. This lack of concern is also true for other types of degeneracy of basic feasible solutions in LP. In this paper we discuss that the specific degeneracy that arises from lack of strict feasibility necessarily causes difficulties in both simplex and interior point methods. In particular, we show that the lack of strict feasibility implies that every basic feasible solution, BFS, is degenerate; thus conversely, the existence of a nondegenerate BFS implies that strict feasibility (regularity) holds. We prove the results using facial reduction and simple linear algebra. In particular, the facially reduced system reveals the implicit non-surjectivity of the linear map of the equality constraint system. As a consequence, we emphasize that facial reduction involves two steps where, the first guarantees strict feasibility, and the second recovers full row rank of the constraint matrix. This illustrates the implicit singularity of problems where strict feasibility fails, and also helps in obtaining new efficient techniques for preproccessing. We include an efficient preprocessing method that can be performed as an extension of phase-I of the two-phase simplex method. We show that this can be used to avoid the loss of precision for many well known problem sets in the literature, e.g., the NETLIB problem set.


title: Distributed Optimal Power Flow for VSC-MTDC Meshed AC/DC Grids Using ALADIN

abstract: The increasing application of voltage source converter (VSC) high voltage direct current (VSC-HVDC) technology in power grids has raised the importance of incorporating DC grids and converters into the existing transmission network. This poses significant challenges in dealing with the resulting optimal power flow (OPF) problem. In this paper, a recently proposed nonconvex distributed optimization algorithm -- Augmented Lagrangian based Alternating Direction Inexact Newton method (ALADIN), is tailored to solve the nonconvex AC/DC OPF problem for emerging voltage source converter (VSC) based multiterminal high voltage direct current (VSC-MTDC) meshed AC/DC hybrid systems. The proposed scheme decomposes this AC/DC hybrid OPF problem and handles it in a fully distributed way. Compared to the existing state-of-art Alternating Direction Method of Multipliers(ADMM), which is in general, not applicable for nonconvex problems, ALADIN has a theoretical convergence guarantee. Applying these two approaches to (VSC-MTDC) coupled with an IEEE benchmark AC power system illustrates that the tailored ALADIN outperforms ADMM in convergence speed and numerical robustness.


title: ClueGraphSum: Let Key Clues Guide the Cross-Lingual Abstractive Summarization

abstract: Cross-Lingual Summarization (CLS) is the task to generate a summary in one language for an article in a different language. Previous studies on CLS mainly take pipeline methods or train the end-to-end model using the translated parallel data. However, the quality of generated cross-lingual summaries needs more further efforts to improve, and the model performance has never been evaluated on the hand-written CLS dataset. Therefore, we first propose a clue-guided cross-lingual abstractive summarization method to improve the quality of cross-lingual summaries, and then construct a novel hand-written CLS dataset for evaluation. Specifically, we extract keywords, named entities, etc. of the input article as key clues for summarization and then design a clue-guided algorithm to transform an article into a graph with less noisy sentences. One Graph encoder is built to learn sentence semantics and article structures and one Clue encoder is built to encode and translate key clues, ensuring the information of important parts are reserved in the generated summary. These two encoders are connected by one decoder to directly learn cross-lingual semantics. Experimental results show that our method has stronger robustness for longer inputs and substantially improves the performance over the strong baseline, achieving an improvement of 8.55 ROUGE-1 (English-to-Chinese summarization) and 2.13 MoverScore (Chinese-to-English summarization) scores over the existing SOTA.


title: pylspack: Parallel algorithms and data structures for sketching, column subset selection, regression and leverage scores

abstract: We present parallel algorithms and data structures for three fundamental operations in Numerical Linear Algebra: (i) Gaussian and CountSketch random projections and their combination, (ii) computation of the Gram matrix and (iii) computation of the squared row norms of the product of two matrices, with a special focus on "tall-and-skinny" matrices, which arise in many applications. We provide a detailed analysis of the ubiquitous CountSketch transform and its combination with Gaussian random projections, accounting for memory requirements, computational complexity and workload balancing. We also demonstrate how these results can be applied to column subset selection, least squares regression and leverage scores computation. These tools have been implemented in pylspack, a publicly available Python package (https://github.com/IBM/pylspack) whose core is written in C++ and parallelized with OpenMP, and which is compatible with standard matrix data structures of SciPy and NumPy. Extensive numerical experiments indicate that the proposed algorithms scale well and significantly outperform existing libraries for tall-and-skinny matrices.


title: Accelerated carrier invoice factoring using predictive freight transport events

abstract: Invoice factoring is an invoice financing process where business organizations sell their invoices to banks or financial institutions at a discount to gain faster access to the invoice amount. In global trade, ocean and land carriers exercise invoice factoring to gain quick access to the money they get paid for the shipment of consignments by shippers. Shippers typically clear the invoice payment within 60-90 days of goods getting delivered. In order to get early access to capital, carriers initiate invoice factoring after the completion of goods delivery to the shipper. In this work, we provide an approach to enable accelerated carrier invoice factoring even before the delivery of goods. Carrier invoice value for a given shipment depends on the actual values of shipment tracking events. We predict the carrier invoice value at different freight transport milestone events as the goods transportation progresses from supplier to shipper using smart contracts on a blockchain network that is operated by the global trade logistics participants. The prediction at a given stage is based on past shipment tracking events and predicted future shipment tracking events. Accurate prediction of invoice value for ongoing shipment enables the carrier organization to initiate invoice factoring on the trade finance network before the completion of goods delivery to the shipper. Further, based on the past accuracy of prediction models, the financial institutions may choose to release the invoice amount in installments at different freight transportation milestone events.


title: Rifts in Rafts

abstract: A particle raft floating on an expanding liquid substrate provides a macroscopic analog for studying material failure. The time scales in this system allow both particle-relaxation dynamics and rift formation to be resolved. In our experiments, a raft, an aggregate of particles, is stretched uniaxially by the expansion of the air-liquid interface on which it floats. Its failure morphology changes continuously with pulling velocity. This can be understood as a competition between two velocity scales: the speed of re-aggregation, in which particles relax towards a low-energy configuration determined by viscous and capillary forces, and the difference of velocity between neighboring particles caused by the expanding fluid. This competition selects the cluster length, i.e., the distance between adjacent rifts. A model based on this competition is consistent with the experimental failure patterns.


title: Object-centric Process Predictive Analytics

abstract: Object-centric processes (a.k.a. Artifact-centric processes) are implementations of a paradigm where an instance of one process is not executed in isolation but interacts with other instances of the same or other processes. Interactions take place through bridging events where instances exchange data. Object-centric processes are recently gaining popularity in academia and industry, because their nature is observed in many application scenarios. This poses significant challenges in predictive analytics due to the complex intricacy of the process instances that relate to each other via many-to-many associations. Existing research is unable to directly exploit the benefits of these interactions, thus limiting the prediction quality. This paper proposes an approach to incorporate the information about the object interactions into the predictive models. The approach is assessed on real-life object-centric process event data, using different KPIs. The results are compared with a naive approach that overlooks the object interactions, thus illustrating the benefits of their use on the prediction quality.


title: On discrepancy, intrinsic Diophantine approximation, and spectral gaps

abstract: In the present paper we establish bounds for the size of the spectral gap for actions of algebraic groups on certain homogeneous spaces. Our approach is based on estimating operator norms of suitable averaging operators, and we develop techniques for establishing both upper and lower bounds for such norms. We shall show that this analytic problem is closely related to the arithmetic problem of establishing bounds on the discrepancy of distribution for rational points on algebraic group varieties. As an application, we show how to establish an effective bound for property $\tau$ of congruence subgroups of arithmetic lattices in algebraic groups which are forms of $SL(2)$, using estimates in intrinsic Diophantine approximation which follow from Heath-Brown's analysis of rational points on 3-dimensional quadratic surfaces.


title: Localized scalar structures around static black holes

abstract: In this work we address a way to capture scalar field solutions on static spacetimes by using BPS formalism and relaxing the general covariance condition. We focus on configurations where the background geometry describes topological black holes and present both analytical and numerical solutions, in addition to discussing the use of conserved charges associated to such field configurations. The obtained solutions are radially stable and the zero-mode arising from the stability equation can be written analytically.


title: A highly efficient tensor network algorithm for multi-asset Fourier options pricing

abstract: Risk assessment and in particular derivatives pricing is one of the core areas in computational finance and accounts for a sizeable fraction of the global computing resources of the financial industry. We outline a quantum-inspired algorithm for multi-asset options pricing. The algorithm is based on tensor networks, which have allowed for major conceptual and numerical breakthroughs in quantum many body physics and quantum computation. In the proof-of-concept example explored, the tensor network approach yields several orders of magnitude speedup over vanilla Monte Carlo simulations. We take this as good evidence that the use of tensor network methods holds great promise for alleviating the computation burden of risk evaluation in the financial and other industries, thus potentially lowering the carbon footprint these simulations incur today.


title: Stability Constrained Characterization of the 23 Myr-old V1298 Tau System: Do Young Planets Form in Mean Motion Resonance Chains?

abstract: A leading theoretical expectation for the final stages of planet formation is that disk migration should naturally drive orbits into chains of mean motion resonances (MMRs). In order to explain the dearth of MMR chains observed at Gyr ages ($<1\%$), this picture requires such configurations to destabilize and scramble period ratios following disk dispersal. Strikingly, the only two known stars with three or more planets younger than $\lesssim 100$ Myrs, HR 8799 and V1298 Tau, have been suggested to be in such MMR chains, given the orbits' near-integer period ratios. We incorporate recent transit and radial velocity observations of the V1298 Tau system, and investigate constraints on the system's orbital architecture imposed by requiring dynamical stability on timescales much shorter than the system's age. We show that the recent radial-velocity mass measurement of V1298 Tau $b$ places it within a factor of two of the instability limit, and that this allows us to set significantly lower limits on the eccentricity ($e_b \leq 0.17$ at $99.7\%$ confidence). Additionally, we rule out a resonant chain configuration for V1298 Tau at $\gtrsim 99\%$ confidence. Thus, if the $\sim 23$ Myr-old V1298 Tau system did form as a resonant chain, it must have undergone instability and rearrangement shortly after disk dispersal. We expect that similar stability constrained characterization of future young multi-planet systems will be valuable in informing planet formation models.


title: Converging approximations of attractors via almost Lyapunov functions and semidefinite programming

abstract: In this paper we combine two existing approaches for approximating attractors. One of them approximates the attractors arbitrarily well by sublevel sets related to solutions of infinite dimensional linear programming problems. A downside there is that these sets are not necessarily positively invariant. On the contrary, the second method provides supersets of the attractor which are positively invariant. Their method on the other hand has the disadvantage that the underlying optimization problem is not computationally tractable without the use of heuristics - and incorporating them comes at the price of losing guaranteed convergence. In this paper we marry both approaches by combining their techniques and we get converging outer approximations of the attractor consisting of positively invariant sets based on convex optimization via sum-of-squares techniques. The method is easy to use and illustrated by numerical examples.


title: Off-Policy Evaluation in Embedded Spaces

abstract: Off-policy evaluation methods are important in recommendation systems and search engines, where data collected under an existing logging policy is used to estimate the performance of a new proposed policy. A common approach to this problem is weighting, where data is weighted by a density ratio between the probability of actions given contexts in the target and logged policies. In practice, two issues often arise. First, many problems have very large action spaces and we may not observe rewards for most actions, and so in finite samples we may encounter a positivity violation. Second, many recommendation systems are not probabilistic and so having access to logging and target policy densities may not be feasible. To address these issues, we introduce the featurized embedded permutation weighting estimator. The estimator computes the density ratio in an action embedding space, which reduces the possibility of positivity violations. The density ratio is computed leveraging recent advances in normalizing flows and density ratio estimation as a classification problem, in order to obtain estimates which are feasible in practice.


title: Unveiling and Manipulating Hidden Symmetries in Graphene Nanoribbons

abstract: Armchair graphene nanoribbons are a highly promising class of semiconductors for all-carbon nanocircuitry. Here, we present a new perspective on their electronic structure from simple model Hamiltonians and $\textit{ab initio}$ calculations. We focus on a specific set of nanoribbons of width $n = 3p+2$, where $n$ is the number of carbon atoms across the nanoribbon axis and $p$ is a positive integer. We demonstrate that the energy-gap opening in these nanoribbons originates from the breaking of a previously unidentified hidden symmetry by long-ranged hopping of $\pi$-electrons and structural distortions occurring at the edges. This hidden symmetry can be restored or manipulated through the application of in-plane lattice strain, which enables continuous energy-gap tuning, the emergence of Dirac points at the Fermi level, and topological quantum phase transitions. Our work establishes an original interpretation of the semiconducting character of armchair graphene nanoribbons and offers guidelines for rationally designing their electronic structure.


title: 28 GHz Phased Array-Based Self-Interference Measurements for Millimeter Wave Full-Duplex

abstract: We present measurements of the 28 GHz self-interference channel for full-duplex sectorized multi-panel millimeter wave (mmWave) systems, such as integrated access and backhaul. We measure the isolation between the input of a transmitting phased array panel and the output of a co-located receiving phased array panel, each of which is electronically steered across a number of directions in azimuth and elevation. In total, nearly 6.5 million measurements were taken in an anechoic chamber to densely inspect the directional nature of the coupling between 256-element phased arrays. We observe that highly directional mmWave beams do not necessarily offer widespread high isolation between transmitting and receiving arrays. Rather, our measurements indicate that steering the transmitter or receiver away from the other tends to offer higher isolation but even slight steering changes can lead to drastic variations in isolation. These measurements can be useful references when developing mmWave full-duplex solutions and can motivate a variety of future topics including beam/user selection and beamforming codebook design.


title: Virtual Reality Digital Twin and Environment for Troubleshooting Lunar-based Infrastructure Assembly Failures

abstract: Humans and robots will need to collaborate in order to create a sustainable human lunar presence by the end of the 2020s. This includes cases in which a human will be required to teleoperate an autonomous rover that has encountered an instrument assembly failure. To aid teleoperators in the troubleshooting process, we propose a virtual reality digital twin placed in a simulated environment. Here, the operator can virtually interact with a digital version of the rover and mechanical arm that uses the same controls and kinematic model. The user can also adopt the egocentric (a first person view through using stereoscopic passthrough) and exocentric (a third person view where the operator can virtually walk around the environment and rover as if they were on site) view. We also discuss our metrics for evaluating the differences between our digital and physical robot, as well as the experimental concept based on real and applicable missions, and future work that would compare our platform to traditional troubleshooting methods.


title: Temperature Distributions in the Protosatellite Disks of Uranus, Jupiter, Neptune and Saturn and Common Structure in Their Satellite Systems

abstract: We observe similar structures in the orbital radii of satellite and ring systems of Uranus, Jupiter and Neptune. This stimulates an investigation that the evolution of these systems and Saturn's system is initiated by a common mechanism that involves the interaction of radiation with their subnebulae (protosatellite disks). A model is presented that is characterized by resonance created through stimulated radiative molecular association (SRMA) reactions. In this model thermal energy is extracted from a protoplanetary disk at specific distances from the protoplanet wherever there is a match between the local thermal energy of the disk and the energy of photons impinging on the disk. Radiation is supplied by a portion of the hydrogen molecule's spectrum for the Uranian, Jovian and Neptunian disks and a portion of the hydrogen atom's spectrum for the Saturnian disk. Findings shed light on the early evolution and structure of satellite systems including the complicated ring system of Saturn. They also link satellite orbital radii size distributions to shapes of temperature distributions (TD's) in protosatellite disks. Theoretically determined TD's in the protosatellite disks of Uranus (Mousis 2004) and Saturn (Mousis Gautier and Bockelee-Moran, 2002) are essential to the present investigation.


title: Partially polaron-transformed quantum master equation for exciton and charge transport dynamics

abstract: Polaron-transformed quantum master equation (PQME) offers a unified framework to describe the dynamics of quantum systems in both limits of weak and strong couplings to environmental degrees of freedom. Thus, PQME serves as an efficient method to describe charge and exciton transfer/transport dynamics for a broad range of parameters in condensed or complex environments. However, in some cases, the polaron transformation (PT) being employed in the formulation invokes an over-relaxation of slow modes and results in premature suppression of important coherence terms. A formal framework to address this issue is developed in the present work by employing a partial PT that has smaller weights for low frequency bath modes. It is shown here that a closed form expression of a 2nd order time-local PQME including all the inhomogeneous terms can be derived for a general form of partial PT, although more complicated than that for the full PT. All the expressions needed for numerical calculation are derived in detail. Applications to a model of two-level system coupled to a bath of harmonic oscillators, with test calculations focused on those due to homogeneous relaxation terms, demonstrate the feasibility and the utility of the present approach.

categories: quant-ph cond-mat.mes-hall physics.bio-ph physics.chem-ph physics.comp-ph


title: Explicit expressions for \vSapovalov elements in Type A

abstract: We give explicit expressions for \vSapovalov elements in Type A Lie algebras and superalgebras. Explicit expressions were already given in arXiv:1710.10528 Section 9, using non-commutative determinants, and in fact our first main results, Theorems 2.3 and 2.6 can be viewed as complete expansions of these determinants. But we give new proofs, which seem easier because they avoid induction and cofactor expansion. We also describe \vSapovalov elements for fgl(m,n)with respect to an arbitrary Borel subalgebra in Theorem 3.7, and interpret \vSapovalov elements in Type A as determinants of Hessenberg matrices in Theorems 4.5 and 4.9.The exact form of the explicit expressions depends on an ordering on the set of positive roots, and Hessenberg matrices are useful in changing the ordering. Having explicit expressions for \vSapovalov elements allows us to give easy proofs of several results in representation theory, see Section 6 and Subsection 4.6.


title: Competitive Exclusion in an Artificial Foraging Ecosystem

abstract: Artificial ecosystems provide an additional experimental tool to support laboratory work, field work, and theoretical development in competitive exclusion research. A novel application of a spatiotemporal agent based model is presented which simulates two foraging species of different intrinsic growth rates competing for the same replenishing resource in stable and seasonal environments. This experimental approach provides precise control over the parameters for the environment, the species, and individual movements. Detailed trajectories of these non-equilibrium populations and their characteristics are produced. Narrow zones of potential coexistence are identified within the environmental and intrinsic growth parameter spaces. An example of commensalism driven by the local spatial dynamics is identified. Results of these experiments are discussed in context of modern coexistence theory and research in movement-mediated community assembly. Constraints on possible origination scenarios are identified.


title: Double Choco is NP-complete

abstract: In the Nikoli pencil-and-paper game Double Choco, a puzzle consists of an m $\times$ n grid of cells of white or gray color, separated by dotted lines where each cell possibly contains an integer. The goal is to partition the grid into blocks by drawing solid lines over the dotted lines, where every block must contain a pair of areas of white and gray cells having the same form (size and shape). An integer indicates the number of cells of that color in the block. A block can contain any number of cells with the integer. We prove this puzzle NP-complete, establishing a Nikoli gap of 2 years.


title: Safely: Safe Stochastic Motion Planning Under Constrained Sensing via Duality

abstract: Consider a robot operating in an uncertain environment with stochastic, dynamic obstacles. Despite the clear benefits for trajectory optimization, it is often hard to keep track of each obstacle at every time step due to sensing and hardware limitations. We introduce the Safely motion planner, a receding-horizon control framework, that simultaneously synthesizes both a trajectory for the robot to follow as well as a sensor selection strategy that prescribes trajectory-relevant obstacles to measure at each time step while respecting the sensing constraints of the robot. We perform the motion planning using sequential quadratic programming, and prescribe obstacles to sense based on the duality information associated with the convex subproblems. We guarantee safety by ensuring that the probability of the robot colliding with any of the obstacles is below a prescribed threshold at every time step of the planned robot trajectory. We demonstrate the efficacy of the Safely motion planner through software and hardware experiments.


title: Asymptotic symmetries and soft charges of fractons

abstract: The asymptotic structure of gauge theories describing fracton interactions is analyzed. Two sets of asymptotic conditions are proposed. Both encompass all known solutions, lead to finite charges and resolve the problem of the divergent energy coming from the monopole contribution. While the first set leads to the expected fracton symmetry algebra, including a dipole charge, the second set provides a soft infinite-dimensional extension of it. These soft charges provide evidence of a rich infrared structure for fracton-like theories and provide one corner of a possible fracton infrared triangle.

categories: hep-th cond-mat.stat-mech cond-mat.str-el gr-qc


title: Fuzzy Forests For Feature Selection in High-Dimensional Survey Data: An Application to the 2020 U.S. Presidential Election

abstract: An increasingly common methodological issue in the field of social science is high-dimensional and highly correlated datasets that are unamenable to the traditional deductive framework of study. Analysis of candidate choice in the 2020 Presidential Election is one area in which this issue presents itself: in order to test the many theories explaining the outcome of the election, it is necessary to use data such as the 2020 Cooperative Election Study Common Content, with hundreds of highly correlated features. We present the Fuzzy Forests algorithm, a variant of the popular Random Forests ensemble method, as an efficient way to reduce the feature space in such cases with minimal bias, while also maintaining predictive performance on par with common algorithms like Random Forests and logit. Using Fuzzy Forests, we isolate the top correlates of candidate choice and find that partisan polarization was the strongest factor driving the 2020 presidential election.


title: Workflow for investigating thermodynamic, structural and energy properties of condensed polymer systems from Molecular Dynamics

abstract: Soft matter materials and polymers are widely used in the controlled delivery of drugs. Simulation and modeling provide insight at the atomic scale enabling a level of control unavailable to experiments. We present a workflow protocol for modeling, simulating, and analyzing structural and thermodynamic response properties of poly-lactic-co-glycolic acid (PLGA), a well-studied and FDA approved material. We concatenate a battery of molecular dynamics, computational chemistry, highly parallel scripting, and analysis tools for generating properties of bulk polymers in the condensed phase. We provide the workflow leading to the glass transition temperature, enthalpy, density, isobaric heat capacity, thermal expansion coefficient, isothermal compressibility, bulk modulus, sonic velocity, cohesive energy, and solubility parameters. Calculated properties agree very well with experiments, when available. This methodology has been extended to a variety of polymer types and environments.


title: Evaluation of Dirichlet Process Gaussian Mixtures for Segmentation on Noisy Hyperspectral Images

abstract: Image segmentation is a fundamental step for the interpretation of Remote Sensing Images. Clustering or segmentation methods usually precede the classification task and are used as support tools for manual labeling. The most common algorithms, such as k-means, mean-shift, and MRS, require an extra manual step to find the scale parameter. The segmentation results are severely affected if the parameters are not correctly tuned and diverge from the optimal values. Additionally, the search for the optimal scale is a costly task, as it requires a comprehensive hyper-parameter search. This paper proposes and evaluates a method for segmentation of Hyperspectral Images using the Dirichlet Process Gaussian Mixture Model. Our model can self-regulate the parameters until it finds the optimal values of scale and the number of clusters in a given dataset. The results demonstrate the potential of our method to find objects in a Hyperspectral Image while bypassing the burden of manual search of the optimal parameters. In addition, our model also produces similar results on noisy datasets, while previous research usually required a pre-processing task for noise reduction and spectral smoothing.


title: Bathymetry Inversion using a Deep-Learning-Based Surrogate for Shallow Water Equations Solvers

abstract: River bathymetry is critical for many aspects of water resources management. We propose and demonstrate a bathymetry inversion method using a deep-learning-based surrogate for shallow water equations solvers. The surrogate uses the convolutional autoencoder with a shared-encoder, separate-decoder architecture. It encodes the input bathymetry and decodes to separate outputs for flow-field variables. A gradient-based optimizer is used to perform bathymetry inversion with the trained surrogate. Two physically-based constraints on both bed elevation value and slope have to be added as inversion loss regularizations to obtain usable inversion results. Using the "L-curve" criterion, a heuristic approach was proposed to determine the regularization parameters. Both the surrogate model and the inversion algorithm show good performance. We found the bathymetry inversion process has two distinctive stages, which resembles the sculptural process of initial broad-brush calving and final detailing. The inversion loss due to flow prediction error reaches its minimum in the first stage and remains almost constant afterward. The bed elevation value and slope regularizations play the dominant role in the second stage in selecting the most probable solution. We also found the surrogate architecture (whether with both velocity and water surface elevation or velocity only as outputs) does not show significant impact on inversion result.


title: Phase trajectories of linear homogeneous autonomous dynamical systems of the third order

abstract: The proposed article is devoted to the study of the problem of constructing phase trajectories in the vicinity of a singular point. This paper presents a more expanded view of this problem in comparison with those previously considered by other authors and suggests a more meaningful and grounded development of this study. Numerous examples of the application of this approach to the study of trajectories at a singular point are given.


title: Improved Hamiltonians for Quantum Simulations

abstract: Quantum simulations of lattice gauge theories for the foreseeable future will be hampered by limited resources. The historical success of improved lattice actions in classical simulations strongly suggests that Hamiltonians with improved discretization errors will reduce quantum resources, i.e. require $\gtrsim 2^d$ fewer qubits in quantum simulations for lattices with $d$ spatial dimensions. In this work, we consider $\mathcal{O}(a^2)$-improved Hamiltonians for pure gauge theories and design the corresponding quantum circuits for its real-time evolution in terms of primitive gates. An explicit demonstration for $\mathbb{Z}_2$ gauge theory is presented including exploratory tests using the ibm_perth device.


title: Distributional Hardness Against Preconditioned Lasso via Erasure-Robust Designs

abstract: Sparse linear regression with ill-conditioned Gaussian random designs is widely believed to exhibit a statistical/computational gap, but there is surprisingly little formal evidence for this belief, even in the form of examples that are hard for restricted classes of algorithms. Recent work has shown that, for certain covariance matrices, the broad class of Preconditioned Lasso programs provably cannot succeed on polylogarithmically sparse signals with a sublinear number of samples. However, this lower bound only shows that for every preconditioner, there exists at least one signal that it fails to recover successfully. This leaves open the possibility that, for example, trying multiple different preconditioners solves every sparse linear regression problem. In this work, we prove a stronger lower bound that overcomes this issue. For an appropriate covariance matrix, we construct a single signal distribution on which any invertibly-preconditioned Lasso program fails with high probability, unless it receives a linear number of samples. Surprisingly, at the heart of our lower bound is a new positive result in compressed sensing. We show that standard sparse random designs are with high probability robust to adversarial measurement erasures, in the sense that if $b$ measurements are erased, then all but $O(b)$ of the coordinates of the signal are still information-theoretically identifiable. To our knowledge, this is the first time that partial recoverability of arbitrary sparse signals under erasures has been studied in compressed sensing.

categories: cs.DS cs.IT cs.LG math.IT math.ST stat.ML stat.TH


title: Almost K\"ahler metrics and pp-wave spacetimes

abstract: We establish a one-to-one correspondence between a class of strictly almost K\"ahler metrics on the one hand, and Lorentzian pp-wave spacetimes on the other; the latter metrics are well known in general relativity, where they model radiation propagating at the speed of light. Specifically, we construct families of complete almost K\"ahler metrics by deforming pp-waves via their propagation wave vector. The almost K\"ahler metrics we obtain exist in all dimensions $2n \geq 4$, and are defined on both $\mathbb{R}^{2n}$ and $\mathbb{S}^1\times\mathbb{S}^1 \times M$, where $M$ is any closed almost K\"ahler manifold; they are not warped products, they include noncompact examples with constant negative scalar curvature, and all of them have the property that their fundamental 2-forms are also co-closed with respect to the Lorentzian pp-wave metric. Finally, we further deepen this relationship between almost K\"ahler and Lorentzian geometry by utilizing Penrose's "plane wave limit," by which every spacetime has, locally, a pp-wave metric as a limit: using Penrose's construction, we show that in all dimensions $2n \geq 4$, every Lorentzian metric admits, locally, an almost K\"ahler metric of this form as a limit.


title: On the Maximum $F_5$-free Subhypergraphs of a Random Hypergraph

abstract: Denote by $F_5$ the $3$-uniform hypergraph on vertex set $\{1,2,3,4,5\}$ with hyperedges $\{123,124,345\}$. Balogh, Butterfield, Hu, and Lenz proved that if $p > K \log n / n$ for some large constant $K$, then every maximum $F_5$-free subhypergraph of $G^3(n,p)$ is tripartite with high probability, and showed that if $p_0 = 0.1\sqrt{\log n} / n$, then with high probability there exists a maximum $F_5$-free subhypergraph of $G^3(n,p_0)$ that is not tripartite. In this paper, we sharpen the upper bound to be best possible up to a constant factor. We prove that if $p > C \sqrt{\log n} / n $ for some large constant $C$, then every maximum $F_5$-free subhypergraph of $G^3(n, p)$ is tripartite with high probability.


title: Data-driven input reconstruction and experimental validation

abstract: This paper addresses a data-driven input reconstruction problem based on Willems' Fundamental Lemma in which unknown input estimators (UIEs) are constructed directly from historical I/O data. Given only output measurements, the inputs are estimated by the UIE, which is shown to asymptotically converge to the true input without knowing the initial conditions. Both open-loop and closed-loop UIEs are developed based on Lyapunov conditions and the Luenberger-observer-type feedback, whose convergence properties are studied. An experimental study is presented demonstrating the efficacy of the closed-loop UIE for estimating the occupancy of a building on the EPFL campus via measured carbon dioxide levels.


title: Universal anomalous turbulent diffusion in quantum fluids

abstract: In classical viscous fluids, turbulent eddies are known to be responsible for the rapid spreading of embedded particles. But in an inviscid quantum fluid where the turbulence is induced by a chaotic tangle of quantized vortices, dispersion of the particles is achieved via a non-classical mechanism, i.e., their binding to the evolving quantized vortices. However, there is limited existing knowledge on how the vortices diffuse and spread in turbulent quantum fluids. Here we report a systematic numerical study of the apparent diffusion of vortices in a random vortex tangle in superfluid helium-4 using full Biot-Savart simulation. We reveal that the vortices in pure superfluid exhibit a universal anomalous diffusion (superdiffusion) at small times, which transits to normal diffusion at large times. This behavior is found to be caused by a generic scaling property of the vortex velocity, which should exist in all quantum fluids where the Biot-Savart law governs the vortex motion. Our simulation at finite temperatures also nicely reproduces recent experimental observations. The knowledge obtained from this study may form the foundation for understanding turbulent transport and universal vortex dynamics in various condensed-matter and cosmic quantum fluids.

categories: physics.flu-dyn cond-mat.mtrl-sci cond-mat.other nlin.CD


title: The limit cycles in a generalized Rayleigh-Li\'enard oscillator

abstract: We compute the cyclicity of open period annuli of the following generalized Rayleigh-Li\'enard equation $$\ddot{x}+ax+bx^3-(\lambda_1+\lambda_2 x^2+\lambda_3\dot{x}^2+\lambda_4 x^4+\lambda_5\dot{x}^4+\lambda_6 x^6)\dot{x}=0$$ and the equivalent planar system $X_\lambda$, where the coefficients of the perturbation $\lambda_j$ are independent small parameters and $a, b$ are fixed nonzero constants. Our main tool is the machinery of the so called higher-order Poincar\'e-Pontryagin-Melnikov functions (Melnikov functions $M_n$ for short), combined with the explicit computation of center conditions and the corresponding Bautin ideal. We consider first arbitrary analytic arcs $\varepsilon \to \lambda(\varepsilon)$ and explicitly compute all possible Melnikov functions $M_n$ related to the deformation $X_{ \lambda(\varepsilon)} $. At a second step we obtain exact bounds for the number of the zeros of the Melnikov functions (complete elliptic integrals depending on parameter) in an appropriate complex domain, using a modification of Petrov's method. To deal with the general case of six-parameter deformations $\lambda \to X_\lambda$, we compute first the related Bautin ideal. To do this we carefully study the Melnikov functions up to order three, and then use Nakayama lemma from Algebraic geometry. The principalization of the Bautin ideal (achieved after a blow up) reduces finally the study of general deformations $X_{ \lambda } $ to the study of one-parameter deformations $X_{ \lambda(\varepsilon)} $.


title: On the surface plasmonic waves excited by a dipole above anisotropic and spatially dispersive two-dimensional surfaces of infinite extent in planarly layered media

abstract: The surface plasmonic waves excited by a vertical or horizontal oriented Hertzian dipole above anisotropic and spatially dispersive two-dimensional surfaces of infinite extent embedded in planarly layered uniaxial media is investigated using the dyadic Green function approach. The spectral-domain transmission line analogy Green function formulation and iso-frequency contours equations are derived. The methods to accurately and efficently evaluate the two-dimensional Fourier integral arisen from the spatial-domain Green function computation are also developed. To resolve the numerical inefficiency due to the highly oscillatory integrand and singularities of surface plasmonic waves possessing large wavenumber, two numerical strategies, the extrapolation of the real-axis integration combined with singularity subtraction, and the deformed vertical integration path, are proposed and applicable to a wide range of observation distance. As a demonstration of the proposed formulation, we compute the scattered fields of a vertical dipole above the graphene biased by drift current which exhibits significant spatial dispersion and show that its light-matter interaction can be significantly reinforced when placed above uniaxially epsilon-near-zero substrates. The proposed formulation may provide methodology for the computational analysis of two-dimensional materials and surface plasmonic waves.


title: The detection matrix as a model-agnostic tool to estimate the number of degrees of freedom in mechanical systems and engineering structures

abstract: Estimating the number of degrees of freedom of a mechanical system or an engineering structure from the time-series of a small set of sensors is a basic problem in diagnostics, which, however, is often overlooked when monitoring their health and integrity. In this work, we demonstrate the applicability of the network-theoretic concept of detection matrix as a tool to solve this problem. From this estimation, we illustrate the possibility to identify damage. The detection matrix, recently introduced by Haehne et al. in the context of network theory, is assembled from the transient response of a few nodes as a result of non-zero initial conditions: its rank offers an estimate of the number of nodes in the network itself. The use of the detection matrix is completely model-agnostic, whereby it does not require any knowledge of the system dynamics. Here, we show that, with a few modifications, this same principle applies to discrete systems, such as spring-mass lattices and trusses. Moreover, we discuss how damage in one or more members causes the appearance of distinct jumps in the singular values of this matrix, thereby opening the door to structural health monitoring applications, without the need for a complete model reconstruction.


title: On the Error of Random Sampling: Uniformly Distributed Random Points on Parametric Curves

abstract: Given a parametric polynomial curve $\gamma:[a,b]\rightarrow \mathbb{R}^n$, how can we sample a random point $\mathfrak{x}\in \mathrm{im}(\gamma)$ in such a way that it is distributed uniformly with respect to the arc-length? Unfortunately, we cannot sample exactly such a point-even assuming we can perform exact arithmetic operations. So we end up with the following question: how does the method we choose affect the quality of the approximate sample we obtain? In practice, there are many answers. However, in theory, there are still gaps in our understanding. In this paper, we address this question from the point of view of complexity theory, providing bounds in terms of the size of the desired error.


title: Tabula: Efficiently Computing Nonlinear Activation Functions for Secure Neural Network Inference

abstract: Multiparty computation approaches to secure neural network inference commonly rely on garbled circuits for securely executing nonlinear activation functions. However, garbled circuits require excessive communication between server and client, impose significant storage overheads, and incur large runtime penalties. To reduce these costs, we propose an alternative to garbled circuits: Tabula, an algorithm based on secure lookup tables. Our approach precomputes lookup tables during an offline phase that contains the result of all possible nonlinear function calls. Because these tables incur exponential storage costs in the number of operands and the precision of the input values, we use quantization to reduce these storage costs to make this approach practical. This enables an online phase where securely computing the result of a nonlinear function requires just a single round of communication, with communication cost equal to twice the number of bits of the input to the nonlinear function. In practice our approach costs 2 bytes of communication per nonlinear function call in the online phase. Compared to garbled circuits with 8-bit quantized inputs, when computing individual nonlinear functions during the online phase, experiments show Tabula with 8-bit activations uses between $280$-$560 \times$ less communication, is over $100\times$ faster, and uses a comparable (within a factor of 2) amount of storage; compared against other state-of-the-art protocols Tabula achieves greater than $40\times$ communication reduction. This leads to significant performance gains over garbled circuits with quantized inputs during the online phase of secure inference of neural networks: Tabula reduces end-to-end inference communication by up to $9 \times$ and achieves an end-to-end inference speedup of up to $50 \times$, while imposing comparable storage and offline preprocessing costs.


title: Future role and economic benefits of hydrogen and synthetic energy carriers in Germany: a systematic review of long-term energy scenarios

abstract: Determining the development of Germany's energy system by taking the energy transition objectives into account is the subject of a series of studies. Since their assumptions and results play a significant role in the political energy debate for understanding the role of hydrogen and synthetic energy carriers, a better discussion is needed. This article provides a comparative assessment of published transition pathways for Germany to assess the role and advantages of hydrogen and synthetic energy carriers. Twelve energy studies were selected and 37 scenarios for the years 2030 and 2050 were evaluated. Despite the variations, the two carriers will play an important future role. While their deployment is expected to have only started by 2030 with a mean demand of 91 TWh/a (4% of the final energy demand) in Germany, they will be an essential part by 2050 with a mean demand of 480 TWh/a (24% of the final energy demand). A moderately positive correlation (0.53) between the decarbonisation targets and the share of hydrogen-based carriers in final energy demand underlines the relevance for reaching the climate targets. Additionally, value creation effects of about 5 bn EUR/a in 2030 can be expected for hydrogen-based carriers. By 2050, these effects will increase to almost 16 bn EUR/a. Hydrogen is expected to be mainly produced domestically while synthetic fuels are projected to be mostly imported. Despite of all the advantages, the construction of the facilities is associated with high costs which should be not neglected in the discussion.


title: Systematic, Lyapunov-Based, Safe and Stabilizing Controller Synthesis for Constrained Nonlinear Systems

abstract: A controller synthesis method for state- and input-constrained nonlinear systems is presented that seeks continuous piecewise affine (CPA) Lyapunov-like functions and controllers simultaneously. Non-convex optimization problems are formulated on triangulated subsets of the admissible states that can be refined to meet primary control objectives, such as stability and safety, alongside secondary performance objectives. A multi-stage design is also given that enlarges the region of attraction (ROA) sequentially while allowing exclusive performance for each stage. A clear boundary for an invariant subset of closed-loop system's ROA is obtained from the resulting Lipschitz Lyapunov function. For control-affine nonlinear systems, the non-convex problem is formulated as a series of conservative, but well-posed, semi-definite programs. These decrease the cost function iteratively until the design objectives are met. Since the resulting CPA Lyapunov-like functions are also Lipschitz control (or barrier) Lyapunov functions, they can be used in online quadratic programming to find minimum-norm control inputs. Numerical examples are provided to demonstrate the effectiveness of the method.