title: On the existence and regularity of local times

abstract: We study the existence and regularity of local times for general $d$-dimensional stochastic processes. We give a general condition for their existence and regularity properties. To emphasize the contribution of our results, we show that they include various prominent examples, among others solutions to stochastic differential equations driven by fractional Brownian motion, where the behavior of the local time was not fully understood up to now and remained as an open problem in the stochastic analysis literature. In particular this completes the picture regarding the local time behavior of such equations, above all includes high dimensions and both large and small Hurst parameters. As other main examples, we also show that by using our general approach, one can quite easily cover and extend some recently obtained results on the local times of the Rosenblatt process and Gaussian quasi-helices.


title: The Small Separation A-Star Companion Population: First Results with CHARA/MIRC-X

abstract: We present preliminary results from our long-baseline interferometry (LBI) survey to constrain the multiplicity properties of intermediate-mass A-type stars within 80pc. Previous multiplicity studies of nearby stars exhibit orbital separation distributions well-fitted with a log-normal with peaks > 15au, increasing with primary mass. The A-star multiplicity survey of De Rosa et al. (2014), sensitive beyond 30au but incomplete below 100 au, found a log-normal peak around 390au. Radial velocity surveys of slowly-rotating, chemically peculiar Am stars identified a significant number of very close companions with periods $\leq$ 5 days, ~ 0.1au, a result similar to surveys of O- and B-type primaries. With the improved performance of LBI techniques, we can probe these close separations for normal A-type stars where other surveys are incomplete. Our initial sample consists of 27 A-type primaries with estimated masses between 1.44-2.49M$_{\odot}$ and ages 10-790Myr, which we observed with the MIRC-X instrument at the CHARA Array. We use the open source software CANDID to detect five companions, three of which are new, and derive a companion frequency of 0.19$^{+0.11}_{-0.06}$ over mass ratios 0.25-1.0 and projected separations 0.288-5.481 au. We find a probability of 10$^{-6}$ that our results are consistent with extrapolations based on previous models of the A-star companion population, over mass ratios and separations sampled. Our results show the need to explore these very close separations to inform our understanding of stellar formation and evolution processes.


title: Inverse Kinematics with Dual-Quaternions, Exponential-Maps, and Joint Limits

abstract: We present a novel approach for solving articulated inverse kinematic problems (e.g., character structures) by means of an iterative dual-quaternion and exponentialmapping approach. As dual-quaternions are a break from the norm and offer a straightforward and computationally efficient technique for representing kinematic transforms (i.e., position and translation). Dual-quaternions are capable of represent both translation and rotation in a unified state space variable with its own set of algebraic equations for concatenation and manipulation. Hence, an articulated structure can be represented by a set of dual-quaternion transforms, which we can manipulate using inverse kinematics (IK) to accomplish specific goals (e.g., moving end-effectors towards targets). We use the projected Gauss-Seidel iterative method to solve the IK problem with joint limits. Our approach is flexible and robust enough for use in interactive applications, such as games. We use numerical examples to demonstrate our approach, which performed successfully in all our test cases and produced pleasing visual results.


title: Generative Entity-to-Entity Stance Detection with Knowledge Graph Augmentation

abstract: Stance detection is typically framed as predicting the sentiment in a given text towards a target entity. However, this setup overlooks the importance of the source entity, i.e., who is expressing the opinion. In this paper, we emphasize the need for studying interactions among entities when inferring stances. We first introduce a new task, entity-to-entity (E2E) stance detection, which primes models to identify entities in their canonical names and discern stances jointly. To support this study, we curate a new dataset with 10,619 annotations labeled at the sentence-level from news articles of different ideological leanings. We present a novel generative framework to allow the generation of canonical names for entities as well as stances among them. We further enhance the model with a graph encoder to summarize entity activities and external knowledge surrounding the entities. Experiments show that our model outperforms strong comparisons by large margins. Further analyses demonstrate the usefulness of E2E stance detection for understanding media quotation and stance landscape, as well as inferring entity ideology.


title: A New Approach to Estimating Effective Resistances and Counting Spanning Trees in Expander Graphs

abstract: We demonstrate that for expander graphs, for all $\epsilon > 0,$ there exists a data structure of size $\widetilde{O}(n\epsilon^{-1})$ which can be used to return $(1 + \epsilon)$-approximations to effective resistances in $\widetilde{O}(1)$ time per query. Short of storing all effective resistances, previous best approaches could achieve $\widetilde{O}(n\epsilon^{-2})$ size and $\widetilde{O}(\epsilon^{-2})$ time per query by storing Johnson-Lindenstrauss vectors for each vertex, or $\widetilde{O}(n\epsilon^{-1})$ size and $\widetilde{O}(n\epsilon^{-1})$ time per query by storing a spectral sketch. Our construction is based on two key ideas: 1) $\epsilon^{-1}$-sparse, $\epsilon$-additive approximations to $DL^+1_u$ for all $u,$ can be used to recover $(1 + \epsilon)$-approximations to the effective resistances, 2) In expander graphs, only $\widetilde{O}(\epsilon^{-1})$ coordinates of a vector similar to $DL^+1_u$ are larger than $\epsilon.$ We give an efficient construction for such a data structure in $\widetilde{O}(m + n\epsilon^{-2})$ time via random walks. This results in an algorithm for computing $(1+\epsilon)$-approximate effective resistances for $s$ vertex pairs in expanders that runs in $\widetilde{O}(m + n\epsilon^{-2} + s)$ time, improving over the previously best known running time of $m^{1 + o(1)} + (n + s)n^{o(1)}\epsilon^{-1.5}$ for $s = \omega(n\epsilon^{-0.5}).$ We employ the above algorithm to compute a $(1+\delta)$-approximation to the number of spanning trees in an expander graph, or equivalently, approximating the (pseudo)determinant of its Laplacian in $\widetilde{O}(m + n^{1.5}\delta^{-1})$ time. This improves on the previously best known result of $m^{1+o(1)} + n^{1.875+o(1)}\delta^{-1.75}$ time, and matches the best known size of determinant sparsifiers.


title: Restricted Secant Varieties of Grassmannians

abstract: Restricted secant varieties of Grassmannians are constructed from sums of points corresponding to $k$-planes with the restriction that their intersection has a prescribed dimension. We study dimensions of restricted secant of Grassmannians and relate them to the analogous question for secants of Grassmannians via an incidence variety construction. We define a notion of expected dimension and give a formula for the dimension of all restricted secant varieties of Grassmannians that holds if the Baur-Draisma-deGraaf Conjecture on non-defectivity of Grassmannians is true. We also demonstrate example calculations in Macaulay2, and point out ways to make these calculations more efficient. We also show a potential application to coding theory.


title: Invisible Charm Exotica

abstract: One possible interpretation of two narrow $\Omega^0_c$ states reported by the LHCb Collaboration at CERN in 2017 is that they are pentaquarks belonging to a $\overline{\boldsymbol{15}}$ exotic SU(3) representation, as predicted by the Chiral Quark--Soliton Model. If so, there must exist a number of other exotic states since the model predicts three different $\overline{\boldsymbol{15}}$ multiplets of heavy baryons. We show, that depending on the soliton spin $J$, these states are either very narrow or very broad. This explains why they might have escaped experimental observation. Furthermore, we show that the lightest members of these multiplets are stable against two body strong decays.


title: Dual Generator Offline Reinforcement Learning

abstract: In offline RL, constraining the learned policy to remain close to the data is essential to prevent the policy from outputting out-of-distribution (OOD) actions with erroneously overestimated values. In principle, generative adversarial networks (GAN) can provide an elegant solution to do so, with the discriminator directly providing a probability that quantifies distributional shift. However, in practice, GAN-based offline RL methods have not performed as well as alternative approaches, perhaps because the generator is trained to both fool the discriminator and maximize return -- two objectives that can be at odds with each other. In this paper, we show that the issue of conflicting objectives can be resolved by training two generators: one that maximizes return, with the other capturing the ``remainder'' of the data distribution in the offline dataset, such that the mixture of the two is close to the behavior policy. We show that not only does having two generators enable an effective GAN-based offline RL method, but also approximates a support constraint, where the policy does not need to match the entire data distribution, but only the slice of the data that leads to high long term performance. We name our method DASCO, for Dual-Generator Adversarial Support Constrained Offline RL. On benchmark tasks that require learning from sub-optimal data, DASCO significantly outperforms prior methods that enforce distribution constraint.


title: The Need for Medically Aware Video Compression in Gastroenterology

abstract: Compression is essential to storing and transmitting medical videos, but the effect of compression on downstream medical tasks is often ignored. Furthermore, systems in practice rely on standard video codecs, which naively allocate bits between medically relevant frames or parts of frames. In this work, we present an empirical study of some deficiencies of classical codecs on gastroenterology videos, and motivate our ongoing work to train a learned compression model for colonoscopy videos. We show that two of the most common classical codecs, H264 and HEVC, compress medically relevant frames statistically significantly worse than medically nonrelevant ones, and that polyp detector performance degrades rapidly as compression increases. We explain how a learned compressor could allocate bits to important regions and allow detection performance to degrade more gracefully. Many of our proposed techniques generalize to medical video domains beyond gastroenterology


title: Fluent APIs in Functional Languages (full version)

abstract: Fluent API is an object-oriented pattern for elegant APIs and embedded DSLs. A smart fluent API can enforce the API protocol or DSL syntax at compile time. As fluent API designs typically rely on function overloading, they are hard to realize in functional programming languages. We show how to write functional fluent APIs using parametric polymorphism and type inference instead of overloading. Our designs support all regular and deterministic context-free API protocols and beyond.


title: Peering into the Young Planetary System AB Pic. Atmosphere, Orbit, Obliquity & Second Planetary Candidate

abstract: We aim to revisit the system AB Pic which has a known companion at the exoplanet/ brown-dwarf boundary. We based this study on a rich set of observations to investigate the companion's orbit and atmosphere. We composed a spectrum of AB Pic b merging archival VLT/SINFONI K-band data, with published spectra at J and H-band (SINFONI) and Lp-band (Magellan-AO), and photometric measurements (HST and Spitzer). We modeled the spectrum with ForMoSA, based on two atmospheric models: ExoREM and BT-SETTL13. We determined the orbital properties of b fitting the astrometric measurements from NaCo (2003 and 2004) and SPHERE (2015). The orbital solutions favor a semi-major axis of $\sim$190au viewed edge-on. With Exo-REM, we derive a T$_{eff}$ of 1700$\pm$50K and surface gravity of 4.5$\pm$0.3dex, consistent with previous works, and we report for the first time a C/O ratio of 0.58$\pm$0.08 ($\sim$solar). The posteriors are sensitive to the wavelength interval and the family of models used. Given the 2.1hr rotation period and our vsin(i) of $\sim$73km/s, we estimate for the first time the true obliquity to be $\sim$45 or $\sim$135deg, indicating a significant misalignment between the planet's spin and orbit orientations. Finally, a proper motion anomaly between the Hipparcos and Gaia eDR3 compared to our SPHERE detection limits and adapted radial velocity limits indicate the existence of a $\sim$6M$_{Jup}$ inner planet orbiting from 2 to 10au (40-200mas). The possible existence of an inner companion and the likely miss-alignment of the spin axis orientation strongly favor a formation path by gravitational instability or core accretion within a disk closer inside followed by dynamical interactions. Confirmation and characterization of planet c and access to a broader wavelength coverage for planet b will be essential to probe the uncertainties associated with the parameters.


title: Insensitizing controls for a fourth order semi-linear parabolic equations

abstract: This paper is concerned with the existence of insensitizing controls for a fourth order semilinear parabolic equation. Here, the initial data is partially unknown, we would like to find controls such that a specific functional is insensitive for small perturbations of the initial data. In general, this kind of problems can be recast as a null controllability problem for a nonlinear cascade system. We will first prove a null controllability result for a linear problem by global Carleman estimates and dual arguments. Then, by virtue of Leray-Schauder's fixed points theorem, we conclude the null controllability for the cascade system in the semi-linear case.


title: Integrated Photonic Tensor Processing Unit for a Matrix Multiply: a Review

abstract: The explosion of artificial intelligence and machine-learning algorithms, connected to the exponential growth of the exchanged data, is driving a search for novel application-specific hardware accelerators. Among the many, the photonics field appears to be in the perfect spotlight for this global data explosion, thanks to its almost infinite bandwidth capacity associated with limited energy consumption. In this review, we will overview the major advantages that photonics has over electronics for hardware accelerators, followed by a comparison between the major architectures implemented on Photonics Integrated Circuits (PIC) for both the linear and nonlinear parts of Neural Networks. By the end, we will highlight the main driving forces for the next generation of photonic accelerators, as well as the main limits that must be overcome.


title: On the practical usefulness of the Hardware Efficient Ansatz

abstract: Variational Quantum Algorithms (VQAs) and Quantum Machine Learning (QML) models train a parametrized quantum circuit to solve a given learning task. The success of these algorithms greatly hinges on appropriately choosing an ansatz for the quantum circuit. Perhaps one of the most famous ansatzes is the one-dimensional layered Hardware Efficient Ansatz (HEA), which seeks to minimize the effect of hardware noise by using native gates and connectives. The use of this HEA has generated a certain ambivalence arising from the fact that while it suffers from barren plateaus at long depths, it can also avoid them at shallow ones. In this work, we attempt to determine whether one should, or should not, use a HEA. We rigorously identify scenarios where shallow HEAs should likely be avoided (e.g., VQA or QML tasks with data satisfying a volume law of entanglement). More importantly, we identify a Goldilocks scenario where shallow HEAs could achieve a quantum speedup: QML tasks with data satisfying an area law of entanglement. We provide examples for such scenario (such as Gaussian diagonal ensemble random Hamiltonian discrimination), and we show that in these cases a shallow HEA is always trainable and that there exists an anti-concentration of loss function values. Our work highlights the crucial role that input states play in the trainability of a parametrized quantum circuit, a phenomenon that is verified in our numerics.


title: A machine learning model to identify corruption in M\'exico's public procurement contracts

abstract: The costs and impacts of government corruption range from impairing a country's economic growth to affecting its citizens' well-being and safety. Public contracting between government dependencies and private sector instances, referred to as public procurement, is a fertile land of opportunity for corrupt practices, generating substantial monetary losses worldwide. Thus, identifying and deterring corrupt activities between the government and the private sector is paramount. However, due to several factors, corruption in public procurement is challenging to identify and track, leading to corrupt practices going unnoticed. This paper proposes a machine learning model based on an ensemble of random forest classifiers, which we call hyper-forest, to identify and predict corrupt contracts in M\'exico's public procurement data. This method's results correctly detect most of the corrupt and non-corrupt contracts evaluated in the dataset. Furthermore, we found that the most critical predictors considered in the model are those related to the relationship between buyers and suppliers rather than those related to features of individual contracts. Also, the method proposed here is general enough to be trained with data from other countries. Overall, our work presents a tool that can help in the decision-making process to identify, predict and analyze corruption in public procurement contracts.


title: Topological Amenability of Semihypergroups

abstract: In this article, we introduce and explore the notion of topological amenability in the broad setting of (locally compact) semihypergroups. We acquire several stationary, ergodic and Banach algebraic characterizations of the same in terms of convergence of certain probability measures, total variation of convolution with probability measures and translation of certain functionals, as well as the F-algebraic properties of the associated measure algebra. We further investigate the interplay between restriction of convolution product and convolution of restrictions of measures on a sub-semihypergroup. Finally, we discuss and characterize topological amenability of sub-semihypergroups in terms of certain invariance properties attained on the corresponding measure algebra of the parent semihypergroup. This in turn provides us with an affirmative answer to an open question posed by J. Wong in 1980.


title: Over-communicate no more: Situated RL agents learn concise communication protocols

abstract: While it is known that communication facilitates cooperation in multi-agent settings, it is unclear how to design artificial agents that can learn to effectively and efficiently communicate with each other. Much research on communication emergence uses reinforcement learning (RL) and explores unsituated communication in one-step referential tasks -- the tasks are not temporally interactive and lack time pressures typically present in natural communication. In these settings, agents may successfully learn to communicate, but they do not learn to exchange information concisely -- they tend towards over-communication and an inefficient encoding. Here, we explore situated communication in a multi-step task, where the acting agent has to forgo an environmental action to communicate. Thus, we impose an opportunity cost on communication and mimic the real-world pressure of passing time. We compare communication emergence under this pressure against learning to communicate with a cost on articulation effort, implemented as a per-message penalty (fixed and progressively increasing). We find that while all tested pressures can disincentivise over-communication, situated communication does it most effectively and, unlike the cost on effort, does not negatively impact emergence. Implementing an opportunity cost on communication in a temporally extended environment is a step towards embodiment, and might be a pre-condition for incentivising efficient, human-like communication.


title: Physics-inspired machine learning for power grid frequency modelling

abstract: The operation of power systems is affected by diverse technical, economic and social factors. Social behaviour determines load patterns, electricity markets regulate the generation and weather-dependent renewables introduce power fluctuations. Thus, power system dynamics must be regarded as a non-autonomous system whose parameters vary strongly with time. However, the external driving factors are usually only available on coarse scales and the actual dependencies of the dynamic system parameters are generally unknown. Here, we propose a physics-inspired machine learning model that bridges the gap between large-scale drivers and short-term dynamics of the power system. Integrating stochastic differential equations and artificial neural networks, we construct a probabilistic model of the power grid frequency dynamics in Continental Europe. Its probabilistic prediction outperforms the daily average profile, which is an important benchmark. Using the integrated model, we identify and explain the parameters of the dynamical system from the data, which reveals their strong time-dependence and their relation to external drivers such as wind power feed-in and fast generation ramps. Finally, we generate synthetic time series from the model, which successfully reproduce central characteristics of the grid frequency such as their heavy-tailed distribution. All in all, our work emphasises the importance of modelling power system dynamics as a stochastic non-autonomous system with both intrinsic dynamics and external drivers.


title: RQUGE: Reference-Free Metric for Evaluating Question Generation by Answering the Question

abstract: Existing metrics for evaluating the quality of automatically generated questions such as BLEU, ROUGE, BERTScore, and BLEURT compare the reference and predicted questions, providing a high score when there is a considerable lexical overlap or semantic similarity between the candidate and the reference questions. This approach has two major shortcomings. First, we need expensive human-provided reference questions. Second, it penalises valid questions that may not have high lexical or semantic similarity to the reference questions. In this paper, we propose a new metric, RQUGE, based on the answerability of the candidate question given the context. The metric consists of a question-answering and a span scorer modules, using pre-trained models from existing literature, thus it can be used without any further training. We demonstrate that RQUGE has a higher correlation with human judgment without relying on the reference question. Additionally, RQUGE is shown to be more robust to several adversarial corruptions. Furthermore, we illustrate that we can significantly improve the performance of QA models on out-of-domain datasets by fine-tuning on synthetic data generated by a question generation model and re-ranked by RQUGE.


title: Unravelling the origin of the peculiar transition in the magnetically ordered phase of the Weyl semimetal Co3Sn2S2

abstract: Recent discovery of topologically non-trivial behavior in Co3Sn2S2 stimulated a notable interest in this itinerant ferromagnet (Tc = 174 K). The exact magnetic state remains ambiguous, with several reports indicating the existence of a second transition in the range 125 -- 130 K, with antiferromagnetic and glassy phases proposed to coexist with the ferromagnetic phase. Using detailed angle-dependent DC and AC magnetization measurements on large, high-quality single crystals we reveal a highly anisotropic behavior of both static and dynamic response of Co3Sn2S2. It is established that many observations related to sharp magnetization changes when B || c are influenced by the demagnetization factor of a sample. On the other hand, a genuine transition has been found at Tp = 128 K, with the magnetic response being strictly perpendicular to the c-axis and several orders of magnitude smaller than for B || c. Calculations using density-functional theory indicate that the ground state magnetic structure consist of magnetic moments canted away from the c-axis by a small angle (~ 1.5deg). We argue that the second transition originates from a small additional canting of moments within the kagome plane, with two equivalent orientations for each spin.


title: Data Level Lottery Ticket Hypothesis for Vision Transformers

abstract: The conventional lottery ticket hypothesis (LTH) claims that there exists a sparse subnetwork within a dense neural network and a proper random initialization method called the winning ticket, such that it can be trained from scratch to almost as good as the dense counterpart. Meanwhile, the research of LTH in vision transformers (ViTs) is scarcely evaluated. In this paper, we first show that the conventional winning ticket is hard to find at the weight level of ViTs by existing methods. Then, we generalize the LTH for ViTs to input data consisting of image patches inspired by the input dependence of ViTs. That is, there exists a subset of input image patches such that a ViT can be trained from scratch by using only this subset of patches and achieve similar accuracy to the ViTs trained by using all image patches. We call this subset of input patches the em winning tickets, which represent a significant amount of information in the input data. We use a ticket selector to generate the winning tickets based on the informativeness of patches for various types of ViT, including DeiT, LV-ViT, and Swin Transformers. The experiments show that there is a clear difference between the performance of models trained with winning tickets and randomly selected subsets, which verifies our proposed theory. We elaborate on the analogical similarity between our proposed Data-LTH-ViTs and the conventional LTH to further verify the integrity of our theory. The Source codes are available at https://github.com/shawnricecake/vit-lottery-ticket-input.


title: Information compression at the turbulent-phase transition in cold atom gases

abstract: The statistical properties of physical systems in thermal equilibrium are blatantly different from their far-from-equilibrium counterparts. In the latter, fluctuations often dominate the dynamics and might cluster in ordered patterns in the form of dissipative coherent structures. Here, we study the transition of a cold atomic cloud, driven close to a sharp electronic resonance, from a stable to a turbulent phase. From the atomic density distribution -- measured using a spatially-resolved pump-probe technique -- we have computed the Shannon entropy on two different basis sets. Information compression, corresponding to a minimum in the Shannon entropy, has been observed at criticality, where the system fluctuations organize into high-order (low-entropy) patterns. Being independent of the representation used, this feature is a property shared by a vast class of physical systems undergoing phase transitions.


title: Assessing Resource-Performance Trade-off of Natural Language Models using Data Envelopment Analysis

abstract: Natural language models are often summarized through a high-dimensional set of descriptive metrics including training corpus size, training time, the number of trainable parameters, inference times, and evaluation statistics that assess performance across tasks. The high dimensional nature of these metrics yields challenges with regard to objectively comparing models; in particular it is challenging to assess the trade-off models make between performance and resources (compute time, memory, etc.). We apply Data Envelopment Analysis (DEA) to this problem of assessing the resource-performance trade-off. DEA is a nonparametric method that measures productive efficiency of abstract units that consume one or more inputs and yield at least one output. We recast natural language models as units suitable for DEA, and we show that DEA can be used to create an effective framework for quantifying model performance and efficiency. A central feature of DEA is that it identifies a subset of models that live on an efficient frontier of performance. DEA is also scalable, having been applied to problems with thousands of units. We report empirical results of DEA applied to 14 different language models that have a variety of architectures, and we show that DEA can be used to identify a subset of models that effectively balance resource demands against performance.


title: Multi-vehicle Conflict Resolution in Highly Constrained Spaces by Merging Optimal Control and Reinforcement Learning

abstract: We present a novel method to address the problem of multi-vehicle conflict resolution in highly constrained spaces. An optimal control problem is formulated to incorporate nonlinear, non-holonomic vehicle dynamics and exact collision avoidance constraints. A solution to the problem can be obtained by first learning configuration strategies with reinforcement learning (RL) in a simplified discrete environment, and then using these strategies to shape the constraint space of the original problem. Simulation results show that our method can explore efficient actions to resolve conflicts in confined space and generate dexterous maneuvers that are both collision-free and kinematically feasible.


title: A study linking patient EHR data to external death data at Stanford Medicine

abstract: This manuscript explores linking real-world patient data with external death data in the context of research Clinical Data Warehouses (r-CDWs). We specifically present the linking of Electronic Health Records (EHR) data for Stanford Health Care (SHC) patients and data from the Social Security Administration (SSA) Limited Access Death Master File (LADMF) made available by the US Department of Commerce's National Technical Information Service (NTIS). The data analysis framework presented in this manuscript extends prior approaches and is generalizable to linking any two cross-organizational real-world patient data sources. Electronic Health Record (EHR) data and NTIS LADMF are heavily used resources at other medical centers and we expect that the methods and learnings presented here will be valuable to others. Our findings suggest that strong linkages are incomplete and weak linkages are noisy i.e., there is no good linkage rule that provides coverage and accuracy. Furthermore, the best linkage rule for any two datasets is different from the best linkage rule for two other datasets i.e., there is no generalization of linkage rules. Finally, LADMF, a commonly used external death data resource for r-CDWs, has a significant gap in death data making it necessary for r-CDWs to seek out more than one external death data source. We anticipate that presentation of multiple linkages will make it hard to present the linkage outcome to the end user. This manuscript is a resource in support of Stanford Medicine STARR (STAnford medicine Research data Repository) r-CDWs. The data are stored and analyzed as PHI in our HIPAA-compliant data center and are used under research and development (R&D) activities of STARR IRB.


title: Chip-scale Simulations in a Quantum-correlated Synthetic Space

abstract: An efficient simulator for quantum systems is one of the original goals for the efforts to develop a quantum computer [1]. In recent years, synthetic dimension in photonics [2] have emerged as a potentially powerful approach for simulation that is free from the constraint of geometric dimensionality. Here we demonstrate a quantum-correlated synthetic crystal, based upon a coherently-controlled broadband quantum frequency comb produced in a chip-scale dynamically modulated lithium niobate microresonator. The time-frequency entanglement inherent with the comb modes significantly extends the dimensionality of the synthetic space, creating a massive nearly 400 x 400 synthetic lattice with electrically-controlled tunability. With such a system, we are able to utilize the evolution of quantum correlations between entangled photons to perform a series of simulations, demonstrating quantum random walks, Bloch oscillations, and multi-level Rabi oscillations in the time and frequency correlation space. The device combines the simplicity of monolithic nanophotonic architecture, high dimensionality of a quantum-correlated synthetic space, and on-chip coherent control, which opens up an avenue towards chip-scale implementation of large-scale analog quantum simulation and computation [1,3,4] in the time-frequency domain.


title: Deep Learning Hamiltonians from Disordered Image Data in Quantum Materials

abstract: The capabilities of image probe experiments are rapidly expanding, providing new information about quantum materials on unprecedented length and time scales. Many such materials feature inhomogeneous electronic properties with intricate pattern formation on the observable surface. This rich spatial structure contains information about interactions, dimensionality, and disorder -- a spatial encoding of the Hamiltonian driving the pattern formation. Image recognition techniques from machine learning are an excellent tool for interpreting information encoded in the spatial relationships in such images. Here, we develop a deep learning framework for using the rich information available in these spatial correlations in order to discover the underlying Hamiltonian driving the patterns. We first vet the method on a known case, scanning near-field optical microscopy on a thin film of VO2. We then apply our trained convolutional neural network architecture to new optical microscope images of a different VO2 film as it goes through the metal-insulator transition. We find that a two-dimensional Hamiltonian with both interactions and random field disorder is required to explain the intricate, fractal intertwining of metal and insulator domains during the transition. This detailed knowledge about the underlying Hamiltonian paves the way to using the model to control the pattern formation via, e.g., tailored hysteresis protocols. We also introduce a distribution-based confidence measure on the results of a multi-label classifier, which does not rely on adversarial training. In addition, we propose a new machine learning based criterion for diagnosing a physical system's proximity to criticality.


title: High Temperature Ferromagnetism in Cr$_{1+x}$Pt$_{5-x}$P

abstract: We present the growth and basic magnetic and transport properties of Cr$_{1+x}$Pt$_{5-x}$P. We show that single crystals can readily be grown from a high-temperature solution created by adding dilute quantities of Cr to Pt-P based melts. Like other 1-5-1 compounds, Cr$_{1+x}$Pt$_{5-x}$P adopts a tetragonal P4/mmm structure composed face-sharing CrPt$_3$ like slabs that are broken up along the c-axis by sheets of P atoms. EDS and X-ray diffraction measurements both suggest Cr$_{1+x}$Pt$_{5-x}$P has mixed occupancy between Cr and Pt atoms, similar to what is found in the closely related compound CrPt$_3$, giving real compositions of Cr$_{1.5}$Pt$_{4.5}$P (x = 0.5). We report that Cr$_{1.5}$Pt$_{4.5}$P orders ferromagnetically at T$_C$ = 464.5 K with a saturated moment of $\approx$ 2.1 $\mu_{\textit{B}}$/Cr at 1.8 K. Likely owing to the strong spin-orbit coupling associated with the large quantity of high Z Pt atoms, Cr$_{1.5}$Pt$_{4.5}$P has exceptionally strong planar anisotropy with estimated anisotropy fields of 345 kOe and 220 kOe at 1.8 K and 300 K respectively. The resistance of Cr$_{1.5}$Pt$_{4.5}$P has a metallic temperature dependence with relatively weak magnetoresistance. Electronic band structure calculations show that CrPt$_5$P has a large peak in the density of states near the Fermi level which is split into spin majority and minority bands in the ferromagnetic state. Furthermore, the calculations suggest substantial hybridization between Cr-3d and Pt-5d states near the Fermi level, in agreement with the experimentally measured anisotropy.


title: Fast, effective, and coherent time series modeling using the sparsity-ranked lasso

abstract: The sparsity-ranked lasso (SRL) has been developed for model selection and estimation in the presence of interactions and polynomials. The main tenet of the SRL is that an algorithm should be more skeptical of higher-order polynomials and interactions *a priori* compared to main effects, and hence the inclusion of these more complex terms should require a higher level of evidence. In time series, the same idea of ranked prior skepticism can be applied to the possibly seasonal autoregressive (AR) structure of the series during the model fitting process, becoming especially useful in settings with uncertain or multiple modes of seasonality. The SRL can naturally incorporate exogenous variables, with streamlined options for inference and/or feature selection. The fitting process is quick even for large series with a high-dimensional feature set. In this work, we discuss both the formulation of this procedure and the software we have developed for its implementation via the **fastTS** R package. We explore the performance of our SRL-based approach in a novel application involving the autoregressive modeling of hourly emergency room arrivals at the University of Iowa Hospitals and Clinics. We find that the SRL is considerably faster than its competitors, while producing more accurate predictions.


title: Magnetic fields and relativistic electrons fill entire galaxy cluster

abstract: The hot plasma within merging galaxy clusters is predicted to be filled with shocks and turbulence that may convert part of their kinetic energy into relativistic electrons and magnetic fields generating synchrotron radiation. Analyzing Low Frequency Array (LOFAR) observations of the galaxy cluster Abell 2255, we show evidence of radio synchrotron emission distributed over very large scales of at least 5 megaparsec. The pervasive radio emission witnesses that shocks and turbulence efficiently transfer kinetic energy into relativistic particles and magnetic fields in a region that extends up to the cluster outskirts. The strength of the emission requires a magnetic field energy density at least 100 times higher than expected from a simple compression of primordial fields, presumably implying that dynamo operates efficiently also in the cluster periphery. It also suggests that nonthermal components may contribute substantially to the pressure of the intracluster medium in the cluster periphery.


title: Regression Compatible Listwise Objectives for Calibrated Ranking with Binary Relevance

abstract: As Learning-to-Rank (LTR) approaches primarily seek to improve ranking quality, their output scores are not scale-calibrated by design. This fundamentally limits LTR usage in score-sensitive applications. Though a simple multi-objective approach that combines a regression and a ranking objective can effectively learn scale-calibrated scores, we argue that the two objectives are not necessarily compatible, which makes the trade-off less ideal for either of them. In this paper, we propose a practical regression compatible ranking (RCR) approach that achieves a better trade-off, where the two ranking and regression components are proved to be mutually aligned. Although the same idea applies to ranking with both binary and graded relevance, we mainly focus on binary labels in this paper. We evaluate the proposed approach on several public LTR benchmarks and show that it consistently achieves either best or competitive result in terms of both regression and ranking metrics, and significantly improves the Pareto frontiers in the context of multi-objective optimization. Furthermore, we evaluated the proposed approach on YouTube Search and found that it not only improved the ranking quality of the production pCTR model, but also brought gains to the click prediction accuracy. The proposed approach has been successfully deployed in the YouTube production system.


title: On the Edge Derivative of the Normalized Laplacian with Applications to Kemeny's Constant

abstract: In a connected graph, Kemeny's constant gives the expected time of a random walk from an arbitrary vertex $x$ to reach a randomly-chosen vertex $y$. Because of this, Kemeny's constant can be interpreted as a measure of how well a graph is connected. It is generally unknown how the addition or removal of edges affects Kemeny's constant. Inspired by the directional derivative of the normalized Laplacian, we derive the directional derivative of Kemeny's constant for several graph families. In addition, we find sharp bounds for the directional derivative of an eigenvalue of the normalized Laplacian and bounds for the directional derivative of Kemeny's constant.


title: Max Markov Chain

abstract: In this paper, we introduce Max Markov Chain (MMC), a novel representation for a useful subset of High-order Markov Chains (HMCs) with sparse correlations among the states. MMC is parsimony while retaining the expressiveness of HMCs. Even though parameter optimization is generally intractable as with HMC approximate models, it has an analytical solution, better sample efficiency, and the desired spatial and computational advantages over HMCs and approximate HMCs. Simultaneously, efficient approximate solutions exist for this type of chains as we show empirically, which allow MMCs to scale to large domains where HMCs and approximate HMCs would struggle to perform. We compare MMC with HMC, first-order Markov chain, and an approximate HMC model in synthetic domains with various data types to demonstrate that MMC is a valuable alternative for modeling stochastic processes and has many potential applications.


title: Optimizing for periodicity: a model-independent approach to flux crosstalk calibration for superconducting circuits

abstract: Flux tunability is an important engineering resource for superconducting circuits. Large-scale quantum computers based on flux-tunable superconducting circuits face the problem of flux crosstalk, which needs to be accurately calibrated to realize high-fidelity quantum operations. Typical calibration methods either assume that circuit elements can be effectively decoupled and simple models can be applied, or require a large amount of data. Such methods become ineffective as the system size increases and circuit interactions become stronger. Here we propose a new method for calibrating flux crosstalk, which is independent of the underlying circuit model. Using the fundamental property that superconducting circuits respond periodically to external fluxes, crosstalk calibration of N flux channels can be treated as N independent optimization problems, with the objective functions being the periodicity of a measured signal depending on the compensation parameters. We demonstrate this method on a small-scale quantum annealing circuit based on superconducting flux qubits, achieving comparable accuracy with previous methods. We also show that the objective function usually has a nearly convex landscape, allowing efficient optimization.


title: On the Safety of Interpretable Machine Learning: A Maximum Deviation Approach

abstract: Interpretable and explainable machine learning has seen a recent surge of interest. We focus on safety as a key motivation behind the surge and make the relationship between interpretability and safety more quantitative. Toward assessing safety, we introduce the concept of maximum deviation via an optimization problem to find the largest deviation of a supervised learning model from a reference model regarded as safe. We then show how interpretability facilitates this safety assessment. For models including decision trees, generalized linear and additive models, the maximum deviation can be computed exactly and efficiently. For tree ensembles, which are not regarded as interpretable, discrete optimization techniques can still provide informative bounds. For a broader class of piecewise Lipschitz functions, we leverage the multi-armed bandit literature to show that interpretability produces tighter (regret) bounds on the maximum deviation. We present case studies, including one on mortgage approval, to illustrate our methods and the insights about models that may be obtained from deviation maximization.


title: Angle-free cluster robust Ritz value bounds for restarted block eigensolvers

abstract: Convergence rates of block iterations for solving eigenvalue problems typically measure errors of Ritz values approximating eigenvalues. The errors of the Ritz values are commonly bounded in terms of principal angles between the initial or iterative subspace and the invariant subspace associated with the target eigenvalues. Such bounds thus cannot be applied repeatedly as needed for restarted block eigensolvers, since the left- and right-hand sides of the bounds use different terms. They must be combined with additional bounds which could cause an overestimation. Alternative repeatable bounds that are angle-free and depend only on the errors of the Ritz values have been pioneered for Hermitian eigenvalue problems in doi:10.1515/rnam.1987.2.5.371 but only for a single extreme Ritz value. We extend this result to all Ritz values and achieve robustness for clustered eigenvalues by utilizing nonconsecutive eigenvalues. Our new bounds cover the restarted block Lanczos method and its modifications with shift-and-invert and deflation, and are numerically advantageous.


title: Learning to Grasp the Ungraspable with Emergent Extrinsic Dexterity

abstract: A simple gripper can solve more complex manipulation tasks if it can utilize the external environment such as pushing the object against the table or a vertical wall, known as "Extrinsic Dexterity." Previous work in extrinsic dexterity usually has careful assumptions about contacts which impose restrictions on robot design, robot motions, and the variations of the physical parameters. In this work, we develop a system based on reinforcement learning (RL) to address these limitations. We study the task of "Occluded Grasping" which aims to grasp the object in configurations that are initially occluded; the robot needs to move the object into a configuration from which these grasps can be achieved. We present a system with model-free RL that successfully achieves this task using a simple gripper with extrinsic dexterity. The policy learns emergent behaviors of pushing the object against the wall to rotate and then grasp it without additional reward terms on extrinsic dexterity. We discuss important components of the system including the design of the RL problem, multi-grasp training and selection, and policy generalization with automatic curriculum. Most importantly, the policy trained in simulation is zero-shot transferred to a physical robot. It demonstrates dynamic and contact-rich motions with a simple gripper that generalizes across objects with various size, density, surface friction, and shape with a 78% success rate. Videos can be found at https://sites.google.com/view/grasp-ungraspable/.


title: A light-induced Weyl semiconductor-to-metal transition mediated by Peierls instability

abstract: Elemental tellurium is a strongly spin-orbit coupled Peierls-distorted semiconductor whose band structure features topologically protected Weyl nodes. Using time-dependent density functional theory calculations, we show that impulsive optical excitation can be used to transiently control the amplitude of the Peierls distortion, realizing a mechanism to switch tellurium between three states: Weyl semiconductor, Weyl metal and non-Weyl metal. Further, we present experimental evidence of this inverse-Peierls distortion using time-resolved optical second harmonic generation measurements. These results provide a pathway to multifunctional ultrafast Weyl devices and introduce Peierls systems as viable hosts of light-induced topological transitions.


title: Detection and characterization of wind-blown charged sand grains on Titan with the DraGMet/EFIELD experiment on Dragonfly

abstract: The EFIELD instrument is part of the geophysics and meteorology sensor package DraGMet on the Dragonfly mission, which will explore the surface of Titan in the mid-2030s. EFIELD consists of two electrodes designed to passively record the AC electric field at each landing site. The exploration zone of Dragonfly will mostly consist of dune fields, covered with sand grains. Little is known on the properties of these grains, although Cassini-Huygens observations suggest they are mostly made of organic material produced by Titan's atmospheric photochemistry and evolved at the surface. Little is known also about dune formation and in general about the transport of sediments by winds. The latter much depends on inter-particle forces and therefore on how grains are charged by friction. We demonstrate here that the EFIELD experiment can bring new insights on these questions. We have developed a hydrodynamic-electrostatic model to simulate the trajectory of a wind-blown charged sand grain in the vicinity of an idealized EFIELD probe and to predict how such a grain flying close to the probe would affect its potential. We show that, in some conditions, the resulting perturbation will be strong enough to be detected by the EFIELD probe. More specifically, we find that the detection of typical charged wind-blown grains (200 microns) on Titan requires an instrument standard deviation noise inferior to 1mV, though occasional larger grains flying close to one electrode could be detected with a higher noise level. Furthermore, we propose a method to retrieve information on the charge and velocity of wind-blown charged grains detected by the EFIELD experiment. This method well applies for cases where the particle trajectory can be regarded as quasi-linear. We validate our inversion approach on both synthetic and experimental data obtained with a laboratory prototype of the EFIELD experiment.


title: Jensen's and Cantelli's Inequalities with Imprecise Previsions

abstract: We investigate how basic probability inequalities can be extended to an imprecise framework, where (precise) probabilities and expectations are replaced by imprecise probabilities and lower/upper previsions. We focus on inequalities giving information on a single bounded random variable $X$, considering either convex/concave functions of $X$ (Jensen's inequalities) or one-sided bounds such as $(X\geq c)$ or $(X\leq c)$ (Markov's and Cantelli's inequalities). As for the consistency of the relevant imprecise uncertainty measures, our analysis considers coherence as well as weaker requirements, notably $2$-coherence, which proves to be often sufficient. Jensen-like inequalities are introduced, as well as a generalisation of a recent improvement to Jensen's inequality. Some of their applications are proposed: extensions of Lyapunov's inequality and inferential problems. After discussing upper and lower Markov's inequalities, Cantelli-like inequalities are proven with different degrees of consistency for the related lower/upper previsions. In the case of coherent imprecise previsions, the corresponding Cantelli's inequalities make use of Walley's lower and upper variances, generally ensuring better bounds.


title: Redundant Field Survey Data of Cyclotron with Imperfect Median Plane

abstract: An accurate and detailed field map is important for cyclotron beam dynamics studies. During the long history of cyclotron studies, many techniques have been developed by cyclotron pioneers for the treatment of median plane field map. In this paper, we take the TRIUMF 500 MeV cyclotron as an example to study the asymmetric field resulting from the imperfect median plane symmetry. The ``Gordon approach'' and a highly accurate compact finite differentiation method are used to investigate the historical field survey data. The redundancy in the survey data is revealed by the expansion method, which also makes it possible to correct the error in the measurement. Finally, both the azimuthal field $B_\theta$ and the axial gradient of the axial field $dB_z/dz$ in the median plane are corrected using the radial field map $B_r$. The influence of the correction is examined by recalculating the equilibrium orbit properties of the TRIUMF cyclotron. The result shows significantly increased vertical centering errors of the closed orbits. Further simulation study suggests that these centering errors can be reduced to below 1.5 cm by adjusting the trim coils' $B_r$ field within the output limits of our trim coils' power supplies. The error in the measurement field data may explain why the calculated trim coils' settings during the cyclotron commissioning in 1974 encountered difficulty.


title: Implicit Neural Representation as a Differentiable Surrogate for Photon Propagation in a Monolithic Neutrino Detector

abstract: Optical photons are used as signal in a wide variety of particle detectors. Modern neutrino experiments employ hundreds to tens of thousands of photon detectors to observe signal from millions to billions of scintillation photons produced from energy deposition of charged particles. These neutrino detectors are typically large, containing kilotons of target volume, with different optical properties. Modeling individual photon propagation in form of look-up table requires huge computational resources. As the size of a table increases with detector volume for a fixed resolution, this method scales poorly for future larger detectors. Alternative approaches such as fitting a polynomial to the model could address the memory issue, but results in poorer performance. Both look-up table and fitting approaches are prone to discrepancies between the detector simulation and the data collected. We propose a new approach using SIREN, an implicit neural representation with periodic activation functions, to model the look-up table as a 3D scene and reproduces the acceptance map with high accuracy. The number of parameters in our SIREN model is orders of magnitude smaller than the number of voxels in the look-up table. As it models an underlying functional shape, SIREN is scalable to a larger detector. Furthermore, SIREN can successfully learn the spatial gradients of the photon library, providing additional information for downstream applications. Finally, as SIREN is a neural network representation, it is differentiable with respect to its parameters, and therefore tunable via gradient descent. We demonstrate the potential of optimizing SIREN directly on real data, which mitigates the concern of data vs. simulation discrepancies. We further present an application for data reconstruction where SIREN is used to form a likelihood function for photon statistics.


title: Real-space imaging of polar and elastic nano-textures in thin films via inversion of diffraction data

abstract: Exploiting the emerging nanoscale periodicities in epitaxial, single-crystal thin films is an exciting direction in quantum materials science: confinement and periodic distortions induce novel properties. The structural motifs of interest are ferroelastic, ferroelectric, multiferroic, and, more recently, topologically protected magnetization and polarization textures. A critical step towards heterostructure engineering is understanding their nanoscale structure, best achieved through real-space imaging. X-ray Bragg coherent diffractive imaging visualizes sub-picometer crystalline displacements with tens of nanometers spatial resolution. Yet, it is limited to objects spatially confined in all three dimensions and requires highly coherent, laser-like x-rays. Here we lift the confinement restriction by developing real-space imaging of periodic lattice distortions: we combine an iterative phase retrieval algorithm with unsupervised machine learning to invert the diffuse scattering in conventional x-ray reciprocal-space mapping into real-space images of polar and elastic textures in thin epitaxial films. We first demonstrate our imaging in PbTiO3/SrTiO3 superlattices to be consistent with published phase-field model calculations. We then visualize strain-induced ferroelastic domains emerging during the metal-insulator transition in Ca2RuO4 thin films. Instead of homogeneously transforming into a low-temperature structure (like in bulk), the strained Mott insulator splits into nanodomains with alternating lattice constants, as confirmed by cryogenic scanning transmission electron microscopy. Our study reveals the type, size, orientation, and crystal displacement field of the nano-textures. The non-destructive imaging of textures promises to improve models for their dynamics and enable advances in quantum materials and microelectronics.


title: Solving the Kerzman's problem on the sup-norm estimate for $\dbar$ on product domains

abstract: In this paper, the author solves the long term open problem of Kerzman on sup-norm estimate for Cauchy-Riemann equation on polydisc in $n$-dimensional complex space. The problem has been open since 1971. He also extends and solves the problem on a bounded product domain $\Omega^n$, where $\Omega$ either is simply connected with $C^{1,\alpha}$ boundary or satisfies a uniform exterior ball condition with piecewise $C^1$ boundary.


title: Partially-Observable Security Games for Automating Attack-Defense Analysis

abstract: Network systems often contain vulnerabilities that remain unfixed in a network for various reasons, such as the lack of a patch or knowledge to fix them. With the presence of such residual vulnerabilities, the network administrator should properly react to the malicious activities or proactively prevent them, by applying suitable countermeasures that minimize the likelihood of an attack by the attacker. In this paper, we propose a stochastic game-theoretic approach for analyzing network security and synthesizing defense strategies to protect a network. To support analysis under partial observation, where some of the attacker's activities are unobservable or undetectable by the defender, we construct a one-sided partially observable security game and transform it into a perfect game for further analysis. We prove that this transformation is sound for a sub-class of security games and a subset of properties specified in the logic rPATL. We implement a prototype that fully automates our approach, and evaluate it by conducting experiments on a real-life network.


title: From Enriques surface to Artin-Mumford counterexample

abstract: After an Introduction to the themes of Enriques surfaces and Rationality questions, the Artin-Mumford counterexample to Lueroth problem is revisited. A construction of it is given, which is related in an explicit way to the geometry of Enriques surfaces, more precisely to the special family of Reye congruences and their classical geometry.


title: Hopfian wreath products and the stable finiteness conjecture

abstract: We study the Hopf property for wreath products of finitely generated groups, focusing on the case of an abelian base group. Our main result establishes a strong connection between this problem and Kaplansky's stable finiteness conjecture. Namely, the latter holds true if and only if for every finitely generated abelian group $A$ and every finitely generated Hopfian group $\Gamma$ the wreath product $A \wr \Gamma$ is Hopfian. In fact, we characterize precisely when $A \wr \Gamma$ is Hopfian, in terms of the existence of one-sided units in certain matrix algebras over $\mathbb{F}_p[\Gamma]$, for every prime factor $p$ occurring as the order of some element in $A$. A tool in our arguments is the fact that fields of positive characteristic locally embed into matrix algebras over $\mathbb{F}_p$ thus reducing the stable finiteness conjecture to the case of $\mathbb{F}_p$. A further application of this result shows that the validity of Kaplansky's stable finiteness conjecture is equivalent to a version of Gottschalk's surjunctivity conjecture for additive cellular automata.


title: The Cosmic-Ray Induced Sputtering Process On Icy Grains

abstract: In molecular cloud cores, the cosmic ray (CR) induced sputtering via CR ion-icy grain collision is one of the desorption processes for ice molecules from mantles around dust grains. The efficiency of this process depends on the incident CR ion properties as well as the physicochemical character of the ice mantle. Our main objective is the examination of the sputtering efficiency for H$_2$O and CO ices found in molecular cloud cores. In the calculation routine, we consider a multi-dimensional parameter space that consists of thirty CR ion types, five different CR ion energy flux distributions, two separate ice mantle components (pure H$_2$O and CO), three ice formation states, and two sputtering regimes (linear and quadratic). We find that the sputtering behavior of H$_2$O and CO ices is dominated by the quadratic regime rather than the linear regime, especially for CO sputtering. The sputtering rate coefficients for H$_2$O and CO ices show distinct variations with respect to the adopted CR ion energy flux as well as the grain size-dependent mantle depth. The maximum radius of the cylindrical latent region is quite sensitive to the effective electronic stopping power. The track radii for CO ice are much bigger than H$_2$O ice values. In contrast to the H$_2$O mantle, even relatively light CR ions ($Z \geq 4$) may lead to a track formation within the CO mantle, depending on ${\rm S_{\rm e,eff}}$. We suggest that the latent track formation threshold can be assumed as a separator between the linear and the quadratic regimes for sputtering.


title: Convergence of the Inexact Langevin Algorithm and Score-based Generative Models in KL Divergence

abstract: We study the Inexact Langevin Dynamics (ILD), Inexact Langevin Algorithm (ILA), and Score-based Generative Modeling (SGM) when utilizing estimated score functions for sampling. Our focus lies in establishing stable biased convergence guarantees in terms of the Kullback-Leibler (KL) divergence. To achieve these guarantees, we impose two key assumptions: 1) the target distribution satisfies the log-Sobolev inequality (LSI), and 2) the score estimator exhibits a bounded Moment Generating Function (MGF) error. Notably, the MGF error assumption we adopt is more lenient compared to the $L^\infty$ error assumption used in existing literature. However, it is stronger than the $L^2$ error assumption utilized in recent works, which often leads to unstable bounds. We explore the question of how to obtain a provably accurate score estimator that satisfies the MGF error assumption. Specifically, we demonstrate that a simple estimator based on kernel density estimation fulfills the MGF error assumption for sub-Gaussian target distribution, at the population level.


title: Optimizing Fiducial Marker Placement for Improved Visual Localization

abstract: Adding fiducial markers to a scene is a well-known strategy for making visual localization algorithms more robust. Traditionally, these marker locations are selected by humans who are familiar with visual localization techniques. This paper explores the problem of automatic marker placement within a scene. Specifically, given a predetermined set of markers and a scene model, we compute optimized marker positions within the scene that can improve accuracy in visual localization. Our main contribution is a novel framework for modeling camera localizability that incorporates both natural scene features and artificial fiducial markers added to the scene. We present optimized marker placement (OMP), a greedy algorithm that is based on the camera localizability framework. We have also designed a simulation framework for testing marker placement algorithms on 3D models and images generated from synthetic scenes. We have evaluated OMP within this testbed and demonstrate an improvement in the localization rate by up to 20 percent on four different scenes.


title: Creating large Fock states and massively squeezed states in optics using systems with nonlinear bound states in the continuum

abstract: The quantization of the electromagnetic field leads directly to the existence of quantum mechanical states, called Fock states, with an exact integer number of photons. Despite these fundamental states being long-understood, and despite their many potential applications, generating them is largely an open problem. For example, at optical frequencies, it is challenging to deterministically generate Fock states of order two and beyond. Here, we predict the existence of a new effect in nonlinear optics, which enables the deterministic generation of large Fock states at arbitrary frequencies. The effect, which we call an n-photon bound state in the continuum, is one in which a photonic resonance (such as a cavity mode) becomes lossless when a precise number of photons n is inside the resonance. Based on analytical theory and numerical simulations, we show that these bound states enable a remarkable phenomenon in which a coherent state of light, when injected into a system supporting this bound state, can spontaneously evolve into a Fock state of a controllable photon number. This effect is also directly applicable for creating (highly) squeezed states of light, whose photon number fluctuations are (far) below the value expected from classical physics (i.e., shot noise). We suggest several examples of systems to experimentally realize the effects predicted here in nonlinear nanophotonic systems, showing examples of generating both optical Fock states with large n (n > 10), as well as more macroscopic photonic states with very large squeezing, with over 90% less noise (10 dB) than the classical value associated with shot noise.


title: MAST: Multiscale Audio Spectrogram Transformers

abstract: We present Multiscale Audio Spectrogram Transformer (MAST) for audio classification, which brings the concept of multiscale feature hierarchies to the Audio Spectrogram Transformer (AST). Given an input audio spectrogram, we first patchify and project it into an initial temporal resolution and embedding dimension, post which the multiple stages in MAST progressively expand the embedding dimension while reducing the temporal resolution of the input. We use a pyramid structure that allows early layers of MAST operating at a high temporal resolution but low embedding space to model simple low-level acoustic information and deeper temporally coarse layers to model high-level acoustic information with high-dimensional embeddings. We also extend our approach to present a new Self-Supervised Learning (SSL) method called SS-MAST, which calculates a symmetric contrastive loss between latent representations from a student and a teacher encoder, leveraging patch-drop, a novel audio augmentation approach that we introduce. In practice, MAST significantly outperforms AST by an average accuracy of 3.4% across 8 speech and non-speech tasks from the LAPE Benchmark, achieving state-of-the-art results on keyword spotting in Speech Commands. Additionally, our proposed SS-MAST achieves an absolute average improvement of 2.6% over the previously proposed SSAST.


title: Afterglow Polarization from Off-Axis GRB Jets

abstract: As we further our studies on Gamma-ray bursts (GRBs), both on theoretical models and observational tools, more and more options begin to open for exploration of its physical properties. As transient events primarily dominated by synchrotron radiation, it is expected that the synchrotron photons emitted by GRBs should present some degree of polarization throughout the evolution of the burst. Whereas observing this polarization can still be challenging due to the constraints on observational tools, especially for short GRBs, it is paramount that the groundwork is laid for the day we have abundant data. In this work, we present a polarization model linked with an off-axis spreading top-hat jet synchrotron scenario in a stratified environment with a density profile $n(r)\propto r^ {-k}$. We present this model's expected temporal polarization evolution for a realistic set of afterglow parameters constrained within the values observed in the GRB literature for four degrees of stratification $k=0,1,1.5 {\rm \, and\,} 2$ and two magnetic field configurations with high extreme anisotropy. We apply this model and predict polarization from a set of GRBs exhibiting off-axis afterglow emission. In particular, for GRB 170817A, we use the available polarimetric upper limits to rule out the possibility of a extremely anisotropic configuration for the magnetic field.


title: Hydration of a side-chain-free n-type semiconducting ladder polymer driven by electrochemical doping

abstract: We study the organic electrochemical transistors (OECTs) performance of the ladder polymer, poly(benzimidazobenzophenanthroline) (BBL) in an attempt to better understand how an apparently hydrophobic side-chain-free polymer is able to operate as an OECT with favorable redox kinetics in an aqueous environment. We examine two BBLs of different molecular masses from different sources. Both BBLs show significant film swelling during the initial reduction step. By combining electrochemical quartz crystal microbalance (eQCM) gravimetry, in-operando atomic force microscopy (AFM), and both ex-situ and in-operando grazing incidence wide-angle x-ray scattering (GIWAXS), we provide a detailed structural picture of the electrochemical charge injection process in BBL in the absence of any hydrophilic side-chains. Compared with ex-situ measurements, in-operando GIWAXS shows both more swelling upon electrochemical doping than has previously been recognized, and less contraction upon dedoping. The data show that BBL films undergo an irreversible hydration driven by the initial electrochemical doping cycle with significant water retention and lamellar expansion that persists across subsequent oxidation/reduction cycles. This swelling creates a hydrophilic environment that facilitates the subsequent fast hydrated ion transport in the absence of the hydrophilic side-chains used in many other polymer systems. Due to its rigid ladder backbone and absence of hydrophilic side-chains, the primary BBL water uptake does not significantly degrade the crystalline order, and the original dehydrated, unswelled state can be recovered after drying. The combination of doping induced hydrophilicity and robust crystalline order leads to efficient ionic transport and good stability.


title: Bayesian Counterfactual Mean Embeddings and Off-Policy Evaluation

abstract: The counterfactual distribution models the effect of the treatment in the untreated group. While most of the work focuses on the expected values of the treatment effect, one may be interested in the whole counterfactual distribution or other quantities associated to it. Building on the framework of Bayesian conditional mean embeddings, we propose a Bayesian approach for modeling the counterfactual distribution, which leads to quantifying the epistemic uncertainty about the distribution. The framework naturally extends to the setting where one observes multiple treatment effects (e.g. an intermediate effect after an interim period, and an ultimate treatment effect which is of main interest) and allows for additionally modelling uncertainty about the relationship of these effects. For such goal, we present three novel Bayesian methods to estimate the expectation of the ultimate treatment effect, when only noisy samples of the dependence between intermediate and ultimate effects are provided. These methods differ on the source of uncertainty considered and allow for combining two sources of data. Moreover, we generalize these ideas to the off-policy evaluation framework, which can be seen as an extension of the counterfactual estimation problem. We empirically explore the calibration of the algorithms in two different experimental settings which require data fusion, and illustrate the value of considering the uncertainty stemming from the two sources of data.


title: SLICER: Learning universal audio representations using low-resource self-supervised pre-training

abstract: We present a new Self-Supervised Learning (SSL) approach to pre-train encoders on unlabeled audio data that reduces the need for large amounts of labeled data for audio and speech classification. Our primary aim is to learn audio representations that can generalize across a large variety of speech and non-speech tasks in a low-resource un-labeled audio pre-training setting. Inspired by the recent success of clustering and contrasting learning paradigms for SSL-based speech representation learning, we propose SLICER (Symmetrical Learning of Instance and Cluster-level Efficient Representations), which brings together the best of both clustering and contrasting learning paradigms. We use a symmetric loss between latent representations from student and teacher encoders and simultaneously solve instance and cluster-level contrastive learning tasks. We obtain cluster representations online by just projecting the input spectrogram into an output subspace with dimensions equal to the number of clusters. In addition, we propose a novel mel-spectrogram augmentation procedure, k-mix, based on mixup, which does not require labels and aids unsupervised representation learning for audio. Overall, SLICER achieves state-of-the-art results on the LAPE Benchmark \cite{9868132}, significantly outperforming DeLoRes-M and other prior approaches, which are pre-trained on $10\times$ larger of unsupervised data. We will make all our codes available on GitHub.


title: Deep Residual Networks for Gravitational Wave Detection

abstract: Traditionally, gravitational waves are detected with techniques such as matched filtering or unmodeled searches based on wavelets. However, in the case of generic black hole binaries with non-aligned spins, if one wants to explore the whole parameter space, matched filtering can become impractical, which sets severe restrictions on the sensitivity and computational efficiency of gravitational-wave searches. Here, we use a novel combination of machine-learning algorithms and arrive at sensitive distances that surpass traditional techniques in a specific setting. Moreover, the computational cost is only a small fraction of the computational cost of matched filtering. The main ingredients are a 54-layer deep residual network (ResNet), a Deep Adaptive Input Normalization (DAIN), a dynamic dataset augmentation, and curriculum learning, based on an empirical relation for the signal-to-noise ratio. We compare the algorithm's sensitivity with two traditional algorithms on a dataset consisting of a large number of injected waveforms of non-aligned binary black hole mergers in real LIGO O3a noise samples. Our machine-learning algorithm can be used in upcoming rapid online searches of gravitational-wave events in a sizeable portion of the astrophysically interesting parameter space. We make our code, AResGW, and detailed results publicly available at https://github.com/vivinousi/gw-detection-deep-learning .


title: Inferring independent sets of Gaussian variables after thresholding correlations

abstract: We consider testing whether a set of Gaussian variables, selected from the data, is independent of the remaining variables. We assume that this set is selected via a very simple approach that is commonly used across scientific disciplines: we select a set of variables for which the correlation with all variables outside the set falls below some threshold. Unlike other settings in selective inference, failure to account for the selection step leads, in this setting, to excessively conservative (as opposed to anti-conservative) results. Our proposed test properly accounts for the fact that the set of variables is selected from the data, and thus is not overly conservative. To develop our test, we condition on the event that the selection resulted in the set of variables in question. To achieve computational tractability, we develop a new characterization of the conditioning event in terms of the canonical correlation between the groups of random variables. In simulation studies and in the analysis of gene co-expression networks, we show that our approach has much higher power than a ``naive'' approach that ignores the effect of selection.


title: Losses Can Be Blessings: Routing Self-Supervised Speech Representations Towards Efficient Multilingual and Multitask Speech Processing

abstract: Self-supervised learning (SSL) for rich speech representations has achieved empirical success in low-resource Automatic Speech Recognition (ASR) and other speech processing tasks, which can mitigate the necessity of a large amount of transcribed speech and thus has driven a growing demand for on-device ASR and other speech processing. However, advanced speech SSL models have become increasingly large, which contradicts the limited on-device resources. This gap could be more severe in multilingual/multitask scenarios requiring simultaneously recognizing multiple languages or executing multiple speech processing tasks. Additionally, strongly overparameterized speech SSL models tend to suffer from overfitting when being finetuned on low-resource speech corpus. This work aims to enhance the practical usage of speech SSL models towards a win-win in both enhanced efficiency and alleviated overfitting via our proposed S$^3$-Router framework, which for the first time discovers that simply discarding no more than 10\% of model weights via only finetuning model connections of speech SSL models can achieve better accuracy over standard weight finetuning on downstream speech processing tasks. More importantly, S$^3$-Router can serve as an all-in-one technique to enable (1) a new finetuning scheme, (2) an efficient multilingual/multitask solution, (3) a state-of-the-art ASR pruning technique, and (4) a new tool to quantitatively analyze the learned speech representation. We believe S$^3$-Router has provided a new perspective for practical deployment of speech SSL models. Our codes are available at: https://github.com/GATECH-EIC/S3-Router.


title: A weak lensing perspective on nonlinear structure formation with fuzzy dark matter

abstract: We investigate nonlinear structure formation in the fuzzy dark matter (FDM) model in comparison to cold dark matter (CDM) models from a weak lensing perspective using perturbative methods. We use Eulerian perturbation theory (PT) up to fourth order to compute the tree-level matter trispectrum and the one-loop matter spectrum and bispectrum from consistently chosen initial conditions. In addition, we predict the non-linear matter power spectra using $N$-body simulations with CDM and FDM initial conditions. We go on to derive the respective lensing spectra, bispectra and trispectra in CDM and FDM in the context of a Euclid-like weak lensing survey. Finally, we compute the attainable cumulative signal-to-noise ratios and an estimate of the attainable $\chi^2$-functionals for distinguishing FDM from CDM at particle masses $m=10^{-21}$ eV, $m = 10^{-22}$ eV and $m = 10^{-23}$ eV. We find that PT predictions cannot be used to reliably distinguish the three models in a weak lensing survey. Assuming that $N$-body simulations overestimate the late-time small-scale power in the FDM model, future weak lensing survey might be used to distinguish between the FDM and CDM cases up to a mass of $m = 10^{-23}$ eV. However, observations probing the local high-$z$ universe are probably more suited to constrain the FDM mass.


title: Addressing interpersonal harm in online gaming communities: the opportunities and challenges for a restorative justice approach

abstract: Most social media platforms implement content moderation to address interpersonal harms such as harassment. Content moderation relies on offender-centered, punitive approaches, e.g., bans and content removal. We consider an alternative justice framework, restorative justice, which aids victims in healing, supports offenders in repairing the harm, and engages community members in addressing the harm collectively. To assess the utility of restorative justice in addressing online harm, we interviewed 23 users from Overwatch gaming communities, including moderators, victims, and offenders; such communities are particularly susceptible to harm, with nearly three quarters of all online game players suffering from some form of online abuse. We study how the communities currently handle harm cases through the lens of restorative justice and examine their attitudes toward implementing restorative justice processes. Our analysis reveals that cultural, technical, and resource-related obstacles hinder implementation of restorative justice within the existing punitive framework despite online community needs and existing structures to support it. We discuss how current content moderation systems can embed restorative justice goals and practices and overcome these challenges.


title: Exploring quantum geometry created by quantum matter

abstract: Exactly soluble models can serve as excellent tools to explore conceptual issues in non-perturbative quantum gravity. In perturbative approaches, it is only the two radiative modes of the linearized gravitational field that are quantized. The goal of this investigation is to probe the `Coulombic' aspects of quantum geometry that are governed entirely by matter sources. Since there are no gravitational waves in 3 dimensions, 3-d gravity coupled to matter provides an ideal arena for this task. Our analysis will reveal novel aspects of quantum gravity that bring out limitations of classical and semi-classical theories in unforeseen regimes: non-linearities of general relativity can magnify small quantum fluctuations in the matter sector to large effects in the gravitational sector. Finally, this analysis leads to thought experiments that bring out rather starkly why understanding of the nature of physical reality depends sensitively on the theoretical lens with which it is probed. As theories becomes richer, new scales emerge, triggering novel effects that could not be imagined before. The model provides a concise realization of this well-known chain.


title: Dust temperature uncertainties hamper the inference of dust and molecular gas masses from the dust continuum emission of quiescent high-redshift galaxies

abstract: Single flux density measurements at observed-frame sub-millimeter and millimeter wavelengths are commonly used to probe dust and gas masses in galaxies. In this Letter, we explore the robustness of this method to infer dust mass, focusing on quiescent galaxies, using a series of controlled experiments on four massive haloes from the Feedback in Realistic Environments (FIRE) project. Our starting point is four star-forming, central galaxies at seven redshifts between z=1.5 and z=4.5. We generate modified quiescent galaxies that have been quenched for 100Myr, 500Myr, or 1Gyr prior to each of the studied redshifts by re-assigning stellar ages. We derive spectral energy distributions for each fiducial and modified galaxy using radiative transfer. We demonstrate that the dust mass inferred is highly dependent on the assumed dust temperature, T_dust, which is often unconstrained observationally. Motivated by recent work on quiescent galaxies that assumed T_dust~25K, we show that the ratio between dust mass and 1.3mm flux density can be higher than inferred by up to an order of magnitude, due to the considerably lower dust temperatures seen in non star-forming galaxies. This can lead to an underestimation of dust mass (and, when sub-mm flux density is used as a proxy for molecular gas content, gas mass). This underestimation is most severe at higher redshifts, where the observed-frame 1.3mm flux density probes rest-frame wavelengths far from the Rayleigh-Jeans regime, and hence depends super-linearly on dust temperature. We fit relations between ratios of rest-frame far-infrared flux densities and mass-weighted dust temperature that can be used to constrain dust temperatures from observations and hence derive more reliable dust and molecular gas masses.


title: Sensor Control for Information Gain in Dynamic, Sparse and Partially Observed Environments

abstract: We present an approach for autonomous sensor control for information gathering under partially observable, dynamic and sparsely sampled environments that maximizes information about entities present in that space. We describe our approach for the task of Radio-Frequency (RF) spectrum monitoring, where the goal is to search for and track unknown, dynamic signals in the environment. To this end, we extend the Deep Anticipatory Network (DAN) Reinforcement Learning (RL) framework by (1) improving exploration in sparse, non-stationary environments using a novel information gain reward, and (2) scaling up the control space and enabling the monitoring of complex, dynamic activity patterns using hybrid convolutional-recurrent neural layers. We also extend this problem to situations in which sampling from the intended RF spectrum/field is limited and propose a model-based version of the original RL algorithm that fine-tunes the controller via a model that is iteratively improved from the limited field sampling. Results in simulated RF environments of differing complexity show that our system outperforms the standard DAN architecture and is more flexible and robust than baseline expert-designed agents. We also show that it is adaptable to non-stationary emission environments.


title: Fair and Optimal Classification via Post-Processing

abstract: To mitigate the bias exhibited by machine learning models, fairness criteria can be integrated into the training process to ensure fair treatment across all demographics, but it often comes at the expense of model performance. Understanding such tradeoffs, therefore, underlies the design of fair algorithms. To this end, this paper provides a complete characterization of the inherent tradeoff of demographic parity on classification problems, under the most general multi-group, multi-class, and noisy setting. Specifically, we show that the minimum error rate achievable by randomized and attribute-aware fair classifiers is given by the optimal value of a Wasserstein-barycenter problem. On the practical side, our findings lead to a simple post-processing algorithm that derives fair classifiers from score functions, which yields the optimal fair classifier when the score is Bayes optimal. We provide suboptimality analysis and sample complexity for our algorithm, and demonstrate its effectiveness on benchmark datasets.


title: Truncated smooth function spaces

abstract: We introduce truncated Besov and Triebel--Lizorkin function spaces and investigate their main properties: embeddings, interpolation, duality, lifting, traces. These new scales allow us to improve several known results in functional analysis and PDE's.


title: Theory of q-commuting contractions : decompositions, dilations and Brehmer's positivity

abstract: In the literature, we have several results associated with canonical decomposition of commuting contractions, dilations and Brehmer's positivity. In this paper, we generalize a few of these results to q-commuting contractions. Here we mainly deal with q-commuting and doubly q-commuting contractions when $|q|=1$ or when $q$ is a unitary operator.


title: Adaptive sparse grid discontinuous Galerkin method: review and software implementation

abstract: This paper reviews the adaptive sparse grid discontinuous Galerkin (aSG-DG) method for computing high dimensional partial differential equations (PDEs) and its software implementation. The C\texttt{++} software package called AdaM-DG, implementing the aSG-DG method, is available on Github at \url{https://github.com/JuntaoHuang/adaptive-multiresolution-DG}. The package is capable of treating a large class of high dimensional linear and nonlinear PDEs. We review the essential components of the algorithm and the functionality of the software, including the multiwavelets used, assembling of bilinear operators, fast matrix-vector product for data with hierarchical structures. We further demonstrate the performance of the package by reporting numerical error and CPU cost for several benchmark test, including linear transport equations, wave equations and Hamilton-Jacobi equations.


title: A 16 Hour Transit of Kepler-167 e Observed by the Ground-based Unistellar Telescope Network

abstract: More than 5,000 exoplanets have been confirmed and among them almost 4,000 were discovered by the transit method. However, few transiting exoplanets have an orbital period greater than 100 days. Here we report a transit detection of Kepler-167 e, a "Jupiter analog" exoplanet orbiting a K4 star with a period of 1,071 days, using the Unistellar ground-based telescope network. From 2021 November 18 to 20, citizen astronomers located in nine different countries gathered 43 observations, covering the 16 hour long transit. Using a nested sampling approach to combine and fit the observations, we detected the mid-transit time to be UTC 2021 November 19 17:20:51 with a 1$\sigma$ uncertainty of 9.8 minutes, making it the longest-period planet to ever have its transit detected from the ground. This is the fourth transit detection of Kepler-167 e, but the first made from the ground. This timing measurement refines the orbit and keeps the ephemeris up to date without requiring space telescopes. Observations like this demonstrate the capabilities of coordinated networks of small telescopes to identify and characterize planets with long orbital periods.


title: The Chern Sectional Curvature of a Hermitian Manifold

abstract: On a Hermitian manifold, the Chern connection can induce a metric connection on the background Riemannian manifold. We call the sectional curvature of the metric connection induced by the Chern connection the Chern sectional curvature of this Hermitian manifold. First, we derive expression of the Chern sectional curvature in local complex coordinates. As an application, we find that a Hermitian metric is K\"ahler if the Riemann sectional curvature and the Chern sectional curvature coincide. As subsequent results, Ricci curvature and scalar curvature of the metric connection induced by the Chern connection are obtained.


title: Constructing the Milky Way Stellar Halo in the Galactic Center by Direct Orbit Integration

abstract: The halo stars on highly radial orbits should inevitably pass the center regions of the Milky Way. Under the assumption that the stellar halo is in dynamical equilibrium and axisymmetric, we integrate the orbits of $\sim 10,000$ halo K-giants at $5\leq r \leq 50$ kpc cross-matched from LAMOST DR5 and $Gaia$ DR3. By carefully considering the selection function, we construct the stellar halo distribution at the entire regions of $r \leq 50$ kpc. We find that a double-broken power-law function well describes the stellar halo density distribution with shallower slopes in the inner regions and the two breaks at $r=10$ kpc and $r=25$ kpc, respectively. The stellar halo becomes flatter from outer to inner regions but has $q\sim 0.5$ at $r \lesssim 5$ kpc. The stellar halo becomes isotropic with a slight prograde rotation in the inner 5 kpc, and reaches velocity dispersions of $\sim 250\rm \ km\ s^{-1}$. We get a weak negative metallicity gradient of $-0.005$ dex kpc$^{-1}$ at $5\leq r \leq 50$ kpc, while there is an excess of relative metal-rich stars with [Fe/H]$>-1$ in the inner 10 kpc. The halo interlopers at $r \leq 5$ kpc from integration of our sample has a mass of $\sim1.2 \times 10^8\ M_{\odot}$ ($\sim 4.7 \times 10^7\ M_{\odot}$ at [Fe/H]$<-1.5$), which can explain 50-100% of the metal-poor stars with [Fe/H]$<-1.5$ directly observed in the Galactic central regions.


title: Reliable Malware Analysis and Detection using Topology Data Analysis

abstract: Increasingly, malwares are becoming complex and they are spreading on networks targeting different infrastructures and personal-end devices to collect, modify, and destroy victim information. Malware behaviors are polymorphic, metamorphic, persistent, able to hide to bypass detectors and adapt to new environments, and even leverage machine learning techniques to better damage targets. Thus, it makes them difficult to analyze and detect with traditional endpoint detection and response, intrusion detection and prevention systems. To defend against malwares, recent work has proposed different techniques based on signatures and machine learning. In this paper, we propose to use an algebraic topological approach called topological-based data analysis (TDA) to efficiently analyze and detect complex malware patterns. Next, we compare the different TDA techniques (i.e., persistence homology, tomato, TDA Mapper) and existing techniques (i.e., PCA, UMAP, t-SNE) using different classifiers including random forest, decision tree, xgboost, and lightgbm. We also propose some recommendations to deploy the best-identified models for malware detection at scale. Results show that TDA Mapper (combined with PCA) is better for clustering and for identifying hidden relationships between malware clusters compared to PCA. Persistent diagrams are better to identify overlapping malware clusters with low execution time compared to UMAP and t-SNE. For malware detection, malware analysts can use Random Forest and Decision Tree with t-SNE and Persistent Diagram to achieve better performance and robustness on noised data.


title: The Andr\'e-Quillen cohomology of commutative monoids

abstract: We observe that Beck modules for a commutative monoid are exactly modules over a graded commutative ring associated to the monoid. Under this identification, the Quillen cohomology of commutative monoids is a special case of Andr\'e-Quillen cohomology for graded commutative rings, generalizing a result of Kurdiani and Pirashvili. To verify this we develop the necessary grading formalism. The partial cochain complex developed by Pierre Grillet for computing Quillen cohomology appears as the start of a modification of the Harrison cochain complex suggested by Michael Barr.


title: Stochastic Treatment Choice with Empirical Welfare Updating

abstract: This paper proposes a novel method to estimate individualised treatment assignment rules. The method is designed to find rules that are stochastic, reflecting uncertainty in estimation of an assignment rule and about its welfare performance. Our approach is to form a prior distribution over assignment rules, not over data generating processes, and to update this prior based upon an empirical welfare criterion, not likelihood. The social planner then assigns treatment by drawing a policy from the resulting posterior. We show analytically a welfare-optimal way of updating the prior using empirical welfare; this posterior is not feasible to compute, so we propose a variational Bayes approximation for the optimal posterior. We characterise the welfare regret convergence of the assignment rule based upon this variational Bayes approximation, showing that it converges to zero at a rate of ln(n)/sqrt(n). We apply our methods to experimental data from the Job Training Partnership Act Study to illustrate the implementation of our methods.


title: $D^2$SLAM: Decentralized and Distributed Collaborative Visual-inertial SLAM System for Aerial Swarm

abstract: Collaborative simultaneous localization and mapping (CSLAM) is essential for autonomous aerial swarms, laying the foundation for downstream algorithms such as planning and control. To address existing CSLAM systems' limitations in relative localization accuracy, crucial for close-range UAV collaboration, this paper introduces $D^2$SLAM-a novel decentralized and distributed CSLAM system. $D^2$SLAM innovatively manages near-field estimation for precise relative state estimation in proximity and far-field estimation for consistent global trajectories. Its adaptable front-end supports both stereo and omnidirectional cameras, catering to various operational needs and overcoming field-of-view challenges in aerial swarms. Experiments demonstrate $D^2$SLAM's effectiveness in accurate ego-motion estimation, relative localization, and global consistency. Enhanced by distributed optimization algorithms, $D^2$SLAM exhibits remarkable scalability and resilience to network delays, making it well-suited for a wide range of real-world aerial swarm applications. The adaptability and proven performance of $D^2$SLAM represent a significant advancement in autonomous aerial swarm technology.


title: Conformal Prediction for STL Runtime Verification

abstract: We are interested in predicting failures of cyber-physical systems during their operation. Particularly, we consider stochastic systems and signal temporal logic specifications, and we want to calculate the probability that the current system trajectory violates the specification. The paper presents two predictive runtime verification algorithms that predict future system states from the current observed system trajectory. As these predictions may not be accurate, we construct prediction regions that quantify prediction uncertainty by using conformal prediction, a statistical tool for uncertainty quantification. Our first algorithm directly constructs a prediction region for the satisfaction measure of the specification so that we can predict specification violations with a desired confidence. The second algorithm constructs prediction regions for future system states first, and uses these to obtain a prediction region for the satisfaction measure. To the best of our knowledge, these are the first formal guarantees for a predictive runtime verification algorithm that applies to widely used trajectory predictors such as RNNs and LSTMs, while being computationally simple and making no assumptions on the underlying distribution. We present numerical experiments of an F-16 aircraft and a self-driving car.