title: The sum of the masses of the Milky Way and M31: a likelihood-free inference approach

abstract: We use Density Estimation Likelihood-Free Inference, $\Lambda$ Cold Dark Matter simulations of $\sim 2M$ galaxy pairs, and data from Gaia and the Hubble Space Telescope to infer the sum of the masses of the Milky Way and Andromeda (M31) galaxies, the two main components of the Local Group. This method overcomes most of the approximations of the traditional timing argument, makes the writing of a theoretical likelihood unnecessary, and allows the non-linear modelling of observational errors that take into account correlations in the data and non-Gaussian distributions. We obtain an $M_{200}$ mass estimate $M_{\rm MW+M31} = 4.6^{+2.3}_{-1.8} \times 10^{12} M_{\odot}$ ($68 \%$ C.L.), in agreement with previous estimates both for the sum of the two masses and for the individual masses. This result is not only one of the most reliable estimates of the sum of the two masses to date, but is also an illustration of likelihood-free inference in a problem with only one parameter and only three data points.


title: Underpinning Theories of Software Engineering: Dynamism in Physical Sources of the Shannon Weaver Communication Model

abstract: This paper aims to contribute to further understanding of dynamism (the dynamic behavior of system models) in the mathematical and conceptual modeling of systems. This study is conducted in the context of the claim that software engineering lacks underpinning scientific theories, both for the software it produces and the processes by which it does so. The research literature proposes that information theory can provide such a benefit for software engineering. We explore the dynamism expressive power of conceptual modeling as a software engineering tool that can represent physical systems in the Shannon Weaver communication model (SWCM). Specifically, the modeled source in the SWCM is a physical phenomenon (a change that can occur in the world, e.g., tossing a coin) resulting in generating observable events and data of unaddressed information. The resultant model reflects the feasibility of extending the SWCM to be applied in conceptual modeling in software engineering.


title: What Can You Learn from Your Muscles? Learning Visual Representation from Human Interactions

abstract: Learning effective representations of visual data that generalize to a variety of downstream tasks has been a long quest for computer vision. Most representation learning approaches rely solely on visual data such as images or videos. In this paper, we explore a novel approach, where we use human interaction and attention cues to investigate whether we can learn better representations compared to visual-only representations. For this study, we collect a dataset of human interactions capturing body part movements and gaze in their daily lives. Our experiments show that our "muscly-supervised" representation that encodes interaction and attention cues outperforms a visual-only state-of-the-art method MoCo (He et al.,2020), on a variety of target tasks: scene classification (semantic), action recognition (temporal), depth estimation (geometric), dynamics prediction (physics) and walkable surface estimation (affordance). Our code and dataset are available at: https://github.com/ehsanik/muscleTorch.


title: Detecting Objectifying Language in Online Professor Reviews

abstract: Student reviews often make reference to professors' physical appearances. Until recently RateMyProfessors.com, the website of this study's focus, used a design feature to encourage a "hot or not" rating of college professors. In the wake of recent #MeToo and #TimesUp movements, social awareness of the inappropriateness of these reviews has grown; however, objectifying comments remain and continue to be posted in this online context. We describe two supervised text classifiers for detecting objectifying commentary in professor reviews. We then ensemble these classifiers and use the resulting model to track objectifying commentary at scale. We measure correlations between objectifying commentary, changes to the review website interface, and teacher gender across a ten-year period.


title: Splendid Morita equivalences for principal blocks with semidihedral defect groups

abstract: We classify principal blocks of finite groups with semidihedral defect groups up to splendid Morita equivalence. This completes the classification of all principal $2$-blocks of tame representation type up to splendid Morita equivalence and shows that Puig's Finiteness Conjecture holds for such blocks.


title: Mischief: A Simple Black-Box Attack Against Transformer Architectures

abstract: We introduce Mischief, a simple and lightweight method to produce a class of human-readable, realistic adversarial examples for language models. We perform exhaustive experimentations of our algorithm on four transformer-based architectures, across a variety of downstream tasks, as well as under varying concentrations of said examples. Our findings show that the presence of Mischief-generated adversarial samples in the test set significantly degrades (by up to $20\%$) the performance of these models with respect to their reported baselines. Nonetheless, we also demonstrate that, by including similar examples in the training set, it is possible to restore the baseline scores on the adversarial test set. Moreover, for certain tasks, the models trained with Mischief set show a modest increase on performance with respect to their original, non-adversarial baseline.


title: Statistics of wide pre-main sequence binaries in the Orion OB1 association

abstract: Statistics of low-mass pre-main sequence binaries in the Orion OB1 association with separations ranging from 0.6" to 20" (220 to 7400 au at 370 pc) are studied using images from the VISTA Orion mini-survey and astrometry from Gaia. The input sample based on the CVSO catalog contains 1137 stars of K and M spectral types (masses between 0.3 and 0.9 Msun), 1021 of which are considered to be association members. There are 135 physical binary companions to these stars with mass ratios above ~0.13. The average companion fraction is 0.09+-0.01 over 1.2 decades in separation, slightly less than, but still consistent with, the field. We found a difference between the Ori OB1a and OB1b groups, the latter being richer in binaries by a factor 1.6+-0.3. No overall dependence of the wide-binary frequency on the observed underlying stellar density is found, although in the Ori OB1a off-cloud population these binaries seem to avoid dense clusters. The multiplicity rates in Ori OB1 and in sparse regions like Taurus differ significantly, hinting that binaries in the field may originate from a mixture of diverse populations.


title: Robust in Practice: Adversarial Attacks on Quantum Machine Learning

abstract: State-of-the-art classical neural networks are observed to be vulnerable to small crafted adversarial perturbations. A more severe vulnerability has been noted for quantum machine learning (QML) models classifying Haar-random pure states. This stems from the concentration of measure phenomenon, a property of the metric space when sampled probabilistically, and is independent of the classification protocol. In order to provide insights into the adversarial robustness of a quantum classifier on real-world classification tasks, we focus on the adversarial robustness in classifying a subset of encoded states that are smoothly generated from a Gaussian latent space. We show that the vulnerability of this task is considerably weaker than that of classifying Haar-random pure states. In particular, we find only mildly polynomially decreasing robustness in the number of qubits, in contrast to the exponentially decreasing robustness when classifying Haar-random pure states and suggesting that QML models can be useful for real-world classification tasks.


title: Sudden collapse of magnetic order in oxygen deficient nickelate films

abstract: Oxygen vacancies play a crucial role in the control of the electronic, magnetic, ionic, and transport properties of functional oxide perovskites. Rare earth nickelates (RENiO$_{3-x}$) have emerged over the years as a rich platform to study the interplay between the lattice, the electronic structure, and ordered magnetism. In this study, we investigate the evolution of the electronic and magnetic structure in thin films of RENiO$_{3-x}$, using a combination of X-ray absorption spectroscopy and imaging, resonant X-ray scattering, and extended multiplet ligand field theory modeling. We find that oxygen vacancies modify the electronic configuration within the Ni-O orbital manifolds, leading to a dramatic evolution of long-range electronic transport pathways despite the absence of nanoscale phase separation. Remarkably, magnetism is robust to substantial levels of carrier doping, and only a moderate weakening of the $(1/4, 1/4, 1/4)_{pc}$ antiferromagnetic order parameter is observed, whereas the magnetic transition temperature is largely unchanged. Only at a certain point long-range magnetism is abruptly erased without an accompanying structural transition. We propose the progressive disruption of the 3D magnetic superexchange pathways upon introduction of point defects as the mechanism behind the sudden collapse of magnetic order in oxygen-deficient nickelates. Our work demonstrates that, unlike most other oxides, ordered magnetism in RENiO$_{3-x}$ is mostly insensitive to carrier doping. The sudden collapse of ordered magnetism upon oxygen removal may provide a new mechanism for solid-state magneto-ionic switching and new applications in antiferromagnetic spintronics.


title: A Generative Model based Adversarial Security of Deep Learning and Linear Classifier Models

abstract: In recent years, machine learning algorithms have been applied widely in various fields such as health, transportation, and the autonomous car. With the rapid developments of deep learning techniques, it is critical to take the security concern into account for the application of the algorithms. While machine learning offers significant advantages in terms of the application of algorithms, the issue of security is ignored. Since it has many applications in the real world, security is a vital part of the algorithms. In this paper, we have proposed a mitigation method for adversarial attacks against machine learning models with an autoencoder model that is one of the generative ones. The main idea behind adversarial attacks against machine learning models is to produce erroneous results by manipulating trained models. We have also presented the performance of autoencoder models to various attack methods from deep neural networks to traditional algorithms by using different methods such as non-targeted and targeted attacks to multi-class logistic regression, a fast gradient sign method, a targeted fast gradient sign method and a basic iterative method attack to neural networks for the MNIST dataset.


title: A Unified Model for Recommendation with Selective Neighborhood Modeling

abstract: Neighborhood-based recommenders are a major class of Collaborative Filtering (CF) models. The intuition is to exploit neighbors with similar preferences for bridging unseen user-item pairs and alleviating data sparseness. Many existing works propose neural attention networks to aggregate neighbors and place higher weights on specific subsets of users for recommendation. However, the neighborhood information is not necessarily always informative, and the noises in the neighborhood can negatively affect the model performance. To address this issue, we propose a novel neighborhood-based recommender, where a hybrid gated network is designed to automatically separate similar neighbors from dissimilar (noisy) ones, and aggregate those similar neighbors to comprise neighborhood representations. The confidence in the neighborhood is also addressed by putting higher weights on the neighborhood representations if we are confident with the neighborhood information, and vice versa. In addition, a user-neighbor component is proposed to explicitly regularize user-neighbor proximity in the latent space. These two components are combined into a unified model to complement each other for the recommendation task. Extensive experiments on three publicly available datasets show that the proposed model consistently outperforms state-of-the-art neighborhood-based recommenders. We also study different variants of the proposed model to justify the underlying intuition of the proposed hybrid gated network and user-neighbor modeling components.


title: Characterizing the Latent Space of Molecular Deep Generative Models with Persistent Homology Metrics

abstract: Deep generative models are increasingly becoming integral parts of the in silico molecule design pipeline and have dual goals of learning the chemical and structural features that render candidate molecules viable while also being flexible enough to generate novel designs. Specifically, Variational Auto Encoders (VAEs) are generative models in which encoder-decoder network pairs are trained to reconstruct training data distributions in such a way that the latent space of the encoder network is smooth. Therefore, novel candidates can be found by sampling from this latent space. However, the scope of architectures and hyperparameters is vast and choosing the best combination for in silico discovery has important implications for downstream success. Therefore, it is important to develop a principled methodology for distinguishing how well a given generative model is able to learn salient molecular features. In this work, we propose a method for measuring how well the latent space of deep generative models is able to encode structural and chemical features of molecular datasets by correlating latent space metrics with metrics from the field of topological data analysis (TDA). We apply our evaluation methodology to a VAE trained on SMILES strings and show that 3D topology information is consistently encoded throughout the latent space of the model.


title: What can we learn from FRB 200428?

abstract: The two radio pulses (TRPs) from SGR 1935+2154 detected by \cite{BKRMHB} and \cite{2020arXiv200510324T} have similar features to that of cosmological fast radio bursts (FRBs). Many authors directly call the TRPs as FRB 200428 without consider two questions carefully. (1) Are the TRPs just two brighter subpulses of a normal radio pulse liking the normal radio pulses seen in other magnetars during their outburst? (2) If the TRPs are two bursts of an FRB, does this make other FRBs difficult to understand? In this paper, we try to clarify these two questions. First, we compare these TRPs with previous observations of normal radio pulses from pulsars and magnetars, and find that the TRPs should be produced by different mechanism from that of normal radio pulses. We then investigate the second question by assuming the TRPs have the same origin as that of periodically repeating FRBs and find that the origin of the periodicity of periodically repeating FRBs should not be induced by precession. Otherwise, we should expect that the TRPs are not an FRB, at least the TRPs don't have the same origin as that of periodically repeating FRBs.


title: Phenomenology of diagrams in Book II of the Elements

abstract: In this paper, we provide an interpretation of Book II of the Elements from the perspective of figures which are represented and not represented on the diagrams. We show that Euclid's reliance on figures not represented on the diagram is a proof technique which enables to turn his diagrams II.11--14 into ideograms of a kind. We also discuss interpretations of Book II developed by J. Baldwin and A. Mueller, L. Corry, D. Fowler, R. Hartshorne, I. Mueller, K. Saito, and the so-called geometric algebraic interpretation in B. van der Waerden's version.


title: Into the UV: The Atmosphere of the Hot Jupiter HAT-P-41b Revealed

abstract: For solar-system objects, ultraviolet spectroscopy has been critical in identifying sources for stratospheric heating and measuring the abundances of a variety of hydrocarbon and sulfur-bearing species, produced via photochemical mechanisms, as well as oxygen and ozone. To date, less than 20 exoplanets have been probed in this critical wavelength range (0.2-0.4 um). Here we use data from Hubble's newly implemented WFC3 UVIS G280 grism to probe the atmosphere of the hot Jupiter HAT-P-41b in the ultraviolet through optical in combination with observations at infrared wavelengths. We analyze and interpret HAT-P-41b's 0.2-5.0 um transmission spectrum using a broad range of methodologies including multiple treatments of data systematics as well as comparisons with atmospheric forward, cloud microphysical, and multiple atmospheric retrieval models. Although some analysis and interpretation methods favor the presence of clouds or potentially a combination of Na, VO, AlO, and CrH to explain the ultraviolet through optical portions of HAT-P-41b's transmission spectrum, we find that the presence of a significant H- opacity provides the most robust explanation. We obtain a constraint for the abundance of H-, log(H-) = -8.65 +/- 0.62 in HAT-P-41b's atmosphere, which is several orders of magnitude larger than predictions from equilibrium chemistry for a 1700 - 1950 K hot Jupiter. We show that a combination of photochemical and collisional processes on hot hydrogen-dominated exoplanets can readily supply the necessary amount of H- and suggest that such processes are at work in HAT-P-41b and many other hot Jupiter atmospheres.


title: Yangian Bootstrap for Massive Feynman Integrals

abstract: We extend the study of the recently discovered Yangian symmetry of massive Feynman integrals and its relation to massive momentum space conformal symmetry. After proving the symmetry statements in detail at one and two loop orders, we employ the conformal and Yangian constraints to bootstrap various one-loop examples of massive Feynman integrals. In particular, we explore the interplay between Yangian symmetry and hypergeometric expressions of the considered integrals. Based on these examples we conjecture single series representations for all dual conformal one-loop integrals in D spacetime dimensions with generic massive propagators.


title: Sleeping beasts: strong toroidal magnetic field in quiescent magnetars explains their large pulsed fraction

abstract: Magnetars are neutron stars (NSs) with extreme magnetic fields of strength $5 \times 10^{13}$ - $10^{15}$ G. They exhibit transient, highly energetic events, such as short X-ray flashes, bursts and giant flares, all of which are powered by their enormous magnetic energy. Quiescent magnetars have X-ray luminosities between $10^{29}$ and $10^{35}$ erg/s, and are further classified as either persistent or transient magnetars. Their X-ray emission is modulated with the rotational period of the NS, with a typical relative amplitude (so-called pulsed fraction) between 10-58 per cent, implying that the surface temperature is significantly non-uniform despite the high thermal conductivity of the star's crust. Here, we present the first 3D magneto-thermal MHD simulations of magnetars with strong toroidal magnetic fields. We show that these models, combined with ray propagation in curved space-time, accurately describe the light-curves of most transient magnetars in quiescence and allow us to further constrain their rotational orientation. We find that the presence of a strong toroidal magnetic field explains the observed asymmetry in the surface temperature, and is the main cause of the strong modulation of thermal X-ray emission in quiescence.


title: A Family Tree of Optical Transients from Narrow-Line Seyfert 1 Galaxies

abstract: The Zwicky Transient Facility (ZTF) has discovered five new events belonging to an emerging class of AGN undergoing smooth flares with large amplitudes and rapid rise times. This sample consists of several transients that were initially classified as supernovae with narrow spectral lines. However, upon closer inspection, all of the host galaxies display resolved Balmer lines characteristic of a narrow-line Seyfert 1 (NLSy1) galaxy. The transient events are long-lived, over 400 days on average. We report UV and X-ray follow-up of the flares and observe persistent UV-bright emission, with two of the five transients detected with luminous X-ray emission, ruling out a supernova interpretation. We compare the properties of this sample to previously reported flaring NLSy1 galaxies, and find that they fall into three spectroscopic categories: Transients with 1) Balmer line profiles and Fe II complexes typical of NLSy1s, 2) strong He II profiles, and 3) He II profiles including Bowen fluorescence features. The latter are members of the growing class of AGN flares attributed to enhanced accretion reported by Trakhtenbrot et al. (2019). We consider physical interpretations in the context of related transients from the literature. For example, two of the sources show high amplitude rebrightening in the optical, ruling out a simple tidal disruption event scenario for those transients. We conclude that three of the sample belong to the Trakhtenbrot et al. (2019) class, and two are TDEs in NLSy1s. We also aim to understand why NLSy1s are preferentially the sites of such rapid enhanced flaring activity.


title: Dynamical evolution of massive perturbers in realistic multi-component galaxy models I: implementation and validation

abstract: Galaxies are self-gravitating structures composed by several components encompassing spherical, axial and triaxial symmetry. Although real systems feature heterogeneous components whose properties are intimately connected, semi-analytical approaches often exploit the linearity of the Poisson's equation to represent the potential and mass distribution of a multi-component galaxy as the sum of the individual components. In this work, we expand the semi-analytical framework developed in Bonetti et al. (2020) by including both a detailed implementation of the gravitational potential of exponential disc (modelled with a ${\rm sech}^2$ and an exponential vertical profile) and an accurate prescription for the dynamical friction experienced by massive perturbers in composite galaxy models featuring rotating disc structures. Such improvements allow us to evolve arbitrary orbits either within or outside the galactic disc plane. We validate the results obtained by our numerical model against public semi-analytical codes as well as full N-body simulations, finding that our model is in excellent agreement to the codes it is compared with. The ability to reproduce the relevant physical processes responsible for the evolution of massive perturber orbits and its computational efficiency make our framework perfectly suited for large parameter-space exploration studies.


title: FPGAs-as-a-Service Toolkit (FaaST)

abstract: Computing needs for high energy physics are already intensive and are expected to increase drastically in the coming years. In this context, heterogeneous computing, specifically as-a-service computing, has the potential for significant gains over traditional computing models. Although previous studies and packages in the field of heterogeneous computing have focused on GPUs as accelerators, FPGAs are an extremely promising option as well. A series of workflows are developed to establish the performance capabilities of FPGAs as a service. Multiple different devices and a range of algorithms for use in high energy physics are studied. For a small, dense network, the throughput can be improved by an order of magnitude with respect to GPUs as a service. For large convolutional networks, the throughput is found to be comparable to GPUs as a service. This work represents the first open-source FPGAs-as-a-service toolkit.

categories: physics.comp-ph cs.DC hep-ex physics.data-an physics.ins-det


title: Can the EHT M87 results be used to test general relativity?

abstract: No. All theoretical predictions for the observational appearance of an accreting supermassive black hole, as measured interferometrically by a sparse Earth-sized array at current observation frequencies, are sensitive to many untested assumptions about accretion flow and emission physics. There is no way to distinguish a violation of general relativity (GR) from the much more likely scenario that the relevant "gastrophysical" assumptions simply do not hold. Tests of GR will become possible with longer interferometric baselines (likely requiring a space mission) that reach the resolution where astrophysics-independent predictions of the theory become observable.


title: Pole skipping away from maximal chaos

abstract: Pole skipping is a recently discovered subtle effect in the thermal energy density retarded two point function at a special point in the complex $(\omega,p)$ planes. We propose that pole skipping is determined by the stress tensor contribution to many-body chaos, and the special point is at $(\omega,p)_\text{p.s.}= i \lambda^{(T)}(1,1/u_B^{(T)})$, where $\lambda^{(T)}=2\pi/\beta$ and $u_B^{(T)}$ are the stress tensor contributions to the Lyapunov exponent and the butterfly velocity respectively. While this proposal is consistent with previous studies conducted for maximally chaotic theories, where the stress tensor dominates chaos, it clarifies that one cannot use pole skipping to extract the Lyapunov exponent of a theory, which obeys $\lambda\leq \lambda^{(T)}$. On the other hand, in a large class of strongly coupled but non-maximally chaotic theories $u_B^{(T)}$ is the true butterfly velocity and we conjecture that $u_B\leq u_B^{(T)}$ is a universal bound. While it remains a challenge to explain pole skipping in a general framework, we provide a stringent test of our proposal in the large-$q$ limit of the SYK chain, where we determine $\lambda,\, u_B,$ and the energy density two point function in closed form for all values of the coupling, interpolating between the free and maximally chaotic limits. Since such an explicit expression for a thermal correlator is one of a kind, we take the opportunity to analyze many of its properties: the coupling dependence of the diffusion constant, the dispersion relations of poles, and the convergence properties of all order hydrodynamics.


title: Leading Nonlinear Tidal Effects and Scattering Amplitudes

abstract: We present the two-body Hamiltonian and associated eikonal phase, to leading post-Minkowskian order, for infinitely many tidal deformations described by operators with arbitrary powers of the curvature tensor. Scattering amplitudes in momentum and position space provide systematic complementary approaches. For the tidal operators quadratic in curvature, which describe the linear response to an external gravitational field, we work out the leading post-Minkowskian contributions using a basis of operators with arbitrary numbers of derivatives which are in one-to-one correspondence with the worldline multipole operators. Explicit examples are used to show that the same techniques apply to both bodies interacting tidally with a spinning particle, for which we find the leading contributions from quadratic in curvature tidal operators with an arbitrary number of derivatives, and to effective field theory extensions of general relativity. We also note that the leading post-Minkowskian order contributions from higher-dimension operators manifest double-copy relations. Finally, we comment on the structure of higher-order corrections.


title: EFT Asymptotics: the Growth of Operator Degeneracy

abstract: We establish formulae for the asymptotic growth (with respect to the scaling dimension) of the number of operators in effective field theory, or equivalently the number of $S$-matrix elements, in arbitrary spacetime dimensions and with generic field content. This we achieve by generalising a theorem due to Meinardus and applying it to Hilbert series -- partition functions for the degeneracy of (subsets of) operators. Although our formulae are asymptotic, numerical experiments reveal remarkable agreement with exact results at very low orders in the EFT expansion, including for complicated phenomenological theories such as the standard model EFT. Our methods also reveal phase transition-like behaviour in Hilbert series. We discuss prospects for tightening the bounds and providing rigorous errors to the growth of operator degeneracy, and of extending the analytic study and utility of Hilbert series to EFT.


title: Differentiable Quantum Architecture Search

abstract: Quantum architecture search (QAS) is the process of automating architecture engineering of quantum circuits. It has been desired to construct a powerful and general QAS platform which can significantly accelerate current efforts to identify quantum advantages of error-prone and depth-limited quantum circuits in the NISQ era. Hereby, we propose a general framework of differentiable quantum architecture search (DQAS), which enables automated designs of quantum circuits in an end-to-end differentiable fashion. We present several examples of circuit design problems to demonstrate the power of DQAS. For instance, unitary operations are decomposed into quantum gates, noisy circuits are re-designed to improve accuracy, and circuit layouts for quantum approximation optimization algorithm are automatically discovered and upgraded for combinatorial optimization problems. These results not only manifest the vast potential of DQAS being an essential tool for the NISQ application developments, but also present an interesting research topic from the theoretical perspective as it draws inspirations from the newly emerging interdisciplinary paradigms of differentiable programming, probabilistic programming, and quantum programming.


title: The impact of inelastic self-interacting dark matter on the dark matter structure of a Milky Way halo

abstract: We study the effects of inelastic dark matter self-interactions on the internal structure of a simulated Milky Way (MW)-size halo. Self-interacting dark matter (SIDM) is an alternative to collisionless cold dark matter (CDM) which offers a unique solution to the problems encountered with CDM on sub-galactic scales. Although previous SIDM simulations have mainly considered elastic collisions, theoretical considerations motivate the existence of multi-state dark matter where transitions from the excited to the ground state are exothermic. In this work, we consider a self-interacting, two-state dark matter model with inelastic collisions, implemented in the Arepo code. We find that energy injection from inelastic self-interactions reduces the central density of the MW halo in a shorter timescale relative to the elastic scale, resulting in a larger core size. Inelastic collisions also isotropize the orbits, resulting in an overall lower velocity anisotropy for the inelastic MW halo. In the inner halo, the inelastic SIDM case (minor-to-major axis ratio $s \equiv c/a \approx 0.65$) is more spherical than the CDM ($s \approx 0.4$), but less spherical than the elastic SIDM case ($s \approx 0.75$). The speed distribution $f(v)$ of dark matter particles at the location of the Sun in the inelastic SIDM model shows a significant departure from the CDM model, with $f(v)$ falling more steeply at high speeds. In addition, the velocity kicks imparted during inelastic collisions produce unbound high-speed particles with velocities up to 500 km s$^{-1}$ throughout the halo. This implies that inelastic SIDM can potentially leave distinct signatures in direct detection experiments, relative to elastic SIDM and CDM.


title: Dark Matter Interpretation of the Fermi-LAT Observations Toward the Outer Halo of M31

abstract: An excess $\gamma$-ray signal toward the outer halo of M31 has recently been reported. Although other explanations are plausible, the possibility that it arises from dark matter (DM) is valid. In this work we interpret the excess in the framework of DM annihilation, using as our representative case WIMP DM annihilating to bottom quarks, and we perform a detailed study of the systematic uncertainty in the $J$-factor for the M31 field. We find that the signal favors a DM particle with a mass of $\sim$46-73 GeV. While the mass is well constrained, the systematic uncertainty in the cross-section spans 2.5 orders of magnitude, ranging from $\sim$8$\times 10^{-27}-4 \times 10^{-24} \ \mathrm{cm^3 \ s^{-1}}$. This high uncertainty is due to two main factors, namely, an uncertainty in the substructure nature and geometry of the DM halos for both M31 and the Milky Way (MW), and correspondingly, an uncertainty in the contribution to the signal from the MW's DM halo along the line of sight. However, under the conditions that the minimum subhalo mass is $\lesssim 10^{-6} \ M_\odot$ and the actual contribution from the MW's DM halo along the line of sight is at least $\sim$30$\%$ of its total value, we show that there is a large overlap with the DM interpretations of both the Galactic center (GC) excess and the antiproton excess, while also being compatible with the limits for the MW dwarf spheroidals. More generally, we summarize the results from numerous complementary DM searches in the energy range 10 GeV $-$ 300 GeV corresponding to the GC excess and identify a region in parameter space that still remains viable for discovery of the DM particle.


title: Review of experiments on the chiral anomaly in Dirac-Weyl semimetals

abstract: We provide a review of recent experimental results on the chiral anomaly in Dirac/Weyl semimetals. After a brief introduction, we trace the steps leading to the prediction of materials that feature protected 3D bulk Dirac nodes. The chiral anomaly is presented in terms of charge pumping between the chiral Landau levels of Weyl fermions in parallel electric and magnetic fields. The related chiral magnetic effect and chiral zero sound are described. Current jetting effects, which present major complications in experiments on the longitudinal magnetoresistance, are carefully analyzed. We describe a recent test that is capable of distinguishing these semiclassical artifacts from intrinsic quantum effects. Turning to experiments, we review critically the longitudinal magnetoresistance experiments in the Dirac/Weyl semimetals Na3Bi, GdPtBi, ZrTe5 and TaAs. Alternate approaches to the chiral anomaly, including experiments on non-local transport, thermopower, thermal conductivity and optical pump-probe response are reviewed. In the Supplement, we provide a brief discussion of the chiral anomaly in the broader context of high energy physics and relativistic quantum field theory, as well the anomaly's starring role at the nexus of quantum physics and differential geometry.


title: WISDOM project -- VI. Exploring the relation between supermassive black hole mass and galaxy rotation with molecular gas

abstract: Empirical correlations between the masses of supermassive black holes (SMBHs) and properties of their host galaxies are well-established. Among these is the correlation with the flat rotation velocity of each galaxy measured either at a large radius in its rotation curve or via a spatially-integrated emission line width. We propose here the use of the de-projected integrated CO emission line width as an alternative tracer of this rotation velocity, that has already been shown useful for the Tully-Fisher (luminosity-rotation velocity) relation. We investigate the correlation between CO line widths and SMBH masses for two samples of galaxies with dynamical SMBH mass measurements, with respectively spatially-resolved and unresolved CO observations. The tightest correlation is found using the resolved sample of 24 galaxies as $\log (M_\mathrm{BH}/\mathrm{M_\odot})=(7.5\pm0.1)+(8.5\pm0.9)[\log(W_\mathrm{50}/\sin i \,\mathrm{km\,s}^{-1})-2.7]$, where $M_\mathrm{BH}$ is the central SMBH mass, $W_{50}$ the full-width at half-maximum of a double-horned emission line profile, and $i$ the inclination of the CO disc. This relation has a total scatter of $0.6\,$dex, comparable to those of other SMBH mass correlations, and dominated by the intrinsic scatter of $0.5\,$dex. A tight correlation is also found between the de-projected CO line widths and the stellar velocity dispersions averaged within one effective radius. We apply our correlation to the COLD GASS sample to estimate the local SMBH mass function.


title: Reflective Decoding: Beyond Unidirectional Generation with Off-the-Shelf Language Models

abstract: Publicly available, large pretrained LanguageModels (LMs) generate text with remarkable quality, but only sequentially from left to right. As a result, they are not immediately applicable to generation tasks that break the unidirectional assumption, such as paraphrasing or text-infilling, necessitating task-specific supervision. In this paper, we present Reflective Decoding, a novel unsupervised algorithm that allows for direct application of unidirectional LMs to non-sequential tasks. Our 2-step approach requires no supervision or even parallel corpora, only two off-the-shelf pretrained LMs in opposite directions: forward and backward. First, in the contextualization step, we use LMs to generate ensembles of past and future contexts which collectively capture the input (e.g. the source sentence for paraphrasing). Second, in the reflection step, we condition on these "context ensembles", generating outputs that are compatible with them. Comprehensive empirical results demonstrate that Reflective Decoding outperforms strong unsupervised baselines on both paraphrasing and abductive text infilling, significantly narrowing the gap between unsupervised and supervised methods. Reflective Decoding surpasses multiple supervised baselines on various metrics including human evaluation.


title: Infinite staircases for Hirzebruch surfaces

abstract: We consider the embedding capacity functions $c_{H_b}(z)$ for symplectic embeddings of ellipsoids of eccentricity $z$ into the family of nontrivial rational Hirzebruch surfaces $H_b$ with symplectic form parametrized by $b\in [0,1)$. This function was known to have an infinite staircase in the monotone cases ($b= 0$ and $ b= 1/3$). It is also known that for each $b$ there is at most one value of $z$ that can be the accumulation point of such a staircase. In this manuscript, we identify three sequences of open, disjoint, blocked $b$-intervals, consisting of $b$-parameters where the embedding capacity function for $H_b$ does not contain an infinite staircase. There is one sequence in each of the intervals $(0,1/5)$, $(1/5,1/3)$, and $(1/3,1)$. We then establish six sequences of associated infinite staircases, one occurring at each endpoint of the blocked $b$-intervals. The staircase numerics are variants of those in the Fibonacci staircase for the projective plane (the case $b=0$). We also show that there is no staircase at the point $b=1/5$, even though this value is not blocked. The focus of this paper is to develop techniques, both graphical and numeric, that allow identification of potential staircases, and then to understand the obstructions well enough to prove that the purported staircases really do have the required properties. A subsequent paper will explore in more depth the set of $b$ that admit infinite staircases.


title: Mining the Geodesic Equation for Scattering Data

abstract: The geodesic equation encodes test-particle dynamics at arbitrary gravitational coupling, hence retaining all orders in the post-Minkowskian (PM) expansion. Here we explore what geodesic motion can tell us about dynamical scattering in the presence of perturbatively small effects such as tidal distortion and higher derivative corrections to general relativity. We derive an algebraic map between the perturbed geodesic equation and the leading PM scattering amplitude at arbitrary mass ratio. As examples, we compute formulas for amplitudes and isotropic gauge Hamiltonians for certain infinite classes of tidal operators such as electric or magnetic Weyl to any power, and for higher derivative corrections to gravitationally interacting bodies with or without electric charge. Finally, we present a general method for calculating closed-form expressions for amplitudes and isotropic gauge Hamiltonians in the test-particle limit at all orders in the PM expansion.


title: Generalizable Machine Learning in Neuroscience using Graph Neural Networks

abstract: Although a number of studies have explored deep learning in neuroscience, the application of these algorithms to neural systems on a microscopic scale, i.e. parameters relevant to lower scales of organization, remains relatively novel. Motivated by advances in whole-brain imaging, we examined the performance of deep learning models on microscopic neural dynamics and resulting emergent behaviors using calcium imaging data from the nematode C. elegans. We show that neural networks perform remarkably well on both neuron-level dynamics prediction, and behavioral state classification. In addition, we compared the performance of structure agnostic neural networks and graph neural networks to investigate if graph structure can be exploited as a favorable inductive bias. To perform this experiment, we designed a graph neural network which explicitly infers relations between neurons from neural activity and leverages the inferred graph structure during computations. In our experiments, we found that graph neural networks generally outperformed structure agnostic models and excel in generalization on unseen organisms, implying a potential path to generalizable machine learning in neuroscience.


title: Generating Fact Checking Summaries for Web Claims

abstract: We present SUMO, a neural attention-based approach that learns to establish the correctness of textual claims based on evidence in the form of text documents (e.g., news articles or Web documents). SUMO further generates an extractive summary by presenting a diversified set of sentences from the documents that explain its decision on the correctness of the textual claim. Prior approaches to address the problem of fact checking and evidence extraction have relied on simple concatenation of claim and document word embeddings as an input to claim driven attention weight computation. This is done so as to extract salient words and sentences from the documents that help establish the correctness of the claim. However, this design of claim-driven attention does not capture the contextual information in documents properly. We improve on the prior art by using improved claim and title guided hierarchical attention to model effective contextual cues. We show the efficacy of our approach on datasets concerning political, healthcare, and environmental issues.


title: Planes of satellites around Milky Way/M31-mass galaxies in the FIRE simulations and comparisons with the Local Group

abstract: We examine the prevalence, longevity, and causes of planes of satellite dwarf galaxies, as observed in the Local Group. We use 14 Milky Way/Andromeda-(MW/M31) mass host galaxies from the FIRE-2 simulations. We select the 14 most massive satellites by stellar mass within 300 kpc of each host and correct for incompleteness from the foreground galactic disc when comparing to the MW. We find that MW-like planes as spatially thin and/or kinematically coherent as observed are uncommon, but they do exist in our simulations. Spatially thin planes occur in 1-2 per cent of snapshots during $z=0-0.2$, and kinematically coherent planes occur in 5 per cent of snapshots. These planes are generally transient, surviving for less than 500 Myr. However, if we select hosts with an LMC-like satellite near first pericentre, the fraction of snapshots with MW-like planes increases dramatically to 7-16 per cent, with lifetimes of 0.7-1 Gyr, likely because of group accretion of satellites. We find that M31's satellite distribution is much more common: M31's satellites lie within about 1 sigma of the simulation median for every plane metric we consider. We find no significant difference in average satellite planarity for isolated hosts versus hosts in LG-like pairs. Baryonic and dark matter-only simulations exhibit similar levels of planarity, even though baryonic subhaloes are less centrally concentrated within their host haloes. We conclude that planes of satellites are not a strong challenge to LCDM cosmology.


title: Horizon-independent Preconditioner Design for Linear Predictive Control

abstract: First-order optimization solvers, such as the Fast Gradient Method, are increasingly being used to solve Model Predictive Control problems in resource-constrained environments. Unfortunately, the convergence rate of these solvers is significantly affected by the conditioning of the problem data, with ill-conditioned problems requiring a large number of iterations. To reduce the number of iterations required, we present a simple method for computing a horizon-independent preconditioning matrix for the Hessian of the condensed problem. The preconditioner is based on the block Toeplitz structure of the Hessian. Horizon-independence allows one to use only the predicted system and cost matrices to compute the preconditioner, instead of the full Hessian. The proposed preconditioner has equivalent performance to an optimal preconditioner in numerical examples, producing speedups between 2x and 9x for the Fast Gradient Method. Additionally, we derive horizon-independent spectral bounds for the Hessian in terms of the transfer function of the predicted system, and show how these can be used to compute a novel horizon-independent bound on the condition number for the preconditioned Hessian.


title: Maximal couplings of the Metropolis-Hastings algorithm

abstract: Couplings play a central role in the analysis of Markov chain Monte Carlo algorithms and appear increasingly often in the algorithms themselves, e.g. in convergence diagnostics, parallelization, and variance reduction techniques. Existing couplings of the Metropolis-Hastings algorithm handle the proposal and acceptance steps separately and fall short of the upper bound on one-step meeting probabilities given by the coupling inequality. This paper introduces maximal couplings which achieve this bound while retaining the practical advantages of current methods. We consider the properties of these couplings and examine their behavior on a selection of numerical examples.


title: Non-intrusive speech intelligibility prediction using automatic speech recognition derived measures

abstract: The estimation of speech intelligibility is still far from being a solved problem. Especially one aspect is problematic: most of the standard models require a clean reference signal in order to estimate intelligibility. This is an issue of some significance, as a reference signal is often unavailable in practice. In this work, therefore a non-intrusive speech intelligibility estimation framework is presented. In it, human listeners' performance in keyword recognition tasks is predicted using intelligibility measures that are derived from models trained for automatic speech recognition (ASR). One such ASR-based and one signal-based measure are combined into a full framework, the proposed NO-Reference Intelligibility (Nori) estimator, which is evaluated in predicting the performance of both normal-hearing and hearing-impaired listeners in multiple noise conditions. It is shown that the Nori framework even outperforms the widely used reference-based (or intrusive) short-term objective intelligibility (STOI) measure in most considered scenarios, while being applicable in fully blind scenarios with no reference signal or transcription, creating perspectives for online and personalized optimization of speech enhancement systems.


title: Search for associations containing young stars (SACY) VIII. An updated census of spectroscopic binary systems showing hints of non-universal multiplicity among these associations

abstract: We seek to update the spectroscopy binary fraction of the SACY (Search for Associations Containing Young stars) sample taking in consideration all possible biases in our identification of binary candidates, such as activity and rotation. Using high-resolution spectroscopic observations we have produced $\sim$1300 cross-correlation functions (CCFs) to disentangle the previously mentioned sources of contamination. The radial velocity values obtained were cross-matched with the literature and were used to revise and update the spectroscopic binary (SB) fraction in each of the SACY association. In order to better describe the CCF profile, we calculated a set of high-order cross-correlation features to determine the origin of the variations in radial velocities. We identified 68 SB candidates from our sample of 410 objects. Our results hint that the youngest associations have a higher SB fraction. Specifically, we found sensitivity-corrected SB fractions of $22 \substack{+15 \\ -11} \%$ for $\epsilon$~Cha , $31 \substack{+16 \\ -14} \%$ for TW Hya and $32 \substack{+9 \\ -8} \%$ for $\beta$~Pictoris, in contrast with the five oldest ($\sim 35-125$ Myr) associations we have sampled which are $\sim 10\%$ or lower. This result seems independent of the methodology used to asses membership to the associations. The new CCF analysis, radial velocity estimates and SB candidates are particularly relevant for membership revision of targets in young stellar associations. These targets would be ideal candidates for follow-up campaigns using high-resolution techniques in order to confirm binarity, resolve the orbits, and ideally calculate dynamical masses. Additionally, if the results on SB fraction in the youngest associations are confirmed, it could hint of non-universal multiplicity among SACY associations.


title: Improving Schroeppel and Shamir's Algorithm for Subset Sum via Orthogonal Vectors

abstract: We present an $\mathcal{O}^\star(2^{0.5n})$ time and $\mathcal{O}^\star(2^{0.249999n})$ space randomized algorithm for solving worst-case Subset Sum instances with $n$ integers. This is the first improvement over the long-standing $\mathcal{O}^\star(2^{n/2})$ time and $\mathcal{O}^\star(2^{n/4})$ space algorithm due to Schroeppel and Shamir (FOCS 1979). We breach this gap in two steps: (1) We present a space efficient reduction to the Orthogonal Vectors Problem (OV), one of the most central problem in Fine-Grained Complexity. The reduction is established via an intricate combination of the method of Schroeppel and Shamir, and the representation technique introduced by Howgrave-Graham and Joux (EUROCRYPT 2010) for designing Subset Sum algorithms for the average case regime. (2) We provide an algorithm for OV that detects an orthogonal pair among $N$ given vectors in $\{0,1\}^d$ with support size $d/4$ in time $\tilde{O}(N\cdot2^d/\binom{d}{d/4})$. Our algorithm for OV is based on and refines the representative families framework developed by Fomin, Lokshtanov, Panolan and Saurabh (J. ACM 2016). Our reduction uncovers a curious tight relation between Subset Sum and OV, because any improvement of our algorithm for OV would imply an improvement over the runtime of Schroeppel and Shamir, which is also a long standing open problem.


title: Direct Numerical Simulation of a Turbulent Boundary Layer with Strong Pressure Gradients

abstract: The turbulent boundary layer over a Gaussian shaped bump is computed by direct numerical simulation (DNS) of the incompressible Navier-Stokes equations. The two-dimensional bump causes a series of strong pressure gradients alternating in rapid succession. At the inflow, the momentum thickness Reynolds number is approximately 1,000 and the boundary layer thickness is 1/8 of the bump height. DNS results show that the strong favorable pressure gradient (FPG) causes the boundary layer to enter a relaminarization process. The near-wall turbulence is significantly weakened and becomes intermittent, however relaminarization does not complete. The streamwise velocity profiles deviate above the standard logarithmic law and the Reynolds shear stress is reduced. The strong acceleration also suppresses the wall-shear normalized turbulent kinetic energy production rate. At the bump peak, where the FPG switches to an adverse gradient (APG), the near-wall turbulence is suddenly enhanced through a partial retransition process. This results in a new highly energized internal layer which is more resilient to the strong APG and only produces incipient flow separation on the downstream side. In the strong FPG and APG regions, the inner and outer layers become largely independent of each other. The near-wall region responds to the pressure gradients and determines the skin friction. The outer layer behaves similarly to a free-shear layer subject to pressure gradients and mean streamline curvature effects. Results from a RANS simulation of the bump are also discussed and clearly show the lack of predictive capacity of the near-wall pressure gradient effects on the mean flow.


title: Positivity and convexity in incomplete cooperative games

abstract: Incomplete cooperative games generalise the classical model of cooperative games by omitting the values of some of the coalitions. This allows to incorporate uncertainty into the model and study the underlying games as well as possible payoff distribution based only on the partial information. In this paper we perform a systematic study of incomplete games, focusing on two important classes of cooperative games: positive and convex games. Regarding positivity, we generalise previous results for a special class of minimal incomplete games to general setting. We characterise non-extendability to a positive game by the existence of a certificate and provide a description of the set of positive extensions using its extreme games. The results are then used to obtain explicit formulas for several classes of incomplete games with special structures. The second part deals with convexity. We begin with considering the case of non-negative minimal incomplete games. Then we survey existing results in the related theory of set functions, namely providing context to the problem of completing partial functions. We provide a characterisation of extendability and a full description of the set of symmetric convex extensions. The set serves as an approximation of the set of convex extensions. Finally, we outline an entirely new perspective on a connection between incomplete cooperative games and cooperative interval games.


title: Effective isotrivial Mordell-Lang in positive characteristic

abstract: The isotrivial Mordell-Lang theorem of Moosa and Scanlon describes the set $X\cap\Gamma$ when $X$ is a subvariety of a semiabelian variety $G$ over a finite field $\mathbb{F}_q$ and $\Gamma$ is a finitely generated subgroup of $G$ that is invariant under the $q$-power Frobenius endomorphism $F$. That description is here made effective, and extended to arbitrary commutative algebraic groups $G$ and arbitrary finitely generated $\mathbb{Z}[F]$-submodules $\Gamma$. The approach is to use finite automata to give a concrete description of $X\cap \Gamma$. These methods and results have new applications even when specialised to the case when $G$ is an abelian variety over a finite field, $X\subseteq G$ a subvariety defined over a function field $K$, and $\Gamma=G(K)$. As an application of the automata-theoretic approach, a dichotomy theorem is established for the growth of the number of points in $X(K)$ of bounded height. As an application of the effective description of $X\cap\Gamma$, decision procedures are given for the following three diophantine problems: Is $X(K)$ nonempty? Is it infinite? Does it contain an infinite coset?


title: Linguistically-Informed Transformations (LIT): A Method for Automatically Generating Contrast Sets

abstract: Although large-scale pretrained language models, such as BERT and RoBERTa, have achieved superhuman performance on in-distribution test sets, their performance suffers on out-of-distribution test sets (e.g., on contrast sets). Building contrast sets often re-quires human-expert annotation, which is expensive and hard to create on a large scale. In this work, we propose a Linguistically-Informed Transformation (LIT) method to automatically generate contrast sets, which enables practitioners to explore linguistic phenomena of interests as well as compose different phenomena. Experimenting with our method on SNLI and MNLI shows that current pretrained language models, although being claimed to contain sufficient linguistic knowledge, struggle on our automatically generated contrast sets. Furthermore, we improve models' performance on the contrast sets by apply-ing LIT to augment the training data, without affecting performance on the original data.


title: A Sampling-Based Method for Tensor Ring Decomposition

abstract: We propose a sampling-based method for computing the tensor ring (TR) decomposition of a data tensor. The method uses leverage score sampled alternating least squares to fit the TR cores in an iterative fashion. By taking advantage of the special structure of TR tensors, we can efficiently estimate the leverage scores and attain a method which has complexity sublinear in the number of input tensor entries. We provide high-probability relative-error guarantees for the sampled least squares problems. We compare our proposal to existing methods in experiments on both synthetic and real data. Our method achieves substantial speedup -- sometimes two or three orders of magnitude -- over competing methods, while maintaining good accuracy. We also provide an example of how our method can be used for rapid feature extraction.


title: CT Image Segmentation for Inflamed and Fibrotic Lungs Using a Multi-Resolution Convolutional Neural Network

abstract: The purpose of this study was to develop a fully-automated segmentation algorithm, robust to various density enhancing lung abnormalities, to facilitate rapid quantitative analysis of computed tomography images. A polymorphic training approach is proposed, in which both specifically labeled left and right lungs of humans with COPD, and nonspecifically labeled lungs of animals with acute lung injury, were incorporated into training a single neural network. The resulting network is intended for predicting left and right lung regions in humans with or without diffuse opacification and consolidation. Performance of the proposed lung segmentation algorithm was extensively evaluated on CT scans of subjects with COPD, confirmed COVID-19, lung cancer, and IPF, despite no labeled training data of the latter three diseases. Lobar segmentations were obtained using the left and right lung segmentation as input to the LobeNet algorithm. Regional lobar analysis was performed using hierarchical clustering to identify radiographic subtypes of COVID-19. The proposed lung segmentation algorithm was quantitatively evaluated using semi-automated and manually-corrected segmentations in 87 COVID-19 CT images, achieving an average symmetric surface distance of $0.495 \pm 0.309$ mm and Dice coefficient of $0.985 \pm 0.011$. Hierarchical clustering identified four radiographical phenotypes of COVID-19 based on lobar fractions of consolidated and poorly aerated tissue. Lower left and lower right lobes were consistently more afflicted with poor aeration and consolidation. However, the most severe cases demonstrated involvement of all lobes. The polymorphic training approach was able to accurately segment COVID-19 cases with diffuse consolidation without requiring COVID-19 cases for training.


title: Spatial distribution of graffiti: A complex network approach

abstract: Despite the great differences among cities, they face similar challenges regarding social inequality, politics and criminality. Urban art express these feelings from the citizen point-of-view. In particular, the drawing and painting of public surfaces may carry rich information about the time and region it was made. Existing studies have explored the spatial distribution of graffiti, but most of them considered graffiti as a whole, with no separation among the types. Also, the analyses rarely take into account the city topology. In this work, we propose to categorize the graffiti into three types: simple scribbles, complex scribbles, canvases. We analyze the spatial distribution and identify the spatial bias of each type. To further analyze the spatial distribution of the types, we apply concepts from complex networks. First, regions (communities) defined by the connectivity profiles of the city network are obtained and the prevalence of each type of graffiti over these regions are analyzed. Next, a measure based on the dynamics of the network (accessibility) is calculated and compared to the distribution of the graffiti types. A case study is performed in the analysis of three different categories of graffiti in the city of Sao Paulo, Brazil. The results showed that the categories present characteristic spatial distributions. The ratio of each type per community of the network, though, does not pose significant deviations. Finally, a small positive correlation was observed between the locations of each graffiti type and the accessibility.


title: Improving significance of binary black hole mergers in Advanced LIGO data using deep learning : Confirmation of GW151216

abstract: We present a novel Machine Learning (ML) based strategy to search for binary black hole (BBH) mergers in data from ground-based gravitational wave (GW) observatories. This is the first ML-based search that not only recovers all the compact binary coalescences (CBCs) in the first GW transients catalog (GWTC-1), but also makes a clean detection of GW151216 by only adding a new coincident ranking statistic (MLStat) to a standard analysis that was used for GWTC-1. In CBC searches, reducing contamination by terrestrial and instrumental transients, which create a loud noise background by triggering numerous false alarms, is crucial to improving the sensitivity for detecting true events. The sheer volume of data and a large number of expected detections also prompts the use of ML techniques. We perform transfer learning to train "InceptionV3", a pre-trained deep neural network, along with curriculum learning to distinguish GW signals from noisy events by analysing their continuous wavelet transform (CWT) maps. MLStat incorporates information from this ML classifier into the coincident search likelihood used by the standard PyCBC search. This leads to at least an order of magnitude improvement in the inverse false-alarm-rate (IFAR) for the previously "low significance" events GW151012, GW170729 and GW151216. We also perform the parameter estimation of GW151216 using SEOBNRv4HM_ROM. We carry out an injection study to show that MLStat brings substantial improvement to the detection sensitivity of Advanced LIGO for all compact binary coalescences. The average improvement in the sensitive volume is ~10% for low chirp masses (0.8-5 Msun), and ~30% for higher masses (5-50 Msun). This work demonstrates the immense potential and readiness of MLStat for finding new sources in current data and the possibility of its adaptation in similar searches.


title: Well-posedness and tamed Euler schemes for McKean-Vlasov equations driven by L\'evy noise

abstract: We prove the well-posedness of solutions to McKean-Vlasov stochastic differential equations driven by L\'evy noise under mild assumptions where, in particular, the L\'evy measure is not required to be finite. The drift, diffusion and jump coefficients are allowed to be random, can grow super-linearly in the state variable, and all may depend on the marginal law of the solution process. We provide a propagation of chaos result under more relaxed conditions than those existing in the literature, and consistent with our well-posedness result. We propose a tamed Euler scheme for the associated interacting particle system and prove that the rate of its strong convergence is arbitrarily close to $1/2$. As a by-product, we also obtain the corresponding results on well-posedness, propagation of chaos and strong convergence of the tamed Euler scheme for McKean-Vlasov stochastic delay differential equations (SDDE) and McKean-Vlasov stochastic differential equations with Markovian switching (SDEwMS), both driven by L\'evy noise. Furthermore, our results on tamed Euler schemes are new even for ordinary SDEs driven by L\'evy noise and with super-linearly growing coefficients.


title: Physical Explanation for the Galaxy Distribution on the $(\lambda_{\rm R}, \varepsilon)$ and $(V/\sigma, \varepsilon)$ Diagrams or for the Limit on Orbital Anisotropy

abstract: In the $(\lambda_{\rm R}, \varepsilon)$ and $(V/\sigma, \varepsilon)$ diagrams for characterizing dynamical states, the fast-rotator galaxies (both early-type and spirals) are distributed within a well-defined leaf-shaped envelope. This was explained as due to an upper limit to the orbital anisotropy increasing with galaxy intrinsic flattening. However, a physical explanation for this empirical trend was missing. Here we construct Jeans Anisotropic Models (JAM), with either cylindrically or spherically aligned velocity ellipsoid (two extreme assumptions), and each with either spatially-constant or -variable anisotropy. We use JAM to build mock samples of axisymmetric galaxies, assuming on average an oblate shape for the velocity ellipsoid (as required to reproduce the rotation of real galaxies), and limiting the radial anisotropy $\beta$ to the range allowed by physical solutions. We find that all four mock samples naturally predict the observed galaxy distribution on the $(\lambda_{\rm R}, \varepsilon)$ and $(V/\sigma, \varepsilon)$ diagrams, without further assumptions. Given the similarity of the results from quite different models, we conclude that the empirical anisotropy upper limit in real galaxies, and the corresponding observed distributions in the $(\lambda_{\rm R}, \varepsilon)$ and $(V/\sigma, \varepsilon)$ diagrams, are due to the lack of physical axisymmetric equilibrium solutions at high $\beta$ anisotropy when the velocity ellipsoid is close to oblate.


title: Learning Dexterous Manipulation from Suboptimal Experts

abstract: Learning dexterous manipulation in high-dimensional state-action spaces is an important open challenge with exploration presenting a major bottleneck. Although in many cases the learning process could be guided by demonstrations or other suboptimal experts, current RL algorithms for continuous action spaces often fail to effectively utilize combinations of highly off-policy expert data and on-policy exploration data. As a solution, we introduce Relative Entropy Q-Learning (REQ), a simple policy iteration algorithm that combines ideas from successful offline and conventional RL algorithms. It represents the optimal policy via importance sampling from a learned prior and is well-suited to take advantage of mixed data distributions. We demonstrate experimentally that REQ outperforms several strong baselines on robotic manipulation tasks for which suboptimal experts are available. We show how suboptimal experts can be constructed effectively by composing simple waypoint tracking controllers, and we also show how learned primitives can be combined with waypoint controllers to obtain reference behaviors to bootstrap a complex manipulation task on a simulated bimanual robot with human-like hands. Finally, we show that REQ is also effective for general off-policy RL, offline RL, and RL from demonstrations. Videos and further materials are available at sites.google.com/view/rlfse.


title: Reinforcement Learning with Neural Networks for Quantum Multiple Hypothesis Testing

abstract: Reinforcement learning with neural networks (RLNN) has recently demonstrated great promise for many problems, including some problems in quantum information theory. In this work, we apply RLNN to quantum hypothesis testing and determine the optimal measurement strategy for distinguishing between multiple quantum states $\{ \rho_{j} \}$ while minimizing the error probability. In the case where the candidate states correspond to a quantum system with many qubit subsystems, implementing the optimal measurement on the entire system is experimentally infeasible. We use RLNN to find locally-adaptive measurement strategies that are experimentally feasible, where only one quantum subsystem is measured in each round. We provide numerical results which demonstrate that RLNN successfully finds the optimal local approach, even for candidate states up to 20 subsystems. We additionally demonstrate that the RLNN strategy meets or exceeds the success probability for a modified locally greedy approach in each random trial. While the use of RLNN is highly successful for designing adaptive local measurement strategies, in general a significant gap can exist between the success probability of the optimal locally-adaptive measurement strategy and the optimal collective measurement. We build on previous work to provide a set of necessary and sufficient conditions for collective protocols to strictly outperform locally adaptive protocols. We also provide a new example which, to our knowledge, is the simplest known state set exhibiting a significant gap between local and collective protocols. This result raises interesting new questions about the gap between theoretically optimal measurement strategies and practically implementable measurement strategies.


title: THz range Faraday rotation in the Weyl Semimetal Candidate $\mathrm{Co_2TiGe}$

abstract: The $\mathrm{Co_2}$ family of ferromagnetic Heusler alloys have attracted interest due to their fully spin-polarized nature, making them ideal for applications in spintronic devices. More recently, the existence of room temperature time-reversal-breaking Weyl nodes near the Fermi level was predicted and confirmed in these systems. As a result of the presence of these Weyl nodes, these systems possess a non-zero momentum space Berry curvature that can dramatically influence transport properties such as the anomalous Hall effect. One of these candidate compounds is $\mathrm{Co_2 Ti Ge}$. Recently, high quality molecular beam epitaxy-grown thin films of $\mathrm{Co_2 Ti Ge}$ have become available. In this work, we present THz-range measurement of MBE-grown $\mathrm{Co_2 Ti Ge}$ films. We measure the THz-range Faraday rotation, which can be understood as a measure of the anomalous Hall effect. We supplement this work with electronic band structure calculations showing that the principal contribution to the anomalous Hall effect in the this material stems from the Berry curvature of the material. Our work shows that this class of Heusler materials shows promise for Weyl semimetal based spintronics.


title: To Lane or Not to Lane? Comparing On-Road Experiences in Developing and Developed Countries using a New Simulator "RoadBird"

abstract: Even though the traffic systems in developed countries have been analyzed with rigor and operated efficiently, the same does not generally hold for developing countries due to inadequate planning, design, and operations of their transportation systems. Because of inherent differences between internal infrastructures, the systems deployed in developed countries may not be amenable to developing ones. Besides, the traffic systems of developing countries are not well-studied in the literature to the best of our knowledge. For example, it is yet to explore how a developed country's lane-based traffic flow would perform in the context of a developing country, which generally experiences non-lane-based traffic. As such, by using our newly developed traffic simulator 'RoadBird', we investigate outcomes of both lane-based and non-lane-based traffic from the contexts of both developing and developed countries. To do so, we run simulations over real road topologies (extracted from the GIS maps of major cities such as Dhaka, Miami, and Riyadh) considering different scenarios such as lane-based or non-lane-based flows, homogeneous or heterogeneous traffic, with or without pedestrians, etc. We also incorporate different car-following and lane-changing models to mimic traffic behaviors and investigate their performances. While the lane changing dilemma remains an open research question, our experimental evidences indicate: (i) lane-based approaches will not necessarily perform better in the case of currently-adopted non-lane-based scenarios; and (ii) non-lane-based strategies may benefit system performance in lane-based scenarios while having heavy mixed traffic. Nonetheless, we reveal several new insights for on-road experiences both in developing and developed countries.


title: A Conglomerate of Multiple OCR Table Detection and Extraction

abstract: Information representation as tables are compact and concise method that eases searching, indexing, and storage requirements. Extracting and cloning tables from parsable documents is easier and widely used, however industry still faces challenge in detecting and extracting tables from OCR documents or images. This paper proposes an algorithm that detects and extracts multiple tables from OCR document. The algorithm uses a combination of image processing techniques, text recognition and procedural coding to identify distinct tables in same image and map the text to appropriate corresponding cell in dataframe which can be stored as Comma-separated values, Database, Excel and multiple other usable formats.


title: The threshold for the square of a Hamilton cycle

abstract: Resolving a conjecture of K\"uhn and Osthus from 2012, we show that $p= 1/\sqrt{n}$ is the threshold for the random graph $G_{n,p}$ to contain the square of a Hamilton cycle.


title: Deep Submodular Networks for Extractive Data Summarization

abstract: Deep Models are increasingly becoming prevalent in summarization problems (e.g. document, video and images) due to their ability to learn complex feature interactions and representations. However, they do not model characteristics such as diversity, representation, and coverage, which are also very important for summarization tasks. On the other hand, submodular functions naturally model these characteristics because of their diminishing returns property. Most approaches for modelling and learning submodular functions rely on very simple models, such as weighted mixtures of submodular functions. Unfortunately, these models only learn the relative importance of the different submodular functions (such as diversity, representation or importance), but cannot learn more complex feature representations, which are often required for state-of-the-art performance. We propose Deep Submodular Networks (DSN), an end-to-end learning framework that facilitates the learning of more complex features and richer functions, crafted for better modelling of all aspects of summarization. The DSN framework can be used to learn features appropriate for summarization from scratch. We demonstrate the utility of DSNs on both generic and query focused image-collection summarization, and show significant improvement over the state-of-the-art. In particular, we show that DSNs outperform simple mixture models using off the shelf features. Secondly, we also show that just using four submodular functions in a DSN with end-to-end learning performs comparably to the state-of-the-art mixture model with a hand-crafted set of 594 components and outperforms other methods for image collection summarization.


title: Geometric cycles and bounded cohomology for a cocompact lattice in $SL_n(\mathbb R)$

abstract: We show there exists a closed locally symmetric manifold $M$ modeled on $SL_n(\mathbb R)/SO(n)$, and a non-trivial homology class in degree $dim(M)-rank(M)$ represented by a totally geodesic submanifold that contains a circle factor. As a result, the comparison map $c^k:H_b^k(M,\mathbb R)\rightarrow H^k(M,\mathbb R)$ is not surjective in degree $k=dim(M)-rank(M)$. This provides a counterpart to a result of Lafont-Wang which states that $c^k$ is always surjective in degree $k\geq dim(M)-rank(M)+2$.


title: Flow-FL: Data-Driven Federated Learning for Spatio-Temporal Predictions in Multi-Robot Systems

abstract: In this paper, we show how the Federated Learning (FL) framework enables learning collectively from distributed data in connected robot teams. This framework typically works with clients collecting data locally, updating neural network weights of their model, and sending updates to a server for aggregation into a global model. We explore the design space of FL by comparing two variants of this concept. The first variant follows the traditional FL approach in which a server aggregates the local models. In the second variant, that we call Flow-FL, the aggregation process is serverless thanks to the use of a gossip-based shared data structure. In both variants, we use a data-driven mechanism to synchronize the learning process in which robots contribute model updates when they collect sufficient data. We validate our approach with an agent trajectory forecasting problem in a multi-agent setting. Using a centralized implementation as a baseline, we study the effects of staggered online data collection, and variations in data flow, number of participating robots, and time delays introduced by the decentralization of the framework in a multi-robot setting.


title: Efficient Robotic Object Search via HIEM: Hierarchical Policy Learning with Intrinsic-Extrinsic Modeling

abstract: Despite the significant success at enabling robots with autonomous behaviors makes deep reinforcement learning a promising approach for robotic object search task, the deep reinforcement learning approach severely suffers from the nature sparse reward setting of the task. To tackle this challenge, we present a novel policy learning paradigm for the object search task, based on hierarchical and interpretable modeling with an intrinsic-extrinsic reward setting. More specifically, we explore the environment efficiently through a proxy low-level policy which is driven by the intrinsic rewarding sub-goals. We further learn our hierarchical policy from the efficient exploration experience where we optimize both of our high-level and low-level policies towards the extrinsic rewarding goal to perform the object search task well. Experiments conducted on the House3D environment validate and show that the robot, trained with our model, can perform the object search task in a more optimal and interpretable way.


title: Regular spinors and fermionic fields

abstract: Bearing in mind the Lounesto spinor classification, we connect the expansion coefficients of well behaved fermionic quantum field, i.e., a local field within a full Lorentz covariant theory, with and only with a given subclass of Type-2 spinors according to Lounesto. We comment on theoretical possibilities as well as physical outputs for the other cases.


title: Evidence for periodic accretion-ejection in LSI+61303

abstract: The stellar binary system LS I +61303, composed of a compact object in an eccentric orbit around a B0 Ve star, emits from radio up to gamma-ray energies. The orbital modulation of radio spectral index, X-ray, and GeV gamma-ray data suggests the presence of two peaks. This two-peaked profile is in line with the accretion theory predicting two accretion-ejection events for LS I +61303 along the 26.5 d orbit. However, the existing multiwavelength data are not simultaneous. In this paper, we report the results of a campaign covering radio, X-ray, and gamma-ray observations of the system along one single orbit. Our results confirm the two predicted events along the orbit and in addition show that the positions of radio and gamma-ray peaks are coincident with X-ray dips as expected for radio and gamma-ray emitting ejections depleting the X-ray emitting accretion flow. We discuss future observing strategies for a systematic study of the accretion-ejection physical processes in LS I +61303.


title: Logical Relations as Types: Proof-Relevant Parametricity for Program Modules

abstract: The theory of program modules is of interest to language designers not only for its practical importance to programming, but also because it lies at the nexus of three fundamental concerns in language design: the phase distinction, computational effects, and type abstraction. We contribute a fresh "synthetic" take on program modules that treats modules as the fundamental constructs, in which the usual suspects of prior module calculi (kinds, constructors, dynamic programs) are rendered as derived notions in terms of a modal type-theoretic account of the phase distinction. We simplify the account of type abstraction (embodied in the generativity of module functors) through a lax modality that encapsulates computational effects. Our main result is a (significant) proof-relevant and phase-sensitive generalization of the Reynolds abstraction theorem for a calculus of program modules, based on a new kind of logical relation called a parametricity structure. Parametricity structures generalize the proof-irrelevant relations of classical parametricity to proof-relevant families, where there may be non-trivial evidence witnessing the relatedness of two programs -- simplifying the metatheory of strong sums over the collection of types, for although there can be no "relation classifying relations", one easily accommodates a "family classifying small families". Using the insight that logical relations/parametricity is itself a form of phase distinction between the syntactic and the semantic, we contribute a new synthetic approach to phase separated parametricity based on the slogan "logical relations as types", iterating our modal account of the phase distinction. Then, to construct a simulation between two implementations of an abstract type, one simply programs a third implementation whose type component carries the representation invariant.


title: Robot Navigation in Constrained Pedestrian Environments using Reinforcement Learning

abstract: Navigating fluently around pedestrians is a necessary capability for mobile robots deployed in human environments, such as buildings and homes. While research on social navigation has focused mainly on the scalability with the number of pedestrians in open spaces, typical indoor environments present the additional challenge of constrained spaces such as corridors and doorways that limit maneuverability and influence patterns of pedestrian interaction. We present an approach based on reinforcement learning (RL) to learn policies capable of dynamic adaptation to the presence of moving pedestrians while navigating between desired locations in constrained environments. The policy network receives guidance from a motion planner that provides waypoints to follow a globally planned trajectory, whereas RL handles the local interactions. We explore a compositional principle for multi-layout training and find that policies trained in a small set of geometrically simple layouts successfully generalize to more complex unseen layouts that exhibit composition of the structural elements available during training. Going beyond walls-world like domains, we show transfer of the learned policy to unseen 3D reconstructions of two real environments. These results support the applicability of the compositional principle to navigation in real-world buildings and indicate promising usage of multi-agent simulation within reconstructed environments for tasks that involve interaction.


title: Information Coefficient as a Performance Measure of Stock Selection Models

abstract: Information coefficient (IC) is a widely used metric for measuring investment managers' skills in selecting stocks. However, its adequacy and effectiveness for evaluating stock selection models has not been clearly understood, as IC from a realistic stock selection model can hardly be materially different from zero and is often accompanies with high volatility. In this paper, we investigate the behavior of IC as a performance measure of stick selection models. Through simulation and simple statistical modeling, we examine the IC behavior both statically and dynamically. The examination helps us propose two practical procedures that one may use for IC-based ongoing performance monitoring of stock selection models.


title: Importance of mirror modes in binary black hole ringdown waveform

abstract: The post-merger signal in binary black hole merger is described by linear, black-hole perturbation theory. Historically, this has been modeled using the dominant positive-frequency (corotating) fundamental mode. Recently, there has been a renewed effort in modeling the post-merger waveform using higher, positive-frequency overtones in an attempt to achieve greater accuracy in describing the waveform at earlier times using linear perturbation theory. It has been shown that the inclusion of higher overtones can shift the linear regime to the peak of $(l,m)=(2,2)$ spherical harmonic mode. In this work, we show that the inclusion of negative-frequency (counterrotating) modes, called 'mirror' modes, extends the validity of linear perturbation theory to even earlier times, with far lower systematic uncertainties in the model in recovering the remnant parameters at these early times. A good description of the signal at early times also enables for a greater signal-to-noise ratio to be accumulated in the ringdown phase, thereby, allowing for a more accurate measurement of remnant parameters and tests of general relativity.


title: A Sequential Framework Towards an Exact SDP Verification of Neural Networks

abstract: Although neural networks have been applied to several systems in recent years, they still cannot be used in safety-critical systems due to the lack of efficient techniques to certify their robustness. A number of techniques based on convex optimization have been proposed in the literature to study the robustness of neural networks, and the semidefinite programming (SDP) approach has emerged as a leading contender for the robust certification of neural networks. The major challenge to the SDP approach is that it is prone to a large relaxation gap. In this work, we address this issue by developing a sequential framework to shrink this gap to zero by adding non-convex cuts to the optimization problem via disjunctive programming. We analyze the performance of this sequential SDP method both theoretically and empirically, and show that it bridges the gap as the number of cuts increases.


title: Revisiting QCD corrections to the forward-backward charge asymmetry of heavy quarks in electron-positron collisions at the Z pole: really a problem?

abstract: We review in some detail the QCD corrections to the measurement of the forward-backward charge asymmetry of heavy quarks in the $\mathrm{e^+e^-\rightarrow Q\overline{Q}(g)}$ process at the Z pole. We show that the size of these corrections can be reduced by an order of magnitude by using simple cuts on jet acollinearity. Such a reduction is expected to lead to systematic uncertainties at the $\Delta \mathrm{A_{FB}^{0,Q}} \approx 10^{-4}$ level, opening up the path to high precision electroweak measurements with heavy flavors at future high luminosity $\mathrm{e^+e^-}$ colliders like the FCC-ee.


title: Predicting Playa Inundation Using a Long Short-Term Memory Neural Network

abstract: In the Great Plains, playas are critical wetland habitats for migratory birds and a source of recharge for the agriculturally-important High Plains aquifer. The temporary wetlands exhibit complex hydrology, filling rapidly via local rain storms and then drying through evaporation and groundwater infiltration. Using a long short-term memory (LSTM) neural network to account for these complex processes, we modeled playa inundation for 71,842 playas in the Great Plains from 1984-2018. At the level of individual playas, the model achieved an F1-score of 0.538 on a withheld test set, displaying the ability to predict complex inundation patterns. When averaging over all the playas in the entire region, the model is able to very closely track inundation trends, even during periods of drought. Our results demonstrate potential for using LSTMs to model complex hydrological dynamics. Our modeling approach could be used to model playa inundation into the future under different climate scenarios to better understand how wetland habitats and groundwater will be impacted by changing climate.


title: Evaluating Attribution Methods using White-Box LSTMs

abstract: Interpretability methods for neural networks are difficult to evaluate because we do not understand the black-box models typically used to test them. This paper proposes a framework in which interpretability methods are evaluated using manually constructed networks, which we call white-box networks, whose behavior is understood a priori. We evaluate five methods for producing attribution heatmaps by applying them to white-box LSTM classifiers for tasks based on formal languages. Although our white-box classifiers solve their tasks perfectly and transparently, we find that all five attribution methods fail to produce the expected model explanations.


title: DeepIntent: ImplicitIntent based Android IDS with E2E Deep Learning architecture

abstract: The Intent in Android plays an important role in inter-process and intra-process communications. The implicit Intent that an application could accept are declared in its manifest and are amongst the easiest feature to extract from an apk. Implicit Intents could even be extracted online and in real-time. So far neither the feasibility of developing an Intrusion Detection System solely on implicit Intent has been explored, nor are any benchmarks available of a malware classifier that is based on implicit Intent alone. We demonstrate that despite Intent is implicit and well declared, it can provide very intuitive insights to distinguish malicious from non-malicious applications. We conducted exhaustive experiments with over 40 different end-to-end Deep Learning configurations of Auto-Encoders and Multi-Layer-Perceptron to create a benchmark for a malware classifier that works exclusively on implicit Intent. Using the results from the experiments we create an intrusion detection system using only the implicit Intents and end-to-end Deep Learning architecture. We obtained an area-under-curve statistic of 0.81, and accuracy of 77.2% along with false-positive-rate of 0.11 on Drebin dataset.


title: DOOM: A Novel Adversarial-DRL-Based Op-Code Level Metamorphic Malware Obfuscator for the Enhancement of IDS

abstract: We designed and developed DOOM (Adversarial-DRL based Opcode level Obfuscator to generate Metamorphic malware), a novel system that uses adversarial deep reinforcement learning to obfuscate malware at the op-code level for the enhancement of IDS. The ultimate goal of DOOM is not to give a potent weapon in the hands of cyber-attackers, but to create defensive-mechanisms against advanced zero-day attacks. Experimental results indicate that the obfuscated malware created by DOOM could effectively mimic multiple-simultaneous zero-day attacks. To the best of our knowledge, DOOM is the first system that could generate obfuscated malware detailed to individual op-code level. DOOM is also the first-ever system to use efficient continuous action control based deep reinforcement learning in the area of malware generation and defense. Experimental results indicate that over 67% of the metamorphic malware generated by DOOM could easily evade detection from even the most potent IDS. This achievement gains significance, as with this, even IDS augment with advanced routing sub-system can be easily evaded by the malware generated by DOOM.


title: Ab initio benchmarks of neutrinoless double beta decay in light nuclei with a chiral Hamiltonian

abstract: We report ab initio benchmark calculations of nuclear matrix elements (NMEs) for neutrinoless double-beta ($0\nu\beta\beta$) decays in light nuclei with mass number ranging from $A=6$ to $A=22$. We use the transition operator derived from light-Majorana neutrino exchange and evaluate the NME with three different methods: two variants of in-medium similarity renormalization group (IMSRG) and importance-truncated no-core shell model (IT-NCSM). The same two-plus-three-nucleon interaction from chiral effective field theory is employed, and both isospin-conserving ($\Delta T=0$) and isospin-changing ($\Delta T=2$) transitions are studied. We compare our resulting ground-state energies and NMEs to those of recent ab initio no-core shell model and coupled-cluster calculations, also with the same inputs. We show that the NMEs of $\Delta T=0$ transitions are in good agreement among all calculations, at the level of 10%. For $\Delta T=2$, relative deviations are more significant in some nuclei. The comparison with the exact IT-NCSM result allows us to analyze these cases in detail, and indicates the next steps towards improving the IMSRG-based approaches. The present study clearly demonstrates the power of consistent cross-checks that are made possible by ab initio methodology. This capability is crucial for providing meaningful many-body uncertainties in the NMEs for the $0\nu\beta\beta$ decays in heavier candidate nuclei, where quasi-exact benchmarks are not available.


title: Szeg\H{o} and Widom theorems for finite codimensional subalgebras of a class of uniform algebras

abstract: We establish versions of Szeg\H{o}'s distance formula and Widom's theorem on invertibility of (a family of) Toeplitz operators in a class of finite codimension subalgebras of uniform algebras, obtained by imposing a finite number of linear constraints. Each such algebra is naturally represented on a family of reproducing kernel Hilbert spaces, which play a central role in the proofs.


title: Minimal enumeration of all possible total effects in a Markov equivalence class

abstract: In observational studies, when a total causal effect of interest is not identified, the set of all possible effects can be reported instead. This typically occurs when the underlying causal DAG is only known up to a Markov equivalence class, or a refinement thereof due to background knowledge. As such, the class of possible causal DAGs is represented by a maximally oriented partially directed acyclic graph (MPDAG), which contains both directed and undirected edges. We characterize the minimal additional edge orientations required to identify a given total effect. A recursive algorithm is then developed to enumerate subclasses of DAGs, such that the total effect in each subclass is identified as a distinct functional of the observed distribution. This resolves an issue with existing methods, which often report possible total effects with duplicates, namely those that are numerically distinct due to sampling variability but are in fact causally identical.


title: Understanding the Role of Intermediaries in Online Social E-commerce: An Exploratory Study of Beidian

abstract: Social e-commerce, as a new form of social computing based marketing platforms, utilizes existing real-world social relationships for promotions and sales of products. It has been growing rapidly in recent years and attracted tens of millions of users in China. A key group of actors who enable market transactions on these platforms are intermediaries who connect producers with consumers by sharing information with and recommending products to their real-world social contacts. Despite their crucial role, the nature and behavior of these intermediaries on these social e-commerce platforms has not been systematically analyzed. Here we address this knowledge gap through a mixed method study. Leveraging 9 months' all-round behavior of about 40 million users on Beidian -- one of the largest social e-commerce sites in China, alongside with qualitative evidence from online forums and interviews, we examine characteristics of intermediaries, identify their behavioral patterns and uncover strategies and mechanisms that make successful intermediaries. We demonstrate that intermediaries on social e-commerce sites act as local trend detectors and "social grocers". Furthermore, successful intermediaries are highly dedicated whenever best sellers appear and broaden items for promotion. To the best of our knowledge, this paper presents the first large-scale analysis on the emerging role of intermediaries in social e-commerce platforms, which provides potential insights for the design and management of social computing marketing platforms.


title: The Horton-Strahler Number of Conditioned Galton-Watson Trees

abstract: The Horton-Strahler number of a tree is a measure of its branching complexity; it is also known in the literature as the register function. We show that for critical Galton-Watson trees with finite variance conditioned to be of size $n$, the Horton-Strahler number grows as $\frac{1}{2}\log_2 n$ in probability. We further define some generalizations of this number. Among these are the rigid Horton-Strahler number and the $k$-ary register function, for which we prove asymptotic results analogous to the standard case.


title: Finite-difference-based simulation and adjoint optimization of gas networks

abstract: The stable operation of gas networks is an important optimization target. While for this task commonly finite volume methods are used, we introduce a new finite difference approach. With a summation by part formulation for the spatial discretization, we get well-defined fluxes between the pipes. This allows a simple and explicit formulation of the coupling conditions at the node. From that, we derive the adjoint equations for the network simply and transparently. The resulting direct and adjoint equations are numerically efficient and easy to implement. The approach is demonstrated by the optimization of two sample gas networks.


title: Decomposability and Parallel Computation of Multi-Agent LQR

abstract: Individual agents in a multi-agent system (MAS) may have decoupled open-loop dynamics, but a cooperative control objective usually results in coupled closed-loop dynamics thereby making the control design computationally expensive. The computation time becomes even higher when a learning strategy such as reinforcement learning (RL) needs to be applied to deal with the situation when the agents dynamics are not known. To resolve this problem, we propose a parallel RL scheme for a linear quadratic regulator (LQR) design in a continuous-time linear MAS. The idea is to exploit the structural properties of two graphs embedded in the $Q$ and $R$ weighting matrices in the LQR objective to define an orthogonal transformation that can convert the original LQR design to multiple decoupled smaller-sized LQR designs. We show that if the MAS is homogeneous then this decomposition retains closed-loop optimality. Conditions for decomposability, an algorithm for constructing the transformation matrix, a parallel RL algorithm, and robustness analysis when the design is applied to non-homogeneous MAS are presented. Simulations show that the proposed approach can guarantee significant speed-up in learning without any loss in the cumulative value of the LQR cost.