title: Magnetic structure and exchange interactions in the layered semiconductor CrPS4

abstract: Compounds with two-dimensional (2D) layers of magnetic ions weakly connected by van der Waals bonding offer routes to enhance quantum behavior, stimulating both fundamental and applied interest. CrPS4 is one such magnetic van der Waals material, however, it has undergone only limited investigation. Here we present a comprehensive series of neutron scattering measurements to determine the magnetic structure and exchange interactions. The observed magnetic excitations allow a high degree of constraint on the model parameters not normally associated with measurements on a powder sample. The results demonstrate the 2D nature of the magnetic interactions, while also revealing the importance of interactions along 1D chains within the layers. The subtle role of competing interactions is observed, which manifest in a non-trivial magnetic transition and a tunable magnetic structure in a small applied magnetic field through a spin-flop transition. Our results on the bulk compound provide insights that can be applied to an understanding of the behavior of reduced layer CrPS4.


title: On spinning loop amplitudes in Anti-de Sitter space

abstract: In this work we present a systematic study of AdS$_{d+1}$ loop amplitudes for gluons and gravitons using momentum space techniques. Inspired by the recent progress in tree level computation, we construct a differential operator that can act on a scalar factor in order to generate gluon and graviton loop integrands: this systematizes the computation for any given loop level Witten diagram. We then give a general prescription in this formalism, and discuss it for bubble, triangle, and box diagrams.


title: Universal location of the Yang-Lee edge singularity in O(N) theories

abstract: We determine a previously unknown universal quantity, the location of the Yang-Lee edge singularity for the O($N$) theories in a wide range of $N$ and various dimensions. At large $N$, we reproduce the $N\to\infty$ analytical result on the location of the singularity and, additionally, we obtain the mean-field result for the location in $d=4$ dimensions. In order to capture the nonperturbative physics for arbitrary $N$, $d$ and complex-valued external fields, we use the functional renormalization group approach.


title: A Multiscale Graph Convolutional Network Using Hierarchical Clustering

abstract: The information contained in hierarchical topology, intrinsic to many networks, is currently underutilised. A novel architecture is explored which exploits this information through a multiscale decomposition. A dendrogram is produced by a Girvan-Newman hierarchical clustering algorithm. It is segmented and fed through graph convolutional layers, allowing the architecture to learn multiple scale latent space representations of the network, from fine to coarse grained. The architecture is tested on a benchmark citation network, demonstrating competitive performance. Given the abundance of hierarchical networks, possible applications include quantum molecular property prediction, protein interface prediction and multiscale computational substrates for partial differential equations.


title: Two-photon exchange from intermediate state resonances in elastic electron-proton scattering

abstract: We use a recently developed dispersive approach to compute the two-photon exchange (TPE) correction to elastic electron-proton scattering, including contributions from hadronic $J^P=1/2^\pm$ and $3/2^\pm$ resonant intermediate states below~1.8~GeV. For the transition amplitudes from the proton ground state to the resonant excited states we employ new exclusive meson electroproduction data from CLAS at $Q^2 \lesssim 5$~GeV$^2$, and we explore the effects of both fixed and dynamic widths for the resonances. Among the resonant states, the $N(1520)~\!3/2^-$ becomes dominant for $Q^2 \gtrsim 2$~GeV$^2$, with a sign opposite to the comparably sized $\Delta(1232)~\!3/2^+$ contribution, leading to an overall increase in the size of the TPE correction to the cross section relative to the nucleon only contribution at higher $Q^2$ values. The results are in good overall agreement with recent $e^+ p$ to $e^- p$ cross section ratio and polarization transfer measurements, and provide compelling evidence for a resolution of the electric to magnetic form factor ratio discrepancy.


title: Asymptotic Analysis of a Two-Phase Model of Solid Tumour Growth

abstract: We investigate avascular tumour growth as a two-phase process consisting of cells and liquid. Based on the one-dimensional continuum moving-boundary model formulated by (Byrne, King, McElwain, Preziosi, Applied Mathematics Letters, 2003, 16, 567-573), we defined boundary conditions for the analogous model of tumour growth in two dimensions. We investigate linear stability of one dimensional time-dependent solution profiles in the moving-boundary formulation of a limit case (with negligible nutrient consumption and cell drag). For this, we obtain an asymptotic limit of the two-dimensional perturbations for large time (in the case where the tumour is growing) by using the method of matched asymptotic approximations. Having characterised an asymptotic limit of the perturbations, we compare it to the time-dependent solution profile in order to analytically obtain a condition for instability. Numerical simulations are mentioned.


title: A variant of Cauchy's argument principle for analytic functions which applies to curves containing zeroes

abstract: It is known that the Cauchy's argument principle, applied to an holomorphic function $f$, requires that $f$ has no zeros on the curve of integration. In this short note, we give a generalization of such a principle which covers the case when $f$ has zeros on the curve, as well as an application.


title: The Riemann Hypothesis is false

abstract: Let $\Theta$ denote the supremum of the real parts of the zeros of the Riemann zeta function. We demonstrate that $\Theta=1$, which entails the existence of infinitely many Riemann zeros off the critical line (thus disproving the Riemann Hypothesis (RH), which asserts that $\Theta = \frac{1}{2}$). The paper is concluded by a brief discussion of why our argument doesn't work for both Weil and Beurling zeta functions whose analogues of the RH are known to be true.


title: The Mondrian Puzzle: A Bound Concerning the $M(n) = 0$ Case

abstract: In response to the Numberphile video regarding the Mondrian Puzzle https://www.youtube.com/watch?v=49KvZrioFB0, we provide a lower bound on how many integers less than a given threshold $x$ satisfy $M(n) \neq 0$ where $M(n)$ is the quantity in which the Mondrian Puzzle is interested, i.e. the minimal difference in area between the largest and smallest rectangle in a set of incongruent, integer-sided rectangles which tile an $n$ by $n$ square.


title: Some Partial Fraction Identities associated with the Cyclotomic Polynomials

abstract: We establish some partial fraction identities for rational functions whose denominators are implicit products of the cyclotomic polynomials. To achieve this, we first develop a general algebraic approach for partial fraction decomposition inspired by the Heaviside's cover-up method. We thus call our method the Extended Cover-Up Method. Using our method we obtain direct formulas for $q$-partial fractions for certain generating functions. As a direct consequence of our formulas one can compute the Sylvester denumerants, the Frobenius number and the Ehrhart polynomials in pseudo-polynomial time. Further, we provide a framework for a generalization of the Fourier-Dedekind sum and their associated Rademacher reciprocity theorem extending the results of Carlitz, Zagier and Gessel. By performing a Fourier analysis we demonstrate that our extended cover-up method explains in simple terms the mechanism behind the reciprocity law.


title: Optimizing Downlink Resource Allocation in Multiuser MIMO Networks via Fractional Programming and the Hungarian Algorithm

abstract: Optimizing the sum-log-utility for the downlink of multi-frequency band, multiuser, multiantenna networks requires joint solutions to the associated beamforming and user scheduling problems through the use of cloud radio access network (CRAN) architecture; optimizing such a network is, however, non-convex and NP-hard. In this paper, we present a novel iterative beamforming and scheduling strategy based on fractional programming and the Hungarian algorithm. The beamforming strategy allows us to iteratively maximize the chosen objective function in a fashion similar to block coordinate ascent. Furthermore, based on the crucial insight that, in the downlink, the interference pattern remains fixed for a given set of beamforming weights, we use the Hungarian algorithm as an efficient approach to optimally schedule users for the given set of beamforming weights. Specifically, this approach allows us to select the best subset of users (amongst the larger set of all available users). Our simulation results show that, in terms of average sum-log-utility, as well as sum-rate, the proposed scheme substantially outperforms both the state-of-the-art multicell weighted minimum mean-squared error (WMMSE) and greedy proportionally fair WMMSE schemes, as well as standard interior-point and sequential quadratic solvers. Importantly, our proposed scheme is also far more computationally efficient than the multicell WMMSE scheme.


title: Modified Gravity (MOG), Cosmology and Black Holes

abstract: A covariant modified gravity (MOG) is formulated by adding to general relativity two new degrees of freedom, a scalar field gravitational coupling strength $G= 1/\chi$ and a gravitational spin 1 vector field $\phi_\mu$. The $G$ is written as $G=G_N(1+\alpha)$ where $G_N$ is Newton's constant, and the gravitational source charge for the vector field is $Q_g=\sqrt{\alpha G_N}M$, where $M$ is the mass of a body. Cosmological solutions of the theory are derived in a homogeneous and isotropic cosmology. Black holes in MOG are stationary as the end product of gravitational collapse and are axisymmetric solutions with spherical topology. It is shown that the scalar field $\chi$ is constant everywhere for an isolated black hole with asymptotic flat boundary condition. A consequence of this is that the scalar field loses its monopole moment radiation.


title: PICO: Primitive Imitation for COntrol

abstract: In this work, we explore a novel framework for control of complex systems called Primitive Imitation for Control PICO. The approach combines ideas from imitation learning, task decomposition, and novel task sequencing to generalize from demonstrations to new behaviors. Demonstrations are automatically decomposed into existing or missing sub-behaviors which allows the framework to identify novel behaviors while not duplicating existing behaviors. Generalization to new tasks is achieved through dynamic blending of behavior primitives. We evaluated the approach using demonstrations from two different robotic platforms. The experimental results show that PICO is able to detect the presence of a novel behavior primitive and build the missing control policy.


title: Role of silicon and carbon on the structural and electrochemical properties of Si-Ni$_{3.4}$Sn$_4$-Al-C anodes for Li-ion batteries

abstract: Varying the amounts of silicon and carbon, different composites have been prepared by ball milling of Si, Ni$_{3.4}$Sn$_4$, Al and C. Silicon and carbon contents are varied from 10 to 30 wt.% Si, and 0 to 20 wt.% C. The microstructural and electrochemical properties of the composites have been investigated by X-Ray Diffraction (XRD), Scanning Electron Microscopy (SEM) and electrochemical galvanostatic cycling up to 1000 cycles. Impact of silicon and carbon contents on the phase occurrence, electrochemical capacity and cycle-life are compared and discussed. For C-content comprised between 9 and 13 wt.% and Si-content >= 20 wt.%, Si nanoparticles are embedded in a Ni$_{3.4}$Sn$_4$-Al-C matrix which is chemically homogeneous at the micrometric scale. For other carbon contents and low Si-amount (10 wt.%), no homogeneous matrix is formed around Si nanoparticles. When homogenous matrix is formed, both Ni$_3$Sn$_4$ and Si participate to the reversible lithiation mechanism, whereas no reaction between Ni$_3$Sn$_4$ and Li is observed for no homogenous matrix. Moreover, best cycle-life performances are obtained when Si nanoparticles are embedded in a homogenous matrix and Si-content is moderate (<= 20 wt.%). Composites with carbon in the 9-13 wt.% range and 20 wt.% silicon lead to the best balance between capacity and life duration upon cycling. This work experimentally demonstrates that embedding Si in an intermetallic/carbon matrix allows to efficiently accommodate Si volume changes on cycling to ensure long cycle-life.


title: Insights into nature of a magnetization plateau of 3$d$-4$f$ coordination polymer [Dy$_2$Cu$_2$]$_n$ from a spin-1/2 Ising-Heisenberg orthogonal-dimer chain

abstract: The ground state and magnetization process of an exactly solved spin-$1/2$ Ising-Heisenberg orthogonal-dimer chain with two different gyromagnetic factors of the Ising and Heisenberg spins are investigated in detail. It is shown that the investigated quantum spin chain exhibits up to seven possible ground states depending on a mutual interplay of the magnetic field, intra- and inter-dimer coupling constants. More specifically, the frustrated and modulated quantum antiferromagnetic phases are responsible in zero-temperature magnetization curves for a zero magnetization plateau. The intermediate 1/11- and 5/11-plateaus emerge due to the frustrated and modulated quantum ferrimagnetic phases, while the intermediate 9/11- and 10/11-plateaus can be attributed to the quantum and classical ferrimagnetic phases. It is conjectured that the magnetization plateau experimentally observed in a high-field magnetization curve of 3$d$-4$f$ heterobimetallic coordination polymer [\{Dy(hfac)$_2$(CH$_3$OH)\}$_2$\{Cu(dmg)(Hdmg)\}$_2$]$_n$ (H$_2$dmg $=$ dimethylglyoxime; Hhfac $=$ 1,1,1,5,5,5-hexafluoropentane-2,4-dione) could be attributed to the classical and quantum ferrimagnetic phases.


title: Scratching the surface: Elastic rotations beneath nanoscratch and nanoindentation tests

abstract: In this paper, we investigate the residual deformation field in the vicinity of nano-scratch tests using two orientations of a Berkovich tip on an (001) Cu single crystal. We compare the deformation with that from indentation, in an attempt to understand the mechanisms of deformation in tangential sliding. The lattice rotation fields are mapped experimentally using high-resolution electron backscatter diffraction (HR-EBSD) on cross-sections prepared using focused ion beam (FIB). A physically-based crystal plasticity finite element model (CPFEM) is used to simulate the lattice rotation fields, and provide insight into the 3D rotation field surrounding nano-scratch experiments, as it transitions from an initial static indentation to a steady-state scratch. The CPFEM simulations capture the experimental rotation fields with good fidelity, and show how the rotations about the scratch direction are reversed as the indenter moves away from the initial indentation.


title: IXmon: Detecting and Analyzing DRDoS Attacks at Internet Exchange Points

abstract: Distributed reflective denial of service (DRDoS) attacks are a popular choice among adversaries. In fact, one of the largest DDoS attacks ever recorded, reaching a peak of 1.3Tbps against GitHub, was a memcached-based DRDoS attack. More recently, a record-breaking 2.3Tbps attack against Amazon AWS was due to a CLDAP-based DRDoS attack. Although reflective attacks have been known for years, DRDoS attacks are unfortunately still popular and largely unmitigated. In this paper, we study in-the-wild DRDoS attacks observed from a large Internet exchange point (IXP) and provide a number of security-relevant measurements and insights. To enable this study, we first developed IXmon, an open-source DRDoS detection system specifically designed for deployment at large IXP-like network connectivity providers and peering hubs. We deployed IXmon at Southern Crossroads (SoX), an IXP-like hub that provides both peering and upstream Internet connectivity services to more than 20 research and education (R&E) networks in the South-East United States. In a period of about 21 months, IXmon detected more than 900 DRDoS attacks towards 31 different victim ASes. An analysis of the real-world DRDoS attacks detected by our system shows that most DRDoS attacks are short lived, lasting only a few minutes, but that large-volume, long-lasting, and highly-distributed attacks against R&E networks are not uncommon. We then use the results of our analysis to discuss possible attack mitigation approaches that can be deployed at the IXP level, before the attack traffic overwhelms the victim's network bandwidth.


title: Frost filtered scale-invariant feature extraction and multilayer perceptron for hyperspectral image classification

abstract: Hyperspectral image (HSI) classification plays a significant in the field of remote sensing due to its ability to provide spatial and spectral information. Due to the rapid development and increasing of hyperspectral remote sensing technology, many methods have been developed for HSI classification but still a lack of achieving the better performance. A Frost Filtered Scale-Invariant Feature Transformation based MultiLayer Perceptron Classification (FFSIFT-MLPC) technique is introduced for classifying the hyperspectral image with higher accuracy and minimum time consumption. The FFSIFT-MLPC technique performs three major processes, namely preprocessing, feature extraction and classification using multiple layers. Initially, the hyperspectral image is divided into number of spectral bands. These bands are given as input in the input layer of perceptron. Then the Frost filter is used in FFSIFT-MLPC technique for preprocessing the input bands which helps to remove the noise from hyper-spectral image at the first hidden layer. After preprocessing task, texture, color and object features of hyper-spectral image are extracted at second hidden layer using Gaussian distributive scale-invariant feature transform. At the third hidden layer, Euclidean distance is measured between the extracted features and testing features. Finally, feature matching is carried out at the output layer for hyper-spectral image classification. The classified outputs are resulted in terms of spectral bands (i.e., different colors). Experimental analysis is performed with PSNR, classification accuracy, false positive rate and classification time with number of spectral bands. The results evident that presented FFSIFT-MLPC technique improves the hyperspectral image classification accuracy, PSNR and minimizes false positive rate as well as classification time than the state-of-the-art methods.


title: Just How Toxic is Data Poisoning? A Unified Benchmark for Backdoor and Data Poisoning Attacks

abstract: Data poisoning and backdoor attacks manipulate training data in order to cause models to fail during inference. A recent survey of industry practitioners found that data poisoning is the number one concern among threats ranging from model stealing to adversarial attacks. However, it remains unclear exactly how dangerous poisoning methods are and which ones are more effective considering that these methods, even ones with identical objectives, have not been tested in consistent or realistic settings. We observe that data poisoning and backdoor attacks are highly sensitive to variations in the testing setup. Moreover, we find that existing methods may not generalize to realistic settings. While these existing works serve as valuable prototypes for data poisoning, we apply rigorous tests to determine the extent to which we should fear them. In order to promote fair comparison in future work, we develop standardized benchmarks for data poisoning and backdoor attacks.


title: Control of microswimmers by spiral nematic vortices: transition from individual to collective motion and contraction, expansion, and stable circulation of bacterial swirls

abstract: Active systems comprised of self-propelled units show fascinating transitions from Brownian-like dynamics to collective coherent motion. Swirling of swimming bacteria is a spectacular example. This study demonstrates that a nematic liquid crystal environment patterned as a spiral vortex controls individual-to-collective transition in bacterial swirls and defines whether they expand or shrink. In dilute dispersions, the bacteria swim along open spiral trajectories, following the pre-imposed molecular orientation. The trajectories are nonpolar. As their concentration exceeds some threshold, the bacteria condense into unipolar circular swirls resembling stable limit cycles. This collective circular motion is controlled by the spiral angle that defines the splay-to-bend ratio of the background director. Vortices with dominating splay shrink the swirls towards the center, while vortices with dominating bend expand them to the periphery. 45o spiraling vortices with splay-bend parity produce the most stable swirls. All the dynamic scenarios are explained by hydrodynamic interactions of bacteria mediated by the patterned passive nematic environment and by the coupling between the concentration and orientation. The acquired knowledge of how to control individual and collective motion of microswimmers by a nematic environment can help in the development of microscopic mechanical systems.


title: Artificial Intelligence-Assisted Energy and Thermal Comfort Control for Sustainable Buildings: An Extended Representation of the Systematic Review

abstract: Different factors such as thermal comfort, humidity, air quality, and noise have significant combined effects on the acceptability and quality of the activities performed by the building occupants who spend most of their times indoors. Among the factors cited, thermal comfort, which contributes to the human well-being because of its connection with the thermoregulation of the human body. Therefore, the creation of thermally comfortable and energy efficient environments is of great importance in the design of the buildings and hence the heating, ventilation and air-conditioning systems. Recent works have been directed towards more advanced control strategies, based mainly on artificial intelligence which has the ability to imitate human behavior. This systematic literature review aims to provide an overview of the intelligent control strategies inside building and to investigate their ability to balance thermal comfort and energy efficiency optimization in indoor environments. Methods. A systematic literature review examined the peer-reviewed research works using ACM Digital Library, Scopus, Google Scholar, IEEE Xplore (IEOL), Web of Science, and Science Direct (SDOL), besides other sources from manual search. With the following string terms: thermal comfort, comfort temperature, preferred temperature, intelligent control, advanced control, artificial intelligence, computational intelligence, building, indoors, and built environment. Inclusion criteria were: English, studies monitoring, mainly, human thermal comfort in buildings and energy efficiency simultaneously based on control strategies using the intelligent approaches. Preferred Reporting Items for Systematic Reviews and Meta-Analysis guidelines were used. Initially, 1,077 articles were yielded, and 120 ultimately met inclusion criteria and were reviewed.


title: Multiverse: Dynamic VM Provisioning for Virtualized High Performance Computing Clusters

abstract: Traditionally, HPC workloads have been deployed in bare-metal clusters; but the advances in virtualization have led the pathway for these workloads to be deployed in virtualized clusters. However, HPC cluster administrators/providers still face challenges in terms of resource elasticity and virtual machine (VM) provisioning at large-scale, due to the lack of coordination between a traditional HPC scheduler and the VM hypervisor (resource management layer). This lack of interaction leads to low cluster utilization and job completion throughput. Furthermore, the VM provisioning delays directly impact the overall performance of jobs in the cluster. Hence, there is a need for effectively provisioning virtualized HPC clusters, which can best-utilize the physical hardware with minimal provisioning overheads. Towards this, we propose Multiverse, a VM provisioning framework, which can dynamically spawn VMs for incoming jobs in a virtualized HPC cluster, by integrating the HPC scheduler along with VM resource manager. We have implemented this framework on the Slurm} scheduler along with the vSphere VM resource manager. In order to reduce the VM provisioning overheads, we use instant cloning which shares both the disk and memory with the parent VM, when compared to full VM cloning which has to boot-up a new VM from scratch. Measurements with real-world HPC workloads demonstrate that, instant cloning is 2.5x faster than full cloning in terms of VM provisioning time. Further, it improves resource utilization by up to 40%, and cluster throughput by up to 1.5x, when compared to full clone for bursty job arrival scenarios.


title: Better approximation algorithms for maximum weight internal spanning trees in cubic graphs and claw-free graphs

abstract: Given a connected vertex-weighted graph $G$, the maximum weight internal spanning tree (MaxwIST) problem asks for a spanning tree of $G$ that maximizes the total weight of internal nodes. This problem is NP-hard and APX-hard, with the currently best known approximation factor $1/2$ (Chen et al., Algorithmica 2019). For the case of claw-free graphs, Chen et al. present an involved approximation algorithm with approximation factor $7/12$. They asked whether it is possible to improve these ratios, in particular for claw-free graphs and cubic graphs. We improve the approximation factors for the MaxwIST problem in cubic graphs and claw-free graphs. For cubic graphs we present an algorithm that computes a spanning tree whose total weight of internal vertices is at least $\frac{3}{4}-\frac{3}{n}$ times the total weight of all vertices, where $n$ is the number of vertices of $G$. This ratio is almost tight for large values of $n$. For claw-free graphs of degree at least three, we present an algorithm that computes a spanning tree whose total internal weight is at least $\frac{3}{5}-\frac{1}{n}$ times the total vertex weight. The degree constraint is necessary as this ratio may not be achievable if we allow vertices of degree less than three. With the above ratios, we immediately obtain better approximation algorithms with factors $\frac{3}{4}-\epsilon$ and $\frac{3}{5}-\epsilon$ for the MaxwIST problem in cubic graphs and claw-free graphs of degree at least three, for any $\epsilon>0$. In addition to improving the approximation factors, the new algorithms are relatively short compared to that of Chen et al.. The new algorithms are fairly simple, and employ a variant of the depth-first search algorithm that selects a relatively-large-weight vertex in every branching step. Moreover, the new algorithms take linear time while previous algorithms for similar problem instances are super-linear.


title: On the independent set sequence of a tree

abstract: Alavi, Malde, Schwenk and Erd\H{o}s asked whether the independent set sequence of every tree is unimodal. Here we make some observations about this question. We show that for the uniformly random (labelled) tree, asymptotically almost surely (a.a.s.) the initial approximately 49.5\% of the sequence is increasing while the terminal approximately 38.8\% is decreasing. Our approach uses the Matrix Tree Theorem, combined with computation. We also present a generalization of a result of Levit and Mandrescu, concerning the final one-third of the independent set sequence of a K\"onig-Egerv\'ary graph.


title: Spin relaxation benchmarks and individual qubit addressability for holes in quantum dots

abstract: We investigate hole spin relaxation in the single- and multi-hole regime in a 2x2 germanium quantum dot array. We use radiofrequency (rf) charge sensing and observe Pauli Spin-Blockade (PSB) for every second interdot transition up to the (1,5)-(0,6) anticrossing, consistent with a standard Fock-Darwin spectrum. We find spin relaxation times $T_1$ as high as 32 ms for a quantum dot with single-hole occupation and 1.2 ms for a quantum dot occupied by five-holes, setting benchmarks for spin relaxation times for hole quantum dots. Furthermore, we investigate the qubit addressability and sensitivity to electric fields by measuring the resonance frequency dependence of each qubit on gate voltages. We are able to tune the resonance frequency over a large range for both the single and multi-hole qubit. Simultaneously, we find that the resonance frequencies are only weakly dependent on neighbouring gates, and in particular the five-hole qubit resonance frequency is more than twenty times as sensitive to its corresponding plunger gate. The excellent individual qubit tunability and long spin relaxation times make holes in germanium promising for addressable and high-fidelity spin qubits in dense two-dimensional quantum dot arrays for large-scale quantum information.


title: A Somewhat Random Walk Through Nuclear and Particle Physics

abstract: These notes are an outgrowth of an advanced undergraduate course taught at the University of Maryland, College Park. They are intended as an introduction to various aspects of particle and nuclear physics with an emphasis on the role of symmetry. The basic philosophy is to introduce many of the fundamental ideas in nuclear and particle physics using relatively sophisticated mathematical tools -- but to do so in as a simplified a context to explain the underlying ideas. Thus, for example, the Higgs mechanism is discussed in terms of an Abelian Higgs model. The emphasis is largely, but not entirely theoretical in orientation. The goal is for readers to develop an understanding of many of the underlying issues in a relatively sophisticated way.


title: The dilogarithm and abelian Chern-Simons

abstract: We construct the (enhanced Rogers) dilogarithm function from the spin Chern-Simons invariant of C*-connections. This leads to geometric proofs of basic dilogarithm identities and a geometric context for other properties, such as the branching structure.


title: Astronomy Paper Seminar Participation Guide & Reading Walkthrough

abstract: Welcome to the wonderful world of scientific inquiry! On this journey you'll be reading many$\times 10^N$ papers in your discipline. Therefore, efficiency in digesting and relaying this information is paramount. In this guide, we'll review how you can participate in your local astronomy seminars. Participation takes many forms, from contributing a recently discovered article to the discussion of a published paper. In this guide, we'll begin by providing some suggested introductory activities for beginner scientists. Then we discuss how to locate papers and assimilate their results. Finally we conclude with a discussion on paper presentation and note storage. This guide is intended for an undergraduate and graduate student audience, and we encourage faculty to read and distribute this guide to students.


title: A Survey on Deep Learning for Localization and Mapping: Towards the Age of Spatial Machine Intelligence

abstract: Deep learning based localization and mapping has recently attracted significant attention. Instead of creating hand-designed algorithms through exploitation of physical models or geometric theories, deep learning based solutions provide an alternative to solve the problem in a data-driven way. Benefiting from ever-increasing volumes of data and computational power, these methods are fast evolving into a new area that offers accurate and robust systems to track motion and estimate scenes and their structure for real-world applications. In this work, we provide a comprehensive survey, and propose a new taxonomy for localization and mapping using deep learning. We also discuss the limitations of current models, and indicate possible future directions. A wide range of topics are covered, from learning odometry estimation, mapping, to global localization and simultaneous localization and mapping (SLAM). We revisit the problem of perceiving self-motion and scene understanding with on-board sensors, and show how to solve it by integrating these modules into a prospective spatial machine intelligence system (SMIS). It is our hope that this work can connect emerging works from robotics, computer vision and machine learning communities, and serve as a guide for future researchers to apply deep learning to tackle localization and mapping problems.


title: Properties of the Hot Ambient Medium of Early-type Galaxies Hosting Powerful Radio Sources

abstract: We present an archival analysis of Chandra X-ray observations for twelve nearby early-type galaxies hosting radio sources with radio power $>10^{23} \, \rm{W}~\rm{Hz}^{-1}$ at 1.4 GHz, similar to the radio power of the radio source in NGC 4261. Previously, in a similar analysis of eight nearby X-ray and optically-bright elliptical galaxies, Werner et al. 2012, found that NGC 4261 exhibited unusually low central gas entropy compared to the full sample. In the central 0.3 kpc of NGC 4261, the ratio of cooling time to freefall time ($t_{\rm{cool}}/t_{\rm ff}$) is less than $10$, indicating that cold clouds may be precipitating out of the hot ambient medium and providing fuel for accretion in the central region. NGC 4261 also hosts the most powerful radio source in the original sample. Because NGC 4261 may represent an important phase during which powerful feedback from a central active galactic nucleus (AGN) is fueled by multiphase condensation in the central kpc, we searched the Chandra archive for analogs to NGC 4261. We present entropy profiles of those galaxies as well as profiles of $t_{\rm{cool}}/t_{\rm ff}$. We find that one of them, IC 4296, exhibits properties similar to NGC 4261, including the presence of only single phase gas outside of $r \sim 2$ kpc and a similar central velocity dispersion. We compare the properties of NGC 4261 and IC 4296 to hydrodynamic simulations of AGN feedback fueled by precipitation. Over the course of those simulations, the single phase galaxy has an entropy gradient that remains similar to the entropy profiles inferred from our observations.


title: Collective radiation from distant emitters

abstract: Waveguides allow for direct coupling of emitters separated by large distances, offering a path to connect remote quantum systems. However, when facing the distances needed for practical applications, retardation effects due to the finite speed of light are often overlooked. Previous works studied the non-Markovian dynamics of emitters with retardation, but the properties of the radiated field remain mostly unexplored. By considering a toy model of two distant two-level atoms coupled through a waveguide, we observe that the spectrum of the radiated field exhibits non-Markovian features such as linewidth broadening beyond standard superradiance, or narrow Fano resonance-like peaks. We also show that the dipole-dipole interaction decays exponentially with distance as a result of retardation, with the range determined by the atomic linewidth. We discuss a proof-of-concept implementation of our results in a superconducting circuit platform.


title: Hybrid Low-Power Wide-Area Mesh Network for IoT Applications

abstract: The recent advancement of the Internet of Things (IoT) enables the possibility of data collection from diverse environments using IoT devices. However, despite the rapid advancement of low-power communication technologies, the deployment of IoT networks still faces many challenges. In this paper, we propose a hybrid, low-power, wide-area network (LPWAN) structure that can achieve wide-area communication coverage and low power consumption on IoT devices by utilizing both sub-GHz long-range radio and 2.4 GHz short-range radio. Specifically, we constructed a low-power mesh network with LoRa, a physical-layer standard that can provide long-range (kilometers) point-to-point communication using custom time-division multiple access (TDMA). Furthermore, we extended the capabilities of the mesh network by enabling ANT, an ultra-low-power, short-range communication protocol to satisfy data collection in dense device deployments. Third, we demonstrate the performance of the hybrid network with two real-world deployments at the Purdue University campus and at the university-owned farm. The results suggest that both networks have superior advantages in terms of cost, coverage, and power consumption vis-\`a-vis other IoT solutions, like LoRaWAN.


title: Linear instability of stationary solutions for the Korteweg-de Vries equation on a star graph

abstract: The aim of this work is to establish a linear instability criterium of stationary solutions for the Korteweg-de Vries model on a star graph with a structure represented by a finite collections of semi-infinite edges. By considering a boundary condition of $\delta$-type interaction at the graph-vertex, we show that the continuous tail and bump profiles are linearly unstable in a balanced star graph. The use of the analytic perturbation theory of operators and the extension theory of symmetric operators is a piece fundamental in our stability analysis. The arguments presented in this investigation has prospects for the study of the instability of stationary waves solutions of other nonlinear evolution equations on star graphs.


title: Opinion Diffusion Software with Strategic Opinion Revelation and Unfriending

abstract: We present a novel software suite for social network modeling and opinion diffusion processes. Much research on social network science has assumed networks with static topologies. More recently, attention has been turned to networks that evolve. Although software for modeling both the topological evolution of networks and diffusion processes are constantly improving, very little attention has been paid to agent modeling. Our software is designed to be robust, modular, and extensible, providing the ability to model dynamic social network topologies and multidimensional diffusion processes, different styles of agent including non-homophilic paradigms, as well as a testing environment for multi-agent reinforcement learning (MARL) experiments with diverse sets of agent types. We also illustrate the value of diverse agent modeling, and environments that allow for strategic unfriending. Our work shows that polarization and consensus dynamics, as well as topological clustering effects, may rely more than previously known on individuals' goals for the composition of their neighborhood's opinions.


title: A Causally Formulated Hazard Ratio Estimation through Backdoor Adjustment on Structural Causal Model

abstract: Identifying causal relationships for a treatment intervention is a fundamental problem in health sciences. Randomized controlled trials (RCTs) are considered the gold standard for identifying causal relationships. However, recent advancements in the theory of causal inference based on the foundations of structural causal models (SCMs) have allowed the identification of causal relationships from observational data, under certain assumptions. Survival analysis provides standard measures, such as the hazard ratio, to quantify the effects of an intervention. While hazard ratios are widely used in clinical and epidemiological studies for RCTs, a principled approach does not exist to compute hazard ratios for observational studies with SCMs. In this work, we review existing approaches to compute hazard ratios as well as their causal interpretation, if it exists. We also propose a novel approach to compute hazard ratios from observational studies using backdoor adjustment through SCMs and do-calculus. Finally, we evaluate the approach using experimental data for Ewing's sarcoma.


title: A dynamo amplifies the magnetic field of a Milky-Way-like galaxy

abstract: The magnetic fields of spiral galaxies are so strong that they cannot be primordial. Their typical values are over one billion times higher than any value predicted for the early Universe. Explaining this immense growth and incorporating it in galaxy evolution theories is one of the long-standing challenges in astrophysics. So far, the most successful theory for the sustained growth of the galactic magnetic field is the alpha-omega dynamo. This theory predicts a characteristic dipolar or quadrupolar morphology for the galactic magnetic field, which has been observed in external galaxies. However, so far, there has been no direct demonstration of a mean-field dynamo operating in direct, multi-physics simulations of spiral galaxies. We do so in this work. We employ numerical models of isolated, star-forming spiral galaxies that include a magnetized gaseous disk, a dark matter halo, stars, and stellar feedback. Naturally, the resulting magnetic field has a complex morphology that includes a strong random component. Using a smoothing of the magnetic field on small scales, we are able to separate the mean from the turbulent component and analyze them individually. We find that a mean-field dynamo naturally occurs as a result of the dynamical evolution of the galaxy and amplifies the magnetic field by an order of magnitude over half a Gyr. Despite the highly dynamical nature of these models, the morphology of the mean component of the field is identical to analytical predictions. This result underlines the importance of the mean-field dynamo in galactic evolution. Moreover, by demonstrating the natural growth of the magnetic field in a complex galactic environment, it brings us a step closer to understanding the cosmic origin of magnetic fields.


title: LAMP: Large Deep Nets with Automated Model Parallelism for Image Segmentation

abstract: Deep Learning (DL) models are becoming larger, because the increase in model size might offer significant accuracy gain. To enable the training of large deep networks, data parallelism and model parallelism are two well-known approaches for parallel training. However, data parallelism does not help reduce memory footprint per device. In this work, we introduce Large deep 3D ConvNets with Automated Model Parallelism (LAMP) and investigate the impact of both input's and deep 3D ConvNets' size on segmentation accuracy. Through automated model parallelism, it is feasible to train large deep 3D ConvNets with a large input patch, even the whole image. Extensive experiments demonstrate that, facilitated by the automated model parallelism, the segmentation accuracy can be improved through increasing model size and input context size, and large input yields significant inference speedup compared with sliding window of small patches in the inference. Code is available\footnote{https://monai.io/research/lamp-automated-model-parallelism}.


title: Graph Neural Networks and Reinforcement Learning for Behavior Generation in Semantic Environments

abstract: Most reinforcement learning approaches used in behavior generation utilize vectorial information as input. However, this requires the network to have a pre-defined input-size -- in semantic environments this means assuming the maximum number of vehicles. Additionally, this vectorial representation is not invariant to the order and number of vehicles. To mitigate the above-stated disadvantages, we propose combining graph neural networks with actor-critic reinforcement learning. As graph neural networks apply the same network to every vehicle and aggregate incoming edge information, they are invariant to the number and order of vehicles. This makes them ideal candidates to be used as networks in semantic environments -- environments consisting of objects lists. Graph neural networks exhibit some other advantages that make them favorable to be used in semantic environments. The relational information is explicitly given and does not have to be inferred. Moreover, graph neural networks propagate information through the network and can gather higher-degree information. We demonstrate our approach using a highway lane-change scenario and compare the performance of graph neural networks to conventional ones. We show that graph neural networks are capable of handling scenarios with a varying number and order of vehicles during training and application.


title: Unsettling physics in the quantum-corrected Schwarzschild black hole

abstract: We study a quantum-corrected Schwarzschild black hole proposed recently in Loop Quantum Gravity. Prompted by the fact that corrections to the innermost stable circular orbit of Schwarzschild diverge, we investigate timelike and null radial geodesics. Massive particles moving radially outwards are confined, while photons make it to infinity with infinite redshift. This unsettling physics, which deviates radically from both Schwarzschild (near the horizon) and Minkowski (at infinity) is due to repulsion by the negative quantum energy density that makes the quasilocal mass vanish as one approaches spatial infinity.


title: An information-theoretic account of semantic interference in word production

abstract: I present a computational-level model of semantic interference effects in word production. Word production is cast as a rate-distortion problem where an agent selects words to minimize a measure of cost while also minimizing the resources used to compute the output word based on perceptual input and behavioral goals. I show that similarity-based interference among words arises naturally in this setup, and I present a series of simulations showing that the model captures some of the key empirical patterns observed in Stroop and Picture-Word Interference paradigms. I argue that the rate-distortion account of interference provides a high-level formalization of computational principles that are instantiated more mechanistically in existing models.


title: Sequentially optimized projections in X-ray imaging

abstract: This work applies Bayesian experimental design to selecting optimal projection geometries in (discretized) parallel beam X-ray tomography assuming the prior and the additive noise are Gaussian. The introduced greedy exhaustive optimization algorithm proceeds sequentially, with the posterior distribution corresponding to the previous projections serving as the prior for determining the design parameters, i.e. the imaging angle and the lateral position of the source-receiver pair, for the next one. The algorithm allows redefining the region of interest after each projection as well as adapting parameters in the (original) prior to the measured data. Both A and D-optimality are considered, with emphasis on efficient evaluation of the corresponding objective functions. Two-dimensional numerical experiments demonstrate the functionality of the approach.


title: Empirical Measures, Geodesic Lengths, and a Variational Formula in First-Passage Percolation

abstract: This monograph resolves - in a dense class of cases - several open problems concerning geodesics in i.i.d. first-passage percolation on $\mathbb{Z}^d$. Our primary interest is in the empirical measures of edge-weights observed along geodesics from $0$ to $n\xi$, where $\xi$ is a fixed unit vector. For various dense families of edge-weight distributions, we prove that these empirical measures converge weakly to a deterministic limit as $n\to\infty$, answering a question of Hoffman. These families include arbitrarily small $L^\infty$-perturbations of any given distribution, almost every finitely supported distribution, uncountable collections of continuous distributions, and certain discrete distributions whose atoms can have any prescribed sequence of probabilities. Moreover, the constructions are explicit enough to guarantee examples possessing certain features, for instance: both continuous and discrete distributions whose support is all of $[0,\infty)$, and distributions given by a density function that is $k$-times differentiable. All results also hold for $\xi$-directed infinite geodesics. In comparison, we show that if $\mathbb{Z}^d$ is replaced by the infinite $d$-ary tree, then any distribution for the weights admits a unique limiting empirical measure along geodesics. In both the lattice and tree cases, our methodology is driven by a new variational formula for the time constant, which requires no assumptions on the edge-weight distribution. Incidentally, this variational approach also allows us to obtain new convergence results for geodesic lengths, which have been unimproved in the subcritical regime since the seminal 1965 manuscript of Hammersley and Welsh.


title: Density-based Stochastic Reachability Computation for Occupancy Prediction in Automated Driving

abstract: We propose a stochastic reachability computation framework for occupancy prediction in automated driving by directly solving the underlying transport partial differential equation governing the advection of the closed-loop joint density functions. The resulting nonparametric gridless computation is based on integration along the characteristic curves, and allows online computation of the time-varying collision probabilities. Numerical simulations highlight the scope the proposed method.


title: Laplacian Mixture Model Point Based Registration

abstract: Point base registration is an important part in many machine VISIOn applications, medical diagnostics, agricultural studies etc. The goal of point set registration is to find correspondences between different data sets and estimate the appropriate transformation that can map one set to another. Here we introduce a novel method for matching of different data sets based on Laplacian distribution. We consider the alignment of two point sets as probability density estimation problem. By using maximum likelihood methods we try to fit the Laplacian mixture model (LMM) centroids (source point set) to the data point set.


title: Exact Support Recovery in Federated Regression with One-shot Communication

abstract: Federated learning provides a framework to address the challenges of distributed computing, data ownership and privacy over a large number of distributed clients with low computational and communication capabilities. In this paper, we study the problem of learning the exact support of sparse linear regression in the federated learning setup. We provide a simple communication efficient algorithm which only needs one-shot communication with the centralized server to compute the exact support. Our method does not require the clients to solve any optimization problem and thus, can be run on devices with low computational capabilities. Our method is naturally robust to the problems of client failure, model poisoning and straggling clients. We formally prove that our method requires a number of samples per client that is polynomial with respect to the support size, but independent of the dimension of the problem. We require the number of distributed clients to be logarithmic in the dimension of the problem. If the predictor variables are mutually independent then the overall sample complexity matches the optimal sample complexity of the non-federated centralized setting. Furthermore, our method is easy to implement and has an overall polynomial time complexity.


title: An evolutionary study of volatile chemistry in protoplanetary disks

abstract: The volatile composition of a planet is determined by the inventory of gas and ice in the parent disk. The volatile chemistry in the disk is expected to evolve over time, though this evolution is poorly constrained observationally. We present ALMA observations of C18O, C2H, and the isotopologues H13CN, HC15N, and DCN towards five Class 0/I disk candidates. Combined with a sample of fourteen Class II disks presented in Bergner et al. (2019b), this data set offers a view of volatile chemical evolution over the disk lifetime. Our estimates of C18O abundances are consistent with a rapid depletion of CO in the first ~0.5-1 Myr of the disk lifetime. We do not see evidence that C2H and HCN formation are enhanced by CO depletion, possibly because the gas is already quite under-abundant in CO. Further CO depletion may actually hinder their production by limiting the gas-phase carbon supply. The embedded sources show several chemical differences compared to the Class II stage, which seem to arise from shielding of radiation by the envelope (impacting C2H formation and HC15N fractionation) and sublimation of ices from infalling material (impacting HCN and C18O abundances). Such chemical differences between Class 0/I and Class II sources may affect the volatile composition of planet-forming material at different stages in the disk lifetime.


title: Semantic Features Aided Multi-Scale Reconstruction of Inter-Modality Magnetic Resonance Images

abstract: Long acquisition time (AQT) due to series acquisition of multi-modality MR images (especially T2 weighted images (T2WI) with longer AQT), though beneficial for disease diagnosis, is practically undesirable. We propose a novel deep network based solution to reconstruct T2W images from T1W images (T1WI) using an encoder-decoder architecture. The proposed learning is aided with semantic features by using multi-channel input with intensity values and gradient of image in two orthogonal directions. A reconstruction module (RM) augmenting the network along with a domain adaptation module (DAM) which is an encoder-decoder model built-in with sharp bottleneck module (SBM) is trained via modular training. The proposed network significantly reduces the total AQT with negligible qualitative artifacts and quantitative loss (reconstructs one volume in approximately 1 second). The testing is done on publicly available dataset with real MR images, and the proposed network shows (approximately 1dB) increase in PSNR over SOTA.


title: Drive-Net: Convolutional Network for Driver Distraction Detection

abstract: To help prevent motor vehicle accidents, there has been significant interest in finding an automated method to recognize signs of driver distraction, such as talking to passengers, fixing hair and makeup, eating and drinking, and using a mobile phone. In this paper, we present an automated supervised learning method called Drive-Net for driver distraction detection. Drive-Net uses a combination of a convolutional neural network (CNN) and a random decision forest for classifying images of a driver. We compare the performance of our proposed Drive-Net to two other popular machine-learning approaches: a recurrent neural network (RNN), and a multi-layer perceptron (MLP). We test the methods on a publicly available database of images acquired under a controlled environment containing about 22425 images manually annotated by an expert. Results show that Drive-Net achieves a detection accuracy of 95%, which is 2% more than the best results obtained on the same database using other methods


title: PipeSim: Trace-driven Simulation of Large-Scale AI Operations Platforms

abstract: Operationalizing AI has become a major endeavor in both research and industry. Automated, operationalized pipelines that manage the AI application lifecycle will form a significant part of tomorrow's infrastructure workloads. To optimize operations of production-grade AI workflow platforms we can leverage existing scheduling approaches, yet it is challenging to fine-tune operational strategies that achieve application-specific cost-benefit tradeoffs while catering to the specific domain characteristics of machine learning (ML) models, such as accuracy, robustness, or fairness. We present a trace-driven simulation-based experimentation and analytics environment that allows researchers and engineers to devise and evaluate such operational strategies for large-scale AI workflow systems. Analytics data from a production-grade AI platform developed at IBM are used to build a comprehensive simulation model. Our simulation model describes the interaction between pipelines and system infrastructure, and how pipeline tasks affect different ML model metrics. We implement the model in a standalone, stochastic, discrete event simulator, and provide a toolkit for running experiments. Synthetic traces are made available for ad-hoc exploration as well as statistical analysis of experiments to test and examine pipeline scheduling, cluster resource allocation, and similar operational mechanisms.


title: The dynamics of a driven harmonic oscillator coupled to independent Ising spins in random fields

abstract: We aim at an understanding of the dynamical properties of a periodically driven damped harmonic oscillator coupled to a Random Field Ising Model (RFIM) at zero temperature, which is capable to show complex hysteresis. The system is a combination of a continuous (harmonic oscillator) and a discrete (RFIM) subsystem, which classifies it as a hybrid system. In this paper we focus on the hybrid nature of the system and consider only independent spins in quenched random local fields, which can already lead to complex dynamics such as chaos and multistability. We study the dynamic behavior of this system by using the theory of piecewise-smooth dynamical systems and discontinuity mappings. Specifically, we present bifurcation diagrams, Lyapunov exponents as well as results for the shape and the dimensions of the attractors and the self-averaging behavior of the attractor dimensions and the magnetization. Furthermore we investigate the dynamical behavior of the system for an increasing number of spins and the transition to the thermodynamic limit, where the system behaves like a driven harmonic oscillator with an additional nonlinear smooth external force.


title: Distributional Individual Fairness in Clustering

abstract: In this paper, we initiate the study of fair clustering that ensures distributional similarity among similar individuals. In response to improving fairness in machine learning, recent papers have investigated fairness in clustering algorithms and have focused on the paradigm of statistical parity/group fairness. These efforts attempt to minimize bias against some protected groups in the population. However, to the best of our knowledge, the alternative viewpoint of individual fairness, introduced by Dwork et al. (ITCS 2012) in the context of classification, has not been considered for clustering so far. Similar to Dwork et al., we adopt the individual fairness notion which mandates that similar individuals should be treated similarly for clustering problems. We use the notion of $f$-divergence as a measure of statistical similarity that significantly generalizes the ones used by Dwork et al. We introduce a framework for assigning individuals, embedded in a metric space, to probability distributions over a bounded number of cluster centers. The objective is to ensure (a) low cost of clustering in expectation and (b) individuals that are close to each other in a given fairness space are mapped to statistically similar distributions. We provide an algorithm for clustering with $p$-norm objective ($k$-center, $k$-means are special cases) and individual fairness constraints with provable approximation guarantee. We extend this framework to include both group fairness and individual fairness inside the protected groups. Finally, we observe conditions under which individual fairness implies group fairness. We present extensive experimental evidence that justifies the effectiveness of our approach.


title: C-SURE: Shrinkage Estimator and Prototype Classifier for Complex-Valued Deep Learning

abstract: The James-Stein (JS) shrinkage estimator is a biased estimator that captures the mean of Gaussian random vectors.While it has a desirable statistical property of dominance over the maximum likelihood estimator (MLE) in terms of mean squared error (MSE), not much progress has been made on extending the estimator onto manifold-valued data. We propose C-SURE, a novel Stein's unbiased risk estimate (SURE) of the JS estimator on the manifold of complex-valued data with a theoretically proven optimum over MLE. Adapting the architecture of the complex-valued SurReal classifier, we further incorporate C-SURE into a prototype convolutional neural network (CNN) classifier. We compare C-SURE with SurReal and a real-valued baseline on complex-valued MSTAR and RadioML datasets. C-SURE is more accurate and robust than SurReal, and the shrinkage estimator is always better than MLE for the same prototype classifier. Like SurReal, C-SURE is much smaller, outperforming the real-valued baseline on MSTAR (RadioML) with less than 1 percent (3 percent) of the baseline size


title: A Survey of $q$-Whittaker polynomials

abstract: Exploiting the fact that the $q$-Whittaker polynomials arise as a specialization of the (modified) Macdonald polynomials, we derive some of their basic properties, and explore interesting identities that they satisfy. We also show how they arise as graded Frobenius characteristics of $\S_n$-modules, and give a combinatorial approach to associated Pieri formulas.


title: An Efficient Smoothing Proximal Gradient Algorithm for Convex Clustering

abstract: Cluster analysis organizes data into sensible groupings and is one of fundamental modes of understanding and learning. The widely used K-means and hierarchical clustering methods can be dramatically suboptimal due to local minima. Recently introduced convex clustering approach formulates clustering as a convex optimization problem and ensures a globally optimal solution. However, the state-of-the-art convex clustering algorithms, based on the alternating direction method of multipliers (ADMM) or the alternating minimization algorithm (AMA), require large computation and memory space, which limits their applications. In this paper, we develop a very efficient smoothing proximal gradient algorithm (Sproga) for convex clustering. Our Sproga is faster than ADMM- or AMA-based convex clustering algorithms by one to two orders of magnitude. The memory space required by Sproga is less than that required by ADMM and AMA by at least one order of magnitude. Computer simulations and real data analysis show that Sproga outperforms several well known clustering algorithms including K-means and hierarchical clustering. The efficiency and superior performance of our algorithm will help convex clustering to find its wide application.


title: Powerfully solvable and powerfully simple groups

abstract: We introduce the notion of a powerfully solvable group. These are powerful groups possessing an abelian series of a special kind. These groups include in particular the class of powerfully nilpotent groups. We will also see that for a certain rich class of powerful groups we can naturally introduce the term powerfully simple group and prove a Jordan-H\"older type theorem that justifies the term.


title: Articulatory-WaveNet: Autoregressive Model For Acoustic-to-Articulatory Inversion

abstract: This paper presents Articulatory-WaveNet, a new approach for acoustic-to-articulator inversion. The proposed system uses the WaveNet speech synthesis architecture, with dilated causal convolutional layers using previous values of the predicted articulatory trajectories conditioned on acoustic features. The system was trained and evaluated on the ElectroMagnetic Articulography corpus of Mandarin Accented English (EMA-MAE),consisting of 39 speakers including both native English speakers and native Mandarin speakers speaking English. Results show significant improvement in both correlation and RMSE between the generated and true articulatory trajectories for the new method, with an average correlation of 0.83, representing a 36% relative improvement over the 0.61 correlation obtained with a baseline Hidden Markov Model (HMM)-Gaussian Mixture Model (GMM) inversion framework. To the best of our knowledge, this paper presents the first application of a point-by-point waveform synthesis approach to the problem of acoustic-to-articulatory inversion and the results show improved performance compared to previous methods for speaker dependent acoustic to articulatory inversion.


title: Locally trimmed least squares: conventional inference in possibly nonstationary models

abstract: A novel IV estimation method, that we term Locally Trimmed LS (LTLS), is developed which yields estimators with (mixed) Gaussian limit distributions in situations where the data may be weakly or strongly persistent. In particular, we allow for nonlinear predictive type of regressions where the regressor can be stationary short/long memory as well as nonstationary long memory process or a nearly integrated array. The resultant t-tests have conventional limit distributions (i.e. N(0; 1)) free of (near to unity and long memory) nuisance parameters. In the case where the regressor is a fractional process, no preliminary estimator for the memory parameter is required. Therefore, the practitioner can conduct inference while being agnostic about the exact dependence structure in the data. The LTLS estimator is obtained by applying certain chronological trimming to the OLS instrument via the utilisation of appropriate kernel functions of time trend variables. The finite sample performance of LTLS based t-tests is investigated with the aid of a simulation experiment. An empirical application to the predictability of stock returns is also provided.


title: Improved Asymptotic Expressions for the Eigenvalues of Laplace's Tidal Equations

abstract: Laplace's tidal equations govern the angular dependence of oscillations in stars when uniform rotation is treated within the so-called traditional approximation. Using a perturbation expansion approach, I derive improved expressions for the eigenvalue associated with these equations, valid in the asymptotic limit of large spin parameter $q$. These expressions have a relative accuracy of order $q^{-3}$ for gravito-inertial modes, and $q^{-1}$ for Rossby and Kelvin modes; the corresponding absolute accuracy is of order $q^{-1}$ for all three mode types. I validate my analysis against numerical calculations, and demonstrate how it can be applied to derive formulae for the periods and eigenfunctions of Rossby modes.


title: Simulating solvation and acidity in complex mixtures with first-principles accuracy: the case of CH$_3$SO$_3$H and H$_2$O$_2$ in phenol

abstract: We present a generally-applicable computational framework for the efficient and accurate characterization of molecular structural patterns and acid properties in explicit solvent using H$_2$O$_2$ and CH$_3$SO$_3$H in phenol as an example. In order to address the challenges posed by the complexity of the problem, we resort to a set of data-driven methods and enhanced sampling algorithms. The synergistic application of these techniques makes the first-principle estimation of the chemical properties feasible without renouncing to the use of explicit solvation, involving extensive statistical sampling. Ensembles of neural network potentials are trained on a set of configurations carefully selected out of preliminary simulations performed at a low-cost density-functional tight-binding (DFTB) level. The energy and forces of these configurations are then recomputed at the hybrid density functional theory (DFT) level and used to train the neural networks. The stability of the NN model is enhanced by using DFTB energetics as a baseline, but the efficiency of the direct NN (i.e., baseline-free) is exploited via a multiple-time step integrator. The neural network potentials are combined with enhanced sampling techniques, such as replica exchange and metadynamics, and used to characterize the relevant protonated species and dominant non-covalent interactions in the mixture, also considering nuclear quantum effects.


title: Meta Learning for Support Recovery in High-dimensional Precision Matrix Estimation

abstract: In this paper, we study meta learning for support (i.e., the set of non-zero entries) recovery in high-dimensional precision matrix estimation where we reduce the sufficient sample complexity in a novel task with the information learned from other auxiliary tasks. In our setup, each task has a different random true precision matrix, each with a possibly different support. We assume that the union of the supports of all the true precision matrices (i.e., the true support union) is small in size. We propose to pool all the samples from different tasks, and \emph{improperly} estimate a single precision matrix by minimizing the $\ell_1$-regularized log-determinant Bregman divergence. We show that with high probability, the support of the \emph{improperly} estimated single precision matrix is equal to the true support union, provided a sufficient number of samples per task $n \in O((\log N)/K)$, for $N$-dimensional vectors and $K$ tasks. That is, one requires less samples per task when more tasks are available. We prove a matching information-theoretic lower bound for the necessary number of samples, which is $n \in \Omega((\log N)/K)$, and thus, our algorithm is minimax optimal. Then for the novel task, we prove that the minimization of the $\ell_1$-regularized log-determinant Bregman divergence with the additional constraint that the support is a subset of the estimated support union could reduce the sufficient sample complexity of successful support recovery to $O(\log(|S_{\text{off}}|))$ where $|S_{\text{off}}|$ is the number of off-diagonal elements in the support union and is much less than $N$ for sparse matrices. We also prove a matching information-theoretic lower bound of $\Omega(\log(|S_{\text{off}}|))$ for the necessary number of samples. Synthetic experiments validate our theory.


title: Phase space learning with neural networks

abstract: This work proposes an autoencoder neural network as a non-linear generalization of projection-based methods for solving Partial Differential Equations (PDEs). The proposed deep learning architecture presented is capable of generating the dynamics of PDEs by integrating them completely in a very reduced latent space without intermediate reconstructions, to then decode the latent solution back to the original space. The learned latent trajectories are represented and their physical plausibility is analyzed. It is shown the reliability of properly regularized neural networks to learn the global characteristics of a dynamical system's phase space from the sample data of a single path, as well as its ability to predict unseen bifurcations.


title: Blow-up for wave equation with the scale-invariant damping and combined nonlinearities

abstract: In this article, we study the blow-up of the damped wave equation in the \textit{scale-invariant case} and in the presence of two nonlinearities. More precisely, we consider the following equation: $$u_{tt}-\Delta u+\frac{\mu}{1+t}u_t=|u_t|^p+|u|^q, \quad \mbox{in}\ \R^N\times[0,\infty), $$ with small initial data.\\ For $\mu < \frac{N(q-1)}{2}$ and $\mu \in (0, \mu_*)$, where $\mu_*>0$ is depending on the nonlinearties' powers and the space dimension ($\mu_*$ satisfies $(q-1)\left((N+2\mu_*-1)p-2\right) = 4$), we prove that the wave equation, in this case, behaves like the one without dissipation ($\mu =0$). Our result completes the previous studies in the case where the dissipation is given by $\frac{\mu}{(1+t)^\beta}u_t; \ \beta >1$ (\cite{LT3}), where, contrary to what we obtain in the present work, the effect of the damping is not significant in the dynamics. Interestingly, in our case, the influence of the damping term $\frac{\mu}{1+t}u_t$ is important.


title: Petz reconstruction in random tensor networks

abstract: We illustrate the ideas of bulk reconstruction in the context of random tensor network toy models of holography. Specifically, we demonstrate how the Petz reconstruction map works to obtain bulk operators from the boundary data by exploiting the replica trick. We also take the opportunity to comment on the differences between coarse-graining and random projections.


title: Analogues of Katona's and Milner's Theorems for two families

abstract: Let $n>s>0$ be integers, $X$ an $n$-element set and $\mathscr{A}, \mathscr{B}\subset 2^X$ two families. If $|A\cup B|\le s$ for all $A\in\mathscr{A}, B\in \mathscr{B}$, then $\mathscr{A}$ and $\mathscr{B}$ are called cross $s$-union. Assuming that neither $\mathscr{A}$ nor $\mathscr{B}$ is empty, we prove several best possible bounds. In particular, we show that $|\mathscr{A}|+|\mathscr{B}|\le 1+\sum\limits_{0\le i\le s}{{n}\choose{i}}$. Supposing $n\ge 2s$ and $\mathscr{A},\mathscr{B}$ are antichains, we show that $|\mathscr{A}|+|\mathscr{B}|\le {{n}\choose{1}}+{{n}\choose{s-1}}$ unless $\mathscr{A}=\{\emptyset\}$ or $\mathscr{B}=\{\emptyset\}$. An analogous result for three families is established as well.


title: Measures of space-time non-separability of electromagnetic pulses

abstract: Electromagnetic pulses are typically treated as space-time (or space-frequency) separable solutions of Maxwell's equations, where spatial and temporal (spectral) dependence can be treated separately. In contrast to this traditional viewpoint, recent advances in structured light and topological optics have highlighted the non-trivial wave-matter interactions of pulses with complex topology and space-time non-separable structure, as well as their potential for energy and information transfer. A characteristic example of such a pulse is the "Flying Doughnut" (FD), a space-time non-separable toroidal few-cycle pulse with links to toroidal and non-radiating (anapole) excitations in matter. Here, we propose a quantum-mechanics-inspired methodology for the characterization of space-time non-separability in structured pulses. In analogy to the non-separability of entangled quantum systems, we introduce the concept of space-spectrum entangled states to describe the space-time non-separability of classical electromagnetic pulses and develop a method to reconstruct the corresponding density matrix by state tomography. We apply our method to the FD pulse and obtain the corresponding fidelity, concurrence, and entanglement of formation. We demonstrate that such properties dug out from quantum mechanics quantitatively characterize the evolution of the general spatiotemporal structured pulse upon propagation.


title: Coherent WaveBurst, a pipeline for unmodeled gravitational-wave data analysis

abstract: coherent WaveBurst (cWB) is a highly configurable pipeline designed to detect a broad range of gravitational-wave (GW) transients in the data of the worldwide network of GW detectors. The algorithmic core of cWB is a time-frequency analysis with the Wilson-Daubechies-Meyer wavelets aimed at the identification of GW events without prior knowledge of the signal waveform. cWB has been in active development since 2003 and it has been used to analyze all scientific data collected by the LIGO-Virgo detectors ever since. On September 14, 2015, the cWB low-latency search detected the first gravitational-wave event, GW150914, a merger of two black holes. In 2019, a public open-source version of cWB has been released with GPLv3 license.


title: Fuzzy field theories and related matrix models

abstract: We review the description of scalar field theories on fuzzy spaces by Hermitian random matrix models. After reminding the reader of the relevant aspects of the random matrix theory and construction of the fuzzy spaces, we summarize the most important results for the scalar fields on such spaces. We then introduce the multi-trace matrix models relevant for the analytical description of scalar field theories on fuzzy spaces and show to what extent they do, and to what extent they do not, describe the know phase structure of $\phi^4$ theory on the fuzzy sphere.


title: Heavy-flavour correlations and jets with ALICE at the LHC

abstract: This contribution summarises the results on heavy-flavour correlations and jets measured with ALICE detector. Studies of D$^{0}$-tagged jets are presented in various collision systems and energies. The measurement of the fraction of jet momentum carried by the D$^{0}$ meson in the direction of the jet axis ($z_{\rm ||}^{\rm ch,jet}$) in pp collisions at $\sqrt{s}$ = 5.02 and 13 TeV as well as the nuclear modification factor in p--Pb and Pb--Pb collisions are presented. The measurements of the production of jets tagged by electrons from heavy-flavour decays in pp and p--Pb are compared to expectations from POWHEG+PYTHIA8 simulations. The $p_{\rm T}$-differential cross section of b-tagged jets in p--Pb collisions is also reported and compared to expectations from POWHEG event generator. The azimuthal correlation of D mesons with charged particles in pp and p--Pb collisions are also presented.


title: Pore-Scale Transport and Two-Phase Fluid Structures in Fibrous Porous Layers: Application to Fuel Cells and Beyond

abstract: We present pore-scale simulations of two-phase flows in a reconstructed fibrous porous layer. The three dimensional microstructure of the material, a fuel cell gas diffusion layer, is acquired via X-ray computed tomography and used as input for lattice Boltzmann simulations. We perform a quantitative analysis of the multiphase pore-scale dynamics and we identify the dominant fluid structures governing mass transport. The results show the existence of three different regimes of transport: a fast inertial dynamics at short times, characterised by a compact uniform front, a viscous-capillary regime at intermediate times, where liquid is transported along a gradually increasing number of preferential flow paths of the size of one-two pores, and a third regime at longer times, where liquid, after having reached the outlet, is exclusively flowing along such flow paths and the two-phase fluid structures are stabilised. We observe that the fibrous layer presents significant variations in its microscopic morphology, which have an important effect on the pore invasion dynamics, and counteract the stabilising viscous force. Liquid transport is indeed affected by the presence of microstructure-induced capillary pressures acting adversely to the flow, leading to capillary fingering transport mechanism and unstable front displacement, even in the absence of hydrophobic treatments of the porous material. We propose a macroscopic model based on an effective contact angle that mimics the effects of the such a dynamic capillary pressure. Finally, we underline the significance of the results for the optimal design of face masks in an effort to mitigate the current COVID-19 pandemic.


title: Similarity Search with Tensor Core Units

abstract: Tensor Core Units (TCUs) are hardware accelerators developed for deep neural networks, which efficiently support the multiplication of two dense $\sqrt{m}\times \sqrt{m}$ matrices, where $m$ is a given hardware parameter. In this paper, we show that TCUs can speed up similarity search problems as well. We propose algorithms for the Johnson-Lindenstrauss dimensionality reduction and for similarity join that, by leveraging TCUs, achieve a $\sqrt{m}$ speedup up with respect to traditional approaches.


title: Unitarizability in Corank Three for Classical p-adic Groups

abstract: Let G be the F-points of a classical group defined over a p-adic field F of characteristic 0. We classify the irreducible unitarizable representation of G that are subquotients of the parabolic induction of cuspidal representations of Levi subgroup of corank at most 3 in G.


title: Dual Topological Insulator Device with Disorder Robustness

abstract: Two-dimensional Na$_3$Bi is a dual topological insulator protected by time-reversal and mirror symmetry, resulting in a promising platform for devices design. However, in reality, the design of topological devices is hindered by a sensitivity against disorder and temperature. We study the topological properties of Na$_3$Bi in the presence of intrinsic defects, investigating the robustness of the edge states and the resulting transport properties. We apply a recursive Green's function technique enabling the study of disordered systems with lengths comparable to experimentally synthesized materials, in the order of micrometers. We combine our findings to propose a topological insulator device, where intrinsic defects are used to filter the response of trivial bulk states. This results in a stable conductance throughout a large range of electronic temperatures, and controllable by a perpendicular electric field. Our proposal is general, enabling the design of various dual topological insulators devices.


title: GW190814: Gravitational Waves from the Coalescence of a 23 M$_\odot$ Black Hole with a 2.6 M$_\odot$ Compact Object

abstract: We report the observation of a compact binary coalescence involving a 22.2 - 24.3 $M_{\odot}$ black hole and a compact object with a mass of 2.50 - 2.67 $M_{\odot}$ (all measurements quoted at the 90$\%$ credible level). The gravitational-wave signal, GW190814, was observed during LIGO's and Virgo's third observing run on August 14, 2019 at 21:10:39 UTC and has a signal-to-noise ratio of 25 in the three-detector network. The source was localized to 18.5 deg$^2$ at a distance of $241^{+41}_{-45}$ Mpc; no electromagnetic counterpart has been confirmed to date. The source has the most unequal mass ratio yet measured with gravitational waves, $0.112^{+0.008}_{-0.009}$, and its secondary component is either the lightest black hole or the heaviest neutron star ever discovered in a double compact-object system. The dimensionless spin of the primary black hole is tightly constrained to $\leq 0.07$. Tests of general relativity reveal no measurable deviations from the theory, and its prediction of higher-multipole emission is confirmed at high confidence. We estimate a merger rate density of 1-23 Gpc$^{-3}$ yr$^{-1}$ for the new class of binary coalescence sources that GW190814 represents. Astrophysical models predict that binaries with mass ratios similar to this event can form through several channels, but are unlikely to have formed in globular clusters. However, the combination of mass ratio, component masses, and the inferred merger rate for this event challenges all current models for the formation and mass distribution of compact-object binaries.


title: In-situ Nanoscale Characterization of Composition and Structure during Formation of Ultrathin Nickel Silicide

abstract: We characterize composition and structure of ultrathin nickel silicide during formation from 3 nm Ni films on Si(100) using in-situ high resolution ion scattering and high resolution transmission electron microscopy. We show the transition to occur in discrete steps, in which an intermediate phase is observed within a narrow range of temperature from 230 oC to 290 oC. The film composition of this intermediate phase is found to be 50% Ni:50% Si, without evidence for long-range structure, indicating the film to be a homogeneous monosilicide NiSi phase. The final phase is resemblant of the cubic disilicide NiSi2, but with slightly off-stoichiometric composition of 38% Ni and 62% Si. Along the [100] axis, the lattices of the film and the substrate are found in perfect alignment. Due to the epitaxial growth of the silicide, a contraction of the c lattice constant of the film by 0.7-1% is detected.


title: Identification of single spectral lines through supervised machine learning in a large HST survey (WISP): a pilot study for Euclid and WFIRST

abstract: Future surveys focusing on understanding the nature of dark energy (e.g., Euclid and WFIRST) will cover large fractions of the extragalactic sky in near-IR slitless spectroscopy. These surveys will detect a large number of galaxies that will have only one emission line in the covered spectral range. In order to maximize the scientific return of these missions, it is imperative that single emission lines are correctly identified. Using a supervised machine-learning approach, we classified a sample of single emission lines extracted from the WFC3 IR Spectroscopic Parallel survey (WISP), one of the closest existing analogs to future slitless surveys. Our automatic software integrates a SED fitting strategy with additional independent sources of information. We calibrated it and tested it on a "gold" sample of securely identified objects with multiple lines detected. The algorithm correctly classifies real emission lines with an accuracy of 82.6%, whereas the accuracy of the SED fitting technique alone is low (~50%) due to the limited amount of photometric data available (<=6 bands). While not specifically designed for the Euclid and WFIRST surveys, the algorithm represents an important precursor of similar algorithms to be used in these future missions.


title: A Critical Look at the Electroweak Phase Transition

abstract: The electroweak phase transition broke the electroweak symmetry. Perturbative methods used to calculate observables related to this phase transition suffer from severe problems such as gauge dependence, infrared divergences, and a breakdown of perturbation theory. In this paper we develop robust perturbative tools for dealing with phase transitions. We argue that gauge and infrared problems are absent in a consistent power-counting. We calculate the finite temperature effective potential to two loops for general gauge-fixing parameters in a generic model. We demonstrate gauge invariance, and perform numerical calculations for the Standard Model in Fermi gauge.


title: HVDC Surface Flashover in Compressed Air for Various Dielectrics

abstract: This study measures the voltage at which flashover occurs in compressed air for a variety of dielectric materials and lengths in a uniform field for DC voltages up to 100 kV. Statistical time lag is recorded and characterized, displaying a roughly exponential dependence on breakdown voltage. Of the materials tested, acrylic is observed to be the most resistant to flashover. These data are intended to facilitate the design of compressed-air insulated high voltage systems as an alternative to SF6 insulated systems.


title: Prediction error-driven memory consolidation for continual learning. On the case of adaptive greenhouse models

abstract: This work presents an adaptive architecture that performs online learning and faces catastrophic forgetting issues by means of episodic memories and prediction-error driven memory consolidation. In line with evidences from the cognitive science and neuroscience, memories are retained depending on their congruency with the prior knowledge stored in the system. This is estimated in terms of prediction error resulting from a generative model. Moreover, this AI system is transferred onto an innovative application in the horticulture industry: the learning and transfer of greenhouse models. This work presents a model trained on data recorded from research facilities and transferred to a production greenhouse.


title: Adaptive County Level COVID-19 Forecast Models: Analysis and Improvement

abstract: Accurately forecasting county level COVID-19 confirmed cases is crucial to optimizing medical resources. Forecasting emerging outbreaks pose a particular challenge because many existing forecasting techniques learn from historical seasons trends. Recurrent neural networks (RNNs) with LSTM-based cells are a logical choice of model due to their ability to learn temporal dynamics. In this paper, we adapt the state and county level influenza model, TDEFSI-LONLY, proposed in Wang et a. [l2020] to national and county level COVID-19 data. We show that this model poorly forecasts the current pandemic. We analyze the two week ahead forecasting capabilities of the TDEFSI-LONLY model with combinations of regularization techniques. Effective training of the TDEFSI-LONLY model requires data augmentation, to overcome this challenge we utilize an SEIR model and present an inter-county mixing extension to this model to simulate sufficient training data. Further, we propose an alternate forecast model, {\it County Level Epidemiological Inference Recurrent Network} (\alg{}) that trains an LSTM backbone on national confirmed cases to learn a low dimensional time pattern and utilizes a time distributed dense layer to learn individual county confirmed case changes each day for a two weeks forecast. We show that the best, worst, and median state forecasts made using CLEIR-Net model are respectively New York, South Carolina, and Montana.


title: A Bayesian incorporated linear non-Gaussian acyclic model for multiple directed graph estimation to study brain emotion circuit development in adolescence

abstract: Emotion perception is essential to affective and cognitive development which involves distributed brain circuits. The ability of emotion identification begins in infancy and continues to develop throughout childhood and adolescence. Understanding the development of brain's emotion circuitry may help us explain the emotional changes observed during adolescence. Our previous study delineated the trajectory of brain functional connectivity (FC) from late childhood to early adulthood during emotion identification tasks. In this work, we endeavour to deepen our understanding from association to causation. We proposed a Bayesian incorporated linear non-Gaussian acyclic model (BiLiNGAM), which incorporated our previous association model into the prior estimation pipeline. In particular, it can jointly estimate multiple directed acyclic graphs (DAGs) for multiple age groups at different developmental stages. Simulation results indicated more stable and accurate performance over various settings, especially when the sample size was small (high-dimensional cases). We then applied to the analysis of real data from the Philadelphia Neurodevelopmental Cohort (PNC). This included 855 individuals aged 8-22 years who were divided into five different adolescent stages. Our network analysis revealed the development of emotion-related intra- and inter- modular connectivity and pinpointed several emotion-related hubs. We further categorized the hubs into two types: in-hubs and out-hubs, as the center of receiving and distributing information. Several unique developmental hub structures and group-specific patterns were also discovered. Our findings help provide a causal understanding of emotion development in the human brain.


title: Mathematical modeling and prediction of COVID-19 in Moscow city and Novosibirsk region

abstract: The paper formulates and solves the problem of identification of unknown parameters of mathematical models of the spread of COVID-19 coronavirus infection, based on SEIR type models, based on additional information about the number of detected cases, mortality, self-isolation coefficient and tests performed for the Moscow city and the Novosibirsk Region from 03.23.2020. Within the framework of the models used, the population is divided into seven (SEIR-HCD) and five (SEIR-D) groups with similar characteristics with transition probabilities between groups depending on a specific region. Identifiability analysis of the SEIR-HCD mathematical model was carried out, which revealed the least sensitive unknown parameters to additional measurements. The tasks of refining the parameters are reduced to minimizing the corresponding target functionals, which were solved using stochastic methods (simulating annealing, differential evolution, genetic algorithm, etc.). For a different amount of tested data, a prognostic scenario for the development of the disease in the city of Moscow and the Novosibirsk region was developed, the peak is predicted the development of the epidemic in Moscow with an error of 2 days and 174 detected cases, and an analysis of the applicability of the developed models was carried out.


title: A maximum-entropy approach to off-policy evaluation in average-reward MDPs

abstract: This work focuses on off-policy evaluation (OPE) with function approximation in infinite-horizon undiscounted Markov decision processes (MDPs). For MDPs that are ergodic and linear (i.e. where rewards and dynamics are linear in some known features), we provide the first finite-sample OPE error bound, extending existing results beyond the episodic and discounted cases. In a more general setting, when the feature dynamics are approximately linear and for arbitrary rewards, we propose a new approach for estimating stationary distributions with function approximation. We formulate this problem as finding the maximum-entropy distribution subject to matching feature expectations under empirical dynamics. We show that this results in an exponential-family distribution whose sufficient statistics are the features, paralleling maximum-entropy approaches in supervised learning. We demonstrate the effectiveness of the proposed OPE approaches in multiple environments.


title: Fairness Through Robustness: Investigating Robustness Disparity in Deep Learning

abstract: Deep neural networks (DNNs) are increasingly used in real-world applications (e.g. facial recognition). This has resulted in concerns about the fairness of decisions made by these models. Various notions and measures of fairness have been proposed to ensure that a decision-making system does not disproportionately harm (or benefit) particular subgroups of the population. In this paper, we argue that traditional notions of fairness that are only based on models' outputs are not sufficient when the model is vulnerable to adversarial attacks. We argue that in some cases, it may be easier for an attacker to target a particular subgroup, resulting in a form of \textit{robustness bias}. We show that measuring robustness bias is a challenging task for DNNs and propose two methods to measure this form of bias. We then conduct an empirical study on state-of-the-art neural networks on commonly used real-world datasets such as CIFAR-10, CIFAR-100, Adience, and UTKFace and show that in almost all cases there are subgroups (in some cases based on sensitive attributes like race, gender, etc) which are less robust and are thus at a disadvantage. We argue that this kind of bias arises due to both the data distribution and the highly complex nature of the learned decision boundary in the case of DNNs, thus making mitigation of such biases a non-trivial task. Our results show that robustness bias is an important criterion to consider while auditing real-world systems that rely on DNNs for decision making. Code to reproduce all our results can be found here: \url{https://github.com/nvedant07/Fairness-Through-Robustness}


title: WD3: Taming the Estimation Bias in Deep Reinforcement Learning

abstract: The overestimation phenomenon caused by function approximation is a well-known issue in value-based reinforcement learning algorithms such as deep Q-networks and DDPG, which could lead to suboptimal policies. To address this issue, TD3 takes the minimum value between a pair of critics. In this paper, we show that the TD3 algorithm introduces underestimation bias in mild assumptions. To obtain a more precise estimation for value function, we unify these two opposites and propose a novel algorithm \underline{W}eighted \underline{D}elayed \underline{D}eep \underline{D}eterministic Policy Gradient (WD3), which can eliminate the estimation bias and further improve the performance by weighting a pair of critics. To demonstrate the effectiveness of WD3, we compare the learning process of value function between DDPG, TD3, and WD3. The results verify that our algorithm does eliminate the estimation error of value functions. Furthermore, we evaluate our algorithm on the continuous control tasks. We observe that in each test task, the performance of WD3 consistently outperforms, or at the very least matches, that of the state-of-the-art algorithms\footnote{Our code is available at~\href{https://sites.google.com/view/ictai20-wd3/}{https://sites.google.com/view/ictai20-wd3/}.}.


title: The Social Welfare Implications of the Zenga Index

abstract: We introduce the social welfare implications of the Zenga index, a recently proposed index of inequality. Our proposal is derived by following the seminal book by Son (2011) and the recent working paper by Kakwani and Son (2019). We compare the Zenga based approach with the classical one, based on the Lorenz curve and the Gini coefficient, as well as the Bonferroni index. We show that the social welfare specification based on the Zenga uniformity curve presents some peculiarities that distinguish it from the other considered indexes. The social welfare specification presented here provides a deeper understanding of how the Zenga index evaluates the inequality in a distribution.


title: Effects of Non-Cognitive Factors on Post-Secondary Persistence of Deaf Students: An Agent-Based Modeling Approach

abstract: Post-secondary education persistence is the likelihood of a student remaining in post-secondary education. Although statistics show that post-secondary persistence for deaf students has increased recently, there are still many obstacles obstructing students from completing their post-secondary degree goals. Therefore, increasing the persistence rate is crucial to increase education and work goals for deaf students. In this work, we present an agent-based model using NetLogo software for the persistence phenomena of deaf students. We consider four non-cognitive factors: having clear goals, social integration, social skills, and academic experience, which influence the departure decision of deaf students. Progress and results of this work suggest that agent-based modeling approaches promise to give better understanding of what will increase persistence.