title: On the $\mathrm{EO}$-orientability of vector bundles

abstract: We study the orientability of vector bundles with respect to a family of cohomology theories called $\mathrm{EO}$-theories. The $\mathrm{EO}$-theories are higher height analogues of real $\mathrm{K}$-theory $\mathrm{KO}$. For each $\mathrm{EO}$-theory, we prove that the direct sum of $i$ copies of any vector bundle is $\mathrm{EO}$-orientable for some specific integer $i$. Using a splitting principal, we reduce to the case of the canonical line bundle over $\mathbb{CP}^{\infty}$. Our method involves understanding the action of an order $p$ subgroup of the Morava stabilizer group on the Morava $\mathrm{E}$-theory of $\mathbb{CP}^{\infty}$. Our calculations have another application: We determine the homotopy type of the $\mathrm{S}^{1}$-Tate spectrum associated to the trivial action of $\mathrm{S}^{1}$ on all $\mathrm{EO}$-theories.


title: Semi-Inclusive Jet Functions and Jet Substructure in $J_{E_T}^{(I)}$ and $J_{E_T}^{(II)}$ Algorithms

abstract: Within the framework of Soft Collinear Effective Theory, we present calculations of semi-inclusive jet functions and fragmenting jet functions at next-to-leading order (NLO) for both quark- and gluon-initiated jets, for jet algorithms of $J_{E_T}^{(I)}$ and $J_{E_T}^{(II)}$ where one maximizes a suitable jet function. We demonstrate the consistency of the obtained results with the standard perturbative QCD calculations for $J_{E_T}^{(I)}$ algorithm, while the results for fragmenting jet functions with the $J_{E_T}^{(II)}$ algorithm are new. The renormalization group (RG) equation for both semi-inclusive jet functions and fragmenting jet functions are derived and shown to follow the time-like DGLAP evolution equations, independent of specific jet algorithms. The RG equation can be used to resum single logarithms of the jet size parameter $\beta$ for highly collimated jets in these algorithms where $\beta \gg 1$.


title: 1D Probabilistic Undersampling Pattern Optimization for MR Image Reconstruction

abstract: Magnetic resonance imaging (MRI) is mainly limited by long scanning time and vulnerable to human tissue motion artifacts, in 3D clinical scenarios. Thus, k-space undersampling is used to accelerate the acquisition of MRI while leading to visually poor MR images. Recently, some studies 1) use effective undersampling patterns, or 2) design deep neural networks to improve the quality of resulting images. However, they are considered as two separate optimization strategies. In this paper, we propose a cross-domain network for MR image reconstruction, in a retrospective data-driven manner, under limited sampling rates. Our method can simultaneously obtain the optimal undersampling pattern (in k-space) and the reconstruction model, which are customized to the type of training data, by using an end-to-end learning strategy. We propose a 1D probabilistic undersampling layer, to obtain the optimal undersampling pattern and its probability distribution in a differentiable way. We propose a 1D inverse Fourier transform layer, which connects the Fourier domain and the image domain during the forward pass and the backpropagation. In addition, by training 3D fully-sampled k-space data and MR images with the traditional Euclidean loss, we discover the universal relationship between the probability distribution of the optimal undersampling pattern and its corresponding sampling rate. Experiments show that the quantitative and qualitative results of recovered MR images by our 1D probabilistic undersampling pattern obviously outperform those of several existing sampling strategies.


title: Efficient mate finding in planktonic copepods swimming in turbulence

abstract: Zooplankton live in dynamic environments where turbulence may challenge their limited swimming abilities. How this interferes with fundamental behavioral processes remains elusive. We reconstruct simultaneously the trajectories of flow tracers and calanoid copepods and we quantify their ability to find mates when ambient flow imposes physical constrains on their motion and impairs their olfactory orientation. We show that copepods achieve high encounter rates in turbulence due to the contribution of advection and vigorous swimming. Males further convert encounters within the perception radius to contacts and then to mating via directed motion toward nearby organisms within the short time frame of the encounter. Inertial effects do not result in preferential concentration, reducing the geometric collision kernel to the clearance rate, which we model accurately by superposing turbulent velocity and organism motion. This behavioral and physical coupling mechanism may account for the ability of copepods to reproduce in turbulent environments.


title: Continuous Experimentation for Automotive Software on the Example of a Heavy Commercial Vehicle in Daily Operation

abstract: As the automotive industry focuses its attention more and more towards the software functionality of vehicles, techniques to deliver new software value at a fast pace are needed. Continuous Experimentation, a practice coming from the web-based systems world, is one of such techniques. It enables researchers and developers to use real-world data to verify their hypothesis and steer the software evolution based on performances and user preferences, reducing the reliance on simulations and guesswork. Several challenges prevent the verbatim adoption of this practice on automotive cyber-physical systems, e.g., safety concerns and limitations from computational resources; nonetheless, the automotive field is starting to take interest in this technique. This work aims at demonstrating and evaluating a prototypical Continuous Experimentation infrastructure, implemented on a distributed computational system housed in a commercial truck tractor that is used in daily operations by a logistic company on public roads. The system comprises computing units and sensors, and software deployment and data retrieval are only possible remotely via a mobile data connection due to the commercial interests of the logistics company. This study shows that the proposed experimentation process resulted in the development team being able to base software development choices on the real-world data collected during the experimental procedure. Additionally, a set of previously identified design criteria to enable Continuous Experimentation on automotive systems was discussed and their validity confirmed in the light of the presented work.


title: One or two small points in thermodynamics

abstract: I present my recollections of what I used to find to be "one or two small points in thermodynamics", following Sommerfeld's famous quote, and review them on the light of present knowledge.


title: Multiset Synchronization with Counting Cuckoo Filters

abstract: Set synchronization is a fundamental task in distributed applications and implementations. Existing methods that synchronize simple sets are mainly based on compact data structures such as Bloom filter and its variants. However, these methods are infeasible to synchronize a pair of multisets which allow an element to appear for multiple times. To this end, in this paper, we propose to leverage the counting cuckoo filter (CCF), a novel variant of cuckoo filter, to represent and thereafter synchronize a pair of multisets. The cuckoo filter (CF) is a minimized hash table that uses cuckoo hashing to resolve collisions. CF has an array of buckets, each of which has multiple slots to store element fingerprints. Based on CF, CCF extends each slot as two fields, the fingerprint field and the counter field. The fingerprint field records the fingerprint of element which is stored by this slot; while the counter field counts the multiplicity of the stored element. With such a design, CCF is competent to represent any multiset. After generating and exchanging the respective CCFs which represent the local multi-sets, we propose the query-based and the decoding-based methods to identify the different elements between the given multisets. The comprehensive evaluation results indicate that CCF outperforms the counting Bloom filter (CBF) when they are used to synchronize multisets, in terms of both synchronization accuracy and the space-efficiency, at the cost of a little higher time-consumption.


title: Probabilistic conformal blocks for Liouville CFT on the torus

abstract: Virasoro conformal blocks are a family of important functions defined as power series via the Virasoro algebra. They are a fundamental input to the conformal bootstrap program for 2D conformal field theory (CFT) and are closely related to four dimensional supersymmetric gauge theory through the Alday-Gaiotto-Tachikawa correspondence. The present work provides a probabilistic construction of the 1-point toric Virasoro conformal block for central change greater than 25. More precisely, we construct an analytic function using a probabilistic tool called Gaussian multiplicative chaos (GMC) and prove that its power series expansion coincides with the 1-point toric Virasoro conformal block. The range $(25,\infty)$ of central charges corresponds to Liouville CFT, an important CFT originating from 2D quantum gravity and bosonic string theory. Our work reveals a new integrable structure underlying GMC and opens the door to the study of non-perturbative properties of Virasoro conformal blocks such as their analytic continuation and modular symmetry. Our proof combines an analysis of GMC with tools from CFT such as Belavin-Polyakov-Zamolodchikov differential equations, operator product expansions, and Dotsenko-Fateev type integrals.


title: Lagrangian schemes for Wasserstein gradient flows

abstract: This paper reviews different numerical methods for specific examples of Wasserstein gradient flows: we focus on nonlinear Fokker-Planck equations,but also discuss discretizations of the parabolic-elliptic Keller-Segel model and of the fourth order thin film equation. The methods under review are of Lagrangian nature, that is, the numerical approximations trace the characteristics of the underlying transport equation rather than solving the evolution equation for the mass density directly. The two main approaches are based on integrating the equation for the Lagrangian maps on the one hand, and on solution of coupled ODEs for individual mass particles on the other hand.


title: Continuous Experimentation and the Cyber-Physical Systems challenge: An overview of the literature and the industrial perspective

abstract: Context: New software development patterns are emerging aiming at accelerating the process of delivering value. One is Continuous Experimentation, which allows to systematically deploy and run instrumented software variants during development phase in order to collect data from the field of application. While currently this practice is used on a daily basis on web-based systems, technical difficulties challenge its adoption in fields where computational resources are constrained, e.g., cyber-physical systems and the automotive industry. Objective: This paper aims at providing an overview of the engagement on the Continuous Experimentation practice in the context of cyber-physical systems. %To provide an understanding of what is the state-of-the-art of the Continuous Experimentation practice in the context of cyber-physical systems, and what is the practitioners' feedback about this practice. Method: A systematic literature review has been conducted to investigate the link between the practice and the field of application. Additionally, an industrial multiple case study is reported. Results: The study presents the current state-of-the-art regarding Continuous Experimentation in the field of cyber-physical systems. The current perspective of Continuous Experimentation in industry is also reported. Conclusions: The field has not reached maturity yet. More conceptual analyses are found than solution proposals and the state-of-practice is yet to be achieved. However it is expected that in time an increasing number of solutions will be proposed and validated.


title: BCOV invariant and blow-up

abstract: Bershadsky, Cecotti, Ooguri and Vafa constructed a real valued invariant for Calabi-Yau manifolds, which is called the BCOV invariant. In this paper, we extend the BCOV invariant to such pairs $(X,D)$, where $X$ is a compact K\"ahler manifold and $D$ is a pluricanonical divisor on X with simple normal crossing support. We also study the behavior of the extended BCOV invariant under blow-up. The results in this paper lead to a joint work with Fu proving that birational Calabi-Yau manifolds have the same BCOV invariant.


title: A model in one-dimensional thermoelasticity

abstract: We study a one-dimensional nonlinear hyperbolic-parabolic initial boundary value problem occurring in the theory of thermoelasticity. We prove existence and uniqueness of the local-in-time strong solution. Also, some global-in-time weak measure valued solutions are proven to exist. To this end we introduce an auxiliary problem with artificial viscosity and prove its global-in-time well-posedness. Next, we show that solutions of the auxiliary problem converge, at some short time interval to the strong solution, and to our measure valued solution for an arbitrary time.


title: Task-Motion Planning for Safe and Efficient Urban Driving

abstract: Autonomous vehicles need to plan at the task level to compute a sequence of symbolic actions, such as merging left and turning right, to fulfill people's service requests, where efficiency is the main concern. At the same time, the vehicles must compute continuous trajectories to perform actions at the motion level, where safety is the most important. Task-motion planning in autonomous driving faces the problem of maximizing task-level efficiency while ensuring motion-level safety. To this end, we develop algorithm Task-Motion Planning for Urban Driving (TMPUD) that, for the first time, enables the task and motion planners to communicate about the safety level of driving behaviors. TMPUD has been evaluated using a realistic urban driving simulation platform. Results suggest that TMPUD performs significantly better than competitive baselines from the literature in efficiency, while ensuring the safety of driving behaviors.


title: PULSE: Self-Supervised Photo Upsampling via Latent Space Exploration of Generative Models

abstract: The primary aim of single-image super-resolution is to construct high-resolution (HR) images from corresponding low-resolution (LR) inputs. In previous approaches, which have generally been supervised, the training objective typically measures a pixel-wise average distance between the super-resolved (SR) and HR images. Optimizing such metrics often leads to blurring, especially in high variance (detailed) regions. We propose an alternative formulation of the super-resolution problem based on creating realistic SR images that downscale correctly. We present an algorithm addressing this problem, PULSE (Photo Upsampling via Latent Space Exploration), which generates high-resolution, realistic images at resolutions previously unseen in the literature. It accomplishes this in an entirely self-supervised fashion and is not confined to a specific degradation operator used during training, unlike previous methods (which require supervised training on databases of LR-HR image pairs). Instead of starting with the LR image and slowly adding detail, PULSE traverses the high-resolution natural image manifold, searching for images that downscale to the original LR image. This is formalized through the "downscaling loss," which guides exploration through the latent space of a generative model. By leveraging properties of high-dimensional Gaussians, we restrict the search space to guarantee realistic outputs. PULSE thereby generates super-resolved images that both are realistic and downscale correctly. We show proof of concept of our approach in the domain of face super-resolution (i.e., face hallucination). We also present a discussion of the limitations and biases of the method as currently implemented with an accompanying model card with relevant metrics. Our method outperforms state-of-the-art methods in perceptual quality at higher resolutions and scale factors than previously possible.


title: Half-space stationary Kardar-Parisi-Zhang equation

abstract: We study the solution of the Kardar-Parisi-Zhang (KPZ) equation for the stochastic growth of an interface of height $h(x,t)$ on the positive half line, equivalently the free energy of the continuum directed polymer in a half space with a wall at $x=0$. The boundary condition $\partial_x h(x,t)|_{x=0}=A$ corresponds to an attractive wall for $A<0$, and leads to the binding of the polymer to the wall below the critical value $A=-1/2$. Here we choose the initial condition $h(x,0)$ to be a Brownian motion in $x>0$ with drift $-(B+1/2)$. When $A+B \to -1$, the solution is stationary, i.e. $h(\cdot,t)$ remains at all times a Brownian motion with the same drift, up to a global height shift $h(0,t)$. We show that the distribution of this height shift is invariant under the exchange of parameters $A$ and $B$. For any $A,B > - 1/2$, we provide an exact formula characterizing the distribution of $h(0,t)$ at any time $t$, using two methods: the replica Bethe ansatz and a discretization called the log-gamma polymer, for which moment formulae were obtained. We analyze its large time asymptotics for various ranges of parameters $A,B$. In particular, when $(A, B) \to (-1/2, -1/2)$, the critical stationary case, the fluctuations of the interface are governed by a universal distribution akin to the Baik-Rains distribution arising in stationary growth on the full-line. It can be expressed in terms of a simple Fredholm determinant, or equivalently in terms of the Painlev\'e II transcendent. This provides an analog for the KPZ equation, of some of the results recently obtained by Betea-Ferrari-Occelli in the context of stationary half-space last-passage-percolation. From universality, we expect that limiting distributions found in both models can be shown to coincide.

categories: cond-mat.stat-mech cond-mat.dis-nn math-ph math.MP math.PR nlin.SI


title: Attacking the DeFi Ecosystem with Flash Loans for Fun and Profit

abstract: Credit allows a lender to loan out surplus capital to a borrower. In the traditional economy, credit bears the risk that the borrower may default on its debt, the lender hence requires upfront collateral from the borrower, plus interest fee payments. Due to the atomicity of blockchain transactions, lenders can offer flash loans, i.e., loans that are only valid within one transaction and must be repaid by the end of that transaction. This concept has lead to a number of interesting attack possibilities, some of which were exploited in February 2020. This paper is the first to explore the implication of transaction atomicity and flash loans for the nascent decentralized finance (DeFi) ecosystem. We show quantitatively how transaction atomicity increases the arbitrage revenue. We moreover analyze two existing attacks with ROIs beyond 500k%. We formulate finding the attack parameters as an optimization problem over the state of the underlying Ethereum blockchain and the state of the DeFi ecosystem. We show how malicious adversaries can efficiently maximize an attack profit and hence damage the DeFi ecosystem further. Specifically, we present how two previously executed attacks can be "boosted" to result in a profit of 829.5k USD and 1.1M USD, respectively, which is a boost of 2.37x and 1.73x, respectively.


title: ASAP-SML: An Antibody Sequence Analysis Pipeline Using Statistical Testing and Machine Learning

abstract: Antibodies are capable of potently and specifically binding individual antigens and, in some cases, disrupting their functions. The key challenge in generating antibody-based inhibitors is the lack of fundamental information relating sequences of antibodies to their unique properties as inhibitors. We develop a pipeline, Antibody Sequence Analysis Pipeline using Statistical testing and Machine Learning (ASAP-SML), to identify features that distinguish one set of antibody sequences from antibody sequences in a reference set. The pipeline extracts feature fingerprints from sequences. The fingerprints represent germline, CDR canonical structure, isoelectric point and frequent positional motifs. Machine learning and statistical significance testing techniques are applied to antibody sequences and extracted feature fingerprints to identify distinguishing feature values and combinations thereof. To demonstrate how it works, we applied the pipeline on sets of antibody sequences known to bind or inhibit the activities of matrix metalloproteinases (MMPs), a family of zinc-dependent enzymes that promote cancer progression and undesired inflammation under pathological conditions, against reference datasets that do not bind or inhibit MMPs. ASAP-SML identifies features and combinations of feature values found in the MMP-targeting sets that are distinct from those in the reference sets.


title: Implicit and Coupled Multi-Fluid Solver for Collisional Low-Temperature Plasma

abstract: We present a new multi-fluid, multi-temperature plasma solver with adaptive Cartesian mesh (ACM) based on a full-Newton (non-linear, implicit) scheme for collisional low-temperature plasma. The particle transport is described using the drift-diffusion approximation for electrons and ions coupled to Poisson equation for electric field. In addition, the electron-energy transport equation is solved to account for electron thermal conductivity, Joule heating, and energy loss of electrons in collisions with neutral species. The rate of electron-induced ionization is a function of electron temperature and could also depend on electron density (important for plasma stratification). The ion and gas temperature are kept constant. The spatial discretization of the transport equations uses non-isothermal Scharfetter-Gummel scheme from semiconductor physics adapted for multi-dimensional ACM framework. We demonstrate the new solver for simulations of direct current (DC) and radio frequency (RF) discharges. The implicit treatment of the coupled equations allows using large time steps, and the full-Newton method enables fast non-linear convergence at each time step, offering greatly improved efficiency of fluid plasma simulations. We discuss the selection of time steps for solving different plasma problems. The new solver enables us to solve several problems we could not solve before with existing software: two- and three-dimensional structures of the entire DC discharges including cathode and anode regions with electric field reversals, normal cathode spot and anode ring, plasma stratification in diffuse and constricted DC discharges, and standing striations in RF discharges.


title: Keeping it simple: Implementation and performance of the proto-principle of adaptation and learning in the language sciences

abstract: In this paper we present the Widrow-Hoff rule and its applications to language data. After contextualizing the rule historically and placing it in the chain of neurally inspired artificial learning models, we explain its rationale and implementational considerations. Using a number of case studies we illustrate how the Widrow-Hoff rule offers unexpected opportunities for the computational simulation of a range of language phenomena that make it possible to approach old problems from a novel perspective.


title: Enhancing Industrial X-ray Tomography by Data-Centric Statistical Methods

abstract: X-ray tomography has applications in various industrial fields such as sawmill industry, oil and gas industry, chemical engineering, and geotechnical engineering. In this article, we study Bayesian methods for the X-ray tomography reconstruction. In Bayesian methods, the inverse problem of tomographic reconstruction is solved with help of a statistical prior distribution which encodes the possible internal structures by assigning probabilities for smoothness and edge distribution of the object. We compare Gaussian random field priors, that favour smoothness, to non-Gaussian total variation, Besov, and Cauchy priors which promote sharp edges and high-contrast and low-contrast areas in the object. We also present computational schemes for solving the resulting high-dimensional Bayesian inverse problem with 100,000-1,000,000 unknowns. In particular, we study the applicability of a no-U-turn variant of Hamiltonian Monte Carlo methods and of a more classical adaptive Metropolis-within-Gibbs algorithm for this purpose. These methods also enable full uncertainty quantification of the reconstructions. For faster computations, we use maximum a posteriori estimates with limited-memory BFGS optimisation algorithm. As the first industrial application, we consider sawmill industry X-ray log tomography. The logs have knots, rotten parts, and even possibly metallic pieces, making them good examples for non-Gaussian priors. Secondly, we study drill-core rock sample tomography, an example from oil and gas industry. We show that Cauchy priors produce smaller number of artefacts than other choices, especially with sparse high-noise measurements, and choosing Hamiltonian Monte Carlo enables systematic uncertainty quantification.


title: Models for knot spaces and Atiyah duality

abstract: Let $\mathrm{Emb}(S^1,M)$ be the space of smooth embeddings from the circle to a closed manifold $M$ of dimension $\geq 4$. We study a cosimplicial model of $\mathrm{Emb}(S^1,M)$ in stable categories, using a spectral version of Poincar\'e-Lefschetz duality called Atiyah duality. We actually deal with a notion of a comodule instead of the cosimplicial model, and prove a comodule version of the duality. As an application, we introduce a new spectral sequence converging to $H^*(\mathrm{Emb}(S^1,M))$ for simply connected $M$ and for major coefficient rings. Using this, we compute $H^*(\mathrm{Emb}(S^1, S^k\times S^l))$ in low degrees with some conditions on $k$, $l$. We also prove the inclusion $\mathrm{Emb}(S^1,M)\to \mathrm{Imm}(S^1,M)$ to the immersions induces an isomorphism on $\pi_1$ for some simply connected $4$-manifolds, related to a question posed by Arone and Szymik. We also prove an equivalence of singular cochain complex of $\mathrm{Emb}(S^1,M)$ and a homotopy colimit of chain complexes of a Thom spectrum of a bundle over a comprehensible space. Our key ingredient is a structured version of the duality due to R. Cohen.


title: Reinterpretation of classic proton charge form factor measurements

abstract: In 1963, a proton radius of $0.805(11)~\mathrm{fm}$ was extracted from electron scattering data and this classic value has been used in the standard dipole parameterization of the form factor. In trying to reproduce this classic result, we discovered that there was a sign error in the original analysis and that the authors should have found a value of $0.851(19)~\mathrm{fm}$. We additionally made use of modern computing power to find a robust function for extracting the radius using this 1963 data's spacing and uncertainty. This optimal function, the Pad\'{e} $(0,1)$ approximant, also gives a result which is consistent with the modern high precision proton radius extractions.


title: On the $ \rho $ and $ \sigma_{tot} $ measurement by the TOTEM Collaboration: in the wake of recent discoveries

abstract: We show that extraction of the quantities} $\rho(s)=\mbox{Re}[T_{N}(s,0)]/\mbox{Im}[T_{N}(s,0)]$ {\it and} $\sigma_{\mbox{tot}}(s)$ {\it from the data on the $pp$ differential cross-section at $\sqrt{s} = 13$} TeV {\it obtained by the} TOTEM {\it Collaboration gives results essentially different from those presented in publication} \cite{Ant} {\it if to use a modified formula for Coulomb-nuclear interference. The physical interpretation of these data changes accordingly.


title: Quantum versus classical approach of dechanneling and incoherent electromagnetic processes in aligned crystals

abstract: Particles traveling in aligned crystals at small angles w.r.t. crystallographic axes or planes are principally steered by the continuous Lindhard potential. This interaction conserves the energy E, the longitudinal momentum p_parallel, the transverse energy of the particle E_perp and is elastic concerning the crystal quantum state. At high enough energy the particle motion is quasi-classical. The time-dependent fluctuations of the positions of the atoms or of the electrons of the crystal create a residual potential, on which the particle can scatter. This interaction does not conserve the previous quantities and is inelastic for the crystal. We compare its treatments with the classical binary collision model and with a phenomenological quantum model. The classical dechanneling rate is estimated to be several ten per cent larger than the quantum one. The influence of correlated vibrations of neighboring atoms is discussed.


title: An energy stable one-field monolithic arbitrary Lagrangian-Eulerian formulation for fluid-structure interaction

abstract: In this article we present a one-field monolithic finite element method in the Arbitrary Lagrangian-Eulerian (ALE) formulation for Fluid-Structure Interaction (FSI) problems. The method only solves for one velocity field in the whole FSI domain, and it solves in a monolithic manner so that the fluid solid interface conditions are satisfied automatically. We prove that the proposed scheme is unconditionally stable, through energy analysis, by utilising a conservative formulation and an exact quadrature rule. We implement the algorithm using both ${\bf F}$-scheme and ${\bf d}$-scheme, and demonstrate that the former has the same formulation in two and three dimensions. Finally several numerical examples are presented to validate this methodology, including combination with remesh techniques to handle the case of very large solid displacement.


title: Classical spectral sum rules and "half-naked" electron effects in radiation from relativistic electrons in external field

abstract: Two properties of the radiation emitted by a relativistic electron in an external field, in the classical approximation, are presented in details: 1) spectral sum rules and their relationship with the sum rules for oscillator strength in atomic physics. Filtered sum rules can be used to remove infrared or ultraviolet divergences. 2) "half-naked electron" effects like the Landau- Pomeranchuk-Migdal and Ternovskij-Shul'ga-Fomin effects. These properties are best understood when introducing the apparent transverse velocity.


title: Dirac-Born-Infeld realization of sound speed resonance mechanism for primordial black holes

abstract: We present a concrete realization of the sound speed resonance (SSR) mechanism for primordial black hole (PBH) formation within a specific model of Dirac-Born-Infeld (DBI) inflation. We perform a perturbative approach to phenomenologically construct such a viable DBI inflation model that involves the nonoscillating stage and the oscillating stage, with a type of specific forms of the warp factor and the potential. We show that the continuous but nonsmooth conjunction of sound speed between two stages does not yield manifest effects on the phenomenology of SSR, and thus, our model gives rise to the same PBH mass spectrum as the original predictions of SSR. Additionally, we also demonstrate that the violation of adiabaticity of the Mukhanov-Sasaki equation does not affect the comoving curvature perturbation after Hubble crossing in the nonresonant region. Making use of observational data, we derive various cosmological constraints on the parameter space. Our analyses show that the predicted tensor-to-scalar ratio is typically small, while the amplitude of primordial non-Gaussianity can meet with cosmic microwave background bounds, and additionally, the consistency relation for single-field slow-roll inflation is softly violated in our case due to the small sound speed variations.


title: Inverse Scattering for Critical Semilinear Wave Equations

abstract: We show that the scattering operator for defocusing energy critical semilinear wave equations \square u+f(u)=0, with f C-infinity and f ~ u^5, in three space dimensions, determines f.


title: On Adiabatic Oscillations of a Stratified Atmosphere on the Flat Earth

abstract: We consider the oscillations of the atmosphere around a stratified back ground density and entropy distribution under the gravitation on the flat Earth. The atmosphere is supposed to be an ideal gas and the motion is supposed to be governed by the compressible Euler equations. The density distribution of the back ground equilibrium is supposed to touch the vacuum at the finite height of the stratosphere. Considering the linearized approximation for small perturbations, we show that time periodic oscillations with a sequence of time periods which accumulate to infinity, say, slow and slow oscillations, so called `g-modes', can appear when the square of the Brunt-Vaisala frequency is positive everywhere for the considered back ground equilibrium.


title: No Surprises: Training Robust Lung Nodule Detection for Low-Dose CT Scans by Augmenting with Adversarial Attacks

abstract: Detecting malignant pulmonary nodules at an early stage can allow medical interventions which may increase the survival rate of lung cancer patients. Using computer vision techniques to detect nodules can improve the sensitivity and the speed of interpreting chest CT for lung cancer screening. Many studies have used CNNs to detect nodule candidates. Though such approaches have been shown to outperform the conventional image processing based methods regarding the detection accuracy, CNNs are also known to be limited to generalize on under-represented samples in the training set and prone to imperceptible noise perturbations. Such limitations can not be easily addressed by scaling up the dataset or the models. In this work, we propose to add adversarial synthetic nodules and adversarial attack samples to the training data to improve the generalization and the robustness of the lung nodule detection systems. To generate hard examples of nodules from a differentiable nodule synthesizer, we use projected gradient descent (PGD) to search the latent code within a bounded neighbourhood that would generate nodules to decrease the detector response. To make the network more robust to unanticipated noise perturbations, we use PGD to search for noise patterns that can trigger the network to give over-confident mistakes. By evaluating on two different benchmark datasets containing consensus annotations from three radiologists, we show that the proposed techniques can improve the detection performance on real CT data. To understand the limitations of both the conventional networks and the proposed augmented networks, we also perform stress-tests on the false positive reduction networks by feeding different types of artificially produced patches. We show that the augmented networks are more robust to both under-represented nodules as well as resistant to noise perturbations.


title: Energy resolution and linearity of XENON1T in the MeV energy range

abstract: Xenon dual-phase time projection chambers designed to search for Weakly Interacting Massive Particles have so far shown a relative energy resolution which degrades with energy above $\sim$200 keV due to the saturation effects. This has limited their sensitivity in the search for rare events like the neutrinoless double-beta decay of $^{136}$Xe at its $Q$-value, $Q_{\beta\beta}\simeq$ 2.46 MeV. For the XENON1T dual-phase time projection chamber, we demonstrate that the relative energy resolution at 1 $\sigma/\mu$ is as low as (0.80$\pm$0.02) % in its one-ton fiducial mass, and for single-site interactions at $Q_{\beta\beta}$. We also present a new signal correction method to rectify the saturation effects of the signal readout system, resulting in more accurate position reconstruction and indirectly improving the energy resolution. The very good result achieved in XENON1T opens up new windows for the xenon dual-phase dark matter detectors to simultaneously search for other rare events.

categories: physics.ins-det astro-ph.IM hep-ex nucl-ex


title: On the existence of balanced metrics on six-manifolds of cohomogeneity one

abstract: We consider balanced metrics on complex manifolds with holomorphically trivial canonical bundle, most commonly known as balanced $\rm{SU}(n)$-structures. Such structures are of interest for both Hermitian geometry and string theory, since they provide the ideal setting for the Hull-Strominger system. In this paper, we provide a non-existence result for balanced non-K\"ahler $\rm{SU}(3)$-structures which are invariant under a cohomogeneity one action on simply connected six-manifolds.


title: The Coming Decade of $h \to \tau\mu$ and $\tau \to \mu\gamma$ Interplay in $\tau$ Flavor Violation Search

abstract: We explore the interplay between $h(125) \to \tau\mu$ search at the LHC and $\tau \to \mu\gamma$ at the up and coming Belle~II experiment, in context of the general two Higgs doublet model with extra Yukawa couplings such as $\rho_{\tau\mu}$. The search for $h \to \tau\mu$ constrains $\rho_{\tau\mu} \cos\gamma$, where $\cos\gamma$ is the $h$--$H$ mixing angle of $h$ with the exotic $CP$-even scalar $H$. For $\tau \to \mu\gamma$, we define the "BSM-benchmark" by setting the extra top Yukawa coupling $\rho_{tt} = 1\ (\cong \lambda_t)$ in two-loop diagrams, and $\cos\gamma = 0$ to decouple $h$. We show that this leading effect due to $H$ and $CP$-odd scalar $A$ can be readily probed by Belle~II, even for the conservative value of $\rho_{\tau\mu} = 0.7\lambda_\tau$. We define the subleading "$h$-benchmark" by setting $\rho_{tt} = 0$ in two-loop diagrams, and take the conservative maximal value of $\cos\gamma = 0.2$. We show that it falls beyond Belle~II reach, but can interfere with the BSM-benchmark effect, which in principle probes the phase of $\rho_{tt}$. We further show that the one-loop $H$, $A$ effect, proportional to $\rho_{\tau\mu}\rho_{\tau\tau}$ in amplitude, is beyond the sensitivity of Belle~II to probe, even for $\rho_{\tau\tau}$ as large as $3\lambda_\tau$. With the working assumption that $\rho_{32}^f, \rho_{33}^f = {\cal O}(\lambda_3^f)$ for all charged fermions $f$, we find good discovery potential for both $\tau$ lepton flavor violation searches in the coming decade.


title: $\Pi-$nets: Deep Polynomial Neural Networks

abstract: Deep Convolutional Neural Networks (DCNNs) is currently the method of choice both for generative, as well as for discriminative learning in computer vision and machine learning. The success of DCNNs can be attributed to the careful selection of their building blocks (e.g., residual blocks, rectifiers, sophisticated normalization schemes, to mention but a few). In this paper, we propose $\Pi$-Nets, a new class of DCNNs. $\Pi$-Nets are polynomial neural networks, i.e., the output is a high-order polynomial of the input. $\Pi$-Nets can be implemented using special kind of skip connections and their parameters can be represented via high-order tensors. We empirically demonstrate that $\Pi$-Nets have better representation power than standard DCNNs and they even produce good results without the use of non-linear activation functions in a large battery of tasks and signals, i.e., images, graphs, and audio. When used in conjunction with activation functions, $\Pi$-Nets produce state-of-the-art results in challenging tasks, such as image generation. Lastly, our framework elucidates why recent generative models, such as StyleGAN, improve upon their predecessors, e.g., ProGAN.


title: RP-CARS reveals molecular spatial order anomalies in myelin of an animal model of Krabbe disease

abstract: Krabbe disease (KD) is a rare demyelinating sphingolipidosis, often fatal in the first years of life. It is caused by the inactivation of the galactocerebrosidase (GALC) enzyme that causes an increase in the cellular levels of psychosine considered to be at the origin of the tissue-level effects. GALC is inactivated also in the Twitcher (TWI) mouse: a genetic model of KD that is providing important insights into the understating of the pathogenetic process and the development of possible treatments. In this article an innovative optical technique, RP-CARS, is proposed as a tool to study the degree of order of the CH2 bonds inside the myelin sheaths of TWI-mice sciatic-nerve fibres. RP-CARS, a recently developed variation of CARS microscopy, is able to combine the intrinsic chemical selectivity of CARS microscopy with molecular-bond-spatial-orientation sensibility. This is the first time RP-CARS is applied to the study of a genetic model of a pathology, leading to the demonstration of a post-onset progressive spatial disorganization of the myelin CH2 bonds. The presented result could be of great interest for a deeper understanding of the pathogenic mechanisms underlying the human KD and, moreover, it is an additional proof of the experimental validity of this microscopy technique.


title: The Fast Loaded Dice Roller: A Near-Optimal Exact Sampler for Discrete Probability Distributions

abstract: This paper introduces a new algorithm for the fundamental problem of generating a random integer from a discrete probability distribution using a source of independent and unbiased random coin flips. We prove that this algorithm, which we call the Fast Loaded Dice Roller (FLDR), is highly efficient in both space and time: (i) the size of the sampler is guaranteed to be linear in the number of bits needed to encode the input distribution; and (ii) the expected number of bits of entropy it consumes per sample is at most 6 bits more than the information-theoretically optimal rate. We present fast implementations of the linear-time preprocessing and near-optimal sampling algorithms using unsigned integer arithmetic. Empirical evaluations on a broad set of probability distributions establish that FLDR is 2x-10x faster in both preprocessing and sampling than multiple baseline algorithms, including the widely-used alias and interval samplers. It also uses up to 10000x less space than the information-theoretically optimal sampler, at the expense of less than 1.5x runtime overhead.

categories: stat.CO cs.DM cs.DS cs.IT math.IT math.PR


title: Why You Cannot (Yet) Write an "Interval Arithmetic" Library in Common Lisp

abstract: "Interval Arithmetic" (IA) appears to be a useful numerical tool to have at hand in several applications. Alas, the current IA descriptions and proposed standards are always formulated in terms of the IEEE-754 standard, and the status of IEEE-754 compliance of most Common Lisp implementations is not up to par. A solution would be for Common Lisp implementations to adhere to the Language Independent Arithmetic} (LIA) IEC standard, which includes IEEE/754. While the LIA standard provides a set of proposed bindings for Common Lisp, the format and depth of the specification documents is not readily usable by a Common Lisp programmer, should an implementation decide to comply with the provisions. Moreover, much latitude is left to each implementation to provide the LIA "environmental" setup. It would be most beneficial if more precision were agreed upon by the Common Lisp community about how to provide LIA compliance in the implementations. In that case, a new set of documentation or manuals in the style of the HyperSpec could be provided, for the benefit of the Common Lisp programmer.


title: A Banach space induced by an almost disjoint family, admitting only few operators and decompositions

abstract: We consider the closed subspace of $\ell_\infty$ generated by $c_0$ and the characteristic functions of elements of an uncountable, almost disjoint family $\mathcal A$ of infinite subsets of $\mathbb N$. This Banach space has the form $C_0(K_{\mathcal A})$ for a locally compact Hausdorff space $K_{\mathcal A}$ that is known under many names, such as $\Psi$-space and Isbell--Mr\'owka space. We construct an uncountable, almost disjoint family ${\mathcal A}$ such that the Banach algebra of all bounded linear operators on $C_0(K_{\mathcal A})$ is as small as possible in the sense that every bounded linear operator on $C_0(K_{\mathcal A})$ is the sum of a scalar multiple of the identity and an operator that factors through $c_0$ (which in this case is equivalent to having separable range). This implies that $C_0(K_{\mathcal A})$ has the fewest possible decompositions: whenever $C_0(K_{\mathcal A})=X\oplus Y$ with $dim({X})=\infty$, $dim({Y})=\infty$, either ${X}$ is isomorphic to $C_0(K_{\mathcal A})$ and ${Y}$ to $c_0$, or vice versa. These results improve previous work of the first named author in which an extra set-theoretic hypothesis was required. We also discuss the consequences of these results for the algebra of all bounded linear operators on our Banach space $C_0(K_{\mathcal A})$ concerning the lattice of closed ideals, characters and automatic continuity of homomorphisms. To exploit the perfect set property for Borel sets as in the classical construction of an almost disjoint family of Mr\'owka we need to deal with $\mathbb N \times \mathbb N$-matrices rather than with the usual partitioners. This noncommutative setting requires new ideas inspired by the theory of compact and weakly compact operators and the use of an extraction principle due to F. van Engelen, K. Kunen and A. Miller concerning Borel subsets of the square.


title: The effect of material defects on resonant spin wave modes in a nanomagnet

abstract: We have theoretically studied how resonant spin wave modes in an elliptical nanomagnet are affected by fabrication defects, such as small local thickness variations. Our results indicate that defects of this nature, which can easily result from the fabrication process, or are sometimes deliberately introduced during the fabrication process, will significantly alter the frequencies, magnetic field dependence of the frequencies, and the power and phase profiles of the resonant spin wave modes. They can also spawn new resonant modes and quench existing ones. All this has important ramifications for multi-device circuits based on spin waves, such as phase locked oscillators for neuromorphic computing, where the device-to-device variability caused by defects can be inhibitory.


title: The shape of the value function under Poisson optimal stopping

abstract: In a classical problem for the stopping of a diffusion process $(X_t)_{t \geq 0}$, where the goal is to maximise the expected discounted value of a function of the stopped process ${\mathbb E}^x[e^{-\beta \tau}g(X_\tau)]$, maximisation takes place over all stopping times $\tau$. In a Poisson optimal stopping problem, stopping is restricted to event times of an independent Poisson process. In this article we consider whether the resulting value function $V_\theta(x) = \sup_{\tau \in {\mathcal T}({\mathbb T}^\theta)}{\mathbb E}^x[e^{-\beta \tau}g(X_\tau)]$ (where the supremum is taken over stopping times taking values in the event times of an inhomogeneous Poisson process with rate $\theta = (\theta(X_t))_{t \geq 0}$) inherits monotonicity and convexity properties from $g$. It turns out that monotonicity (respectively convexity) of $V_\theta$ in $x$ depends on the monotonicity (respectively convexity) of the quantity $\frac{\theta(x) g(x)}{\theta(x) + \beta}$ rather than $g$. Our main technique is stochastic coupling.


title: Multivariate Boosted Trees and Applications to Forecasting and Control

abstract: Gradient boosted trees are competition-winning, general-purpose, non-parametric regressors, which exploit sequential model fitting and gradient descent to minimize a specific loss function. The most popular implementations are tailored to univariate regression and classification tasks, precluding the possibility of capturing multivariate target cross-correlations and applying structured penalties to the predictions. In this paper, we present a computationally efficient algorithm for fitting multivariate boosted trees. We show that multivariate trees can outperform their univariate counterpart when the predictions are correlated. Furthermore, the algorithm allows to arbitrarily regularize the predictions, so that properties like smoothness, consistency and functional relations can be enforced. We present applications and numerical results related to forecasting and control.


title: Fine-Grained Visual Classification via Progressive Multi-Granularity Training of Jigsaw Patches

abstract: Fine-grained visual classification (FGVC) is much more challenging than traditional classification tasks due to the inherently subtle intra-class object variations. Recent works mainly tackle this problem by focusing on how to locate the most discriminative parts, more complementary parts, and parts of various granularities. However, less effort has been placed to which granularities are the most discriminative and how to fuse information cross multi-granularity. In this work, we propose a novel framework for fine-grained visual classification to tackle these problems. In particular, we propose: (i) a progressive training strategy that effectively fuses features from different granularities, and (ii) a random jigsaw patch generator that encourages the network to learn features at specific granularities. We obtain state-of-the-art performances on several standard FGVC benchmark datasets, where the proposed method consistently outperforms existing methods or delivers competitive results. The code will be available at https://github.com/PRIS-CV/PMG-Progressive-Multi-Granularity-Training.


title: Hardware Architecture Proposal for TEDA algorithm to Data Streaming Anomaly Detection

abstract: The amount of data in real-time, such as time series and streaming data, available today continues to grow. Being able to analyze this data the moment it arrives can bring an immense added value. However, it also requires a lot of computational effort and new acceleration techniques. As a possible solution to this problem, this paper proposes a hardware architecture for Typicality and Eccentricity Data Analytic (TEDA) algorithm implemented on Field Programmable Gate Arrays (FPGA) for use in data streaming anomaly detection. TEDA is based on a new approach to outlier detection in the data stream context. In order to validate the proposals, results of the occupation and throughput of the proposed hardware are presented. Besides, the bit accurate simulation results are also presented. The project aims to Xilinx Virtex-6 xc6vlx240t-1ff1156 as the target FPGA.


title: Amplitude hysteresis and the synchronization region: Prediction of vortex-induced vibration using a freely forced van der Pol oscillator

abstract: Vortex-induced vibration is a nonlinear phenomenon that can damage buildings or produce energy. Here, the range of synchronization between fluid and structure is an important parameter. Using a coupled van der Pol and linear oscillator system, we derived an analytic, closed-form expression for the synchronization range as a function of fluid shedding frequency, or lock-in as it's also called. Our new definition is an improvement to an earlier synchronization definition as it includes the effect of dissipation; without it hysteresis is missing. The fluid speed range with hysteresis is found to be non-monotonic with forcing strength while the region of absolute synchronization increases monotonically. Our proposed synchronization definition can be a useful tool when designing structural dampers or energy harvesters but is a general solution for all systems that can be modeled with similarly coupled equations.


title: Quantum Einstein equations

abstract: We derive the quantum Einstein equations (which are the quantum generalisation of the Einstein equations of classical gravity) from Bohmian quantum gravity. Bohmian quantum gravity is a non-classical geometrodynamics (in the ADM formalism) which describes the time evolution of a 3-geometry and of a matter field (or other matter degrees of freedom) on a three manifold. The evolution is determined by a velocity law which is defined by the wave function. The wave function itself satisfies the Wheeler-DeWitt equation. We cast the Bohmian dynamics into the form of the Einstein field equations, where the interesting novelty is a contribution to the energy-momentum tensor that depends on the quantum potential.


title: On the geometry of nearly orthogonal lattices

abstract: Nearly orthogonal lattices were formally defined in [4], where their applications to image compression were also discussed. The idea of ``near orthogonality" in $2$-dimensions goes back to the work of Gauss. In this paper, we focus on well-rounded nearly orthogonal lattices in~$\mathbb R^n$ and investigate their geometric and optimization properties. Specifically, we prove that the sphere packing density function on the space of well-rounded lattices in dimension $n\geq 3$ does not have any local maxima on the nearly orthogonal set and has only one local minimum there: at the integer lattice~$\mathbb Z^n$. Further, we show that the nearly orthogonal set cannot contain any perfect lattices for~$n \geq 3$, although it contains multiple eutactic (and even strongly eutactic) lattices in every dimension. This implies that eutactic lattices, while always critical points of the packing density function, are not necessarily local maxima or minima even among the well-rounded lattices. We also prove that a (weakly) nearly orthogonal lattice in~$\mathbb R^n$ contains no more than~$4n-2$ minimal vectors (with any smaller even number possible) and establish some bounds on coherence of these lattices.


title: MAGIC Highlights

abstract: MAGIC has been exploring the sky at Very High Energy gamma-rays (50 GeV - 50 TeV) since 2004, operating first with a single telescope and from 2009 with two telescopes in stereoscopic mode. MAGIC has carried out a observational program involving fundamental physics and astrophysics topics. In this paper we present some of the most important results obtained by MAGIC, with a special emphasis in the most recent ones.


title: Upper bounds for roots of B-functions, following Kashiwara and Lichtin

abstract: By building on a method introduced by Kashiwara and refined by Lichtin, we give upper bounds for the roots of certain b-functions associated to a regular function f in terms of a log resolution of singularities. As applications, we recover with more elementary methods a result of Budur and Saito describing the multiplier ideals of f in terms of the V-filtration of f and a result of the second named author with Popa giving a lower bound for the minimal exponent of f in terms of a log resolution.


title: Software-testing education: A systematic literature mapping

abstract: Context: With the rising complexity and scale of software systems, there is an ever-increasing demand for sophisticated and cost-effective software testing. To meet such a demand, there is a need for a highly-skilled software testing work-force (test engineers) in the industry. To address that need, many university educators worldwide have included software-testing education in their software engineering (SE) or computer science (CS) programs. Objective: Our objective in this paper is to summarize the body of experience and knowledge in the area of software-testing education to benefit the readers (both educators and researchers) in designing and delivering software testing courses in university settings, and to also conduct further education research in this area. Method: To address the above need, we conducted a systematic literature mapping (SLM) to synthesize what the community of educators have published on this topic. After compiling a candidate pool of 307 papers, and applying a set of inclusion/exclusion criteria, our final pool included 204 papers published between 1992 and 2019. Results: The topic of software-testing education is becoming more active, as we can see by the increasing number of papers. Many pedagogical approaches (how to best teach testing), course-ware, and specific tools for testing education have been proposed. Many challenges in testing education and insights on how to overcome those challenges have been proposed. Conclusion: This paper provides educators and researchers with a classification of existing studies within software-testing education. We further synthesize challenges and insights reported when teaching software testing. The paper also provides a reference ("index") to the vast body of knowledge and experience on teaching software testing.


title: Discrete Morse Theory, Persistent Homology and Forman-Ricci Curvature

abstract: Using Banchoff's discrete Morse Theory, in tandem with Bloch's result on the strong connection between the former and Forman's Morse Theory, and our own previous algorithm based on the later, we show that there exists a curvature-based, efficient Persistent Homology scheme for networks and hypernetworks. We also broaden the proposed method to include more general types of networks, by using Bloch's extension of Banchoff's work. Moreover, we show the connection between defect and Forman's Ricci curvature that exists in the combinatorial setting, thus explaining previous empirical results showing very strong correlation between Persistent Homology results obtained using Forman's Morse Theory on the one hand, and Forman's Ricci curvature, on the other.


title: Cross-tier web programming for curated databases: A case study

abstract: Curated databases have become important sources of information across scientific disciplines, and due to the manual work of experts, often become important reference works. Features such as provenance tracking, archiving, and data citation are widely regarded as important features for curated databases, but implementing such features is challenging, and small database projects often lack the resources to do so. A scientific database application is not just the database itself, but also an ecosystem of web applications to display the data, and applications supporting data curation. Supporting advanced curation features requires changing all of these components, and there is currently no way to provide such capabilities in a reusable way. Cross-tier programming languages have been proposed to simplify the creation of web applications, where developers write an application in a single, uniform language. Consequently, database queries and updates can be written in the same language as the rest of the program, and at least in principle, it should be possible to provide curation features reusably via program transformations. As a first step, it is important to establish that realistic curated databases can be implemented in a cross-tier programming language. In this paper, we describe such a case study: reimplementing the web frontend of a real-world scientific database, the IUPHAR/BPS Guide to Pharmacology (GtoPdb), in the Links programming language. We show how features such as language-integrated query simplify the development process, and rule out common errors. We show that the Links implementation performs fewer database queries, while the time needed to handle the queries is comparable to the Java version. While there is some overhead to using Links because of its comparative immaturity compared to Java, the Links version is viable as a proof-of-concept case study.


title: Normal density and moment of inertia of a moving superfluid

abstract: In this work, the normal density $\rho_n$ and moment of inertia of a moving superfluid are investigated. We find that, even at zero temperature, there exists a finite normal density for the moving superfluid. When the velocity of superfluid reaches sound velocity, the normal density becomes total mass density $\rho$, which indicates that the system losses superfluidity. At the same time, the Landau's critical velocity also becomes zero. The existence of the non-zero normal density is attributed to the coupling between the motion of superflow and density fluctuation in transverse directions. With Josephson relation, the superfluid density $\rho_s$ is also calculated and the identity $\rho_s+\rho_n=\rho$ holds. Further more, we find that the finite normal density also results in a quantized moment of inertia in a moving superfluid trapped by a ring. The normal density and moment of inertia at zero temperature could be verified experimentally by measuring the angular momentum of a moving superfluid in a ring trap.


title: Nonlinear least-squares spline fitting with variable knots

abstract: In this paper, we present a nonlinear least-squares fitting algorithm using B-splines with free knots. Since its performance strongly depends on the initial estimation of the free parameters (i.e. the knots), we also propose a fast and efficient knot-prediction algorithm that utilizes numerical properties of first-order B-splines. Using $\ell_p\;(p=1,2,\infty)$ norm solutions, we also provide three different strategies for properly selecting the free knots. Our initial predictions are then iteratively refined by means of a gradient-based variable projection optimization. Our method is general in nature and can be used to estimate the optimal number of knots in cases in which no a-priori information is available. To evaluate the performance of our method, we approximated a one-dimensional discrete time series and conducted an extensive comparative study using both synthetic and real-world data. We chose the problem of electrocardiogram (ECG) signal compression as a real-world case study. Our experiments on the well-known PhysioNet MIT-BIH Arrhythmia database show that the proposed method outperforms other knot-prediction techniques in terms of accuracy while requiring much lower computational complexity.


title: Coronavirus Perceptions And Economic Anxiety

abstract: We provide one of the first systematic assessments of the development and determinants of economic anxiety at the onset of the coronavirus pandemic. Using a global dataset on internet searches and two representative surveys from the US, we document a substantial increase in economic anxiety during and after the arrival of the coronavirus. We also document a large dispersion in beliefs about the pandemic risk factors of the coronavirus, and demonstrate that these beliefs causally affect individuals' economic anxieties. Finally, we show that individuals' mental models of infectious disease spread understate non-linear growth and shape the extent of economic anxiety.


title: Active Fine-Tuning from gMAD Examples Improves Blind Image Quality Assessment

abstract: The research in image quality assessment (IQA) has a long history, and significant progress has been made by leveraging recent advances in deep neural networks (DNNs). Despite high correlation numbers on existing IQA datasets, DNN-based models may be easily falsified in the group maximum differentiation (gMAD) competition with strong counterexamples being identified. Here we show that gMAD examples can be used to improve blind IQA (BIQA) methods. Specifically, we first pre-train a DNN-based BIQA model using multiple noisy annotators, and fine-tune it on multiple subject-rated databases of synthetically distorted images, resulting in a top-performing baseline model. We then seek pairs of images by comparing the baseline model with a set of full-reference IQA methods in gMAD. The resulting gMAD examples are most likely to reveal the relative weaknesses of the baseline, and suggest potential ways for refinement. We query ground truth quality annotations for the selected images in a well controlled laboratory environment, and further fine-tune the baseline on the combination of human-rated images from gMAD and existing databases. This process may be iterated, enabling active and progressive fine-tuning from gMAD examples for BIQA. We demonstrate the feasibility of our active learning scheme on a large-scale unlabeled image set, and show that the fine-tuned method achieves improved generalizability in gMAD, without destroying performance on previously trained databases.


title: A Compiler Assisted Scheduler for Detecting and Mitigating Cache-Based Side Channel Attacks

abstract: Side channel attacks steal secret keys by cleverly leveraging information leakages and can, therefore, break encryption. Thus, detection and mitigation of side channel attacks is a very important problem, but the solutions proposed in the literature have limitations in that they do not work in a real-world multi-tenancy setting on servers, have high false positives, or have high overheads. In this work, we demonstrate a compiler guided scheduler, Biscuit, that detects cache-based side channel attacks for processes scheduled on multi-tenancy server farms. A key element of this solution involves the use of a cache-miss model which is inserted by the compiler at the entrances of loop nests to predict the cache misses of the corresponding loop. Such inserted library calls, or beacons, convey the cache miss information to the scheduler at run time, which uses it to co-schedule processes such that their combined cache footprint does not exceed the maximum capacity of the last level cache. The scheduled processes are then monitored for actual vs predicted cache misses, and when an anomaly is detected, the scheduler performs a search to isolate the attacker. We show that Biscuit is able to detect and mitigate Prime+Probe, Flush+Reload, and Flush+Flush attacks on OpenSSL cryptography algorithms with an F-score of 1, and also to detect and mitigate degradation of service on a vision application suite with an F-score of 0.9375. Under a no-attack scenario, the scheme poses low overheads (up to a maximum of 6 percent). In the case of an attack, the scheme ends up with less than 11 percent overhead and is able to reduce the degradation of service in some cases by 40 percent. With these many desirable features such as an ability to deal with multi-tenancy, its ability to detect attacks early, its ability to mitigate those attacks, and low runtime overheads, Biscuit is a practical solution.


title: On solutions of a partial integro-differential equation in Bessel potential spaces with applications in option pricing models

abstract: In this paper we focus on qualitative properties of solutions to a nonlocal nonlinear partial integro-differential equation (PIDE). Using the theory of abstract semilinear parabolic equations we prove existence and uniqueness of a solution in the scale of Bessel potential spaces. Our aim is to generalize known existence results for a wide class of L\'evy measures including with a strong singular kernel. As an application we consider a class of PIDEs arising in the financial mathematics. The classical linear Black-Scholes model relies on several restrictive assumptions such as liquidity and completeness of the market. Relaxing the complete market hypothesis and assuming a Levy stochastic process dynamics for the underlying stock price process we obtain a model for pricing options by means of a PIDE. We investigate a model for pricing call and put options on underlying assets following a Levy stochastic process with jumps. We prove existence and uniqueness of solutions to the penalized PIDE representing approximation of the linear complementarity problem arising in pricing American style of options under Levy stochastic processes. We also present numerical results and comparison of option prices for various Levy stochastic processes modelling underlying asset dynamics.


title: Low Precision Floating-point Arithmetic for High Performance FPGA-based CNN Acceleration

abstract: Low precision data representation is important to reduce storage size and memory access for convolutional neural networks (CNNs). Yet, existing methods have two major limitations: (1) requiring re-training to maintain accuracy for deep CNNs, and (2) needing 16-bit floating-point or 8-bit fixed-point for a good accuracy. In this paper, we propose a low precision (8-bit) floating-point (LPFP) quantization method for FPGA-based acceleration to overcome the above limitations. Without any re-training, LPFP finds an optimal 8-bit data representation with negligible top-1/top-5 accuracy loss (within 0.5%/0.3% in our experiments, respectively, and significantly better than existing methods for deep CNNs). Furthermore, we implement one 8-bit LPFP multiplication by one 4-bit multiply-adder (MAC) and one 3-bit adder, and therefore implement four 8-bit LPFP multiplications using one DSP slice of Xilinx Kintex 7 family (KC705 in this paper) while one DSP can implement only two 8-bit fixed-point multiplications. Experiments on six typical CNNs for inference show that on average, we improve throughput by 64.5x over Intel i9 CPU and by 1.5x over existing FPGA accelerators. Particularly for VGG16 and YOLO, compared to six recent FPGA accelerators, we improve average throughput by 3.5x and 27.5x and improve average throughput per DSP by 4.1x and 5x, respectively. To the best of our knowledge, this is the first in-depth study to simplify one multiplication for CNN inference to one 4-bit MAC and implement four multiplications within one DSP while maintaining comparable accuracy without any re-training.


title: H-infinity Strong Stabilization via HIFOO, a Package for Fixed-Order Controller Design

abstract: We report on our experience with strong stabilization using HIFOO, a toolbox for H-infinity fixed-order controller design. We applied HIFOO to 21 fixed-order stable H-infinity controller design problems in the literature, comparing the results with those published for other methods. The results show that HIFOO often achieves good H-infinity performance with low-order stable controllers, unlike other methods in the literature.


title: Twisted submanifolds of R^n

abstract: We propose a general procedure to construct noncommutative deformations of an embedded submanifold $M$ of $\mathbb{R}^n$ determined by a set of smooth equations $f^a(x)=0$. We use the framework of Drinfel'd twist deformation of differential geometry of [Aschieri et al., Class. Quantum Gravity 23 (2006), 1883]; the commutative pointwise product is replaced by a (generally noncommutative) $\star$-product determined by a Drinfel'd twist. The twists we employ are based on the Lie algebra $\Xi_t$ of vector fields that are tangent to all the submanifolds that are level sets of the $f^a$; the twisted Cartan calculus is automatically equivariant under twisted tangent infinitesimal diffeomorphisms. We can consistently project a connection from the twisted $\mathbb{R}^n$ to the twisted $M$ if the twist is based on a suitable Lie subalgebra $\mathfrak{e}\subset\Xi_t$. If we endow $\mathbb{R}^n$ with a metric then twisting and projecting to the normal and tangent vector fields commute, and we can project the Levi-Civita connection consistently to the twisted $M$, provided the twist is based on the Lie subalgebra $\mathfrak{k}\subset\mathfrak{e}$ of the Killing vector fields of the metric; a twisted Gauss theorem follows, in particular. Twisted algebraic manifolds can be characterized in terms of generators and polynomial relations. We present in some detail twisted cylinders embedded in twisted Euclidean $\mathbb{R}^3$ and twisted hyperboloids embedded in twisted Minkowski $\mathbb{R}^3$ [these are twisted (anti-)de Sitter spaces $dS_2,AdS_2$].

categories: math-ph gr-qc hep-th math.DG math.MP math.QA


title: Atomic physics studies at the Gamma Factory at CERN

abstract: The Gamma Factory initiative proposes to develop novel research tools at CERN by producing, accelerating and storing highly relativistic, partially stripped ion beams in the SPS and LHC storage rings. By exciting the electronic degrees of freedom of the stored ions with lasers, high-energy narrow-band photon beams will be produced by properly collimating the secondary radiation that is peaked in the direction of ions' propagation. Their intensities, up to $10^{17}$ photons per second, will be several orders of magnitude higher than those of the presently operating light sources in the particularly interesting $\gamma$--ray energy domain reaching up to 400 MeV. This article reviews opportunities that may be afforded by utilizing the primary beams for spectroscopy of partially stripped ions circulating in the storage ring, as well as the atomic-physics opportunities afforded by the use of the secondary high-energy photon beams. The Gamma Factory will enable ground breaking experiments in spectroscopy and novel ways of testing fundamental symmetries of nature.

categories: physics.atom-ph hep-ex hep-ph physics.acc-ph


title: A Tracking System For Baseball Game Reconstruction

abstract: The baseball game is often seen as many contests that are performed between individuals. The duel between the pitcher and the batter, for example, is considered the engine that drives the sport. The pitchers use a variety of strategies to gain competitive advantage against the batter, who does his best to figure out the ball trajectory and react in time for a hit. In this work, we propose a system that captures the movements of the pitcher, the batter, and the ball in a high level of detail, and discuss several ways how this information may be processed to compute interesting statistics. We demonstrate on a large database of videos that our methods achieve comparable results as previous systems, while operating solely on video material. In addition, state-of-the-art AI techniques are incorporated to augment the amount of information that is made available for players, coaches, teams, and fans.


title: Limiting distributions for eigenvalues of sample correlation matrices from heavy-tailed populations

abstract: Consider a $p$-dimensional population ${\mathbf x} \in\mathbb{R}^p$ with iid coordinates in the domain of attraction of a stable distribution with index $\alpha\in (0,2)$. Since the variance of ${\mathbf x}$ is infinite, the sample covariance matrix ${\mathbf S}_n=n^{-1}\sum_{i=1}^n {{\mathbf x}_i}{\mathbf x}'_i$ based on a sample ${\mathbf x}_1,\ldots,{\mathbf x}_n$ from the population is not well behaved and it is of interest to use instead the sample correlation matrix ${\mathbf R}_n= \{\operatorname{diag}({\mathbf S}_n)\}^{-1/2}\, {\mathbf S}_n \{\operatorname{diag}({\mathbf S}_n)\}^{-1/2}$. This paper finds the limiting distributions of the eigenvalues of ${\mathbf R}_n$ when both the dimension $p$ and the sample size $n$ grow to infinity such that $p/n\to \gamma \in (0,\infty)$. The family of limiting distributions $\{H_{\alpha,\gamma}\}$ is new and depends on the two parameters $\alpha$ and $\gamma$. The moments of $H_{\alpha,\gamma}$ are fully identified as sum of two contributions: the first from the classical Mar\v{c}enko-Pastur law and a second due to heavy tails. Moreover, the family $\{H_{\alpha,\gamma}\}$ has continuous extensions at the boundaries $\alpha=2$ and $\alpha=0$ leading to the Mar\v{c}enko-Pastur law and a modified Poisson distribution, respectively. Our proofs use the method of moments, the path-shortening algorithm developed in [18] and some novel graph counting combinatorics. As a consequence, the moments of $H_{\alpha,\gamma}$ are expressed in terms of combinatorial objects such as Stirling numbers of the second kind. A simulation study on these limiting distributions $H_{\alpha,\gamma}$ is also provided for comparison with the Mar\v{c}enko-Pastur law.


title: K-theory for semigroup C*-algebras and partial crossed products

abstract: Using the Baum-Connes conjecture with coefficients, we develop a K-theory formula for reduced C*-algebras of strongly $0$-$E$-unitary inverse semigroups, or equivalently, for certain reduced partial crossed products. In the case of semigroup C*-algebras, we obtain a generalization of previous K-theory results of Cuntz, Echterhoff and the author without having to assume the Toeplitz condition. As applications, we discuss semigroup C*-algebras of Artin monoids, Baumslag-Solitar monoids, one-relator monoids, C*-algebras generated by right regular representations of semigroups from number theory, and C*-algebras of inverse semigroups arising in the context of tilings.


title: Semi-device-independent certification of independent quantum state and measurement devices

abstract: Certifying that quantum devices behave as intended is crucial for quantum information science. Here, methods are developed for certification of both state preparation devices and measurement devices based on prepare-and-measure experiments with independent devices. The experimenter assumes the independence of the devices and knowledge of the Hilbert space dimension. Thus no precise characterisation of any part of the experiment is required. The certification is based on a randomised version of unambiguous state discrimination and targets the class of state ensembles corresponding to quantum $t$-designs of any size and any dimension. These quantum designs are sets of states over which the average of any $t$-degree polynomial equals its average over all pure states, and they accommodate many of the most useful discrete structures in quantum information processing. Furthermore, it is shown that the same experiments also certify the detection efficiency of the measurement devices, as well as their non-projective nature. The presented methods can readily be implemented in experiments.


title: A Modular Small-Signal Analysis Framework for Inverter Penetrated Power Grids: Measurement, Assembling, Aggregation, and Stability Assessment

abstract: Unprecedented dynamic phenomena may appear in power grids due to higher and higher penetration of inverter-based resources (IBR), e.g., wind and solar photovoltaic (PV). A major challenge in dynamic modeling and analysis is that unlike synchronous generators, whose analytical models are well studied and known to system planners, inverter models are proprietary information with black box models provided to utilities. Thus, measurement based characterization of IBR is a popular approach to find frequency-domain response of an IBR. The resulting admittances are essentially small-signal current/voltage relationship in frequency domain. Integrating admittances for grid dynamic modeling and analysis requires a new framework, namely modular small-signal analysis framework. In this visionary paper, we examine the current state-of-the-art of dynamic modeling and analysis of power grids with IBR, including inverter admittance characterization, the procedure of component assembling and aggregation, and stability assessment. We push forward a computing efficient modular modeling and analysis framework via four visions: (i) efficient and accurate admittance model characterization via model building and time-domain responses, (ii) accurate assembling of components, (iii) efficient aggregation, and (iv) stability assessment relying on network admittance matrices. Challenges of admittance-based modular analysis are demonstrated using examples and techniques to tackle those challenges are pointed out in this visionary paper.


title: Binomial ideals and congruences on $\mathbb{N}^n$

abstract: A \emph{congruence} on $\mathbb{N}^n$ is an equivalence relation on $\mathbb{N}^n$ that is compatible with the additive structure. If $\Bbbk$ is a field, and $I$ is a \emph{binomial ideal} in $\Bbbk[X_1,\dots,X_n]$ (that is, an ideal generated by polynomials with at most two terms), then $I$ induces a congruence on $\mathbb{N}^n$ by declaring $\mathbf{u}$ and $\mathbf{v}$ to be equivalent if there is a linear combination with nonzero coefficients of $\mathbf{X}^{\mathbf{u}}$ and $\mathbf{X}^{\mathbf{v}}$ that belongs to $I$. While every congruence on $\mathbb{N}^n$ arises this way, this is not a one-to-one correspondence, as many binomial ideals may induce the same congruence. Nevertheless, the link between a binomial ideal and its corresponding congruence is strong, and one may think of congruences as the underlying combinatorial structures of binomial ideals. In the current literature, the theories of binomial ideals and congruences on $\mathbb{N}^n$ are developed separately. The aim of this survey paper is to provide a detailed parallel exposition, that provides algebraic intuition for the combinatorial analysis of congruences. For the elaboration of this survey paper, we followed mainly [Kahle and Miller, Algebra Number Theory 8(6):1297-1364, 2014] with an eye on [Eisenbud and Sturmfels. Duke Math J 84(1):1-45, 1996] and [Ojeda and Piedra S\'anchez, J. Symbolic Comput 30(4):383-400, 2000].


title: Improving Training on Noisy Stuctured Labels

abstract: Fine-grained annotations---e.g. dense image labels, image segmentation and text tagging---are useful in many ML applications but they are labor-intensive to generate. Moreover there are often systematic, structured errors in these fine-grained annotations. For example, a car might be entirely unannotated in the image, or the boundary between a car and street might only be coarsely annotated. Standard ML training on data with such structured errors produces models with biases and poor performance. In this work, we propose a novel framework of Error-Correcting Networks (ECN) to address the challenge of learning in the presence structured error in fine-grained annotations. Given a large noisy dataset with commonly occurring structured errors, and a much smaller dataset with more accurate annotations, ECN is able to substantially improve the prediction of fine-grained annotations compared to standard approaches for training on noisy data. It does so by learning to leverage the structures in the annotations and in the noisy labels. Systematic experiments on image segmentation and text tagging demonstrate the strong performance of ECN in improving training on noisy structured labels.


title: The rencontre problem

abstract: Let $\left\{X^{1}_k\right\}_{k=1}^{\infty}, \left\{X^{2}_k\right\}_{k=1}^{\infty}, \cdots, \left\{X^{d}_k\right\}_{k=1}^{\infty}$ be $d$ independent sequences of Bernoulli random variables with success-parameters $p_1, p_2, \cdots, p_d$ respectively, where $d \geq 2$ is a positive integer, and $ 0<p_j<1$ for all $j=1,2,\cdots,d.$ Let \begin{equation*} S^{j}(n) = \sum_{i=1}^{n} X^{j}_{i} = X^{j}_{1} + X^{j}_{2} + \cdots + X^{j}_{n}, \quad n =1,2 , \cdots. \end{equation*} We declare a "rencontre" at time $n$, or, equivalently, say that $n$ is a "rencontre-time," if \begin{equation*} S^{1}(n) = S^{2}(n) = \cdots = S^{d}(n). \end{equation*} We motivate and study the distribution of the first (provided it is finite) rencontre time.


title: Low-noise octave-spanning mid-infrared supercontinuum generation in a multimode chalcogenide fiber

abstract: We demonstrate the generation of a low-noise, octave-spanning mid-infrared supercontinuum from 1700 to 4800 nm by injecting femtosecond pulses into the normal dispersion regime of a multimode step-index chalcogenide fiber with 100 $\mu$m core diameter. We conduct a systematic study of the intensity noise across the supercontinuum spectrum and show that the initial fluctuations of the pump laser are at most amplified by a factor of three. We also perform a comparison with the noise characteristics of an octave-spanning supercontinuum generated in the anomalous dispersion regime of a multimode fluoride fiber with similar core size and show that the all-normal dispersion supercontinuum in the multimode chalcogenide fiber has superior noise characteristics. Our results open up novel perspective for many practical applications such as long-distance remote sensing where high power and low noise are paramount.


title: Response solutions for strongly dissipative quasi-periodically forced systems with arbitrary nonlinearities and frequencies

abstract: We consider quasi-periodically systems in the presence of dissipation and study the existence of response solutions, i.e. quasi-periodic solutions with the same frequency vector as the forcing term. When the dissipation is large enough and a suitable function involving the forcing has a simple zero, response solutions are known to exist without assuming any non resonance condition on the frequency vector. We analyse the case of non-simple zeroes and, in order to deal with the small divisors problem, we confine ourselves to two-dimensional frequency vectors, so as to use the properties of continued fractions. We show that, if the order of the zero is odd (if it is even, in general no response solution exists), a response solution still exists provided the inverse of the parameter measuring the dissipation belongs to a set given by the union of infinite intervals depending on the convergents of the ratio of the two components of the frequency vector. The intervals may be disjoint and as a consequence we obtain the existence of response solutions in a set with "holes". If we want the set to be connected we have to require some non-resonance condition on the frequency: in fact, we need a condition weaker than the Bryuno condition usually considered in small divisors problems.


title: Online Data-Enabled Predictive Control

abstract: We develop an online data-enabled predictive (ODeePC) control method for optimal control of unknown systems, building on the recently proposed DeePC [1]. Our proposed ODeePC method leverages a primal-dual algorithm with real-time measurement feedback to iteratively compute the corresponding real-time optimal control policy as system conditions change. The proposed ODeePC conceptual-wise resembles standard adaptive system identification and model predictive control (MPC), but it provides a new alternative for the standard methods. ODeePC is enabled by computationally efficient methods that exploit the special structure of the Hankel matrices in the context of DeePC with Fast Fourier Transform (FFT) and primal-dual algorithm. We provide theoretical guarantees regarding the asymptotic behavior of ODeePC, and we demonstrate its performance through numerical examples.


title: Strategic Abilities of Asynchronous Agents: Semantic Side Effects and How to Tame Them

abstract: Recently, we have proposed a framework for verification of agents' abilities in asynchronous multi-agent systems, together with an algorithm for automated reduction of models. The semantics was built on the modeling tradition of distributed systems. As we show here, this can sometimes lead to counterintuitive interpretation of formulas when reasoning about the outcome of strategies. First, the semantics disregards finite paths, and thus yields unnatural evaluation of strategies with deadlocks. Secondly, the semantic representations do not allow to capture the asymmetry between proactive agents and the recipients of their choices. We propose how to avoid the problems by a suitable extension of the representations and change of the execution semantics for asynchronous MAS. We also prove that the model reduction scheme still works in the modified framework.


title: CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations

abstract: CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular and biological systems. It is especially aimed at massively-parallel and linear-scaling electronic structure methods and state-of-the-art ab-initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2k to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension.


title: Optical Thermometry with Quantum Emitters in Hexagonal Boron Nitride

abstract: Nanoscale optical thermometry is a promising non-contact route for measuring local temperature with both high sensitivity and spatial resolution. In this work, we present a deterministic optical thermometry technique based on quantum emitters in nanoscale hexagonal boron-nitride. We show that these nanothermometers exhibit better performance than that of homologous, all-optical nanothermometers both in sensitivity and range of working temperature. We demonstrate their effectiveness as nanothermometers by monitoring the local temperature at specific locations in a variety of custom-built micro-circuits. This work opens new avenues for nanoscale temperature measurements and heat flow studies in miniaturized, integrated devices.


title: 3-Symmetric Graphs

abstract: An intuitive property of a random graph is that its subgraphs should also appear randomly distributed. We consider graphs whose subgraph densities exactly match their expected values. We call graphs with this property for all subgraphs with $k$ vertices to be $k$-symmetric. We discuss some properties and examples of such graphs. We construct 3-symmetric graphs and provide some statistics.


title: Accurate and efficient description of interacting carriers in quantum nanostructures by selected configuration interaction and perturbation theory

abstract: We present a method to calculate many-body states of interacting carriers in million atom quantum nanostructures based on atomistic tight-binding calculations and a combination of iterative selection of configurations and perturbation theory. This method enables investigations of large excitonic complexes and multi-electron systems with near full configuration interaction accuracy, even though only a small subspace of the full many-body Hilbert space is sampled, thus saving orders of magnitudes in computational resources. Important advantages of this method are that the convergence is controlled by a single parameter, the threshold, and that ground and excited states can be treated on an equal footing. We demonstrate the extreme efficiency of the method by numerical studies of complexes composed of up to 13 excitons, which requires filling of states up to the fourth electronic shell. We find that the method generally converges fast as a function of the threshold, profiting from a significant enhancement due to the perturbative corrections. The role of the choice of single-particle basis states is discussed. It is found that the algorithm converges faster in the Hartree-Fock basis only for highly charged systems, where Coulomb repulsion dominates. Finally, based on the observation that second order perturbative energy corrections only depend on off-diagonal elements of the many-body Hamiltonian, we present a way to accurately calculate many-body states that requires only a relatively small number of Coulomb matrix elements.


title: Graph Clustering Via QUBO and Digital Annealing

abstract: This article empirically examines the computational cost of solving a known hard problem, graph clustering, using novel purpose-built computer hardware. We express the graph clustering problem as an intra-cluster distance or dissimilarity minimization problem. We formulate our poblem as a quadratic unconstrained binary optimization problem and employ a novel computer architecture to obtain a numerical solution. Our starting point is a clustering formulation from the literature. This formulation is then converted to a quadratic unconstrained binary optimization formulation. Finally, we use a novel purpose-built computer architecture to obtain numerical solutions. For benchmarking purposes, we also compare computational performances to those obtained using a commercial solver, Gurobi, running on conventional hardware. Our initial results indicate the purpose-built hardware provides equivalent solutions to the commercial solver, but in a very small fraction of the time required.


title: Minimum Wiener Index of Triangulations and Quadrangulations

abstract: The Wiener index of a connected graph is the sum of the distances between all unordered pairs of vertices. We provide formulae for the minimum Wiener index of simple triangulations and quadrangulations with connectivity at least $c$, and provide the extremal structures, which attain those values. Our main tool is setting upper bounds for the maximum degree in highly connected triangulations and quadrangulations.


title: Dynamical, value-based decision making among $N$ options

abstract: Decision making is a fundamental capability of autonomous systems. As decision making is a process which happens over time, it can be well modeled by dynamical systems. Often, decisions are made on the basis of perceived values of the underlying options and the desired outcome is to select the option with the highest value. This can be encoded as a bifurcation which produces a stable equilibrium corresponding to the high-value option. When some options have identical values, it is natural to design the decision-making model to be indifferent among the equally-valued options, leading to symmetries in the underlying dynamical system. For example, when all $N$ options have identical values, the dynamical system should have $S_N$ symmetry. Unfortunately, constructing a dynamical system that unfolds the $S_N$-symmetric pitchfork bifurcation is non-trivial. In this paper, we develop a method to construct an unfolding of the pitchfork bifurcation with a symmetry group that is a significant subgroup of $S_N$. The construction begins by parsing the decision among $N$ options into a hierarchical set of $N-1$ binary decisions encoded in a binary tree. By associating the unfolding of a standard $S_2$-symmetric pitchfork bifurcation with each of these binary decisions, we develop an unfolding of the pitchfork bifurcation with symmetries corresponding to isomorphisms of the underlying binary tree.


title: Overview of the CCKS 2019 Knowledge Graph Evaluation Track: Entity, Relation, Event and QA

abstract: Knowledge graph models world knowledge as concepts, entities, and the relationships between them, which has been widely used in many real-world tasks. CCKS 2019 held an evaluation track with 6 tasks and attracted more than 1,600 teams. In this paper, we give an overview of the knowledge graph evaluation tract at CCKS 2019. By reviewing the task definition, successful methods, useful resources, good strategies and research challenges associated with each task in CCKS 2019, this paper can provide a helpful reference for developing knowledge graph applications and conducting future knowledge graph researches.


title: How much is your Strangle worth? On the relative value of the $\delta-$Symmetric Strangle under the Black-Scholes model

abstract: Trading option strangles is a highly popular strategy often used by market participants to mitigate volatility risks in their portfolios. In this paper we propose a measure of the relative value of a delta-Symmetric Strangle and compute it under the standard Black-Scholes option pricing model. This new measure accounts for the price of the strangle, relative to the Present Value of the spread between the two strikes, all expressed, after a natural re-parameterization, in terms of delta and a volatility parameter. We show that under the standard BS option pricing model, this measure of relative value is bounded by a simple function of delta only and is independent of the time to expiry, the price of the underlying security or the prevailing volatility used in the pricing model. We demonstrate how this bound can be used as a quick {\it benchmark} to assess, regardless the market volatility, the duration of the contract or the price of the underlying security, the market (relative) value of the $\delta-$strangle in comparison to its BS (relative) price. In fact, the explicit and simple expression for this measure and bound allows us to also study in detail the strangle's exit strategy and the corresponding {\it optimal} choice for a value of delta.


title: FoCL: Feature-Oriented Continual Learning for Generative Models

abstract: In this paper, we propose a general framework in continual learning for generative models: Feature-oriented Continual Learning (FoCL). Unlike previous works that aim to solve the catastrophic forgetting problem by introducing regularization in the parameter space or image space, FoCL imposes regularization in the feature space. We show in our experiments that FoCL has faster adaptation to distributional changes in sequentially arriving tasks, and achieves the state-of-the-art performance for generative models in task incremental learning. We discuss choices of combined regularization spaces towards different use case scenarios for boosted performance, e.g., tasks that have high variability in the background. Finally, we introduce a forgetfulness measure that fairly evaluates the degree to which a model suffers from forgetting. Interestingly, the analysis of our proposed forgetfulness score also implies that FoCL tends to have a mitigated forgetting for future tasks.


title: Distinguishing between Flaring and Non-Flaring Active Regions: A Machine Learning Perspective

abstract: Large scale solar eruptions significantly impact space weather and damages space-based human infrastructures. It is necessary to predict large scale solar eruptions, which will enable us to protect our vulnerable infrastructures of modern society. We aim to investigate the difference between flaring and non-flaring active regions. We use photospheric vector magnetogram data from Solar Dynamic Observatory's Helioseismic Magnetic Imager to study the time evolution of photospheric magnetic parameters on the solar surface. We build a database of flaring and non-flaring active region observed on the solar surface from the years 2010 to 2017. We train the machine learning algorithm by the time evolution of these active region parameters. Finally, we estimate the performance obtained from this machine learning algorithm. We find the strength of some magnetic parameters namely total unsigned magnetic flux, total unsigned magnetic helicity, total unsigned vertical current and total photospheric magnetic energy density in flaring active regions are much higher compared to the non-flaring ones. These magnetic parameters in the flaring active region are highly evolving and complex. We are able to obtain good forecasting capability with a relatively high value of true skill statistic (TSS). We also find that time evolution of total unsigned magnetic helicity and total unsigned magnetic flux have very high ability to distinguish flaring and non-flaring active regions. It is possible to distinguish flaring active region from the non-flaring one with good accuracy. We confirm that there is no single common parameter which can distinguish all flaring active regions from the non flaring one. However, time evolution of top few magnetic parameters namely total unsigned magnetic flux and total unsigned magnetic helicity have very high distinguishing capability.