title: Covariance-Aware Private Mean Estimation Without Private Covariance Estimation

abstract: We present two sample-efficient differentially private mean estimators for $d$-dimensional (sub)Gaussian distributions with unknown covariance. Informally, given $n \gtrsim d/\alpha^2$ samples from such a distribution with mean $\mu$ and covariance $\Sigma$, our estimators output $\tilde\mu$ such that $\| \tilde\mu - \mu \|_{\Sigma} \leq \alpha$, where $\| \cdot \|_{\Sigma}$ is the Mahalanobis distance. All previous estimators with the same guarantee either require strong a priori bounds on the covariance matrix or require $\Omega(d^{3/2})$ samples. Each of our estimators is based on a simple, general approach to designing differentially private mechanisms, but with novel technical steps to make the estimator private and sample-efficient. Our first estimator samples a point with approximately maximum Tukey depth using the exponential mechanism, but restricted to the set of points of large Tukey depth. Its accuracy guarantees hold even for data sets that have a small amount of adversarial corruption. Proving that this mechanism is private requires a novel analysis. Our second estimator perturbs the empirical mean of the data set with noise calibrated to the empirical covariance, without releasing the covariance itself. Its sample complexity guarantees hold more generally for subgaussian distributions, albeit with a slightly worse dependence on the privacy parameter. For both estimators, careful preprocessing of the data is required to satisfy differential privacy.


title: Borel combinatorics fail in HYP

abstract: We characterize the completely determined Borel subsets of HYP as exactly the omega_1^{ck} subsets of HYP. As a result, HYP believes there is a Borel well-ordering of the reals, that the Borel Dual Ramsey Theorem fails, and that every Borel d-regular bipartite graph has a Borel perfect matching, among other examples. Therefore, the Borel Dual Ramsey Theorem and several theorems of descriptive combinatorics are not theories of hyperarithmetic analysis. In the case of the Borel Dual Ramsey Theorem, this answers a question of Astor, Dzhafarov, Montalban, Solomon & the third author.


title: Linear Multifractional Stable Sheets in the Broad Sense: Existence and Joint Continuity of Local Times

abstract: We introduce the notion of linear multifractional stable sheets in the broad sense (LMSS) with $\alpha\in(0,2]$, to include both linear multifractional Brownian sheets ($\alpha=2$) and linear multifractional stable sheets ($\alpha<2$). The purpose of the present paper is to study the existence and joint continuity of the local times of LMSS, and also the local H\"older condition of the local times in the set variable. Among the main results of this paper, Theorem 2.4 provides a sufficient and necessary condition for the existence of local times of LMSS; Theorem 3.1 shows a sufficient condition for the joint continuity of local times; and Theorem 4.1 proves a sharp local H\"older condition for the local times in the set variable. All these theorems improve significantly the existing results for the local times of multifractional Brownian sheets and linear multifractional stable sheets in the literature.


title: Imaging arbitrary incoherent source distributions with near quantum-limited resolution

abstract: We demonstrate an approach to obtaining near quantum-limited far-field imaging resolution of incoherent sources with arbitrary distributions. Our method assumes no prior knowledge of the source distribution, but rather uses an adaptive approach to imaging via spatial mode demultiplexing that iteratively updates both the form of the spatial imaging modes and the estimate of the source distribution. The optimal imaging modes are determined by minimizing the estimated Cram\'er-Rao bound over the manifold of all possible sets of orthogonal imaging modes. We have observed through Monte Carlo simulations that the manifold-optimized spatial mode demultiplexing measurement consistently outperforms standard imaging techniques in the accuracy of source reconstructions and comes within a factor of 2 of the absolute quantum limit as set by the quantum Cram\'er-Rao bound. The adaptive framework presented here allows for a consistent approach to achieving near quantum-limited imaging resolution of arbitrarily distributed sources through spatial mode imaging techniques.


title: Large-data equicontinuity for the derivative NLS

abstract: We consider the derivative NLS equation in one spatial dimension, which is known to be completely integrable. We prove that the orbits of $L^2$ bounded and equicontinuous sets of initial data remain bounded and equicontinuous, not only under this flow, but under the entire hierarchy. This allows us to remove the small-data restriction from prior conservation laws and global well-posedness results.


title: Advanced calculations of x-ray spectroscopies with FEFF10 and Corvus

abstract: The real-space Green's function code FEFF has been extensively developed and used for calculations of x-ray and related spectra, including x-ray absorption (XAS), x-ray emission (XES), inelastic x-ray scattering, and electron energy loss spectra (EELS). The code is particularly useful for the analysis and interpretation of the XAS fine-structure (EXAFS) and the near-edge structure (XANES) in materials throughout the periodic table. Nevertheless, many applications, such as non-equilibrium systems, and the analysis of ultra-fast pump-probe experiments, require extensions of the code including finite-temperature and auxiliary calculations of structure and vibrational properties. To enable these extensions, we have developed in tandem, a new version FEFF10, and new FEFF based workflows for the Corvus workflow manager, which allow users to easily augment the capabilities of FEFF10 via auxiliary codes. This coupling facilitates simplified input and automated calculations of spectra based on advanced theoretical techniques. The approach is illustrated with examples of high temperature behavior, vibrational properties, many-body excitations in XAS, super-heavy materials, and fits of calculated spectra to experiment.


title: A False Vacuum Skyrme Model for Nuclear Matter

abstract: The low energy regime of Quantum Chromodynamics (QCD) presents enormous challenges due to its large coupling. Effective field theories, like the Skyrme model, are useful approaches to study properties of strong interaction at hadronic scales. We propose a Skyrme-type model with a self-dual sector and that treats the density of the baryonic charge as a self-interacting fluid. The dynamics reduces to Coleman's false vacuum problem for a scalar field that is a fractional power of that density. The main result is that such a Skyrme-type model is the first one to reproduce, with good accuracy, the experimental values of radii and binding energies for a very wide range of the mass number. The robust and simple properties of the model lead to many possible generalizations with implications not only in nuclear physics but also in other areas of Physics.


title: Persistent chlorophyll maxima in the Eastern Tropical North Pacific

abstract: This study aims to describe the response of two persistent chlorophyll-a maxima to physical processes affecting the thermocline and nitracline position in the Eastern Tropical North Pacific (ETNP). We focused on Long Rossby Waves given their relevance to the ETNP circulation and its potential as a mechanism introducing nutrients into the euphotic zone. We found the shallower chlorophyll-a maximum in oxygenated waters became more intense when denser waters (more nutrients) moved toward the surface. It suggests that isopycnals and nitracline displacements modify the nutrient supply in the euphotic zone, which produces changes in phytoplankton growth. The suboxic and deeper chlorophyll a maximum showed a strong association with the 26 kg m^{-3} isopycnal being only mechanically displaced, and its chlorophyll-a content does not seem to covary with irradiance or nutrients. The different responses of the chlorophyll-a maxima could be explained if different phytoplankton groups are associated with them. Long Rossby Waves can affect the position of the thermocline, nitracline and isopycnals in an annual cycle, but it seems to be a background signal modulated by higher frequency processes such as mesoscale eddies and other Rossby waves. The co-occurrence of processes can control the nitracline depth, and hence the input of nutrients into the euphotic zone that can cause sporadic enhancements of the chlorophyll-a concentration of one maximum.


title: On Zero Point Calibration of the Red Giant Branch Tip in the Magellanic Clouds

abstract: A zero point calibration of the Red Giant Branch Tip (TRGB) in the $I$-band is determined from OGLE photometry of the Magellanic Clouds (MCs). It is shown that TRGB measurements made in star-forming regions, with concomitantly high quantities of gas and dust, are less precise and biased to fainter magnitudes, as compared to the same measurements made in quiescent regions. Once these low accuracy fields are excluded from consideration, the TRGB can be used for the first time to constrain the three-dimensional plane geometry of the LMC. Composite CMDs are constructed for the SMC and LMC from only those fields with well-defined TRGB features, and the highest accuracy TRGB zero point calibration to date is presented. The $I$-band TRGB magnitude is measured to be flat over the color range $ 1.45 < (V-I)_0 < 1.95$ mag, with a modest slope introduced when including metal-rich (up to $(V-I)_0 = 2.2$ mag) Tip stars into the fit. Both the flat, blue zero point and the shallow slope calibration are consistent with the canonical value of $-4.05$ mag for the old, metal-poor TRGB, and would appear to resolve a recent debate in the literature over the method's absolute calibration.


title: Distributed IDA-PBC for a Class of Nonholonomic Mechanical Systems

abstract: Nonholonomic mechanical systems encompass a large class of practically interesting robotic structures, such as wheeled mobile robots, space manipulators, and multi-fingered robot hands. However, few results exist on the cooperative control of such systems in a generic, distributed approach. In this work we extend a recently developed distributed Interconnection and Damping Assignment Passivity-Based Control (IDA-PBC) method to such systems. More specifically, relying on port-Hamiltonian system modelling for networks of mechanical systems, we propose a full-state stabilization control law for a class of nonholonomic systems within the framework of distributed IDA-PBC. This enables the cooperative control of heterogeneous, underactuated and nonholonomic systems with a unified control law. This control law primarily relies on the notion of Passive Configuration Decomposition (PCD) and a novel, non-smooth desired potential energy function proposed here. A low-level collision avoidance protocol is also implemented in order to achieve dynamic inter-agent collision avoidance, enhancing the practical relevance of this work. Theoretical results are tested in different simulation scenarios in order to highlight the applicability of the derived method.


title: Blockchain-based Security Framework for Critical Industry 4.0 Cyber-physical System

abstract: There has been an intense concern for security alternatives because of the recent rise of cyber attacks, mainly targeting critical systems such as industry, medical, or energy ecosystem. Though the latest industry infrastructures largely depend on AI-driven maintenance, the prediction based on corrupted data undoubtedly results in loss of life and capital. Admittedly, an inadequate data-protection mechanism can readily challenge the security and reliability of the network. The shortcomings of the conventional cloud or trusted certificate-driven techniques have motivated us to exhibit a unique Blockchain-based framework for a secure and efficient industry 4.0 system. The demonstrated framework obviates the long-established certificate authority after enhancing the consortium Blockchain that reduces the data processing delay, and increases cost-effective throughput. Nonetheless, the distributed industry 4.0 security model entails cooperative trust than depending on a single party, which in essence indulges the costs and threat of the single point of failure. Therefore, multi-signature technique of the proposed framework accomplishes the multi-party authentication, which confirms its applicability for the real-time and collaborative cyber-physical system.


title: Subgradient Ellipsoid Method for Nonsmooth Convex Problems

abstract: In this paper, we present a new ellipsoid-type algorithm for solving nonsmooth problems with convex structure. Examples of such problems include nonsmooth convex minimization problems, convex-concave saddle-point problems and variational inequalities with monotone operator. Our algorithm can be seen as a combination of the standard Subgradient and Ellipsoid methods. However, in contrast to the latter one, the proposed method has a reasonable convergence rate even when the dimensionality of the problem is sufficiently large. For generating accuracy certificates in our algorithm, we propose an efficient technique, which ameliorates the previously known recipes.


title: Guessing Based on Compressed Side Information

abstract: A source sequence is to be guessed with some fidelity based on a rate-limited description of an observed sequence with which it is correlated. The trade-off between the description rate and the exponential growth rate of the least power mean of the number of guesses is characterized.


title: The Complexity of Boolean Conjunctive Queries with Intersection Joins

abstract: Intersection joins over interval data are relevant in spatial and temporal data settings. A set of intervals join if their intersection is non-empty. In case of point intervals, the intersection join becomes the standard equality join. We establish the complexity of Boolean conjunctive queries with intersection joins by a many-one equivalence to disjunctions of Boolean conjunctive queries with equality joins. The complexity of any query with intersection joins is that of the hardest query with equality joins in the disjunction exhibited by our equivalence. This is captured by a new width measure called the IJ-width. We also introduce a new syntactic notion of acyclicity called iota-acyclicity to characterise the class of Boolean queries with intersection joins that admit linear time computation modulo a poly-logarithmic factor in the data size. Iota-acyclicity is for intersection joins what alpha-acyclicity is for equality joins. It strictly sits between gamma-acyclicity and Berge-acyclicity. The intersection join queries that are not iota-acyclic are at least as hard as the Boolean triangle query with equality joins, which is widely considered not computable in linear time.


title: Second-order estimates for collapsed limits of Ricci-flat K\"ahler metrics

abstract: We show that the singularities of the twisted K\"ahler--Einstein metric arising as the long-time solution of the K\"ahler--Ricci flow or in the collapsed limit of Ricci-flat K\"ahler metrics is intimately related to the holomorphic sectional curvature of the reference conical geometry. This provides an alternative proof of the second-order estimate obtained by Gross--Tosatti--Zhang with explicit constants appearing in the divisorial pole.


title: First light observations of the solar wind in the outer corona with the Metis coronagraph

abstract: The investigation of the wind in the solar corona initiated with the observations of the resonantly scattered UV emission of the coronal plasma obtained with UVCS-SOHO, designed to measure the wind outflow speed by applying the Doppler dimming diagnostics. Metis on Solar Orbiter complements the UVCS spectroscopic observations, performed during solar activity cycle 23, by simultaneously imaging the polarized visible light and the HI Ly-alpha corona in order to obtain high-spatial and temporal resolution maps of the outward velocity of the continuously expanding solar atmosphere. The Metis observations, on May 15, 2020, provide the first HI Ly-alpha images of the extended corona and the first instantaneous map of the speed of the coronal plasma outflows during the minimum of solar activity and allow us to identify the layer where the slow wind flow is observed. The polarized visible light (580-640 nm), and the UV HI Ly-alpha (121.6 nm) coronal emissions, obtained with the two Metis channels, are combined in order to measure the dimming of the UV emission relative to a static corona. This effect is caused by the outward motion of the coronal plasma along the direction of incidence of the chromospheric photons on the coronal neutral hydrogen. The plasma outflow velocity is then derived as a function of the measured Doppler dimming. The static corona UV emission is simulated on the basis of the plasma electron density inferred from the polarized visible light. This study leads to the identification, in the velocity maps of the solar corona, of the high-density layer about +/-10 deg wide, centered on the extension of a quiet equatorial streamer present at the East limb where the slowest wind flows at about (160 +/- 18) km/s from 4 Rs to 6 Rs. Beyond the boundaries of the high-density layer, the wind velocity rapidly increases, marking the transition between slow and fast wind in the corona.


title: The Hanson-Wright Inequality for Random Tensors

abstract: We provide moment bounds for expressions of the type $(X^{(1)} \otimes \dots \otimes X^{(d)})^T A (X^{(1)} \otimes \dots \otimes X^{(d)})$ where $\otimes$ denotes the Kronecker product and $X^{(1)}, \dots, X^{(d)}$ are random vectors with independent, mean 0, variance 1, subgaussian entries. The bounds are tight up to constants depending on $d$ for the case of Gaussian random vectors. Our proof also provides a decoupling inequality for expressions of this type. Using these bounds, we obtain new, improved concentration inequalities for expressions of the form $\|B (X^{(1)} \otimes \dots \otimes X^{(d)})\|_2$.


title: What will it take to generate fairness-preserving explanations?

abstract: In situations where explanations of black-box models may be useful, the fairness of the black-box is also often a relevant concern. However, the link between the fairness of the black-box model and the behavior of explanations for the black-box is unclear. We focus on explanations applied to tabular datasets, suggesting that explanations do not necessarily preserve the fairness properties of the black-box algorithm. In other words, explanation algorithms can ignore or obscure critical relevant properties, creating incorrect or misleading explanations. More broadly, we propose future research directions for evaluating and generating explanations such that they are informative and relevant from a fairness perspective.


title: Relationship between Cultural Values, Sense of Community and Trust and the Effect of Trust in Workplace

abstract: This paper provides a general overview of different perspectives and studies on trust, offers a definition of trust, and provides factors that play a substantial role in developing social trust, and shows from which perspectives it can be fostered. The results showed that trust is playing an important role in success for organizations involved in cross-national strategic partnerships. Trust can reduce transaction costs, promotes inter-organizational relationships, and improve subordinate relationships between managers.


title: Implications of a search for intergalactic civilizations on prior estimates of human survival and travel speed

abstract: We present a model where some proportion of extraterrestrial civilizations expand uniformly over time to reach a cosmological scale. We then ask what humanity could infer if a sky survey were to find zero, one, or more such civilizations. We show how the results of this survey, in combination with an approach to anthropics called the Self Indication Assumption (SIA), would shift any prior estimates of two quantities: 1) The chance that a technological civilization like ours survives to embark on such expansion, and 2) the maximum feasible speed at which it could expand. The SIA gives pessimistic estimates for both, but survey results (even null results) can reverse some of the effect. The SIA-induced shift in expectations is strong enough to be regarded as a falsifiable technological prediction, and a test of the SIA in cosmological reasoning.


title: Johnson-Lindenstrauss Embeddings with Kronecker Structure

abstract: We prove the Johnson-Lindenstrauss property for matrices $\Phi D_\xi$ where $\Phi$ has the restricted isometry property and $D_\xi$ is a diagonal matrix containing the entries of a Kronecker product $\xi = \xi^{(1)} \otimes \dots \otimes \xi^{(d)}$ of $d$ independent Rademacher vectors. Such embeddings have been proposed in recent works for a number of applications concerning compression of tensor structured data, including the oblivious sketching procedure by Ahle et al. for approximate tensor computations. For preserving the norms of $p$ points simultaneously, our result requires $\Phi$ to have the restricted isometry property for sparsity $C(d) (\log p)^d$. In the case of subsampled Hadamard matrices, this can improve the dependence of the embedding dimension on $p$ to $(\log p)^d$ while the best previously known result required $(\log p)^{d + 1}$. That is, for the case of $d=2$ at the core of the oblivious sketching procedure by Ahle et al., the scaling improves from cubic to quadratic. We provide a counterexample to prove that the scaling established in our result is optimal under mild assumptions.


title: Reference Dependence and Random Attention

abstract: We explore the ways that a reference point may direct attention. Utilizing a stochastic choice framework, we provide behavioral foundations for the Reference-Dependent Random Attention Model (RD-RAM). Our characterization result shows that preferences may be uniquely identified even when the attention process depends arbitrarily on both the menu and the reference point. The RD-RAM is able to capture rich behavioral patterns, including frequency reversals among non-status quo alternatives and choice overload. We also analyze specific attention processes, characterizing reference-dependent versions of several prominent models of stochastic consideration.


title: Yield stress and other flow and wall slip parameters of viscoplastic fluids from steady torsional flow

abstract: The ubiquitous wall slip behavior of viscoplastic fluids renders the analysis of their steady torsional flow data to determine their yield stress and other parameters of their shear viscosity material function challenging. Roughened surfaces to prevent wall slip can be used but such methods frequently result in the fracturing of the sample. Instead, it is shown here that the yield stress of a viscoplastic fluid can be determined from the torque versus apparent shear rate data of steady torsional flow using a single gap. The step change in the slope of the torque versus apparent shear rate data of a viscoplastic fluid marks the yield stress of the fluid. For demonstration, the earlier-characterized parameters of the shear viscosity and apparent wall slip behavior of a viscoplastic suspension [He et al., J. Rheology, 63, 19 (2019)] were used for obtaining the shear stress distribution and the torque at each rotational speed. The yield stress, determined from the torque versus apparent shear rate data, could be verified by using a second method relying on the comparison of the wall slip velocities with wall velocities, while at the same time enabling the determination of the other parameters of the shear viscosity material function of the viscoplastic fluid.


title: Local drag of a slender rod parallel to a plane wall in a viscous fluid

abstract: The viscous drag on a slender rod by a wall is important to many biological and industrial systems. This drag critically depends on the separation between the rod and the wall and can be approximated asymptotically in specific regimes, namely far from, or very close to, the wall, but is typically determined numerically for general separations. In this note we determine an asymptotic representation of the local drag for a slender rod parallel to a wall which is valid for all separations. This is possible through matching the behaviour of a rod close to the wall and a rod far from the wall. We show that the leading order drag in both these regimes has been known since 1981 and that they can used to produce a composite representation of the drag which is valid for all separations. This is in contrast to a sphere above a wall, where no simple uniformly valid representation exists. We estimate the error on this composite representation as the separation increases, discuss how the results could be used as resistive-force theory and demonstrate their use on a two-hinged swimmer above a wall.


title: Cutting Down on Prompts and Parameters: Simple Few-Shot Learning with Language Models

abstract: Prompting language models (LMs) with training examples and task descriptions has been seen as critical to recent successes in few-shot learning. In this work, we show that finetuning LMs in the few-shot setting can considerably reduce the need for prompt engineering. In fact, one can use null prompts, prompts that contain neither task-specific templates nor training examples, and achieve competitive accuracy to manually-tuned prompts across a wide range of tasks. While finetuning LMs does introduce new parameters for each downstream task, we show that this memory overhead can be substantially reduced: finetuning only the bias terms can achieve comparable or better accuracy than standard finetuning while only updating 0.1% of the parameters. All in all, we recommend finetuning LMs for few-shot learning as it is more accurate, robust to different prompts, and can be made nearly as efficient as using frozen LMs.


title: Speckle Observations of TESS Exoplanet Host Stars. II. Stellar Companions at 1-1000 AU and Implications for Small Planet Detection

abstract: We present high angular resolution imaging observations of 517 host stars of TESS exoplanet candidates using the `Alopeke and Zorro speckle cameras at Gemini North and South. The sample consists mainly of bright F, G, K stars at distances of less than 500 pc. Our speckle observations span angular resolutions of ~20 mas out to 1.2 arcsec, yielding spatial resolutions of <10 to 500 AU for most stars, and our contrast limits can detect companion stars 5-9 magnitudes fainter than the primary at optical wavelengths. We detect 102 close stellar companions and determine the separation, magnitude difference, mass ratio, and estimated orbital period for each system. Our observations of exoplanet host star binaries reveal that they have wider separations than field binaries, with a mean orbital semi-major axis near 100 AU. Other imaging studies have suggested this dearth of very closely separated binaries in systems which host exoplanets, but incompleteness at small separations makes it difficult to disentangle unobserved companions from a true lack of companions. With our improved angular resolution and sensitivity, we confirm that this lack of close exoplanet host binaries is indeed real. We also search for a correlation between planetary orbital radii vs. binary star separation, but given the very short orbital periods of the TESS planets, we do not find any clear trend. We do note that in exoplanet systems containing binary host stars, there is an observational bias against detecting Earth-size planet transits due to transit depth dilution caused by the companion star.


title: Cohomology ring of tree braid groups and exterior face rings

abstract: For a tree $T$ and a positive integer $n$, let $B_nT$ denote the $n$-strand braid group on $T$. We use discrete Morse theory techniques to show that the cohomology ring $H^*(B_nT)$ is encoded by an explicit abstract simplicial complex $K_nT$ that measures $n$-local interactions among essential vertices of $T$. We show that, in many cases (for instance when $T$ is a binary tree), $H^*(B_nT)$ is the exterior face ring determined by $K_nT$.


title: Monte Carlo Modeling of Spin-polarized Photoemission from p-doped GaAs Activated to Negative Electron Affinity

abstract: The anticorrelation between quantum efficiency (QE) and electron spin polarization (ESP) from a p-doped GaAs activated to negative electron affinity (NEA) is studied in detail using an ensemble Monte Carlo approach. The photoabsorption, momentum and spin relaxation during transport, and tunnelling of electrons through the surface potential barrier are modeled to identify fundamental mechanisms, which limit the efficiency of GaAs spin-polarized electron sources. In particular, we study the response of QE and ESP to various parameters such as the photoexcitation energy, doping density, and electron affinity level. Our modeling results for various transport and emission characteristics are in a good agreement with available experimental data. Our findings show that the behaviour of both QE and ESP at room temperature can be fully explained by the bulk relaxation mechanisms and the time which electrons spend in the material before being emitted.


title: On the In-Out-Proper Orientations of Graphs

abstract: An orientation of a graph $G$ is {\it in-out-proper} if any two adjacent vertices have different in-out-degrees, where the in-out-degree of each vertex is equal to the in-degree minus the out-degree of that vertex. The {\it in-out-proper orientation number} of a graph $G$, denoted by $\overleftrightarrow{\chi}(G)$, is $ \min_{D\in \Gamma}\max_{v\in V(G)} |d_D^{\pm}(v)|$, where $\Gamma$ is the set of in-out-proper orientations of $G$ and $d_D^{\pm}(v)$ is the in-out-degree of the vertex $v$ in the orientation $D$.


title: Scalable Perception-Action-Communication Loops with Convolutional and Graph Neural Networks

abstract: In this paper, we present a perception-action-communication loop design using Vision-based Graph Aggregation and Inference (VGAI). This multi-agent decentralized learning-to-control framework maps raw visual observations to agent actions, aided by local communication among neighboring agents. Our framework is implemented by a cascade of a convolutional and a graph neural network (CNN / GNN), addressing agent-level visual perception and feature learning, as well as swarm-level communication, local information aggregation and agent action inference, respectively. By jointly training the CNN and GNN, image features and communication messages are learned in conjunction to better address the specific task. We use imitation learning to train the VGAI controller in an offline phase, relying on a centralized expert controller. This results in a learned VGAI controller that can be deployed in a distributed manner for online execution. Additionally, the controller exhibits good scaling properties, with training in smaller teams and application in larger teams. Through a multi-agent flocking application, we demonstrate that VGAI yields performance comparable to or better than other decentralized controllers, using only the visual input modality and without accessing precise location or motion state information.


title: A numerically stable online implementation and exploration of WAIC through variations of the predictive density, using NIMBLE

abstract: We go through the process of crafting a robust and numerically stable online algorithm for the computation of the Watanabe-Akaike information criteria (WAIC). We implement this algorithm in the NIMBLE software. The implementation is performed in an online manner and does not require the storage in memory of the complete samples from the posterior distribution. This algorithm allows the user to specify a specific form of the predictive density to be used in the computation of WAIC, in order to cater to specific prediction goals. We then comment and explore via simulations the use of different forms of the predictive density in the context of different predictive goals. We find that when using marginalized predictive densities, WAIC is sensitive to the grouping of the observations into a joint density.


title: An atomic compass -- detecting 3D magnetic field alignment with vector vortex light

abstract: We describe and demonstrate how 3D magnetic field alignment can be inferred from single absorption images of an atomic cloud. While optically pumped magnetometers conventionally rely on temporal measurement of the Larmor precession of atomic dipoles, here a cold atomic vapour provides a spatial interface between vector light and external magnetic fields. Using a vector vortex beam, we inscribe structured atomic spin polarisation in a cloud of cold rubidium atoms, and record images of the resulting absorption patterns. The polar angle of an external magnetic field can be deduced with spatial Fourier analysis. This effect presents an alternative concept for detecting magnetic vector fields, and demonstrates, more generally, how introducing spatial phases between atomic energy levels can translate transient effects to the spatial domain.


title: Multifidelity Modeling for Physics-Informed Neural Networks (PINNs)

abstract: Multifidelity simulation methodologies are often used in an attempt to judiciously combine low-fidelity and high-fidelity simulation results in an accuracy-increasing, cost-saving way. Candidates for this approach are simulation methodologies for which there are fidelity differences connected with significant computational cost differences. Physics-informed Neural Networks (PINNs) are candidates for these types of approaches due to the significant difference in training times required when different fidelities (expressed in terms of architecture width and depth as well as optimization criteria) are employed. In this paper, we propose a particular multifidelity approach applied to PINNs that exploits low-rank structure. We demonstrate that width, depth, and optimization criteria can be used as parameters related to model fidelity, and show numerical justification of cost differences in training due to fidelity parameter choices. We test our multifidelity scheme on various canonical forward PDE models that have been presented in the emerging PINNs literature.


title: Magnetic Field Reconstruction for a Realistic Multi-Point, Multi-Scale Spacecraft Observatory

abstract: Future in situ space plasma investigations will likely involve spatially distributed observatories comprised of multiple spacecraft, beyond the four and five spacecraft configurations currently in operation. Inferring the magnetic field structure across the observatory, and not simply at the observation points, is a necessary step towards characterizing fundamental plasma processes using these unique multi-point, multi-scale data sets. We propose improvements upon the classic first-order reconstruction method, as well as a second-order method, utilizing magnetometer measurements from a realistic nine-spacecraft observatory. The improved first-order method, which averages over select ensembles of four spacecraft, reconstructs the magnetic field associated with simple current sheets and numerical simulations of turbulence accurately over larger volumes compared to second-order methods or first-order methods using a single regular tetrahedron. Using this averaging method on data sets with fewer than nine measurement points, the volume of accurate reconstruction compared to a known magnetic vector field improves approximately linearly with the number of measurement points.


title: Hardy's inequality and (almost) the Landau equation

abstract: In this manuscript we establish an $L^\infty$ estimate for the isotropic analogue of the homogeneous Landau equation. This is done for values of the interaction exponent $\gamma$ in (a part of) the range of very soft potentials. The main observation in our proof is that the classical weighted Hardy inequality leads to a weighted Poincar\'e inequality, which in turn implies the propagation of some $L^p$ norms of solutions. From here, the $L^\infty$ estimate follows from certain weighted Sobolev inequalities and De Giorgi-Nash-Moser theory.


title: CausalCity: Complex Simulations with Agency for Causal Discovery and Reasoning

abstract: The ability to perform causal and counterfactual reasoning are central properties of human intelligence. Decision-making systems that can perform these types of reasoning have the potential to be more generalizable and interpretable. Simulations have helped advance the state-of-the-art in this domain, by providing the ability to systematically vary parameters (e.g., confounders) and generate examples of the outcomes in the case of counterfactual scenarios. However, simulating complex temporal causal events in multi-agent scenarios, such as those that exist in driving and vehicle navigation, is challenging. To help address this, we present a high-fidelity simulation environment that is designed for developing algorithms for causal discovery and counterfactual reasoning in the safety-critical context. A core component of our work is to introduce \textit{agency}, such that it is simple to define and create complex scenarios using high-level definitions. The vehicles then operate with agency to complete these objectives, meaning low-level behaviors need only be controlled if necessary. We perform experiments with three state-of-the-art methods to create baselines and highlight the affordances of this environment. Finally, we highlight challenges and opportunities for future work.


title: RSN: Range Sparse Net for Efficient, Accurate LiDAR 3D Object Detection

abstract: The detection of 3D objects from LiDAR data is a critical component in most autonomous driving systems. Safe, high speed driving needs larger detection ranges, which are enabled by new LiDARs. These larger detection ranges require more efficient and accurate detection models. Towards this goal, we propose Range Sparse Net (RSN), a simple, efficient, and accurate 3D object detector in order to tackle real time 3D object detection in this extended detection regime. RSN predicts foreground points from range images and applies sparse convolutions on the selected foreground points to detect objects. The lightweight 2D convolutions on dense range images results in significantly fewer selected foreground points, thus enabling the later sparse convolutions in RSN to efficiently operate. Combining features from the range image further enhance detection accuracy. RSN runs at more than 60 frames per second on a 150m x 150m detection region on Waymo Open Dataset (WOD) while being more accurate than previously published detectors. As of 11/2020, RSN is ranked first in the WOD leaderboard based on the APH/LEVEL 1 metrics for LiDAR-based pedestrian and vehicle detection, while being several times faster than alternatives.


title: PHANGS-HST: New Methods for Star Cluster Identification in Nearby Galaxies

abstract: We present an innovative and widely applicable approach for the detection and classification of stellar clusters, developed for the PHANGS-HST Treasury Program, an $NUV$-to-$I$ band imaging campaign of 38 spiral galaxies. Our pipeline first generates a unified master source list for stars and candidate clusters, to enable a self-consistent inventory of all star formation products. To distinguish cluster candidates from stars, we introduce the Multiple Concentration Index (MCI) parameter, and measure inner and outer MCIs to probe morphology in more detail than with a single, standard concentration index (CI). We improve upon cluster candidate selection, jointly basing our criteria on expectations for MCI derived from synthetic cluster populations and existing cluster catalogues, yielding model and semi-empirical selection regions (respectively). Selection purity (confirmed clusters versus candidates, assessed via human-based classification) is high (up to 70$\%$) for moderately luminous sources in the semi-empirical selection region, and somewhat lower overall (outside the region or fainter). The number of candidates rises steeply with decreasing luminosity, but pipeline-integrated Machine Learning (ML) classification prevents this from being problematic. We quantify the performance of our PHANGS-HST methods in comparison to LEGUS for a sample of four galaxies in common to both surveys, finding overall agreement with $50-75\%$ of human verified star clusters appearing in both catalogues, but also subtle differences attributable to specific choices adopted by each project. The PHANGS-HST ML-classified Class 1 or 2 catalogues reach $\sim1$ magnitude fainter, $\sim2\times$ lower stellar mass, and are $2{-}5\times$ larger in number, than attained in the human classified samples.


title: SeaNet -- Towards A Knowledge Graph Based Autonomic Management of Software Defined Networks

abstract: Automatic network management driven by Artificial Intelligent technologies has been heatedly discussed over decades. However, current reports mainly focus on theoretic proposals and architecture designs, works on practical implementations on real-life networks are yet to appear. This paper proposes our effort toward the implementation of knowledge graph driven approach for autonomic network management in software defined networks (SDNs), termed as SeaNet. Driven by the ToCo ontology, SeaNet is reprogrammed based on Mininet (a SDN emulator). It consists three core components, a knowledge graph generator, a SPARQL engine, and a network management API. The knowledge graph generator represents the knowledge in the telecommunication network management tasks into formally represented ontology driven model. Expert experience and network management rules can be formalized into knowledge graph and by automatically inferenced by SPARQL engine, Network management API is able to packet technology-specific details and expose technology-independent interfaces to users. The Experiments are carried out to evaluate proposed work by comparing with a commercial SDN controller Ryu implemented by the same language Python. The evaluation results show that SeaNet is considerably faster in most circumstances than Ryu and the SeaNet code is significantly more compact. Benefit from RDF reasoning, SeaNet is able to achieve O(1) time complexity on different scales of the knowledge graph while the traditional database can achieve O(nlogn) at its best. With the developed network management API, SeaNet enables researchers to develop semantic-intelligent applications on their own SDNs.


title: On Kaczmarz method with oblique projection for solving large overdetermined linear systems

abstract: In this paper, an extension of Kaczmarz method, the Kaczmarz method with oblique projection (KO), is introduced and analyzed. Using this method, a number of iteration steps to solve the over-determined systems of linear equations are significantly reduced, and the the computing time is much saved, especially for those problems that contain some linear equations with near-linear correlation. Simultaneously, a randomized version--randomized Kaczmarz method with oblique projection (RKO) is established. The convergence proofs of these two methods are given and numerical experiments show the effectiveness of the two methods for uniformly distributed random data. Especially when the system has correlated rows, the improvement of experimental results is very prominent.


title: Distributed Nash Equilibrium Seeking Algorithm Design for Multi-Cluster Games with High-Order Players

abstract: In this paper, a multi-cluster game with high-order players is investigated. Different from the well-known multi-cluster games, the dynamics of players are taken into account in our problem. Due to the high-order dynamics of players, existing algorithms for multi-cluster games cannot solve the problem. For purpose of seeking the Nash equilibrium of the game, we design a distributed algorithm based on gradient descent and state feedback, where a distributed estimator is embedded for the players to estimate the decisions of other players. Furthermore, we analyze the exponential convergence of the algorithm via variational analysis and Lyapunov stability theory. Finally, a numerical simulation verifies the effectiveness of our method.


title: $K_1$-injectivity of the Paschke dual algebra, and uniqueness

abstract: We prove that a large class of Paschke dual algebras of simple unital C*-algebras are $K_1$-injective. As a consequence, we obtain interesting $KK$-uniqueness theorems which generalize the Brown--Douglas--Fillmore essential codimension property.


title: Blistering Failure of Elastic Coatings with Applications to Corrosion Resistance

abstract: A variety of polymeric surfaces, such as anti-corrosion coatings and polymer-modified asphalts, are prone to blistering when exposed to moisture and air. As water and oxygen diffuse through the material, dissolved species are produced, which generate osmotic pressure that deforms and debonds the coating.These mechanisms are experimentally well-supported; however, comprehensive macroscopic models capable of predicting the formation osmotic blisters, without extensive data-fitting, is scant. Here, we develop a general mathematical theory of blistering and apply it to the failure of anti-corrosion coatings on carbon steel. The model is able to predict the irreversible, nonlinear blister growth dynamics, which eventually reaches a stable state, ruptures, or undergoes runaway delamination, depending on the mechanical and adhesion properties of the coating. For runaway delamination, the theory predicts a critical delamination length, beyond which unstable corrosion-driven growth occurs. The model is able to fit multiple sets of blister growth data with no fitting parameters. Corrosion experiments are also performed to observe undercoat rusting on carbon steel, which yielded trends comparable with model predictions. The theory is used to define three dimensionless numbers which can be used for engineering design of elastic coatings capable of resisting visible deformation, rupture, and delamination.


title: Graphs with Many Hamiltonian Paths

abstract: A graph is \emph{hamiltonian-connected} if every pair of vertices can be connected by a hamiltonian path, and it is \emph{hamiltonian} if it contains a hamiltonian cycle. We construct families of non-hamiltonian graphs for which the ratio of pairs of vertices connected by hamiltonian paths to all pairs of vertices approaches 1. We then consider minimal graphs that are hamiltonian-connected. It is known that any order-$n$ graph that is hamiltonian-connected must have $\geq 3n/2$ edges. We construct an infinite family of graphs realizing this minimum.


title: Temperature Control of PDE Constrained Optimization Problems Governed by Kobayashi--Warren--Carter Type Models of Grain Boundary Motions

abstract: In this paper, we consider a class of optimal control problems governed by state-equations of Kobayashi--Warren--Carter type. The control is given by physical temperature. The focus is on problems in dimensions less than equal to 4. The results are divided in four Main Theorems, concerned with: solvability and parameter-dependence of state-equations and optimal control problems; the first order necessary optimality conditions for these regularized optimal control problems. Subsequently, we derive the limiting systems and optimality conditions and study their well-posedness.


title: Nonequilibrium Stationary Process and Fluctuation-Dissipation Relations

abstract: A stochastic dynamics has a natural decomposition into a drift capturing mean rate of change and a martingale increment capturing randomness. They are two statistically uncorrelated, but not necessarily independent mechanisms contributing to the overall fluctuations of the dynamics, representing the uncertainties in the past and in the future. A generalized Einstein relation is a consequence solely because the dynamics being stationary; and the Green-Kubo formula reflects a balance between the two mechanisms. Equilibrium with reversibility is characterized by a novel covariance symmetry.


title: Domain-Specific Pretraining for Vertical Search: Case Study on Biomedical Literature

abstract: Information overload is a prevalent challenge in many high-value domains. A prominent case in point is the explosion of the biomedical literature on COVID-19, which swelled to hundreds of thousands of papers in a matter of months. In general, biomedical literature expands by two papers every minute, totalling over a million new papers every year. Search in the biomedical realm, and many other vertical domains is challenging due to the scarcity of direct supervision from click logs. Self-supervised learning has emerged as a promising direction to overcome the annotation bottleneck. We propose a general approach for vertical search based on domain-specific pretraining and present a case study for the biomedical domain. Despite being substantially simpler and not using any relevance labels for training or development, our method performs comparably or better than the best systems in the official TREC-COVID evaluation, a COVID-related biomedical search competition. Using distributed computing in modern cloud infrastructure, our system can scale to tens of millions of articles on PubMed and has been deployed as Microsoft Biomedical Search, a new search experience for biomedical literature: https://aka.ms/biomedsearch.


title: Exact solutions for scalars and spinors on quantized (Euclidean) $AdS_2$ and the correspondence principle

abstract: We obtain the exact solutions to the field equations for massless scalar and massless spinor fields on quantized two-dimensional anti-de Sitter space. We then apply the AdS/CFT correspondence principle to get exact answers for the two point correlation functions of the associated operators on the boundary. The results support the conclusion that conformal symmetry on the boundary is not spoiled by quantization of the bulk. Moreover, quantization of the bulk has no effect on the spinor two-point correlation function, while it induces an overall re-scaling in the scalar two-point correlation function.


title: Minimal quadrangulations of surfaces

abstract: A quadrangular embedding of a graph in a surface $\Sigma$, also known as a quadrangulation of $\Sigma$, is a cellular embedding in which every face is bounded by a $4$-cycle. A quadrangulation of $\Sigma$ is minimal if there is no quadrangular embedding of a (simple) graph of smaller order in $\Sigma$. In this paper we determine $n(\Sigma)$, the order of a minimal quadrangulation of a surface $\Sigma$, for all surfaces, both orientable and nonorientable. Letting $S_0$ denote the sphere and $N_2$ the Klein bottle, we prove that $n(S_0)=4, n(N_2)=6$, and $n(\Sigma)=\lceil (5+\sqrt{25-16\chi(\Sigma)})/2\rceil$ for all other surfaces $\Sigma$, where $\chi(\Sigma)$ is the Euler characteristic. Our proofs use a `diagonal technique', introduced by Hartsfield in 1994. We explain the general features of this method.


title: Schubert polynomials, the inhomogeneous TASEP, and evil-avoiding permutations

abstract: Consider a lattice of n sites arranged around a ring, with the $n$ sites occupied by particles of weights $\{1,2,\dots,n\}$; the possible arrangements of particles in sites thus corresponds to the $n!$ permutations in $S_n$. The inhomogeneous totally asymmetric simple exclusion process (or TASEP) is a Markov chain on the set of permutations, in which two adjacent particles of weights $i<j$ swap places at rate $x_i - y_{n+1-j}$ if the particle of weight $j$ is to the right of the particle of weight $i$. (Otherwise nothing happens.) In the case that $y_i=0$ for all $i$, the stationary distribution was conjecturally linked to Schubert polynomials by Lam-Williams, and explicit formulas for steady state probabilities were subsequently given in terms of multiline queues by Ayyer-Linusson and Arita-Mallick. In the case of general $y_i$, Cantini showed that $n$ of the $n!$ states have probabilities proportional to products of double Schubert polynomials. In this paper we introduce the class of evil-avoiding permutations, which are the permutations avoiding the patterns $2413, 4132, 4213$ and $3214$. We show that there are $\frac{(2+\sqrt{2})^{n-1}+(2-\sqrt{2})^{n-1}}{2}$ evil-avoiding permutations in $S_n$, and for each evil-avoiding permutation $w$, we give an explicit formula for the steady state probability $\psi_w$ as a product of double Schubert polynomials. We also show that the Schubert polynomials that arise in these formulas are flagged Schur functions, and give a bijection in this case between multiline queues and semistandard Young tableaux.


title: Bayesian Inference in High-Dimensional Time-Serieswith the Orthogonal Stochastic Linear Mixing Model

abstract: Many modern time-series datasets contain large numbers of output response variables sampled for prolonged periods of time. For example, in neuroscience, the activities of 100s-1000's of neurons are recorded during behaviors and in response to sensory stimuli. Multi-output Gaussian process models leverage the nonparametric nature of Gaussian processes to capture structure across multiple outputs. However, this class of models typically assumes that the correlations between the output response variables are invariant in the input space. Stochastic linear mixing models (SLMM) assume the mixture coefficients depend on input, making them more flexible and effective to capture complex output dependence. However, currently, the inference for SLMMs is intractable for large datasets, making them inapplicable to several modern time-series problems. In this paper, we propose a new regression framework, the orthogonal stochastic linear mixing model (OSLMM) that introduces an orthogonal constraint amongst the mixing coefficients. This constraint reduces the computational burden of inference while retaining the capability to handle complex output dependence. We provide Markov chain Monte Carlo inference procedures for both SLMM and OSLMM and demonstrate superior model scalability and reduced prediction error of OSLMM compared with state-of-the-art methods on several real-world applications. In neurophysiology recordings, we use the inferred latent functions for compact visualization of population responses to auditory stimuli, and demonstrate superior results compared to a competing method (GPFA). Together, these results demonstrate that OSLMM will be useful for the analysis of diverse, large-scale time-series datasets.


title: Asymptotic behavior of solutions to the Yamabe equation with an asymptotically flat metric

abstract: We prove that any positive solution of the Yamabe equation on an asymptotically flat $n$-dimensional manifold of flatness order at least $\frac{n-2}{2}$ and $n\le 24$ must converge at infinity either to a fundamental solution of the Laplace operator on the Euclidean space or to a radial Fowler solution defined on the entire Euclidean space. The flatness order $\frac{n-2}{2}$ is the minimal flatness order required to define ADM mass in general relativity; the dimension $24$ is the dividing dimension of the validity of compactness of solutions to the Yamabe problem. We also prove such alternatives for bounded solutions when $n>24$. We prove these results by establishing appropriate asymptotic behavior near an isolated singularity of solutions to the Yamabe equation when the metric has a flatness order of at least $\frac{n-2}{2}$ at the singularity and $n\le 24$, also when $n>24$ and the solution grows no faster than the fundamental solution of the flat metric Laplacian at the singularity. These results extend earlier results of L. Caffarelli, B. Gidas and J. Spruck, also of N. Korevaar, R. Mazzeo, F. Pacard and R. Schoen, when the metric is conformally flat, and work of C.C. Chen and C. S. Lin when the scalar curvature is a non-constant function with appropriate flatness at the singular point, also work of F. Marques when the metric is not necessarily conformally flat but smooth, and the dimension of the manifold is three, four, or five, as well as recent similar results by the second and third authors in dimension six.


title: To the Point: Efficient 3D Object Detection in the Range Image with Graph Convolution Kernels

abstract: 3D object detection is vital for many robotics applications. For tasks where a 2D perspective range image exists, we propose to learn a 3D representation directly from this range image view. To this end, we designed a 2D convolutional network architecture that carries the 3D spherical coordinates of each pixel throughout the network. Its layers can consume any arbitrary convolution kernel in place of the default inner product kernel and exploit the underlying local geometry around each pixel. We outline four such kernels: a dense kernel according to the bag-of-words paradigm, and three graph kernels inspired by recent graph neural network advances: the Transformer, the PointNet, and the Edge Convolution. We also explore cross-modality fusion with the camera image, facilitated by operating in the perspective range image view. Our method performs competitively on the Waymo Open Dataset and improves the state-of-the-art AP for pedestrian detection from 69.7% to 75.5%. It is also efficient in that our smallest model, which still outperforms the popular PointPillars in quality, requires 180 times fewer FLOPS and model parameters


title: A Source-Criticism Debiasing Method for GloVe Embeddings

abstract: It is well-documented that word embeddings trained on large public corpora consistently exhibit known human social biases. Although many methods for debiasing exist, almost all fixate on completely eliminating biased information from the embeddings and often diminish training set size in the process. In this paper, we present a simple yet effective method for debiasing GloVe word embeddings (Pennington et al., 2014) which works by incorporating explicit information about training set bias rather than removing biased data outright. Our method runs quickly and efficiently with the help of a fast bias gradient approximation method from Brunet et al. (2019). As our approach is akin to the notion of 'source criticism' in the humanities, we term our method Source-Critical GloVe (SC-GloVe). We show that SC-GloVe reduces the effect size on Word Embedding Association Test (WEAT) sets without sacrificing training data or TOP-1 performance.


title: Highly r-process enhanced stars in ultra-faint dwarf galaxies

abstract: Highly r-process enhanced metal-poor stars (MP r-II, $\rm [Eu/Fe]>1$ and $\rm [Fe/H]\lesssim-1.5$) have been observed in ultra-faint dwarf (UFD) galaxy, specifically in Reticulum~II (Ret~II). The fact that only a few UFDs contain such stars implies that the r-process site may reflect very rare, but individually prolific events, such as neutron star mergers (NSMs). Considering the relatively short star formation history (SFH) of UFDs, it is puzzling how they could experience such a rare phenomenon. In this work, we show the results of cosmological hydrodynamic zoom-in simulations of isolated UFDs ($M_{vir}\approx10^7-10^8$ solar mass and $M_{\ast}\approx10^3-10^4$ solar mass at $z=0$) to explain the formation of MP r-II stars in UFDs. We employ a simple toy model for NSM events, adopting parameters consistent with observations, such as the NSM rate (1 per $M_{\ast}\approx10^5$ solar mass) and europium (Eu) mass ($M_{Eu}\approx10^{-5}$ solar mass). We identify only one simulated galaxy ($ M_{vir}\approx4.6\times10^7$, $M_{\ast}\approx 3.4\times 10^3$ solar mass at $z=0$) with abundances similar to Ret~II in a simulation volume that hosts $\sim30$ UFD analogs, indicating that such abundances are possible but rare. By exploring a range of key parameters, we demonstrate that the most important factor in determining the formation of MP r-II stars in UFDs is how quickly subsequent stars can be formed out of r-process enriched gas. We find that it takes between 10 to 100~Myr to form the first and second burst of MP r-II stars. Over this period, Eu-polluted gas maintains the required high abundance ratios of $\rm [Eu/Fe]>1$.


title: Finite elements for div and divdiv conforming symmetric tensors in arbitrary dimension

abstract: Several div-conforming and divdiv-conforming finite elements for symmetric tensors on simplexes in arbitrary dimension are constructed in this work. The shape function space is first split as the trace space and the bubble space. The later is further decomposed into the null space of the differential operator and its orthogonal complement. Instead of characterization of these subspaces of the shape function space, characterization of the dual spaces are provided. Vector div-conforming finite elements are firstly constructed as an introductory example. Then new symmetric div-conforming finite elements are constructed. The dual subspaces are then used as build blocks to construct divdiv conforming finite elements.


title: Trends, Politics, Sentiments, and Misinformation: Understanding People's Reactions to COVID-19 During its Early Stages

abstract: The sudden outbreak of COVID-19 resulted in large volumes of data shared on different social media platforms. Analyzing and visualizing these data is doubtlessly essential to having a deep understanding of the pandemic's impacts on people's lives and their reactions to them. In this work, we conduct a large-scale spatiotemporal data analytic study to understand peoples' reactions to the COVID-19 pandemic during its early stages. In particular, we analyze a JSON-based dataset that is collected from news/messages/boards/blogs in English about COVID-19 over a period of 4 months, for a total of 5.2M posts. The data are collected from December 2019 to March 2020 from several social media platforms such as Facebook, LinkedIn, Pinterest, StumbleUpon and VK. Our study aims mainly to understand which implications of COVID-19 have interested social media users the most and how did they vary over time, the spatiotemporal distribution of misinformation, and the public opinion toward public figures during the pandemic. Our results can be used by many parties (e.g., governments, psychologists, etc.) to make more informative decisions, taking into account the actual interests and opinions of the people.


title: Balancing Accuracy and Fairness for Interactive Recommendation with Reinforcement Learning

abstract: Fairness in recommendation has attracted increasing attention due to bias and discrimination possibly caused by traditional recommenders. In Interactive Recommender Systems (IRS), user preferences and the system's fairness status are constantly changing over time. Existing fairness-aware recommenders mainly consider fairness in static settings. Directly applying existing methods to IRS will result in poor recommendation. To resolve this problem, we propose a reinforcement learning based framework, FairRec, to dynamically maintain a long-term balance between accuracy and fairness in IRS. User preferences and the system's fairness status are jointly compressed into the state representation to generate recommendations. FairRec aims at maximizing our designed cumulative reward that combines accuracy and fairness. Extensive experiments validate that FairRec can improve fairness, while preserving good recommendation quality.


title: Bayesian Eye Tracking

abstract: Model-based eye tracking has been a dominant approach for eye gaze tracking because of its ability to generalize to different subjects, without the need of any training data and eye gaze annotations. Model-based eye tracking, however, is susceptible to eye feature detection errors, in particular for eye tracking in the wild. To address this issue, we propose a Bayesian framework for model-based eye tracking. The proposed system consists of a cascade-Bayesian Convolutional Neural Network (c-BCNN) to capture the probabilistic relationships between eye appearance and its landmarks, and a geometric eye model to estimate eye gaze from the eye landmarks. Given a testing eye image, the Bayesian framework can generate, through Bayesian inference, the eye gaze distribution without explicit landmark detection and model training, based on which it not only estimates the most likely eye gaze but also its uncertainty. Furthermore, with Bayesian inference instead of point-based inference, our model can not only generalize better to different sub-jects, head poses, and environments but also is robust to image noise and landmark detection errors. Finally, with the estimated gaze uncertainty, we can construct a cascade architecture that allows us to progressively improve gaze estimation accuracy. Compared to state-of-the-art model-based and learning-based methods, the proposed Bayesian framework demonstrates significant improvement in generalization capability across several benchmark datasets and in accuracy and robustness under challenging real-world conditions.


title: Influences on Drivers' Understandings of Systems by Presenting Image Recognition Results

abstract: It is essential to help drivers have appropriate understandings of level 2 automated driving systems for keeping driving safety. A human machine interface (HMI) was proposed to present real time results of image recognition by the automated driving systems to drivers. It was expected that drivers could better understand the capabilities of the systems by observing the proposed HMI. Driving simulator experiments with 18 participants were preformed to evaluate the effectiveness of the proposed system. Experimental results indicated that the proposed HMI could effectively inform drivers of potential risks continuously and help drivers better understand the level 2 automated driving systems.


title: Energy-Based Generative Cooperative Saliency Prediction

abstract: Conventional saliency prediction models typically learn a deterministic mapping from an image to its saliency map, and thus fail to explain the subjective nature of human attention. In this paper, to model the uncertainty of visual saliency, we study the saliency prediction problem from the perspective of generative models by learning a conditional probability distribution over the saliency map given an input image, and treating the saliency prediction as a sampling process from the learned distribution. Specifically, we propose a generative cooperative saliency prediction framework, where a conditional latent variable model (LVM) and a conditional energy-based model (EBM) are jointly trained to predict salient objects in a cooperative manner. The LVM serves as a fast but coarse predictor to efficiently produce an initial saliency map, which is then refined by the iterative Langevin revision of the EBM that serves as a slow but fine predictor. Such a coarse-to-fine cooperative saliency prediction strategy offers the best of both worlds. Moreover, we propose a "cooperative learning while recovering" strategy and apply it to weakly supervised saliency prediction, where saliency annotations of training images are partially observed. Lastly, we find that the learned energy function in the EBM can serve as a refinement module that can refine the results of other pre-trained saliency prediction models. Experimental results show that our model can produce a set of diverse and plausible saliency maps of an image, and obtain state-of-the-art performance in both fully supervised and weakly supervised saliency prediction tasks.


title: Implementation of an alternative method for assessing competing risks: restricted mean time lost

abstract: In clinical and epidemiological studies, hazard ratios are often applied to compare treatment effects between two groups for survival data. For competing risks data, the corresponding quantities of interest are cause-specific hazard ratios (cHRs) and subdistribution hazard ratios (sHRs). However, they both have some limitations related to model assumptions and clinical interpretation. Therefore, we recommend restricted mean time lost (RMTL) as an alternative that is easy to interpret in a competing risks framework. Based on the difference in restricted mean time lost (RMTLd), we propose a new estimator, hypothetical test and sample size formula. The simulation results show that the estimation of the RMTLd is accurate and that the RMTLd test has robust statistical performance (both type I error and power). The results of three example analyses also verify the performance of the RMTLd test. From the perspectives of clinical interpretation, application conditions and statistical performance, we recommend that the RMTLd be reported with the HR in the analysis of competing risks data and that the RMTLd even be regarded as the primary outcome when the proportional hazard assumption fails. The R code (crRMTL) is publicly available from Github (https://github.com/chenzgz/crRMTL.1). Keywords: survival analysis, competing risks, hazard ratio, restricted mean time lost, sample size, hypothesis test


title: HAN: An Efficient Hierarchical Self-Attention Network for Skeleton-Based Gesture Recognition

abstract: Previous methods for skeleton-based gesture recognition mostly arrange the skeleton sequence into a pseudo picture or spatial-temporal graph and apply deep Convolutional Neural Network (CNN) or Graph Convolutional Network (GCN) for feature extraction. Although achieving superior results, these methods have inherent limitations in dynamically capturing local features of interactive hand parts, and the computing efficiency still remains a serious issue. In this work, the self-attention mechanism is introduced to alleviate this problem. Considering the hierarchical structure of hand joints, we propose an efficient hierarchical self-attention network (HAN) for skeleton-based gesture recognition, which is based on pure self-attention without any CNN, RNN or GCN operators. Specifically, the joint self-attention module is used to capture spatial features of fingers, the finger self-attention module is designed to aggregate features of the whole hand. In terms of temporal features, the temporal self-attention module is utilized to capture the temporal dynamics of the fingers and the entire hand. Finally, these features are fused by the fusion self-attention module for gesture classification. Experiments show that our method achieves competitive results on three gesture recognition datasets with much lower computational complexity.


title: Topological cavity based on slow light topological edge mode for broadband Purcell enhancement

abstract: Slow light in topological valley photonic crystal structures offers new possibilities to enhance light-matter interaction. We report a topological cavity based on slow light topological edge mode for broadband Purcell enhancement. The topological edge modes with large group indices over 100 can be realized with a bearded interface between two topologically distinct valley photonic crystals, featuring the greatly enhanced Purcell factor because of the increased local density of states. In the slow light regime, the topological cavity supports much more cavity modes with higher quality factor than that in the fast light regime, which is both demonstrated theoretically and experimentally. We demonstrate the cavity enables the broadband Purcell enhancement together with substantial Purcell factor, benefiting from dense cavity modes with high quality factor in a wide spectral range. It has great benefit to the realization of high-efficiency quantum-dot-based single-photon sources and entangled-photon sources with less restriction on spectral match. Such topological cavity could serve as a significant building block toward the development of photonic integrated circuits with embedded quantum emitters.


title: Interpreting Depression From Question-wise Long-term Video Recording of SDS Evaluation

abstract: Self-Rating Depression Scale (SDS) questionnaire has frequently been used for efficient depression preliminary screening. However, the uncontrollable self-administered measure can be easily affected by insouciantly or deceptively answering, and producing the different results with the clinician-administered Hamilton Depression Rating Scale (HDRS) and the final diagnosis. Clinically, facial expression (FE) and actions play a vital role in clinician-administered evaluation, while FE and action are underexplored for self-administered evaluations. In this work, we collect a novel dataset of 200 subjects to evidence the validity of self-rating questionnaires with their corresponding question-wise video recording. To automatically interpret depression from the SDS evaluation and the paired video, we propose an end-to-end hierarchical framework for the long-term variable-length video, which is also conditioned on the questionnaire results and the answering time. Specifically, we resort to a hierarchical model which utilizes a 3D CNN for local temporal pattern exploration and a redundancy-aware self-attention (RAS) scheme for question-wise global feature aggregation. Targeting for the redundant long-term FE video processing, our RAS is able to effectively exploit the correlations of each video clip within a question set to emphasize the discriminative information and eliminate the redundancy based on feature pair-wise affinity. Then, the question-wise video feature is concatenated with the questionnaire scores for final depression detection. Our thorough evaluations also show the validity of fusing SDS evaluation and its video recording, and the superiority of our framework to the conventional state-of-the-art temporal modeling methods.


title: Countering Adversarial Examples: Combining Input Transformation and Noisy Training

abstract: Recent studies have shown that neural network (NN) based image classifiers are highly vulnerable to adversarial examples, which poses a threat to security-sensitive image recognition task. Prior work has shown that JPEG compression can combat the drop in classification accuracy on adversarial examples to some extent. But, as the compression ratio increases, traditional JPEG compression is insufficient to defend those attacks but can cause an abrupt accuracy decline to the benign images. In this paper, with the aim of fully filtering the adversarial perturbations, we firstly make modifications to traditional JPEG compression algorithm which becomes more favorable for NN. Specifically, based on an analysis of the frequency coefficient, we design a NN-favored quantization table for compression. Considering compression as a data augmentation strategy, we then combine our model-agnostic preprocess with noisy training. We fine-tune the pre-trained model by training with images encoded at different compression levels, thus generating multiple classifiers. Finally, since lower (higher) compression ratio can remove both perturbations and original features slightly (aggressively), we use these trained multiple models for model ensemble. The majority vote of the ensemble of models is adopted as final predictions. Experiments results show our method can improve defense efficiency while maintaining original accuracy.


title: Volume and Monge-Amp\`ere energy on polarized affine varieties

abstract: Let $(X, \xi)$ be a polarized affine variety, i.e. an affine variety $X$ with a (possibly irrational) Reeb vector field $\xi$. We define the volume of a filtration of the coordinate ring of $X$ in terms of the asymptotics of the average of jumping numbers. When the filtration is finitely generated, it induces a Fubini-Study function $\varphi$ on the Berkovich analytification of $X$. In this case, we define the Monge-Amp\`ere energy for $\varphi$ using the theory of forms and currents on Berkovich spaces developed by Chambert-Loir and Ducros, and show that it agrees with the volume of the filtration. In the special case when the filtration comes from a test configuration, we recover the functional defined by Collins-Sz\'ekelyhidi and Li-Xu.


title: Explicit formula of radiation fields of free waves with applications on channel of energy

abstract: In this work we give a few explicit formulas regarding the radiation fields of linear free waves. We then apply these formulas on the channel of energy theory. We characterize all the radial weakly non-radiative solutions in all dimensions and give a few new exterior energy estimates.


title: Pheno-Mapper: An Interactive Toolbox for the Visual Exploration of Phenomics Data

abstract: High-throughput technologies to collect field data have made observations possible at scale in several branches of life sciences. The data collected can range from the molecular level (genotypes) to physiological (phenotypic traits) and environmental observations (e.g., weather, soil conditions). These vast swathes of data, collectively referred to as phenomics data, represent a treasure trove of key scientific knowledge on the dynamics of the underlying biological system. However, extracting information and insights from these complex datasets remains a significant challenge owing to their multidimensionality and lack of prior knowledge about their complex structure. In this paper, we present Pheno-Mapper, an interactive toolbox for the exploratory analysis and visualization of large-scale phenomics data. Our approach uses the mapper framework to perform a topological analysis of the data, and subsequently render visual representations with built-in data analysis and machine learning capabilities. We demonstrate the utility of this new tool on real-world plant (e.g., maize) phenomics datasets. In comparison to existing approaches, the main advantage of Pheno-Mapper is that it provides rich, interactive capabilities in the exploratory analysis of phenomics data, and it integrates visual analytics with data analysis and machine learning in an easily extensible way. In particular, Pheno-Mapper allows the interactive selection of subpopulations guided by a topological summary of the data and applies data mining and machine learning to these selected subpopulations for in-depth exploration.


title: Code-Verification Techniques for the Method-of-Moments Implementation of the Electric-Field Integral Equation

abstract: The method-of-moments implementation of the electric-field integral equation yields many code-verification challenges due to the various sources of numerical error and their possible interactions. Matters are further complicated by singular integrals, which arise from the presence of a Green's function. In this paper, we provide approaches to separately assess the numerical errors arising from the use of basis functions to approximate the solution and the use of quadrature to approximate the integration. Through these approaches, we are able to verify the code and compare the error from different quadrature options.


title: On the Resolution of Reductive Monoids and Multiplicativity of $\gamma$-factors

abstract: In this article, we give a proof of multiplicativity for $\gamma$-factors, an equality of parabolically induced and inducing factors, in the context of the Braverman-Kazhdan/Ngo program, under the assumption of commutativity of the corresponding Fourier transforms and a certain generalized Harish-Chandra transform. We also discuss the resolution of singularities and their rationality for reductive monoids, which are among the basic objects in the program.


title: Emergent gapless topological Luttinger liquid

abstract: Gapless Luttinger liquid is conventionally viewed as topologically trivial, unless it hosts degenerate ground states and or entanglement spectrum, which necessitates partial bulk degree of freedom to be gapped out. Here we predict an emergent gapless topological Luttinger liquid which is beyond the conventional scenarios and is characterized by the nontrivial many-body bulk spin texture, and propose feasible scheme for experimental observation. We consider a one-dimensional spin-orbit coupled Fermi-Hubbard model with fractional filling, whose low-energy physics is effectively described by a spinless Luttinger liquid and is trivial in the conventional characterization. We show that, as being tuned by the filling factor and interaction strength, the many-body ground state may exhibit nontrivial winding in its bulk spin texture in the projected momentum space, manifesting an emergent topological phase. A topological transition occurs when the projected spin-state at a high symmetry momentum becomes fully mixed one, resulting from competing processes of particle scattering to the lower and higher subbands, for which the spin texture at such momentum point is ill-defined, but the Luttinger liquid keeps gapless through the transition. Surprisingly, at relatively small filling the Luttinger liquid remains topologically nontrivial even at infinitely strong interaction. The results can be generalized to finite temperature which facilitates the real experimental detection. This work shows a novel gapless topological Luttinger liquid whose characterization is beyond the low-energy effective theory, and can be verified based on current experiments.


title: Decomposed Mutual Information Estimation for Contrastive Representation Learning

abstract: Recent contrastive representation learning methods rely on estimating mutual information (MI) between multiple views of an underlying context. E.g., we can derive multiple views of a given image by applying data augmentation, or we can split a sequence into views comprising the past and future of some step in the sequence. Contrastive lower bounds on MI are easy to optimize, but have a strong underestimation bias when estimating large amounts of MI. We propose decomposing the full MI estimation problem into a sum of smaller estimation problems by splitting one of the views into progressively more informed subviews and by applying the chain rule on MI between the decomposed views. This expression contains a sum of unconditional and conditional MI terms, each measuring modest chunks of the total MI, which facilitates approximation via contrastive bounds. To maximize the sum, we formulate a contrastive lower bound on the conditional MI which can be approximated efficiently. We refer to our general approach as Decomposed Estimation of Mutual Information (DEMI). We show that DEMI can capture a larger amount of MI than standard non-decomposed contrastive bounds in a synthetic setting, and learns better representations in a vision domain and for dialogue generation.


title: Efficient algorithms for computing rank-revealing factorizations on a GPU

abstract: Standard rank-revealing factorizations such as the singular value decomposition and column pivoted QR factorization are challenging to implement efficiently on a GPU. A major difficulty in this regard is the inability of standard algorithms to cast most operations in terms of the Level-3 BLAS. This paper presents two alternative algorithms for computing a rank-revealing factorization of the form $A = U T V^*$, where $U$ and $V$ are orthogonal and $T$ is triangular. Both algorithms use randomized projection techniques to cast most of the flops in terms of matrix-matrix multiplication, which is exceptionally efficient on the GPU. Numerical experiments illustrate that these algorithms achieve an order of magnitude acceleration over finely tuned GPU implementations of the SVD while providing low-rank approximation errors close to that of the SVD.


title: ParaLaw Nets -- Cross-lingual Sentence-level Pretraining for Legal Text Processing

abstract: Ambiguity is a characteristic of natural language, which makes expression ideas flexible. However, in a domain that requires accurate statements, it becomes a barrier. Specifically, a single word can have many meanings and multiple words can have the same meaning. When translating a text into a foreign language, the translator needs to determine the exact meaning of each element in the original sentence to produce the correct translation sentence. From that observation, in this paper, we propose ParaLaw Nets, a pretrained model family using sentence-level cross-lingual information to reduce ambiguity and increase the performance in legal text processing. This approach achieved the best result in the Question Answering task of COLIEE-2021.


title: Viscoelastic response of impact process on dense suspensions

abstract: We numerically study impact processes on dense suspensions using the lattice Boltzmann method to elucidate the connection between the elastic rebound of an impactor and relations among the impact speed $u_0$, maximum force acting on the impactor $F_{\rm max}$, and elapsed time $t_{\rm max}$ to reach $F_{\rm max}$. We find that $t_{\rm max}$ emerges in the early stage of the impact, while the rebound process takes place in the late stage. We find a crossover of $F_{\rm max}$ from $u_0$ independent regime for low $u_0$ to a power law regime satisfying $F_{\rm max}\propto u_0^\alpha$ with $\alpha\approx 1.5$ for high $u_0$. Similarly, $t_{\rm max}$ satisfies $t_{\rm max}\propto u_0^{\beta}$ with $\beta\approx -0.5$ for high $u_0$. Both power-law relations for $F_{\rm max}$ and $t_{\rm max}$ versus $u_0$ for high $u_0$ are independent of the system size, but the rebound phenomenon strongly depends on the depth of the container for suspensions. Thus, we indicate that the rebound phenomenon is not directly related to the relations among $u_0$, $F_{\rm max}$ and $t_{\rm max}$. We propose a floating + force chain model, where the rebound process is caused by an elastic term that is proportional to the number of the connected force chains from the impactor to the bottom plate. On the other hand, there are no elastic contributions in the relations for $F_{\rm max}$ and $t_{\rm max}$ against $u_0$ because of the absence of percolated force chains in the early stage. This phenomenology predicts $F_{\rm max}\propto u_0^{3/2}$ and $t_{\rm max}\propto u_0^{-1/2}$ for high $u_0$ and also recovers the behavior of the impactor quantitatively even if there is the rebound process.


title: JNLP Team: Deep Learning Approaches for Legal Processing Tasks in COLIEE 2021

abstract: COLIEE is an annual competition in automatic computerized legal text processing. Automatic legal document processing is an ambitious goal, and the structure and semantics of the law are often far more complex than everyday language. In this article, we survey and report our methods and experimental results in using deep learning in legal document processing. The results show the difficulties as well as potentials in this family of approaches.


title: Bounded geometry with no bounded pants decomposition

abstract: We construct a quasiconformally homogeneous hyperbolic Riemann surface-other than the hyperbolic plane-that does not admit a bounded pants decomposition. Also, given a connected orientable topological surface of infinite type with compact boundary components, we construct a complete hyperbolic metric on the surface that has bounded geometry but does not admit a bounded pants decomposition.


title: $L^2$-Extension of Adjoint bundles and Koll\'ar's Conjecture

abstract: We give a new proof of Koll\'ar's conjecture on the pushforward of the dualizing sheaf twisted by a variation of Hodge structure. This conjecture was settled by M. Saito via mixed Hodge modules and has applications in the investigation of Albanese maps. Our technique is the $L^2$-method and we give a concrete construction and proofs of the conjecture. The $L^2$ point of view allows us to generalize Koll\'ar's conjecture to the context of non-abelian Hodge theory.


title: The Enigmatic Brown Dwarf WISEA J153429.75-104303.3 (aka "The Accident")

abstract: Continued follow-up of WISEA J153429.75-104303.3, announced in Meisner et al (2020), has proven it to have an unusual set of properties. New imaging data from Keck/MOSFIRE and HST/WFC3 show that this object is one of the few faint proper motion sources known with J-ch2 > 8 mag, indicating a very cold temperature consistent with the latest known Y dwarfs. Despite this, it has W1-W2 and ch1-ch2 colors ~1.6 mag bluer than a typical Y dwarf. A new trigonometric parallax measurement from a combination of WISE, Spitzer, and HST astrometry confirms a nearby distance of $16.3^{+1.4}_{-1.2}$ pc and a large transverse velocity of $207.4{\pm}15.9$ km/s. The absolute J, W2, and ch2 magnitudes are in line with the coldest known Y dwarfs, despite the highly discrepant W1-W2 and ch1-ch2 colors. We explore possible reasons for the unique traits of this object and conclude that it is most likely an old, metal-poor brown dwarf and possibly the first Y subdwarf. Given that the object has an HST F110W magnitude of 24.7 mag, broad-band spectroscopy and photometry from JWST are the best options for testing this hypothesis.


title: Generative Modeling for Multi-task Visual Learning

abstract: Generative modeling has recently shown great promise in computer vision, but it has mostly focused on synthesizing visually realistic images. In this paper, motivated by multi-task learning of shareable feature representations, we consider a novel problem of learning a shared generative model that is useful across various visual perception tasks. Correspondingly, we propose a general multi-task oriented generative modeling (MGM) framework, by coupling a discriminative multi-task network with a generative network. While it is challenging to synthesize both RGB images and pixel-level annotations in multi-task scenarios, our framework enables us to use synthesized images paired with only weak annotations (i.e., image-level scene labels) to facilitate multiple visual tasks. Experimental evaluation on challenging multi-task benchmarks, including NYUv2 and Taskonomy, demonstrates that our MGM framework improves the performance of all the tasks by large margins, consistently outperforming state-of-the-art multi-task approaches.


title: Invariant Data-Driven Subgrid Stress Modeling in the Strain-Rate Eigenframe for Large Eddy Simulation

abstract: We present a new approach for constructing data-driven subgrid stress models for large eddy simulation of turbulent flows. The key to our approach is representation of model input and output tensors in the filtered strain rate eigenframe. Provided inputs and outputs are selected and non-dimensionalized in a suitable manner, this yields a model form that is symmetric, Galilean invariant, rotationally invariant, reflectionally invariant, and unit invariant. We use this model form to train a simple and efficient neural network model using only one time step of filtered direct numerical simulation data from a forced homogeneous isotropic turbulence simulation. We demonstrate the accuracy of this model as well as the model's ability to generalize to previously unseen filter widths, Reynolds numbers, and flow physics using a priori and a posteriori tests.


title: Fine-grained Geolocation Prediction of Tweets with Human Machine Collaboration

abstract: Twitter is a useful resource to analyze peoples' opinions on various topics. Often these topics are correlated or associated with locations from where these Tweet posts are made. For example, restaurant owners may need to know where their target customers eat with respect to the sentiment of the posts made related to food, policy planners may need to analyze citizens' opinion on relevant issues such as crime, safety, congestion, etc. with respect to specific parts of the city, or county or state. As promising as this is, less than $1\%$ of the crawled Tweet posts come with geolocation tags. That makes accurate prediction of Tweet posts for the non geo-tagged tweets very critical to analyze data in various domains. In this research, we utilized millions of Twitter posts and end-users domain expertise to build a set of deep neural network models using natural language processing (NLP) techniques, that predicts the geolocation of non geo-tagged Tweet posts at various level of granularities such as neighborhood, zipcode, and longitude with latitudes. With multiple neural architecture experiments, and a collaborative human-machine workflow design, our ongoing work on geolocation detection shows promising results that empower end-users to correlate relationship between variables of choice with the location information.


title: Broad-Band High-Energy Emission of the Gamma-ray Binary System LS 5039: Spectral and Temporal Features using NuSTAR and Fermi Observations

abstract: We report on detailed analysis of the hard X-ray and GeV gamma-ray spectra of LS 5039, one of the brightest gamma-ray binary system in the Galaxy. The NuSTAR observation covering its entire orbit in 2016 allowed us for the first time to study the orbital variability of the spectrum above 10 keV. The hard X-ray spectrum is well described with a single power-law component up to 78 keV. The X-ray flux showed a slight deviation from those observed previously with Suzaku in 2007. The fast X-ray brightening observed with Suzaku, around the inferior conjunction, was not observed in this observation. We also analyzed 11 years of Fermi Large Area Telescope data of LS 5039. The GeV spectrum around the inferior conjunction was well described with two non-thermal components; a power law with a photon index of $\sim 3$ and a cut-off power law with a cutoff energy of $\sim 2$ GeV. The orbital flux variability also changed gradually around a few GeV. These results indicate that there are two emission components in the GeV band, and the dominant component above $\sim 1$ GeV does not depend on the orbital phase. By combining these results, we update the spectral energy distribution of LS 5039 with the highest available statistics. Theoretical models proposed so far cannot explain the obtained multi-wavelength spectrum, especially the emission from $\sim$ 1 MeV to $\sim$ 400 MeV, and we discuss a possibility that particle acceleration in LS 5039 is different from the shock acceleration.


title: Is the free locally convex space $L(X)$ nuclear?

abstract: Given a class $\mathcal P$ of Banach spaces, a locally convex space (LCS) $E$ is called {\em multi-$\mathcal P$} if $E$ can be isomorphically embedded into a product of spaces that belong to $\mathcal P$. We investigate the question whether the free locally convex space $L(X)$ is strongly nuclear, nuclear, Schwartz, multi-Hilbert or multi-reflexive. If $X$ is a Tychonoff space containing an infinite compact subset then, as it follows from the results of \cite{Aus}, $L(X)$ is not nuclear. We prove that for such $X$ the free LCS $L(X)$ has the stronger property of not being multi-Hilbert. We deduce that if $X$ is a $k$-space, then the following properties are equivalent: (1) $L(X)$ is strongly nuclear; (2) $L(X)$ is nuclear; (3) $L(X)$ is multi-Hilbert; (4) $X$ is countable and discrete. On the other hand, we show that $L(X)$ is strongly nuclear for every projectively countable $P$-space (in particular, for every Lindel\"of $P$-space) $X$. We observe that every Schwartz LCS is multi-reflexive. It is known that if $X$ is a $k_\omega$-space, then $L(X)$ is a Schwartz LCS \cite{Chasco}, hence $L(X)$ is multi-reflexive. We show that for any first-countable paracompact (in particular, metrizable) space $X$ the converse is true, so $L(X)$ is multi-reflexive if and only if $X$ is a $k_\omega$-space, equivalently, if $X$ is a locally compact and $\sigma$-compact space. Similarly, we show that for any first-countable paracompact space $X$ the free abelian topological group $A(X)$ is a Schwartz group if and only if $X$ is a locally compact space such that the set $X^{(1)}$ of all non-isolated points of $X$ is $\sigma$-compact.


title: The Price of Tolerance in Distribution Testing

abstract: We revisit the problem of tolerant distribution testing. That is, given samples from an unknown distribution $p$ over $\{1, \dots, n\}$, is it $\varepsilon_1$-close to or $\varepsilon_2$-far from a reference distribution $q$ (in total variation distance)? Despite significant interest over the past decade, this problem is well understood only in the extreme cases. In the noiseless setting (i.e., $\varepsilon_1 = 0$) the sample complexity is $\Theta(\sqrt{n})$, strongly sublinear in the domain size. At the other end of the spectrum, when $\varepsilon_1 = \varepsilon_2/2$, the sample complexity jumps to the barely sublinear $\Theta(n/\log n)$. However, very little is known about the intermediate regime. We fully characterize the price of tolerance in distribution testing as a function of $n$, $\varepsilon_1$, $\varepsilon_2$, up to a single $\log n$ factor. Specifically, we show the sample complexity to be \[\tilde \Theta\left(\frac{\sqrt{n}}{\varepsilon_2^{2}} + \frac{n}{\log n} \cdot \max \left\{\frac{\varepsilon_1}{\varepsilon_2^2},\left(\frac{\varepsilon_1}{\varepsilon_2^2}\right)^{\!\!2}\right\}\right),\] providing a smooth tradeoff between the two previously known cases. We also provide a similar characterization for the problem of tolerant equivalence testing, where both $p$ and $q$ are unknown. Surprisingly, in both cases, the main quantity dictating the sample complexity is the ratio $\varepsilon_1/\varepsilon_2^2$, and not the more intuitive $\varepsilon_1/\varepsilon_2$. Of particular technical interest is our lower bound framework, which involves novel approximation-theoretic tools required to handle the asymmetry between $\varepsilon_1$ and $\varepsilon_2$, a challenge absent from previous works.

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