title: Thermal history of matrix forsterite grains from Murchison based on high-resolution tomography

abstract: Protoplanetary disks are dust- and gas-rich structures surrounding protostars. Depending on the distance from the protostar, this dust is thermally processed to different degrees and accreted to form bodies of varying chemical compositions. The primordial accretion processes occurring in the early protoplanetary disk such as chondrule formation and metal segregation are not well understood. One way to constrain them is to study the morphology and composition of forsteritic grains from the matrix of carbonaceous chondrites. Here, we present high-resolution ptychographic X-ray nanotomography and multimodal chemical micro-tomography (X-ray diffraction and X-ray fluorescence) to reveal the early history of forsteritic grains extracted from the matrix of the Murchison CM2.5 chondrite. The 3D electron density maps revealed, at unprecedented resolution (64~nm), spherical inclusions containing Fe-Ni, very little silica-rich glass and void caps (i.e., volumes where the electron density is consistent with conditions close to vacuum) trapped in forsterite. The presence of the voids along with the overall composition, petrological textures, and shrinkage calculations is consistent with the grains experiencing one or more heating events with peak temperatures close to the melting point of forsterite ($\sim$2100~K) and subsequently cooled and contracted, in agreement with chondrule-forming conditions.

categories: astro-ph.EP astro-ph.GA astro-ph.IM physics.geo-ph


title: Hyperbolic domains in real Euclidean spaces

abstract: The second named author and David Kalaj introduced a pseudometric on any domain in the real Euclidean space $\mathbb R^n$, $n\ge 3$, defined in terms of conformal harmonic discs, by analogy with Kobayashi's pseudometric on complex manifolds, which is defined in terms of holomorphic discs. They showed that on the unit ball of $\mathbb R^n$, the minimal metric coincides with the classical Beltrami-Cayley-Klein metric, one of the models of hyperbolic geometry. In the present paper we investigate properties of the minimal pseudometric and give sufficient conditions for a domain to be (complete) hyperbolic, meaning that the minimal pseudometric is a (complete) metric. We show in particular that a convex domain is complete hyperbolic if and only if it does not contain any affine 2-planes. One of our main results is that a domain with a negative minimal plurisubharmonic exhaustion function is hyperbolic, and a bounded strongly minimally convex domain is complete hyperbolic. We also prove a localization theorem for the minimal pseudometric.


title: Minimum Path Star Topology Algorithms for Weighted Regions and Obstacles

abstract: Shortest path algorithms have played a key role in the past century, paving the way for modern day GPS systems to find optimal routes along static systems in fractions of a second. One application of these algorithms includes optimizing the total distance of power lines (specifically in star topological configurations). Due to the relevancy of discovering well-connected electrical systems in certain areas, finding a minimum path that is able to account for geological features would have far-reaching consequences in lowering the cost of electric power transmission. We initialize our research by proving the convex hull as an effective bounding mechanism for star topological minimum path algorithms. Building off this bounding, we propose novel algorithms to manage certain cases that lack existing methods (weighted regions and obstacles) by discretizing Euclidean space into squares and combining pre-existing algorithms that calculate local minimums that we believe have a possibility of being the absolute minimum. We further designate ways to evaluate iterations necessary to reach some level of accuracy. Both of these novel algorithms fulfill certain niches that past literature does not cover.


title: The supersymmetric extension of the replica model

abstract: We perform a $\mathcal{N}=1$ supersymmetric extension of the replica model quantized in the Landau gauge and compute the gluon and gluino propagators at tree-level, such results display a supersymmetric confined model very similar to the supersymmetric version of the Gribov-Zwanziger approach.


title: Inverse magnetocaloric effect in synthetic antiferromagnets

abstract: The magnetocaloric effect in exchange-coupled synthetic-antiferromagnet multilayers is investigated experimentally and theoretically. We observe a temperature-controlled inversion of the effect, where the entropy increases on switching the individual ferromagnetic layers from anti-parallel to parallel alignment near their Curie point. Using a microscopic analytical model as well as numerical atomistic-spin simulations of the system, we explain the observed effect as due to the interplay between the intra- and inter-layer exchange interactions, which either add up or counteract to effectively modulate the Curie temperature of the dilute ferromagnetic layers. The proposed method of designing tunable, strongly magneto-caloric materials should be of interest for such applications as heat-assisted spintronics and magnetic refrigeration.


title: Exploring the role of binarity in the origin of the bimodal rotational velocity distribution in stellar clusters

abstract: Many young and intermediate age massive stellar clusters host bimodal distributions in the rotation rates of their stellar populations, with a dominant peak of rapidly rotating stars and a secondary peak of slow rotators. The origin of this bimodal rotational distribution is currently debated and two main theories have been put forward in the literature. The first is that all/most stars are born as rapid rotators and that interacting binaries brake a fraction of the stars, resulting in two populations. The second is that the rotational distribution is a reflection of the early evolution of pre-main sequence stars, in particular, whether they are able to retain or lose their protoplanetary discs during the first few Myr. Here, we test the binary channel by exploiting multi-epoch VLT/MUSE observations of NGC 1850, a 100Myr massive cluster in the LMC, to search for differences in the binary fractions of the slow and fast rotating populations. If binarity is the cause of the rotational bimodality, we would expect that the slowly rotating population should have a much larger binary fraction than the rapid rotators. However, in our data we detect similar fractions of binary stars in the slow and rapidly rotating populations (5.9+/-1.1% and 4.5+/-0.6%, respectively).Hence, we conclude that binarity is not a dominant mechanism in the formation of the observed bimodal rotational distributions.


title: Generating diffusions with fractional Brownian motion

abstract: We study fast / slow systems driven by a fractional Brownian motion $B$ with Hurst parameter $H\in (\frac 13, 1]$. Surprisingly, the slow dynamic converges on suitable timescales to a limiting Markov process and we describe its generator. More precisely, if $Y^\varepsilon$ denotes a Markov process with sufficiently good mixing properties evolving on a fast timescale $\varepsilon \ll 1$, the solutions of the equation $$ dX^\varepsilon = \varepsilon^{\frac 12-H} F(X^\varepsilon,Y^\varepsilon)\,dB+F_0(X^\varepsilon,Y^\varepsilon)\,dt\; $$ converge to a regular diffusion without having to assume that $F$ averages to $0$, provided that $H< \frac 12$. For $H > \frac 12$, a similar result holds, but this time it does require $F$ to average to $0$. We also prove that the $n$-point motions converge to those of a Kunita type SDE. One nice interpretation of this result is that it provides a continuous interpolation between the homogenisation theorem for random ODEs with rapidly oscillating right-hand sides ($H=1$) and the averaging of diffusion processes ($H= \frac 12$).


title: Targeted Cross-Validation

abstract: In many applications, we have access to the complete dataset but are only interested in the prediction of a particular region of predictor variables. A standard approach is to find the globally best modeling method from a set of candidate methods. However, it is perhaps rare in reality that one candidate method is uniformly better than the others. A natural approach for this scenario is to apply a weighted $L_2$ loss in performance assessment to reflect the region-specific interest. We propose a targeted cross-validation (TCV) to select models or procedures based on a general weighted $L_2$ loss. We show that the TCV is consistent in selecting the best performing candidate under the weighted $L_2$ loss. Experimental studies are used to demonstrate the use of TCV and its potential advantage over the global CV or the approach of using only local data for modeling a local region. Previous investigations on CV have relied on the condition that when the sample size is large enough, the ranking of two candidates stays the same. However, in many applications with the setup of changing data-generating processes or highly adaptive modeling methods, the relative performance of the methods is not static as the sample size varies. Even with a fixed data-generating process, it is possible that the ranking of two methods switches infinitely many times. In this work, we broaden the concept of the selection consistency by allowing the best candidate to switch as the sample size varies, and then establish the consistency of the TCV. This flexible framework can be applied to high-dimensional and complex machine learning scenarios where the relative performances of modeling procedures are dynamic.


title: Automatically Exposing Problems with Neural Dialog Models

abstract: Neural dialog models are known to suffer from problems such as generating unsafe and inconsistent responses. Even though these problems are crucial and prevalent, they are mostly manually identified by model designers through interactions. Recently, some research instructs crowdworkers to goad the bots into triggering such problems. However, humans leverage superficial clues such as hate speech, while leaving systematic problems undercover. In this paper, we propose two methods including reinforcement learning to automatically trigger a dialog model into generating problematic responses. We show the effect of our methods in exposing safety and contradiction issues with state-of-the-art dialog models.


title: Understanding two same-sign and three leptons with $b$-jets in four top quark events at the LHC

abstract: The top quark is the heaviest known elementary particle of the Standard Model (SM) of particle physics and, therefore, it is expected to have large couplings to hypothetical new physics in many models beyond the SM (BSM). Various studies have predicted the presence of multi-lepton anomalies at the LHC. One of those anomalies is the excess production of two same-sign leptons and three isolated leptons in association with $b$-jets. These are reasonably well described by a 2HDM+$S$ model, where $S$ is a singlet scalar. Both the ATLAS and CMS experiments have reported sustained excesses in these final states. This includes corners of the phase-space where production of top quark pairs in association with a $W$ boson contributes to. Here, we investigate the production of two same-sign and three leptons from the production of four top quark final states. Our focus is on understanding the differences between the SM and BSM production mechanisms of four top quarks from $t\overline{t} A$ ($A \rightarrow t\overline{t}$) using Machine Leaning techniques with twelve discriminating kinematic variables.


title: Residual Adapters for Parameter-Efficient ASR Adaptation to Atypical and Accented Speech

abstract: Automatic Speech Recognition (ASR) systems are often optimized to work best for speakers with canonical speech patterns. Unfortunately, these systems perform poorly when tested on atypical speech and heavily accented speech. It has previously been shown that personalization through model fine-tuning substantially improves performance. However, maintaining such large models per speaker is costly and difficult to scale. We show that by adding a relatively small number of extra parameters to the encoder layers via so-called residual adapter, we can achieve similar adaptation gains compared to model fine-tuning, while only updating a tiny fraction (less than 0.5%) of the model parameters. We demonstrate this on two speech adaptation tasks (atypical and accented speech) and for two state-of-the-art ASR architectures.


title: The converse of Sturm's Separation Theorem

abstract: We show that Sturm's classical separation theorem on the interlacing of the zeros of linearly independent solutions of real second order two-term ordinary differential equations necessarily fails in the presence of a unique turning point in the principal part of the equation. Related results are discussed. The last section contains an extension of the main result to a finite number of turning points.


title: Simultaneous Successive Twinning Captured by Atomic Electron Tomography

abstract: Shape-controlled synthesis of multiply twinned nanostructures is heavily emphasized in nanoscience, in large part due to the desire to control the size, shape, and terminating facets of metal nanoparticles for applications in catalysis. Direct control of the size and shape of solution-grown nanoparticles relies on an understanding of how synthetic parameters alter nanoparticle structures during synthesis. However, while outcome populations can be effectively studied with standard electron microscopy methods, transient structures that appear during some synthetic routes are difficult to study using conventional high resolution imaging methods due to the high complexity of the 3D nanostructures. Here, we have studied the prevalence of transient structures during growth of multiply twinned particles and employed atomic electron tomography to reveal the atomic-scale three dimensional structure of a Pd nanoparticle undergoing a shape transition. By identifying over 20,000 atoms within the structure and classifying them according to their local crystallographic environment, we observe a multiply-twinned structure consistent with a simultaneous successive twinning from a decahedral to icosahedral structure. We also observe a high degree of structural disorder in the nanoparticle, and a disordered crystal structure on the particle surface. Our results shed light on the transition mechanism for formation of icosahedral nanoparticle structures.


title: Learning trends of COVID-19 using semi-supervised clustering

abstract: A finite mixture model is used to learn trends from the currently available data on coronavirus (COVID-19). Data on the number of confirmed COVID-19 related cases and deaths for European countries and the United States (US) are explored. A semi-supervised clustering approach with positive equivalence constraints is used to incorporate country and state information into the model. The analysis of trends in the rates of cases and deaths is carried out jointly using a mixture of multivariate Gaussian non-linear regression models with a mean trend specified using a generalized logistic function. The optimal number of clusters is chosen using the Bayesian information criterion. The resulting clusters provide insight into different mitigation strategies adopted by US states and European countries. The obtained results help identify the current relative standing of individual states and show a possible future if they continue with the chosen mitigation technique


title: A fast, high-order numerical method for the simulation of single-excitation states in quantum optics

abstract: We consider the numerical solution of a nonlocal partial differential equation which models the process of collective spontaneous emission in a two-level atomic system containing a single photon. We reformulate the problem as an integro-differential equation for the atomic degrees of freedom, and describe an efficient solver for the case of a Gaussian atomic density. The problem of history dependence arising from the integral formulation is addressed using sum-of-exponentials history compression. We demonstrate the solver on two systems of physical interest: in the first, an initially-excited atom decays into a photon by spontaneous emission, and in the second, a photon pulse is used to an excite an atom, which then decays.


title: Critical Points in Quantum Generative Models

abstract: One of the most important properties of neural networks is the clustering of local minima of the loss function near the global minimum, enabling efficient training. Though generative models implemented on quantum computers are known to be more expressive than their traditional counterparts, it has empirically been observed that these models experience a transition in the quality of their local minima. Namely, below some critical number of parameters, all local minima are far from the global minimum in function value; above this critical parameter count, all local minima are good approximators of the global minimum. Furthermore, for a certain class of quantum generative models, this transition has empirically been observed to occur at parameter counts exponentially large in the problem size, meaning practical training of these models is out of reach. Here, we give the first proof of this transition in trainability, specializing to this latter class of quantum generative model. We use techniques inspired by those used to study the loss landscapes of classical neural networks. We also verify that our analytic results hold experimentally even at modest model sizes.


title: Probabilistic Analysis of Euclidean Capacitated Vehicle Routing

abstract: We give a probabilistic analysis of the unit-demand Euclidean capacitated vehicle routing problem in the random setting, where the input distribution consists of $n$ unit-demand customers modeled as independent, identically distributed uniform random points in the two-dimensional plane. The objective is to visit every customer using a set of routes of minimum total length, such that each route visits at most $k$ customers, where $k$ is the capacity of a vehicle. All of the following results are in the random setting and hold asymptotically almost surely. The best known polynomial-time approximation for this problem is the iterated tour partitioning (ITP) algorithm, introduced in 1985 by Haimovich and Rinnooy Kan. They showed that the ITP algorithm is near-optimal when $k$ is either $o(\sqrt{n})$ or $\omega(\sqrt{n})$, and they asked whether the ITP algorithm was also effective in the intermediate range. In this work, we show that when $k=\sqrt{n}$, the ITP algorithm is at best a $(1+c_0)$-approximation for some positive constant $c_0$. On the other hand, the approximation ratio of the ITP algorithm was known to be at most $0.995+\alpha$ due to Bompadre, Dror, and Orlin, where $\alpha$ is the approximation ratio of an algorithm for the traveling salesman problem. In this work, we improve the upper bound on the approximation ratio of the ITP algorithm to $0.915+\alpha$. Our analysis is based on a new lower bound on the optimal cost for the metric capacitated vehicle routing problem, which may be of independent interest.


title: Topological invariants beyond symmetry indicators: Boundary diagnostics for twofold rotationally symmetric superconductors

abstract: Topological crystalline superconductors are known to have possible higher-order topology, which results in Majorana modes on $d-2$ or lower-dimensional boundaries. Given the rich possibilities of boundary signatures, it is desirable to have topological invariants that can predict the type of Majorana modes from band structures. Although symmetry indicators, a type of invariant that depends only on the band data at high-symmetry points, have been proposed for certain crystalline superconductors, there exist symmetry classes in which symmetry indicators fail to distinguish superconductors with different Majorana boundaries. Here, we systematically obtain topological invariants for an example of this kind, two-dimensional time-reversal symmetric superconductors with twofold rotational symmetry $C_2$. First, we show that the nontrivial topology is independent of band data on the high-symmetry points by conducting a momentum-space classification study. Then, from the resulting K groups, we derive calculable expressions for four $\mathbb{Z}_2$ invariants defined on high-symmetry lines or general points in the Brillouin zone. Finally, together with a real-space classification study, we establish the bulk-boundary correspondence and show that the four $\mathbb{Z}_2$ invariants can predict Majorana boundary types from band structures. Our proposed invariants can fuel practical material searches for $C_2$-symmetric topological superconductors featuring Majorana edge and corner modes.


title: OSPREI: A Coupled Approach to Modeling CME-Driven Space Weather with Automatically-Generated, User-Friendly Outputs

abstract: Coronal Mass Ejections (CMEs) drive space weather activity at Earth and throughout the solar system. Current CME-related space weather predictions rely on information reconstructed from coronagraphs, sometimes from only a single viewpoint, to drive a simple interplanetary propagation model, which only gives the arrival time or limited additional information. We present the coupling of three established models into OSPREI (Open Solar Physics Rapid Ensemble Information), a new tool that describes Sun-to-Earth CME behavior, including the location, orientation, size, shape, speed, arrival time, and internal thermal and magnetic properties, on the timescale needed for forecasts. First, ForeCAT describes the trajectory that a CME takes through the solar corona. Second, ANTEATR simulates the propagation, including expansion and deformation, of a CME in interplanetary space and determines the evolution of internal properties via conservation laws. Finally, FIDO produces in situ profiles for a CME's interaction with a synthetic spacecraft. OSPREI includes ensemble modeling by varying each input parameter to probe any uncertainty in their values, yielding probabilities for all outputs. Standardized visualizations are automatically generated, providing easily-accessible, essential information for space weather forecasting. We show OSPREI results for CMEs observed in the corona on 2021 April 22 and 2021 May 09. We approach these CME as a forecasting proof-of-concept, using information analogous to what would be available in real time rather than fine-tuning input parameters to achieve a best fit for a detailed scientific study. The OSPREI prediction shows good agreement with the arrival time and in situ properties.


title: Multihop: Leveraging Complex Models to Learn Accurate Simple Models

abstract: Knowledge transfer from a complex high performing model to a simpler and potentially low performing one in order to enhance its performance has been of great interest over the last few years as it finds applications in important problems such as explainable artificial intelligence, model compression, robust model building and learning from small data. Known approaches to this problem (viz. Knowledge Distillation, Model compression, ProfWeight, etc.) typically transfer information directly (i.e. in a single/one hop) from the complex model to the chosen simple model through schemes that modify the target or reweight training examples on which the simple model is trained. In this paper, we propose a meta-approach where we transfer information from the complex model to the simple model by dynamically selecting and/or constructing a sequence of intermediate models of decreasing complexity that are less intricate than the original complex model. Our approach can transfer information between consecutive models in the sequence using any of the previously mentioned approaches as well as work in 1-hop fashion, thus generalizing these approaches. In the experiments on real data, we observe that we get consistent gains for different choices of models over 1-hop, which on average is more than 2\% and reaches up to 8\% in a particular case. We also empirically analyze conditions under which the multi-hop approach is likely to be beneficial over the traditional 1-hop approach, and report other interesting insights. To the best of our knowledge, this is the first work that proposes such a multi-hop approach to perform knowledge transfer given a single high performing complex model, making it in our opinion, an important methodological contribution.


title: A doubly monotone flow for constant width bodies in $\mathbb{R}^3$

abstract: We introduce a flow in the space of constant width bodies in three-dimensional Euclidean space that simultaneously increases the volume and decreases the circumradius of the shape as time increases. Starting from any initial constant width figure, we show that the flow exists for all positive times and converges to a closed ball as time tends to plus infinity. We also anticipate that this flow is interesting to study for negative times and that it would provide a mechanism to decrease the volume and increase the circumradius of a constant width body.


title: Exosphere Modeling of Proxima b: A Case Study of Photochemical Escape with a Venus-like Atmosphere

abstract: Exoplanets orbiting M-dwarfs within habitable zones are exposed to stellar environments more extreme than that terrestrial planets experience in our Solar System, which can significantly impact the atmospheres of the exoplanets and affect their habitability and sustainability. This study provides the first prediction of hot oxygen corona structure and the associated photochemical loss from a 1 bar CO2-dominated atmosphere of a Venus-like rocky exoplanet, where dissociative recombination of O2+ ions is assumed to be the major source reaction for the escape of neutral O atoms and formation of the hot O corona (or exospheres) as on Mars and Venus. We employ a 3D Monte Carlo code to simulate the exosphere of Proxima Centauri b (PCb) based on the ionosphere simulated by a 3D magnetohydrodynamic model. Our simulation results show that variability of the stellar wind dynamic pressure over one orbital period of PCb does not affect the overall spatial structure of the hot O corona but contributes to the change in the global hot O escape rate that varies by an order of magnitude. The escape increases dramatically when the planet possesses its intrinsic magnetic fields as the ionosphere becomes more extended with the presence of a global magnetic field. The extended hot O corona may lead to a more extended H exosphere through collisions between thermal H and hot O, which exemplifies the importance of considering nonthermal populations in exospheres to interpret future observations.


title: Photo-memristive sensing with charge storing 2D carbon nitrides

abstract: We report the charge storing 2D carbon nitride potassium poly(heptazine imide), K-PHI, as a direct memristive (bio)sensing platform. Memristive devices have the potential to innovate current (bio)electronic systems such as photo-electrochemical sensors by incorporating new sensing capabilities including non-invasive, wireless remote and time-delayed (memory) readout. We demonstrate a direct photomemristive sensing platform that capitalizes on K PHI's visible light bandgap, large oxidation potential and intrinsic optoionic light energy storage properties. Our system simultaneously enables analyte concentration information storage as well as potentiometric, impedimetric and coulo-metric readouts on the same material, with no additional reagents required. Utilizing the light-induced charge storage function of K-PHI, we demonstrate analyte sensing via charge accumulation and present various methods to write/erase this information from the material. Additionally, fully wireless colorimetric and fluorometric detection of the charged state of K-PHI is demonstrated and could facilitate its use as particle-based in-situ sensing probe. The various readout options of the K PHI's response enable us to adapt the sensitivities and dynamic ranges without modifying the sensor. We demonstrate these features using glucose as an example analyte over a wide range of concentrations (50 $\mu$M to 50 mM). Moreover, due to the strong oxidative power of K-PHI, this sensing platform is able to detect a large variety of organic or biologically relevant analytes. Since PHI is easily synthesized, based on earth abundant precursors, biocompatible, chemically robust and responsive to visible light, we anticipate that the sensing platform presented herein opens up novel memristive and neuromorphic functions.


title: Plasma-fixated nitrogen as fertilizer for turf grass

abstract: We investigated the use of plasma-fixated nitrogen, which produces nitrates (NO3-) in water, as a possible nitrogen fertilizer for recreational turf such as ryegrass and bentgrass. Experiments were carried out to study the effects of nitrate concentration on growth, the further effects of adding phosphorous (P) and potassium (K) to the plasma nitrated solution to make an N-P-K complete fertilizer, and to compare the efficacy of plasma-fixated nitrogen to sodium nitrate (NaNO3) and potassium nitrate (KNO3). The results indicate that the growth and biomass of the plants were strongly dependent on the concentration of the plasma-fixated nitrogen. Adding P-K to the plasma-fixated nitrogen improved grass growth. Grass that was supplied plasma-fixated nitrogen had improved growth compared to those supplied with equal amounts of NaNO3and KNO3. This work highlights the potential use of plasma-fixated nitrogen as a fertilizer source for commonly used turf grass.


title: Searching for More Efficient Dynamic Programs

abstract: Computational models of human language often involve combinatorial problems. For instance, a probabilistic parser may marginalize over exponentially many trees to make predictions. Algorithms for such problems often employ dynamic programming and are not always unique. Finding one with optimal asymptotic runtime can be unintuitive, time-consuming, and error-prone. Our work aims to automate this laborious process. Given an initial correct declarative program, we search for a sequence of semantics-preserving transformations to improve its running time as much as possible. To this end, we describe a set of program transformations, a simple metric for assessing the efficiency of a transformed program, and a heuristic search procedure to improve this metric. We show that in practice, automated search -- like the mental search performed by human programmers -- can find substantial improvements to the initial program. Empirically, we show that many common speed-ups described in the NLP literature could have been discovered automatically by our system.


title: Grounding-aware RRT* for Path Planning and Safe Navigation of Marine Crafts in Confined Waters

abstract: The paper presents a path planning algorithm based on RRT* that addresses the risk of grounding during evasive manoeuvres to avoid collision. The planner achieves this objective by integrating a collective navigation experience with the systematic use of water depth information from the electronic navigational chart. Multivariate kernel density estimation is applied to historical AIS data to generate a probabilistic model describing seafarer's best practices while sailing in confined waters. This knowledge is then encoded into the RRT* cost function to penalize path deviations that would lead own ship to sail in shallow waters. Depth contours satisfying the own ship draught define the actual navigable area, and triangulation of this non-convex region is adopted to enable uniform sampling. This ensures the optimal path deviation.


title: Differential equations for the recurrence coefficients of semi-classical orthogonal polynomials and their relation to the Painlev\'e equations via the geometric approach

abstract: In this paper we present a general scheme for how to relate differential equations for the recurrence coefficients of semi-classical orthogonal polynomials to the Painlev\'e equations using the geometric framework of the Okamoto Space of Initial Conditions. We demonstrate this procedure in two examples. For semi-classical Laguerre polynomials appearing in \cite{HC17}, we show how the recurrence coefficients are connected to the fourth Painlev\'e equation. For discrete orthogonal polynomials associated with the hypergeometric weight appearing in \cite{FVA18} we discuss the relation of the recurrence coefficients to the sixth Painlev\'e equation, extending the results of \cite{DFS19}, where a similar approach was used for a discrete system for the same recurrence coefficients. Though the discrete and differential systems here share the same geometry, the construction of the space of initial conditions from the differential system is different and reveals extra considerations that must be made. We also discuss a number of related topics in the context of the geometric approach, such as Hamiltonian forms of the differential equations for the recurrence coefficients, Riccati solutions for special parameter values, and associated discrete Painlev\'e equations.


title: Multi-modal Wound Classification using Wound Image and Location by Deep Neural Network

abstract: Wound classification is an essential step of wound diagnosis. An efficient classifier can assist wound specialists in classifying wound types with less financial and time costs and help them decide an optimal treatment procedure. This study developed a deep neural network-based multi-modal classifier using wound images and their corresponding locations to categorize wound images into multiple classes, including diabetic, pressure, surgical, and venous ulcers. A body map is also developed to prepare the location data, which can help wound specialists tag wound locations more efficiently. Three datasets containing images and their corresponding location information are designed with the help of wound specialists. The multi-modal network is developed by concatenating the image-based and location-based classifier's outputs with some other modifications. The maximum accuracy on mixed-class classifications (containing background and normal skin) varies from 77.33% to 100% on different experiments. The maximum accuracy on wound-class classifications (containing only diabetic, pressure, surgical, and venous) varies from 72.95% to 98.08% on different experiments. The proposed multi-modal network also shows a significant improvement in results from the previous works of literature.


title: Optical Observations and Modeling of the Superluminous Supernova 2018lfe

abstract: We present optical imaging and spectroscopy of SN\,2018lfe, which we classify as a Type I superluminous supernova (SLSN-I) at a redshift of $z = 0.3501$ with a peak absolute magnitude of $M_r\approx -22.1$ mag, one of the brightest SLSNe discovered. SN\,2018lfe was identified for follow-up using our FLEET machine learning pipeline. Both the light curve and the spectra of SN\,2018lfe are consistent with the broad population of SLSNe. We fit the light curve with a magnetar central engine model and find an ejecta mass of $M_{\rm ej}\approx 3.8$ M$_\odot$, a magnetar spin period of $P\approx 2.9$ ms and a magnetic field strength of $B_{\perp}\approx 2.8\times 10^{14}$ G. The magnetic field strength is near the top of the distribution for SLSNe, while the spin period and ejecta mass are near the median values of the distribution for SLSNe. From late-time imaging and spectroscopy we find that the host galaxy of SN\,2018lfe has an absolute magnitude of $M_r\approx -17.85$, ($L_B \approx 0.029$ $L^*$), and an inferred metallicity of $Z\approx 0.3$ Z$_\odot$, and star formation rate of $\approx 0.8$ M$_\odot$ yr$^{-1}$.


title: Contactless Series Resistance Imaging of Perovskite Solar Cells via Inhomogeneous Illumination

abstract: A contactless effective series resistance imaging method for large area perovskite solar cells that is based on photoluminescence imaging with non-uniform illumination is introduced and demonstrated experimentally. The proposed technique is applicable to partially and fully processed perovskite solar cells if laterally conductive layers are present. The capability of the proposed contactless method to detect features with high effective series resistance is validated by comparison with various contacted mode luminescence imaging techniques. The method can reliably provide information regarding the severeness of the detected series resistance through photo-excitation pattern manipulation. Application of the method to sub-cells in monolithic tandem devices, without the need for electrical contacting the terminals, appears feasible.


title: Combining GEDI and Sentinel-2 for wall-to-wall mapping of tall and short crops

abstract: High resolution crop type maps are an important tool for improving food security, and remote sensing is increasingly used to create such maps in regions that possess ground truth labels for model training. However, these labels are absent in many regions, and models trained in other regions on typical satellite features, such as those from optical sensors, often exhibit low performance when transferred. Here we explore the use of NASA's Global Ecosystem Dynamics Investigation (GEDI) spaceborne lidar instrument, combined with Sentinel-2 optical data, for crop type mapping. Using data from three major cropped regions (in China, France, and the United States) we first demonstrate that GEDI energy profiles are capable of reliably distinguishing maize, a crop typically above 2m in height, from crops like rice and soybean that are shorter. We further show that these GEDI profiles provide much more invariant features across geographies compared to spectral and phenological features detected by passive optical sensors. GEDI is able to distinguish maize from other crops within each region with accuracies higher than 84%, and able to transfer across regions with accuracies higher than 82% compared to 64% for transfer of optical features. Finally, we show that GEDI profiles can be used to generate training labels for models based on optical imagery from Sentinel-2, thereby enabling the creation of 10m wall-to-wall maps of tall versus short crops in label-scarce regions. As maize is the second most widely grown crop in the world and often the only tall crop grown within a landscape, we conclude that GEDI offers great promise for improving global crop type maps.


title: Simulation of the Kinetic Energy Spectrum of Primary Electron Generated by a 6 keV X-Ray Source in a Gas Volume

abstract: With the intention of creating a complete simulation of the Thick GEM detector, this work presents a partial simulation of the interaction of 6 keV X-Ray radiation with a gas volume. For the gas, it was used pure Argon, Xenon and CO2 as well as a mixture of each noble gas with CO2 in a 70\%-30\% proportion. For this range of energy, the dominant physical effect is the photoelectric, therefore knowing the binding energies of the atoms we predict the kinetic energy of the primary electron. A simulated energy spectrum of the primary electron kinetic energy was obtained, which contains peaks corresponding to the inner shell's ionization of the atoms. The spectrum shows very little Compton effect, transferring only small energies (in the order of eVs). We concluded that energy spectrum for the mixture is a linear combination of the spectrum of its individual components.


title: Algorithmic Auditing and Social Justice: Lessons from the History of Audit Studies

abstract: Algorithmic audits have been embraced as tools to investigate the functioning and consequences of sociotechnical systems. Though the term is used somewhat loosely in the algorithmic context and encompasses a variety of methods, it maintains a close connection to audit studies in the social sciences--which have, for decades, used experimental methods to measure the prevalence of discrimination across domains like housing and employment. In the social sciences, audit studies originated in a strong tradition of social justice and participatory action, often involving collaboration between researchers and communities; but scholars have argued that, over time, social science audits have become somewhat distanced from these original goals and priorities. We draw from this history in order to highlight difficult tensions that have shaped the development of social science audits, and to assess their implications in the context of algorithmic auditing. In doing so, we put forth considerations to assist in the development of robust and engaged assessments of sociotechnical systems that draw from auditing's roots in racial equity and social justice.


title: A Speckle-Imaging Search for Close and Very Faint Companions to the Nearest and Brightest Wolf-Rayet Stars

abstract: Gravitationally bound companions to stars enable determinations of their masses, and offer clues to their formation, evolution and dynamical histories. So motivated, we have carried out a speckle imaging survey of eight of the nearest and brightest Wolf-Rayet (WR) stars to directly measure the frequency of their resolvable companions, and to search for much fainter companions than hitherto possible. We found one new, close companion to each of WR 113, WR 115 and WR 120 in the separation range 0.2" - 1.2". Our results provide more evidence that similar-brightness, close companions to WR stars are common. More remarkably, they also demonstrate that the predicted, but much fainter and thus elusive companions to WR stars are now within reach of modern speckle cameras on 8m class telescopes by finding the first example. The new companion to WR 113 is just 1.16" distant from it, and is 8 magnitudes fainter than the WR star. The empirical probability of a chance line-of-sight of the faint companion at the position of WR 113 is < 0.5%, though we cannot yet prove or disprove if the two stars are gravitationally bound. If these three new detections are physical companions we suggest, based on their narrowband magnitudes, colors, reddenings and GAIA distances that the companions to WR113, WR 115 and WR 120 are an F-type dwarf, an early B-type dwarf, and a WNE-type WR star, respectively.


title: GRiD: GPU-Accelerated Rigid Body Dynamics with Analytical Gradients

abstract: We introduce GRiD: a GPU-accelerated library for computing rigid body dynamics with analytical gradients. GRiD was designed to accelerate the nonlinear trajectory optimization subproblem used in state-of-the-art robotic planning, control, and machine learning, which requires tens to hundreds of naturally parallel computations of rigid body dynamics and their gradients at each iteration. GRiD leverages URDF parsing and code generation to deliver optimized dynamics kernels that not only expose GPU-friendly computational patterns, but also take advantage of both fine-grained parallelism within each computation and coarse-grained parallelism between computations. Through this approach, when performing multiple computations of rigid body dynamics algorithms, GRiD provides as much as a 7.2x speedup over a state-of-the-art, multi-threaded CPU implementation, and maintains as much as a 2.5x speedup when accounting for I/O overhead. We release GRiD as an open-source library for use by the wider robotics community.


title: Design Guidelines for Prompt Engineering Text-to-Image Generative Models

abstract: Text-to-image generative models are a new and powerful way to generate visual artwork. However, the open-ended nature of text as interaction is double-edged; while users can input anything and have access to an infinite range of generations, they also must engage in brute-force trial and error with the text prompt when the result quality is poor. We conduct a study exploring what prompt keywords and model hyperparameters can help produce coherent outputs. In particular, we study prompts structured to include subject and style keywords and investigate success and failure modes of these prompts. Our evaluation of 5493 generations over the course of five experiments spans 51 abstract and concrete subjects as well as 51 abstract and figurative styles. From this evaluation, we present design guidelines that can help people produce better outcomes from text-to-image generative models.


title: Event-Triggered Distributed Stabilization of Interconnected Multiagent Systems with Abnormal Agent and Control Layers: Theoretical Analysis

abstract: A graph theoretic framework recently has been proposed to stabilize interconnected multiagent systems in a distributed fashion, while systematically capturing the architectural aspect of cyber-physical systems with separate agent or physical layer and control or cyber layer. Based on that development, in addition to the modeling uncertainties over the agent layer, we consider a scenario where the control layer is subject to the denial of service attacks. We propose a step-by-step procedure to design a control layer that, in the presence of the aforementioned abnormalities, guarantees a level of robustness and resiliency for the final two-layer interconnected multiagent system. The incorporation of an event-triggered strategy further ensures an effective use of the limited energy and communication resources over the control layer. We theoretically prove the resilient, robust, and Zeno-free convergence of all state trajectories to the origin and, via a simulation study, discuss the feasibility of the proposed ideas.


title: CORNET 2.0: A Co-Simulation Middleware for Robot Networks

abstract: We present a networked co-simulation framework for multi-robot systems applications. We require a simulation framework that captures both physical interactions and communications aspects to effectively design such complex systems. This is necessary to co-design the multi-robots' autonomy logic and the communication protocols. The proposed framework extends existing tools to simulate the robot's autonomy and network-related aspects. We have used Gazebo with ROS/ROS2 to develop the autonomy logic for robots and mininet-WiFi as the network simulator to capture the cyber-physical systems properties of the multi-robot system. This framework addresses the need to seamlessly integrate the two simulation environments by synchronizing mobility and time, allowing for easy migration of the algorithms to real platforms. The framework supports container-based virtualization and extends a generic robotic framework by decoupling the data plane and control plane.


title: Explainable Identification of Dementia from Transcripts using Transformer Networks

abstract: Alzheimer's disease (AD) is the main cause of dementia which is accompanied by loss of memory and may lead to severe consequences in peoples' everyday life if not diagnosed on time. Very few works have exploited transformer-based networks and despite the high accuracy achieved, little work has been done in terms of model interpretability. In addition, although Mini-Mental State Exam (MMSE) scores are inextricably linked with the identification of dementia, research works face the task of dementia identification and the task of the prediction of MMSE scores as two separate tasks. In order to address these limitations, we employ several transformer-based models, with BERT achieving the highest accuracy accounting for 87.50%. Concurrently, we propose an interpretable method to detect AD patients based on siamese networks reaching accuracy up to 83.75%. Next, we introduce two multi-task learning models, where the main task refers to the identification of dementia (binary classification), while the auxiliary one corresponds to the identification of the severity of dementia (multiclass classification). Our model obtains accuracy equal to 86.25% on the detection of AD patients in the multi-task learning setting. Finally, we present some new methods to identify the linguistic patterns used by AD patients and non-AD ones, including text statistics, vocabulary uniqueness, word usage, correlations via a detailed linguistic analysis, and explainability techniques (LIME). Findings indicate significant differences in language between AD and non-AD patients.


title: Convex-cocompact subgroups of the Goeritz group

abstract: We show that finitely-generated, purely pseudo-Anosov subgroups of the genus-2 Goeritz group are convex cocompact in the genus-2 mapping class group.


title: Generalized Multivariable Grid-Forming Control Design for Power Converters

abstract: The grid-forming converter is an important unit in the future power system with more inverter-interfaced generators. However, improving its performance is still a key challenge. This paper proposes a generalized architecture of the grid-forming converter from the view of multivariable feedback control. As a result, many of the existing popular control strategies, i.e., droop control, power synchronization control, virtual synchronous generator control, matching control, dispatchable virtual oscillator control, and their improved forms are unified into a multivariable feedback control transfer matrix working on several linear and nonlinear error signals. Meanwhile, unlike the traditional assumptions of decoupling between AC and DC control, active power and reactive power control, the proposed configuration simultaneously takes all of them into consideration, which therefore can provide better performance. As an example, a new multi-input-multi-output-based grid-forming (MIMO-GFM) control is proposed based on the generalized configuration. To cope with the multivariable feedback, an optimal and structured $H_{\infty}$ synthesis is used to design the control parameters. At last, simulation and experimental results show superior performance and robustness of the proposed configuration and control.


title: Dirac lines and loop at the Fermi level in the Time-Reversal Symmetry Breaking Superconductor LaNiGa$_2$

abstract: Unconventional superconductors have Cooper pairs with lower symmetries than in conventional superconductors. In most unconventional superconductors, the additional symmetry breaking occurs in relation to typical ingredients such as strongly correlated Fermi liquid phases, magnetic fluctuations, or strong spin-orbit coupling in noncentrosymmetric structures. In this article, we show that the time-reversal symmetry breaking in the superconductor LaNiGa$_2$ is enabled by its previously unknown topological electronic band structure. Our single crystal diffraction experiments indicate a nonsymmorphic crystal structure, in contrast to the previously reported symmorphic structure. The nonsymmorphic symmetries transform the $k_z=\pi/c$ plane of the Brillouin zone boundary into a node-surface. Band-structure calculations reveal that distinct Fermi surfaces become degenerate on the node-surface and form Dirac lines and a Dirac loop at the Fermi level. Two symmetry related Dirac points remain degenerate under spin-orbit coupling. ARPES measurements confirm the calculations and provide evidence for the Fermi surface degeneracies on the node-surface. These unique topological features enable an unconventional superconducting gap in which time-reversal symmetry can be broken in the absence of other typical ingredients. LaNiGa$_2$ is therefore a topological crystalline superconductor that breaks time-reversal symmetry without any overlapping magnetic ordering or fluctuations. Our findings will enable future discoveries of additional topological superconductors.


title: Monotonicity in the averaging process

abstract: We investigate an averaging process that describes how interacting agents approach consensus through binary interactions. In each elementary step, two agents are selected at random and they reach compromise by adopting their opinion average. We show that the fraction of agents with a monotonically decreasing opinion decays as $e^{-\alpha t}$, and that the exponent $\alpha=\tfrac{1}{2}-\tfrac{1+\ln \ln 2}{4\ln 2}$ is selected as the extremum from a continuous spectrum of possible values. The opinion distribution of monotonic agents is asymmetric, and it becomes self-similar at large times. Furthermore, the tails of the opinion distribution are algebraic, and they are characterized by two distinct and nontrivial exponents. We also explore statistical properties of agents with an opinion strictly above average.


title: Compact Quantum Metric Spaces from Free Graph Algebras

abstract: Starting with a vertex-weighted pointed graph $(\Gamma,\mu,v_0)$, we form the free loop algebra $\mathcal{S}_0$ defined in Hartglass-Penneys' article on canonical $\rm C^*$-algebras associated to a planar algebra. Under mild conditions, $\mathcal{S}_0$ is a non-nuclear simple $\rm C^*$-algebra with unique tracial state. There is a canonical polynomial subalgebra $A\subset \mathcal{S}_0$ together with a Dirac number operator $N$ such that $(A, L^2A,N)$ is a spectral triple. We prove the Haagerup-type bound of Ozawa-Rieffel to verify $(\mathcal{S}_0, A, N)$ yields a compact quantum metric space in the sense of Rieffel. We give a weighted analog of Benjamini-Schramm convergence for vertex-weighted pointed graphs. As our $\rm C^*$-algebras are non-nuclear, we adjust the Lip-norm coming from $N$ to utilize the finite dimensional filtration of $A$. We then prove that convergence of vertex-weighted pointed graphs leads to quantum Gromov-Hausdorff convergence of the associated adjusted compact quantum metric spaces. As an application, we apply our construction to the Guionnet-Jones-Shyakhtenko (GJS) $\rm C^*$-algebra associated to a planar algebra. We conclude that the compact quantum metric spaces coming from the GJS $\rm C^*$-algebras of many infinite families of planar algebras converge in quantum Gromov-Hausdorff distance.


title: Finite rigid sets and the separating curve complex

abstract: For each closed surface of genus $g\ge3$, we find a finite subcomplex of the separating curve complex that is rigid with respect to incidence-preserving maps.


title: NOPE: A Corpus of Naturally-Occurring Presuppositions in English

abstract: Understanding language requires grasping not only the overtly stated content, but also making inferences about things that were left unsaid. These inferences include presuppositions, a phenomenon by which a listener learns about new information through reasoning about what a speaker takes as given. Presuppositions require complex understanding of the lexical and syntactic properties that trigger them as well as the broader conversational context. In this work, we introduce the Naturally-Occurring Presuppositions in English (NOPE) Corpus to investigate the context-sensitivity of 10 different types of presupposition triggers and to evaluate machine learning models' ability to predict human inferences. We find that most of the triggers we investigate exhibit moderate variability. We further find that transformer-based models draw correct inferences in simple cases involving presuppositions, but they fail to capture the minority of exceptional cases in which human judgments reveal complex interactions between context and triggers.


title: Investigating the relation between gravitational wave tests of general relativity

abstract: Gravitational wave observations of compact binary coalescences provide precision probes of strong-field gravity. There is thus now a standard set of null tests of general relativity (GR) applied to LIGO-Virgo detections and many more such tests proposed. However, the relation between all these tests is not yet well understood. We start to investigate this by applying a set of standard tests to simulated observations of binary black holes in GR and with phenomenological deviations from GR. We consider four types of tests: residuals, inspiral-merger-ringdown consistency, parameterized phasing (two varieties), and modified dispersion relation. We also check the consistency of the unmodeled reconstruction of the waveforms with the waveform recovered using GR templates. These tests are applied to simulated observations similar to GW150914 with both large and small deviations from GR and similar to GW170608 just with small deviations from GR. We find that while very large deviations from GR are picked up with high significance by almost all tests, more moderate deviations are picked up by only a few tests, and some deviations are not recognized as GR violations by any test at the moderate signal-to-noise ratios we consider. Moreover, the tests that identify various deviations with high significance are not necessarily the expected ones. In particular, the parameterized tests recover PN test parameters much closer to zero than their true values in some cases. Additionally, we find that of the GR deviations we consider, the residuals test is only able to detect extreme deviations from GR that no longer look like binary black hole coalescences in GR. The reconstruction comparison shows more promise for detecting relatively small GR deviations in an unmodeled framework, at least for high-mass systems.


title: A short proof of Kneser's theorem via transversals

abstract: A very short proof of Kneser's theorem via transversal is given.


title: Personalization, Privacy, and Me

abstract: News recommendation and personalization is not a solved problem. People are growing concerned of their data being collected in excess in the name of personalization and the usage of it for purposes other than the ones they would think reasonable. Our experience in building personalization products for publishers while adhering to safeguard user privacy led us to investigate more on the user perspective of privacy and personalization. We conducted a survey to explore people's experience with personalization and privacy and the viewpoints of different age groups. In this paper, we share our major findings with publishers and the community that can inform algorithmic design and implementation of the next generation of news recommender systems, which must put the human at its core and reach a balance between personalization experiences and privacy to reap the benefits of both.


title: Self-Similarity of $k$-Nearest Neighbor Distributions in Scale-Free Simulations

abstract: We use the $k$-nearest neighbor probability distribution function ($k$NN-PDF, Banerjee & Abel 2021) to assess convergence in a scale-free $N$-body simulation. Compared to our previous two-point analysis, the $k$NN-PDF allows us to quantify our results in the language of halos and numbers of particles, while also incorporating non-Gaussian information. We find good convergence for 32 particles and greater at densities typical of halos, while 16 particles and fewer appears unconverged. Halving the softening length extends convergence to higher densities, but not to fewer particles. Our analysis is less sensitive to voids, but we analyze a limited range of underdensities and find evidence for convergence at 16 particles and greater even in sparse voids.


title: ML-aided power allocation for Tactical MIMO

abstract: We study the problem of optimal power allocation in single-hop multi-antenna ad-hoc wireless networks. A standard technique to solve this problem involves optimizing a tri-convex function under power constraints using a block-coordinate-descent based iterative algorithm. This approach, termed WMMSE, tends to be computationally complex and time consuming. Several learning-based approaches have been proposed to speed up the power allocation process. A recent work, UWMMSE, learns an affine transformation of a WMMSE parameter in an unfolded structure to accelerate convergence. In spite of achieving promising results, its application is limited to single-antenna wireless networks. In this work, we present a UWMMSE framework for power allocation in (multiple-input multiple-output) MIMO interference networks. A major advantage of this method lies in its use of low-complexity learnable systems in which the number of parameters scales linearly with respect to the hidden layer size of embedded neural architectures and the product of the number of transmitter and receiver antennas only, fully independent of the number of transceivers in the network. We illustrate the superiority of our method through an empirical study of our approach in comparison to WMMSE and also analyze its robustness to changes in channel conditions and network size.


title: On non-normal subgroup perfect codes

abstract: Let $X = (V,E)$ be a graph. A subset $C \subseteq V(X)$ is a \emph{perfect code} of $X$ if $C$ is a coclique of $X$ with the property that any vertex in $V(X)\setminus C$ is adjacent to exactly one vertex in $C$. Given a finite group $G$ with identity element $e$ and $H\leq G$, $H$ is a \emph{subgroup perfect code} of $G$ if there exists an inverse-closed subset $S \subseteq G\setminus \{e\}$ such that $H$ is a perfect code of the Cayley graph $\operatorname{Cay}(G,S)$ of $G$ with connection set $S$. In this short note, we give an infinite family of finite groups $G$ admitting a non-normal subgroup perfect code $H$ such that there exists $ g\in G$ with $g^2\in H$ but $(gh)^2 \neq e$, for all $h \in H$; thus, answering a question raised by Wang, Xia, and Zhou in [Perfect sets in Cayley graphs. {\it arXiv preprint} arXiv:2006.05100, 2020].


title: Remarks on the preservation and breaking of translational symmetry for a class of ODEs

abstract: In this paper, we provide both a preservation and breaking of symmetry theorem for $2\pi$-periodic problems of the form \begin{align*} \begin{cases} -u''(t) + g(u(t)) = f(t)\cr u(0) - u(2\pi) = u'(0) - u'(2\pi) = 0 \end{cases} \end{align*} where $g: \mathbb{R} \to \mathbb{R}$ is a given $C^1$ function and $f: [0,2\pi] \to \mathbb{R}$ is continuous. We provide a preservation of symmetry result that is analogous to one given by Willem (Willem, 1989) and a generalization of the theorem given by Costa-Fang (Costa and Fang, 2019). Both of these theorems use group actions that are not normally considered in the literature.


title: Strategic Inventories in a Supply Chain with Downstream Cournot Duopoly

abstract: The inventories carried in a supply chain as a strategic tool to influence the competing firms are considered to be strategic inventories (SI). We present a two-period game-theoretic supply chain model, in which a singular manufacturer supplies products to a pair of identical Cournot duopolistic retailers. We show that the SI carried by the retailers under dynamic contract is Pareto-dominating for the manufacturer, retailers, consumers, the channel, and the society as well. We also find that the retailer's SI, however, can be eliminated when the manufacturer commits wholesale contract or inventory holding cost is too high. In comparing the cases with and without downstream competition, we also show that the downstream Cournot duopoly undermines the retailers in profits, but benefits all others.


title: Scalable Average Consensus with Compressed Communications

abstract: We propose a new decentralized average consensus algorithm with compressed communication that scales linearly with the network size n. We prove that the proposed method converges to the average of the initial values held locally by the agents of a network when agents are allowed to communicate with compressed messages. The proposed algorithm works for a broad class of compression operators (possibly biased), where agents interact over arbitrary static, undirected, and connected networks. We further present numerical experiments that confirm our theoretical results and illustrate the scalability and communication efficiency of our algorithm.


title: Phase Stability and Sodium-Vacancy Orderings in a NaSICON Electrode

abstract: We elucidate the thermodynamics of sodium (Na) intercalation into the sodium super-ionic conductor (NaSICON)-type electrode, Na$_x$V$_2$(PO$_4$)$_3$, for promising Na-ion batteries with high-power density. This is the first report of a computational temperature-composition phase diagram of the NaSICON-type electrode Na$_x$V$_2$(PO$_4$)$_3$. We identify two thermodynamically stable phases at the compositions Na$_2$V$_2$(PO$_4$)$_3$ and Na$_{3.5}$V$_2$(PO$_4$)$_3$, and their structural features are described for the first time based on our computational analysis. We unveil the crystal-structure and the electronic-structure origins of the ground-state compositions associated with specific Na/vacancy arrangements, which are driven by charge orderings on the vanadium sites. These results are significant for the optimization of high-energy and power densities electrodes for sustainable Na-ion batteries


title: $\mathcal{L}_1$ Adaptive Augmentation for Geometric Tracking Control of Quadrotors

abstract: This paper introduces an $\mathcal{L}_1$ adaptive control augmentation for geometric tracking control of quadrotors. In the proposed design, the $\mathcal{L}_1$ augmentation handles nonlinear (time- and state-dependent) uncertainties in the quadrotor dynamics without assuming or enforcing parametric structures, while the baseline geometric controller achieves stabilization of the known nonlinear model of the system dynamics. The $\mathcal{L}_1$ augmentation applies to both the rotational and the translational dynamics. Experimental results demonstrate that the augmented geometric controller shows consistent and (on average five times) smaller trajectory tracking errors compared with the geometric controller alone when tested for different trajectories and under various types of uncertainties/disturbances.


title: Behavior of k-NN as an Instance-Based Explanation Method

abstract: Adoption of DL models in critical areas has led to an escalating demand for sound explanation methods. Instance-based explanation methods are a popular type that return selective instances from the training set to explain the predictions for a test sample. One way to connect these explanations with prediction is to ask the following counterfactual question - how does the loss and prediction for a test sample change when explanations are removed from the training set? Our paper answers this question for k-NNs which are natural contenders for an instance-based explanation method. We first demonstrate empirically that the representation space induced by last layer of a neural network is the best to perform k-NN in. Using this layer, we conduct our experiments and compare them to influence functions (IFs) ~\cite{koh2017understanding} which try to answer a similar question. Our evaluations do indicate change in loss and predictions when explanations are removed but we do not find a trend between $k$ and loss or prediction change. We find significant stability in the predictions and loss of MNIST vs. CIFAR-10. Surprisingly, we do not observe much difference in the behavior of k-NNs vs. IFs on this question. We attribute this to training set subsampling for IFs.


title: Koopman Linearization for Data-Driven Batch State Estimation of Control-Affine Systems

abstract: We present the Koopman State Estimator (KoopSE), a framework for model-free batch state estimation of control-affine systems that makes no linearization assumptions, requires no problem-specific feature selections, and has an inference computational cost that is independent of the number of training points. We lift the original nonlinear system into a higher-dimensional Reproducing Kernel Hilbert Space (RKHS), where the system becomes bilinear. The time-invariant model matrices can be learned by solving a least-squares problem on training trajectories. At test time, the system is algebraically manipulated into a linear time-varying system, where standard batch linear state estimation techniques can be used to efficiently compute state means and covariances. Random Fourier Features (RFF) are used to combine the computational efficiency of Koopman-based methods and the generality of kernel-embedding methods. KoopSE is validated experimentally on a localization task involving a mobile robot equipped with ultra-wideband receivers and wheel odometry. KoopSE estimates are more accurate and consistent than the standard model-based extended Rauch-Tung-Striebel (RTS) smoother, despite KoopSE having no prior knowledge of the system's motion or measurement models.


title: ZFlow: Gated Appearance Flow-based Virtual Try-on with 3D Priors

abstract: Image-based virtual try-on involves synthesizing perceptually convincing images of a model wearing a particular garment and has garnered significant research interest due to its immense practical applicability. Recent methods involve a two stage process: i) warping of the garment to align with the model ii) texture fusion of the warped garment and target model to generate the try-on output. Issues arise due to the non-rigid nature of garments and the lack of geometric information about the model or the garment. It often results in improper rendering of granular details. We propose ZFlow, an end-to-end framework, which seeks to alleviate these concerns regarding geometric and textural integrity (such as pose, depth-ordering, skin and neckline reproduction) through a combination of gated aggregation of hierarchical flow estimates termed Gated Appearance Flow, and dense structural priors at various stage of the network. ZFlow achieves state-of-the-art results as observed qualitatively, and on quantitative benchmarks of image quality (PSNR, SSIM, and FID). The paper presents extensive comparisons with other existing solutions including a detailed user study and ablation studies to gauge the effect of each of our contributions on multiple datasets.


title: Experimental Demonstration of a New Oscillator Concept Based on Degenerate Band Edge in Microstrip Circuit

abstract: The first experimental demonstration of an oscillator based on a periodic, resonant microstrip circuit with a degenerate band edge (DBE) is presented. The DBE is a fourth-order exceptional degeneracy of the eigenmodes in a lossless periodic waveguide that is visible in the wavenumber-frequency dispersion diagram, and the periodic microstrip behaves as a frequency selective medium. The presence of the DBE condition and the associated DBE resonance allow for a stable, single-mode oscillation as well as stability with varying the load resistance.


title: Search for and Identification of Young Compact Galactic Supernova Remnants Using THOR

abstract: Young Supernova remnants (SNRs) with smaller angular sizes are likely missing from existing radio SNR catalogues, caused by observational constraints and selection effects. In order to find new compact radio SNR candidates, we searched the high angular resolution (25") THOR radio survey of the first quadrant of the galaxy. We selected sources with non-thermal radio spectra. HI absorption spectra and channel maps were used to identify which sources are galactic and to estimate their distances. Two new compact SNRs were found: G31.299$-$0.493 and G18.760$-$0.072, of which the latter was a previously suggested SNR candidate. The distances to these SNRs are 5.0 $\pm$ 0.3 kpc and 4.7 $\pm$ 0.2 kpc, respectively. Based on the SN rate in the galaxy or on the statistics of known SNRs, we estimate that there are 15$-$20 not yet detected compact SNRs in the galaxy and that the THOR survey area should contain three or four. Our detection of two SNRs (half the expected number) is consistent with the THOR sensitivity limit compared with the distribution of integrated flux densities of SNRs.


title: Vibrational density of states capture the role of dynamic allostery in protein evolution

abstract: Previous studies of the flexibilities of ancestral proteins suggests that proteins evolve their function by altering their native state ensemble. Here we propose a more direct method of visualizing this by measuring the changes in the vibrational density of states (VDOS) of proteins as they evolve. Through analysis of VDOS profiles of ancestral and extant proteins we observe that $\beta$-lactamase and thioredoxins evolve by altering their density of states in the terahertz region. Particularly, the shift in VDOS profiles between ancestral and extant proteins suggests that nature utilize dynamic allostery for functional evolution. Moreover, we also show that VDOS profile of individual position can be used to describe the flexibility changes, particularly those without any amino acid substitution.


title: WaveCorr: Correlation-savvy Deep Reinforcement Learning for Portfolio Management

abstract: The problem of portfolio management represents an important and challenging class of dynamic decision making problems, where rebalancing decisions need to be made over time with the consideration of many factors such as investors preferences, trading environments, and market conditions. In this paper, we present a new portfolio policy network architecture for deep reinforcement learning (DRL)that can exploit more effectively cross-asset dependency information and achieve better performance than state-of-the-art architectures. In particular, we introduce a new property, referred to as \textit{asset permutation invariance}, for portfolio policy networks that exploit multi-asset time series data, and design the first portfolio policy network, named WaveCorr, that preserves this invariance property when treating asset correlation information. At the core of our design is an innovative permutation invariant correlation processing layer. An extensive set of experiments are conducted using data from both Canadian (TSX) and American stock markets (S&P 500), and WaveCorr consistently outperforms other architectures with an impressive 3%-25% absolute improvement in terms of average annual return, and up to more than 200% relative improvement in average Sharpe ratio. We also measured an improvement of a factor of up to 5 in the stability of performance under random choices of initial asset ordering and weights. The stability of the network has been found as particularly valuable by our industrial partner.


title: A Three Step Training Approach with Data Augmentation for Morphological Inflection

abstract: We present the BME submission for the SIGMORPHON 2021 Task 0 Part 1, Generalization Across Typologically Diverse Languages shared task. We use an LSTM encoder-decoder model with three step training that is first trained on all languages, then fine-tuned on each language families and finally finetuned on individual languages. We use a different type of data augmentation technique in the first two steps. Our system outperformed the only other submission. Although it remains worse than the Transformer baseline released by the organizers, our model is simpler and our data augmentation techniques are easily applicable to new languages. We perform ablation studies and show that the augmentation techniques and the three training steps often help but sometimes have a negative effect.


title: The Empirical Content of Bayesianism

abstract: This paper characterizes the conditions under which the observed beliefs of a group of agents are consistent with Bayesian updating. Beliefs are consistent with Bayesianism if they arise from the application of Bayes' rule given some subjective distribution for the state and the signals agents observe between periods. The paper's main finding is that beliefs are consistent with Bayesianism if and only if the mean of the distribution of posteriors is uniformly absolutely continuous with respect to the prior. Furthermore, the paper shows that the existing results on the empirical content of Bayesianism rely on additional restrictions on permissible subjective distributions, such as the requirement that agents have correct beliefs about the distribution of signals.


title: HeMI: Multi-view Embedding in Heterogeneous Graphs

abstract: Many real-world graphs involve different types of nodes and relations between nodes, being heterogeneous by nature. The representation learning of heterogeneous graphs (HGs) embeds the rich structure and semantics of such graphs into a low-dimensional space and facilitates various data mining tasks, such as node classification, node clustering, and link prediction. In this paper, we propose a self-supervised method that learns HG representations by relying on knowledge exchange and discovery among different HG structural semantics (meta-paths). Specifically, by maximizing the mutual information of meta-path representations, we promote meta-path information fusion and consensus, and ensure that globally shared semantics are encoded. By extensive experiments on node classification, node clustering, and link prediction tasks, we show that the proposed self-supervision both outperforms and improves competing methods by 1% and up to 10% for all tasks.


title: Will this Question be Answered? Question Filtering via Answer Model Distillation for Efficient Question Answering

abstract: In this paper we propose a novel approach towards improving the efficiency of Question Answering (QA) systems by filtering out questions that will not be answered by them. This is based on an interesting new finding: the answer confidence scores of state-of-the-art QA systems can be approximated well by models solely using the input question text. This enables preemptive filtering of questions that are not answered by the system due to their answer confidence scores being lower than the system threshold. Specifically, we learn Transformer-based question models by distilling Transformer-based answering models. Our experiments on three popular QA datasets and one industrial QA benchmark demonstrate the ability of our question models to approximate the Precision/Recall curves of the target QA system well. These question models, when used as filters, can effectively trade off lower computation cost of QA systems for lower Recall, e.g., reducing computation by ~60%, while only losing ~3-4% of Recall.


title: Exact eigenstates of extended SU($N$) Hubbard models: Generalizations of $\eta$-pairing states with $N$-particle off-diagonal long-range order

abstract: We consider $N$-particle generalizations of $\eta$-pairing states in a chain of $N$-component fermions and show that these states are exact (high-energy) eigenstates of an extended SU($N$) Hubbard model. We compute the singlet correlation function of the states and find that its behavior is qualitatively different for even and odd $N$. When $N$ is even, these states exhibit off-diagonal long-range order in $N$-particle reduced density matrix. On the other hand, when $N$ is odd, the correlations decay exponentially with distance in the bulk, but end-to-end correlations do not vanish in the thermodynamic limit. Finally, we prove that these states are the unique ground states of suitably tailored Hamiltonians.

categories: cond-mat.str-el cond-mat.quant-gas cond-mat.stat-mech cond-mat.supr-con math-ph math.MP


title: Testing Self-Organized Criticality Across the Main Sequence using Stellar Flares from TESS

abstract: Self-organized criticality describes a class of dynamical systems that maintain themselves in an attractor state with no intrinsic length or time scale. Fundamentally, this theoretical construct requires a mechanism for instability that may trigger additional instabilities locally via dissipative processes. This concept has been invoked to explain nonlinear dynamical phenomena such as featureless energy spectra that have been observed empirically for earthquakes, avalanches, and solar flares. If this interpretation proves correct, it implies that the solar coronal magnetic field maintains itself in a critical state via a delicate balance between the dynamo-driven injection of magnetic energy and the release of that energy via flaring events. All-sky high-cadence surveys like the Transiting Exoplanet Survey Satellite (TESS) provide the necessary data to compare the energy distribution of flaring events in stars of different spectral types to that observed in the Sun. We identified $\sim 10^6$ flaring events on $\sim 10^5$ stars observed by TESS at 2-minute cadence. By fitting the flare frequency distribution for different mass bins, we find that all main sequence stars exhibit distributions of flaring events similar to that observed in the Sun, independent of their mass or age. This may suggest that stars universally maintain a critical state in their coronal topologies via magnetic reconnection events. If this interpretation proves correct, we may be able to infer properties of magnetic fields, interior structure, and dynamo mechanisms for stars that are otherwise unresolved point sources.


title: Searching for Representation: A sociotechnical audit of googling for members of U.S. Congress

abstract: High-quality online civic infrastructure is increasingly critical for the success of democratic processes. There is a pervasive reliance on search engines to find facts and information necessary for political participation and oversight. We find that approximately 10\% of the top Google search results are likely to mislead California information seekers who use search to identify their congressional representatives. 70\% of the misleading results appear in featured snippets above the organic search results. We use both qualitative and quantitative methods to understand what aspects of the information ecosystem lead to this sociotechnical breakdown. Factors identified include Google's heavy reliance on Wikipedia, the lack of authoritative, machine parsable, high accuracy data about the identity of elected officials based on geographic location, and the search engine's treatment of under-specified queries. We recommend steps that Google can take to meet its stated commitment to providing high quality civic information, and steps that information providers can take to improve the legibility and quality of information about congressional representatives available to search algorithms.


title: The relativistic Pythagorean three-body problem

abstract: We study the influence of relativity on the chaotic properties and dynamical outcomes of an unstable triple system; the Pythagorean three-body problem. To this end, we extend the Brutus N-body code to include Post-Newtonian pairwise terms up to 2.5 order, and the first order Taylor expansion to the Einstein-Infeld-Hoffmann equations of motion. The degree to which our system is relativistic depends on the scaling of the total mass (the unit size was 1 parsec). Using the Brutus method of convergence, we test for time-reversibility in the conservative regime, and demonstrate that we are able to obtain definitive solutions to the relativistic three-body problem. It is also confirmed that the minimal required numerical accuracy for a successful time-reversibility test correlates with the amplification factor of an initial perturbation. When we take into account dissipative effects through gravitational wave emission, we find that the duration of the resonance, and the amount of exponential growth of small perturbations depend on the mass scaling. For a unit mass <= 10 MSun, the system behavior is indistinguishable from the Newtonian case, and the resonance always ends in a binary and one escaping body. For a mass scaling up to 1e7 MSun, relativity gradually becomes more prominent, but the majority of the systems still dissolve. The first mergers start to appear for a mass of ~1e5 MSun, and between 1e7 MSun and 1e9 MSun all systems end prematurely in a merger. These mergers are preceded by a gravitational wave driven in-spiral. For a mass scaling >= 1e9 MSun, all systems result in a gravitational wave merger upon the first close encounter. Relativistic three-body encounters thus provide an efficient pathway for resolving the final parsec problem. The onset of mergers at the characteristic mass scale of 1e7 MSun potentially leaves an imprint in the mass function of supermassive black holes.


title: Smooth solutions to the heat equation which are nowhere analytic in time

abstract: The existence of smooth but nowhere analytic functions is well-known (du Bois-Reymond, Math. Ann., 21(1):109-117, 1883). However, smooth solutions to the heat equation are usually analytic in the space variable. It is also well-known (Kowalevsky, Crelle, 80:1-32, 1875) that a solution to the heat equation may not be time-analytic at $t=0$ even if the initial function is real analytic. Recently, it was shown in \cite{Zha20, DZ20, DP20} that solutions to the heat equation in the whole space, or half space with zero boundary value, are analytic in time under essentially optimal conditions. In this paper, we show that time analyticity is not always true in domains with general boundary conditions or without suitable growth conditions. More precisely, we construct two bounded solutions to the heat equation in the half plane which are nowhere analytic in time. In addition, for any $\delta>0$, we find a solution to the heat equation on the whole plane, with exponential growth of order $2+\delta$, which is nowhere analytic in time.


title: Non-uniqueness in law of three-dimensional magnetohydrodynamics system forced by random noise

abstract: We prove non-uniqueness in law of the three-dimensional magnetohydrodynamics system that is forced by random noise of an additive and a linear multiplicative type and has viscous and magnetic diffusion, both of which are weaker than a full Laplacian. We apply convex integration to both equations of velocity and magnetic fields in order to obtain the non-uniqueness in law in the class of probabilistically strong solutions.


title: Graph Embedding via Diffusion-Wavelets-Based Node Feature Distribution Characterization

abstract: Recent years have seen a rise in the development of representational learning methods for graph data. Most of these methods, however, focus on node-level representation learning at various scales (e.g., microscopic, mesoscopic, and macroscopic node embedding). In comparison, methods for representation learning on whole graphs are currently relatively sparse. In this paper, we propose a novel unsupervised whole graph embedding method. Our method uses spectral graph wavelets to capture topological similarities on each k-hop sub-graph between nodes and uses them to learn embeddings for the whole graph. We evaluate our method against 12 well-known baselines on 4 real-world datasets and show that our method achieves the best performance across all experiments, outperforming the current state-of-the-art by a considerable margin.


title: Written Justifications are Key to Aggregate Crowdsourced Forecasts

abstract: This paper demonstrates that aggregating crowdsourced forecasts benefits from modeling the written justifications provided by forecasters. Our experiments show that the majority and weighted vote baselines are competitive, and that the written justifications are beneficial to call a question throughout its life except in the last quarter. We also conduct an error analysis shedding light into the characteristics that make a justification unreliable.


title: Discretization-independent surrogate modeling over complex geometries using hypernetworks and implicit representations

abstract: Numerical solutions of partial differential equations (PDEs) require expensive simulations, limiting their application in design optimization, model-based control, and large-scale inverse problems. Surrogate modeling techniques seek to decrease the computational expense while retaining dominant solution features and behavior. Traditional Convolutional Neural Network-based frameworks for surrogate modeling require lossy pixelization and data-preprocessing, and generally are not effective in realistic engineering applications. We propose alternative deep-learning based surrogate models for discretization-independent, continuous representations of PDE solutions, which can be used for learning and prediction over domains with complex, variable geometry and mesh topology. Three methods are proposed and compared; design-variable-coded multi-layer perceptron (DV-MLP), design-variable hypernetworks (DV-Hnet), and non-linear independent dual system (NIDS). Each method utilizes a main network which consumes pointwise spatial information to provide a continuous representation, allowing predictions at any location in the domain. Input features include a minimum-distance function evaluation to implicitly encode the problem geometry. The geometric design variables, which define and distinguish problem instances, are used differently by each method, appearing as additional main-network input features (DV-MLP), or as hypernetwork inputs (DV-Hnet and NIDS). The methods are applied to predict solutions around complex, parametrically-defined geometries on non-parametrically-defined meshes with model predictions obtained many orders of magnitude faster than the full order models. Test cases include a vehicle-aerodynamics problem with complex geometry and limited training data, with a design-variable hypernetwork performing best, with a competitive time-to-best-model despite a much greater parameter count.


title: Coarse-Graining of Observables

abstract: We first define the coarse-graining of probability measures in terms of stochastic kernels. We define when a probability measure is part of another probability measure and say that two probability measures coexist if they are both parts of a single probability measure. We then show that any two probability measures coexist. We extend these concepts to observables and instruments and mention that two observables need not coexist. We define the discretization of an observable as a special case of coarse-graining and show that these have \zeroone stochastic kernels. We next consider finite observables and instruments and show that in these cases, stochastic kernels are replaced by stochastic matrices. We also show that coarse-graining is the same as post-processing in this finite case. We then consider sequential products of observables and discuss the sequential product of a post-processed observable with another observable. We briefly discuss SIC observables and the example of qubit observables.


title: Frequency Effects on Syntactic Rule Learning in Transformers

abstract: Pre-trained language models perform well on a variety of linguistic tasks that require symbolic reasoning, raising the question of whether such models implicitly represent abstract symbols and rules. We investigate this question using the case study of BERT's performance on English subject-verb agreement. Unlike prior work, we train multiple instances of BERT from scratch, allowing us to perform a series of controlled interventions at pre-training time. We show that BERT often generalizes well to subject-verb pairs that never occurred in training, suggesting a degree of rule-governed behavior. We also find, however, that performance is heavily influenced by word frequency, with experiments showing that both the absolute frequency of a verb form, as well as the frequency relative to the alternate inflection, are causally implicated in the predictions BERT makes at inference time. Closer analysis of these frequency effects reveals that BERT's behavior is consistent with a system that correctly applies the SVA rule in general but struggles to overcome strong training priors and to estimate agreement features (singular vs. plural) on infrequent lexical items.


title: Pure quantum extension of the semiclassical Boltzmann-Uehling-Uhlenbeck equation

abstract: The Boltzmann equation is the traditional framework in which one extends the time-dependent mean field classical description of a many-body system to include the effect of particle-particle collisions in an approximate manner. A semiclassical extension of this approach to quantum many body systems was suggested by Uehling and Uhlenbeck in 1933 for both Fermi and Bose statistics, and many further developments of this approach are known as the BoltzmannUehling-Uhlenbeck (BUU) equations. Here I introduce a pure quantum version of the BUU type of equations, which is mathematically equivalent to a generalized Time-Dependent Density Functional Theory extended to superfluid systems. As expected, during non-equilibrium processes the quantum Boltzmann one body entropy increases during evolution.


title: How Does Counterfactually Augmented Data Impact Models for Social Computing Constructs?

abstract: As NLP models are increasingly deployed in socially situated settings such as online abusive content detection, it is crucial to ensure that these models are robust. One way of improving model robustness is to generate counterfactually augmented data (CAD) for training models that can better learn to distinguish between core features and data artifacts. While models trained on this type of data have shown promising out-of-domain generalizability, it is still unclear what the sources of such improvements are. We investigate the benefits of CAD for social NLP models by focusing on three social computing constructs -- sentiment, sexism, and hate speech. Assessing the performance of models trained with and without CAD across different types of datasets, we find that while models trained on CAD show lower in-domain performance, they generalize better out-of-domain. We unpack this apparent discrepancy using machine explanations and find that CAD reduces model reliance on spurious features. Leveraging a novel typology of CAD to analyze their relationship with model performance, we find that CAD which acts on the construct directly or a diverse set of CAD leads to higher performance.


title: Embedding Node Structural Role Identity Using Stress Majorization

abstract: Nodes in networks may have one or more functions that determine their role in the system. As opposed to local proximity, which captures the local context of nodes, the role identity captures the functional "role" that nodes play in a network, such as being the center of a group, or the bridge between two groups. This means that nodes far apart in a network can have similar structural role identities. Several recent works have explored methods for embedding the roles of nodes in networks. However, these methods all rely on either approximating or indirect modeling of structural equivalence. In this paper, we present a novel and flexible framework using stress majorization, to transform the high-dimensional role identities in networks directly (without approximation or indirect modeling) to a low-dimensional embedding space. Our method is also flexible, in that it does not rely on specific structural similarity definitions. We evaluated our method on the tasks of node classification, clustering, and visualization, using three real-world and five synthetic networks. Our experiments show that our framework achieves superior results than existing methods in learning node role representations.


title: DPMPC-Planner: A real-time UAV trajectory planning framework for complex static environments with dynamic obstacles

abstract: Safe UAV navigation is challenging due to the complex environment structures, dynamic obstacles, and uncertainties from measurement noises and unpredictable moving obstacle behaviors. Although plenty of recent works achieve safe navigation in complex static environments with sophisticated mapping algorithms, such as occupancy map and ESDF map, these methods cannot reliably handle dynamic environments due to the mapping limitation from moving obstacles. To address the limitation, this paper proposes a trajectory planning framework to achieve safe navigation considering complex static environments with dynamic obstacles. To reliably handle dynamic obstacles, we divide the environment representation into static mapping and dynamic object representation, which can be obtained from computer vision methods. Our framework first generates a static trajectory based on the proposed iterative corridor shrinking algorithm. Then, reactive chance-constrained model predictive control with temporal goal tracking is applied to avoid dynamic obstacles with uncertainties. The simulation results in various environments demonstrate the ability of our algorithm to navigate safely in complex static environments with dynamic obstacles.


title: Globally-Attractive Logarithmic Geometric Control of a Quadrotor for Aggressive Trajectory Tracking

abstract: We present a new quadrotor geometric control scheme that is capable of tracking highly aggressive trajectories. Our geometric controller uses the logarithmic map of SO(3) to express rotational error in the Lie algebra, and we show that it is globally attractive without requiring a complicated hybrid switching scheme. We show the performance of our controller against highly aggressive trajectories in simulation experiments. Additionally, we present an adaptation of this controller that allows us to interface effectively with the angular rate controllers on an onboard flight control unit and show the ability of this adapted control scheme to track aggressive trajectories on a quadrotor hardware platform.


title: Varying the medium surrounding an optical resonator: An efficient and rigorous way to calculate its spectral changes

abstract: Finding reliably and efficiently the spectrum of the resonant states of an optical system under varying parameters of the medium surrounding it is a technologically important task, primarily due to various sensing applications. Computationally, it presents, however, a fundamental challenge owing to the nature of the eigenstates of an open system lacking completeness outside it. We solve this challenge by making a linear transformation of Maxwell's equations which maps perturbations of the surrounding medium onto effective perturbations within the system where the resonant states are complete. By treating such perturbations with the rigorous resonant state expansion, we find the modified modes of the system for arbitrary perturbations of the medium with any required accuracy. Numerically efficient single and few mode approximations are shown to be precise in practically important cases of, respectively, plasmonic nanoparticles and dielectric micro resonators.


title: Guaranteed Safe Spacecraft Docking with Control Barrier Functions

abstract: This paper presents a strategy for control of a spacecraft docking with a non-maneuvering target in the presence of safety constraints and bounded disturbances. The presence of disturbances prevents convergence to a unique docking state, so in our formulation, docking is defined as occurring within a set constructed using prescribed tolerances. Safety is ensured via application of Robust Control Barrier Functions to render a designated safe set forward invariant for any allowable disturbance. However, this safety strategy necessarily presumes a worst-case disturbance, and thus restricts trajectories to a subset of the safe set when a worst-case disturbance is not present. The presented controller accounts for this restriction, and guarantees that the spacecraft both remains safe and achieves docking in finite time for any allowable disturbance. The controller is then validated in simulation for a spacecraft landing on an asteroid, and two spacecraft docking in low Earth orbit.