title: Optical brain imaging using a semi-transparent organic light-emitting diode

abstract: We report optical brain imaging using a semi-transparent organic light-emitting diode (OLED) based on the orange light-emitting polymer (LEP) Livilux PDO-124. The OLED serves as a compact, extended light source which is capable of uniformly illuminating the cortical surface when placed across a burr hole in the skull. Since all layers of the OLED are substantially transparent to photons with energies below the optical gap of the LEP, light emitted or reflected by the cortical surface may be efficiently transmitted through the OLED and into the objective lens of a low magnification microscope ('macroscope'). The OLED may be placed close to the cortical surface, providing efficient coupling of incident light into the brain cavity; furthermore, the macroscope may be placed close to the upper surface of the OLED, enabling efficient collection of reflected/emitted light from the cortical surface. Hence the use of a semi-transparent OLED simplifies the optical setup, while at the same time maintaining high sensitivity. The OLED is applied here to one of the most demanding forms of optical brain imaging, namely extrinsic optical imaging involving a voltage sensitive dye (VSD). Specifically, we carry out functional imaging of the primary visual cortex (V1) of a rat, using the voltage sensitive dye RH-1691 as a reporter. Imaging through the OLED light-source, we are able to resolve small (~ 0.1 %) changes in the fluorescence intensity of the dye due to changes in the neuronal membrane potential following a visual stimulus. Results are obtained on a single trial basis -- i.e. without averaging over multiple measurements -- with a time-resolution of ten milliseconds.


title: Rotational state-changing collisions of C$_2$H$^-$ and C$_2$N$^-$ anions with He under interstellar and cold ion trap conditions: a computational comparison

abstract: We present an extensive range of quantum calculations for the state-changing rotational dynamics involving two simple molecular anions which are expected to play some role in evolutionary analysis of chemical networks in the Interstellar environments, C$_2$H$^-$($X^1\Sigma^+$) and C$_2$N$^-$ ($X^3 \Sigma^-$) but for which inelastic rates are only known for C$_2$H$^-$. The same systems are also of direct interest in modelling selective photo-detachment (PD) experiments in cold ion traps where the He atoms function as the chief buffer gas at the low trap temperatures. This study employs accurate, \textit{ab initio} calculations of the interaction potential energy surfaces (PESs) for these anions, treated as Rigid Rotors (RR) and the He atom to obtain a wide range of state-changing quantum cross sections and rates at temperatures up to about 100 K. The results are analysed and compared for the two systems, to show differences and similarities between their rates of state-changing dynamics.


title: Affordance as general value function: A computational model

abstract: General value functions (GVFs) in the reinforcement learning (RL) literature are long-term predictive summaries of the outcomes of agents following specific policies in the environment. Affordances as perceived action possibilities with specific valence may be cast into predicted policy-relative goodness and modelled as GVFs. A systematic explication of this connection shows that GVFs and especially their deep learning embodiments (1) realize affordance prediction as a form of direct perception, (2) illuminate the fundamental connection between action and perception in affordance, and (3) offer a scalable way to learn affordances using RL methods. Through an extensive review of existing literature on GVF applications and representative affordance research in robotics, we demonstrate that GVFs provide the right framework for learning affordances in real-world applications. In addition, we highlight a few new avenues of research opened up by the perspective of "affordance as GVF", including using GVFs for orchestrating complex behaviors.


title: Post Training Uncertainty Calibration of Deep Networks For Medical Image Segmentation

abstract: Neural networks for automated image segmentation are typically trained to achieve maximum accuracy, while less attention has been given to the calibration of their confidence scores. However, well-calibrated confidence scores provide valuable information towards the user. We investigate several post hoc calibration methods that are straightforward to implement, some of which are novel. They are compared to Monte Carlo (MC) dropout. They are applied to neural networks trained with cross-entropy (CE) and soft Dice (SD) losses on BraTS 2018 and ISLES 2018. Surprisingly, models trained on SD loss are not necessarily less calibrated than those trained on CE loss. In all cases, at least one post hoc method improves the calibration. There is limited consistency across the results, so we can't conclude on one method being superior. In all cases, post hoc calibration is competitive with MC dropout. Although average calibration improves compared to the base model, subject-level variance of the calibration remains similar.


title: Fast Local Attack: Generating Local Adversarial Examples for Object Detectors

abstract: The deep neural network is vulnerable to adversarial examples. Adding imperceptible adversarial perturbations to images is enough to make them fail. Most existing research focuses on attacking image classifiers or anchor-based object detectors, but they generate globally perturbation on the whole image, which is unnecessary. In our work, we leverage higher-level semantic information to generate high aggressive local perturbations for anchor-free object detectors. As a result, it is less computationally intensive and achieves a higher black-box attack as well as transferring attack performance. The adversarial examples generated by our method are not only capable of attacking anchor-free object detectors, but also able to be transferred to attack anchor-based object detector.


title: Counterfactual Ghost Imaging

abstract: We give a protocol for ghost imaging in a way that is always counterfactual - while imaging an object, no light interacts with that object. This extends the idea of counterfactuality beyond communication, showing how this interesting phenomenon can be leveraged for metrology. Given, in the infinite limit, no photons ever go to the imaged object, it presents a method of imaging even the most light-sensitive of objects without damaging them. Even when not in the infinite limit, it still provides a many-fold improvement in visibility and signal-to-noise ratio over previous protocols, with over an order of magnitude reduction in absorbed intensity.


title: Search for squarks and gluinos in final states with jets and missing transverse momentum using 139 fb$^{-1}$ of $\sqrt{s}$ =13 TeV $pp$ collision data with the ATLAS detector

abstract: A search for the supersymmetric partners of quarks and gluons (squarks and gluinos) in final states containing jets and missing transverse momentum, but no electrons or muons, is presented. The data used in this search were recorded by the ATLAS experiment in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}$ = 13 TeV during Run 2 of the Large Hadron Collider, corresponding to an integrated luminosity of 139 fb$^{-1}$. The results are interpreted in the context of various $R$-parity-conserving models where squarks and gluinos are produced in pairs or in association and a neutralino is the lightest supersymmetric particle. An exclusion limit at the 95% confidence level on the mass of the gluino is set at 2.30 TeV for a simplified model containing only a gluino and the lightest neutralino, assuming the latter is massless. For a simplified model involving the strong production of mass-degenerate first- and second-generation squarks, squark masses below 1.85 TeV are excluded if the lightest neutralino is massless. These limits extend substantially beyond the region of supersymmetric parameter space excluded previously by similar searches with the ATLAS detector.


title: Robust Odometry and Mapping for Multi-LiDAR Systems with Online Extrinsic Calibration

abstract: Combining multiple LiDARs enables a robot to maximize its perceptual awareness of environments and obtain sufficient measurements, which is promising for simultaneous localization and mapping (SLAM). This paper proposes a system to achieve robust and simultaneous extrinsic calibration, odometry, and mapping for multiple LiDARs. Our approach starts with measurement preprocessing to extract edge and planar features from raw measurements. After a motion and extrinsic initialization procedure, a sliding window-based multi-LiDAR odometry runs onboard to estimate poses with online calibration refinement and convergence identification. We further develop a mapping algorithm to construct a global map and optimize poses with sufficient features together with a method to model and reduce data uncertainty. We validate our approach's performance with extensive experiments on ten sequences (4.60km total length) for the calibration and SLAM and compare them against the state-of-the-art. We demonstrate that the proposed work is a complete, robust, and extensible system for various multi-LiDAR setups. The source code, datasets, and demonstrations are available at https://ram-lab.com/file/site/m-loam.


title: The flaring X-ray corona in the quasar PDS 456

abstract: New Swift monitoring observations of the variable, radio-quiet quasar, PDS 456, are presented. A bright X-ray flare was captured in September 2018, the flux increasing by a factor of 4 and with a doubling time-scale of 2 days. From the light crossing argument, the coronal size is inferred to be about 30 gravitational radii for a black hole mass of $10^{9} {\rm M}_{\odot}$ and the total flare energy exceeds $10^{51}$ erg. A hardening of the X-ray emission accompanied the flare, with the photon index decreasing from $\Gamma=2.2$ to $\Gamma=1.7$ and back again. The flare is produced in the X-ray corona, the lack of any optical or UV variability being consistent with a constant accretion rate. Simultaneous XMM-Newton and NuSTAR observations were performed, $1-3$ days after the flare peak and during the decline phase. These caught PDS 456 in a bright, bare state, where no disc wind absorption features are apparent. The hard X-ray spectrum shows a high energy roll-over, with an e-folding energy of $E_{\rm fold}=51^{+11}_{-8}$ keV. The deduced coronal temperature, of $kT=13$ keV, is one of the coolest measured in any AGN and PDS 456 lies well below the predicted pair annihilation line in X-ray corona. The spectral variability, becoming softer when fainter following the flare, is consistent with models of cooling X-ray coronae. Alternatively, an increase in a non-thermal component could contribute towards the hard X-ray flare spectrum.


title: Fit to Measure: Reasoning about Sizes for Robust Object Recognition

abstract: Service robots can help with many of our daily tasks, especially in those cases where it is inconvenient or unsafe for us to intervene: e.g., under extreme weather conditions or when social distance needs to be maintained. However, before we can successfully delegate complex tasks to robots, we need to enhance their ability to make sense of dynamic, real world environments. In this context, the first prerequisite to improving the Visual Intelligence of a robot is building robust and reliable object recognition systems. While object recognition solutions are traditionally based on Machine Learning methods, augmenting them with knowledge based reasoners has been shown to improve their performance. In particular, based on our prior work on identifying the epistemic requirements of Visual Intelligence, we hypothesise that knowledge of the typical size of objects could significantly improve the accuracy of an object recognition system. To verify this hypothesis, in this paper we present an approach to integrating knowledge about object sizes in a ML based architecture. Our experiments in a real world robotic scenario show that this combined approach ensures a significant performance increase over state of the art Machine Learning methods.


title: Field theory generalizations of two-body Calogero-Moser models in the form of Landau-Lifshitz equations

abstract: We give detailed description for continuous version of the classical IRF-Vertex relation, where on the IRF side we deal with the Calogero-Moser-Sutherland models. Our study is based on constructing modifications of the Higgs bundles of infinite rank over elliptic curve and its degenerations. In this way the previously predicted gauge equivalence between L-A pairs of the Landau-Lifshitz type equations and 1+1 field theory generalization of the Calogero-Moser-Sutherland models is described. In this paper the ${\rm sl}_2$ case is studied. Explicit changes of variables are obtained between the rational, trigonometric and elliptic models.


title: A Statistical Framework for Low-bitwidth Training of Deep Neural Networks

abstract: Fully quantized training (FQT), which uses low-bitwidth hardware by quantizing the activations, weights, and gradients of a neural network model, is a promising approach to accelerate the training of deep neural networks. One major challenge with FQT is the lack of theoretical understanding, in particular of how gradient quantization impacts convergence properties. In this paper, we address this problem by presenting a statistical framework for analyzing FQT algorithms. We view the quantized gradient of FQT as a stochastic estimator of its full precision counterpart, a procedure known as quantization-aware training (QAT). We show that the FQT gradient is an unbiased estimator of the QAT gradient, and we discuss the impact of gradient quantization on its variance. Inspired by these theoretical results, we develop two novel gradient quantizers, and we show that these have smaller variance than the existing per-tensor quantizer. For training ResNet-50 on ImageNet, our 5-bit block Householder quantizer achieves only 0.5% validation accuracy loss relative to QAT, comparable to the existing INT8 baseline.


title: Simulations and the Lamplighter group

abstract: We introduce a notion of "simulation" for labelled graphs, in which edges of the simulated graph are realized by regular expressions in the simulating graph, and prove that the tiling problem (aka "domino problem") for the simulating graph is at least as difficult as that for the simulated graph. We apply this to the Cayley graph of the "lamplighter group" $L=\mathbb Z/2\wr\mathbb Z$, and more generally to "Diestel-Leader graphs". We prove that these graphs simulate the plane, and thus deduce that the seeded tiling problem is unsolvable on the group $L$. We note that $L$ does not contain any plane in its Cayley graph, so our undecidability criterion by simulation covers cases not covered by Jeandel's criterion based on translation-like action of a product of finitely generated infinite groups. Our approach to tiling problems is strongly based on categorical constructions in graph theory.


title: Ice Monitoring in Swiss Lakes from Optical Satellites and Webcams using Machine Learning

abstract: Continuous observation of climate indicators, such as trends in lake freezing, is important to understand the dynamics of the local and global climate system. Consequently, lake ice has been included among the Essential Climate Variables (ECVs) of the Global Climate Observing System (GCOS), and there is a need to set up operational monitoring capabilities. Multi-temporal satellite images and publicly available webcam streams are among the viable data sources to monitor lake ice. In this work we investigate machine learning-based image analysis as a tool to determine the spatio-temporal extent of ice on Swiss Alpine lakes as well as the ice-on and ice-off dates, from both multispectral optical satellite images (VIIRS and MODIS) and RGB webcam images. We model lake ice monitoring as a pixel-wise semantic segmentation problem, i.e., each pixel on the lake surface is classified to obtain a spatially explicit map of ice cover. We show experimentally that the proposed system produces consistently good results when tested on data from multiple winters and lakes. Our satellite-based method obtains mean Intersection-over-Union (mIoU) scores >93%, for both sensors. It also generalises well across lakes and winters with mIoU scores >78% and >80% respectively. On average, our webcam approach achieves mIoU values of 87% (approx.) and generalisation scores of 71% (approx.) and 69% (approx.) across different cameras and winters respectively. Additionally, we put forward a new benchmark dataset of webcam images (Photi-LakeIce) which includes data from two winters and three cameras.


title: SIRI: Spatial Relation Induced Network For Spatial Description Resolution

abstract: Spatial Description Resolution, as a language-guided localization task, is proposed for target location in a panoramic street view, given corresponding language descriptions. Explicitly characterizing an object-level relationship while distilling spatial relationships are currently absent but crucial to this task. Mimicking humans, who sequentially traverse spatial relationship words and objects with a first-person view to locate their target, we propose a novel spatial relationship induced (SIRI) network. Specifically, visual features are firstly correlated at an implicit object-level in a projected latent space; then they are distilled by each spatial relationship word, resulting in each differently activated feature representing each spatial relationship. Further, we introduce global position priors to fix the absence of positional information, which may result in global positional reasoning ambiguities. Both the linguistic and visual features are concatenated to finalize the target localization. Experimental results on the Touchdown show that our method is around 24\% better than the state-of-the-art method in terms of accuracy, measured by an 80-pixel radius. Our method also generalizes well on our proposed extended dataset collected using the same settings as Touchdown.


title: From frieze patterns to cluster categories

abstract: Motivated by Conway and Coxeter's combinatorial results concerning frieze patterns, we sketch an introduction to the theory of cluster algebras and cluster categories for acyclic quivers. The goal is to show how these more abstract theories provide a conceptual explanation for phenomena concerning friezes, principally integrality and periodicity.


title: A local-realistic model for quantum theory

abstract: We provide a rigorous definition of local realism. We show that the universal wave function cannot be a complete description of a local reality. Finally, we construct a local-realistic model for quantum theory.


title: Optimization of Robot Grasping Forces and Worst Case Loading

abstract: We consider the optimization of the vector of grasping forces that support a known generalized force acting on the grasped object---a rigid body or a mechanism. Working in the framework of finite-dimensional normed vector spaces and their dual spaces, the cost function to be minimized is assumed to be a norm on the space of grasping forces. We present an expression for the optimum which depends on the external force and the kinematics of the grasping system. Next, assuming that optimal grasping forces are applied using force control, and assuming that there is a bound on the norm of the admissible grasping forces, we characterize the largest norm of an external force that the grasping system may support, that is, the norm of the worst-case loading that may be applied and still be supported. A few simple examples are given for the sake of illustration.


title: No evidence of phosphine in the atmosphere of Venus by independent analyses

abstract: The detection of phosphine (PH3) in the atmosphere of Venus has been recently reported based on millimeter-wave radio observations (Greaves et al. 2020), and its re-analyses (Greaves et al. 2021a/b). In this Matters Arising we perform an independent reanalysis, identifying several issues in the interpretation of the spectroscopic data. As a result, we determine sensitive upper-limits for PH3 in Venus' atmosphere (>75 km, above the cloud decks) that are discrepant with the findings in G2020 and G2021a/b. The measurements target the fundamental first rotational transition of PH3 (J=1-0) at 266.944513 GHz, which was observed with the James Clerk Maxwell Telescope (JCMT) in June 2017 and with the Atacama Large Millimeter/submillimeter Array (ALMA) in March 2019. This line's center is near the SO2 (J=309,21-318,24) transition at 266.943329 GHz (only 1.3 km/s away from the PH3 line) which represents a potential source of contamination. The JCMT and ALMA data, as presented in G2020, are at spectral resolutions comparable to the frequency separation of the two lines. Moreover, the spectral features identified are several km/s in width, and therefore do not permit distinct spectroscopic separation of the candidate spectral lines of PH3 and SO2. We present the radiative transfer modelling we have performed and then discuss the ALMA and JCMT analyses in turn.


title: Adsorption-controlled growth of MnTe(Bi2Te3)n by molecular beam epitaxy exhibiting stoichiometry-controlled magnetism

abstract: We report the growth of the intrinsic magnetic topological system MnTe(Bi2Te3)n by molecular beam epitaxy. By mapping the temperature and the Bi:Mn flux ratio, it is shown that there is a narrow growth window for the n=1 phase MnBi2Te4 with 2.0<Bi:Mn<2.6 at 225 {\deg}C. Here the films are stoichiometric and excess Bi and Te is not incorporated. At higher flux ratios (Bi:Mn>4.5) it is found that the n = 2 MnBi4Te7 phase is stabilized. Transport measurements indicate that the MnBi2Te4 and MnBi4Te7 undergo magnetic transitions around 25 K, and 10 K, respectively, consistent with antiferromagnetic phases found in the bulk. Further, for Mn-rich conditions (Bi:Mn<2), ferromagnetism emerges that exhibits a clear hysteretic state in the Hall effect, which likely indicates Mn-doped MnBi2Te4. Understanding how to grow ternary chalcogenide phases is the key to synthesizing new materials and to interface magnetism and topology, which together are routes to realize and control exotic quantum phenomena.


title: What Color is this? Explaining Art Restoration Research Methods using Interactive Museum Installations

abstract: This case study describes an approach to designing interactive museum installations as a student project with the aim of presenting the research results of the restoration process of paintings to a wide range of visitors. During one and a half years, the Chair of Media Design created five interactive media stations in two lectures to enrich the special exhibition "Veronese: The Cuccina Cycle. The Restored Masterpiece". The project was realised in close communication with the conservators of the Dresden State Art Collections and the employees of the Science and Archaeometric Laboratory of the Dresden University of Fine Arts. The students had to learn about the foreign content and how to translate it into a media-related environment. With suitable teaching methods, we pushed the students towards a deeper understanding of the matter.


title: A constrained Cosserat shell model up to order $O(h^5)$: Modelling, existence of minimizers, relations to classical shell models and scaling invariance of the bending tensor

abstract: We consider a recently introduced geometrically nonlinear elastic Cosserat shell model incorporating effects up to order $O(h^5)$ in the shell thickness $h$. We develop the corresponding geometrically nonlinear constrained Cosserat shell model, we show the existence of minimizers for the $O(h^5)$ and $O(h^3)$ case and we draw some connections to existing models and classical shell strain measures. Notably, the role of the appearing new bending tensor is highlighted and investigated with respect to an invariance condition of Acharya [Int. J. Solids and Struct., 2000] which will be further strengthened.


title: Microcanonical Tunneling Rates from Density-of-States Instanton Theory

abstract: Semiclassical instanton theory is a form of quantum transition-state theory which can be applied to computing thermal reaction rates for complex molecular systems including quantum tunneling effects. There have been a number of attempts to extend the theory to treat microcanonical rates. However, the previous formulations are either computationally unfeasible for large systems due to an explicit sum over states or they involve extra approximations which make them less reliable. We propose a robust and practical microcanonical formulation called density-of-states instanton theory, which avoids the sum over states altogether. In line with the semiclassical approximations inherent to the instanton approach, we employ the stationary-phase approximation to the inverse Laplace transform to obtain the densities of states. This can be evaluated using only post-processing of the data available from a small set of instanton calculations, such that our approach remains computationally efficient. We show that the new formulation predicts results that agree well with quantum scattering theory for an atom-diatom reaction and with experiments for a photoexcited unimolecular hydrogen transfer in a Criegee intermediate. When the thermal rate is evaluated from a Boltzmann average over our new microcanonical formalism, it can overcome some problems of conventional instanton theory. In particular, it predicts a smooth transition at the crossover temperature and is able to describe bimolecular reactions with pre-reactive complexes such as CH3OH + OH.


title: A normalized solitary wave solution of the Maxwell-Dirac equations

abstract: We prove the existence of a $L^2$-normalized solitary wave solution for the Maxwell-Dirac equations in (3+1)-Minkowski space. In addition, for the Coulomb-Dirac model, describing fermions with attractive Coulomb interactions in the mean-field limit, we prove the existence of the (positive) energy minimizer.


title: The large-scale magnetic field of Proxima Centauri near activity maximum

abstract: We report the detection of a large-scale magnetic field at the surface of the slowly-rotating fully-convective M dwarf Proxima Centauri. Ten circular polarization spectra, collected from April to July 2017 with the HARPS-Pol spectropolarimeter, exhibit rotationally-modulated Zeeman signatures suggesting a stellar rotation period of $89.8 \pm 4.0$ d. Using Zeeman-Doppler Imaging, we invert the circular polarization spectra into a surface distribution of the large-scale magnetic field. We find that Proxima Cen hosts a large-scale magnetic field of typical strength 200 G, whose topology is mainly poloidal, and moderately axisymmetric, featuring, in particular, a dipole component of 135 G tilted at 51$^{\circ}$ to the rotation axis. The large-scale magnetic flux is roughly 3 times smaller than the flux measured from the Zeeman broadening of unpolarized lines, which suggests that the underlying dynamo is efficient at generating a magnetic field at the largest spatial scales. Our observations occur $\sim$1 yr after the maximum of the reported 7 yr-activity cycle of Proxima Cen, which opens the door for the first long-term study of how the large-scale field evolves with the magnetic cycle in a fully-convective very-low-mass star. Finally, we find that Proxima Cen's habitable zone planet, Proxima-b, is likely orbiting outside the Alfv\`en surface, where no direct magnetic star-planet interactions occur.


title: The twinning operation on graphs does not always preserve $e$-positivity

abstract: Motivated by Stanley's $\mathbf{(3+1)}$-free conjecture on chromatic symmetric functions, Foley, Ho\`{a}ng and Merkel introduced the concept of strong $e$-positivity and conjectured that a graph is strongly $e$-positive if and only if it is (claw, net)-free. In order to study strongly $e$-positive graphs, they further introduced the twinning operation on a graph $G$ with respect to a vertex $v$, which adds a vertex $v'$ to $G$ such that $v$ and $v'$ are adjacent and any other vertex is adjacent to both of them or neither of them. Foley, Ho\`{a}ng and Merkel conjectured that if $G$ is $e$-positive, then so is the resulting twin graph $G_v$ for any vertex $v$. Based on the theory of chromatic symmetric functions in non-commuting variables developed by Gebhard and Sagan, we establish the $e$-positivity of a class of graphs called tadpole graphs. By considering the twinning operation on a subclass of these graphs with respect to certain vertices we disprove the latter conjecture of Foley, Ho\`{a}ng and Merkel. We further show that if $G$ is $e$-positive, the twin graph $G_v$ and more generally the clan graphs $G^{(k)}_v$ ($k \ge 1$) may not even be $s$-positive, where $G^{(k)}_v$ is obtained from $G$ by applying $k$ twinning operations to $v$.


title: On homotopy exponentials in path categories

abstract: In the first part of this paper we show that path categories are enriched over groupoids, in a way that is compatible with a suitable 2-category of path categories. In the second part we introduce a new notion of homotopy exponential and homotopy Pi-type for path categories, the strong homotopy exponential and the strong homotopy Pi-type. We prove some of their basic properties, along with two important technical propositions. The results in the second part of this paper are mostly obtained using the groupoid enrichment of path categories established in the first part.


title: Faster Lagrangian-Based Methods in Convex Optimization

abstract: In this paper, we aim at unifying, simplifying and improving the convergence rate analysis of Lagrangian-based methods for convex optimization problems. We first introduce the notion of nice primal algorithmic map, which plays a central role in the unification and in the simplification of the analysis of most Lagrangian-based methods. Equipped with a nice primal algorithmic map, we then introduce a versatile generic scheme, which allows for the design and analysis of Faster LAGrangian (FLAG) methods with new provably sublinear rate of convergence expressed in terms of function values and feasibility violation of the original (non-ergodic) generated sequence. To demonstrate the power and versatility of our approach and results, we show that most well-known iconic Lagrangian-based schemes admit a nice primal algorithmic map, and hence share the new faster rate of convergence results within their corresponding FLAG.


title: A proof of the mod $4$ unimodal sequence conjectures and related mock theta functions

abstract: In 2012 Bryson, Ono, Pitman and Rhoades showed how the generating functions for certain strongly unimodal sequences are related to quantum modular and mock modular forms. They proved some parity results and conjectured some mod 4 congruences for the coefficients of these generating functions. In 2016 Kim, Lim and Lovejoy obtained similar results for odd-balanced unimodal sequences and made similar mod 4 conjectures. We prove all of these mod 4 conjectures and similar congruences for the Andrews spt-function and related mock theta functions. Our method of proof involves new Hecke-Rogers type identities for indefinite binary quadratic forms and the Hurwitz class number.


title: Computation of Large Asymptotics of 3-Manifold Quantum Invariants

abstract: Quantum topological invariants have played an important role in computational topology, and they are at the heart of major modern mathematical conjectures. In this article, we study the experimental problem of computing large $r$ values of Turaev-Viro invariants $\mathrm{TV}_r$. We base our approach on an optimized backtracking algorithm, consisting of enumerating combinatorial data on a triangulation of a 3-manifold. We design an easily computable parameter to estimate the complexity of the enumeration space, based on lattice point counting in polytopes, and show experimentally its accuracy. We apply this parameter to a preprocessing strategy on the triangulation, and combine it with multi-precision arithmetics in order to compute the Turaev-Viro invariants. We finally study the improvements brought by these optimizations compared to state-of-the-art implementations, and verify experimentally Chen and Yang's volume conjecture on a census of closed 3-manifolds.


title: Observation of Coulomb-assisted nuclear bound state of $\Xi^-$-$^{14}$N system

abstract: In an emulsion-counter hybrid experiment performed at J-PARC, a $\Xi^-$ absorption event was observed which decayed into twin single-$\Lambda$ hypernuclei. Kinematic calculations enabled a unique identification of the reaction process as $\Xi^{-} + ^{14}$N$\ \rightarrow\ ^{10}_\Lambda$Be + $^5_\Lambda$He. For the binding energy of the $\Xi^{-}$ hyperon in the $\Xi^-$-$^{14}$N system a value of $1.27 \pm 0.21$ MeV was deduced. The energy level of $\Xi^-$ is likely a nuclear $1p$ state which indicates a weak ${\Xi}N$-$\Lambda\Lambda$ coupling.


title: Multitask Training with Text Data for End-to-End Speech Recognition

abstract: We propose a multitask training method for attention-based end-to-end speech recognition models. We regularize the decoder in a listen, attend, and spell model by multitask training it on both audio-text and text-only data. Trained on the 100-hour subset of LibriSpeech, the proposed method, without requiring an additional language model, leads to an 11% relative performance improvement over the baseline and approaches the performance of language model shallow fusion on the test-clean evaluation set. We observe a similar trend on the whole 960-hour LibriSpeech training set. Analyses of different types of errors and sample output sentences demonstrate that the proposed method can incorporate language level information, suggesting its effectiveness in real-world applications.


title: Atmosphere Models of Brown Dwarfs Irradiated by White Dwarfs: Analogues for Hot and Ultra-Hot Jupiters

abstract: Irradiated brown dwarfs (BDs) provide natural laboratories to test our understanding of substellar and irradiated atmospheres. A handful of short-period BDs around white dwarfs (WDs) have been observed, but the uniquely intense UV-dominated irradiation presents a modeling challenge. Here, we present the first fully self-consistent 1D atmosphere models that take into account the UV irradiation's effect on the object's temperature structure. We explore two BD-WD systems, namely WD-0137-349 and EPIC-212235321. WD-0137-349B has an equilibrium temperature that would place it in the transition between hot and ultra-hot Jupiters, while EPIC-212235321B has an equilibrium temperature higher than all ultra-hot Jupiters except KELT-9b. We explore some peculiar aspects of irradiated BD atmospheres and show that existing photometry can be well-fit with our models. Additionally, the detections of atomic emission lines from these BDs can be explained by a strong irradiation-induced temperature inversion, similar to inversions recently explored in ultra-hot Jupiters. Our models of WD-0137-349B can reproduce the observed equivalent width of many but not all of these atomic lines. We use the observed photometry of these objects to retrieve the temperature structure using the PHOENIX ExoplaneT Retrieval Algorithm (PETRA) and demonstrate that the structures are consistent with our models, albeit somewhat cooler at low pressures. We then discuss the similarities and differences between this class of irradiated brown dwarf and the lower-mass ultra-hot Jupiters. Lastly, we describe the behavior of irradiated BDs in color-magnitude space to show the difficulty in classifying irradiated BDs using otherwise well-tested methods for isolated objects.


title: Dynamics of zeroes under repeated differentiation

abstract: Consider a random polynomial $P_n$ of degree $n$ whose roots are independent random variables sampled according to some probability distribution $\mu_0$ on the complex plane $\mathbb C$. It is natural to conjecture that, for a fixed $t\in [0,1)$ and as $n\to\infty$, the zeroes of the $[tn]$-th derivative of $P_n$ are distributed according to some measure $\mu_t$ on $\mathbb C$. Assuming either that $\mu_0$ is concentrated on the real line or that it is rotationally invariant, Steinerberger [Proc. AMS, 2019] and O'Rourke and Steinerberger [arXiv:1910.12161] derived nonlocal transport equations for the density of roots. We introduce a different method to treat such problems. In the rotationally invariant case, we obtain a closed formula for $\psi(x,t)$, the asymptotic density of the radial parts of the roots of the $[tn]$-th derivative of $P_n$. Although its derivation is non-rigorous, we provide numerical evidence for its correctness and prove that it solves the PDE of O'Rourke and Steinerberger. Moreover, we present several examples in which the solution is fully explicit (including the special case in which the initial condition $\psi(x,0)$ is an arbitrary convex combination of delta functions) and analyze some properties of the solutions such as the behavior of void annuli and circles of zeroes. As an additional support for the correctness of the method, we show that a similar method, applied to the case when $\mu_0$ is concentrated on the real line, gives a correct result which is known to have an interpretation in terms of free probability.


title: Summing a family of generalized Pell numbers

abstract: A new family of generalized Pell numbers was recently introduced and studied by Br\'od \cite{Dorota}. These number possess, as Fibonacci numbers, a Binet formula. Using this, partial sums of arbitrary powers of generalized Pell numbers can be summed explicitly. For this, as a first step, a power $P_n^l$ is expressed as a linear combination of $P_{mn}$. The summation of such expressions is then manageable using generating functions. Since the new family contains a parameter $R=2^r$, the relevant manipulations are quite involved, and computer algebra produced huge expressions that where not trivial to handle at times.


title: An efficient nonconvex reformulation of stagewise convex optimization problems

abstract: Convex optimization problems with staged structure appear in several contexts, including optimal control, verification of deep neural networks, and isotonic regression. Off-the-shelf solvers can solve these problems but may scale poorly. We develop a nonconvex reformulation designed to exploit this staged structure. Our reformulation has only simple bound constraints, enabling solution via projected gradient methods and their accelerated variants. The method automatically generates a sequence of primal and dual feasible solutions to the original convex problem, making optimality certification easy. We establish theoretical properties of the nonconvex formulation, showing that it is (almost) free of spurious local minima and has the same global optimum as the convex problem. We modify PGD to avoid spurious local minimizers so it always converges to the global minimizer. For neural network verification, our approach obtains small duality gaps in only a few gradient steps. Consequently, it can quickly solve large-scale verification problems faster than both off-the-shelf and specialized solvers.


title: Sub-sampling for Efficient Non-Parametric Bandit Exploration

abstract: In this paper we propose the first multi-armed bandit algorithm based on re-sampling that achieves asymptotically optimal regret simultaneously for different families of arms (namely Bernoulli, Gaussian and Poisson distributions). Unlike Thompson Sampling which requires to specify a different prior to be optimal in each case, our proposal RB-SDA does not need any distribution-dependent tuning. RB-SDA belongs to the family of Sub-sampling Duelling Algorithms (SDA) which combines the sub-sampling idea first used by the BESA [1] and SSMC [2] algorithms with different sub-sampling schemes. In particular, RB-SDA uses Random Block sampling. We perform an experimental study assessing the flexibility and robustness of this promising novel approach for exploration in bandit models.


title: Probabilistic learning on manifolds constrained by nonlinear partial differential equations for small datasets

abstract: A novel extension of the Probabilistic Learning on Manifolds (PLoM) is presented. It makes it possible to synthesize solutions to a wide range of nonlinear stochastic boundary value problems described by partial differential equations (PDEs) for which a stochastic computational model (SCM) is available and depends on a vector-valued random control parameter. The cost of a single numerical evaluation of this SCM is assumed to be such that only a limited number of points can be computed for constructing the training dataset (small data). Each point of the training dataset is made up realizations from a vector-valued stochastic process (the stochastic solution) and the associated random control parameter on which it depends. The presented PLoM constrained by PDE allows for generating a large number of learned realizations of the stochastic process and its corresponding random control parameter. These learned realizations are generated so as to minimize the vector-valued random residual of the PDE in the mean-square sense. Appropriate novel methods are developed to solve this challenging problem. Three applications are presented. The first one is a simple uncertain nonlinear dynamical system with a nonstationary stochastic excitation. The second one concerns the 2D nonlinear unsteady Navier-Stokes equations for incompressible flows in which the Reynolds number is the random control parameter. The last one deals with the nonlinear dynamics of a 3D elastic structure with uncertainties. The results obtained make it possible to validate the PLoM constrained by stochastic PDE but also provide further validation of the PLoM without constraint.


title: Distributed Optimization with Quantization for Computing Wasserstein Barycenters

abstract: We study the problem of the decentralized computation of entropy-regularized semi-discrete Wasserstein barycenters over a network. Building upon recent primal-dual approaches, we propose a sampling gradient quantization scheme that allows efficient communication and computation of approximate barycenters where the factor distributions are stored distributedly on arbitrary networks. The communication and algorithmic complexity of the proposed algorithm are shown, with explicit dependency on the size of the support, the number of distributions, and the desired accuracy. Numerical results validate our algorithmic analysis.


title: Charge Density Waves and Coplanar Magnetism in Gd2PdSi3

abstract: The intermetallic Gd2PdSi3 has recently generated a lot of excitement after reports that it hosts chiral magnetic nano-skyrmions despite its centrosymmetric crystal structure. Using magnetic-field-dependent polarized resonant elastic x-ray scattering (REXS), we find than an unexpected incommensurate charge density wave (CDW) appears below the ordering transition with a wavevector equal to that of the magnetic textures. Furthermore, we show these incommensurate magnetic textures in Gd2PdSi3 are highly anisotropic, with the vast majority of the spin modulation lying within the hexagonal ab-plane. This observation is not compatible with the previously suggested non-coplanar magnetic textures and coplanar alternatives are discussed. Our results thus refute the interpretation of the observed large anomalous Hall and Nernst effects in Gd2PdSi3 as arising from topologically-nontrivial magnetic skyrmions.


title: Two-component galaxy models with a central BH -- II. The ellipsoidal case

abstract: Recently, two-component spherical galaxy models have been presented, where the stellar profile is described by a Jaffe law, and the total density by another Jaffe law, or by an $r^{-3}$ law at large radii. We extend these two families to their ellipsoidal axisymmetric counterparts: the JJe and J3e models. The total and stellar density distributions can have different flattenings and scale lengths, and the dark matter halo is defined by difference. First, the analytical conditions required to have a nowhere negative dark matter halo density are derived. The Jeans equations for the stellar component are then solved analytically, in the limit of small flattenings, also in presence of a central BH. The azimuthal velocity dispersion anisotropy is described by the Satoh $k$-decomposition. Finally, we present the analytical formulae for velocity fields near the center and at large radii, together with the various terms entering the Virial Theorem. The JJe and J3e models can be useful in a number of theoretical applications, e.g. to explore the role of the various parameters (flattening, relative scale lengths, mass ratios, rotational support) in determining the behavior of the stellar kinematical fields before performing more time-expensive integrations with specific galaxy models, to test codes of stellar dynamics, and in numerical simulations of gas flows in galaxies.


title: Graphene on noncommutative plane and the Seiberg-Witten map

abstract: Graphene on two dimensional (2D) noncommutative (NC) plane in the presence of a constant background magnetic field has been studied. To handel the gauge-invariance issue we start our analysis by a effective massles NC Dirac field theory where we incorporate the Seiberg-Witten (SW) map along with the Moyal star ($\star$) product. The gauge-invariant Hamiltonian of a massless Dirac particle is then computed which is used to study the relativistic Landau problem of graphene on NC plane. Specifically we study the quantum dynamics of a massless relativistic electron moves on monolayer graphene, in the presence of a constant background magnetic field, on NC plane. We also compute the energy spectrum of the NC Landau system in graphene. The results obtained are corrected by the spatial NC parameter $\theta$. Finally we visit the Weyl equation for electron in graphene on NC plane. Interestingly, in this case helicity is found to be $\theta$ modified.


title: Dynamical Analysis of Screening in Scalar-Tensor Theory

abstract: We investigate linear and non-linear dynamics of spherically symmetric perturbations on a static configuration in scalar-tensor theories focusing on the chameleon screening mechanism. We particularly address two questions: how much the perturbations can source the fifth force when the static background is well screened, and whether the resultant fifth force can change the stability and structure of the background configuration. For linear perturbations, we derive a lower bound for the square of the Fourier mode frequency $\omega^2$ using the adiabatic approximation. There may be unstable modes if this lower bound is negative, and we find that the condition of the instability can be changed by the fifth force although this effect is suppressed by the screening parameter. For non-linear perturbations, because we are mainly interested in short wavelength modes for which the fifth force may become stronger, we perform numerical simulations under the planar approximation. For a sufficiently large initial amplitude of the density perturbation, we find that the magnitude of the fifth force can be comparable to that of Newtonian gravity even when the model parameters are chosen so that the static background is well screened. It is also shown that if the screening is effective for the static background, the fluid dynamics is mostly governed by the pressure gradient and is not significantly affected by the fifth force.


title: Effect of increased stability of peptide-based coatings in the Casimir regime via nanoparticle doping

abstract: We find that thin peptide films and coatings doped with metallic nanoparticles are more stable due to the role of electromagnetic fluctuations. It is shown that for the doped freestanding in vacuum peptide film the Casimir attraction becomes larger in magnitude. For dielectric substrates coated with peptide films, the nanoparticle doping leads to a wider range of film thicknesses where the Casimir pressure is attractive and to larger pressure magnitudes at the points of extremum. The doping of peptide coatings with magnetic nanoparticles preserves all the advantages of nonmagnetic ones and simultaneously imparts superparamagnetic properties to the coating which could extend significantly the application areas of bioelectronics.


title: Are Multi-language Design Smells Fault-prone? An Empirical Study

abstract: Nowadays, modern applications are developed using components written in different programming languages. These systems introduce several advantages. However, as the number of languages increases, so does the challenges related to the development and maintenance of these systems. In such situations, developers may introduce design smells (i.e., anti-patterns and code smells) which are symptoms of poor design and implementation choices. Design smells are defined as poor design and coding choices that can negatively impact the quality of a software program despite satisfying functional requirements. Studies on mono-language systems suggest that the presence of design smells affects code comprehension, thus making systems harder to maintain. However, these studies target only mono-language systems and do not consider the interaction between different programming languages. In this paper, we present an approach to detect multi-language design smells in the context of JNI systems. We then investigate the prevalence of those design smells. Specifically, we detect 15 design smells in 98 releases of nine open-source JNI projects. Our results show that the design smells are prevalent in the selected projects and persist throughout the releases of the systems. We observe that in the analyzed systems, 33.95% of the files involving communications between Java and C/C++ contains occurrences of multi-language design smells. Some kinds of smells are more prevalent than others, e.g., Unused Parameters, Too Much Scattering, Unused Method Declaration. Our results suggest that files with multi-language design smells can often be more associated with bugs than files without these smells, and that specific smells are more correlated to fault-proneness than others.


title: Generalized permutahedra and Schubert calculus

abstract: We connect generalized permutahedra with Schubert calculus. Thereby, we give sufficient vanishing criteria for Schubert intersection numbers of the flag variety. Our argument utilizes recent developments in the study of Schubitopes, which are Newton polytopes of Schubert polynomials. The resulting tableau test executes in polynomial time.


title: Real-time, Software-Defined, GPU-Based Receiver Field Trial

abstract: We demonstrate stable real-time operation of a software-defined, GPU-based receiver over a metropolitan network. Massive parallelization is exploited for implementing direct-detection and coherent Kramers-Kronig detection in real time at 2 and 1 GBaud, respectively.


title: Implications of the lowest frequency detection of the persistent counterpart of FRB121102

abstract: Context. The repeating FRB121102 is so far the only extra-galactic Fast Radio Burst found to be associated with a counterpart, a steady radio source with a nearly flat spectral energy distribution (SED) in centimeter wavelengths. Aims. Previous observations of the persistent source down to $1.6$~GHz has shown no sign of a spectral turn-over. Absorption is expected to eventually cause a turn-over at lower frequencies. Better constraints on the physical parameters of the emitting medium can be derived by detecting the self-absorption frequency. Methods. We used the Giant Metre-Wave Radio Telescope (GMRT) during the period of July to December 2017 to observe the source at low radio frequencies down to $400$~MHz. Results. The spectral energy distribution of the source remains optically thin even at $400$~MHz, with a spectral index of $\nu^{-(0.07 \pm 0.03)}$ similar to what is seen in Galactic plerions. Using a generic synchrotron radiation model, we obtain constraints on properties of the non-thermal plasma and the central engine powering it. Conclusions. We present low-frequency detections of the persistent source associated with FRB121102. Its characteristic flat SED extends down to $400$~MHz. Like Galactic plerions, the energy in the persistent source is carried predominantly by leptons. The emitting plasma has a $B< 0.01$~G, and its age is $> 524 \left(\frac{B}{0.01 {\rm G}} \right)^{-3/2}$. We show that the energetics of the persistent source requires an initial spin period shorter than 36~ms, and the magnetic field of the neutron star must exceed $4.5\times 10^{12}$~G. This implies that the persistent source does not necessarily require energetic input from a magnetar.


title: Blowing-up Nonautonomous Vector Fields: Infinite Delay Equations and Invariant Manifolds

abstract: We show the existence of nonautonomous invariant manifolds for planar, asymptotically autonomous differential equations, that have equilibrium solutions with zero Lyapunov spectrum. These invariant manifolds correspond to the stable and unstable manifold of a desingularized equation, that we obtain by using the blow-up method. More precisely, the blow-up method is extended to the nonautonomous setting and transforms the original finite-dimensional ordinary differential equation into an infinite-dimensional delay equation with infinite delay, but hyperbolic structure. In the technical construction of the invariant manifolds for the delay equation, we have to carefully study the effect of the time reparametrization used for desingularization in the blown-up space to guarantee sufficient regularity. This allows us to employ a Lyapunov-Perron argument to obtain existence of an invariant manifold. We combine the last step with an implicit function argument to study differentiability of the manifold. Finally, we reverse the blow-up transformation of space and time variables obtaining invariant manifolds of the initially considered equation.


title: Corrections of order O(E^2_e/m^2_N), caused by weak magnetism and proton recoil, to the neutron lifetime and correlation coefficients of the neutron beta decay

abstract: We calculate the contributions of weak magnetism and proton recoil of order O(E^2_e/m^2_N)~10^{-5}, i.e. to next-to-next-to-leading order in the large nucleon mass expansion, to the neutron lifetime and correlation coefficients of the neutron beta decay, where E_e and m_N are the electron energy and the nucleon mass, respectively. We analyze the electron-energy and angular distribution for the neutron beta decay with a polarized neutron, a polarized electron and an unpolarized proton. Together with Wilkinson's corrections (Nucl. Phys. A 377, 474 (1982) and radiative corrections of order O(alpha E_e/m_N) ~ 10^{-5} (Phys. Rev. D 99, 093006 (2019)), calculated as next--to--leading order corrections in the large nucleon mass $m_N$ expansion to Sirlin's corrections of order O(alpha/pi) (Phys. Rev. 164, 1767 (1967)), the corrections of order O(E^2_e/m^2_N) ~ 10^{-5} provide an improved level of precision of the theoretical background of the neutron beta decay, calculated in the Standard Model, for experimental searches of contributions of interactions beyond the Standard Model.


title: Acoustic echo cancellation with the dual-signal transformation LSTM network

abstract: This paper applies the dual-signal transformation LSTM network (DTLN) to the task of real-time acoustic echo cancellation (AEC). The DTLN combines a short-time Fourier transformation and a learned feature representation in a stacked network approach, which enables robust information processing in the time-frequency and in the time domain, which also includes phase information. The model is only trained on 60~h of real and synthetic echo scenarios. The training setup includes multi-lingual speech, data augmentation, additional noise and reverberation to create a model that should generalize well to a large variety of real-world conditions. The DTLN approach produces state-of-the-art performance on clean and noisy echo conditions reducing acoustic echo and additional noise robustly. The method outperforms the AEC-Challenge baseline by 0.30 in terms of Mean Opinion Score (MOS).


title: On Minimum Generalized Manhattan Connections

abstract: We consider minimum-cardinality Manhattan connected sets with arbitrary demands: Given a collection of points $P$ in the plane, together with a subset of pairs of points in $P$ (which we call demands), find a minimum-cardinality superset of $P$ such that every demand pair is connected by a path whose length is the $\ell_1$-distance of the pair. This problem is a variant of three well-studied problems that have arisen in computational geometry, data structures, and network design: (i) It is a node-cost variant of the classical Manhattan network problem, (ii) it is an extension of the binary search tree problem to arbitrary demands, and (iii) it is a special case of the directed Steiner forest problem. Since the problem inherits basic structural properties from the context of binary search trees, an $O(\log n)$-approximation is trivial. We show that the problem is NP-hard and present an $O(\sqrt{\log n})$-approximation algorithm. Moreover, we provide an $O(\log\log n)$-approximation algorithm for complete bipartite demands as well as improved results for unit-disk demands and several generalizations. Our results crucially rely on a new lower bound on the optimal cost that could potentially be useful in the context of BSTs.


title: Geometric quantization and quantum moment maps on coadjoint orbits and K\"ahler-Einstein manifolds

abstract: Deformation quantization and geometric quantization on K\"ahler manifolds give the mathematical description of the algebra of quantum observables and the Hilbert spaces respectively, where the later forms a representation of quantum observables asymptotically via Toeplitz operators. When there is a Hamiltonian $G$-action on a K\"ahler manifold, there are associated symmetries on both the quantum algebra and representation aspects. We show that in nice cases of coadjoint orbits and K\"ahler-Einstein manifolds, these symmetries are strictly compatible (not only asymptotically).


title: Nonparametric estimation of highest density regions for COVID-19

abstract: Highest density regions refer to level sets containing points of relatively high density. Their estimation from a random sample, generated from the underlying density, allows to determine the clusters of the corresponding distribution. This task can be accomplished considering different nonparametric perspectives. From a practical point of view, reconstructing highest density regions can be interpreted as a way of determining hot-spots, a crucial task for understanding COVID-19 space-time evolution. In this work, we compare the behavior of classical plug-in methods and a recently proposed hybrid algorithm for highest density regions estimation through an extensive simulation study. Both methodologies are applied to analyze a real data set about COVID-19 cases in the United States.


title: Linear Stochastic Dyadic model

abstract: We discuss a stochastic interacting particles' system connected to dyadic models of turbulence, defining suitable classes of solutions and proving their existence and uniqueness. We investigate the regularity of a particular family of solutions, called moderate, and we conclude with existence and uniqueness of invariant measures associated with such moderate solutions.


title: Listener's Social Identity Matters in Personalised Response Generation

abstract: Personalised response generation enables generating human-like responses by means of assigning the generator a social identity. However, pragmatics theory suggests that human beings adjust the way of speaking based on not only who they are but also whom they are talking to. In other words, when modelling personalised dialogues, it might be favourable if we also take the listener's social identity into consideration. To validate this idea, we use gender as a typical example of a social variable to investigate how the listener's identity influences the language used in Chinese dialogues on social media. Also, we build personalised generators. The experiment results demonstrate that the listener's identity indeed matters in the language use of responses and that the response generator can capture such differences in language use. More interestingly, by additionally modelling the listener's identity, the personalised response generator performs better in its own identity.


title: Learning to Infer Unseen Attribute-Object Compositions

abstract: The composition recognition of unseen attribute-object is critical to make machines learn to decompose and compose complex concepts like people. Most of the existing methods are limited to the composition recognition of single-attribute-object, and can hardly distinguish the compositions with similar appearances. In this paper, a graph-based model is proposed that can flexibly recognize both single- and multi-attribute-object compositions. The model maps the visual features of images and the attribute-object category labels represented by word embedding vectors into a latent space. Then, according to the constraints of the attribute-object semantic association, distances are calculated between visual features and the corresponding label semantic features in the latent space. During the inference, the composition that is closest to the given image feature among all compositions is used as the reasoning result. In addition, we build a large-scale Multi-Attribute Dataset (MAD) with 116,099 images and 8,030 composition categories. Experiments on MAD and two other single-attribute-object benchmark datasets demonstrate the effectiveness of our approach.


title: Fast computation of angular power spectra and covariances of high-resolution cosmic microwave background maps using the Toeplitz approximation

abstract: We present a simple approximation that can speed up the computation of the mode-coupling matrices, which are usually the bottleneck for computing unbiased angular power spectra, as well as their associated covariance matrices, of the cosmic microwave background temperature and polarization anisotropies. The approximation results in the speed up of the MASTER algorithm by more than an order of magnitude with very little loss of precision. We demonstrate the performance on simulations of forthcoming cosmic microwave background surveys such as the Simons Observatory and CMB-S4 for a wide variety of survey window functions.


title: Supreme Torsion Forms

abstract: We study formally real, non-pythagorean fields which have an anisotropic torsion form that contains every anisotropic torsion form as a subform. We obtain consequences for certain invariants and the Witt ring of such fields and construct examples. We obtain a theory analogous to the theory of supreme Pfister forms introduced by Karim Becher and see examples in which the Pythagoras number for formally real fields behaves like the level for nonreal fields.


title: A multi-technique approach to understanding delithiation damage in LiCoO2 thin films

abstract: We report on the delithiation of LiCoO2 thin films using oxalic acid (C2H2O4) with the goal of understanding the structural degradation of an insertion oxide associated with Li chemical extraction. Using a multi-technique approach that includes synchrotron radiation x-ray diffraction, scanning electron microscopy, micro Raman spectroscopy, photoelectron spectroscopy and conductive atomic force microscopy we reveal the balance between selective Li extraction and structural damage. We identify three different delithiation regimes, related to surface processes, bulk delithiation and damage generation. We find that only a fraction of the grains is affected by the delithiation process, which may create local inhomogeneities. The chemical route to Li extraction provides additional opportunities to investigate delithiation while avoiding the complications associated with electrolyte breakdown and could simplify in situ measurements.


title: Gravitational Dyonic Amplitude at One-Loop and its Inconsistency with the Classical Impulse

abstract: The recent proposal [1911.06318, 2010.07861] of implementing electric-magnetic duality rotation at the level of perturbative scattering amplitudes and its generalisation to gravitational context where usual gravitational mass is rotated to the NUT parameter of the Taub-NUT spacetime opens up an interesting avenue for studying NUT-charged objects as dynamical entities, in contrast to the usual approach where NUT-charged objects are considered as a static background. We extend the tree-order analysis to one-loop order and find a disagreement between geodesic motion on Taub-NUT background and impulse computation of scattering amplitudes. As a by-product of our analysis, we find a relation between tidal response parameters and resonance excitation parameters in the language of quantum field theory scattering amplitudes.


title: Numerical scheme for evaluating the collision integrals for triple interactions in relativistic plasma

abstract: Binary interactions in relativistic plasma, such as Coulomb and Compton scattering as well as pair creation and annihilation are well known and studied in details. Triple interactions, namely relativistic bremsstrahlung, double Compton scattering, radiative pair production, triple pair production/annihilation and their inverse processes, are usually considered as emission processes in astrophysical problems, as well as in laboratory plasmas. Their role in plasma kinetics is fundamental [1]. We present a new conservative scheme for computation of Uehling-Uhlenbeck collision integral for all triple interactions in relativistic plasma based on direct integration of exact QED matrix elements. Reaction rates for thermal distributions are compared, where possible, with the corresponding analytic expressions, showing good agreement. Our results are relevant for quantitative description of relativistic plasmas out of equilibrium, both in astrophysical and laboratory conditions.


title: Star edge-coloring of some special graphs

abstract: The star chromatic index of a multigraph $G$, denoted by $\chi_{\mathrm{star}}'(G)$, is the minimum number of colors needed to properly color the edges of $G$ such that no path or cycle of length $4$ is bicolored. In this paper, we study the star edge-coloring of Halin graphs, $k$-power graphs and the generalized Petersen graphs $P(3n, n)$.


title: CT Reconstruction with PDF: Parameter-Dependent Framework for Multiple Scanning Geometries and Dose Levels

abstract: Current mainstream of CT reconstruction methods based on deep learning usually needs to fix the scanning geometry and dose level, which will significantly aggravate the training cost and need more training data for clinical application. In this paper, we propose a parameter-dependent framework (PDF) which trains data with multiple scanning geometries and dose levels simultaneously. In the proposed PDF, the geometry and dose level are parameterized and fed into two multi-layer perceptrons (MLPs). The MLPs are leveraged to modulate the feature maps of CT reconstruction network, which condition the network outputs on different scanning geometries and dose levels. The experiments show that our proposed method can obtain competing performance similar to the original network trained with specific geometry and dose level, which can efficiently save the extra training cost for multiple scanning geometries and dose levels.


title: Real-space observation of vibrational strong coupling between propagating phonon polaritons and organic molecules

abstract: Phonon polaritons (PPs) in van der Waals (vdW) materials can strongly enhance light-matter interactions at mid-infrared frequencies, owing to their extreme infrared field confinement and long lifetimes. PPs thus bear potential for achieving vibrational strong coupling (VSC) with molecules. Although the onset of VSC has recently been observed spectroscopically with PP nanoresonators, no experiments so far have resolved VSC in real space and with propagating modes in unstructured layers. Here, we demonstrate by real-space nanoimaging that VSC can be achieved between propagating PPs in thin vdW crystals (specifically h-BN) and molecular vibrations in adjacent thin molecular layers. To that end, we performed near-field polariton interferometry, showing that VSC leads to the formation of a propagating hybrid mode with a pronounced anti-crossing region in its dispersion, in which propagation with negative group velocity is found. Numerical calculations predict VSC for nanometer-thin molecular layers and PPs in few-layer vdW materials, which could make propagating PPs a promising platform for ultra-sensitive on-chip spectroscopy and strong coupling experiments.


title: Certifying the adiabatic preparation of ultracold lattice bosons in the vicinity of the Mott transition

abstract: We present a joint experimental and theoretical analysis to assess the adiabatic experimental preparation of ultracold bosons in optical lattices aimed at simulating the three-dimensional Bose-Hubbard model. Thermometry of lattice gases is realized from the superfluid to the Mott regime by combining the measurement of three-dimensional momentum-space densities with ab-initio quantum Monte Carlo (QMC) calculations of the same quantity. The measured temperatures across the superfluid-to-Mott transition are in agreement with isentropic lines reconstructed via QMC for the experimental parameters of interest, with a conserved entropy per particle of $S/N=0.8(1) k_{B}$. In addition, the Fisher information associated with this thermometry method shows that the latter is most accurate in the critical regime close to the Mott transition, as confirmed in the experiment. These results prove that equilibrium states of the Bose-Hubbard model - including those in the quantum-critical regime above the Mott transition - can be adiabatically prepared in cold-atom apparatus.


title: Implementing efficient balanced networks with mixed-signal spike-based learning circuits

abstract: Efficient Balanced Networks (EBNs) are networks of spiking neurons in which excitatory and inhibitory synaptic currents are balanced on a short timescale, leading to desirable coding properties such as high encoding precision, low firing rates, and distributed information representation. It is for these benefits that it would be desirable to implement such networks in low-power neuromorphic processors. However, the degree of device mismatch in analog mixed-signal neuromorphic circuits renders the use of pre-trained EBNs challenging, if not impossible. To overcome this issue, we developed a novel local learning rule suitable for on-chip implementation that drives a randomly connected network of spiking neurons into a tightly balanced regime. Here we present the integrated circuits that implement this rule and demonstrate their expected behaviour in low-level circuit simulations. Our proposed method paves the way towards a system-level implementation of tightly balanced networks on analog mixed-signal neuromorphic hardware. Thanks to their coding properties and sparse activity, neuromorphic electronic EBNs will be ideally suited for extreme-edge computing applications that require low-latency, ultra-low power consumption and which cannot rely on cloud computing for data processing.


title: On the Cauchy problem for the wave equation in a two-dimensional domain with data on the boundary

abstract: The subject of the paper is the Cauchy problem for the wave equation in a space-time cylinder $\Omega\times{\mathbb R}$, $\Omega\subset{\mathbb R}^2$, with the data on the surface $\partial\Omega\times I$, where $I$ is a finite time interval. The algorithm for solving the Cauchy problem with the data on $S\times I$, $S\subset\partial\Omega$, was obtained previously. Here we adapt this algorithm to the special case $S=\partial\Omega$ and show that in this situation, the solution is determined with higher stability in comparison with the case $S\subsetneqq\partial\Omega$.


title: Random walks and community detection in hypergraphs

abstract: We propose a one parameter family of random walk processes on hypergraphs, where a parameter biases the dynamics of the walker towards hyperedges of low or high cardinality. We show that for each value of the parameter the resulting process defines its own hypergraph projection on a weighted network. We then explore the differences between them by considering the community structure associated to each random walk process. To do so, we generalise the Markov stability framework to hypergraphs and test it on artificial and real-world hypergraphs.

categories: cond-mat.stat-mech cs.SI math.DS physics.soc-ph


title: Upsampling artifacts in neural audio synthesis

abstract: A number of recent advances in neural audio synthesis rely on upsampling layers, which can introduce undesired artifacts. In computer vision, upsampling artifacts have been studied and are known as checkerboard artifacts (due to their characteristic visual pattern). However, their effect has been overlooked so far in audio processing. Here, we address this gap by studying this problem from the audio signal processing perspective. We first show that the main sources of upsampling artifacts are: (i) the tonal and filtering artifacts introduced by problematic upsampling operators, and (ii) the spectral replicas that emerge while upsampling. We then compare different upsampling layers, showing that nearest neighbor upsamplers can be an alternative to the problematic (but state-of-the-art) transposed and subpixel convolutions which are prone to introduce tonal artifacts.


title: Carnahan Starling type equations of state for stable hard disk and hard sphere fluids

abstract: The well-known Carnahan-Starling (CS) equation of state (EoS) [1] for the hard sphere (HS) fluid was derived from a quadratic relation between the integer portions of the virial coefficients, Bn, and their orders, n. Here we extend the method to the full virial coefficients Bn for the general D-dimensional case. We assume a polynomial function of (D-1)th order for the virial coefficients starting from n=4 and EoS are derived from it. For the hard rob (D=1) case, the exact solution is obtained. For the stable hard disk fluid (D=2), the most recent virial coefficients up to the 10th [2] and accurate compressibility data[3,4] are employed to construct and test the EoS. For the stable hard sphere (D=3) fluid, a new CS-type EoS is constructed and tested with the most recent virial coefficients [5,2] up to the 11th and with the highly-accurate simulation data for compressibility [6-8]. The simple new EoS turn out to be as accurate as the highest-level Pade approximations based on all available virial coefficients, and significantly improve the CS-type EoS in the hard sphere case. We also shown that as long as the virial coefficients obey a polynomial function any EoS derived from it will diverge at the non-physical packing fraction=1.


title: A Data-Driven Sparse Polynomial Chaos Expansion Method to Assess Probabilistic Total Transfer Capability for Power Systems with Renewables

abstract: The increasing uncertainty level caused by growing renewable energy sources (RES) and aging transmission networks poses a great challenge in the assessment of total transfer capability (TTC) and available transfer capability (ATC). In this paper, a novel data-driven sparse polynomial chaos expansion (DDSPCE) method is proposed for estimating the probabilistic characteristics (e.g., mean, variance, probability distribution) of probabilistic TTC (PTTC). Specifically, the proposed method, requiring no pre-assumed probabilistic distributions of random inputs, exploits data sets directly in estimating the PTTC. Besides, a sparse scheme is integrated to improve the computational efficiency. Numerical studies on the modified IEEE 118-bus system demonstrate that the proposed DDSPCE method can achieve accurate estimation for the probabilistic characteristics of PTTC with a high efficiency. Moreover, numerical results reveal the great significance of incorporating discrete random inputs in PTTC and ATC assessment, which nevertheless was not given sufficient attention.


title: Improved Inference of Gaussian Mixture Copula Model for Clustering and Reproducibility Analysis using Automatic Differentiation

abstract: Copulas provide a modular parameterization of multivariate distributions that decouples the modeling of marginals from the dependencies between them. Gaussian Mixture Copula Model (GMCM) is a highly flexible copula that can model many kinds of multi-modal dependencies, as well as asymmetric and tail dependencies. They have been effectively used in clustering non-Gaussian data and in Reproducibility Analysis, a meta-analysis method designed to verify the reliability and consistency of multiple high-throughput experiments. Parameter estimation for GMCM is challenging due to its intractable likelihood. The best previous methods have maximized a proxy-likelihood through a Pseudo Expectation Maximization (PEM) algorithm. They have no guarantees of convergence or convergence to the correct parameters. In this paper, we use Automatic Differentiation (AD) tools to develop a method, called AD-GMCM, that can maximize the exact GMCM likelihood. In our simulation studies and experiments with real data, AD-GMCM finds more accurate parameter estimates than PEM and yields better performance in clustering and Reproducibility Analysis. We discuss the advantages of an AD-based approach, to address problems related to monotonic increase of likelihood and parameter identifiability in GMCM. We also analyze, for GMCM, two well-known cases of degeneracy of maximum likelihood in GMM that can lead to spurious clustering solutions. Our analysis shows that, unlike GMM, GMCM is not affected in one of the cases.


title: Identification of strongly interacting organic semimetals

abstract: Dirac- and Weyl point- and line-node semimetals are characterized by a zero band gap with simultaneously vanishing density of states. Given a sufficient interaction strength, such materials can undergo an interaction instability, e.g., into an excitonic insulator phase. Due to generically flat bands, organic crystals represent a promising materials class in this regard. We combine machine learning, density functional theory, and effective models to identify specific example materials. Without taking into account the effect of many-body interactions, we found the organic charge transfer salts (EDT-TTF-I$_2$)$_2$(DDQ)$\cdot($CH$_3$CN) and TSeF-TCNQ and a bis-1,2,3-dithiazolyl radical conductor to exhibit a semimetallic phase in our ab initio calculations. Adding the effect of strong particle-hole interactions for (EDT-TTF-I$_2$)$_2$(DDQ)$\cdot($CH$_3$CN) and TSeF-TCNQ opens an excitonic gap in the order of 60 meV and 100 meV, which is in good agreement with previous experiments on these materials.


title: Fourth-Order Nonlocal Tensor Decomposition Model for Spectral Computed Tomography

abstract: Spectral computed tomography (CT) can reconstruct spectral images from different energy bins using photon counting detectors (PCDs). However, due to the limited photons and counting rate in the corresponding spectral fraction, the reconstructed spectral images usually suffer from severe noise. In this paper, a fourth-order nonlocal tensor decomposition model for spectral CT image reconstruction (FONT-SIR) method is proposed. Similar patches are collected in both spatial and spectral dimensions simultaneously to form the basic tensor unit. Additionally, principal component analysis (PCA) is applied to extract latent features from the patches for a robust and efficient similarity measure. Then, low-rank and sparsity decomposition is performed on the produced fourth-order tensor unit, and the weighted nuclear norm and total variation (TV) norm are used to enforce the low-rank and sparsity constraints, respectively. The alternating direction method of multipliers (ADMM) is adopted to optimize the objective function. The experimental results with our proposed FONT-SIR demonstrates a superior qualitative and quantitative performance for both simulated and real data sets relative to several state-of-the-art methods, in terms of noise suppression and detail preservation.


title: Redesign of the ATLAS Tile Calorimeter link Daughterboard for the HL-LHC

abstract: The Phase-2 ATLAS upgrade for the High Luminosity Large Hadron Collider (HL-LHC) has motivated progressive redesigns of the ATLAS Tile Calorimeter (TileCal) read-out link and control board (Daughterboard). The Daughterboard (DB) communicates with the off-detector electronics via two 4.6 Gbps downlinks and two pairs of 9.6 Gbps uplinks. Configuration commands and LHC timing is received through the downlinks by two CERN radiation hard GBTx ASICs and propagated through Ultrascale+ FPGAs to the front-end. Simultaneously, the FPGAs send continuous high-speed readout of digitized PMT samples, slow control and monitoring data through the uplink. The design minimizes single points of failure and reduces sensitivity to SEUs and radiation damage by employing a double-redundant scheme, using Triple Mode Redundancy (TMR) and Xilinx Soft Error Mitigation (SEM) in the FPGAs, adopting Cyclic Redundancy Check (CRC) error verification in the uplinks and Forward Error Correction (FEC) in the downlinks. We present a DB redesign that brings an enhanced timing scheme, and improved radiation tolerance by mitigating Single Event Latch-up (SEL) induced errors and implementing a more robust power-up and current monitoring scheme.


title: Classical simulation of Gaussian quantum circuits with non-Gaussian input states

abstract: We consider Gaussian quantum circuits supplemented with non-Gaussian input states and derive sufficient conditions for efficient classical strong simulation of these circuits. In particular, we generalise the stellar representation of continuous-variable quantum states to the multimode setting and relate the stellar rank of the input non-Gaussian states, a recently introduced measure of non-Gaussianity, to the cost of evaluating classically the output probability densities of these circuits. Our results have consequences for the strong simulability of a large class of near-term continuous-variable quantum circuits.


title: Study on $\eta_{c2}(\eta_{b2})$ electromagnetic decay into double photons

abstract: Within the framework of nonrelativistic QCD (NRQCD) factorization formalism, we compute the helicity amplitude as well as the decay width of $\eta_{Q2}$ ($Q=c,b$) electromagnetic decay into two photons up to next-to-next-to-leading order (NNLO) in $\alpha_s$ expansion. For the first time, we verify the validity of NRQCD factorization for the D-wave quarkonium decay at NNLO. We find that the $\mathcal{O}(\alpha_s)$ and $\mathcal{O}(\alpha_s^2)$ corrections to the helicity amplitude are negative and moderate, nevertheless both corrections combine to suppress the leading-order prediction for the decay width significantly. By approximating the total decay width of $\eta_{Q2}$ as the sum of those for the hadronic decay and the electric $E1$ transition, we obtain the branching ratios ${\rm Br}(\eta_{c2}\to 2\gamma)\approx 5\times10^{-6}$ and ${\rm Br}(\eta_{b2}\to 2\gamma)\approx 4\times10^{-7}$. To explore the potential measurement on $\eta_{Q2}$, we further evaluate the production cross section of $\eta_{Q2}$ at LHCb at the lowest order in $\alpha_s$ expansion. With the kinematic constraint on the longitudinal rapidity $4.5>y>2$ and transverse momentum $P_T>(2-4)m_Q$ for $\eta_{Q2}$, we find the cross section can reach $2-50$ nb for $\eta_{c2}$, and $1-22$ pb for $\eta_{b2}$. Considering the integrated luminosity $\mathcal{L}=10\, {\rm fb}^{-1}$ at $\sqrt{s}=7$ TeV and $\sqrt{s}=13$ TeV, we estimate that there are several hundreds events of $pp\to \eta_{c2}\to 2\gamma$. Since the background is relatively clean, it is promising to reconstruct $\eta_{c2}$ through its electromagnetic decay. On the contrary, due to small branching ratio and production cross section, it is quite challenging to detect $\eta_{b2}\to 2\gamma$ at LHCb.


title: A Poisson limit theorem for Gibbs-Markov maps

abstract: We prove for Gibbs-Markov maps that the number of visits to a sequence of shrinking sets with bounded cylindrical lengths converges in distribution to a Poisson law. Applying to continued fractions, this result extends Doeblin's Poisson limit theorem.


title: Modeling and control of malaria dynamics in fish farming regions

abstract: In this work we propose a model that represents the relation between fish ponds, the mosquito population and the transmission of malaria. It has been observed that in the Amazonic region of Acre, in the North of Brazil, fish farming is correlated to the transmission of malaria when carried out in artificial ponds that become breeding sites. Evidence has been found indicating that cleaning the vegetation from the edges of the crop tanks helps to control the size of the mosquito population. We use our model to determine the effective contribution of fish farming practices on malaria transmission dynamics. The model consists of a nonlinear system of ordinary differential equations with jumps at the cleaning time, which act as impulsive controls. We study the asymptotic behaviour of the system in function of the intensity and periodicity of the cleaning, and the value of the parameters. In particular, we state sufficient conditions under which the mosquito population is eliminated or persists, and under which the malaria is eliminated or becomes endemic. We prove our conditions by applying results for cooperative systems with concave nonlinearities.


title: Simultaenous Sieves: A Deterministic Streaming Algorithm for Non-Monotone Submodular Maximization

abstract: In this work, we present a combinatorial, deterministic single-pass streaming algorithm for the problem of maximizing a submodular function, not necessarily monotone, with respect to a cardinality constraint (SMCC). In the case the function is monotone, our algorithm reduces to the optimal streaming algorithm of Badanidiyuru et al. (2014). In general, our algorithm achieves ratio $\alpha / (1 + \alpha) - \varepsilon$, for any $\varepsilon > 0$, where $\alpha$ is the ratio of an offline (deterministic) algorithm for SMCC used for post-processing. Thus, if exponential computation time is allowed, our algorithm deterministically achieves nearly the optimal $1/2$ ratio. These results nearly match those of a recently proposed, randomized streaming algorithm that achieves the same ratios in expectation. For a deterministic, single-pass streaming algorithm, our algorithm achieves in polynomial time an improvement of the best approximation factor from $1/9$ of previous literature to $\approx 0.2689$.


title: Destruction of eccentric planetesimals by ram pressure and erosion

abstract: Small, pebble-sized objects and large bodies of planetesimal size both play important roles in planet formation. They form the evolutionary steps of dust growth in their own respect. However, at later times, they are also thought to provide background populations of mass that larger bodies might feed upon. What we suggest in this work is that starting at times of viscous stirring, planetesimals on eccentric orbits could simply explode as they become supersonic in comparison to small, porous planetary bodies entering Earth's atmosphere. We present a toy model of planetesimal motion and destruction to show the key aspects of this process. The consequences are quite severe. At all times, it is shown that only planetesimals on more or less circular orbits exist in the inner disk. After the destruction of a planetesimal, the remaining matter is continuously redistributed to the pebble reservoir of the protoplanetary disk. Since destruction typically occurs at small stellar distances due to supersonic speeds, it is expected to boost pebble accretion in the inner protoplanetary disk as one of its main effects.


title: Constraining the density dependence of the symmetry energy: the isospin transport ratio revisited

abstract: The isospin diffusion of the quasi-projectile formed in the $^{64,58}Ni$ on $^{64,58}Ni$ reactions in the Fermi energy domain is investigated in the framework of the Boltzmann-Uehling-Uhlenbeck transport model. The well known isospin transport ratio observable is revisited, with the aim of insuring an optimal comparison between experimental data and theoretical calculations and reducing the present uncertainties in the extraction of empirical equation of state parameters. We show that isospin transport ratios are sensitive to all the low order isovector parameters ($E_{sym}$, $L_{sym}$ and $K_{sym}$). We demonstrate that realistic models of the equation of state, covering the uncertainty that presently affects the theoretical description of neutron stars static observables, can be effectively discriminated by isospin diffusion experiments, provided the neutron to proton ratio of the projectile remnant is precisely measured as a function of centrality.


title: Quantum Size Effects in the Magnetic Susceptibility of a Metallic Nanoparticle

abstract: We theoretically study quantum size effects in the magnetic response of a spherical metallic nanoparticle (e.g. gold). Using the Jellium model in spherical coordinates, we compute the induced magnetic moment and the magnetic susceptibility for a nanoparticle in the presence of a static external magnetic field. Below a critical magnetic field the magnetic response is diamagnetic, whereas above such field the magnetization is characterized by sharp, step-like increases of several tenths of Bohr magnetons, associated with the Zeeman crossing of energy levels above and below the Fermi sea. We quantify the robustness of these regimes against thermal excitations and finite linewidth of the electronic levels. Finally, we propose two methods for experimental detection of the quantum size effects based on the coupling to superconducting quantum interference devices.


title: A rigid, not infinitesimally rigid surface with K ample

abstract: We produce an example of a rigid, but not infinitesimally rigid smooth compact complex surface with ample canonical bundle using results about arrangements of lines inspired by work of Hirzebruch, Kapovich and Millson, Manetti and Vakil.


title: Measuring Dark Matter in Galaxies: The Mass Fraction Within 5 Effective Radii

abstract: Large galaxies may contain an "atmosphere" of hot interstellar X-ray gas, and the temperature and radial density profile of this gas can be used to measure the total mass of the galaxy contained within a given radius r. We use this technique for 102 early-type galaxies (ETGs) with stellar masses M_* > 10^10 M_Sun, to evaluate the mass fraction of dark matter (DM) within the fiducial radius r = 5 r_e, denoted f_5 = f_{DM}(5r_e). On average, these systems have a median f_5 = 0.8 - 0.9 with a typical galaxy-to-galaxy scatter +-0.15. Comparisons with mass estimates made through the alternative techniques of satellite dynamics (e.g. velocity distributions of globular clusters, planetary nebulae, satellite dwarfs) as well as strong lensing show encouraging consistency over the same range of stellar mass. We find that many of the disk galaxies (S0/SA0/SB0) have a significantly higher mean $f_5$ than do the pure ellipticals, by Delta f_5 = 0.1. We suggest that this higher level may be a consequence of sparse stellar haloes and quieter histories with fewer major episodes of feedback or mergers. Comparisons are made with the Magneticum Pathfinder suite of simulations for both normal and centrally dominant "Brightest Cluster" galaxies. Though the observed data exhibit somewhat larger scatter at a given galaxy mass than do the simulations, the mean level of DM mass fraction for all classes of galaxies is in good first-order agreement with the simulations. Lastly, we find that the group galaxies with stellar masses near M_* ~ 10^11 M_Sun have relatively more outliers at low $f_5$ than in other mass ranges, possibly the result of especially effective AGN feedback in that mass range leading to expansion of their dark matter halos.


title: Matrix Engines for High Performance Computing:A Paragon of Performance or Grasping at Straws?

abstract: Matrix engines or units, in different forms and affinities, are becoming a reality in modern processors; CPUs and otherwise. The current and dominant algorithmic approach to Deep Learning merits the commercial investments in these units, and deduced from the No.1 benchmark in supercomputing, namely High Performance Linpack, one would expect an awakened enthusiasm by the HPC community, too. Hence, our goal is to identify the practical added benefits for HPC and machine learning applications by having access to matrix engines. For this purpose, we perform an in-depth survey of software stacks, proxy applications and benchmarks, and historical batch job records. We provide a cost-benefit analysis of matrix engines, both asymptotically and in conjunction with state-of-the-art processors. While our empirical data will temper the enthusiasm, we also outline opportunities to misuse these dense matrix-multiplication engines if they come for free.


title: Explainable Machine Learning for Public Policy: Use Cases, Gaps, and Research Directions

abstract: Explainability is highly-desired in Machine Learning (ML) systems supporting high-stakes policy decisions in areas such as health, criminal justice, education, and employment. While the field of explainable ML has expanded in recent years, much of this work has not taken real-world needs into account. A majority of proposed methods are designed with \textit{generic} explainability goals without well-defined use-cases or intended end-users and evaluated on simplified tasks, benchmark problems/datasets, or with proxy users (e.g., AMT). We argue that these simplified evaluation settings do not capture the nuances and complexities of real-world applications. As a result, the applicability and effectiveness of this large body of theoretical and methodological work in real-world applications are unclear. In this work, we take steps toward addressing this gap for the domain of public policy. First, we identify the primary use-cases of explainable ML within public policy problems. For each use case, we define the end-users of explanations and the specific goals the explanations have to fulfill. Finally, we map existing work in explainable ML to these use-cases, identify gaps in established capabilities, and propose research directions to fill those gaps to have a practical societal impact through ML. The contribution is 1) a methodology for explainable ML researchers to identify use cases and develop methods targeted at them and 2) using that methodology for the domain of public policy and giving an example for the researchers on developing explainable ML methods that result in real-world impact.


title: E-Commerce Delivery Demand Modeling Framework for An Agent-Based Simulation Platform

abstract: The e-commerce delivery demand has grown rapidly in the past two decades and such trend has accelerated tremendously due to the ongoing coronavirus pandemic. Given the situation, the need for predicting e-commerce delivery demand and evaluating relevant logistics solutions is increasing. However, the existing simulation models for e-commerce delivery demand are still limited and do not consider the delivery options and their attributes that shoppers face on e-commerce order placements. We propose a novel modeling framework which jointly predicts the average total value of e-commerce purchase, the purchase amount per transaction, and delivery option choices. The proposed framework can simulate the changes in e-commerce delivery demand attributable to the changes in delivery options. We assume the model parameters based on various sources of relevant information and conduct a demonstrative sensitivity analysis. Furthermore, we have applied the model to the simulation for the Auto-Innovative Prototype city. While the calibration of the model using real-world survey data is required, the result of the analysis highlights the applicability of the proposed framework.


title: The DigitalTwin from an Artificial Intelligence Perspective

abstract: Services for Cyber-Physical Systems based on Artificial Intelligence and Machine Learning require a virtual representation of the physical. To reduce modeling efforts and to synchronize results, for each system, a common and unique virtual representation used by all services during the whole system life-cycle is needed, i.e. a DigitalTwin. In this paper such a DigitalTwin, namely the AI reference model AITwin, is defined. This reference model is verified by using a running example from process industry and by analyzing the work done in recent projects.


title: Designing optimal networks for multi-commodity transport problem

abstract: Designing and optimizing different flows in networks is a relevant problem in many contexts. While a number of methods have been proposed in the physics and optimal transport literature for the one-commodity case, we lack similar results for the multi-commodity scenario. In this paper we present a model based on optimal transport theory for finding optimal multi-commodity flow configurations on networks. This model introduces a dynamics that regulates the edge conductivities to achieve, at infinite times, a minimum of a Lyapunov functional given by the sum of a convex transport cost and a concave infrastructure cost. We show that the long time asymptotics of this dynamics are the solutions of a standard constrained optimization problem that generalizes the one-commodity framework. Our results provide new insights into the nature and properties of optimal network topologies. In particular, they show that loops can arise as a consequence of distinguishing different flow types, complementing previous results where loops, in the one-commodity case, were obtained as a consequence of imposing dynamical rules to the sources and sinks or when enforcing robustness to damage. Finally, we provide an efficient implementation of our model which convergences faster than standard optimization methods based on gradient descent.

categories: physics.soc-ph cs.SI cs.SY eess.SY nlin.AO