title: Augmentations, annuli, and Alexander polynomials

abstract: The augmentation variety of a knot is the locus, in the 3-dimensional coefficient space of the knot contact homology dg-algebra, where the algebra admits a unital chain map to the complex numbers. We explain how to express the Alexander polynomial of a knot in terms of the augmentation variety: it is the exponential of the integral of a ratio of two partial derivatives. The expression is derived from a description of the Alexander polynomial as a count of Floer strips and holomorphic annuli, in the cotangent bundle of Euclidean 3-space, stretching between a Lagrangian with the topology of the knot complement and the zero-section, and from a description of the boundary of the moduli space of such annuli with one positive puncture.


title: SPARCQ: A new approach for fat fraction mapping using asymmetries in the phase-cycled bSSFP signal profile

abstract: Purpose: To develop SPARCQ (Signal Profile Asymmetries for Rapid Compartment Quantification), a novel approach to quantify fat fraction using asymmetries in the phase-cycled bSSFP profile. Methods: SPARCQ uses phase-cycling to obtain bSSFP frequency profiles, which display asymmetries in the presence of fat and water at certain repetition times. For each voxel, the measured signal profile is decomposed into a weighted sum of simulated profiles via multi-compartment dictionary matching. Each dictionary entry represents a single-compartment bSSFP profile with a specific off-resonance frequency and relaxation time ratio. Using the results of dictionary matching, the fractions of the different off-resonance components are extracted for each voxel, generating quantitative maps of water and fat fraction and banding-artifact-free images for the entire image volume. SPARCQ was validated using simulations, experiments in a water-fat phantom and in knees of healthy volunteers. Experimental results were compared with reference proton density fat fractions obtained with 1H-MRS (phantoms) and with multiecho gradient-echo MRI (phantoms and volunteers). SPARCQ repeatability was evaluated in six scan-rescan experiments. Results: Simulations showed that fat fraction quantification is accurate and robust for signal-to-noise ratios greater than 20. Phantom experiments demonstrated good agreement between SPARCQ and gold standard fat fractions. In volunteers, banding-artifact-free quantitative maps and water-fat-separated images obtained with SPARCQ and ME-GRE demonstrated the expected contrast between fatty and non-fatty tissues. The coefficient of repeatability of SPARCQ fat fraction was 0.0512. Conclusion: SPARCQ demonstrates potential for fat quantification using asymmetries in bSSFP profiles, and may be a promising alternative to conventional fat fraction quantification techniques.


title: Self-Excited Ising Game

abstract: Effects of dynamical activity spillover in a noisy binary choice game (Ising game) on a complete graph are studied. Binary choice games are very important for both economics and statistical physics playing a role of the bridge between these two fields. In this paper we investigate the effects of self-excited activity induced by activity spillover on relaxation to equilibria and transitions between metastable equilibria at finite times. Using the formalism of master equations we show that both relaxation and interequilibria transitions at finite time are accelerated by the effects of activity spillover.


title: Hardware-Accelerated SAR Simulation with NVIDIA-RTX Technology

abstract: Synthetic Aperture Radar (SAR) is a critical sensing technology that is notably independent of the sensor-to-target distance and has numerous cross-cutting applications, e.g., target recognition, mapping, surveillance, oceanography, geology, forestry (biomass, deforestation), disaster monitoring (volcano eruptions, oil spills, flooding), and infrastructure tracking (urban growth, structure mapping). SAR uses a high-power antenna to illuminate target locations with electromagnetic radiation, e.g., 10GHz radio waves, and illuminated surface backscatter is sensed by the antenna which is then used to generate images of structures. Real SAR data is difficult and costly to produce and, for research, lacks a reliable source ground truth. This article proposes a open source SAR simulator to compute phase histories for arbitrary 3D scenes using newly available ray-tracing hardware made available commercially through the NVIDIA's RTX graphics cards series. The OptiX GPU ray tracing library for NVIDIA GPUs is used to calculate SAR phase histories at unprecedented computational speeds. The simulation results are validated against existing SAR simulation code for spotlight SAR illumination of point targets. The computational performance of this approach provides orders of magnitude speed increases over CPU simulation. An additional order of magnitude of GPU acceleration when simulations are run on RTX GPUs which include hardware specifically to accelerate OptiX ray tracing. The article describes the OptiX simulator structure, processing framework and calculations that afford execution on massively parallel GPU computation device. The shortcoming of the OptiX library's restriction to single precision float representation is discussed and modifications of sensitive calculations are proposed to reduce truncation error thereby increasing the simulation accuracy under this constraint.


title: Laminar separated flows over finite-aspect-ratio swept wings

abstract: We perform direct numerical simulations of laminar separated flows over finite-aspect-ratio swept wings at a chord-base Reynolds number of $Re = 400$ to reveal a variety of wake structures generated for a range of aspect ratios ($sAR=0.5-4$), angles of attack ($\alpha=16^{\circ}-30^{\circ}$), and sweep angles ($\Lambda=0^{\circ}-45^{\circ}$). Flows behind swept wings exhibit increased complexity in their dynamical features compared to unswept-wing wakes. For unswept wings, the wake dynamics are predominantly influenced by the tip effects. Steady wakes are mainly limited to low-aspect-ratio wings. Unsteady vortex shedding takes place near the midspan of higher-$AR$ wings due to weakened downwash induced by the tip vortices. With increasing sweep angle, the source of three dimensionality transitions from the tip to the midspan. A pair of symmetric vortical structures forms along the two sides of the midspan, imposing downward velocity upon each other in a manner similar to the downwash induced by tip vortices. Such stabilizing midspan effects not only expand the steady wake region to higher aspect ratios, but also enhance lift. At higher aspect ratios, the midspan effects of swept wings diminish at outboard region, allowing unsteady vortex shedding to develop near the tip. In the wakes of highly swept wings, streamwise finger-like structures form repetitively along the wing span, providing a stabilizing effect. The insights revealed from this study can aid the design of high-lift devices and serve as a stepping stone for understanding the complex wake dynamics at higher Reynolds numbers and those generated by unsteady wing maneuvers.


title: Matching methods for obtaining survival functions to estimate the effect of a time-dependent treatment

abstract: In observational studies of survival time featuring a binary time-dependent treatment, the hazard ratio (an instantaneous measure) is often used to represent the treatment effect. However, investigators are often more interested in the difference in survival functions. We propose semiparametric methods to estimate the causal effect of treatment among the treated with respect to survival probability. The objective is to compare post-treatment survival with the survival function that would have been observed in the absence of treatment. For each patient, we compute a prognostic score (based on the pre-treatment death hazard) and a propensity score (based on the treatment hazard). Each treated patient is then matched with an alive, uncensored and not-yet-treated patient with similar prognostic and/or propensity scores. The experience of each treated and matched patient is weighted using a variant of Inverse Probability of Censoring Weighting to account for the impact of censoring. We propose estimators of the treatment-specific survival functions (and their difference), computed through weighted Nelson-Aalen estimators. Closed-form variance estimators are proposed which take into consideration the potential replication of subjects across matched sets. The proposed methods are evaluated through simulation, then applied to estimate the effect of kidney transplantation on survival among end-stage renal disease patients using data from a national organ failure registry.


title: The Possibility of a Non-Lagrangian Theory of Gravity

abstract: General Relativity resembles a very elegant crystal glass: If we touch its principles, that is, its Lagrangian, there is a risk of breaking everything. Or, if we will, it is like a short blanket: Curing some problems creates new problems. This paper is devoted to bring to light the reasons why we pursue the possibility of a non-Lagrangian theory of gravity under the hypothesis of an extension of the original general relativity with an ansatz inspired in the fundamental principles of classical and quantum physics.


title: GLEAKE: Global and Local Embedding Automatic Keyphrase Extraction

abstract: Automated methods for granular categorization of large corpora of text documents have become increasingly more important with the rate scientific, news, medical, and web documents are growing in the last few years. Automatic keyphrase extraction (AKE) aims to automatically detect a small set of single or multi-words from within a single textual document that captures the main topics of the document. AKE plays an important role in various NLP and information retrieval tasks such as document summarization and categorization, full-text indexing, and article recommendation. Due to the lack of sufficient human-labeled data in different textual contents, supervised learning approaches are not ideal for automatic detection of keyphrases from the content of textual bodies. With the state-of-the-art advances in text embedding techniques, NLP researchers have focused on developing unsupervised methods to obtain meaningful insights from raw datasets. In this work, we introduce Global and Local Embedding Automatic Keyphrase Extractor (GLEAKE) for the task of AKE. GLEAKE utilizes single and multi-word embedding techniques to explore the syntactic and semantic aspects of the candidate phrases and then combines them into a series of embedding-based graphs. Moreover, GLEAKE applies network analysis techniques on each embedding-based graph to refine the most significant phrases as a final set of keyphrases. We demonstrate the high performance of GLEAKE by evaluating its results on five standard AKE datasets from different domains and writing styles and by showing its superiority with regards to other state-of-the-art methods.


title: Data-driven feedback stabilization of nonlinear systems: Koopman-based model predictive control

abstract: In this work, a predictive control framework is presented for feedback stabilization of nonlinear systems. To achieve this, we integrate Koopman operator theory with Lyapunov-based model predictive control (LMPC). The main idea is to transform nonlinear dynamics from state-space to function space using Koopman eigenfunctions - for control affine systems this results in a bilinear model in the (lifted) function space. Then, a predictive controller is formulated in Koopman eigenfunction coordinates which uses an auxiliary Control Lyapunov Function (CLF) based bounded controller as a constraint to ensure stability of the Koopman system in the function space. Provided there exists a continuously differentiable inverse mapping between the original state-space and (lifted) function space, we show that the designed controller is capable of translating the feedback stabilizability of the Koopman bilinear system to the original nonlinear system. Remarkably, the feedback control design proposed in this work remains completely data-driven and does not require any explicit knowledge of the original system. Furthermore, due to the bilinear structure of the Koopman model, seeking a CLF is no longer a bottleneck for LMPC. Benchmark numerical examples demonstrate the utility of the proposed feedback control design.


title: Characterization of the lengths of binary circular words containing no squares other than 00, 11, and 0101

abstract: We characterize exactly the lengths of binary circular words containing no squares other than 00, 11, and 0101. Key words: combinatorics on words, circular words, necklaces, square-free words, non-repetitive sequences


title: On Evaluating Weakly Supervised Action Segmentation Methods

abstract: Action segmentation is the task of temporally segmenting every frame of an untrimmed video. Weakly supervised approaches to action segmentation, especially from transcripts have been of considerable interest to the computer vision community. In this work, we focus on two aspects of the use and evaluation of weakly supervised action segmentation approaches that are often overlooked: the performance variance over multiple training runs and the impact of selecting feature extractors for this task. To tackle the first problem, we train each method on the Breakfast dataset 5 times and provide average and standard deviation of the results. Our experiments show that the standard deviation over these repetitions is between 1 and 2.5% and significantly affects the comparison between different approaches. Furthermore, our investigation on feature extraction shows that, for the studied weakly-supervised action segmentation methods, higher-level I3D features perform worse than classical IDT features.


title: On a family of discrete log-symmetric distributions

abstract: The use of continuous probability distributions has been widespread in problems with purely discrete nature. In general, such distributions are not appropriate in this scenario. In this paper, we introduce a class of discrete and asymmetric distributions based on the family of continuous log-symmetric distributions. Some properties are discussed as well as estimation by the maximum likelihood method. A Monte Carlo simulation study is carried out to evaluate the performance of the estimators, and censored and uncensored data sets are used to illustrate the proposed methodology.


title: Optimal Resource Allocation for Elastic and Inelastic Jobs

abstract: Modern data centers are tasked with processing heterogeneous workloads consisting of various classes of jobs. These classes differ in their arrival rates, size distributions, and job parallelizability. With respect to paralellizability, some jobs are elastic, meaning they can parallelize linearly across many servers. Other jobs are inelastic, meaning they can only run on a single server. Although job classes can differ drastically, they are typically forced to share a single cluster. When sharing a cluster among heterogeneous jobs, one must decide how to allocate servers to each job at every moment in time. In this paper, we design and analyze allocation policies which aim to minimize the mean response time across jobs, where a job's response time is the time from when it arrives until it completes. We model this problem in a stochastic setting where each job may be elastic or inelastic. Job sizes are drawn from exponential distributions, but are unknown to the system. We show that, in the common case where elastic jobs are larger on average than inelastic jobs, the optimal allocation policy is Inelastic-First, giving inelastic jobs preemptive priority over elastic jobs. We obtain this result by introducing a novel sample path argument. We also show that there exist cases where Elastic-First (giving priority to elastic jobs) performs better than Inelastic-First. We then provide the first analysis of mean response time under both Elastic-First and Inelastic-First by leveraging recent techniques for solving high-dimensional Markov chains.


title: Bridging physics-based and equivalent circuit models for lithium-ion batteries

abstract: In this article, a novel implementation of a widely used pseudo-two-dimensional (P2D) model for lithium-ion battery simulation is presented with a transmission line circuit structure. This implementation represents an interplay between physical and equivalent circuit models. The discharge processes of an NMC-graphite lithium-ion battery under different currents are simulated, and it is seen the results from the circuit model agree well with the results obtained from a physical simulation carried out in COMSOL Multiphysics, including both terminal voltage and concentration distributions. Finally we demonstrated how the circuit model can contribute to the understanding of the cell electrochemistry, exemplified by an analysis of the overpotential contributions by various processes.


title: An Efficient Machine-Learning Approach for PDF Tabulation in Turbulent Combustion Closure

abstract: Probability density function (PDF) based turbulent combustion modelling is limited by the need to store multi-dimensional PDF tables that can take up large amounts of memory. A significant saving in storage can be achieved by using various machine-learning techniques that represent the thermo-chemical quantities of a PDF table using mathematical functions. These functions can be computationally more expensive than the existing interpolation methods used for thermo-chemical quantities. More importantly, the training time can amount to a considerable portion of the simulation time. In this work, we address these issues by introducing an adaptive training algorithm that relies on multi-layer perception (MLP) neural networks for regression and self-organizing maps (SOMs) for clustering data to tabulate using different networks. The algorithm is designed to address both the multi-dimensionality of the PDF table as well as the computational efficiency of the proposed algorithm. SOM clustering divides the PDF table into several parts based on similarities in data. Each cluster of data is trained using an MLP algorithm on simple network architectures to generate local functions for thermo-chemical quantities. The algorithm is validated for the so-called DLR-A turbulent jet diffusion flame using both RANS and LES simulations and the results of the PDF tabulation are compared to the standard linear interpolation method. The comparison yields a very good agreement between the two tabulation techniques and establishes the MLP-SOM approach as a viable method for PDF tabulation.


title: The Virtual Block Interface: A Flexible Alternative to the Conventional Virtual Memory Framework

abstract: Computers continue to diversify with respect to system designs, emerging memory technologies, and application memory demands. Unfortunately, continually adapting the conventional virtual memory framework to each possible system configuration is challenging, and often results in performance loss or requires non-trivial workarounds. To address these challenges, we propose a new virtual memory framework, the Virtual Block Interface (VBI). We design VBI based on the key idea that delegating memory management duties to hardware can reduce the overheads and software complexity associated with virtual memory. VBI introduces a set of variable-sized virtual blocks (VBs) to applications. Each VB is a contiguous region of the globally-visible VBI address space, and an application can allocate each semantically meaningful unit of information (e.g., a data structure) in a separate VB. VBI decouples access protection from memory allocation and address translation. While the OS controls which programs have access to which VBs, dedicated hardware in the memory controller manages the physical memory allocation and address translation of the VBs. This approach enables several architectural optimizations to (1) efficiently and flexibly cater to different and increasingly diverse system configurations, and (2) eliminate key inefficiencies of conventional virtual memory. We demonstrate the benefits of VBI with two important use cases: (1) reducing the overheads of address translation (for both native execution and virtual machine environments), as VBI reduces the number of translation requests and associated memory accesses; and (2) two heterogeneous main memory architectures, where VBI increases the effectiveness of managing fast memory regions. For both cases, VBI significanttly improves performance over conventional virtual memory.


title: Ring Laser Gyro: Noise Improves Performance

abstract: We present a systematic formulation of the performance of a ring laser gyro (RLG) under the influence of a randomized sinusoidal dither. We develop the Fokker Planck equation of such a system and find its steady state solution that we use to find exact equations for the average output rate and its variance. The former determines the RLG's scale factor and the latter its angular random walk. We find a key parameter, called the residual lock-in rate, which affects both the scale factor and the random walk of the gyro. We show that when the noise amplitude is much larger than the residual lock-in rate the gyro rate output approaches that of an ideal gyro. However, this is at the expense of generating an angular random walk that is proportional to the residual lock-in rate. We present a formulation of the statistics of gyro output count and relate the average count and its standard deviation to gyro parameters. For example, we show that the angular random walk, in addition to being proportional to the residual lock-in rate, is also proportional to the inverse square root of the dither frequency.


title: Spectrum-free integrated photonic remote molecular identification and sensing

abstract: Absorption spectroscopy is widely used in sensing and astronomy to understand molecular compositions on microscopic to cosmological scales. However, typical dispersive spectroscopic techniques require multichannel detection, fundamentally limiting the ability to detect extremely weak signals when compared to direct photometric methods. We report the realization of direct spectral molecular detection using a silicon nanophotonic waveguide resonator, obviating dispersive spectral acquisition. We use a thermally tunable silicon ring resonator with a transmission spectrum matched and cross-correlated to the quasi-periodic vibronic absorption lines of hydrogen cyanide. We show that the correlation peak amplitude is proportional to the number of overlapping ring resonances and gas lines, and that molecular specificity is obtained from the phase of the correlation signal in a single detection channel. Our results demonstrate on-chip correlation spectroscopy that is less restricted by the signal-to-noise penalty of other spectroscopic approaches, enabling the detection of faint spectral signatures.


title: Phase transitions and control measures for network epidemics caused by infections with presymptomatic, asymptomatic,and symptomatic stages

abstract: We investigate phase transitions associated with three control methods for epidemics on small world networks. Motivated by the behavior of SARS-CoV-2, we construct a theoretical SIR model of a virus that exhibits presymptomatic, asymptomatic, and symptomatic stages in two possible pathways. Using agent-based simulations on small world networks, we observe phase transitions for epidemic spread related to: 1) Global social distancing with a fixed probability of adherence. 2) Individually initiated social isolation when a threshold number of contacts are infected. 3) Viral shedding rate. The primary driver of total number of infections is the viral shedding rate, with probability of social distancing being the next critical factor. Individually initiated social isolation was effective when initiated in response to a single infected contact. For each of these control measures, the total number of infections exhibits a sharp phase transition as the strength of the measure is varied.


title: Learning Representations using Spectral-Biased Random Walks on Graphs

abstract: Several state-of-the-art neural graph embedding methods are based on short random walks (stochastic processes) because of their ease of computation, simplicity in capturing complex local graph properties, scalability, and interpretibility. In this work, we are interested in studying how much a probabilistic bias in this stochastic process affects the quality of the nodes picked by the process. In particular, our biased walk, with a certain probability, favors movement towards nodes whose neighborhoods bear a structural resemblance to the current node's neighborhood. We succinctly capture this neighborhood as a probability measure based on the spectrum of the node's neighborhood subgraph represented as a normalized laplacian matrix. We propose the use of a paragraph vector model with a novel Wasserstein regularization term. We empirically evaluate our approach against several state-of-the-art node embedding techniques on a wide variety of real-world datasets and demonstrate that our proposed method significantly improves upon existing methods on both link prediction and node classification tasks.


title: High-Contrast Integral Field Spectrograph (HCIFS): multi-spectral wavefront control and reduced-dimensional system identification

abstract: Any high-contrast imaging instrument in a future large space-based telescope will include an integral field spectrograph (IFS) for measuring broadband starlight residuals and characterizing the exoplanet's atmospheric spectrum. In this paper, we report the development of a high-contrast integral field spectrograph (HCIFS) at Princeton University and demonstrate its application in multi-spectral wavefront control. Moreover, we propose and experimentally validate a new reduced-dimensional system identification algorithm for an IFS imaging system, which improves the system's wavefront control speed, contrast and computational and data storage efficiency.


title: Non-twist tori in conformally symplectic systems

abstract: Dissipative mechanical systems on the torus with a friction that is proportional to the velocity are modeled by conformally symplectic maps on the annulus, which are maps that transport the symplectic form into a multiple of itself (with a conformal factor smaller than 1). It is important to understand the structure and the dynamics on the attractors. With the aid of parameters, and under suitable non-degeneracy conditions, one can obtain that, by adjusting parameters, there is an attractor that is an invariant torus whose internal dynamics is conjugate to a rotation [CCdlL13]. By analogy with symplectic dynamics, there have been some debate in establishing appropriate definitions for twist and non-twist invariant tori (or systems). The purpose of this paper is two-fold: (a) to establish proper definitions of twist and non-twist invariant tori in families of conformally symplectic systems; (b) to derive algorithms of computation of non-twist invariant tori. The last part of the paper is devoted to implementations of the algorithms, illustrating the definitions presented in this paper, and exploring the mechanisms of breakdown of non-twist tori. For the sake of simplicity we have considered here 2D systems, i.e. defined in the 2D annulus, but generalization to higher dimensions is straightforward.


title: Adapting a Kidney Exchange Algorithm to Align with Human Values

abstract: The efficient and fair allocation of limited resources is a classical problem in economics and computer science. In kidney exchanges, a central market maker allocates living kidney donors to patients in need of an organ. Patients and donors in kidney exchanges are prioritized using ad-hoc weights decided on by committee and then fed into an allocation algorithm that determines who gets what--and who does not. In this paper, we provide an end-to-end methodology for estimating weights of individual participant profiles in a kidney exchange. We first elicit from human subjects a list of patient attributes they consider acceptable for the purpose of prioritizing patients (e.g., medical characteristics, lifestyle choices, and so on). Then, we ask subjects comparison queries between patient profiles and estimate weights in a principled way from their responses. We show how to use these weights in kidney exchange market clearing algorithms. We then evaluate the impact of the weights in simulations and find that the precise numerical values of the weights we computed matter little, other than the ordering of profiles that they imply. However, compared to not prioritizing patients at all, there is a significant effect, with certain classes of patients being (de)prioritized based on the human-elicited value judgments.


title: Improving Proper Noun Recognition in End-to-End ASR By Customization of the MWER Loss Criterion

abstract: Proper nouns present a challenge for end-to-end (E2E) automatic speech recognition (ASR) systems in that a particular name may appear only rarely during training, and may have a pronunciation similar to that of a more common word. Unlike conventional ASR models, E2E systems lack an explicit pronounciation model that can be specifically trained with proper noun pronounciations and a language model that can be trained on a large text-only corpus. Past work has addressed this issue by incorporating additional training data or additional models. In this paper, we instead build on recent advances in minimum word error rate (MWER) training to develop two new loss criteria that specifically emphasize proper noun recognition. Unlike past work on this problem, this method requires no new data during training or external models during inference. We see improvements ranging from 2% to 7% relative on several relevant benchmarks.


title: Studying Type II supernovae as cosmological standard candles using the Dark Energy Survey

abstract: Despite vast improvements in the measurement of the cosmological parameters, the nature of dark energy and an accurate value of the Hubble constant (H$_0$) in the Hubble-Lema\^itre law remain unknown. To break the current impasse, it is necessary to develop as many independent techniques as possible, such as the use of Type II supernovae (SNe II). The goal of this paper is to demonstrate the utility of SNe II for deriving accurate extragalactic distances, which will be an asset for the next generation of telescopes where more-distant SNe II will be discovered. More specifically, we present a sample from the Dark Energy Survey Supernova Program (DES-SN) consisting of 15 SNe II with photometric and spectroscopic information spanning a redshift range up to 0.35. Combining our DES SNe with publicly available samples, and using the standard candle method (SCM), we construct the largest available Hubble diagram with SNe II in the Hubble flow (70 SNe II) and find an observed dispersion of 0.27 mag. We demonstrate that adding a colour term to the SN II standardisation does not reduce the scatter in the Hubble diagram. Although SNe II are viable as distance indicators, this work points out important issues for improving their utility as independent extragalactic beacons: find new correlations, define a more standard subclass of SNe II, construct new SN II templates, and dedicate more observing time to high-redshift SNe II. Finally, for the first time, we perform simulations to estimate the redshift-dependent distance-modulus bias due to selection effects.


title: Predictor Antenna Systems: Exploiting Channel State Information for Vehicle Communications

abstract: Vehicle communication is one of the most important use cases in the fifth generation of wireless networks (5G). The growing demand for quality of service (QoS) characterized by performance metrics, such as spectrum efficiency, peak data rate, and outage probability, is mainly limited by inaccurate prediction/estimation of channel state information (CSI) of the rapidly changing environment around moving vehicles. One way to increase the prediction horizon of CSI in order to improve the QoS is deploying predictor antennas (PAs). A PA system consists of two sets of antennas typically mounted on the roof of a vehicle, where the PAs positioned at the front of the vehicle are used to predict the CSI observed by the receive antennas (RAs) that are aligned behind the PAs. In realistic PA systems, however, the actual benefit is affected by a variety of factors, including spatial mismatch, antenna utilization, temporal correlation of scattering environment, and CSI estimation error. This thesis investigates different resource allocation schemes for the PA systems under practical constraints.


title: Observational appearance of rapidly rotating neutron stars: X-ray bursts, cooling tail method, and radius determination

abstract: Neutron stars (NSs) in low-mass X-ray binaries rotate at frequencies high enough to significantly deviate from sphericity ($\nu_* \sim$ 200--600 Hz). We investigate the effects of rapid rotation on the observational appearance of a NS. We propose analytical formulae relating gravitational mass and equatorial radius of the rapidly rotating NS to the mass $M$ and radius $R$ of a non-rotating NS of the same baryonic mass using accurate fully relativistic computations. We compute spectra from an oblate rotating NS observed at different inclination angles using the modified oblate Schwarzschild (MOS) approximation, where light bending is computed in Schwarzschild metric, but frame dragging and quadrupole moment of a NS are approximately accounted for in the photon redshift calculations. We generalize the cooling tail method to the case of a rapidly rotating NS to obtain the most probable values of $M$ and $R$ of the corresponding non-rotating NS with the same baryonic mass. We approximate the local spectra from the NS surface by a diluted blackbody using previously computed NS atmosphere models. We show that the NS radius could be overestimated by 3--3.5 km for face-on stars of $R\approx 11$ km rotating at $\nu_* =$ 700 Hz if the version of the cooling tail method for a non-rotating NS is used. We apply the method to an X-ray burst observed from the NS rotating at $\nu_* \approx$ 532 Hz in SAX J1810.8$-$2609. The resulting radius of the non-rotating NS (assuming $M=1.5 M_\odot$) becomes $11.8\pm0.5$ km if it is viewed at inclination i=60 deg and $R=11.2\pm0.5$ km for a face-on view, which are smaller by 0.6 and 1.2 km than the radius obtained using standard cooling tail method ignoring rotation. The corresponding equatorial radii of these rapidly rotating NSs are 12.3$\pm 0.6$ km (for i=60 deg) and 11.6$\pm 0.6$\,km (for i=0 deg).


title: Some remarks on a coupling method for the practical computation of homogenized coefficients

abstract: We numerically investigate, and improve upon, a computational approach originally introduced in [Cottereau, IJNME 2013] which aims at evaluating the effective coefficient of a medium modelled by a highly oscillatory coefficient. This computational approach is based on a Arlequin type coupling. It combines the original fine-scale description of the medium (modelled by an oscillatory coefficient) with an effective description (modelled by a constant coefficient) and optimizes upon the coefficient of the effective medium to best fit the response of the actual heterogeneous medium using a purely homogeneous medium. We present here a mathematical formalization of the approach along with various improvements of the algorithms, in order to obtain a procedure as efficient as possible. Representative numerical results demonstrate the added value of our approach in comparison to the original approach.


title: Analytical correspondence between shadow radius and black hole quasinormal frequencies

abstract: We consider the equivalence of quasinormal modes and geodesic quantities recently brought back due to the black hole shadow observation by Event Horizon Telescope. Using WKB method we found an analytical relation between the real part of quasinormal frequencies at the eikonal limit and black hole shadow radius. We verify this correspondence with two black hole families in $4$ and $D$ dimensions, respectively.


title: Digraph Signal Processing with Generalized Boundary Conditions

abstract: Signal processing on directed graphs (digraphs) is problematic, since the graph shift, and thus associated filters, are in general not diagonalizable. Furthermore, the Fourier transform in this case is now obtained from the Jordan decomposition, which may not be computable at all for large graphs. We propose a novel and general solution for this problem based on matrix perturbation theory: We design an algorithm that adds a small number of edges to a given digraph to destroy nontrivial Jordan blocks. The obtained digraph is then diagonalizable and yields, as we show, an approximate eigenbasis and Fourier transform for the original digraph. We explain why and how this construction can be viewed as generalized form of boundary conditions, a common practice in signal processing. Our experiments with random and real world graphs show that we can scale to graphs with a few thousands nodes, and obtain Fourier transforms that are close to orthogonal while still diagonalizing an intuitive notion of convolution. Our method works with adjacency and Laplacian shift and can be used as preprocessing step to enable further processing as we show with a prototypical Wiener filter application.


title: Interplay between coherent and dissipative dynamics of bosonic doublons in an optical lattice

abstract: We observe the dissipative dynamics of a dense, strongly interacting gas of bosonic atom pairs in an optical lattice, controlling the strength of the two-body interactions over a wide parameter regime. We study how three-body losses contribute to the lattice dynamics, addressing a number of open questions related to the effects of strong dissipation in a many-body system, including the relationship to the continuous quantum Zeno effect. We observe rapid break-up of bound pairs for weak interactions, and for stronger interactions we see doublon decay rates that are asymmetric when changing from attractive and repulsive interactions, and which strongly depend on the interactions and on-site loss rates. By comparing our experimental data with a theoretical analysis of few-body dynamics, we show that these features originate from a non-trivial combination of dissipative dynamics described by a lattice model beyond a standard Bose-Hubbard Hamiltonian, and the modification of three-atom dynamics on a single site, which is generated alongside strong three-body loss. Our results open new possibilities for investigating bosonic atoms with strong three-body loss features, and allow for the better understanding of the parameter regimes that are required to realize strong effective three-body interactions.


title: Debye screening for the stationary Vlasov-Poisson equation in interaction with a point charge

abstract: We prove that the Debye screening length emerges in a spatially homogeneous plasma described by the nonlinear Vlasov-Poisson equation interacting with a point charge. While screened stationary states as well as the time-dependent problem are well-studied for the linearized equation, there are few rigorous results on screening in the nonlinear setting. As such, the results presented here cover the stationary case under the assumptions predicted by the linearized theory.


title: Vortices in Maxwell-Chern-Simons-Higgs models with nonminimal coupling

abstract: We investigate the presence of vortex configurations in generalized Maxwell-Chern-Simons models with nonminimal coupling, in which we introduce a function that modifies the dynamical term of the scalar field in the Lagrangian. We first follow a route already considered in previous works to develop the Bogomol'nyi procedure, and, in this context, we use the first order equations to obtain a vortex with a novel behavior at its core. We then go further and introduce a novel procedure to develop the Bogomol'nyi methodology. It supports distinct first order equations, and we then investigate another model, in which the vortex may engender inversion of the magnetic flux, an effect with no precedents in the study of vortices within the nonminimal context.


title: Self-supervised Dynamic CT Perfusion Image Denoising with Deep Neural Networks

abstract: Dynamic computed tomography perfusion (CTP) imaging is a promising approach for acute ischemic stroke diagnosis and evaluation. Hemodynamic parametric maps of cerebral parenchyma are calculated from repeated CT scans of the first pass of iodinated contrast through the brain. It is necessary to reduce the dose of CTP for routine applications due to the high radiation exposure from the repeated scans, where image denoising is necessary to achieve a reliable diagnosis. In this paper, we proposed a self-supervised deep learning method for CTP denoising, which did not require any high-dose reference images for training. The network was trained by mapping each frame of CTP to an estimation from its adjacent frames. Because the noise in the source and target was independent, this approach could effectively remove the noise. Being free from high-dose training images granted the proposed method easier adaptation to different scanning protocols. The method was validated on both simulation and a public real dataset. The proposed method achieved improved image quality compared to conventional denoising methods. On the real data, the proposed method also had improved spatial resolution and contrast-to-noise ratio compared to supervised learning which was trained on the simulation data


title: Many flows in the group connectivity setting

abstract: Two well-known results in the world of nowhere-zero flows are Jaeger's 4-flow theorem asserting that every 4-edge-connected graph has a nowhere-zero $\mathbb{Z}_2 \times \mathbb{Z}_2$-flow and Seymour's 6-flow theorem asserting that every 2-edge-connected graph has a nowhere-zero $\mathbb{Z}_6$-flow. Dvo\v{r}\'ak and the last two authors of this paper extended these results by proving the existence of exponentially many nowhere-zero flows under the same assumptions. We revisit this setting and provide extensions and simpler proofs of these results. The concept of a nowhere-zero flow was extended in a significant paper of Jaeger, Linial, Payan, and Tarsi to a choosability-type setting. For a fixed abelian group $\Gamma$, an oriented graph $G = (V,E)$ is called $\Gamma$-connected if for every function $f : E \rightarrow \Gamma$ there is a flow $\phi : E \rightarrow \Gamma$ with $\phi(e) \neq f(e)$ for every $e \in E$ (note that taking $f = 0$ forces $\phi$ to be nowhere-zero). Jaeger et al. proved that every oriented 3-edge-connected graph is $\Gamma$-connected whenever $|\Gamma| \ge 6$. We prove that there are exponentially many solutions whenever $|\Gamma| \ge 8$. For the group $\mathbb{Z}_6$ we prove that for every oriented 3-edge-connected $G = (V,E)$ with $\ell = |E| - |V| \ge 11$ and every $f: E \rightarrow \mathbb{Z}_6$, there are at least $2^{ \sqrt{\ell} / \log \ell}$ flows $\phi$ with $\phi(e) \neq f(e)$ for every $e \in E$.


title: Perceptual similarity between piano notes: Simulations with a template-based perception model

abstract: In this paper the auditory model developed by Dau et al. [J. Acoust. Soc. Am. 102, 2892-2905 (1997)] was used to simulate the perceptual similarity between complex sounds. For this purpose, a central processor stage was developed and attached as a back-end module to the auditory model. As complex sounds, a set of recordings of one note played on seven different pianos was used, whose similarity has been recently measured by Osses et al. [J. Acoust. Soc. Am. 146, 1024-1035 (2019)] using a 3-AFC discrimination task in noise. The auditory model has several processing stages that are to a greater or lesser extent inspired by physiological aspects of the human normal-hearing system. A set of configurable parameters in each stage affects directly the sound (internal) representations that are further processed in the developed central processor. Therefore, a comprehensive review of the model parameters is given, indicating the configuration we chose. This includes an in-depth description of the auditory adaptation stage, the adaptation loops. Simulations of the similarity task were compared with (1) existing experimental data, where they had a moderate to high correlation, and with (2) simulations using an alternative but similar background noise to that of the experiments, which were used to obtain further information about how the participants' responses were weighted as a function of frequency.


title: Controlling the Hidden Growth of COVID-19

abstract: The COVID-19 pandemic has plagued the world for months. The U.S. has taken measures to counter it. On a daily basis, newly confirmed cases have been reported. In the early days, these numbers showed an increasing trend. Recently, the numbers have been generally flattened out. This report tries to estimate the hidden number of currently alive infections in the population by using the confirmed cases. A major result indicates an existing infections estimate at about 10-50 times the daily confirmed new cases, with the stringent social distancing policy tipping to the upper end of this range. It clarifies the relationship between the infection rate and the test rate to put the epidemic under control, which says that the test rate shall keep up at the same pace as infection rate to prevent an outbreak. This relationship is meaningful in the wake of business re-opening in the U.S. and the world. The report also reveals the connections of all the measures taken to the epidemic spread. A stratified sampling method is proposed to add to the current tool kits of epidemic control. Again, this report is a summary of some straight observations and thoughts, not through a thorough study backed with field data. The results appear obvious and suitable for general education to interested policymakers and the public.


title: Solubility of Additive Sextic Forms over $\mathbb{Q}_2(\sqrt{-1})$ and $\mathbb{Q}_2(\sqrt{-5})$

abstract: Michael Knapp, in a previous work, conjectured that every additive sextic form over $\mathbb{Q}_2(\sqrt{-1})$ and $\mathbb{Q}_2(\sqrt{-5})$ in seven variables has a nontrivial zero. In this paper, we show that this conjecture is true, establishing that $\Gamma^*(6, \mathbb{Q}_2(\sqrt{-1})) = \Gamma^*(6, \mathbb{Q}_2(\sqrt{-5})) = 7 $.


title: Left invariant special K\"ahler structures

abstract: We construct left invariant special K\"ahler structures on the cotangent bundle of a flat pseudo-Riemannian Lie group. We introduce the twisted cartesian product of two special K\"ahler Lie algebras according to two linear representations by infinitesimal K\"ahler transformations. We also exhibit a double extension process of a special K\"ahler Lie algebra which allows us to get all simply connected special K\"ahler Lie groups with bi-invariant symplectic connections. All Lie groups constructed by performing this double extension process can be identified with a subgroup of symplectic (or K\"ahler) affine transformations of its Lie algebra containing a nontrivial $1$-parameter subgroup formed by central translations. We show a characterization of left invariant flat special K\"ahler structures using \'etale K\"ahler affine representations, exhibit some immediate consequences of the constructions mentioned above, and give several non-trivial examples.


title: A CMV connection between orthogonal polynomials on the unit circle and the real line

abstract: M. Derevyagin, L. Vinet and A. Zhedanov introduced in Constr. Approx. 36 (2012) 513-535 a new connection between orthogonal polynomials on the unit circle and the real line. It maps any real CMV matrix into a Jacobi one depending on a real parameter $\lambda$. These authors prove that this map yields a natural link between the Jacobi polynomials on the unit circle and the little and big $-1$ Jacobi polynomials on the real line. They also provide explicit expressions for the measure and orthogonal polynomials associated with the Jacobi matrix in terms of those related to the CMV matrix, but only for the value $\lambda=1$ which simplifies the connection -basic DVZ connection-. However, similar explicit expressions for an arbitrary value of $\lambda$ -- (general) DVZ connection -- are missing. This is the main problem overcome in this paper. This work introduces a new approach to the DVZ connection which formulates it as a two-dimensional eigenproblem by using known properties of CMV matrices. This allows us to go further, providing explicit relations between the measures and orthogonal polynomials for the general DVZ connection. It turns out that this connection maps a measure on the unit circle into a rational perturbation of an even measure supported on two symmetric intervals of the real line, which reduce to a single interval for the basic DVZ connection, while the perturbation becomes a degree one polynomial. Some instances of the DVZ connection are shown to give new one-parameter families of orthogonal polynomials on the real line.


title: Polytropic Behavior of Solar Wind Protons Observed by Parker Solar Probe

abstract: A polytropic process describes the transition of a fluid from one state to another through a specific relationship between the fluid density and temperature. The value of the polytropic index that governs this relationship determines the heat transfer and the effective degrees of freedom during the process. In this study, we analyze solar wind proton plasma measurements, obtained by the Faraday cup instrument on-board Parker Solar Probe. We examine the large-scale variations of the proton plasma density and temperature within the inner heliosphere explored by the spacecraft. We also address a polytropic behavior in the density and temperature fluctuations in short-time intervals, which we analyze in order to derive the effective polytropic index of small time-scale processes. The large-scale variations of the solar wind proton density and temperature which are associated with the plasma expansion through the heliosphere, follow a polytropic model with a polytropic index ~5/3. On the other hand, the short time-scale fluctuations which may be associated with turbulence, follow a model with a larger polytropic index. We investigate possible correlations between the polytropic index of short time-scale fluctuations and the plasma speed, plasma beta, and the magnetic field direction. We discuss the scenario of mechanisms including energy transfer or mechanisms that restrict the particle effective degrees of freedom.


title: Weak and Semi-Contraction for Network Systems and Diffusively-Coupled Oscillators

abstract: We develop two generalizations of contraction theory, namely, semi-contraction and weak-contraction theory. First, using the notion of semi-norm, we propose a geometric framework for semi-contraction theory. We introduce matrix semi-measures and characterize their properties. We show that the spectral abscissa of a matrix is the infimum over weighted semi-measures. For dynamical systems, we use the semi-measure of their Jacobian to characterize the contractivity properties of their trajectories. Second, for weakly contracting systems, we prove a dichotomy for the asymptotic behavior of their trajectories and novel sufficient conditions for convergence to an equilibrium. Third, we show that every trajectory of a doubly-contracting system, i.e., a system that is both weakly and semi-contracting, converges to an equilibrium point. Finally, we apply our results to various important network systems including affine averaging and affine flow systems, continuous-time distributed primal-dual algorithms, and networks of diffusively-coupled dynamical systems. For diffusively-coupled systems, the semi-contraction theory leads to a sufficient condition for synchronization that is sharper, in general, than previously-known tests.


title: Random Toeplitz Matrices: The Condition Number under High Stochastic Dependence

abstract: In this paper, we study the condition number of a random Toeplitz matrix. Since a Toeplitz matrix is a diagonal constant matrix, its rows or columns cannot be stochastically independent. This situation does not permit us to use the classic strategy to analyze its minimum singular value when all the entries of a random matrix are stochastically independent. Using a circulant embedding, we can break the stochastic dependence of the structure of the Toeplitz matrix and reduce the problem to analyze the extreme singular values of a random circulant matrix. Among our results, we show the condition number of non--symmetric random circulant matrix of dimension $n$ under the existence of moment generating function of the random entries is $\kappa\left(\mathcal{C}_n\right) = \mbox{O}\left( \frac{1}{\varepsilon} n^{\rho+1/2} \left(\log n\right)^{1/2} \right)$ with probability $1-\mbox{O}\left((\varepsilon^2 + \varepsilon) n^{-2\rho} + n^{-1/2+\scriptstyle{o}(1)}\right)$ for any $\varepsilon >0$, $\rho\in(0,1/4)$. Moreover, if the random entries only have the second moment, we have $\kappa\left(\mathcal{C}_n\right) = \mbox{O}\left( \frac{1}{\varepsilon} n^{\rho+1/2} \log n\right)$ with probability $1-\mbox{O}\left((\varepsilon^2 + \varepsilon) n^{-2\rho} + \left(\log n\right)^{-1/2}\right)$. For the condition number of a random (non--symmetric or symmetric) Toeplitz matrix $\mathcal{T}_n$ we establish $\kappa\left(\mathcal{T}_n\right) \leq \kappa\left(\mathcal{C}_{2n}\right) \left(\sigma_{\min}\left( C_{2n} \right)\sigma_{\min}\left( S_n \right)\right)^{-1}$, where $\sigma_{\min}(A)$ is the minimum singular value of the matrix $A$. The matrix $C_{2n}$ is a random circulant matrix and $S_n:=F^*_{2,n} D_{1,n}^{-1}F_{2,n} + F^*_{4,n} D^{-1}_2 F_{4,n}$, where $F_{2,n},F_{4,n}$ are deterministic matrices and $D_{1,n}, D_{2,n}$ are random diagonal matrices. We conjeture $S_n$ is well conditioned.


title: Exploiting machine learning to efficiently predict multidimensional optical spectra in complex environments

abstract: The excited state dynamics of chromophores in complex environments determine a range of vital biological and energy capture processes. Time-resolved, multidimensional optical spectroscopies provide a key tool to investigate these processes. Although theory has the potential to decode these spectra in terms of the electronic and atomistic dynamics, the need for large numbers of excited state electronic structure calculations severely limits first principles predictions of multidimensional optical spectra for chromophores in the condensed phase. Here, we leverage the locality of chromophore excitations to develop machine learning models to predict the excited state energy gap of chromophores in complex environments for efficiently constructing linear and multidimensional optical spectra. By analyzing the performance of these models, which span a hierarchy of physical approximations, across a range of chromophore-environment interaction strengths, we provide strategies for the construction of ML models that greatly accelerate the calculation of multidimensional optical spectra from first principles.


title: Polarity dependent heating at the phase interface in metal-insulator transitions

abstract: Current-driven insulator-metal transitions are in many cases driven by Joule heating proportional to the square of the applied current. Recent nano-imaging experiments in Ca$_2$RuO$_4$ reveal a metal-insulator phase boundary that depends on the direction of an applied current, suggesting an important non-heating effect. Motivated by these results, we study the effects of an electric current in a system containing interfaces between metallic and insulating phases. Derivation of a heat balance equation from general macroscopic Onsager transport theory, reveals a heating term proportional to the product of the current across the interface and the discontinuity in the Seebeck coefficient, so that heat can either be generated or removed at an interface, depending on the direction of the current relative to the change in material properties. For parameters appropriate to Ca$_2$RuO$_4$, this heating can be comparable to or larger than Joule heating. A simplified model of the relevant experimental geometry is shown to provide results consistent with the experiments. Extension of the results to other inhomogeneous metal-insulator transition systems is discussed.


title: Ultracold collisions of the lithium monoxide radical

abstract: Ultracold collisions of LiO molecules in the $^{2}\Pi_{3/2}$ ground state are considered, under the influence of either an external magnetic or electric field. Inelastic collisions are shown to be suppressed in the presence of modest laboratory strength magnetic and electric fields. The rate of elastic collisions that rethermalize the thermal distribution, and the corresponding low rate of heating state-changing collisions, suggest that quantum degeneracy or even molecular Bose-Einstein condensation of LiO gas may be attainable, provided that the initial temperatures in the milliKelvin range are achievable.


title: The Beurling-Wintner problem for characteristic functions

abstract: This paper concerns a long-standing problem raised by Beurling and Wintner on completeness of the dilation system $\{\varphi(kx):k=1,2,\cdots\}$ generated by the odd periodic extension on $\mathbb{R}$ of any $\varphi\in L^2[0,1]$. Up to now there has been no explicit description of solutions of the Beurling-Wintner problem even for characteristic functions. We focus on characteristic function $\mathbf{1}_V$ of an open subset $V$ of $(0,1)$ where $V$ is the union of finitely many intervals with rational endpoints. Using substantially techniques from analytic number theory, we fully solved the Beurling-Wintner problem in most interesting situations and exhibit the explicit form of such $V$. As a consequence, it yields a complete solution for the rational version of Kozlov's problem. Moreover, we find that the Beurling-Wintner problem is closely related to the Twin Prime Conjecture and the Sophie Germain Prime Conjecture.


title: The Impact of Unmeasured Within- and Between-Cluster Confounding on the Bias of Effect Estimators from Fixed Effect, Mixed effect and Instrumental Variable Models

abstract: Instrumental variable methods are popular choices in combating unmeasured confounding to obtain less biased effect estimates. However, we demonstrate that alternative methods may give less biased estimates depending on the nature of unmeasured confounding. Treatment preferences of clusters (e.g., physician practices) are the most f6requently used instruments in instrumental variable analyses (IVA). These preference-based IVAs are usually conducted on data clustered by region, hospital/facility, or physician, where unmeasured confounding often occurs within or between clusters. We aim to quantify the impact of unmeasured confounding on the bias of effect estimators in IVA, as well as alternative methods including ordinary least squares regression, linear mixed models (LMM) and fixed effect models (FE) to study the effect of a continuous exposure (e.g., treatment dose). We derive bias formulae of estimators from these four methods in the presence of unmeasured within- and/or between-cluster confounders. We show that IVAs can provide consistent estimates when unmeasured within-cluster confounding exists, but not when between-cluster confounding exists. On the other hand, FEs and LMMs can provide consistent estimates when unmeasured between-cluster confounding exits, but not for within-cluster confounding. Whether IVAs are advantageous in reducing bias over FEs and LMMs depends on the extent of unmeasured within-cluster confounding relative to between-cluster confounding. Furthermore, the impact of unmeasured between-cluster confounding on IVA estimates is larger than the impact of unmeasured within-cluster confounding on FE and LMM estimates. We illustrate these methods through data applications. Our findings provide guidance for choosing appropriate methods to combat the dominant types of unmeasured confounders and help interpret statistical results in the context of unmeasured confounding


title: SwarmControl: An Automated Distributed Control Framework for Self-Optimizing Drone Networks

abstract: Networks of Unmanned Aerial Vehicles (UAVs), composed of hundreds, possibly thousands of highly mobile and wirelessly connected flying drones will play a vital role in future Internet of Things (IoT) and 5G networks. However, how to control UAV networks in an automated and scalable fashion in distributed, interference-prone, and potentially adversarial environments is still an open research problem. This article introduces SwarmControl, a new software-defined control framework for UAV wireless networks based on distributed optimization principles. In essence, SwarmControl provides the Network Operator (NO) with a unified centralized abstraction of the networking and flight control functionalities. High-level control directives are then automatically decomposed and converted into distributed network control actions that are executed through programmable software-radio protocol stacks. SwarmControl (i) constructs a network control problem representation of the directives of the NO; (ii) decomposes it into a set of distributed sub-problems; and (iii) automatically generates numerical solution algorithms to be executed at individual UAVs. We present a prototype of an SDR-based, fully reconfigurable UAV network platform that implements the proposed control framework, based on which we assess the effectiveness and flexibility of SwarmControl with extensive flight experiments. Results indicate that the SwarmControl framework enables swift reconfiguration of the network control functionalities, and it can achieve an average throughput gain of 159% compared to the state-of-the-art solutions.


title: On Lipschitz-free spaces over spheres of Banach spaces

abstract: We prove that, for each Banach space $X$ which is isomorphic to its hyperplanes, the Lipschitz-free spaces over $X$ and over its sphere are isomorphic.


title: Classification of full exceptional collections on smooth toric Fano varieties with Picard rank two

abstract: In this paper, we give the complete classification of full exceptional collections on smooth toric Fano threefolds and fourfolds with Picard rank two. To be precise, we give a partial answer to the conjecture in \cite{Kuz} and \cite{LYY}: we provide all the exceptional collections of maximal length and prove that they are in fact full exceptional collections.


title: Images and Misinformation in Political Groups: Evidence from WhatsApp in India

abstract: WhatsApp is a key medium for the spread of news and rumors, often shared as images. We study a large collection of politically-oriented WhatsApp groups in India, focusing on the period leading up to the 2019 Indian national elections. By labeling samples of random and popular images, we find that around 13% of shared images are known misinformation and most fall into three types of images. Machine learning methods can be used to predict whether a viral image is misinformation, but are brittle to shifts in content over time.


title: Approximation properties in Lipschitz-free spaces over groups

abstract: We study Lipschitz-free spaces over compact and uniformly discrete metric spaces enjoying certain high regularity properties - having group structure with left-invariant metric. Using methods of harmonic analysis we show that, given a compact metrizable group $G$ equipped with an arbitrary compatible left-invariant metric $d$, the Lipschitz-free space over $G$, $\mathcal{F}(G,d)$, satisfies the metric approximation property. We show also that, given a finitely generated group $G$, with its word metric $d$, from a class of groups admitting a certain special type of combing, which includes all hyperbolic groups and Artin groups of large type, $\mathcal{F}(G,d)$ has a Schauder basis. Examples and applications are discussed. In particular, for any net $N$ in a real hyperbolic $n$-space $\mathbb{H}^n$, $\mathcal{F}(N)$ has a Schauder basis.


title: Optimal quarantine strategies for the COVID-19 pandemic in a population with a discrete age structure

abstract: The goal of this work is to study the optimal controls for the COVID-19 epidemic in Brazil. We consider an age-structured SEIRQ model with quarantine compartment, where the controls are the quarantine entrance parameters. We then compare the optimal controls for different quarantine lengths and distribution of the total control cost by assessing their respective reductions in deaths in comparison to the same period without quarantine. The best strategy provides a calendar of when to relax the isolation measures for each age group. Finally, we analyse how a delay in the beginning of the quarantine affects this calendar by changing the initial conditions.


abstract: Many environments currently employ machine learning models for data processing and analytics that were built using a limited number of training data points. Once deployed, the models are exposed to significant amounts of previously-unseen data, not all of which is representative of the original, limited training data. However, updating these deployed models can be difficult due to logistical, bandwidth, time, hardware, and/or data sensitivity constraints. We propose a framework, Self-Updating Models with Error Remediation (SUMER), in which a deployed model updates itself as new data becomes available. SUMER uses techniques from semi-supervised learning and noise remediation to iteratively retrain a deployed model using intelligently-chosen predictions from the model as the labels for new training iterations. A key component of SUMER is the notion of error remediation as self-labeled data can be susceptible to the propagation of errors. We investigate the use of SUMER across various data sets and iterations. We find that self-updating models (SUMs) generally perform better than models that do not attempt to self-update when presented with additional previously-unseen data. This performance gap is accentuated in cases where there is only limited amounts of initial training data. We also find that the performance of SUMER is generally better than the performance of SUMs, demonstrating a benefit in applying error remediation. Consequently, SUMER can autonomously enhance the operational capabilities of existing data processing systems by intelligently updating models in dynamic environments.


title: Non-destructive study of collision cascade damage in self-ion irradiated tungsten using HR-EBSD and ECCI

abstract: Understanding defect production and evolution under irradiation is a long-standing multi-scale problem. Conventionally, experimental examination of irradiation-induced defects (IIDs) has mainly relied on transmission electron microscopy (TEM), which offers high spatial resolution but requires destructive sample preparation. Furthermore, limited field of view and low strain sensitivity make multi-scale characterisation and quantitative strain measurements difficult. Here we explore the potential of using advanced techniques in the scanning electron microscope (SEM) to non-destructively probe irradiation damage at the surface of bulk materials. Electron channelling contrast imaging (ECCI) is used to image nano-scale irradiation-induced defects in 20 MeV self-ion irradiated tungsten, the main candidate material for fusion reactor armour. The results show an evolution of the damage microstructure from uniformly and randomly distributed nano-scale defects at 0.01 dpa (displacement per atom) to string structures extending over hundreds of nanometres at 1 dpa. Cross-correlation based high-resolution EBSD (HR-EBSD) is used to probe the lattice strain fields associated with IIDs. While there is little strain fluctuation at 0.01 dpa, significant heterogeneity in the lattice strains is observed at 0.1 dpa, increasing with dose until saturation at 0.32 dpa. The characteristic length scale of strain fluctuations is ~500 nm. Together, ECCI and HR-EBSD reveal a transition from a structure where defects are disordered to a structure with long-range order driven by elastic interactions between pre-existing defects and new cascade damage. This study demonstrates that SEM provides an attractive tool for rapid throughput, non-destructive, multi-scale and multi-aspect characterisation of irradiation damage.


title: Adjunction in 2-categories

abstract: We give an overview of the parts of arXiv:2004.04279 that deal with 2-categories, up to and including adjunction, and explain how the Segal-type approach to 2-categories adopted there is related to the more standard approaches. As an application, we construct the derived Morita and the Fourier-Mukai 2-categories over a Noetherian ring, and show how to embed one into the other.


abstract: The Layer 2 Atomic Cross-Blockchain Function Calls protocol allows composable programming across Ethereum blockchains. It allows for inter-contract and inter-blockchain function calls that are both synchronous and atomic: if one part fails, the whole call graph of function calls is rolled back. Existing atomic cross-blockchain function call protocols are Blockchain Layer 1 protocols, which require changes to the blockchain platform software to operate. Blockchain Layer 2 technologies such as the one described in this paper require no such changes. They operate on top of the infrastructure provided by the blockchain platform software. This paper introduces the protocol and a more scalable variant, provides an initial safety and liveness analysis, and presents the expected overhead of using this technology when compared to using multiple non-atomic single blockchain transactions. The overhead is analysed for three scenarios involving multiple blockchains: the Hotel and Train problem, Supply Chain with Provenance, and an Oracle. The protocol is shown to provide 93.8 or 186 cross-blockchain function calls per second for the Hotel and Train scenario when there are many travel agencies, for the standard and scalable variant of the protocol respectively, given the Ethereum client, Hyperledger Besu's performance of 375 tps, assuming a block period of one second, and assuming all transactions take the same amount of time to execute as the benchmark transactions.


title: Hydrostatic and uniaxial pressure tuning of iron-based superconductors: Insights into superconductivity, magnetism, nematicity and collapsed tetragonal transitions

abstract: Iron-based superconductors are well-known for their intriguing phase diagrams, which manifest a complex interplay of electronic, magnetic and structural degrees of freedom. Among the phase transitions observed are superconducting, magnetic, and several types of structural transitions, including a tetragonal-to-orthorhombic and a collapsed-tetragonal transition. In particular, the widely-observed tetragonal-to-orthorhombic transition is believed to be a result of an electronic order that is coupled to the crystalline lattice and is, thus, referred to as nematic transition. Nematicity is therefore a prominent feature of these materials, which signals the importance of the coupling of electronic and lattice properties. Correspondingly, these systems are particularly susceptible to tuning via pressure (hydrostatic, uniaxial, or some combination). We review efforts to probe the phase diagrams of pressure-tuned iron-based superconductors, with a strong focus on our own recent insights into the phase diagrams of several members of this material class under hydrostatic pressure. These studies on FeSe, Ba(Fe$_{1-x}$Co$_x$)$_2$As$_2$, Ca(Fe$_{1-x}$Co$_x$)$_2$As$_2$ and CaK(Fe$_{1-x}$Ni$_x$)$_4$As$_4$ were, to a significant extent, made possible by advances of what measurements can be adapted to the use under differing pressure environments. We point out the potential impact of these tools for the study of the wider class of strongly correlated electron systems.


title: Lie algebra structure of fitness and replicator control

abstract: For over fifty years, the dynamical systems perspective has had a prominent role in evolutionary biology and economics, through the lens of game theory. In particular, the study of replicator differential equations on the standard (probability) simplex, specified by fitness maps or payoff functions, has yielded insights into the temporal behavior of such systems. However behavior is influenced by context and environmental factors with a game-changing quality (i.e., fitness maps are manipulated). This paper develops a principled geometric approach to model and understand such influences by incorporating replicator dynamics into a broader control-theoretic framework. Central to our approach is the construction of a Lie algebra structure on the space of fitness maps, mapping homomorphically to the Lie algebra of replicator vector fields. This is akin to classical mechanics, where the Poisson bracket Lie algebra of functions maps to associated Hamiltonian vector fields. We show, extending the work of Svirezhev in 1972, that a trajectory of a replicator vector field is the base integral curve of a solution to a Hamiltonian system defined on the cotangent bundle of the simplex. Further, we exploit the Lie algebraic structure of fitness maps to determine controllability properties of a class of replicator systems.


title: Quark-Novae in the outskirts of galaxies: An explanation of the Fast Radio Burst phenomenon

abstract: We show that old isolated neutron stars in groups and clusters of galaxies experiencing a Quark-Nova phase (QN: an explosive transition to a quark star) may be the source of FRBs. Each of the millions of fragments of the ultra-relativistic QN ejecta provides a collisionless plasma for which the ambient medium (galactic/halo, the intra-group/intra-cluster medium) acts as a relativistic plasma beam. The Buneman and the Weibel instabilities, successively induced by the beam in the fragment, generate particle bunching and observed coherent emission at GHz frequency with a corresponding fluence in the Jy ms range. The duration, frequency drift and the rate are in agreement with observed properties of FRBs. Repeats (on timescales of minutes to months) are due to seeing multiple fragments each beaming at a different direction and coming in at different times. Single (non-repeating) FRBs, occur when only emission from the primary fragment is within the detector's sensitivity. Key properties of FRB 121102 (its years of activity) and of FRB 180916.J0158+65 (its 16 day period) are recovered. The spatial and temporal coincidence between SGR 1935+2154 and FRB 200428 finds an explanation in our model. We give testable predictions.


title: Pairs Trading with Nonlinear and Non-Gaussian State Space Models

abstract: This paper studies pairs trading using a nonlinear and non-Gaussian state-space model framework. We model the spread between the prices of two assets as an unobservable state variable and assume that it follows a mean-reverting process. This new model has two distinctive features: (1) The innovations to the spread is non-Gaussianity and heteroskedastic. (2) The mean reversion of the spread is nonlinear. We show how to use the filtered spread as the trading indicator to carry out statistical arbitrage. We also propose a new trading strategy and present a Monte Carlo based approach to select the optimal trading rule. As the first empirical application, we apply the new model and the new trading strategy to two examples: PEP vs KO and EWT vs EWH. The results show that the new approach can achieve a 21.86% annualized return for the PEP/KO pair and a 31.84% annualized return for the EWT/EWH pair. As the second empirical application, we consider all the possible pairs among the largest and the smallest five US banks listed on the NYSE. For these pairs, we compare the performance of the proposed approach with that of the existing popular approaches, both in-sample and out-of-sample. Interestingly, we find that our approach can significantly improve the return and the Sharpe ratio in almost all the cases considered.


title: Chirp Mass and Spin of Binary Black Holes from First Star Remnants

abstract: We performed Population III (Pop III) binary evolution by using population synthesis simulations for seven different models. We found that Pop III binaries tend to be binary black holes (BBHs) with chirp mass $M_{\rm chirp} \sim 30~M_{\odot}$ and they can merge at present day due to long merger time. The merger rate densities of Pop III BBHs at $z=0$ ranges 3.34--21.2 $\rm /yr/Gpc^3$ which is consistent with the aLIGO/aVIRGO result of 9.7--101 $\rm /yr/Gpc^3$. These Pop III binaries might contribute to some part of the massive BBH gravitational wave (GW) sources detected by aLIGO/aVIRGO. We also calculated the redshift dependence of Pop III BBH mergers. We found that Pop III low spin BBHs tend to merge at low redshift, while Pop III high spin BBHs do at high redshift, which can be confirmed by future GW detectors such as ET, CE, and DECIGO. These detectors can also check the redshift dependence of BBH merger rate and spin distribution. Our results show that except for one model, the mean effective spin $\left\langle \chi_{\rm eff} \right\rangle$ at $z=0$ ranges $0.02$--$0.3$ while at $z=10$ it does $0.16$--$0.64$. Therefore, massive stellar-mass BBH detection by GWs will be a key for the stellar evolution study in the early universe.


title: List Decodable Mean Estimation in Nearly Linear Time

abstract: Learning from data in the presence of outliers is a fundamental problem in statistics. Until recently, no computationally efficient algorithms were known to compute the mean of a high dimensional distribution under natural assumptions in the presence of even a small fraction of outliers. In this paper, we consider robust statistics in the presence of overwhelming outliers where the majority of the dataset is introduced adversarially. With only an $\alpha < 1/2$ fraction of "inliers" (clean data) the mean of a distribution is unidentifiable. However, in their influential work, [CSV17] introduces a polynomial time algorithm recovering the mean of distributions with bounded covariance by outputting a succinct list of $O(1/\alpha)$ candidate solutions, one of which is guaranteed to be close to the true distributional mean; a direct analog of 'List Decoding' in the theory of error correcting codes. In this work, we develop an algorithm for list decodable mean estimation in the same setting achieving up to constants the information theoretically optimal recovery, optimal sample complexity, and in nearly linear time up to polylogarithmic factors in dimension. Our conceptual innovation is to design a descent style algorithm on a nonconvex landscape, iteratively removing minima to generate a succinct list of solutions. Our runtime bottleneck is a saddle-point optimization for which we design custom primal dual solvers for generalized packing and covering SDP's under Ky-Fan norms, which may be of independent interest.


title: B\'ezoutians and injectivity of polynomial maps

abstract: We prove that an endomorphism $f$ of affine space is injective on rational points if its B\'ezoutian is constant. Similarly, $f$ is injective at a given rational point if its reduced B\'ezoutian is constant. We also show that if the Jacobian determinant of $f$ is invertible, then $f$ is injective at a given rational point if and only if its reduced B\'ezoutian is constant.


title: Groups for which it is easy to detect graphical regular representations

abstract: We say that a finite group G is "DRR-detecting" if, for every subset S of G, either the Cayley digraph Cay(G,S) is a digraphical regular representation (that is, its automorphism group acts regularly on its vertex set) or there is a nontrivial group automorphism phi of G such that phi(S) = S. We show that every nilpotent DRR-detecting group is a p-group, but that the wreath product of two cyclic groups of order p is not DRR-detecting, for every odd prime p. We also show that if G and H are nontrivial groups that admit a digraphical regular representation and either gcd(|G|,|H|) = 1, or H is not DRR-detecting, then the direct product G x H is not DRR-detecting. Some of these results also have analogues for graphical regular representations.


title: Dimension formula of the affine Deligne-Lusztig variety $X(\mu, b)$

abstract: The study of certain union $X(\mu, b)$ of affine Deligne-Lusztig varieties in the affine flag varieties arose from the study of Shimura varieties with Iwahori level structure. In this paper, we give an explicit formula of $\dim X(\mu, b)$ for sufficiently large dominant coweight $\mu$.


title: Fingerprinting Encrypted Voice Traffic on Smart Speakers with Deep Learning

abstract: This paper investigates the privacy leakage of smart speakers under an encrypted traffic analysis attack, referred to as voice command fingerprinting. In this attack, an adversary can eavesdrop both outgoing and incoming encrypted voice traffic of a smart speaker, and infers which voice command a user says over encrypted traffic. We first built an automatic voice traffic collection tool and collected two large-scale datasets on two smart speakers, Amazon Echo and Google Home. Then, we implemented proof-of-concept attacks by leveraging deep learning. Our experimental results over the two datasets indicate disturbing privacy concerns. Specifically, compared to 1% accuracy with random guess, our attacks can correctly infer voice commands over encrypted traffic with 92.89\% accuracy on Amazon Echo. Despite variances that human voices may cause on outgoing traffic, our proof-of-concept attacks remain effective even only leveraging incoming traffic (i.e., the traffic from the server). This is because the AI-based voice services running on the server side response commands in the same voice and with a deterministic or predictable manner in text, which leaves distinguishable pattern over encrypted traffic. We also built a proof-of-concept defense to obfuscate encrypted traffic. Our results show that the defense can effectively mitigate attack accuracy on Amazon Echo to 32.18%.


title: FashionBERT: Text and Image Matching with Adaptive Loss for Cross-modal Retrieval

abstract: In this paper, we address the text and image matching in cross-modal retrieval of the fashion industry. Different from the matching in the general domain, the fashion matching is required to pay much more attention to the fine-grained information in the fashion images and texts. Pioneer approaches detect the region of interests (i.e., RoIs) from images and use the RoI embeddings as image representations. In general, RoIs tend to represent the "object-level" information in the fashion images, while fashion texts are prone to describe more detailed information, e.g. styles, attributes. RoIs are thus not fine-grained enough for fashion text and image matching. To this end, we propose FashionBERT, which leverages patches as image features. With the pre-trained BERT model as the backbone network, FashionBERT learns high level representations of texts and images. Meanwhile, we propose an adaptive loss to trade off multitask learning in the FashionBERT modeling. Two tasks (i.e., text and image matching and cross-modal retrieval) are incorporated to evaluate FashionBERT. On the public dataset, experiments demonstrate FashionBERT achieves significant improvements in performances than the baseline and state-of-the-art approaches. In practice, FashionBERT is applied in a concrete cross-modal retrieval application. We provide the detailed matching performance and inference efficiency analysis.


title: A central limit theorem for descents of a Mallows permutation and its inverse

abstract: This paper studies the asymptotic distribution of descents $\des(w)$ in a permutation $w$, and its inverse, distributed according to the Mallows measure. The Mallows measure is a non-uniform probability measure on permutations introduced to study ranked data. Under this measure, permutations are weighted according to the number of inversions they contain, with the weighting controlled by a parameter $q$. The main results are a Berry-Esseen theorem for $\des(w)+\des(w^{-1})$ as well as a joint central limit theorem for $(\des(w),\des(w^{-1}))$ to a bivariate normal with a non-trivial correlation depending on $q$. The proof uses Stein's method with size-bias coupling along with a regenerative process associated to the Mallows measure.


title: A Computational Analysis of Polarization on Indian and Pakistani Social Media

abstract: Between February 14, 2019 and March 4, 2019, a terrorist attack in Pulwama, Kashmir followed by retaliatory airstrikes led to rising tensions between India and Pakistan, two nuclear-armed countries. In this work, we examine polarizing messaging on Twitter during these events, particularly focusing on the positions of Indian and Pakistani politicians. We use a label propagation technique focused on hashtag co-occurrences to find polarizing tweets and users. Our analysis reveals that politicians in the ruling political party in India (BJP) used polarized hashtags and called for escalation of conflict more so than politicians from other parties. Our work offers the first analysis of how escalating tensions between India and Pakistan manifest on Twitter and provides a framework for studying polarizing messages.


title: Dessins d'enfants and some holomorphic structures on the Loch Ness Monster

abstract: The classical theory of dessin d'enfants, which are bipartite maps on compact orientable surfaces, are combinatorial objects used to study branched covers between compact Riemann surfaces and the absolute Galois group of the field of rational numbers. In this paper, we show how this theory is naturally extended to non-compact orientable surfaces and, in particular, we observe that the Loch Ness monster (the surface of infinite genus with exactly one end) admits infinitely many regular dessins d'enfants (either chiral or reflexive). In addition, we study different holomorphic structures on the Loch Ness monster, which come from homology covers of compact Riemann surfaces, infinite hyperelliptic and infinite superelliptic curves.


title: Molecular Characterization of Ebselen Binding Activity to SARS-CoV-2 Main Protease

abstract: The Coronavirus Disease (COVID-19) pandemic caused by the SARS-coronavirus 2 (SARS-CoV-2) urgently calls for the design of drugs directed against this new virus. Given its essential role in proteolytic processing, the main protease Mpro has been identified as an attractive candidate for drugs against SARS-CoV-2 and similar coronaviruses. Recent high-throughput screening studies have identified a set of existing, small-molecule drugs as potent Mpro inhibitors. Amongst these, Ebselen (2-Phenyl-1,2-benzoselenazol-3-one), a glutathione peroxidase mimetic seleno-organic compound, is particularly attractive. Recent experiments suggest that its effectiveness is higher than that of other molecules that also act at the enzyme catalytic site. By relying on extensive simulations with all-atom models, in this study we examine at a molecular level the potential of Ebselen to decrease Mpro catalytic activity. Our results indicate that Ebselen exhibits a distinct affinity for the catalytic site cavity of Mpro. In addition, our molecular models reveal a second, previously unkown binding site for Ebselen in the dimerization region localized between the II and III domains of the protein. A detailed analysis of the free energy of binding indicates that the affinity of Ebselen to this second binding site is in fact significantly larger than that to the catalytic site. A strain analysis indicates that Ebselen bound between the II-III domains exerts a pronounced allosteric effect that regulates catalytic site access through surface loop interactions, and induces a displacement and reconfiguration of water hotspots, including the catalytic water, which could interfere with normal enzymatic function. Taken together, these findings provide a framework for the future design of more potent and specific Mpro inhibitors, based on the Ebselen scaffold, that could lead to new therapeutic strategies for COVID-19.


title: Coverage Analysis of Net Inscriptions in Coloured Petri Net Models

abstract: High-level Petri net such as Coloured Petri Nets (CPNs) are characterised by the combination of Petri nets and a high-level programming language. In the context of CPNs and CPN Tools, the inscriptions (e.g., arc expressions and guards) are specified using Standard ML (SML). The application of simulation and state space exploration (SSE) for validating CPN models traditionally focuses on behavioural properties related to net structure, i.e., places and transitions. This means that the net inscriptions are only implicitly validated, and the extent to which these have been covered is not made explicit. The contribution of this paper is an approach that establishes a link between coverage analysis known from programming languages and net inscriptions of CPN models. Specifically, we consider Modified Condition/Decision Coverage (MC/DC) which generalises branch coverage of SML decisions. We have implemented our approach in a library for CPN Tools comprised of an annotation and instrumentation mechanism that transparently intercept and collect evaluations of Boolean conditions, and a post-processing tool that determines whether each decision is MC/DC-covered by a set of models executions (runs). We evaluate our approach on four larger public-available CPN models.


title: Neural Ordinary Differential Equation based Recurrent Neural Network Model

abstract: Neural differential equations are a promising new member in the neural network family. They show the potential of differential equations for time series data analysis. In this paper, the strength of the ordinary differential equation (ODE) is explored with a new extension. The main goal of this work is to answer the following questions: (i)~can ODE be used to redefine the existing neural network model? (ii)~can Neural ODEs solve the irregular sampling rate challenge of existing neural network models for a continuous time series, i.e., length and dynamic nature, (iii)~how to reduce the training and evaluation time of existing Neural ODE systems? This work leverages the mathematical foundation of ODEs to redesign traditional RNNs such as Long Short-Term Memory (LSTM) and Gated Recurrent Unit (GRU). The main contribution of this paper is to illustrate the design of two new ODE-based RNN models (GRU-ODE model and LSTM-ODE) which can compute the hidden state and cell state at any point of time using an ODE solver. These models reduce the computation overhead of hidden state and cell state by a vast amount. The performance evaluation of these two new models for learning continuous time series with irregular sampling rate is then demonstrated. Experiments show that these new ODE based RNN models require less training time than Latent ODEs and conventional Neural ODEs. They can achieve higher accuracy quickly, and the design of the neural network is simpler than, previous neural ODE systems.


title: Bernoulli generalized likelihood ratio test for signal detection from photon counting images

abstract: Because exoplanets are extremely dim, an Electron Multiplying Charged Coupled Device (EMCCD) operating in photon counting (PC) mode is necessary to reduce the detector noise level and enable their detection. Typically, PC images are added together as a co-added image before processing. We present here a signal detection and estimation technique that works directly with individual PC images. The method is based on the generalized likelihood ratio test (GLRT) and uses a Bernoulli distribution between PC images. The Bernoulli distribution is derived from a stochastic model for the detector, which accurately represents its noise characteristics. We show that our technique outperforms a previously used GLRT method that relies on co-added images under a Gaussian noise assumption and two detection algorithms based on signal-to-noise ratio (SNR). Furthermore, our method provides the maximum likelihood estimate of exoplanet intensity and background intensity while doing detection. It can be applied online, so it is possible to stop observations once a specified threshold is reached, providing confidence for the existence (or absence) of planets. As a result, the observation time is efficiently used. Besides the observation time, the analysis of detection performance introduced in the paper also gives quantitative guidance on the choice of imaging parameters, such as the threshold. Lastly, though this work focuses on the example of detecting point source, the framework is widely applicable.


title: A Semicircle Law for Derivatives of Random Polynomials

abstract: Let $x_1, \dots, x_n$ be $n$ independent and identically distributed random variables with mean zero, unit variance, and finite moments of all remaining orders. We study the random polynomial $p_n$ having roots at $x_1, \dots, x_n$. We prove that for $\ell \in \mathbb{N}$ fixed as $n \rightarrow \infty$, the $(n-\ell)-$th derivative of $p_n^{}$ behaves like a Hermite polynomial: for $x$ in a compact interval,$${n^{\ell/2}} \frac{\ell!}{n!} \cdot p_n^{(n-\ell)}\left( \frac{x}{\sqrt{n}}\right) \rightarrow He_{\ell}(x + \gamma_n),$$ where $He_{\ell}$ is the $\ell-$th probabilists' Hermite polynomial and $\gamma_n$ is a random variable converging to the standard $\mathcal{N}(0,1)$ Gaussian as $n \rightarrow \infty$. Thus, there is a universality phenomenon when differentiating a random polynomial many times: the remaining roots follow a Wigner semicircle distribution.


title: $p$-Norm Flow Diffusion for Local Graph Clustering

abstract: Local graph clustering and the closely related seed set expansion problem are primitives on graphs that are central to a wide range of analytic and learning tasks such as local clustering, community detection, nodes ranking and feature inference. Prior work on local graph clustering mostly falls into two categories with numerical and combinatorial roots respectively. In this work, we draw inspiration from both fields and propose a family of convex optimization formulations based on the idea of diffusion with p-norm network flow for $p\in (1,\infty)$. In the context of local clustering, we characterize the optimal solutions for these optimization problems and show their usefulness in finding low conductance cuts around input seed set. In particular, we achieve quadratic approximation of conductance in the case of $p=2$ similar to the Cheeger-type bounds of spectral methods, constant factor approximation when $p\rightarrow\infty$ similar to max-flow based methods, and a smooth transition for general $p$ values in between. Thus, our optimization formulation can be viewed as bridging the numerical and combinatorial approaches, and we can achieve the best of both worlds in terms of speed and noise robustness. We show that the proposed problem can be solved in strongly local running time for $p\ge 2$ and conduct empirical evaluations on both synthetic and real-world graphs to illustrate our approach compares favorably with existing methods.


title: A magneto-optical trap with millimeter ball lenses

abstract: We present a magneto-optical trap (MOT) design based on millimeter ball lenses, contained within a metal cube of 0.75$^{\prime \prime}$ side length. We present evidence of trapping approximately $4.2\times 10^5$ of $^{85}$Rb atoms with a number density of $3.2\times 10^9$ atoms/cm$^{3}$ and a loading time of 1.3 s. Measurement and a kinetic laser-cooling model are used to characterize the atom trap design. The design provides several advantages over other types of MOTs: the laser power requirement is low, the small lens and cube sizes allow for miniaturization of MOT applications, and the lack of large-diameter optical beam pathways prevents external blackbody radiation from entering the trapping region.


title: Active Speakers in Context

abstract: Current methods for active speak er detection focus on modeling short-term audiovisual information from a single speaker. Although this strategy can be enough for addressing single-speaker scenarios, it prevents accurate detection when the task is to identify who of many candidate speakers are talking. This paper introduces the Active Speaker Context, a novel representation that models relationships between multiple speakers over long time horizons. Our Active Speaker Context is designed to learn pairwise and temporal relations from an structured ensemble of audio-visual observations. Our experiments show that a structured feature ensemble already benefits the active speaker detection performance. Moreover, we find that the proposed Active Speaker Context improves the state-of-the-art on the AVA-ActiveSpeaker dataset achieving a mAP of 87.1%. We present ablation studies that verify that this result is a direct consequence of our long-term multi-speaker analysis.


title: Polygon Coherent Modes in aWeakly Perturbed Whispering Gallery Microresonator for Efficient Second Harmonic, Optomechanic and Frequency Comb Generations

abstract: We observe high optical quality factor (Q) polygonal and star coherent optical modes in a lithium niobate microdisk. In contrast to the previous polygon modes achieved by deformed microcavities at lower mechanical and optical Q, we adopted weak perturbation from a tapered fiber for the polygon mode formation. The resulting high intracavity optical power of the polygon modes triggers second harmonic generation at high efficiency. With the combined advantage of high mechanical Q cavity optomechanical oscillation was observed for the first time. Finally, we observe frequency microcomb generation from the polygon modes with an ultra stable taper-on-disk coupling mechanism.