title: Improved Surrogate Modeling using Machine Learning for Industrial Civil Aircraft Aerodynamics

abstract: Predicting and simulating aerodynamic fields for civil aircraft over wide flight envelopes represent a real challenge mainly due to significant numerical costs and complex flows. Surrogate models and reduced-order models help to estimate aerodynamic fields from a few well-selected simulations. However, their accuracy dramatically decreases when different physical regimes are involved. Therefore, a method of local non-intrusive reduced-order models using machine learning, called Local Decomposition Method, has been developed to mitigate this issue. This paper introduces several enhancements to this method and presents a complex application to an industrial-like three-dimensional aircraft configuration over a full flight envelope. The enhancements of the method cover several aspects: choosing the best number of models, estimating apriori errors, improving the adaptive sampling for parallel issues, and better handling the borders between local models. The application is supported by an analysis of the model behavior, with a focus on the machine learning methods and the local properties. The model achieves strong levels of accuracy, in particular with two sub-models: one for the subsonic regime and one for the transonic regime. These results highlight that local models and machine learning represent very promising solutions to deal with surrogate models for aerodynamics.


title: Oscillatory asymptotics for the Airy kernel determinant on two intervals

abstract: We obtain asymptotics for the Airy kernel Fredholm determinant on two intervals. We give explicit formulas for all the terms up to and including the oscillations of order $1$, which are expressed in terms of Jacobi $\theta$-functions.


title: Distribution of rare saddles in the $p$-spin energy landscape

abstract: We compute the statistical distribution of index-1 saddles surrounding a given local minimum of the $p$-spin energy landscape, as a function of their distance to the minimum in configuration space and of the energy of the latter. We identify the saddles also in the region of configuration space in which they are subdominant in number (i.e., rare) with respect to local minima, by computing large deviation probabilities of the extremal eigenvalues of their Hessian. As an independent result, we determine the joint large deviation probability of the smallest eigenvalue and eigenvector of a GOE matrix perturbed with both an additive and multiplicative finite-rank perturbation.


title: Measurement Scheduling for Cooperative Localization in Resource-Constrained Conditions

abstract: This paper studies the measurement scheduling problem for a group of N mobile robots moving on a flat surface that are preforming cooperative localization (CL). We consider a scenario in which due to the limited on-board resources such as battery life and communication bandwidth only a given number of relative measurements per robot are allowed at observation and update stage. Optimal selection of which teammates a robot should take a relative measurement from such that the updated joint localization uncertainty of the team is minimized is an NP-hard problem. In this paper, we propose a suboptimal greedy approach that allows each robot to choose its landmark robots locally in polynomial time. Our method, unlike the known results in the literature, does not assume full-observability of CL algorithm. Moreover, it does not require inter-robot communication at scheduling stage. That is, there is no need for the robots to collaborate to carry out the landmark robot selections. We discuss the application of our method in the context of an state-of-the-art decentralized CL algorithm and demonstrate its effectiveness through numerical simulations. Even though our solution does not come with rigorous performance guarantees, its low computational cost along with no communication requirement makes it an appealing solution for operatins with resource constrained robots.


title: Unitarity Flow in 2+1 Dimensional Massive Gravity

abstract: We analyze the most general case of third-order Chern-Simons-like theories of massive 3D gravity. Results show the conditions for finding the unitary regions on the parameter space. There exists $(n-1)$th order theories on the boundary of a unitary $n$th order model on parameter space under certain conditions, as the first example recently demonstrated from the bimetric generalization of exotic massive gravity. We investigated the mechanism that causes this type of transition for third-order models. Hamiltonian analysis of the theory also presents that ghost and no-ghost regions can be separated by Chern-Simons theories.


title: End-to-end facial and physiological model for Affective Computing and applications

abstract: In recent years, Affective Computing and its applications have become a fast-growing research topic. Furthermore, the rise of Deep Learning has introduced significant improvements in the emotion recognition system compared to classical methods. In this work, we propose a multi-modal emotion recognition model based on deep learning techniques using the combination of peripheral physiological signals and facial expressions. Moreover, we present an improvement to proposed models by introducing latent features extracted from our internal Bio Auto-Encoder (BAE). Both models are trained and evaluated on AMIGOS datasets reporting valence, arousal, and emotion state classification. Finally, to demonstrate a possible medical application in affective computing using deep learning techniques, we applied the proposed method to the assessment of anxiety therapy. To this purpose, a reduced multi-modal database has been collected by recording facial expressions and peripheral signals such as Electrocardiogram (ECG) and Galvanic Skin Response (GSR) of each patient. Valence and arousal estimation was extracted using the proposed model from the beginning until the end of the therapy, with successful evaluation to the different emotional changes in the temporal domain.


title: Shuttling-Based Trapped-Ion Quantum Information Processing

abstract: Moving trapped-ion qubits in a microstructured array of radiofrequency traps offers a route towards realizing scalable quantum processing nodes. Establishing such nodes, providing sufficient functionality to represent a building block for emerging quantum technologies, e.g. a quantum computer or quantum repeater, remains a formidable technological challenge. In this review, we present a holistic view on such an architecture, including the relevant components, their characterization and their impact on the overall system performance. We present a hardware architecture based on a uniform linear segmented multilayer trap, controlled by a custom-made fast multi-channel arbitrary waveform generator. The latter allows for conducting a set of different ion shuttling operations at sufficient speed and quality. We describe the relevant parameters and performance specifications for microstructured ion traps, waveform generators and additional circuitry, along with suitable measurement schemes to verify the system performance. Furthermore, a set of different basic shuttling operations for dynamic qubit register reconfiguration is described and characterized in detail.


title: Neural-IR-Explorer: A Content-Focused Tool to Explore Neural Re-Ranking Results

abstract: In this paper we look beyond metrics-based evaluation of Information Retrieval systems, to explore the reasons behind ranking results. We present the content-focused Neural-IR-Explorer, which empowers users to browse through retrieval results and inspect the inner workings and fine-grained results of neural re-ranking models. The explorer includes a categorized overview of the available queries, as well as an individual query result view with various options to highlight semantic connections between query-document pairs. The Neural-IR-Explorer is available at: https://neural-ir-explorer.ec.tuwien.ac.at/


title: Rare event asymptotics for exploration processes for random graphs

abstract: Much work in the study of large deviations for random graph models is focused on the dense regime where the theory of graphons has emerged as a principal tool. These tools do not give a good approach to large deviation problems for random graph models in the sparse regime. The aim of this paper is to study an approach for large deviation problems in this regime by establishing Large Deviation Principles (LDP) on suitable path spaces for certain exploration processes of the associated random graph sequence. Our work focuses on the study of one particular class of random graph models, namely the configuration model; however the general approach of using exploration processes for studying large deviation properties of sparse random graph models has broader applicability. The goal is to study asymptotics of probabilities of non-typical behavior in the large network limit. The first key step for this is to establish a LDP for an exploration process associated with the configuration model. A suitable exploration process here turns out to be an infinite dimensional Markov process with transition probability rates that diminish to zero in certain parts of the state space. Large deviation properties of such Markovian models is challenging due to poor regularity behavior of the associated local rate functions. Next, using the rate function in the LDP for the exploration process we formulate a calculus of variations problem associated with the asymptotics of component degree distributions. The second key ingredient in our study is a careful analysis of the infinite dimensional Euler-Lagrange equations associated with this calculus of variations problem. Exact solutions are identified which then provide explicit formulas for decay rates of probabilities of non-typical component degree distributions and related quantities. Please see the paper for the complete abstract.


title: Functional central limit theorems for random vectors under sub-linear expectations

abstract: The central limit theorem of martingales is the fundamental tool for studying the convergence of stochastic processes. The central limit theorem and functional central limit theorem are obtained for martingale like random variables under the sub-linear expectation by Zhang (2019). In this paper, we consider the multi-dimensional martingale like random vectors and establish a functional central limit theorem. As applications, the Lindeberg central limit theorem for independent random vectors is established, and the sufficient and necessary conditions of the central limit theorem for independent and identically distributed random vectors are obtained.


title: Spectrum Prediction and Interference Detection for Satellite Communications

abstract: Spectrum monitoring and interference detection are crucial for the satellite service performance and the revenue of SatCom operators. Interference is one of the major causes of service degradation and deficient operational efficiency. Moreover, the satellite spectrum is becoming more crowded, as more satellites are being launched for different applications. This increases the risk of interference, which causes anomalies in the received signal, and mandates the adoption of techniques that can enable the automatic and real-time detection of such anomalies as a first step towards interference mitigation and suppression. In this paper, we present a Machine Learning (ML)-based approach able to guarantee a real-time and automatic detection of both short-term and long-term interference in the spectrum of the received signal at the base station. The proposed approach can localize the interference both in time and in frequency and is universally applicable across a discrete set of different signal spectra. We present experimental results obtained by applying our method to real spectrum data from the Swedish Space Corporation. We also compare our ML-based approach to a model-based approach applied to the same spectrum data and used as a realistic baseline. Experimental results show that our method is a more reliable interference detector.


title: Feeding of active galactic nuclei by dynamical perturbations

abstract: There possibly was an AGN episode in the Galactic Centre about 6 Myr ago, powerful enough to produce the Fermi bubbles. We present numerical simulations of a possible scenario giving rise to an activity episode: a collision between a central gas ring surrounding the supermassive black hole (SMBH) and an infalling molecular cloud. We investigate different initial collision angles between the cloud and the ring. We follow the hydrodynamical evolution of the system following the collision using Gadget-3 hybrid N-body/SPH code and calculate the feeding rate of the SMBH accretion disc. This rate is then used as an input for a 1D thin $\alpha$-disc model in order to calculate the AGN luminosity. By varying the disc feeding radii we determine the limiting values for possible AGN accretion disc luminosity. Small angle collisions do not result in significant mass transport to the centre of the system, while models with highest collision angles transport close to $40\%$ of the initial matter to the accretion disc. Even with ring and cloud masses equal to $10^4 \, M_\odot$, which is the lower limit of present-day mass of the Circumnuclear ring in the Galactic Centre, the energy released over an interval of 1.5 Myr can produce $\sim 10\%$ of that required to inflate the Fermi bubbles. If the gas ring in the Galactic Centre 6 Myr ago had a mass of at least $10^5 \, M_\odot$, our proposed scenario can explain the formation of the Fermi bubbles. We estimate that such high-impact collisions might occur once every $\sim 10^8$ yr in our Galaxy.


title: Duality of sum of nonnegative circuit polynomials and optimal SONC bounds

abstract: Circuit polynomials are polynomials satisfying a number of conditions that make it easy to compute sharp and certifiable global lower bounds for them. Consequently, one may use them to find certifiable lower bounds for any polynomial by writing it as a sum of circuit polynomials with known lower bounds (if possible), in a fashion similar to the better-known sum-of-squares polynomials. Seidler and de Wolff recently showed that sums of nonnegative circuit (SONC) polynomials can be used to compute global lower bounds (called SONC bounds) for polynomials in this manner in polynomial time, as long as the polynomial is bounded from below and its support satisfies a nondegeneracy assumption. The quality of the SONC bound depends on the circuits used in the computation, but finding the set of circuits that yield the best attainable SONC bound among the astronomical number of candidate circuits is a non-trivial task that has not been addressed so far. In this paper we propose an efficient method to compute the optimal SONC lower bound by iteratively identifying the optimal circuits to use in the SONC bounding process. The method is based on a new proof of a recent result by Wang which states that (under the same nondegeneracy assumption) every SONC polynomial decomposes into SONC polynomials on the same support. Our proof, based on convex programming duality, removes the nondegeneracy assumption and motivates an algorithm that generates an optimal set of circuits and computes the corresponding SONC bound in a manner that is particularly attractive for sparse polynomials. The method is implemented and tested on a large set of sparse polynomial optimization problems with up to 40 unknowns, of degree up to 60, and up to 3000 monomials in the support. The results indicate that the method is efficient in practice. [abstract truncated for arXiv character limits]


title: PLIERS: A Process that Integrates User-Centered Methods into Programming Language Design

abstract: Programming language design requires making many usability-related design decisions. However, existing HCI methods can be impractical to apply to programming languages: they have high iteration costs, programmers require significant learning time, and user performance has high variance. To address these problems, we adapted both formative and summative HCI methods to make them more suitable for programming language design. We integrated these methods into a new process, PLIERS, for designing programming languages in a user-centered way. We evaluated PLIERS by using it to design two new programming languages. Glacier extends Java to enable programmers to express immutability properties effectively and easily. Obsidian is a language for blockchains that includes verification of critical safety properties. Summative usability studies showed that programmers were able to program effectively in both languages after short training periods.


title: Calculating the static gravitational two-body potential to fifth post-Newtonian order with Feynman diagrams

abstract: We discuss the first-time calculation of the static gravitational two-body potential up to fifth post-Newtonian(PN) order. The results are achieved through a manifest factorization property of the odd PN diagrams. The factorization property is illustrated also at first and third PN order.


title: Arbitrary optical wave evolution with Fourier transforms and phase masks

abstract: A large number of applications in classical and quantum photonics require the capability of implementing arbitrary linear unitary transformations on a set of optical modes. In a seminal work by Reck et al. it was shown how to build such multiport universal interferometers with a mesh of beam splitters and phase shifters, and this design became the basis for most experimental implementations in the last decades. However, the design of Reck et al. is difficult to scale up to a large number of modes, which would be required for many applications. Here we present a constructive proof that it is possible to realize a multiport universal interferometer on N modes with a succession of 6N Fourier transforms and 6N+1 phase masks, for any even integer N. Furthermore, we provide an algorithm to find the correct succesion of Fourier transforms and phase masks to realize a given arbitrary unitary transformation. Since Fourier transforms and phase masks are routinely implemented in several optical setups and they do not suffer from the scalability issues associated with building extensive meshes of beam splitters, we believe that our design can be useful for many applications in photonics.


title: Subleading corrections to the free energy in a theory with $N^{5/3}$ scaling

abstract: We numerically investigate the sphere partition function of a Chern-Simons-matter theory with $SU(N)$ gauge group at level $k$ coupled to three adjoint chiral multiplets that is dual to massive IIA theory. Beyond the leading order $N^{5/3}$ behavior of the free energy, we find numerical evidence for a term of the form $(2/9)\log N-(1/18)\log k$. We conjecture that the $(2/9)\log N$ term may be universal in theories with $N^{5/3}$ scaling in the large-$N$ limit with the Chern-Simons level $k$ held fixed.


title: Interfacial Load Monitoring and Failure Detection in Total Joint Replacements via Piezoresistive Bone Cement and Electrical Impedance Tomography

abstract: Aseptic loosening, or loss of implant fixation, is a common complication following total joint replacement. Revision surgeries cost the healthcare system over $8 billion annually in the US. Despite the prevalence of aseptic loosening, timely and accurate detection remains a challenge because traditional imaging modalities such as plain radiographs struggle to reliably detect the early stages of implant loosening. Motivated by this challenge, we present a novel approach for in vivo monitoring and failure detection of cemented joint replacements. Poly(methyl methacrylate) (PMMA) bone cement is modified with low volume fractions of chopped carbon fiber (CF) to impart piezoresistive-based self-sensing. Electrical impedance tomography (EIT) is then used to detect and monitor load-induced deformation and fracture of CF/PMMA in a phantom tank. We therefore show that EIT indeed is able to adeptly detect loading force on a prosthetic surrogate, distinguish between increasing load magnitudes, detect failure of implant fixation, and even distinguish between cement cracking and cement de-bonding. Because EIT is a low-cost, physiologically benign, and potentially real-time imaging modality, the feasibility study herein presented has potential to transform the standard of care for post-operative monitoring of implant conditions. Beyond clinical relevance, this technique could positively impact orthopedic researchers by providing, via in vivo monitoring, insight into the factors that initiate aseptic loosening.


title: A review of mathematical representations of biomolecules

abstract: Recently, machine learning (ML) has established itself in various worldwide benchmarking competitions in computational biology, including Critical Assessment of Structure Prediction (CASP) and Drug Design Data Resource (D3R) Grand Challenges. However, the intricate structural complexity and high ML dimensionality of biomolecular datasets obstruct the efficient application of ML algorithms in the field. In addition to data and algorithm, an efficient ML machinery for biomolecular predictions must include structural representation as an indispensable component. Mathematical representations that simplify the biomolecular structural complexity and reduce ML dimensionality have emerged as a prime winner in D3R Grand Challenges. This review is devoted to the recent advances in developing low-dimensional and scalable mathematical representations of biomolecules in our laboratory. We discuss three classes of mathematical approaches, including algebraic topology, differential geometry, and graph theory. We elucidate how the physical and biological challenges have guided the evolution and development of these mathematical apparatuses for massive and diverse biomolecular data. We focus the performance analysis on the protein-ligand binding predictions in this review although these methods have had tremendous success in many other applications, such as protein classification, virtual screening, and the predictions of solubility, solvation free energy, toxicity, partition coefficient, protein folding stability changes upon mutation, etc.


title: Some combinatorial results on smooth permutations

abstract: We show that any smooth permutation $\sigma\in S_n$ is characterized by the set ${\mathbf{C}}(\sigma)$ of transpositions and $3$-cycles in the Bruhat interval $(S_n)_{\leq\sigma}$, and that $\sigma$ is the product (in a certain order) of the transpositions in ${\mathbf{C}}(\sigma)$. We also characterize the image of the map $\sigma\mapsto{\mathbf{C}}(\sigma)$. As an application, we show that $\sigma$ is smooth if and only if the intersection of $(S_n)_{\leq\sigma}$ with every conjugate of a parabolic subgroup of $S_n$ admits a maximum. This also gives another approach for enumerating smooth permutations and subclasses thereof. Finally, we relate covexillary permutations to smooth ones and rephrase the results in terms of the (co)essential set in the sense of Fulton.


title: A Write-Friendly and Fast-Recovery Scheme for Security Metadata in NVM

abstract: Non-Volatile Memories (NVMs) have attracted the attentions of academia and industry, which is expected to become the next-generation memory. However, due to the nonvolatile property, NVMs become vulnerable to attacks and require security mechanisms, e.g., counter mode encryption and integrity tree, which introduce the security metadata. NVMs promise to recover these security metadata after a system crash, including the counter and integrity tree. However, unlike merkle tree reconstructed from user data, recovering SGX integrity tree (SIT) has to address the challenges from unique top-down hierarchical dependency. Moreover, writing overhead and recovery time are important metrics for evaluating persistent memory system due to the high costs of NVM writes and IT downtime. How to recover the security metadata, i.e., counter blocks and integrity tree nodes, with low write overhead and short recovery time, becomes much important. To provide a fast recovery scheme with low write overhead, we propose STAR, a cost-efficient scheme for recovering counter blocks and SGX integrity tree nodes after crashes. For fast recovery and verification, STAR synergizes the MAC and correct data, uses bitmap lines in ADR to indicate the location of stale node and constructs a cached merkle tree to verify the correctness of the recovery process. Moreover, STAR uses a multi-layer index to speed up the recovery process. STAR also allows different configurations to meet adaptive requirements for write overhead and recovery time. Our evaluation results show that the proposed STAR reduces the number of memory writes by up to 87\% compared with state-of-the-art work, Anubis, which needs extra 1x memory writes. For a 4MB security metadata cache, STAR needs 0.039s/0.023s/0.004s in three different configurations to recover the metadata cache while Anubis needs 0.020s.


title: Cosmology and Dark Matter

abstract: Cosmology and astroparticle physics give strongest possible evidence for the incompleteness of the Standard Model of particle physics. Leaving aside misterious dark energy, which may or may not be just the cosmological constant, two properties of the Universe cannot be explained by the Standard Model: dark matter and matter-antimatter asymmtery. Dark matter particles may well be discovered in foreseeable future; this issue is under intense experimental investigation. Theoretical hypotheses on the nature of the dark matter particles are numerous, so we concentrate on several well motivated candidates, such as WIMPs, axions and sterile neutrinos, and also give examples of less motivated and more elusive candidates such as fuzzy dark matter. This gives an idea of the spectrum of conceivable dark matter candidates, while certainly not exhausting it. We then consider the matter-antimatter asymmetry and discuss whether it may result from physics at 100~GeV -- TeV scale. Finally, we turn to the earliest epoch of the cosmological evolution. Although the latter topic does not appear immediately related to contemporary particle physics, it is of great interest due to its fundamental nature. We emphasize that the cosmological data, notably, on CMB anisotropies, unequivocally show that the well understood hot stage was not the earliest one. The best guess for the earlier stage is inflation, which is consistent with everything known to date; however, there are alternative scenarios. We discuss the ways to study the earliest epoch, with emphasis on future cosmological observations.


title: Primitivoids and inversions of plane curves

abstract: The pedal of a curve in the Euclidean plane is a classical subject which has a singular point at the inflection point of the original curve. The primitive of a curve is a curve given by the inverse construction for making the pedal. We consider relatives of the primitive of a plane curve which we call primitivoids. We investigate the relationship of primitivoids and pedals of plane curves.


title: Accurate OH Maser Positions From The SPLASH Survey III: The Final 96 Square Degrees

abstract: We present high spatial resolution observations of ground-state OH masers achieved with the Australia Telescope Compact Array (ATCA). These observations targeted 253 pointing centres containing OH maser candidates at all four ground-state OH transitions identified in the Southern Parkes Large-Area Survey in Hydroxyl (SPLASH) across 96 square degrees of the Southern Galactic plane (332\degree$<l<$334\degree and $-$2\degree$<b<+$2\degree, 344\degree$<l<$355\degree and $-$2\degree$<b<+$2\degree, 358\degree$<l<$4\degree and $+$2\degree$<b<+$6\degree, 5\degree$<l<$10\degree and $-$2\degree$<b<+$2\degree). We detect maser emission towards 236 fields and suggest that 7 out of 17 non-detections are due to the slightly lower sensitivity of the ATCA observations, combined with some temporal variability. The superior resolution provided by the ATCA data has allowed us to identify 362 OH maser sites in the 236 target fields. Almost half (160 of 362) of these masers have been detected for the first time. Comparison between these 362 maser sites with information presented in the literature allowed us to categorize 238 sites as evolved star sites (66\%), 63 as star formation (17\%), eight as supernova remnants and 53 unknown maser sites (15\%). We present analysis of the OH masers across the full SPLASH survey range (176 square degrees) and find that the detection rate of 1.7 GHz radio continuum sources (18\%) is lower than that previously found at 8.2 and 9.2 GHz (38\%). We also find that the velocity separations of evolved star sites with symmetric 1612 MHz maser profiles are generally smaller than those with asymmetric profiles.


title: On the oriented Thompson subgroup $\vec{F}_3$ and its relatives in higher Brown-Thompson groups

abstract: A few years ago the so-called oriented subgroup $\vec F$ of the Thompson group $F$ was introduced by V. Jones while investigating the connections between subfactors and conformal field theories. In the coding of links and knots by elements of $F$ it corresponds exactly to the oriented ones. Thanks to the work of Golan and Sapir, $\vec F$ provided the first example of a maximal subgroup of infinite index in $F$ different from the parabolic subgroups that fix a point in $(0,1)$. In this paper we investigate possible analogues of $\vec F$ in higher Thompson groups $F_k, k\geq 2$, with $F=F_2$, introduced by Brown. Most notably, we study algebraic properties of the oriented subgroup $\vec{F}_3$ of $F_3$, as described recently by Jones, and prove in particular that it gives rise to a non-parabolic maximal subgroup of infinite index in $F_3$ and that the corresponding quasi-regular representation is irreducible.


title: On a question of Dikranjan and Zava

abstract: Let $G$ be a non-discrete countable metrizable abelian topological group endowed with the coarse structure $ \mathcal{C} $ generated by compact subsets of $G$. We prove that $asdim (G, \mathcal{C} ) = \infty$. For an infinite cyclic subgroup $G$ of the circle, this answers a question of Dikranjan and Zava [3].


title: Charge Constraints of Macroscopic Dark Matter

abstract: Macroscopic dark matter (macros) refers to a broad class of alternative candidates to particle dark matter with still unprobed regions of parameter space. Prior work on macros has considered elastic scattering to be the dominant energy transfer mechanism in deriving constraints on the abundance of macros for some range of masses $M_x$ and (geometric) cross-sections $\sigma_x$ However, macros with a significant amount of electric charge would, through Coulomb interactions, interact strongly enough to have produced observable signals on terrestrial, galactic and cosmological scales. We determine the expected signals and constrain the corresponding regions of parameter space, based on the lack of these signals in observations.


title: Off-Grid Aware Spatial Covariance Estimation in mmWave Communications

abstract: This work investigates the problem of spatial covariance matrix estimation in a millimeter-wave (mmWave) hybrid multiple-input multiple-output (MIMO) system with an emphasis on the basis-mismatch effect. The basis mismatch is prevalent in the compressed sensing (CS) schemes which adopt discretization procedure. In such an approach, the algorithm yields a finite discrete point which is an approximation to the continuous parametric space. The quality of this approximation depends on the number of discretized points in the dictionary. Instead of increasing the number of discretized points to combat this off-grid effect, we propose an efficient parameter perturbed framework which uses a controlled perturbation mechanism in conjunction with the orthogonal matching pursuit (OMP) algorithm. Numerical results verify the performance improvement through our proposed algorithm in terms of relative efficiency metric, which is basically due to taking care of the off-grid effect carefully that is ignored in the conventional CS algorithms.


title: Transformed Subspace Clustering

abstract: Subspace clustering assumes that the data is sepa-rable into separate subspaces. Such a simple as-sumption, does not always hold. We assume that, even if the raw data is not separable into subspac-es, one can learn a representation (transform coef-ficients) such that the learnt representation is sep-arable into subspaces. To achieve the intended goal, we embed subspace clustering techniques (locally linear manifold clustering, sparse sub-space clustering and low rank representation) into transform learning. The entire formulation is jointly learnt; giving rise to a new class of meth-ods called transformed subspace clustering (TSC). In order to account for non-linearity, ker-nelized extensions of TSC are also proposed. To test the performance of the proposed techniques, benchmarking is performed on image clustering and document clustering datasets. Comparison with state-of-the-art clustering techniques shows that our formulation improves upon them.


title: Covert Channel-Based Transmitter Authentication in Controller Area Networks

abstract: In recent years, the security of automotive Cyber-Physical Systems (CPSs) is facing urgent threats due to the widespread use of legacy in-vehicle communication systems. As a representative legacy bus system, the Controller Area Network (CAN) hosts Electronic Control Units (ECUs) that are crucial vehicle functioning. In this scenario, malicious actors can exploit CAN vulnerabilities, such as the lack of built-in authentication and encryption schemes, to launch CAN bus attacks with life-threatening consequences (e.g., disabling brakes). In this paper, we present TACAN (Transmitter Authentication in CAN), which provides secure authentication of ECUs on the legacy CAN bus by exploiting the covert channels, without introducing CAN protocol modifications or traffic overheads. TACAN turns upside-down the originally malicious concept of covert channels and exploits it to build an effective defensive technique that facilitates transmitter authentication via a centralized, trusted Monitor Node. TACAN consists of three different covert channels for ECU authentication: 1) the Inter-Arrival Time (IAT)-based; 2) the Least Significant Bit (LSB)-based; and 3) a hybrid covert channel, exploiting the combination of the first two. In order to validate TACAN, we implement the covert channels on the University of Washington (UW) EcoCAR (Chevrolet Camaro 2016) testbed. We further evaluate the bit error, throughput, and detection performance of TACAN through extensive experiments using the EcoCAR testbed and a publicly available dataset collected from Toyota Camry 2010. We demonstrate the feasibility of TACAN and the effectiveness of detecting CAN bus attacks, highlighting no traffic overheads and attesting the regular functionality of ECUs.


title: Theoretical status of $\varepsilon'/\varepsilon$

abstract: We briefly overview the historical controversy around Standard Model predictions of $\varepsilon'/\varepsilon$ and clarify the underlying physics. A full update of this important observable is presented, with all known short- and long-distance contributions, including isospin-breaking corrections. The current Standard Model prediction, $\mathrm{Re}(\varepsilon'/\varepsilon) = (14\pm 5)\cdot 10^{-4}$, is in excellent agreement with the experimentally measured value.


title: Solving Partial Differential Equations with Neural Networks

abstract: Many scientific and industrial applications require solving Partial Differential Equations (PDEs) to describe the physical phenomena of interest. Some examples can be found in the fields of aerodynamics, astrodynamics, combustion and many others. In some exceptional cases an analytical solution to the PDEs exists, but in the vast majority of the applications some kind of numerical approximation has to be computed. In this work, an alternative approach is proposed using neural networks (NNs) as the approximation function for the PDEs. Unlike traditional numerical methods, NNs have the property to be able to approximate any function given enough parameters. Moreover, these solutions are continuous and derivable over the entire domain removing the need for discretization. Another advantage that NNs as function approximations provide is the ability to include the free-parameters in the process of finding the solution. As a result, the solution can generalize to a range of situations instead of a particular case, avoiding the need of performing new calculations every time a parameter is changed dramatically decreasing the optimization time. We believe that the presented method has the potential to disrupt the physics simulation field enabling real-time physics simulation and geometry optimization without the need of big supercomputers to perform expensive and time consuming simulations


title: Histogram Transform Ensembles for Large-scale Regression

abstract: We propose a novel algorithm for large-scale regression problems named histogram transform ensembles (HTE), composed of random rotations, stretchings, and translations. First of all, we investigate the theoretical properties of HTE when the regression function lies in the H\"{o}lder space $C^{k,\alpha}$, $k \in \mathbb{N}_0$, $\alpha \in (0,1]$. In the case that $k=0, 1$, we adopt the constant regressors and develop the na\"{i}ve histogram transforms (NHT). Within the space $C^{0,\alpha}$, although almost optimal convergence rates can be derived for both single and ensemble NHT, we fail to show the benefits of ensembles over single estimators theoretically. In contrast, in the subspace $C^{1,\alpha}$, we prove that if $d \geq 2(1+\alpha)/\alpha$, the lower bound of the convergence rates for single NHT turns out to be worse than the upper bound of the convergence rates for ensemble NHT. In the other case when $k \geq 2$, the NHT may no longer be appropriate in predicting smoother regression functions. Instead, we apply kernel histogram transforms (KHT) equipped with smoother regressors such as support vector machines (SVMs), and it turns out that both single and ensemble KHT enjoy almost optimal convergence rates. Then we validate the above theoretical results by numerical experiments. On the one hand, simulations are conducted to elucidate that ensemble NHT outperform single NHT. On the other hand, the effects of bin sizes on accuracy of both NHT and KHT also accord with theoretical analysis. Last but not least, in the real-data experiments, comparisons between the ensemble KHT, equipped with adaptive histogram transforms, and other state-of-the-art large-scale regression estimators verify the effectiveness and accuracy of our algorithm.


title: The confining color field in SU(3) gauge theory

abstract: We extend a previous numerical study of SU(3) Yang-Mills theory in which we measured the spatial distribution of all components of the color fields surrounding a static quark-antiquark pair for a wide range of quark-antiquark separations, and provided evidence that the simulated gauge invariant chromoelectric field can be separated into a Coulomb-like 'perturbative' field and a 'non-perturbative' field, identified as the confining part of the SU(3) flux tube field. In this paper we hypothesize that the fluctuating color fields not measured in our simulations do not contribute to the string tension. Under this assumption the string tension is determined by the color fields we measure, which form a tensor $F_{\mu \nu}$ pointing in a single direction in color space. We call this the Maxwell mechanism of confinement. We provide an additional procedure to isolate the non-perturbative (confining) field. We then extract the string tension from a stress energy-momentum tensor $T_{\mu \nu}$ having the Maxwell form, constructed from the non-perturbative part of the tensor $F_{\mu \nu}$ obtained from our simulations. To test our hypothesis we calculate the string tension from our simulations of the color fields for ten values of the quark-antiquark separation ranging from 0.37 fm to 1.2 fm. We also calculate the spatial distributions of the energy-momentum tensor $T_{\mu \nu}$ surrounding static quarks for this range of separations, and we compare these distributions with those obtained from direct simulations of the energy-momentum tensor in SU(3) Yang-Mills theory.


title: A note on the spectra and eigenspaces of the universal adjacency matrices of arbitrary lifts of graphs

abstract: The universal adjacency matrix $U$ of a graph $\Gamma$, with adjacency matrix $A$, is a linear combination of $A$, the diagonal matrix $D$ of vertex degrees, the identity matrix $I$, and the all-1 matrix $J$ with real coefficients, that is, $U=c_1 A+c_2 D+c_3 I+ c_4 J$, with $c_i\in {\mathbb R}$ and $c_1\neq 0$. Thus, as particular cases, $U$ may be the adjacency matrix, the Laplacian, the signless Laplacian, and the Seidel matrix. In this note, we show that basically the same method introduced before by the authors can be applied for determining the spectra and bases of all the corresponding eigenspaces of arbitrary lifts of graphs (regular or not).


title: Sequential collision-free optimal motion planning algorithms in punctured Euclidean spaces

abstract: In robotics, a topological theory of motion planning was initiated by M. Farber. The multitasking motion planning problem is new and its theoretical part via topological complexity has hardly been developed, but the concrete implementations are still non-existent, and in fact this work takes the first step in this last direction (producing explicit algorithms.) We present optimal motion planning algorithms which can be used in designing practical systems controlling objects moving in Euclidean space without collisions between them and avoiding obstacles. Furthermore, we present the multitasking version of the algorithms.


title: Ensemble Estimation of Network Parameters: A Tool to Improve the Real-time Estimation of GICs in the South African Power Network

abstract: It has long been known that large grounded conducting networks on the surface of Earth are affected by solar activity and geomagnetic storms. Power networks are such extensive grounded conductors and are susceptible to geomagnetically induced currents (GICs). GICs at any specific node in a power network are assumed to be linearly related to the horizontal vector components of an induced plane-wave geoelectric field by a pair of network parameters. These network parameters are not easily measured in the network, but may be estimated empirically. In this work, we present a new approach of using an ensemble of network parameters estimates. The ensembles include a huge number of parameter pair estimates calculated from simultaneously solving pairs of time instances of the governing GIC equation. Each individual estimate is not the true state of the system, but a possible state. Taking the ensemble as a whole though gives the most probable parameter estimate. The most probable parameter estimate for both network parameters, as defined by their respective ensembles, is used directly in the modelling of GICs. The ensembles themselves however allow for further analysis into the nature of GICs. An improvement is seen when comparing the out-of-sample performance of the ensemble estimates with previous GIC modelling in the South African power network during the Halloween Storm of 2003. For the first time, it is shown that errors in the GIC modelling chain are absorbed into the network parameter estimates. Using a range of estimates from the ensemble, a GIC prediction band is produced. This band corresponds to an error estimate for predicted GIC. Furthermore, it has been explicitly shown that estimated network parameters vary with GIC magnitude during an event. This behaviour is then used to refine the parameter estimation further and allow for real time dynamic network parameter estimation.


title: dSECS: Including Variometers in Geomagnetic Field Interpolation

abstract: One of the best tools we have in geophysical modelling is the ability to interpolate the horizontal geomagnetic field at the surface of the Earth. This is especially useful in regions, such as southern Africa, where there is a sparse array of absolute magnetometers available for geomagnetic field measurements. In terms of geomagnetic field interpolation, the spherical elementary current systems (SECS) spatial interpolation scheme has shown to be very successful, and the planar approximation of this method adequate for modelling at mid-latitudes. The SECS interpolation scheme is physics based, making use of the Biot-Savart law and equivalent ionospheric currents to interpolate measured geomagnetic field data. As with most interpolation methods, more data points result in lower error. Therefore, we adapt the SECS method to work with variometers. These instruments measure variations in magnetic field and are more abundant in southern Africa. Merging the two resulting interpolated datasets, the initial absolute geomagnetic field interpolation can be significantly improved. This improved interpolation scheme is not only incredibly useful locally, where a sparse magnetometer array is a challenge, but can also be applied just as effectively in other cases across the globe where there are numerous magnetometers and variometers available.


title: Near-Optimal Rapid MPC using Neural Networks: A Primal-Dual Policy Learning Framework

abstract: In this paper, we propose a novel framework for approximating the explicit MPC policy for linear parameter-varying systems using supervised learning. Our learning scheme guarantees feasibility and near-optimality of the approximated MPC policy with high probability. Furthermore, in contrast to most existing approaches that only learn the MPC policy, we also learn the "dual policy", which enables us to keep a check on the approximated MPC's optimality online during the control process. If the check deems the control input from the approximated MPC policy safe and near-optimal, then it is applied to the plant, otherwise a backup controller is invoked, thus filtering out (severely) suboptimal control inputs. The backup controller is only invoked with a bounded (low) probability, where the exact probability level can be chosen by the user. Since our framework does not require solving any optimization problem during the control process, it enables the deployment of MPC on resource-constrained systems. Specifically, we illustrate the utility of the proposed framework on a vehicle dynamics control problem. Compared to online optimization methods, we demonstrate a speedup of up to 62x on a desktop computer and 10x on an automotive-grade electronic control unit, while maintaining a high control performance.


title: Infinite families of $2$-designs from a class of non-binary Kasami cyclic codes

abstract: Combinatorial $t$-designs have been an important research subject for many years, as they have wide applications in coding theory, cryptography, communications and statistics. The interplay between coding theory and $t$-designs has been attracted a lot of attention for both directions. It is well known that a linear code over any finite field can be derived from the incidence matrix of a $t$-design, meanwhile, that the supports of all codewords with a fixed weight in a code also may hold a $t$-design. In this paper, by determining the weight distribution of a class of linear codes derived from non-binary Kasami cyclic codes, we obtain infinite families of $2$-designs from the supports of all codewords with a fixed weight in these codes, and calculate their parameters explicitly.


title: On the Security of Pixel-Based Image Encryption for Privacy-Preserving Deep Neural Networks

abstract: This paper aims to evaluate the safety of a pixel-based image encryption method, which has been proposed to apply images with no visual information to deep neural networks (DNN), in terms of robustness against ciphertext-only attacks (COA). In addition, we propose a novel DNN-based COA that aims to reconstruct the visual information of encrypted images. The effectiveness of the proposed attack is evaluated under two encryption key conditions: same encryption key, and different encryption keys. The results show that the proposed attack can recover the visual information of the encrypted images if images are encrypted under same encryption key. Otherwise, the pixel-based image encryption method has robustness against COA.


title: Oversampling Log Messages Using a Sequence Generative Adversarial Network for Anomaly Detection and Classification

abstract: Dealing with imbalanced data is one of the main challenges in machine/deep learning algorithms for classification. This issue is more important with log message data as it is typically very imbalanced and negative logs are rare. In this paper, a model is proposed to generate text log messages using a SeqGAN network. Then features are extracted using an Autoencoder and anomaly detection is done using a GRU network. The proposed model is evaluated with two imbalanced log data sets, namely BGL and Openstack. Results are presented which show that oversampling and balancing data increases the accuracy of anomaly detection and classification.


title: Fraud detection in telephone conversations for financial services using linguistic features

abstract: Detecting the elements of deception in a conversation is one of the most challenging problems for the AI community. It becomes even more difficult to design a transparent system, which is fully explainable and satisfies the need for financial and legal services to be deployed. This paper presents an approach for fraud detection in transcribed telephone conversations using linguistic features. The proposed approach exploits the syntactic and semantic information of the transcription to extract both the linguistic markers and the sentiment of the customer's response. We demonstrate the results on real-world financial services data using simple, robust and explainable classifiers such as Naive Bayes, Decision Tree, Nearest Neighbours, and Support Vector Machines.


title: Efficient Differentiable Neural Architecture Search with Meta Kernels

abstract: The searching procedure of neural architecture search (NAS) is notoriously time consuming and cost prohibitive.To make the search space continuous, most existing gradient-based NAS methods relax the categorical choice of a particular operation to a softmax over all possible operations and calculate the weighted sum of multiple features, resulting in a large memory requirement and a huge computation burden. In this work, we propose an efficient and novel search strategy with meta kernels. We directly encode the supernet from the perspective on convolution kernels and "shrink" multiple convolution kernel candidates into a single one before these candidates operate on the input feature. In this way, only a single feature is generated between two intermediate nodes. The memory for storing intermediate features and the resource budget for conducting convolution operations are both reduced remarkably. Despite high efficiency, our search strategy can search in a more fine-grained way than existing works and increases the capacity for representing possible networks. We demonstrate the effectiveness of our search strategy by conducting extensive experiments. Specifically, our method achieves 77.0% top-1 accuracy on ImageNet benchmark dataset with merely 357M FLOPs, outperforming both EfficientNet and MobileNetV3 under the same FLOPs constraints. Compared to models discovered by the start-of-the-art NAS method, our method achieves the same (sometimes even better) performance, while faster by three orders of magnitude.


title: Two-color phase-of-the-phase spectroscopy applied to nonperturbative electron-positron pair production in strong oscillating electric fields

abstract: Production of electron-positron pairs from vacuum in strong bichromatic electric fields, oscillating in time with a fundamental frequency and its second harmonic, is studied. Strong-field processes occuring in such field configurations are generally known to be sensitive to the relative phase between the field modes. Phase-of-the-phase spectroscopy has recently been introduced in the context of strong-field photoionization as a systematic means to analyze these coherence effects. We apply this method to field-induced pair production by calculating the phase dependence of the momentum-resolved particle yields. We show that asymmetric checkerboard patterns arise in the phase-of-the-phase spectra, similarly to those found in strong-field photoionization. The physical origin of these characteristic structures, which differ between the created electron and positron, are discussed.


title: Characterisation of the ground layer of turbulence at Paranal using a robotic SLODAR system

abstract: We describe the implementation of a robotic SLODAR instrument at the Cerro Paranal observatory. The instrument measures the vertical profile of the optical atmospheric turbulence strength, in 8 resolution elements, to a maximum altitude ranging between 100 m and 500 m. We present statistical results of measurements of the turbulence profile on a total of 875 nights between 2014 and 2018. The vertical profile of the ground layer of turbulence is very varied, but in the median case most of the turbulence strength in the ground layer is concentrated within the first 50 m altitude, with relatively weak turbulence at higher altitudes up to 500 m. We find good agreement between measurements of the seeing angle from the SLODAR and from the Paranal DIMM seeing monitor, and also for seeing values extracted from the Shack-Hartmann active optics sensor of VLT UT1, adjusting for the height of each instrument above ground level. The SLODAR data suggest that a median improvement in the seeing angle from 0.689 arcsec to 0.481 arcsec at wavelength 500 nm would be obtained by fully correcting the ground-layer turbulence between the height of the UTs (taken as 10 m) and altitude 500 m.


title: Closed conformal Killing-Yano initial data

abstract: Through an exhaustive search, we produce a 5-parameter family of propagation identities for the closed conformal Killing-Yano equation on 2-forms, which hold on an Einstein cosmological vacuum spacetime in any dimension $n>4$. It is well-known that spacetimes admitting a non-degenerate 2-form of this type are exhausted by the Kerr-NUT-(A)dS family of exact higher dimensional black hole solutions. As a consequence, we identify a set of necessary and sufficient conditions ensuring that the cosmological vacuum development of an initial data set for Einstein's field equations admits a closed conformal Killing-Yano 2-form. We refer to these conditions as \emph{closed conformal Killing-Yano initial data} (cCYKID) equations. The 4-dimensional case is special and is treated separately, where we can also handle the conformal Killing-Yano equation without the closed condition.


title: Optimizing and accelerating space-time Ripley's K function based on Apache Spark for distributed spatiotemporal point pattern analysis

abstract: With increasing point of interest (POI) datasets available with fine-grained spatial and temporal attributes, space-time Ripley's K function has been regarded as a powerful approach to analyze spatiotemporal point process. However, space-time Ripley's K function is computationally intensive for point-wise distance comparisons, edge correction and simulations for significance testing. Parallel computing technologies like OpenMP, MPI and CUDA have been leveraged to accelerate the K function, and related experiments have demonstrated the substantial acceleration. Nevertheless, previous works have not extended optimization of Ripley's K function from space dimension to space-time dimension. Without sophisticated spatiotemporal query and partitioning mechanisms, extra computational overhead can be problematic. Meanwhile, these researches were limited by the restricted scalability and relative expensive programming cost of parallel frameworks and impeded their applications for large POI dataset and Ripley's K function variations. This paper presents a distributed computing method to accelerate space-time Ripley's K function upon state-of-the-art distributed computing framework Apache Spark, and four strategies are adopted to simplify calculation procedures and accelerate distributed computing respectively. Based on the optimized method, a web-based visual analytics framework prototype has been developed. Experiments prove the feasibility and time efficiency of the proposed method, and also demonstrate its value on promoting applications of space-time Ripley's K function in ecology, geography, sociology, economics, urban transportation and other fields.


title: Deep Latent Factor Model for Collaborative Filtering

abstract: Latent factor models have been used widely in collaborative filtering based recommender systems. In recent years, deep learning has been successful in solving a wide variety of machine learning problems. Motivated by the success of deep learning, we propose a deeper version of latent factor model. Experiments on benchmark datasets shows that our proposed technique significantly outperforms all state-of-the-art collaborative filtering techniques.


title: Additive manufacturing for energy storage: Methods, designs and materials selection for customizable 3D printed batteries and supercapacitors

abstract: Additive manufacturing and 3D printing in particular have the potential to revolutionize existing fabrication processes where objects with complex structures and shapes can be built with multifunctional material systems. For electrochemical energy storage devices such as batteries and supercapacitors, 3D printing methods allows alternative form factors to be conceived based on the end use application need in mind at the design stage. Additively manufactured energy storage devices require active materials and composites that are printable and this is influenced by performance requirements and the basic electrochemistry. The interplay between electrochemical response, stability, material type, object complexity and end use application are key to realising 3D printing for electrochemical energy storage. Here, we summarise recent advances and highlight the important role of methods, designs and material selection for energy storage devices made by 3D printing, which is general to the majority of methods in use currently.


title: Electron-Acoustic Solitons in Magnetized Collisional Nonthermal Plasmas

abstract: The properties of obliquely propagating dissipative electron-acoustic solitary waves (OPdEASWs) have been investigated in a magnetized collisional superthermal plasma consisting of inertial cold electrons, inertialess hot electrons featuring kappa-distribution and static ions via the fluid dynamical approach. By using the reductive perturbation technique, a nonlinear Schamel equation that governs the nonlinear features of OPdEASWS is obtained. The solitary wave solution of the Schamel equation is used to examine the basic features of small, but finite amplitude OPdEASWs in such a magnetized collisional superthermal plasma. The basic features (width, amplitude, speed, etc.) of OPdEASWs are found to be modified by the different plasma configuration parameters, such as the superthermality index, obliquity angle, collisional parameter, trapping parameter, and external magnetic field. The nature of electrostatic disturbances, that may propagate in different realistic space and laboratory plasma systems (e.g., in Saturn ring), are briefly discussed.


title: Sign-tunable anomalous Hall effect induced by two-dimensional symmetry-protected nodal structures in ferromagnetic perovskite oxide thin films

abstract: Magnetism and spin-orbit coupling (SOC) are two quintessential ingredients underlying novel topological transport phenomena in itinerant ferromagnets. When spin-polarized bands support nodal points/lines with band degeneracy that can be lifted by SOC, the nodal structures become a source of Berry curvature; this leads to a large anomalous Hall effect (AHE). Contrary to three-dimensional systems that naturally host nodal points/lines, two-dimensional (2D) systems can possess stable nodal structures only when proper crystalline symmetry exists. Here we show that 2D spin-polarized band structures of perovskite oxides generally support symmetry-protected nodal lines and points that govern both the sign and the magnitude of the AHE. To demonstrate this, we performed angle-resolved photoemission studies of ultrathin films of SrRuO$_3$, a representative metallic ferromagnet with SOC. We show that the sign-changing AHE upon variation in the film thickness, magnetization, and chemical potential can be well explained by theoretical models. Our study is the first to directly characterize the topological band structure of 2D spin-polarized bands and the corresponding AHE, which could facilitate new switchable devices based on ferromagnetic ultrathin films.


title: Generalised Network Autoregressive Processes and the GNAR package

abstract: This article introduces the GNAR package, which fits, predicts, and simulates from a powerful new class of generalised network autoregressive processes. Such processes consist of a multivariate time series along with a real, or inferred, network that provides information about inter-variable relationships. The GNAR model relates values of a time series for a given variable and time to earlier values of the same variable and of neighbouring variables, with inclusion controlled by the network structure. The GNAR package is designed to fit this new model, while working with standard ts objects and the igraph package for ease of use.


title: Pore Network and Medial Axis simultaneous extraction through Maximal Ball Algorithm

abstract: Network Extraction algorithms from X-ray microcomputed tomography have become a routine method to obtain pore connectivity and pore morphology information from porous media. The main approaches for this extraction are either based on Max Ball Algorithm or Medial Axis Extraction. The first is a robust method to separate the pore space into pore and throats, which provides the pore and throat sizes distributions directly. The second simplifies the medium by thinning the volume into a one voxel wide centerline that preserves the volume's topological information. Since each method has drastically different implementations, it is not usual to use both to characterize the porous structure. This work presents a method to extract a simplified centerline of porous materials by adding few more steps to the well-established Max Ball algorithm. Results for sandstones and carbonate rock samples show that this Medial-Axis network can be used for single phase flow simulations, as well as preserves the mediums morphology.

categories: physics.data-an physics.app-ph physics.comp-ph physics.flu-dyn physics.geo-ph


title: Is Your Smartband Smart Enough to Know Who You Are: Continuous Physiological Authentication in The Wild

abstract: The use of cloud services that process privacy-sensitive information such as digital banking, pervasive healthcare, smart home applications requires an implicit continuous authentication solution which will make these systems less vulnerable to the spoofing attacks. Physiological signals can be used for continuous authentication due to their personal uniqueness. Ubiquitous wrist-worn wearable devices are equipped with photoplethysmogram sensors which enable to extract heart rate variability (HRV) features. In this study, we show that these devices can be used for continuous physiological authentication, for enhancing the security of the cloud, edge services, and IoT devices. A system that is suitable for the smartband framework comes with new challenges such as relatively low signal quality and artifacts due to placement which were not encountered in full lead electrocardiogram systems. After the artifact removal, cleaned physiological signals are fed to the machine learning algorithms. In order to train our machine learning models, we collected physiological data using off-the-shelf smartbands and smartwatches in a real-life event. Performance evaluation of selected machine learning algorithms shows that HRV is a strong candidate for continuous unobtrusive implicit physiological authentication.


title: Sound Event Detection of Weakly Labelled Data with CNN-Transformer and Automatic Threshold Optimization

abstract: Sound event detection (SED) is a task to detect sound events in an audio recording. One challenge of the SED task is that many datasets such as the Detection and Classification of Acoustic Scenes and Events (DCASE) datasets are weakly labelled. That is, there are only audio tags for each audio clip without the onset and offset times of sound events. \qk{We compare segment-wise and clip-wise training for SED that is lacking in previous works. We propose a convolutional neural network transformer (CNN-Transfomer) for audio tagging and SED, and show that CNN-Transformer performs similarly to a convolutional recurrent neural network (CRNN)}. Another challenge of SED is that thresholds are required for detecting sound events. Previous works set thresholds empirically, and are not an optimal approaches. To solve this problem, we propose an automatic threshold optimization method. The first stage is to optimize the system with respect to metrics that do not depend on thresholds, such as mean average precision (mAP). The second stage is to optimize the thresholds with respect to metrics that depends on those thresholds. Our proposed automatic threshold optimization system achieves a state-of-the-art audio tagging F1 of 0.646, outperforming that without threshold optimization of 0.629, and a sound event detection F1 of 0.584, outperforming that without threshold optimization of 0.564.


title: Dispersive Two-Loop Calculations: Methodology and Applications

abstract: As the new-generation precision experiments such as MOLLER and P2 look for physics beyond Standard Model, it is becoming increasingly important to evaluate the higher-order electroweak radiative corrections to a sub-percent level of uncertainty. However, due to propagators with different masses and higher-order tensor Feynman integrals, the two-loop calculations involving thousands of Feynman graphs become a demanding task requiring novel computational approaches. In this paper, we describe our dispersive sub-loop insertion approach and develop two-loop integrals using two-point functions basis which is applicable to wide range of processes.


title: First evidence of enhanced recombination in astrophysical environments and the implications for plasma diagnostics

abstract: We report the first unambiguous observational evidence of Rydberg Enhanced Recombination (RER), a potentially important recombination mechanism that has hitherto been unexplored in low-temperature photoionized plasmas. RER shares similarities with dielectronic recombination, with the difference that the electron is captured into a highly excited state below the ionization threshold -- rather than above the threshold -- of the recombining ion. We predict transitions of carbon and oxygen ions that are formed via the RER process, and their relative strengths with collisional-radiative spectral models. Optical C II RER features are detected in published high-resolution spectra of eight planetary nebulae, and a C III transition has been found in the ultraviolet spectrum in a symbiotic star system. The relative intensities of these lines are consistent with their production by this recombination mechanism. Because RER has not previously been accounted for in photoionized plasmas, its inclusion in models can significantly impact the predicted ionization balance and hence abundance calculations of important astrophysical species. Calculations for C+ suggest that the enhancement in the total recombination rate can amount to a factor of 2.2 at 8100 K, increasing to 7.5 at T = 3500 K. These results demonstrate the importance of including RER in models of photoionized astrophysical plasmas and in elemental abundance determinations.


title: Competition between tenants with strategic resource allocation in 5G network slicing

abstract: We propose and analyze a business model for 5G operators. Each operator is entitled to a share of a network operated by an Infrastructure Provider (InP) and use network slicing mechanisms to request network resources as needed for service provision. The network operators become Network Slice Tenants (NSTs). The InP performs the resource allocation based on a vector of weights chosen strategically by each NST. The weights distribute the NST's share of resources between its subscribers in each cell. We propose a strategy profile in which the NST chooses weights equal to the product of its share by the ratio between the total number of subscribers in the cell and the total number of subscribers in the network. We characterize the proposed solution in terms of subscription ratios and fractions of subscribers, for different cell capacities and user sensitivities. The proposed solution provides the exact values for the Nash equilibrium if the cells are homogeneous in terms of normalized capacity, which is a measure of the total amount of resources available in the cell. Otherwise, if the cells are heterogeneous, it provides an accurate approximation. We quantify the deviation from the equilibrium and conclude that it is highly accurate.


title: Formation of moon systems around giant planets: Capture and ablation of planetesimals as foundation for a pebble accretion scenario

abstract: The four major satellites of Jupiter, known as the Galilean moons, and Saturn's most massive satellite, Titan, are believed to have formed in a predominantly gaseous circum-planetary disk, during the last stages of formation of their parent planet. Pebbles from the protoplanetary disk are blocked from flowing into the circumplanetary disk by the positive pressure gradient at the outer edge of the planetary gap, so the gas drag assisted capture of planetesimals should be the main contributor to the delivery of solids onto circum-planetary disks. However, a consistent framework for the subsequent accretion of the moons remains to be built. Here we use numerical integrations to show that most planetesimals being captured within a circum-planetary disk are strongly ablated due to the frictional heating they experience, thus supplying the disk with small dust grains, whereas only a small fraction 'survives' their capture. We then construct a simple model of a circum-planetary disk supplied by ablation, where the flux of solids through the disk is at equilibrium with the ablation supply rate, and investigate the formation of moons in such disks. We show that the growth of satellites is driven mainly by accretion of the pebbles that coagulate from the ablated material. The pebble-accreting protosatellites rapidly migrate inward and pile up in resonant chains at the inner edge of the circum-planetary disk. We propose that dynamical instabilities in these resonant chains are at the origin of the different architectures of Jupiter's and Saturn's moon systems. The assembly of moon systems through pebble accretion can therefore be seen as a down-scaled manifestation of the same process that forms systems of super-Earths and terrestrial-mass planets around solar-type stars and M-dwarfs.


title: On the monotonicity of additive representation functions

abstract: We study the monotonicity behavior of three slightly differently defined additive representation functions (as initiated by Erd\H{o}s, S\'{a}rk\"{o}zy and S\'{o}s), answering one open question and another one partially, and give a slightly simpler proof for a result of Chen and Tang.


title: Dynamic Structure of Hadrons in ChPT

abstract: The Chiral Perturbation Theory (ChPT) has been very successful in describing low-energy hadronic properties in the non-perturbative regime of Quantum Chromodynamics. The results of ChPT, many of which are currently under active experimental investigation, provide stringent predictions of many fundamental properties of hadrons, including quantities such as electromagnetic polarizabilities. The paper outlines our semi-automated calculations in ChPT, the corresponding results for the electric and magnetic polarizabilities of the proton and our predictions for Compton differential cross sections.


title: Doublet-Triplet Splitting in Fertile Left-Right Symmetric Heterotic String Vacua

abstract: Classification of Left-Right Symmetric (LRS) heterotic-string vacua in the free fermionic formulation, using random generation of generalised GSO (GGSO) projection coefficients, produced phenomenologically viable models with probability $4\times 10^{-11}$. Extracting substantial number of phenomenologically viable models requires modification of the classification method. This is achieved by identifying phenomenologically amenable conditions on the Generalised GSO projection coefficients that are randomly generated at the $SO(10)$ level. Around each of these fertile cores we perform a complete LRS classification, generating viable models with probabilility $1.4\times 10^{-2}$, hence increasing the probability of generating phenomenologically viable models by nine orders of magnitude, and producing some $1.4\times 10^5$ such models. In the process we identify a doublet-triplet selection mechanism that operates in twisted sectors of the string models that break the $SO(10)$ symmetry to the Pati-Salam subgroup. This mechanism therefore operates as well in free fermionic models with Pati-Salam and Standard-like Model $SO(10)$ subgroups.


title: Post-quantum Zero Knowledge in Constant Rounds

abstract: We construct a constant-round zero-knowledge classical argument for NP secure against quantum attacks. We assume the existence of Quantum Fully-Homomorphic Encryption and other standard primitives, known based on the Learning with Errors Assumption for quantum algorithms. As a corollary, we also obtain a constant-round zero-knowledge quantum argument for QMA. At the heart of our protocol is a new no-cloning non-black-box simulation technique.


title: Spin splitting and strain in epitaxial monolayer WSe$_2$ on graphene

abstract: We present the electronic and structural properties of monolayer WSe$_{2}$ grown by pulsed-laser deposition on monolayer graphene (MLG) on SiC. The spin splitting in the WSe$_{2}$ valence band at $\overline{\mathrm{K}}$ was $\Delta_\mathrm{SO}=0.469\pm0.008$ eV by angle-resolved photoemission spectroscopy (ARPES). Synchrotron-based grazing-incidence in-plane X-ray diffraction (XRD) revealed the in-plane lattice constant of monolayer WSe$_{2}$ to be $a_\mathrm{WSe_2}=3.2757\pm0.0008 \mathrm{\r{A}}$. This indicates a lattice compression of -0.19 % from bulk WSe$_{2}$. By using experimentally determined graphene lattice constant ($a_\mathrm{MLG}=2.4575\pm0.0007 \mathrm{\r{A}}$), we found that a 3$\times$3 unit cell of the slightly compressed WSe$_{2}$ is perfectly commensurate with a 4$\times$4 graphene lattice with a mismatch below 0.03 %, which could explain why the monolayer WSe$_{2}$ is compressed on MLG. From XRD and first-principles calculations, however, we conclude that the observed size of strain is negligibly small to account for a discrepancy in $\Delta_\mathrm{SO}$ found between exfoliated and epitaxial monolayers in earlier ARPES. In addition, angle-resolved, ultraviolet and X-ray photoelectron spectroscopy shed light on the band alignment between WSe$_{2}$ and MLG/SiC and indicate electron transfer from graphene to the WSe$_{2}$ monolayer. As further revealed by atomic force microscopy, the WSe$_{2}$ island size depends on the number of carbon layers on top of the SiC substrate. This suggests that the epitaxy of WSe$_{2}$ favors the weak van der Waals interactions with graphene while it is perturbed by the influence of the SiC substrate and its carbon buffer layer.


title: Optimally Resilient Strategies in Pushdown Safety Games

abstract: Infinite-duration games with disturbances extend the classical framework of infinite-duration games, which captures the reactive synthesis problem, with a discrete measure of resilience against non-antagonistic external influence. This concerns events where the observed system behavior differs from the intended one prescribed by the controller. For games played on finite arenas it is known that computing optimally resilient strategies only incurs a polynomial overhead over solving classical games. This paper studies safety games with disturbances played on infinite arenas induced by pushdown systems. We show how to compute optimally resilient strategies in triply-exponential time. For the subclass of safety games played on one-counter configuration graphs, we show that determining the degree of resilience of the initial configuration is PSPACE-complete and that optimally resilient strategies can be computed in doubly-exponential time.


title: Asymptotic performance of the Grimmett-McDiarmid heuristic

abstract: Grimmett and McDiarmid suggested a simple heuristic for finding stable sets in random graphs. They showed that the heuristic finds a stable set of size $\sim\log_2 n$ (with high probability) on a $G(n, 1/2)$ random graph. We determine the asymptotic distribution of the size of the stable set found by the algorithm.


title: Manifestly Lorentz invariant chiral boson action

abstract: A manifestly Lorentz invariant action is found for the Floreanini-Jackiw chiral boson. The method involves a novel chiral reduction of the phase-space action for a string, and can be adapted to describe chiral bosons on the heterotic string worldsheet. A similar manifestly Lorentz invariant action is found for an entire class of conformal chiral 2k-form electrodynamics in (4k+2) dimensions which includes the Floreanini-Jackiw theory as the k=0$ case.


title: Voluntary Disclosure and Personalized Pricing

abstract: Central to privacy concerns is that firms may use consumer data to price discriminate. A common policy response is that consumers should be given control over which firms access their data and how. Since firms learn about a consumer's preferences based on the data seen and the consumer's disclosure choices, the equilibrium implications of consumer control are unclear. We study whether such measures improve consumer welfare in monopolistic and competitive markets. We find that consumer control can improve consumer welfare relative to both perfect price discrimination and no personalized pricing. First, consumers can use disclosure to amplify competitive forces. Second, consumers can disclose information to induce even a monopolist to lower prices. Whether consumer control improves welfare depends on the disclosure technology and market competitiveness. Simple disclosure technologies suffice in competitive markets. When facing a monopolist, a consumer needs partial disclosure possibilities to obtain any welfare gains.


title: Context-Aware Dynamic Feature Extraction for 3D Object Detection in Point Clouds

abstract: Varying density of point clouds increases the difficulty of 3D detection. In this paper, we present a context-aware dynamic network (CADNet) to capture the variance of density by considering both point context and semantic context. Point-level contexts are generated from original point clouds to enlarge the effective receptive filed. They are extracted around the voxelized pillars based on our extended voxelization method and processed with the context encoder in parallel with the pillar features. With a large perception range, we are able to capture the variance of features for potential objects and generate attentive spatial guidance to help adjust the strengths for different regions. In the region proposal network, considering the limited representation ability of traditional convolution where same kernels are shared among different samples and positions, we propose a decomposable dynamic convolutional layer to adapt to the variance of input features by learning from local semantic context. It adaptively generates the position-dependent coefficients for multiple fixed kernels and combines them to convolve with local feature windows. Based on our dynamic convolution, we design a dual-path convolution block to further improve the representation ability. We conduct experiments with our Network on KITTI dataset and achieve good performance on 3D detection task for both precision and speed. Our one-stage detector outperforms SECOND and PointPillars by a large margin and achieves the speed of 30 FPS.


title: Search for a narrow resonance lighter than 200 GeV decaying to a pair of muons in proton-proton collisions at $\sqrt{s} =$ 13 TeV

abstract: A search is presented for a narrow resonance decaying to a pair of oppositely charged muons using $\sqrt{s} =$ 13 TeV proton-proton collision data recorded at the LHC. In the 45-75 and 110-200 GeV resonance mass ranges, the search is based on conventional triggering and event reconstruction techniques. In the 11.5-45 GeV mass range, the search uses data collected with dimuon triggers with low transverse momentum thresholds, recorded at high rate by storing a reduced amount of trigger-level information. The data correspond to integrated luminosities of 137 fb$^{-1}$ and 96.6 fb$^{-1}$ for conventional and high-rate triggering, respectively. No significant resonant peaks are observed in the probed mass ranges. The search sets the most stringent constraints to date on a dark photon in the $\sim$30-75 and 110-200 GeV mass ranges.


title: Some results about Mittag-Leffler function's integral representations

abstract: In many articles on the integral expressions of Mittag-Leffler functions, we have found that whether the integral expression can be used at the origin is still unresolved. In this article we give the applicable conditions and proof. And we also generalize some important conclusions about Mittag-Leffler function.


title: GeBioToolkit: Automatic Extraction of Gender-Balanced Multilingual Corpus of Wikipedia Biographies

abstract: We introduce GeBioToolkit, a tool for extracting multilingual parallel corpora at sentence level, with document and gender information from Wikipedia biographies. Despite thegender inequalitiespresent in Wikipedia, the toolkit has been designed to extract corpus balanced in gender. While our toolkit is customizable to any number of languages (and different domains), in this work we present a corpus of 2,000 sentences in English, Spanish and Catalan, which has been post-edited by native speakers to become a high-quality dataset for machinetranslation evaluation. While GeBioCorpus aims at being one of the first non-synthetic gender-balanced test datasets, GeBioToolkit aims at paving the path to standardize procedures to produce gender-balanced datasets


title: Wrapped Branes and Punctured Horizons

abstract: Large classes of AdS$_p$ supergravity backgrounds describing the IR dynamics of $p$-branes wrapped on a Riemann surface are determined by a solution to the Liouville equation. The regular solutions of this equation lead to the well-known wrapped brane supergravity solutions associated with the constant curvature metric on a compact Riemann surface. We show that some singular solutions of the Liouville equation have a physical interpretation as explicit point-like brane sources on the Riemann surface. We uncover the details of this picture by focusing on $\mathcal{N}=1$ theories of class $\mathcal{S}$ arising from M5-branes on a punctured Riemann surface. We present explicit AdS$_5$ solutions dual to these SCFTs and check the holographic duality by showing the non-trivial agreement of 't Hooft anomalies.


title: Testing Smart Contracts Gets Smarter

abstract: Smart contracts are immutable, verifiable, and autonomous pieces of code that can be deployed and ran on blockchain networks like Ethereum. Due to the immutability nature of blockchain, no change is possible on a deployed smart contract or a verified transaction. On the other hand, there are millions of dollars carried by smart contracts in Ethereum blockchain, and hence, a faulty smart contract can lead to a huge monetary loss. Therefore, it is important for smart contract developers to fully test and check the correctness of their code before deploying it on the blockchain. In this paper, we propose a testing mechanism for smart contracts in Solidity language, based on mutation testing. We analyzed a comprehensive list of known bugs in Solidity smart contracts, and designed 10 classes of mutation operators inspired by the real faults. Our experimental results show that our proposed mutation operators can regenerate 10 of 15 famous faulty smart contracts, which have resulted in millions of dollars loss. The results show the effectiveness of our proposed mutation operators in detecting real faults in Solidity smart contracts. We have also extended {\em Universal Mutator } tool with our mutation operators, so that it can automatically generate mutants for smart contracts written in Solidity.


title: Multiple character of non-monotonic size-dependence for relaxation dynamics in polymer-particle and binary mixtures

abstract: Adding plasticizers is a well-known procedure to reduce the glass transition temperature in polymers. It has been recently shown that this effect shows a non-monotonic dependence on the size of additive molecules [The Journal of Chemical Physics 150 (2019) 024903]. In this work, we demonstrate that, as the size of the additive molecules is changed at fixed concentration, multiple extrema emerge in the dependence of the system's relaxation time on the size ratio. The effect occurs on all relevant length scales including single monomer dynamics, decay of Rouse modes and relaxation of the chain's end-to-end vector. A qualitatively similar trend is found within mode-coupling theoretical results for a binary hard-sphere (HS) mixture. An interpretation of the effect in terms of local packing efficiency and coupling between the dynamics of minority and majority species is provided.


title: On-demand generation of entangled photon pairs in the telecom C-band for fiber-based quantum networks

abstract: On-demand sources of entangled photons for the transmission of quantum information in the telecom C-band are required to realize fiber-based quantum networks. So far, non-deterministic sources of quantum states of light were used for long distance entanglement distribution in this lowest loss wavelength range. However, they are fundamentally limited in either efficiency or security due to their Poissonian emission statistics. Here, we show on-demand generation of entangled photon pairs in the telecom C-band by an InAs/GaAs semiconductor quantum dot. Using a robust phonon-assisted excitation scheme we measure a concurrence of $91.4\,\%$ and a fidelity of $95.2\,\%$ to $\Phi^+$. On-demand generation of polarization entangled photons will enable secure quantum communication in fiber-based networks. Furthermore, applying this excitation scheme to several remote quantum dots tuned into resonance will enable first on-demand entanglement distribution over large distances for scalable real-life quantum applications.


title: Frivolous Units: Wider Networks Are Not Really That Wide

abstract: A remarkable characteristic of overparameterized deep neural networks (DNNs) is that their accuracy does not degrade when the network's width is increased. Recent evidence suggests that developing compressible representations is key for adjusting the complexity of large networks to the learning task at hand. However, these compressible representations are poorly understood. A promising strand of research inspired from biology is understanding representations at the unit level as it offers a more granular and intuitive interpretation of the neural mechanisms. In order to better understand what facilitates increases in width without decreases in accuracy, we ask: Are there mechanisms at the unit level by which networks control their effective complexity as their width is increased? If so, how do these depend on the architecture, dataset, and training parameters? We identify two distinct types of "frivolous" units that proliferate when the network's width is increased: prunable units which can be dropped out of the network without significant change to the output and redundant units whose activities can be expressed as a linear combination of others. These units imply complexity constraints as the function the network represents could be expressed by a network without them. We also identify how the development of these units can be influenced by architecture and a number of training factors. Together, these results help to explain why the accuracy of DNNs does not degrade when width is increased and highlight the importance of frivolous units toward understanding implicit regularization in DNNs.


title: A Novel Topology for End-to-end Temporal Classification and Segmentation with Recurrent Neural Network

abstract: Connectionist temporal classification (CTC) has matured as an alignment free to sequence transduction and shows competitive for end-to-end speech recognition. In the CTC topology, the blank symbol occupies more than half of the state trellis, which results the spike phenomenon of the non-blank symbols. For classification task, the spikes work quite well, but as to the segmentation task it does not provide boundaries information. In this paper, a novel topology is introduced to combine the temporal classification and segmentation ability in one framework.


title: On the Convex Properties of Wireless Power Transfer with Nonlinear Energy Harvesting

abstract: The convex property of a nonlinear wireless power transfer (WPT) is characterized in this work. Following a nonlinear energy harvesting model, we express the relationship between the harvested direct current (DC) power and the power of the received radio-frequency signal via an implicit function, based on which the convex property is further proved. In particular, for a predefined rectifier's input signal distribution, we show that the harvested DC power of the nonlinear model is convex in the reciprocal of the rectifier's input signal power. Finally, we provide an example to show the advantages of applying the convex property in WPT network designs.


title: edBB: Biometrics and Behavior for Assessing Remote Education

abstract: We present a platform for student monitoring in remote education consisting of a collection of sensors and software that capture biometric and behavioral data. We define a collection of tasks to acquire behavioral data that can be useful for facing the existing challenges in automatic student monitoring during remote evaluation. Additionally, we release an initial database including data from 20 different users completing these tasks with a set of basic sensors: webcam, microphone, mouse, and keyboard; and also from more advanced sensors: NIR camera, smartwatch, additional RGB cameras, and an EEG band. Information from the computer (e.g. system logs, MAC, IP, or web browsing history) is also stored. During each acquisition session each user completed three different types of tasks generating data of different nature: mouse and keystroke dynamics, face data, and audio data among others. The tasks have been designed with two main goals in mind: i) analyse the capacity of such biometric and behavioral data for detecting anomalies during remote evaluation, and ii) study the capability of these data, i.e. EEG, ECG, or NIR video, for estimating other information about the users such as their attention level, the presence of stress, or their pulse rate.


title: Deformations of representations of fundamental groups of complex varieties

abstract: We describe locally the representation varieties of fundamental groups for smooth complex varieties at representations coming from the monodromy of a variation of mixed Hodge structure. Given such a manifold $X$ and such a linear representation $\rho$ of its fundamental group $\pi_1(X,x)$, we use the theory of Goldman-Millson and pursue our previous work that combines mixed Hodge theory with derived deformation theory to construct a mixed Hodge structure on the formal local ring $\widehat{\mathcal{O}}_\rho$ to the representation variety of $\pi_1(X,x)$ at $\rho$. Then we show how a weighted-homogeneous presentation of $\widehat{\mathcal{O}}_\rho$ is induced directly from a splitting of the weight filtration of its mixed Hodge structure. In this way we recover and generalize theorems of Eyssidieux-Simpson ($X$ compact) and of Kapovich-Millson ($\rho$ finite).


title: The trace of the local $\mathbf{A}^1$-degree

abstract: We prove that the local $\mathbb{A}^1$-degree of a polynomial function at an isolated zero with finite separable residue field is given by the trace of the local $\mathbb{A}^1$-degree over the residue field. This fact was originally suggested by Morel's work on motivic transfers and by Kass and Wickelgren's work on the Scheja-Storch bilinear form. As a corollary, we generalize a result of Kass and Wickelgren's relating the Scheja-Storch form and the local $\mathbb{A}^1$-degree.


title: Feedback control of surface roughness in a one-dimensional KPZ growth process

abstract: Control of generically scale-invariant systems, i.e., targeting specific cooperative features in non-linear stochastic interacting systems with many degrees of freedom subject to strong fluctuations and correlations that are characterized by power laws, remains an important open problem. We study the control of surface roughness during a growth process described by the Kardar--Parisi--Zhang (KPZ) equation in $(1+1)$ dimensions. We achieve the saturation of the mean surface roughness to a prescribed value using non-linear feedback control. Numerical integration is performed by means of the pseudospectral method, and the results are used to investigate the coupling effects of controlled (linear) and uncontrolled (non-linear) KPZ dynamics during the control process. While the intermediate time kinetics is governed by KPZ scaling, at later times a linear regime prevails, namely the relaxation towards the desired surface roughness. The temporal crossover region between these two distinct regimes displays intriguing scaling behavior that is characterized by non-trivial exponents and involves the number of controlled Fourier modes. Due to the control, the height probability distribution becomes negatively skewed, which affects the value of the saturation width.


title: Space-time waves from a collapse with a time dependent cosmological parameter

abstract: We study the emission of space-time waves produced by back-reaction effects during a collapse of a spherically symmetric universe with a time dependent cosmological parameter, which is driven by a scalar field. As in a previous work the final state avoids the final singularity due to the fact the co-moving relativistic observer never reaches the center, because the physical time evolution $d\tau=U_{0}\,dx^0$, decelerates for a co-moving observer which falls with the collapse. The equation of state of the system depends on the rate of the collapse, but always is positive: $0 < \omega(p) < 0.25$.


title: Efficient analytical implementation of the DOT Riemann solver for the de Saint Venant-Exner morphodynamic model

abstract: Within the framework of the de Saint Venant equations coupled with the Exner equation for morphodynamic evolution, this work presents a new efficient implementation of the Dumbser-Osher-Toro (DOT) scheme for non-conservative problems. The DOT path-conservative scheme is a robust upwind method based on a complete Riemann solver, but it has the drawback of requiring expensive numerical computations. Indeed, to compute the non-linear time evolution in each time step, the DOT scheme requires numerical computation of the flux matrix eigenstructure (the totality of eigenvalues and eigenvectors) several times at each cell edge. In this work, an analytical and compact formulation of the eigenstructure for the de Saint Venant-Exner (dSVE) model is introduced and tested in terms of numerical efficiency and stability. Using the original DOT and PRICE-C (a very efficient FORCE-type method) as reference methods, we present a convergence analysis (error against CPU time) to study the performance of the DOT method with our new analytical implementation of eigenstructure calculations (A-DOT). In particular, the numerical performance of the three methods is tested in three test cases: a movable bed Riemann problem with an analytical solution; a problem with smooth analytical solution; a test in which the water flow is characterised by subcritical and supercritical regions. For a given target error, the A-DOT method is always the most efficient choice. Finally, two experimental data sets and different transport formulae are considered to test the A-DOT model in more practical case studies.


title: Training Provably Robust Models by Polyhedral Envelope Regularization

abstract: Training certifiable neural networks enables one to obtain models with robustness guarantees against adversarial attacks. In this work, we introduce a framework to bound the adversary-free region in the neighborhood of the input data by a polyhedral envelope, which yields finer-grained certified robustness. We further introduce polyhedral envelope regularization (PER) to encourage larger polyhedral envelopes and thus improve the provable robustness of the models. We demonstrate the flexibility and effectiveness of our framework on standard benchmarks; it applies to networks of different architectures and general activation functions. Compared with the state-of-the-art methods, PER has very little computational overhead and better robustness guarantees without over-regularizing the model.