title: Diverging quantum speed limits: a herald of classicality

abstract: When is the quantum speed limit (QSL) really quantum? While vanishing QSL times often indicate emergent classical behavior, it is still not entirely understood what precise aspects of classicality are at the origin of this dynamical feature. Here, we show that vanishing QSL times (or, equivalently, diverging quantum speeds) can be traced back to reduced uncertainty in quantum observables and thus can be understood as a consequence of emerging classicality for these particular observables. We illustrate this mechanism by developing a QSL formalism for continuous variable quantum systems undergoing general Gaussian dynamics. For these systems, we show that three typical scenarios leading to vanishing QSL times, namely large squeezing, small effective Planck's constant, and large particle number, can be fundamentally connected to each other. In contrast, by studying the dynamics of open quantum systems and mixed states, we show that the classicality that emerges due to incoherent mixing of states from the addition of classical noise typically increases the QSL time.


title: On the Performance Analysis of the Adversarial System Variant Approximation Method to Quantify Process Model Generalization

abstract: Process mining algorithms discover a process model from an event log. The resulting process model is supposed to describe all possible event sequences of the underlying system. Generalization is a process model quality dimension of interest. A generalization metric should quantify the extent to which a process model represents the observed event sequences contained in the event log and the unobserved event sequences of the system. Most of the available metrics in the literature cannot properly quantify the generalization of a process model. A recently published method [1] called Adversarial System Variant Approximation leverages Generative Adversarial Networks to approximate the underlying event sequence distribution of a system from an event log. While this method demonstrated performance gains over existing methods in measuring the generalization of process models, its experimental evaluations have been performed under ideal conditions. This paper experimentally investigates the performance of Adversarial System Variant Approximation under non-ideal conditions such as biased and limited event logs. Moreover, experiments are performed to investigate the originally proposed sampling hyperparameter value of the method on its performance to measure the generalization. The results confirm the need to raise awareness about the working conditions of the Adversarial System Variant Approximation method. The outcomes of this paper also serve to initiate future research directions. [1] Theis, Julian, and Houshang Darabi. "Adversarial System Variant Approximation to Quantify Process Model Generalization." IEEE Access 8 (2020): 194410-194427.


title: Dyson-Schwinger equations in scalar electrodynamics

abstract: Counterparts of the Dyson-Schwinger equations for scalar QED in an external electromagnetic field are derived. Exact structure and diagrammatic interpretation of the corresponding mass and polarization operators are obtained. It is shown that in the presence of an external field the final equations include an exact three-photon interaction vertex.


title: A New Multipoint Symmetric Secant Method with a Dense Initial Matrix

abstract: In large-scale optimization, when either forming or storing Hessian matrices are prohibitively expensive, quasi-Newton methods are often used in lieu of Newton's method because they only require first-order information to approximate the true Hessian. Multipoint symmetric secant (MSS) methods can be thought of as generalizations of quasi-Newton methods in that they attempt to impose additional requirements on their approximation of the Hessian. Given an initial Hessian approximation, MSS methods generate a sequence of possibly-indefinite matrices using rank-2 updates to solve nonconvex unconstrained optimization problems. For practical reasons, up to now, the initialization has been a constant multiple of the identity matrix. In this paper, we propose a new limited-memory MSS method for large-scale nonconvex optimization that allows for dense initializations. Numerical results on the CUTEst test problems suggest that the MSS method using a dense initialization outperforms the standard initialization. Numerical results also suggest that this approach is competitive with both a basic L-SR1 trust-region method and an L-PSB method.


title: Canonical bases arising from $\imath$quantum covering groups

abstract: For $\imath$quantum covering groups $(\mathbf{U}, \mathbf{U}^\imath)$ of super Kac-Moody type, we construct $\imath$-canonical bases for the highest weight integrable $\mathbf{U}$-modules and their tensor products regarded as $\mathbf{U}^\imath$-modules, as well as a canonical basis for the modified form of the $\imath$quantum group $\mathbf{U}^\imath$, using the $\imath^\pi$-divided powers, rank one canonical basis for $\mathbf{U}^\imath$.


title: Observational Signatures of Planets in Protoplanetary Disks: Temperature structures in spiral arms

abstract: High-resolution imaging of protoplanetary disks has unveiled a rich diversity of spiral structure, some of which may arise from disk-planet interaction. Using 3D hydrodynamics with $\beta$-cooling to a vertically-stratified background, as well as radiative-transfer modeling, we investigate the temperature rise in planet-driven spirals. In rapidly cooling disks, the temperature rise is dominated by a contribution from stellar irradiation, 0.3-3% inside the planet radius but always <0.5% outside. When cooling time equals or exceeds dynamical time, however, this is overwhelmed by hydrodynamic PdV work, which introduces a 10-20% perturbation within a factor of 2 from the planet's orbital radius. We devise an empirical fit of the spiral amplitude $\Delta (T)$ to take into account both effects. Where cooling is slow, we find also that temperature perturbations from buoyancy spirals -- a strictly 3D, non-isothermal phenomenon -- become nearly as strong as those from Lindblad spirals, which are amenable to 2D and isothermal studies. Our findings may help explain observed thermal features in disks like TW Hydrae and CQ Tauri, and underscore that 3D effects have a qualitatively important effect on disk structure.


title: Expansion of harmonic functions near the boundary of Dini domains

abstract: Let $u$ be a harmonic function in a $C^1$-Dini domain, such that $u$ vanishes on an open set of the boundary. We show that near every point in the open set, $u$ can be written uniquely as the sum of a non-trivial homogeneous harmonic polynomial and an error term of higher degree (depending on the Dini parameter). In particular, this implies that $u$ has a unique tangent function at every such point, and that the convergence rate to the tangent function can be estimated. We also study the relationship of tangent functions at nearby points in a special case.


title: Graphhopper: Multi-Hop Scene Graph Reasoning for Visual Question Answering

abstract: Visual Question Answering (VQA) is concerned with answering free-form questions about an image. Since it requires a deep semantic and linguistic understanding of the question and the ability to associate it with various objects that are present in the image, it is an ambitious task and requires multi-modal reasoning from both computer vision and natural language processing. We propose Graphhopper, a novel method that approaches the task by integrating knowledge graph reasoning, computer vision, and natural language processing techniques. Concretely, our method is based on performing context-driven, sequential reasoning based on the scene entities and their semantic and spatial relationships. As a first step, we derive a scene graph that describes the objects in the image, as well as their attributes and their mutual relationships. Subsequently, a reinforcement learning agent is trained to autonomously navigate in a multi-hop manner over the extracted scene graph to generate reasoning paths, which are the basis for deriving answers. We conduct an experimental study on the challenging dataset GQA, based on both manually curated and automatically generated scene graphs. Our results show that we keep up with a human performance on manually curated scene graphs. Moreover, we find that Graphhopper outperforms another state-of-the-art scene graph reasoning model on both manually curated and automatically generated scene graphs by a significant margin.


title: Supercritical percolation on graphs of polynomial growth

abstract: We consider Bernoulli percolation on transitive graphs of polynomial growth. In the subcritical regime ($p<p_c$), it is well known that the connection probabilities decay exponentially fast. In the present paper, we study the supercritical phase ($p>p_c$) and prove the exponential decay of the truncated connection probabilities (probabilities that two points are connected by an open path, but not to infinity). This sharpness result was established by [CCN87] on $\mathbb Z^d$ and uses the difficult slab result of Grimmett and Marstrand. However, the techniques used there are very specific to the hypercubic lattices and do not extend to more general geometries. In this paper, we develop new robust techniques based on the recent progress in the theory of sharp thresholds and the sprinkling method of Benjamini and Tassion. On $\mathbb Z^d$, these methods can be used to produce a new proof of the slab result of Grimmett and Marstrand.


title: Contextual Games: Multi-Agent Learning with Side Information

abstract: We formulate the novel class of contextual games, a type of repeated games driven by contextual information at each round. By means of kernel-based regularity assumptions, we model the correlation between different contexts and game outcomes and propose a novel online (meta) algorithm that exploits such correlations to minimize the contextual regret of individual players. We define game-theoretic notions of contextual Coarse Correlated Equilibria (c-CCE) and optimal contextual welfare for this new class of games and show that c-CCEs and optimal welfare can be approached whenever players' contextual regrets vanish. Finally, we empirically validate our results in a traffic routing experiment, where our algorithm leads to better performance and higher welfare compared to baselines that do not exploit the available contextual information or the correlations present in the game.


title: Searching for Rigidity in Algebraic Starscapes

abstract: We create plots of algebraic integers in the complex plane, exploring the effect of sizing the integers according to various arithmetic invariants. We focus on Galois theoretic invariants, in particular creating plots which emphasize algebraic integers whose Galois group is not the full symmetric group--these integers we call rigid. We then give some analysis of the resulting images, suggesting avenues for future research about the geometry of so-called rigid algebraic integers.


title: Efficient exact computation of the conjunctive and disjunctive decompositions of D-S Theory for information fusion: Translation and extension

abstract: Dempster-Shafer Theory (DST) generalizes Bayesian probability theory, offering useful additional information, but suffers from a high computational burden. A lot of work has been done to reduce the complexity of computations used in information fusion with Dempster's rule. Yet, few research had been conducted to reduce the complexity of computations for the conjunctive and disjunctive decompositions of evidence, which are at the core of other important methods of information fusion. In this paper, we propose a method designed to exploit the actual evidence (information) contained in these decompositions in order to compute them. It is based on a new notion that we call focal point, derived from the notion of focal set. With it, we are able to reduce these computations up to a linear complexity in the number of focal sets in some cases. In a broader perspective, our formulas have the potential to be tractable when the size of the frame of discernment exceeds a few dozen possible states, contrary to the existing litterature. This article extends (and translates) our work published at the french conference GRETSI in 2019.


title: Multi-variance replica exchange stochastic gradient MCMC for inverse and forward Bayesian physics-informed neural network

abstract: Physics-informed neural network (PINN) has been successfully applied in solving a variety of nonlinear non-convex forward and inverse problems. However, the training is challenging because of the non-convex loss functions and the multiple optima in the Bayesian inverse problem. In this work, we propose a multi-variance replica exchange stochastic gradient Langevin diffusion method to tackle the challenge of the multiple local optima in the optimization and the challenge of the multiple modal posterior distribution in the inverse problem. Replica exchange methods are capable of escaping from the local traps and accelerating the convergence. However, it may not be efficient to solve mathematical inversion problems by using the vanilla replica method directly since the method doubles the computational cost in evaluating the forward solvers (likelihood functions) in the two chains. To address this issue, we propose to make different assumptions on the energy function estimation and this facilities one to use solvers of different fidelities in the likelihood function evaluation. We give an unbiased estimate of the swapping rate and give an estimation of the discretization error of the scheme. To verify our idea, we design and solve four inverse problems which have multiple modes. The proposed method is also employed to train the Bayesian PINN to solve the forward and inverse problems; faster and more accurate convergence has been observed when compared to the stochastic gradient Langevin diffusion (SGLD) method and vanila replica exchange methods.


title: The Unintended Consequences of Minimizing the Price of Anarchy in Congestion Games

abstract: This work focuses on the design of taxes in atomic congestion games, a commonly studied model for competitive resource sharing. While most related studies focus on optimizing either the worst- or best-case performance (i.e., Price of Anarchy (PoA) or Price of Stability (PoS)), we investigate whether optimizing for the PoA has consequences on the PoS. Perhaps surprisingly, our results reveal a fundamental trade-off between the two performance metrics. Our main result demonstrates that the taxation rule that optimizes the PoA inherits a matching PoS, implying that the best outcome is no better than the worst outcome under such a design choice. We then study this trade-off in terms of the Pareto frontier between the PoA and PoS. Our results also establish that any taxes with PoS equal to 1 incur a much higher PoA, and that, in several well-studied cases, the untaxed setting lies strictly above the Pareto frontier.


title: Measurement of Smartphone Sensor Efficiency to Cosmic Ray Muons

abstract: A measurement of the efficiency of CMOS sensors in smartphone cameras to cosmic ray muons is presented. A coincidence in external scintillators indicates the passage of a cosmic ray muon, allowing the measurement of the efficiency of the CMOS sensor. The observed flux is consistent with well-established values, and efficiencies are presented as a function of the number of photo-electrons collected from the CMOS silicon photodiode pixels. These efficiencies are vital to understanding the feasibility of large-scale smartphone networks operating as air-shower observatories.


title: Broadband variability study of MAXI J1631-479 in its Hard-Intermediate State observed with Insight-HXMT

abstract: We report the energy-resolved broadband timing analysis of the black hole X-ray transient MAXI J1631-479 during its 2019 outburst from February 11 to April 9, using data from the Insight-Hard X-ray Modulation Telescope (Insight-HXMT), which caught the source from its hard intermediate state to the soft state. Thanks to the large effective area of Insight-HXMT at high energies, we are able to present the energy dependence of fast variability up to ~100 keV. Type-C quasi-periodic oscillations (QPOs) with frequency varying between 4.9 Hz and 6.5 Hz are observed in the 1-100 keV energy band. While the QPO fractional rms increases with photon energy from 1 keV to ~10 keV and remains more or less constant from ~10 keV to ~100 keV, the rms of the flat-top noise first increases from 1 keV to ~8 keV then drops to less than 0.1% above ~30 keV. We suggest that the disappearance of the broadband variability above 30 keV could be caused by the non-thermal acceleration in the Comptonizing plasma. At the same time, the QPOs could be produced by the precession of either a small-scale jet or a hot inner flow model.


title: A Diffuse Interface Model for Cell Blebbing Including Membrane-Cortex Coupling with Linker Dynamics

abstract: The aim of this paper is to develop suitable models for the phenomenon of cell blebbing, which allow for computational predictions of mechanical effects including the crucial interaction of the cell membrane and the actin cortex. For this sake we resort to a two phase-field model that uses diffuse descriptions of both the membrane and the cortex, which in particular allows for a suitable description of the interaction via linker protein densities. Besides the detailed modelling we discuss some energetic aspects of the models and present a numerical scheme, which allows to carry out several computational studies. In those we demonstrate that several effects found in experiments can be reproduced, in particular bleb formation by cortex rupture, which was not possible by previous models without the linker dynamics.


title: Bosonization of Majorana modes and edge states

abstract: We present a bosonization procedure which replaces fermions with generalized spin variables subject to local constraints. It requires that the number of Majorana modes per lattice site matches the coordination number modulo two. If this condition is not obeyed, then bosonization introduces additional fermionic excitations not present in the original model. In the case of one Majorana mode per site on a honeycomb lattice, we recover a sector of Kitaev's model. We discuss also decagonal and rectangular geometries and present bosonization of the Hubbard model. For geometries with a boundary we find that certain fermionic edge modes naturally emerge. They are of different nature than edge modes encountered in topological phases of matter. Euclidean representation for the unconstrained version of a spin system of the type arising in our construction is derived and briefly studied by computing some exact averages for small volumes.


title: Improving Cooperative Game Theory-based Data Valuation via Data Utility Learning

abstract: The Shapley value (SV) and Least core (LC) are classic methods in cooperative game theory for cost/profit sharing problems. Both methods have recently been proposed as a principled solution for data valuation tasks, i.e., quantifying the contribution of individual datum in machine learning. However, both SV and LC suffer computational challenges due to the need for retraining models on combinatorially many data subsets. In this work, we propose to boost the efficiency in computing Shapley value or Least core by learning to estimate the performance of a learning algorithm on unseen data combinations. Theoretically, we derive bounds relating the error in the predicted learning performance to the approximation error in SV and LC. Empirically, we show that the proposed method can significantly improve the accuracy of SV and LC estimation.


title: Anomalous measurements: recent results deepen flavour puzzle

abstract: The crown jewel of particle physics, the Standard Model (SM), has withstood numerous experimental trials. However, there are still some observations it cannot explain. Examples such as dark matter and the matter-antimatter imbalance in the Universe come to mind. The SM may be extended, by including additional particles and interactions, so as to explain such phenomena. These new particles and interactions are collectively referred to as "new physics" (NP), and the results covered by this article provide a promising lead for their discovery.


title: Spectral Recovery of Binary Censored Block Models

abstract: Community detection is the problem of identifying community structure in graphs. Often the graph is modeled as a sample from the Stochastic Block Model, in which each vertex belongs to a community. The probability that two vertices are connected by an edge depends on the communities of those vertices. In this paper, we consider a model of {\em censored} community detection with two communities, where most of the data is missing as the status of only a small fraction of the potential edges is revealed. In this model, vertices in the same community are connected with probability $p$ while vertices in opposite communities are connected with probability $q$. The connectivity status of a given pair of vertices $\{u,v\}$ is revealed with probability $\alpha$, independently across all pairs, where $\alpha = \frac{t \log(n)}{n}$. We establish the information-theoretic threshold $t_c(p,q)$, such that no algorithm succeeds in recovering the communities exactly when $t < t_c(p,q)$. We show that when $t > t_c(p,q)$, a simple spectral algorithm based on a weighted, signed adjacency matrix succeeds in recovering the communities exactly. While spectral algorithms are shown to have near-optimal performance in the symmetric case, we show that they may fail in the asymmetric case where the connection probabilities inside the two communities are allowed to be different. In particular, we show the existence of a parameter regime where a simple two-phase algorithm succeeds but any algorithm based on the top two eigenvectors of the weighted, signed adjacency matrix fails.


title: Generation of Photon Pairs by Stimulated Emission in Ring Resonators

abstract: Third-order parametric down-conversion (TOPDC) describes a class of nonlinear interactions in which a pump photon is converted into a photon triplet. This process can occur spontaneously, or it can be stimulated by seeding fields. In the former case, one typically has the generation of non-Gaussian states of light. In the latter, the situation is more variegated, for stimulated TOPDC (StTOPDC) can be implemented in many ways, depending on the number and properties of the seeding fields. Here we show that StTOPDC can be exploited for the generation of quantum correlated photon pairs. We examine the peculiar features of this approach when compared with second-order spontaneous parametric down-conversion and spontaneous four-wave mixing. We model StTOPDC in a microring resonator, predicting observable generation rates in a microring engineered for third-harmonic generation. We conclude that if the experimental difficulties associated with implementing StTOPDC can be overcome, it may soon be possible to demonstrate this process in resonant integrated devices.


title: EC 22536-5304: a lead-rich and metal-poor long-period binary

abstract: Helium-burning hot subdwarf stars of spectral types O and B (sdO/B) are thought to be produced through various types of binary interaction. The helium-rich hot subdwarf star EC22536-5304 was recently found to be extremely enriched in lead. Here, we show that EC22536-5304 is a binary star with a metal-poor subdwarf F-type (sdF) companion. We have performed a detailed analysis of high-resolution SALT/HRS and VLT/UVES spectra, deriving metal abundances for the hot subdwarf, as well as atmospheric parameters for both components. Because we consider the contribution of the sdF star, the derived lead abundance for the sdOB, +6.3 +/- 0.3 dex relative to solar, is even higher than previously thought. We derive Teff = 6210 +/- 70 K, log g = 4.64 +/- 0.10, [Fe/H] = -1.95 +/- 0.04, and [$\alpha$/Fe] = +0.40 +/- 0.04 for the sdF component. Radial velocity variations, although poorly sampled at present, indicate that the binary system has a long orbital period of about 457 days. This suggests that the system was likely formed through stable Roche lobe overflow (RLOF). A kinematic analysis shows that EC22536-5304 is on an eccentric orbit around the Galactic centre. This, as well as the low metallicity and strong alpha-enhancement of the sdF-type companion, indicate that EC22536-5304 is part of the Galactic halo or metal-weak thick disc. As the first long-period hot subdwarf binary below about [Fe/H] = -1, EC22536-5304 may help to constrain the RLOF mechanism for mass transfer from low-mass, low-metallicity red giant branch (RGB) stars to main-sequence companions.


title: Hybrid A Posteriori Error Estimators for Conforming Finite Element Approximations to Stationary Convection-Diffusion-Reaction equations

abstract: We consider the a posteriori error estimation for convection-diffusion-reaction equations in both diffusion-dominated and convection/reaction-dominated regimes. We present an explicit hybrid estimator, which, in each regime, is proved to be reliable and efficient with constants independent of the parameters in the underlying problem. For convection-dominated problems, the norm introduced by Verf{\"u}rth \cite{verf2005confusion} is used to measure the approximation error. Various numerical experiments are performed to (1) demonstrate the robustness of the hybrid estimator; (2) show that the hybrid estimator is more accurate than the explicit residual estimator and is less sensitive to the size of reaction, even though both of them are robust.


title: Design, simulation and characterization of integrated photonic spectrographs for Astronomy I: Generation-I AWG devices based on canonical layouts

abstract: We present an experimental study on our first generation of custom-developed arrayed waveguide gratings (AWG) on silica platform for spectroscopic applications in near-infrared astronomy. We provide a comprehensive description of the design, numerical simulation and characterization of several AWG devices aimed at spectral resolving powers of 15,000 - 60,000 in the astronomical H-band. We evaluate the spectral characteristics of the fabricated devices in terms of insertion loss and estimated spectral resolving power and compare the results with numerical simulations. We estimate resolving powers of up to 18,900 from the output channel 3-dB transmission bandwidth. Based on the first characterization results, we select two candidate AWGs for further processing by removal of the output waveguide array and polishing the output facet to optical quality with the goal of integration as the primary diffractive element in a cross-dispersed spectrograph. We further study the imaging properties of the processed AWGs with regards to spectral resolution in direct imaging mode, geometry-related defocus aberration, and polarization sensitivity of the spectral image. We identify phase error control, birefringence control, and aberration suppression as the three key areas of future research and development in the field of high-resolution AWG-based spectroscopy in astronomy.


title: Design of Adaptive Backstepping Control for Direct Power Control of Three-Phase PWM Rectifier)

abstract: In this paper, we focused on the design of an adaptive backstepping controller (adaptive-BSC) for direct power control (DPC) of a three-phase PWM rectifier. In the proposed system, it is desired to control both the output DC voltage of the rectifier and the reactive power simultaneously by making them track desired respective values. This was done by having independent virtual control signals for both the output voltage and the reactive power. The adaptive control signals were gotten from the dynamic equations of the three-phase system. For comparison, both the BSC and adaptive-BSC equations were developed. Numerical simulations were performed on both of them on a 5kW system. The proposed adaptive-BSC was designed to work under more challenging system variations as compared with the BSC as it has to estimate the value of the unknown system load. Despite this, it still performed better than its BSC counterpart.


title: Inverse Reinforcement Learning Based Stochastic Driver Behavior Learning

abstract: Drivers have unique and rich driving behaviors when operating vehicles in traffic. This paper presents a novel driver behavior learning approach that captures the uniqueness and richness of human driver behavior in realistic driving scenarios. A stochastic inverse reinforcement learning (SIRL) approach is proposed to learn a distribution of cost function, which represents the richness of the human driver behavior with a given set of driver-specific demonstrations. Evaluations are conducted on the realistic driving data collected from the 3D driver-in-the-loop driving simulation. The results show that the learned stochastic driver model is capable of expressing the richness of the human driving strategies under different realistic driving scenarios. Compared to the deterministic baseline driver behavior model, the results reveal that the proposed stochastic driver behavior model can better replicate the driver's unique and rich driving strategies in a variety of traffic conditions.


title: Extended L-ensembles: a new representation for Determinantal Point Processes

abstract: Determinantal point processes (DPPs) are a class of repulsive point processes, popular for their relative simplicity. They are traditionally defined via their marginal distributions, but a subset of DPPs called "L-ensembles" have tractable likelihoods and are thus particularly easy to work with. Indeed, in many applications, DPPs are more naturally defined based on the L-ensemble formulation rather than through the marginal kernel. The fact that not all DPPs are L-ensembles is unfortunate, but there is a unifying description. We introduce here extended L-ensembles, and show that all DPPs are extended L-ensembles (and vice-versa). Extended L-ensembles have very simple likelihood functions, contain L-ensembles and projection DPPs as special cases. From a theoretical standpoint, they fix some pathologies in the usual formalism of DPPs, for instance the fact that projection DPPs are not L-ensembles. From a practical standpoint, they extend the set of kernel functions that may be used to define DPPs: we show that conditional positive definite kernels are good candidates for defining DPPs, including DPPs that need no spatial scale parameter. Finally, extended L-ensembles are based on so-called ``saddle-point matrices'', and we prove an extension of the Cauchy-Binet theorem for such matrices that may be of independent interest.


title: Driven dissipative preparation of few-body Laughlin states of Rydberg polaritons in twisted cavities

abstract: We present a driven dissipative protocol for creating an optical analog of the Laughlin state in a system of Rydberg polaritons in a twisted optical cavity. We envision resonantly driving the system into a 4-polariton state by injecting photons in carefully selected modes. The dissipative nature of the polariton-polariton interactions leads to a decay into a two-polariton analog of the Laughlin state. Generalizations of this technique could be used to explore fractional statistics and anyon based quantum information processing. We also model recent experiments that attempt to coherently drive into this same state.


title: Enhancing Resilience of Distribution Networks by Coordinating Microgrids and Demand Response Programs in Service Restoration

abstract: In case of high impact low probability events, in order to restore the critical loads of the distribution network as much as possible, it is necessary to employ all available resources such as microgrids and distributed generations. This paper presents a two-stage method for critical load restoration after an extreme event utilizing the coordination of all available distributed generations, microgrids, and demand response programs. In the first stage, the post-disaster reconfiguration of the network is determined, and in the second stage, the demand response outputs of the responsive loads and the restoration status of all critical loads are ascertained. In the first stage, an algorithm is proposed to determine electrical islands, and in the second stage, a new model based on the generalized Benders decomposition algorithm is proposed, which considers demand response programs and all operational constraints and aims to maximize the restored critical loads. The effectiveness of the proposed method is verified by simulating the 33-bus and 69-bus test systems.


title: Vertical compositional variations of liquid hydrocarbons in Titan's alkanofers

abstract: According to clues left by the Cassini mission, Titan, one of the two Solar System bodies with a hydrologic cycle, may harbor liquid hydrocarbon-based analogs of our terrestrial aquifers, referred to as "alkanofers". On the Earth, petroleum and natural gas reservoirs show a vertical gradient in chemical composition, established over geological timescales. In this work, we aim to investigate the conditions under which Titan's processes could lead to similar situations. We built numerical models including barodiffusion and thermodiffusion (Soret's effect) in N_2+CH_4+C_2H_6 liquid mixtures, which are relevant for Titan's possible alkanofers. Our main assumption is the existence of reservoirs of liquids trapped in a porous matrix with low permeability. Due to the small size of the molecule, nitrogen seems to be more sensitive to gravity than ethane, even if the latter has a slightly larger mass. This behavior, noticed for an isothermal crust, is reinforced by the presence of a geothermal gradient. Vertical composition gradients, formed over timescales of between a fraction of a mega-year to several tens of mega-years, are not influenced by molecular diffusion coefficients. We find that ethane does not accumulate at the bottom of the alkanofers under diffusion, leaving the question of why ethane is not observed on Titan's surface unresolved. If the alkanofer liquid was in contact with water-ice, we checked that N_2 did not, in general, impede the clathration of C_2H_6, except in some layers. Interestingly, we found that noble gases could easily accumulate at the bottom of an alkanofer.


title: Arbitrage-free pricing of CVA for cross-currency swap with wrong-way risk under stochastic correlation modeling framework

abstract: A positive correlation between exposure and counterparty credit risk gives rise to the so-called Wrong-Way Risk (WWR). Even after a decade of the financial crisis, addressing WWR in both sound and tractable ways remains challenging. Academicians have proposed arbitrage-free set-ups through copula methods but those are computationally expensive and hard to use in practice. Resampling methods are proposed by the industry but they lack mathematical foundations. The purpose of this article is to bridge this gap between the approaches used by academicians and industry. To this end, we propose a stochastic correlation approach to asses WWR. The methods based on constant correlation to model the dependency between exposure and counterparty credit risk assume a linear dependency, thus fail to capture the tail dependence. Using a stochastic correlation we move further away from the Gaussian copula and can capture the tail risk. This effect is reflected in the results where the impact of stochastic correlation on calculated CVA is substantial when compared to the case when a high constant correlation is assumed between exposure and credit. Given the uncertainty inherent to CVA, the proposed method is believed to provide a promising way to model WWR.


title: Shift-Invariance of the Colored TASEP and Finishing Times of the Oriented Swap Process

abstract: We prove a new shift-invariance property of the colored TASEP. From the shift-invariance of the colored stochastic six-vertex model (proved in Borodin-Gorin-Wheeler or Galashin), one can get a shift-invariance property of the colored TASEP at one time, and our result generalizes this to multiple times. Our proof takes the single-time shift-invariance as an input, and uses analyticity of the probability functions and induction arguments. We apply our shift-invariance to prove a distributional identity between the finishing times of the oriented swap process and the point-to-line passage times in exponential last-passage percolation, which is conjectured by Bisi-Cunden-Gibbons-Romik and Bufetov-Gorin-Romik, and is also equivalent to a purely combinatorial identity related to the Edelman-Greene correspondence. With known results from last-passage percolation, we also get new asymptotic results on the colored TASEP and the finishing times of the oriented swap process.


title: BRIMA: low-overhead BRowser-only IMage Annotation tool (Preprint)

abstract: Image annotation and large annotated datasets are crucial parts within the Computer Vision and Artificial Intelligence fields.At the same time, it is well-known and acknowledged by the research community that the image annotation process is challenging, time-consuming and hard to scale. Therefore, the researchers and practitioners are always seeking ways to perform the annotations easier, faster, and at higher quality. Even though several widely used tools exist and the tools' landscape evolved considerably, most of the tools still require intricate technical setups and high levels of technical savviness from its operators and crowdsource contributors. In order to address such challenges, we develop and present BRIMA -- a flexible and open-source browser extension that allows BRowser-only IMage Annotation at considerably lower overheads. Once added to the browser, it instantly allows the user to annotate images easily and efficiently directly from the browser without any installation or setup on the client-side. It also features cross-browser and cross-platform functionality thus presenting itself as a neat tool for researchers within the Computer Vision, Artificial Intelligence, and privacy-related fields.


title: Pairing symmetries in the Zeeman-coupled extended attractive Hubbard model

abstract: By introducing the possibility of equal- and opposite-spin pairings concurrently, we show that the extended attractive Hubbard model (EAHM) exhibits rich ground state phase diagrams with a variety of singlet, triplet, and mixed parity superconducting orders. We study the competition between these superconducting pairing symmetries invoking an unrestricted Hartree-Fock- Bogoliubov-de Gennes (HFBdG) mean-field approach, and we use the d-vector formalism to characterize the nature of the stabilized superconducting orders. We discover that, while all other types of orders are suppressed, a non-unitary triplet order dominates the phase space in the presence of an in-plane external magnetic field. We also find a transition between a non-unitary to unitary superconducting phase driven by the change in average electron density. Our results serve as a reference for identifying and understanding the nature of superconductivity based on the symmetries of the pairing correlations. The results further highlight that EAHM is a suitable effective model for describing most of the pairing symmetries discovered in different materials.


title: Distributionally Robust Policy Learning via Adversarial Environment Generation

abstract: Our goal is to train control policies that generalize well to unseen environments. Inspired by the Distributionally Robust Optimization (DRO) framework, we propose DRAGEN - Distributionally Robust policy learning via Adversarial Generation of ENvironments - for iteratively improving robustness of policies to realistic distribution shifts by generating adversarial environments. The key idea is to learn a generative model for environments whose latent variables capture cost-predictive and realistic variations in environments. We perform DRO with respect to a Wasserstein ball around the empirical distribution of environments by generating realistic adversarial environments via gradient ascent on the latent space. We demonstrate strong Out-of-Distribution (OoD) generalization in simulation for (i) swinging up a pendulum with onboard vision and (ii) grasping realistic 3D objects. Grasping experiments on hardware demonstrate better sim2real performance compared to domain randomization.


title: Mid-infrared Single-photon Detection Using Superconducting NbTiN Nanowires with Sub-15 ps Time Resolution in a Gifford-McMahon Cryocooler

abstract: Shortly after their inception, superconducting nanowire single-photon detectors (SNSPDs) became the leading quantum light detection technology. With the capability of detecting single-photons with near-unity efficiency, high time resolution, low dark count rate, and fast recovery time, SNSPDs outperform conventional single-photon detection techniques. However, detecting lower energy single-photons (<0.8 eV) with high efficiency and low timing jitter has remained a challenge. To achieve unity internal efficiency at mid-infrared wavelengths, previous works used amorphous superconducting materials with low energy gaps at the expense of reduced time resolution (close to a nanosecond), and by operating them in complex mK dilution refrigerators. In this work, we provide an alternative approach with SNSPDs fabricated from 5-9.5 nm thick NbTiN superconducting films and devices operated in conventional Gifford-McMahon (GM) cryocoolers. By optimizing the superconducting film deposition process, film thickness and nanowire design, our fiber-coupled devices achieved > 70% system detection efficiency (SDE) at 2 {\mu}m and sub-15 ps timing jitter. Furthermore, detectors from the same batch demonstrated unity internal detection efficiency at 3 {\mu}m and 80% internal efficiency at 4 {\mu}m, paving the road for an efficient mid-infrared single-photon detection technology with unparalleled time resolution and without mK cooling requirements.


title: Shaping the equation of state to improve numerical accuracy and stability of the pseudopotential lattice Boltzmann method

abstract: Recently it was discovered that altering the shape of the meta stable and unstable branches of an equation of state (EOS) can greatly improve the numerical accuracy of liquid and gas densities in the pseudopotential method. Inspired by this approach we develop an improved approach that is benchmarked for both equilibrium and non-equilibrium situations. We show here that the original approach reduces the method stability in non-equilibrium situations. Here we propose a new procedure to replace the metastable and unstable regions of these EOS by alternative functions. Our approach does not affects the coexistence densities or the speed of sound of the liquid phase while maintaining continuity of the sound speed in the pressure-density curve. Using this approach we were able to reduce the relative error of the planar interface vapor density compared to the thermodynamic consistent value by increasing the vapor phase sound speed. To allow for the benchmarking of dynamic results we also developed a finite difference method (FD) that solves the same macroscopic conservation equation as the pseudopotential lattice Boltzmann method (LBM). With this FD scheme we are able to perform mesh refinement and obtain reference solutions for the dynamic tests. We observed excellent agreement between the FD solutions and our proposed scheme. We also performed a detailed study of the stability of the methods using simulations of a droplet impacting on a liquid film for reduced temperatures down to 0.35 with Reynolds number of 300. Our approach remains stable for a density ratio up to $3.38\cdot10^{4}$.


title: Real-Time Pothole Detection Using Deep Learning

abstract: Roads are connecting line between different places, and used daily. Roads' periodic maintenance keeps them safe and functional. Detecting and reporting the existence of potholes to responsible departments can help in eliminating them. This study deployed and tested on different deep learning architecture to detect potholes. The images used for training were collected by cellphone mounted on the windshield of the car, in addition to many images downloaded from the internet to increase the size and variability of the database. Second, various object detection algorithms are employed and compared to detect potholes in real-time like SDD-TensorFlow, YOLOv3Darknet53 and YOLOv4Darknet53. YOLOv4 achieved the best performance with 81% recall, 85% precision and 85.39% mean Average Precision (mAP). The speed of processing was 20 frame per second. The system was able to detect potholes from a range on 100 meters away from the camera. The system can increase the safety of drivers and improve the performance of self-driving cars by detecting pothole time ahead.


title: Purcell radiative rate enhancement of label-free proteins with ultraviolet aluminum plasmonics

abstract: The vast majority of proteins are intrinsically fluorescent in the ultraviolet, thanks to the emission from their tryptophan and tyrosine amino-acid constituents. However, the protein autofluorescence quantum yields are generally very low due to the prevailing quenching mechanisms by other amino acids inside the protein. This motivates the interest to enhance the radiative emission rate of proteins using nanophotonic structures. Although there have been numerous reports of Purcell effect and local density of optical states (LDOS) control in the visible range using single dipole quantum emitters, the question remains open to apply these concepts in the UV on real proteins containing several tryptophan and tyrosine amino acids arranged in a highly complex manner. Here, we report the first complete characterization of the Purcell effect and radiative rate enhancement for the UV intrinsic fluorescence of label-free \b{eta}-galactosidase and streptavidin proteins in plasmonic aluminum nanoapertures. We find an excellent agreement with a calibration performed using a high quantum yield UV fluorescent dye. Demonstrating and intensifying the Purcell effect is essential for the applications of UV plasmonics and the label-free detection of single proteins.


title: Ramification loci of non-archimedean cubic rational functions

abstract: For a cubic rational function with coefficients in a non-archimedean field $K$ whose residue characteristic is $0$ or greater than $3$, there are $2$ possibilities for the shape of its Berkovich ramification locus, considered as an endomorphism of the Berkovich projective line: one is the connected hull of all the critical points, and the other is consisting of $2$ disjoint segments. In this paper, we list up all the possible forms of cubic rational functions and calculate their ramification loci.


title: Strong bounds and exact solutions to the minimum broadcast time problem

abstract: Given a graph and a subset of its nodes, referred to as source nodes, the minimum broadcast problem asks for the minimum number of steps in which a signal can be transmitted from the sources to all other nodes in the graph. In each step, the sources and the nodes that already have received the signal can forward it to at most one of their neighbour nodes. The problem has previously been proved to be NP-hard. In the current work, we develop a compact integer programming model for the problem. We also devise procedures for computing lower bounds on the minimum number of steps required, along with methods for constructing near-optimal solutions. Computational experiments demonstrate that in a wide range of instances, in particular instances with sufficiently dense graphs, the lower and upper bounds under study collapse. In instances where this is not the case, the integer programming model proves strong capabilities in closing the remaining gap, and proves to be considerably more efficient than previously studied models.


title: Developmental Stage Classification of Embryos Using Two-Stream Neural Network with Linear-Chain Conditional Random Field

abstract: The developmental process of embryos follows a monotonic order. An embryo can progressively cleave from one cell to multiple cells and finally transform to morula and blastocyst. For time-lapse videos of embryos, most existing developmental stage classification methods conduct per-frame predictions using an image frame at each time step. However, classification using only images suffers from overlapping between cells and imbalance between stages. Temporal information can be valuable in addressing this problem by capturing movements between neighboring frames. In this work, we propose a two-stream model for developmental stage classification. Unlike previous methods, our two-stream model accepts both temporal and image information. We develop a linear-chain conditional random field (CRF) on top of neural network features extracted from the temporal and image streams to make use of both modalities. The linear-chain CRF formulation enables tractable training of global sequential models over multiple frames while also making it possible to inject monotonic development order constraints into the learning process explicitly. We demonstrate our algorithm on two time-lapse embryo video datasets: i) mouse and ii) human embryo datasets. Our method achieves 98.1 % and 80.6 % for mouse and human embryo stage classification, respectively. Our approach will enable more profound clinical and biological studies and suggests a new direction for developmental stage classification by utilizing temporal information.


title: Quantum spins and hybridization in artificially-constructed chains of magnetic adatoms on a superconductor

abstract: Magnetic adatom chains on surfaces constitute fascinating quantum spin systems. Superconducting substrates suppress interactions with bulk electronic excitations but couple the adatom spins to a chain of subgap Yu-Shiba-Rusinov (YSR) quasiparticles. Using a scanning tunneling microscope, we investigate such correlated spin-fermion systems by constructing Fe chains adatom by adatom on superconducting NbSe$_2$. The adatoms couple entirely via the substrate, retaining their quantum spin nature. In dimers, we observe that the deepest YSR state undergoes a quantum phase transition due to Ruderman-Kittel-Kasuya-Yosida interactions, a distinct signature of quantum spins. Chains exhibit coherent hybridization and band formation of the YSR excitations, indicating ferromagnetic coupling. Longer chains develop separate domains due to coexisting charge-density-wave order of NbSe$_2$. Despite the spin-orbit-coupled substrate, we find no signatures of Majoranas, possibly because quantum spins reduce the parameter range for topological superconductivity. We suggest that adatom chains are versatile systems for investigating correlated-electron physics and its interplay with topological superconductivity.


title: Information Spread with Error Correction

abstract: We study the process of information dispersal in a network with communication errors and local error-correction. Specifically we consider a simple model where a single bit of information initially known to a single source is dispersed through the network, and communication errors lead to differences in the agents' opinions on this information. Naturally, such errors can very quickly make the communication completely unreliable, and in this work we study to what extent this unreliability can be mitigated by local error-correction, where nodes periodically correct their opinion based on the opinion of (some subset of) their neighbors. We analyze how the error spreads in the "early stages" of information dispersal by monitoring the average opinion, i.e., the fraction of agents that have the correct information among all nodes that hold an opinion at a given time. Our main results show that even with significant effort in error-correction, tiny amounts of noise can lead the average opinion to be nearly uncorrelated with the truth in early stages. We also propose some local methods to help agents gauge when the information they have has stabilized.


title: Lattices in Tate modules

abstract: Refining a theorem of Zarhin, we prove that given a $g$-dimensional abelian variety $X$ and an endomorphism $u$ of $X$, there exists a matrix $A \in \operatorname{M}_{2g}(\mathbb{Z})$ such that each Tate module $T_\ell X$ has a $\mathbb{Z}_\ell$-basis on which the action of $u$ is given by $A$.


title: Scaling Theory of Few-Particle Delocalization

abstract: We develop a scaling theory of interaction-induced delocalization of few-particle states in disordered quantum systems. In the absence of interactions, all single-particle states are localized in $d<3$, while in $d \geq 3$ there is a critical disorder below which states are delocalized. We hypothesize that such a delocalization transition occurs for $n$-particle bound states in $d$ dimensions when $d+n\geq 4$. Exact calculations of disorder-averaged $n$-particle Greens functions support our hypothesis. In particular, we show that $3$-particle states in $d=1$ with nearest-neighbor repulsion will delocalize with $W_c \approx 1.4t$ and with localization length critical exponent $\nu = 1.5 \pm 0.3$. The delocalization transition can be understood by means of a mapping onto a non-interacting problem with symplectic symmetry. We discuss the importance of this result for many-body delocalization, and how few-body delocalization can be probed in cold atom experiments.


title: The yielding of amorphous solids at finite temperatures

abstract: We analyze the effect of temperature on the yielding transition of amorphous solids using different coarse-grained model approaches. On one hand we use an elasto-plastic model, with temperature introduced in the form of an Arrhenius activation law over energy barriers. On the other hand, we implement a Hamiltonian model with a relaxational dynamics, where temperature is introduced in the form of a Langevin stochastic force. In both cases, temperature transforms the sharp transition of the athermal case in a smooth crossover. We show that this thermally smoothed transition follows a simple scaling form that can be fully explained using a one-particle system driven in a potential under the combined action of a mechanical and a thermal noise, the stochastically-driven Prandtl-Tomlinson model. Our work harmonizes the results of simple models for amorphous solids with the phenomenological $\sim T^{2/3}$ law proposed by Johnson and Samwer [Phys. Rev. Lett. 95, 195501 (2005)] in the framework of experimental metallic glasses yield observations, and extend it to a generic case. Finally, our results strengthen the interpretation of the yielding transition as an effective mean-field phenomenon.


title: Dynamical Surface Imaging of lambda Andromedae

abstract: We present temperature maps of RS CVn star lambda Andromedae, reconstructed from interferometric data acquired in 2010 and 2011 by the MIRC instrument at the Center for High Angular Resolution Astronomy Array. To constrain the stellar parameters required for this imaging task, we first modeled the star using our GPU-accelerated code SIMTOI. The stellar surface was then imaged using our open source interferometric imaging code ROTIR, in the process further refining the estimation of stellar parameters. We report that the measured angular diameter is 2.742 +/- 0.010 mas with a limb-darkening coefficient of 0.231 +/- 0.024. While our images are consistent with those of prior works, we provide updated physical parameters for lambda Andromedae (R_star = 7.78 +/- 0.05 R_odot, M_star = 1.24 +/- 0.72 M_odot, log L/L_odot = 1.46 +/- 0.04).


title: Optical Charge Injection and Full Coherent Control of Spin-Qubit in the Telecom C-band Emitting Quantum Dot

abstract: Solid-state quantum emitters with manipulable spin-qubits are promising platforms for quantum communication applications. Although such light-matter interfaces could be realized in many systems only a few allow for light emission in the telecom bands necessary for long-distance quantum networks. Here, we propose and implement a new optically active solid-state spin-qubit based on a hole confined in a single InAs/GaAs quantum dot grown on an InGaAs metamorphic buffer layer emitting photons in the C-band. We lift the hole spin-degeneracy using an external magnetic field and demonstrate hole injection, initialization, read-out and complete coherent control using picosecond optical pulses. These results showcase a new solid-state spin-qubit platform compatible with preexisting optical fibre networks.


title: Fault-tolerant parity readout on a shuttling-based trapped-ion quantum computer

abstract: Quantum error correction requires the detection of errors by reliable measurements of suitable multi-qubit correlation operators. Here, we experimentally demonstrate a fault-tolerant weight-4 parity check measurement scheme. An additional 'flag' qubit serves to detect errors occurring throughout the parity measurement, which would otherwise proliferate into uncorrectable weight-2 errors on the qubit register. We achieve a flag-conditioned parity measurement single-shot fidelity of 93.2(2)\%. Deliberately injecting bit and phase-flip errors, we show that the fault-tolerant protocol is capable of reliably intercepting such faults. For holistic benchmarking of the parity measurement scheme, we use entanglement witnessing to show that the implemented circuit generates genuine six-qubit multi-partite entanglement. The fault-tolerant parity measurement scheme is an essential building block in a broad class of stabilizer quantum error correction protocols, including topological color codes. Our hardware platform is based on atomic ions stored in a segmented microchip ion trap. The qubit register is dynamically reconfigured via shuttling operations, enabling effective full connectivity without operational cross-talk, which provides key capabilities for scalable fault-tolerant quantum computing.


title: Exploring DMD-type Algorithms for Modeling Signalised Intersections

abstract: This paper explores a novel data-driven approach based on recent developments in Koopman operator theory and dynamic mode decomposition (DMD) for modeling signalized intersections. Vehicular flow and queue formation on signalized intersections have complex nonlinear dynamics, making system identification, modeling, and controller design tasks challenging. We employ a Koopman theoretic approach to transform the original nonlinear dynamics into locally linear infinite-dimensional dynamics. The data-driven approach relies entirely on spatio-temporal snapshots of the traffic data. We investigate several key aspects of the approach and provide insights into the usage of DMD-type algorithms for application in adaptive signalized intersections. To demonstrate the utility of the obtained linearized dynamics, we perform prediction of the queue lengths at the intersection; and compare the results with the state-of-the-art long short term memory (LSTM) method. The case study involves the morning peak vehicle movements and queue lengths at two Orlando area signalized intersections. It is observed that DMD-based algorithms are able to capture complex dynamics with a linear approximation to a reasonable extent.


title: Tame and relatively elliptic $\mathbb{CP}^1$-structures on the thrice-punctured sphere

abstract: Suppose a relatively elliptic representation $\rho$ of the fundamental group of the thrice-punctured sphere $S$ is given. We prove that all projective structures on $S$ with holonomy $\rho$ and satisfying a tameness condition at the punctures can be obtained by grafting certain circular triangles. The specific collection of triangles is determined by a natural framing of $\rho$. In the process, we show that (on a general surface $\Sigma$ of negative Euler characteristics) structures satisfying these conditions can be characterized in terms of their M\"obius completion, and in terms of certain meromorphic quadratic differentials.


title: Intersection Problems in Extremal Combinatorics: Theorems, Techniques and Questions Old and New

abstract: The study of intersection problems in Extremal Combinatorics dates back perhaps to 1938, when Paul Erd\H{o}s, Chao Ko and Richard Rado proved the (first) `Erd\H{o}s-Ko-Rado theorem' on the maximum possible size of an intersecting family of $k$-element subsets of a finite set. Since then, a plethora of results of a similar flavour have been proved, for a range of different mathematical structures, using a wide variety of different methods. Structures studied in this context have included families of vector subspaces, families of graphs, subsets of finite groups with given group actions, and of course uniform hypergraphs with stronger or weaker intersection conditions imposed. The methods used have included purely combinatorial ones such as shifting/compressions, algebraic methods (including linear-algebraic, Fourier analytic and representation-theoretic), and more recently, analytic, probabilistic and regularity-type methods. As well as being natural problems in their own right, intersection problems have connections with many other parts of Combinatorics and with Theoretical Computer Science (and indeed with many other parts of Mathematics), both through the results themselves, and the methods used. In this survey paper, we discuss both old and new results (and both old and new methods), in the field of intersection problems. Many interesting open problems remain; we will discuss several. For expositional and pedagogical purposes, we also take this opportunity to give slightly streamlined versions of proofs (due to others) of several classical results in the area. This survey is intended to be useful to PhD students, as well as to more established researchers. It is a personal perspective on the field, and is not intended to be exhaustive; we apologise for any omissions. It is an expanded version of a paper that will appear in the Proceedings of the 29th British Combinatorial Conference.


title: On the Analysis of MUD-Files' Interactions, Conflicts, and Configuration Requirements Before Deployment

abstract: Manufacturer Usage Description (MUD) is an Internet Engineering Task Force (IETF) standard designed to protect IoT devices and networks by creating an out-of-the-box access control list for an IoT device. %The protocol defines a conceptually straightforward method to implement an isolation-based defensive mechanism based on the rules that are introduced by the manufacturer of the device. However, in practice, the access control list of each device is defined in its MUD-File and may contain possibly hundreds of access control rules. As a result, reading and validating these files is a challenge; and determining how multiple IoT devices interact is difficult for the developer and infeasible for the consumer. To address this we introduce the MUD-Visualizer to provide a visualization of any number of MUD-Files. MUD-Visualizer is designed to enable developers to produce correct MUD-Files by providing format correction, integrating them with other MUD-Files, and identifying conflicts through visualization. MUD-Visualizer is scalable and its core task is to merge and illustrate ACEs for multiple devices; both within and beyond the local area network. MUD-Visualizer is made publicly available and can be found on GitHub.


title: The Initial-Final Mass Relation for Hydrogen-Deficient White Dwarfs

abstract: The initial-final mass relation (IFMR) represents the total mass lost by a star during the entirety of its evolution from the zero age main sequence to the white dwarf cooling track. The semi-empirical IFMR is largely based on observations of DA white dwarfs, the most common spectral type of white dwarf and the simplest atmosphere to model. We present a first derivation of the semi-empirical IFMR for hydrogen deficient white dwarfs (non-DA) in open star clusters. We identify a possible discrepancy between the DA and non-DA IFMRs, with non-DA white dwarfs $\approx 0.07 M_\odot$ less massive at a given initial mass. Such a discrepancy is unexpected based on theoretical models of non-DA formation and observations of field white dwarf mass distributions. If real, the discrepancy is likely due to enhanced mass loss during the final thermal pulse and renewed post-AGB evolution of the star. However, we are dubious that the mass discrepancy is physical and instead is due to the small sample size, to systematic issues in model atmospheres of non-DAs, and to the uncertain evolutionary history of Procyon B (spectral type DQZ). A significantly larger sample size is needed to test these assertions. In addition, we also present Monte Carlo models of the correlated errors for DA and non-DA white dwarfs in the initial-final mass plane. We find the uncertainties in initial-final mass determinations for individual white dwarfs can be significantly asymmetric, but the recovered functional form of the IFMR is grossly unaffected by the correlated errors.


title: Bilinear Control of Convection-Cooling: From Open-Loop to Closed-Loop

abstract: This paper is concerned with a bilinear control problem for enhancing convection-cooling via an incompressible velocity field. Both optimal open-loop control and closed-loop feedback control designs are addressed. First and second order optimality conditions for characterizing the optimal solution are discussed. In particular, the method of instantaneous control is applied to establish the feedback laws. Moreover, the construction of feedback laws is also investigated by directly utilizing the optimality system with appropriate numerical discretization schemes. Computationally, it is much easier to implement the closed-loop feedback control than the optimal open-loop control, as the latter requires to solve the state equations forward in time, coupled with the adjoint equations backward in time together with a nonlinear optimality condition. Rigorous analysis and numerical experiments are presented to demonstrate our ideas and validate the efficacy of the control designs.


title: Exoplanets with ELT-METIS I: Estimating the direct imaging exoplanet yield around stars within 6.5 parsecs

abstract: Direct imaging is a powerful exoplanet discovery technique that is complementary to other techniques and offers great promise in the era of 30 meter class telescopes. Space-based transit surveys have revolutionized our understanding of the frequency of planets at small orbital radii around Sun-like stars. The next generation of extremely large ground-based telescopes will have the angular resolution and sensitivity to directly image planets with $R < 4R_\oplus$ around the very nearest stars. Here, we predict yields from a direct imaging survey of a volume-limited sample of Sun-like stars with the Mid-Infrared ELT Imager and Spectrograph (METIS) instrument, planned for the 39 m European Southern Observatory (ESO) Extremely Large Telescope (ELT) that is expected to be operational towards the end of the decade. Using Kepler occurrence rates, a sample of stars with spectral types A-K within 6.5 pc, and simulated contrast curves based on an advanced model of what is achievable from coronagraphic imaging with adaptive optics, we estimated the expected yield from METIS using Monte Carlo simulations. We find the METIS expected yield of planets in the N2 band (10.10 - 12.40 $\mu$m) is 1.14 planets, which is greater than comparable observations in the L (3.70 - 3.95 $\mu$m) and M (4.70 - 4.90 $\mu$m) bands. We also determined a 24.6% chance of detecting at least one Jovian planet in the background limited regime assuming a 1 hour integration. We calculated the yield per star and estimate optimal observing revisit times to increase the yield. We also analyzed a northern hemisphere version of this survey and found there are additional targets worth considering. In conclusion, we present an observing strategy aimed to maximize the possible yield for limited telescope time, resulting in 1.48 expected planets in the N2 band.


title: CARMENES input catalog of M dwarfs VI. A time-resolved Ca II H&K catalog from archival data

abstract: Radial-velocity (RV) jitter caused by stellar magnetic activity is an important factor in state-of-the-art exoplanet discovery surveys such as CARMENES. Stellar rotation, along with heterogeneities in the photosphere and chromosphere caused by activity, can result in false-positive planet detections. Hence, it is necessary to determine the stellar rotation period and compare it to any putative planetary RV signature. Long-term measurements of activity indicators such as the chromospheric emission in the Ca II H&K lines enable the identification of magnetic activity cycles. In order to determine stellar rotation periods and study the long-term behavior of magnetic activity of the CARMENES guaranteed time observations (GTO) sample, it is advantageous to extract Ca II H&K time series from archival data, since the CARMENES spectrograph does not cover the blue range of the stellar spectrum containing the Ca II H&K lines. We have assembled a catalog of 11634 archival spectra of 186 M dwarfs acquired by seven different instruments covering the Ca II H&K regime: ESPADONS, FEROS, HARPS, HIRES, NARVAL, TIGRE, and UVES. The relative chromospheric flux in these lines was directly extracted from the spectra by rectification with PHOENIX synthetic spectra via narrow passbands around the Ca ii H&K line cores. The combination of archival spectra from various instruments results in time series for 186 stars from the CARMENES GTO sample. As an example of the use of the catalog, we report the tentative discovery of three previously unknown activity cycles of M dwarfs. We conclude that the method of extracting Ca II H&K fluxes with the use of model spectra yields consistent results for different instruments and that the compilation of this catalog will enable the analysis of long-term activity time series for a large number of M dwarfs.


title: The Blazar Hadronic Code Comparison Project

abstract: Blazar hadronic models have been developed in the past decades as an alternative to leptonic ones. In hadronic models the gamma-ray emission is associated with synchrotron emission by protons, and/or secondary leptons produced in proton-photon interactions. Together with photons, hadronic emission models predict the emission of neutrinos that are therefore the smoking gun for acceleration of relativistic hadrons in blazar jets. The simulation of proton-photon interactions and all associated radiative processes is a complex numerical task, and different approaches to the problem have been adopted in the literature. So far, no systematic comparison between the different codes has been performed, preventing a clear understanding of the underlying uncertainties in the numerical simulations. To fill this gap, we have undertaken the first comprehensive comparison of blazar hadronic codes, and the results from this effort will be presented in this contribution.


title: Study of using Quantum Computer to Solve Poisson Equation in Gate Insulators

abstract: In this paper, the application of quantum computing (QC) in solving gate insulator Poisson equation is studied, through QC simulator and hardware in IBM. Various gate insulator stacks with and without fixed charges are studied. It is found that by increasing the number of clock bits and by choosing appropriate evolution time, accurate solutions can be obtained in QC simulation. However, when the real quantum hardware is used, the accuracy is substantially reduced. Therefore, a more robust quantum circuit or error correction should be employed and developed.


title: Separation of Learning and Control for Cyber-Physical Systems

abstract: Most cyber-physical systems (CPS) encounter a large volume of data which is added to the system gradually in real time and not altogether in advance. In this paper, we provide a theoretical framework that yields optimal control strategies for such CPS at the intersection of control theory and learning. In the proposed framework, we use the actual CPS, i.e., the "true" system that we seek to optimally control online, in parallel with a model of the CPS that is available. We then institute an information state for the system which does not depend on the control strategy. An important consequence of this independence is that for any given choice of a control strategy and a realization of the system's variables until time t, the information states at future times do not depend on the choice of the control strategy at time t but only on the realization of the decision at time t, and thus they are related to the concept of separation between estimation of the state and control. Namely, the future information states are separated from the choice of the current control strategy. Such control strategies are called separated control strategies. Hence, we can derive offline the optimal control strategy of the system with respect to the information state, which might not be precisely known due to model uncertainties or complexity of the system, and then use standard learning approaches to learn the information state online while data are added gradually to the system in real time. We show that after the information state becomes known, the separated control strategy of the CPS model derived offline is optimal for the actual system. We illustrate the proposed framework in a dynamic system consisting of two subsystems with a delayed sharing information structure.


title: Bivariate Lagrange interpolation at the checkerboard nodes

abstract: In this paper, we derive an explicit formula for the bivariate Lagrange basis polynomials of a general set of checkerboard nodes. This formula generalizes existing results of bivariate Lagrange basis polynomials at the Padua nodes, Chebyshev nodes, Morrow-Patterson nodes, and Geronimus nodes. We also construct a subspace spanned by linearly independent bivariate vanishing polynomials that vanish at the checkerboard nodes and prove the uniqueness of the set of bivariate Lagrange basis polynomials in the quotient space defined as the space of bivariate polynomials with a certain degree by the subspace of bivariate vanishing polynomials.


title: Fast Parallel-in-Time Quasi-Boundary Value Methods for Backward Heat Conduction Problems

abstract: In this paper we proposed two new quasi-boundary value methods for regularizing the ill-posed backward heat conduction problems. With a standard finite difference discretization in space and time, the obtained all-at-once nonsymmetric sparse linear systems have the desired block $\omega$-circulant structure, which can be utilized to design an efficient parallel-in-time (PinT) direct solver that built upon an explicit FFT-based diagonalization of the time discretization matrix. Convergence analysis is presented to justify the optimal choice of the regularization parameter. Numerical examples are reported to validate our analysis and illustrate the superior computational efficiency of our proposed PinT methods.


title: Stochastic Geometry based Interference Analysis of Multiuser mmWave Networks with RIS

abstract: In this paper, we utilize tools from stochastic geometry to estimate the interference propagation via reconfigurable intelligent surface (RIS) in the millimeter wave (mmWave, 30-300 GHz) band and specifically on the D band (110-170 GHz). The RISs have been of great interest lately to maximize the channel gains in non-line-of-sight (NLOS) communication situations. We derive expressions for stochastic interference level in RIS powered systems and validate those with simulations. It will be shown that the interference levels via RIS link are rather small compared to the designed RIS link or the LOS interference as the random interference loses significant part of the RIS gain. We also analyse the validity of far field channel and antenna gains in the near field of a large array. It is shown that, while the high frequency systems require large arrays that push the far field far away from the antenna, the far field equations are very accurate up to about half way of the near field.


title: How Much Can CLIP Benefit Vision-and-Language Tasks?

abstract: Most existing Vision-and-Language (V&L) models rely on pre-trained visual encoders, using a relatively small set of manually-annotated data (as compared to web-crawled data), to perceive the visual world. However, it has been observed that large-scale pretraining usually can result in better generalization performance, e.g., CLIP (Contrastive Language-Image Pre-training), trained on a massive amount of image-caption pairs, has shown a strong zero-shot capability on various vision tasks. To further study the advantage brought by CLIP, we propose to use CLIP as the visual encoder in various V&L models in two typical scenarios: 1) plugging CLIP into task-specific fine-tuning; 2) combining CLIP with V&L pre-training and transferring to downstream tasks. We show that CLIP significantly outperforms widely-used visual encoders trained with in-domain annotated data, such as BottomUp-TopDown. We achieve competitive or better results on diverse V&L tasks, while establishing new state-of-the-art results on Visual Question Answering, Visual Entailment, and V&L Navigation tasks. We release our code at https://github.com/clip-vil/CLIP-ViL.


title: Bellman function method for general operators on martingales

abstract: It is shown that the Bellman function method can be applied to study the $L^p$-norms of general operators on martingales, i.e., of operators that are not necessarily martingale transforms. Informally, we provide a single Bellman-type function that "encodes" the $L^p$-boundedness of "almost all" operators from Gundy's extrapolation theorem. As examples of such operators, we consider the Haar transforms and the operator whose $L^p$-boundedness underlies Rubio de Francia's inequality for the Walsh system.


title: Automorphisms of simple quotients of the Poisson and universal enveloping algebras of $\mathrm{sl}_2$

abstract: Let $P(\mathrm{sl}_2(K))$ be the Poisson enveloping algebra of the Lie algebra $\mathrm{sl}_2(K)$ over an algebraically closed field $K$ of characteristic zero. The quotient algebras $ $ $P(\mathrm{sl}_2(K))/(C_P-\lambda)$, where $C_P$ is the standard Casimir element of $\mathrm{sl}_2(K)$ in $P(\mathrm{sl}_2(K))$ and $0\neq \lambda\in K$, are proven to be simple in \cite{UZh}. Using a result by L. Makar-Limanov \cite{ML90}, we describe generators of the automorphism group of $P(\mathrm{sl}_2(K))/(C_P-\lambda)$ and represent this group as an amalgamated product of its subgroups. Moreover, using similar results by J. Dixmier \cite{Dixmier73} and O. Fleury \cite{Fleury} for the quotient algebras $U(\mathrm{sl}_2(K))/(C_U-\lambda)$, where $C_U$ is the standard Casimir element of $\mathrm{sl}_2(K)$ in the universal enveloping algebra $U(\mathrm{sl}_2(K))$, we prove that the automorphism groups of $P(\mathrm{sl}_2(K))/(C_P-\lambda)$ and $U(\mathrm{sl}_2(K))/(C_U-\lambda)$ are isomorphic.


title: Geometry and Generalization: Eigenvalues as predictors of where a network will fail to generalize

abstract: We study the deformation of the input space by a trained autoencoder via the Jacobians of the trained weight matrices. In doing so, we prove bounds for the mean squared errors for points in the input space, under assumptions regarding the orthogonality of the eigenvectors. We also show that the trace and the product of the eigenvalues of the Jacobian matrices is a good predictor of the MSE on test points. This is a dataset independent means of testing an autoencoder's ability to generalize on new input. Namely, no knowledge of the dataset on which the network was trained is needed, only the parameters of the trained model.


title: Finite memory time and anisotropy effects for initial magnetic energy growth in random flow of conducting media

abstract: The dynamo mechanism is a process of magnetic field self-excitation in a moving electrically conducting fluid. One of the most interesting applications of this mechanism related to the astrophysical systems is the case of a random motion of plasma. For the very first stage of the process, the governing dynamo equation can be reduced to a system of first order ordinary differential equations. For this case we suggest a regular method to calculate the growth rate of magnetic energy. Based on this method we calculate the growth rate for random flow with finite memory time and anisotropic statistical distribution of the stretching matrix and compare the results with corresponding ones for isotropic case and for short-correlated approximation. We find that for moderate Strouhal numbers and moderate anisotropy the analytical results reproduce the numerically estimated growth rates reasonably well, while for larger governing parameters the quantitative difference becomes substantial. In particular, analytical approximation is applicable for the Strouhal numbers s<0.6 and we find some numerical models and observational examples for which this region might be relevant. Rather unexpectedly, we find that the mirror asymmetry does not contribute to the growth rates obtained, although the mirror asymmetry effects are known to be crucial for later stages of dynamo action.


title: Whiteout: when do fixed-X knockoffs fail?

abstract: A core strength of knockoff methods is their virtually limitless customizability, allowing an analyst to exploit machine learning algorithms and domain knowledge without threatening the method's robust finite-sample false discovery rate control guarantee. While several previous works have investigated regimes where specific implementations of knockoffs are provably powerful, general negative results are more difficult to obtain for such a flexible method. In this work we recast the fixed-$X$ knockoff filter for the Gaussian linear model as a conditional post-selection inference method. It adds user-generated Gaussian noise to the ordinary least squares estimator $\hat\beta$ to obtain a "whitened" estimator $\widetilde\beta$ with uncorrelated entries, and performs inference using $\text{sgn}(\widetilde\beta_j)$ as the test statistic for $H_j:\; \beta_j = 0$. We prove equivalence between our whitening formulation and the more standard formulation involving negative control predictor variables, showing how the fixed-$X$ knockoffs framework can be used for multiple testing on any problem with (asymptotically) multivariate Gaussian parameter estimates. Relying on this perspective, we obtain the first negative results that universally upper-bound the power of all fixed-$X$ knockoff methods, without regard to choices made by the analyst. Our results show roughly that, if the leading eigenvalues of $\text{Var}(\hat\beta)$ are large with dense leading eigenvectors, then there is no way to whiten $\hat\beta$ without irreparably erasing nearly all of the signal, rendering $\text{sgn}(\widetilde\beta_j)$ too uninformative for accurate inference. We give conditions under which the true positive rate (TPR) for any fixed-$X$ knockoff method must converge to zero even while the TPR of Bonferroni-corrected multiple testing tends to one, and we explore several examples illustrating this phenomenon.


title: Prime power variations of higher $Lie_n$ modules

abstract: We define, for each subset $S$ of the set $\mathcal{P}$ of primes, an $S_n$-module $Lie_n^S$ with interesting properties. $Lie_n^\emptyset$ is the well-known representation $Lie_n$ of $S_n$ afforded by the free Lie algebra, while $Lie_n^\mathcal{P}$ is the module $C\!onj_n$ of the conjugacy action of $S_n$ on $n$-cycles. For arbitrary $S$ the module $Lie_n^{S}$ interpolates between the representations $Lie_n$ and $C\!onj_n.$ We consider the symmetric and exterior powers of $Lie_n^S.$ These are the analogues of the higher Lie modules of Thrall. We show that the Frobenius characteristic of these higher $Lie_n^S$ modules can be elegantly expressed as a multiplicity-free sum of power sums. In particular this establishes the Schur positivity of new classes of sums of power sums. More generally, for each nonempty subset $T$ of positive integers we define a sequence of symmetric functions $f_n^T$ of homogeneous degree $n.$ We show that the series $\sum_{\lambda, \lambda_i\in T} p_\lambda$ can be expressed as symmetrised powers of the functions $f_n^T$, analogous to the higher Lie modules first defined by Thrall. This in turn allows us to unify previous results on the Schur positivity of multiplicity-free sums of power sums, as well as investigate new ones. We also uncover some curious plethystic relationships between $f_n^T$, the conjugacy action and the Lie representation.


title: Spectral diffusion of phosphorus donors in silicon at high magnetic field

abstract: We characterize the phase memory time of phosphorus donor electron spins in lightly-doped natural silicon at high magnetic field (8.58 T) in the dark and under low-power optical excitation. The spin echo decays are dominated by spectral diffusion due to the presence of the 4.7% abundant spin-1/2 silicon-29 nuclei. At 4.2 K, the spectral diffusion time (T$_{SD}$) measured in the dark is $124 \pm 7$ $\mu$s, a factor of 2 smaller than that measured at low magnetic fields (0.35 T). Using a tunable laser we also measured the echo decay as the wavelength of the optical excitation is swept across the band edge from 1050 nm to 1090 nm. Above-bandgap optical excitation is seen to increase the spectral diffusion time of the donor electron spin to $201 \pm 11$ $\mu$s. The physical mechanism underlying both the decrease of T$_{SD}$ at high field and the subsequent increase under optical excitation remains unclear.


title: Observation of quasi-ballistic thermal transport of surface phonon-polaritons over hundreds of micrometres

abstract: Long-distance propagation of heat carriers is essential for efficient heat dissipation in microelectronics. However, in dielectric nanomaterials, the primary heat carriers - phonons - can propagate ballistically only for hundreds of nanometres, which limits their heat conduction efficiency. Theory predicts that surface phonon-polaritons (SPhPs) can overcome this limitation and conduct heat without dissipation for hundreds of micrometres. In this work, we experimentally demonstrate such long-distance heat transport by SPhPs. Using the 3$\omega$ technique, we measure the in-plane thermal conductivity of SiN nanomembranes for different heater-sensor distances (100 and 200 $\mu$m), membrane thicknesses (30 - 200 nm), and temperatures (300 - 400 K). We find that in contrast with thick membranes, thin nanomembranes support heat conduction by SPhPs, as evidenced by an increase in the thermal conductivity with temperature. Remarkably, the thermal conductivity measured 200 $\mu$m away from the heater are consistently higher than that measured 100 $\mu$m closer. This result suggests that heat conduction by SPhPs is quasi-ballistic over at least hundreds of micrometres. Thus, our findings show that SPhPs can enhance heat dissipation in polar nanomembranes and find applications in thermal management, near-field radiation, and polaritonics.


title: Active Crowds

abstract: This chapter focuses on the mathematical modelling of active particles (or agents) in crowded environments. We discuss several microscopic models found in literature and the derivation of the respective macroscopic partial differential equations for the particle density. The macroscopic models share common features, such as cross diffusion or degenerate mobilities. We then take the diversity of macroscopic models to a uniform structure and work out potential similarities and differences. Moreover, we discuss boundary effects and possible applications in life and social sciences. This is complemented by numerical simulations that highlight the effects of different boundary conditions.


title: Hybrid Memoised Wake-Sleep: Approximate Inference at the Discrete-Continuous Interface

abstract: Modeling complex phenomena typically involves the use of both discrete and continuous variables. Such a setting applies across a wide range of problems, from identifying trends in time-series data to performing effective compositional scene understanding in images. Here, we propose Hybrid Memoised Wake-Sleep (HMWS), an algorithm for effective inference in such hybrid discrete-continuous models. Prior approaches to learning suffer as they need to perform repeated expensive inner-loop discrete inference. We build on a recent approach, Memoised Wake-Sleep (MWS), which alleviates part of the problem by memoising discrete variables, and extend it to allow for a principled and effective way to handle continuous variables by learning a separate recognition model used for importance-sampling based approximate inference and marginalization. We evaluate HMWS in the GP-kernel learning and 3D scene understanding domains, and show that it outperforms current state-of-the-art inference methods.


title: Compressive Representations of Weather Scenes for Strategic Air Traffic Flow Management

abstract: Terse representation of high-dimensional weather scene data is explored, in support of strategic air traffic flow management objectives. Specifically, we consider whether aviation-relevant weather scenes are compressible, in the sense that each scene admits a possibly-different sparse representation in a basis of interest. Here, compression of weather scenes extracted from METAR data (including temperature, flight categories, and visibility profiles for the contiguous United States) is examined, for the graph-spectral basis. The scenes are found to be compressible, with 75-95% of the scene content captured using 0.5-4% of the basis vectors. Further, the dominant basis vectors for each scene are seen to identify time-varying spatial characteristics of the weather, and reconstruction from the compressed representation is demonstrated. Finally, potential uses of the compressive representations in strategic TFM design are briefly scoped.


title: Pulling Radiation Force State on Active Spherical Carriers via Varying Phase Incident Wave Field

abstract: On the way of investigating a complete maneuverability technique and based on the emergence of negative acoustic radiation force on an active carrier in a progressive plane wave field, we propose a new technique in which not only the attractive force is experienced but also the average of acoustic radiation force over the location of the object with respect to the sound source is negative. Keeping in mind that the produced pulling force upon a carrier activated on only one mode always have a positive average, along with the point that the positive average force means that the net motion will be pushing in lieu of a pulling effect associated with a negative force, made us to put forward a technique to generate the acoustic radiation force with negative average over a complete wavelength to ensure that the net motion is attractive. The idea here is to excite the object in two modes (zero and first) synchronously with a varying difference between the phases of excitations on two modes with the hope of generating acoustic radiation force with a negative average over a complete wavelength.


title: Forecasting Thermoacoustic Instabilities in Liquid Propellant Rocket Engines Using Multimodal Bayesian Deep Learning

abstract: The 100 MW cryogenic liquid oxygen/hydrogen multi-injector combustor BKD operated by the DLR Institute of Space Propulsion is a research platform that allows the study of thermoacoustic instabilities under realistic conditions, representative of small upper stage rocket engines. We use data from BKD experimental campaigns in which the static chamber pressure and fuel-oxidizer ratio are varied such that the first tangential mode of the combustor is excited under some conditions. We train an autoregressive Bayesian neural network model to forecast the amplitude of the dynamic pressure time series, inputting multiple sensor measurements (injector pressure/ temperature measurements, static chamber pressure, high-frequency dynamic pressure measurements, high-frequency OH* chemiluminescence measurements) and future flow rate control signals. The Bayesian nature of our algorithms allows us to work with a dataset whose size is restricted by the expense of each experimental run, without making overconfident extrapolations. We find that the networks are able to accurately forecast the evolution of the pressure amplitude and anticipate instability events on unseen experimental runs 500 milliseconds in advance. We compare the predictive accuracy of multiple models using different combinations of sensor inputs. We find that the high-frequency dynamic pressure signal is particularly informative. We also use the technique of integrated gradients to interpret the influence of different sensor inputs on the model prediction. The negative log-likelihood of data points in the test dataset indicates that predictive uncertainties are well-characterized by our Bayesian model and simulating a sensor failure event results as expected in a dramatic increase in the epistemic component of the uncertainty.


title: SurgeonAssist-Net: Towards Context-Aware Head-Mounted Display-Based Augmented Reality for Surgical Guidance

abstract: We present SurgeonAssist-Net: a lightweight framework making action-and-workflow-driven virtual assistance, for a set of predefined surgical tasks, accessible to commercially available optical see-through head-mounted displays (OST-HMDs). On a widely used benchmark dataset for laparoscopic surgical workflow, our implementation competes with state-of-the-art approaches in prediction accuracy for automated task recognition, and yet requires 7.4x fewer parameters, 10.2x fewer floating point operations per second (FLOPS), is 7.0x faster for inference on a CPU, and is capable of near real-time performance on the Microsoft HoloLens 2 OST-HMD. To achieve this, we make use of an efficient convolutional neural network (CNN) backbone to extract discriminative features from image data, and a low-parameter recurrent neural network (RNN) architecture to learn long-term temporal dependencies. To demonstrate the feasibility of our approach for inference on the HoloLens 2 we created a sample dataset that included video of several surgical tasks recorded from a user-centric point-of-view. After training, we deployed our model and cataloged its performance in an online simulated surgical scenario for the prediction of the current surgical task. The utility of our approach is explored in the discussion of several relevant clinical use-cases. Our code is publicly available at https://github.com/doughtmw/surgeon-assist-net.


title: Planning a method for covariate adjustment in individually-randomised trials: a practical guide

abstract: Background: It has long been advised to account for baseline covariates in the analysis of confirmatory randomised trials, with the main statistical justifications being that this increases power and, when a randomisation scheme balanced covariates, permits a valid estimate of experimental error. There are various methods available to account for covariates but it is not clear how to choose among them. Methods: Taking the perspective of writing a statistical analysis plan, we consider how to choose between the three most promising broad approaches: direct adjustment, standardisation and inverse-probability-of-treatment weighting. Results: The three approaches are similar in being asymptotically efficient, in losing efficiency with mis-specified covariate functions, and in handling designed balance. If a marginal estimand is targeted (for example, a risk difference or survival difference), then direct adjustment should be avoided because it involves fitting non-standard models that are subject to convergence issues. Convergence is most likely with IPTW. Robust standard errors used by IPTW are anti-conservative at small sample sizes. All approaches can use similar methods to handle missing covariate data. With missing outcome data, each method has its own way to estimate a treatment effect in the all-randomised population. We illustrate some issues in a reanalysis of GetTested, a randomised trial designed to assess the effectiveness of an electonic sexually-transmitted-infection testing and results service. Conclusions: No single approach is always best: the choice will depend on the trial context. We encourage trialists to consider all three methods more routinely.


title: The Perfect Matching Cut Problem Revisited

abstract: In a graph, a perfect matching cut is an edge cut that is a perfect matching. Perfect Matching Cut (PMC) is the problem of deciding whether a given graph has a perfect matching cut, and is known to be NP-complete. We revisit the problem and show that PMC remains NP-complete when restricted to bipartite graphs of maximum degree 3 and arbitrarily large girth. Complementing this hardness result, we give two graph classes in which PMC is polynomial time solvable. The first one includes claw-free graphs and graphs without an induced path on five vertices, the second one properly contains all chordal graphs. Assuming the Exponential Time Hypothesis, we show there is no $O^*(2^{o(n)})$-time algorithm for PMC even when restricted to $n$-vertex bipartite graphs, and also show that PMC can be solved in $O^*(1.2721^n)$ time by means of an exact branching algorithm.