title: Random walks on mated-CRT planar maps and Liouville Brownian motion

abstract: We prove a scaling limit result for random walk on certain random planar maps with its natural time parametrization. In particular, we show that for $\gamma \in (0,2)$, the random walk on the mated-CRT map with parameter $\gamma$ converges to $\gamma$-Liouville Brownian motion, the natural quantum time parametrization of Brownian motion on a $\gamma$-Liouville quantum gravity (LQG) surface. Our result applies if the mated-CRT map is embedded into the plane via the embedding which comes from SLE / LQG theory or via the Tutte embedding (a.k.a. the harmonic or barycentric embedding). In both cases, the convergence is with respect to the local uniform topology on curves and it holds in the quenched sense, i.e., the conditional law of the walk given the map converges. Previous work by Gwynne, Miller, and Sheffield (2017) showed that the random walk on the mated-CRT map converges to Brownian motion modulo time parametrization. This is the first work to show the convergence of the parametrized walk. As an intermediate result of independent interest, we derive an axiomatic characterisation of Liouville Brownian motion, for which the notion of Revuz measure of a Markov process plays a crucial role.


title: Convergence of a Godunov scheme to an Audusse-Perthame adapted entropy solution for conservation laws with BV spatial flux

abstract: In this article we consider the initial value problem for a scalar conservation law in one space dimension with a spatially discontinuous flux. There may be infinitely many flux discontinuities, and the set of discontinuities may have accumulation points. Thus the existence of traces cannot be assumed. In [6] Audusse and Perthame proved a uniqueness result that does not require the existence of traces, using adapted entropies. We generalize the Godunov-type scheme of Adimurthi, Jaffre and Gowda [2] for this problem with the following assumptions on the flux function, (i) the flux is BV in the spatial variable and (ii) the critical point of the flux is BV as a function of the space variable. We prove that the Godunov approximations converge to an adapted entropy solution, thus providing an existence result, and extending the convergence result of Adimurthi, Jaffre and Gowda.


title: Catalogue with visual morphological classification of 32,616 radio galaxies with optical hosts

abstract: We present the catalogue of Radio sources associated with Optical Galaxies and having Unresolved or Extended morphologies I (ROGUE I). It was generated by cross-matching galaxies from the Sloan Digital Sky Survey Data Release 7 (SDSS DR 7) as well as radio sources from the First Images of Radio Sky at Twenty Centimetre (FIRST) and the National Radio Astronomical Observatory VLA Sky Survey (NVSS) catalogues. We created the largest handmade catalogue of visually classified radio objects and associated with them optical host galaxies, containing 32,616 galaxies with a FIRST core within 3 arcsec of the optical position. All listed objects possess the good quality SDSS DR 7 spectra with the signal-to-noise ratio $>$10 and spectroscopic redshifts up to $z=0.6$. The radio morphology classification was performed by a visual examination of the FIRST and the NVSS contour maps overlaid on a DSS image, while an optical morphology classification was based on the 120 arcsec snapshot images from SDSS DR 7. The majority of radio galaxies in ROGUE I, i.e. $\sim$ 93%, are unresolved (compact or elongated), while the rest of them exhibit extended morphologies, such as Fanaroff-Riley (FR) type I, II, and hybrid, wide-angle tail, narrow-angle tail, head-tail sources, and sources with intermittent or reoriented jet activity, i.e. double-double, X-shaped, and Z-shaped. Most of FR IIs have low radio luminosities, comparable to the luminosities of FR Is. Moreover, due to visual check of all radio maps and optical images, we were able to discover or reclassify a number of radio objects as giant, double-double, X-shaped, and Z-shaped radio galaxies. The presented sample can serve as a database for training automatic methods of identification and classification of optical and radio galaxies.


title: Monte Carlo integration of non-differentiable functions on $[0,1]^\iota$, $\iota=1,\dots,d$, using a single determinantal point pattern defined on $[0,1]^d$

abstract: This paper concerns the use of a particular class of determinantal point processes (DPP), a class of repulsive spatial point processes, for Monte Carlo integration. Let $d\ge 1$, $I\subseteq \overline d=\{1,\dots,d\}$ with $\iota=|I|$. Using a single set of $N$ quadrature points $\{u_1,\dots,u_N\}$ defined, once for all, in dimension $d$ from the realization of the DPP model, we investigate "minimal" assumptions on the integrand in order to obtain unbiased Monte Carlo estimates of $\mu(f_I)=\int_{[0,1]^\iota} f_I(u) \mathrm{d} u$ for any known $\iota$-dimensional integrable function on $[0,1]^\iota$. In particular, we show that the resulting estimator has variance with order $N^{-1-(2s\wedge 1)/d}$ when the integrand belongs to some Sobolev space with regularity $s > 0$. When $s>1/2$ (which includes a large class of non-differentiable functions), the variance is asymptotically explicit and the estimator is shown to satisfy a Central Limit Theorem.

categories: stat.CO cs.NA math.CA math.NA math.ST stat.TH


title: Fast Alternating Projections on Manifolds Based on Tangent Spaces

abstract: In this paper, we study alternating projections on nontangential manifolds based on the tangent spaces. The main motivation is that the projection of a point onto a manifold can be computational expensive. We propose to use the tangent space of the point in the manifold to approximate the projection onto the manifold in order to reduce the computational cost. We show that the sequence generated by alternating projections on two nontangential manifolds based on tangent spaces, converges linearly to a point in the intersection of the two manifolds where the convergent point is close to the optimal solution. Numerical examples for nonnegative low rank matrix approximation and low rank image quaternion matrix (color image) approximation, are given to demonstrate that the performance of the proposed method is better than that of the classical alternating projection method in terms of computational time.


title: DYSAN: Dynamically sanitizing motion sensor data against sensitive inferences through adversarial networks

abstract: With the widespread adoption of the quantified self movement, an increasing number of users rely on mobile applications to monitor their physical activity through their smartphones. Granting to applications a direct access to sensor data expose users to privacy risks. Indeed, usually these motion sensor data are transmitted to analytics applications hosted on the cloud leveraging machine learning models to provide feedback on their health to users. However, nothing prevents the service provider to infer private and sensitive information about a user such as health or demographic attributes.In this paper, we present DySan, a privacy-preserving framework to sanitize motion sensor data against unwanted sensitive inferences (i.e., improving privacy) while limiting the loss of accuracy on the physical activity monitoring (i.e., maintaining data utility). To ensure a good trade-off between utility and privacy, DySan leverages on the framework of Generative Adversarial Network (GAN) to sanitize the sensor data. More precisely, by learning in a competitive manner several networks, DySan is able to build models that sanitize motion data against inferences on a specified sensitive attribute (e.g., gender) while maintaining a high accuracy on activity recognition. In addition, DySan dynamically selects the sanitizing model which maximize the privacy according to the incoming data. Experiments conducted on real datasets demonstrate that DySan can drasticallylimit the gender inference to 47% while only reducing the accuracy of activity recognition by 3%.


title: The possible explanation for the mysterious hot solar corona

abstract: The coronal heating has remained a puzzle for decades since its observation in 1940. Several space missions have been planned to resolve this conundrum in the past and many more intend to target this issue in the future. The unfolding of this issue will not only advance the fundamentals of astrophysical science but also promise an improvement in space weather prediction. Its origin has been debated without complete convergence; the acoustic waves, magneto-hydrodynamic waves, and micro/nano-flares are the strongest candidates so far to explain the mystery of coronal heating. However, none of these processes significantly justifies the million-degree temperature of the solar corona and the problem remains unsolved to date. Here, we propose a new physical mechanism to explain the observed heating of the solar corona. The statistical energy created during the interaction of the spin magnetic moment of the plasma particles with the turbulent and continuously evolving coronal magnetic field could substantiate the observed million-degree coronal temperature.


title: Error bounds of regularized gap functions for polynomial variational inequalities

abstract: This paper is devoted to present new error bounds of regularized gap functions for polynomial variational inequalities with exponents explicitly determined by the dimension of the underlying space and the number/degree of the involved polynomials. The developed techniques are largely based on variational analysis and semialgebraic geometry, which allow us to establish a nonsmooth extension of the seminal \L ojasiewicz's gradient inequality to regularized gap functions with explicitly calculated exponents.


title: A Physically Universal Turing Machine

abstract: We construct a two-dimensional Turing machine that is physically universal in both the moving tape and moving head model. In particular, it is mixing of all finite orders in both models. We also provide a variant that is physically universal in the moving tape model, but not in the moving head model.


title: Critical exponent for the wave equation with a time-dependent scale invariant damping and a cubic convolution

abstract: In the present paper, we study the Cauchy problem for the wave equation with a time-dependent scale invariant damping $\frac{2}{1+t}\partial_t v$ and a cubic convolution $(|x|^{-\gamma}*v^2)v$ with $\gamma\in \left(-\frac{1}{2},3\right)$ in three spatial dimension for initial data $\left(v(x,0),\partial_tv(x,0)\right)\in C^2(\mathbb{R}^3)\times C^1(\mathbb{R}^3)$ with a compact support, where $v=v(x,t)$ is an unknown function to the problem on $\mathbb{R}^3\times[0,T)$. Here $T$ denotes a maximal existence time of $v$. The first aim of the present paper is to prove unique global existence of the solution to the problem and asymptotic behavior of the solution in the supercritical case $\gamma\in (0,3)$, and show a lower estimate of the lifespan in the critical or subcritical case $\gamma\in \left(-\frac{1}{2},0\right]$. The essential part for their proofs is to derive a weaker estimate under the weaker condition than the case without damping and to recover the weakness by the effect of the dissipative term. The second aim of the present paper is to prove a small data blow-up and the almost sharp upper estimate of the lifespan for positive data with a compact support in the subcritical case $\gamma\in \left(-\frac{1}{2},0\right)$. The essential part for the proof is to refine the argument for the proof of Theorem 6.1 in \cite{H20} to obtain the upper estimate of the lifespan. Our two results determine that a critical exponent $\gamma_c$ which divides global existence and blow-up for small solutions is $0$, namely $\gamma_c=0$. As the result, we can see that the critical exponent shift from $2$ to $0$ due to the effect of the scale invariant damping term.


title: Latent Model Extreme Value Index Estimation

abstract: We propose a novel strategy for multivariate extreme value index estimation. In applications such as finance, volatility and risk present in the components of a multivariate time series are often driven by the same underlying factors, such as the subprime crisis in the US. To estimate the latent risk, we apply a two-stage procedure. First, a set of independent latent series is estimated using a method of latent variable analysis. Then, univariate risk measures are estimated individually for the latent series to assess their contribution to the overall risk. As our main theoretical contribution, we derive conditions under which the effect of the first step to the asymptotic behavior of the risk estimators is negligible. Simulations demonstrate the theory under both i.i.d. and dependent data, and an application into financial data illustrates the usefulness of the method in extracting joint sources of risk in practice.


title: Varying Newton constant and black hole to white hole quantum tunneling

abstract: The thermodynamics of black holes is discussed for the case, when the Newton constant $G$ is not a constant, but is the thermodynamic variable. This gives for the first law of the Schwarzschild black hole thermodynamics: $dS_\text{BH}= -AdK + \frac{dM}{T_\text{BH}}$, where the gravitational coupling $K=1/4G$, $M$ is the black hole mass, $A$ is the area of horizon, and $T_\text{BH}$ is Hawking temperature. From this first law it follows that the dimensionless quantity $M^2/K$ is the adiabatic invariant, which in principle can be quantized if to follow the Bekenstein conjecture. From the Euclidean action for the black hole it follows that $K$ and $A$ serve as dynamically conjugate variables. This allows us to calculate the quantum tunneling from the black hole to the white hole, and determine the temperature and entropy of the white hole.


title: Stochastic Event-based Sensor Schedules for Remote State Estimation in Cognitive Radio Sensor Networks

abstract: We consider the problem of communication allocation for remote state estimation in a cognitive radio sensor network~(CRSN). A sensor collects measurements of a physical plant, and transmits the data to a remote estimator as a secondary user (SU) in the shared network. The existence of the primal users (PUs) brings exogenous uncertainties into the transmission scheduling process, and how to design an event-based scheduling scheme considering these uncertainties has not been addressed in the literature. In this work, we start from the formulation of a discrete-time remote estimation process in the CRSN, and then analyze the hidden information contained in the absence of data transmission. In order to achieve a better tradeoff between estimation performance and communication consumption, we propose both open-loop and closed-loop schedules using the hidden information under a Bayesian setting. The open-loop schedule does not rely on any feedback signal but only works for stable plants. For unstable plants, a closed-loop schedule is designed based on feedback signals. The parameter design problems in both schedules are efficiently solved by convex programming. Numerical simulations are included to illustrate the theoretical results.


title: Neural Contours: Learning to Draw Lines from 3D Shapes

abstract: This paper introduces a method for learning to generate line drawings from 3D models. Our architecture incorporates a differentiable module operating on geometric features of the 3D model, and an image-based module operating on view-based shape representations. At test time, geometric and view-based reasoning are combined with the help of a neural module to create a line drawing. The model is trained on a large number of crowdsourced comparisons of line drawings. Experiments demonstrate that our method achieves significant improvements in line drawing over the state-of-the-art when evaluated on standard benchmarks, resulting in drawings that are comparable to those produced by experienced human artists.


title: Two-path interference for enantiomer-selective state transfer of chiral molecules

abstract: With a microwave-regime cyclic three-state configuration, an enantiomer-selective state transfer~(ESST) is carried out through the two-path interference between a direct one-photon coupling and an effective two-photon coupling. The $\pi$-phase difference in the one-photon process between two enantiomers makes the interference constructive for one enantiomer but destructive for the other. Therefore only one enantiomer is excited into a higher rotational state while the other remains in the ground state. The scheme is of flexibility in the pulse waveforms and the time order of two paths. We simulate the scheme in a sample of cyclohexylmethanol~(C$_7$H$_{14}$O) molecules. Simulative results show the robust and high-fidelity ESST can be obtained when experimental concerns are considered. Finally, we propose to employ the finished ESST in implementing enantio-separation and determining enantiomeric excess.


title: Heptagon relation in a direct sum

abstract: An ansatz is proposed for heptagon relation, that is, algebraic imitation of five-dimensional Pachner move 4--3. Our relation is realized in terms of matrices acting in a direct sum of one-dimensional linear spaces corresponding to 4-faces.


title: Efficient Bayesian Inference of General Gaussian Models on Large Phylogenetic Trees

abstract: Phylogenetic comparative methods correct for shared evolutionary history among a set of non-independent organisms by modeling sample traits as arising from a diffusion process along on the branches of a possibly unknown history. To incorporate such uncertainty, we present a scalable Bayesian inference framework under a general Gaussian trait evolution model that exploits Hamiltonian Monte Carlo (HMC). HMC enables efficient sampling of the constrained model parameters and takes advantage of the tree structure for fast likelihood and gradient computations, yielding algorithmic complexity linear in the number of observations. This approach encompasses a wide family of stochastic processes, including the general Ornstein-Uhlenbeck (OU) process, with possible missing data and measurement errors. We implement inference tools for a biologically relevant subset of all these models into the BEAST phylogenetic software package and develop model comparison through marginal likelihood estimation. We apply our approach to study the morphological evolution in the superfamilly of Musteloidea (including weasels and allies) as well as the heritability of HIV virulence. This second problem furnishes a new measure of evolutionary heritability that demonstrates its utility through a targeted simulation study.


title: Small Weight Code Words of Projective Geometric Codes

abstract: We investigate small weight code words of the $p$-ary linear code $\mathcal C_{j,k}(n,q)$ generated by the incidence matrix of $k$-spaces and $j$-spaces of PG$(n,q)$ and its dual, with $q$ a prime power and $0 \leq j < k < n$. Firstly, we prove that all code words of $\mathcal C_{j,k}(n,q)$ up to weight $\left(3 - \mathcal{O}\left(\frac 1 q \right) \right) \genfrac{[}{]}{0pt}{}{k+1}{j+1}_q$ are linear combinations of at most two $k$-spaces (i.e. two rows of the incidence matrix). As for the dual code $\mathcal C_{j,k}(n,q)^\perp$, we manage to reduce both problems of determining its minimum weight (1) and characterising its minimum weight code words (2) to the case $\mathcal C_{0,1}(n,q)^\perp$. This implies the solution to both problem (1) and (2) if $q$ is prime and the solution to problem (1) if $q$ is even.


title: Quantum critical thermal transport in the unitary Fermi gas

abstract: Strongly correlated systems are often associated with an underlying quantum critical point which governs their behavior in the finite temperature phase diagram. Their thermodynamical and transport properties arise from critical fluctuations and follow universal scaling laws. Here, we develop a microscopic theory of thermal transport in the quantum critical regime expressed in terms of a thermal sum rule and an effective scattering time. We explicitly compute the characteristic scaling functions in a quantum critical model system, the unitary Fermi gas. Moreover, we derive an exact thermal sum rule for heat and energy currents and evaluate it numerically using the nonperturbative Luttinger-Ward approach. For the thermal scattering times we find a simple quantum critical scaling form. Together, the sum rule and the scattering time determine the heat conductivity, thermal diffusivity, Prandtl number and sound diffusivity from high temperatures down into the quantum critical regime. The results provide a quantitative description of recent sound attenuation measurements in ultracold Fermi gases.


title: Diffusion-based Deep Active Learning

abstract: The remarkable performance of deep neural networks depends on the availability of massive labeled data. To alleviate the load of data annotation, active deep learning aims to select a minimal set of training points to be labelled which yields maximal model accuracy. Most existing approaches implement either an `exploration'-type selection criterion, which aims at exploring the joint distribution of data and labels, or a `refinement'-type criterion which aims at localizing the detected decision boundaries. We propose a versatile and efficient criterion that automatically switches from exploration to refinement when the distribution has been sufficiently mapped. Our criterion relies on a process of diffusing the existing label information over a graph constructed from the hidden representation of the data set as provided by the neural network. This graph representation captures the intrinsic geometry of the approximated labeling function. The diffusion-based criterion is shown to be advantageous as it outperforms existing criteria for deep active learning.


title: Entropy as a measure of attractiveness and socioeconomic complexity in Rio de Janeiro metropolitan area

abstract: Defining and measuring spatial inequalities across the urban environment remains a complex and elusive task that has been facilitated by the increasing availability of large geolocated databases. In this study, we rely on a mobile phone dataset and an entropy-based metric to measure the attractiveness of a location in the Rio de Janeiro Metropolitan Area (Brazil) as the diversity of visitors' location of residence. The results show that the attractiveness of a given location measured by entropy is an important descriptor of the socioeconomic status of the location, and can thus be used as a proxy for complex socioeconomic indicators.


title: Insights into the "cross-world" independence assumption of causal mediation analysis

abstract: Causal mediation analysis is a useful tool for epidemiological research, but it has been criticized for relying on a "cross-world" independence assumption that is empirically difficult to verify and problematic to justify based on background knowledge. In the present article we aim to assist the applied researcher in understanding this assumption. Synthesizing what is known about the cross-world independence assumption, we discuss the relationship between assumptions for causal mediation analyses, causal models, and non-parametric identification of natural direct and indirect effects. In particular we give a practical example of an applied setting where the cross-world independence assumption is violated even without any post-treatment confounding. Further, we review possible alternatives to the cross-world independence assumption, including the use of computation of bounds that avoid the assumption altogether. Finally, we carry out a numerical study in which the cross-world independence assumption is violated to assess the ensuing bias in estimating natural direct and indirect effects. We conclude with recommendations for carrying out causal mediation analyses.


title: Rate of Convergence for Distributed Optimization with Uncertain Communications

abstract: We consider the distributed optimization problem for the sum of convex functions where the underlying communications network connecting agents at each time is drawn at random from a collection of directed graphs. Building on an earlier work [15], where a modified version of the subgradient-push algorithm is shown to be almost surely convergent to an optimizer on sequences of random directed graphs, we find an upper bound of the order of $\sim O(\frac{1}{\sqrt{t}})$ on the convergence rate of our proposed algorithm, establishing the first convergence bound in such random settings.


title: Spin-wave gap collapse in Rh-doped Sr2IrO4

abstract: We use resonant inelastic x-ray scattering (RIXS) at the Ir L3 edge to study the effect of hole doping upon the Jeff=1/2 Mott-insulating state in Sr2IrO4, via Rh replacement of the Ir site. The spin-wave gap, associated with XY-type spin-exchange anisotropy, collapses with increasing Rh content, prior to the suppression of the Mott-insulating state and in contrast to electron doping via La substitution of the Sr site. At the same time, despite heavy damping, the d-d excitation spectra retain their overall amplitude and dispersion character. A careful study of the spin-wave spectrum reveals that deviations from the J1-J2-J3 Heisenberg used to model the pristine system disappear at intermediate doping levels. These findings are interpreted in terms of a modulation of Ir-Ir correlations due to the influence of Rh impurities upon nearby Ir wave functions, even as the single-band Jeff=1/2 model remains valid up to full carrier delocalization. They underline the importance of the transition metal site symmetry when doping pseudospin systems such as Sr2IrO4.


title: Purely inseparable coverings of rational double points in positive characteristic

abstract: We classify purely inseparable morphisms of degree $p$ between rational double points (RDPs) in characteristic $p > 0$. Using such morphisms, we refine a result of Artin that any RDP admits a finite smooth covering.


title: Berezin-Toeplitz quantization and the least unsharpness principle

abstract: We show that compatible almost-complex structures on symplectic manifolds correspond to optimal quantizations.


title: Averaging 2d stochastic wave equation

abstract: We consider a 2D stochastic wave equation driven by a Gaussian noise, which is temporally white and spatially colored described by the Riesz kernel. Our first main result is the functional central limit theorem for the spatial average of the solution. And we also establish a quantitative central limit theorem for the marginal and the rate of convergence is described by the total-variation distance. A fundamental ingredient in our proofs is the pointwise $L^p$-estimate of Malliavin derivative, which is of independent interest.


title: Distributed Set-Based Observers Using Diffusion Strategies

abstract: We propose two distributed set-based observers using strip-based and set-propagation approaches for linear discrete-time dynamical systems with bounded modeling and measurement uncertainties. Both algorithms utilize a set-based diffusion step, which decreases the estimation errors and the size of estimated sets, and can be seen as a lightweight approach to achieve partial consensus between the distributed estimated sets. Every node shares its measurement with its neighbor in the measurement update step. In the diffusion step, the neighbors intersect their estimated sets using our novel lightweight zonotope intersection technique. A localization example demonstrates the applicability of our algorithms.


title: Distributed Discontinuous Coupling for Convergence in Networks of Heterogeneous Nonlinear Systems

abstract: Synchronization is a crucial phenomenon in many natural and artificial complex network systems. Applications include neuronal networks, formation control and coordination in robotics, and frequency synchronization in electrical power grids. In this paper, we propose the use of a distributed discontinuous coupling protocol to achieve convergence and synchronization in networks of non-identical nonlinear dynamical systems. We show that the synchronous dynamics is a solution to the average of the nodes' vector fields, and derive analytical estimates of the critical coupling gains required to achieve convergence. Numerical simulations are used to illustrate and validate the theoretical results.


title: Hydrodynamic-colloidal interactions of an oil droplet and a membrane surface

abstract: Understanding the influence of hydrodynamic stresses and colloidal interactions between droplets and membranes will inform better engineered membrane surfaces and process conditions. Using a long-wave hydrodynamic description that incorporates wettability and electrostatic repulsion, we numerically investigate the influence of droplet shapes and membrane characteristics on the overall repulsion-attraction towards the membrane surface. A stability phase diagram classifies stable-unstable drop shapes and investigate whether there is reversible deposition of droplets onto the membrane.


title: Weakly Supervised 3D Human Pose and Shape Reconstruction with Normalizing Flows

abstract: Monocular 3D human pose and shape estimation is challenging due to the many degrees of freedom of the human body and thedifficulty to acquire training data for large-scale supervised learning in complex visual scenes. In this paper we present practical semi-supervised and self-supervised models that support training and good generalization in real-world images and video. Our formulation is based on kinematic latent normalizing flow representations and dynamics, as well as differentiable, semantic body part alignment loss functions that support self-supervised learning. In extensive experiments using 3D motion capture datasets like CMU, Human3.6M, 3DPW, or AMASS, as well as image repositories like COCO, we show that the proposed methods outperform the state of the art, supporting the practical construction of an accurate family of models based on large-scale training with diverse and incompletely labeled image and video data.


title: A modular phantom and software to characterize 3D geometric distortion in MRI

abstract: MRI offers outstanding soft tissue contrast that may reduce uncertainties in target and organ-at-risk delineation and enable online adaptive image-guided treatment. Spatial distortions resulting from non-linearities in the gradient fields and non-uniformity in the main magnetic field must be accounted for across the imaging field-of-view to prevent systematic errors during treatment delivery. This work presents a modular phantom and software application to characterize geometric distortion (GD) within the large field-of-view MRI images required for radiation therapy simulation. The modular phantom is assembled from a series of rectangular foam blocks containing high-contrast fiducial markers in a known configuration. The modular phantom design facilitates transportation of the phantom between different MR scanners and MR-guided linear accelerators and allows the phantom to be adapted to fit different sized bores or coils. The phantom was evaluated using a 1.5T MR-guided linear accelerator (MR-Linac) and 1.5T and 3.0T diagnostic scanners. Performance was assessed by varying acquisition parameters to induce image distortions in a known manner. Imaging was performed using T1 and T2 weighted pulse sequences with 2D and 3D distortion correction algorithms and the receiver bandwidth (BW) varied as 250-815 Hz/pixel. Phantom set-up reproducibility was evaluated across independent set-ups. The software was validated by comparison with a non-modular phantom. Average geometric distortion was 0.94+/-0.58 mm for the MR-Linac, 0.90+/-0.53 mm for the 1.5 T scanner, and 1.15+/-0.62 mm for the 3.0T scanner, for a 400 mm diameter volume-of-interest. GD increased, as expected, with decreasing BW, and with the 2D versus 3D correction algorithm. Differences in GD attributed to phantom set-up were 0.13 mm or less. Differences in GD for the two software applications were less than 0.07 mm.


title: Timing stability of black widow pulsars

abstract: We study the timing stability of three black widow pulsars, both in terms of their long-term spin evolution and their shorter-term orbital stability. The erratic timing behaviour and radio eclipses of the first two black widow pulsar systems discovered (PSRs B1957+20 and J2051$-$0827) was assumed to be representative for this class of pulsars. With several new black widow systems added to this population in the last decade, there are now several systems known that do not show these typical orbital variations or radio eclipses. We present timing solutions using 7$-$8 yrs of observations from four of the European Pulsar Timing Array telescopes for PSRs\,J0023+0923, J2214+3000 and J2234+0944, and confirm that two of these systems do not show any significant orbital variability over our observing time span, both in terms of secular or orbital parameters. The third pulsar PSR J0023+0923 shows orbital variability and we discuss the implications for the timing solution. Our results from the long-term timing of these pulsars provide several new or improved parameters compared to earlier works. We discuss our results regarding the stability of these pulsars, and the stability of the class of black widow pulsars in general, in the context of the binary parameters, and discuss the potential of the Roche-lobe filling factor of the companion star being an indicator for stability of these systems.


title: Effects of MiFID II on stock price formation

abstract: This paper examines effects of MiFID II on European stock markets. We study the effects of the new tick size regime, both intraday and in the closing auction. An increase (decrease) in tick size is associated with a decrease (increase) in intraday liquidity, but a more (less) stable market. In the closing auction an increase in tick size has a positive effect on liquidity. Moreover, we report a positive relationship between tick size and transacted volume, in particular in the closing auction. Finally, closing auction volumes increased heavily since MiFID II and price formation in closing auctions became more efficient.


title: Fairway: A Way to Build Fair ML Software

abstract: Machine learning software is increasingly being used to make decisions that affect people's lives. But sometimes, the core part of this software (the learned model), behaves in a biased manner that gives undue advantages to a specific group of people (where those groups are determined by sex, race, etc.). This "algorithmic discrimination" in the AI software systems has become a matter of serious concern in the machine learning and software engineering community. There have been works done to find "algorithmic bias" or "ethical bias" in the software system. Once the bias is detected in the AI software system, the mitigation of bias is extremely important. In this work, we a)explain how ground-truth bias in training data affects machine learning model fairness and how to find that bias in AI software,b)propose a methodFairwaywhich combines pre-processing and in-processing approach to remove ethical bias from training data and trained model. Our results show that we can find bias and mitigate bias in a learned model, without much damaging the predictive performance of that model. We propose that (1) test-ing for bias and (2) bias mitigation should be a routine part of the machine learning software development life cycle. Fairway offers much support for these two purposes.


title: Exciton localization in tubular molecular aggregates: size effects and optical response

abstract: We study the exciton localization and resulting optical response for disordered tubular aggregates of optically active molecules. It has been shown previously that such tubular structures allow for excitons delocalized over more than a thousand molecules, owing to the combined effects of long-range dipole-dipole interactions and the higher-dimensional (not truly one-dimensional) nature of the aggregate. Such large delocalization sizes prompt the question to what extent in experimental systems the delocalization may still be determined by the aggregate size (diameter and length) and how this affects the aggregate's optical response and dynamics. We perform a systematic study of the size effects on the localization properties, using numerical simulations of the exciton states in a cylindrical model structure inspired by the previously derived geometry of a cylindrical aggregate of cyanine dye molecules (C8S3). To characterize the exciton localization, we calculate the participation ratio and the autocorrelation function of the exciton wave function. Also, we calculate the density of states and absorption spectrum. We find strong effects of the tube's radius on the localization and optical properties in the range of parameters relevant to experiment. In addition, surprisingly, we find that even for tubes as long as 750 nm, the localization size is limited by the tube's length for disorder values that are relevant to experimental circumstances, while observable effects of the tube's length in the absorption spectrum still occur for tube lengths up to about 150 nm. The latter may explain changes in the optical spectra observed during the aging process of bromine-substituted C8S3 aggregates. For weak disorder, the exciton wave functions exhibit a scattered, fractal-like nature, similar to the quasi-particles in two-dimensional disordered systems.


title: Dispersive estimates for 2D-wave equations with critical potentials

abstract: We study the 2D-wave equation with a scaling-critical electromagnetic potential. This problem is doubly critical, because of the scaling invariance of the model and the singularities of the potentials, which are not locally integrable. In particular, the diamagnetic phenomenon allows to consider negative electric potential which can be singular in the same fashion as the inverse-square potential. We prove sharp time-decay estimates in the purely magnetic case, and Strichartz estimates for the complete model, involving a critical electromagnetic field.


title: Projective toric codes

abstract: Any integral convex polytope $P$ in $\mathbb{R}^N$ provides a $N$-dimensional toric variety $X_P$ and an ample divisor $D_P$ on this variety. This paper gives an explicit construction of the algebraic geometric error-correcting code on $X_P$ , obtained by evaluating global section of $\mathcal{L}(D_P)$ on every rational point of $X_P$. This work presents an extension of toric codes analogous to the one of Reed-Muller codes into projective ones, by evaluating on the whole variety instead of considering only points with non-zero coordinates. The dimension of the code is given in terms of the number of integral points in the polytope $P$ and an algorithmic technique to get a lowerbound on the minimum distance is described.


title: Atomistic perspective of long lifetimes of small skyrmions at room temperature

abstract: The current development to employ magnetic skyrmions in novel spintronic device designs has led to a demand for room temperature-stable skyrmions of ever smaller size. We present extensive studies on skyrmion stability in atomistic magnetic systems in two- and three-dimensional geometries. We show that for materials described by the same micromagnetic parameters, the variation of the atomistic exchange between different neighbors, the stacking order, and the number of layers of the atomic lattice can significantly influence the rate of the thermally activated decay of a skyrmion. These factors alone are important considerations, but it is shown that their combination can open up novel avenues of materials design in the search for sub-10nm skyrmions, as their lifetime can be extended by several orders of magnitude.


title: Understanding the perception of COVID-19 policies by mining a multilanguage Twitter dataset

abstract: The objective of this work is to explore popular discourse about the COVID-19 pandemic and policies implemented to manage it. Using Natural Language Processing, Text Mining, and Network Analysis to analyze corpus of tweets that relate to the COVID-19 pandemic, we identify common responses to the pandemic and how these responses differ across time. Moreover, insights as to how information and misinformation were transmitted via Twitter, starting at the early stages of this pandemic, are presented. Finally, this work introduces a dataset of tweets collected from all over the world, in multiple languages, dating back to January 22nd, when the total cases of reported COVID-19 were below 600 worldwide. The insights presented in this work could help inform decision makers in the face of future pandemics, and the dataset introduced can be used to acquire valuable knowledge to help mitigate the COVID-19 pandemic.


title: Bayesian Models Applied to Cyber Security Anomaly Detection Problems

abstract: Cyber security is an important concern for all individuals, organisations and governments globally. Cyber attacks have become more sophisticated, frequent and dangerous than ever, and traditional anomaly detection methods have been proved to be less effective when dealing with these new classes of cyber threats. In order to address this, both classical and Bayesian models offer a valid and innovative alternative to the traditional signature-based methods, motivating the increasing interest in statistical research that it has been observed in recent years. In this review we provide a description of some typical cyber security challenges, typical types of data and statistical methods, paying special attention to Bayesian approaches for these problems.


title: Nodal Semimetals: A Survey on Optical Conductivity

abstract: Among different topological and related phases of condensed matter, nodal semimetals occupy a special place - the electronic band topology in these materials is related to three-dimensional bulk, rather than to surface, states. A great variety of different realizations of electronic band crossings (the nodes) leads to a plethora of different electronic properties, ranging from the chiral anomaly to solid-state realizations of a black-hole horizon. The different nodal phases have similar low-energy band structure and quasiparticle dynamics, which both can be accessed experimentally by a number of methods. Optical measurements with their large penetration depth and high energy resolution are ideally suited as such a bulk probe; especially at low energies where other spectroscopic methods often lack the required resolution. In this contribution, we review recent optical-conductivity studies of different nodal semimetals, discuss possible limitations of such measurements, and provide a comparison between the experimental results, simple theoretical models, and band-structure-based calculations.


title: Maintaining Hard Infection Caps in Epidemics via the Theory of Barriers

abstract: Research in epidemiology often focusses on designing interventions that result in the number of infected individuals asymptotically approaching zero, without considering that this number may peak at high values during transients. Recent research has shown that a set-based approach could be used to address the problem, and we build on this idea by applying the theory of barriers to construct admissible and invariant sets for an epidemic model. We describe how these sets may be used to choose intervention strategies that maintain infection caps during epidemics. We also derive algebraic conditions of the model parameters that classify a system as being either comfortable, comfortable-viable, viable, or desperate.


title: Experimental identification of the $T = 1$, $J^{\pi} = 6^+$ state of $^{54}$Co and isospin symmetry in $A = 54$ studied via one-nucleon knockout reactions

abstract: New experimental data obtained from $\gamma$-ray tagged one-neutron and one-proton knockout from $^{55}$Co is presented. A candidate for the sought-after $T=1, T_z = 0, J^{\pi} = 6^+$ state in $^{54}$Co is proposed based on a comparison to the new data on $^{54}$Fe, the corresponding observables predicted by large-scale-shell-model (LSSM) calculations in the full $fp$-model space employing charge-dependent contributions, and isospin-symmetry arguments. Furthermore, possible isospin-symmetry breaking in the $A=54$, $T=1$ triplet is studied by calculating the experimental $c$ coefficients of the isobaric mass multiplet equation (IMME) up to the maximum possible spin $J=6$ expected for the $(1f_{7/2})^{-2}$ two-hole configuration relative to the doubly-magic nucleus $^{56}$Ni. The experimental quantities are compared to the theoretically predicted $c$ coefficients from LSSM calculations using two-body matrix elements obtained from a realistic chiral effective field theory potential at next-to-next-to-next-to-leading order (N$^3$LO).


title: Optical-Field Driven Charge-Transfer Modulations near Composite Nanostructures

abstract: Optical activation of material properties illustrates the potentials held by tuning light-matter interactions with impacts ranging from basic science to technological applications. Here, we demonstrate for the first time that composite nanostructures providing nonlocal environments can be engineered to optically trigger photoinduced charge transfer dynamic (CTD) modulations in the solid state. The nanostructures herein explored lead to unprecedented out-of-phase behaviour between charge separation and recombination dynamics, along with linear CTD variations with the optical-field amplitude. Using transient absorption spectroscopy, up to 270 % increase in charge separation rate is obtained in organic semiconductor thin films. We provide evidences that composite nanostructures allow for surface photovoltages to be created, which kinetics vary with the composite architecture and last beyond optical pulse temporal characteristics. Furthermore, by generalizing Marcus theory framework, we explain why CTD modulations can only be unveiled when optic field effects are enhanced by nonlocal image dipole interactions. Demonstrating that composite nanostructures can be designed to use optical fields as CTD remote actuators opens the path for their use in practical and original applications ranging from photochemistry to optoelectronics.


title: On Interactive Machine Learning and the Potential of Cognitive Feedback

abstract: In order to increase productivity, capability, and data exploitation, numerous defense applications are experiencing an integration of state-of-the-art machine learning and AI into their architectures. Especially for defense applications, having a human analyst in the loop is of high interest due to quality control, accountability, and complex subject matter expertise not readily automated or replicated by AI. However, many applications are suffering from a very slow transition. This may be in large part due to lack of trust, usability, and productivity, especially when adapting to unforeseen classes and changes in mission context. Interactive machine learning is a newly emerging field in which machine learning implementations are trained, optimized, evaluated, and exploited through an intuitive human-computer interface. In this paper, we introduce interactive machine learning and explain its advantages and limitations within the context of defense applications. Furthermore, we address several of the shortcomings of interactive machine learning by discussing how cognitive feedback may inform features, data, and results in the state of the art. We define the three techniques by which cognitive feedback may be employed: self reporting, implicit cognitive feedback, and modeled cognitive feedback. The advantages and disadvantages of each technique are discussed.


title: Morita invariance for infinitesimal deformations

abstract: Let $A$ and $B$ be two Morita equivalent finite dimensional associative algebras over a field $\Bbbk$. It is well known that Hochschild cohomology is invariant under Morita equivalence. Since infinitesimal deformations are connected with the second Hochschild cohomology group, we explicitly describe the transfer map connecting $\mathsf{HH}^2(A)$ with $\mathsf{HH}^2(B)$. This allows us to transfer Morita equivalence between $A$ and $B$ to that between infinitesimal deformations of them. As an application, when $\Bbbk$ is algebraically closed, we consider the quotient path algebra associated to $A$ and describe the presentation by quiver and relations of the infinitesimal deformations of $A$.


title: Entropic singularities give rise to quantum transmission

abstract: When can noiseless quantum information be sent across noisy quantum devices? And at what maximum rate? These questions lie at the heart of quantum technology, but remain unanswered because of non-additivity -- a fundamental synergy which allows quantum devices (aka quantum channels) to send more information than expected. Previously, non-additivity was known to occur in very noisy channels with coherent information much smaller than that of a perfect channel; but, our work shows non-additivity in a simple low-noise channel. Our results extend even further. We prove a general theorem concerning positivity of a channel's coherent information. A corollary of this theorem gives a simple dimensional test for a channel's capacity. Applying this corollary solves an open problem by characterizing all qubit channels whose complement has non-zero capacity. Another application shows a wide class of zero quantum capacity qubit channels can assist an incomplete erasure channel in sending quantum information. These results arise from introducing and linking logarithmic singularities in the von-Neumann entropy with quantum transmission: changes in entropy caused by this singularity are a mechanism responsible for both positivity and non-additivity of the coherent information. Analysis of such singularities may be useful in other physics problems.


title: Closed 1-Forms and Twisted Cohomology

abstract: We show that the first twisted cohomology group associated to closed 1-forms on compact manifolds is related to certain 2-dimensional representations of the fundamental group. In particular, we construct examples of nowhere-vanishing 1-forms with non-trivial twisted cohomology.


title: Low Latency End-to-End Streaming Speech Recognition with a Scout Network

abstract: The attention-based Transformer model has achieved promising results for speech recognition (SR) in the offline mode. However, in the streaming mode, the Transformer model usually incurs significant latency to maintain its recognition accuracy when applying a fixed-length look-ahead window in each encoder layer. In this paper, we propose a novel low-latency streaming approach for Transformer models, which consists of a scout network and a recognition network. The scout network detects the whole word boundary without seeing any future frames, while the recognition network predicts the next subword by utilizing the information from all the frames before the predicted boundary. Our model achieves the best performance (2.7/6.4 WER) with only 639 ms latency on the test-clean and test-other data sets of Librispeech.


title: Local duality symmetry in gauge field theories

abstract: We develop the idea of local duality symmetry (LDS) in gauge field theories. Using Clifford algebra techniques we construct dually invariant scalar Lagrangian of electrodynamics in the presence of sources and demonstrate that in tensor formalism it is exactly the same as the usual one. Then we localize the duality symmetry with two possible options for the appearing pseudovector field - massive and massless. The first option might be interpreted as a candidate for dark matter. Perspectives of the application of LDS in QCD are briefly discussed.


title: Nonperturbative Dyson-Schwinger equation approach to strongly interacting Dirac fermion systems

abstract: Studying the strong correlation effects in interacting Dirac fermion systems is one of the most challenging problems in modern condensed matter physics. The long-range Coulomb interaction and the fermion-phonon interaction can lead to a variety of intriguing properties. In the strong-coupling regime, weak-coupling perturbation theory breaks down. The validity of $1/N$ expansion with $N$ being the fermion flavor is also in doubt since $N$ equals to $2$ or $4$ in realistic systems. Here, we investigate the interaction between (1+2)- and (1+3)-dimensional massless Dirac fermions and a generic scalar boson, and develop an efficient non-perturbative approach to access the strong-coupling regime. We first derive a number of self-consistently coupled Ward-Takahashi identities based on a careful symmetry analysis and then use these identities to show that the full fermion-boson vertex function is solely determined by the full fermion propagator. Making use of this result, we rigorously prove that the full fermion propagator satisfies an exact and self-closed Dyson-Schwinger integral equation, which can be solved by employing numerical methods. A major advantage of our non-perturbative approach is that there is no need to employ any small expansion parameter. Our approach provides a unified theoretical framework for studying strong Coulomb and fermion-phonon interactions. It may also be used to approximately handle the Yukawa coupling between fermions and order-parameter fluctuations around continuous quantum critical points. Our approach is applied to treat the Coulomb interaction in undoped graphene. We find that the renormalized fermion velocity exhibits a logarithmic momentum-dependence but is nearly energy independent, and that no excitonic gap is generated by the Coulomb interaction. These theoretical results are consistent with experiments in graphene.


title: Comparison of different spike train synchrony measures regarding their robustness to erroneous data from bicuculline induced epileptiform activity

abstract: As synchronized activity is associated with basic brain functions and pathological states, spike train synchrony has become an important measure to analyze experimental neuronal data. Many different measures of spike train synchrony have been proposed, but there is no gold standard allowing for comparison of results between different experiments. This work aims to provide guidance on which synchrony measure is best suitable to quantify the effect of epileptiform inducing substances (e.g. bicuculline (BIC)) in in vitro neuronal spike train data. Spike train data from recordings are likely to suffer from erroneous spike detection, such as missed spikes (false negative) or noise (false positive). Therefore, different time-scale dependent (cross-correlation, mutual information, spike time tiling coefficient) and time-scale independent (Spike-contrast, phase synchronization, A-SPIKE-synchronization, A-ISI-distance, ARI-SPIKE-distance) synchrony measures were compared in terms of their robustness to erroneous spike trains. As a result of the in silico manipulated data, Spike-contrast was the only measure being robust to false negative as well as false-positive spikes. Analyzing the experimental data set revealed that all measures were able to capture the effect of BIC in a statistically significant way, with Spike-contrast showing the highest statistical significance even at low spike detection thresholds. In summary, we suggest the usage of Spike-contrast to complement established synchrony measures, as it is time scale independent and robust to erroneous spike trains.


title: BusTime: Which is the Right Prediction Model for My Bus Arrival Time?

abstract: With the rise of big data technologies, many smart transportation applications have been rapidly developed in recent years including bus arrival time predictions. This type of applications help passengers to plan trips more efficiently without wasting unpredictable amount of waiting time at bus stops. Many studies focus on improving the prediction accuracy of various machine learning and statistical models, while much less work demonstrate their applicability of being deployed and used in realistic urban settings. This paper tries to fill this gap by proposing a general and practical evaluation framework for analysing various widely used prediction models (i.e. delay, k-nearest-neighbour, kernel regression, additive model, and recurrent neural network using long short term memory) for bus arrival time. In particular, this framework contains a raw bus GPS data pre-processing method that needs much less number of input data points while still maintain satisfactory prediction results. This pre-processing method enables various models to predict arrival time at bus stops only, by using a KD-tree based nearest point search method. Based on this framework, using raw bus GPS dataset in different scales from the city of Dublin, Ireland, we also present preliminary results for city managers by analysing the practical strengths and weaknesses in both training and predicting stages of commonly used prediction models.


title: Markovian Score Climbing: Variational Inference with KL(p||q)

abstract: Modern variational inference (VI) uses stochastic gradients to avoid intractable expectations, enabling large-scale probabilistic inference in complex models. VI posits a family of approximating distributions q and then finds the member of that family that is closest to the exact posterior p. Traditionally, VI algorithms minimize the "exclusive Kullback-Leibler (KL)" KL(q || p), often for computational convenience. Recent research, however, has also focused on the "inclusive KL" KL(p || q), which has good statistical properties that makes it more appropriate for certain inference problems. This paper develops a simple algorithm for reliably minimizing the inclusive KL using stochastic gradients with vanishing bias. This method, which we call Markovian score climbing (MSC), converges to a local optimum of the inclusive KL. It does not suffer from the systematic errors inherent in existing methods, such as Reweighted Wake-Sleep and Neural Adaptive Sequential Monte Carlo, which lead to bias in their final estimates. We illustrate convergence on a toy model and demonstrate the utility of MSC on Bayesian probit regression for classification as well as a stochastic volatility model for financial data.


title: FTT-NAS: Discovering Fault-Tolerant Convolutional Neural Architecture

abstract: With the fast evolvement of embedded deep-learning computing systems, applications powered by deep learning are moving from the cloud to the edge. When deploying neural networks (NNs) onto the devices under complex environments, there are various types of possible faults: soft errors caused by cosmic radiation and radioactive impurities, voltage instability, aging, temperature variations, and malicious attackers. Thus the safety risk of deploying NNs is now drawing much attention. In this paper, after the analysis of the possible faults in various types of NN accelerators, we formalize and implement various fault models from the algorithmic perspective. We propose Fault-Tolerant Neural Architecture Search (FT-NAS) to automatically discover convolutional neural network (CNN) architectures that are reliable to various faults in nowadays devices. Then we incorporate fault-tolerant training (FTT) in the search process to achieve better results, which is referred to as FTT-NAS. Experiments on CIFAR-10 show that the discovered architectures outperform other manually designed baseline architectures significantly, with comparable or fewer floating-point operations (FLOPs) and parameters. Specifically, with the same fault settings, F-FTT-Net discovered under the feature fault model achieves an accuracy of 86.2% (VS. 68.1% achieved by MobileNet-V2), and W-FTT-Net discovered under the weight fault model achieves an accuracy of 69.6% (VS. 60.8% achieved by ResNet-20). By inspecting the discovered architectures, we find that the operation primitives, the weight quantization range, the capacity of the model, and the connection pattern have influences on the fault resilience capability of NN models.


title: COVID-19 attack rate increases with city size

abstract: The current outbreak of novel coronavirus disease 2019 (COVID-19) poses an unprecedented global health and economic threat to interconnected human societies. Until a vaccine is developed, strategies for controlling the outbreak rely on aggressive social distancing. These measures largely disconnect the social network fabric of human societies, especially in urban areas. Here, we estimate the growth rates and reproductive numbers of COVID-19 in US cities from March 14th through March 19th to reveal a power-law scaling relationship to city population size. This means that COVID-19 is spreading faster on average in larger cities with the additional implication that, in an uncontrolled outbreak, larger fractions of the population are expected to become infected in more populous urban areas. We discuss the implications of these observations for controlling the COVID-19 outbreak, emphasizing the need to implement more aggressive distancing policies in larger cities while also preserving socioeconomic activity.


title: Extended lens reconstructions with Grale: exploiting time domain, substructural and weak-lensing information

abstract: The information about the mass density of galaxy clusters provided by the gravitational lens effect has inspired many inversion techniques. In this article, updates to the previously introduced method in Grale are described, and explored in a number of examples. The first looks into a different way of incorporating time delay information, not requiring the unknown source position. It is found that this avoids a possible bias that leads to "over-focusing" the images, i.e. providing source position estimates that lie in a considerably smaller region than the true positions. The second is inspired by previous reconstructions of the cluster of galaxies MACS J1149.6+2223, where a multiply-imaged background galaxy contained a supernova, SN Refsdal, of which four additional images were produced by the presence of a smaller cluster galaxy. The inversion for the cluster as a whole, was not able to recover sufficient detail interior to this quad. We show how constraints on such different scales, from the entire cluster to a single member galaxy, can now be used, allowing such small scale substructures to be resolved. Finally, the addition of weak lensing information to this method is investigated. While this clearly helps recover the environment around the strong lensing region, the mass sheet degeneracy may make a full strong and weak inversion difficult, depending on the quality of the ellipticity information at hand. We encounter ring-like structure at the boundary of the two regimes, argued to be the result of combining strong and weak lensing constraints, possibly affected by degeneracies.


title: Weighting NTBEA for Game AI Optimisation

abstract: The N-Tuple Bandit Evolutionary Algorithm (NTBEA) has proven very effective in optimising algorithm parameters in Game AI. A potential weakness is the use of a simple average of all component Tuples in the model. This study investigates a refinement to the N-Tuple model used in NTBEA by weighting these component Tuples by their level of information and specificity of match. We introduce weighting functions to the model to obtain Weighted- NTBEA and test this on four benchmark functions and two game environments. These tests show that vanilla NTBEA is the most reliable and performant of the algorithms tested. Furthermore we show that given an iteration budget it is better to execute several independent NTBEA runs, and use part of the budget to find the best recommendation from these runs.


title: Moment State Dynamical Systems for Nonlinear Chance-Constrained Motion Planning

abstract: Chance-constrained motion planning requires uncertainty in dynamics to be propagated into uncertainty in state. When nonlinear models are used, Gaussian assumptions on the state distribution do not necessarily apply since almost all random variables propagated through nonlinear dynamics results in non-Gaussian state distributions. To address this, recent works have developed moment-based approaches for enforcing chance-constraints on non-Gaussian state distributions. However, there still lacks fast and accurate moment propagation methods to determine the necessary statistical moments of these state distributions. To address this gap, we present a framework that, given a stochastic dynamical system, can algorithmically search for a new dynamical system in terms of moment state that can be used to propagate moments of disturbance random variables into moments of the state distribution. The key algorithm, TreeRing, can be applied to a large class of nonlinear systems which we refer to as trigonometric polynomial systems. As an example application, we present a distributionally robust RRT (DR-RRT) algorithm that propagates uncertainty through the nonlinear Dubin's car model without linearization.


title: Elliptic Equations With Degenerate weights

abstract: We obtain new local Calderon-Zygmund estimates for elliptic equations with matrix-valued weights for linear as well as non-linear equations. We introduce a novel log-BMO condition on the weight M. In particular, we assume smallness of the logarithm of the matrix-valued weight in BMO. This allows to include degenerate, discontinuous weights. We provide examples that show the sharpness of the estimates in terms of the log-BMO-norm.


title: Ambiguity in Sequential Data: Predicting Uncertain Futures with Recurrent Models

abstract: Ambiguity is inherently present in many machine learning tasks, but especially for sequential models seldom accounted for, as most only output a single prediction. In this work we propose an extension of the Multiple Hypothesis Prediction (MHP) model to handle ambiguous predictions with sequential data, which is of special importance, as often multiple futures are equally likely. Our approach can be applied to the most common recurrent architectures and can be used with any loss function. Additionally, we introduce a novel metric for ambiguous problems, which is better suited to account for uncertainties and coincides with our intuitive understanding of correctness in the presence of multiple labels. We test our method on several experiments and across diverse tasks dealing with time series data, such as trajectory forecasting and maneuver prediction, achieving promising results.


title: Thermal equilibration on the edges of topological liquids

abstract: Thermal conductance has emerged as a powerful probe of topological order in the quantum Hall effect and beyond. The interpretation of experiments crucially depends on the ratio of the sample size and the equilibration length, on which energy exchange among contra-propagating chiral modes becomes significant. We show that at low temperatures the equilibration length diverges as $1/T^2$ for almost all Abelian and non-Abelian topological orders. A faster $1/T^4$ divergence is present on the edges of the non-Abelian PH-Pfaffian and negative-flux Read-Rezayi liquids. We address experimental consequences of the $1/T^2$ and $1/T^4$ laws in a sample, shorter than the equilibration length.


title: Eigenvectors From Eigenvalues: The Case of One-Dimensional Schr\"odinger Operators

abstract: We revisit an archive submission by P. B. Denton, S. J. Parke, T. Tao, and X. Zhang, arXiv:1908.03795, on $n \times n$ self-adjoint matrices from the point of view of self-adjoint Dirichlet Schr\"odinger operators on a compact interval.


title: Implicit parametrizations in shape optimization: boundary observation

abstract: We present first a brief review of the existing literature on shape optimization, stressing the recent use of Hamiltonian systems in topology optimization. In the second section, we collect some preliminaries on the implicit parametrization theorem, especially in dimension two, which is a case of interest in shape optimization. The formulation of the problem is also discussed. The approximation via penalization and its differentiability properties are analyzed in Section 3. Next, we investigate the discretization process in Section 4. The last section is devoted to numerical experiments.


title: On the collapse of states and the uncertainty principle

abstract: There had been previous successful explanations of Quantum Mechanics using the popular Copenhagen interpretation. In this paper,we build an equivalent mathematical structure of Copenhagen Interpretation from the Electromagnetic Field Theory and show how un-certainty principle is associated with the collapse of states. We also try to provide a more intuitive method towards our understanding of Quantum Mechanics.


title: Incorporating Relational Background Knowledge into Reinforcement Learning via Differentiable Inductive Logic Programming

abstract: Relational Reinforcement Learning (RRL) can offers various desirable features. Most importantly, it allows for incorporating expert knowledge into the learning, and hence leading to much faster learning and better generalization compared to the standard deep reinforcement learning. However, most of the existing RRL approaches are either incapable of incorporating expert background knowledge (e.g., in the form of explicit predicate language) or are not able to learn directly from non-relational data such as image. In this paper, we propose a novel deep RRL based on a differentiable Inductive Logic Programming (ILP) that can effectively learn relational information from image and present the state of the environment as first order logic predicates. Additionally, it can take the expert background knowledge and incorporate it into the learning problem using appropriate predicates. The differentiable ILP allows an end to end optimization of the entire framework for learning the policy in RRL. We show the efficacy of this novel RRL framework using environments such as BoxWorld, GridWorld as well as relational reasoning for the Sort-of-CLEVR dataset.


title: Chaos-assisted tunneling resonances in a synthetic Floquet superlattice

abstract: The field of quantum simulation, which aims at using a tunable quantum system to simulate another, has been developing fast in the past years as an alternative to the all-purpose quantum computer. In particular, the use of temporal driving has attracted a huge interest recently as it was shown that certain fast drivings can create new topological effects, while a strong driving leads to e.g. Anderson localization physics. In this work, we focus on the intermediate regime to observe a quantum chaos transport mechanism called chaos-assisted tunneling which provides new possibilities of control for quantum simulation. Indeed, this regime generates a rich classical phase space where stable trajectories form islands surrounded by a large sea of unstable chaotic orbits. This mimics an effective superlattice for the quantum states localized in the regular islands, with new controllable tunneling properties. Besides the standard textbook tunneling through a potential barrier, chaos-assisted tunneling corresponds to a much richer tunneling process where the coupling between quantum states located in neighboring regular islands is mediated by other states spread over the chaotic sea. This process induces sharp resonances where the tunneling rate varies by orders of magnitude over a short range of parameters. We experimentally demonstrate and characterize these resonances for the first time in a quantum system. This opens the way to new kinds of quantum simulations with long-range transport and new types of control of quantum systems through complexity.


title: Generating Natural Language Adversarial Examples on a Large Scale with Generative Models

abstract: Today text classification models have been widely used. However, these classifiers are found to be easily fooled by adversarial examples. Fortunately, standard attacking methods generate adversarial texts in a pair-wise way, that is, an adversarial text can only be created from a real-world text by replacing a few words. In many applications, these texts are limited in numbers, therefore their corresponding adversarial examples are often not diverse enough and sometimes hard to read, thus can be easily detected by humans and cannot create chaos at a large scale. In this paper, we propose an end to end solution to efficiently generate adversarial texts from scratch using generative models, which are not restricted to perturbing the given texts. We call it unrestricted adversarial text generation. Specifically, we train a conditional variational autoencoder (VAE) with an additional adversarial loss to guide the generation of adversarial examples. Moreover, to improve the validity of adversarial texts, we utilize discrimators and the training framework of generative adversarial networks (GANs) to make adversarial texts consistent with real data. Experimental results on sentiment analysis demonstrate the scalability and efficiency of our method. It can attack text classification models with a higher success rate than existing methods, and provide acceptable quality for humans in the meantime.


title: Integration by Parts Formulae for the Laws of Bessel Bridges via Hypergeometric Functions

abstract: In this article, we extend the integration by parts formulae (IbPF) for the laws of Bessel bridges obtained in a recent work by Elad Altman and Zambotti to linear functionals. Our proof relies on properties of hypergeometric functions, thus providing a new interpretation of these formulae.


title: Generalized $f$-Sum Rules and Kohn formulas on Non-linear Conductivities

abstract: The $f$-sum rule and the Kohn formula are well-established general constraints on the electric conductivity in quantum many-body systems. We present their generalization to non-linear conductivities at all orders of the response in a unified manner, by considering two limiting quantum time-evolution processes: a quench process and an adiabatic process. Our generalized formulas are valid in any stationary state, including the ground state and finite temperature Gibbs states, regardless of the details of the system such as the specific form of the kinetic term, the strength of the many-body interactions, or the presence of disorders.

categories: cond-mat.str-el cond-mat.other cond-mat.quant-gas cond-mat.stat-mech


title: The heavier the better: how to constrain mass ratios and spins of high-mass neutron-star mergers

abstract: The first binary neutron-star merger event, GW170817, and its bright electromagnetic counterpart have provided a remarkable amount of information. By contrast, the second event, GW190425, with $M_{\rm tot}=3.4^{+0.3}_{-0.1}\,M_{\odot}$ and the lack of an electromagnetic counterpart, has hardly improved our understanding of neutron-star physics. While GW190425 is compatible with a scenario in which the merger has lead to a prompt collapse to a black hole and little ejected matter to power a counterpart, determining the mass ratio and the effective spin $\tilde{\chi}$ of the binary remains difficult. This is because gravitational waveforms cannot yet well constrain the component spins of the binary. However, since the mass of GW190425 is significantly larger than the maximum mass for nonrotating neutron stars, $M_{_{\rm TOV}}$, the mass ratio $q$ cannot be too small, as the heavier star would not be gravitationally stable. Making use of universal relations and a large number of equations of state, we provide limits in the $(\tilde{\chi},q)$ plane for GW190425, namely: $q_{\rm min} \geq 0.38$ and $ \tilde{\chi}_{\rm max}\leq 0.20$, assuming $M_\mathrm{tot} \simeq 3.4\, M_\odot$. Finally, we show how future observations of high-mass binaries can provide a lower bound on $M_{_{\rm TOV}}$.


title: Steepest Descent Neural Architecture Optimization: Escaping Local Optimum with Signed Neural Splitting

abstract: Developing efficient and principled neural architecture optimization methods is a critical challenge of modern deep learning. Recently, Liu et al.[19] proposed a splitting steepest descent (S2D) method that jointly optimizes the neural parameters and architectures based on progressively growing network structures by splitting neurons into multiple copies in a steepest descent fashion. However, S2D suffers from a local optimality issue when all the neurons become "splitting stable", a concept akin to local stability in parametric optimization. In this work, we develop a significant and surprising extension of the splitting descent framework that addresses the local optimality issue. The idea is to observe that the original S2D is unnecessarily restricted to splitting neurons into positive weighted copies. By simply allowing both positive and negative weights during splitting, we can eliminate the appearance of splitting stability in S2D and hence escape the local optima to obtain better performance. By incorporating signed splittings, we significantly extend the optimization power of splitting steepest descent both theoretically and empirically. We verify our method on various challenging benchmarks such as CIFAR-100, ImageNet and ModelNet40, on which we outperform S2D and other advanced methods on learning accurate and energy-efficient neural networks.


title: Think Globally, Act Locally: On the Optimal Seeding for Nonsubmodular Influence Maximization

abstract: We study the $r$-complex contagion influence maximization problem. In the influence maximization problem, one chooses a fixed number of initial seeds in a social network to maximize the spread of their influence. In the $r$-complex contagion model, each uninfected vertex in the network becomes infected if it has at least $r$ infected neighbors. In this paper, we focus on a random graph model named the stochastic hierarchical blockmodel, which is a special case of the well-studied stochastic blockmodel. When the graph is not exceptionally sparse, in particular, when each edge appears with probability $\omega(n^{-(1+1/r)})$, under certain mild assumptions, we prove that the optimal seeding strategy is to put all the seeds in a single community. This matches the intuition that in a nonsubmodular cascade model placing seeds near each other creates synergy. However, it sharply contrasts with the intuition for submodular cascade models (e.g., the independent cascade model and the linear threshold model) in which nearby seeds tend to erode each others' effects. Our key technique is a novel time-asynchronized coupling of four cascade processes. Finally, we show that this observation yields a polynomial time dynamic programming algorithm which outputs optimal seeds if each edge appears with a probability either in $\omega(n^{-(1+1/r)})$ or in $o(n^{-2})$.


title: Orientational Distribution of an Active Brownian Particle: an analytical study

abstract: We use the Fokker Planck equation as a starting point for studying the orientational probability distribution of an Active Brownian Particle (ABP) in $(d+1)$ dimensions. This Fokker Planck equation admits an exact solution in series form which is, however, unwieldly to use because of poor convergence for short and intermediate times. A truncated version of this series is a reasonable approximation for long times. In this paper, we present an analytical closed form expression, which gives a good approximate orientational probability distribution, which is derived using saddle point methods for short times. However, it works well even for intermediate times. Thus, we have simple analytical forms for the ${\it entire}$ range of time scales for the orientational probability distribution of an ABP. Our predictions can be tested against future experiments and simulations probing orientational probability distribution of an ABP.


title: Inverse heat source problem and experimental design for determining iron loss distribution

abstract: Iron loss determination in the magnetic core of an electrical machine, such as a motor or a transformer, is formulated as an inverse heat source problem. The sensor positions inside the object are optimized in order to minimize the uncertainty in the reconstruction in the sense of the A-optimality of Bayesian experimental design. This paper focuses on the problem formulation and an efficient numerical solution of the discretized sensor optimization and source reconstruction problems. A semirealistic linear model is discretized by finite elements and studied numerically.


title: Evaluating complexity and resilience trade-offs in emerging memory inference machines

abstract: Neuromorphic-style inference only works well if limited hardware resources are maximized properly, e.g. accuracy continues to scale with parameters and complexity in the face of potential disturbance. In this work, we use realistic crossbar simulations to highlight that compact implementations of deep neural networks are unexpectedly susceptible to collapse from multiple system disturbances. Our work proposes a middle path towards high performance and strong resilience utilizing the Mosaics framework, and specifically by re-using synaptic connections in a recurrent neural network implementation that possesses a natural form of noise-immunity.


title: Critical Point-Finding Methods Reveal Gradient-Flat Regions of Deep Network Losses

abstract: Despite the fact that the loss functions of deep neural networks are highly non-convex, gradient-based optimization algorithms converge to approximately the same performance from many random initial points. One thread of work has focused on explaining this phenomenon by characterizing the local curvature near critical points of the loss function, where the gradients are near zero, and demonstrating that neural network losses enjoy a no-bad-local-minima property and an abundance of saddle points. We report here that the methods used to find these putative critical points suffer from a bad local minima problem of their own: they often converge to or pass through regions where the gradient norm has a stationary point. We call these gradient-flat regions, since they arise when the gradient is approximately in the kernel of the Hessian, such that the loss is locally approximately linear, or flat, in the direction of the gradient. We describe how the presence of these regions necessitates care in both interpreting past results that claimed to find critical points of neural network losses and in designing second-order methods for optimizing neural networks.


title: Clocking Auger Electrons

abstract: Intense X-ray free-electron lasers (XFELs) can rapidly excite matter, leaving it in inherently unstable states that decay on femtosecond timescales. As the relaxation occurs primarily via Auger emission, excited state observations are constrained by Auger decay. In situ measurement of this process is therefore crucial, yet it has thus far remained elusive at XFELs due to inherent timing and phase jitter, which can be orders of magnitude larger than the timescale of Auger decay. Here, we develop a new approach termed self-referenced attosecond streaking, based upon simultaneous measurements of streaked photo- and Auger electrons. Our technique enables sub-femtosecond resolution in spite of jitter. We exploit this method to make the first XFEL time-domain measurement of the Auger decay lifetime in atomic neon, and, by using a fully quantum-mechanical description, retrieve a lifetime of $2.2^{ + 0.2}_{ - 0.3}$ fs for the KLL decay channel. Importantly, our technique can be generalised to permit the extension of attosecond time-resolved experiments to all current and future FEL facilities.


title: Inherent Adversarial Robustness of Deep Spiking Neural Networks: Effects of Discrete Input Encoding and Non-Linear Activations

abstract: In the recent quest for trustworthy neural networks, we present Spiking Neural Network (SNN) as a potential candidate for inherent robustness against adversarial attacks. In this work, we demonstrate that adversarial accuracy of SNNs under gradient-based attacks is higher than their non-spiking counterparts for CIFAR datasets on deep VGG and ResNet architectures, particularly in blackbox attack scenario. We attribute this robustness to two fundamental characteristics of SNNs and analyze their effects. First, we exhibit that input discretization introduced by the Poisson encoder improves adversarial robustness with reduced number of timesteps. Second, we quantify the amount of adversarial accuracy with increased leak rate in Leaky-Integrate-Fire (LIF) neurons. Our results suggest that SNNs trained with LIF neurons and smaller number of timesteps are more robust than the ones with IF (Integrate-Fire) neurons and larger number of timesteps. Also we overcome the bottleneck of creating gradient-based adversarial inputs in temporal domain by proposing a technique for crafting attacks from SNN


title: Transforming Commercial Contracts through Computable Contracting

abstract: Contracts are an essential and fundamental component of commerce and society, serving to clarify agreement between multiple parties. While digital technologies have helped to automate many activities associated with contracting, the contracts themselves continue, in the main, to be in the form of unstructured, natural-language text. This limits the scope for improvements in productivity and automation, as well as the emergence of new business models. To this end, this paper examines the concept of computable contracts as objects that are understandable by both humans and computers, and goes on to present a framework that unifies a range of technologies and approaches that collectively will help to make computable contracting a reality.


title: Learning Dynamic Routing for Semantic Segmentation

abstract: Recently, numerous handcrafted and searched networks have been applied for semantic segmentation. However, previous works intend to handle inputs with various scales in pre-defined static architectures, such as FCN, U-Net, and DeepLab series. This paper studies a conceptually new method to alleviate the scale variance in semantic representation, named dynamic routing. The proposed framework generates data-dependent routes, adapting to the scale distribution of each image. To this end, a differentiable gating function, called soft conditional gate, is proposed to select scale transform paths on the fly. In addition, the computational cost can be further reduced in an end-to-end manner by giving budget constraints to the gating function. We further relax the network level routing space to support multi-path propagations and skip-connections in each forward, bringing substantial network capacity. To demonstrate the superiority of the dynamic property, we compare with several static architectures, which can be modeled as special cases in the routing space. Extensive experiments are conducted on Cityscapes and PASCAL VOC 2012 to illustrate the effectiveness of the dynamic framework. Code is available at https://github.com/yanwei-li/DynamicRouting.


title: DevSecOps in Robotics

abstract: Quality in software is often understood as "execution according to design purpose" whereas security means that "software will not put data or computing systems at risk of unauthorized access." There seems to be a connection between these two aspects but, how do we integrate both of them in the robotics development cycle? In this article we introduce DevSecOps in Robotics, a set of best practices designed to help roboticists implant security deep in the heart of their development and operations processes. First, we briefly describe DevOps, introduce the value added with DevSecOps and describe and illustrate how these practices may be implemented in the robotics field. We finalize with a discussion on the relationship between security, quality and safety, open problems and future research questions.


title: Stress-driven spin-down of a viscous fluid within a spherical shell

abstract: We investigate the linear properties of the steady and axisymmetric stress-driven spin-down flow of a viscous fluid inside a spherical shell, both within the incompressible and anelastic approximations, and in the asymptotic limit of small viscosities. From boundary layer analysis, we derive an analytical geostrophic solution for the 3D incompressible steady flow, inside and outside the cylinder $\mathcal{C}$ that is tangent to the inner shell. The Stewartson layer that lies on $\mathcal{C}$ is composed of two nested shear layers of thickness $O(E^{2/7})$ and $O(E^{1/3})$. We derive the lowest order solution for the $E^{2/7}$-layer. A simple analysis of the $E^{1/3}$-layer laying along the tangent cylinder, reveals it to be the site of an upwelling flow of amplitude $O(E^{1/3})$. Despite its narrowness, this shear layer concentrates most of the global meridional kinetic energy of the spin-down flow. Furthermore, a stable stratification does not perturb the spin-down flow provided the Prandtl number is small enough. If this is not the case, the Stewartson layer disappears and meridional circulation is confined within the thermal layers. The scalings for the amplitude of the anelastic secondary flow have been found to be the same as for the incompressible flow in all three regions, at the lowest order. However, because the velocity no longer conforms the Taylor-Proudman theorem, its shape differs outside the tangent cylinder $\mathcal{C}$, that is, where differential rotation takes place. Finally, we find the settling of the steady-state to be reached on a viscous time for the weakly, strongly and thermally unstratified incompressible flows. Large density variations relevant to astro- and geophysical systems, tend to slightly shorten the transient.


title: Spatial Modulation for Joint Radar-Communications Systems: Design, Analysis, and Hardware Prototype

abstract: Dual-function radar-communications (DFRC) systems implement radar and communication functionalities on a single platform. Jointly designing these subsystems can lead to substantial gains in performance as well as size, cost, and power consumption. In this paper, we propose a DFRC system, which utilizes generalized spatial modulation (GSM) to realize coexisting radar and communications waveforms. Our proposed GSM-based scheme, referred to as spatial modulation based communication-radar (SpaCoR) system, allocates antenna elements among the subsystems based on the transmitted message, thus achieving increased communication rates by embedding additional data bits in the antenna selection. We formulate the resulting signal models, and present a dedicated radar processing scheme. To evaluate the radar performance, we characterize the statistical properties of the transmit beam pattern. Then, we present a hardware prototype of the proposed DFRC system, demonstrating the feasibility of the scheme. Our results show that the proposed GSM system achieves improved communication performance compared to techniques utilizing fixed allocations operating at the same data rate. For the radar subsystem, our experiments show that the spatial agility induced by the GSM transmission improves the angular resolution and reduces the sidelobe level in the transmit beam pattern compared to using fixed antenna allocations.


title: Mechanical Ventilator Milano (MVM): A Novel Mechanical Ventilator Designed for Mass Scale Production in Response to the COVID-19 Pandemic

abstract: Presented here is the design of the Mechanical Ventilator Milano (MVM), a novel mechanical ventilator designed for rapid mass production in response to the COVID-19 pandemic to address the urgent shortage of intensive therapy ventilators in many countries, and the growing difficulty in procuring these devices through normal supply chains across borders. This ventilator is an electro-mechanical equivalent of the old and reliable Manley Ventilator, and is able to operate in both pressure-controlled and pressure-supported ventilation modes. MVM is optimized for the COVID-19 emergency, thanks to the collaboration with medical doctors in the front line. MVM is designed for large-scale production in a short amount of time and at a limited cost, as it relays on off-the-shelf components, readily available worldwide. Operation of the MVM requires only a source of compressed oxygen (or compressed medical air) and electrical power. Initial tests of a prototype device with a breathing simulator are also presented. Further tests and developments are underway. At this stage the MVM is not yet a certified medical device but certification is in progress.


title: Dynamic load balancing with enhanced shared-memory parallelism for particle-in-cell codes

abstract: Furthering our understanding of many of today's interesting problems in plasma physics---including plasma based acceleration and magnetic reconnection with pair production due to quantum electrodynamic effects---requires large-scale kinetic simulations using particle-in-cell (PIC) codes. However, these simulations are extremely demanding, requiring that contemporary PIC codes be designed to efficiently use a new fleet of exascale computing architectures. To this end, the key issue of parallel load balance across computational nodes must be addressed. We discuss the implementation of dynamic load balancing by dividing the simulation space into many small, self-contained regions or "tiles," along with shared-memory (e.g., OpenMP) parallelism both over many tiles and within single tiles. The load balancing algorithm can be used with three different topologies, including two space-filling curves. We tested this implementation in the code OSIRIS and show low overhead and improved scalability with OpenMP thread number on simulations with both uniform load and severe load imbalance. Compared to other load-balancing techniques, our algorithm gives order-of-magnitude improvement in parallel scalability for simulations with severe load imbalance issues.


title: Asymptotic horizon formation, spacetime stretching and causality

abstract: In this work we analyse asymptotically flat, spherically symmetric spacetimes in which an event horizon is present without any trapped surfaces. We identify two types of such spacetimes, each related to the asymptotic behaviour (in time) of one of the two degrees of freedom of the metric. We study the causal structure of both types, showing that one almost always has a Cauchy horizon beyond which it is extendable, while the other is inextendable but has two separate future null infinity regions on either side of the horizon. We also study what energy conditions can be satisfied by the matter around the horizon. Some of these spacetimes were first introduced in an earlier work in which semiclassical effects near black-hole horizons were analysed. Here we generalise this analysis to a larger family of geometries.


title: Countable graphs are majority 3-choosable

abstract: The Unfriendly Partition Conjecture posits that every countable graph admits a 2-colouring in which for each vertex there are at least as many bichromatic edges containing that vertex as monochromatic ones. This is not known in general, but it is known that a 3-colouring with this property always exists. Anholcer, Bosek and Grytczuk recently gave a list-colouring version of this conjecture, and proved that such a colouring exists for lists of size 4. We improve their result to lists of size 3; the proof extends to directed acyclic graphs. We also discuss some generalisations.


title: Online stochastic gradient descent on non-convex losses from high-dimensional inference

abstract: Stochastic gradient descent (SGD) is a popular algorithm for optimization problems arising in high-dimensional inference tasks. Here one produces an estimator of an unknown parameter from independent samples of data by iteratively optimizing a loss function. This loss function is random and often non-convex. We study the performance of the simplest version of SGD, namely online SGD, from a random start in the setting where the parameter space is high-dimensional. We develop nearly sharp thresholds for the number of samples needed for consistent estimation as one varies the dimension. Our thresholds depend only on an intrinsic property of the population loss which we call the information exponent. In particular, our results do not assume uniform control on the loss itself, such as convexity or uniform derivative bounds. The thresholds we obtain are polynomial in the dimension and the precise exponent depends explicitly on the information exponent. As a consequence of our results, we find that except for the simplest tasks, almost all of the data is used simply in the initial search phase to obtain non-trivial correlation with the ground truth. Upon attaining non-trivial correlation, the descent is rapid and exhibits law of large numbers type behavior. We illustrate our approach by applying it to a wide set of inference tasks such as phase retrieval, and parameter estimation for generalized linear models, online PCA, and spiked tensor models, as well as to supervised learning for single-layer networks with general activation functions.


title: A Critique of Covariant Emergent Gravity

abstract: I address some problems encountered in the formulation of relativistic models encompassing the MOND phenomenology of radial acceleration. I explore scalar and vector theories with fractional kinetic terms and $f(R)$-type gravity, demanding that the energy density be bounded from below and that superluminal modes be absent, but also that some consistency constraints with observational results hold. I identify configurations whose energy is unbounded from below and formulate some no-go statements for vector field models and modified gravity theories. Finally, I discuss superfluid dark matter as a hybrid theory lying between CDM and modified gravity, highlighting some difficulties present also in this case, which appears preferable to the others.