title: Urban Analytics: History, Trajectory, and Critique

abstract: Urban analytics combines spatial analysis, statistics, computer science, and urban planning to understand and shape city futures. While it promises better policymaking insights, concerns exist around its epistemological scope and impacts on privacy, ethics, and social control. This chapter reflects on the history and trajectory of urban analytics as a scholarly and professional discipline. In particular, it considers the direction in which this field is going and whether it improves our collective and individual welfare. It first introduces early theories, models, and deductive methods from which the field originated before shifting toward induction. It then explores urban network analytics that enrich traditional representations of spatial interaction and structure. Next it discusses urban applications of spatiotemporal big data and machine learning. Finally, it argues that privacy and ethical concerns are too often ignored as ubiquitous monitoring and analytics can empower social repression. It concludes with a call for a more critical urban analytics that recognizes its epistemological limits, emphasizes human dignity, and learns from and supports marginalized communities.


title: Modeling of quantum dot based CNOT and Toffoli gates in a noisy environment

abstract: Quantum dot-based system is one of the most promising systems owing to its integrability with classical computation hardware and its versatility in realizing qubits and quantum gates. In this work, we investigate the functionality of the two-qubit CNOT gate and three-qubit Toffoli gate in the presence of hyperfine fluctuation noise and phononic noise. We model two and three-qubit gates using the Lindblad Master Equation to estimate the operating range of the external static magnetic field. In these ranges, we observe a successful gate operation under decoherence.


title: On single crystal total scattering data reduction and correction protocols for analysis in direct space

abstract: We explore data reduction and correction steps and processed data reproducibility in the emerging single crystal total scattering based technique of three-dimensional differential atomic pair distribution function (3D-$\Delta$PDF) analysis. All steps from sample measurement to data-processing are outlined in detail using a CuIr$_2$S$_4$ example crystal studied in a setup equipped with a high-energy x-ray beam and a flat panel area detector. Computational overhead as it pertains to data-sampling and the associated data processing steps is also discussed. Various aspects of the final 3D-$\Delta$PDF reproducibility are explicitly tested by varying data-processing order and included steps, and by carrying out a crystal-to-crystal data comparison. We identify situations in which the 3D-$\Delta$PDF is robust, and caution against a few particular cases which can lead to inconsistent 3D-$\Delta$PDFs. Although not all the approaches applied here-in will be valid across all systems, and a more in-depth analysis of some of the effects of the data processing steps may still needed, the methods collected here-in represent the start of a more systematic discussion about data processing and corrections in this field.


title: The Dense Warm Ionized Medium in the Inner Galaxy

abstract: Ionized interstellar gas is an important component of the interstellar medium and its lifecycle. The recent evidence for a widely distributed highly ionized warm interstellar gas with a density intermediate between the warm ionized medium (WIM) and compact HII regions suggests that there is a major gap in our understanding of the interstellar gas. Here we investigate the properties of the dense warm ionized medium (D-WIM) in the Milky Way using spectrally resolved SOFIA GREAT [NII] 205 micron line emission and Green Bank Telescope hydrogen radio recombination lines (RRL) data, supplemented by Herschel PACS [NII] 122 micron data, and spectrally resolved 12CO. We observed eight lines of sight in the 20deg <l < 30deg region in the Galactic plane. We derived the kinetic temperature, and the thermal and turbulent velocity dispersions from the [NII] and RRL linewidths. The regions with [NII] 205 micron emission are characterized by electron densities, n(e) ~ 10 to 35 cm(-3), temperatures from 3400 to 8500 K, and column densities N(N+) ~ 7e16 to 3e17 cm(-2). The ionized hydrogen column densities range from 6e20 to 1.7e21 cm(-2) and the fractional nitrogen ion abundance x(N+) ~1 to 3e-4, implying an enhanced nitrogen abundance at ~ 4.3 kpc from the Galactic Center. The [NII] 205 micron emission coincides with CO emission, although often with an offset in velocity, which suggests that the D-WIM gas is located in, or near, star-forming regions, which themselves are associated with molecular gas. These dense ionized regions are found to contribute > 50% of the observed [CII] intensity along these LOS. The kinetic temperatures we derive are too low to explain the presence of N+ resulting from electron collisional ionization and/or proton charge transfer of atomic nitrogen. Rather, these regions most likely are ionized by extreme ultraviolet radiation.


title: A hyperparameter-tuning approach to automated inverse planning

abstract: Radiotherapy inverse planning often requires planners to modify parameters in the treatment planning system's objective function to produce clinically acceptable plans. Due to the manual steps in this process, plan quality can vary depending on the planning time available and the planner's skills. This study investigates two hyperparameter-tuning methods for automated inverse planning. Because this framework does not train a model on previously-optimized plans, it can be readily adapted to practice pattern changes, and plan quality is not limited by that of a training cohort. We selected 10 patients who received lung SBRT using manually-generated clinical plans. We used random sampling (RS) and Bayesian optimization (BO) to tune parameters using linear-quadratic utility functions based on 11 clinical goals. Normalizing all plans to have PTV D95 equal to 48 Gy, we compared plan quality for the automatically-generated and manually-generated plans. We also investigated the impact of iteration count on the automatically-generated plans, comparing planning time and plan utility for RS and BO plans with and without stopping criteria. Without stopping criteria, the median planning time was 1.9 and 2.3 hours for RS and BO plans. The OAR doses in the RS and BO plans had a median percent difference (MPD) of 48.7% and 60.4% below clinical dose limits and an MPD of 2.8% and 3.3% below clinical plan doses. With stopping criteria, the utility decreased by an MPD of 5.3% and 3.9% for RS and BO plans, but the median planning time was reduced to 0.5 and 0.7 hours, and the OAR doses still had an MPD of 42.9% and 49.7% below clinical dose limits and an MPD of 0.3% and 1.8% below clinical plan doses. This study demonstrates that hyperparameter-tuning approaches to automated inverse planning can reduce active planning time with plan quality that is similar to or better than manually-generated plans.


title: Minimal Cycle Representatives in Persistent Homology using Linear Programming: an Empirical Study with User's Guide

abstract: Cycle representatives of persistent homology classes can be used to provide descriptions of topological features in data. However, the non-uniqueness of these representatives creates ambiguity and can lead to many different interpretations of the same set of classes. One approach to solving this problem is to optimize the choice of representative against some measure that is meaningful in the context of the data. In this work, we provide a study of the effectiveness and computational cost of several $\ell_1$-minimization optimization procedures for constructing homological cycle bases for persistent homology with rational coefficients in dimension one, including uniform-weighted and length-weighted edge-loss algorithms as well as uniform-weighted and area-weighted triangle-loss algorithms. We conduct these optimizations via standard linear programming methods, applying general-purpose solvers to optimize over column bases of simplicial boundary matrices. Our key findings are: (i) optimization is effective in reducing the size of cycle representatives, (ii) the computational cost of optimizing a basis of cycle representatives exceeds the cost of computing such a basis in most data sets we consider, (iii) the choice of linear solvers matters a lot to the computation time of optimizing cycles, (iv) the computation time of solving an integer program is not significantly longer than the computation time of solving a linear program for most of the cycle representatives, using the Gurobi linear solver, (v) strikingly, whether requiring integer solutions or not, we almost always obtain a solution with the same cost and almost all solutions found have entries in {-1, 0, 1} and therefore, are also solutions to a restricted $\ell_0$ optimization problem, and (vi) we obtain qualitatively different results for generators in Erd\H{o}s-R\'enyi random clique complexes.


title: A Monotone Approximate Dynamic Programming Approach for the Stochastic Scheduling, Allocation, and Inventory Replenishment Problem: Applications to Drone and Electric Vehicle Battery Swap Stations

abstract: There is a growing interest in using electric vehicles (EVs) and drones for many applications. However, battery-oriented issues, including range anxiety and battery degradation, impede adoption. Battery swap stations are one alternative to reduce these concerns that allow the swap of depleted for full batteries in minutes. We consider the problem of deriving actions at a battery swap station when explicitly considering the uncertain arrival of swap demand, battery degradation, and replacement. We model the operations at a battery swap station using a finite horizon Markov Decision Process model for the stochastic scheduling, allocation, and inventory replenishment problem (SAIRP), which determines when and how many batteries are charged, discharged, and replaced over time. We present theoretical proofs for the monotonicity of the value function and monotone structure of an optimal policy for special SAIRP cases. Due to the curses of dimensionality, we develop a new monotone approximate dynamic programming (ADP) method, which intelligently initializes a value function approximation using regression. In computational tests, we demonstrate the superior performance of the new regression-based monotone ADP method as compared to exact methods and other monotone ADP methods. Further, with the tests, we deduce policy insights for drone swap stations.


title: Bayesian inference under model misspecification using transport-Lagrangian distances: an application to seismic inversion

abstract: Model misspecification constitutes a major obstacle to reliable inference in many inverse problems. Inverse problems in seismology, for example, are particularly affected by misspecification of wave propagation velocities. In this paper, we focus on a specific seismic inverse problem - full-waveform moment tensor inversion - and develop a Bayesian framework that seeks robustness to velocity misspecification. A novel element of our framework is the use of transport-Lagrangian (TL) distances between observed and model predicted waveforms to specify a loss function, and the use of this loss to define a generalized belief update via a Gibbs posterior. The TL distance naturally disregards certain features of the data that are more sensitive to model misspecification, and therefore produces less biased or dispersed posterior distributions in this setting. To make the latter notion precise, we use several diagnostics to assess the quality of inference and uncertainty quantification, i.e., continuous rank probability scores and rank histograms. We interpret these diagnostics in the Bayesian setting and compare the results to those obtained using more typical Gaussian noise models and squared-error loss, under various scenarios of misspecification. Finally, we discuss potential generalizability of the proposed framework to a broader class of inverse problems affected by model misspecification.


title: Methods Included: Standardizing Computational Reuse and Portability with the Common Workflow Language

abstract: Computational Workflows are widely used in data analysis, enabling innovation and decision-making. In many domains (bioinformatics, image analysis, & radio astronomy) the analysis components are numerous and written in multiple different computer languages by third parties. However, many competing workflow systems exist, severely limiting portability of such workflows, thereby hindering the transfer of workflows between different systems, between different projects and different settings, leading to vendor lock-ins and limiting their generic re-usability. Here we present the Common Workflow Language (CWL) project which produces free and open standards for describing command-line tool based workflows. The CWL standards provide a common but reduced set of abstractions that are both used in practice and implemented in many popular workflow systems. The CWL language is declarative, which allows expressing computational workflows constructed from diverse software tools, executed each through their command-line interface. Being explicit about the runtime environment and any use of software containers enables portability and reuse. Workflows written according to the CWL standards are a reusable description of that analysis that are runnable on a diverse set of computing environments. These descriptions contain enough information for advanced optimization without additional input from workflow authors. The CWL standards support polylingual workflows, enabling portability and reuse of such workflows, easing for example scholarly publication, fulfilling regulatory requirements, collaboration in/between academic research and industry, while reducing implementation costs. CWL has been taken up by a wide variety of domains, and industries and support has been implemented in many major workflow systems.


title: Learning a Universal Template for Few-shot Dataset Generalization

abstract: Few-shot dataset generalization is a challenging variant of the well-studied few-shot classification problem where a diverse training set of several datasets is given, for the purpose of training an adaptable model that can then learn classes from new datasets using only a few examples. To this end, we propose to utilize the diverse training set to construct a universal template: a partial model that can define a wide array of dataset-specialized models, by plugging in appropriate components. For each new few-shot classification problem, our approach therefore only requires inferring a small number of parameters to insert into the universal template. We design a separate network that produces an initialization of those parameters for each given task, and we then fine-tune its proposed initialization via a few steps of gradient descent. Our approach is more parameter-efficient, scalable and adaptable compared to previous methods, and achieves the state-of-the-art on the challenging Meta-Dataset benchmark.


title: A numerical study of long-term multi-wavelength blazar variability

abstract: Decade-long monitoring of blazars at optical and infrared (OIR) wavelengths with the Small and Moderate Aperture Research Telescope System (SMARTS) in Chile and in $\gamma$-rays with the Fermi Large Area Telescope (LAT) has enabled the systematic study of their multi-wavelength long-term variability. In this work we investigate, from a theoretical perspective, the long-term variability properties of blazar emission by introducing an observationally motivated time-dependence to four main parameters of the one-zone leptonic model: injection luminosity of relativistic electrons, strength of magnetic field, Doppler factor, and external photon field luminosity. For the first time, we use both the probability density function and the power spectral density of the 10 year-long Fermi-LAT light curves to create variation patterns for the model parameters. Using as test beds two bright blazars from the SMARTS sample (PKS 2155-304 and 3C 273), we compute 10 year-long OIR, X-ray, and $\gamma$-ray model light curves for different varying parameters. We compare the findings of our theoretical investigation with multi-wavelength observations using various measures of variability. While no single-varying parameter simulation can explain all multi-wavelength variability properties, changes in the electron luminosity and external radiation field in PKS 2155-304 and 3C 273, respectively, can account for most of them. Our results motivate future time-dependent studies with coupling between two or more physical parameters to describe the multi-wavelength long-term blazar variability.


title: The Benefit Of Temporally-Strong Labels In Audio Event Classification

abstract: To reveal the importance of temporal precision in ground truth audio event labels, we collected precise (~0.1 sec resolution) "strong" labels for a portion of the AudioSet dataset. We devised a temporally strong evaluation set (including explicit negatives of varying difficulty) and a small strong-labeled training subset of 67k clips (compared to the original dataset's 1.8M clips labeled at 10 sec resolution). We show that fine-tuning with a mix of weak and strongly labeled data can substantially improve classifier performance, even when evaluated using only the original weak labels. For a ResNet50 architecture, d' on the strong evaluation data including explicit negatives improves from 1.13 to 1.41. The new labels are available as an update to AudioSet.


title: Variable Reduction For Quadratic Unconstrained Binary Optimization

abstract: Quadratic Unconstrained Binary Optimization models are useful for solving a diverse range of optimization problems. Constraints can be added by incorporating quadratic penalty terms into the objective, often with the introduction of slack variables needed for conversion of inequalities. This transformation can lead to a significant increase in the size and density of the problem. Herein, we propose an efficient approach for recasting inequality constraints that reduces the number of linear and quadratic variables. Experimental results illustrate the efficacy.


title: Cause and Effect: Hierarchical Concept-based Explanation of Neural Networks

abstract: In many scenarios, human decisions are explained based on some high-level concepts. In this work, we take a step in the interpretability of neural networks by examining their internal representation or neuron's activations against concepts. A concept is characterized by a set of samples that have specific features in common. We propose a framework to check the existence of a causal relationship between a concept (or its negation) and task classes. While the previous methods focus on the importance of a concept to a task class, we go further and introduce four measures to quantitatively determine the order of causality. Moreover, we propose a method for constructing a hierarchy of concepts in the form of a concept-based decision tree which can shed light on how various concepts interact inside a neural network towards predicting output classes. Through experiments, we demonstrate the effectiveness of the proposed method in explaining the causal relationship between a concept and the predictive behaviour of a neural network as well as determining the interactions between different concepts through constructing a concept hierarchy.


title: Subgraph Games in the Semi-Random Graph Process and Its Generalization to Hypergraphs

abstract: The semi-random graph process is a single-player game that begins with an empty graph on $n$ vertices. In each round, a vertex $u$ is presented to the player independently and uniformly at random. The player then adaptively selects a vertex $v$ and adds the edge $uv$ to the graph. For a fixed monotone graph property, the objective of the player is to force the graph to satisfy this property with high probability in as few rounds as possible. We focus on the problem of constructing a subgraph isomorphic to an arbitrary, fixed graph $G$. Let $\omega = \omega(n)$ be any function tending to infinity as $n \to \infty$. In (Omri Ben-Eliezer et al. "Semi-random graph process". In: Random Structures & Algorithms 56.3 (2020), pp. 648-675) it was proved that asymptotically almost surely one can construct $G$ in less than $n^{(d-1)/d} \omega$ rounds where $d \ge 2$ is the degeneracy of $G$. It was also proved that the result is sharp for $G = K_{d+1}$, that is, asymptotically almost surely it takes at least $n^{(d-1)/d} / \omega$ rounds to create $K_{d+1}$. Moreover, the authors conjectured that their general upper bound is sharp for all graphs $G$. We prove this conjecture here. We also consider a natural generalization of the process to $s$-uniform hypergraphs, the semi-random hypergraph process in which $r \ge 1$ vertices are presented at random, and the player then selects $s-r \ge 1$ vertices to form an edge of size~$s$. Our results for graphs easily generalize to hypergraphs when $r=1$; the threshold for constructing a fixed $s$-uniform hypergraph $G$ is, again, determined by the degeneracy of $G$. However, new challenges are mounting when $r \ge 2$; thresholds are not even known for complete hypergraphs. We provide bounds for this family and determine thresholds for some sparser hypergraphs.


title: Continuum Dark Matter

abstract: We initiate the study of dark matter models based on a gapped continuum. Dark matter consists of a mixture of states with a continuous mass distribution, which evolves as the universe expands. We present an effective field theory describing the gapped continuum, outline the structure of the Hilbert space and show how to deal with the thermodynamics of such a system. This formalism enables us to study the cosmological evolution and phenomenology of gapped continuum DM in detail. As a concrete example, we consider a weakly-interacting continuum (WIC) model, a gapped continuum counterpart of the familiar WIMP. The DM interacts with the SM via a Z-portal. The model successfully reproduces the observed relic density, while direct detection constraints are avoided due to the effect of continuum kinematics. The model has striking observational consequences, including continuous decays of DM states throughout cosmological history, as well as cascade decays of DM states produced at colliders. We also describe how the WIC theory can arise from a local, unitary scalar QFT propagating on a five-dimensional warped background with a soft wall.


title: Quantum Langevin dynamics of a charged particle in a magnetic field : Response function, position-velocity and velocity autocorrelation functions

abstract: We use the Quantum Langevin equation as a starting point to study the response function, the position-velocity correlation function and the velocity autocorrelation function of a charged Quantum Brownian particle in the presence of a magnetic field and linearly coupled to a heat bath via position coordinate. We study two bath models -- the Ohmic bath model and the Drude bath model -- and make a detailed comparison in various time-temperature regimes. For both bath models there is a competition between the cyclotron frequency and the viscous damping rate giving rise to a transition from an oscillatory to a monotonic behaviour as the damping rate is increased. In the zero point fluctuation dominated low temperature regime, non-trivial noise correlations lead to some interesting features in this transition. We study the role of the memory time scale which comes into play in the Drude model and study the effect of this additional time scale. We discuss the experimental implications of our analysis in the context of experiments in cold ions.


title: Information Theoretic Key Agreement Protocol based on ECG signals

abstract: Wireless body area networks (WBANs) are becoming increasingly popular as they allow individuals to continuously monitor their vitals and physiological parameters remotely from the hospital. With the spread of the SARS-CoV-2 pandemic, the availability of portable pulse-oximeters and wearable heart rate detectors has boomed in the market. At the same time, in 2020 we assisted to an unprecedented increase of healthcare breaches, revealing the extreme vulnerability of the current generation of WBANs. Therefore, the development of new security protocols to ensure data protection, authentication, integrity and privacy within WBANs are highly needed. Here, we targeted a WBAN collecting ECG signals from different sensor nodes on the individual's body, we extracted the inter-pulse interval (i.e., R-R interval) sequence from each of them, and we developed a new information theoretic key agreement protocol that exploits the inherent randomness of ECG to ensure authentication between sensor pairs within the WBAN. After proper pre-processing, we provide an analytical solution that ensures robust authentication; we provide a unique information reconciliation matrix, which gives good performance for all ECG sensor pairs; and we can show that a relationship between information reconciliation and privacy amplification matrices can be found. Finally, we show the trade-off between the level of security, in terms of key generation rate, and the complexity of the error correction scheme implemented in the system.


title: Low Diameter Monochromatic Covers of Complete Multipartite Graphs

abstract: Let the diameter cover number, $D^t_r(G)$, denote the least integer $d$ such that under any $r$-coloring of the edges of the graph $G$, there exists a collection of $t$ monochromatic subgraphs of diameter at most $d$ such that every vertex of $G$ is contained in at least one of the subgraphs. We explore the diameter cover number with two colors and two subgraphs when $G$ is a complete multipartite graph with at least three parts. We determine exactly the value of $D_2^2(G)$ for all complete tripartite graphs $G$, and almost all complete multipartite graphs with more than three parts.


title: Dark matter from dark photons

abstract: In this talk, I explained how the observed dark matter (DM) relic abundance can be accounted for in models composed of three sectors (the DM, the Standard Model (SM) and a light mediator) connected to each other. This scenario is explored in the context of the massive dark photon model in which the DM is connected to the SM through the kinetic mixing between the hypercharge gauge boson and the vector gauge boson associated to a new $U(1)'$ gauge symmetry. In such portal models, the DM relic abundance can be produced in no less than nine regimes along five dynamical mechanisms. In particular, I emphasize the sequential freeze-in dynamical mechanism which can occurs when the dark photon is massive and consists in two successive freeze-in mechanisms.


title: Speckle Interferometry at SOAR in 2020

abstract: The results of speckle interferometric observations at the 4.1 m Southern Astrophysical Research Telescope (SOAR) in 2020, as well as earlier unpublished data, are given, totaling 1735 measurements of 1288 resolved pairs and non-resolutions of 1177 targets. We resolved for the first time 59 new pairs or subsystems in known binaries, mostly among nearby dwarf stars. This work continues our long-term speckle program. Its main goal is to monitor orbital motion of close binaries, including members of high-order hierarchies and Hipparcos pairs in the solar neighborhood. We also report observations of 892 members of young moving groups and associations, where we resolved 103 new pairs.


title: A local Cauchy integral formula for slice-regular functions

abstract: We prove a Cauchy-type integral formula for slice-regular functions where the integration is performed on the boundary of an open subset of the quaternionic space, with no requirement of axial symmetry. In particular, we get a local Cauchy-type integral formula. As a step towards the proof, we provide a decomposition of a slice-regular function as a combination of two axially monogenic functions.


title: The Entropic Dynamics of Relativistic Quantum Fields in Curved Space-time

abstract: It has often been the case in history that the laws of physics have been used as the framework for understanding and implementing information processing. The tacit assumption is that the laws of physics are fundamental and that the notion of information is derived from these laws. Here we take the opposite view: the laws of physics are an application of the rules for processing information. In this dissertation we apply the Entropic Dynamics (ED) framework to construct a quantum dynamics for scalar fields in space-time. We begin by considering a toy model consisting of many interacting particles, resulting in the familiar Schrodinger equation for non-relativistic particles. Using a similar methodology, we construct a theory of quantum scalar fields in flat space-time that is relativistic, but not manifestly so. Here we also discuss a novel way in which the ED of quantum scalar fields appears to evade the so-called Wallstrom objection. To go further towards constructing a manifestly covariant quantum ED of fields on a curved space-time, both fixed and dynamical, we borrow from the "many-time" approaches of P. Weiss, P. Dirac, K. Kuchar, and C. Teitelboim. For a fixed background the result is a manifestly covariant ED of scalar fields that is in the spirit of the covariant quantum theories proposed by S. Tomonaga and J. Schwinger. However, the formalism is sufficiently flexible so as to allow for the possibility of modeling the back reaction of the quantum matter fields on a fully dynamical classical background. The simplest realization of this classical-quantum interaction shares some formal similarity to semi-classical gravity models, and the semi-classical Einstein equations, in particular. We consider such a theory and discuss its plausibility as a candidate for a quantum gravity theory.


title: Post-processing Multi-Model Medium-Term Precipitation Forecasts Using Convolutional Neural Networks

abstract: The goal of this study was to improve the post-processing of precipitation forecasts using convolutional neural networks (CNNs). Instead of post-processing forecasts on a per-pixel basis, as is usually done when employing machine learning in meteorological post-processing, input forecast images were combined and transformed into probabilistic output forecast images using fully convolutional neural networks. CNNs did not outperform regularized logistic regression. Additionally, an ablation analysis was performed. Combining input forecasts from a global low-resolution weather model and a regional high-resolution weather model improved performance over either one.


title: SA-GAN: Structure-Aware GAN for Organ-Preserving Synthetic CT Generation

abstract: In medical image synthesis, model training could be challenging due to the inconsistencies between images of different modalities even with the same patient, typically caused by internal status/tissue changes as different modalities are usually obtained at a different time. This paper proposes a novel deep learning method, Structure-aware Generative Adversarial Network (SA-GAN), that preserves the shapes and locations of in-consistent structures when generating medical images. SA-GAN is employed to generate synthetic computed tomography (synCT) images from magnetic resonance imaging (MRI) with two parallel streams: the global stream translates the input from the MRI to the CT domain while the local stream automatically segments the inconsistent organs, maintains their locations and shapes in MRI, and translates the organ intensities to CT. Through extensive experiments on a pelvic dataset, we demonstrate that SA-GAN provides clinically acceptable accuracy on both synCTs and organ segmentation and supports MR-only treatment planning in disease sites with internal organ status changes.


title: Hybrid Schr\"odinger-Feynman Simulation of Quantum Circuits With Decision Diagrams

abstract: Classical simulations of quantum computations are vital for the future development of this emerging technology. To this end, decision diagrams have been proposed as a complementary technique which frequently allows to tackle the inherent exponential complexity of these simulations. In the worst case, however, they still cannot escape this complexity. Additionally, while other techniques make use of all the available processing power, decision diagram-based simulation to date cannot exploit the many processing units of today's systems. In this work, we show that both problems can be tackled together by employing a hybrid Schr\"odinger-Feynman scheme for the simulation. More precisely, we show that realizing such a scheme with decision diagrams is indeed possible, we discuss the resulting problems in its realization, and propose solutions how they can be handled. Experimental evaluations confirm that this significantly advances the state of the art in decision diagram-based simulation -- allowing to simulate certain hard circuits within minutes that could not be simulated in a whole day thus far.


title: Time Evolution of Quantum Effects

abstract: For quantum effects $a$ and $b$ we define the $a$-evolution of $b$ at time $t$ denoted by $b(t\mid a)$. We interpret $b(t\mid a)$ as the influence that $a$ has on $b$ at time $t$ when $a$ occurs, but is not measured at time $t=0$. Using $b(t\mid a)$ we define the time-dependent sequential product $a[t]b$. This is interpreted as an effect that results from first measuring $a$ and then measuring $b$ after a time delay $t$. Various properties of $a[t]b$ are derived and it is shown that $a[t]b$ is constant in time if and only if $a$ and $b$ commute or $a$ is a multiple of a projection. These concepts are extended to observables for a quantum system. The ideas are illustrated with some examples.


title: On Performance of Energy Harvested Cooperative NOMA Under Imperfect CSI and Imperfect SIC

abstract: With the advent of 5G and the need for energy-efficient massively connected wireless networks, in this work, we consider an energy harvesting (EH) based multi-relay downlink cooperative non-orthogonal multiple access (NOMA) system with practical constraints. The base station serves NOMA users with the help of decode-and-forward based multiple EH relays, where relays harvest the energy from the base station's radio frequency. A relay is selected from the multiple K-relay by using a partial relay selection protocol. The system is considered to operate in half-duplex mode over a generalized independent and identical Nakagami$-m$ fading channel. The closed-form expression of outage probability and ergodic rate are derived for users, under the assumption of imperfect channel state information (CSI) and imperfect successive interference cancellation (SIC) at the receiver node. Expression of outage probability and ergodic rate for both users under the assumption of perfect CSI and perfect SIC are also presented. Further, the asymptotic expression for the outage probability is also shown. The derived analytical expressions are verified through Monte-Carlo simulations.


title: Stripe order enhanced superconductivity in the Hubbard model

abstract: Unidirectional ("stripe") charge-density-wave order has now been established as a ubiquitous feature in the phase diagram of the cuprate high temperature (HT) superconductors, where it generally competes with superconductivity (SC). None-the-less, on theoretical grounds it has been conjectured that stripe order (or other forms of "optimal" inhomogeneities) may play an essential positive role in the mechanism of HTSC. Here we report density matrix renormalization group studies of the Hubbard model on long 4 and 6 leg cylinders where the hopping matrix elements transverse to the long direction are periodically modulated - mimicing the effect of putative period-2 stripe order. We find even modest amplitude modulations can enhance the long-distance SC correlations by many orders of magnitude, and drive the system into a phase with a substantial spin gap and SC quasi-long-range-order with a Luttinger exponent, $K_{sc} \sim 1$.


title: Handlebody diagram algebras

abstract: In this paper we study handlebody versions of some classical diagram algebras, most prominently, handlebody versions of Temperley-Lieb, blob, Brauer, BMW, Hecke and Ariki-Koike algebras. Moreover, motivated by Green-Kazhdan-Lusztig's theory of cells, we reformulate the notion of (sandwich, inflated or affine) cellular algebras. We explain this reformulation and how all of the above algebras are part of this theory.


title: A generalized Haus master equation model for mode-locked class-B lasers

abstract: Using an asymptotic technique we develop a generalized version of class-B Haus partial differential equation mode-locking model that accounts for both the slow gain response to the averaged value of the field intensity and the fast gain dynamics on the scale comparable to the pulse duration. We show that unlike the conventional class-B Haus mode-locked model, our model is able to describe not only Q-switched instability of the fundamental mode-locked regime, but also the leading edge instability leading to harmonic mode-locked regimes with the increase of the pump power.


title: A pivoting Mechanically Controllable Break junction setup enabling partially wet phase MCB-junctions

abstract: A technique is presented, which creates MCB junctions that can be pivoted to any desirable angle. The MCB junction equipped with a specific glass liquid cell can be used to produce a MCB junction, of which the electrodes are covered with a microscopic layer of fluid, thus producing a partially wet phase MCB junction.


title: Slicing-Based AI Service Provisioning on Network Edge

abstract: Edge intelligence leverages computing resources on network edge to provide artificial intelligence (AI) services close to network users. As it enables fast inference and distributed learning, edge intelligence is envisioned to be an important component of 6G networks. In this article, we investigate AI service provisioning for supporting edge intelligence. First, we present the features and requirements of AI services. Then, we introduce AI service data management, and customize network slicing for AI services. Specifically, we propose a novel resource pooling method to jointly manage service data and network resources for AI services. A trace-driven case study demonstrates the effectiveness of the proposed resource pooling method. Through this study, we illustrate the necessity, challenge, and potential of AI service provisioning on network edge.


title: Palatini gravity with nonmetricity, torsion, and boundaries in metric and connection variables

abstract: We prove the equivalence in the covariant phase space of the metric and connection formulations for Palatini gravity, with nonmetricity and torsion, on a spacetime manifold with boundary. To this end, we will rely on the cohomological approach provided by the relative bicomplex framework. Finally, we discuss some of the physical implications derived from this equivalence in the context of singularity identification through curvature invariants.


title: Face Attributes as Cues for Deep Face Recognition Understanding

abstract: Deeply learned representations are the state-of-the-art descriptors for face recognition methods. These representations encode latent features that are difficult to explain, compromising the confidence and interpretability of their predictions. Most attempts to explain deep features are visualization techniques that are often open to interpretation. Instead of relying only on visualizations, we use the outputs of hidden layers to predict face attributes. The obtained performance is an indicator of how well the attribute is implicitly learned in that layer of the network. Using a variable selection technique, we also analyze how these semantic concepts are distributed inside each layer, establishing the precise location of relevant neurons for each attribute. According to our experiments, gender, eyeglasses and hat usage can be predicted with over 96% accuracy even when only a single neural output is used to predict each attribute. These performances are less than 3 percentage points lower than the ones achieved by deep supervised face attribute networks. In summary, our experiments show that, inside DCNNs optimized for face identification, there exists latent neurons encoding face attributes almost as accurately as DCNNs optimized for these attributes.


title: 3D Two-Hop Cellular Networks with Wireless Backhauled UAVs: Modeling and Fundamentals

abstract: In this paper, we characterize the performance of a three-dimensional (3D) two-hop cellular network in which terrestrial base stations (BSs) coexist with unmanned aerial vehicles (UAVs) to serve a set of ground user equipment (UE). In particular, a UE connects either directly to its serving terrestrial BS by an access link or connects first to its serving UAV which is then wirelessly backhauled to a terrestrial BS (joint access and backhaul). We consider realistic antenna radiation patterns for both BSs and UAVs using practical models developed by the third generation partnership project (3GPP). We assume a probabilistic channel model for the air-to-ground transmission, which incorporates both line-of-sight (LoS) and non-line-of-sight (NLoS) links. Assuming the max-power association policy, we study the performance of the network in both amplify-and-forward (AF) and decode-and-forward (DF) relaying protocols. Using tools from stochastic geometry, we analyze the joint distribution of distance and zenith angle of the closest (and serving) UAV to the origin in a 3D setting. Further, we identify and extensively study key mathematical constructs as the building blocks of characterizing the received signal-to-interference-plus-noise ratio (SINR) distribution. Using these results, we obtain exact mathematical expressions for the coverage probability in both AF and DF relaying protocols. Furthermore, considering the fact that backhaul links could be quite weak because of the downtilted antennas at the BSs, we propose and analyze the addition of a directional uptilted antenna at the BS that is solely used for backhaul purposes. The superiority of having directional antennas with wirelessly backhauled UAVs is further demonstrated via simulation.


title: Convergence of the boundary integral method for interfacial Stokes flow

abstract: Boundary integral numerical methods are among the most accurate methods for interfacial Stokes flow, and are widely applied. They have the advantage that only the boundary of the domain must be discretized, which reduces the number of discretization points and allows the treatment of complicated interfaces. Despite their popularity, there is no analysis of the convergence of these methods for interfacial Stokes flow. In practice, the stability of discretizations of the boundary integral formulation can depend sensitively on details of the discretization and on the application of numerical filters. We present a convergence analysis of the boundary integral method for Stokes flow, focusing on a rather general method for computing the evolution of an elastic capsule, viscous drop, or inviscid bubble in 2D strain and shear flows. The analysis clarifies the role of numerical filters in practical computations.


title: Attainment of a High Concentration of Salt Ions Near a Metallic or Dielectric Wall in a Salt Solution as a result of Electrical Image Forces Near the Wall

abstract: Electrical image potentials near a metallic or a dielectric wall of higher dielectric constant than that of the solution are attractive, and therefore, can concentrate salt ions near the wall. It has recently been observed that near a metallic surface, ions in room temperature ionic liquids precipitate (but not near a nonmetallic surface). It will be argued that a likely reason for why precipitation of ions in salt water, as a result of electrical image forces, has not as yet been observed is the existence of an energy barrier near a solid surface, resulting from the decrease of ion solvation as a result of the large decrease of the dielectric constant of water normal to a solid wall within a short distance from the wall. We will explore the conditions under which ions are able to get past this barrier and concentrate at a solid wall, either as a result of a reduction of this barrier caused by screening at high ion concentration or as a result of thermal activation over this solvation energy potential barrier.


title: Valley Isospin Controlled Fractional Quantum Hall States in Bilayer Graphene

abstract: A two-dimensional electron system placed in a magnetic field develops Landau levels, where strong Coulomb interactions lead to the appearance of many-body correlated ground states. Quantum numbers similar to the electron spin enable the understanding and control of complex ground state order and collective excitations. Owing to its spin, valley and orbital degrees of freedom, Bernal-stacked bilayer graphene offers a rich platform to pursue correlated phenomena in two dimensions. In this work, we fabricate dual-gated Bernal-stacked bilayer graphene devices and demonstrate unprecedented fine control over its valley isospin degrees of freedom using a perpendicular electric field. Higher sample quality enables us to probe regimes obscured by disorder in previous studies. We present evidence for a new even-denominator fractional quantum Hall state at filling factor {\nu} = 5/2. The 5/2 state is found to be spontaneously valley polarized in the limit of vanishing valley Zeeman splitting, consistent with a theoretical prediction made regarding the spin polarization of the Moore-Read state. In the vicinity of the even-denominator fractional quantum Hall states, we observe the appearance of the predicted Levin-Halperin daughter states of the Moore-Read Pfaffian wave function at {\nu}= 3/2, 7/2, and -1/2 and of the anti-Pfaffian at {\nu}= 5/2. These observations suggest the breaking of particle-hole symmetry in bilayer graphene. We construct a comprehensive valley polarization phase diagram for the Jain sequence fractional states surrounding filling factor 3/2. These results are well explained by a two-component composite fermion model, further demonstrating the SU(2) nature of the valley isospin in bilayer graphene. Our experiment paves the path for future efforts of manipulating the valley isospin in bilayer graphene to engineer exotic topological orders and quantum information processes.


title: Momentum Contrastive Voxel-wise Representation Learning for Semi-supervised Volumetric Medical Image Segmentation

abstract: Contrastive learning (CL) aims to learn useful representation without relying on expert annotations in the context of medical image segmentation. Existing approaches mainly contrast a single positive vector (i.e., an augmentation of the same image) against a set of negatives within the entire remainder of the batch by simply mapping all input features into the same constant vector. Despite the impressive empirical performance, those methods have the following shortcomings: (1) it remains a formidable challenge to prevent the collapsing problems to trivial solutions; and (2) we argue that not all voxels within the same image are equally positive since there exist the dissimilar anatomical structures with the same image. In this work, we present a novel Contrastive Voxel-wise Representation Learning (CVRL) method to effectively learn low-level and high-level features by capturing 3D spatial context and rich anatomical information along both the feature and the batch dimensions. Specifically, we first introduce a novel CL strategy to ensure feature diversity promotion among the 3D representation dimensions. We train the framework through bi-level contrastive optimization (i.e., low-level and high-level) on 3D images. Experiments on two benchmark datasets and different labeled settings demonstrate the superiority of our proposed framework. More importantly, we also prove that our method inherits the benefit of hardness-aware property from the standard CL approaches.


title: Trimmed Match Design for Randomized Paired Geo Experiments

abstract: How to measure the incremental Return On Ad Spend (iROAS) is a fundamental problem for the online advertising industry. A standard modern tool is to run randomized geo experiments, where experimental units are non-overlapping ad-targetable geographical areas (Vaver & Koehler 2011). However, how to design a reliable and cost-effective geo experiment can be complicated, for example: 1) the number of geos is often small, 2) the response metric (e.g. revenue) across geos can be very heavy-tailed due to geo heterogeneity, and furthermore 3) the response metric can vary dramatically over time. To address these issues, we propose a robust nonparametric method for the design, called Trimmed Match Design (TMD), which extends the idea of Trimmed Match (Chen & Au 2019) and furthermore integrates the techniques of optimal subset pairing and sample splitting in a novel and systematic manner. Some simulation and real case studies are presented. We also point out a few open problems for future research.


title: Efficient Least Squares Monte-Carlo Technique for PFE/EE Calculations

abstract: We describe a regression-based method, generally referred to as the Least Squares Monte Carlo (LSMC) method, to speed up exposure calculations of a portfolio. We assume that the portfolio contains several exotic derivatives that are priced using Monte-Carlo on each real world scenario and time step. Such a setting is often referred to as a Monte Carlo over a Monte Carlo or a Nested Monte Carlo method.


title: Measuring the User Satisfaction in a Recommendation Interface with Multiple Carousels

abstract: It is common for video-on-demand and music streaming services to adopt a user interface composed of several recommendation lists, i.e. widgets or swipeable carousels, each generated according to a specific criterion or algorithm (e.g. most recent, top popular, recommended for you, editors' choice, etc.). Selecting the appropriate combination of carousel has significant impact on user satisfaction. A crucial aspect of this user interface is that to measure the relevance a new carousel for the user it is not sufficient to account solely for its individual quality. Instead, it should be considered that other carousels will already be present in the interface. This is not considered by traditional evaluation protocols for recommenders systems, in which each carousel is evaluated in isolation, regardless of (i) which other carousels are displayed to the user and (ii) the relative position of the carousel with respect to other carousels. Hence, we propose a two-dimensional evaluation protocol for a carousel setting that will measure the quality of a recommendation carousel based on how much it improves upon the quality of an already available set of carousels. Our evaluation protocol takes into account also the position bias, i.e. users do not explore the carousels sequentially, but rather concentrate on the top-left corner of the screen. We report experiments on the movie domain and notice that under a carousel setting the definition of which criteria has to be preferred to generate a list of recommended items changes with respect to what is commonly understood.


title: On some properties of weak solutions to the Maxwell equations

abstract: This paper is concerned with weak solutions {e,h} in L^2 x L^2 of the time-dependent Maxwell equations. We show that these solutions obey an energy equality. Our method of proof is based on the approximation of {e,h} by its Steklov mean with respect to time t. This approximation technique is well-known for establishing integral estimates for weak solutions of parabolic equations. In addition we prove the uniqueness of {e,h}.


title: Simulation of particle identification with the cluster counting technique

abstract: In this paper we show the potential of the cluster counting technique for particle identification. Simulations based on Garfield++ software prove that this technique improves the particle separation capabilities with respect to the ones obtained with the traditional method of dE/dx. Moreover three different algorithms to reproduce the clusters number and the cluster size distribution with Geant4 software are discussed.


title: Visual analogy: Deep learning versus compositional models

abstract: Is analogical reasoning a task that must be learned to solve from scratch by applying deep learning models to massive numbers of reasoning problems? Or are analogies solved by computing similarities between structured representations of analogs? We address this question by comparing human performance on visual analogies created using images of familiar three-dimensional objects (cars and their subregions) with the performance of alternative computational models. Human reasoners achieved above-chance accuracy for all problem types, but made more errors in several conditions (e.g., when relevant subregions were occluded). We compared human performance to that of two recent deep learning models (Siamese Network and Relation Network) directly trained to solve these analogy problems, as well as to that of a compositional model that assesses relational similarity between part-based representations. The compositional model based on part representations, but not the deep learning models, generated qualitative performance similar to that of human reasoners.


title: Node Selection Toward Faster Convergence for Federated Learning on Non-IID Data

abstract: Federated Learning (FL) is a distributed learning paradigm that enables a large number of resource-limited nodes to collaboratively train a model without data sharing. The non-independent-and-identically-distributed (non-i.i.d.) data samples invoke discrepancies between the global and local objectives, making the FL model slow to converge. In this paper, we proposed Optimal Aggregation algorithm for better aggregation, which finds out the optimal subset of local updates of participating nodes in each global round, by identifying and excluding the adverse local updates via checking the relationship between the local gradient and the global gradient. Then, we proposed a Probabilistic Node Selection framework (FedPNS) to dynamically change the probability for each node to be selected based on the output of Optimal Aggregation. FedPNS can preferentially select nodes that propel faster model convergence. The unbiasedness of the proposed FedPNS design is illustrated and the convergence rate improvement of FedPNS over the commonly adopted Federated Averaging (FedAvg) algorithm is analyzed theoretically. Experimental results demonstrate the effectiveness of FedPNS in accelerating the FL convergence rate, as compared to FedAvg with random node selection.


title: Geometric Approaches to Quantum Fields and Strings at Strong Couplings

abstract: Geometric structures and dualities arise naturally in quantum field theories and string theory. In fact, these tools become very useful when studying strong coupling effects, where standard perturbative techniques can no longer be used. In this thesis we look at several conformal field theories in various dimensions. We first discuss the structure of the nilpotent networks stemming from T-brane deformations in 4D $\mathcal{N}=1$ theories and then go to the stringy origins of 6D superconformal field theories to realize deformations associated with T-branes in terms of simple combinatorial data. We then analyze non-perturbative generalizations of orientifold 3-planes (i.e. S-folds) in order to produce different 4D $\mathcal{N}=2$ theories. Afterwards, we turn our attention towards a few dualities found at strong coupling. For instance, abelian T-duality is known to be a full duality in string theory between type IIA and type IIB. Its nonabelian generalization, Poisson-Lie T-duality, has only been conjectured to be so. We show that Poisson-Lie symmetric $\sigma$-models are at least two-loop renormalizable and their $\beta$-functions are invariant under Poisson-Lie T-duality. Finally, we review recent progress leading to phenomenologically relevant dualities between M-theory on local $G_2$ spaces and F-theory on locally elliptically fibered Calabi-Yau fourfolds. In particular, we find that the 3D $\mathcal{N}=1$ effective field theory defined by M-theory on a local $Spin(7)$ space unifies the Higgs bundle data associated with 4D $\mathcal{N}=1$ M-theory and F-theory vacua. We finish with some comments on 3D interfaces at strong coupling.


title: Global ionospheric response to a periodic sequence of HSS/CIR events during the 2007-2008 solar minimum

abstract: In this study, we investigate the global ionospheric impact of high-speed solar wind streams/corotating interaction regions (HSS/CIR). A series of ten such events are identified between December 1st 2007 and April 29th 2008, characterized in the frequency domain by the main spectral peaks corresponding to 27, 13.5, 9 and 6.75 days. The spectra of solar wind magnetic field, speed and proton density, as well as those of the geomagnetic indices AE and SYM-H are solely dominated by these features. By contrast, the ionospheric NmF2 and to a lesser extent the hmF2 spectra have a much more complex structure, with secondary peaks adding to or replacing the main ones. We argue that this is evidence of the nonlinear nature of the magnetosphere-ionosphere coupling, highlighted particularly in the NmF2 ionospheric response. Additionally, we show that hmF2 is more closely correlated than NmF2 to all parameters describing the solar wind and geomagnetic activity. Finally, the ionospheric response shows higher correlation with Bz than any other solar wind parameter, and higher with SYM-H than AE, indicating that for the low-frequency part of the spectrum, high-latitude Joule heating and particle precipitation play a secondary role to that of prompt penetration electric fields in dictating the ionospheric response to geomagnetic activity, in the case of this sequence of HSS/CIR events.


title: Edges in Brain Networks: Contributions to Models of Structure and Function

abstract: Network models describe the brain as sets of nodes and edges that represent its distributed organization. So far, most discoveries in network neuroscience have prioritized insights that highlight distinct groupings and specialized functional contributions of network nodes. Importantly, these functional contributions are determined and expressed by the web of their interrelationships, formed by network edges. Here, we underscore the important contributions made by brain network edges for understanding distributed brain organization. Different types of edges represent different types of relationships, including connectivity and similarity among nodes. Adopting a specific definition of edges can fundamentally alter how we analyze and interpret a brain network. Furthermore, edges can associate into collectives and higher-order arrangements, describe time series, and form edge communities that provide insights into brain network topology complementary to the traditional node-centric perspective. Focusing on the edges, and the higher-order or dynamic information they can provide, discloses previously underappreciated aspects of structural and functional network organization.


title: The Multivariate Theory of Functional Connections: An n-Dimensional Constraint Embedding Technique Applied to Partial Differential Equations

abstract: The Theory of Functional Connections (TFC) is a functional interpolation framework founded upon the so-called constrained expression: a functional that expresses the family of all possible functions that satisfy some user-specified, linear constraints. These constrained expressions can be utilized to transform constrained problems into unconstrained ones. The benefits of doing so include faster solution times, more accurate solutions, and more robust convergence. This dissertation contains a comprehensive, self-contained presentation of the TFC theory beginning with simple univariate point constraints and ending with general linear constraints in $n$-dimensions; relevant mathematical theorems and clarifying examples are included throughout the presentation to expand and solidify the reader's understanding. Furthermore, this dissertation describes how TFC can be applied to estimate differential equations' solutions, its primary application to date. In addition, comparisons with other state-of-the-art algorithms that estimate differential equations' solutions are included to showcase the advantages and disadvantages of the TFC approach. Lastly, the aforementioned concepts are leveraged to estimate solutions of differential equations from the field of flexible body dynamics.


title: Listen with Intent: Improving Speech Recognition with Audio-to-Intent Front-End

abstract: Comprehending the overall intent of an utterance helps a listener recognize the individual words spoken. Inspired by this fact, we perform a novel study of the impact of explicitly incorporating intent representations as additional information to improve a recurrent neural network-transducer (RNN-T) based automatic speech recognition (ASR) system. An audio-to-intent (A2I) model encodes the intent of the utterance in the form of embeddings or posteriors, and these are used as auxiliary inputs for RNN-T training and inference. Experimenting with a 50k-hour far-field English speech corpus, this study shows that when running the system in non-streaming mode, where intent representation is extracted from the entire utterance and then used to bias streaming RNN-T search from the start, it provides a 5.56% relative word error rate reduction (WERR). On the other hand, a streaming system using per-frame intent posteriors as extra inputs for the RNN-T ASR system yields a 3.33% relative WERR. A further detailed analysis of the streaming system indicates that our proposed method brings especially good gain on media-playing related intents (e.g. 9.12% relative WERR on PlayMusicIntent).


title: I3T: Intensity Interferometry Imaging Telescope

abstract: We propose a new approach, based on the Hanbury Brown and Twiss intensity interferometry, to transform a Cherenkov telescope to its equivalent optical telescope. We show that, based on the use of photonics components borrowed from quantum-optical applications, we can recover spatial details of the observed source down to the diffraction limit of the Cherenkov telescope, set by its diameter at the mean wavelength of observation. For this, we propose to apply aperture synthesis techniques from pairwise and triple correlation of sub-pupil intensities, in order to reconstruct the image of a celestial source from its Fourier moduli and phase information, despite atmospheric turbulence. We examine the sensitivity of the method, i.e. limiting magnitude, and its implementation on existing or future high energy arrays of Cherenkov telescopes. We show that despite its poor optical quality compared to extremely large optical telescopes under construction, a Cherenkov telescope can provide diffraction limited imaging of celestial sources, in particular at the visible, down to violet wavelengths.


title: N-ary Huffman Encoding Using High-Degree Trees -- A Performance Comparison

abstract: In this paper we implement an n-ary Huffman Encoding and Decoding application using different degrees of tree structures. Our goal is to compare the performance of the algorithm in terms of compression ratio, decompression speed and weighted path length when using higher degree trees, compared to the 2-ary Huffman Code. The Huffman tree degrees that we compare are 2-ary, 3-ary, 4-ary, 5-ary, 6-ary, 7-ary, 8-ary and 16-mal. We also present the impact that branch prediction has on the performance of the n-ary Huffman Decoding.


title: Closed-form Characterization of the MGF of AoI in Energy Harvesting Status Update Systems

abstract: This paper considers a real-time status update system in which an energy harvesting (EH)-powered transmitter node observes some physical process, and sends its sensed measurements in the form of status updates to a destination node. The status update and harvested energy packets are assumed to arrive at the transmitter according to independent Poisson processes, and the service time of each status update is assumed to be exponentially distributed. We quantify the freshness of status updates when they reach the destination using the concept of Age of Information (AoI). Unlike most of the existing analyses of AoI focusing on the evaluation of its average value when the transmitter is not subject to energy constraints, our analysis is focused on understanding the distributional properties of AoI through the characterization of its moment generating function (MGF). In particular, we use the stochastic hybrid systems (SHS) framework to derive closed-form expressions of the MGF of AoI under several queueing disciplines at the transmitter, including non-preemptive and preemptive in service/waiting strategies. Using these MGF results, we further obtain closed-form expressions for the first and second moments of AoI in each queueing discipline. We demonstrate the generality of this analysis by recovering several existing results for the corresponding system with no energy constraints as special cases of the new results. Our numerical results verify the analytical findings, and demonstrate the necessity of incorporating the higher moments of AoI in the implementation/optimization of real-time status update systems rather than just relying on its average value.


title: Strong gravitational lens image of the M87 black hole with a simple accreting matter model

abstract: We study simulated images generated from an accretion disk surrounding the supermassive black hole hosted in the nearby galaxy M87. We approach the problem employing very simple accreting models inspired from magnetohydrodynamical simulations and introducing a new recipe for dealing with the combined integration of the geodesic and geodesic deviation equations in Kerr spacetime, which allows for a convenient and efficient way to manage the system of equations. The geometry of the basic emission model is given by a two temperature thin disk in the equatorial plane of the black hole supplemented by an asymmetric bar structure. We show that this configuration permits to generate the most salient features appearing in the EHT Collaboration images of M87 with impressive fidelity.


title: Efficient Algorithms for Constructing an Interpolative Decomposition

abstract: Low-rank approximations are essential in modern data science. The interpolative decomposition provides one such approximation. Its distinguishing feature is that it reuses columns from the original matrix. This enables it to preserve matrix properties such as sparsity and non-negativity. It also helps save space in memory. In this work, we introduce two optimized algorithms to construct an interpolative decomposition along with numerical evidence that they outperform the current state of the art.


title: FORESEE: FORward Experiment SEnsitivity Estimator for the LHC and future hadron colliders

abstract: We introduce a numerical package FORward Experiment SEnsitivity Estimator, or FORESEE, that can be used to simulate the expected sensitivity reach of experiments placed in the far-forward direction from the proton-proton interaction point. The simulations can be performed for $14$ TeV collision energy characteristic for the LHC, as well as for larger energies: $27$ and $100$ TeV. In the package, a comprehensive list of validated forward spectra of various SM species is also provided. The capabilities of FORESEE are illustrated for the popular dark photon and dark Higgs boson models, as well as for the search for light up-philic scalars. For the dark photon portal, we also comment on the complementarity between such searches and dark matter direct detection bounds. Additionally, for the first time, we discuss the prospects for the LLP searches in the proposed future hadron colliders: High-Energy LHC (HE-LHC), Super proton-proton Collider (SppC), and Future Circular Collider (FCC-hh).


title: High-Robustness, Low-Transferability Fingerprinting of Neural Networks

abstract: This paper proposes Characteristic Examples for effectively fingerprinting deep neural networks, featuring high-robustness to the base model against model pruning as well as low-transferability to unassociated models. This is the first work taking both robustness and transferability into consideration for generating realistic fingerprints, whereas current methods lack practical assumptions and may incur large false positive rates. To achieve better trade-off between robustness and transferability, we propose three kinds of characteristic examples: vanilla C-examples, RC-examples, and LTRC-example, to derive fingerprints from the original base model. To fairly characterize the trade-off between robustness and transferability, we propose Uniqueness Score, a comprehensive metric that measures the difference between robustness and transferability, which also serves as an indicator to the false alarm problem.


title: Dynamic network analysis improves protein 3D structural classification

abstract: Protein structural classification (PSC) is a supervised problem of assigning proteins into pre-defined structural (e.g., CATH or SCOPe) classes based on the proteins' sequence or 3D structural features. We recently proposed PSC approaches that model protein 3D structures as protein structure networks (PSNs) and analyze PSN-based protein features, which performed better than or comparable to state-of-the-art sequence or other 3D structure-based approaches in the task of PSC. However, existing PSN-based PSC approaches model the whole 3D structure of a protein as a static PSN. Because folding of a protein is a dynamic process, where some parts of a protein fold before others, modeling the 3D structure of a protein as a dynamic PSN might further help improve the existing PSC performance. Here, we propose for the first time a way to model 3D structures of proteins as dynamic PSNs, with the hypothesis that this will improve upon the current state-of-the-art PSC approaches that are based on static PSNs (and thus upon the existing state-of-the-art sequence and other 3D structural approaches). Indeed, we confirm this on 71 datasets spanning ~44,000 protein domains from CATH and SCOPe


title: Iterative Algorithms for Assessing Network Resilience Against Structured Perturbations

abstract: This paper studies network resilience against structured additive perturbations to its topology. We consider dynamic networks modeled as linear time-invariant systems subject to perturbations of bounded energy satisfying specific sparsity and entry-wise constraints. Given an energy level, the structured pseudospectral abscissa captures the worst-possible perturbation an adversary could employ to de-stabilize the network, and the structured stability radius is the maximum energy in the structured perturbation that the network can withstand without becoming unstable. Building on a novel characterization of the worst-case structured perturbation, we propose iterative algorithms that efficiently compute the structured pseudospectral abscissa and structured stability radius. We provide theoretical guarantees of the local convergence of the algorithms and illustrate their efficacy and accuracy on several network examples.


title: Rare-Earth Molecular Crystals with Ultra-narrow Optical Linewidths for Photonic Quantum Technologies

abstract: Rare-earth ions are promising solid state systems to build light-matter interfaces at the quantum level. This relies on their potential to show narrow optical homogeneous linewidths or, equivalently, long-lived optical quantum states. In this letter, we report on europium molecular crystals that exhibit linewidths in the 10s of kHz range, orders of magnitude narrower than other molecular centers. We harness this property to demonstrate efficient optical spin initialization, coherent storage of light using an atomic frequency comb, and optical control of ion-ion interactions towards implementation of quantum gates. These results illustrate the utility of rare-earth molecular crystals as a new platform for photonic quantum technologies that combines highly coherent emitters with the unmatched versatility in composition, structure, and integration capability of molecular materials.


title: Interpretable Drug Synergy Prediction with Graph Neural Networks for Human-AI Collaboration in Healthcare

abstract: We investigate molecular mechanisms of resistant or sensitive response of cancer drug combination therapies in an inductive and interpretable manner. Though deep learning algorithms are widely used in the drug synergy prediction problem, it is still an open problem to formulate the prediction model with biological meaning to investigate the mysterious mechanisms of synergy (MoS) for the human-AI collaboration in healthcare systems. To address the challenges, we propose a deep graph neural network, IDSP (Interpretable Deep Signaling Pathways), to incorporate the gene-gene as well as gene-drug regulatory relationships in synergic drug combination predictions. IDSP automatically learns weights of edges based on the gene and drug node relations, i.e., signaling interactions, by a multi-layer perceptron (MLP) and aggregates information in an inductive manner. The proposed architecture generates interpretable drug synergy prediction by detecting important signaling interactions, and can be implemented when the underlying molecular mechanism encounters unseen genes or signaling pathways. We test IDWSP on signaling networks formulated by genes from 46 core cancer signaling pathways and drug combinations from NCI ALMANAC drug combination screening data. The experimental results demonstrated that 1) IDSP can learn from the underlying molecular mechanism to make prediction without additional drug chemical information while achieving highly comparable performance with current state-of-art methods; 2) IDSP show superior generality and flexibility to implement the synergy prediction task on both transductive tasks and inductive tasks. 3) IDSP can generate interpretable results by detecting different salient signaling patterns (i.e. MoS) for different cell lines.


title: Assessing the detectability of the secondary spin in extreme mass-ratio inspirals with fully-relativistic numerical waveforms

abstract: Extreme mass-ratio inspirals~(EMRIs) detectable by the Laser Inteferometric Space Antenna~(LISA) are unique probes of astrophysics and fundamental physics. Parameter estimation for these sources is challenging, especially because the waveforms are long, complicated, known only numerically, and slow to compute in the most relevant regime, where the dynamics is relativistic. We perform a time-consuming Fisher-matrix error analysis of the EMRI parameters using fully-relativistic numerical waveforms to leading order in an adiabatic expansion on a Kerr background, taking into account the motion of the LISA constellation, higher harmonics, and also including the leading correction from the spin of the secondary in the post-adiabatic approximation. We pay particular attention to the convergence of the numerical derivatives in the Fisher matrix and to the numerical stability of the covariance matrix, which for some systems requires computing the numerical waveforms with approximately $90$-digit precision. Our analysis confirms previous results (obtained with approximated but much more computationally efficient waveforms) for the measurement errors on the binary's parameters. We also show that the inclusion of higher harmonics improves the errors on the luminosity distance and on the orbital angular momentum angles by one order and two orders of magnitude, respectively, which might be useful to identify the environments where EMRIs live. We particularly focus on the measurability of the spin of the secondary, confirming that it cannot be measured with sufficient accuracy. However, due to correlations, its inclusion in the waveform model can deteriorate the accuracy on the measurements of other parameters by orders of magnitude, unless a physically-motivated prior on the secondary spin is imposed.


title: Monodromies of Projective Structures on Surface of Finite-type

abstract: We characterize the monodromies of projective structures with fuchsian-type singularities. Namely, any representation from the fundamental group of a Riemann surface of finite-type in $PSL_2(\mathbb{C})$ can be represented as the holonomy of branched projective structure with fuchsian-type singularities over the cusps. We made a geometrical/topological study of all local conical projective structures whose Schwarzian derivative admits a simple pole at the cusp. Finally, we explore isomonodromic deformations of such projective structures and the problem of minimizing angles.


title: MutualNet: Adaptive ConvNet via Mutual Learning from Different Model Configurations

abstract: Most existing deep neural networks are static, which means they can only do inference at a fixed complexity. But the resource budget can vary substantially across different devices. Even on a single device, the affordable budget can change with different scenarios, and repeatedly training networks for each required budget would be incredibly expensive. Therefore, in this work, we propose a general method called MutualNet to train a single network that can run at a diverse set of resource constraints. Our method trains a cohort of model configurations with various network widths and input resolutions. This mutual learning scheme not only allows the model to run at different width-resolution configurations but also transfers the unique knowledge among these configurations, helping the model to learn stronger representations overall. MutualNet is a general training methodology that can be applied to various network structures (e.g., 2D networks: MobileNets, ResNet, 3D networks: SlowFast, X3D) and various tasks (e.g., image classification, object detection, segmentation, and action recognition), and is demonstrated to achieve consistent improvements on a variety of datasets. Since we only train the model once, it also greatly reduces the training cost compared to independently training several models. Surprisingly, MutualNet can also be used to significantly boost the performance of a single network, if dynamic resource constraint is not a concern. In summary, MutualNet is a unified method for both static and adaptive, 2D and 3D networks. Codes and pre-trained models are available at \url{https://github.com/taoyang1122/MutualNet}.


title: Divergence Estimation in Message Passing algorithms

abstract: Many modern imaging applications can be modeled as compressed sensing linear inverse problems. When the measurement operator involved in the inverse problem is sufficiently random, denoising Scalable Message Passing (SMP) algorithms have a potential to demonstrate high efficiency in recovering compressed data. One of the key components enabling SMP to achieve fast convergence, stability and predictable dynamics is the Onsager correction that must be updated at each iteration of the algorithm. This correction involves the denoiser's divergence that is traditionally estimated via the Black-Box Monte Carlo (BB-MC) method \cite{MC-divergence}. While the BB-MC method demonstrates satisfying accuracy of estimation, it requires executing the denoiser additional times at each iteration and might lead to a substantial increase in computational cost of the SMP algorithms. In this work we develop two Large System Limit models of the Onsager correction for denoisers operating within SMP algorithms and use these models to propose two practical classes of divergence estimators that require no additional executions of the denoiser and demonstrate similar or superior correction compared to the BB-MC method.


title: Effective description of bistability and irreversibility in apoptosis

abstract: Apoptosis is a mechanism of programmed cell death in which cells engage in a controlled demolition and prepare to be digested without damaging their environment. In normal conditions apoptosis is repressed, until it is irreversibly induced by an appropriate signal. In adult organisms apoptosis is a natural way to dispose of damaged cells, and its disruption or excess is associated with cancer and autoimmune diseases. Apoptosis is regulated by a complex signaling network controlled by caspases, specialized enzymes that digest essential cellular components and promote the degradation of genomic DNA. In this work we propose an effective description of the signaling network focused on caspase-3 as a readout of cell fate. We integrate intermediate network interactions into a nonlinear feedback function acting on caspase-3 and introduce the effect of pro-apoptotic stimuli and regulatory elements as a saturating activation function. We find that the theory has a robust bistable regime where two possible states coexist, representing survival and cell death fates. For a broad range of parameters, strong stimuli can induce an irreversible switch to the death fate. We use the theory to explore dynamical stimulation conditions and determine how cell fate depends on different stimuli patterns. This analysis reveals a critical relation between transient stimuli intensity and duration to trigger irreversible apoptosis.


title: On Achilles Heel of Some Optical Network Designs and Performance Comparisons

abstract: This non-conventional paper represents the first attempt to uncover a possible vulnerability in some proposals for optical network designs and performance comparisons. While optical network designs and planning lie at the heart of achieving fiber capacity efficiency and/or operational efficiency, its combinatorial nature makes it computationally hard to reach optimal solutions for realistic scenarios. Therefore, the well-established way that have been taken for granted by not-so-small number of research papers is that an optimization model based on mixed integer linear programming (MILP) is first proposed and then due to the intractability of such combinatorial model, an heuristic algorithm is offered as an approximation. The solution-quality comparison between the MILP and heuristic is then carried out on small-scale instances including topologies and traffic tests to verify the efficacy of the proposed heuristic and the next step is to use such allegedly verified heuristic for optical network designs of realistic scenarios. This approach may nevertheless leave a critical vulnerability as there is no guarantee that one performs well in small tests will generalize adequately for large-scale cases, a common pitfall widely referred as the peril of extrapolation and/or overfitting. Besides, it is not uncommon that in some research works, for benchmarking purpose, the comparison between a new design proposal whose performance is obtained from on one heuristic and a reference design based on another heuristic is carried out. As the result of missing solution quality check, such performance comparison relied merely on heuristic solutions may be equally vulnerable as its results can be distorted and thus, be far from the possibly achieved zones. In this work, we pinpoint those issues and provide a realistic case study to highlight and demonstrate the impact of such vulnerabilities.


title: On the Stability Region of Intermittent Interference Networks

abstract: Recent information-theoretic studies have resulted in several interference management (IM) techniques that promise significant capacity improvements over interference avoidance techniques. However, in practice, the stable throughput region is a more relevant metric compared to the capacity region. In this work, we focus on the stable throughput region of a two-pair intermittent interference network with distributed transmitters and propose a queue-based transmission protocol in different regimes to handle the data between queues. In this context, we translate physical-layer IM protocols to accommodate stochastic message arrivals. To evaluate our proposed techniques, we compare the stable throughput region to the capacity region and show, through simulations, that the stable throughput region matches the capacity region when the latter is known. We show that in order to achieve the optimal stable throughput region, new ingredients are needed when compared to prior results. We quantify the trade-off between the encoding/decoding complexity of the proposed scheme (in terms of number of required algebraic operations), and the achievable rates. Finally, we study the lifetime of messages (i.e. the duration from arrival to successful delivery) vis-a-vis the total communication time, and we observe that the average lifetime scales as the square root of the total communication time.


title: Quadratic transformation and matrix biorthogonal polynomials: an $\mathcal{LU}$ factorization approach

abstract: The manuscript presents the $LU$ approach to matrix biorthogonal polynomials when all the even ordered entries in the Gram matrix are zero. This arises in case of a quadratic transformation which is briefly discussed. Further, the main diagonal of the Gram matrix is a zero diagonal and we present the theory that follows from this fact. Precisely, we discuss the Christoffel transformation and matrix representations of the kernel polynomials, usually called the ABC Theorem. Finally, we provide an illustration of our results assuming the Gram matrix has Hankel symmetry.


title: Verification of Image-based Neural Network Controllers Using Generative Models

abstract: Neural networks are often used to process information from image-based sensors to produce control actions. While they are effective for this task, the complex nature of neural networks makes their output difficult to verify and predict, limiting their use in safety-critical systems. For this reason, recent work has focused on combining techniques in formal methods and reachability analysis to obtain guarantees on the closed-loop performance of neural network controllers. However, these techniques do not scale to the high-dimensional and complicated input space of image-based neural network controllers. In this work, we propose a method to address these challenges by training a generative adversarial network (GAN) to map states to plausible input images. By concatenating the generator network with the control network, we obtain a network with a low-dimensional input space. This insight allows us to use existing closed-loop verification tools to obtain formal guarantees on the performance of image-based controllers. We apply our approach to provide safety guarantees for an image-based neural network controller for an autonomous aircraft taxi problem. We guarantee that the controller will keep the aircraft on the runway and guide the aircraft towards the center of the runway. The guarantees we provide are with respect to the set of input images modeled by our generator network, so we provide a recall metric to evaluate how well the generator captures the space of plausible images.


title: Influential groups for seeding and sustaining nonlinear contagion in heterogeneous hypergraphs

abstract: Several biological and social contagion phenomena, such as superspreading events or social reinforcement, are the results of multi-body interactions, for which hypergraphs offer a natural mathematical description. In this paper, we develop a novel mathematical framework based on approximate master equations to study contagions on random hypergraphs with a heterogeneous structure, both in terms of group size (hyperedge cardinality) and of membership of nodes to groups (hyperdegree). The characterization of the inner dynamics of groups provides an accurate description of the contagion process, without losing the analytical tractability. Using a contagion model where multi-body interactions are mapped onto a nonlinear infection rate, our two main results show how large groups are influential, in the sense that they drive both the early spread of a contagion and its endemic state (i.e., its stationary state). First, we provide a detailed characterization of the phase transition, which can be continuous or discontinuous with a bistable regime, and derive analytical expressions for the critical and tricritical points. We find that large values of the third moment of the membership distribution suppress the emergence of a discontinuous phase transition. Furthermore, the combination of heterogeneous group sizes and nonlinear contagion facilitates the onset of a mesoscopic localization phase, where contagion is sustained only by the largest groups, thereby inhibiting bistability as well. Second, we formulate the problem of optimal seeding for hypergraph contagion, and we compare two strategies: tuning the allocation of seeds according to either individual node or group properties. We find that, when the contagion is sufficiently nonlinear, groups are more effective seeds of contagion than individual nodes.


title: Origins of electronic bands in antiferromagnetic topological insulator MnBi$_2$Te$_4$

abstract: Despite the rapid progress in understanding the first intrinsic magnetic topological insulator MnBi$_2$Te$_4$, its electronic structure remains a topic under debates. Here we perform a thorough spectroscopic investigation into the electronic structure of MnBi$_2$Te$_4$ via laser-based angle-resolved photoemission spectroscopy. Through quantitative analysis, we estimate an upper bound of 3 meV for the gap size of the topological surface state. Furthermore, our circular dichroism measurements reveal band chiralities for both the topological surface state and quasi-2D bands, which can be well reproduced in a band hybridization model. A numerical simulation of energy-momentum dispersions based on a four-band model with an additional step potential near the surface provides a promising explanation for the origin of the quasi-2D bands. Our study represents a solid step forward in reconciling the existing controversies in the electronic structure of MnBi$_2$Te$_4$, and provides an important framework to understand the electronic structures of other relevant topological materials MnBi$_{2n}$Te$_{3n+1}$.


title: Two examples based on the properties of discrete measures

abstract: In the paper we represent two examples which are based on the properties of discrete measures. In the first part of the paper we prove that for each probability measure $\mu$, $\operatorname{supp}{\mu}=[-1,1]$, which logarithmic potential is a continuous function on $[-1,1]$ there exists a (discrete) measure $\sigma=\sigma(\mu)$, $\operatorname{supp}{\sigma}=[-1,1]$, with the following property. Let $\{P_n(x;\sigma)\}$ be the sequence of polynomials orthogonal with respect to $\sigma$. Then $\dfrac1n\chi(P_n(\cdot;\sigma))\overset{*}\to\mu$, $n\to\infty$, where $\chi(\cdot)$ is zero counting measure for the corresponding polynomial. The construction of the measure $\sigma$ is based of the properties of weighted Leja points. In the second part we give an example of a compact set and a sequence of discrete measures supported on that compact set with the following property. The sequence of measures converges in weak-$*$ topology to the equilibrium measure for the compact set but the corresponding sequence of the logarithmic potentials does not converge to the equilibrium potential in any neighbourhood of the compact set.


title: Topological Rigidity and Positive scalar curvature

abstract: In this paper, we study the topological rigidity and its relationship with the positivity of scalar curvature. Precisely, we show that any complete contractible $3$-manifold with non-negative scalar curvature is homeomorphic to $\mathbf{R}^3$.


title: Thompson-like groups, Reidemeister numbers, and fixed points

abstract: We investigate fixed-point properties of automorphisms of groups similar to R. Thompson's group $F$. Revisiting work of Gon\c{c}alves-Kochloukova, we deduce a cohomological criterion to detect infinite fixed-point sets in the abelianization, implying the so-called property $R_\infty$. Using the BNS $\Sigma$-invariant and drawing from works of Gon\c{c}alves-Sankaran-Strebel and Zaremsky, we show that our tool applies to many $F$-like groups, including Stein's $F_{2,3}$, Cleary's $F_\tau$, the Lodha-Moore groups, and the braided version of $F$.


title: Identifying Wisdom (of the Crowd): A Regression Approach

abstract: Experts in a population hold (a) beliefs over a state (call these state beliefs), as well as (b) beliefs over the distribution of beliefs in the population (call these hypothetical beliefs). If these are generated via updating a common prior using a fixed information structure, then the information structure can (generically) be derived by regressing hypothetical beliefs on state beliefs, provided there are at least as many signals as states. In addition, the prior solves an eigenvector equation derived from a matrix determined by the state beliefs and the hypothetical beliefs. Thus, the ex-ante informational environment (i.e., how signals are generated) can be determined using ex-post data (i.e., the beliefs in the population). I discuss implications of this finding, as well as what is identified when there are more states than signals.


title: Stability of the Phase Separation State for Compressible Navier-Stokes/Allen-Cahn System

abstract: This paper is concerned with the large time behavior of the Cauchy problem for Navier-Stokes/Allen-Cahn system describing the interface motion of immiscible two-phase flow in 3-D. The existence and uniqueness of global solutions and the stability of the phase separation state is proved under the small initial perturbations. Moreover, the optimal time decay rates are obtained for higher-order spatial derivatives of density, velocity and phase. Our results implies that if the immiscible two-phase flow is initially located near the phase separation state, then under small perturbation conditions, the solution exists globally and decays algebraically to the complete separation state of the two-phase flow, that is, there will be no interface fracture, vacuum, shock wave, mass concentration at any time, and the interface thickness tends to zero as the time $t\rightarrow+\infty$.


title: Feature-Based Interpretable Reinforcement Learning based on State-Transition Models

abstract: Growing concerns regarding the operational usage of AI models in the real-world has caused a surge of interest in explaining AI models' decisions to humans. Reinforcement Learning is not an exception in this regard. In this work, we propose a method for offering local explanations on risk in reinforcement learning. Our method only requires a log of previous interactions between the agent and the environment to create a state-transition model. It is designed to work on RL environments with either continuous or discrete state and action spaces. After creating the model, actions of any agent can be explained in terms of the features most influential in increasing or decreasing risk or any other desirable objective function in the locality of the agent. Through experiments, we demonstrate the effectiveness of the proposed method in providing such explanations.


title: Convergence of the Scalar- and Vector-Valued Allen-Cahn Equation to Mean Curvature Flow with $90${\deg}-Contact Angle in Higher Dimensions

abstract: We consider the sharp interface limit for the scalar-valued and vector-valued Allen-Cahn equation with homogeneous Neumann boundary condition in a bounded smooth domain $\Omega$ of arbitrary dimension $N\geq 2$ in the situation when a two-phase diffuse interface has developed and intersects the boundary $\partial\Omega$. The limit problem is mean curvature flow with $90${\deg}-contact angle and we show convergence in strong norms for well-prepared initial data as long as a smooth solution to the limit problem exists. To this end we assume that the limit problem has a smooth solution on $[0,T]$ for some time $T>0$. Based on the latter we construct suitable curvilinear coordinates and set up an asymptotic expansion for the scalar-valued and the vector-valued Allen-Cahn equation. Finally, we prove a spectral estimate for the linearized Allen-Cahn operator in both cases in order to estimate the difference of the exact and approximate solutions with a Gronwall-type argument.


title: Calculating black hole shadows: Review of analytical studies

abstract: In this article, we provide a review of the current state of the research of the black hole shadow, focusing on analytical (as opposed to numerical and observational) studies. We start with particular attention to the definition of the shadow and its relation to the often used concepts of escape cone, critical impact parameter and particle cross-section. For methodological purposes, we present the derivation of the angular size of the shadow for an arbitrary spherically symmetric and static space-time, which allows one to calculate the shadow for an observer at arbitrary distance from the center. Then we discuss the calculation of the shadow of a Kerr black hole, for an observer anywhere outside of the black hole. For observers at large distances we present and compare two methods used in the literature. Special attention is given to calculating the shadow in space-times which are not asymptotically flat. Shadows of wormholes and other black-hole impostors are reviewed. Then we discuss the calculation of the black hole shadow in an expanding universe as seen by a comoving observer. The influence of a plasma on the shadow of a black hole is also considered.