title: Gaussian Process Latent Class Choice Models

abstract: We present a Gaussian Process - Latent Class Choice Model (GP-LCCM) to integrate a non-parametric class of probabilistic machine learning within discrete choice models (DCMs). Gaussian Processes (GPs) are kernel-based algorithms that incorporate expert knowledge by assuming priors over latent functions rather than priors over parameters, which makes them more flexible in addressing nonlinear problems. By integrating a Gaussian Process within a LCCM structure, we aim at improving discrete representations of unobserved heterogeneity. The proposed model would assign individuals probabilistically to behaviorally homogeneous clusters (latent classes) using GPs and simultaneously estimate class-specific choice models by relying on random utility models. Furthermore, we derive and implement an Expectation-Maximization (EM) algorithm to jointly estimate/infer the hyperparameters of the GP kernel function and the class-specific choice parameters by relying on a Laplace approximation and gradient-based numerical optimization methods, respectively. The model is tested on two different mode choice applications and compared against different LCCM benchmarks. Results show that GP-LCCM allows for a more complex and flexible representation of heterogeneity and improves both in-sample fit and out-of-sample predictive power. Moreover, behavioral and economic interpretability is maintained at the class-specific choice model level while local interpretation of the latent classes can still be achieved, although the non-parametric characteristic of GPs lessens the transparency of the model.


title: Adiabatic evolution of Hayward black hole

abstract: In this letter we use the Carath\'{e}odory's approach to thermodynamics, to construct the thermodynamic manifold of the Hayward black hole. The Pfaffian form representing the infinitesimal heat exchange reversibly is considered to be $\delta Q_{\mathrm{rev}}\equiv {\rm d}r_s-\mathcal{F}_H {\rm d}l$, previously obtained by Molina \& Villanueva \cite{fmv20}, where $r_s$ is the Schwarzschild radius, $l$ is the Hayward's parameter responsible for the possible regularization of the Schwarzschild black hole, and $\mathcal{F}_H$ is the intensive variable called the Hayward's force. By solving the associated Cauchy problem, the adiabatic paths are confined to the non-extremal manifold, and therefore, the status of the second and third laws are preserved. Consequently, the extremal sub-manifold corresponds to the {adiabatically disconnected} boundary of the manifold. In addition, the merger of two extremal Hayward black holes is analyzed.


title: Reliable COVID-19 Detection Using Chest X-ray Images

abstract: Coronavirus disease 2019 (COVID-19) has emerged the need for computer-aided diagnosis with automatic, accurate, and fast algorithms. Recent studies have applied Machine Learning algorithms for COVID-19 diagnosis over chest X-ray (CXR) images. However, the data scarcity in these studies prevents a reliable evaluation with the potential of overfitting and limits the performance of deep networks. Moreover, these networks can discriminate COVID-19 pneumonia usually from healthy subjects only or occasionally, from limited pneumonia types. Thus, there is a need for a robust and accurate COVID-19 detector evaluated over a large CXR dataset. To address this need, in this study, we propose a reliable COVID-19 detection network: ReCovNet, which can discriminate COVID-19 pneumonia from 14 different thoracic diseases and healthy subjects. To accomplish this, we have compiled the largest COVID-19 CXR dataset: QaTa-COV19 with 124,616 images including 4603 COVID-19 samples. The proposed ReCovNet achieved a detection performance with 98.57% sensitivity and 99.77% specificity.


title: Hybrid leg compliance enables robots to operate with sensorimotor delays and low control update frequencies

abstract: Animals locomote robustly and agile, albeit significant sensorimotor delays of their nervous system. The sensorimotor control of legged robots is implemented with much higher frequencies-often in the kilohertz range-and sensor and actuator delays in the low millisecond range. But especially at harsh impacts with unknown touch-down timing, legged robots show unstable controller behaviors, while animals are seemingly not impacted. Here we examine this discrepancy and suggest a hybrid robotic leg and controller design. We implemented a physical, parallel joint compliance dimensioned in combination with an active, virtual leg length controller. We present an extensive set of systematic experiments both in computer simulation and hardware. Our hybrid leg and controller design shows previously unseen robustness, in the presence of sensorimotor delays up to 60 ms, or control frequencies as low as 20 Hz, for a drop landing task from 1.3 leg lengths high and with a passive compliance ratio of 0.7. In computer simulations, we report successful drop-landings of the hybrid compliant leg from 3.8 leg lengths (1.2 m) for a 2 kg quadruped robot with 100 Hz control frequency and a sensorimotor delay of 35 ms. The results of our presented hybrid leg design and control provide a further explanation for the performance robustness of animals, and the resulting discrepancy between animals and legged robots.


title: Generation possibility of gamma-ray glows induced by photonuclear reactions

abstract: Relativistic runaway electron avalanches (RREAs) imply a large multiplication of high energy electrons (~1 MeV). Two factors are necessary for this phenomenon: a high electric field sustained over a large distance and an energetic particle to serve as a seed. The former sustains particle energies as they keep colliding and lose energy randomly, and the latter serves as a multiplication starting point that promotes avalanches. RREA is usually connected to both terrestrial gamma-ray flashes (TGFs) and gamma-ray glows (also known as Thunderstorm Ground Enhancement (TGE) when detected at ground level) as possible generation mechanism of both events, but the current knowledge does not provide a clear relationship between these events (TGF and TGE), beyond their possible common source mechanism, still as they have different characteristics. In particular, their timescales differ by several orders of magnitude. This work shows that chain reactions by TGF byproducts can continue for the timescale of gamma-ray glows and even provide energetic particles as seeds for RREAs of gamma-ray glows.


title: Integrals of Motion in Time-periodic Hamiltonian Systems: The Case of the Mathieu Equation

abstract: We present an algorithm for constructing analytically approximate integrals of motion in simple time periodic Hamiltonians of the form $H=H_0+ \varepsilon H_i$, where $\varepsilon$ is a perturbation parameter. We apply our algorithm in a Hamiltonian system whose dynamics is governed by the Mathieu equation and examine in detail the orbits and their stroboscopic invariant curves for different values of $\varepsilon$. We find the values of $\varepsilon_{crit}$ beyond which the orbits escape to infinity and construct integrals which are expressed as series in the perturbation $\varepsilon$ and converge up to $\varepsilon_{crit}$. In the absence of resonances the invariant curves are concentric ellipses which are approximated very well by our integrals. Finally we construct an integral of motion which describes the hyperbolic stroboscopic invariant curve of a resonant case.


title: Microswimmers learning chemotaxis with genetic algorithms

abstract: Various microorganisms and some mammalian cells are able to swim in viscous fluids by performing nonreciprocal body deformations, such as rotating attached flagella or by distorting their entire body. In order to perform chemotaxis, i.e. to move towards and to stay at high concentrations of nutrients, they adapt their swimming gaits in a nontrivial manner. We propose a model how microswimmers are able to autonomously adapt their shape in order to swim in one dimension towards high field concentrations using an internal decision making machinery modeled by an artificial neural network. We present two methods to measure chemical gradients, spatial and temporal sensing, as known for swimming mammalian cells and bacteria, respectively. Using the NEAT genetic algorithm surprisingly simple neural networks evolve which control the shape deformations of the microswimmer and allow them to navigate in static and complex time-dependent chemical environments. By introducing noisy signal transmission in the neural network the well-known biased run-and-tumble motion emerges. Our work demonstrates that the evolution of a simple internal decision-making machinery, which we can fully interpret and is coupled to the environment, allows navigation in diverse chemical landscapes. These findings are of relevance for intracellular biochemical sensing mechanisms of single cells, or for the simple nervous system of small multicellular organisms such as C. elegans.

categories: physics.bio-ph cond-mat.dis-nn cond-mat.soft physics.flu-dyn


title: Exceptional surgeries in 3-manifolds

abstract: Myers shows that every compact, connected, orientable $3$--manifold with no $2$--sphere boundary components contains a hyperbolic knot. We use work of Ikeda with an observation of Adams-Reid to show that every $3$--manifold subject to the above conditions contains a hyperbolic knot which admits a non-trivial non-hyperbolic surgery, a toroidal surgery in particular. We conclude with a question and a conjecture about reducible surgeries.


title: Feedback in Forming Star Clusters: The Mass-Radius Relation and Mass Function of Molecular Clumps in the Large Magellanic Cloud

abstract: We derive the mass-radius relation and mass function of molecular clumps in the Large Magellanic Cloud (LMC) and interpret them in terms of the simple feedback model proposed by Fall, Krumholz, and Matzner (FKM). Our work utilizes the dendrogram-based catalog of clumps compiled by Wong et al. from $^{12}$CO and $^{13}$CO maps of six giant molecular clouds in the LMC observed with the Atacama Large Millimeter Array (ALMA). The Magellanic Clouds are the only external galaxies for which this type of analysis is possible at the necessary spatial resolution ($\sim1$ pc). We find that the mass-radius relation and mass function of LMC clumps have power-law forms, $R \propto M^{\alpha}$ and $dN/dM \propto M^{\beta}$, with indices $\alpha = 0.36 \pm 0.03$ and $\beta= -1.8 \pm 0.1 $ over the mass ranges $10^2 M_\odot \lesssim M \lesssim 10^5 M_\odot$ and $10^2 M_\odot \lesssim M \lesssim 10^4 M_\odot$, respectively. With these values of $\alpha$ and $\beta$ for the clumps (i.e., protoclusters), the predicted index for the mass function of young LMC clusters from the FKM model is $\beta \approx 1.7$, in good agreement with the observed index. The situation portrayed here for clumps and clusters in the LMC replicates that in the Milky Way.


title: Five-percent measurements of the growth rate from simulation-based modelling of redshift-space clustering in BOSS LOWZ

abstract: We use a simulation-based modelling approach to analyse the anisotropic clustering of the BOSS LOWZ sample over the radial range $0.4 \, h^{-1} \, \mathrm{Mpc}$ to $63 \, h^{-1} \, \mathrm{Mpc}$, significantly extending what is possible with a purely analytic modelling framework. Our full-scale analysis yields constraints on the growth of structure that are a factor of two more stringent than any other study on large scales at similar redshifts. We infer $f \sigma_8 = 0.471 \pm 0.024$ at $z \approx 0.25$, and $f \sigma_8 = 0.431 \pm 0.025$ at $z \approx 0.40$; the corresponding $\Lambda$CDM predictions of the Planck CMB analysis are $0.470 \pm 0.006$ and $0.476 \pm 0.005$, respectively. Our results are thus consistent with Planck, but also follow the trend seen in previous low-redshift measurements of $f \sigma_8$ falling slightly below the $\Lambda$CDM+CMB prediction. We find that small and large radial scales yield mutually consistent values of $f \sigma_8$, but there are $1-2.5 \sigma$ hints of small scales ($< 10 \, h^{-1} \, \mathrm{Mpc}$) preferring lower values for $f \sigma_8$ relative to larger scales. We analyse the constraining power of the full range of radial scales, finding that most of the multipole information about $f\sigma_8$ is contained in the scales $2 \, h^{-1} \, \mathrm{Mpc} \lesssim s \lesssim 20 \, h^{-1} \, \mathrm{Mpc}$. Evidently, once the cosmological information of the quasi-to-nonlinear regime has been harvested, large-scale modes contain only modest additional information about structure growth. Finally, we compare predictions for the galaxy-galaxy lensing amplitude of the two samples against measurements from SDSS and assess the lensing-is-low effect in light of our findings.


title: Measuring non-exchangeable tail dependence using tail copulas

abstract: Quantifying tail dependence is an important issue in insurance and risk management. The prevalent tail dependence coefficient (TDC), however, is known to underestimate the degree of tail dependence and it does not capture non-exchangeable tail dependence since it evaluates the limiting tail probability only along the main diagonal. To overcome these issues, two novel tail dependence measures called the maximal tail concordance measure (MTCM) and the average tail concordance measure (ATCM) are proposed. Both measures are constructed based on tail copulas and possess clear probabilistic interpretations in that the MTCM evaluates the largest limiting probability among all comparable rectangles in the tail, and the ATCM is a normalized average of these limiting probabilities. In contrast to the TDC, the proposed measures can capture non-exchangeable tail dependence. Analytical forms of the proposed measures are also derived for various copulas. A real data analysis reveals striking tail dependence and tail non-exchangeability of the return series of stock indices, particularly in periods of financial distress.


title: Explicit zero density for the Riemann zeta function

abstract: Let $N(\sigma,T)$ denote the number of nontrivial zeros of the Riemann zeta function with real part greater than $\sigma$ and imaginary part between $0$ and $T$. We provide explicit upper bounds for $N(\sigma,T)$ commonly referred to as a zero density result. In 1937, Ingham showed the following asymptotic result $N(\sigma,T)=\mathcal{O} ( T^{\frac83(1-\sigma)} (\log T)^5 )$. Ramar\'{e} recently proved an explicit version of this estimate. We discuss a generalization of the method used in these two results which yields an explicit bound of a similar shape while also improving the constants.


title: On the Kodaira dimension of Hurwitz spaces

abstract: We show that the canonical bundle of the Hurwitz stack classifying covers of genus g>1 and degree k>2 of the projective line is big. We show that all coarse moduli spaces of trigonal curves of genus g>1 are of general type.


title: Speckle Statistics of Biological Tissues in Optical Coherence Tomography

abstract: The speckle statistics of optical coherence tomography images of biological tissue have been studied using several historical probability density functions. A recent hypothesis implies that underlying power-law distributions in the medium structure, such as the fractal branching vasculature, will contribute to power-law probability distributions of speckle statistics. Specifically, these are the Burr type XII distribution for speckle amplitude, the Lomax distribution for intensity, and the generalized logistic distribution for log amplitude. In this study, these three distributions are fitted to histogram data from nine optical coherence tomography scans of various biological tissues and samples. The distributions are also compared with conventional distributions such as the Rayleigh, K, and gamma distributions. The results indicate that these newer distributions based on power laws are, in general, more appropriate models and support the plausibility of their use for characterizing biological tissue. Potentially, the governing power-law parameter of these distributions could be used as a biomarker for tissue disease or pathology.


title: Detecting violations of macrorealism when the original Leggett-Garg inequalities are satisfied

abstract: Macroscopic realism (MR) is the notion that a time-evolving system possesses definite properties, irrespective of past or future measurements. Quantum mechanical theories can, however, produce violations of MR. Most research to date has focused on a single set of conditions for MR, the Leggett-Garg inequalities (LGIs), and on a single data set, the "standard data set", which consists of single-time averages and second-order correlators of a dichotomic variable Q for three times. However, if such conditions are all satisfied, then where is the quantum behaviour? In this paper, we provide an answer to this question by considering expanded data sets obtained from finer-grained measurements and MR conditions on those sets. We consider three different situations in which there are violations of MR that go undetected by the standard LGIs. First, we explore higher-order LGIs on a data set involving third- and fourth-order correlators, using a spin-1/2 and spin-1 system. Second, we explore the pentagon inequalities (PIs) and a data set consisting of all possible averages and second-order correlators for measurements of Q at five times. Third, we explore the LGIs for a trichotomic variable and measurements made with a trichotomic operator to, again, identify violations for a spin-1 system beyond those seen with a single dichotomic variable. We also explore the regimes in which combinations of two and three-time LGIs can be satisfied and violated in a spin-1 system, extending recent work. We discuss the possible experimental implementation of all the above results.


title: A Bayesian Model of Cash Bail Decisions

abstract: The use of cash bail as a mechanism for detaining defendants pre-trial is an often-criticized system that many have argued violates the presumption of "innocent until proven guilty." Many studies have sought to understand both the long-term effects of cash bail's use and the disparate rate of cash bail assignments along demographic lines (race, gender, etc). However, such work is often susceptible to problems of infra-marginality -- that the data we observe can only describe average outcomes, and not the outcomes associated with the marginal decision. In this work, we address this problem by creating a hierarchical Bayesian model of cash bail assignments. Specifically, our approach models cash bail decisions as a probabilistic process whereby judges balance the relative costs of assigning cash bail with the cost of defendants potentially skipping court dates, and where these skip probabilities are estimated based upon features of the individual case. We then use Monte Carlo inference to sample the distribution over these costs for different magistrates and across different races. We fit this model to a data set we have collected of over 50,000 court cases in the Allegheny and Philadelphia counties in Pennsylvania. Our analysis of 50 separate judges shows that they are uniformly more likely to assign cash bail to black defendants than to white defendants, even given identical likelihood of skipping a court appearance. This analysis raises further questions about the equity of the practice of cash bail, irrespective of its underlying legal justification.


title: Trace estimates of Toeplitz operators on Bergman spaces and applications to composition operators

abstract: Let $\Omega$ be a subdomain of $\mathbb{C}$ and let $\mu$ be a positive Borel measure on $\Omega$. In this paper, we study the asymptotic behavior of the eigenvalues of compact Toeplitz operator $T_\mu$ acting on Bergman spaces on $\Omega$. Let $(\lambda_n(T_\mu))$ be the decreasing sequence of the eigenvalues of $T_\mu$ and let $\rho$ be an increasing function such that $\rho (n)/n^A$ is decreasing for some $A>0$. We give an explicit necessary and sufficient geometric condition on $\mu$ in order to have $\lambda_n(T_\mu)\asymp 1/\rho (n)$. As applications, we consider composition operators $C_\varphi$, acting on some standard analytic spaces on the unit disc $\mathbb{D}$. First, we give a general criterion ensuring that the singular values of $C_\varphi$ satisfy $s_n(C_\varphi ) \asymp 1/\rho(n)$. Next, we focus our attention on composition operators with univalent symbols, where we express our general criterion in terms of the harmonic measure of $\varphi \mathbb{D})$. We finally study the case where $\partial \varphi (\mathbb{D})$ meets the unit circle in one point and give several concrete examples. Our method is based on upper and lower estimates of the trace of $h(T_\mu)$, where $h$ is suitable concave or convex functions.


title: Prospects for detecting gravitational waves from eccentric subsolar mass compact binaries

abstract: Due to their small mass, subsolar mass black hole binaries would have to be primordial in origin instead of the result of stellar evolution. Soon after formation in the early universe, primordial black holes can form binaries after decoupling from the cosmic expansion. Alternatively, primordial black holes as dark matter could also form binaries in the late universe due to dynamical encounters and gravitational-wave braking. A significant feature for this channel is the possibility that some sources retain nonzero eccentricity in the LIGO/Virgo band. Assuming all dark matter is primordial black holes with a delta function mass distribution, $1M_\odot-1M_\odot$ binaries formed in this late-universe channel can be detected by Advanced LIGO and Virgo with their design sensitivities at a rate of $\mathcal{O}(1)$/year, where $12\%(3\%)$ of events have eccentricity at a gravitational-wave frequency of 10 Hz, $e^\mathrm{10Hz}\geq0.01(0.1)$, and nondetection can constrain the binary formation rate within this model. Third generation detectors would be expected to detect subsolar mass eccentric binaries as light as $0.01 M_\odot$ within this channel, if they accounted for the majority of the dark matter. Furthermore, we use simulated gravitational-wave data to study the ability to search for eccentric gravitational-wave signals using a quasi-circular waveform template bank with Advanced LIGO design sensitivity. For a match-filtering targeted search, assuming binaries with a delta function mass of $0.1(1)M_\odot$ and the eccentricity distribution derived from this late-universe formation channel, $41\%(6\%)$ of the signals would be missed compared to the ideal detection rate due to the mismatch in the gravitational-wave signal from eccentricity.


title: An Experimental Analysis of Attack Classification Using Machine Learning in IoT Networks

abstract: In recent years, there has been a massive increase in the amount of Internet of Things (IoT) devices as well as the data generated by such devices. The participating devices in IoT networks can be problematic due to their resource-constrained nature, and integrating security on these devices is often overlooked. This has resulted in attackers having an increased incentive to target IoT devices. As the number of attacks possible on a network increases, it becomes more difficult for traditional intrusion detection systems (IDS) to cope with these attacks efficiently. In this paper, we highlight several machine learning (ML) methods such as k-nearest neighbour (KNN), support vector machine (SVM), decision tree (DT), naive Bayes (NB), random forest (RF), artificial neural network (ANN), and logistic regression (LR) that can be used in IDS. In this work, ML algorithms are compared for both binary and multi-class classification on Bot-IoT dataset. Based on several parameters such as accuracy, precision, recall, F1 score, and log loss, we experimentally compared the aforementioned ML algorithms. In the case of HTTP distributed denial-of-service (DDoS) attack, the accuracy of RF is 99%. Furthermore, other simulation results-based precision, recall, F1 score, and log loss metric reveal that RF outperforms on all types of attacks in binary classification. However, in multi-class classification, KNN outperforms other ML algorithms with an accuracy of 99%, which is 4% higher than RF.


title: Fixpoints and relative precompleteness

abstract: We study relative precompleteness in the context of the theory of numberings, and relate this to a notion of lowness. We introduce a notion of divisibility for numberings, and use it to show that for the class of divisible numberings, lowness and relative precompleteness coincide with being computable. We also study the complexity of Skolem functions arising from Arslanov's completeness criterion with parameters. We show that for suitably divisible numberings, these Skolem functions have the maximal possible Turing degree. In particular this holds for the standard numberings of the partial computable functions and the c.e. sets.


title: Skin marker-based subject-specific spinal alignment modeling: A feasibility study

abstract: Musculoskeletal models have the potential to improve diagnosis and optimize clinical treatment by predicting accurate outcomes on an individual basis. However, the subject-specific modeling of spinal alignment is often strongly simplified or is based on radiographic assessments, exposing subjects to unnecessary radiation. We therefore developed a novel skin marker-based approach for modeling subject-specific spinal alignment and evaluated its feasibility by comparing the predicted with the actual intervertebral joint (IVJ) locations/orientations (ground truth) using lateral-view radiographic images. Moreover, the predictive performance of the subject-specific models was evaluated by comparing the predicted L1/L2 spinal loads during various functional activities with in vivo measured data obtained from the OrthoLoad database. IVJ locations/orientations were predicted closer to ground truth as opposed to standard model scaling, with average location prediction errors of 0.99+/-0.68 cm on the frontal and 1.21+/-0.97 cm on the transverse axis as well as an average orientation prediction error of 4.74{\deg}+/-2.80{\deg}. Simulated spinal loads showed similar curve patterns but considerably larger values as compared to in vivo measured data. Differences in spinal loads between generic and subject-specific models become only apparent on an individual subject level. These results underline the feasibility of the proposed method and associated workflow for inter- and intra-subject investigations using musculoskeletal simulations. When implemented into standard model scaling workflows, it is expected to improve the accuracy of muscle activity and joint loading simulations, which is crucial for investigations of treatment effects or pathology-dependent deviations.


title: Teaching Turn-Taking Skills to Children with Autism using a Parrot-Like Robot

abstract: Robot Assisted Therapy is a new paradigm in many therapies such as the therapy of children with autism spectrum disorder. In this paper we present the use of a parrot-like robot as an assistive tool in turn taking therapy. The therapy is designed in the form of a card game between a child with autism and a therapist or the robot. The intervention was implemented in a single subject study format and the effect sizes for different turn taking variables are calculated. The results show that the child robot interaction had larger effect size than the child trainer effect size in most of the turn taking variables. Furthermore the therapist point of view on the proposed Robot Assisted Therapy is evaluated using a questionnaire. The therapist believes that the robot is appealing to children which may ease the therapy process. The therapist suggested to add other functionalities and games to let children with autism to learn more turn taking tasks and better generalize the learned tasks


title: On the well-posedness problem for the derivative nonlinear Schr\"odinger equation

abstract: We consider the derivative nonlinear Schr\"odinger equation in one space dimension, posed both on the line and on the circle. This model is known to be completely integrable and $L^2$-critical with respect to scaling. The first question we discuss is whether ensembles of orbits with $L^2$-equicontinuous initial data remain equicontinuous under evolution. We prove that this is true under the restriction $M(q)=\int |q|^2 < 4\pi$. We conjecture that this restriction is unnecessary. Further, we prove that the problem is globally well-posed for initial data in $H^{1/6}$ under the same restriction on $M$. Moreover, we show that this restriction would be removed by a successful resolution of our equicontinuity conjecture.


title: Stochastic social behavior coupled to COVID-19 dynamics leads to waves, plateaus and an endemic state

abstract: It is well recognized that population heterogeneity plays an important role in the spread of epidemics. While individual variations in social activity are often assumed to be persistent, i.e. constant in time, here we discuss the consequences of dynamic heterogeneity. By integrating the stochastic dynamics of social activity into traditional epidemiological models we demonstrate the emergence of a new long timescale governing the epidemic in broad agreement with empirical data. Our model captures multiple features of real-life epidemics such as COVID-19, including prolonged plateaus and multiple waves, which are transiently suppressed due to the dynamic nature of social activity. The existence of the long timescale due to the interplay between epidemic and social dynamics provides a unifying picture of how a fast-paced epidemic typically will transition to the endemic state.


title: A causal convolutional neural network for multi-subject motion modeling and generation

abstract: Inspired by the success of WaveNet in multi-subject speech synthesis, we propose a novel neural network based on causal convolutions for multi-subject motion modeling and generation. The network can capture the intrinsic characteristics of the motion of different subjects, such as the influence of skeleton scale variation on motion style. Moreover, after fine-tuning the network using a small motion dataset for a novel skeleton that is not included in the training dataset, it is able to synthesize high-quality motions with a personalized style for the novel skeleton. The experimental results demonstrate that our network can model the intrinsic characteristics of motions well and can be applied to various motion modeling and synthesis tasks.


title: Stellar Evolution and Tidal Dissipation in REBOUNDx

abstract: To study the post-main sequence evolution of the Solar system and exoplanetary systems more accurately and efficiently, we introduce two new features to REBOUNDx, an extended library for the N-body integrator REBOUND. The first is a convenient parameter interpolator for coupling different physics and integrators using numerical splitting schemes. The second implements a constant time lag model for tides without evolving spins. We demonstrate various uses of these features using stellar evolution data from MESA (Modules for Experiments in Stellar Astrophysics) as an example. The results of our tests agree with several studies in the literature on post-main sequence orbital evolution, and our convergence and performance studies respectively demonstrate our implementations' accuracy and efficiency. These additional effects are publicly available as of REBOUNDx's latest release.


title: Markov dynamics on the cone of discrete Radon measures

abstract: We start with a brief overview of the known facts about the spaces of discrete Radon measures those may be considered as generalizations of configuration spaces. Then we study three Markov dynamics on the spaces of discrete Radon measures: analogues of the contact model, of the Bolker--Dieckmann--Law--Pacala model, and of the Glauber-type dynamics. We show how the results obtained previously for the configuration spaces can be modified for the case of the spaces of discrete Radon measures.


title: GF-Flush: A GF(2) Algebraic Attack on Secure Scan Chains

abstract: Scan chains provide increased controllability and observability for testing digital circuits. The increased testability, however, can also be a source of information leakage for sensitive designs. The state-of-the-art defenses to secure scan chains apply dynamic keys to pseudo-randomly invert the scan vectors. In this paper, we pinpoint an algebraic vulnerability of these dynamic defenses that involves creating and solving a system of linear equations over the finite field GF(2). In particular, we propose a novel GF(2)-based flush attack that breaks even the most rigorous version of state-of-the-art dynamic defenses. Our experimental results demonstrate that our attack recovers the key as long as 500 bits in less than 7 seconds, the attack times are about one hundredth of state-of-the-art SAT based attacks on the same defenses. We then demonstrate how our attacks can be extended to scan chains compressed with Multiple-Input Signature Registers (MISRs).


title: On the Degeneracy of Spin Ice Graphs, and its Estimate via the Bethe Permanent

abstract: The concept of spin ice can be extended to a general graph. We study the degeneracy of spin ice graph on arbitrary interaction structures via graph theory. Via the mapping of spin ices to the Ising model, we clarify whether the inverse mapping is possible via a modified Krausz construction. From the gauge freedom of frustrated Ising systems, we derive exact, general results about frustration and degeneracy. We demonstrate for the first time that every spin ice graph, with the exception of the 1D Ising model, is degenerate. We then study how degeneracy scales in size, using the mapping between Eulerian trails and spin ice manifolds, and a permanental identity for the number of Eulerian orientations. We show that the Bethe permanent technique provides both an estimate and a lower bound to the frustration of spin ices on arbitrary graphs of even degree. While such technique can be used also to obtain an upper bound, we find that in all the examples we studied but one, another upper bound based on Schrijver inequality is tighter.


title: Voltage Controlled Spin-Orbit Torque Switching in W/CoFeB/MgO

abstract: Voltage control of magnetism and spintronics have been highly desirable, but rarely realized. In this work, we show voltage-controlled spin-orbit torque (SOT) switching in W/CoFeB/MgO films with perpendicular magnetic anisotropy (PMA) with voltage administered through SrTiO3 with a high dielectric constant. We show that a DC voltage can significantly lower PMA by 45%, reduce switching current by 23%, and increase the damping-like torque as revealed by the first and second-harmonic measurements. These are characteristics that are prerequisites for voltage-controlled and voltage-select SOT switching spintronic devices.


title: Simple Adaptive Estimation of Quadratic Functionals in Nonparametric IV Models

abstract: This paper considers adaptive, minimax estimation of a quadratic functional in a nonparametric instrumental variables (NPIV) model, which is an important problem in optimal estimation of a nonlinear functional of an ill-posed inverse regression with an unknown operator. We first show that a leave-one-out, sieve NPIV estimator of the quadratic functional can attain a convergence rate that coincides with the lower bound previously derived in Chen and Christensen [2018]. The minimax rate is achieved by the optimal choice of the sieve dimension (a key tuning parameter) that depends on the smoothness of the NPIV function and the degree of ill-posedness, both are unknown in practice. We next propose a Lepski-type data-driven choice of the key sieve dimension adaptive to the unknown NPIV model features. The adaptive estimator of the quadratic functional is shown to attain the minimax optimal rate in the severely ill-posed case and in the regular mildly ill-posed case, but up to a multiplicative $\sqrt{\log n}$ factor in the irregular mildly ill-posed case.


title: Optically tunable Mie-resonance VO2 nanoantennas for metasurfaces in the visible

abstract: Metasurfaces are ultrathin nanostructured surfaces that can allow arbitrary manipulation of light. Implementing dynamic tunability into their design could allow the optical functions of metasurfaces to be rapidly modified at will. The most pronounced and robust tunability of optical properties is provided by phase-change materials such as vanadium dioxide (VO2) and germanium antimony telluride (GST), but their implementations have been limited only to near-infrared wavelengths. Here, we demonstrate that VO2 nanoantennas with widely tunable Mie resonances can be utilized for designing tunable metasurfaces in the visible range. In contrast to the dielectric-metallic phase transition-induced tunability in previous demonstrations, we show that dielectric Mie resonances in VO2 nanoantennas offer remarkable scattering and extinction modulation depths (5-8 dB and 1-3 dB, respectively) for tunability in the visible. Moreover, these strong resonances are optically switchable using a continuous-wave laser. Our results establish VO2 nanostructures as low-loss building blocks of optically tunable metasurfaces.


title: AffectiveSpotlight: Facilitating the Communication of Affective Responses from Audience Members during Online Presentations

abstract: The ability to monitor audience reactions is critical when delivering presentations. However, current videoconferencing platforms offer limited solutions to support this. This work leverages recent advances in affect sensing to capture and facilitate communication of relevant audience signals. Using an exploratory survey (N = 175), we assessed the most relevant audience responses such as confusion, engagement, and head-nods. We then implemented AffectiveSpotlight, a Microsoft Teams bot that analyzes facial responses and head gestures of audience members and dynamically spotlights the most expressive ones. In a within-subjects study with 14 groups (N = 117), we observed that the system made presenters significantly more aware of their audience, speak for a longer period of time, and self-assess the quality of their talk more similarly to the audience members, compared to two control conditions (randomly-selected spotlight and default platform UI). We provide design recommendations for future affective interfaces for online presentations based on feedback from the study.


title: Semiparametric point process modeling of blinking artifacts in PALM

abstract: Photoactivated localization microscopy (PALM) is a powerful imaging technique for characterization of protein organization in biological cells. Due to the stochastic blinking of fluorescent probes, and camera discretization effects, each protein gives rise to a cluster of artificial observations. These blinking artifacts are an obstacle for quantitative analysis of PALM data, and tools for their correction are in high demand. We develop the Independent Blinking Cluster point process (IBCpp) family of models, which is suited for modeling of data from single-molecule localization microscopy modalities, and we present results on the mark correlation function. We then construct the PALM-IBCpp - a semiparametric IBCpp tailored for PALM data, and we describe a procedure for estimation of parameters, which can be used without parametric assumptions on the spatial organization of proteins. Our model is validated on nuclear pore complex reference data, where the ground truth was accurately recovered, and we demonstrate how the estimated blinking parameters can be used to perform a blinking corrected test for protein clustering in a cell expressing the adaptor protein LAT. Finally, we consider simulations with varying degrees of blinking and protein clustering to shed light on the expected performance in a range of realistic settings.


title: Measures of pion and kaon structure from generalised parton distributions

abstract: Pion and kaon structural properties provide insights into the emergence of mass within the Standard Model and attendant modulations by the Higgs boson. Novel expressions of these effects, in impact parameter space and in mass and pressure profiles, are exposed via $\pi$ and $K$ generalised parton distributions, built using the overlap representation from light-front wave functions constrained by one-dimensional valence distribution functions that describe available data. Notably, e.g. $K$ pressure profiles are spatially more compact than $\pi$ profiles and both achieve near-core pressures of similar magnitude to that found in neutron stars.


title: Support varieties over skew complete intersections via derived braided Hochschild cohomology

abstract: In this article we study a theory of support varieties over a skew complete intersection $R$, i.e. a skew polynomial ring modulo an ideal generated by a sequence of regular normal elements. We compute the derived braided Hochschild cohomology of $R$ relative to the skew polynomial ring and show its action on $\mathrm{Ext}_R(M,N)$ is noetherian for finitely generated $R$-modules $M$ and $N$ respecting the braiding of $R$. When the parameters defining the skew polynomial ring are roots of unity we use this action to define a support theory. In this setting applications include a proof of the Generalized Auslander-Reiten Conjecture and that $R$ possesses symmetric complexity.


title: From Geometry to Topology: Inverse Theorems for Distributed Persistence

abstract: What is the "right" topological invariant of a large point cloud X? Prior research has focused on estimating the full persistence diagram of X, a quantity that is very expensive to compute, unstable to outliers, and far from a sufficient statistic. We therefore propose that the correct invariant is not the persistence diagram of X, but rather the collection of persistence diagrams of many small subsets. This invariant, which we call "distributed persistence," is perfectly parallelizable, more stable to outliers, and has a rich inverse theory. The map from the space of point clouds (with the quasi-isometry metric) to the space of distributed persistence invariants (with the Hausdorff-Bottleneck distance) is a global quasi-isometry. This is a much stronger property than simply being injective, as it implies that the inverse of a small neighborhood is a small neighborhood, and is to our knowledge the only result of its kind in the TDA literature. Moreover, the quasi-isometry bounds depend on the size of the subsets taken, so that as the size of these subsets goes from small to large, the invariant interpolates between a purely geometric one and a topological one. Lastly, we note that our inverse results do not actually require considering all subsets of a fixed size (an enormous collection), but a relatively small collection satisfying certain covering properties that arise with high probability when randomly sampling subsets. These theoretical results are complemented by two synthetic experiments demonstrating the use of distributed persistence in practice.


title: Probabilistic Data with Continuous Distributions

abstract: Statistical models of real world data typically involve continuous probability distributions such as normal, Laplace, or exponential distributions. Such distributions are supported by many probabilistic modelling formalisms, including probabilistic database systems. Yet, the traditional theoretical framework of probabilistic databases focusses entirely on finite probabilistic databases. Only recently, we set out to develop the mathematical theory of infinite probabilistic databases. The present paper is an exposition of two recent papers which are cornerstones of this theory. In (Grohe, Lindner; ICDT 2020) we propose a very general framework for probabilistic databases, possibly involving continuous probability distributions, and show that queries have a well-defined semantics in this framework. In (Grohe, Kaminski, Katoen, Lindner; PODS 2020) we extend the declarative probabilistic programming language Generative Datalog, proposed by (B\'ar\'any et al.~2017) for discrete probability distributions, to continuous probability distributions and show that such programs yield generative models of continuous probabilistic databases.


title: Independent Hyperplanes in Oriented Paving Matroids

abstract: In 1993, Csima and Sawyer proved that in a non-pencil arrangement of n pseudolines, there are at least $\frac{6}{13}n$ simple points of intersection. Since pseudoline arrangements are the topological representations of reorientation classes of oriented matroids of rank $3$, in this paper, we will use this result to prove by induction that an oriented paving matroid of rank $r \ge 3$ on $n$ elements, where $n \geq 5+ r$, has at least $\frac{12}{13(r-1)} \binom{n}{r-2}$ independent hyperplanes, yielding a new necessary condition for a paving matroid to be orientable.


title: Magic Conditions for Multiple Rotational States of Bialkali Molecules in Optical Lattices

abstract: We investigate magic-wavelength trapping of ultracold bialkali molecules in the vicinity of weak optical transitions from the vibrational ground state of the X$^1\Sigma^+$ potential to low-lying rovibrational states of the b$^3\Pi_0$ potential, focussing our discussion on the $^{87}$Rb$^{133}$Cs molecule in a magnetic field of $B=181\,$G. We show that a frequency window exists between two nearest neighbor vibrational poles in the dynamic polarizability where the trapping potential is "near magic" for multiple rotational states simultaneously. We show that the addition of a modest DC electric field of $E=0.13\,\text{kV}/\text{cm}$ leads to an exact magic-wavelength trap for the lowest three rotational states at a angular-frequency detuning of $\Delta_{v'=0} = 2\pi\times 218.22$\,GHz from the X$^1\Sigma^+ (v=0, J=0)\rightarrow$ b$^3\Pi_0 (v'=0, J=1)$ transition. We derive a set of analytical criteria that must be fulfilled to ensure the existence of such magic frequency windows and present an analytic expression for the position of the frequency window in terms of a set of experimentally measurable parameters. These results should inform future experiments requiring long coherence times on multiple rotational transitions in ultracold polar molecules.


title: Universal field dependence of magnetic resonance near zero frequency

abstract: Magnetic resonance is a widely-established phenomenon that probes magnetic properties such as magnetic damping and anisotropy. Even though the typical resonance frequency of a magnet ranges from gigahertz to terahertz, experiments also report the resonance near zero frequency in a large class of magnets. Here we revisit this phenomenon by analyzing the symmetry of the system and find that the resonance frequency ($\omega$) follows a universal power law $\omega \varpropto |H-H_c|^p$, where $H_c$ is the critical field at which the resonance frequency is zero. When the magnet preserves the rotational symmetry around the external field ($H$), $p = 1$. Otherwise, $p=1/2$. The magnon excitations are gapped above $H_c$, gapless at $H_c$ and gapped again below $H_c$. The zero frequency is often accompanied by a reorientation transition in the magnetization. For the case that $p=1/2$, this transition is described by a Landau theory for second-order phase transitions. We further show that the spin current driven by thermal gradient and spin-orbit effects can be significantly enhanced when the resonance frequency is close to zero, which can be measured electrically by converting the spin current into electric signals. This may provide an experimentally accessible way to characterize the critical field. Our findings provide a unified understanding of the magnetization dynamics near the critical field, and may, furthermore, inspire the study of magnon transport near magnetic transitions.


title: Jones Matrix Characterization of Optical Elements via Evolutionary Algorithms

abstract: Jones calculus provides a robust and straightforward method to characterize polarized light and polarizing optical systems using two-element vectors (Jones vectors) and $2 \times 2$ matrices (Jones matrices). Jones matrices are used to determine the retardance and diattenuation introduced by an optical element or a sequence of elements. Moreover, they are the tool of choice to study optical geometric phases. However, the current sampling method for characterizing the Jones matrix of an optical element is inefficient, since the search space of the problem is in the realm of the real numbers and so applying a general sampling method is time-consuming. In this study, we present an initial approach for solving the problem of finding the eigenvectors that characterize the Jones matrix of a homogeneous optical element through Evolutionary Algorithms (EAs). We evaluate the analytical performance of an EA with a Polynomial Mutation operator and a Genetic Algorithm (GA) with a Simulated Binary crossover operator and a Polynomial Mutation operator, and compare the results with those obtained through a general sampling method. The results show that both the EA and the GA out-performed a general sampling method of 6,000 measurements, by requiring in average 103 and 188 fitness functions measurements respectively, while having a perfect rate of convergence.


title: Combining pre-trained language models and structured knowledge

abstract: In recent years, transformer-based language models have achieved state of the art performance in various NLP benchmarks. These models are able to extract mostly distributional information with some semantics from unstructured text, however it has proven challenging to integrate structured information, such as knowledge graphs into these models. We examine a variety of approaches to integrate structured knowledge into current language models and determine challenges, and possible opportunities to leverage both structured and unstructured information sources. From our survey, we find that there are still opportunities at exploiting adapter-based injections and that it may be possible to further combine various of the explored approaches into one system.


title: Spectral zeta-Functions and zeta-Regularized Functional Determinants for Regular Sturm-Liouville Operators

abstract: The principal aim in this paper is to employ a recently developed unified approach to the computation of traces of resolvents and $\zeta$-functions to efficiently compute values of spectral $\zeta$-functions at positive integers associated to regular (three-coefficient) self-adjoint Sturm--Liouville differential expressions $\tau$. Depending on the underlying boundary conditions, we express the $\zeta$-function values in terms of a fundamental system of solutions of $\tau y = z y$ and their expansions about the spectral point $z=0$. Furthermore, we give the full analytic continuation of the $\zeta$-function through a Liouville transformation and provide an explicit expression for the $\zeta$-regularized functional determinant in terms of a particular set of this fundamental system of solutions. An array of examples illustrating the applicability of these methods is provided, including regular Schr\"{o}dinger operators with zero, piecewise constant, and a linear potential on a compact interval.


title: Asymptotics of $p$-torsion subgroup sizes in class groups of monogenized cubic fields

abstract: Bhargava, Hanke, and Shankar have recently shown that the asymptotic average $2$-torsion subgroup size of the family of class groups of monogenized cubic fields with positive and negative discriminants is $3/2$ and $2$, respectively. In this paper, we provide strong computational evidence for these asymptotes. We then develop a pair of novel conjectures that predicts, for $p$ prime, the asymptotic average $p$-torsion subgroup size in class groups of monogenized cubic fields.


title: Growth and Collapse of an Isolated Bubble Driven by a Single Negative Histotripsy Cycle in Agarose Gel: Stress, Strain, and Strain Rate Fields

abstract: Histotripsy relies on cavitation to mechanically homogenize soft tissue. There is strong evidence that the high stresses, strains, and strain rates developed as bubbles grow and collapse contribute to this tissue homogenization. While such stresses and strains have been examined computationally in model systems with assumed constitutive models (e.g., finite-deformation Neo-Hookean model) and viscoelastic properties determined under quasi-static conditions, recent studies proposed that the Quadratic Law Kelvin-Voigt (QLKV) constitutive model, which additionally accounts for strain stiffening, more accurately represents the viscoelastic response of soft materials subjected to cavitation; this model has also been used to infer viscoelastic properties at high rates. In this work, we use the QLKV model and these properties to calculate the time-dependent stress, strain, and strain rate fields produced during the growth and collapse of individual bubbles subjected to a histotripsy-relevant pressure waveform in agarose gels of 0.3~\% and 1.0~\% concentration and corresponding to actual (past) experiments. We find that, as the gel concentration is increased, strain stiffening manifests in larger elastic stresses and compressive stresses extending into the collapse phase, particularly for the 1.0~\% concentration gel. As a result, the duration of the collapse phase also increases. In comparison with the conventional Neo-Hookean model, the compressive stress has a larger magnitude, extends farther into the surrounding medium, and shows an increased departure from growth/collapse symmetry close to the bubble; all of these effects are magnified in the stiffer gel.


title: Atomic Clocks in Space: A Search for Rubidium and Cesium Masers in M- and L-Dwarfs

abstract: I searched for the ground state 6.8 and 9.2 GHz hyperfine transitions of rubidium and cesium toward M- and L-dwarfs that show Rb and Cs optical resonance lines. The optical lines can pump the hyperfine transitions, potentially forming masers. These spin-flip transitions of Rb and Cs are the principal transitions used in atomic clocks (the $^{133}$Cs hyperfine transition defines the second). If they are detected in stellar atmospheres, these transitions would provide exceptionally precise clocks that can be used as accelerometers, as exoplanet detectors, as probes of the predictions of general relativity, as probes of light propagation effects, and as a means to do fundamental physics with telescopes. Observations of 21 M- and L-dwarfs, however, show no evidence for Rb or Cs maser action, and a previous survey of giant stars made no Rb maser detections.


title: Putting gradual types to work

abstract: In this paper, we describe our experience incorporating gradual types in a statically typed functional language with Hindley-Milner style type inference. Where most gradually typed systems aim to improve static checking in a dynamically typed language, we approach it from the opposite perspective and promote dynamic checking in a statically typed language. Our approach provides a glimpse into how languages like SML and OCaml might handle gradual typing. We discuss our implementation and challenges faced -- specifically how gradual typing rules apply to our representation of composite and recursive types. We review the various implementations that add dynamic typing to a statically typed language in order to highlight the different ways of mixing static and dynamic typing and examine possible inspirations while maintaining the gradual nature of our type system. This paper also discusses our motivation for adding gradual types to our language, and the practical benefits of doing so in our industrial setting.


title: A sub-Neptune and a non-transiting Neptune-mass companion unveiled by ESPRESSO around the bright late-F dwarf HD 5278 (TOI-130)

abstract: [Abridged] We exploit the extreme radial velocity (RV) precision of the ultra-stable echelle spectrograph ESPRESSO on the VLT to unveil the physical properties of the transiting sub-Neptune TOI-130 b, uncovered by TESS orbiting the nearby, bright, late F-type star HD 5278 (TOI-130) with a period $P_{\rm b}=14.3$. We use 43 ESPRESSO high-resolution spectra and broad-band photometry information to derive accurate stellar atmospheric and physical parameters of HD 5278. We exploit the TESS light curve (LC) and spectroscopic diagnostics to gauge the impact of stellar activity on the ESPRESSO RVs. We perform a joint ESPRESSO RVs + TESS LC analysis using fully Bayesian frameworks to determine the system parameters. The updated stellar parameters of HD 5278 are T$_\mathrm{eff}=6203\pm64$ K, $\log g =4.50\pm0.11$ dex, [Fe/H]=$-0.12\pm0.04$ dex, M$_\star=1.126_{-0.035}^{+0.036}$ M$_\odot$ and R$_\star=1.194_{-0.016}^{+0.017}$ R$_\odot$. We determine HD 5278 b's mass and radius to be $M_{\rm b} = 7.8_{-1.4}^{+1.5}$ M$_\oplus$ and $R_{\rm b} = 2.45\pm0.05$ R$_\oplus$. The derived mean density, $\varrho_{\rm b} = 2.9_{-0.5}^{+0.6}$ g cm$^{-3}$, is consistent with a bulk composition with a substantial ($\sim30\%$) water mass fraction and a gas envelope comprising $\sim17\%$ of the measured radius. Given the host brightness and irradiation levels, HD 5278 b is one of the best targets orbiting G-F primaries for follow-up atmospheric characterization measurements with HST and JWST. We discover a second, non-transiting companion in the system, with a period $P_{\rm c}=40.87_{-0.17}^{+0.18}$ days and a minimum mass $M_{\rm c}\sin i_{\rm c} =18.4_{-1.9}^{+1.8}$ M$_\oplus$. We study emerging trends in the growing population of transiting sub-Neptunes, and provide statistical evidence for a low occurrence of close-in, $10-15$ M$_\oplus$ companions around G-F primaries with $T_\mathrm{eff}\gtrsim5500$ K.


title: Strong Field Ionization of Water: Nuclear Dynamics Revealed by Varying the Pulse Duration

abstract: Polyatomic molecules in strong laser fields can undergo substantial nuclear motion within tens of femtoseconds. Ion imaging methods based on dissociation or Coulomb explosion therefore have difficulty faithfully recording the geometry dependence of the field ionization that initiates the dissociation process. Here we compare the strong-field double ionization and subsequent dissociation of water (both H$_2$O and D$_2$O) in 10-fs and 40-fs 800-nm laser pulses. We find that 10-fs pulses turn off before substantial internuclear motion occurs, whereas rapid internuclear motion can take place during the double ionization process for 40-fs pulses. The short-pulse measurements are consistent with a simple tunnel ionization picture, whose predictions help interpret the motion observed in the long-pulse measurements.


title: Existence and uniqueness for non-Markovian triangular quadratic BSDEs

abstract: We prove the existence and uniqueness of solutions to a class of quadratic BSDE systems which we call triangular quadratic. Our results generalize several existing results about diagonally quadratic BSDEs in the non-Markovian setting. As part of our analysis, we obtain new results about linear BSDEs with unbounded coefficients, which may be of independent interest. Through a non-uniqueness example, we answer a "crucial open question" raised by Harter and Richou by showing that the stochastic exponential of an n x n matrix-valued BMO martingale need not satisfy a reverse H\"older inequality.


title: Transparent and flexible high-power supercapacitor based on carbon nanotube fibre aerogels

abstract: In this work, we report on the fabrication of continuous transparent and flexible supercapacitors by depositing a CNT network onto a polymer electrolyte membrane directly from an aerogel of ultra-long CNTs produced floating in the gas phase. The supercapacitors combine record power density of $1370 kW kg^{-1}$ at high transmittance ($ca. 70%$), high electrochemical stability during extended cycling ($>94%$ capacitance retention over $20 000 cycles$) as well as against repeated $180{\deg}$ flexural deformation. They represent a significant enhancement of 1-3 orders of magnitude compared to the prior-art transparent supercapacitors based on graphene, CNTs, and rGO. These features mainly arise from the exceptionally long length of the CNTs, which makes the material behave as a bulk conductor instead of an aspect ratio-limited percolating network, even for electrodes with $>90%$ transparency. The electrical and capacitive figures-of-merit for the transparent conductor are $FoMe = 2.7$, and $FoMc = 0.46 F S^{-1} cm^{-2}$ respectively


title: Search for $K^+$ decays to a muon and invisible particles

abstract: The NA62 experiment at CERN reports searches for $K^+\to\mu^+N$ and $K^+\to\mu^+\nu X$ decays, where $N$ and $X$ are massive invisible particles, using the 2016-2018 data set. The $N$ particle is assumed to be a heavy neutral lepton, and the results are expressed as upper limits of ${\cal O}(10^{-8})$ of the neutrino mixing parameter $|U_{\mu4}|^2$ for $N$ masses in the range 200-384 MeV/$c^2$ and lifetime exceeding 50 ns. The $X$ particle is considered a scalar or vector hidden sector mediator decaying to an invisible final state, and upper limits of the decay branching fraction for $X$ masses in the range 10-370 MeV/$c^2$ are reported for the first time, ranging from ${\cal O}(10^{-5})$ to ${\cal O}(10^{-7})$. An improved upper limit of $1.0\times 10^{-6}$ is established at 90% CL on the $K^+\to\mu^+\nu\nu\bar\nu$ branching fraction.


title: Evaluating Complex Queries on Streaming Graphs

abstract: We study the problem of evaluating persistent queries over streaming graphs in a principled fashion. These queries need to be evaluated over unbounded and very high speed graph streams. We define a streaming graph data model and query model incorporating navigational queries, subgraph queries and paths as first-class citizens. To support this full-fledged query model we develop a streaming graph algebra that describes the precise semantics of persistent graph queries with their complex constructs. We present transformation rules and describe query formulation and plan generation for persistent graph queries over streaming graphs. Our implementation of a streaming graph query processor shows the feasibility of our approach and allows us to gauge the high performance gains obtained for query processing over streaming graphs.


title: New Formulations of Ambiguous Volatility with an Application to Optimal Dynamic Contracting

abstract: I introduce novel preference formulations which capture aversion to ambiguity about unknown and potentially time-varying volatility. I compare these preferences with Gilboa and Schmeidler's maxmin expected utility as well as variational formulations of ambiguity aversion. The impact of ambiguity aversion is illustrated in a simple static model of portfolio choice, as well as a dynamic model of optimal contracting under repeated moral hazard. Implications for investor beliefs, optimal design of corporate securities, and asset pricing are explored.


title: Magnetoresistive Sensor Detectivity: A Comparative Analysis

abstract: We report on the noise performance characteristics of magnetic sensors using both magnetic tunnel junction (MTJ) and giant magnetoresistance (GMR) elements. Each sensor studied has a notably different noise and detectivity. Of the sensors we measured, those based on GMR multilayers have the lowest noise and detectivity. However, the GMR sensor also has a significantly smaller linear range. To make a direct comparison between sensors we scale the linear operating ranges of each sensor to be the same. This is the phenomenological equivalent of modifying the flux concentration. Upon scaling the low frequency detectivity of the TMR sensors becomes essentially equal to that of the GMR sensor. Using the scaling approach we are able to place the detectivity in the context of other key parameters, namely size and power consumption. Lastly, we use this technique to examine the upper limit for magnetoresistive sensor performance based on a notional MTJ sensor using present record setting TMR values.


title: The Initial Degree of Symbolic Powers of Ideals of Fermat Configuration of Points

abstract: Let $n \ge 2$ be an integer and consider the defining ideal of the Fermat configuration of points in $\mathbb{P}^2$: $I_n=(x(y^n-z^n),y(z^n-x^n),z(x^n-y^n)) \subset R=\mathbb{C}[x,y,z]$. In this paper, we compute explicitly the least degree of generators of its symbolic powers in all unknown cases. As direct applications, we easily verify Chudnovsky's Conjecture, Demailly's Conjecture and Harbourne-Huneke Containment problem as well as calculating explicitly the Waldschmidt constant and (asymptotic) resurgence number.


title: Observation of the Decrease of Larmor Tunneling Times with Lower Incident Energy

abstract: How much time does a tunneling particle spend in a barrier? A Larmor clock, one proposal to answer this question, measures the interaction between the particle and the barrier region using an auxiliary degree of freedom of the particle to clock the dwell time inside the barrier. We report on precise Larmor time measurements of ultra-cold $^{87}$Rb atoms tunneling through an optical barrier, which confirm longstanding predictions of tunneling times. We observe that atoms generally spend less time tunneling through higher barriers and that this time decreases for lower energy particles. For the lowest measured incident energy, at least $90\%$ of transmitted atoms tunneled through the barrier, spending an average of $0.59(2)$ms inside. This is $0.11(3)$ms faster than atoms traversing the same barrier with energy close to the barrier's peak and $0.21(3)$ms faster than when the atoms traverse a barrier with $23\%$ less energy.


title: Flow-Sensory Contact Electrification of Graphene

abstract: All-electronic interrogation of biofluid flow velocity by sensors incorporated in ultra-low-power or self-sustained systems offers the promise of enabling multifarious emerging research and applications. Electrical sensors based on nanomaterials are of high spatiotemporal resolution and exceptional sensitivity to external flow stimulus and easily integrated and fabricated using scalable techniques. But existing nano-based electrical flow-sensing technologies remain lacking in precision and stability and are typically only applicable to simple aqueous solutions or liquid/gas dual-phase mixtures, making them unsuitable for monitoring low-flow (~micrometer/second) yet important characteristics of continuous biofluids (e.g., hemorheological behaviors in microcirculation). Here we show that monolayer-graphene single microelectrodes harvesting charge from continuous aqueous flow provide an ideal flow sensing strategy: Our devices deliver over six months stability and sub-micrometer/second resolution in real-time quantification of whole-blood flows with multiscale amplitude-temporal characteristics in a microfluidic chip. The flow transduction is enabled by low-noise charge transfer at graphene/water interface in response to flow-sensitive rearrangement of the interfacial electrical double layer. Our results demonstrate the feasibility of using a graphene-based self-powered strategy for monitoring biofluid flow velocity with key performance metrics orders of magnitude higher than other electrical approaches.


title: Hamiltonian Gotay-Nester-Hinds analysis of the parametrized unimodular extension of the Holst action

abstract: We give a detailed account of the Hamiltonian GNH analysis of the parametrized unimodular extension of the Holst action. The purpose of the paper is to derive, through the clear geometric picture furnished by the GNH method, a simple Hamiltonian formulation for this model and explain why it is difficult to arrive at it in other approaches. We will also show how to take advantage of the field equations to anticipate the simple form of the constraints that we find in the paper.


title: The Bootstrap for Network Dependent Processes

abstract: This paper focuses on the bootstrap for network dependent processes under the conditional $\psi$-weak dependence. Such processes are distinct from other forms of random fields studied in the statistics and econometrics literature so that the existing bootstrap methods cannot be applied directly. We propose a block-based approach and a modification of the dependent wild bootstrap for constructing confidence sets for the mean of a network dependent process. In addition, we establish the consistency of these methods for the smooth function model and provide the bootstrap alternatives to the network heteroskedasticity-autocorrelation consistent (HAC) variance estimator. We find that the modified dependent wild bootstrap and the corresponding variance estimator are consistent under weaker conditions relative to the block-based method, which makes the former approach preferable for practical implementation.


title: Third-order ladder operators, generalized Okamoto and exceptional orthogonal polynomials

abstract: We extend and generalize the construction of Sturm-Liouville problems for a family of Hamiltonians constrained to fulfill a third-order shape-invariance condition and focusing on the "$-2x/3$" hierarchy of solutions to the fourth Painlev\'e transcendent. Such a construction has been previously addressed in the literature for some particular cases but we realize it here in the most general case. The corresponding potential in the Hamiltonian operator is a rationally extended oscillator defined in terms of the conventional Okamoto polynomials, from which we identify three different zero-modes constructed in terms of the generalized Okamoto polynomials. The third-order ladder operators of the system reveal that the complete set of eigenfunctions is decomposed as a union of three disjoint sequences of solutions, generated from a set of three-term recurrence relations. We also identify a link between the eigenfunctions of the Hamiltonian operator and a special family of exceptional Hermite polynomial.


title: Fourier multipliers for Triebel-Lizorkin spaces on compact Lie groups

abstract: We investigate the boundedness of Fourier multipliers on a compact Lie group when acting on Triebel-Lizorkin spaces. Criteria are given in terms of the H\"ormander-Mihlin-Marcinkiewicz condition. In our analysis, we use the difference structure of the unitary dual of a compact Lie group. Our results cover the sharp H\"ormander-Mihlin theorem on Lebesgue spaces and also other historical results on the subject.


title: Phase transitions in GLSMs and defects

abstract: In this paper, we construct defects (domain walls) that connect different phases of two-dimensional gauged linear sigma models (GLSMs), as well as defects that embed those phases into the GLSMs. Via their action on boundary conditions these defects give rise to functors between the D-brane categories, which respectively describe the transport of D-branes between different phases, and embed the D-brane categories of the phases into the category of D-branes of the GLSMs.


title: Byzantine Fault-Tolerance in Peer-to-Peer Distributed Gradient-Descent

abstract: We consider the problem of Byzantine fault-tolerance in the peer-to-peer (P2P) distributed gradient-descent method -- a prominent algorithm for distributed optimization in a P2P system. In this problem, the system comprises of multiple agents, and each agent has a local cost function. In the fault-free case, when all the agents are honest, the P2P distributed gradient-descent method allows all the agents to reach a consensus on a solution that minimizes their aggregate cost. However, we consider a scenario where a certain number of agents may be Byzantine faulty. Such faulty agents may not follow an algorithm correctly, and may share arbitrary incorrect information to prevent other non-faulty agents from solving the optimization problem. In the presence of Byzantine faulty agents, a more reasonable goal is to allow all the non-faulty agents to reach a consensus on a solution that minimizes the aggregate cost of all the non-faulty agents. We refer to this fault-tolerance goal as $f$-resilience where $f$ is the maximum number of Byzantine faulty agents in a system of $n$ agents, with $f < n$. Most prior work on fault-tolerance in P2P distributed optimization only consider approximate fault-tolerance wherein, unlike $f$-resilience, all the non-faulty agents' compute a minimum point of a non-uniformly weighted aggregate of their cost functions. We propose a fault-tolerance mechanism that confers provable $f$-resilience to the P2P distributed gradient-descent method, provided the non-faulty agents satisfy the necessary condition of $2f$-redundancy, defined later in the paper. Moreover, compared to prior work, our algorithm is applicable to a larger class of high-dimensional convex distributed optimization problems.


title: Lippmann-Schwinger-Lanczos algorithm for inverse scattering problems

abstract: Data-driven reduced order models (ROMs) are combined with the Lippmann-Schwinger integral equation to produce a direct nonlinear inversion method. The ROM is viewed as a Galerkin projection and is sparse due to Lanczos orthogonalization. Embedding into the continuous problem, a data-driven internal solution is produced. This internal solution is then used in the Lippmann-Schwinger equation, thus making further iterative updates unnecessary. We show numerical experiments for spectral domain domain data for which our inversion is far superior to the Born inversion and works as well as when the true internal solution is known.


title: Experimentation for Homogenous Policy Change

abstract: When the Stable Unit Treatment Value Assumption (SUTVA) is violated and there is interference among units, there is not a uniquely defined Average Treatment Effect (ATE), and alternative estimands may be of interest, among them average unit-level differences in outcomes under different homogeneous treatment policies. We term this target the Homogeneous Assignment Average Treatment Effect (HAATE). We consider approaches to experimental design with multiple treatment conditions under partial interference and, given the estimand of interest, we show that difference-in-means estimators may perform better than correctly specified regression models in finite samples on root mean squared error (RMSE). With errors correlated at the cluster level, we demonstrate that two-stage randomization procedures with intra-cluster correlation of treatment strictly between zero and one may dominate one-stage randomization designs on the same metric. Simulations demonstrate performance of this approach; an application to online experiments at Facebook is discussed.


title: General conditions for universality of Quantum Hamiltonians

abstract: Recent work has demonstrated the existence of universal Hamiltonians - simple spin lattice models that can simulate any other quantum many body system to any desired level of accuracy. Until now proofs of universality have relied on explicit constructions, tailored to each specific family of universal Hamiltonians. In this work we go beyond this approach, and completely classify the simulation ability of quantum Hamiltonians by their complexity classes. We do this by deriving necessary and sufficient complexity theoretic conditions characterising universal quantum Hamiltonians. Although the result concerns the theory of analogue Hamiltonian simulation - a promising application of near-term quantum technology - the proof relies on abstract complexity theoretic concepts and the theory of quantum computation. As well as providing simplified proofs of previous Hamiltonian universality results, and offering a route to new universal constructions, the results in this paper give insight into the origins of universality. For example, finally explaining the previously noted coincidences between families of universal Hamiltonian and classes of Hamiltonians appearing in complexity theory.


title: Observation of ultra-slow shock waves in a tunable magnetic lattice

abstract: The combination of fast propagation speeds and highly localized nature has hindered the direct observation of the evolution of shock waves at the molecular scale. To address this limitation, an experimental system is designed by tuning a one-dimensional magnetic lattice to evolve benign wave forms into shock waves at observable spatial and temporal scales, thus serving as a 'magnifying glass' to illuminate shock processes. An accompanying analysis confirms that the formation of strong shocks is fully captured. The exhibited lack of a steady state induced by indefinite expansion of a disordered transition zone points to the absence of local thermodynamic equilibrium, and resurfaces lingering questions on the validity of continuum assumptions in presence of strong shocks.


title: Topological transitions during grain growth on a finite element mesh

abstract: The topological transitions that occur to the grain boundary network during grain growth in a material with uniform grain boundary energies are believed to be known. The same is not true for more realistic materials, since more general grain boundary energies in principle allow many more viable grain boundary configurations. A simulation of grain growth in such a material therefore requires a procedure to enumerate all possible topological transitions and select the most energetically favorable one. Such a procedure is developed and implemented here for a microstructure represented by a volumetric finite element mesh. As a specific example, all possible transitions for a typical configuration with five grains around a junction point are enumerated, and some exceptional transitions are found to be energetically similar to the conventional ones even for a uniform boundary energy. A general discrete formulation to calculate grain boundary velocities is used to simulate grain growth for an example microstructure. The method is implemented as a C++ library based on SCOREC, an open source massively parallelizable library for finite element simulations with adaptive meshing.


title: Position, Padding and Predictions: A Deeper Look at Position Information in CNNs

abstract: In contrast to fully connected networks, Convolutional Neural Networks (CNNs) achieve efficiency by learning weights associated with local filters with a finite spatial extent. An implication of this is that a filter may know what it is looking at, but not where it is positioned in the image. In this paper, we first test this hypothesis and reveal that a surprising degree of absolute position information is encoded in commonly used CNNs. We show that zero padding drives CNNs to encode position information in their internal representations, while a lack of padding precludes position encoding. This gives rise to deeper questions about the role of position information in CNNs: (i) What boundary heuristics enable optimal position encoding for downstream tasks?; (ii) Does position encoding affect the learning of semantic representations?; (iii) Does position encoding always improve performance? To provide answers, we perform the largest case study to date on the role that padding and border heuristics play in CNNs. We design novel tasks which allow us to quantify boundary effects as a function of the distance to the border. Numerous semantic objectives reveal the effect of the border on semantic representations. Finally, we demonstrate the implications of these findings on multiple real-world tasks to show that position information can both help or hurt performance.


title: The simplicity of the C*-algebras associated to arbitrary labeled spaces

abstract: In this paper, we consider the simplicity of the C*-algebra associated to an arbitrary weakly left-resolving labeled space (E, L, E), where E is the smallest non-degenerate accommodating set. We classify all gauge-invariant ideals of C*(E, L, E) and characterize minimality of (E, L, E) in terms of ideal structure of C*(E, L, E). Using these results, we prove simplicity results for C*-algebras associated to arbitrary labeled spaces.


title: Geodesic length and shifted weights in first-passage percolation

abstract: We study first-passage percolation through related optimization problems over paths of restricted length. The path length variable is in duality with a shift of the weights. This puts into a convex duality framework old observations about the convergence of the normalized Euclidean length of geodesics due to Hammersley and Welsh, Smythe and Wierman, and Kesten, and leads to new results about geodesic length and the regularity of the shape function as a function of the weight shift. For points far enough away from the origin, the ratio of the geodesic length and the $\ell^1$ distance to the endpoint is uniformly bounded away from one. The shape function is a strictly concave function of the weight shift. Atoms of the weight distribution generate singularities, that is, points of nondifferentiability, in this function. We generalize to all distributions, directions and dimensions an old singularity result of Steele and Zhang for the planar Bernoulli case. When the weight distribution has two or more atoms, a dense set of shifts produce singularities. The results come from a combination of the convex duality, the shape theorems of the different first-passage optimization problems, and modification arguments.


title: Hyperpolarisation of external nuclear spins using nitrogen-vacancy centre ensembles

abstract: The nitrogen-vacancy (NV) centre in diamond has emerged as a candidate to non-invasively hyperpolarise nuclear spins in molecular systems to improve the sensitivity of nuclear magnetic resonance (NMR) experiments. Several promising proof of principle experiments have demonstrated small-scale polarisation transfer from single NVs to hydrogen spins outside the diamond. However, the scaling up of these results to the use of a dense NV ensemble, which is a necessary prerequisite for achieving realistic NMR sensitivity enhancement, has not yet been demonstrated. In this work, we present evidence for a polarising interaction between a shallow NV ensemble and external nuclear targets over a micrometre scale, and characterise the challenges in achieving useful polarisation enhancement. In the most favourable example of the interaction with hydrogen in a solid state target, a maximum polarisation transfer rate of $\approx 7500$ spins per second per NV is measured, averaged over an area containing order $10^6$ NVs. Reduced levels of polarisation efficiency are found for liquid state targets, where molecular diffusion limits the transfer. Through analysis via a theoretical model, we find that our results suggest implementation of this technique for NMR sensitivity enhancement is feasible following realistic diamond material improvements.


title: The Optimal Dynamic Treatment Rule SuperLearner: Considerations, Performance, and Application

abstract: The optimal dynamic treatment rule (ODTR) framework offers an approach for understanding which kinds of patients respond best to specific treatments -- in other words, treatment effect heterogeneity. Recently, there has been a proliferation of methods for estimating the ODTR. One such method is an extension of the SuperLearner algorithm -- an ensemble method to optimally combine candidate algorithms extensively used in prediction problems -- to ODTRs. Following the "causal roadmap," we causally and statistically define the ODTR and provide an introduction to estimating it using the ODTR SuperLearner. Additionally, we highlight practical choices when implementing the algorithm, including choice of candidate algorithms, metalearners to combine the candidates, and risk functions to select the best combination of algorithms. Using simulations, we illustrate how estimating the ODTR using this SuperLearner approach can uncover treatment effect heterogeneity more effectively than traditional approaches based on fitting a parametric regression of the outcome on the treatment, covariates and treatment-covariate interactions. We investigate the implications of choices in implementing an ODTR SuperLearner at various sample sizes. Our results show the advantages of: (1) including a combination of both flexible machine learning algorithms and simple parametric estimators in the library of candidate algorithms; (2) using an ensemble metalearner to combine candidates rather than selecting only the best-performing candidate; (3) using the mean outcome under the rule as a risk function. Finally, we apply the ODTR SuperLearner to the "Interventions" study, an ongoing randomized controlled trial, to identify which justice-involved adults with mental illness benefit most from cognitive behavioral therapy (CBT) to reduce criminal re-offending.


title: Counting orientations of graphs with no strongly connected tournaments

abstract: Let $S_k(n)$ be the maximum number of orientations of an $n$-vertex graph $G$ in which no copy of $K_k$ is strongly connected. For all integers $n$, $k\geq 4$ where $n\geq 5$ or $k\geq 5$, we prove that $S_k(n) = 2^{t_{k-1}(n)}$, where $t_{k-1}(n)$ is the number of edges of the $n$-vertex $(k-1)$-partite Tur\'an graph $T_{k-1}(n)$, and that $T_{k-1}(n)$ is the only $n$-vertex graph with this number of orientations. Furthermore, $S_4(4) = 40$ and this maximality is achieved only by $K_4$.


title: Coalition Game Based Full-duplex Popular Content Distribution in mmWave Vehicular Networks

abstract: The millimeter wave (mmWave) communication has drawn intensive attention with abundant band resources. In this paper, we consider the popular content distribution (PCD) problem in the mmWave vehicular network. In order to offload the communication burden of base stations (BSs), vehicle-to-vehicle (V2V) communication is introduced into the PCD problem to transmit contents between on-board units (OBUs) and improve the transmission efficiency. We propose a full-duplex (FD) cooperative scheme based on coalition formation game, and the utility function is provided based on the maximization of the number of received contents. The contribution of each member in the coalition can be transferable to its individual profit. While maximizing the number of received contents in the fixed time, the cooperative scheme also ensures the individual profit of each OBU in the coalition. We evaluate the proposed scheme by extensive simulations in mmWave vehicular networks. Compared with other existing schemes, the proposed scheme has superior performances on the number of possessed contents and system fairness. Besides, the low complexity of the proposed algorithm is demonstrated by the switch operation number and CPU time.


title: Dynamically reversible and strong circular dichroism based on Babinet-invertible chiral metasurfaces

abstract: We propose a Babinet-invertible chiral metasurface for achieving dynamically reversible and strong circular dichroism (CD). The proposed metasurface is composed of VO$_2$-metal hybrid structure, and when VO$_2$ transits between the dielectric state and the metallic state, the metasurface unit cell switches between complementary structures that are designed according to the Babinet principle. This leads to a large and reversible CD tuning range between $\pm 0.5$ at 0.97~THz, which is larger than the literature. We attribute the CD effect to extrinsic chirality of the proposed metasurface. We envision that the Babinet-invertible chiral metasurface proposed here will advance the engineering of active and tunable chiro-optical devices and promote their applications.


title: Multi-Dimensional Parameter Study of Double Detonation Type Ia Supernovae Originating from Thin-Helium-Shell White Dwarfs

abstract: Despite the importance of Type Ia supernovae (SNe Ia) throughout astronomy, the precise progenitor systems and explosion mechanisms that drive SNe Ia are still unknown. An explosion scenario that has gained traction recently is the double detonation in which an accreted shell of He detonates and triggers a secondary detonation in the underlying white dwarf. Our research presents a number of high resolution, multi-dimensional, full star simulations of thin-He-shell, sub-Chandrasekhar-mass white dwarf progenitors that undergo a double detonation. This suite of thin-shell progenitors incorporates He shells that are thinner than those in previous multi-dimensional studies. We confirm the viability of the double detonation across a range of He shell parameter space as well as present bulk yields and ejecta profiles for each progenitor. The yields obtained are generally consistent with previous works and indicate the likelihood of producing observables that resemble SNe Ia. The dimensionality of our simulations allow us to examine features of the double detonation more closely, including the details of the off-center secondary ignition and asymmetric ejecta. We find considerable differences in the high-velocity extent of post-detonation products across different lines of sight. The data from this work will be used to generate predicted observables and may further support the viability of the double detonation scenario as a SNe Ia channel as well as show how properties of the progenitor or viewing angle may influence trends in observable characteristics.


title: A Pub-Sub Architecture to Promote Blockchain Interoperability

abstract: The maturing of blockchain technology leads to heterogeneity, where multiple solutions specialize in a particular use case. While the development of different blockchain networks shows great potential for blockchains, the isolated networks have led to data and asset silos, limiting the applications of this technology. Blockchain interoperability solutions are essential to enable distributed ledgers to reach their full potential. Such solutions allow blockchains to support asset and data transfer, resulting in the development of innovative applications. This paper proposes a novel blockchain interoperability solution for permissioned blockchains based on the publish/subscribe architecture. We implemented a prototype of this platform to show the feasibility of our design. We evaluate our solution by implementing examples of the different publisher and subscriber networks, such as Hyperledger Besu, which is an Ethereum client, and two different versions of Hyperledger Fabric. We present a performance analysis of the whole network that indicates its limits and bottlenecks. Finally, we discuss the extensibility and scalability of the platform in different scenarios. Our evaluation shows that our system can handle a throughput in the order of the hundreds of transactions per second.


title: Atomic Swaps between Bitcoin and Monero

abstract: Due to the evergrowing blockchain ecosystem, interoperability has become a matter of great importance. Atomic swaps allow connecting otherwise isolated blockchains while adhering to the core principles of censorship resistance and permissionlessnes. Up until recently, atomic swap protocols have mostly relied on complex script support, excluding certain types of blockchains. With advances in cryptography, it is now possible to build a bridge between almost any two blockchains. In this work, we give an explanation of one such protocol which applies adaptor signatures on Bitcoin to procure atomic swaps between Monero and Bitcoin. We dive into the cryptographic details, discuss its limitations and give an outlook on our current work where we use adaptor signatures on the Monero signature scheme.


title: Performance and Application of Estimators for the Value of an Optimal Dynamic Treatment Rule

abstract: Given an (optimal) dynamic treatment rule, it may be of interest to evaluate that rule -- that is, to ask the causal question: what is the expected outcome had every subject received treatment according to that rule? In this paper, we study the performance of estimators that approximate the true value of: 1) an $a$ $priori$ known dynamic treatment rule 2) the true, unknown optimal dynamic treatment rule (ODTR); 3) an estimated ODTR, a so-called "data-adaptive parameter," whose true value depends on the sample. Using simulations of point-treatment data, we specifically investigate: 1) the impact of increasingly data-adaptive estimation of nuisance parameters and/or of the ODTR on performance; 2) the potential for improved efficiency and bias reduction through the use of semiparametric efficient estimators; and, 3) the importance of sample splitting based on CV-TMLE for accurate inference. In the simulations considered, there was very little cost and many benefits to using the cross-validated targeted maximum likelihood estimator (CV-TMLE) to estimate the value of the true and estimated ODTR; importantly, and in contrast to non cross-validated estimators, the performance of CV-TMLE was maintained even when highly data-adaptive algorithms were used to estimate both nuisance parameters and the ODTR. In addition, we apply these estimators for the value of the rule to the "Interventions" Study, an ongoing randomized controlled trial, to identify whether assigning cognitive behavioral therapy (CBT) to criminal justice-involved adults with mental illness using an ODTR significantly reduces the probability of recidivism, compared to assigning CBT in a non-individualized way.


title: sGrapp: Butterfly Approximation in Streaming Graphs

abstract: We study the fundamental problem of butterfly (i.e. (2,2)-bicliques) counting in bipartite streaming graphs. Similar to triangles in unipartite graphs, enumerating butterflies is crucial in understanding the structure of bipartite graphs. This benefits many applications where studying the cohesion in a graph shaped data is of particular interest. Examples include investigating the structure of computational graphs or input graphs to the algorithms, as well as dynamic phenomena and analytic tasks over complex real graphs. Butterfly counting is computationally expensive, and known techniques do not scale to large graphs; the problem is even harder in streaming graphs. In this paper, following a data-driven methodology, we first conduct an empirical analysis to uncover temporal organizing principles of butterflies in real streaming graphs and then we introduce an approximate adaptive window-based algorithm, sGrapp, for counting butterflies as well as its optimized version sGrapp-x. sGrapp is designed to operate efficiently and effectively over any graph stream with any temporal behavior. Experimental studies of sGrapp and sGrapp-x show superior performance in terms of both accuracy and efficiency.


title: Personalization and Recommendation Technologies for MaaS

abstract: Over the last few years, MaaS has been extensively studied and evolved into offering a multitude of mobility services that continuously increase, from alternative car or bike-sharing modes to autonomous vehicles, that aspire to become a part of this novel ecosystem. MaaS provides end-users with multimodal, integrated, and digital mobility solutions, including a multitude of different choices able to cover users specific needs in a personalized manner. This practically leads to a range of novel MaaS products, that may have complex structures and the challenge of matching them to user preferences and needs so that suitable products can be provided to end-users. Moreover, in the everyday use of MaaS, travelers require support to identify routes to reach their destination that adhere to their personal preferences and are aligned to the MaaS product they have purchased. This paper tackles these two user-centric challenges by exploiting state-of-the-art techniques from the field of Personalization and Recommendation systems and integrating them in MaaS platforms and route planning applications.


title: Community Detection in the Stochastic Block Model by Mixed Integer Programming

abstract: The Degree-Corrected Stochastic Block Model (DCSBM) is a popular model to generate random graphs with community structure given an expected degree sequence. The standard approach of community detection based on the DCSBM is to search for the model parameters that are the most likely to have produced the observed network data through maximum likelihood estimation (MLE). Current techniques for the MLE problem are heuristics, and therefore do not guarantee convergence to the optimum. We present mathematical programming formulations and exact solution methods that can provably find the model parameters and community assignments of maximum likelihood given an observed graph. We compare these exact methods with classical heuristic algorithms based on expectation-maximization (EM). The solutions given by exact methods give us a principled way of measuring the experimental performance of classical heuristics and comparing different variations thereof.


title: The geometrical origin of the Doppler factor in the Lienard-Wiechert potentials

abstract: We present an in depth analysis and a new derivation of the Doppler factor in the Lienard-Wiechert potentials, based on geometrical considerations in Minkowski space. We argue that, contrary to a common assumption, the methods used for deriving the Doppler factor in the case of an electrically charged extended particle are not applicable in the case of a point particle. A geometrical interpretation of the Doppler factor is nonetheless found in the later case, based on the electromagnetic interaction model of Fokker. In this model the interaction takes place between infinitesimal worldline segments with end points connected by light signals, just like the points where light pulses are emitted and detected in general. This analogy reveals that the relativistic Doppler effect is the missing link between the classical Doppler effect and the Doppler factor in the Lienard-Wiechert potentials.