title: The Rashba Scale: Emergence of Band Anti-Crossing as a Design Principle for Materials with Large Rashba coefficient

abstract: The spin-orbit -induced spin splitting of energy bands in low symmetry compounds (the Rashba Effect) has a long-standing relevance to spintronic applications and to the fundamental understanding of symmetry breaking in solids, yet the knowledge of what controls its magnitude in different materials is difficult to anticipate. Indeed, rare discoveries of compounds with large Rashba coefficients are invariably greeted as pleasant surprises. We advance the understanding of the "Rashba Scale" using the "inverse design" approach by formulating theoretically the relevant design principle and then identifying compounds that satisfy it. We show that the presence of energy band anti-crossing provides a causal design principle of compounds with large Rashba coefficients, leading to the identification via first-principles calculations of 34 rationally designed strong-Rashba compounds. Since topological insulators must have band anti crossing, this leads us to establish an interesting cross functionality of "Topological Rashba Insulators" (TRI) that may provide a platform for the simultaneous control of spin splitting and spin-polarization.


title: Automatic Generation of Algorithms for Black-Box Robust Optimisation Problems

abstract: We develop algorithms capable of tackling robust black-box optimisation problems, where the number of model runs is limited. When a desired solution cannot be implemented exactly the aim is to find a robust one, where the worst case in an uncertainty neighbourhood around a solution still performs well. This requires a local maximisation within a global minimisation. To investigate improved optimisation methods for robust problems, and remove the need to manually determine an effective heuristic and parameter settings, we employ an automatic generation of algorithms approach: Grammar-Guided Genetic Programming. We develop algorithmic building blocks to be implemented in a Particle Swarm Optimisation framework, define the rules for constructing heuristics from these components, and evolve populations of search algorithms. Our algorithmic building blocks combine elements of existing techniques and new features, resulting in the investigation of a novel heuristic solution space. As a result of this evolutionary process we obtain algorithms which improve upon the current state of the art. We also analyse the component level breakdowns of the populations of algorithms developed against their performance, to identify high-performing heuristic components for robust problems.


title: Understanding high pressure hydrogen with a hierarchical machine-learned potential

abstract: The hydrogen phase diagram has a number of unusual features which are generally well reproduced by density functional calculations. Unfortunately, these calculations fail to provide good physical insights into why those features occur. In this paper, we parameterize a model potential for molecular hydrogen which permits long and large simulations. The model shows excellent reproduction of the phase diagram, including the broken-symmetry Phase II, an efficiently-packed phase III and the maximum in the melt curve. It also gives an excellent reproduction of the vibrational frequencies, including the maximum in the vibrational frequency $\nu(P)$ and negative thermal expansion. By detailed study of lengthy molecular dynamics, we give intuitive explanations for observed and calculated properties. All solid structures approximate to hexagonal close packed, with symmetry broken by molecular orientation. At high pressure, Phase I shows significant short-ranged correlations between molecular orientations. The turnover in Raman frequency is due to increased coupling between neighboring molecules, rather than weakening of the bond. The liquid is denser than the close-packed solid because, at molecular separations below 2.3\AA, the favoured relative orientation switches from quadrupole-energy-minimising to steric-repulsion-minimising. The latter allows molecules to get closer together, without atoms getting closer but this cannot be achieved within the constraints of a close-packed layer.


title: Clustering Time Series Data through Autoencoder-based Deep Learning Models

abstract: Machine learning and in particular deep learning algorithms are the emerging approaches to data analysis. These techniques have transformed traditional data mining-based analysis radically into a learning-based model in which existing data sets along with their cluster labels (i.e., train set) are learned to build a supervised learning model and predict the cluster labels of unseen data (i.e., test set). In particular, deep learning techniques are capable of capturing and learning hidden features in a given data sets and thus building a more accurate prediction model for clustering and labeling problem. However, the major problem is that time series data are often unlabeled and thus supervised learning-based deep learning algorithms cannot be directly adapted to solve the clustering problems for these special and complex types of data sets. To address this problem, this paper introduces a two-stage method for clustering time series data. First, a novel technique is introduced to utilize the characteristics (e.g., volatility) of given time series data in order to create labels and thus be able to transform the problem from unsupervised learning into supervised learning. Second, an autoencoder-based deep learning model is built to learn and model both known and hidden features of time series data along with their created labels to predict the labels of unseen time series data. The paper reports a case study in which financial and stock time series data of selected 70 stock indices are clustered into distinct groups using the introduced two-stage procedure. The results show that the proposed procedure is capable of achieving 87.5\% accuracy in clustering and predicting the labels for unseen time series data.


title: Secure protocol to protect location privacy in distance calculation

abstract: Several applications require computing distances between different people. For example,this is required if we want to obtain the close contacts of people in case of and epidemic,or when restraining orders are imposed by a judge. However, periodically revealing location might pose a privacy threat to the involved parties. Continuous location data may be used to infer personal information about the owner, like behaviors, religious beliefs, buying habits, routines, etc. In this paper, we show that it is possible to calculate distance between two parties without disclosing their latitude and longitude data. For this purpose, we design a secure protocol based on the ElGamal cryptosystem and its homomorphic properties. The proposed protocol allows the calculation of distances while preserving location privacy. The protocol is analyzed in terms of security and performance. The security analysis shows that no involved party can learn any information about location.


title: Mixing properties of generalized $T, T^{-1}$ transformations

abstract: We study mixing properties of generalized $T, T^{-1}$ transformation. We discuss two mixing mechanisms. In the case the fiber dynamics is mixing, it is sufficient that the driving cocycle is small with small probability. In the case the fiber dynamics is only assumed to be ergodic, one needs to use the shearing properties of the cocycle. Applications include the Central Limit Theorem for sufficiently fast mixing systems and the estimates on deviations of ergodic averages.


title: Asymptotic realization of the super-BMS algebra at spatial infinity

abstract: Explicit boundary conditions are given at spatial infinity for four-dimensional supergravity, which provide a realization of the super-BMS algebra of Awada, Gibbons and Shaw. The results are then generalized to the $N$- extended super-BMS algebras.


title: Gumbel-softmax-based Optimization: A Simple General Framework for Optimization Problems on Graphs

abstract: In computer science, there exist a large number of optimization problems defined on graphs, that is to find a best node state configuration or a network structure such that the designed objective function is optimized under some constraints. However, these problems are notorious for their hardness to solve because most of them are NP-hard or NP-complete. Although traditional general methods such as simulated annealing (SA), genetic algorithms (GA) and so forth have been devised to these hard problems, their accuracy and time consumption are not satisfying in practice. In this work, we proposed a simple, fast, and general algorithm framework based on advanced automatic differentiation technique empowered by deep learning frameworks. By introducing Gumbel-softmax technique, we can optimize the objective function directly by gradient descent algorithm regardless of the discrete nature of variables. We also introduce evolution strategy to parallel version of our algorithm. We test our algorithm on three representative optimization problems on graph including modularity optimization from network science, Sherrington-Kirkpatrick (SK) model from statistical physics, maximum independent set (MIS) and minimum vertex cover (MVC) problem from combinatorial optimization on graph. High-quality solutions can be obtained with much less time consuming compared to traditional approaches.


title: ESResNet: Environmental Sound Classification Based on Visual Domain Models

abstract: Environmental Sound Classification (ESC) is an active research area in the audio domain and has seen a lot of progress in the past years. However, many of the existing approaches achieve high accuracy by relying on domain-specific features and architectures, making it harder to benefit from advances in other fields (e.g., the image domain). Additionally, some of the past successes have been attributed to a discrepancy of how results are evaluated (i.e., on unofficial splits of the UrbanSound8K (US8K) dataset), distorting the overall progression of the field. The contribution of this paper is twofold. First, we present a model that is inherently compatible with mono and stereo sound inputs. Our model is based on simple log-power Short-Time Fourier Transform (STFT) spectrograms and combines them with several well-known approaches from the image domain (i.e., ResNet, Siamese-like networks and attention). We investigate the influence of cross-domain pre-training, architectural changes, and evaluate our model on standard datasets. We find that our model out-performs all previously known approaches in a fair comparison by achieving accuracies of 97.0 % (ESC-10), 91.5 % (ESC-50) and 84.2 % / 85.4 % (US8K mono / stereo). Second, we provide a comprehensive overview of the actual state of the field, by differentiating several previously reported results on the US8K dataset between official or unofficial splits. For better reproducibility, our code (including any re-implementations) is made available.


title: Uniqueness Criteria for the Oseen Vortex in the 3d Navier-Stokes Equations

abstract: In this paper, we consider the uniqueness of solutions to the 3d Navier-Stokes equations with initial vorticity given by $\omega_0 = \alpha e_z \delta_{x = y = 0}$, where $\delta_{x=y= 0}$ is the one dimensional Hausdorff measure of an infinite, vertical line and $\alpha \in \mathbb R$ is an arbitrary circulation. This initial data corresponds to an idealized, infinite vortex filament. One smooth, mild solution is given by the self-similar Oseen vortex column, which coincides with the heat evolution. Previous work by Germain, Harrop-Griffiths, and the first author implies that this solution is unique within a class of mild solutions that converge to the Oseen vortex in suitable self-similar weighted spaces. In this paper, the uniqueness class of the Oseen vortex is expanded to include any solution that converges to the initial data in a sufficiently strong sense. This gives further evidence in support of the expectation that the Oseen vortex is the only possible mild solution that is identifiable as a vortex filament. The proof is a 3d variation of a 2d compactness/rigidity argument in $t \searrow 0$ originally due to Gallagher and Gallay.


title: Surface Drinfeld Torsors I : Higher Genus Associators

abstract: We develop a higher genus version of Drinfeld associators by means of operad theory. We start by introducing a framed version of rational associators and Grothendieck-Teichm\"uller groups and show that their definition is independent of the framing data. Next, we define a framed version of the universal KZ connection and we use it to show that over the complex numbers, the rational framed Drinfeld torsor is not empty. Next, we concentrate on the higher genus version of this story. We define an operad module of framed parenthesized higher genus braidings in prounipotent groupoids and we define its chord diagram counterpart. We then use these operadic modules to operadicly define higher genus associators and Grothendieck-Teichm\"uller groups, which again do not depend on the framing data. Finally, we compare our results in the genus $1$ case with those appearing in the litterature.


title: $^{12}$C(e,e'pN) Measurements of Short Range Correlations in the Tensor-to-Scalar Interaction Transition Region

abstract: High-momentum configurations of nucleon pairs at short-distance are probed using measurements of the $^{12}$C$(e,e'p)$ and $^{12}$C$(e,e'pN)$ reactions (where $N$ is either $n$ or $p$), at high-$Q^2$ and $x_B>1.1$. The data span a missing-momentum range of 300--1000 MeV/c and are predominantly sensitive to the transition region of the strong nuclear interaction from a Tensor to Scalar interaction. The data are well reproduced by theoretical calculations using the Generalized Contact Formalism with both chiral and phenomenological nucleon-nucleon ($NN$) interaction models. This agreement suggests that the measured high missing-momentum protons up to $1000$ MeV/c predominantly belong to short-ranged correlated (SRC) pairs. The measured $^{12}$C$(e,e'pN)$ / $^{12}$C$(e,e'p)$ and $^{12}$C$(e,e'pp)$ / $^{12}$C$(e,e'pn)$ cross-section ratios are consistent with a decrease in the fraction of proton-neutron SRC pairs and increase in the fraction of proton-proton SRC pairs with increasing missing momentum. This confirms the transition from an isospin-dependent tensor $NN$ interaction at $\sim 400$ MeV/c to an isospin-independent scalar interaction at high-momentum around $\sim 800$ MeV/c as predicted by theoretical calculation.


title: Internal Kinematics of M10 and M71

abstract: We present a study of the internal kinematics of two globular clusters, M10 (NGC 6254) and M71 (NGC 6838), using individual radial velocity (RV) measurements obtained from observations using the Hydra multi-object spectrograph on the WIYN 3.5 m telescope. We measured 120 RVs for stars in M10, of which 107 were determined to be cluster members. In M71, we measured 82 RVs and determined 78 of those measurements belonged to cluster members. Using the cluster members, we determine a mean RV of $75.9 \pm 4.0$ (s.d.) km s$^{-1}$ and $-22.9 \pm 2.2$ (s.d.) km s$^{-1}$ for M10 and M71, respectively. We combined the Hydra RV measurements with literature samples and performed a line-of-sight rotational analysis on both clusters. Our analysis has not revealed a statistically significant rotation in either of these clusters with the exception of the inner region (10\arcsec - 117\arcsec) of M10 for which we find hints of a marginally significant rotation with amplitude V$_{rot} = 1.14 \pm 0.18$ km s$^{-1}$. For M10, we calculate a central velocity dispersion of $\sigma_{0} = 5.44 \pm 0.61$ km s$^{-1}$, which gives a ratio of the amplitude of rotation to the central velocity dispersion $V_{rot}/\sigma_{0} = 0.21 \pm 0.04$. We also explored the rotation of the multiple stellar populations identified in M10 and M71 and found rotation (or lack thereof) in each population consistent with each other and the cluster global rotation signatures.


title: A new approach to mod 2 decompositions of BSU(2) and BSO(3)

abstract: Dwyer, Miller and Wilkerson proved that at the prime 2, the classifying spaces of SU(2) and SO(3) can be obtained as a homotopy pushout of the classifying spaces of certain subgroups. In this paper we show explicitly how these decompositions arise from the fusion systems of SU(2) and SO(3) over maximal discrete 2-toral subgroups.


title: The arrow of time across five centuries of classical music

abstract: The concept of time series irreversibility -- the degree by which the statistics of signals are not invariant under time reversal -- naturally appears in non-equilibrium physics in stationary systems which operate away from equilibrium and produce entropy. This concept has not been explored to date in the realm of musical scores as these are typically short sequences whose time reversibility estimation could suffer from strong finite size effects which preclude interpretability. Here we show that the so-called Horizontal Visibility Graph method -- which recently was shown to quantify such statistical property even in non-stationary signals -- is a method that can estimate time reversibility of short symbolic sequences, thus unlocking the possibility of exploring such properties in the context of musical compositions. Accordingly, we analyse over 8000 musical pieces ranging from the Renaissance to the early Modern period and certify that, indeed, most of them display clear signatures of time irreversibility. Since by construction stochastic processes with a linear correlation structure (such as 1/f noise) are time reversible, we conclude that musical compositions have a considerably richer structure, that goes beyond the traditional properties retrieved by the power spectrum or similar approaches. We also show that musical compositions display strong signs of nonlinear correlations, that nonlinearity is correlated to irreversibility, and that these are also related to asymmetries in the abundance of musical intervals, which we associate to the narrative underpinning a musical composition. These findings provide tools for the study of musical periods and composers, as well as criteria related to music appreciation and cognition.


title: Hopf monoids of ordered simplicial complexes

abstract: We study ordered matroids and generalized permutohedra from a Hopf theoretic point of view. Our main object is a Hopf monoid in the vector species of extended generalized permutahedra equipped with an order of the coordinates; this monoid extends the Hopf monoid of generalized permutahedra studied by Aguiar and Ardila. Our formula for the antipode is cancellation-free and multiplicity-free, and is supported only on terms that are compatible with the local geometry of the polyhedron. Our result is part of a larger program to understand orderings on ground sets of simplicial complexes (for instance, on shifted and matroid independence complexes). In this vein, we show that shifted simplicial complexes and broken circuit complexes generate Hopf monoids that are expected to exhibit similar behavior.


title: Algebras with matchings and link Floer homology

abstract: We explain how to use bordered algebras to compute a version of link Floer homology. As a corollary, we also give a fast computation of the Thuston polytope for links in the three-sphere.


title: Stability of doubly-intractable distributions

abstract: Doubly-intractable distributions appear naturally as posterior distributions in Bayesian inference frameworks whenever the likelihood contains a normalizing function $Z$. Having two such functions $Z$ and $\widetilde Z$ we provide estimates of the total variation and Wasserstein distance of the resulting posterior probability measures. As a consequence this leads to local Lipschitz continuity w.r.t. $Z$. In the more general framework of a random function $\widetilde Z$ we derive bounds on the expected total variation and expected Wasserstein distance. The applicability of the estimates is illustrated within the setting of two representative Monte Carlo recovery scenarios.


title: Edge Computing For Smart Health: Context-aware Approaches, Opportunities, and Challenges

abstract: Improving efficiency of healthcare systems is a top national interest worldwide. However, the need of delivering scalable healthcare services to the patients while reducing costs is a challenging issue. Among the most promising approaches for enabling smart healthcare (s-health) are edge-computing capabilities and next-generation wireless networking technologies that can provide real-time and cost-effective patient remote monitoring. In this paper, we present our vision of exploiting multi-access edge computing (MEC) for s-health applications. We envision a MEC-based architecture and discuss the benefits that it can bring to realize in-network and context-aware processing so that the s-health requirements are met. We then present two main functionalities that can be implemented leveraging such an architecture to provide efficient data delivery, namely, multimodal data compression and edge-based feature extraction for event detection. The former allows efficient and low distortion compression, while the latter ensures high-reliability and fast response in case of emergency applications. Finally, we discuss the main challenges and opportunities that edge computing could provide and possible directions for future research.


title: RescueNet: Joint Building Segmentation and Damage Assessment from Satellite Imagery

abstract: Accurate and fine-grained information about the extent of damage to buildings is essential for directing Humanitarian Aid and Disaster Response (HADR) operations in the immediate aftermath of any natural calamity. In recent years, satellite and UAV (drone) imagery has been used for this purpose, sometimes aided by computer vision algorithms. Existing Computer Vision approaches for building damage assessment typically rely on a two stage approach, consisting of building detection using an object detection model, followed by damage assessment through classification of the detected building tiles. These multi-stage methods are not end-to-end trainable, and suffer from poor overall results. We propose RescueNet, a unified model that can simultaneously segment buildings and assess the damage levels to individual buildings and can be trained end-toend. In order to to model the composite nature of this problem, we propose a novel localization aware loss function, which consists of a Binary Cross Entropy loss for building segmentation, and a foreground only selective Categorical Cross-Entropy loss for damage classification, and show significant improvement over the widely used Cross-Entropy loss. RescueNet is tested on the large scale and diverse xBD dataset and achieves significantly better building segmentation and damage classification performance than previous methods and achieves generalization across varied geographical regions and disaster types.


title: Evaluation of Generalizability of Neural Program Analyzers under Semantic-Preserving Transformations

abstract: The abundance of publicly available source code repositories, in conjunction with the advances in neural networks, has enabled data-driven approaches to program analysis. These approaches, called neural program analyzers, use neural networks to extract patterns in the programs for tasks ranging from development productivity to program reasoning. Despite the growing popularity of neural program analyzers, the extent to which their results are generalizable is unknown. In this paper, we perform a large-scale evaluation of the generalizability of two popular neural program analyzers using seven semantically-equivalent transformations of programs. Our results caution that in many cases the neural program analyzers fail to generalize well, sometimes to programs with negligible textual differences. The results provide the initial stepping stones for quantifying robustness in neural program analyzers.


title: Associative submanifolds and gradient cycles

abstract: We discuss a model for associative submanifolds in $G_{2}$-manifolds with K3 fibrations, in the adiabatic limit. The model involves graphs in a 3-manifold whose edges are locally gradient flow lines. We show that this model produces analogues of known singularity formation phenomena for associative submanifolds.


title: Hard Synchrotron Spectra from Magnetically Dominated Plasma Turbulence

abstract: Synchrotron emission from astrophysical nonthermal sources usually assumes that the emitting particles are isotropic. By means of large-scale two- and three-dimensional particle-in-cell simulations, we demonstrate that the dissipation of magnetically-dominated ($\sigma_0\gg1$) turbulence in pair plasmas leads to strongly anisotropic particle distributions. At Lorentz factors $\sim \sigma_0 \gamma_{th0}$ (here, $\gamma_{th0}$ is the initial Lorentz factor), the particle velocity is preferentially aligned with the local magnetic field; instead, the highest energy particles are roughly isotropic. This energy-dependent anisotropy leads to a synchrotron spectral flux $\nu F_\nu\propto \nu^s$ that is much harder than for isotropic particles. Remarkably, for $\sigma_0\gg1$ we find that the solid-angle-averaged spectral slope in the slow cooling regime is $s\sim 0.5-0.7$ for a wide range of turbulence fluctuations, $0.25\lesssim \delta B_{\rm rms0}^2/B_0^2\lesssim 4$, despite significant variations in the power-law energy spectrum of nonthermal particles. This is because weaker turbulence levels imprint a stronger degree of anisotropy, thereby counteracting the effect of the steeper particle spectrum. The synchrotron spectral slope may be even harder, $s\gtrsim 0.7$, if the observer is in the plane perpendicular to the mean magnetic field. Our results are independent of domain size and dimensionality. Our findings may help explaining the origin of hard synchrotron spectra of astrophysical nonthermal sources, most notably the radio spectrum of Pulsar Wind Nebulae.


title: Modelling with Volna-OP2: Towards tsunami threat reduction

abstract: Accurate and efficient tsunami modelling is essential for providing tsunami forecasts and hazard assessments. Volna-OP2 is a finite volume solver of the nonlinear shallow water equations and its capabilities of producing both faster than real time ensembles and high resolution inundation studies are presented here. The code is massively parallelised and can utilise various high performance computing architectures. When an earthquake is detected there is always some uncertainty on the source parameters. Generating a faster than real time ensemble for maximum wave heights which captures this uncertainty would be of great benefit to tsunami warning centres. The 2003 Boumerdes earthquake (Algeria) acts as a test case for showing Volna-OP2's ability at rapidly forecasting regional maximum wave heights. Drawing on various earthquake sources proposed in the literature and scaling the magnitudes to mimic uncertainty on the source, 20 separate earthquake realisations are simulated for 4 hours real time in 97s on two Nvidia V100 GPUs. Further a reduced ensemble of the Lisbon 1755 tsunami with an emphasis on the effects to the Irish coastline is presented. Where again various earthquake sources have been drawn from the literature and simulated on a regional scale. Finally, a pilot study which builds upon the reduced ensemble results investigates the inundation of a Lisbon tsunami on key sections of the Irish coastline. The results of this pilot study highlight that the inundation is constrained to low-lying areas with maximum run-up heights of $\approx 3.4m$ being found.


title: An Evaluation of DNN Architectures for Page Segmentation of Historical Newspapers

abstract: One important and particularly challenging step in the optical character recognition (OCR) of historical documents with complex layouts, such as newspapers, is the separation of text from non-text content (e.g. page borders or illustrations). This step is commonly referred to as page segmentation. While various rule-based algorithms have been proposed, the applicability of Deep Neural Networks (DNNs) for this task recently has gained a lot of attention. In this paper, we perform a systematic evaluation of 11 different published DNN backbone architectures and 9 different tiling and scaling configurations for separating text, tables or table column lines. We also show the influence of the number of labels and the number of training pages on the segmentation quality, which we measure using the Matthews Correlation Coefficient. Our results show that (depending on the task) Inception-ResNet-v2 and EfficientNet backbones work best, vertical tiling is generally preferable to other tiling approaches, and training data that comprises 30 to 40 pages will be sufficient most of the time.


title: Scaling of high-energy elastic scattering and the observation of Odderon

abstract: We provide a statistically significant observation of the elusive Odderon exchange, based on novel and model-independent analysis of the scaling properties of the differential cross sections of elastic $pp$ and $p\bar p$ scattering in the TeV energy range. We report the statistical significance of the observed Odderon signal at the level of 6.26$\sigma$.


title: The Impact of Heterogeneity and Geometry on the Proof Complexity of Random Satisfiability

abstract: Satisfiability is considered the canonical NP-complete problem and is used as a starting point for hardness reductions in theory, while in practice heuristic SAT solving algorithms can solve large-scale industrial SAT instances very efficiently. This disparity between theory and practice is believed to be a result of inherent properties of industrial SAT instances that make them tractable. Two characteristic properties seem to be prevalent in the majority of real-world SAT instances, heterogeneous degree distribution and locality. To understand the impact of these two properties on SAT, we study the proof complexity of random k-SAT models that allow to control heterogeneity and locality. Our findings show that heterogeneity alone does not make SAT easy as heterogeneous random k-SAT instances have superpolynomial resolution size. This implies intractability of these instances for modern SAT-solvers. On the other hand, modeling locality with an underlying geometry leads to small unsatisfiable subformulas, which can be found within polynomial time. A key ingredient for the result on geometric random k-SAT can be found in the complexity of higher-order Voronoi diagrams. As an additional technical contribution, we show a linear upper bound on the number of non-empty Voronoi regions, that holds for points with random positions in a very general setting. In particular, it covers arbitrary p-norms, higher dimensions, and weights affecting the area of influence of each point multiplicatively. This is in stark contrast to quadratic lower bounds for the worst case.


title: Training with Quantization Noise for Extreme Model Compression

abstract: We tackle the problem of producing compact models, maximizing their accuracy for a given model size. A standard solution is to train networks with Quantization Aware Training, where the weights are quantized during training and the gradients approximated with the Straight-Through Estimator. In this paper, we extend this approach to work beyond int8 fixed-point quantization with extreme compression methods where the approximations introduced by STE are severe, such as Product Quantization. Our proposal is to only quantize a different random subset of weights during each forward, allowing for unbiased gradients to flow through the other weights. Controlling the amount of noise and its form allows for extreme compression rates while maintaining the performance of the original model. As a result we establish new state-of-the-art compromises between accuracy and model size both in natural language processing and image classification. For example, applying our method to state-of-the-art Transformer and ConvNet architectures, we can achieve 82.5% accuracy on MNLI by compressing RoBERTa to 14MB and 80.0 top-1 accuracy on ImageNet by compressing an EfficientNet-B3 to 3.3MB.


title: Generating infinite monoids of cellular automata

abstract: For a group $G$ and a set $A$, let $\text{End}(A^G)$ be the monoid of all cellular automata over $A^G$, and let $\text{Aut}(A^G)$ be its group of units. By establishing a characterisation of surjunctuve groups in terms of the monoid $\text{End}(A^G)$, we prove that the rank of $\text{End}(A^G)$ (i.e. the smallest cardinality of a generating set) is equal to the rank of $\text{Aut}(A^G)$ plus the relative rank of $\text{Aut}(A^G)$ in $\text{End}(A^G)$, and that the latter is infinite when $G$ has an infinite decreasing chain of normal subgroups of finite index, condition which is satisfied, for example, for any infinite residually finite group. Moreover, when $A=V$ is a vector space over a field $\mathbb{F}$, we study the monoid $\text{End}_{\mathbb{F}}(V^G)$ of all linear cellular automata over $V^G$ and its group of units $\text{Aut}_{\mathbb{F}}(V^G)$. We show that if $G$ is an indicable group and $V$ is finite-dimensional, then $\text{End}_{\mathbb{F}}(V^G)$ is not finitely generated; however, for any finitely generated indicable group $G$, the group $\text{Aut}_{\mathbb{F}}(\mathbb{F}^G)$ is finitely generated if and only if $\mathbb{F}$ is finite.


title: Regularity for $C^{1,\alpha}$ interface transmission problems

abstract: We study existence, uniqueness, and optimal regularity of solutions to transmission problems for harmonic functions with $C^{1,\alpha}$ interfaces. For this, we develop a novel geometric stability argument based on the mean value property.


title: The Maximum Distance Problem and Minimal Spanning Trees

abstract: Given a compact $E \subset \mathbb{R}^n$ and $s > 0$, the maximum distance problem seeks a compact and connected subset of $\mathbb{R}^n$ of smallest one dimensional Hausdorff measure whose $s$-neighborhood covers $E$. For $E\subset \mathbb{R}^2$, we prove that minimizing over minimum spanning trees that connect the centers of balls of radius $s$, which cover $E$, solves the maximum distance problem. The main difficulty in proving this result is overcome by the proof of Lemma 3.5, which states that one is able to cover the $s$-neighborhood of a Lipschitz curve $\Gamma$ in $\mathbb{R}^2$ with a finite number of balls of radius $s$, and connect their centers with another Lipschitz curve $\Gamma_\ast$, where $\mathcal{H}^1(\Gamma_\ast)$ is arbitrarily close to $\mathcal{H}^1(\Gamma)$. We also present an open source package for computational exploration of the maximum distance problem using minimum spanning trees, available at https://github.com/mtdaydream/MDP_MST.


title: Building a Multi-domain Neural Machine Translation Model using Knowledge Distillation

abstract: Lack of specialized data makes building a multi-domain neural machine translation tool challenging. Although emerging literature dealing with low resource languages starts to show promising results, most state-of-the-art models used millions of sentences. Today, the majority of multi-domain adaptation techniques are based on complex and sophisticated architectures that are not adapted for real-world applications. So far, no scalable method is performing better than the simple yet effective mixed-finetuning, i.e finetuning a generic model with a mix of all specialized data and generic data. In this paper, we propose a new training pipeline where knowledge distillation and multiple specialized teachers allow us to efficiently finetune a model without adding new costs at inference time. Our experiments demonstrated that our training pipeline allows improving the performance of multi-domain translation over finetuning in configurations with 2, 3, and 4 domains by up to 2 points in BLEU.


title: Experimental observation of magnetic dimers in diluted Yb:YAlO$_3$

abstract: We present a comprehensive experimental investigation of Yb magnetic dimers in Yb$_{0.04}$Y$_{0.96}$AlO$_3$, an Yb-doped Yttrium Aluminum Perovskite (YAP) YAlO$_3$ by means of specific heat, magnetization and high-resolution inelastic neutron scattering (INS) measurements. In our sample, the Yb ions are randomly distributed over the lattice and $\sim 7$\% of Yb ions form quantum dimers due to nearest-neighbor antiferromagnetic coupling along the $c$-axis. At zero field, the dimer formation manifests itself in an appearance of an inelastic peak at $\Delta \approx 0.2$~meV in the INS spectrum and a Schottky-like anomaly in the specific heat. The structure factor of the INS peak exhibits a cosine modulation along the $L$ direction, in agreement with the $c$-axis nearest-neighbor intra-dimer coupling. A careful fitting of the low-temperature specific heat shows that the excited state is a degenerate triplet, which indicates a surprisingly small anisotropy of the effective Yb-Yb exchange interaction despite the low crystal symmetry and anisotropic magnetic dipole contribution, in agreement with previous reports for the Yb parent compound, YbAlO$_3$ [arXiv:1904.11513, arXiv:1902.04112], and in contrast to Yb$_2$Pt$_2$Pb [arXiv:1606.01309, arXiv:1907.01067]. The obtained results are precisely reproduced by analytical calculations for the Yb dimers.


title: Modeling Hall viscosity in magnetic-skyrmion systems

abstract: Magnetic skyrmions are topologically stable objects that are made with a bunch of spins tightly arranged in a smooth fashion. Their topological nature provides unusual and complex transport properties, such as the skyrmion Hall effect. Extensive Hall data have further revealed asymmetry between skyrmion and antiskyrmion Hall angles, which cannot be accounted by known mechanisms. Here, we explain this asymmetry by utilizing another universal transport coefficient called `Hall viscosity,' extensively studied in quantum Hall systems. Hall viscosity is modeled in steady-state skyrmions motion by generalizing the Thiele equation with a transverse velocity component and is independent of the skyrmion charge. Our analyses, based on available asymmetric Hall angle data, reveal this transverse force amounts 3\% - 5.4\% of the force due to the skyrmion Hall effect. Further clarification of Hall viscosity will be essential for designing next generation storage devices properly, not to mention for our deeper understanding of fundamental properties of nature.

categories: cond-mat.mtrl-sci cond-mat.mes-hall cond-mat.str-el quant-ph


title: Astronomical masers and Dicke's superradiance

abstract: We consider the radiation properties and processes of a gas with a population inversion using the formalism based on the Maxwell-Bloch equations. We focus on the maser action and Dicke's superradiance to establish their relationship in the overall radiation process during the temporal evolution of the system as a function of position. We show that the maser action and superradiance are not competing phenomena but are rather complementary, and define two distinct limits for the intensity of radiation. Masers characterise the quasi-steady state limit, when the population inversion density and the polarisation amplitude vary on time-scales longer than those of non-coherent processes affecting their evolution (e.g., collisions), while superradiance defines the fast transient regime taking place when these conditions are reversed. We show how a transition from a maser regime to superradiance will take place whenever a critical threshold for the column density of the population inversion is reached, at which point a strong level of coherence is established in the system and a powerful burst of radiation can ensue during the transient regime. This critical level also determines the spatial region where a transition from the unsaturated to the saturated maser regimes will take place; superradiance can thus be seen as the intermediary between the two. We also quantify the gain in radiation intensity attained during the superradiance phase relative to the two maser regimes, and show how the strong coherence level during superradiance is well suited to explain observations that reveal intense and fast radiation flares in maser-hosting regions.


title: The Reduced Proper Motion selected halo: methods and description of the catalogue

abstract: The Gaia mission has provided the largest ever astrometric chart of the Milky Way. Using it to map the Galactic halo is helpful for disentangling its merger history. The identification of halo stars in Gaia DR2 with reliable distance estimates requires special methods because such stars are typically farther away and scarce. We apply the reduced proper motion (RPM) method to identify halo main sequence stars on the basis of Gaia photometry and proper motions. Using the colour-absolute-magnitude relation for this type of stars, we calculate photometric distances. Our selection results in a set of $\sim10^7$ tentative main sequence halo stars with typical distance uncertainties of $7\%$ and with median velocity errors of 20 km/s. The median distance of our sample is $\sim 4.4$ kpc, with the faintest stars located at $\sim 16$ kpc. The spatial distribution of the stars in our sample is centrally concentrated. Visual inspection of the mean velocities of stars on the sky reveals large-scale patterns as well as clear imprints of the GD-1 stream and tentative hints of the Jhelum and Leiptr streams. Incompleteness and selection effects limit our ability to interpret the patterns reliably as well as to identify new substructures. We define a pseudo-velocity space by setting to zero the line-of-sight velocities of our sample stars. In this space, we recover several known structures such as the footprint of Gaia-Enceladus (i.e. the Gaia-Sausage) as well as the Helmi streams and some other retrograde substructures (Sequoia, Thamnos). We show that the two-point velocity correlation function reveals significant clustering on scales smaller than 100 km/s, of similar amplitude as found for the 6D Gaia halo sample. This clustering indicates the presence of nearby streams that are predominantly phase-mixed.


title: Augmented Lagrangian preconditioners for the Oseen-Frank model of nematic and cholesteric liquid crystals

abstract: We propose a robust and efficient augmented Lagrangian-type preconditioner for solving linearizations of the Oseen-Frank model arising in cholesteric liquid crystals. By applying the augmented Lagrangian method, the Schur complement of the director block can be better approximated by the weighted mass matrix of the Lagrange multiplier, at the cost of making the augmented director block harder to solve. In order to solve the augmented director block, we develop a robust multigrid algorithm which includes an additive Schwarz relaxation that captures a pointwise version of the kernel of the semi-definite term. Furthermore, we prove that the augmented Lagrangian term improves the discrete enforcement of the unit-length constraint. Numerical experiments verify the efficiency of the algorithm and its robustness with respect to problem-related parameters (Frank constants and cholesteric pitch) and the mesh size.


title: A New Constrained Optimization Model for Solving the Nonsymmetric Stochastic Inverse Eigenvalue Problem

abstract: The stochastic inverse eigenvalue problem aims to reconstruct a stochastic matrix from its spectrum. While there exists a large literature on the existence of solutions for special settings, there are only few numerical solution methods available so far. Recently, Zhao et al. (2016) proposed a constrained optimization model on the manifold of so-called isospectral matrices and adapted a modified Polak-Ribi\`ere-Polyak conjugate gradient method to the geometry of this manifold. However, not every stochastic matrix is an isospectral one and the model from Zhao et al. is based on the assumption that for each stochastic matrix there exists a (possibly different) isospectral, stochastic matrix with the same spectrum. We are not aware of such a result in the literature, but will see that the claim is at least true for $3 \times 3$ matrices. In this paper, we suggest to extend the above model by considering matrices which differ from isospectral ones only by multiplication with a block diagonal matrix with $2 \times 2$ blocks from the special linear group $SL(2)$, where the number of blocks is given by the number of pairs of complex-conjugate eigenvalues. Every stochastic matrix can be written in such a form, which was not the case for the form of the isospectral matrices. We prove that our model has a minimizer and show how the Polak-Ribi\`ere-Polyak conjugate gradient method works on the corresponding more general manifold. We demonstrate by numerical examples that the new, more general method performs similarly as the one from Zhao et al.


title: Shear-induced gelation of self-yielding active networks

abstract: Molecular-motor generated active stresses drive the cytoskeleton away from equilibrium, endowing it with tunable mechanical properties that are essential for diverse functions such as cell division and motility[1-5]. Designing analogous biomimetic systems is a key prerequisite for creating active matter that can emulate cellular functions[6-7]. These long-term goals requires understanding of how motor-generated stresses tune the mechanics of filamentous networks[8-11]. In microtubule-based active matter, kinesin motors generate extensile motion that leads to persistent breaking and reforming of the network links[12]. We study how such microscopic dynamics modifies the network's mechanical properties, uncovering that the network viscosity first increases with the imposed shear rate before transitioning back to a low-viscosity state. The non-monotonic shear-dependent viscosity can be controlled by tuning the speed of molecular motors. A two-state phenomenological model that incorporates liquid- and solid-like elements quantitatively relates the non-monotonic shear-rate-dependent viscosity to locally-measured flows. These studies show that rheology of extensile networks are different from previously studied active gels[13], where contractility enhances mechanical stiffness. Moreover, the flow induced gelation is not captured by continuum models of hydrodynamically interacting swimmers[14-21]. Observation of activity-dependent viscoelasticity necessitates the development of models for self-yielding of soft active solids whose intrinsic active stresses fluidize or stiffen the network.


title: Tests for multivariate normality -- a critical review with emphasis on weighted $L^2$-statistics

abstract: This article gives a synopsis on new developments in affine invariant tests for multivariate normality in an i.i.d.-setting, with special emphasis on asymptotic properties of several classes of weighted $L^2$-statistics. Since weighted $L^2$-statistics typically have limit normal distributions under fixed alternatives to normality, they open ground for a neighborhood of model validation for normality. The paper also reviews several other invariant tests for this problem, notably the energy test, and it presents the results of a large-scale simulation study. All tests under study are implemented in the accompanying R-package mnt.


title: Reinforcement Learning in a Physics-Inspired Semi-Markov Environment

abstract: Reinforcement learning (RL) has been demonstrated to have great potential in many applications of scientific discovery and design. Recent work includes, for example, the design of new structures and compositions of molecules for therapeutic drugs. Much of the existing work related to the application of RL to scientific domains, however, assumes that the available state representation obeys the Markov property. For reasons associated with time, cost, sensor accuracy, and gaps in scientific knowledge, many scientific design and discovery problems do not satisfy the Markov property. Thus, something other than a Markov decision process (MDP) should be used to plan / find the optimal policy. In this paper, we present a physics-inspired semi-Markov RL environment, namely the phase change environment. In addition, we evaluate the performance of value-based RL algorithms for both MDPs and partially observable MDPs (POMDPs) on the proposed environment. Our results demonstrate deep recurrent Q-networks (DRQN) significantly outperform deep Q-networks (DQN), and that DRQNs benefit from training with hindsight experience replay. Implications for the use of semi-Markovian RL and POMDPs for scientific laboratories are also discussed.


title: Visual explanation of country specific differences in Covid-19 dynamics

abstract: This report provides a visual examination of Covid-19 case and death data. In particular, it shows that country specific differences can too a large extend be explained by two easily interpreted parameters. Namely, the delay between reported cases and deaths and the fraction of cases observed. Furthermore, this allows to lower bound the actual total number of people already infected.


title: Nontrivial upper bounds for the least common multiple of an arithmetic progression

abstract: In this paper, we establish some nontrivial and effective upper bounds for the least common multiple of consecutive terms of a finite arithmetic progression. Precisely, we prove that for any two coprime positive integers $a$ and $b$, with $b\geq 2$, we have \[\mathrm{lcm}\left(a,a+b,\dots,a+nb\right) \leq \left(c_1\cdot b\log b\right)^{n+\left\lfloor \frac{a}{b}\right\rfloor}~~~~(\forall n\geq b+1),\] where $c_1=41.30142$. If in addition $b$ is a prime number and $a<b$, then we prove that for any $n\geq b+1$, we have $\mathrm{lcm}\left(a,a+b,\dots,a+nb\right) \leq \left(c_2\cdot b^{\frac{b}{b-1}}\right)^n$, where $c_2=12.30641$. Finally, we apply those inequalities to estimate the arithmetic function $M$ defined by $M(n):=\frac{1}{\varphi(n)}\sum_{\substack{1\leq\ell\leq n \\ \ell \wedge n=1}}\frac{1}{\ell}$ ($\forall n \geq 1$), as well as some values of the generalized Chebyshev function $\theta(x;k,\ell)$.


title: Parasitic neutron beam monitoring: proof of concept on gamma monitoring of neutron chopper phases

abstract: Neutron beam monitors are an essential diagnostic component of neutron scattering facilities. They are used to measure neutron flux, calibrating experiments performed on the instruments, allowing measurement of facility performance, understanding of the effect on the neutrons of beam-line components (such as choppers), calibration of detectors and tracking of beam stability. Ideally beam monitors {should} not perturb the beam. Previous work shows commercial beam monitors attenuate the beam by a few percent in the worst case due to the 1-2 mm thick Aluminium entrance and exit windows and the material inside. Parasitic methods of neutron beam diagnostics, where there is no beam monitor directly in the beam, would be preferable. This paper presents the concept of a parasitic method of monitoring the beam which can be used for neutron chopper phasing. This is achieved by placing a gamma detector close to a rotating chopper and measures a signal proportional to the flux absorbed by the chopper. Neutrons interact with the Boron absorber on the chopper disc lead to gamma emission at 480 keV. Detection of these gamma rays is used to determine the chopper phasing and timing. Potentially information on the flux of the beamline can be extracted. Results from a proof of concept implementation show that diagnosis of neutron chopper phases is feasible.


title: Nonequilibrium Solvent Polarization Effects in Real-Time Electronic Dynamics of Solute Molecules Subject to Time-Dependent Electric Fields: A New Feature of the Polarizable Continuum Model

abstract: We develop an extension of the time-dependent equation-of-motion formulation of the polarizable continuum model (EOM-TDPCM) to introduce nonequilibrium cavity field effects in quantum mechanical calculations of solvated molecules subject to time-dependent electric fields. This method has been implemented in Octopus, a state-of-the-art code for real-space, real-time time-dependent density functional theory (RT-TDDFT) calculations. To show the potential of our methodology, we perform EOM-TDPCM/RT-TDDFT calculations of trans-azobenzene in water and in other model solvents with shorter relaxation times. Our results for the optical absorption spectrum of trans-azobenzene show (i) that cavity field effects have a clear impact in the overall spectral shape and (ii) that an accurate description of the solute shape (as the one provided within PCM) is key to correctly account for cavity field effects.


title: Effects of in-medium $k_T$ broadening on di-jet observables

abstract: We study the influence of in-medium transverse momentum broadening on observables of jet pairs produced in nuclear collisions. These dijets were obtained numerically via the Monte-Carlo simulations of nuclear collisions, where the partons within the colliding nuclei were described via unintegrated parton densities, and their subsequent propagation within the medium of a quark-gluon plasma (QGP). For jet-evolution in the QGP both, kicks transverse to the momentum of the incoming particle as well as medium induced radiation were considered. After favorable comparison of our results with experimental LHC-data on jet-quenching we make predictions for the decorrelation of dijets. In particular, we study deviations from a transverse momentum broadening that follows a Gaussian distribution.


title: An Adaptive Intelligence Algorithm for Undersampled Knee MRI Reconstruction

abstract: Adaptive intelligence aims at empowering machine learning techniques with the additional use of domain knowledge. In this work, we present the application of adaptive intelligence to accelerate MR acquisition. Starting from undersampled k-space data, an iterative learning-based reconstruction scheme inspired by compressed sensing theory is used to reconstruct the images. We adopt deep neural networks to refine and correct prior reconstruction assumptions given the training data. The network was trained and tested on a knee MRI dataset from the 2019 fastMRI challenge organized by Facebook AI Research and NYU Langone Health. All submissions to the challenge were initially ranked based on similarity with a known groundtruth, after which the top 4 submissions were evaluated radiologically. Our method was evaluated by the fastMRI organizers on an independent challenge dataset. It ranked #1, shared #1, and #3 on respectively the 8x accelerated multi-coil, the 4x multi-coil, and the 4x single-coil track. This demonstrates the superior performance and wide applicability of the method.


title: Magic DIAMOND: Multi-Fascicle Diffusion Compartment Imaging with Tensor Distribution Modeling and Tensor-Valued Diffusion Encoding

abstract: Diffusion tensor imaging provides increased sensitivity to microstructural tissue changes compared to conventional anatomical imaging but also presents limited specificity. To tackle this problem, the DIAMOND model subdivides the voxel content into diffusion compartments and draws from diffusion-weighted data to estimate compartmental non-central matrix-variate Gamma distribution of diffusion tensors, thereby resolving crossing fascicles while accounting for their respective heterogeneity. Alternatively, tensor-valued diffusion encoding defines new acquisition schemes tagging specific features of the intra-voxel diffusion tensor distribution directly from the outcome of the measurement. However, the impact of such schemes on estimating brain microstructural features has only been studied in a handful of parametric single-fascicle models. In this work, we derive a general Laplace transform for the non-central matrix-variate Gamma distribution, which enables the extension of DIAMOND to tensor-valued encoded data. We then evaluate this "Magic DIAMOND" model in silico and in vivo on various combinations of tensor-valued encoded data. Assessing uncertainty on parameter estimation via stratified bootstrap, we investigate both voxel-based and fixel-based metrics by carrying out multi-peak tractography. We show that our estimated metrics can be mapped along tracks robustly across regions of fiber crossing, which opens new perspectives for tractometry and microstructure mapping along specific white-matter tracts.


title: Drug-Drug Interaction Prediction with Wasserstein Adversarial Autoencoder-based Knowledge Graph Embeddings

abstract: Interaction between pharmacological agents can trigger unexpected adverse events. Capturing richer and more comprehensive information about drug-drug interactions (DDI) is one of the key tasks in public health and drug development. Recently, several knowledge graph embedding approaches have received increasing attention in the DDI domain due to their capability of projecting drugs and interactions into a low-dimensional feature space for predicting links and classifying triplets. However, existing methods only apply a uniformly random mode to construct negative samples. As a consequence, these samples are often too simplistic to train an effective model. In this paper, we propose a new knowledge graph embedding framework by introducing adversarial autoencoders (AAE) based on Wasserstein distances and Gumbel-Softmax relaxation for drug-drug interactions tasks. In our framework, the autoencoder is employed to generate high-quality negative samples and the hidden vector of the autoencoder is regarded as a plausible drug candidate. Afterwards, the discriminator learns the embeddings of drugs and interactions based on both positive and negative triplets. Meanwhile, in order to solve vanishing gradient problems on the discrete representation--an inherent flaw in traditional generative models--we utilize the Gumbel-Softmax relaxation and the Wasserstein distance to train the embedding model steadily. We empirically evaluate our method on two tasks, link prediction and DDI classification. The experimental results show that our framework can attain significant improvements and noticeably outperform competitive baselines.


title: Mathematical modeling of blood flow through arterial bifurcation

abstract: The blood vascular system consists of blood vessels such as arteries, arterioles, capillaries and veins that convey blood throughout the body. The pressure difference which exists between the ends of the vessels provides the living force required to the flow of blood. In this study, we have used a model that is appropriate for use with blood flow simulation applications. We study the effect of decomposition of plaques and or tumour cells on the velocity profile, pressure of the incompressible, Newtonian fluid. We consider two different cases; vessels with bifurcation and without bifurcations and we display all the results through using Commercial software COMSOL Multiphysics 5.2.


title: Microscopic Mechanisms of Glass-Like Lattice Thermal Transport in Cubic Cu$_{12}$Sb$_{4}$S$_{13}$ Tetrahedrites

abstract: Materials based on cubic tetrahedrites (Cu$_{12}$Sb$_{4}$S$_{13}$) are useful thermoelectrics with unusual thermal and electrical transport properties, such as very low and nearly temperature-independent lattice thermal conductivity ($\kappa_{L}$). We explain the microscopic origin of the glass-like $\kappa_{L}$ in Cu$_{12}$Sb$_{4}$S$_{13}$ by explicitly treating anharmonicity up to quartic terms for both phonon energies and phonon scattering rates. We show that the strongly unstable phonon modes associated with trigonally coordinated Cu atoms are anharmonically stabilized above approximately $100$ K and continue hardening with increasing temperature, in accord with experimental data. This temperature induced hardening effect reduces scattering of heat carrying acoustic modes by reducing the available phase space for three-phonon processes, thereby balancing the conventional $\propto T$ increase in scattering due to phonon population and yielding nearly temperature-independent $\kappa_{L}$. Furthermore, we find that very strong phonon broadening lead to a qualitative breakdown of the conventional phonon-gas model and modify the dominant heat transport mechanism from the particle-like phonon wave packet propagation to incoherent tunneling described by off-diagonal terms in the heat-flux operator, which are typically prevailing in glasses and disordered crystals. Our work paves the way to a deeper understanding of glass-like thermal conductivity in complex crystals with strong anharmonicity.


title: From quasi-incompressible to semi-compressible fluids

abstract: A new concept of semi-compressible fluids is introduced for slightly compressible visco-elastic fluids (typically rather liquids than gasses) where mass density variations are negligible in some sense, while being directly controlled by pressure which is very small in comparison with the elastic bulk modulus. The physically consistent fully Eulerian models with specific dispersion of pressure-wave speed are devised. This contrasts to the so-called quasi-incompressible fluids which are described not physically consistently and, in fact, only approximate ideally incompressible ones in the limit. After surveying and modifying models for the quasi-incompressible fluids, we eventually devise some fully convective models complying with energy conservation and capturing phenomena as pressure-wave propagation with wave-length (and possibly also pressure) dependent velocity.


title: Active-Region Tilt Angles from White-Light Images and Magnetograms: The Role of Magnetic Tongues

abstract: The presence of elongations in active region (AR) polarities, called magnetic tongues, are mostly visible during their emergence phase. AR tilts have been measured thoroughly using long-term white-light (WL) databases, sometimes combined with magnetic field information. Since the influence of magnetic tongues on WL tilt measurements has not been taken into account before, we aim to investigate their role in tilt-angle values and to compare them with those derived from LOS magnetograms. We apply four methods to compute the tilt angle of generally bipolar ARs: one applies the k-means algorithm to WL data, a second one includes the magnetic field sign of the polarities to WL data, and a third one uses the magnetic flux-weighted center of each polarity. The tilt values computed in any of these ways are affected by the presence of magnetic tongues. Therefore, we apply the newly developed Core Field Fit Estimator (CoFFE) method to separate the magnetic flux in the tongues from that in the AR core. We compare the four computed tilt-angle values, as well as these with the ones reported in long-term WL databases. For ARs with low magnetic flux tongues the different methods report consistent tilt-angle values. But for ARs with high flux tongues there are noticeable discrepancies between all methods indicating that magnetic tongues affect differently WL and magnetic data. However, in general, CoFFE achieves a better estimation of the main bipole tilt because it removes both the effect of tongues as well as the emergence of secondary bipoles when it occurs in between the main bipole magnetic polarities.


title: Online Multiserver Convex Chasing and Optimization

abstract: We introduce the problem of $k$-chasing of convex functions, a simultaneous generalization of both the famous k-server problem in $R^d$, and of the problem of chasing convex bodies and functions. Aside from fundamental interest in this general form, it has natural applications to online $k$-clustering problems with objectives such as $k$-median or $k$-means. We show that this problem exhibits a rich landscape of behavior. In general, if both $k > 1$ and $d > 1$ there does not exist any online algorithm with bounded competitiveness. By contrast, we exhibit a class of nicely behaved functions (which include in particular the above-mentioned clustering problems), for which we show that competitive online algorithms exist, and moreover with dimension-free competitive ratio. We also introduce a parallel question of top-$k$ action regret minimization in the realm of online convex optimization. There, too, a much rougher landscape emerges for $k > 1$. While it is possible to achieve vanishing regret, unlike the top-one action case the rate of vanishing does not speed up for strongly convex functions. Moreover, vanishing regret necessitates both intractable computations and randomness. Finally we leave open whether almost dimension-free regret is achievable for $k > 1$ and general convex losses. As evidence that it might be possible, we prove dimension-free regret for linear losses via an information-theoretic argument.


title: HybridQA: A Dataset of Multi-Hop Question Answering over Tabular and Textual Data

abstract: Existing question answering datasets focus on dealing with homogeneous information, based either only on text or KB/Table information alone. However, as human knowledge is distributed over heterogeneous forms, using homogeneous information alone might lead to severe coverage problems. To fill in the gap, we present HybridQA https://github.com/wenhuchen/HybridQA, a new large-scale question-answering dataset that requires reasoning on heterogeneous information. Each question is aligned with a Wikipedia table and multiple free-form corpora linked with the entities in the table. The questions are designed to aggregate both tabular information and text information, i.e., lack of either form would render the question unanswerable. We test with three different models: 1) a table-only model. 2) text-only model. 3) a hybrid model that combines heterogeneous information to find the answer. The experimental results show that the EM scores obtained by two baselines are below 20\%, while the hybrid model can achieve an EM over 40\%. This gap suggests the necessity to aggregate heterogeneous information in HybridQA. However, the hybrid model's score is still far behind human performance. Hence, HybridQA can serve as a challenging benchmark to study question answering with heterogeneous information.


title: Learning 1-Dimensional Submanifolds for Subsequent Inference on Random Dot Product Graphs

abstract: A random dot product graph (RDPG) is a generative model for networks in which vertices correspond to positions in a latent Euclidean space and edge probabilities are determined by the dot products of the latent positions. We consider RDPGs for which the latent positions are randomly sampled from an unknown $1$-dimensional submanifold of the latent space. In principle, restricted inference, i.e., procedures that exploit the structure of the submanifold, should be more effective than unrestricted inference; however, it is not clear how to conduct restricted inference when the submanifold is unknown. We submit that techniques for manifold learning can be used to learn the unknown submanifold well enough to realize benefit from restricted inference. To illustrate, we test $1$- and $2$-sample hypotheses about the Fr\'{e}chet means of small communities of vertices, using the complete set of vertices to infer latent structure. We propose test statistics that deploy the Isomap procedure for manifold learning, using shortest path distances on neighborhood graphs constructed from estimated latent positions to estimate arc lengths on the unknown $1$-dimensional submanifold. Unlike conventional applications of Isomap, the estimated latent positions do not lie on the submanifold of interest. We extend existing convergence results for Isomap to this setting and use them to demonstrate that, as the number of auxiliary vertices increases, the power of our test converges to the power of the corresponding test when the submanifold is known. Finally, we apply our methods to an inference problem that arises in studying the connectome of the Drosophila larval mushroom body. The univariate learnt manifold test rejects ($p<0.05$), while the multivariate ambient space test does not ($p\gg0.05$), illustrating the value of identifying and exploiting low-dimensional structure for subsequent inference.


title: Regularization of backward time-fractional parabolic equations by Sobolev-type equations

abstract: The problem of determining the initial condition from noisy final observations in time-fractional parabolic equations is considered. This problem is well-known to be ill-posed and it is regularized by backward Sobolev-type equations. Error estimates of Holder type are obtained with a priori and a posteriori regularization parameter choice rules. The proposed regularization method results in a stable noniterative numerical scheme. The theoretical error estimates are confirmed by numerical tests for one- and two-dimensional equations


title: General Mapping between Complex Spatial and Temporal Frequencies by Analytical Continuation

abstract: This paper introduces a general technique for inter-mapping the complex spatial frequency (or propagation constant) $\gamma=\alpha+j\beta$ and the temporal frequency $\omega = \omega_\text{r}+j\omega_\text{i}$ of an arbitrary electromagnetic structure. This technique, based on the analytic property of complex functions describing physical phenomena, invokes the analytic continuity theorem to assert the unicity of the mapping function, and find this function from known data within a restricted domain of its analycity by curve-fitting to a generic polynomial expansion. It is not only applicable to canonical problems admitting an analytical solution, but to \emph{any} problems, from eigen-mode or driven-mode full-wave simulation results. The proposed technique is demonstrated for several systems, namely an unbounded lossy medium, a dielectric-filled rectangular waveguide, a periodically-loaded transmission line, a one-dimensional photonic crystal and a series-fed patch (SFP) leaky-wave antenna (LWA), and it is validated either by analytical results or by full-wave simulated results.


title: Communication Efficient Federated Learning with Energy Awareness over Wireless Networks

abstract: In federated learning (FL), reducing the communication overhead is one of the most critical challenges since the parameter server and the mobile devices share the training parameters over wireless links. With such consideration, we adopt the idea of SignSGD in which only the signs of the gradients are exchanged. Moreover, most of the existing works assume Channel State Information (CSI) available at both the mobile devices and the parameter server, and thus the mobile devices can adopt fixed transmission rates dictated by the channel capacity. In this work, only the parameter server side CSI is assumed, and channel capacity with outage is considered. In this case, an essential problem for the mobile devices is to select appropriate local processing and communication parameters (including the transmission rates) to achieve a desired balance between the overall learning performance and their energy consumption. Two optimization problems are formulated and solved, which optimize the learning performance given the energy consumption requirement, and vice versa. Furthermore, considering that the data may be distributed across the mobile devices in a highly uneven fashion in FL, a stochastic sign-based algorithm is proposed. Extensive simulations are performed to demonstrate the effectiveness of the proposed methods.


title: Ownership at Large -- Open Problems and Challenges in Ownership Management

abstract: Software-intensive organizations rely on large numbers of software assets of different types, e.g., source-code files, tables in the data warehouse, and software configurations. Who is the most suitable owner of a given asset changes over time, e.g., due to reorganization and individual function changes. New forms of automation can help suggest more suitable owners for any given asset at a given point in time. By such efforts on ownership health, accountability of ownership is increased. The problem of finding the most suitable owners for an asset is essentially a program comprehension problem: how do we automatically determine who would be best placed to understand, maintain, evolve (and thereby assume ownership of) a given asset. This paper introduces the Facebook Ownesty system, which uses a combination of ultra large scale data mining and machine learning and has been deployed at Facebook as part of the company's ownership management approach. Ownesty processes many millions of software assets (e.g., source-code files) and it takes into account workflow and organizational aspects. The paper sets out open problems and challenges on ownership for the research community with advances expected from the fields of software engineering, programming languages, and machine learning.


title: Nucleus I: Adjunction spectra in recommender systems and descent

abstract: Recommender systems build user profiles using concept analysis of usage matrices. The concepts are mined as spectra and form Galois connections. Descent is a general method for spectral decomposition in algebraic geometry and topology which also leads to generalized Galois connections. Both recommender systems and descent theory are vast research areas, separated by a technical gap so large that trying to establish a link would seem foolish. Yet a formal link emerged, all on its own, bottom-up, against authors' intentions and better judgment. Familiar problems of data analysis led to a novel solution in category theory. The present paper arose from a series of earlier efforts to provide a top-down account of these developments.


title: Efficient Algorithms for Battleship

abstract: We consider an algorithmic problem inspired by the Battleship game. In the variant of the problem that we investigate, there is a unique ship of shape $S \subset Z^2$ which has been translated in the lattice $Z^2$. We assume that a player has already hit the ship with a first shot and the goal is to sink the ship using as few shots as possible, that is, by minimizing the number of missed shots. While the player knows the shape $S$, which position of $S$ has been hit is not known. Given a shape $S$ of $n$ lattice points, the minimum number of misses that can be achieved in the worst case by any algorithm is called the Battleship complexity of the shape $S$ and denoted $c(S)$. We prove three bounds on $c(S)$, each considering a different class of shapes. First, we have $c(S) \leq n-1$ for arbitrary shapes and the bound is tight for parallelogram-free shapes. Second, we provide an algorithm that shows that $c(S) = O(\log n)$ if $S$ is an HV-convex polyomino. Third, we provide an algorithm that shows that $c(S) = O(\log \log n)$ if $S$ is a digital convex set. This last result is obtained through a novel discrete version of the Blaschke-Lebesgue inequality relating the area and the width of any convex body.


title: Time-domain Wavefield Reconstruction Inversion in a TTI medium

abstract: We introduce a generalization of time-domain wavefield reconstruction inversion to anisotropic acoustic modeling. Wavefield reconstruction inversion has been extensively researched in recent years for its ability to mitigate cycle skipping. The original method was formulated in the frequency domain with acoustic isotropic physics. However, frequency-domain modeling requires sophisticated iterative solvers that are difficult to scale to industrial-size problems and more realistic physical assumptions, such as tilted transverse isotropy, object of this study. The work presented here is based on a recently proposed dual formulation of wavefield reconstruction inversion, which allows time-domain propagator that are suitable to both large scales and more accurate physics.


title: A practical Response Adaptive Block Randomization (RABR) design with analytic type I error protection

abstract: Response adaptive randomization (RAR) is appealing from methodological, ethical, and pragmatic perspectives in the sense that subjects are more likely to be randomized to better performing treatment groups based on accumulating data. However, applications of RAR in confirmatory drug clinical trials with multiple active arms are limited largely due to its complexity, and lack of control of randomization ratios to different treatment groups. To address the aforementioned issues, we propose a Response Adaptive Block Randomization (RABR) design allowing arbitrarily pre-specified randomization ratios for the control and high-performing groups to meet clinical trial objectives. We show the validity of the conventional unweighted test in RABR with a controlled type I error rate based on the weighted combination test for sample size adaptive design invoking no large sample approximation. The advantages of the proposed RABR in terms of robustly reaching target final sample size to meet regulatory requirements and increasing statistical power as compared with the popular Doubly Adaptive Biased Coin Design (DBCD) are demonstrated by statistical simulations and a practical clinical trial design example.


title: Bounded Orbits of Quadratic Collatz-type Recursions

abstract: We characterize all bounded orbits of two similar Collatz-type quadratic mappings of the set of non-negative integers. In one case, where cycles of all possible lengths may occur, an orbit is bounded if and only if it reaches a cycle. For the other map we prove that every bounded orbit must reach 0 (in particular, there are no cycles).


title: Interband tunneling effects on materials transport properties using the first principles Wigner distribution

abstract: Electronic transport in narrow gap semiconductors is characterized by spontaneous vertical transitions between carriers in the valence and conduction bands, a phenomenon also known as Zener tunneling. However, this effect is not captured by existing models based on the Boltzmann transport equation. In this work, we propose a new fully first principles model for electronic transport using the Wigner distribution function and implement it to solve the equations of motion for electrons. The formalism generalizes the Boltzmann equation to materials with strong interband coupling and include transport contributions from off-diagonal components of the charge current operator. We illustrate the method with a study of Bi$_2$Se$_3$, showing that interband tunneling dominates the electron transport dynamics at experimentally relevant small doping concentrations, a behavior that is likely shared with other semiconductors, including topological insulators. Surprisingly, Zener tunneling occurs also between band subvalleys separated by energy much larger than the band gap.


title: Usable, Acceptable, Appropriable: Towards Practicable Privacy

abstract: A majority of the work on digital privacy and security has focused on users from developed countries who account for only around 20\% of the global population. Moreover, the privacy needs for population that is already marginalized and vulnerable differ from users who have privilege to access a greater social support system. We reflect on our experiences of introducing computers and the Internet to a group of sex-trafficking survivors in Nepal and highlight a few socio-political factors that have influenced the design space around digital privacy. These factors include the population's limited digital and text literacy skills and the fear of stigma against trafficked persons widely prevalent in Nepali society. We underscore the need to widen our perspective by focusing on practicable privacy, that is, privacy practices that are (1) usable, (2) acceptable, and (3) appropriable.


title: Reconfiguring and ramping-up ventilator production in the face of COVID-19: Can robots help?

abstract: As the COVID-19 pandemic expands, the shortening of medical equipment is swelling. A key piece of equipment getting far-out attention has been ventilators. The difference between supply and demand is substantial to be handled with normal production techniques, especially under social distancing measures in place. The study explores the rationale of human-robot teams to ramp up production using advantages of both the ease of integration and maintaining social distancing. The paper presents a model for faster integration of collaborative robots and design guidelines for workstation. The scenarios are evaluated for an open source ventilator through continuous human-robot simulation and amplification of results in a discrete event simulation.


title: Renormalization of group field theories for quantum gravity: new scaling results and some suggestions

abstract: We discuss motivation and goals of renormalization analyses of group field theory models of simplicial 4d quantum gravity, and review briefly the status of this research area. We present some new computations of perturbative GFT (spin foam) amplitudes, concerning in particular the scaling behaviour of radiative corrections to N-point functions. Finally, we single out key open issues and suggest a number of research directions for further progress in this area.


title: Hodge decompositions and Poincare duality models

abstract: We extend a CDGA $V$ with a perfect pairing of degree $n$ on cohomology to a CDGA $\hat V$ with a pairing of degree $n$ on chain level such that $\hat V$ admits a Hodge decomposition and retracts onto $V$ preserving the pairing on cohomology; here we suppose that $V$ is either 1-connected, or that $V$ is connected, of finite type, and $n$ is odd. We show that a Hodge decomposition of $\hat V$ induces a differential Poincar\'e duality model of $V$ in a natural way. Assuming that $H(V)$ is 1-connected, we apply our extension to a Sullivan model of $V$ in the proof of the existence and "uniqueness" of a 1-connected differential Poincar\'e duality model of $V$ by Lambrechts & Stanley; we eliminate their extra assumptions in the uniqueness statement, including $H^2(V)=0$ if $n$ is odd.


title: A detailed and direct proof of Skorohod-Wichura's theorem

abstract: The representation Skorohod theorem of weak convergence of random variables on a metric space goes back to Skorohod (1956) in the case where the metric space is the class of real-valued functions defined on [0,1] which are right-continuous and have left-hand limits when endowed with the Skorohod metric. Among the extensions of that to metric spaces, the version by Wichura (1970) seems to be the most fundamental. But the proof of Wichura seems to be destined to a very restricted public. We propose a more detailed proof to make it more accessible at the graduate level. However we do far more by simplifying it since important steps in the original proof are dropped, which leads to a direct proof that we hope to be more understandable to a larger spectrum of readers. The current version is more appropriate for different kinds of generalizations.


title: Microlocally accurate solution of the inverse source problem of thermoacoustic tomography

abstract: We consider the inverse source problem of thermo- and photoacoustic tomography, with data registered on an open surface partially surrounding the source of acoustic waves. Under the assumption of constant speed of sound we develop an explicit non-iterative reconstruction procedure that recovers the Radon transform of the sought source, up to an infinitely smooth additive error term. The source then can be found by inverting the Radon transform. Our analysis is microlocal in nature and does not provide a norm estimate on the error in the so obtained image. However, numerical simulations show that this error is quite small in practical terms. We also present an asymptotically fast implementation of this procedure for the case when the data are given on a circular arc in 2D.


title: Featureless adaptive optimization accelerates functional electronic materials design

abstract: Electronic materials exhibiting phase transitions between metastable states (e.g., metal-insulator transition materials with abrupt electrical resistivity transformations) are challenging to decode. For these materials, conventional machine learning methods display limited predictive capability due to data scarcity and the absence of features impeding model training. In this article, we demonstrate a discovery strategy based on multi-objective Bayesian optimization to directly circumvent these bottlenecks by utilizing latent variable Gaussian processes combined with high-fidelity electronic structure calculations for validation in the chalcogenide lacunar spinel family. We directly and simultaneously learn phase stability and band gap tunability from chemical composition alone to efficiently discover all superior compositions on the design Pareto front. Previously unidentified electronic transitions also emerge from our featureless adaptive optimization engine. Our methodology readily generalizes to optimization of multiple properties, enabling co-design of complex multifunctional materials, especially where prior data is sparse.


title: Spatially Resolved Chandra Spectroscopy of the Large Magellanic Cloud Supernova Remnant N132D

abstract: We perform detailed spectroscopy of the X-ray brightest supernova remnant (SNR) in the Large Magellanic Cloud (LMC), N132D, using Chandra archival observations. By analyzing the spectra of the entire well-defined rim, we determine the mean abundances for O, Ne, Mg, Si, S and Fe for the local LMC environment. We find evidence of enhanced O on the north-western and S on the north-eastern blast wave. By analyzing spectra interior to the remnant, we confirm the presence of a Si-rich relatively hot plasma (> 1.5 kev) that is also responsible for the Fe K emission. Chandra images show that the Fe K emission is distributed throughout the interior of the southern half of the remnant but does not extend out to the blast wave. We estimate the progenitor mass to be $15\pm5\,M_{\odot}$ using abundance ratios in different regions that collectively cover a large fraction of the remnant, as well as from the radius of the forward shock compared with models of an explosion in a cavity created by stellar winds. We fit ionizing and recombining plasma models to the Fe K emission and find that the current data cannot distinguish between the two, hence the origin of the high-temperature plasma remains uncertain. Our analysis is consistent with N132D being the result of a core-collapse supernova in a cavity created by its intermediate mass progenitor.


title: Ahlfors regularity and fractal dimension of Smale spaces

abstract: We prove that, up to topological conjugacy, every Smale space admits an Ahlfors regular Bowen measure. Bowen's construction of Markov partitions implies that Smale spaces are factors of topological Markov chains. The latter are equipped with Parry's measure which is Ahlfors regular. By extending Bowen's construction we create a tool for transferring the Ahlfors regularity of the Parry measure down to the Bowen measure of the Smale space. An essential part of our method uses a refined notion of approximation graphs over compact metric spaces. Moreover, we obtain new estimates for the Hausdorff, box-counting and Assouad dimensions of a large class of Smale spaces.


title: A Study on the Challenges of Using Robotics Simulators for Testing

abstract: Robotics simulation plays an important role in the design, development, and verification and validation of robotic systems. Recent studies have shown that simulation may be used as a cheaper, safer, and more reliable alternative to manual, and widely used, process of field testing. This is particularly important in the context of continuous integration pipelines, where integrated automated testing is key to reducing costs while maintaining system safety. However, simulation and automated testing are not seeing the degree of widespread adoption in practice that their potential would motivate. Our goal in this paper is to develop a principled understanding of the ways developers use simulation in their process, and the challenges they face in doing so. This type of understanding can guide the development of more effective simulators and testing techniques for modern robotics development. To that end, we conduct a survey of 82 robotics developers from a diversity of backgrounds that addresses the current capabilities and limits of simulation technology in practice. We find that simulation is used by 85% of our participants for testing, and that many participants desire to use simulation as part of their test automation. We identify 10 high-level challenges that impede developers from using simulation for manual and automated testing, and general purposes. These challenges include the gap between simulation and reality, a lack of reproducibility, and considerable resource costs associated with using simulators. Finally, we outline avenues for improvement in the development of new simulators that can help simulation reach its potential as a means of verification and validation.


title: Raman Spectrum and Polarizability of Liquid Water from Deep Neural Networks

abstract: We introduce a scheme based on machine learning and deep neural networks to model the environmental dependence of the electronic polarizability in insulating materials. Application to liquid water shows that training the network with a relatively small number of molecular configurations is sufficient to predict the polarizability of arbitrary liquid configurations in close agreement with ab initio density functional theory calculations. In combination with a neural network representation of the interatomic potential energy surface,the scheme allows us to calculate the Raman spectra along 2-nanosecond classical trajectories at different temperatures for H2O and D2O. The vast gains in efficiency provided by the machine learning approach enable longer trajectories and larger system sizes relative to ab initio methods, reducing the statistical error and improving the resolution of the low-frequency Raman spectra. Decomposing the spectra into intramolecular and intermolecular contributions elucidates the mechanisms behind the temperature dependence of the low-frequency and stretch modes.


title: F0-consistent many-to-many non-parallel voice conversion via conditional autoencoder

abstract: Non-parallel many-to-many voice conversion remains an interesting but challenging speech processing task. Many style-transfer-inspired methods such as generative adversarial networks (GANs) and variational autoencoders (VAEs) have been proposed. Recently, AutoVC, a conditional autoencoders (CAEs) based method achieved state-of-the-art results by disentangling the speaker identity and speech content using information-constraining bottlenecks, and it achieves zero-shot conversion by swapping in a different speaker's identity embedding to synthesize a new voice. However, we found that while speaker identity is disentangled from speech content, a significant amount of prosodic information, such as source F0, leaks through the bottleneck, causing target F0 to fluctuate unnaturally. Furthermore, AutoVC has no control of the converted F0 and thus unsuitable for many applications. In the paper, we modified and improved autoencoder-based voice conversion to disentangle content, F0, and speaker identity at the same time. Therefore, we can control the F0 contour, generate speech with F0 consistent with the target speaker, and significantly improve quality and similarity. We support our improvement through quantitative and qualitative analysis.


title: Joint Traffic-Aware UAV Placement and Predictive Routing for Aerial Networks

abstract: Aerial networks, composed of Unmanned Aerial Vehicles (UAVs) acting as Wi-Fi access points or cellular base stations, are emerging as an interesting solution to provide on-demand wireless connectivity to users, when there is no network infrastructure available, or to enhance the network capacity. This article proposes a traffic-aware topology control solution for aerial networks that holistically combines the placement of UAVs with a predictive and centralized routing protocol. The synergy created by the combination of the UAV placement and routing solutions allows the aerial network to seamlessly update its topology according to the users' traffic demand, whilst minimizing the disruption caused by the movement of the UAVs. As a result, the Quality of Service (QoS) provided to the users is improved. The components of the proposed solution are described and evaluated individually in this article by means of simulation and an experimental testbed. The results show that all the components improve the QoS provided to the users when compared to the corresponding baseline solutions.


title: Quantum Phases and Correlations Drive the Dynamics of Macroscopic Quantum Tunneling Escape in Quantum Simulators

abstract: Quantum tunneling remains unexplored in many regimes of many-body quantum physics, including the effect of quantum phase transitions on tunneling dynamics. In general, the quantum phase is a statement about the ground state and has no relation to far-from-equilibrium dynamics. Although tunneling is a highly dynamical process involving many excited states, we find that the quantum phase of the Bose-Hubbard model determines phase-dependent tunneling outcomes for the quantum tunneling escape, or quasi-bound problem. Superfluid and Mott insulator correlations lead to a new quantum tunneling rate, the quantum fluctuation rate. This rate shows surprising and highly dynamical features, such as oscillatory interference between trapped and escaped atoms and a completely different macroscopic quantum tunneling behavior for superfluid and Mott insulator phases. In the superfluid phase we find that escape dynamics are wave-like and coherent, leading to interference patterns in the density with a rapid decay process which is non-exponential. Quantum entropy production peaks when about half the atoms have escaped. In the Mott phase, despite stronger repulsive interactions, tunneling is significantly slowed by the presence of a Mott gap, creating an effective extra barrier to overcome. Only one atom can tunnel at a time, yet the decay process is nearly linear, completely defying the single-particle exponential model. Moreover, quantum entropy peaks when only about one quarter of the atoms have escaped. These and many other such effects go beyond the usual notions of single-particle quantum tunneling, quantum statistical effects on tunneling, and well-known semi-classical approaches from WKB to instanton theory. These results thus open up a new regime of exploration of far-from-equilibrium dynamics for quantum simulators and quantum dynamics.


title: Improv practices in Mathematics Active Teaching

abstract: This article explores the parallels between improvisational theater (Improv) and teaching in an Active Learning environment. It presents the notions of Active Teaching as a natural complement to Active Learning, and discusses how unexpected situations give rise to valuable Teaching Moments. These Teaching Moments can be strategically utilized following the rules of Improv. This article presents some examples of this in the Mathematics classroom, as well as the implementation of an Improv Seminar in the Mathematics Department at the University of California, Santa Cruz.


title: Homological mirror symmetry for Milnor fibers of simple singularities

abstract: We prove homological mirror symmetry for Milnor fibers of simple singularities in dimensions greater than one, which are among the log Fano cases of Conjecture 1.5 in arXiv:1806.04345. The proof is based on a relation between matrix factorizations and Calabi--Yau completions. As an application, we give an explicit computation of the Hochschild cohomology group of the derived $n$-preprojective algebra of a Dynkin quiver for any $n \geq 1$, and the symplectic cohomology group of the Milnor fiber of any simple singularity in any dimension greater than one.


title: A Practical Introduction to Bayesian Estimation of Causal Effects: Parametric and Nonparametric Approaches

abstract: Substantial advances in Bayesian methods for causal inference have been developed in recent years. We provide an introduction to Bayesian inference for causal effects for practicing statisticians who have some familiarity with Bayesian models and would like an overview of what it can add to causal estimation in practical settings. In the paper, we demonstrate how priors can induce shrinkage and sparsity on parametric models and be used to perform probabilistic sensitivity analyses around causal assumptions. We provide an overview of nonparametric Bayesian estimation and survey their applications in the causal inference literature. Inference in the point-treatment and time-varying treatment settings are considered. For the latter, we explore both static and dynamic treatment regimes. Throughout, we illustrate implementation using off-the-shelf open source software. We hope the reader will walk away with implementation-level knowledge of Bayesian causal inference using both parametric and nonparametric models. All synthetic examples and code used in the paper are publicly available on a companion GitHub repository.